  In the 21st century, several emergent viruses have posed a global threat.
Each pathogen has emphasized the value of rapid and scalable vaccine
development programs. The ongoing SARS-CoV-2 pandemic has made the importance
of such efforts especially clear. New biotechnological advances in vaccinology
allow for recent advances that provide only the nucleic acid building blocks of
an antigen, eliminating many safety concerns. During the COVID-19 pandemic,
these DNA and RNA vaccines have facilitated the development and deployment of
vaccines at an unprecedented pace. This success was attributable at least in
part to broader shifts in scientific research relative to prior epidemics; the
genome of SARS-CoV-2 was available as early as January 2020, facilitating
global efforts in the development of DNA and RNA vaccines within two weeks of
the international community becoming aware of the new viral threat.
Additionally, these technologies that were previously only theoretical are not
only safe but also highly efficacious. Although historically a slow process,
the rapid development of vaccines during the COVID-19 crisis reveals a major
shift in vaccine technologies. Here, we provide historical context for the
emergence of these paradigm-shifting vaccines. We describe several DNA and RNA
vaccines and in terms of their efficacy, safety, and approval status. We also
discuss patterns in worldwide distribution. The advances made since early 2020
provide an exceptional illustration of how rapidly vaccine development
technology has advanced in the last two decades in particular and suggest a new
era in vaccines against emerging pathogens.


  Models of counts-of-counts data have been extensively used in the biological
sciences, for example in cancer, population genetics, sampling theory and
ecology. In this paper we explore properties of one model that is embedded into
a continuous-time process and can describe the appearance of certain biological
data such as covid DNA sequences in a database. More specifically, we consider
an evolving model of counts-of-counts data that arises as the family size
counts of samples taken sequentially from a Birth process with Immigration
(BI). Here, each family represents a type or species, and the family size
counts represent the type or species frequency spectrum in the population. We
study the correlation of $S(a,b)$ and $S(c,d)$, the number of families observed
in two disjoint time intervals $(a,b)$ and $(c,d)$. We find the expected sample
variance and its asymptotics for $p$ consecutive sequential samples
$\mathbf{S}_p:=(S(t_0,t_1),\dots, S(t_{p-1},t_p))$, for any given
$0=t_0<t_1<\dots<t_p$. By conditioning on the sizes of the samples, we provide
a connection between $\mathbf{S}_p$ and $p$ sequential samples of sizes
$n_1,n_2,\dots,n_p$, drawn from a single run of a Chinese Restaurant Process.
The properties of the latter were studied in da Silva et al. (2022). We show
how the continuous-time framework helps to make asymptotic calculations easier
than its discrete-time counterpart. As an application, for a specific choice of
$t_1,t_2,\dots, t_p$, we revisit Fisher's 1943 multi-sampling problem and give
another explanation of what Fisher's model could have meant in the world of
sequential samples drawn from a BI process.


  The crisis caused by the COVID-19 outbreak around the globe raised an
increasing concern about the ongoing emergence of variants of SARS-CoV-2 that
may evade the immune response provided by vaccines. New variants appear due to
mutation, and as the cases accumulate, the probability of the emergence of a
variant of concern increases. In this article, we propose a modified SIR model
with waning immunity that captures the competition of two strain classes of an
infectious disease under the effect of vaccination with a highly contagious and
deadly strain class emerging from a prior strain due to mutation. When these
strains compete for a limited supply of susceptible individuals, changes in the
efficiency of vaccines may affect the behaviour of the disease in a non-trivial
way, resulting in complex outcomes. We characterise the parameter space
including intrinsic parameters of the disease, and using the vaccine
efficiencies as control variables. We find different types of transcritical
bifurcations between endemic fixed points and a disease-free equilibrium and
identify a region of strain competition where the two strain classes coexist
during a transient period. We show that a strain can be extinguished either due
to strain competition or vaccination, and we obtain the critical values of the
efficiency of vaccines to eradicate the disease. Numerical studies using
parameters estimated from publicly reported data agree with our theoretical
results. Our mathematical model could be a tool to assess quantitatively the
vaccination policies of competing and emerging strains using the dynamics in
epidemics of infectious diseases.


  Hypergraph neural networks can model multi-way connections among nodes of the
graphs, which are common in real-world applications such as genetic medicine.
In particular, genetic pathways or gene sets encode molecular functions driven
by multiple genes, naturally represented as hyperedges. Thus, hypergraph-guided
embedding can capture functional relations in learned representations. Existing
hypergraph neural network models often focus on node-level or graph-level
inference. There is an unmet need in learning powerful representations of
subgraphs of hypergraphs in real-world applications. For example, a cancer
patient can be viewed as a subgraph of genes harboring mutations in the
patient, while all the genes are connected by hyperedges that correspond to
pathways representing specific molecular functions. For accurate inductive
subgraph prediction, we propose SubHypergraph Inductive Neural nEtwork (SHINE).
SHINE uses informative genetic pathways that encode molecular functions as
hyperedges to connect genes as nodes. SHINE jointly optimizes the objectives of
end-to-end subgraph classification and hypergraph nodes' similarity
regularization. SHINE simultaneously learns representations for both genes and
pathways using strongly dual attention message passing. The learned
representations are aggregated via a subgraph attention layer and used to train
a multilayer perceptron for inductive subgraph inferencing. We evaluated SHINE
against a wide array of state-of-the-art (hyper)graph neural networks, XGBoost,
NMF and polygenic risk score models, using large scale NGS and curated
datasets. SHINE outperformed all comparison models significantly, and yielded
interpretable disease models with functional insights.


  The Susceptible-Infected (SI) and Susceptible-Infected-Recovered (SIR) models
provide two distinct representations of epidemic evolution, distinguished by
the lack of spontaneous recovery in the SI model. Here we introduce a new
active matter epidemic model, the ``Susceptible-Cleric-Zombie-Recovered''
(SCZR) model, in which spontaneous recovery is absent but zombies can recover
with probability $\gamma$ via interaction with a cleric. Upon interacting with
a zombie, both susceptibles and clerics can enter the zombie state with
probability $\beta$ and $\alpha$, respectively. By changing the intial fraction
of clerics or their healing ability rate $\gamma$, we can tune the SCZR model
between SI dynamics, in which no susceptibles or clerics remain at long times,
and SIR dynamics, in which no zombies remain at long times. The model is
relevant to certain real world diseases such as HIV where spontaneous recovery
is impossible but where medical interventions by a limited number of caregivers
can reduce or eliminate the spread of infection.


  High throughput screening of small molecules and natural products is costly,
requiring significant amounts of time, reagents, and operating space. Although
microarrays have proven effective in the miniaturization of screening for
certain biochemical assays, such as nucleic acid hybridization or antibody
binding, they are not widely used for drug discovery in cell culture due to the
need for cells to internalize lipophilic drug candidates. Lipid droplet
microarrays are a promising solution to this problem as they are capable of
delivering lipophilic drugs to cells at dosages comparable to solution
delivery. However, the scalablility of the array fabrication, assay validation,
and screening steps has limited the utility of this approach. Here we
demonstrate a scalable process for lipid droplet array fabrication, assay
validation in cell culture, and drug screening. A nanointaglio printing process
has been adapted for use with a printing press. The arrays are stabilized for
immersion into aqueous solution using a vapor coating process. In addition to
delivery of lipophilic compounds, we found that we are also able to encapsulate
and deliver a water-soluble compound in this way. The arrays can be
functionalized by extracellular matrix proteins such as collagen prior to cell
culture as the mechanism for uptake is based on direct contact with the lipid
delivery vehicles rather than diffusion of the drug out of the microarray
spots. We demonstrate this method for delivery to 3 different cell types and
the screening of 90 natural product extracts on a microarray covering an area
of less than 0.1 cm2. The arrays are suitable for miniaturized screening, for
instance in BSL-3 conditions where space is limited and for applications where
cell numbers are limited, such as in functional precision medicine.


  Cryo-electron microscopy (cryo-EM) is an imaging modality that provides
unique insights into the dynamics of proteins and other building blocks of
life. The algorithmic challenge of jointly estimating the poses, 3D structure,
and conformational heterogeneity of a biomolecule from millions of noisy and
randomly oriented 2D projections in a computationally efficient manner,
however, remains unsolved. Our method, cryoFIRE, performs ab initio
heterogeneous reconstruction with unknown poses in an amortized framework,
thereby avoiding the computationally expensive step of pose search while
enabling the analysis of conformational heterogeneity. Poses and conformation
are jointly estimated by an encoder while a physics-based decoder aggregates
the images into an implicit neural representation of the conformational space.
We show that our method can provide one order of magnitude speedup on datasets
containing millions of images without any loss of accuracy. We validate that
the joint estimation of poses and conformations can be amortized over the size
of the dataset. For the first time, we prove that an amortized method can
extract interpretable dynamic information from experimental datasets.


  Rheumatoid arthritis (RA) is an autoimmune condition caused when patients'
immune system mistakenly targets their own tissue. Machine learning (ML) has
the potential to identify patterns in patient electronic health records (EHR)
to forecast the best clinical treatment to improve patient outcomes. This study
introduced a Drug Response Prediction (DRP) framework with two main goals: 1)
design a data processing pipeline to extract information from tabular clinical
data, and then preprocess it for functional use, and 2) predict RA patient's
responses to drugs and evaluate classification models' performance. We propose
a novel two-stage ML framework based on European Alliance of Associations for
Rheumatology (EULAR) criteria cutoffs to model drug effectiveness. Our model
Stacked-Ensemble DRP was developed and cross-validated using data from 425 RA
patients. The evaluation used a subset of 124 patients (30%) from the same data
source. In the evaluation of the test set, two-stage DRP leads to improved
classification accuracy over other end-to-end classification models for binary
classification. Our proposed method provides a complete pipeline to predict
disease activity scores and identify the group that does not respond well to
anti-TNF treatments, thus showing promise in supporting clinical decisions
based on EHR information.


  We lay the foundation of a circuit theory for chemical reaction networks.
Chemical reactions are grouped into chemical modules solely characterized by
their current-concentration characteristic, as electrical devices by their
current-voltage (I-V) curve in electronic circuit theory. This, combined with
the chemical analog of Kirchhoff's current and voltage laws, provides a
powerful tool to predict reaction currents and dissipation across complex
chemical networks. The theory can serve to build accurate reduced models of
complex networks as well as to design networks performing desired tasks.


  Motivation: Protein-to-genome alignment is critical to annotating genes in
non-model organisms. While there are a few tools for this purpose, all of them
were developed over ten years ago and did not incorporate the latest advances
in alignment algorithms. They are inefficient and could not keep up with the
rapid production of new genomes and quickly growing protein databases.
  Results: Here we describe miniprot, a new aligner for mapping protein
sequences to a complete genome. Miniprot integrates recent techniques such as
k-mer sketch and SIMD-based dynamic programming. It is tens of times faster
than existing tools while achieving comparable accuracy on real data.
  Availability and implementation: https://github.com/lh3/miniprot


  In protein biophysics, the separation between the functionally important
residues (forming the active site or binding surface) and those that create the
overall structure (the fold) is a well-established and fundamental concept.
Identifying and modifying those functional sites is critical for protein
engineering but computationally non-trivial, and requires significant domain
knowledge. To automate this process from a data-driven perspective, we propose
a disentangled Wasserstein autoencoder with an auxiliary classifier, which
isolates the function-related patterns from the rest with theoretical
guarantees. This enables one-pass protein sequence editing and improves the
understanding of the resulting sequences and editing actions involved. To
demonstrate its effectiveness, we apply it to T-cell receptors (TCRs), a
well-studied structure-function case. We show that our method can be used to
alter the function of TCRs without changing the structural backbone,
outperforming several competing methods in generation quality and efficiency,
and requiring only 10% of the running time needed by baseline models. To our
knowledge, this is the first approach that utilizes disentangled
representations for TCR engineering.


  Inference of gene regulatory networks has been an active area of research for
around 20 years, leading to the development of sophisticated inference
algorithms based on a variety of assumptions and approaches. With the always
increasing demand for more accurate and powerful models, the inference problem
remains of broad scientific interest. The abstract representation of biological
systems through gene regulatory networks represents a powerful method to study
such systems, encoding different amounts and types of information. In this
review, we summarize the different types of inference algorithms specifically
based on time-series transcriptomics, giving an overview of the main
applications of gene regulatory networks in computational biology. This review
is intended to give an updated overview of regulatory networks inference tools
to biologists and researchers new to the topic and guide them in selecting the
appropriate inference method that best fits their questions, aims and
experimental data.


  In the scope of drug discovery, the molecular design aims to identify novel
compounds from the chemical space where the potential drug-like molecules are
estimated to be in the order of 10^60 - 10^100. Since this search task is
computationally intractable due to the unbounded search space, deep learning
draws a lot of attention as a new way of generating unseen molecules. As we
seek compounds with specific target proteins, we propose a Transformer-based
deep model for de novo target-specific molecular design. The proposed method is
capable of generating both drug-like compounds (without specified targets) and
target-specific compounds. The latter are generated by enforcing different keys
and values of the multi-head attention for each target. In this way, we allow
the generation of SMILES strings to be conditional on the specified target.
Experimental results demonstrate that our method is capable of generating both
valid drug-like compounds and target-specific compounds. Moreover, the sampled
compounds from conditional model largely occupy the real target-specific
molecules' chemical space and also cover a significant fraction of novel
compounds.


  Proteins are macromolecules that perform essential functions in all living
organisms. Designing novel proteins with specific structures and desired
functions has been a long-standing challenge in the field of bioengineering.
Existing approaches generate both protein sequence and structure using either
autoregressive models or diffusion models, both of which suffer from high
inference costs. In this paper, we propose a new approach capable of protein
sequence and structure co-design, which iteratively translates both protein
sequence and structure into the desired state from random initialization, based
on context features given a priori. Our model consists of a trigonometry-aware
encoder that reasons geometrical constraints and interactions from context
features, and a roto-translation equivariant decoder that translates protein
sequence and structure interdependently. Notably, all protein amino acids are
updated in one shot in each translation step, which significantly accelerates
the inference process. Experimental results across multiple tasks show that our
model outperforms previous state-of-the-art baselines by a large margin, and is
able to design proteins of high fidelity as regards both sequence and
structure, with running time orders of magnitude less than sampling-based
methods.


  In their recent article, Madej et al. 1 proposed an original way to solve the
recurrent issue of controlling for the false discovery rate (FDR) in
peptide-spectrum-match (PSM) validation. Briefly, they proposed to derive a
single precise distribution of decoy matches termed the Common Decoy
Distribution (CDD) and to use it to control for FDR during a target-only
search. Conceptually, this approach is appealing as it takes the best of two
worlds, i.e., decoy-based approaches (which leverage a large-scale collection
of empirical mismatches) and decoy-free approaches (which are not subject to
the randomness of decoy generation while sparing an additional database
search). Interestingly, CDD also corresponds to a middle-of-the-road approach
in statistics with respect to the two main families of FDR control procedures:
Although historically based on estimating the falsepositive distribution, FDR
control has recently been demonstrated to be possible thanks to competition
between the original variables (in proteomics, target sequences) and their
fictional counterparts (in proteomics, decoys). Discriminating between these
two theoretical trends is of prime importance for computational proteomics. In
addition to highlighting why proteomics was a source of inspiration for
theoretical biostatistics, it provides practical insights into the improvements
that can be made to FDR control methods used in proteomics, including CDD.


  Repeat asymptomatic testing in order to identify and quarantine infectious
individuals has become a widely-used intervention to control SARS-CoV-2
transmission. In some workplaces, and in particular health and social care
settings with vulnerable patients, regular asymptomatic testing has been
deployed to staff to reduce the likelihood of workplace outbreaks. We have
developed a model based on data available in the literature to predict the
potential impact of repeat asymptomatic testing on SARS-CoV-2 transmission. The
results highlight features that are important to consider when modelling
testing interventions, including population heterogeneity of infectiousness and
correlation with test-positive probability, as well as adherence behaviours in
response to policy. Furthermore, the model based on the reduction in
transmission potential presented here can be used to parameterise existing
epidemiological models without them having to explicitly simulate the testing
process. Overall, we find that even with different model paramterisations, in
theory, regular asymptomatic testing is likely to be a highly effective measure
to reduce transmission in workplaces, subject to adherence.


  Oscillatory behaviors are ubiquitous in nature and the human society.
However, most previous works fail to reproduce them in the two-strategy
game-theoretical models. Here we show that oscillatory behaviors naturally
emerge if incomplete information is incorporated into the cooperation evolution
of a non-Markov model. Specifically, we consider a population playing
prisoner's dilemma game, where each individual can only probabilistically get
access to their neighbors' payoff information and store them within their
memory with a given length. They make their decisions based upon these
memories. Interestingly, we find that the level of cooperation generally cannot
stabilize but render quasi-periodic oscillation, and this observation is
strengthened for a longer memory and a smaller information acquisition
probability. The mechanism uncovered shows that there are misperceived payoffs
about the player's neighborhood, facilitating the growth of cooperators and
defectors at different stages that leads to oscillatory behaviors as a result.
Our findings are robust to the underlying structure of the population. Given
the omnipresence of incomplete information, our findings may provide a
plausible explanation for the phenomenon of oscillatory behaviors in the real
world.


  Direct reciprocity facilitates the evolution of cooperation when individuals
interact repeatedly. Most previous studies on direct reciprocity implicitly
assume compulsory interactions. Yet, interactions are often voluntary in human
societies. Here, we consider repeated optional games, where individuals can
freely opt out of each interaction and rejoin later. We find that voluntary
participation greatly promotes cooperation in repeated interactions, even in
harsh situations where repeated compulsory games and one-shot optional games
yield low cooperation rates. Moreover, we theoretically characterize all Nash
equilibria that support cooperation among reactive strategies, and identify
three novel classes of strategies that are error-robust, readily become
equilibria, and dominate in the evolutionary dynamics. The success of these
strategies hinges on the effect of opt-out: it not only avoids trapping in
mutual defection but also poses additional threats to intentional defectors.
Our work highlights that voluntary participation is a simple and effective
mechanism to enhance cooperation in repeated interactions.


  It is recognized that social heterogeneities in terms of the contact
distribution have a strong influence on the spread of infectious diseases.
Nevertheless, few data are available on the group composition of social
contacts, and their statistical description does not possess universal patterns
and may vary spatially and temporally. It is therefore essential to design
robust control strategies, mimicking the effects of non-pharmaceutical
interventions, to limit efficiently the number of infected cases. In this work,
starting from a recently introduced kinetic model for epidemiological dynamics
that takes into account the impact of social contacts of individuals, we
consider an uncertain contact formation dynamics leading to slim-tailed as well
as fat-tailed distributions of contacts. Hence, we analyse the effects of an
optimally robust control strategy of the system of agents. Thanks to classical
methods of kinetic theory, we couple uncertainty quantification methods with
the introduced mathematical model to assess the effects of social limitations.
Finally, using the proposed modelling approach and starting from available
data, we show the effectiveness of the proposed selective measures to dampen
uncertainties together with the epidemic trends.


  We propose an interacting particle system with a moving boundary to model the
spread of an epidemic. Each individual in a given population starts from some
level of shielding, given by a non-negative real number, and this level then
evolves over time according to diffusive dynamics driven by independent
Brownian motions. If the level of shielding gets too low, individuals will find
themselves in `at-risk' situations, as captured by proximity to a lower moving
boundary, which we call the `advancing front' of the epidemic. Specifically,
local time is accumulated along this boundary and, while individuals are
initially reflected, infection may eventually occur with a likelihood depending
on the accumulated local time as well as the intrinsic infectiousness of the
disease and the current contagiousness within the population. Our main
technical contribution is to give a rigorous formulation of this model in the
form of a well-posed interacting particle system. Moreover, we exploit the
precise construction of the system to establish an important martingale
property of the infected proportion and present a result on its limiting
behaviour if the size of the population is sent to infinity.


  When complex systems with nonlinear dynamics achieve an output performance
objective, only a fraction of the state dynamics significantly impacts that
output. Those minimal state dynamics can be identified using the differential
geometric approach to the observability of nonlinear systems, but the theory is
limited to only analytical systems. In this paper, we extend the notion of
nonlinear observable decomposition to the more general class of data-informed
systems. We employ Koopman operator theory, which encapsulates nonlinear
dynamics in linear models, allowing us to bridge the gap between linear and
nonlinear observability notions. We propose a new algorithm to learn Koopman
operator representations that capture the system dynamics while ensuring that
the output performance measure is in the span of its observables. We show that
a transformation of this linear, output-inclusive Koopman model renders a new
minimum Koopman representation. This representation embodies only the
observable portion of the nonlinear observable decomposition of the original
system. A prime application of this theory is to identify genes in biological
systems that correspond to specific phenotypes, the performance measure. We
simulate two biological gene networks and demonstrate that the observability of
Koopman operators can successfully identify genes that drive each phenotype. We
anticipate our novel system identification tool will effectively discover
reduced gene networks that drive complex behaviors in biological systems.


  Aligning electron density maps from Cryogenic electron microscopy (cryo-EM)
is a first key step for studying multiple conformations of a biomolecule. As
this step remains costly and challenging, with standard alignment tools being
potentially stuck in local minima, we propose here a new procedure, called
AlignOT, which relies on the use of computational optimal transport (OT) to
align EM maps in 3D space. By embedding a fast estimation of OT maps within a
stochastic gradient descent algorithm, our method searches for a rotation that
minimizes the Wasserstein distance between two maps, represented as point
clouds. We quantify the impact of various parameters on the precision and
accuracy of the alignment, and show that AlignOT can outperform the standard
local alignment methods, with an increased range of rotation angles leading to
proper alignment. We further benchmark AlignOT on various pairs of experimental
maps, which account for different types of conformational heterogeneities and
geometric properties. As our experiments show good performance, we anticipate
that our method can be broadly applied to align 3D EM maps.


  For a biological system to grow, the biomass must be incorporated,
transferred, and accumulated into the underlying reaction network. There are
two perspectives for studying growth dynamics of reaction networks: one way is
to focus on each node in the networks and study its associated influxes and
effluxes. The other way is to focus on a fraction of biomass and study its
trajectory along the reaction pathways. The former perspective (analogous to
the "Eulerian representation" in fluid mechanics) has been studied extensively,
while the latter perspective (analogous to the "Lagrangian representation" in
fluid mechanics) has not been systematically explored. In this work, I
characterized the biomass transfer process for autocatalytic, growing systems
with scalable reaction fluxes. Under balanced growth, the long-term growth
dynamics of the systems are described by delay differential equations (DDEs).
The kernel function of the DDE serves as a unique pattern for the catalytic
delay for a reaction network, and in frequency domain the delay spectrum
provides a geometric interpretation for long-term growth rate. The DDE
formulation provides a clear intuition on how autocatalytic reaction pathways
lead to system growth, it also enables us to classify and compare reaction
networks with different network structures.


  Accurate and reliable forecasting of emerging dominant severe acute
respiratory syndrome coronavirus 2 (SARS-CoV-2) variants enables policymakers
and vaccine makers to get prepared for future waves of infections. The last
three waves of SARS-CoV-2 infections caused by dominant variants Omicron
(BA.1), BA.2, and BA.4/BA.5 were accurately foretold by our artificial
intelligence (AI) models built with biophysics, genotyping of viral genomes,
experimental data, algebraic topology, and deep learning. Based on newly
available experimental data, we analyzed the impacts of all possible viral
spike (S) protein receptor-binding domain (RBD) mutations on the SARS-CoV-2
infectivity. Our analysis sheds light on viral evolutionary mechanisms, i.e.,
natural selection through infectivity strengthening and antibody resistance. We
forecast that BA.2.10.4, BA.2.75, BQ.1.1, and particularly, BA.2.75+R346T, have
high potential to become new dominant variants to drive the next surge.


  Rule-based models have been successfully used to represent different aspects
of the COVID-19 pandemic, including age, testing, hospitalisation, lockdowns,
immunity, infectivity, behaviour, mobility and vaccination of individuals.
These rule-based approaches are motivated by chemical reaction rules which are
traditionally solved numerically with the standard Gillespie algorithm proposed
in the context of molecular dynamics. When applying reaction system type of
approaches to epidemiology, generalisations of the Gillespie algorithm are
required due to the time-dependency of the problems. In this article, we
present different generalisations of the standard Gillespie algorithm which
address discrete subtypes (e.g., incorporating the age structure of the
population), time-discrete updates (e.g., incorporating daily imposed change of
rates for lockdowns) and deterministic delays (e.g., given waiting time until a
specific change in types such as release from isolation occurs). These
algorithms are complemented by relevant examples in the context of the COVID-19
pandemic and numerical results.


  Background: RNA guanine-7 methyltransferase (RNMT) is one of the main
regulators of N7-methylguanosine, and the deregulation of RNMT correlated with
tumor development and immune metabolism. However, the specific function of RNMT
in pan-cancer remains unclear.
  Methods: RNMT expression in different cancers was analyzed using multiple
databases, including Cancer Cell Line Encyclopedia (CCLE), Genotype-Tissue
Expression Project (GTEx), and The Cancer Genome Atlas (TCGA). Cox regression
analysis and Kaplan-Meier analysis were used to estimate the correlation of
RNMT expression to prognosis. The data was also used to research the
relationship between RNMT expression and common immunoregulators, tumor
mutation burden (TMB), microsatellite instability (MSI), mismatch repair (MMR),
and DNA methyltransferase (DNMT). Additionally, the cBioPortal website was used
to evaluate the characteristics of RNMT alteration. The TISDB database was used
to obtain the expression of different subtypes. The Tumor Immune Estimation
Resource (TIMER) database was used to analyze the association between RNMT and
tumor immune infiltration. Gene set enrichment analysis (GSEA) was used to
identify the relevant pathways.
  Results: RNMT was ubiquitously highly expressed across cancers and survival
analysis revealed that its expression was highly associated with the clinical
prognosis of various cancer types. Remarkably, RNMT participates in immune
regulation and plays a crucial part in the tumor microenvironment. A positive
association was found between RNMT expression and six immune cell types
expression in colon adenocarcinoma, kidney renal clear cell carcinoma, and
liver hepatocellular carcinoma. Moreover, RNMT expression was highly associated
with immunoregulators in most cancer types, and correlated to TMB, MSI, MMR,
and DNMT. Finally, GSEA indicated that RNMT may correlate with tumor immunity.


  The genetic code refers to a rule that maps 64 codons to 20 amino acids.
Nearly all organisms, with few exceptions, share the same genetic code, the
standard genetic code (SGC). While it remains unclear why this universal code
has arisen and been maintained during evolution, it may have been preserved
under selection pressure. Theoretical studies comparing the SGC and numerically
created hypothetical random genetic codes have suggested that the SGC has been
subject to strong selection pressure for being robust against translation
errors. However, these prior studies have searched for random genetic codes in
only a small subspace of the possible code space due to limitations in
computation time. Thus, how the genetic code has evolved, and the
characteristics of the genetic code fitness landscape, remain unclear. By
applying multicanonical Monte Carlo, an efficient rare-event sampling method,
we efficiently sampled random codes from a much broader random ensemble of
genetic codes than in previous studies, estimating that only one out of every
$10^{20}$ random codes is more robust than the SGC. This estimate is
significantly smaller than the previous estimate, one in a million. We also
characterized the fitness landscape of the genetic code that has four major
fitness peaks, one of which includes the SGC. Furthermore, genetic algorithm
analysis revealed that evolution under such a multi-peaked fitness landscape
could be strongly biased toward a narrow peak, in an evolutionary
path-dependent manner.


  Background: DNA, RNA, and protein sequence motifs can be recognition sites
for biological functions such as regulation, DNA base modification, and
molecular binding in general. The gain and loss of such motifs can carry
important consequences. When comparing sequences, the analysis of individual
variants does not impart an understanding of the impact on these sites. Rather,
only when these variants considered together with their neighbors and the
original sequence context does this become possible.
  Results: The motif-variant probe (mvp) makes this consideration, counting
instances of user-specified sequence motifs before and after mutation and
reports those that result in motif gain or loss. mvp can perform a similar
analysis for proteins with amino acid variant data. The software is freely
available at https://lpcdrp.gitlab.io/mvp and also installable with the conda
package manager.
  Conclusions: The ability to easily search for variants affecting any motif,
together with the simultaneous consideration of neighboring variants makes mvp
a versatile tool to aid in a less-frequented dimension of comparative genomics.


  We study an epidemic model for a constant population by taking into account
four compartments of the individuals characterizing their states of health.
Each individual is in one of the compartments susceptible (S); incubated -
infected yet not infectious (C), infected and infectious (I), and recovered -
immune (R). An infection is 'visible' only when an individual is in state I.
Upon infection, an individual performs the transition pathway S to C to I to R
to S remaining in each compartments C, I, and R for certain random waiting
times, respectively. The waiting times for each compartment are independent and
drawn from specific probability density functions (PDFs) introducing memory
into the model. We derive memory evolution equations involving convolutions
(time derivatives of general fractional type). We obtain formulae for the
endemic equilibrium and a condition of its existence for cases when the waiting
time PDFs have existing means. We analyze the stability of healthy and endemic
equilibria and derive conditions for which the endemic state becomes
oscillatory (Hopf) unstable. We implement a simple multiple random walker's
approach (microscopic model of Brownian motion of Z independent walkers) with
random SCIRS waiting times into computer simulations. Infections occur with a
certain probability by collisions of walkers in compartments I and S. We
compare the endemic states predicted in the macroscopic model with the
numerical results of the simulations and find accordance of high accuracy. We
conclude that a simple random walker's approach offers an appropriate
microscopic description for the macroscopic model.


  Genome assembly, the process of reconstructing a long genetic sequence by
aligning and merging short fragments, or reads, is known to be NP-hard, either
as a version of the shortest common superstring problem or in a
Hamiltonian-cycle formulation. That is, the computing time is believed to grow
exponentially with the the problem size in the worst case. Despite this fact,
high-throughput technologies and modern algorithms currently allow
bioinformaticians to produce and assemble datasets of billions of reads. Using
methods from statistical mechanics, we address this conundrum by demonstrating
the existence of a phase transition in the computational complexity of the
problem and showing that practical instances always fall in the `easy' phase
(solvable by polynomial-time algorithms). In addition, we propose a
Markov-chain Monte Carlo method that outperforms common deterministic
algorithms in the hard regime.


  Motivated by bacterial chemotaxis and multi-species ecological interactions
in heterogeneous environments, we study a general one-dimensional
reaction-cross-diffusion system in the presence of spatial heterogeneity in
both transport and reaction terms. Under a suitable asymptotic assumption that
the transport is slow over the domain, while gradients in the reaction
heterogeneity are not too sharp, we study the stability of a heterogeneous
steady state approximated by the system in the absence of transport. Using a
WKB ansatz, we find that this steady state can undergo a Turing-type
instability in subsets of the domain, leading to the formation of localized
patterns. The boundaries of the pattern-forming regions are given
asymptotically by `local' Turing conditions corresponding to a spatially
homogeneous analysis parameterized by the spatial variable. We developed a
general open-source code which is freely available, and show numerical examples
of this localized pattern formation in a Schnakenberg cross-diffusion system, a
Keller-Segel chemotaxis model, and the Shigesada-Kawasaki-Teramoto model with
heterogeneous parameters. We numerically show that the patterns may undergo
secondary instabilities leading to spatiotemporal movement of spikes, though
these remain approximately within the asymptotically predicted localized
regions. This theory can elegantly differentiate between spatial structure due
to background heterogeneity, from spatial patterns emergent from Turing-type
instabilities.


  When a dynamical system can be modeled as a sequence of observations, Granger
causality is a powerful approach for detecting predictive interactions between
its variables. However, traditional Granger causal inference has limited
utility in domains where the dynamics need to be represented as directed
acyclic graphs (DAGs) rather than as a linear sequence, such as with cell
differentiation trajectories. Here, we present GrID-Net, a framework based on
graph neural networks with lagged message passing for Granger causal inference
on DAG-structured systems. Our motivating application is the analysis of
single-cell multimodal data to identify genomic loci that mediate the
regulation of specific genes. To our knowledge, GrID-Net is the first
single-cell analysis tool that accounts for the temporal lag between a genomic
locus becoming accessible and its downstream effect on a target gene's
expression. We applied GrID-Net on multimodal single-cell assays that profile
chromatin accessibility (ATAC-seq) and gene expression (RNA-seq) in the same
cell and show that it dramatically outperforms existing methods for inferring
regulatory locus-gene links, achieving up to 71% greater agreement with
independent population genetics-based estimates. By extending Granger causality
to DAG-structured dynamical systems, our work unlocks new domains for causal
analyses and, more specifically, opens a path towards elucidating gene
regulatory interactions relevant to cellular differentiation and complex human
diseases at unprecedented scale and resolution.


  Computing observables from conditioned dynamics is typically computationally
hard, because, although obtaining independent samples efficiently from the
unconditioned dynamics is usually feasible, generally most of the samples must
be discarded (in a form of importance sampling) because they do not satisfy the
imposed conditions. Sampling directly from the conditioned distribution is
non-trivial, as conditioning breaks the causal properties of the dynamics which
ultimately renders the sampling procedure efficient. One standard way of
achieving it is through a Metropolis Monte-Carlo procedure, but this procedure
is normally slow and a very large number of Monte-Carlo steps is needed to
obtain a small number of statistically independent samples. In this work, we
propose an alternative method to produce independent samples from a conditioned
distribution. The method learns the parameters of a generalized dynamical model
that optimally describe the conditioned distribution in a variational sense.
The outcome is an effective, unconditioned, dynamical model, from which one can
trivially obtain independent samples, effectively restoring causality of the
conditioned distribution. The consequences are twofold: on the one hand, it
allows us to efficiently compute observables from the conditioned dynamics by
simply averaging over independent samples. On the other hand, the method gives
an effective unconditioned distribution which is easier to interpret. The
method is flexible and can be applied virtually to any dynamics. We discuss an
important application of the method, namely the problem of epidemic risk
assessment from (imperfect) clinical tests, for a large family of
time-continuous epidemic models endowed with a Gillespie-like sampler. We show
that the method compares favorably against the state of the art, including the
soft-margin approach and mean-field methods.


  As the growth of genomics samples rapidly expands due to increased access to
high resolution DNA sequencing technology, the need for a scalable platform to
aggregate dispersed datasets enable easy access to the vast wealth of DNA
sequences available is paramount. In this work, we introduce and demonstrate a
novel way to use Named Data Networking (NDN) in conjunction with a Kubernetes
cluster to design a flexible and scalable genomics Data Lake in the cloud. In
addition, the use of the NDN Data Plane Development Kit (DPDK) provides an
efficient and accessible distribution of the datasets to researchers anywhere.
This report will explain the need to deploy a Data Lake for genomics data, what
is necessary to deploy successfully, and detailed instructions to replicate the
proposed design. Finally, this technical report outlines future enhancement
options for further improvements.


  Non-consumptive effects such as fear of depredation, can strongly influence
predator-prey dynamics. These effects have not been as well studied in the case
of purely competitive systems, despite ecological and social motivations for
the same. In this work we consider the classic two species ODE and PDE
Lokta-Volterra competition models, where \emph{one} of the competitors is
"fearful" of the other. We find that the presence of fear can have several
interesting dynamical effects on the classical scenarios of weak and strong
competition, and competitive exclusion. Notably, for fear levels in certain
regimes, we show bi-stability between interior equilibrium and boundary
equilibrium is possible - contrary to the classical strong competition
situation where bi-stability is only possible between boundary equilibrium.
Furthermore, in the spatially explicit setting, the effects of several
spatially heterogeneous fear functions are investigated. In particular, we show
that under certain $\mathbb{L}^{1}$ restrictions on the fear function, a weak
competition type situation can change to competitive exclusion. Applications of
these results to ecological as well as sociopolitical settings are discussed,
that connect to the "landscape of fear" (LOF) concept in ecology.


  Eco-evolutionary dynamics is crucial to understand how individuals' behaviors
and the surrounding environment interplay with each other. Typically, it is
assumed that individuals update their behaviors via linear imitation function,
i.e., the replicator dynamics. It has been proved that there cannot be limit
circles in such eco-evolutionary dynamics. It suggests that eco-evolutionary
dynamics alone is not sufficient to explain the widely observed fluctuating
behavior in both the biological and social systems. We extrapolate from the
linear imitation function to general imitation function, which can be
non-linear. It is shown that the general imitation does not change the internal
equilibrium and its local stability. It however leads to limit cycles, which
are never present in classical eco-evolutionary dynamics. Moreover, the average
cooperation level and average environment state agree with the first and the
second component of the internal fixed point. Furthermore, we estimate the
location of the emergent limit cycle. On the one hand, our results provide an
alternative mechanism to make cooperators and defectors coexist in a periodic
way in the "tragedy of the commons" and indicates the global eco-evolutionary
dynamics is sensitive to the imitation function. Our work indicates that the
way of imitation, i.e., the imitation function, is crucial in eco-evolutionary
dynamics, which is not true for evolutionary dynamics with a static
environment.


  The COVID-19 pandemic has affected countries across the world, demanding
drastic public health policies to mitigate the spread of infection, leading to
economic crisis as a collateral damage. In this work, we investigated the
impact of human mobility (described via international commercial flights) on
COVID-19 infection dynamics at the global scale. For this, we developed a graph
neural network-based framework referred to as Dynamic Connectivity GraphSAGE
(DCSAGE), which operates over spatiotemporal graphs and is well-suited for
dynamically changing adjacency information. To obtain insights on the relative
impact of different geographical locations, due to their associated air
traffic, on the evolution of the pandemic, we conducted local sensitivity
analysis on our model through node perturbation experiments. From our analyses,
we identified Western Europe, North America, and Middle East as the leading
geographical locations fueling the pandemic, attributed to the enormity of air
traffic originating or transiting through these regions. We used these
observations to identify tangible air traffic reduction strategies that can
have a high impact on controlling the pandemic, with minimal interference to
human mobility. Our work provides a robust deep learning-based tool to study
global pandemics and is of key relevance to policy makers to take informed
decisions regarding air traffic restrictions during future outbreaks.


  We explored the Protein DataBank (PDB) to collect protein-ssDNA structures
and create a multiconformational docking benchmark including both bound and
unbound protein structures. Due to ssDNA high flexibility when not bound, no
ssDNA unbound structure is included in the benchmark. For the 91
sequence-identity groups identified as bound-unbound structures of the same
protein, we studied the conformational changes in the protein induced by the
ssDNA binding. Moreover, based on several bound or unbound protein structures
in some groups, we also assessed the intrinsic conformational variability in
either bound or unbound conditions, and compared it to the supposedly
binding-induced modifications. To illustrate a use case of this benchmark, we
performed docking experiments using ATTRACT docking software. This benchmark
is, to our knowledge, the first one made to peruse available structures of
ssDNA-protein interactions to such an extent, aiming to improve computational
docking tools dedicated to this kind of molecular interactions.


  This article describes a method to estimate the mortality rate ratio R from
current status data with duration in a chronic condition in case the general
mortality of the overall population is known. Apart from the general mortality,
the method requires four pieces of information from the study participants: age
and time at the survey/interview, whether the chronic condition is present
(current status) and if so, for how long the condition is present (duration).
The method uses a differential equation that relates prevalence, incidence and
mortality to estimate R of the people with the chronic condition compared to
those without the condition.
  To demonstrate feasibility, a simulation based on the illness-death model
(multi-state model) with transition rates motivated from long-term care is run.
It is found that the method requires a large number of study participants
(100000 or more) to estimate R with a reasonably low relative error. Despite
the large sample size, the method can be useful in settings when
cross-sectional information are easily available, e.g., in claims data, and
national age-specific general mortality rates are accessible from vital
statistics.


  In order to grow in any given environment, bacteria need to collect
information about the medium composition and implement suitable growth
strategies by adjusting their regulatory and metabolic degrees of freedom. In
the standard sense, optimal strategy selection is achieved when bacteria grow
at the fastest rate possible in that medium. While this view of optimality is
well suited for cells that have perfect knowledge about their surroundings
(e.g. nutrient levels), things are more involved in uncertain or fluctuating
conditions, especially when changes occur over timescales comparable to (or
faster than) those required to organize a response. Information theory however
provides recipes for how cells can choose the optimal growth strategy under
uncertainty about the stress levels they will face. Here we analyse the
theoretically optimal scenarios for a coarse-grained, experiment-inspired model
of bacterial metabolism for growth in a medium described by the (static)
probability density of a single variable (the `stress level'). We show that
heterogeneity in growth rates consistently emerges as the optimal response when
the environment is sufficiently complex and/or when perfect adjustment of
metabolic degrees of freedom is not possible (e.g. due to limited resources).
In addition, outcomes close to those achievable with unlimited resources are
often attained effectively with a modest amount of fine tuning. In other terms,
heterogeneous population structures in complex media may be rather robust with
respect to the resources available to probe the environment and adjust reaction
rates.


  In this paper we elaborate on homogeneous and heterogeneous SIR-type
epidemiological models. We find an unexpected correspondence between the
epidemic trajectory of a transmissible disease in a homogeneous SIR-type model
and radial null geodesics in the Schwarzschild spacetime. We also discuss
modeling of population heterogeneity effects by considering both a one- and
two-parameter gamma-distributed function for the initial susceptibility
distribution, and deriving the associated herd immunity threshold. We
furthermore describe how mitigation measures can be taken into account by model
fitting.


  The coexistence of multiple phytoplankton species despite their reliance on
similar resources is often explained with mean-field models assuming mixed
populations. In reality, observations of phytoplankton indicate spatial
aggregation at all scales, including at the scale of a few individuals. Local
spatial aggregation can hinder competitive exclusion since individuals then
interact mostly with other individuals of their own species, rather than
competitors from different species. To evaluate how microscale spatial
aggregation might explain phytoplankton diversity maintenance, an
individual-based, multispecies representation of cells in a hydrodynamic
environment is required. We formulate a three-dimensional and multispecies
individual-based model of phytoplankton population dynamics at the Kolmogorov
scale. The model is studied through both simulations and the derivation of
spatial moment equations, in connection with point process theory. The spatial
moment equations show a good match between theory and simulations. We
parameterized the model based on phytoplankters' ecological and physical
characteristics, for both large and small phytoplankton. Defining a zone of
potential interactions as the overlap between nutrient depletion volumes, we
show that local species composition-within the range of possible
interactions-depends on the size class of phytoplankton. In large
phytoplankton, individuals are surrounded by cells from other species, while in
small phytoplankton, individuals remain in mostly monospecific clusters.
Spatial structure therefore favours intra- over inter-specific interactions for
small phytoplankton, which likely contributes to coexistence mechanisms. Other
factors behind diversity maintenance must be examined for large phytoplankton.


  Molecular docking, which aims to find the most stable interacting
configuration of a set of molecules, is of critical importance to drug
discovery. Although a considerable number of classical algorithms have been
developed to carry out molecular docking, most focus on the limiting case of
docking two molecules. Since the number of possible configurations of N
molecules is exponential in N, those exceptions which permit docking of more
than two molecules scale poorly, requiring exponential resources to find
high-quality solutions. Here, we introduce a one-hot encoded quadratic
unconstrained binary optimization formulation (QUBO) of the multibody molecular
docking problem, which is suitable for solution by quantum annealer. Our
approach involves a classical pre-computation of pairwise interactions, which
scales only quadratically in the number of bodies while permitting well-vetted
scoring functions like the Rosetta REF2015 energy function to be used. In a
second step, we use the quantum annealer to sample low-energy docked
configurations efficiently, considering all possible docked configurations
simultaneously through quantum superposition. We show that we are able to
minimize the time needed to find diverse low-energy docked configurations by
tuning the strength of the penalty used to enforce the one-hot encoding,
demonstrating a 3-4 fold improvement in solution quality and diversity over
performance achieved with conventional penalty strengths. By mapping the
configurational search to a form compatible with current- and future-generation
quantum annealers, this work provides an alternative means of solving multibody
docking problems that may prove to have performance advantages for large
problems, potentially circumventing the exponential scaling of classical
approaches and permitting a much more efficient solution to a problem central
to drug discovery and validation pipelines.


  We consider a modified nonlocal FKPP equation to describe the spatial
organization of a single-species population in a refuge (high-quality region
immersed in a hostile environment). In addition to the spatial-dependence of
the growth rate, which is the main distinguishing feature of a refuge, it is
realistic to consider that also the mobility can be heterogeneous. We focus on
two cases where the diffusivity responds to the population density, either
decreasing or increasing with population density. This state dependency can
produce a feedback that is missing in homogeneous environments. For comparison,
we also address space-dependent diffusivity with different values inside and
outside the refuge. The objective is to characterize how the diffusion
variability affects the spatial organization of the population inside and
outside the refuge. We note that the conditions for pattern formation within
the refuge are not strongly affected by the heterogeneity in the diffusivity
profile, but this affects the shape and stability of the patterns.


  Methods from artificial intelligence (AI), in general, and machine learning,
in particular, have kept conquering new territories in numerous areas of
science. Most of the applications of these techniques are restricted to the
classification of large data sets, but new scientific knowledge can seldom be
inferred from these tools. Here we show that an AI-based amyloidogenecity
predictor can strongly differentiate the border- and the internal hexamers of
$\beta$-pleated sheets when screening all the Protein Data Bank-deposited
homology-filtered protein structures. Our main result shows that more than 30\%
of internal hexamers of $\beta$ sheets are predicted to be amyloidogenic, while
just outside the border regions, only 3\% are predicted as such. This result
may elucidate a general protection mechanism of proteins against turning into
amyloids: if the borders of $\beta$-sheets were amyloidogenic, then the whole
$\beta$ sheet could turn more easily into an insoluble amyloid-structure,
characterized by periodically repeated parallel $\beta$-sheets. We also present
that no analogous phenomenon exists on the borders of $\alpha$-helices or
randomly chosen subsequences of the studied protein structures.


  Eco-evolutionary dynamics, or eco-evolution for short, are thought to involve
rapid demography (ecology) and equally rapid phenotypic changes (evolution)
leading to novel, emergent system behaviours. This focus on contemporary
dynamics is likely due to accumulating evidence for rapid evolution, from
classical laboratory microcosms and natural populations, including the iconic
Trinidadian guppies. We argue that this view is too narrow, preventing the
successful integration of ecology and evolution. While maintaining that
eco-evolution involves emergence, we highlight that this may also be true for
slow ecology and evolution which unfold over thousands or millions of years,
such as the feedbacks between riverine geomorphology and plant evolution. We
thereby integrate geomorphology and biome-level feedbacks into eco-evolution,
significantly extending its scope. Most importantly, we emphasize that
eco-evolutionary systems need not be frozen in state-space: We identify
modulators of ecological and evolutionary rates, like temperature or
sensitivity to mutation, which can synchronize or desynchronize ecology and
evolution. We speculate that global change may increase the occurrence of
eco-evolution and emergent system behaviours which represents substantial
challenges for prediction. Our perspective represents an attempt to integrate
ecology and evolution across disciplines, from gene-regulatory networks to
geomorphology and across timescales, from contemporary dynamics to deep time.


  Adaptive therapy is a dynamic cancer treatment protocol that updates (or
"adapts") treatment decisions in anticipation of evolving tumor dynamics. This
broad term encompasses many possible dynamic treatment protocols of
patient-specific dose modulation or dose timing. Adaptive therapy maintains
high levels of tumor burden to benefit from the competitive suppression of
treatment-sensitive subpopulations on treatment-resistant subpopulations. This
evolution-based approach to cancer treatment has been integrated into several
ongoing or planned clinical trials, including treatment of metastatic castrate
resistant prostate cancer, ovarian cancer, and BRAF-mutant melanoma. In the
previous few decades, experimental and clinical investigation of adaptive
therapy has progressed synergistically with mathematical and computational
modeling. In this work, we discuss 11 open questions in cancer adaptive therapy
mathematical modeling. The questions are split into three sections: 1) the
necessary components of mathematical models of adaptive therapy 2) design and
validation of dosing protocols, and 3) challenges and opportunities in clinical
translation.


  There exist several methods for simulating biological and physical systems as
represented by chemical reaction networks. Systems with low numbers of
particles are frequently modelled as discrete-state Markov jump processes and
are typically simulated via a stochastic simulation algorithm (SSA). An SSA,
while accurate, is often unsuitable for systems with large numbers of
individuals, and can become prohibitively expensive with increasing reaction
frequency. Large systems are often modelled deterministically using ordinary
differential equations, sacrificing accuracy and stochasticity for
computational efficiency and analytical tractability. In this paper, we present
a novel hybrid technique for the accurate and efficient simulation of large
chemical reaction networks. This technique, which we name the mass-conversion
method, couples a discrete-state Markov jump process to a system of ordinary
differential equations by simulating a reaction network using both techniques
simultaneously. Individual molecules in the network are represented by exactly
one regime at any given time, and may switch their governing regime depending
on particle density. In this manner, we model high copy-number species using
the cheaper continuum method and low copy-number species using the more
expensive, discrete-state stochastic method to preserve the impact of
stochastic fluctuations at low copy number. The motivation, as with similar
methods, is to retain the advantages while mitigating the shortfalls of each
method. We demonstrate the performance and accuracy of our method for several
test problems that exhibit varying degrees of inter-connectivity and complexity
by comparing averaged trajectories obtained from both our method and from exact
stochastic simulation.


  Understanding genetic variation, e.g., through mutations, in organisms is
crucial to unravel their effects on the environment and human health. A
fundamental characterization can be obtained by solving the haplotype assembly
problem, which yields the variation across multiple copies of chromosomes.
Variations among fast evolving viruses that lead to different strains (called
quasispecies) are also deciphered with similar approaches. In both these cases,
high-throughput sequencing technologies that provide oversampled mixtures of
large noisy fragments (reads) of genomes, are used to infer constituent
components (haplotypes or quasispecies). The problem is harder for polyploid
species where there are more than two copies of chromosomes. State-of-the-art
neural approaches to solve this NP-hard problem do not adequately model
relations among the reads that are important for deconvolving the input signal.
We address this problem by developing a new method, called NeurHap, that
combines graph representation learning with combinatorial optimization. Our
experiments demonstrate substantially better performance of NeurHap in real and
synthetic datasets compared to competing approaches.


  Membrane curvature sensing is essential for a diverse range of biological
processes. Recent experiments have revealed that a single nanometer-sized
septin protein can distinguish between membrane-coated glass beads of one
micron and three micron diameters, even though the septin is orders of
magnitude smaller than the beads. This sensing ability is especially surprising
since curvature-sensing proteins must deal with persistent thermal fluctuations
of the membrane, leading to discrepancies between the bead's curvature and the
local membrane curvature sensed instantaneously by a protein. Using continuum
models of fluctuating membranes, we investigate whether it is feasible for a
protein acting as a perfect observer of the membrane to sense micron-scale
curvature either by measuring local membrane curvature or by using bilayer
lipid densities as a proxy. To do this, we develop algorithms to simulate lipid
density and membrane shape fluctuations. We derive physical limits to the
sensing efficacy of a protein in terms of protein size, membrane thickness,
membrane bending modulus, membrane-substrate adhesion strength, and bead size.
To explain the experimental protein-bead association rates, we develop two
classes of predictive models: i) for proteins that maximally associate to a
preferred curvature, and ii) for proteins with enhanced association rates above
a threshold curvature. We find that the experimentally observed sensing
efficacy is close to the theoretical sensing limits imposed on a septin-sized
protein, and that the variance in membrane curvature fluctuations sensed by a
protein determines how sharply the association rate depends on the curvature of
the bead.


  In genetically admixed populations, admixed individuals possess ancestry from
multiple source groups. Studies of human genetic admixture frequently estimate
ancestry components corresponding to fractions of individual genomes that trace
to specific ancestral populations. However, the same numerical ancestry
fraction can represent a wide array of admixture scenarios. Using a mechanistic
model of admixture, we characterize admixture genealogically: how many distinct
ancestors from the source populations does the admixture represent? We consider
African Americans, for whom continent-level estimates produce a 75-85% value
for African ancestry on average and 15-25% for European ancestry. Genetic
studies together with key features of African-American demographic history
suggest ranges for model parameters. Using the model, we infer that if
genealogical lineages of a random African American born during 1960-1965 are
traced back until they reach members of source populations, the expected number
of genealogical lines terminating with African individuals is 314, and the
expected number terminating in Europeans is 51. Across discrete generations,
the peak number of African genealogical ancestors occurs for birth cohorts from
the early 1700s. The probability exceeds 50% that at least one European
ancestor was born more recently than 1835. Our genealogical perspective can
contribute to further understanding the admixture processes that underlie
admixed populations. For African Americans, the results provide insight both on
how many of the ancestors of a typical African American might have been
forcibly displaced in the Transatlantic Slave Trade and on how many separate
European admixture events might exist in a typical African-American genealogy.


  The mobility patterns of individuals in China during the early outbreak of
the COVID-19 pandemic exhibit reversible changes -- in many regions, the
mobility first decreased significantly and later restored. Based on this
observation, here we study the classical SIR model on a particular type of
time-varying network where the links undergo a freeze-recovery process. We
first focus on an isolated network and find that the recovery mechanism could
lead to the resurgence of an epidemic. The influence of link freezing on
epidemic dynamics is subtle. In particular, we show that there is an optimal
value of the freezing rate for links which corresponds to the lowest prevalence
of the epidemic. This result challenges our conventional idea that stricter
prevention measures (corresponding to a larger freezing rate) could always have
a better inhibitory effect on epidemic spreading. We further investigate an
open system where a small fraction of nodes in the network may acquire the
disease from the "environment" (the outside infected nodes). In this case, the
second wave would appear even if the number of infected nodes has declined to
zero, which can not be explained by the isolated network model.


  Consider the diffusion process defined by the forward equation $u_t(t, x) =
\tfrac{1}{2}\{x u(t, x)\}_{xx} - \alpha \{x u(t, x)\}_{x}$ for $t, x \ge 0$ and
$-\infty < \alpha < \infty$, with an initial condition $u(0, x) = \delta(x -
x_0)$. This equation was introduced and solved by Feller to model the growth of
a population of independently reproducing individuals. We explore important
coalescent processes related to Feller's solution. For any $\alpha$ and $x_0 >
0$ we calculate the distribution of the random variable $A_n(s; t)$, defined as
the finite number of ancestors at a time $s$ in the past of a sample of size
$n$ taken from the infinite population of a Feller diffusion at a time $t$
since since its initiation. In a subcritical diffusion we find the distribution
of population and sample coalescent trees from time $t$ back, conditional on
non-extinction as $t \to \infty$. In a supercritical diffusion we construct a
coalescent tree which has a single founder and derive the distribution of
coalescent times.


  The COVID-19 epidemic has been spreading around the world for nearly three
years, and asymptomatic infections have exacerbated the spread of the epidemic.
To evaluate the role of asymptomatic infections in the spread of the epidemic,
we develop mathematical models to assess the proportion of asymptomatic
infections caused by different strains of the main covid-19 variants. The
analysis shows that when the control reproduction number is less than 1, the
disease-free equilibrium point of the model is globally asymptotically stable;
and when the control reproduction number is greater than 1, the endemic
equilibrium point exists and is unique, and is locally asymptotically stable.
We fit the epidemic data in the four time periods corresponding to the selected
614G, Alpha, Delta and Omicron variants. The fitting results show that, from
the comparison of the four time periods, the proportion of asymptomatic persons
among the infected persons gradually increased. We also predict the peak time
and peak value for the four time periods, and the results indicate that the
transmission speed and transmission intensity of the variant strains increased
to some extent. Finally, we discuss the impact of the detection ratio of
symptomatic infections on the spread of the epidemic. The results show that
with the increase of the detection ratio, the cumulative number of cases has
dropped significantly, but the decline in the proportion of asymptomatic
infections is not obvious. Therefore, in view of the hidden transmission of
asymptomatic infections, the cooperation between various epidemic prevention
and control policies is required to effectively curb the spread of the
epidemic.


  The number of confirmed COVID-19 cases reached over 1.3 million in Ontario,
Canada by June 4, 2022. The continued spread of the virus underlying COVID-19
has been spurred by the emergence of variants since the initial outbreak in
December, 2019. Much attention has thus been devoted to tracking and modelling
the transmission of COVID-19. Compartmental models are commonly used to mimic
epidemic transmission mechanisms and are easy to understand. Their performance
in real-world settings, however, needs to be more thoroughly assessed. In this
comparative study, we examine five compartmental models -- four existing ones
and an extended model that we propose -- and analyze their ability to describe
COVID-19 transmission in Ontario from January 2022 to June 2022.


  Cancer cells possess a nearly universal set of characteristics termed the
hallmarks of cancer, including replicative immortality and resisting cell
death. Dysregulated differentiation is present in virtually all cancers yet has
not yet been described as a cancer hallmark. Like other hallmarks, dysregulated
differentiation involves a breakdown of the cellular cooperation that typically
makes multicellularity possible - in this case disrupting the division of labor
among the cells of a body. At the time that the original hallmarks of cancer
were described, it was not known that dysregulated differentiation was
mechanistically distinct from growth inhibition, but now that this is known, it
is a further reason to consider dysregulated differentiation a hallmark of
cancer. Dysregulated differentiation also has clinical utility, as it forms the
basis of pathological grading, predicts clinical outcomes, and is a viable
target for therapies aimed at inducing differentiation. Here we argue that
hallmarks of cancer should be near universal, mechanistically distinct, and
have clinical utility for prognosis and/or therapy. Dysregulated
differentiation meets all of these criteria.


  Estimating a causal query from observational data is an essential task in the
analysis of biomolecular networks. Estimation takes as input a network
topology, a query estimation method, and observational measurements on the
network variables. However, estimations involving many variables can be
experimentally expensive, and computationally intractable. Moreover, using the
full set of variables can be detrimental, leading to bias, or increasing the
variance in the estimation. Therefore, designing an experiment based on a
well-chosen subset of network components can increase estimation accuracy, and
reduce experimental and computational costs. We propose a simulation-based
algorithm for selecting sub-networks that support unbiased estimators of the
causal query under a constraint of cost, ranked with respect to the variance of
the estimators. The simulations are constructed based on historical
experimental data, or based on known properties of the biological system. Three
case studies demonstrated the effectiveness of well-chosen network subsets for
estimating causal queries from observational data. All the case studies are
reproducible and available at https://github.com/srtaheri/Simplified_LVM.


  There has been a recent surge of interest in what causes aging. This has been
matched by unprecedented research investment in the field from tech companies.
But, despite considerable effort from a broad range of researchers, we do not
have a rigorous mathematical theory of programmed aging. To address this, we
recently derived a mortality equation that governs the transition matrix of an
evolving population with a given maximum age. Here, we characterize the
spectrum of eigenvalues of the solution to this equation. The eigenvalues fall
into two classes. The complex and negative real eigenvalues, which we call the
flower, are always contained in the unit circle in the complex plane. They play
a negligible role in controlling the dynamics of an aging population. The
positive real eigenvalues, which we call the stem, are the only eigenvalues
which can exceed the unit circle. They control the most important properties of
the dynamics. In particular, the spectral radius increases with the maximum
allowed age. This suggests that programmed aging confers no advantage in a
constant environment. However, the spectral gap, which governs the rate of
convergence to equilibrium, decreases with the maximum allowed age. This opens
the door to an evolutionary advantage in a changing environment.


  Structure-based drug design (SBDD) aims to design small-molecule ligands that
bind with high affinity and specificity to pre-determined protein targets. In
this paper, we formulate SBDD as a 3D-conditional generation problem and
present DiffSBDD, an SE(3)-equivariant 3D-conditional diffusion model that
generates novel ligands conditioned on protein pockets. Comprehensive in silico
experiments demonstrate the efficiency and effectiveness of DiffSBDD in
generating novel and diverse drug-like ligands with competitive docking scores.
We further explore the flexibility of the diffusion framework for a broader
range of tasks in drug design campaigns, such as off-the-shelf property
optimization and partial molecular design with inpainting.


  In many organisms, homologous (or repetitive) chromosomal regions can
associate or/and undergo concerted epigenetic changes in the absence of DNA
breakage and recombination. The direct specific pairing of DNA duplexes with
similar nucleotide sequences represents an attractive mechanism for recognizing
such regions. Whereas the pairing of B-DNA duplexes may involve a large energy
barrier, C-DNA duplexes are expected to pair much more readily. This unique
feature of C-DNA is largely due to the fact that its major groove is wide and
very shallow, permitting almost perfect initial homologous contacts between two
duplexes without clashing. Overall, the conjectured role of C-DNA in
recombination-independent pairing should revive the efforts to understand its
structure and function in the cell.


  Understanding the temporal spread of gene drive alleles -- alleles that bias
their own transmission -- through modeling is essential before any field
experiments. In this paper, we present a deterministic reaction-diffusion model
describing the interplay between demographic and allelic dynamics, in a
one-dimensional spatial context. We focused on the traveling wave solutions,
and more specifically, on the speed of gene drive invasion (if successful). We
considered various timings of gene conversion (in the zygote or in the
germline) and different probabilities of gene conversion (instead of assuming
100$\%$ conversion as done in a previous work). We compared the types of
propagation when the intrinsic growth rate of the population takes extreme
values, either very large or very low. When it is infinitely large, the wave
can be either successful or not, and, if successful, it can be either pulled or
pushed, in agreement with previous studies (extended here to the case of
partial conversion). In contrast, it cannot be pushed when the intrinsic growth
rate is vanishing. In this case, analytical results are obtained through an
insightful connection with an epidemiological SI model. We conducted extensive
numerical simulations to bridge the gap between the two regimes of large and
low growth rate. We conjecture that, if it is pulled in the two extreme
regimes, then the wave is always pulled, and the wave speed is independent of
the growth rate. This occurs for instance when the fitness cost is small
enough, or when there is stable coexistence of the drive and the wild-type in
the population after successful drive invasion. Our model helps delineate the
conditions under which demographic dynamics can affect the spread of a gene
drive.


  Antimicrobial resistance is a silent pandemic that is being exacerbated by
the uncontrolled use of antibiotics. Since the discovery of penicillin, we have
been largely dependent on microbe-derived small molecules to treat bacterial
infections. However, the golden era of antibiotics is coming to an end as the
emergence and spread of antimicrobial resistance against these antibacterial
compounds is outpacing the discovery and development of new antibiotics. The
current antibiotic market suffers from various shortcomings, including the
absence of profitability and investment. However, the most important underlying
issue of traditional antibiotics arises from the inherent properties of these
small molecules being mostly broad-spectrum and non-programmable. As the
scientific knowledge of microbes progresses, the scientific community is
starting to explore entirely novel approaches to tackling antimicrobial
resistance. One of the most prominent approaches is to develop next-generation
antibiotics. In this review, we discuss three innovations of next-generation
antibiotics compared to traditional antibiotics as specificity, evolvability,
and non-immunogenicity. We present a number of potential antimicrobial agents,
including bacteriophage-based therapy, CRISPR-Cas-based antimicrobials, and
microbiome-derived antimicrobial compounds. These alternative antimicrobial
agents possess innovative properties that may overcome the inherent
shortcomings of traditional antibiotics, and some of these next-generation
antibiotics are not merely far-fetched ideas but are currently in clinical
development. We further discuss some related issues and challenges such as
infection diagnostics and regulatory frameworks that still need to be addressed
to bring these next-generation antibiotics to the antibiotic market as viable
products to combat antimicrobial resistance using a diversified set of
strategies.


  Uncertainty can be classified as either aleatoric (intrinsic randomness) or
epistemic (imperfect knowledge of parameters). The majority of frameworks
assessing infectious disease risk consider only epistemic uncertainty. We only
ever observe a single epidemic, and therefore cannot empirically determine
aleatoric uncertainty. Here, we characterise both epistemic and aleatoric
uncertainty using a time-varying general branching process. Our framework
explicitly decomposes aleatoric variance into mechanistic components,
quantifying the contribution to uncertainty produced by each factor in the
epidemic process, and how these contributions vary over time. The aleatoric
variance of an outbreak is itself a renewal equation where past variance
affects future variance. We find that, superspreading is not necessary for
substantial uncertainty, and profound variation in outbreak size can occur even
without overdispersion in the offspring distribution (i.e. the distribution of
the number of secondary infections an infected person produces). Aleatoric
forecasting uncertainty grows dynamically and rapidly, and so forecasting using
only epistemic uncertainty is a significant underestimate. Therefore, failure
to account for aleatoric uncertainty will ensure that policymakers are misled
about the substantially higher true extent of potential risk. We demonstrate
our method, and the extent to which potential risk is underestimated, using two
historical examples.


  The search for life beyond the Earth is the overarching goal of the NASA
Astrobiology Program, and it underpins the science of missions that explore the
environments of Solar System planets and exoplanets. However, the detection of
extraterrestrial life, in our Solar System and beyond, is sufficiently
challenging that it is likely that multiple measurements and approaches,
spanning disciplines and missions, will be needed to make a convincing claim.
Life detection will therefore not be an instantaneous process, and it is
unlikely to be unambiguous-yet it is a high-stakes scientific achievement that
will garner an enormous amount of public interest. Current and upcoming
research efforts and missions aimed at detecting past and extant life could be
supported by a consensus framework to plan for, assess and discuss life
detection claims (c.f. Green et al., 2021). Such a framework could help
increase the robustness of biosignature detection and interpretation, and
improve communication with the scientific community and the public. In response
to this need, and the call to the community to develop a confidence scale for
standards of evidence for biosignature detection (Green et al., 2021), a
community-organized workshop was held on July 19-22, 2021. The meeting was
designed in a fully virtual (flipped) format. Preparatory materials including
readings, instructional videos and activities were made available prior to the
workshop, allowing the workshop schedule to be fully dedicated to active
community discussion and prompted writing sessions. To maximize global
interaction, the discussion components of the workshop were held during
business hours in three different time zones, Asia/Pacific, European and US,
with daily information hand-off between group organizers.


  Gene regulation is a dynamic process that connects genotype and phenotype.
Given the difficulty of physically mapping mammalian gene circuitry, we require
new computational methods to learn regulatory rules. Natural language is a
valuable analogy to the communication of regulatory control. Machine learning
systems model natural language by explicitly learning context dependencies
between words. We propose a similar system applied to single-cell RNA
expression profiles to learn context dependencies between genes. Our model,
Exceiver, is trained across a diversity of cell types using a self-supervised
task formulated for discrete count data, accounting for feature sparsity. We
found agreement between the similarity profiles of latent sample
representations and learned gene embeddings with respect to biological
annotations. We evaluated Exceiver on a new dataset and a downstream prediction
task and found that pretraining supports transfer learning. Our work provides a
framework to model gene regulation on a single-cell level and transfer
knowledge to downstream tasks.


  Predicting cellular metabolic states is a central problem in biophysics.
Conventional approaches, however, sensitively depend on the microscopic details
of individual metabolic systems. In this Letter, we derived a universal linear
relationship between the metabolic responses against nutrient conditions and
metabolic inhibition, with the aid of a microeconomic theory. The relationship
holds in arbitrary metabolic systems as long as the law of mass conservation
stands, as supported by extensive numerical calculations. It offers
quantitative predictions without prior knowledge of systems.


  The cultivated strawberry Fragaria ananassa Duch. is widely planted in
greenhouses in China. Its production heavily depends on pollination services.
Compared with artificial pollination, bee pollination can significantly improve
fruit quality and save considerable labor requirement. Multiple factors such as
bee foraging behavior, planting pattern and the spatial complexity of the
greenhouse environment interacting over time and space are major obstacles to
understanding of bee pollination dynamics. We propose a spatially-explicit
agent-based simulation model which allows users to explore how various factors
including bee foraging behavior and strawberry phenology conditions as well as
the greenhouse environment influence pollination efficiency and fruit quality.
Simulation experiments allowed us to compare pollination efficiencies in
different conditions. Especially, the cause of bee pollination advantage,
optimal bee density and bee hive location were discussed based on sensitivity
analysis. In addition, simulation results provide some insights for strawberry
planting in a greenhouse. The firmly validated open-source model is a useful
tool for hypothesis testing and theory development for strawberry pollination
research.


  Machine learning models have found numerous successful applications in
computational drug discovery. A large body of these models represents molecules
as sequences since molecular sequences are easily available, simple, and
informative. The sequence-based models often segment molecular sequences into
pieces called chemical words (analogous to the words that make up sentences in
human languages) and then apply advanced natural language processing techniques
for tasks such as $\textit{de novo}$ drug design, property prediction, and
binding affinity prediction. However, the chemical characteristics and
significance of these building blocks, chemical words, remain unexplored. This
study aims to investigate the chemical vocabularies generated by popular
subword tokenization algorithms, namely Byte Pair Encoding (BPE), WordPiece,
and Unigram, and identify key chemical words associated with protein-ligand
binding. To this end, we build a language-inspired pipeline that treats high
affinity ligands of protein targets as documents and selects key chemical words
making up those ligands based on tf-idf weighting. Further, we conduct case
studies on a number of protein families to analyze the impact of key chemical
words on binding. Through our analysis, we find that these key chemical words
are specific to protein targets and correspond to known pharmacophores and
functional groups. Our findings will help shed light on the chemistry captured
by the chemical words, and by machine learning models for drug discovery at
large.


  All organisms, fundamentally, are made from the same raw material, namely the
elements of the periodic table. Biochemical diversity is achieved with how
these elements are utilized, for what purpose and in which physical location.
Determining elemental distributions, especially those of trace elements that
facilitate metabolism as cofactors in the active centers of essential enzymes,
can determine the state of metabolism, the nutritional status or the
developmental stage of an organism. Photosynthetic eukaryotes, especially
algae, are excellent subjects for quantitative analysis of elemental
distribution. These microbes utilize unique metabolic pathways that require
various trace nutrients at their core to enable its operation. Photosynthetic
microbes also have important environmental roles as primary producers in
habitats with limited nutrient supply or toxin contaminations. Accordingly,
photosynthetic eukaryotes are of great interest for biotechnological
exploitation, carbon sequestration and bioremediation, with many of the
applications involving various trace elements and consequently affecting their
quota and intracellular distribution. A number of diverse applications were
developed for elemental imaging allowing subcellular resolution, with X-ray
fluorescence microscopy (XFM) being at the forefront, enabling quantitative
descriptions of intact cells in a non-destructive method. This Tutorial Review
summarizes the workflow of a quantitative, single-cell elemental distribution
analysis of a eukaryotic alga using XFM.


  Many ncRNAs function through RNA-RNA interactions. Fast and reliable RNA
structure prediction with consideration of RNA-RNA interaction is useful. Some
existing tools are less accurate due to omitting the competing of
intermolecular and intramolecular base pairs, or focus more on predicting the
binding region rather than predicting the complete secondary structure of two
interacting strands. Vienna RNAcofold, which reduces the problem into the
classical single sequence folding by concatenating two strands, scales in cubic
time against the combined sequence length, and is slow for long sequences. To
address these issues, we present LinearCoFold, which predicts the complete
minimum free energy structure of two strands in linear runtime, and
LinearCoPartition, which calculates the cofolding partition function and base
pairing probabilities in linear runtime. LinearCoFold and LinearCoPartition
follows the concatenation strategy of RNAcofold, but are orders of magnitude
faster than RNAcofold. For example, on a sequence pair with combined length of
26,190 nt, LinearCoFold is 86.8x faster than RNAcofold MFE mode (0.6 minutes
vs. 52.1 minutes), and LinearCoPartition is 642.3x faster than RNAcofold
partition function mode (1.8 minutes vs. 1156.2 minutes). Different from the
local algorithms, LinearCoFold and LinearCoPartition are global cofolding
algorithms without restriction on base pair length. Surprisingly, LinearCoFold
and LinearCoPartition's predictions have higher PPV and sensitivity of
intermolecular base pairs. Furthermore, we apply LinearCoFold to predict the
RNA-RNA interaction between SARS-CoV-2 gRNA and human U4 snRNA, which has been
experimentally studied, and observe that LinearCoFold's prediction correlates
better to the wet lab results.


  We constructed a Susceptible-Addicted-Reformed model and explored the
dynamics of nonlinear relapse in the Reformed population. The transition from
susceptible considered {\it at-risk} is modeled using a strictly decreasing
general function, mimicking an influential factor that reduces the flow into
the addicted class. The {\it basic reproductive number} is computed.
Furthermore, $R_0$ determines the local asymptotically stability of the
addicted-free equilibrium. Conditions for a forward-backward bifurcation were
established using $R_0$ and other threshold quantities. A stochastic version of
the model is presented, and some numerical examples are shown. Results showed
that the influence of the temporarily reformed individuals is highly sensitive
to the initial addicted population.


  The ubiquitous existence of microbial communities marks the importance of
understanding how species interact within the community to coexist and their
spatial organization. We study a two-species mutualistic cross-feeding model
through a stochastic cellular automaton on a square lattice using kinetic Monte
Carlo simulation. Our model encapsulates the essential dynamic processes such
as cell growth, and nutrient excretion, diffusion and uptake. Focusing on the
interplay among nutrient diffusion and individual cell division, we discover
three general classes of colony morphology: co-existing sectors, co-existing
spirals, and engulfment. When the cross-feeding nutrient is widely available,
either through high excretion or fast diffusion, a stable circular colony with
alternating species sector emerges. When the consumer cells rely on being
spatially close to the producers, we observe a stable spiral. We also see one
species being engulfed by the other when species interfaces merge due to
stochastic fluctuation. By tuning the diffusion rate and the growth rate, we
are able to gain quantitative insights into the structures of the sectors and
the spirals.


  Social dilemmas are situations in which collective welfare is at odds with
individual gain. One widely studied example, due to the conflict it poses
between human behaviour and game theoretic reasoning, is the Traveler's
Dilemma. The dilemma relies on the players' incentive to undercut their
opponent at the expense of losing a collective high payoff. Such individual
incentive leads players to a systematic mutual undercutting until the lowest
possible payoff is reached, which is the game's unique Nash equilibrium.
However, if players were satisfied with a high payoff -- that is not
necessarily higher than their opponent's -- they would both be better off
individually and collectively. Here, we explain how it is possible to converge
to this cooperative high payoff equilibrium. Our analysis focuses on
decomposing the dilemma into a local and a global game. We show that players
need to escape the local maximisation and jump to the global game, in order to
reach the cooperative equilibrium. Using a dynamic approach, based on
evolutionary game theory and learning theory models, we find that diversity,
understood as the presence of suboptimal strategies, is the general mechanism
that enables the jump towards cooperation.


  Binding affinity prediction of three-dimensional (3D) protein ligand
complexes is critical for drug repositioning and virtual drug screening.
Existing approaches transform a 3D protein-ligand complex to a two-dimensional
(2D) graph, and then use graph neural networks (GNNs) to predict its binding
affinity. However, the node and edge features of the 2D graph are extracted
based on invariant local coordinate systems of the 3D complex. As a result, the
method can not fully learn the global information of the complex, such as, the
physical symmetry and the topological information of bonds. To address these
issues, we propose a novel Equivariant Line Graph Network (ELGN) for affinity
prediction of 3D protein ligand complexes. The proposed ELGN firstly adds a
super node to the 3D complex, and then builds a line graph based on the 3D
complex. After that, ELGN uses a new E(3)-equivariant network layer to pass the
messages between nodes and edges based on the global coordinate system of the
3D complex. Experimental results on two real datasets demonstrate the
effectiveness of ELGN over several state-of-the-art baselines.


  Molecule optimization is an important problem in chemical discovery and has
been approached using many techniques, including generative modeling,
reinforcement learning, genetic algorithms, and much more. Recent work has also
applied zeroth-order (ZO) optimization, a subset of gradient-free optimization
that solves problems similarly to gradient-based methods, for optimizing latent
vector representations from an autoencoder. In this paper, we study the
effectiveness of various ZO optimization methods for optimizing molecular
objectives, which are characterized by variable smoothness, infrequent optima,
and other challenges. We provide insights on the robustness of various ZO
optimizers in this setting, show the advantages of ZO sign-based gradient
descent (ZO-signGD), discuss how ZO optimization can be used practically in
realistic discovery tasks, and demonstrate the potential effectiveness of ZO
optimization methods on widely used benchmark tasks from the Guacamol suite.
Code is available at: https://github.com/IBM/QMO-bench.


  Disordered proteins and nucleic acids play key roles in cellular function and
disease. Here we review recent advances in the computational exploration of the
conformational dynamics of flexible biomolecules. We focus on hierarchical
chain growth (HCG) from fragment libraries built with atomistic molecular
dynamics simulations. HCG combines chain fragments in a statistically
reproducible manner into ensembles of full-length atomically detailed
biomolecular structures. The input fragment structures are typically collected
from molecular dynamics simulations, but could also come from structural
databases. Experimental data can be integrated during and after chain assembly.
Applications to the neurodegeneration-linked proteins $\alpha$-synuclein, tau,
and TDP-43, including as condensate, illustrate the use of HCG. We conclude by
highlighting the emerging connections to AI-based structural modeling.