  Pathway enrichment analysis has become a widely used knowledge-based approach
for the interpretation of biomedical data. Its popularity has led to an
explosion of both enrichment methods and pathway databases. While the elegance
of pathway enrichment lies in its simplicity, multiple factors can impact the
results of such an analysis which may not be accounted for. Researchers may
fail to give influential aspects their due, resorting instead to popular
methods and gene set collections, or default settings. Despite ongoing efforts
to establish set guidelines, meaningful results are still hampered by a lack of
consensus or gold standards around how enrichment analysis should be conducted.
Nonetheless, such concerns have prompted a series of benchmark studies
specifically focused on evaluating the influence of various factors on pathway
enrichment results. In this review, we organize and summarize the findings of
these benchmarks to provide a comprehensive overview on the influence of these
factors. Our work covers a broad spectrum of factors, spanning from
methodological assumptions to those related to prior biological knowledge, such
as pathway definitions and database choice. In doing so, we aim to shed light
on how these aspects can lead to insignificant, uninteresting, or even
contradictory results. Finally, we conclude the review by proposing future
benchmarks as well as solutions to overcome some of the challenges which
originate from the outlined factors.


  DNA is strong polyelectrolyte macromolecule making metal ions (counterions)
condense to a cloud around the double helix. The counterions may be localized
outside the macromolecule and inside the minor and major grooves of the double
helix. In the present work, the distribution of condensed counterions between
inner and outer regions of DNA has been studied using the approaches of
counterion condensation theory. The results have shown that the number of
counterions trapped inside the macromolecule should be greater than 0.16 per
one phosphate group. The maximal number of counterions that may be localized
inside the DNA double helix is limited to about 0.4 per one phosphate group and
it is much lower than the total number of condensed counterions. To analyze the
structure of counterion cloud the molecular dynamics simulations of
\emph{B}-DNA with K$^{+}$ counterions have been performed. The obtained number
of the counterions trapped inside the grooves of the double helix is about
0.22$\pm$0.06 per one phosphate group that agree with the model estimations.
The developed model describes general features of the structure of counterion
cloud around DNA and is able to predict the number of counterions inside the
grooves of the double helix.


  Increasing interest in the acquisition of biotic and abiotic resources from
within the deep sea (e.g. fisheries, oil-gas extraction, and mining) urgently
imposes the development of novel monitoring technologies, beyond the
traditional vessel-assisted, time-consuming, high-cost sampling surveys. The
implementation of permanent networks of seabed and water-column cabled (fixed)
and docked mobile platforms is presently enforced, to cooperatively measure
biological features and environmental (physico-chemical) parameters. Video and
acoustic (i.e. optoacoustic) imaging are becoming central approaches for
studying benthic fauna (e.g. quantifying species presence, behaviour, and
trophic interactions) in a remote, continuous, and prolonged fashion. Imaging
is also being complemented by in situ environmental-DNA sequencing
technologies, allowing the traceability of a wide range of organisms (including
prokaryotes) beyond the reach of optoacoustic tools. Here, we describe the
different fixed and mobile platforms of those benthic and pelagic monitoring
networks, proposing at the same time an innovative roadmap for the automated
computing of hierarchical ecological information of deep-sea ecosystems (i.e.
from single species abundance and life traits, to community composition, and
overall biodiversity)


  Phenomenological and deterministic models are often used for the estimation
of transmission parameters in an epidemic and for the prediction of its growth
trajectory. Such analyses are usually based on single peak outbreak dynamics.
In light of the present COVID-19 pandemic, there is a pressing need to better
understand observed epidemic growth with multiple peak structures, preferably
using first-principles methods. Along the lines of our previous work [Physica A
574, 126014 (2021)], here we apply 2D random-walk Monte Carlo calculations to
better understand COVID-19 spread through contact interactions. Lockdown
scenarios and all other control interventions are imposed through mobility
restrictions and a regulation of the infection rate within the stochastically
interacting population. The susceptible, infected and recovered populations are
tracked over time, with daily infection rates obtained without recourse to the
solution of differential equations.
  The simulations were carried out for population densities corresponding to
four countries, India, Serbia, South Africa and USA. In all cases our results
capture the observed infection growth rates. More importantly, the simulation
model is shown to predict secondary and tertiary waves of infections with
reasonable accuracy. This predictive nature of multiple wave structures
provides a simple and effective tool that may be useful in planning mitigation
strategies during the present pandemic.


  We propose a new deterministic methodology to predict RNA sequence and
protein folding. Is stem enough for structure prediction? The main idea is to
consider all possible stem formation in the given sequence. With the stem loop
energy and the strength of stem, we explore how to deterministically utilize
stem information for RNA sequence and protein folding structure prediction. We
use graph notation, where all possible stems are represented as vertices, and
co-existence as edges. This full Stem-graph presents all possible folding
structure, and we pick sub-graph(s) which give the best matching energy for
folding structure prediction. We introduce a Stem-Loop score to add structure
information and to speed up the computation. The proposed method can handle
secondary structure prediction as well as protein folding with pseudo knots.
Numerical experiments are done using a laptop and results take only a few
minutes or seconds. One of the strengths of this approach is in the simplicity
and flexibility of the algorithm, and it gives deterministic answer. We explore
protein sequences from Protein Data Bank, rRNA 5S sequences, and tRNA sequences
from the Gutell Lab. Various experiments and comparisons are included to
validate the propose method.


  A genome read data set can be quickly and efficiently remapped from one
reference to another similar reference (e.g., between two reference versions or
two similar species) using a variety of tools, e.g., the commonly-used CrossMap
tool. With the explosion of available genomic data sets and references,
high-performance remapping tools will be even more important for keeping up
with the computational demands of genome assembly and analysis.
  We provide FastRemap, a fast and efficient tool for remapping reads between
genome assemblies. FastRemap provides up to a 7.82$\times$ speedup
(6.47$\times$, on average) and uses as low as 61.7% (80.7%, on average) of the
peak memory consumption compared to the state-of-the-art remapping tool,
CrossMap.
  FastRemap is written in C++. The source code and user manual are freely
available at: github.com/CMU-SAFARI/FastRemap. Docker image available at:
https://hub.docker.com/r/alkanlab/fast. Also available in Bioconda at:
https://anaconda.org/bioconda/fastremap-bio.


  Since December 2020, variants of COVID-19 (especially Delta and Omicron)
appeared with different characteristics that influenced death and
transmissibility emerged around the world. To address the novel dynamics of the
disease, we propose a dynamical model of two strains, namely native and mutant,
transmission dynamics with mutation and imperfect vaccination. It is also
assumed that the recuperated individuals from the native strain can be infected
with mutant strain through the direct contact with individual or contaminated
surfaces or aerosols. We compute the basic reproduction number for each strain
independently and take the maximum for $R_0$. We prove the nonexistence of
backward bifurcation using the center manifold theory, and global stability of
disease-free equilibrium when the basic reproduction number $R_0<1. An
intermediate mutation rate $\nu_1$ leads to oscillations. When $\nu_1$
increases over a threshold, the system regains its stability and exhibits an
interesting dynamics called endemic bubble. An analytical expression for
vaccine-induced herd immunity is derived. The model is parameterized using the
Indian data of the cumulative number of confirmed cases and deaths of COVID-19
from March 1 to September 27 in 2021, using MCMC method. The cumulative cases
and deaths can be reduced by increasing the vaccine efficacies to both native
and mutant strains. We observe that by considering the vaccine efficacy to
native strain as 90\%, the cumulative cases and deaths would be reduced by
3.27\% and 5.2\%, respectively; and by considering the vaccine efficacy to
mutant strain as 90\%, the cumulative cases and deaths would be reduced by
0.9\% and 2.5\%, respectively. Our study demonstrates that the COVID-19
pandemic may be worse due to the occurrence of oscillations for certain
mutation rates but better due to stability at a lower infection level with a
larger mutation rate.


  We assess costs and efficiency of state-of-the-art high performance cloud
computing compared to a traditional on-premises compute cluster. Our use case
are atomistic simulations carried out with the GROMACS molecular dynamics (MD)
toolkit with a focus on alchemical protein-ligand binding free energy
calculations.
  We set up a compute cluster in the Amazon Web Services (AWS) cloud that
incorporates various different instances with Intel, AMD, and ARM CPUs, some
with GPU acceleration. Using representative biomolecular simulation systems we
benchmark how GROMACS performs on individual instances and across multiple
instances. Thereby we assess which instances deliver the highest performance
and which are the most cost-efficient ones for our use case.
  We find that, in terms of total costs including hardware, personnel, room,
energy and cooling, producing MD trajectories in the cloud can be as
cost-efficient as an on-premises cluster given that optimal cloud instances are
chosen. Further, we find that high-throughput ligand-screening for
protein-ligand binding affinity estimation can be accelerated dramatically by
using global cloud resources. For a ligand screening study consisting of 19,872
independent simulations, we used all hardware that was available in the cloud
at the time of the study. The computations scaled-up to reach peak performances
using more than 4,000 instances, 140,000 cores, and 3,000 GPUs simultaneously
around the globe. Our simulation ensemble finished in about two days in the
cloud, while weeks would be required to complete the task on a typical
on-premises cluster consisting of several hundred nodes. We demonstrate that
the costs of such and similar studies can be drastically reduced with a
checkpoint-restart protocol that allows to use cheap Spot pricing and by using
instance types with optimal cost-efficiency.


  How do species coexist? A framework known as Modern Coexistence Theory
measures mechanisms of coexistence by comparing a species perturbed to low
density (the invader) to other species that remain at their typical densities
(the residents); this invader-resident comparison measures a rare-species
advantage that results from specialization. However, there are several
reasonable ways (i.e., methods) to compare invaders and residents, each
differing in practicality and biological interpretation. Here, using
theoretical arguments and case studies, we compare four such methods for
calculating coexistence mechanisms: 1) Scaling factors, the traditional
approach where resident growth rates are scaled by a measure of relative
sensitivity to competition, obtained by solving a system of linear equations;
2) The simple comparison, which gives equal weight to all resident species; 3)
Speed conversion factors, a novel method in which resident growth rates are
scaled by a ratio of generation times, and; 4) The invader-invader comparison,
another novel method in which a focal species is compared to itself at high vs.
low density. We conclude that the conventional scaling factors can be useful in
some theoretical research, but are not recommended for empirical applications,
i.e., determining the mechanisms of coexistence in real communities. Instead,
we recommend the simple comparison and speed conversion factor methods. The
speed conversion factors are most useful when comparing species with dissimilar
generation times. However, ecologists often study coexistence in guilds of
species with similar life-histories, and therefore, similar generation times.
In such scenarios, the easier-to-use simple comparison method is reasonable.


  The storage effect is a well-known explanation for coexistence in temporally
varying environments. Like many complex ecological theories, the storage effect
is often used as an explanation for observed coexistence on the basis of
heuristic understanding, rather than careful application of a detailed model.
One interpretation states that species coexist by specializing on specific
environmental states, and therefore must have a robust life-stage (e.g.,
long-lived adults, a seed-bank) in order to "wait it out" for favorable
conditions. Here we show that this widely employed interpretation can be
misleading. Multiple models show that stage-structure, long lifespans, and
overlapping generations are neither necessary nor sufficient for the storage
effect. In models where a robust life-stage does engender a storage effect, it
does not do so by preventing stochastic extinction or by improving relative
bet-hedging. A robust life-stage is best understood as one of many ways to
fulfill an abstract condition for the storage effect: an interaction effect of
environment and competition on per capita growth rates. Using a dataset of
annual plants from a Mediterranean grassland in Spain, we show that such
interaction effects occur between water availability and the number of
germinant competitors, leading to storage in the absence of a persistent seed
bank. Empiricists hoping to uncover the storage effect should look for
interaction effects between environmental conditions and competition -- easily
identifiable with multiple regression -- at all stages of a species'
life-cycle.


  The storage effect is a general explanation for coexistence in a variable
environment. The generality of the storage effect is both a strength - it can
be quantified in many systems - and a challenge - there is not a clear
relationship between the abstract conditions for storage effect and species'
life-history traits (e.g., dormancy, stage-structure, non-overlapping
generations), thus precluding a simple ecological interpretation of the storage
effect. Our goal here is to provide a clearer understanding of the conditions
for the storage effect as a step towards a better general explanation for
coexistence in a variable environment. Our approach focuses on dividing one of
the key conditions for the storage effect, covariance between environment and
competition, into two pieces, namely that there must be a causal relationship
between environment and competition, and that the effects of the environment do
not change too quickly. This finer-grained definition can explain a number of
previous results, including 1) that the storage effect promotes annual plant
coexistence when the germination rate fluctuates, but not when the seed yield
fluctuates, 2) that the storage effect is more likely to be induced by resource
competition than apparent competition, and 3) that the spatial storage effect
is more probable than the temporal storage effect. Additionally, our expanded
definition suggests two novel mechanisms by which the temporal storage effect
can arise: transgenerational plasticity, and causal chains of environmental
variables. These mechanisms produce coexistence via the storage effect without
any need for stage structure or a temporally autocorrelated environment.


  This paper focuses on and analyzes realistic SIR models that take
stochasticity into account. The proposed systems are applicable to most
incidence rates that are used in the literature including the bilinear
incidence rate, the Beddington-DeAngelis incidence rate, and a Holling type II
functional response. Given that many diseases can lead to asymptomatic
infections, we look at a system of stochastic differential equations that also
includes a class of hidden state individuals, for which the infection status is
unknown. We assume that the direct observation of the percentage of hidden
state individuals that are infected, $\alpha(t)$, is not given and only a
noise-corrupted observation process is available. Using the nonlinear filtering
techniques in conjunction with an invasion type analysis (or analysis using
Lyapunov exponents from the dynamical system point of view), this paper proves
that the long-term behavior of the disease is governed by a threshold
$\lambda\in \mathbb{R}$ that depends on the model parameters. It turns out that
if $\lambda<0$ the number $I(t)$ of infected individuals converges to zero
exponentially fast, or the extinction happens. In contrast, if $\lambda>0$, the
infection is endemic and the system is permanent. We showcase our results by
applying them in specific illuminating examples. Numerical simulations are also
given to illustrate our results.


  Recently developed technologies to generate single-cell genomic data have
made a revolutionary impact in the field of biology. Multi-omics assays offer
even greater opportunities to understand cellular states and biological
processes. However, the problem of integrating different -omics data with very
different dimensionality and statistical properties remains quite challenging.
A growing body of computational tools are being developed for this task,
leveraging ideas ranging from machine translation to the theory of networks and
representing a new frontier on the interface of biology and data science. Our
goal in this review paper is to provide a comprehensive, up-to-date survey of
computational techniques for the integration of multi-omics and alignment of
multiple modalities of genomics data in the single cell research field.


  Indirect reciprocity is a mechanism by which individuals cooperate with those
who have cooperated with others. This creates a regime in which repeated
interactions are not necessary to incent cooperation (as would be required for
direct reciprocity). However, indirect reciprocity creates a new problem: how
do agents know who has cooperated with others? To know this, agents would need
to access some form of reputation information. Perhaps there is a communication
system to disseminate reputation information, but how does it remain truthful
and informative? Most papers assume the existence of a truthful, forthcoming,
and informative communication system; in this paper, we seek to explain how
such a communication system could remain evolutionarily stable in the absence
of exogenous pressures. Specifically, we present three conditions that together
maintain both the truthfulness of the communication system and the prevalence
of cooperation: individuals (1) use a norm that rewards the behaviors that it
prescribes (an aligned norm), (2) can signal not only about the actions of
other agents, but also about their truthfulness (by acting as third party
observers to an interaction), and (3) make occasional mistakes, demonstrating
how error can create stability by introducing diversity.


  Bivalency confers several concentration-dependent phenomena, including
avidity, competitive exchange and multi-site competitive exchange. Since these
concepts are crucial for a wide variety of topics in cell and molecular
biology, their extension, modification and/or re-purposing is also increasingly
important for the design and construction of de-novo synthetic systems at the
nanoscale. In this context, we draw upon classical techniques of statistical
physics to revisit bivalency, highlighting that receptor site geometry offers a
design modality independent of the chemistry of the individual binding
interfaces themselves. Recasting the problem in terms of many-body
coordination, we explore extended, translationally-invariant chains and
lattices of receptor sites. This not only brings clarity to behaviours
associated with simpler motifs, but also enables us to distil core principles
for the rational design of concentration-dependent kinetics in synthetic
soft-systems, which centre on the notion of geometric frustration. In doing so,
we also reveal the possibility of other tunable spatio-temporal features, such
as correlation lengths, mean-squared displacements and percolation-like
transitions.


  SIRS models capture transmission dynamics of infectious diseases for which
immunity is not lifelong. Extending these models by a W compartment for
individuals with waning immunity, the boosting of the immune system upon
repeated exposure may be incorporated. Previous analyses assumed identical
waning rates from R to W and from W to S. This implicitly assumes equal length
for the period of full immunity and of waned immunity. We relax this
restriction, and allow an asymmetric partitioning of the total immune period.
Stability switches of the endemic equilibrium are investigated with a
combination of analytic and numerical tools. Then, continuation methods are
applied to track bifurcations along the equilibrium branch. We find rich
dynamics: Hopf bifurcations, endemic double bubbles, and regions of
bistability. Our results highlight that the length of the period in which
waning immunity can be boosted is a crucial parameter significantly influencing
long term epidemiological dynamics.


  Inferring dependencies between complex biological traits while accounting for
evolutionary relationships between specimens is of great scientific interest
yet remains infeasible when trait and specimen counts grow large. The
state-of-the-art approach uses a phylogenetic multivariate probit model to
accommodate binary and continuous traits via a latent variable framework, and
utilizes an efficient bouncy particle sampler (BPS) to tackle the computational
bottleneck -- integrating many latent variables from a high-dimensional
truncated normal distribution. This approach breaks down as the number of
specimens grows and fails to reliably characterize conditional dependencies
between traits. Here, we propose an inference pipeline for phylogenetic probit
models that greatly outperforms BPS. The novelty lies in 1) a combination of
the recent Zigzag Hamiltonian Monte Carlo (Zigzag-HMC) with linear-time
gradient evaluations and 2) a joint sampling scheme for highly correlated
latent variables and correlation matrix elements. In an application exploring
HIV-1 evolution from 535 viruses, the inference requires joint sampling from an
11,235-dimensional truncated normal and a 24-dimensional covariance matrix. Our
method yields a 5-fold speedup compared to BPS and makes it possible to learn
partial correlations between candidate viral mutations and virulence.
Computational speedup now enables us to tackle even larger problems: we study
the evolution of influenza H1N1 glycosylations on around 900 viruses. For
broader applicability, we extend the phylogenetic probit model to incorporate
categorical traits, and demonstrate its use to study Aquilegia flower and
pollinator co-evolution.


  The 21st century is presenting humankind with unprecedented environmental and
medical challenges. The ability to design novel proteins tailored for specific
purposes could transform our ability to respond timely to these issues. Recent
advances in the field of artificial intelligence are now setting the stage to
make this goal achievable. Protein sequences are inherently similar to natural
languages: Amino acids arrange in a multitude of combinations to form
structures that carry function, the same way as letters form words and
sentences that carry meaning. Therefore, it is not surprising that throughout
the history of Natural Language Processing (NLP), many of its techniques have
been applied to protein research problems. In the last few years, we have
witnessed revolutionary breakthroughs in the field of NLP. The implementation
of Transformer pre-trained models has enabled text generation with human-like
capabilities, including texts with specific properties such as style or
subject. Motivated by its considerable success in NLP tasks, we expect
dedicated Transformers to dominate custom protein sequence generation in the
near future. Finetuning pre-trained models on protein families will enable the
extension of their repertoires with novel sequences that could be highly
divergent but still potentially functional. The combination of control tags
such as cellular compartment or function will further enable the controllable
design of novel protein functions. Moreover, recent model interpretability
methods will allow us to open the 'black box' and thus enhance our
understanding of folding principles. While early initiatives show the enormous
potential of generative language models to design functional sequences, the
field is still in its infancy. We believe that protein language models are a
promising and largely unexplored field and discuss their foreseeable impact on
protein design.


  Chimeric antigen receptor (CAR) T-cell based immunotherapy has shown its
potential in treating blood cancers, and its application to solid tumors is
currently being extensively investigated. For glioma brain tumors, various CAR
T-cell targets include IL13Ra2, EGFRvIII, HER2, EphA2, GD2, B7-H3, and
chlorotoxin. In this work, we are interested in developing a mathematical model
of IL13Ra2 targeting CAR T-cells for treating glioma. We focus on extending the
work of Kuznetsov et al. (1994) by considering binding of multiple CAR T-cells
to a single glioma cell, and the dynamics of these multi-cellular conjugates.
Our model more accurately describes experimentally observed CAR T-cell killing
assay data than a model which does not consider cell binding. Moreover, we
derive conditions in the CAR T-cell expansion rate that determines treatment
success or failure. Finally, we show that our model captures distinct CAR
T-cell killing dynamics at low, medium, and high antigen receptor densities in
patient-derived brain tumor cells.


  During the first months, the Covid-19 pandemic has required most countries to
implement complex sequences of non-pharmaceutical interventions, with the aim
of controlling the transmission of the virus in the population. To be able to
take rapid decisions, a detailed understanding of the current situation is
necessary. Estimates of time-varying, instantaneous reproduction numbers
represent a way to quantify the viral transmission in real time. They are often
defined through a mathematical compartmental model of the epidemic, like a
stochastic SEIR model, whose parameters must be estimated from multiple time
series of epidemiological data. Because of very high dimensional parameter
spaces (partly due to the stochasticity in the spread models) and incomplete
and delayed data, inference is very challenging. We propose a state space
formalisation of the model and a sequential Monte Carlo approach which allow to
estimate a daily-varying reproduction number for the Covid-19 epidemic in
Norway with sufficient precision, on the basis of daily hospitalisation and
positive test incidences. The method is in regular use in Norway and is a
powerful instrument for epidemic monitoring and management.


  Quantification and classification of protein structures, such as knotted
proteins, often requires noise-free and complete data. Here we develop a
mathematical pipeline that systematically analyzes protein structures. We
showcase this geometric framework on proteins forming open-ended trefoil knots,
and we demonstrate that the mathematical tool, persistent homology, faithfully
represents their structural homology. This topological pipeline identifies
important geometric features of protein entanglement and clusters the space of
trefoil proteins according to their depth. Persistence landscapes quantify the
topological difference between a family of knotted and unknotted proteins in
the same structural homology class. This difference is localized and
interpreted geometrically with recent advancements in systematic computation of
homology generators. The topological and geometric quantification we find is
robust to noisy input data, which demonstrates the potential of this approach
in contexts where standard knot theoretic tools fail.


  Proteins can exhibit dynamic structural flexibility as they carry out their
functions, especially in binding regions that interact with other molecules.
For the key SARS-CoV-2 spike protein that facilitates COVID-19 infection,
studies have previously identified several such highly flexible regions with
therapeutic importance. However, protein structures available from the Protein
Data Bank are presented as static snapshots that may not adequately depict this
flexibility, and furthermore these cannot keep pace with new mutations and
variants. In this paper we present a sequential Monte Carlo method for broadly
sampling the 3-D conformational space of protein structure, according to the
Boltzmann distribution of a given energy function. Our approach is distinct
from previous sampling methods that focus on finding the lowest-energy
conformation for predicting a single stable structure. We exemplify our method
on the SARS-CoV-2 omicron variant as an application of timely interest. Our
results identify sequence positions 495-508 as a key region where omicron
mutations have the most impact on the space of possible conformations, which
coincides with the findings of other preliminary studies on the binding
properties of the omicron variant.


  In this paper, we discuss the conceptual underpinnings of Modern Coexistence
Theory (MCT), a quantitative framework for understanding ecological
coexistence. In order to use MCT to infer how species are coexisting, one must
relate a complex model (which simulates coexistence in the real world) to
simple models in which previously proposed explanations for coexistence have
been codified. This can be accomplished in three steps: 1) relating the
construct of coexistence to invasion growth rates, 2) mathematically
partitioning the invasion growth rates into coexistence mechanisms (i.e.,
classes of explanations for coexistence), and 3) relating coexistence
mechanisms to simple explanations for coexistence. Previous research has
primarily focused on step 2. Here, we discuss the other crucial steps and their
implications for inferring the mechanisms of coexistence in real communities.
  Our discussion of step 3 -- relating coexistence mechanisms to simple
explanations for coexistence -- serves a heuristic guide for hypothesizing
about the causes of coexistence in new models; but also addresses
misconceptions about coexistence mechanisms. For example, the storage effect
has little to do with bet-hedging or "storage" via a robust life-history stage;
relative nonlinearity is more likely to promote coexistence than originally
thought; and fitness-density covariance is an amalgam of a large number of
previously proposed explanations for coexistence (e.g., the
competition-colonization trade-off, heteromyopia, spatially-varying resource
supply ratios). Additionally, we review a number of topics in MCT, including
the role of "scaling factors"; whether coexistence mechanisms are
approximations; whether the magnitude or sign of invasion growth rates matters
more; whether Hutchinson solved the paradox of the plankton; the
scale-dependence of coexistence mechanisms; and much more.


  Ecologists have put forward many explanations for coexistence, but these are
only partial explanations; nature is complex, so it is reasonable to assume
that in any given ecological community, multiple mechanisms of coexistence are
operating at the same time. Here, we present a methodology for quantifying the
relative importance of different explanations for coexistence, based on an
extension of Modern Coexistence Theory. Current versions of Modern Coexistence
Theory only allow for the analysis of communities that are affected by spatial
or temporal environmental variation, but not both. We show how to analyze
communities with spatiotemporal fluctuations, how to parse the importance of
spatial variation and temporal variation, and how to measure everything with
either mathematical expressions or simulation experiments. Our extension of
Modern Coexistence Theory allows empiricists to use realistic models and more
data to better infer the mechanisms of coexistence in real communities.


  Machine learning potentials have emerged as a means to enhance the accuracy
of biomolecular simulations. However, their application is constrained by the
significant computational cost arising from the vast number of parameters
compared to traditional molecular mechanics. To tackle this issue, we introduce
an optimized implementation of the hybrid method (NNP/MM), which combines
neural network potentials (NNP) and molecular mechanics (MM). This approach
models a portion of the system, such as a small molecule, using NNP while
employing MM for the remaining system to boost efficiency. By conducting
molecular dynamics (MD) simulations on various protein-ligand complexes and
metadynamics (MTD) simulations on a ligand, we showcase the capabilities of our
implementation of NNP/MM. It has enabled us to increase the simulation speed by
5 times and achieve a combined sampling of one microsecond for each complex,
marking the longest simulations ever reported for this class of simulation.


  The voter process is a classic stochastic process that models the invasion of
a mutant trait $A$ (e.g., a new opinion, belief, legend, genetic mutation,
magnetic spin) in a population of agents (e.g., people, genes, particles) who
share a resident trait $B$, spread over the nodes of a graph. An agent may
adopt the trait of one of its neighbors at any time, while the invasion bias
$r\in(0,\infty)$ quantifies the stochastic preference towards ($r>1$) or
against ($r<1$) adopting $A$ over $B$. Success is measured in terms of the
fixation probability, i.e., the probability that eventually all agents have
adopted the mutant trait $A$. In this paper we study the problem of fixation
probability maximization under this model: given a budget $k$, find a set of
$k$ agents to initiate the invasion that maximizes the fixation probability. We
show that the problem is NP-hard for both $r>1$ and $r<1$, while the latter
case is also inapproximable within any multiplicative factor. On the positive
side, we show that when $r>1$, the optimization function is submodular and thus
can be greedily approximated within a factor $1-1/e$. An experimental
evaluation of some proposed heuristics corroborates our results.


  It is well established that people with diabetes are more likely to have
serious complications from COVID-19. Nearly 1 in 5 COVID-19 deaths in the
African region are linked to diabetes. World Health Organization (WHO) finds
that 18.3% of COVID-19 deaths in Africa are among people with diabetes. In this
paper, we have formulated and analysed a mathematical comorbidity model of
diabetes - COVID-19 of the deterministic type. The basic properties of the
model were explored. The basic reproductive number, equilibrium points and
stability of the equilibrium points were examined. Sensitivity analysis of the
model was carried on to determine the impact of the model parameters on the
basic reproduction number of the model. The model had a unique endemic
equilibrium point, which was stable for R_0>1. Time-dependent optimal controls
were incorporated into the model with the sole aim of determining the best
strategy for curtailing the spread of the disease. COVID-19 cases from March to
September 2020 in Ghana were used to validate the model. Results of the
numerical simulation suggest a greater number of individuals deceased when the
infected individual had an underlying condition of diabetes. More so COVID-19
is endemic in Ghana with the basic reproduction number found to be R_0=1.4722.
The numerical simulation of the optimal control model reveals the lockdown
control minimized the rate of decay of the susceptible individuals whereas the
vaccination led to a number of susceptible individuals becoming immune to
COVID-19 infections. In all the two preventive control measures were both
effective in curbing the spread of the COVID-19 disease as the number of
COVID-19 infections was greatly reduced. We conclude that more attention should
be paid to COVID-19 patients with an underlying condition of diabetes as the
probability of death in this population was significantly higher.


  Networks with absolute concentration robustness (ACR) have the property that
a translation of a coordinate hyperplane either contains all steady states
(static ACR) or attracts all trajectories (dynamic ACR). The implication for
the underlying biological system is robustness in the concentration of one of
the species independent of the initial conditions as well as independent of the
concentration of all other species. Identifying network conditions for dynamic
ACR is a challenging problem. We lay the groundwork in this paper by studying
small reaction networks, those with 2 reactions and 2 species. We give a
complete classification by ACR properties of these minimal reaction networks.
The dynamics is rich even within this simple setting. Insights obtained from
this work will help illuminate the properties of more complex networks with
dynamic ACR.


  The tetanus neurotoxin (TeNT) is one of the most toxic proteins known to man,
which prior to the use of the vaccine against the TeNT producing bacteria
Clostridium tetani, resulted in a 20 % mortality rate upon infection. The
clinical detrimental effects of tetanus have decreased immensely since the
introduction of global vaccination programs, which depend on sustainable
vaccine production. One of the major critical points in the manufacturing of
these vaccines is the stable and reproducible production of high levels of
toxin by the bacterial seed strains. In order to minimize time loss, the amount
of TeNT is often monitored during and at the end of the bacterial culturing.
The different methods that are currently available to assess the amount of TeNT
in the bacterial medium suffer from variability, lack of sensitivity, and/or
require specific antibodies. In accordance with the consistency approach and
the three Rs (3Rs), both aiming to reduce the use of animals for testing,
in-process monitoring of TeNT production could benefit from animal and
antibody-free analytical tools. In this paper, we describe the development and
validation of a new and reliable antibody free targeted LC-MS/MS method that is
able to identify and quantify the amount of TeNT present in the bacterial
medium during the different production time points up to the harvesting of the
TeNT just prior to further upstream purification and detoxification. The
quantitation method, validated according to ICH guidelines and by the
application of the total error approach, was utilized to assess the amount of
TeNT present in the cell culture medium of two TeNT production batches during
different steps in the vaccine production process prior to the generation of
the toxoid. The amount of TeNT generated under different physical stress
conditions applied during bacterial culture was also monitored.


  Research focus on optimal control problems brought on by fractional
differential equations has been extensively applied in practice. However,
because they are still open-ended and challenging, a number of problems with
fractional mathematical modeling and problems with optimal control require
additional study. Using fractional-order derivatives defined in the Atangana
Baleanu Caputo sense, we alter the integer-order model that has been proposed
in the literature. We prove the solution's existence, uniqueness, equilibrium
points, fundamental reproduction number, and local stability of the equilibrium
points. The operator's numerical approach was put into practice to obtain a
numerical simulation to back up the analytical conclusions. Fractional optimum
controls were incorporated into the model to identify the most efficient
intervention strategies for controlling the disease. Utilizing actual data from
Ghana for the months of March 2020 to March 2021, the model is validated. The
simulation's results show that the fractional operator significantly affected
each compartment and that the incidence rate of the population rose when v>0.6.
The examination of the most effective control technique discovered that social
exclusion and vaccination were both very effective methods for halting the
development of the illness.


  Straight-forward conformation generation models, which generate 3-D
structures directly from input molecular graphs, play an important role in
various molecular tasks with machine learning, such as 3D-QSAR and virtual
screening in drug design. However, existing loss functions in these models
either cost overmuch time or fail to guarantee the equivalence during
optimization, which means treating different items unfairly, resulting in poor
local geometry in generated conformation. So, we propose Equivalent Distance
Geometry Error (EDGE) to calculate the differential discrepancy between
conformations where the essential factors of three kinds in conformation
geometry (i.e. bond lengths, bond angles and dihedral angles) are equivalently
optimized with certain weights. And in the improved version of our method, the
optimization features minimizing linear transformations of atom-pair distances
within 3-hop. Extensive experiments show that, compared with existing loss
functions, EDGE performs effectively and efficiently in two tasks under the
same backbones.


  We study a dynamic infection spread model, inspired by the discrete time SIR
model, where infections are spread via non-isolated infected individuals. While
infection keeps spreading over time, a limited capacity testing is performed at
each time instance as well. In contrast to the classical, static, group testing
problem, the objective in our setup is not to find the minimum number of
required tests to identify the infection status of every individual in the
population, but to control the infection spread by detecting and isolating the
infections over time by using the given, limited number of tests. In order to
analyze the performance of the proposed algorithms, we focus on the mean-sense
analysis of the number of individuals that remain non-infected throughout the
process of controlling the infection. We propose two dynamic algorithms that
both use given limited number of tests to identify and isolate the infections
over time, while the infection spreads. While the first algorithm is a dynamic
randomized individual testing algorithm, in the second algorithm we employ the
group testing approach similar to the original work of Dorfman. By considering
weak versions of our algorithms, we obtain lower bounds for the performance of
our algorithms. Finally, we implement our algorithms and run simulations to
gather numerical results and compare our algorithms and theoretical
approximation results under different sets of system parameters.


  Returning universities to full on-campus operations while the COVID-19
pandemic is ongoing has been a controversial discussion in many countries. The
risk of large outbreaks in dense course settings is contrasted by the benefits
of in-person teaching. Transmission risk depends on a range of parameters, such
as vaccination coverage and efficacy, number of contacts and adoption of
non-pharmaceutical intervention measures (NPIs). Due to the generalised
academic freedom in Europe, many universities are asked to autonomously decide
on and implement intervention measures and regulate on-campus operations. In
the context of rapidly changing vaccination coverage and parameters of the
virus, universities often lack sufficient scientific insight to base these
decisions on. To address this problem, we analyse a calibrated, data-driven
agent-based simulation of transmission dynamics of 10755 students and 974
faculty members in a medium-sized European university. We use a co-location
network reconstructed from student enrollment data and calibrate transmission
risk based on outbreak size distributions in education institutions. We focus
on actionable interventions that are part of the already existing
decision-making process of universities to provide guidance for concrete policy
decisions. Here we show that, with the Omicron variant of the SARS-CoV-2 virus,
even a reduction to 25% occupancy and universal mask mandates are not enough to
prevent large outbreaks given the vaccination coverage of about 80% recently
reported for students in Austria. Our results show that controlling the spread
of the virus with available vaccines in combination with NPIs is not feasible
in the university setting if presence of students and faculty on campus is
required.


  Many multi-genic systemic diseases such as neurological disorders,
inflammatory diseases, and the majority of cancers do not have effective
treatments yet. Reinforcement learning powered systems pharmacology is a
potentially effective approach to design personalized therapies for untreatable
complex diseases. In this survey, state-of-the-art reinforcement learning
methods and their latest applications to drug design are reviewed. The
challenges on harnessing reinforcement learning for systems pharmacology and
personalized medicine are discussed. Potential solutions to overcome the
challenges are proposed. In spite of successful application of advanced
reinforcement learning techniques to target-based drug discovery, new
reinforcement learning strategies are needed to address systems
pharmacology-oriented personalized de novo drug design.


  Recently, symbolic computation and computer algebra systems have been
successfully applied in systems biology, especially in chemical reaction
network theory. One advantage of symbolic computation is its potential for
qualitative answers to biological questions. Qualitative methods analyze
dynamical input systems as formal objects, in contrast to investigating only
part of the state space, as is the case with numerical simulation. However,
symbolic computation tools and libraries have a different set of requirements
for their input data than their numerical counterparts. A common format used in
mathematical modeling of biological processes is SBML. We illustrate that the
use of SBML data in symbolic computation requires significant pre-processing,
incorporating external biological and mathematical expertise. ODEbase provides
high quality symbolic computation input data derived from established existing
biomodels, covering in particular the BioModels database.


  The inclusion of cognitive processes, such as perception, learning and
memory, are inevitable in mechanistic animal movement modelling. Cognition is
the unique feature that distinguishes animal movement from mere particle
movement in chemistry or physics. Hence, it is essential to incorporate such
knowledge-based processes into animal movement models. Here, we summarize
popular deterministic mathematical models derived from first principles that
begin to incorporate such influences on movement behaviour mechanisms. Most
generally, these models take the form of nonlocal reaction-diffusion-advection
equations, where the nonlocality may appear in the spatial domain, the temporal
domain, or both. Mathematical rules of thumb are provided to judge the model
rationality, to aid in model development or interpretation, and to streamline
an understanding of the range of difficulty in possible model conceptions. To
emphasize the importance of biological conclusions drawn from these models, we
briefly present available mathematical techniques and introduce some existing
"measures of success" to compare and contrast the possible predictions and
outcomes. Throughout the review, we propose numerous open problems relevant to
this relatively new area, ranging from precise technical mathematical
challenges to broader conceptual challenges at the cross-section between
mathematics and ecology. This review paper is expected to act as a synthesis of
existing efforts while also pushing the boundaries of current modelling
perspectives to better understand the influence of cognitive movement
mechanisms on movement behaviours and space use outcomes.


  In this paper, we use SEIR equations to make predictions for the number of
mortality due to COVID-19 in \.Istanbul. Using excess mortality method, we find
the number of mortality for the previous three waves in 2020 and 2021. We show
that the predictions of our model is consistent with number of moralities for
each wave. Furthermore, we predict the number of mortality for the second wave
of 2021. We also extend our analysis for Germany, Italy and Turkey to compare
the basic reproduction number $R_0$ for Istanbul. Finally, we calculate the
number of infected people in Istanbul for herd immunity.


  We developed an integrated recurrent neural network and nonlinear regression
spatio-temporal model for vector-borne disease evolution. We take into account
climate data and seasonality as external factors that correlate with disease
transmitting insects (e.g. flies), also spill-over infections from neighboring
regions surrounding a region of interest. The climate data is encoded to the
model through a quadratic embedding scheme motivated by recommendation systems.
The neighboring regions' influence is modeled by a long short-term memory
neural network. The integrated model is trained by stochastic gradient descent
and tested on leish-maniasis data in Sri Lanka from 2013-2018 where infection
outbreaks occurred. Our model outperformed ARIMA models across a number of
regions with high infections, and an associated ablation study renders support
to our modeling hypothesis and ideas.


  We introduce a derangement model of ligand-receptor binding that allows us to
quantitatively frame the question "How can ligands seek out and bind to their
optimal receptor sites in a sea of other competing ligands and suboptimal
receptor sites?" To answer the question, we first derive a formula to count the
number of partial generalized derangements in a list, thus extending the
derangement result of Gillis and Even. We then compute the general partition
function for the ligand-receptor system and derive the equilibrium expressions
for the average number of bound ligands and the average number of optimally
bound ligands. A visual model of squares assembling onto a grid allows us to
easily identify fully optimal bound states. Equilibrium simulations of the
system reveal its extremes to be one of two types, qualitatively distinguished
by whether optimal ligand-receptor binding is the dominant form of binding at
all temperatures and quantitatively distinguished by the relative values of two
critical temperatures. One of those system types (termed "search-limited," as
it was in previous work) does not exhibit kinetic traps and we thus infer that
biomolecular systems where optimal ligand-receptor binding is functionally
important are likely to be search-limited.


  In this paper we show that the dynamics of a class of kinetic compartmental
models with bounded capacities, monotone reaction rates and a strongly
connected interconnection structure is persistent. The result is based on the
chemical reaction network (CRN) and the corresponding Petri net representation
of the system. For the persistence analysis, it is shown that all siphons in
the Petri net of the studied model class can be characterized efficiently.
Additionally, the existence and stability of equilibria are also analyzed
building on the persistence and the theory of general compartmental systems.
The obtained results can be applied in the analysis of general kinetic models
based on the simple exclusion principle.


  The AlphaFold computer program predicted protein structures for the whole
human genome, which has been considered as a remarkable breakthrough both in
artificial intelligence (AI) application and structural biology. Despite the
varying confidence level, these predicted structures still could significantly
contribute to structure-based drug design of novel targets, especially the ones
with no or limited structural information. In this work, we successfully
applied AlphaFold in our end-to-end AI-powered drug discovery engines
constituted of a biocomputational platform PandaOmics and a generative
chemistry platform Chemistry42, to identify a first-in-class hit molecule of a
novel target without an experimental structure starting from target selection
towards hit identification in a cost- and time-efficient manner. PandaOmics
provided the targets of interest and Chemistry42 generated the molecules based
on the AlphaFold predicted structure, and the selected molecules were
synthesized and tested in biological assays. Through this approach, we
identified a small molecule hit compound for CDK20 with a Kd value of 8.9 +/-
1.6 uM (n = 4) within 30 days from target selection and after only synthesizing
7 compounds. Based on the available data, the second round of AI-powered
compound generation was conducted and through which, a more potent hit
molecule, ISM042-2 048, was discovered with a Kd value of 210.0 +/- 42.4 nM (n
= 2), within 30 days and after synthesizing 6 compounds from the discovery of
the first hit ISM042-2-001. To the best of our knowledge, this is the first
reported small molecule targeting CDK20 and more importantly, this work is the
first demonstration of AlphaFold application in the hit identification process
in early drug discovery.


  While understanding of periodic recurrent waves of Covid-19 epidemics would
aid to combat the pandemics, quantitative analysis of data over a two years
period from the outbreak, is lacking. The complexity of Covid-19 recurrent
waves is related with the concurrent role of i) the containment measures
enforced to mitigate the epidemics spreading ii) the rate of viral gene
mutations, and iii) the variable immune response of the host implemented by
vaccination. This work focuses on the effect of massive vaccination and gene
variants on the recurrent waves in a representative case of countries enforcing
mitigation and vaccination strategy. The spreading rate is measured by the
ratio between the reproductive number Rt(t) and the doubling time Td(t) called
RIC-index and the daily fatalities number. The dynamics of the Covid-19
epidemics has been studied by wavelet analysis and represented by a non-linear
helicoid vortex in a 3D space where both RIC-index and fatalities change with
time. The onset of periodic recurrent waves has been identified by the
transition from convergent to divergent trajectories on the helicoid vortex. We
report a main period of recurrent waves of 120 days and the elongation of this
period after the vaccination campaign


  There are many scales at which to quantify stability in spatial and
ecological networks. Local-scale analyses focus on specific nodes of the
spatial network, while regional-scale analyses consider the whole network.
Similarly, species- and community-level analyses either account for single
species or for the whole community. Furthermore, stability itself can be
defined in multiple ways, including resistance (the inverse of the relative
displacement caused by a perturbation), initial resilience (the rate of return
after a perturbation), and invariability (the inverse of the relative amplitude
of the population fluctuations). Here, we analyze the scale-dependence of these
stability properties. More specifically, we ask how spatial scale (local vs
regional) and ecological scale (species vs community) influence these stability
properties. We find that regional initial resilience is the weighted arithmetic
mean of the local initial resiliences. The regional resistance is the harmonic
mean of local resistances, which makes regional resistance particularly
vulnerable to nodes with low stability, unlike regional initial resilience.
Analogous results hold for the relationship between community- and
species-level initial resilience and resistance. Both resistance and initial
resilience are ``scale-free'' properties: regional and community values are
simply the biomass-weighted means of the local and species values,
respectively. Thus, one can easily estimate both stability metrics of whole
networks from partial sampling. In contrast, invariability generally is greater
at the regional and community-level than at the local and species-level,
respectively. Hence, estimating the invariability of spatial or ecological
networks from measurements at the local or species level is more complicated,
requiring an unbiased estimate of the network (i.e. region or community) size.


  We present a stochastic epidemic model to study the effect of various
preventive measures, such as uniform reduction of contacts and transmission,
vaccination, isolation, screening and contact tracing, on a disease outbreak in
a homogeneously mixing community. The model is based on an infectivity process,
which we define through stochastic contact and infectiousness processes, so
that each individual has an independent infectivity profile. In particular, we
monitor variations of the reproduction number and of the distribution of
generation times. We show that some interventions, i.e. uniform reduction and
vaccination, affect the former while leaving the latter unchanged, whereas
other interventions, i.e. isolation, screening and contact tracing, affect both
quantities. We provide a theoretical analysis of the variation of these
quantities, and we show that, in practice, the variation of the generation time
distribution can be significant and that it can cause biases in the estimation
of reproduction numbers. The framework, because of its general nature, captures
the properties of many infectious diseases, but particular emphasis is on
COVID-19, for which numerical results are provided.


  Infectious disease transmission dynamics are particularly sensitive to social
contact patterns, and the precautions people take to limit disease
transmission. It depends on the age distribution of the community. Thus,
knowing the age$-$specific prevalence and incidence of infectious diseases is
critical for predicting future disease burden and the efficacy of interventions
like immunization. This study uses an SEIR age$-$structured multi-group
epidemic model to understand how social contact affects disease control. We
created location$-$specific social contact matrices in the community to see how
social mixing has affected illness spread. We estimated the basic reproduction
number $(R_0)$ of the system and plotted its global behavior in terms of $R_0$.
Optimal control for the problem has also been established quantitatively. The
proposed model's transmission rate for India from September 1, 2020, to
December 31, 2020, has also been estimated. We replicated the lifting of
non-pharmaceutical therapies by allowing participants to return to work in
phases and studied the impact of this. Our findings imply that identifying
symptomatic sick people aged $20-49$ can help lower the number of infected
people when schools are closed. When some schools are partially open, awareness
of symptomatic infected persons helps reduce cases. The simulation results also
suggest that limiting contact at school and other meeting areas could
significantly lower the number of instances. Using the least square approach,
it was discovered that the time$-$dependent transmission rate is more realistic
than the constant spread rate for COVID$-19$ in India. To reduce the COVID$-19$
in burden, we found that gradually loosening control methods could flatten and
lower other peaks. Our findings may help health policymakers decide on timely
age$-$based immunization distribution strategies and hence control the disease.


  We will study the relation between two well-known theories, genetic evolution
and random matrix theory in the context of many-body systems and chaos theory.
We show that the time evolution of certain chaotic quantum mechanical models is
similar to the evolution of a living cell. It is also suggested that natural
selection can be described by a random matrix theory with statistical
distribution in which the genetic evolution acts as a Gross-Witten-Wadia phase
transition.


  A new strategy was developed to investigate the effect of volatile fatty
acids (VFAs) on the efficiency of biogas production with a focus on improving
bio-H$_2$. The inoculum used, anaerobic granular sludge obtained from a UASB
reactor treating poultry slaughterhouse wastewater, was pretreated with five
different pretreatments. The relationship between VFAs and biogas compounds was
studied as time-dependent components. In time-dependent processes with small
sample size data, regression models may not be good enough at estimating
responses. Therefore, a deep learning neural network (DNN) model was developed
to estimate the biogas compounds based on the VFAs. The accuracy of this model
to predict the biogas compounds was higher than that of multivariate regression
models. Further, it could predict the effect of time changes on biogas
compounds. Analysis showed that all the pretreatments were able to increase the
ratio of butyric acid / acetic acid successfully, decrease propionic acid
drastically, and increase the efficiency of bio-H$_2$ production. As
discovered, butyric acid had the greatest effect on bio-H$_2$, and propionic
acid had the greatest effect on CH$_4$ production. The best amounts of the VFAs
were determined using an optimization method, integrated DNN and desirability
analysis, dynamically retrained based on digestion time. Accordingly, optimal
ranges of acetic, propionic, and butyric acids were 823.2 - 1534.3, 36.3 -
47.4, and 1522 - 1822 mg/L, respectively, determined for digestion time of
25.23 - 123.63 h. These values resulted in the production of bio-H$_2$, N$_2$,
CO$_2$, and CH$_4$ in ranges of 6.4 - 26.2, 12.2 - 43.2, 5 - 25.3, and 0 - 1.4
mmol/L, respectively. The optimum ranges of VFAs are relatively wide ranges and
practically can be used in biogas plants.


  In many social mammals, early life social adversity and social integration
largely predict individual health, lifespan and reproductive success. Efforts
in identifying the physiological mechanisms mediating the relationship between
the social environment and individual fitness have so far concentrated on
socially-induced stress, mediated by alterations in neuroendocrine signaling
and immune function. Here, we propose a much-needed alternative mechanism
relying on microbially-mediated effects: social relationships with
conspecifics, both in early life and adulthood, might strongly contribute both
to the transmission of beneficial microbes and to diversifying host
microbiomes. In turn, more valuable and diverse microbiomes would promote
pathogen resistance and optimal health and thus translate into positive fitness
outcomes. This mechanism relies on two emerging findings from empirical
studies, namely that microbiomes (i) are largely socially transmitted via
vertical and horizontal routes, and (ii) play a pervasive role in host
development, physiology, metabolism, and susceptibility to pathogens. We
suggest that the social transmission of microbiomes has the potential to
explain the sociality-fitness nexus, to a similar - or even higher - extent
than chronic social stress, in ways that have yet to be studied empirically in
social mammals.


  Performing full-resolution atomistic simulations of nucleic acid folding has
remained a challenge for biomolecular modeling. Understanding how nucleic acids
fold and how they transition between different folded structures as they unfold
and refold has important implications for biology. This paper reports a
theoretical model and computer simulation of the ab initio folding of DNA
inverted repeat sequences. The formulation is based on an all-atom
conformational model of the sugar-phosphate backbone via chain closure, and it
incorporates three major molecular-level driving forces - base stacking,
counterion-induced backbone self-interactions and base pairing - via separate
analytical theories designed to capture and reproduce the effects of the
solvent without requiring explicit water and ions in the simulation. To
accelerate computational throughput, a mixed numerical/analytical algorithm for
the calculation of the backbone conformational volume is incorporated into the
Monte Carlo simulation, and special stochastic sampling techniques were
employed to achieve the computational efficiency needed to fold nucleic acids
from scratch. This paper describes implementation details, benchmark results
and the advantages and technical challenges with this approach.


  A theoretical study of cell evolution is presented here. By using a toolbox
containing an intracellular catalytic reaction network model and a
mutation-selection process, four distinct phases of self-organization were
unveiled. First, the nutrients prevail as the central substrate of the chemical
reactions. Second, the cell becomes a small-world. Third, a highly connected
core component emerges, concurrently with the nutrient carriers becoming the
central product of reactions. Finally, the cell reaches a steady configuration
where the concentrations of the core chemical species are described by Zipf's
law.


  Deep learning has contributed to major advances in the prediction of protein
structure from sequence, a fundamental problem in structural bioinformatics.
With predictions now approaching the accuracy of crystallographic resolution in
some cases, and with accelerators like GPUs and TPUs making inference using
large models rapid, fast genome-level structure prediction becomes an obvious
aim. Leadership-class computing resources can be used to perform genome-scale
protein structure prediction using state-of-the-art deep learning models,
providing a wealth of new data for systems biology applications. Here we
describe our efforts to efficiently deploy the AlphaFold2 program, for
full-proteome structure prediction, at scale on the Oak Ridge Leadership
Computing Facility's resources, including the Summit supercomputer. We
performed inference to produce the predicted structures for 35,634 protein
sequences, corresponding to three prokaryotic proteomes and one plant proteome,
using under 4,000 total Summit node hours, equivalent to using the majority of
the supercomputer for one hour. We also designed an optimized structure
refinement that reduced the time for the relaxation stage of the AlphaFold
pipeline by over 10X for longer sequences. We demonstrate the types of analyses
that can be performed on proteome-scale collections of sequences, including a
search for novel quaternary structures and implications for functional
annotation.


  A hallmark of living systems is the ability to employ a common set of
versatile building blocks that can self-organize into a multitude of different
structures, in a way that can be controlled with minimal cost. This capability
can only be afforded in non-equilibrium conditions, as evident from the
energy-consuming nature of the plethora of such dynamical processes. In the
last three decades, synthetic self-assembly has experienced a significant boost
with the development of tools to design specific interactions at different
scales, from nucleic acids and peptides to proteins and colloids. To achieve
automated dynamical control of such self-assembled structures and transitions
between them, we need to identify the relevant fundamental aspects of
non-equilibrium dynamics that can enable such processes. Here, we identify
programmable non-reciprocal interactions as a potential paradigm using which
such functionalities can be achieved. In particular, we propose a model that
enables a system to learn and retrieve predetermined desired structures and
transition between them, thus behaving as a shape-shifter. The learning rule is
composed of reciprocal interactions that lead to the equilibrium assembly of
the structures, and non-reciprocal interactions that give rise to
non-equilibrium dynamical transitions between the structures.


  Decomposing a network flow into weighted paths has numerous applications.
Some applications require any decomposition that is optimal w.r.t. some
property such as number of paths, robustness, or length. Many bioinformatic
applications require a specific decomposition where the paths correspond to
some underlying data that generated the flow. For real inputs, no optimization
criteria guarantees to uniquely identify the correct decomposition. Therefore,
we propose to report safe paths, i.e., subpaths of at least one path in every
flow decomposition.
  Ma, Zheng, and Kingsford [WABI 2020] addressed the existence of multiple
optimal solutions in a probabilistic framework, i.e., non-identifiability.
Later [RECOMB 2021], they gave a quadratic-time algorithm based on a global
criterion for solving a problem called AND-Quant, which generalizes the problem
of reporting whether a given path is safe.
  We give the first local characterization of safe paths for flow
decompositions in directed acyclic graphs (DAGs), leading to a practical
algorithm for finding the complete set of safe paths. We evaluated our
algorithms against the trivial safe algorithms (unitigs, extended unitigs) and
the popularly used heuristic (greedy-width) for flow decomposition on RNA
transcripts datasets. Despite maintaining perfect precision our algorithm
reports significantly higher coverage ($\approx 50\%$ more) than trivial safe
algorithms. The greedy-width algorithm though reporting a better coverage, has
significantly lower precision on complex graphs. Overall, our algorithm
outperforms (by $\approx 20\%$) greedy-width on a unified metric (F-Score) when
the dataset has significant number of complex graphs. Moreover, it has superior
time ($3-5\times$) and space efficiency ($1.2-2.2\times$), resulting in a
better and more practical approach for bioinformatics applications of flow
decomposition.


  In recent years, the monarch butterfly's iconic migration patterns have come
under threat from a number of factors, from climate change to pesticide use. To
track trends in their populations, scientists as well as citizen scientists
must identify individuals accurately. This is uniquely key for the study of
monarch butterflies because there exist other species of butterfly, such as
viceroy butterflies, that are "look-alikes" (coined by the Convention on
International Trade in Endangered Species of Wild Fauna and Flora), having
similar phenotypes. To tackle this problem and to aid in more efficient
identification, we present MonarchNet, the first comprehensive dataset
consisting of butterfly imagery for monarchs and five look-alike species. We
train a baseline deep-learning classification model to serve as a tool for
differentiating monarch butterflies and its various look-alikes. We seek to
contribute to the study of biodiversity and butterfly ecology by providing a
novel method for computational classification of these particular butterfly
species. The ultimate aim is to help scientists track monarch butterfly
population and migration trends in the most precise and efficient manner
possible.


  During the COVID pandemic, periods of exponential growth of the disease have
been mitigated by containment measures that in different occasions have
resulted in a power-law growth of the number of cases. The first observation of
such behaviour has been obtained from 2020 late spring data coming from China
by Ziff and Ziff in Ref. [1]. After this important observation the power-law
scaling (albeit with different exponents) has also been observed in other
countries during periods of containment of the spread. Early interpretations of
these results suggest that this phenomenon might be due to spatial effects of
the spread. Here we show that temporal modulations of infectivity of
individuals due to containment measures can also cause power-law growth of the
number of cases over time. To this end we propose a stochastic well-mixed
Susceptible-Infected-Removed (SIR) model of epidemic spreading in presence of
containment measures resulting in time-dependent infectivity and we explore the
statistical properties of the resulting branching process at criticality. We
show that at criticality it is possible to observe power-law growth of the
number of cases with exponents ranging between one and two. Our asymptotic
analytical results are confirmed by extensive Monte Carlo simulations. Although
these results do not exclude that spatial effects might be important in
modulating the power-law growth of the number of cases at criticality, this
work shows that even well-mixed populations may already feature non-trivial
power-law exponents at criticality.


  Phase separation has emerged as an essential concept for the spatial
organization inside biological cells. However, despite the clear relevance to
virtually all physiological functions, we understand surprisingly little about
what phases form in a system of many interacting components, like in cells.
Here, we introduce a new numerical method based on physical relaxation dynamics
to study the coexisting phases in such systems. We use our approach to optimize
interactions between components, similar to how evolution might have optimized
the interactions of proteins. These evolved interactions robustly lead to a
defined number of phases, despite substantial uncertainties in the initial
composition, while random or designed interactions perform much worse.
Moreover, the optimized interactions are robust to perturbations and they allow
fast adaption to new target phase counts. We thus show that genetically encoded
interactions of proteins provide versatile control of phase behavior. The
phases forming in our system are also a concrete example of a robust emergent
property that does not rely on fine-tuning the parameters of individual
constituents.


  Minimum flow decomposition (MFD) (the problem of finding a minimum set of
paths that perfectly decomposes a flow) is a classical problem in Computer
Science, and variants of it are powerful models in multiassembly problems in
Bioinformatics (e.g. RNA assembly). However, because this problem and its
variants are NP-hard, practical multiassembly tools either use heuristics or
solve simpler, polynomial-time solvable versions of the problem, which may
yield solutions that are not mini-mal or do not perfectly decompose the flow.
Many RNA assemblers also use integer linear programming(ILP) formulations of
such practical variants, having the major limitation they need to encode all
the potentially exponentially many solution paths. Moreover, the only exact
solver for MFD does not scale to large instances and cannot be efficiently
generalized to practical MFD variants. In this work, we provide the first
practical ILP formulation for MFD (and thus the first fast and exact solver for
MFD), based on encoding all of the exponentially many solution paths using only
a quadratic number of variables. On both simulated and real flow graphs, our
approach solves any instance in under 13 seconds. We also show that our ILP
formulation can be easily and efficiently adapted for many practical variants,
such as incorporating longer or paired-end reads or minimizing flow errors. We
hope that our results can remove the current tradeoff between the complexity of
a multi assembly model and its tractability and can lie at the core of future
practical RNA assembly tools.


  The Sackin and Colless indices are two widely-used metrics for measuring the
balance of trees and for testing evolutionary models in phylogenetics. This
short paper contributes two results about the Sackin and Colless indices of
trees. One result is the asymptotic analysis of the expected Sackin and Colless
indices of a tree shape (which are full binary rooted unlabelled trees) under
the uniform model where tree shapes are sampled with equal probability. Another
is a short elementary proof of the closed formula for the expected Sackin index
of phylogenetic trees (which are full binary rooted trees with leaves being
labelled with taxa) under the uniform model.


  Self-supervised protein language models have proved their effectiveness in
learning the proteins representations. With the increasing computational power,
current protein language models pre-trained with millions of diverse sequences
can advance the parameter scale from million-level to billion-level and achieve
remarkable improvement. However, those prevailing approaches rarely consider
incorporating knowledge graphs (KGs), which can provide rich structured
knowledge facts for better protein representations. We argue that informative
biology knowledge in KGs can enhance protein representation with external
knowledge. In this work, we propose OntoProtein, the first general framework
that makes use of structure in GO (Gene Ontology) into protein pre-training
models. We construct a novel large-scale knowledge graph that consists of GO
and its related proteins, and gene annotation texts or protein sequences
describe all nodes in the graph. We propose novel contrastive learning with
knowledge-aware negative sampling to jointly optimize the knowledge graph and
protein embedding during pre-training. Experimental results show that
OntoProtein can surpass state-of-the-art methods with pre-trained protein
language models in TAPE benchmark and yield better performance compared with
baselines in protein-protein interaction and protein function prediction. Code
and datasets are available in https://github.com/zjunlp/OntoProtein.


  Focal cortical dysplasias (FCDs) are the most common cause of treatment
resistant epilepsy affecting the pediatric population. Most individuals with
FCD have seizure onset during the first five years of life and the majority
will have seizures by the age of sixteen. Many cases of FCD are postulated to
be the result of abnormal brain development in utero by germline or somatic
gene mutations regulating neuronal growth and migration during corticogenesis.
Other cases of FCD are thought to be related to infections during brain
development, or even other causes still unable to be fully determined. Typical
anti-seizure medications are oftentimes ineffective in FCD as well as surgery
is unable to be successfully performed due to the involvement of eloquent areas
of the brain or insufficient resection of the epileptogenic focus, posing a
challenge for physicians. The genetic nature of FCD provides an avenue for drug
development with several genetic and molecular targets undergoing study over
the last two decades.


  Post-translational modifications (PTMs) in proteins occur after the process
of translation. PTMs account for many cellular processes such as
deoxyribonucleic acid (DNA) repair, cell signaling and cell death. One of the
recent PTMs is succinylation. Succinylation modifies lysine residue from $-1$
to $+1$. Locating succinylation sites using experimental methods, such as mass
spectrometry is very laborious. Hence, computational methods are favored using
machine learning techniques. This paper proposes a deep learning architecture
to predict succinylation sites. The performance of the proposed architecture is
compared to the state-of-the-art deep learning architecture and other
traditional machine learning techniques for succinylation. It is shown from the
performance metrics that the proposed architecture provides a good trade-off
between speed of computation and classification accuracy.


  In mathematical modeling, several different functional forms can often be
used to fit a data set equally well, especially if the data is sparse. In such
cases, these mathematically different but similar looking functional forms are
typically considered interchangeable. Recent work, however, shows that similar
functional responses may nonetheless result in significantly different
bifurcation points for the Rosenzweig-MacArthur predator-prey system. Since the
bifurcation behaviours include destabilising oscillations, predicting the
occurrence of such behaviours is clearly important. Ecologically, different
bifurcation behaviours mean that different predictions may be obtained from the
models. These predictions can range from stable coexistence to the extinction
of both species, so obtaining more accurate predictions is also clearly
important for conservationists. Mathematically, this difference in bifurcation
structure given similar functional responses is called structural sensitivity.
We extend the existing work to find that the Leslie-Gower-May predator-prey
system is also structurally sensitive to the functional response. Using the
Rosenzweig-MacArthur and Leslie-Gower-May models, we then aim to determine if
there is some way to obtain a functional description of data such that our
model is not structurally sensitive. We first add stochasticity to the
functional responses and find that better similarity of the resulting
bifurcation structures is achieved. Then, we analyze the functional responses
using two different methods to determine which part of each function
contributes most to the observed bifurcation behaviour. We find that prey
densities around the coexistence steady state are most important in defining
the functional response. Lastly, we propose a procedure for ecologists and
mathematical modelers to increase the accuracy of model predictions in
predator-prey systems.


  Physical mechanisms of phase separation in living systems can play key
physiological roles and have recently been the focus of intensive studies. The
strongly heterogeneous and disordered nature of such phenomena in the
biological domain poses difficult modeling challenges that require going beyond
mean field approaches based on postulating a free energy landscape. The
alternative pathway we take in this work is to tackle the full statistical
mechanics problem of calculating the partition function in these systems,
starting from microscopic interactions, by means of cavity methods. We
illustrate the procedure first on the simple binary case, and we then apply it
successfully to ternary systems, in which the naive mean field approximations
are proved inadequate. We then demonstrate the agreement with lattice model
simulations, to finally contrast our theory also with experiments of coacervate
formation by associative de-mixing of nucleotides and poly-lysine in aqueous
solution. In this way, different types of evidence are provided to support
cavity methods as ideal tools for quantitative modeling of biomolecular
condensation, giving an optimal balance between the accurate consideration of
spatial aspects of the microscopic dynamics and the fast computational results
rooted in their analytical tractability.


  Microscale oxygenation plays a prominent role in tumour progression.
Spatiotemporal variability of oxygen distribution in the tumour
microenvironment contributes to cellular heterogeneity and to the emergence of
normoxic and hypoxic populations. Local levels of oxygen strongly affect the
response of tumours to the administration of different therapeutic modalities
and, more generally, to the phenomenon of resistance to treatments. Several
interventions have been proposed to improve tumour oxygenation, being the
elevation of the local temperature (hyperthermia) an important one. While other
factors such as the metabolic activity have to be considered, the proficiency
of the tumour vascular system is a key factor both for the tissue oxygenation
and for its temperature maps. Consequently, the interplay of these factors has
attracted considerable attention from the mathematical modelling perspective.
Here we put forward a transport-based system of partial differential equations
aimed at describing the dynamics of healthy and tumour cell subpopulations at
the microscale in a region placed between two blood vessels. By using this
model with diverse flow conditions, we analyse the oxygen and temperature
profiles that arise in different scenarios of vascular status, both during free
progression and under thermal therapy. We find that low oxygen levels are
associated to elevations of temperature in locations preferentially populated
by hypoxic cells, and hyperthermia-induced cell death, being strongly dependent
on blood flow, would only appear under highly disrupted conditions of the local
vasculature. This results in a noticeable effect of heat on hypoxic cells.
Additionally, when pronounced cell death occurs, it is followed by a
significant increase in the oxygen levels.


  Studies have shown that some people with underlying conditions such as
cancer, heart failure, diabetes and hypertension are more likely to get
COVID-19 and have worse outcomes. In this paper, a fractional-order derivative
is proposed to study the transmission dynamics of COVID-19 taking into
consideration population having an underlying condition. The fractional
derivative is defined in the Atangana Beleanu and Caputo (ABC) sense. For the
proposed model, we find the basic reproductive number, the equilibrium points
and determine the stability of these equilibrium points. The existence and the
uniqueness of the solution are established along with Hyers Ulam Stability. The
numerical scheme for the operator was carried out to obtain a numerical
simulation to support the analytical results. COVID-19 cases from March to June
2020 of Ghana were used to validate the model. The numerical simulation
revealed a decline in infections as the fractional operator was increased from
0.6 within the 120 days. Time-dependent optimal control was incorporated into
the model. The numerical simulation of the optimal control revealed,
vaccination reduces the number of individuals susceptible to the COVID-19,
exposed to the COVID-19 and Covid-19 patients with and without an underlying
health condition.


  Phase separation of several different overall neutral polyampholyte species
(with zero net charge) is studied in solution with two oppositely charged ion
species that can form ion-pairs through an association reaction. A field theory
description of the system, that treats polyampholyte charge sequence dependent
electrostatic interactions as well as excluded volume effects, is hereby given.
Interestingly, analysis of the model using random phase approximation and field
theoretic simulation consistently show evidence of a re-entrant polyampholyte
phase separation at high ion concentrations when there is an overall decrease
of volume upon ion-association. As an illustration of the ramifications of our
theoretical framework, several polyampholyte concentration vs ion concentration
phase diagrams under constant temperature conditions are presented to elucidate
the dependence of phase separation behavior on polyampholyte sequence charge
pattern as well as ion-pair dissociation constant, volumetric effects on ion
association, solvent quality, and temperature.


  Nuclear Magnetic Resonance (NMR) spectroscopy is one of the major techniques
in structural biology with over 11,800 protein structures deposited in the
Protein Data Bank. NMR can elucidate structures and dynamics of small and
medium size proteins in solution, living cells, and solids, but has been
limited by the tedious data analysis process. It typically requires weeks or
months of manual work of a trained expert to turn NMR measurements into a
protein structure. Automation of this process is an open problem, formulated in
the field over 30 years ago. Here, we present a solution to this challenge that
enables the completely automated analysis of protein NMR data within hours
after completing the measurements. Using only NMR spectra and the protein
sequence as input, our machine learning-based method, ARTINA, delivers signal
positions, resonance assignments, and structures strictly without any human
intervention. Tested on a 100-protein benchmark comprising 1329
multidimensional NMR spectra, ARTINA demonstrated its ability to solve
structures with 1.44 {\AA} median RMSD to the PDB reference and to identify
91.36% correct NMR resonance assignments. ARTINA can be used by non-experts,
reducing the effort for a protein assignment or structure determination by NMR
essentially to the preparation of the sample and the spectra measurements.


  Significant efforts have been done in last two decades to develop nanoscale
spectroscopy techniques owning to their great potential for single-molecule
structural detection and in addition, to resolve open questions in
heterogeneous biological systems, such as protein-DNA complexes. Applying
AFM-IR technique has become a powerful leverage for obtaining simultaneous
absorption spectra with a nanoscale spatial resolution for studied proteins,
however the IR-AFM investigation of DNA molecules on surface, as a benchmark
for a nucleoprotein complexes nanocharacterization, has remained elusive.
Herein, we demonstrate methodological approach for acquisition of IR-AFM
mapping modalities with corresponding absorption spectra based on two different
DNA deposition protocols on spermidine and Ni2+ pretreated mica surface. The
nanoscale IR absorbance of distinctly formed DNA morphologies on mica are
demonstrated through series of IR-AFM absorption maps with corresponding IR
spectrum. Our results thus demonstrate the sensitivity of IR-AFM
nanospectroscopy for a nucleic acid research with an open potential to be
employed in further investigation of nucleoprotein complexes.


  With the increasing availability and size of multi-omics datasets,
investigating the casual relationships between molecular phenotypes has become
an important aspect of exploring underlying biology and genetics. This paper
aims to introduce and review the available methods for building large-scale
causal molecular networks that have been developed in the past decade. Existing
methods have their own strengths and limitations so there is no one best
approach, and it is instead down to the discretion of the researcher. This
review also aims to discuss some of the current limitations to biological
interpretation of these networks, and important factors to consider for future
studies on molecular networks.


  When considering large sets of molecules, it is helpful to place them in the
context of a "chemical space" - a multidimensional space defined by a set of
descriptors that can be used to visualize and analyze compound grouping as well
as identify regions that might be void of valid structures. The chemical space
of all possible molecules in a given biological or environmental sample can be
vast and largely unexplored, mainly due to current limitations in processing of
'big data' by brute force methods (e.g., enumeration of all possible compounds
in a space). Recent advances in artificial intelligence (AI) have led to
multiple new cheminformatics tools that incorporate AI techniques to
characterize and learn the structure and properties of molecules in order to
generate plausible compounds, thereby contributing to more accessible and
explorable regions of chemical space without the need for brute force methods.
We have used one such tool, a deep-learning software called DarkChem, which
learns a representation of the molecular structure of compounds by compressing
them into a latent space. With DarkChem's design, distance in this latent space
is often associated with compound similarity, making sparse regions interesting
targets for compound generation due to the possibility of generating novel
compounds. In this study, we used 1 million small molecules (less than 1000 Da)
to create a representative chemical space (defined by calculated molecular
properties) of all small molecules. We identified regions with few or no
compounds and investigated their location in DarkChem's latent space. From
these spaces, we generated 694,645 valid molecules, all of which represent
molecules not found in any chemical database to date. These molecules filled
50.8% of the probed empty spaces in molecular property space. Generated
molecules are provided in the supporting information.


  Consider a branching Markov process with values in some general type space.
Conditional on survival up to generation $N$, the genealogy of the extant
population defines a random marked metric measure space, where individuals are
marked by their type and pairwise distances are measured by the time to the
most recent common ancestor. In the present manuscript, we devise a general
method of moments to prove convergence of such genealogies in the Gromov-weak
topology when $N \to \infty$.
  Informally, the moment of order $k$ of the population is obtained by
observing the genealogy of $k$ individuals chosen uniformly at random after
size-biasing the population at time $N$ by its $k$-th factorial moment. We show
that the sampled genealogy can be expressed in terms of a $k$-spine
decomposition of the original branching process, and that convergence reduces
to the convergence of the underlying $k$-spines.
  As an illustration of our framework, we analyse the large-time behavior of a
branching approximation of the biparental Wright-Fisher model with
recombination. The model exhibits some interesting mathematical features. It
starts in a supercritical state but is naturally driven to criticality. We show
that the limiting behavior exhibits both critical and supercritical
characteristics.


  Antibodies are canonically Y-shaped multimeric proteins capable of highly
specific molecular recognition. The CDRH3 region located at the tip of variable
chains of an antibody dominates antigen-binding specificity. Therefore, it is a
priority to design optimal antigen-specific CDRH3 regions to develop
therapeutic antibodies. However, the combinatorial nature of CDRH3 sequence
space makes it impossible to search for an optimal binding sequence
exhaustively and efficiently using computational approaches. Here, we present
\texttt{AntBO}: a combinatorial Bayesian optimisation framework enabling
efficient \textit{in silico} design of the CDRH3 region. Ideally, antibodies
are expected to have high target specificity and developability. We introduce a
CDRH3 trust region that restricts the search to sequences with favourable
developability scores to achieve this goal. For benchmarking, \texttt{AntBO}
uses the \texttt{Absolut!} software suite as a black-box oracle to score the
target specificity and affinity of designed antibodies \textit{in silico} in an
unconstrained fashion~\citep{robert2021one}. The experiments performed for
$159$ discretised antigens used in \texttt{Absolut!} demonstrate the benefit of
\texttt{AntBO} in designing CDRH3 regions with diverse biophysical properties.
In under $200$ calls to black-box oracle, \texttt{AntBO} can suggest antibody
sequences that outperform the best binding sequence drawn from 6.9 million
experimentally obtained CDRH3s and a commonly used genetic algorithm baseline.
Additionally, \texttt{AntBO} finds very-high affinity CDRH3 sequences in only
38 protein designs whilst requiring no domain knowledge. We conclude
\texttt{AntBO} brings automated antibody design methods closer to what is
practically viable for in vitro experimentation.


  Motivation: A branching processes model yields an unevenly stochastically
distributed dataset that consists of sparse and dense regions. This work
addresses the problem of precisely evaluating parameters for such a model.
Applying a branching processes model to an area such as cancer cell evolution
faces a number of obstacles, including high dimensionality and the rare
appearance of a result of interest. We take on the ambitious task of obtaining
the coefficients of a model that reflects the relationship of driver gene
mutations and cancer hallmarks on the basis of personal data regarding variant
allele frequencies. Results: An approximate Bayesian computation method based
on Isolation Kernel is developed. The method involves the transformation of row
data to a Hilbert space (mapping) and the measurement of the similarity between
simulated points and maxima weighted Isolation Kernel mapping related to the
observation point. We also design a heuristic algorithm for parameter
estimation that requires no calculation and is dimension independent. The
advantages of the proposed machine learning method are illustrated using
multidimensional test data as well as a specific example focused on cancer cell
evolution.


  With the recent advance of geometric deep learning, neural networks have been
extensively used for data in non-Euclidean domains. In particular, hyperbolic
neural networks have proved successful in processing hierarchical information
of data. However, many hyperbolic neural networks are numerically unstable
during training, which precludes using complex architectures. This crucial
problem makes it difficult to build hyperbolic generative models for real and
complex data. In this work, we propose a hyperbolic generative network in which
we design novel architecture and layers to improve stability in training. Our
proposed network contains three parts: first, a hyperbolic autoencoder (AE)
that produces hyperbolic embedding for input data; second, a hyperbolic
generative adversarial network (GAN) for generating the hyperbolic latent
embedding of the AE from simple noise; third, a generator that inherits the
decoder from the AE and the generator from the GAN. We call this network the
hyperbolic AE-GAN, or HAEGAN for short. The architecture of HAEGAN fosters
expressive representation in the hyperbolic space, and the specific design of
layers ensures numerical stability. Experiments show that HAEGAN is able to
generate complex data with state-of-the-art structure-related performance.