  Contact matrices are an important ingredient in age-structured epidemic
models to inform the simulated spread of the disease between sub-groups of the
population. These matrices are generally derived using resource-intensive
diary-based surveys and few exist in the Global South or tailored to vulnerable
populations. In particular, no contact matrices exist for refugee settlements -
locations under-served by epidemic models in general. In this paper we present
a novel, mixed-method approach, for deriving contact matrices in populations
which combines a lightweight, rapidly deployable, survey with an agent-based
model of the population informed by census and behavioural data. We use this
method to derive the first set of contact matrices for the Cox's Bazar refugee
settlement in Bangladesh. The matrices from the refugee settlement show strong
banding effects due to different age cut-offs in attendance at certain venues,
such as distribution centres and religious sites, as well as the important
contribution of the demographic profile of the settlement which was encoded in
the model. These can have significant implications to the modelled disease
dynamics. To validate our approach, we also apply our method to the population
of the UK and compare our derived matrices against well-known contact matrices
previously collected using traditional approaches. Overall, our findings
demonstrate that our mixed-method approach can address some of the challenges
of both the traditional and previously proposed agent-based approaches to
deriving contact matrices, and has the potential to be rolled-out in other
resource-constrained environments. This work therefore contributes to a broader
aim of developing new methods and mechanisms of data collection for modelling
disease spread in refugee and IDP settlements and better serving these
vulnerable communities.


  De novo molecule generation can suffer from data inefficiency; requiring
large amounts of training data or many sampled data points to conduct objective
optimization. The latter is a particular disadvantage when combining deep
generative models with computationally expensive molecule scoring functions
(a.k.a. oracles) commonly used in computer-aided drug design. Recent works have
therefore focused on methods to improve sample efficiency in the context of de
novo molecule drug design, or to benchmark it. In this work, we discuss and
adapt a recent sample efficiency benchmark to better reflect realistic goals
also with respect to the quality of chemistry generated, which must always be
considered in the context of small-molecule drug design; we then re-evaluate
all benchmarked generative models. We find that accounting for molecular weight
and LogP with respect to the training data, and the diversity of chemistry
proposed, re-orders the ranking of generative models. In addition, we benchmark
a recently proposed method to improve sample efficiency (Augmented Hill-Climb)
and found it ranked top when considering both the sample efficiency and
chemistry of molecules generated. Continual improvements in sample efficiency
and chemical desirability enable more routine integration of computationally
expensive scoring functions on a more realistic timescale.


  Motivation: Enhancers are important cis-regulatory elements that regulate a
wide range of biological functions and enhance the transcription of target
genes. Although many feature extraction methods have been proposed to improve
the performance of enhancer identification, they cannot learn position-related
multiscale contextual information from raw DNA sequences.
  Results: In this article, we propose a novel enhancer identification method
(iEnhancer-ELM) based on BERT-like enhancer language models. iEnhancer-ELM
tokenizes DNA sequences with multi-scale k-mers and extracts contextual
information of different scale k-mers related with their positions via an
multi-head attention mechanism. We first evaluate the performance of different
scale k-mers, then ensemble them to improve the performance of enhancer
identification. The experimental results on two popular benchmark datasets show
that our model outperforms stateof-the-art methods. We further illustrate the
interpretability of iEnhancer-ELM. For a case study, we discover 30 enhancer
motifs via a 3-mer-based model, where 12 of motifs are verified by STREME and
JASPAR, demonstrating our model has a potential ability to unveil the
biological mechanism of enhancer.
  Availability and implementation: The models and associated code are available
at https://github.com/chen-bioinfo/iEnhancer-ELM
  Contact: junjiechen@hit.edu.cn
  Supplementary information: Supplementary data are available at Bioinformatics
Advances online.


  Drug addiction or drug overdose is a global public health crisis, and the
design of anti-addiction drugs remains a major challenge due to intricate
mechanisms. Since experimental drug screening and optimization are too
time-consuming and expensive, there is urgent need to develop innovative
artificial intelligence (AI) methods for addressing the challenge. We tackle
this challenge by topology-inferred drug addiction learning (TIDAL) built from
integrating topological Laplacian, deep bidirectional transformer, and
ensemble-assisted neural networks (EANNs). The topological Laplacian is a novel
algebraic topology tool that embeds molecular topological invariants and
algebraic invariants into its harmonic spectra and non-harmonic spectra,
respectively. These invariants complement sequence information extracted from a
bidirectional transformer. We validate the proposed TIDAL framework on 22 drug
addiction related, 4 hERG, and 12 DAT datasets, showing that TIDAL is a
state-of-the-art framework for the modeling and analysis of drug addiction
data. We carry out cross-target analysis of the current drug addiction
candidates to alert their side effects and identify their repurposing
potentials, revealing drugmediated linear and bilinear target correlations.
Finally, TIDAL is applied to shed light on relative efficacy, repurposing
potential, and potential side effects of 12 existing anti-addiction
medications. Our results suggest that TIDAL provides a new computational
strategy for pressingly-needed anti-substance addiction drug development.


  In this work, we propose MEDICO, a Multi-viEw Deep generative model for
molecule generation, structural optimization, and the SARS-CoV-2 Inhibitor
disCOvery. To the best of our knowledge, MEDICO is the first-of-this-kind graph
generative model that can generate molecular graphs similar to the structure of
targeted molecules, with a multi-view representation learning framework to
sufficiently and adaptively learn comprehensive structural semantics from
targeted molecular topology and geometry. We show that our MEDICO significantly
outperforms the state-of-the-art methods in generating valid, unique, and novel
molecules under benchmarking comparisons. In particular, we showcase the
multi-view deep learning model enables us to generate not only the molecules
structurally similar to the targeted molecules but also the molecules with
desired chemical properties, demonstrating the strong capability of our model
in exploring the chemical space deeply. Moreover, case study results on
targeted molecule generation for the SARS-CoV-2 main protease (Mpro) show that
by integrating molecule docking into our model as chemical priori, we
successfully generate new small molecules with desired drug-like properties for
the Mpro, potentially accelerating the de novo design of Covid-19 drugs.
Further, we apply MEDICO to the structural optimization of three well-known
Mpro inhibitors (N3, 11a, and GC376) and achieve ~88% improvement in their
binding affinity to Mpro, demonstrating the application value of our model for
the development of therapeutics for SARS-CoV-2 infection.


  In this work, we outline a mathematical description of biodiversity evolution
throughout the Phanerozoic based on a simple coupled system of two differential
equations and on the division of genera in two classes - Short-lived and
Long-lived types, as used by Rohde and Muller. We show that while the division
in only two classes cannot capture the complexity of biodiversity evolution in
short scales, broad trends can be described well. We also compute the systems'
Green functions as a way to qualitatively describe the possible effects of
intense perturbations of astrophyisical origin on the subsequent biodiversity
evolution.


  We propose a general framework for the study of the genealogy of neutral
discrete-time populations. We remove the standard assumption of exchangeability
of offspring distributions appearing in Cannings' models, and replace it by a
less restrictive condition of non-heritability of reproductive success. We
provide a general criterion for the weak convergence of their genealogies to
$\Xi$-coalescents, and apply it to a simple parametrization of our scenario
(which, under mild conditions, we also prove to essentially include the general
case). We provide examples for such populations, including models with
highly-asymmetric offspring distributions and populations undergoing random but
recurrent bottlenecks. Finally we study the limit genealogy of a new
exponential model which, as previously shown for related models and in spite of
its built in (fitness) inheritance mechanism, can be brought into our setting.


  Evolutionary scenarios describing the evolution of a family of genes within a
collection of species comprise the mapping of the vertices of a gene tree $T$
to vertices and edges of a species tree $S$. The relative timing of the last
common ancestors of two extant genes (leaves of $T$) and the last common
ancestors of the two species (leaves of $S$) in which they reside is indicative
of horizontal gene transfers (HGT) and ancient duplications. Orthologous gene
pairs, on the other hand, require that their last common ancestors coincides
with a corresponding speciation event. The relative timing information of gene
and species divergences is captured by three colored graphs that have the
extant genes as vertices and the species in which the genes are found as vertex
colors: the equal-divergence-time (EDT) graph, the later-divergence-time (LDT)
graph and the prior-divergence-time (PDT) graph, which together form an edge
partition of the complete graph.
  Here we give a complete characterization in terms of informative and
forbidden triples that can be read off the three graphs and provide a
polynomial time algorithm for constructing an evolutionary scenario that
explains the graphs, provided such a scenario exists. We show that every EDT
graph is perfect. While the information about LDT and PDT graphs is necessary
to recognize EDT graphs in polynomial-time for general scenarios, this extra
information can be dropped in the HGT-free case. However, recognition of EDT
graphs without knowledge of putative LDT and PDT graphs is NP-complete for
general scenarios. In contrast, PDT graphs can be recognized in
polynomial-time. We finally connect the EDT graph to the alternative
definitions of orthology that have been proposed for scenarios with horizontal
gene transfer. With one exception, the corresponding graphs are shown to be
colored cographs.


  Single-cell RNA sequencing (scRNA-seq) is powerful technology that allows
researchers to understand gene expression patterns at the single-cell level.
However, analysing scRNA-seq data is challenging due to issues and biases in
data collection. In this work, we construct an integrated Bayesian model that
simultaneously addresses normalization, imputation and batch effects and also
nonparametrically clusters cells into groups across multiple datasets. A Gibbs
sampler based on a finite-dimensional approximation of the HDP is developed for
posterior inference.


  We present a computational algebra solution to reverse engineering the
network structure of discrete dynamical systems from data. We use monomial
ideals to determine dependencies between variables that encode constraints on
the possible wiring diagrams underlying the process generating the
discrete-time, continuous-space data. Our work assumes that each variable is
either monotone increasing or decreasing. We prove that with enough data, even
in the presence of small noise, our method can reconstruct the correct unique
wiring diagram.


  Molecular simulations enable the study of biomolecules and their dynamics on
an atomistic scale. A common task is to compare several simulation conditions -
like mutations or different ligands - to find significant differences and
interrelations between them. However, the large amount of data produced for
ever larger and more complex systems often renders it difficult to identify the
structural features that are relevant for a particular phenomenon. We present a
flexible software package named PENSA that enables a comprehensive and thorough
investigation into biomolecular conformational ensembles. It provides a wide
variety of featurizations and feature transformations that allow for a complete
representation of biomolecules like proteins and nucleic acids, including water
and ion cavities within the biomolecular structure, thus avoiding bias that
would come with manual selection of features. PENSA implements various methods
to systematically compare the distributions of these features across ensembles
to find the significant differences between them and identify regions of
interest. It also includes a novel approach to quantify the state-specific
information between two regions of a biomolecule which allows, e.g., the
tracing of information flow to identify signaling pathways. PENSA also comes
with convenient tools for loading data and visualizing results in ways that
make them quick to process and easy to interpret. PENSA is an open-source
Python library maintained at https://github.com/drorlab/pensa along with an
example workflow and a tutorial. Here we demonstrate its usefulness in
real-world examples by showing how it helps to determine molecular mechanisms
efficiently.


  Mechanisms of nucleic acid accumulation were likely critical to the emergence
of life in the ferruginous oceans of the early Earth. How exactly prebiotic
geological settings accumulated nucleic acids from dilute aqueous solutions, is
poorly understood. As a possible solution to this concentration problem, we
simulated the conditions of prebiotic low-temperature alkaline hydrothermal
vents in co-precipitation experiments to investigate the potential of
ferruginous chemical gardens to accumulate nucleic acids via sorption. The
injection of an alkaline solution into an artificial ferruginous solution under
anoxic conditions (O2 <0.01% of present atmospheric levels) and at ambient
temperatures, caused the precipitation of amakinite (white rust), which quickly
converted to chloride-containing fougerite (green rust). RNA was only
extractable from the ferruginous solution in the presence of a phosphate
buffer, suggesting RNA in solution was bound to Fe2+ ions. During chimney
formation, this iron-bound RNA rapidly accumulated in the white and green rust
chimney structure, as it was depleted from the surrounding solution. Our
findings reveal that in the oceans of the early Earth, white and green rust
chimneys were likely key geochemical features that can rapidly sequester and
accumulate RNA. This represents a new mechanism for nucleic acid accumulation,
in addition to wet dry cycles, and may have promoted RNA survival in a dilute
prebiotic ocean.


  The self-protection of alliances against external invaders is a key concept
behind the maintenance of biodiversity in the face of natural selection. But
since these alliances, which can be formed by different numbers of competitors,
can also compete against each other, it is important to identify their
strengths and weaknesses. Here, we therefore compare the vitalities of two
two-species alliances whose members either beat each other mutually via a
bidirectional invasion or they exchange their positions during an inner
dynamics. The resulting four-species model shows rich behavior in dependence on
the model parameter $p$, which characterizes the inner invasions, and $\beta$,
which determines the intensity of site exchanges. In the low $p$ and the large
$p$ limit, when the inner invasion becomes biased, three-member
rock-scissors-paper-type solutions emerge, where one of the members is
oppressed by having the smallest average concentration due to heterogeneous
inner invasion rates. Interestingly, however, if we allow a more intensive site
exchange between the oppressed species, they can morph into a winning pair and
dominate the full parameter plane. We show that their victory utilizes the
vulnerability of the rival alliance based on cyclic dominance, where a species
can easily fixate a limited-size domain.


  Motivation: A Genomic Dictionary, i.e., the set of the k-mers appearing in a
genome, is a fundamental source of genomic information: its collection is the
first step in strategic computational methods ranging from assembly to sequence
comparison and phylogeny. Unfortunately, it is costly to store. This motivates
some recent studies regarding the compression of those k-mer sets. However,
such an area does not have the maturity of genomic compression, lacking an
homogeneous and methodologically sound experimental foundation that allows to
fairly compare the relative merits of the available solutions, and that takes
into account also the rich choices of compression methods that can be used.
  Results: We provide such a foundation here, supporting it with an extensive
set of experiments that use reference datasets and a carefully selected set of
representative data compressors. Our results highlight the spectrum of
compressor choices one has in terms of Pareto Optimality of compression vs.
post-processing, this latter being important when the Dictionary needs to be
decompressed many times. In addition to the useful indications, not available
elsewhere, that this study offers to the researchers interested in storing
k-mer dictionaries in compressed form, a software system that can be readily
used to explore the Pareto Optimal solutions available r a given Dictionary is
also provided.
  Availability: The software system is available at
https://github.com/GenGrim76/Pareto-Optimal-GDC, together with user manuals and
installation instructions.
  Contact: raffaele.giancarlo@unipa.it
  Supplementary information: Additional data are available in the Supplementary
Material.


  Homophily and heterophobia, the tendency for people with similar
characteristics to preferentially interact with (or avoid) each other are
pervasive in human social networks. Here, we develop an extension of the
mathematical theory of urban scaling which describes the effects of homophily
and heterophobia on social interactions and resulting economic outputs of
cities. Empirical tests of our model show that increased residential racial
heterophobia and segregation in U.S. cities are associated with reduced
economic outputs and that the strength of this relationship increased
throughout the 2010s. Our findings provide the means for the formal
incorporation of general homophilic and heterophobic effects into theories of
modern urban science and suggest that racial segregation is increasingly and
adversely impacting the economic performance and connectivity of urban
societies in the U.S.


  By demonstrating the role that historical population replacements and waves
of admixture have played around the world, the genetics work of Reich and
colleagues has provided a paradigm for understanding human history [Reich et
al. 2009; Reich et al. 2012; Patterson et al. 2012]. Although we show in
Ioannidis et al. [2021] that the peopling of Polynesia was a range expansion,
and not, as suggested by Huang et al. [2022], yet another example of waves of
admixture and large-scale gene flow between populations, we believe that our
result in this recently settled oceanic expanse is the exception that proves
the rule.


  Predicting the physical interaction of proteins is a cornerstone problem in
computational biology. New classes of learning-based algorithms are actively
being developed, and are typically trained end-to-end on protein complex
structures extracted from the Protein Data Bank. These training datasets tend
to be large and difficult to use for prototyping and, unlike image or natural
language datasets, they are not easily interpretable by non-experts. We present
Dock2D-IP and Dock2D-IF, two "toy" datasets that can be used to select
algorithms predicting protein-protein interactions$\unicode{x2014}$or any other
type of molecular interactions. Using two-dimensional shapes as input, each
example from Dock2D-IP ("interaction pose") describes the interaction pose of
two shapes known to interact and each example from Dock2D-IF ("interaction
fact") describes whether two shapes form a stable complex or not. We propose a
number of baseline solutions to the problem and show that the same underlying
energy function can be learned either by solving the interaction pose task
(formulated as an energy-minimization "docking" problem) or the
fact-of-interaction task (formulated as a binding free energy estimation
problem).


  There is considerable interest in predicting the pathogenicity of protein
variants in human genes. Due to the sparsity of high quality labels, recent
approaches turn to \textit{unsupervised} learning, using Multiple Sequence
Alignments (MSAs) to train generative models of natural sequence variation
within each gene. These generative models then predict variant likelihood as a
proxy to evolutionary fitness. In this work we instead combine this
evolutionary principle with pretrained protein language models (LMs), which
have already shown promising results in predicting protein structure and
function. Instead of training separate models per-gene, we find that a single
protein LM trained on broad sequence datasets can score pathogenicity for any
gene variant zero-shot, without MSAs or finetuning. We call this unsupervised
approach \textbf{VELM} (Variant Effect via Language Models), and show that it
achieves scoring performance comparable to the state of the art when evaluated
on clinically labeled variants of disease-related genes.


  In stochastic modeling of infectious diseases, it has been established that
variations in infectivity affect the probability of a major outbreak, but not
the shape of the curves during a major outbreak, which is predicted by
deterministic models [Diekmann et al., 2012]. However, such conclusions are
derived under idealized assumptions such as the population size tending to
infinity, and the individual degree of infectivity only depending on variations
in the infectiousness period. In this paper we show that the same conclusions
hold true in a finite population representing a medium size city, where the
degree of infectivity is determined by the offspring distribution, which we try
to make as realistic as possible for SARS-CoV-2. In particular, we consider
distributions with fat tails, to incorporate the existence of super-spreaders.
We also provide new theoretical results on convergence of stochastic models
which allows to incorporate any offspring distribution with a finite variance.


  Epiabm is a fully tested, open-source software package for epidemiological
agent-based modelling, re-implementing the well-known CovidSim model from the
MRC Centre for Global Infectious Disease Analysis at Imperial College London.
It has been developed as part of the first-year training programme in the EPSRC
SABS:R3 Centre for Doctoral Training at the University of Oxford. The model
builds an age-stratified, spatially heterogeneous population and offers a
modular approach to configure and run epidemic scenarios, allowing for a broad
scope of investigative and comparative studies. Two simulation backends are
provided: a pedagogical Python backend (with full functionality) and a high
performance C++ backend for use with larger population simulations. Both are
highly modular, with comprehensive testing and documentation for ease of
understanding and extensibility. Epiabm is publicly available through GitHub at
https://github.com/SABS-R3-Epidemiology/epiabm.


  Basecalling is an essential step in nanopore sequencing analysis where the
raw signals of nanopore sequencers are converted into nucleotide sequences,
i.e., reads. State-of-the-art basecallers employ complex deep learning models
to achieve high basecalling accuracy. This makes basecalling
computationally-inefficient and memory-hungry; bottlenecking the entire genome
analysis pipeline. However, for many applications, the majority of reads do no
match the reference genome of interest (i.e., target reference) and thus are
discarded in later steps in the genomics pipeline, wasting the basecalling
computation. To overcome this issue, we propose TargetCall, the first
pre-basecalling filter to eliminate the wasted computation in basecalling.
TargetCall's key idea is to discard reads that will not match the target
reference (i.e., off-target reads) prior to basecalling. TargetCall consists of
two main components: (1) LightCall, a lightweight neural network basecaller
that produces noisy reads; and (2) Similarity Check, which labels each of these
noisy reads as on-target or off-target by matching them to the target
reference. TargetCall aims to filter out all off-target reads before
basecalling. The highly-accurate but slow basecalling is performed only on the
raw signals whose noisy reads are labeled as on-target. Our thorough
experimental evaluations using both real and simulated data show that
TargetCall 1) improves the end-to-end basecalling performance while maintaining
high sensitivity in keeping on-target reads, 2) maintains high accuracy in
downstream analysis, 3) precisely filters out up to 94.71% of off-target reads,
and 4) achieves better performance, throughput, sensitivity, precision, and
generality compared to prior works. We open-source TargetCall at
https://github.com/CMU-SAFARI/TargetCall


  Fragment molecular orbital (FMO) calculations were performed in aqueous media
which allowed us to obtain the interaction energy between the human estrogen
receptor $\alpha$ ligand-binding domain (ER) and the selected ligands (L):
17$\beta$-estradiol (E2), 17$\alpha$-estradiol (17$\alpha$-E2), estriol (E3),
genistein (GNT), diethylstilbestrol (DES), bisphenol A (BPA), bisphenol AF
(BPAF), hydroxychlor (HPTE) and methoxychlor (DMDT). These calculations were
carried out on representative structures of L-ER complexes obtained from
molecular dynamics simulations. The MP2/6-31G(d) L-ER FMO interaction energy in
kcal/mol is as follows:
  E3 (-100.1) < GNT (-95.8) < E2 (-88.5) < BPA (-84.7) < DES (-82.6) < BPAF
(-80.6) < 17$\alpha$-E2 (-78.7) < HPTE (-75.9) < DMDT (-46.3)
  The central hydrophobic core of the ligands interacts attractively with
several apolar amino acid residues of ER. Glu 353 and His 524 interacts
strongly with most ligands through a hydrogen bond with the hydroxyl group of
the phenol A-ring and the terminal hydroxylated ring, respectively. Water
molecules were found at the binding site of receptor. In our model systems we
have demonstrated what is generally observed in ligand-receptor complexes: the
steric and chemical complementarity of the groups on the ligand and binding
site surfaces.


  Building prediction models for outcomes of clinical relevance when only a
limited number of mutational features are available causes considerable
challenges due to the sparseness and low-dimensionality of the data. In this
article, we present a method to augment the predictive power of these features
by leveraging multi-modal associative relationships between an individual's
mutational profile and their corresponding gene expression or knockout effect
profiles. We can thus reconstruct expression or effect scores for genes of
interest from the available mutation features and then use this reconstructed
representation directly to model and predict clinical outcomes. We show that
our method produces significant improvements in predictive accuracy compared to
models utilizing only the raw mutational data, and results in conclusions
comparable to those obtained using real expression or effect profiles.


  Single-cell sequencing assay for transposase-accessible chromatin
(scATAC-seq) assesses genome-wide chromatin accessibility in thousands of cells
to reveal regulatory landscapes in high resolutions. However, the analysis
presents challenges due to the high dimensionality and sparsity of the data.
Several methods have been developed, including transformation techniques of
term-frequency inverse-document frequency (TF-IDF), dimension reduction methods
such as singular value decomposition (SVD), factor analysis, and autoencoders.
Yet, a comprehensive study on the mentioned methods has not been fully
performed. It is not clear what is the best practice when analyzing scATAC-seq
data. We compared several scenarios for transformation and dimension reduction
as well as the SVD-based feature analysis to investigate potential enhancements
in scATAC-seq information retrieval. Additionally, we investigate if
autoencoders benefit from the TF-IDF transformation. Our results reveal that
the TF-IDF transformation generally leads to improved clustering and
biologically relevant feature extraction.


  It is well-known that interactions between species determine the population
composition in an ecosystem. Conventional studies have focused on fixed
population structures to reveal how interactions shape population compositions.
However, interaction structures are not fixed, but change over time due to
invasions. Thus, invasion and interaction play an important role in shaping
communities. Despite its importance, however, the interplay between invasion
and interaction has not been well explored. Here, we investigate how invasion
affects the population composition with interactions in open evolving systems
considering generalized Lotka-Volterra-type dynamics. Our results show that the
system has two distinct regimes. One is characterized by low diversity with
abrupt changes of dominant species in time, appearing when the interaction
between species is strong and invasion slowly occurs. On the other hand,
frequent invasions can induce higher diversity with slow changes in abundances
despite strong interactions. It is because invasion happens before the system
reaches its equilibrium, which drags the system from its equilibrium all the
time. All species have similar abundances in this regime, which implies that
fast invasion induces regime shift. Therefore, whether invasion or interaction
dominates determines the population composition.


  The bacterium Xylella fastidiosa (Xf) is mainly transmitted by the
spittlebug, Philaenus spumarius, in Europe, where it has caused significant
economic damage to olive and almond trees. Understanding the factors that
determine disease dynamics in pathosystems that share similarities can help
design control strategies focused on minimizing transmission chains. Here we
introduce a compartmental model for Xf-caused diseases in Europe that accounts
for the main relevant epidemiological processes, including the seasonal
dynamics of P. spumarius. The model was confronted with epidemiological data
from the two major outbreaks of Xf in Europe, the olive quick disease syndrome
(OQDS) in Apulia, Italy, caused by the subspecies pauca, and the almond leaf
scorch disease (ALSD) in Majorca, Spain, caused by subspecies multiplex and
fastidiosa. Using a Bayesian inference framework, we show how the model
successfully reproduces the general field data in both diseases. In a global
sensitivity analysis, the vector-plant and plant-vector transmission rates,
together with the vector removal rate, were the most influential parameters in
determining the time of the infected host population peak, the incidence peak
and the final number of dead hosts. We also used our model to check different
vector-based control strategies, showing that a joint strategy focused on
increasing the rate of vector removal while lowering the number of annual
newborn vectors is optimal for disease control.


  Cryo-electron microscopy (cryo-EM) has recently become a premier method for
obtaining high-resolution structures of biological macromolecules. However, it
is limited to biomolecular samples with low conformational heterogeneity, where
all the conformations can be well-sampled at many projection angles. While
cryo-EM technically provides single-molecule data for heterogeneous molecules,
most existing reconstruction tools cannot extract the full distribution of
possible molecular configurations. To overcome these limitations, we build on a
prior Bayesian approach and develop an ensemble refinement framework that
estimates the ensemble density from a set of cryo-EM particles by reweighting a
prior ensemble of conformations, e.g., from molecular dynamics simulations or
structure prediction tools. Our work is a general approach to recovering the
equilibrium probability density of the biomolecule directly in conformational
space from single-molecule data. To validate the framework, we study the
extraction of state populations and free energies for a simple toy model and
from synthetic cryo-EM images of a simulated protein that explores multiple
folded and unfolded conformations.


  We have presented some arguments to substantiate the usefulness of neural
replicator analysis (NRA) for constructing variants of the natural binomial
classification of virus genomes based only on knowledge of their complete
genomic sequences, without involving other data on the phenotype, functions,
encoded proteins, etc., and also without the need of genomic sequences
alignment. Perhaps this will make sense when processing metagenomic data. This
makes it possible to construct the binomial classification accepted for the
viruses themselves. We restrict ourselves to three families of viruses having
dsDNA circular genomes (Papillomaviridae, Polyomaviridae and Caulimoviridae)
and partly to the family Geminiviridae having ssDNA genomes though the approach
presented can be also applied to genomes of other dsDNA, ssDNA and ssRNA
viruses, including linear ones (some results for Mitoviridae are also
presented). It is argued that binomial classification of virus genomes which is
difficult to apply in all cases can nevertheless be informative tool of
revealing virus properties, areal of hosts, forms of diseases and can also show
the connections of the viruses belonging to different families and even to
different kingdoms.


  The Leinster and Cobbold diversity index possesses a number of merits; in
particular, it generalises many existing indices and defines an effective
number. We present a scheme to quantify the contribution of richness, evenness,
and taxonomic similarity to this index. Compared to the work of van Dam (2019),
our approach gives unbiased estimates of both evenness and similarity in a
non-homogeneous community. We also introduce a notion of taxonomic tree
equilibration which should be of use in the description of community structure.


  Our Microbiome Network Alignment Algorithm (MiNAA) aligns two microbial
networks using a combination of the GRAph ALigner (GRAAL) algorithm and the
Hungarian algorithm. Network alignment algorithms find pairs of nodes (one node
from the first network and the other node from the second network) that have
the highest similarity. Traditionally, similarity has been defined as
topological similarity such that the neighborhoods around the two nodes are
similar. Recent implementations of network alignment methods such as NETAL and
L-GRAAL also include measures of biological similarity, yet these methods are
restricted to one specific type of biological similarity (e.g. sequence
similarity in L-GRAAL). Our work extends existing network alignment
implementations by allowing any type of biological similarity to be input by
the user. This flexibility allows the user to choose whatever measure of
biological similarity is suitable for the study at hand. In addition, unlike
most existing network alignment methods that are tailored for protein or gene
interaction networks, our work is the first one suited for microbiome networks.


  Using theory and simulations, we carried out a first systematic
characterization of DNA unzipping via nanopore translocation. Starting from
partially unzipped states, we found three dynamical regimes depending on the
applied force, f: (i) heterogeneous DNA retraction and rezipping (f < 17pN),
(ii) normal (17pN < f < 60pN) and (iii) anomalous (f > 60pN) drift-diffusive
behavior. We show that the normal drift-diffusion regime can be effectively
modelled as a one-dimensional stochastic process in a tilted periodic
potential. We use the theory of stochastic processes to recover the potential
from nonequilibrium unzipping trajectories and show that it corresponds to the
free-energy landscape for single base-pairs unzipping. Applying this general
approach to other single-molecule systems with periodic potentials ought to
yield detailed free-energy landscapes from out-of-equilibrium trajectories.


  Protein subcellular localization is an important factor in normal cellular
processes and disease. While many protein localization resources treat it as
static, protein localization is dynamic and heavily influenced by biological
context. Biological pathways are graphs that represent a specific biological
context and can be inferred from large-scale data. We develop graph algorithms
to predict the localization of all interactions in a biological pathway as an
edge-labeling task. We compare a variety of models including graph neural
networks, probabilistic graphical models, and discriminative classifiers for
predicting localization annotations from curated pathway databases. We also
perform a case study where we construct biological pathways and predict
localizations of human fibroblasts undergoing viral infection. Pathway
localization prediction is a promising approach for integrating publicly
available localization data into the analysis of large-scale biological data.


  BEF studies aim to understand how ecosystems respond to a gradient of species
diversity. Diversity-Interactions (DI) models are suitable for analysing the
BEF relationship. These models relate an ecosystem function response of a
community to the identity of the species in the community, their evenness
(proportions) and interactions. The number of species in the community
(richness) is also implicitly modelled through this approach. It is common in
BEF studies to model an ecosystem function as a function of richness; while
this can uncover trends in the BEF relationship, by definition, species
diversity is much broader than richness alone, and important patterns in the
BEF relationship may remain hidden. In this paper, we introduce the DImodels R
package for implementing DI models. We also compare DI models to traditional
modelling approaches to highlight the advantages of using a multi-dimensional
definition of species diversity. We show that using DI models can lead to
considerably improved model fit over other methods; it does this by
incorporating variation due to the multiple facets of species diversity.
Predicting from a DI model is not limited to the study design points, the model
can extrapolate to predict for any species composition and proportions
(assuming there is sufficient coverage of this space in the study design).
Expressing the BEF relationship as a function of richness alone can be useful
to capture overall trends. However, collapsing the multiple dimensions of
species diversity to a single dimension (such as richness) can result in
valuable ecological information being lost. DI modelling provides a framework
to test the multiple components of species diversity in the BEF relationship.
It facilitates uncovering a deeper ecological understanding of the BEF
relationship and can lead to enhanced inference.


  Ecosystems represent archetypal complex dynamical systems, often modelled by
coupled differential equations of the form $$ \frac{d x_i}{d t} = x_i
\varphi_i(x_1,\cdots, x_N)\ , $$ where $N$ represents the number of species and
$x_i$, the abundance of species $i$. Among these families of coupled diffential
equations, Lotka-Volterra (LV) equations $$ \frac{d x_i}{d t} = x_i ( r_i - x_i
+(\Gamma \mathbf{x})_i)\ , $$ play a privileged role, as the LV model
represents an acceptable trade-off between complexity and tractability. Here,
$r_i$ represents the intrinsic growth of species $i$ and $\Gamma$ stands for
the interaction matrix: $\Gamma_{ij}$ represents the effect of species $j$ over
species $i$. For large $N$, estimating matrix $\Gamma$ is often an overwhelming
task and an alternative is to draw $\Gamma$ at random, parametrizing its
statistical distribution by a limited number of model features. Dealing with
large random matrices, we naturally rely on Random Matrix Theory (RMT).
  The aim of this review article is to present an overview of the work at the
junction of theoretical ecology and large random matrix theory. It is intended
to an interdisciplinary audience spanning theoretical ecology, complex systems,
statistical physics and mathematical biology.


  Even though deep neural networks (DNNs) achieve state-of-the-art results for
a number of problems involving genomic data, getting DNNs to explain their
decision-making process has been a major challenge due to their black-box
nature. One way to get DNNs to explain their reasoning for prediction is via
attribution methods which are assumed to highlight the parts of the input that
contribute to the prediction the most. Given the existence of numerous
attribution methods and a lack of quantitative results on the fidelity of those
methods, selection of an attribution method for sequence-based tasks has been
mostly done qualitatively. In this work, we take a step towards identifying the
most faithful attribution method by proposing a computational approach that
utilizes point mutations. Providing quantitative results on seven popular
attribution methods, we find Layerwise Relevance Propagation (LRP) to be the
most appropriate one for translation initiation, with LRP identifying two
important biological features for translation: the integrity of Kozak sequence
as well as the detrimental effects of premature stop codons.


  Mixtures of several macromolecular species can lead to the formation of
higher-order structures that often display non-ideal mixing behavior. In this
work, we propose a minimal model of a quaternary system which considers the
formation of a complex via a chemical reaction involving two macromolecular
species; the complex may then phase separate from the buffer and undergo a
further transition into a gel-like state over time. First, a ternary phase
diagram that captures the volume fraction of each species and phases at
equilibrium is constructed. Specifically, we investigate how physical
parameters such as stoichiometric coefficients, molecular sizes and interaction
parameters affect LLPS and aging. Finally, we analyze the thermodynamic
stability of the two-phase system and identify the spinodal regions, and
outline the generalization of our approach to reactive biomolecular systems
with an arbitrary number of components.


  Aims To propose and analyze a general, dynamic, process-oriented theory of
the area of distribution. Methods The area of distribution is modelled by
combining (by multiplication) three matrices: one matrix represents movements,
another niche tolerances, and a third, biotic interactions. Results are derived
from general properties of this product and from simulation of a cellular
automaton defined in terms of the matrix operations. Everything is implemented
practically in an R package. Results Results are obtained by simulation and by
mathematical analysis. We show that the mid-domain effect is a direct
consequence of dispersal; that to include movements to Ecological Niche
Modeling significantly affects results, but cannot be done without choosing an
ancestral area of distribution. We discuss ways of estimating such ancestral
areas. We show that, in our approach, movements and niche effects are mixed in
ways almost impossible to disentangle, and show this is a consequence of the
singularity of a matrix. We introduce a tool (the
Connectivity-Suitability-Dispersal plot) to extend the results of simple niche
modeling to understand the effects of dispersal. Main conclusions The
conceptually straightforward scheme we present for the area of distribution
integrates, in a mathematically sound and computationally feasible way, several
key ideas in biogeography: the geographic and environmental matrix, the
Grinnellian niche, dispersal capacity and the ancestral area of origin of
groups of species. We show that although full simulations are indispensable to
obtain the dynamics of an area of distribution, interesting results can be
derived simply by analyzing the matrices representing the dynamics.


  Spatial transcriptomics has the potential to transform our understanding of
RNA expression in tissues. Classical array-based technologies produce
multiple-cell-scale measurements requiring deconvolution to recover single cell
information. However, rapid advances in subcellular measurement of RNA
expression at whole-transcriptome depth necessitate a fundamentally different
approach. To integrate single-cell RNA-seq data with nanoscale spatial
transcriptomics, we present a topological method for automatic cell type
identification (TopACT). Unlike popular decomposition approaches to
multicellular resolution data, TopACT is able to pinpoint the spatial locations
of individual sparsely dispersed cells without prior knowledge of cell
boundaries. Pairing TopACT with multiparameter persistent homology landscapes
predicts immune cells forming a peripheral ring structure within kidney
glomeruli in a murine model of lupus nephritis, which we experimentally
validate with immunofluorescent imaging. The proposed topological data analysis
unifies multiple biological scales, from subcellular gene expression to
multicellular tissue organization.


  We present an exploratory population dynamics approach, described by
Lotka-Volterra (LV) generalized equations, to explain/predict the dynamics and
competition between land use/land cover (LULC) classes over vulnerable regions
in Uruguay. We use the Mapbiomas-Pampa dataset composed by 20 annual LULC maps
from 2000-2019. From these LULC maps we extract the main LULC classes, their
spatial distribution and the time series of areas covered for each class. The
interaction coefficients between species are inferred through the pairwise
maximum entropy (PME) method from the spatial covariance matrices for different
training periods. The main finding is that this LVPME method globally
outperforms the more traditional Markov chains approach at predicting the
trajectories of areas of LULC classes.


  In a wide variety of natural systems, closely-related microbial strains
coexist stably, resulting in high levels of fine-scale biodiversity. However,
the mechanisms that stabilize this coexistence are not fully understood.
Spatial heterogeneity is one common stabilizing mechanism, but the rate at
which organisms disperse throughout the heterogeneous environment may strongly
impact the stabilizing effect that heterogeneity can provide. An intriguing
example is the gut microbiome, where active mechanisms exist to control the
movement of microbes and potentially maintain diversity. We investigate how
biodiversity is affected by migration rate using a simple evolutionary model
with heterogeneous selection pressure. We find that the biodiversity-migration
rate relationship is shaped by multiple phase transitions, including a
reentrant phase transition to coexistence. At each transition, an ecotype goes
extinct and dynamics exhibit critical slowing down (CSD). CSD is encoded in the
statistics of fluctuations due to demographic noise -- this may provide an
experimental means for detecting and altering impending extinction.


  We performed a thorough sensitivity analysis of the herd immunity threshold
for discrete-time SIR compartmental models with a static network structure. We
find unexpectedly that these models violate classical intuition which holds
that the herd immunity threshold should monotonically increase with the
transmission parameter. We find the existence of bifurcations in the herd
immunity threshold in the high transmission probability regime. The extent of
these bifurcations is modulated by the graph heterogeneity, the recovery
parameter, and the network size. In the limit of large, well-mixed networks,
the behavior approaches that of difference equation models, suggesting this
behavior is a universal feature of all discrete-time SIR models. These results
suggest careful attention is needed in both selecting the assumptions on how to
model time and heterogeneity in epidemiological models and the subsequent
conclusions that can be drawn.


  Protein-ligand (un)binding simulations are a recent focus of biased molecular
dynamics simulations. Such binding and unbinding can occur via different
pathways in and out of a binding site. We here present a theoretical framework
how to compute kinetics along separate paths and to combine the path-specific
rates into global binding and unbinding rates for comparison with experiment.
Using dissipation-corrected targeted molecular dynamics in combination with
temperature-boosted Langevin equation simulations [Nat. Commun. \textbf{11},
2918 (2020)] applied to a two-dimensional model and the trypsin-benzamidine
complex as test systems, we assess the robustness of the procedure and discuss
aspects of its practical applicability to predict multisecond kinetics of
complex biomolecular systems.


  The transmission dynamics of infectious diseases in animal production are
driven by several propagation routes. Contaminated vehicles traveling between
farms have been associated with indirect disease transmission. In this study,
we used transportation vehicle data to analyze the magnitude of farm visits by
different vehicles and to propose a methodology to reconstruct vehicle contact
networks considering pathogen stability and cleaning and disinfection
effectiveness. Here, we collected information from 6,363 farms and Global
Positioning System (GPS) records from 567 vehicles used to transport feed,
animals, and people. We reconstructed vehicle contacts among the farms,
conserving pathogen stability decay and different probabilities of cleaning and
disinfection. Results showed that vehicle movement networks were densely
connected, with up to 86% of farms connected by these movements. Movements of
vehicle transporting feed and pig among farms showed the highest network
connectivity. The cleaning effectiveness of was variable among the different
vehicle types and highly influenced by the frequency of vehicles stopping at
clean stations. A large number of between-farm contacts with a pathogen
stability >0.8 were present in the vehicle network even with 100% cleaning
effectiveness. Finally, we identified that vehicles contacted farms from
different companies. Thus, our results suggest the vehicle network is a
potential mechanism for spreading pathogens among farms. Moreover, even with
scenarios with high effectiveness of cleaning and disinfection, the risk of
vehicles spreading diseases was not completely eliminated.


  A generalized understanding of protein dynamics is an unsolved scientific
problem, the solution of which is critical to the interpretation of the
structure-function relationships that govern essential biological processes.
Here, we approach this problem by constructing coarse-grained molecular
potentials based on artificial neural networks and grounded in statistical
mechanics. For training, we build a unique dataset of unbiased all-atom
molecular dynamics simulations of approximately 9 ms for twelve different
proteins with multiple secondary structure arrangements. The coarse-grained
models are capable of accelerating the dynamics by more than three orders of
magnitude while preserving the thermodynamics of the systems. Coarse-grained
simulations identify relevant structural states in the ensemble with comparable
energetics to the all-atom systems. Furthermore, we show that a single
coarse-grained potential can integrate all twelve proteins and can capture
experimental structural features of mutated proteins. These results indicate
that machine learning coarse-grained potentials could provide a feasible
approach to simulate and understand protein dynamics.


  In recent years it has become evident the need of understanding how failure
of coordination imposes constraints on the size of stable groups that highly
social mammals can live in. We examine here the forces that keep animals
together as a herd and others that drive them apart. Different phenotypes (e.g.
genders) have different rates of gut fill, causing them to spend different
amounts of time performing activities. By modeling a group as a set of
semi-coupled oscillators on a disc, we show that the members of the group may
become less and less coupled until the group dissolves and breaks apart. We
show that when social bonding creates a stickiness, or gravitational pull,
between pairs of individuals, fragmentation is reduced.


  Genome sequencing technology has improved significantly in few last years and
resulted in abundance genetic data. Artificial intelligence has been employed
to analyze genetic data in response to its sheer size and variability. Gene
prediction on single DNA has been conducted using various deep learning
architectures to discover splice sites and therefore discover intron and exon
region. Recent predictions are carried out with models trained on sequence with
fixed splice site location which eliminates possibility of multiple splice
sites existence in single sequence. This paper proposes sequential labelling to
predict splice sites regardless their position in sequence. Sequential
labelling is carried out on DNA to determine intron and exon region and thus
discover splice sites. Sequential labelling models used are based on pretrained
DNABERT-3 which has been trained on human genome. Both fine-tuning and
feature-based approach are tested. Proposed model is benchmarked against latest
sequential labelling model designed for mutation type and location prediction.
While achieving high F1 scores on validation data, both baseline and proposed
model perform poorly on test data. Error and test results analysis reveal that
model experience overfitting and therefore, model is deemed not suitable for
splice site prediction.


  In Dec 2020, the results of AlphaFold2 were presented at CASP14, sparking a
revolution in the field of protein structure predictions. For the first time, a
purely computational method could challenge experimental accuracy for structure
prediction of single protein domains. The code of AlphaFold2 was released in
the summer of 2021, and since then, it has been shown that it can be used to
accurately predict the structure of most (ordered) proteins and many
protein-protein interactions. It has also sparked an explosion of development
in the field, improving AI-based methods to predict protein complexes,
disordered regions, and protein design. Here I will review some of the
inventions sparked by the release of AlphaFold.


  We propose an agent-based opinion formation model characterised by a two-fold
novelty. First, we realistically assume that each agent cannot measure the
opinion of its neighbours with infinite resolution and accuracy, and hence it
can only classify the opinion of others as agreeing much more, or more, or
comparably, or less, or much less (than itself) with a given statement. This
leads to a classification-based rule for opinion update. Second, we consider
three complementary agent traits suggested by significant sociological and
psychological research: conformism, radicalism and stubbornness. We rely on
World Values Survey data to show that the proposed model has the potential to
predict the evolution of opinions in real life: the classification-based
approach and complementary agent traits produce rich collective behaviours,
such as polarisation, consensus, and clustering, which can yield predicted
opinions similar to survey results.


  Recombination is a fundamental evolutionary force, but it is difficult to
quantify because the effect of a recombination event on patterns of variation
in a sample of genetic data can be hard to discern. Estimators for the
recombination rate, which are usually based on the idea of integrating over the
unobserved possible evolutionary histories of a sample, can therefore be noisy.
Here we consider a related question: how would an estimator behave if the
evolutionary history actually was observed? This would offer an upper bound on
the performance of estimators used in practice. In this paper we derive an
expression for the maximum likelihood estimator for the recombination rate
based on a continuously observed, multi-locus, Wright--Fisher diffusion of
haplotype frequencies, complementing existing work for an estimator of
selection. We show that, contrary to selection, the estimator has unusual
properties because the observed information matrix can explode in finite time
whereupon the recombination parameter is learned without error. We also show
that the recombination estimator is robust to the presence of selection in the
sense that incorporating selection into the model leaves the estimator
unchanged. We study the properties of the estimator by simulation and show that
its distribution can be quite sensitive to the underlying mutation rates.


  In molecular research, simulation \& design of molecules are key areas with
significant implications for drug development, material science, and other
fields. Current classical computational power falls inadequate to simulate any
more than small molecules, let alone protein chains on hundreds of peptide.
Therefore these experiment are done physically in wet-lab, but it takes a lot
of time \& not possible to examine every molecule due to the size of the search
area, tens of billions of dollars are spent every year in these research
experiments. Molecule simulation \& design has lately advanced significantly by
machine learning models, A fresh perspective on the issue of chemical synthesis
is provided by deep generative models for graph-structured data. By optimising
differentiable models that produce molecular graphs directly, it is feasible to
avoid costly search techniques in the discrete and huge space of chemical
structures. But these models also suffer from computational limitations when
dimensions become huge and consume huge amount of resources. Quantum Generative
machine learning in recent years have shown some empirical results promising
significant advantages over classical counterparts.


  Phylogenetic networks are a special type of graph that generalize
phylogenetic trees and that are used to model non-treelike evolutionary
processes such as recombination and hybridization. In this paper, we consider
unrooted phylogenetic networks, i.e. simple, connected graphs
$\mathcal{N}=(V,E)$ with leaf set $X$, for $X$ some set of species, in which
every internal vertex in $\mathcal{N}$ has degree three. One approach used to
construct such phylogenetic networks is to take as input a collection
$\mathcal{P}$ of phylogenetic trees and to look for a network $\mathcal{N}$
that contains each tree in $\mathcal{P}$ and that minimizes the quantity
$r(\mathcal{N}) = |E|-(|V|-1)$ over all such networks. Such a network always
exists, and the quantity $r(\mathcal{N})$ for an optimal network $\mathcal{N}$
is called the hybrid number of $\mathcal{P}$. In this paper, we give a new
characterization for the hybrid number in case $\mathcal{P}$ consists of two
trees. This characterization is given in terms of a cherry picking sequence for
the two trees, although to prove that our characterization holds we need to
define the sequence more generally for two forests. Cherry picking sequences
have been intensively studied for collections of rooted phylogenetic trees, but
our new sequences are the first variant of this concept that can be applied in
the unrooted setting. Since the hybrid number of two trees is equal to the
well-known tree bisection and reconnect distance between the two trees, our new
characterization also provides an alternative way to understand this important
tree distance.


  The space of possible human cultures is vast, but some cultural
configurations are more consistent with cognitive and social constraints than
others. This leads to a ``landscape'' of possibilities that our species has
explored over millennia of cultural evolution. But what does this fitness
landscape, which constrains and guides cultural evolution, look like? The
machine-learning algorithms that can answer these questions are typically
developed for large-scale datasets. Applications to the sparse, inconsistent,
and incomplete data found in the historical record have received less
attention, and standard recommendations can lead to bias against marginalized,
under-studied, or minority cultures. We show how to adapt the Minimum
Probability Flow algorithm and the Inverse Ising model, a physics-inspired
workhorse of machine learning, to the challenge. A series of natural extensions
-- including dynamical estimation of missing data, and cross-validation with
regularization -- enables reliable reconstruction of the underlying
constraints. We demonstrate our methods on a curated subset of the Database of
Religious History: records from 407 religious groups throughout human history,
ranging from the Bronze Age to the present day. This reveals a complex, rugged,
landscape, with both sharp, well-defined peaks where state-endorsed religions
tend to concentrate, and diffuse cultural floodplains where evangelical
religions, non-state spiritual practices, and mystery religions can be found.


  Dengue fever is a virulent disease spreading over 100 tropical and
subtropical countries in Africa, the Americas, and Asia. This arboviral disease
affects around 400 million people globally, severely distressing the healthcare
systems. The unavailability of a specific drug and ready-to-use vaccine makes
the situation worse. Hence, policymakers must rely on early warning systems to
control intervention-related decisions. Forecasts routinely provide critical
information for dangerous epidemic events. However, the available forecasting
models (e.g., weather-driven mechanistic, statistical time series, and machine
learning models) lack a clear understanding of different components to improve
prediction accuracy and often provide unstable and unreliable forecasts. This
study proposes an ensemble wavelet neural network with exogenous factor(s)
(XEWNet) model that can produce reliable estimates for dengue outbreak
prediction for three geographical regions, namely San Juan, Iquitos, and
Ahmedabad. The proposed XEWNet model is flexible and can easily incorporate
exogenous climate variable(s) confirmed by statistical causality tests in its
scalable framework. The proposed model is an integrated approach that uses
wavelet transformation into an ensemble neural network framework that helps in
generating more reliable long-term forecasts. The proposed XEWNet allows
complex non-linear relationships between the dengue incidence cases and
rainfall; however, mathematically interpretable, fast in execution, and easily
comprehensible. The proposal's competitiveness is measured using computational
experiments based on various statistical metrics and several statistical
comparison tests. In comparison with statistical, machine learning, and deep
learning methods, our proposed XEWNet performs better in 75% of the cases for
short-term and long-term forecasting of dengue incidence.


  Conformational changes are observed in many enzymes, but their role in
catalysis is highly controversial. Here we present a theoretical model that
illustrates how rigid catalysts can be fundamentally limited and how a
conformational change induced by substrate binding can overcome this
limitation, ultimately enabling barrier-free catalysis. The model is
deliberately minimal, but the principle it illustrates is general and
consistent with unique features of proteins as well as with previous informal
proposals to explain the superiority of enzymes over other classes of
catalysts. Implementing the discriminative switch suggested by the model could
help overcome limitations currently encountered in the design of artificial
catalysts.


  With the development of gene sequencing technology, an explosive growth of
gene data has been witnessed. And the storage of gene data has become an
important issue. Traditional gene data compression methods rely on general
software like G-zip, which fails to utilize the interrelation of nucleotide
sequence. Recently, many researchers begin to investigate deep learning based
gene data compression method. In this paper, we propose a transformer-based
gene compression method named GeneFormer. Specifically, we first introduce a
modified transformer structure to fully explore the nucleotide sequence
dependency. Then, we propose fixed-length parallel grouping to accelerate the
decoding speed of our autoregressive model. Experimental results on real-world
datasets show that our method saves 29.7% bit rate compared with the
state-of-the-art method, and the decoding speed is significantly faster than
all existing learning-based gene compression methods.


  To better understand DNA's 3D folding in cell nuclei, researchers developed
chromosome capture methods such as Hi-C that measure the contact frequencies
between all DNA segment pairs across the genome. As Hi-C data sets often are
massive, it is common to use bioinformatics methods to group DNA segments into
3D regions with correlated contact patterns, such as Topologically Associated
Domains (TADs) and A/B compartments. Recently, another research direction
emerged that treats the Hi-C data as a network of 3D contacts. In this
representation, one can use community detection algorithms from complex network
theory that group nodes into tightly connected mesoscale communities. However,
because Hi-C networks are so densely connected, several node partitions may
represent feasible solutions to the community detection problem but are
indistinguishable unless including other data. Because this limitation is a
fundamental property of the network, this problem persists regardless of the
community-finding or data-clustering method. To help remedy this problem, we
developed a method that charts the solution landscape of network partitions in
Hi-C data from human cells. Our approach allows us to scan seamlessly through
the scales of the network and determine regimes where we can expect reliable
community structures. We find that some scales are more robust than others and
that strong clusters may differ significantly. Our work highlights that finding
a robust community structure hinges on thoughtful algorithm design or method
cross-evaluation.


  Compartmental models are a tool commonly used in epidemiology for the
mathematical modelling of the spread of infectious diseases, with their most
popular representative being the Susceptible-Infected-Removed (SIR) model and
its derivatives. However, current SIR models are bounded in their capabilities
to model government policies in the form of non-pharmaceutical interventions
(NPIs) and weather effects and offer limited predictive power. More capable
alternatives such as agent based models (ABMs) are computationally expensive
and require specialized hardware. We introduce a neural network augmented SIR
model that can be run on commodity hardware, takes NPIs and weather effects
into account and offers improved predictive power as well as counterfactual
analysis capabilities. We demonstrate our models improvement of the
state-of-the-art modeling COVID-19 in Austria during the 03.2020 to 03.2021
period and provide an outlook for the future up to 01.2024.


  In this paper, we rigorously study the problem of cost optimisation of hybrid
(mixed) institutional incentives, which are a plan of actions involving the use
of reward and punishment by an external decision-maker, for maximising the
level (or guaranteeing at least a certain level) of cooperative behaviour in a
well-mixed, finite population of self-regarding individuals who interact via
cooperation dilemmas (Donation Game or Public Goods Game). We show that a mixed
incentive scheme can offer a more cost-efficient approach for providing
incentives while ensuring the same level or standard of cooperation in the
long-run. We establish the asymptotic behaviour (namely neutral drift, strong
selection, and infinite-population limits). We prove the existence of a phase
transition, obtaining the critical threshold of the strength of selection at
which the monotonicity of the cost function changes and providing an algorithm
for finding the optimal value of the individual incentive cost. Our analytical
results are illustrated with numerical investigations. Overall, our analysis
provides novel theoretical insights into the design of cost-efficient
institutional incentive mechanisms for promoting the evolution of cooperation
in stochastic systems.


  Machine learning methods have been used to accelerate the molecule
optimization process. However, efficient search for optimized molecules
satisfying several properties with scarce labeled data remains a challenge for
machine learning molecule optimization. In this study, we propose MOMO, a
multi-objective molecule optimization framework to address the challenge by
combining learning of chemical knowledge with Pareto-based multi-objective
evolutionary search. To learn chemistry, it employs a self-supervised codec to
construct an implicit chemical space and acquire the continues representation
of molecules. To explore the established chemical space, MOMO uses
multi-objective evolution to comprehensively and efficiently search for similar
molecules with multiple desirable properties. We demonstrate the high
performance of MOMO on four multi-objective property and similarity
optimization tasks, and illustrate the search capability of MOMO through case
studies. Remarkably, our approach significantly outperforms previous approaches
in optimizing three objectives simultaneously. The results show the
optimization capability of MOMO, suggesting to improve the success rate of lead
molecule optimization.


  Understanding the impact of collective social phenomena in epidemic dynamics
is a crucial task to effectively contain the disease spread. In this work we
build a mathematical description for assessing the interplay between opinion
polarization and the evolution of a disease. The proposed kinetic approach
describes the evolution of aggregate quantities characterizing the agents
belonging to epidemiologically relevant states, and will show that the spread
of the disease is closely related to consensus dynamics distribution in which
opinion polarization may emerge. In the present modelling framework,
microscopic consensus formation dynamics can be linked to macroscopic epidemic
trends to trigger the collective adherence to protective measures. We conduct
numerical investigations which confirm the ability of the model to describe
different phenomena related to the spread of an epidemic.


  We investigate the evolutionary dynamics of an age-structured population
under weak frequency-dependent selection. It turns out that the weak selection
is affected in a non-trivial way by the life-history trait. We can disentangle
the dynamics, based on the appearance of different time scales. These time
scales, which seem to form a universal structure in the interplay of weak
selection and life-history traits, allow us to reduce the infinite dimensional
model to a one-dimensional modified replicator equation. The modified
replicator equation is then used to investigate cooperation (the prisoner's
dilemma) by means of adaptive dynamics. We identify conditions under which age
structure is able to promote cooperation. At the end we discuss the relevance
of our findings.


  Antimicrobial peptides (AMPs) are promising therapeutic approaches against
drug-resistant pathogens. Recently, deep generative models are used to discover
new AMPs. However, previous studies mainly focus on peptide sequence attributes
and do not consider crucial structure information. In this paper, we propose a
latent sequence-structure model for designing AMPs (LSSAMP). LSSAMP exploits
multi-scale vector quantization in the latent space to represent secondary
structures (e.g. alpha helix and beta sheet). By sampling in the latent space,
LSSAMP can simultaneously generate peptides with ideal sequence attributes and
secondary structures. Experimental results show that the peptides generated by
LSSAMP have a high probability of antimicrobial activity. Our wet laboratory
experiments verified that two of the 21 candidates exhibit strong antimicrobial
activity. The code is released at https://github.com/dqwang122/LSSAMP.


  Aging is a multifactorial process and a key factor of morbidity and
mortality. Alzheimer's disease (AD) is an age-related disorder and a main cause
of worldwide disability. Both aging and AD can be characterized by metabolic
dysfunction. Metabolomics can quantify the complete set of metabolites in a
studied sample and is helpful for studying metabolic alterations in aging and
AD. In this review, we summarize the metabolomic changes regarding aging and
AD, discuss their biological functions, and highlight their potential
application as diagnostic biomarkers or therapeutic targets. Recent advances in
multi-omics approaches for understanding the metabolic mechanism of aging and
AD are also reviewed.


  Environmental variations can significantly influence how populations compete
for resources, and hence shape their evolution. Here, we study population
dynamics subject to a fluctuating environment modeled by a varying carrying
capacity changing continuously in time according to either binary random
switches, or by being driven by a noise of continuous range. We consider a
prototypical example of two competing strains, one growing slightly slower than
the other, and consider also the scenario where the slow strain is a public
goods producer. By systematically comparing the effect of binary- versus
continuously-varying environment, we study how different noise statistics
(mean, variance) influence the population size and fixation properties. We show
that the slow strain fixation probability can be greatly enhanced for a
continuously-varying environment compared to binary switches, even when the
first two moments of the carrying capacity coincide.


  Biomolecular condensates have been shown to play a fundamental role in
localizing biochemistry in a cell. RNA is a common constituent of condensates,
and can determine their biophysical properties. Functions of biomolecular
condensates are varied including activating, inhibiting, and localizing
reactions. Recent theoretical work has shown that the phase separation of
proteins into droplets can diminish cell to cell variability in protein
abundance. However, the extent to which phase separation involving mRNAs may
also buffer noise has yet to be explored. In this paper, we introduce a
phenomenological model for the phase separation of mRNAs into RNP condensates,
and quantify noise suppression as a function of gene expression kinetic
parameters. Through stochastic simulations, we highlight the ability for
condensates formed from just a handful of mRNAs to regulate the abundance and
suppress the fluctuations of proteins. We place particular emphasis on how this
mechanism can facilitate efficient transcription by reducing noise even in the
situation of infrequent bursts of transcription by exploiting the physics of a
concentration-dependent, deterministic phase separation threshold. We
investigate two biologically relevant models in which phase separation acts to
either "buffer" noise by storing mRNA in inert droplets, or "filter" phase
separated mRNAs by accelerating their decay, and quantify expression noise as a
function of kinetic parameters. In either case the most efficient expression
occurs when bursts produce mRNAs close the phase separation threshold, which we
find to be broadly consistent with observations of an RNP-droplet forming
cyclinin multinucleate Ashbya gossypii cells. We finally consider the
contribution of noise in the phase separation threshold, and show that protein
copy number noise can be suppressed by phase separation threshold fluctuations
in certain conditions.


  A long-standing goal of machine-learning-based protein engineering is to
accelerate the discovery of novel mutations that improve the function of a
known protein. We introduce a sampling framework for evolving proteins in
silico that supports mixing and matching a variety of unsupervised models, such
as protein language models, and supervised models that predict protein function
from sequence. By composing these models, we aim to improve our ability to
evaluate unseen mutations and constrain search to regions of sequence space
likely to contain functional proteins. Our framework achieves this without any
model fine-tuning or re-training by constructing a product of experts
distribution directly in discrete protein space. Instead of resorting to brute
force search or random sampling, which is typical of classic directed
evolution, we introduce a fast MCMC sampler that uses gradients to propose
promising mutations. We conduct in silico directed evolution experiments on
wide fitness landscapes and across a range of different pre-trained
unsupervised models, including a 650M parameter protein language model. Our
results demonstrate an ability to efficiently discover variants with high
evolutionary likelihood as well as estimated activity multiple mutations away
from a wild type protein, suggesting our sampler provides a practical and
effective new paradigm for machine-learning-based protein engineering.


  We develop computational methods that incorporate shear into fluctuating
hydrodynamics methods. We are motivated by the rheological responses of complex
fluids and soft materials. Our approach is based on continuum stochastic
hydrodynamic equations that are subject to shear boundary conditions on the
unit periodic cell in a manner similar to the Lees-Edwards conditions of
molecular dynamics. Our methods take into account consistently the
microstructure elastic mechanics, fluid-structure hydrodynamic coupling, and
thermal fluctuations. For practical simulations, we develop numerical methods
for efficient stochastic field generation that handle the sheared generalized
periodic boundary conditions. We show that our numerical methods are consistent
with fluctuation dissipation balance and near-equilibrium statistical
mechanics. As a demonstration in practice, we present several prototype
rheological response studies. These include (i) shear thinning of a polymeric
fluid, (ii) complex moduli for the oscillatory responses of a polymerized lipid
vesicle, and (iii) aging under shear of a gel-like material.


  Ordinary differential equation (ODE) models are widely used to describe
chemical or biological processes. This article considers the estimation and
assessment of such models on the basis of time-course data. Due to experimental
limitations, time-course data are often noisy and some components of the system
may not be observed. Furthermore, the computational demands of numerical
integration have hindered the widespread adoption of time-course analysis using
ODEs. To address these challenges, we explore the efficacy of the recently
developed MAGI (MAnifold-constrained Gaussian process Inference) method for ODE
inference. First, via a range of examples we show that MAGI is capable of
inferring the parameters and system trajectories, including unobserved
components, with appropriate uncertainty quantification. Second, we illustrate
how MAGI can be used to assess and select different ODE models with time-course
data based on MAGI's efficient computation of model predictions. Overall, we
believe MAGI is a useful method for the analysis of time-course data in the
context of ODE models, which bypasses the need for any numerical integration.


  The stochastic kinetics of BRN are described by a chemical master equation
(CME) and the underlying laws of mass action. The CME must be usually solved
numerically by generating enough traces of random reaction events. The
resulting event-time series can be evaluated statistically to identify, for
example, the reaction clusters, rare reaction events, and the periods of
increased or steady-state activity. The aim of this paper is to newly exploit
the empirical statistics of the reaction events in order to obtain causally and
anti-causally related sub-sequences of reactions. This allows discovering some
of the causal dynamics of the reaction networks as well as uncovering their
more deterministic behaviors. In particular, it is proposed that the reaction
sub-sequences that are conditionally nearly certain or nearly uncertain can be
considered as being causally related or unrelated, respectively. Moreover,
since time-ordering of reactions is locally irrelevant, the reaction
sub-sequences can be transformed into the reaction event sets or multi-sets.
The appropriately defined distance metrics can be then used to define
equivalences between the reaction sub-sequences. The proposed framework for
identifying causally associated reaction sub-sequences has been implemented as
a computationally efficient query-response mechanism. The framework was
evaluated assuming five selected models of genetic reaction networks in seven
defined numerical experiments. The models were simulated in BioNetGen using
NFsim, which had to be modified to allow recording of the traces of reaction
events. The generated event time-series were analyzed by Python and Matlab
scripts. The whole process of data generation, analysis and visualization has
been nearly fully automated using shell scripts.


  Generalized mass-action systems are power-law dynamical systems arising from
chemical reaction networks. Essentially, every nonnegative ODE model used in
chemistry and biology (for example, in ecology and epidemiology) and even in
economics and engineering can be written in this form. Previous results have
focused on existence and uniqueness of special steady states (complex-balanced
equilibria) for all rate constants, thereby ruling out multiple (special)
steady states. Recently, necessary conditions for linear stability have been
obtained.
  In this work, we provide sufficient conditions for the linear stability of
complex-balanced equilibria for all rate constants (and also for the
non-existence of other steady states). In particular, via sign-vector
conditions (on the stoichiometric coefficients and kinetic orders), we
guarantee that the Jacobian matrix is a $P$-matrix. Technically, we use a new
decomposition of the graph Laplacian which allows to consider orders of
(generalized) monomials. Alternatively, we use cycle decomposition which allows
a linear parametrization of all Jacobian matrices. In any case, we guarantee
stability without explicit computation of steady states. We illustrate our
results in examples from chemistry and biology: generalized Lotka-Volterra
systems and SIR models, a two-component signaling system, and an enzymatic
futile cycle.


  This paper deals with the existence, multiplicity, minimal complexity and
global structure of the subharmonic solutions to a class of planar Hamiltonian
systems with periodic coefficients, being the classical predator-prey model of
V. Volterra its most paradigmatic example. By means of a topological approach
based on techniques from global bifurcation theory, the first part of the paper
ascertains their nature, multiplicity and minimal complexity, as well as their
global minimal structure, in terms of the configuration of the function
coefficients in the setting of the model. The second part of the paper
introduces a dynamical system approach based on the theory of topological
horseshoes that permits to detect, besides subharmonic solutions,
``chaotic-type'' solutions. As a byproduct of our analysis, the simplest
predator-prey prototype models in periodic environments can provoke chaotic
dynamics. This cannot occur in cooperative and quasi-cooperative dynamics, as a
consequence of the ordering imposed by the maximum principle.


  Machine learning methods have seen increased application to geospatial
environmental problems, such as precipitation nowcasting, haze forecasting, and
crop yield prediction. However, many of the machine learning methods applied to
mosquito population and disease forecasting do not inherently take into account
the underlying spatial structure of the given data. In our work, we apply a
spatially aware graph neural network model consisting of GraphSAGE layers to
forecast the presence of West Nile virus in Illinois, to aid mosquito
surveillance and abatement efforts within the state. More generally, we show
that graph neural networks applied to irregularly sampled geospatial data can
exceed the performance of a range of baseline methods including logistic
regression, XGBoost, and fully-connected neural networks.


  In this letter we consider a single parameter generalization of the standard
three species Rock-Paper-Scissors (RPS) model allowing for predator-prey
reversal. This model, which shall be referred to as $\kappa$RPS model,
incorporates bidirectional predator-prey interactions between all the species
in addition to the unidirectional predator-prey interactions of the standard
RPS model. We study the dynamics of a May-Leonard formulation of the
$\kappa$RPS model using lattice based spatial stochastic simulations with
random initial conditions. We find that if the simulation lattices are
sufficiently large for the coexistence of all three species to be maintained,
the model asymptotically leads to the formation of spiral patterns whose
evolution is qualitatively similar to that of the standard RPS model, albeit
with larger characteristic length and time scales. We show that there are in
general two distinct scaling regimes: one transient curvature dominated regime
in which the characteristic length of the population network grows with time
and another where it becomes a constant. We also estimate the dependence of the
asymptotic value of the characteristic length of the population network on the
likelihood of predator-prey reversal and show that if the simulation lattices
are not sufficiently large then predator-prey reversal can have a significant
negative impact on coexistence. Finally, we interpret these results by
considering the much simpler dynamics of circular domains.


  Retrosynthesis is the task of breaking down a chemical compound recursively
step-by-step into molecular precursors until a set of commercially available
molecules is found. Consequently, the goal is to provide a valid synthesis
route for a molecule. As more single-step models develop, we see increasing
accuracy in the prediction of molecular disconnections, potentially improving
the creation of synthetic paths. Multi-step approaches repeatedly apply the
chemical information stored in single-step retrosynthesis models. However, this
connection is not reflected in contemporary research, fixing either the
single-step model or the multi-step algorithm in the process. In this work, we
establish a bridge between both tasks by benchmarking the performance and
transfer of different single-step retrosynthesis models to the multi-step
domain by leveraging two common search algorithms, Monte Carlo Tree Search and
Retro*. We show that models designed for single-step retrosynthesis, when
extended to multi-step, can have a tremendous impact on the route finding
capabilities of current multi-step methods, improving performance by up to +30%
compared to the most widely used model. Furthermore, we observe no clear link
between contemporary single-step and multi-step evaluation metrics, showing
that single-step models need to be developed and tested for the multi-step
domain and not as an isolated task to find synthesis routes for molecules of
interest.


  We present a novel method for the automatic detection of pockets on protein
molecular surfaces. The algorithm is based on an ad hoc hierarchical clustering
of virtual SES probe spheres obtained from the geometrical primitives generated
by the NanoShaper software. The final ranking of putative pockets is based on
the Isolation Forest method, an unsupervised learning approach originally
developed for anomaly detection. A detailed importance analysis of pocket
features provides insight on which geometrical (clustering) and chemical
(residues) properties characterize a good binding site. The method also
provides a segmentation of pockets into smaller subpockets. We prove that
subpockets are a reliable representation that pinpoint the binding site with
greater precision. Site Ferret is outstanding in its versatility, accurately
predicting a wide range of binding sites, from small molecules to peptides and
difficult shallow sites.


  B-cell repertoires are characterized by a diverse set of receptors of
distinct specificities generated through two processes of somatic
diversification: V(D)J recombination and somatic hypermutations. B cell clonal
families stem from the same V(D)J recombination event, but differ in their
hypermutations. Clonal families identification is key to understanding B-cell
repertoire function, evolution and dynamics. We present HILARy (High-precision
Inference of Lineages in Antibody Repertoires), an efficient, fast and precise
method to identify clonal families from high-throughput sequencing datasets.
HILARy combines probabilistic models that capture the receptor generation and
selection statistics with adapted clustering methods to achieve consistently
high inference accuracy. It automatically leverages the phylogenetic signal of
shared mutations in difficult repertoire subsets. Exploiting the high
sensitivity of the method, we find the statistics of evolutionary properties
such as the site frequency spectrum and dN/dS ratio do not depend on the
junction length. We also identify a broad range of selection pressures scanning
two orders of magnitude.


  Applying deep learning concepts from image detection and graph theory has
greatly advanced protein-ligand binding affinity prediction, a challenge with
enormous ramifications for both drug discovery and protein engineering. We
build upon these advances by designing a novel deep learning architecture
consisting of a 3-dimensional convolutional neural network utilizing
channel-wise attention and two graph convolutional networks utilizing
attention-based aggregation of node features. HAC-Net (Hybrid Attention-Based
Convolutional Neural Network) obtains state-of-the-art results on the PDBbind
v.2016 core set, the most widely recognized benchmark in the field. We
extensively assess the generalizability of our model using multiple train-test
splits, each of which maximizes differences between either protein structures,
protein sequences, or ligand extended-connectivity fingerprints of complexes in
the training and test sets. Furthermore, we perform 10-fold cross-validation
with a similarity cutoff between SMILES strings of ligands in the training and
test sets, and also evaluate the performance of HAC-Net on lower-quality data.
We envision that this model can be extended to a broad range of supervised
learning problems related to structure-based biomolecular property prediction.
All of our software is available as open source at
https://github.com/gregory-kyro/HAC-Net/, and the HACNet Python package is
available through PyPI.


  Information of relevance to the organism often is represented by the
concentrations of particular molecules inside a cell. As outside observers we
can now measure these concentrations precisely, but the cell's own mechanisms
must be noisier. We show that, in certain limits, there is a universal tradeoff
between the information capacity of these concentration measurements and the
amount of relevant information that is captured, a version of the information
bottleneck problem. This universal tradeoff is confirmed, quantitatively, in an
analysis of the positional information encoded by the gap genes in the
developing fly embryo.


  Over the past decade, neural networks have been successful at making
predictions from biological sequences, especially in the context of regulatory
genomics. As in other fields of deep learning, tools have been devised to
extract features such as sequence motifs that can explain the predictions made
by a trained network. Here we intend to go beyond explainable machine learning
and introduce SEISM, a selective inference procedure to test the association
between these extracted features and the predicted phenotype. In particular, we
discuss how training a one-layer convolutional network is formally equivalent
to selecting motifs maximizing some association score. We adapt existing
sampling-based selective inference procedures by quantizing this selection over
an infinite set to a large but finite grid. Finally, we show that sampling
under a specific choice of parameters is sufficient to characterize the
composite null hypothesis typically used for selective inference-a result that
goes well beyond our particular framework. We illustrate the behavior of our
method in terms of calibration, power and speed and discuss its power/speed
trade-off with a simpler data-split strategy. SEISM paves the way to an easier
analysis of neural networks used in regulatory genomics, and to more powerful
methods for genome wide association studies (GWAS).


  Regulatory networks describe the interactions between molecular or cellular
regulators, like transcription factors and genes in gene regulatory networks,
kinases and their receptors in signalling networks, or neurons in neural
networks. A long-standing aim of quantitative biology is to reconstruct such
networks on the basis of large-scale data. Our aim is to leverage fluctuations
around the non-equilibrium steady state for network inference. To this end, we
use a stochastic model of gene regulation or neural dynamics and solve it
approximately within a Gaussian mean-field theory. We develop a likelihood
estimate based on this stochastic theory to infer regulatory interactions from
perturbation data on the network nodes. We apply this approach to artificial
perturbation data as well as to phospho-proteomic data from cell-line
experiments and compare our results to inference schemes restricted to mean
activities in the steady state.


  Neoplasms (NPs) and neurological diseases and disorders (NDDs) are amongst
the major classes of diseases underlying deaths of a disproportionate number of
people worldwide. To determine if there exist some distinctive features in the
local wiring patterns of protein interactions emerging at the onset of a
disease belonging to either of these two classes, we examined 112 and 175
protein interaction networks belonging to NPs and NDDs, respectively. Orbit
usage profiles (OUPs) for each of these networks were enumerated by
investigating the networks' local topology. 56 non-redundant OUPs (nrOUPs) were
derived and used as network features for classification between these two
disease classes. Four machine learning classifiers, namely, k-nearest neighbour
(KNN), support vector machine (SVM), deep neural network (DNN), random forest
(RF) were trained on these data. DNN obtained the greatest average AUPRC
(0.988) among these classifiers. DNNs developed on node2vec and the proposed
nrOUPs embeddings were compared using 5-fold cross validation on the basis of
average values of the six of performance measures, viz., AUPRC, Accuracy,
Sensitivity, Specificity, Precision and MCC. It was found that nrOUPs based
classifier performed better in all of these six performance measures.


  Artificial intelligence (AI) in the form of deep learning bears promise for
drug discovery and chemical biology, $\textit{e.g.}$, to predict protein
structure and molecular bioactivity, plan organic synthesis, and design
molecules $\textit{de novo}$. While most of the deep learning efforts in drug
discovery have focused on ligand-based approaches, structure-based drug
discovery has the potential to tackle unsolved challenges, such as affinity
prediction for unexplored protein targets, binding-mechanism elucidation, and
the rationalization of related chemical kinetic properties. Advances in deep
learning methodologies and the availability of accurate predictions for protein
tertiary structure advocate for a $\textit{renaissance}$ in structure-based
approaches for drug discovery guided by AI. This review summarizes the most
prominent algorithmic concepts in structure-based deep learning for drug
discovery, and forecasts opportunities, applications, and challenges ahead.


  In this paper, we develop a discrete time stochastic model under partial
information to explain the evolution of Covid-19 pandemic. Our model is a
modification of the well-known SIR model for epidemics, which accounts for some
peculiar features of Covid-19. In particular, we work with a random
transmission rate and we assume that the true number of infectious people at
any observation time is random and not directly observable, to account for
asymptomatic and non-tested people. We elaborate a nested particle filtering
approach to estimate the reproduction rate and the model parameters. We apply
our methodology to Austrian Covid-19 infection data in the period from May 2020
to June 2022. Finally, we discuss forecasts and model tests.