  The infectivity of SARS-CoV-2 depends on the binding affinity of the
receptor-binding domain (RBD) of the spike protein with the angiotensin
converting enzyme 2 (ACE2) receptor. The calculated RBD-ACE2 binding energies
indicate that the difference in transmission efficiency of SARS-CoV-2 variants
cannot be fully explained by electrostatic interactions, hydrogen-bond
interactions, van der Waals interactions, internal energy, and nonpolar
solvation energies. Here, we demonstrate that low-entropy regions of hydration
shells around proteins drive hydrophobic attraction between shape-matched
low-entropy regions of the hydration shells, which essentially coordinates
protein-protein binding in rotational-configurational space of mutual
orientations and determines the binding affinity. An innovative method was used
to identify the low-entropy regions of the hydration shells of the RBDs of
multiple SARS-CoV-2 variants and the ACE2. We observed integral low-entropy
regions of hydration shells covering the binding sites of the RBDs and matching
in shape to the low-entropy region of hydration shell at the binding site of
the ACE2. The RBD-ACE2 binding is thus found to be guided by hydrophobic
collapse between the shape-matched low-entropy regions of the hydration shells.
A measure of the low-entropy of the hydration shells can be obtained by
counting the number of hydrophilic groups expressing hydrophilicity within the
binding sites. The low-entropy level of hydration shells at the binding site of
a spike protein is found to be an important indicator of the contagiousness of
the coronavirus.


  We study the distribution of the unobserved states of two measure-valued
diffusions of Fleming-Viot and Dawson-Watanabe type, conditional on
observations from the underlying populations collected at past, present and
future times. If seen as nonparametric hidden Markov models, this amounts to
finding the smoothing distributions of these processes, which we show can be
explicitly described in recursive form as finite mixtures of laws of Dirichlet
and gamma random measures respectively. We characterize the time-dependent
weights of these mixtures, accounting for potentially different time intervals
between data collection times, and fully describe the implications of assuming
a discrete or a nonatomic distribution for the underlying process that drives
mutations. In particular, we show that with a nonatomic mutation offspring
distribution, the inference automatically upweights mixture components that
carry, as atoms, observed types shared at different collection times. The
predictive distributions for further samples from the population conditional on
the data are also identified and shown to be mixtures of generalized Polya
urns, conditionally on a latent variable in the Dawson-Watanabe case.


  The paper deals with two interrelated topics, identifying transient
amplifiers in an iterative process and analyzing the process by its spectral
dynamics, which is the change in the graph spectra by edge manipulations.
Transient amplifiers are networks representing population structures which
shift the balance between natural selection and random drift. Thus, amplifiers
are highly relevant for understanding the relationships between spatial
structures and evolutionary dynamics. We study an iterative procedure to
identify transient amplifiers for death-Birth updating. The algorithm starts
with a regular input graph and iteratively removes edges until desired
structures are achieved. Thus, a sequence of candidate graphs is obtained. The
edge removals are guided by quantities derived from the sequence of candidate
graphs. Moreover, we are interested in the Laplacian spectra of the candidate
graphs and analyze the iterative process by its spectral dynamics. The results
show that although transient amplifiers for death-Birth updating are rare, a
substantial number of them can be obtained by the proposed procedure. The
graphs identified share structural properties and have some similarity to
dumbbell and barbell graphs. Also, the spectral dynamics possesses
characteristic features useful for deducing links between structural and
spectral properties and for distinguishing transient amplifiers among
evolutionary graphs in general.


  Computational protein design has the potential to deliver novel molecular
structures, binders, and catalysts for myriad applications. Recent neural
graph-based models that use backbone coordinate-derived features show
exceptional performance on native sequence recovery tasks and are promising
frameworks for design. A statistical framework for modeling protein sequence
landscapes using Tertiary Motifs (TERMs), compact units of recurring structure
in proteins, has also demonstrated good performance on protein design tasks. In
this work, we investigate the use of TERM-derived data as features in neural
protein design frameworks. Our graph-based architecture, TERMinator,
incorporates TERM-based and coordinate-based information and outputs a Potts
model over sequence space. TERMinator outperforms state-of-the-art models on
native sequence recovery tasks, suggesting that utilizing TERM-based and
coordinate-based features together is beneficial for protein design.


  Epistasis occurs when the effect of a mutation depends on its carrier's
genetic background. Despite increasing evidence that epistasis for fitness is
common, its role during evolution is contentious. Fitness landscapes, mappings
of genotype or phenotype to fitness, capture the full extent and complexity of
epistasis. Fitness landscape theory has shown how epistasis affects the course
and the outcome of evolution. Moreover, by measuring the competitive fitness of
sets of tens to thousands of connected genotypes, empirical fitness landscapes
have shown that epistasis is frequent and depends on the fitness measure, the
choice of mutations for the landscape, and the environment in which it was
measured. Here, I review fitness landscape theory and experiments and their
implications for the role of epistasis in adaptation. I discuss theoretical
expectations in the light of empirical fitness landscapes and highlight open
challenges and future directions towards integrating theory and data, and
incorporating ecological factors.


  Rooted acyclic graphs appear naturally when the phylogenetic relationship of
a set $X$ of taxa involves not only speciations but also recombination,
horizontal transfer, or hybridization, that cannot be captured by trees. A
variety of classes of such networks have been discussed in the literature,
including phylogenetic, level-1, tree-child, tree-based, galled tree, regular,
or normal networks as models of different types of evolutionary processes.
Clusters arise in models of phylogeny as the sets $\mathtt{C}(v)$ of descendant
taxa of a vertex $v$. The clustering system $\mathscr{C}_N$ comprising the
clusters of a network $N$ conveys key information on $N$ itself. In the special
case of rooted phylogenetic trees, $T$ is uniquely determined by its clustering
system $\mathscr{C}_T$. Although this is no longer true for networks in
general, it is of interest to relate properties of $N$ and $\mathscr{C}_N$.
Here, we systematically investigate the relationships of several well-studied
classes of networks and their clustering systems. The main results are
correspondences of classes of networks and clustering system of the following
form: If $N$ is a network of type $\mathbb{X}$, then $\mathcal{C}_N$ satisfies
$\mathbb{Y}$, and conversely if $\mathscr{C}$ is a clustering system satisfying
$\mathbb{Y}$ then there is network $N$ of type $\mathbb{X}$ such that
$\mathscr{C}\subseteq\mathscr{C}_N$.This, in turn, allows us to investigate the
mutual dependencies between the distinct types of networks in much detail.


  Single-cell transcriptomics enabled the study of cellular heterogeneity in
response to perturbations at the resolution of individual cells. However,
scaling high-throughput screens (HTSs) to measure cellular responses for many
drugs remains a challenge due to technical limitations and, more importantly,
the cost of such multiplexed experiments. Thus, transferring information from
routinely performed bulk RNA HTS is required to enrich single-cell data
meaningfully. We introduce chemCPA, a new encoder-decoder architecture to study
the perturbational effects of unseen drugs. We combine the model with an
architecture surgery for transfer learning and demonstrate how training on
existing bulk RNA HTS datasets can improve generalisation performance. Better
generalisation reduces the need for extensive and costly screens at single-cell
resolution. We envision that our proposed method will facilitate more efficient
experiment designs through its ability to generate in-silico hypotheses,
ultimately accelerating drug discovery.


  Genomics data such as RNA gene expression, methylation and micro RNA
expression are valuable sources of information for various clinical predictive
tasks. For example, predicting survival outcomes, cancer histology type and
other patients' related information is possible using not only clinical data
but molecular data as well. Moreover, using these data sources together, for
example in multitask learning, can boost the performance. However, in practice,
there are many missing data points which leads to significantly lower patient
numbers when analysing full cases, which in our setting refers to all
modalities being present.
  In this paper we investigate how imputing data with missing values using deep
learning coupled with multitask learning can help to reach state-of-the-art
performance results using combined genomics modalities, RNA, micro RNA and
methylation. We propose a generalised deep imputation method to impute values
where a patient has all modalities present except one. Interestingly enough,
deep imputation alone outperforms multitask learning alone for the
classification and regression tasks across most combinations of modalities. In
contrast, when using all modalities for survival prediction we observe that
multitask learning alone outperforms deep imputation alone with statistical
significance (adjusted p-value 0.03). Thus, both approaches are complementary
when optimising performance for downstream predictive tasks.


  Protecting interventions of many types (both pharmaceutical and
non-pharmaceutical) can be deployed against the spreading of a communicable
disease, as the worldwide COVID-19 pandemic has dramatically shown. Here we
investigate in detail the effects at the population level of interventions that
provide an asymmetric protection between the people involved in a single
interaction. Masks of different filtration types, either protecting mainly the
wearer or the contacts of the wearer, are a prominent example of these
interventions. By means of analytical calculations and extensive simulations of
simple epidemic models on networks, we show that interventions protecting more
efficiently the adopter (e.g the mask wearer) are more effective than
interventions protecting primarily the contacts of the adopter in reducing the
prevalence of the disease and the number of concurrently infected individuals
("flattening the curve"). This observation is backed up by the study of a more
realistic epidemic model on an empirical network representing the patterns of
contacts in the city of Portland. Our results point out that promoting
wearer-protecting face masks and other self-protecting interventions, though
deemed selfish and inefficient, can actually be a better strategy to
efficiently curtail pandemic spreading.


  Purpose: In this study we propose to identify significant economic,
political, and social factors that helped contribute to the incidence and
mortality of SARS-CoV-2 around the world during the first seven months of the
COVID-19 pandemic. Methods: Retrospective analyses were conducted using data
acquired by the World Health Organization regarding the number of new cases and
deaths from 76 countries from January through July 2020. Hierarchical cluster
and principal component analyses were used to categorize different countries
based on their ability to control the outbreak based on number of cases and
deaths throughout these initial months of the pandemic. Furthermore, using
non-parametrical analyses, different variables were correlated with multiple
economic, political, and social indexes. Results: Multiple significant
correlations were found, including the control of corruption, government
effectiveness, and the human life indicator, with indexes containing the higher
correlations demonstrating the capability of individual countries to help drive
the course of this pandemic. Geographical hotspots were identified that
required further attention to aid in outbreak control. Conclusion:
Collectively, the findings in this study indicate that multiple factors beyond
the biological should be considered in the implementation of measures to
contain this pandemic


  The absence of a conventional association between the cell-cell cohabitation
and its emergent dynamics into cliques during development has hindered our
understanding of how cell populations proliferate, differentiate, and compete,
i.e. the cell ecology. With the recent advancement of the single-cell
RNA-sequencing (RNA-seq), we can potentially describe such a link by
constructing network graphs that characterize the similarity of the gene
expression profiles of the cell-specific transcriptional programs, and
analyzing these graphs systematically using the summary statistics informed by
the algebraic topology. We propose the single-cell topological simplicial
analysis (scTSA). Applying this approach to the single-cell gene expression
profiles from local networks of cells in different developmental stages with
different outcomes reveals a previously unseen topology of cellular ecology.
These networks contain an abundance of cliques of single-cell profiles bound
into cavities that guide the emergence of more complicated habitation forms. We
visualize these ecological patterns with topological simplicial architectures
of these networks, compared with the null models. Benchmarked on the
single-cell RNA-seq data of zebrafish embryogenesis spanning 38,731 cells, 25
cell types and 12 time steps, our approach highlights the gastrulation as the
most critical stage, consistent with consensus in developmental biology. As a
nonlinear, model-independent, and unsupervised framework, our approach can also
be applied to tracing multi-scale cell lineage, identifying critical stages, or
creating pseudo-time series.


  We investigate solutions to the TAP equation, a phenomenological
implementation of the Theory of the Adjacent Possible. Several implementations
of TAP are studied, with potential applications in a range of topics including
economics, social sciences, environmental change, evolutionary biological
systems, and the nature of physical laws. The generic behaviour is an extended
plateau followed by a sharp explosive divergence. We find accurate analytic
approximations for the blow-up time that we validate against numerical
simulations, and explore the properties of the equation in the vicinity of
equilibrium between innovation and extinction. A particular variant, the
two-scale TAP model, replaces the initial plateau with a phase of exponential
growth, a widening of the TAP equation phenomenology that may enable it to be
applied in a wider range of contexts.


  Fermat-Weber points with respect to an asymmetric tropical distance function
are studied. It turns out that they correspond to the optimal solutions of a
transportation problem. The results are applied to obtain a new method for
computing consensus trees in phylogenetics. This method has several desirable
properties; e.g., it is Pareto and co-Pareto on rooted triplets.


  Understanding the mechanisms of species coexistence has always been a
fundamental topic in ecology. Classical theory predicts that interspecific
competition may select for traits that stabilize niche differences, although
recent work shows that this is not strictly necessary. Here we ask whether
adaptive phenotypic plasticity could allow species coexistence (i.e., some
stability at an equilibrium point) without ecological differentiation in
habitat use. We used individual-based stochastic simulations defining a
landscape composed of spatially uncorrelated or autocorrelated environmental
patches, where two species with the same competitive strategies, not able to
coexist without some form of phenotypic plasticity, expanded their ranges in
the absence of a competition-colonization trade-off (a well-studied mechanism
for species diversity). Each patch is characterized by a random environmental
value that determines the optimal phenotype of its occupants. In such a
scenario, only local adaptation and gene flow (migration) may interact to
promote genetic variation and coexistence in the metapopulation. Results show
that a competitively inferior species with adaptive phenotypic plasticity can
coexist in a same patch with a competitively superior, non-plastic species,
provided the migration rates and variances of the patches' environmental values
are sufficiently large.


  Ecosystems are formed by networks of species and their interactions.
Traditional models of such interactions assume a constant interaction strength
between a given pair of species. However, there is often significant trait
variation among individual organisms even within the same species, causing
heterogeneity in their interaction strengths with other species. The
consequences of such heterogeneous interactions for the ecosystem have not been
studied systematically. As a theoretical exploration, we analyze a simple
ecosystem with trophic interactions between two predators and a shared prey,
which would exhibit competitive exclusion in models with homogeneous
interactions. We consider several scenarios where individuals of the prey
species differentiate into subpopulations with different interaction strengths.
We show that in all these cases, whether the heterogeneity is inherent,
reversible, or adaptive, the ecosystem can stabilize at a new equilibrium where
all three species coexist. Moreover, the prey population that has heterogeneous
interactions with its predators reaches a higher density than it would without
heterogeneity, and can even reach a higher density in the presence of two
predators than with just one. Our results suggest that heterogeneity may be a
naturally selected feature of ecological interactions that have important
consequences for the stability and diversity of ecosystems.


  Due to its high transmissibility, Omicron BA.1 ousted the Delta variant to
become a dominating variant in late 2021 and was replaced by more transmissible
Omicron BA.2 in March 2022. An important question is which new variants will
dominate in the future. Topology-based deep learning models have had tremendous
success in forecasting emerging variants in the past. However, topology is
insensitive to homotopic shape variations in virus-human protein-protein
binding, which are crucial to viral evolution and transmission. This challenge
is tackled with persistent Laplacian, which is able to capture both the
topology and shape of data. Persistent Laplacian-based deep learning models are
developed to systematically evaluate variant infectivity. Our comparative
analysis of Alpha, Beta, Gamma, Delta, Lambda, Mu, and Omicron BA.1, BA.1.1,
BA.2, BA.2.11, BA.2.12.1, BA.3, BA.4, and BA.5 unveils that Omicron BA.2.11,
BA.2.12.1, BA.3, BA.4, and BA.5 are more contagious than BA.2. In particular,
BA.4 and BA.5 are about 36\% more infectious than BA.2 and are projected to
become new dominating variants by natural selection. Moreover, the proposed
models outperform the state-of-the-art methods on three major benchmark
datasets for mutation-induced protein-protein binding free energy changes.


  Lotka-Volterra (LV) equations play a key role in the mathematical modeling of
various ecological, biological and chemical systems. When the number of species
(or, depending on the viewpoint, chemical components) becomes large, basic but
fundamental questions such as computing the number of surviving species still
lack theoretical answers. In this paper, we consider a large system of LV
equations where the interactions between the various species are a realization
of a random matrix. We provide conditions to have a unique equilibrium and
present a heuristics to compute the number of surviving species. This
heuristics combines arguments from Random Matrix Theory, mathematical
optimization (LCP), and standard extreme value theory. Numerical simulations,
together with an empirical study where the strength of interactions evolves
with time, illustrate the accuracy and scope of the results.


  We consider a stochastic model for the evolution of a discrete population
structured by a trait with values on a finite grid of the torus, and with
mutation and selection. Traits are vertically inherited unless a mutation
occurs, and influence the birth and death rates. We focus on a parameter
scaling where population is large, individual mutations are small but not rare,
and the grid mesh for the trait values is much smaller than the size of
mutation steps. When considering the evolution of the population in a long time
scale, the contribution of small sub-populations may strongly influence the
dynamics. Our main result quantifies the asymptotic dynamics of sub-population
sizes on a logarithmic scale. We establish that under the parameter scaling the
logarithm of the stochastic population size process, conveniently normalized,
converges to the unique viscosity solution of a Hamilton-Jacobi equation. Such
Hamilton-Jacobi equations have already been derived from parabolic
integro-differential equations and have been widely developed in the study of
adaptation of quantitative traits. Our work provides a justification of this
framework directly from a stochastic individual based model, leading to a
better understanding of the results obtained within this approach. The proof
makes use of almost sure maximum principles and careful controls of the
martingale parts.


  Life is a nonequilibrium phenomenon: metabolism provides a continuous supply
of energy that drives nearly all cellular processes. However, very little is
known about how much energy different cellular processes use, i.e. their
energetic costs. The most direct experimental measurements of these costs
involve modulating the activity of cellular processes and determining the
resulting changes in energetic fluxes. In this review, we present a flux
balance framework to aid in the design and interpretation of such experiments,
and discuss the challenges associated with measuring the relevant metabolic
fluxes. We then describe selected techniques that enable measurement of these
fluxes. Finally, we review prior experimental and theoretical work that has
employed techniques from biochemistry and nonequilibrium physics to determine
the energetic costs of cellular processes.


  In this paper we study intra-host viral adaptation by antigenic cooperation -
a mechanism of immune escape that serves as an alternative to the standard
mechanism of escape by continuous genomic diversification and allows to explain
a number of experimental observations associated with the establishment of
chronic infections by highly mutable viruses. Within this mechanism, the
topology of a cross-immunoreactivity network forces intra-host viral variants
to specialize for complementary roles and adapt to host's immune response as a
quasi-social ecosystem. Here we study dynamical changes in immune adaptation
caused by evolutionary and epidemiological events. First, we show that the
emergence of a viral variant with altered antigenic features may result in a
rapid re-arrangement of the viral ecosystem and a change in the roles played by
existing viral variants. In particular, it may push the population under immune
escape by genomic diversification towards the stable state of adaptation by
antigenic cooperation. Next, we study the effect of a viral transmission
between two chronically infected hosts, which results in merging of two
intra-host viral populations in the state of stable immune-adapted equilibrium.
In this case, we also describe how the newly formed viral population adapts to
the host's environment by changing the functions of its members. The results
are obtained analytically for minimal cross-immunoreactivity networks and
numerically for larger populations.


  The cellular composition of the tumor microenvironment can directly impact
cancer progression and the efficacy of therapeutics. Understanding immune cell
activity, the body's natural defense mechanism, in the vicinity of cancerous
cells is essential for developing beneficial treatments. Single cell RNA
sequencing (scRNA-seq) enables the examination of gene expression on an
individual cell basis, providing crucial information regarding both the
disturbances in cell functioning caused by cancer and cell-cell communication
in the tumor microenvironment. This novel technique generates large amounts of
data, which require proper processing. Various tools exist to facilitate this
processing but need to be organized to standardize the workflow from data
wrangling to visualization, cell type identification, and analysis of changes
in cellular activity, both from the standpoint of malignant cells and immune
stromal cells that eliminate them. We aimed to develop a standardized pipeline
(GeoTyper, https://github.com/celineyayifeng/GeoTyper) that integrates multiple
scRNA-seq tools for processing raw sequence data extracted from NCBI GEO,
visualization of results, statistical analysis, and cell type identification.
This pipeline leverages existing tools, such as Cellranger from 10X Genomics,
Alevin, and Seurat, to cluster cells and identify cell types based on gene
expression profiles. We successfully tested and validated the pipeline on
several publicly available scRNA-seq datasets, resulting in clusters
corresponding to distinct cell types. By determining the cell types and their
respective frequencies in the tumor microenvironment across multiple cancers,
this workflow will help quantify changes in gene expression related to
cell-cell communication and identify possible therapeutic targets.


  Most amino acids and sugars molecules occur in mirror, or chiral, images of
each other, knowns as enantiomers. However, life on Earth is mostly homochiral:
proteins contain almost exclusively L-amino acids, while only D-sugars appear
in RNA and DNA. The mechanism behind this fundamental asymmetry of life remains
unknown, despite much progress in the theoretical and experimental
understanding of homochirality in the past decades. We review three potential
mechanisms for the emergence of biological homochirality on primal Earth and
explore their implications for astrobiology: the first, that biological
homochirality is a stochastic process driven by local environmental
fluctuations; the second, that it is driven by circularly-polarized ultraviolet
radiation in star-forming regions; and the third, that it is driven by parity
violation at the elementary particle level. We argue that each of these
mechanisms leads to different observational consequences for the existence of
enantiomeric excesses in our solar system and in exoplanets, pointing to the
possibility that the search for life elsewhere will help elucidate the origins
of homochirality on Earth.


  Epidemics like Covid-19 and Ebola have impacted people's lives significantly.
The impact of mobility of people across the countries or states in the spread
of epidemics has been significant. The spread of disease due to factors local
to the population under consideration is termed the endogenous spread. The
spread due to external factors like migration, mobility, etc. is called the
exogenous spread. In this paper, we introduce the Exo-SIR model, an extension
of the popular SIR model and a few variants of the model. The novelty in our
model is that it captures both the exogenous and endogenous spread of the
virus. First, we present an analytical study. Second, we simulate the Exo-SIR
model with and without assuming contact network for the population. Third, we
implement the Exo-SIR model on real datasets regarding Covid-19 and Ebola. We
found that endogenous infection is influenced by exogenous infection.
Furthermore, we found that the Exo-SIR model predicts the peak time better than
the SIR model. Hence, the Exo-SIR model would be helpful for governments to
plan policy interventions at the time of a pandemic.


  Deep learning has significantly advanced and accelerated de novo molecular
generation. Generative networks, namely Variational Autoencoders (VAEs) can not
only randomly generate new molecules, but also alter molecular structures to
optimize specific chemical properties which are pivotal for drug-discovery.
While VAEs have been proposed and researched in the past for pharmaceutical
applications, they possess deficiencies which limit their ability to both
optimize properties and decode syntactically valid molecules. We present a
recurrent, conditional $\beta$-VAE which disentangles the latent space to
enhance post hoc molecule optimization. We create a mutual information driven
training protocol and data augmentations to both increase molecular validity
and promote longer sequence generation. We demonstrate the efficacy of our
framework on the ZINC-250k dataset, achieving SOTA unconstrained optimization
results on the penalized LogP (pLogP) and QED scores, while also matching
current SOTA results for validity, novelty and uniqueness scores for random
generation. We match the current SOTA on QED for top-3 molecules at 0.948,
while setting a new SOTA for pLogP optimization at 104.29, 90.12, 69.68 and
demonstrating improved results on the constrained optimization task.


  Plasmodium (P.) falciparum and P. vivax are the two most common causes of
malaria. While the majority of deaths and severe morbidity are due to P.
falciparum, P. vivax poses a greater challenge to eliminating malaria outside
of Africa due to its ability to form latent liver stage parasites
(hypnozoites), which can cause relapsing episodes within an individual patient.
In areas where P. falciparum and P. vivax are co-endemic, individuals can carry
parasites of both species simultaneously. These mixed infections complicate
dynamics in several ways; treatment of mixed infections will simultaneously
affect both species, P. falciparum can mask the detection of P. vivax, and it
has been hypothesised that clearing P. falciparum may trigger a relapse of
dormant P. vivax. When mixed infections are treated for only blood-stage
parasites, patients are at risk of relapse infections due to P. vivax
hypnozoites.
  We present a stochastic mathematical model that captures interactions between
P. falciparum and P. vivax, and incorporates both standard schizontocidal
treatment (which targets blood-stage parasites) and radical treatment (which
additionally targets liver-stage parasites). We apply this model to assess the
implications of different treatment coverage of radical cure for mixed and P.
vivax infections and a so-called "unified radical cure" treatment strategy for
P. falciparum, P. vivax and mixed infections. We find that a unified radical
cure strategy, with G6PD screening, leads to a substantially lower incidence of
malaria cases and deaths overall. We perform a one-way sensitivity analysis to
highlight important model parameters.


  In this paper we create a compartmental, stochastic process model of
SARS-CoV-2 transmission, where the process's mean and variance have distinct
dynamics. The model is fit to time series data from Washington from January
2020 to March 2021 using a deterministic, biologically-motivated signal
processing approach, and we show that the model's hidden states, like
population prevalence, agree with survey and other estimates. Then, in the
paper's second half, we demonstrate that the same model can be reframed as a
branching process with a dynamic degree distribution. This perspective allows
us to generate approximate transmission trees and estimate some higher order
statistics, like the clustering of cases as outbreaks, which we find to be
consistent with related observations from contact tracing and phylogenetics.


  Learning-to-rank, a machine learning technique widely used in information
retrieval, has recently been applied to the problem of ligand-based virtual
screening, to accelerate the early stages of new drug development. Ranking
prediction models learn based on ordinal relationships, making them suitable
for integrating assay data from various environments. Existing studies of rank
prediction in compound screening have generally used a learning-to-rank method
called RankSVM. However, they have not been compared with or validated against
the gradient boosting decision tree (GBDT)-based learning-to-rank methods that
have gained popularity recently. Furthermore, although the ranking metric
called Normalized Discounted Cumulative Gain (NDCG) is widely used in
information retrieval, it only determines whether the predictions are better
than those of other models. In other words, NDCG is incapable of recognizing
when a prediction model produces worse than random results. Nevertheless, NDCG
is still used in the performance evaluation of compound screening using
learning-to-rank. This study used the GBDT model with ranking loss functions,
called lambdarank and lambdaloss, for ligand-based virtual screening; results
were compared with existing RankSVM methods and GBDT models using regression.
We also proposed a new ranking metric, Normalized Enrichment Discounted
Cumulative Gain (NEDCG), which aims to properly evaluate the goodness of
ranking predictions. Results showed that the GBDT model with learning-to-rank
outperformed existing regression methods using GBDT and RankSVM on diverse
datasets. Moreover, NEDCG showed that predictions by regression were comparable
to random predictions in multi-assay, multi-family datasets, demonstrating its
usefulness for a more direct assessment of compound screening performance.


  Both natural ecosystems and biochemical reaction networks involve populations
of heterogeneous agents whose cooperative and competitive interactions lead to
a rich dynamics of species' abundances, albeit at vastly different scales. The
maintenance of diversity in large ecosystems is a longstanding puzzle, towards
which recent progress has been made by the derivation of dynamical mean-field
theories of random models. In particular, it has recently been shown that these
random models have a chaotic phase in which abundances display wild
fluctuations. When modest spatial structure is included, these fluctuations are
stabilized and diversity is maintained. If and how these phenomena have
parallels in biochemical reaction networks is currently unknown. Making this
connection is of interest since life requires cooperation among a large number
of molecular species. In this work, we find a reaction network whose
large-scale behavior recovers the random Lotka-Volterra model recently
considered in theoretical ecology. We clarify the assumptions necessary to
derive its large-scale description, and reveal the underlying assumptions made
on the noise to recover previous dynamical mean-field theories. Then, we show
how local detailed balance and the positivity of reaction rates, which are key
physical requirements of chemical reaction networks, provide obstructions
towards the construction of an associated dynamical mean-field theory of
biochemical reaction networks. Finally, we outline prospects and challenges for
the future.


  The algebraic properties of flattenings and subflattenings provide direct
methods for identifying edges in the true phylogeny -- and by extension the
complete tree -- using pattern counts from a sequence alignment. The relatively
small number of possible internal edges among a set of taxa (compared to the
number of binary trees) makes these methods attractive, however more could be
done to evaluate their effectiveness for inferring phylogenetic trees. This is
the case particularly for subflattenings, and our work makes progress in this
area. We introduce software for constructing and evaluating subflattenings for
splits, utilising a number of methods to make computing subflattenings more
tractable. We then present the results of simulations we have performed in
order to compare the effectiveness of subflattenings to that of flattenings in
terms of split score distributions, and susceptibility to possible biases. We
find that subflattenings perform similarly to flattenings in terms of the
distribution of split scores on the trees we examined, but may be less affected
by bias arising from both split size/balance and long branch attraction. These
insights are useful for developing effective algorithms to utilise these tools
for the purpose of inferring phylogenetic trees.


  Structure, composition and stability of ecological populations are shaped by
the inter- and intra-species interactions within these communities. It remains
to be fully understood how the interplay of these interactions with other
factors, such as immigration, control the structure, diversity and the long
term stability of ecological systems in the presence of noise and fluctuations.
We address this problem using a minimal model of interacting multi-species
ecological communities that incorporates competition, immigration and
demographic noise. We find that the complete phase diagram exhibits rich
behavior with multiple regimes that go beyond the classical 'niche' and
'neutral' regimes, extending and modifying the 'neutral-like' or 'niche-like'
dichotomy. In particular, we observe novel regimes that cannot be characterized
as either 'niche' or 'neutral' where a multimodal species abundance
distribution is observed. We characterize the transitions between the different
regimes and show how they arise from the underlying kinetics of the species
turnover, extinction and invasion. Our model serves as a minimal null model of
noisy competitive ecological systems, against which more complex models that
include factors such as mutations and environmental noise can be compared.


  Cell signaling networks are complex and often incompletely characterized,
making it difficult to obtain a comprehensive picture of the mechanisms they
encode. Mathematical modeling of these networks provides important clues, but
the models themselves are often complex, and it is not always clear how to
extract falsifiable predictions. Here we take an inverse approach, using
experimental data at the cell level to {deduce} the minimal signaling network.
We focus on cells' response to multiple cues, specifically on the surprising
case in which the response is antagonistic: the response to multiple cues is
weaker than the response to the individual cues. We systematically build
candidate signaling networks one node at a time, using the ubiquitous
ingredients of (i) up- or down-regulation, (ii) molecular conversion, or (iii)
reversible binding. In each case, our method reveals a minimal, interpretable
signaling mechanism that explains the antagonistic response. Our work provides
a systematic way to {deduce} molecular mechanisms from cell-level data.


  After more than 6 million deaths worldwide, the ongoing vaccination to
conquer the COVID-19 disease is now competing with the emergence of
increasingly contagious mutations, repeatedly supplanting earlier strains.
Following the near-absence of historical examples of the long-time evolution of
infectious diseases under similar circumstances, models are crucial to
exemplify possible scenarios. Accordingly, in the present work we
systematically generalize the popular susceptible-infected-recovered model to
account for mutations leading to repeatedly occurring new strains, which we
coarse grain based on tools from statistical mechanics to derive a model
predicting the most likely outcomes. The model predicts that mutations can
induce a super-exponential growth of infection numbers at early times, which
self-amplify to giant infection waves which are caused by a positive feedback
loop between infection numbers and mutations and lead to a simultaneous
infection of the majority of the population. At later stage -- if vaccination
progresses too slowly -- mutations can interrupt an ongoing decrease of
infection numbers and can cause infection revivals which occur as single waves
or even as whole wave trains featuring alternative periods of decreasing and
increasing infection numbers. This panorama of possible mutation-induced
scenarios should be tested in more detailed models to explore their concrete
significance for specific infectious diseases. Further, our results might be
useful for discussions regarding the importance of a release of vaccine-patents
to reduce the risk of mutation-induced infection revivals but also to
coordinate the release of measures following a downwards trend of infection
numbers.


  Machine learning (ML) has emerged as a pervasive tool in science,
engineering, and beyond. Its success has also led to several synergies with
molecular dynamics (MD) simulations, which we use to identify and characterize
the major metastable states of molecular systems. Typically, we aim to
determine the relative stabilities of these states and how rapidly they
interchange. This information allows mechanistic descriptions of molecular
mechanisms, enables a quantitative comparison with experiments, and facilitates
their rational design. ML impacts all aspects of MD simulations -- from
analyzing the data and accelerating sampling to defining more efficient or more
accurate simulation models.


  Predicting structural and energetic properties of a molecular system is one
of the fundamental tasks in molecular simulations, and it has use cases in
chemistry, biology, and medicine. In the past decade, the advent of machine
learning algorithms has impacted on molecular simulations for various tasks,
including property prediction of atomistic systems. In this paper, we propose a
novel methodology for transferring knowledge obtained from simple molecular
systems to a more complex one, possessing a significantly larger number of
atoms and degrees of freedom. In particular, we focus on the classification of
high and low free-energy states. Our approach relies on utilizing (i) a novel
hypergraph representation of molecules, encoding all relevant information for
characterizing the potential energy of a conformation, and (ii) novel message
passing and pooling layers for processing and making predictions on such
hypergraph-structured data. Despite the complexity of the problem, our results
show a remarkable AUC of 0.92 for transfer learning from tri-alanine to the
deca-alanine system. Moreover, we show that the very same transfer learning
approach can be used to group, in an unsupervised way, various secondary
structures of deca-alanine in clusters having similar free-energy values. Our
study represents a proof of concept that reliable transfer learning models for
molecular systems can be designed paving the way to unexplored routes in
prediction of structural and energetic properties of biologically relevant
systems.


  Many studies have shown that the protection of the host
$\it{Acyrthosiphon~pisum}$ (Hemiptera, Aphididae) against the parasitoid
$\it{Aphidius~ervi}$ (Hymenoptera, Braconidae) is conferred by the interaction
between the secondary endosymbiont $\it{Hamiltonella~defensa}$ and the
bacteriophage $\it{APSE}$ ($\it{Acyrthosiphon~pisum}$ secondary endosymbiont).
This interaction consists of the production of toxins by the endosymbiont's
molecular machinery, which is encoded by the inserted $\it{APSE}$ genes. The
toxins prevent the development of the parasitoid's egg, conferring protection
for the host. However, the effects of this microscopic interaction on
host-parasitoid dynamics are still an open question. We presented a new
mathematical model based on the bacteriophage effect on parasitism resistance.
We identified that the vertical transmission of the bacteriophage and the host
survival after the parasitoid attack are potential drivers of coexistence.
Also, we showed that the vertical transmission of $\it{H.~defensa}$ is
proportional to the time that the protected population became extinct. Our
results showed that the protected and unprotected hosts' survival after the
parasitoid attack is fundamental to understanding the equilibrium of long
host-parasitoid dynamics. Finally, we illustrated our model considering its
parameters based on experiments performed with $\it{A.~pisum}$ biotypes
$\it{Genista~tinctoria}$ and $\it{Medicago~sativa}$.


  Ontology Matching (OM) plays an important role in many domains such as
bioinformatics and the Semantic Web, and its research is becoming increasingly
popular, especially with the application of machine learning (ML) techniques.
Although the Ontology Alignment Evaluation Initiative (OAEI) represents an
impressive effort for the systematic evaluation of OM systems, it still suffers
from several limitations including limited evaluation of subsumption mappings,
suboptimal reference mappings, and limited support for the evaluation of
ML-based systems. To tackle these limitations, we introduce five new biomedical
OM tasks involving ontologies extracted from Mondo and UMLS. Each task includes
both equivalence and subsumption matching; the quality of reference mappings is
ensured by human curation, ontology pruning, etc.; and a comprehensive
evaluation framework is proposed to measure OM performance from various
perspectives for both ML-based and non-ML-based OM systems. We report
evaluation results for OM systems of different types to demonstrate the usage
of these resources, all of which are publicly available as part of the new
BioML track at OAEI 2022.


  Inference of species networks from genomic data under the Network
Multispecies Coalescent Model is currently severely limited by heavy
computational demands. It also remains unclear how complicated networks can be
for consistent inference to be possible. As a step toward inferring a general
species network, this work considers its tree of blobs, in which non-cut edges
are contracted to nodes, so only tree-like relationships between the taxa are
shown. An identifiability theorem, that most features of the unrooted tree of
blobs can be determined from the distribution of gene quartet topologies, is
established. This depends upon an analysis of gene quartet concordance factors
under the model, together with a new combinatorial inference rule. The
arguments for this theoretical result suggest a practical algorithm for tree of
blobs inference, to be fully developed in a subsequent work.


  Protein corona formation and nanoparticle aggregation have been heavily
discussed over the last years since the lack of fine-mapping of these two
combined effects has hindered the targeted delivery evolution and the
personalized nanomedicine development. We present a multi-technique approach
that combines Dynamic Light and Small-Angle X-ray Scattering techniques with
cryo-Transmission Electron Microscopy in a given fashion that efficiently
distinguishes protein corona from aggregates formation. This methodology was
tested using 25-nm model silica nanoparticles incubated with either model
proteins or biologically relevant proteomes (such as fetal bovine serum and
human plasma) in buffers of low and high ionic strengths to precisely tune
particle-to-protein interactions. In this work, we were able to differentiate
protein corona, small aggregates formation, and massive aggregation, as well as
obtain fractal information of the aggregates reliably and straightforwardly.
The strategy presented here can be expanded to other particle-to-protein
mixtures and might be employed as a quality control platform for samples that
undergo biological tests.


  Our paper presents three new classes of models: SIR-PH, SIR-PH-FA, and
SIR-PH-IA, and states two problems we would like to solve about them. Recall
that deterministic mathematical epidemiology has one basic general law, the R0
alternative" of [52, 51], which states that the local stability condition of
the disease free equilibrium may be expressed as R0 < 1, where R0 is the famous
basic reproduction number, which plays also a major role in the theory of
branching processes. The literature suggests that it is impossible to find
general laws concerning the endemic points. However, it is quite common that
  1. When R0 > 1, there exists a unique fixed endemic point, and
  2. the endemic point is locally stable when R0 > 1.
  One would like to establish these properties for a large class of realistic
epidemic models (and we do not include here epidemics without casualties). We
have introduced in [7, 5] a "simple", but broad class of "SIR-PH models" with
varying population, with the express purpose of establishing for these
processes the two properties above. Since that seemed still hard, we have
introduced a further class of "SIR-PH-FA" models, which may be interpreted as
approximations for the SIR-PH models, and which includes simpler models
typically studied in the literature (with constant population, without loss of
immunity, etc). The goal of our paper is to draw attention to the two open
problems above, for the SIR-PH, SIR-PH-FA, and also for a second, more refined
"intermediate approximation" SIR-PH-IA. We illustrate the current status-quo by
presenting new results on a generalization of the SAIRS epidemic model of [44,
40].


  Dynamic exploration for a predator-prey bio-system of two species with
ratio-dependent functional response is carried out, where the capability to
predate in both the stages of the predator, the juvenile and the matured, is
taken into account. But, only the matured predators are inferred to be
efficient in killing the prey without any negative repercussions. The mortality
risks for the juvenile predators are attributable to the inefficiency rate of
juveniles coupled with habitat complexity which is either in the form of
anti-predator behavior of the prey taken with the aid of their habitat or in
the form of a territorial generalist mesopredator. So as to avoid extinction of
either of the species and to preserve the food chain of the ecological system,
the results pertaining to the existence and stability of all the equilibrium
points of the bio-system along with permanence, transcritical and Hopf
bifurcation has been thoroughly studied. Corroboration of the results along
with the dependence of the biosystem on some crucial parameters is done through
numerical simulation. It is found that juvenile predators' inefficiency
relative to the resistance confronted, plays a crucial role to control each
species density of the ecosystem, as an intriguing limit cycle between the
trivial and axial equilibriums of the proposed system along with the
co-existing periodic point, because of some ineffeciency parametric value of
the juvenile predator has been witnessed.


  It is known that the trajectory of an endoreversibly driven system with
minimal dissipation is a geodesic on the equilibrium state space. Thereby, the
state space is equipped with the Riemannian metric given by the Hessian of the
free energy function, known as Fisher information metric. However, the
derivations given until now require both the system and the driving reservoir
to be in local equilibrium. In the present work, we rederive the framework for
chemical reaction networks and thereby enhance its scope of applicability to
the nonequilibrium situation. Moreover, because our results are derived without
restrictive assumptions, we are able to discuss phenomena that could not been
seen previously. We introduce a suitable weighted Fisher information metric on
the space of chemical concentrations and show that it characterizes the
dissipation caused by diffusive driving, with arbitrary diffusion rate
constants. This allows us to consider driving far from equilibrium. As the main
result, we show that the isometric embedding of a steady state manifold into
the concentration space yields a lower bound for the dissipation when the
system is driven along the manifold. We give an analytic expression for this
bound and for the corresponding geodesic, and thereby are able to dissect the
contributions from the driving kinetics and from thermodynamics. Finally, we
discuss in detail the application to quasi-thermostatic steady states.


  Deep learning is an important method for molecular design and exhibits
considerable ability to predict molecular properties, including
physicochemical, bioactive, and ADME/T (absorption, distribution, metabolism,
excretion, and toxicity) properties. In this study, we advanced a novel deep
learning architecture, termed FP-GNN, which combined and simultaneously learned
information from molecular graphs and fingerprints. To evaluate the FP-GNN
model, we conducted experiments on 13 public datasets, an unbiased LIT-PCBA
dataset, and 14 phenotypic screening datasets for breast cell lines. Extensive
evaluation results showed that compared to advanced deep learning and
conventional machine learning algorithms, the FP-GNN algorithm achieved
state-of-the-art performance on these datasets. In addition, we analyzed the
influence of different molecular fingerprints, and the effects of molecular
graphs and molecular fingerprints on the performance of the FP-GNN model.
Analysis of the anti-noise ability and interpretation ability also indicated
that FP-GNN was competitive in real-world situations.


  The automatic assignment of species information to the corresponding genes in
a research article is a critically important step in the gene normalization
task, whereby a gene mention is normalized and linked to a database record or
identifier by a text-mining algorithm. Existing methods typically rely on
heuristic rules based on gene and species co-occurrence in the article, but
their accuracy is suboptimal. We therefore developed a high-performance method,
using a novel deep learning-based framework, to classify whether there is a
relation between a gene and a species. Instead of the traditional binary
classification framework in which all possible pairs of genes and species in
the same article are evaluated, we treat the problem as a sequence-labeling
task such that only a fraction of the pairs needs to be considered. Our
benchmarking results show that our approach obtains significantly higher
performance compared to that of the rule-based baseline method for the species
assignment task (from 65.8% to 81.3% in accuracy). The source code and data for
species assignment are freely available at
https://github.com/ncbi/SpeciesAssignment.


  This paper considers a single-population model with age structure and
psychological effects in a polluted environment. We divide the single
population into two stages of larval and adult structure. The model uses
Logistic input, and the larvae are converted into adult bodies by constant
ratio. We only consider adulthood. The role of psychological effects makes the
contact between adult and environmental toxins a functional form, while the
contact between larvae and environmental toxins is linear.
  For the deterministic model embodied as a nonlinear time-varying system, we
discuss the asymptotic stability of the system by Lyapunov one-time
approximation theory, and give a sufficient condition for stability to be
established.
  Considering that the contact rate between biological and environmental toxins
in nature is not always constant, we make the contact rate interfere with white
noise, and then modify the contact rate into a stochastic process, thus
establishing a corresponding random single-population model. According to It\^o
formula and Lyapunov in the function method, we first prove the existence of
globally unique positive solutions for stochastic models under arbitrary
initial conditions, and then give sufficient conditions for weak average
long-term survival and random long-term survival for single populations in the
expected sense.


  In this article, we construct an age-structured model for COVID-19 with
vaccination and analyze it from multiple perspectives. We derive the unique
disease-free equilibrium point and the basic reproduction number $
\mathscr{R}_0 $, then we show that the disease-free equilibrium is locally
asymptotically stable when $ \mathscr{R}_0 < 1 $, while is unstable when $
\mathscr{R}_0 > 1 $. We also work out endemic equilibrium points and reveal the
stability. We use sensitivity analysis to explore how parameters influence $
\mathscr{R}_0 $. Sensitivity analysis helps us develop more targeted strategies
to control epidemics. Finally, this model is used to discuss the cases in
Shijiazhuang, Hebei Province at the beginning of 2021. We compare reported
cases with the simulation to evaluate the measures taken by Shijiazhuang
government. Our study shows how age structure, vaccination and drastic
containment measures can affect the epidemic.


  Deep generative modeling for biological sequences presents a unique challenge
in reconciling the bias-variance trade-off between explicit biological insight
and model flexibility. The deep manifold sampler was recently proposed as a
means to iteratively sample variable-length protein sequences by exploiting the
gradients from a function predictor. We introduce an alternative approach to
this guided sampling procedure, multi-segment preserving sampling, that enables
the direct inclusion of domain-specific knowledge by designating preserved and
non-preserved segments along the input sequence, thereby restricting variation
to only select regions. We present its effectiveness in the context of antibody
design by training two models: a deep manifold sampler and a GPT-2 language
model on nearly six million heavy chain sequences annotated with the IGHV1-18
gene. During sampling, we restrict variation to only the
complementarity-determining region 3 (CDR3) of the input. We obtain log
probability scores from a GPT-2 model for each sampled CDR3 and demonstrate
that multi-segment preserving sampling generates reasonable designs while
maintaining the desired, preserved regions.


  We investigate the large population dynamics of a family of stochastic
particle systems with three-state cyclic individual behaviour and
parameter-dependent transition rates. On short time scales, the dynamics turns
out to be approximated by an integrable Hamiltonian system whose phase space is
foliated by periodic trajectories. This feature suggests to consider the
effective dynamics of the long-term process that results from averaging over
the rapid oscillations. We establish the convergence of this process in the
large population limit to the solutions of an explicit stochastic differential
equation. Remarkably, this averaging phenomenon is complemented by the
convergence of stationary measures. The proof of averaging follows the
Stroock-Varadhan approach to martingale problems and relies on a fine analysis
of the system's dynamical features.


  We human beings show remarkable adaptability in response to complex
surroundings, we adopt different behavioral modes at different occasions, such
response multimodality is critical to our survival. Yet, how this behavioral
multimodality affects the evolution of cooperation remains largely unknown.
Here we build a toy model to address this issue by considering a population
with bimodal response behaviors, or specifically, with the Fermi and
Tit-for-tat updating rules. While the former rule tends to imitate the
strategies of those neighbors who are doing well, the latter repeats what their
neighbors did to them. In a structural mixing implementation, where the
updating rule is fixed for each individual, we find that a moderate mode
mixture unexpectedly boosts the overall cooperation level of the population.
The boost is even more pronounced in the probabilistic mixing, where each
individual randomly chooses one of the two modes at each step, and full
cooperation is seen in a wide range. These findings are robust to the
underlying topology of the population. Our mean-field treatment reveals that
the cooperation prevalence within the players with the Fermi rule linearly
increases with the fraction of TFT players and explains the non-monotonic
dependence in the structural mixing. Our study shows that the diversity in
response behaviors may help to explain the emergence of cooperation in
realistic contexts.


  The outbreak of COVID-19, beginning in 2019 and continuing through the time
of writing, has led to renewed interest in the mathematical modeling of
infectious disease. Recent works have focused on partial differential equation
(PDE) models, particularly reaction-diffusion models, able to describe the
progression of an epidemic in both space and time. These studies have shown
generally promising results in describing and predicting COVID-19 progression.
However, people often travel long distances in short periods of time, leading
to nonlocal transmission of the disease. Such contagion dynamics are not
well-represented by diffusion alone. In contrast, ordinary differential
equation (ODE) models may easily account for this behavior by considering
disparate regions as nodes in a network, with the edges defining nonlocal
transmission. In this work, we attempt to combine these modeling paradigms via
the introduction of a network structure within a reaction-diffusion PDE system.
This is achieved through the definition of a population-transfer operator,
which couples disjoint and potentially distant geographic regions, facilitating
nonlocal population movement between them. We provide analytical results
demonstrating that this operator does not disrupt the physical consistency or
mathematical well-posedness of the system, and verify these results through
numerical experiments. We then use this technique to simulate the COVID-19
epidemic in the Brazilian region of Rio de Janeiro, showcasing its ability to
capture important nonlocal behaviors, while maintaining the advantages of a
reaction-diffusion model for describing local dynamics.


  The tragedy of the commons (TOC) is a ubiquitous social dilemma witnessed in
interactions between a population of living entities and shared resources
available to them: The individuals in the population tend to selfishly
overexploit a common resource as it is arguably the rational choice, or in case
of non-human beings, it may be an evolutionarily uninvadable action. How to
avert the TOC is a significant problem related to the conservation of
resources. It is not hard to envisage situations where the resource could be
self-renewing and the size of the population may be dependent on the state of
the resource through the fractions of the population employing different
exploitation rates. If the self-renewal rate of the resource lies between the
maximum and the minimum exploitation rates, it is not a priori obvious under
what conditions the TOC can be averted. In this paper, we address this question
analytically and numerically using the setup of an evolutionary game
theoretical replicator equation that models the Darwinian tenet of natural
selection. Through the replicator equation, while we investigate how a
population of replicators exploit the shared resource, the latter's dynamical
feedback on the former is also not ignored. We also present a transparent
bottom-up derivation of the game-resource feedback model to facilitate future
studies on the stochastic effects on the findings presented herein.


  We collected individualized (unidentifiable) and aggregated openly available
data from various sources related to suspected/confirmed SARS-CoV-2 infections,
vaccinations, non-pharmaceutical government interventions, human mobility, and
levels of population inequality in Brazil. In addition, a data structure
allowing real-time data collection, curation, integration, and
extract-transform-load processes for different objectives was developed. The
granularity of this dataset (state- and municipality-wide) enables its
application to individualized and ecological epidemiological studies,
statistical, mathematical, and computational modeling, data visualization as
well as the scientific dissemination of information on the COVID-19 pandemic in
Brazil.


  Electrical detection of messenger ribonucleic acid (mRNA) is a promising
approach to enhancing transcriptomics and disease diagnostics because of its
sensitivity, rapidity, and modularity. Reported here is a fast SARS-CoV-2 mRNA
biosensor (<1 minute) with a limit of detection of 1 aM, and dynamic range of 4
orders of magnitude and a linear sensitivity of 22 mV per molar decade. These
figures of merit were obtained on photoresistlessly patterned monolayer
graphene field-effect transistors (FETs) derived from commercial four-inch
graphene on 90 nm of silicon dioxide on p-type silicon. Then, to facilitate
mRNA hybridization, graphene sensing mesa were coated with an ultrathin
sub-percolation threshold gold film for bonding 3'-thiolated single-stranded
deoxyribonucleic acid (ssDNA) probes complementary to SARS-CoV-2 nucleocapsid
phosphoprotein (N) gene. Sub-percolated gold was used to minimize the distance
between the graphene material and surface hybridization events. The
liquid-transfer characteristics of the graphene FETs repeatedly shows
correlation between the Dirac voltage and the copy number of polynucleotide.
Ultrathin percolated gold films on graphene FETs facilitate two-dimensional
electron gas (2DEG) mRNA biosensors for transcriptomic profiling.


  There is a crescent use of enzymes in multiple industries and sciences,
ranging from materials and fuel synthesis to pharmaceutical and food
production. Their applicability in this variety of fields depends not only on
their biochemical function but also on their physicochemical properties. In the
present work, we describe how the coincidence methodology can be employed to
construct similarity networks of seventy well-studied enzymes of the Glycoside
Hydrolase Family 13 and to identify communities of physicochemically related
enzymes. More specifically, each of the selected enzymes is mapped into a
network node, while the links between pairs of enzymes are determined by the
coincidence similarity between selected physicochemical features of interest.
The obtained networks have modularity and number of isolated nodes optimized
respectively to two parameters involved in the coincidence methodology,
resulting in highly modular networks. In order to investigate the effect of the
considered physicochemical features on the enzymes relationships, the
coincidence-based method also is applied to create a meta-network, in which the
enzymes similarity networks obtained by the combination of every possible
feature becomes nodes of a feature combination network, and the coincidence
similarity between those networks defines the respective links. The obtained
feature combination network systematically and comprehensively indicates the
impact of the selected physicochemical features on enzyme similarity. Several
interesting results are reported and discussed, including the identification of
subgroups of enzymes with similar physicochemical features within catalytical
classes, providing important information for the selection and design of
enzymes for targeted biotechnological applications.


  In this article, we investigate the importance of demographic and contact
patterns in determining the spread of COVID-19 and to the effectiveness of
social distancing policies. We investigate these questions proposing an
augmented epidemiological model with an age-structured model, with the
population divided into susceptible (S), exposed (E), infected and asymptomatic
(A), hospitalized (H), infected and symptomatic (I), and recovered individuals
(R), to simulate COVID-19 dissemination. The simulations were carried out using
six combinations of four types of isolation policies (work restrictions,
isolation of the elderly, community distancing and school closures) and four
representative fictitious countries generated over alternative demographic
transition stage patterns (aged developed, developed, developing and least
developed countries). We concluded that the basic reproduction number depends
on the age profile and the contact patterns. The aged developed country had the
lowest basic reproduction number ($R0=1.74$) due to the low contact rate among
individuals, followed by the least developed country ($R0=2.00$), the
developing country ($R0=2.43$) and the developed country ($R0=2.64$). Because
of these differences in the basic reproduction numbers, the same intervention
policies had higher efficiencies in the aged and least developed countries. Of
all intervention policies, the reduction in work contacts and community
distancing were the ones which produced the highest decrease in the $R0$ value,
prevalence, maximum hospitalization demand and fatality rate. The isolation of
the elderly was more effective in the developed and aged developed countries.
The school closure was the less effective intervention policy, though its
effects were not negligible in the least developed and developing countries.


  Recent studies have shown that packings of cells, both eukaryotic cellular
tissues and growing or swarming bacterial colonies, can often be understood as
active nematic fluids. A key property of volume-conserving active nematic model
systems is chaotic self-mixing characterized by motile topological defects.
However, for active nematics driven by growth rather than motility, less is
understood about mixing and defect motion. Mixing could affect evolutionary
outcomes in bacterial colonies by counteracting the tendency to spatially
segregate into monoclonal sectors, which reduces the local genetic diversity
and confines competition between subpopulations to the boundaries between
neighboring sectors. To examine whether growth-driven active nematic physics
could influence this genetic demixing process, we conduct agent-based
simulations of growing, dividing, and sterically repelling rod-like bacteria of
various aspect ratios, and we analyze colony morphology using tools from both
soft matter physics and population genetics. We find that despite measurable
defect self-propulsion in growth-driven active nematics, the radial expansion
flow prevents chaotic mixing. Even so, at biologically relevant cell aspect
ratios, self-mixing is more effective in growing active nematics of rod-like
cells compared to growing isotropic colonies of round cells. This suggests
potential evolutionary consequences associated with active nematic dynamics.


  Density dependence is important in the ecology and evolution of microbial and
cancer cells. Typically, we can only measure net growth rates, but the
underlying density-dependent mechanisms that give rise to the observed dynamics
can manifest in birth processes, death processes, or both. Therefore, we
utilize the mean and variance of cell number fluctuations to separately
identify birth and death rates from time series that follow stochastic
birth-death processes with logistic growth. Our method provides a novel
perspective on stochastic parameter identifiability, which we validate by
analyzing the accuracy in terms of the discretization bin size. We apply our
method to the scenario where a homogeneous cell population goes through three
stages: (1) grows naturally to its carrying capacity, (2) is treated with a
drug that reduces its carrying capacity, and (3) overcomes the drug effect to
restore its original carrying capacity. In each stage, we disambiguate whether
it happens through the birth process, death process, or some combination of the
two, which contributes to understanding drug resistance mechanisms. In the case
of limited data sets, we provide an alternative method based on maximum
likelihood and solve a constrained nonlinear optimization problem to identify
the most likely density dependence parameter for a given cell number time
series. Our methods can be applied to other biological systems at different
scales to disambiguate density-dependent mechanisms underlying the same net
growth rate.


  Although lessons have been learned from previous severe acute respiratory
syndrome (SARS) and Middle East respiratory syndrome (MERS) outbreaks, the
rapid evolution of the viruses means that future outbreaks of a much larger
scale are possible, as shown by the current coronavirus disease 2019 (COVID-19)
outbreak. Therefore, it is necessary to better understand the evolution of
coronaviruses as well as viruses in general. This study reports a comparative
analysis of the amino acid usage within several key viral families and genera
that are prone to triggering outbreaks, including coronavirus (SARS-CoV-2,
SARS-CoV, MERS-CoV, HCoV-HKU1, HCoV-OC43, HCoV-NL63, HCoV-229E), influenza A
(H1N1, H3N2), flavivirus (dengue virus serotypes 1-4, Zika) and ebolavirus
(Zaire, Sudan, Bundibugyo ebolavirus). Our analysis reveals that the
distribution of amino acid usage in the viral genome is constrained to follow a
linear order, and the distribution remains closely related to the viral species
within the family or genus. This constraint can be adapted to predict viral
mutations and future variants of concern. By studying previous SARS and MERS
outbreaks, we have adapted this naturally occurring pattern to determine that
although pangolin plays a role in the outbreak of COVID-19, it may not be the
sole agent as an intermediate animal. In addition to this study, our findings
contribute to the understanding of viral mutations for subsequent development
of vaccines and toward developing a model to determine the source of the
outbreak.


  Under certain conditions, the dynamics of coarse-grained models of solvated
proteins can be described using a Markov state model, which tracks the
evolution of populations of configurations. The transition rates among states
that appear in the Markov model can be determined by computing the relative
entropy of states and their mean first passage times. In this paper, we present
an adaptive method to evaluate the configurational entropy and the mean first
passage times for linear chain models with discontinuous potentials. The
approach is based on event-driven dynamical sampling in a massively parallel
architecture. Using the fact that the transition rate matrix can be calculated
for any choice of interaction energies at any temperature, it is demonstrated
how each state's energy can be chosen such that the average time to transition
between any two states is minimized. The methods are used to analyze the
optimization of the folding process of two protein systems: the crambin
protein, and a model with frustration and misfolding. It is shown that the
folding pathways for both systems are comprised of two regimes: first, the
rapid establishment of local bonds, followed by the subsequent formation of
more distant contacts. The state energies that lead to the most rapid folding
encourage multiple pathways, and either penalize folding pathways through
kinetic traps by raising the energies of trapping states, or establish an
escape route from the trapping states by lowering free energy barriers to other
states that rapidly reach the native state.


  A critical step of genome sequence analysis is the mapping of sequenced DNA
fragments (i.e., reads) collected from an individual to a known linear
reference genome sequence (i.e., sequence-to-sequence mapping). Recent works
replace the linear reference sequence with a graph-based representation of the
reference genome, which captures the genetic variations and diversity across
many individuals in a population. Mapping reads to the graph-based reference
genome (i.e., sequence-to-graph mapping) results in notable quality
improvements in genome analysis. Unfortunately, while sequence-to-sequence
mapping is well studied with many available tools and accelerators,
sequence-to-graph mapping is a more difficult computational problem, with a
much smaller number of practical software tools currently available.
  We analyze two state-of-the-art sequence-to-graph mapping tools and reveal
four key issues. We find that there is a pressing need to have a specialized,
high-performance, scalable, and low-cost algorithm/hardware co-design that
alleviates bottlenecks in both the seeding and alignment steps of
sequence-to-graph mapping.
  To this end, we propose SeGraM, a universal algorithm/hardware co-designed
genomic mapping accelerator that can effectively and efficiently support both
sequence-to-graph mapping and sequence-to-sequence mapping, for both short and
long reads. To our knowledge, SeGraM is the first algorithm/hardware co-design
for accelerating sequence-to-graph mapping. SeGraM consists of two main
components: (1) MinSeed, the first minimizer-based seeding accelerator; and (2)
BitAlign, the first bitvector-based sequence-to-graph alignment accelerator.
  We demonstrate that SeGraM provides significant improvements for multiple
steps of the sequence-to-graph and sequence-to-sequence mapping pipelines.


  The Critical Assessment of Genome Interpretation (CAGI) aims to advance the
state of the art for computational prediction of genetic variant impact,
particularly those relevant to disease. The five complete editions of the CAGI
community experiment comprised 50 challenges, in which participants made blind
predictions of phenotypes from genetic data, and these were evaluated by
independent assessors. Overall, results show that while current methods are
imperfect, they have major utility for research and clinical applications.
Missense variant interpretation methods are able to estimate biochemical
effects with increasing accuracy. Performance is particularly strong for
clinical pathogenic variants, including some difficult-to-diagnose cases, and
extends to interpretation of cancer-related variants. Assessment of methods for
regulatory variants and complex trait disease risk is less definitive, and
indicates performance potentially suitable for auxiliary use in the clinic.
Emerging methods and increasingly large, robust datasets for training and
assessment promise further progress ahead.


  Analysis of bond rupture data from single-molecule force spectroscopy
experiments commonly relies on the strong assumption that the bond dissociation
process is irreversible. However, with increased spatiotemporal resolution of
instruments it is now possible to observe multiple unbinding-rebinding events
in a single pulling experiment. Here, we augment the theory of force-induced
unbinding by explicitly taking into account rebinding kinetics, and provide
approximate analytic solutions of the resulting rate equations. Furthermore, we
use a short-time expansion of the exact kinetics to construct numerically
efficient maximum likelihood estimators for the parameters of the
force-dependent unbinding and rebinding rates, which pair well with and
complement established methods, such as the analysis of rate maps. We provide
an open-source implementation of the theory, evaluated for Bell-like rates,
which we apply to synthetic data generated by a Gillespie stochastic simulation
algorithm for time-dependent rates.


  Quantum annealing is a promising approach for obtaining good approximate
solutions to difficult optimization problems. Folding a protein sequence into
its minimum-energy structure represents such a problem. For testing new
algorithms and technologies for this task, the minimal lattice-based HP model
is well suited, as it represents a considerable challenge despite its
simplicity. The HP model has favorable interactions between adjacent, not
directly bound hydrophobic residues. Here, we develop a novel spin
representation for lattice protein folding tailored for quantum annealing. With
a distributed encoding onto the lattice, it differs from earlier attempts to
fold lattice proteins on quantum annealers, which were based upon chain growth
techniques. With our encoding, the Hamiltonian by design has the quadratic
structure required for calculations on an Ising-type annealer, without having
to introduce any auxiliary spin variables. This property greatly facilitates
the study of long chains. The approach is robust to changes in the parameters
required to constrain the spin system to chain-like configurations, and
performs very well in terms of solution quality. The results are evaluated
against existing exact results for HP chains with up to $N=30$ beads with 100%
hit rate, thereby also outperforming classical simulated annealing. In
addition, the method allows us to recover the lowest known energies for $N=48$
and $N=64$ HP chains, with similar hit rates. These results are obtained by the
commonly used hybrid quantum-classical approach. For pure quantum annealing,
our method successfully folds an $N=14$ HP chain. The calculations were
performed on a D-Wave Advantage quantum annealer.


  Biochemical reactions inside living cells often occur in the presence of
crowders -- molecules that do not participate in the reactions but influence
the reaction rates through excluded volume effects. However the standard
approach to modelling stochastic intracellular reaction kinetics is based on
the chemical master equation (CME) whose propensities are derived assuming no
crowding effects. Here, we propose a machine learning strategy based on
Bayesian Optimisation utilising synthetic data obtained from spatial cellular
automata (CA) simulations (that explicitly model volume-exclusion effects) to
learn effective propensity functions for CMEs. The predictions from a small CA
training data set can then be extended to the whole range of parameter space
describing physiologically relevant levels of crowding by means of Gaussian
Process regression. We demonstrate the method on an enzyme-catalyzed reaction
and a genetic feedback loop, showing good agreement between the time-dependent
distributions of molecule numbers predicted by the effective CME and CA
simulations.


  Can a micron sized sack of interacting molecules understand, and adapt to a
constantly-fluctuating environment? Cellular life provides an existence proof
in the affirmative, but the principles that allow for life's existence are far
from being proven. One challenge in engineering and understanding biochemical
computation is the intrinsic noise due to chemical fluctuations. In this paper,
we draw insights from machine learning theory, chemical reaction network
theory, and statistical physics to show that the broad and biologically
relevant class of detailed balanced chemical reaction networks is capable of
representing and conditioning complex distributions. These results illustrate
how a biochemical computer can use intrinsic chemical noise to perform complex
computations. Furthermore, we use our explicit physical model to derive
thermodynamic costs of inference.


  We present SODA, a lightweight and open-source visualization library for
biological sequence annotations that enables straightforward development of
flexible, dynamic, and interactive web graphics. SODA is implemented in
TypeScript and can be used as a library within TypeScript and JavaScript.


  In this work we derive new analytic expressions for fixation time in
Wright-Fisher model with selection. The three standard cases for fixation are
considered: fixation to zero, to one or both. Second order differential
equations for fixation time are obtained by a simplified approach using only
the law of total probability and Taylor expansions. The obtained solutions are
given by a combination of exponential integral functions with elementary
functions. We then state approximate formulas involving only elementary
functions valid for small selection effects. The quality of our results are
explored throughout an extensive simulation study. We show that our results
approximate the discrete problem very accurately even for small population size
(a few hundreds) and large selection coefficients.


  Mutating variants of COVID-19 have been reported across many US states since
2021. In the fight against COVID-19, it has become imperative to study the
heterogeneity in the time-varying transmission rates for each variant in the
presence of pharmaceutical and non-pharmaceutical mitigation measures. We
develop a Susceptible-Exposed-Infected-Recovered mathematical model to
highlight the differences in the transmission of the B.1.617.2 delta variant
and the original SARS-CoV-2. Theoretical results for the well-posedness of the
model are discussed. A Deep neural network is utilized and a deep learning
algorithm is developed to learn the time-varying heterogeneous transmission
rates for each variant. The accuracy of the algorithm for the model is shown
using error metrics in the data-driven simulation for COVID-19 variants in the
US states of Florida, Alabama, Tennessee, and Missouri. Short-term forecasting
of daily cases is demonstrated using long short term memory neural network and
an adaptive neuro-fuzzy inference system.


  The conditional intensity (CI) of a counting process $Y_t$ is based on the
minimal knowledge $\mathcal{F}_t^Y$, i.e., on the observation of $Y_t$ alone.
Prominently, the mutual information rate of a signal and its Poisson channel
output is a difference functional between the CI and the intensity that has
full knowledge about the input. While the CI of Markov-modulated Poisson
processes evolves according to Snyder's filter, self-exciting processes, e.g.,
Hawkes processes, specify the CI via the history of $Y_t$. The emergence of the
CI as a self-contained stochastic process prompts us to bring its statistical
ensemble into focus. We investigate the asymptotic conditional intensity
distribution (ACID) and emphasize its rich information content. We assume the
case in which the CI is determined from a sufficient statistic that progresses
as a Markov process. We present a simulation-free method to compute the ACID
when the dimension of the sufficient statistic is low. The method is made
possible by introducing a backward recurrence time parametrization, which has
the advantage to align all probability inflow in a boundary condition for the
master equation. Case studies illustrate the usage of ACID for three primary
examples: 1) the Poisson channels with binary Markovian input (as an example of
a Markov-modulated Poisson process), 2) the standard Hawkes process with
exponential kernel (as an example of a self-exciting counting process) and 3)
the Gamma filter (as an example of an approximate filter to a Markov-modulated
Poisson process).


  Deep generative models have achieved tremendous success in designing novel
drug molecules in recent years. A new thread of works have shown the great
potential in advancing the specificity and success rate of in silico drug
design by considering the structure of protein pockets. This setting posts
fundamental computational challenges in sampling new chemical compounds that
could satisfy multiple geometrical constraints imposed by pockets. Previous
sampling algorithms either sample in the graph space or only consider the 3D
coordinates of atoms while ignoring other detailed chemical structures such as
bond types and functional groups. To address the challenge, we develop
Pocket2Mol, an E(3)-equivariant generative network composed of two modules: 1)
a new graph neural network capturing both spatial and bonding relationships
between atoms of the binding pockets and 2) a new efficient algorithm which
samples new drug candidates conditioned on the pocket representations from a
tractable distribution without relying on MCMC. Experimental results
demonstrate that molecules sampled from Pocket2Mol achieve significantly better
binding affinity and other drug properties such as druglikeness and synthetic
accessibility.


  Recent evidences show that individuals who recovered from COVID-19 can be
reinfected. However, this phenomenon has rarely been studied using mathematical
models. In this paper, we propose a SEIRE epidemic model to describe the spread
of the epidemic with reinfection. We obtain the important thresholds $R_0$ (the
basic reproduction number) and Rc (a threshold less than one). Our
investigations show that when $R_0 > 1$, the system has an endemic equilibrium,
which is globally asymptotically stable. When $R_c < R_0 < 1$, the epidemic
system exhibits bistable dynamics. That is, the system has backward bifurcation
and the disease cannot be eradicated. In order to eradicate the disease, we
must ensure that the basic reproduction number $R_0$ is less than $R_c$. The
basic reinfection number is obtained to measure the reinfection force, which
turns out to be a new tipping point for disease dynamics. We also give
definition of robustness, a new concept to measure the difficulty of completely
eliminating the disease for a bistable epidemic system. Numerical simulations
are carried out to verify the conclusions.


  Deep learning has achieved tremendous success in designing novel chemical
compounds with desirable pharmaceutical properties. In this work, we focus on a
new type of drug design problem -- generating a small "linker" to physically
attach two independent molecules with their distinct functions. The main
computational challenges include: 1) the generation of linkers is conditional
on the two given molecules, in contrast to generating full molecules from
scratch in previous works; 2) linkers heavily depend on the anchor atoms of the
two molecules to be connected, which are not known beforehand; 3) 3D structures
and orientations of the molecules need to be considered to avoid atom clashes,
for which equivariance to E(3) group are necessary. To address these problems,
we propose a conditional generative model, named 3DLinker, which is able to
predict anchor atoms and jointly generate linker graphs and their 3D structures
based on an E(3) equivariant graph variational autoencoder. So far as we know,
there are no previous models that could achieve this task. We compare our model
with multiple conditional generative models modified from other molecular
design tasks and find that our model has a significantly higher rate in
recovering molecular graphs, and more importantly, accurately predicting the 3D
coordinates of all the atoms.


  We discuss a method to describe the qualitative dynamics of chemical reaction
networks in terms of symbolic dynamics. The method, that can be applied to
mass-action reaction networks with separated timescales, uses solutions of the
partial tropical equilibration problem as proxies for symbolic states. The
partial tropical equilibration solutions are found algorithmically. These
solutions also provide the scaling needed for slow-fast decomposition and model
reduction. Any trace of the model can thus be represented as a sequence of
local approximations of the full model. We illustrate the method using as case
study a biochemical model of the cell cycle.


  Understanding the process of protein unfolding plays a crucial role in
various applications such as design of folding-based protein engines. Using the
well-established kinetostatic compliance (KCM)-based method for modeling of
protein conformation dynamics and a recent nonlinear control theoretic approach
to KCM-based protein folding, this paper formulates protein unfolding as a
destabilizing control analysis/synthesis problem. In light of this formulation,
it is shown that the Chetaev instability framework can be used to investigate
the KCM-based unfolding dynamics. In particular, a Chetaev function for
analysis of unfolding dynamics under the effect of optical tweezers and a class
of control Chetaev functions for synthesizing control inputs that elongate
protein strands from their folded conformations are presented. Based on the
presented control Chetaev function, an unfolding input is derived from the
Artstein-Sontag universal formula and the results are compared against optical
tweezer-based unfolding.


  Despite the availability of effective vaccines, the persistence of SARS-CoV-2
suggests that co-circulation with other pathogens and resulting multi-epidemics
may become increasingly frequent. To better forecast and control the risk of
such multi-epidemics, it is essential to elucidate the potential interactions
of SARS-CoV-2 with other pathogens; these interactions, however, remain poorly
defined. Here, we aimed to review the current body of evidence about SARS-CoV-2
interactions. To study pathogen interactions in a systematic way, we first
developed a general framework to capture their major components: sign,
strength, symmetry, duration, and mechanism. We then reviewed the experimental
evidence from animal models about SARS-CoV-2 interactions. Of the 14 studies
identified, 11 focused on the outcomes of co-infection with non-attenuated
influenza A viruses and generally demonstrated that co-infection increased
disease severity compared with either mono-infection. By contrast, the effect
of co-infection on the viral load of either virus was variable and inconsistent
across studies. Next, we reviewed the epidemiological evidence about SARS-CoV-2
interactions in human populations. Although numerous studies were identified,
only few were specifically designed to infer interaction and many were prone to
multiple biases, including confounding. Nevertheless, their results suggested
that influenza and pneumococcal conjugate vaccinations were associated with
reduced riskof SARS-CoV-2 infection. Finally, we formulated simple transmission
models of SARS-CoV-2 co-circulation with a viral or a bacterial pathogen,
showing how they can naturally incorporate the proposed framework. More
generally, we argue that such models, when designed with an integrative and
multidisciplinary perspective, will be invaluable tools to resolve the
substantial uncertainties that remain about SARS-CoV-2 interactions.


  Transformer models have been developed in molecular science with excellent
performance in applications including quantitative structure-activity
relationship (QSAR) and virtual screening (VS). Compared with other types of
models, however, they are large, which results in a high hardware requirement
to abridge time for both training and inference processes. In this work,
cross-layer parameter sharing (CLPS), and knowledge distillation (KD) are used
to reduce the sizes of transformers in molecular science. Both methods not only
have competitive QSAR predictive performance as compared to the original BERT
model, but also are more parameter efficient. Furthermore, by integrating CLPS
and KD into a two-state chemical network, we introduce a new deep lite chemical
transformer model, DeLiCaTe. DeLiCaTe captures general-domains as well as
task-specific knowledge, which lead to a 4x faster rate of both training and
inference due to a 10- and 3-times reduction of the number of parameters and
layers, respectively. Meanwhile, it achieves comparable performance in QSAR and
VS modeling. Moreover, we anticipate that the model compression strategy
provides a pathway to the creation of effective generative transformer models
for organic drug and material design.


  In more than two years since the COVID-19 virus was first detected in China,
hundreds of millions of individuals have been infected, and millions have died.
Aside from the immediate need for medical solutions (such as vaccines and
medications) to treat the epidemic, the Corona pandemic has strengthened the
demand for mathematical models that can predict the spread of the pandemic in
an ever-changing reality. Here, we present a novel, dynamic particle model
based on the basic principles of statistical physics that enables the
prediction of the spreading of Covid-19 in the presence of effective vaccines.
This particle model enables us to accurately examine the effects of the vaccine
on different subgroups of the vaccinated population and the entire population
and to identify the vaccine waning. Furthermore, a particle model can predict
the prevalence of two competing variants over time and their associated
morbidity.


  The rapid growth in the number of experimental and predicted protein
structures and more complicated protein structures challenge users in
computational biology for utilizing the structural information and protein
surface property representation. Recently, AlphaFold2 released the
comprehensive proteome of various species, and protein surface property
representation plays a crucial role in protein-molecule interaction prediction
such as protein-protein interaction, protein-nucleic acid interaction, and
protein-compound interaction. Here, we proposed the first comprehensive
database, namely ProNet DB, which incorporates multiple protein surface
representations and RNA-binding landscape for more than 326,175 protein
structures covering 16 model organism proteomes from AlphaFold Protein
Structure Database (AlphaFold DB) and experimentally validated protein
structures deposited in Protein Data Bank (PDB). For each protein, we provided
the original protein structure, surface property representation including
hydrophobicity, charge distribution, hydrogen bond, interacting face, and
RNA-binding landscape such as RNA binding sites and RNA binding preference. To
interpret protein surface property representation and RNA binding landscape
intuitively, we also integrate Mol* and Online 3D Viewer to visualize the
representation on the protein surface. The pre-computed features are available
for the users instantaneously and boost computational biology development
including molecular mechanism exploration, geometry-based drug discovery and
novel therapeutics development. The server is now available on
https://proj.cse.cuhk.edu.hk/aihlab/pronet/.


  We now need more than ever to make genome analysis more intelligent. We need
to read, analyze, and interpret our genomes not only quickly, but also
accurately and efficiently enough to scale the analysis to population level.
There currently exist major computational bottlenecks and inefficiencies
throughout the entire genome analysis pipeline, because state-of-the-art genome
sequencing technologies are still not able to read a genome in its entirety. We
describe the ongoing journey in significantly improving the performance,
accuracy, and efficiency of genome analysis using intelligent algorithms and
hardware architectures. We explain state-of-the-art algorithmic methods and
hardware-based acceleration approaches for each step of the genome analysis
pipeline and provide experimental evaluations. Algorithmic approaches exploit
the structure of the genome as well as the structure of the underlying
hardware. Hardware-based acceleration approaches exploit specialized
microarchitectures or various execution paradigms (e.g., processing inside or
near memory) along with algorithmic changes, leading to new hardware/software
co-designed systems. We conclude with a foreshadowing of future challenges,
benefits, and research directions triggered by the development of both very low
cost yet highly error prone new sequencing technologies and specialized
hardware chips for genomics. We hope that these efforts and the challenges we
discuss provide a foundation for future work in making genome analysis more
intelligent. The analysis script and data used in our experimental evaluation
are available at: https://github.com/CMU-SAFARI/Molecules2Variations