  While large language models (LLMs) have been successfully applied to various
tasks, they still face challenges with hallucinations. Augmenting LLMs with
domain-specific tools such as database utilities can facilitate easier and more
precise access to specialized knowledge. In this paper, we present GeneGPT, a
novel method for teaching LLMs to use the Web APIs of the National Center for
Biotechnology Information (NCBI) for answering genomics questions.
Specifically, we prompt Codex to solve the GeneTuring tests with NCBI Web APIs
by in-context learning and an augmented decoding algorithm that can detect and
execute API calls. Experimental results show that GeneGPT achieves
state-of-the-art performance on eight tasks in the GeneTuring benchmark with an
average score of 0.83, largely surpassing retrieval-augmented LLMs such as the
new Bing (0.44), biomedical LLMs such as BioMedLM (0.08) and BioGPT (0.04), as
well as GPT-3 (0.16) and ChatGPT (0.12). Our further analyses suggest that: (1)
API demonstrations have good cross-task generalizability and are more useful
than documentations for in-context learning; (2) GeneGPT can generalize to
longer chains of API calls and answer multi-hop questions in GeneHop, a novel
dataset introduced in this work; (3) Different types of errors are enriched in
different tasks, providing valuable insights for future improvements.


  We propose a general population dynamics model for two seagrass species
growing and interacting in two spatial dimensions. The model includes spatial
terms accounting for the clonal growth characteristics of seagrasses, and
coupling between species through the net mortality rate. We consider both
intraspecies and interspecies facilitative and competitive interactions,
allowing density-dependent interaction mechanisms. Here we study the case of
very similar species with reciprocal interactions, which allows reducing the
number of the model parameters to just four, and whose bifurcation structure
can be considered the backbone of the complete system. We find that the
parameter space can be divided into ten regions with qualitatively different
bifurcation diagrams. These regimes can be further grouped into just five
regimes with different ecological interpretations. Our analysis allows the
classifying of all possible density distributions and dynamical behaviors of
meadows with two coexisting species.


  Satellite DNA are long tandemly repeating sequences in a genome and may be
organized as high-order repeats (HORs). They are enriched in centromeres and
are challenging to assemble. Existing algorithms for identifying satellite
repeats either require the complete assembly of satellites or only work for
simple repeat structures without HORs. Here we describe Satellite Repeat Finder
(SRF), a new algorithm for reconstructing satellite repeat units and HORs from
accurate reads or assemblies without prior knowledge on repeat structures.
Applying SRF to real sequence data, we showed that SRF could reconstruct known
satellites in human and well-studied model organisms. We also found satellite
repeats are pervasive in various other species, accounting for up to 12% of
their genome contents but are often underrepresented in assemblies. With the
rapid progress on genome sequencing, SRF will help the annotation of new
genomes and the study of satellite DNA evolution even if such repeats are not
fully assembled.


  Evolutionary models are used to study the self-organisation of collective
action, often incorporating population structure due to its ubiquitous presence
and long-known impact on emerging phenomena. We investigate the evolution of
multiplayer cooperation in mobile structured populations, where individuals
move strategically on networks and interact with those they meet in groups of
variable size. We find that the evolution of multiplayer cooperation primarily
depends on the network topology and movement cost while using different
stochastic update rules seldom influences evolutionary outcomes. Cooperation
robustly co-evolves with movement on complete networks and structure has a
partially detrimental effect on it. These findings contrast an established
wisdom in evolutionary graph theory that cooperation can only emerge under some
update rules and if the average degree is low. We find that group-dependent
movement erases the locality of interactions, suppresses the impact of
evolutionary structural viscosity on the fitness of individuals, and leads to
assortative behaviour that is much more powerful than viscosity in promoting
cooperation. We analyse the differences remaining between update rules through
a comparison of evolutionary outcomes and fixation probabilities.


  Efficient coverage for newly developed vaccines requires knowing which groups
of individuals will accept the vaccine immediately and which will take longer
to accept or never accept. Of those who may eventually accept the vaccine,
there are two main types: success-based learners, basing their decisions on
others' satisfaction, and myopic rationalists, attending to their own immediate
perceived benefit. We used COVID-19 vaccination data to fit a mechanistic model
capturing the distinct effects of the two types on the vaccination progress. We
estimated that 47 percent of Americans behaved as myopic rationalist with a
high variations across the jurisdictions, from 31 percent in Mississippi to 76
percent in Vermont. The proportion was correlated with the vaccination
coverage, proportion of votes in favor of Democrats in 2020 presidential
election, and education score.


  This work brings the leading accuracy, sample efficiency, and robustness of
deep equivariant neural networks to the extreme computational scale. This is
achieved through a combination of innovative model architecture, massive
parallelization, and models and implementations optimized for efficient GPU
utilization. The resulting Allegro architecture bridges the accuracy-speed
tradeoff of atomistic simulations and enables description of dynamics in
structures of unprecedented complexity at quantum fidelity. To illustrate the
scalability of Allegro, we perform nanoseconds-long stable simulations of
protein dynamics and scale up to a 44-million atom structure of a complete,
all-atom, explicitly solvated HIV capsid on the Perlmutter supercomputer. We
demonstrate excellent strong scaling up to 100 million atoms and 70% weak
scaling to 5120 A100 GPUs.


  Understanding the emergent behavior of chemical reaction networks (CRNs) is a
fundamental aspect of biology and its origin from inanimate matter. A closed
CRN monotonically tends to thermal equilibrium, but when it is opened to
external reservoirs, a range of behaviors is possible, including transition to
a new equilibrium state, a non-equilibrium state, or indefinite growth. This
study shows that slowly driven CRNs are governed by the conserved quantities of
the closed system, which are generally far fewer in number than the species.
Considering both deterministic and stochastic dynamics, a universal slow
dynamics equation is derived with singular perturbation methods, and is shown
to be thermodynamically consistent. The slow dynamics is highly robust against
microscopic details of the network, which may be unknown in practical
situations. In particular, non-equilibrium states of realistic large CRNs can
be sought without knowledge of bulk reaction rates. The framework is
successfully tested against a suite of networks of increasing complexity and
argued to be relevant in the treatment of open CRNs as chemical machines.


  Motivated by long-range dispersal in ecological systems, we formulate and
apply a general strong-disorder renormalization group (SDRG) framework to
describe one-dimensional disordered contact processes with heavy-tailed, such
as power law, stretched exponential, and log-normal dispersal kernels, widely
used in ecology. The focus is on the close-to-critical scaling of the order
parameters, including the commonly used density, as well as the less known
persistence, which is non-zero in the inactive phase. Our analytic and
numerical results obtained by SDRG schemes at different levels of approximation
reveal that the more slowly decaying dispersal kernels lead to more smoothly
vanishing densities as the critical point is approached. The persistence,
however, shows an opposite tendency: the broadening of the dispersal makes its
decline more singular at the critical point, becoming discontinuous for the
extreme case of power-law dispersal. The SDRG schemes presented here also
describe the quantum phase transition of random transverse-field Ising chains
with ferromagnetic long-range interactions, the density corresponding to the
magnetization of this model.


  In this article, a model, discrete in space and time, is developed to
describe the growth and spread of the Scottish populations of red squirrels
(Sciurus vulgaris) and grey squirrel (Sciurus carolinensis). The initial state
for the model is designed using a large dataset of records of sightings of
individuals of both species reported by members of the public. Choices of
parameters involved in the model and their values are informed by the analysis
of this dataset for the period 2011-2016, and model predictions are compared to
records for the years 2006-2019.


  Conformation Generation is a fundamental problem in drug discovery and
cheminformatics. And organic molecule conformation generation, particularly in
vacuum and protein pocket environments, is most relevant to drug design.
Recently, with the development of geometric neural networks, the data-driven
schemes have been successfully applied in this field, both for molecular
conformation generation (in vacuum) and binding pose generation (in protein
pocket). The former beats the traditional ETKDG method, while the latter
achieves similar accuracy compared with the widely used molecular docking
software. Although these methods have shown promising results, some researchers
have recently questioned whether deep learning (DL) methods perform better in
molecular conformation generation via a parameter-free method. To our surprise,
what they have designed is some kind analogous to the famous infinite monkey
theorem, the monkeys that are even equipped with physics education. To discuss
the feasibility of their proving, we constructed a real infinite stochastic
monkey for molecular conformation generation, showing that even with a more
stochastic sampler for geometry generation, the coverage of the benchmark
QM-computed conformations are higher than those of most DL-based methods. By
extending their physical monkey algorithm for binding pose prediction, we also
discover that the successful docking rate also achieves near-best performance
among existing DL-based docking models. Thus, though their conclusions are
right, their proof process needs more concern.


  We explore the inclusion of vaccination in compartmental epidemiological
models concerning the delta and omicron variants of the SARS-CoV-2 virus that
caused the COVID-19 pandemic. We expand on our earlier compartmental-model work
by incorporating vaccinated populations. We present two classes of models that
differ depending on the immunological properties of the variant. The first one
is for the delta variant, where we do not follow the dynamics of the vaccinated
individuals since infections of vaccinated individuals were rare. The second
one for the far more contagious omicron variant incorporates the evolution of
the infections within the vaccinated cohort. We explore comparisons with
available data involving two possible classes of counts, fatalities and
hospitalizations. We present our results for two regions, Andalusia and
Switzerland (including the Principality of Liechtenstein), where the necessary
data are available. In the majority of the considered cases, the models are
found to yield good agreement with the data and have a reasonable predictive
capability beyond their training window, rendering them potentially useful
tools for the interpretation of the COVID-19 and further pandemic waves, and
for the design of intervention strategies during these waves.


  A fundamental problem in ecology is to understand how competition shapes
biodiversity and species coexistence. Historically, one important approach for
addressing this question has been to analyze consumer resource models using
geometric arguments. This has led to broadly applicable principles such as
Tilman's $R^*$ and species coexistence cones. Here, we extend these arguments
by constructing a novel geometric framework for understanding species
coexistence based on convex polytopes in the space of consumer preferences. We
show how the geometry of consumer preferences can be used to predict species
which may coexist and enumerate ecologically-stable steady states and
transitions between them. Collectively, these results provide a framework for
understanding the role of species traits within niche theory.


  Phylogenetic diversity is a popular measure for quantifying the biodiversity
of a collection $Y$ of species, while phylogenetic diversity indices provide a
way to apportion phylogenetic diversity to individual species. Typically, for
some specific diversity index, the phylogenetic diversity of $Y$ is not equal
to the sum of the diversity indices of the species in $Y.$ In this paper, we
investigate the extent of this difference for two commonly-used indices: Fair
Proportion and Equal Splits. In particular, we determine the maximum value of
this difference under various instances including when the associated rooted
phylogenetic tree is allowed to vary across all root phylogenetic trees with
the same leaf set and whose edge lengths are constrained by either their total
sum or their maximum value.


  Computationally generating novel synthetically accessible compounds with high
affinity and low toxicity is a great challenge in drug design. Machine-learning
models beyond conventional pharmacophoric methods have shown promise in
generating novel small molecule compounds, but require significant tuning for a
specific protein target. Here, we introduce a method called selective iterative
latent variable refinement (SILVR) for conditioning an existing diffusion-based
equivariant generative model without retraining. The model allows the
generation of new molecules that fit into a binding site of a protein based on
fragment hits. We use the SARS-CoV-2 Main protease fragments from Diamond
X-Chem that form part of the COVID Moonshot project as a reference dataset for
conditioning the molecule generation. The SILVR rate controls the extent of
conditioning and we show that moderate SILVR rates make it possible to generate
new molecules of similar shape to the original fragments, meaning that the new
molecules fit the binding site without knowledge of the protein. We can also
merge up to 3 fragments into a new molecule without affecting the quality of
molecules generated by the underlying generative model. Our method is
generalizable to any protein target with known fragments and any
diffusion-based model for molecule generation.


  Large pre-trained language models (LLMs) have been shown to have significant
potential in few-shot learning across various fields, even with minimal
training data. However, their ability to generalize to unseen tasks in more
complex fields, such as biology, has yet to be fully evaluated. LLMs can offer
a promising alternative approach for biological inference, particularly in
cases where structured data and sample size are limited, by extracting prior
knowledge from text corpora. Our proposed few-shot learning approach uses LLMs
to predict the synergy of drug pairs in rare tissues that lack structured data
and features. Our experiments, which involved seven rare tissues from different
cancer types, demonstrated that the LLM-based prediction model achieved
significant accuracy with very few or zero samples. Our proposed model, the
CancerGPT (with $\sim$ 124M parameters), was even comparable to the larger
fine-tuned GPT-3 model (with $\sim$ 175B parameters). Our research is the first
to tackle drug pair synergy prediction in rare tissues with limited data. We
are also the first to utilize an LLM-based prediction model for biological
reaction prediction tasks.


  Antigenic variation is the main immune escape mechanism for RNA viruses like
influenza or SARS-CoV-2. While high mutation rates promote antigenic escape,
they also induce large mutational loads and reduced fitness. It remains unclear
how this cost-benefit trade-off selects the mutation rate of viruses. Using a
traveling wave model for the co-evolution of viruses and host immune systems in
a finite population, we investigate how immunity affects the evolution of the
mutation rate and other non-antigenic traits, such as virulence. We first show
that the nature of the wave depends on how cross-reactive immune systems are,
reconciling previous approaches. The immune-virus system behaves like a Fisher
wave at low cross-reactivities, and like a fitness wave at high
cross-reactivities. These regimes predict different outcomes for the evolution
of non-antigenic traits. At low cross-reactivities, the evolutionarily stable
strategy is to maximize the speed of the wave, implying a higher mutation rate
and increased virulence. At large cross-reactivities, where our estimates place
H3N2 influenza, the stable strategy is to increase the basic reproductive
number, keeping the mutation rate to a minimum and virulence low.


  Many researchers have been encouraged to describe the dynamics of chromosomal
loci in chromatin using the classical Rouse model of polymer dynamics by the
agreement between the measured mean square displacement (MSD) versus time of
fluorescently-labelled loci and the Rouse-model predictions. However, the
discovery of intermediate-scale chromatin organization, known as topologically
associating domains (TADs), together with the proposed explanation of TADs in
terms of chromatin loops and loop extrusion, is at odds with the classical
Rouse model, which does not contain loops. Accordingly, we introduce an
extended Rouse model that incorporates chromatin loop configurations from
loop-extrusion-factor-model simulations. Specifically, we extend the classical
Rouse model by modifying the polymer's dynamical matrix to incorporate extra
springs that represent loop bases. We also theoretically generalize the
friction coefficient matrix so that the Rouse beads with non-uniform friction
coefficients are compatible with our Rouse model simulation method. This
extended Rouse model allowes us to investigate the impact of loops and loop
extrusion on the dynamics of chromatin. We show that loops significantly
suppress the averaged MSD of a chromosomal locus, consistent with recent
experiments that track fluorescently-labelled chromatin loci in fission yeast
[M. L. P. Bailey, I. Surovtsev, J. F. Williams, H. Yan, T. Yuan, S. G. Mochrie,
and M. C. King, Mol. Biol. Cell (in press)]. We also find that loops slightly
reduce the MSD's stretching exponent from the classical Rouse-model value of
0.5 to a loop-density-dependent value in the 0.45-0.40 range. Remarkably,
stretching exponent values in this range have also been reported in recent
experiments [S. C. Weber, A. J. Spakowitz, and J. A. Theriot, Phys. Rev. Lett.
104, 238102 (2010) and Bailey et al., Mol. Biol. Cell (in press)].


  Biological systems must be robust for stable function against perturbations,
but robustness alone is not sufficient. The ability to switch between
appropriate states (phenotypes) in response to different conditions is
essential for biological functions. How are robustness and plasticity
simultaneously acquired through evolution? We examine the evolution of
genotypes that realize plastic switching between two phenotypes upon external
inputs, as well as stationary expressions of phenotypes. We introduce a
statistical physics model consisting of spins, with active and regulatory
sites, which are distinct from each other. We represent the phenotype and
genotype as spin configurations and the spin-spin interactions, respectively.
The fitness for selection is given so that it takes a higher value as more of
the active sites take two requested spin configurations depending on the
regulation. We numerically evolve the interaction matrix by changing them with
mutations and selection of those with higher fitness. Our numerical simulations
show that characteristic genotypes evolve slightly above the transition
temperature between replica symmetric and replica symmetry breaking phase.
These genotypes shape two spin configurations separately depending on the
regulation, where the two phenotypes are dominantly represented by the
genotypes' first and second eigenmodes, and smooth switching of two phenotypes
are achieved by following a one-dimensional path connecting the two phenotypes.
Upon changes in regulations, spin configurations are attracted to this path,
which allows for robust and plastic switching between the two phenotypes. The
statistical-physics analysis show that the free energy landscape has a valley
along the switching path. Our finding indicates that the compatibility of the
robustness and plasticity is acquired by the evolution of the
low-dimensionality in the phenotype space.


  In exponential population growth, variability in the timing of individual
division events and environmental factors (including stochastic inoculation)
compound to produce variable growth trajectories. In several stochastic models
of exponential growth we show power-law relationships that relate variability
in the time required to reach a threshold population size to growth rate and
inoculum size. Population-growth experiments in E. coli and S. aureus with
inoculum sizes ranging between 1 and 100 are consistent with these
relationships. We quantify how noise accumulates over time, finding that it
encodes -- and can be used to deduce -- information about the early growth rate
of a population.


  Most genes are part of larger families of evolutionary related genes. The
history of gene families typically involves duplications and losses of genes as
well as horizontal transfers into other organisms. The reconstruction of
detailed gene family histories, i.e., the precise dating of evolutionary events
relative to phylogenetic tree of the underlying species has remained a
challenging topic despite their importance as a basis for detailed
investigations into adaptation and functional evolution of individual members
of the gene family. The identification of orthologs, moreover, is a
particularly important subproblem of the more general setting considered here.
In the last few years, an extensive body of mathematical results has appeared
that tightly links orthology, a formal notion of best matches among genes, and
horizontal gene transfer. The purpose of this chapter is the broadly outline
some of the key mathematical insights and to discuss their implication for
practical applications. In particular, we focus on tree-free methods, i.e.,
methods to infer orthology or horizontal gene transfer as well as gene trees,
species trees and reconciliations between them without using \emph{a priori}
knowledge of the underlying trees or statistical models for the inference of
phylogenetic trees. Instead, the initial step aims to extract binary relations
among genes.


  Epidemic forecasts are only as good as the accuracy of epidemic measurements.
Is epidemic data, particularly COVID-19 epidemic data, clean and devoid of
noise? Common sense implies the negative answer. While we cannot evaluate the
cleanliness of the COVID-19 epidemic data in a holistic fashion, we can assess
the data for the presence of reporting delays. In our work, through the
analysis of the first COVID-19 wave, we find substantial reporting delays in
the published epidemic data. Motivated by the desire to enhance epidemic
forecasts, we develop a statistical framework to detect, uncover, and remove
reporting delays in the infectious, recovered, and deceased epidemic time
series. Our framework can uncover and analyze reporting delays in 8 regions
significantly affected by the first COVID-19 wave. Further, we demonstrate that
removing reporting delays from epidemic data using our statistical framework
may decrease the error in epidemic forecasts. While our statistical framework
can be used in combination with any epidemic forecast method that intakes
infectious, recovered, and deceased data, to make a basic assessment, we
employed the classical SIRD epidemic model. Our results indicate that the
removal of reporting delays from the epidemic data may decrease the forecast
error by up to 50. We anticipate that our framework will be indispensable in
the analysis of novel COVID-19 strains and other existing or novel infectious
diseases.


  Disease propagation between countries strongly depends on their effective
distance, a measure derived from the world air transportation network (WAN). It
reduces the complex spreading patterns of a pandemic to a wave-like propagation
from the outbreak country, i.e. a linear relationship to the arrival time of
the unmitigated spread of a disease. However, in the early stage of an
outbreak, what matters to countries' decision makers is knowledge about the
relative risk of arrival of active cases, i.e. how likely it is that an active
case that boarded at the outbreak location will arrive in their country. As
accurate mechanistic models to estimate such risks are still lacking, we
propose here the ``import risk'' model that defines an import probability by
means of the effective-distance framework. The model assumes that airline
passengers are distributed along the shortest path tree that starts at the
outbreak's origin. In combination with a random walk, we account for all
possible paths, thus inferring predominant connecting flights. Our model
outperforms other mobility models, such as the radiation and gravity model with
varying distance types, and it improves further if additional geographic
information is included. The import risk model's precision increases for
countries that are more connected within the WAN, and recovers a geographic
distance-dependence that suggests a pull- rather than a push- dynamic of the
distribution process.


  Cats use the connection governing parallel transport in the space of shapes
to land safely on their feet. Here we argue that this connection also explains
the impressive performance of molecular motors by enabling molecules to evade
conclusions of Feynman's ratchet-and-pawl analysis. We first demonstrate, using
simple molecular models, how directed rotational motion can emerge from shape
changes even without angular momentum. We then computationally design knotted
polyalanine molecules and show how their shape space connection organizes
individual atom thermal vibrations into collective rotational motion,
independently of angular momentum. Our simulations show that rotational motion
arises effortlessly even in ambient water, making the molecule an effective
theory time crystal. Our findings have potential for practical molecular motor
design and engineering and can be verified through high-precision nuclear
magnetic resonance measurements.


  Recently deep learning based quantitative structure-activity relationship
(QSAR) models has shown surpassing performance than traditional methods for
property prediction tasks in drug discovery. However, most DL based QSAR models
are restricted to limited labeled data to achieve better performance, and also
are sensitive to model scale and hyper-parameters. In this paper, we propose
Uni-QSAR, a powerful Auto-ML tool for molecule property prediction tasks.
Uni-QSAR combines molecular representation learning (MRL) of 1D sequential
tokens, 2D topology graphs, and 3D conformers with pretraining models to
leverage rich representation from large-scale unlabeled data. Without any
manual fine-tuning or model selection, Uni-QSAR outperforms SOTA in 21/22 tasks
of the Therapeutic Data Commons (TDC) benchmark under designed parallel
workflow, with an average performance improvement of 6.09\%. Furthermore, we
demonstrate the practical usefulness of Uni-QSAR in drug discovery domains.


  Understanding the host-specificity of different families of viruses sheds
light on the origin of, e.g., SARS-CoV-2, rabies, and other such zoonotic
pathogens in humans. It enables epidemiologists, medical professionals, and
policymakers to curb existing epidemics and prevent future ones promptly. In
the family Coronaviridae (of which SARS-CoV-2 is a member), it is well-known
that the spike protein is the point of contact between the virus and the host
cell membrane. On the other hand, the two traditional mammalian orders,
Carnivora (carnivores) and Chiroptera (bats) are recognized to be responsible
for maintaining and spreading the Rabies Lyssavirus (RABV). We propose
Virus2Vec, a feature-vector representation for viral (nucleotide or amino acid)
sequences that enable vector-space-based machine learning models to identify
viral hosts. Virus2Vec generates numerical feature vectors for unaligned
sequences, allowing us to forego the computationally expensive sequence
alignment step from the pipeline. Virus2Vec leverages the power of both the
\emph{minimizer} and position weight matrix (PWM) to generate compact feature
vectors. Using several classifiers, we empirically evaluate Virus2Vec on
real-world spike sequences of Coronaviridae and rabies virus sequence data to
predict the host (identifying the reservoirs of infection). Our results
demonstrate that Virus2Vec outperforms the predictive accuracies of baseline
and state-of-the-art methods.


  The aim of this manuscript is to contain the arguments and define the
theoretical objects for building a general framework to model population
dynamics from the ground up, relying mainly on the probabilistic landscapes
defining the dynamics instead of the context-dependent physical specification
of systems. I intend to keep updating and correcting this manuscript. The goal
is for all the different parts to be able to communicate with each other and
for models to be directly comparable and to maintain an explicit connection to
the first principles sustaining them. This modeling paradigm will stem from a
Bayesian perspective on model definition and interpretation and will be
primarily concerned with ecological and evolutionary processes. Populations are
considered to be abstract collections of elements that relate in the same ways,
and the laws of motion ultimately depend on relational properties of elements,
at first irrespective of their constitution. The states of populations are
taken to be their spatial densities, the fundamental quantities shaping the
dynamics of their interactions.


  Designing molecules with desirable physiochemical properties and
functionalities is a long-standing challenge in chemistry, material science,
and drug discovery. Recently, machine learning-based generative models have
emerged as promising approaches for \emph{de novo} molecule design. However,
further refinement of methodology is highly desired as most existing methods
lack unified modeling of 2D topology and 3D geometry information and fail to
effectively learn the structure-property relationship for molecule design. Here
we present MolCode, a roto-translation equivariant generative framework for
\underline{Mol}ecular graph-structure \underline{Co-de}sign. In MolCode, 3D
geometric information empowers the molecular 2D graph generation, which in turn
helps guide the prediction of molecular 3D structure. Extensive experimental
results show that MolCode outperforms previous methods on a series of
challenging tasks including \emph{de novo} molecule design, targeted molecule
discovery, and structure-based drug design. Particularly, MolCode not only
consistently generates valid (99.95$\%$ Validity) and diverse (98.75$\%$
Uniqueness) molecular graphs/structures with desirable properties, but also
generate drug-like molecules with high affinity to target proteins (61.8$\%$
high-affinity ratio), which demonstrates MolCode's potential applications in
material design and drug discovery. Our extensive investigation reveals that
the 2D topology and 3D geometry contain intrinsically complementary information
in molecule design, and provide new insights into machine learning-based
molecule representation and generation.


  Binary phylogenetic trees inferred from biological data are central to
understanding the shared evolutionary history of organisms. Inferring the
placement of latent nodes in a tree by any optimality criterion (e.g., maximum
likelihood) is an NP-hard problem, propelling the development of myriad
heuristic approaches. Yet, these heuristics often lack a systematic means of
uniformly sampling random trees or effectively exploring a tree space that
grows factorially, which are crucial to optimisation problems such as machine
learning. Accordingly, we present Phylo2Vec, a new parsimonious representation
of a phylogenetic tree. Phylo2Vec maps any binary tree with $n$ leaves to an
integer vector of length $n$. We prove that Phylo2Vec is both well-defined and
bijective to the space of phylogenetic trees. The advantages of Phylo2Vec are
twofold: i) easy uniform sampling of binary trees and ii) systematic ability to
traverse tree space in very large or small jumps. As a proof of concept, we use
Phylo2Vec for maximum likelihood inference on five real-world datasets and show
that a simple hill climbing-based optimisation efficiently traverses the
vastness of tree space from a random to an optimal tree.


  Designing molecules that bind to specific target proteins is a fundamental
task in drug discovery. Recent models leverage geometric constraints to
generate ligand molecules that bind cohesively with specific protein pockets.
However, these models cannot effectively generate 3D molecules with 2D skeletal
curtailments and property constraints, which are pivotal to drug potency and
development. To tackle this challenge, we propose GraphVF, a variational
flow-based framework that combines 2D topology and 3D geometry, for
controllable generation of binding 3D molecules. Empirically, our method
achieves state-of-the-art binding affinity and realistic sub-structural layouts
for protein-specific generation. In particular, GraphVF represents the first
controllable geometry-aware, protein-specific molecule generation method, which
can generate binding 3D molecules with tailored sub-structures and
physio-chemical properties. Our code is available at
https://github.com/Franco-Solis/GraphVF-code.


  This study uses the Lyapunov method, the Poincar\'e-Bendixson theorem, and
the Dulac criterion to analyze the stability of SIR and SIRS with non-monotone
incidence and different mortality rates.


  In recent years, several machine learning approaches have been proposed to
predict gene expression and epigenetic signals from the DNA sequence alone.
These models are often used to deduce, and, to some extent, assess putative new
biological insights about gene regulation, and they have led to very
interesting advances in regulatory genomics. This article reviews a selection
of these methods, ranging from linear models to random forests, kernel methods,
and more advanced deep learning models. Specifically, we detail the different
techniques and strategies that can be used to extract new gene-regulation
hypotheses from these models. Furthermore, because these putative insights need
to be validated with wet-lab experiments, we emphasize that it is important to
have a measure of confidence associated with the extracted hypotheses. We
review the procedures that have been proposed to measure this confidence for
the different types of machine learning models, and we discuss the fact that
they do not provide the same kind of information.


  Single-cell analysis is currently one of the most high-resolution techniques
to study biology. The large complex datasets that have been generated have
spurred numerous developments in computational biology, in particular the use
of advanced statistics and machine learning. This review attempts to explain
the deeper theoretical concepts that underpin current state-of-the-art analysis
methods. Single-cell analysis is covered from cell, through instruments, to
current and upcoming models. A minimum of mathematics and statistics has been
used, but the reader is assumed to either have basic knowledge of single-cell
analysis workflows, or have a solid knowledge of statistics. The aim of this
review is to spread concepts which are not yet in common use, especially from
topology and generative processes, and how new statistical models can be
developed to capture more of biology. This opens epistemological questions
regarding our ontology and models, and some pointers will be given to how
natural language processing (NLP) may help overcome our cognitive limitations
for understanding single-cell data.


  There are two main life cycles in plants, annual and perennial. These life
cycles are associated with different traits that determine ecosystem function.
Although life cycles are textbook examples of plant adaptation to different
environments, we lack comprehensive knowledge regarding their global
distributional patterns. Here, we assembled an extensive database of plant life
cycle assignments of 235,000 plant species coupled with millions of
georeferenced data points to map the worldwide biogeography thereof. We found
that annuals are half as common as initially thought, accounting for only 6% of
plant species. Our analyses indicate annuals are favored in hot and dry
regions. However, a more accurate model shows annual species' prevalence is
driven by temperature and precipitation in the driest quarter (rather than
yearly means), explaining, for example, why some Mediterranean systems have
more annuals than deserts. Furthermore, this pattern remains consistent among
different families, indicating convergent evolution. Finally, we demonstrate
that increasing climate variability and anthropogenic disturbance increase
annual favorability. Considering future climate change, we predict an increase
in annual prevalence for 69% of the world's ecoregions by 2060. Overall, our
analyses raise concerns for ecosystem services provided by perennials, as
ongoing changes are leading to a more annuals-dominated world.


  The spatial positioning of chromosomes relative to functional nuclear bodies
is intertwined with genome functions such as transcription. However, the
sequence patterns and epigenomic features that collectively influence chromatin
spatial positioning in a genome-wide manner are not well understood. Here, we
develop a new transformer-based deep learning model called UNADON, which
predicts the genome-wide cytological distance to a specific type of nuclear
body, as measured by TSA-seq, using both sequence features and epigenomic
signals. Evaluations of UNADON in four cell lines (K562, H1, HFFc6, HCT116)
show high accuracy in predicting chromatin spatial positioning to nuclear
bodies when trained on a single cell line. UNADON also performed well in an
unseen cell type. Importantly, we reveal potential sequence and epigenomic
factors that affect large-scale chromatin compartmentalization to nuclear
bodies. Together, UNADON provides new insights into the principles between
sequence features and large-scale chromatin spatial localization, which has
important implications for understanding nuclear structure and function.


  Genetic algorithm (GA) is inspired by biological evolution of genetic
organisms by optimizing the genotypic combinations encoded within each
individual with the help of evolutionary operators, suggesting that GA may be a
suitable model for studying real-life evolutionary processes. This paper
describes the design of a Python library for artificial life simulation,
Digital Organism Simulation Environment (DOSE), based on GA and biological
hierarchy starting from genetic sequence to population. A 3-character
instruction set that does not take any operand is introduced as genetic code
for digital organism. This mimics the 3-nucleotide codon structure in naturally
occurring DNA. In addition, the context of a 3-dimensional world composing of
ecological cells is introduced to simulate a physical ecosystem. Using DOSE, an
experiment to examine the changes in genetic sequences with respect to mutation
rates is presented.


  We formulate an integrated pest management model to control natural pests of
the crop through the periodic application of biopesticide and chemical
pesticides. In a theoretical analysis of the system pest eradication, a
periodic solution is found and established. All the system variables are proved
to be bounded. Our main goal is then to ensure that pesticides are optimized,
in terms of pesticide concentration and pesticide application frequency, and
that the optimum combination of pesticides is found to provide the most benefit
to the crop. By using Floquet theory and the small amplitude perturbation
method, we prove that the pest eradication periodic solution is locally and
globally stable. The acquired results establish a threshold time limit for the
impulsive release of various controls as well as some valid theoretical
conclusions for effective pest management. Furthermore, after a numerical
comparison, we conclude that integrated pest management is more effective than
single biological or chemical controls. Finally, we illustrate the analytical
results through numerical simulations.


  In this paper, we discover that the class of random polynomials arising from
the equilibrium analysis of random asymmetric evolutionary games is
\textit{exactly} the Kostlan-Shub-Smale system of random polynomials, revealing
an intriguing connection between evolutionary game theory and the theory of
random polynomials. Through this connection, we analytically characterize the
statistics of the number of internal equilibria of random asymmetric
evolutionary games, namely its mean value, probability distribution, central
limit theorem and universality phenomena. Moreover, comparing symmetric and
asymmetric random games, we establish that symmetry in the group interaction
increases the expected number of internal equilibria. Our research establishes
new theoretical understanding of asymmetric evolutionary games and why
symmetry/asymmetry of group interactions matters.


  Consider a birth and death process started from one individual in which each
individual gives birth at rate $\lambda$ and dies at rate $\mu$, so that the
population size grows at rate $r = \lambda - \mu$. Lambert and Harris,
Johnston, and Roberts came up with methods for constructing the exact genealogy
of a sample of size $n$ taken from this population at time $T$. We use the
construction of Lambert, which is based on the coalescent point process, to
obtain asymptotic results for the site frequency spectrum associated with this
sample. In the supercritical case $r > 0$, our results extend results of
Durrett for exponentially growing populations. In the critical case $r = 0$,
our results parallel those that Dahmer and Kersting obtained for Kingman's
coalescent.


  Modelling host behavioral change in response to epidemics is important to
describe disease dynamics and many previous studies proposed mathematical
models describing it. Indeed, the epidemic of COVID-19 clearly demonstrated
that people changed their activity in response to the epidemic, which
subsequently modified the disease dynamics. To predict the behavioral change
relevant to the disease dynamics, we need to know the epidemic situation (e.g.,
the number of reported cases) at the moment of decision to change behavior.
However, it is difficult to identify the timing of decision-making. In this
study, we analyzed travel accommodation reservation data in four prefectures of
Japan to observe decision-making timings and how it responded to the changing
epidemic situation during Japan's Coronavirus Disease 2019 (eight waves until
February 2023). To this end, we defined 'mobility avoidance index' to indicate
people's decision of mobility avoidance and quantified it using the time-series
of the accommodation booking/cancellation data. Our analysis revealed
semi-quantitative rules for day-to-day decision-making of human mobility under
a given epidemic situation. We observed matches of the peak dates of the index
and the number of reported cases. Additionally, we found that mobility
avoidance index increased/decreased linearly with the logarithmic number of
reported cases during the first epidemic wave. This pattern agrees with
Weber-Fechner law in psychophysics. We also found that the slope of the
mobility avoidance index against the change of the logarithmic number of
reported cases were similar among the waves, while the intercept of that was
much reduced as the first epidemic wave passed by. It suggests that the
people's response became weakened after the first experience, as if the number
of reported cases were multiplied by a constant small factor.


  Molecule generation is a very important practical problem, with uses in drug
discovery and material design, and AI methods promise to provide useful
solutions. However, existing methods for molecule generation focus either on 2D
graph structure or on 3D geometric structure, which is not sufficient to
represent a complete molecule as 2D graph captures mainly topology while 3D
geometry captures mainly spatial atom arrangements. Combining these
representations is essential to better represent a molecule. In this paper, we
present a new model for generating a comprehensive representation of molecules,
including atom features, 2D discrete molecule structures, and 3D continuous
molecule coordinates, by combining discrete and continuous diffusion processes.
The use of diffusion processes allows for capturing the probabilistic nature of
molecular processes and exploring the effect of different factors on molecular
structures. Additionally, we propose a novel graph transformer architecture to
denoise the diffusion process. The transformer adheres to 3D roto-translation
equivariance constraints, allowing it to learn invariant atom and edge
representations while preserving the equivariance of atom coordinates. This
transformer can be used to learn molecular representations robust to geometric
transformations. We evaluate the performance of our model through experiments
and comparisons with existing methods, showing its ability to generate more
stable and valid molecules. Our model is a promising approach for designing
stable and diverse molecules and can be applied to a wide range of tasks in
molecular modeling.


  The origin of cellular life can be described in terms of the transition from
inorganic matter: solids, liquids and gases, to the emergence of cooperative
assemblies of organic matter, DNA and proteins,capable of replication and
metabolism. Directionality Theory is a mathematical model of the collective
behavior of populations of organic matter: cells and higher organisms.
Evolutionary entropy, the cornerstone of the theory, is a statistical measure
of the cooperativity of the interacting components that comprise the
population. The main tenet of Directionality Theory is the Entropic Principle
of Collective Behavior: The collective behavior of aggregates of organic matter
is contingent on the population size and the external energy source, and
characterized by extremal states of evolutionary entropy. This article invokes
Directionality Theory to provide an evolutionary rationale for the following
sequence of transformations which define the emergence of cellular life: 1. The
self-assembly of activated macromolecules from inorganic matter 2. The
emergence of an RNA world, defined by RNA molecules with catalytic and
replicative properties 3. The origin of cellular life, the integration of the
three carbon-based polymers: DNA, proteins and lipids, to generate a metabolic
and replicative unit.


  Self-organization is the autonomous assembly of a network of interacting
components into a stable, organized pattern. This article shows that the
process of self-assembly can be encoded in terms of evolutionary entropy, a
statistical measure of the cooperativity of the interacting components.
Evolutionary entropy describes the rate at which a network of interacting
metabolic units convert an external energy source into mechanical energy and
work. We invoke Directionality Theory, an analytic model of Darwinian evolution
to analyze self-assembly as a variation-selection process, and to derive a
general tenet, namely, the Entropic Principle of Self-Organization: The
equilibrium states of a self-organizing process are states which maximize
evolutionary entropy, contingent on the production rate of the external energy
source. This principle is a universal rule, applicable to the self-assembly of
structures ranging from the folding of proteins, to branching morphogenesis,
and the emergence of social organization. The principle also elucidates the
origin of cellular life: the transition from inorganic matter to the emergence
of cells, capable of replication and metabolism.


  The stable cooperation ratio of spatial evolutionary games has been widely
studied using simulations or approximate analysis methods. However, sometimes
such ``stable'' cooperation ratios obtained via approximate methods might not
be actually stable, but correspond to quasi-equilibriums instead. We find that
various classic game models, like the evolutionary snowdrift game, evolutionary
prisoner's dilemma, and spatial public goods game on square lattices and
scale-free networks, exhibit the phase transition in convergence time to the
equilibrium state. Moreover, mathematical principles are provided to explain
the phase transition of convergence time and quasi-equilibrium of cooperation
ratio. The findings explain why and when cooperation and defection have a
long-term coexistence.


  In biology, predicting RNA secondary structures plays a vital role in
determining its physical and chemical properties. Although we have powerful
energy models to predict them as well as parametric analysis to understand the
models themselves, the large number of parameters involved makes exploring the
parameter space and effective fine-tuning complicated at best. The literature
describes an approach via so-called RNA polytopes and several attempts to
compute them entirely, but computing explicitly the polytopes is both
practically and theoretically intractable. In this thesis, we demonstrate how
to further modify the dynamic programming algorithms used in RNA secondary
structure prediction, and more generally how to use only supporting functions
to gather some information about the polytopes without explicit construction.
We provide the mathematical frameworks with proofs or sketch thereof whenever
necessary, and carry out some numerical experiments to show that our proposed
methods are practical even when the number of parameters is large. As it turns
out, one of our methods provides a solution to another problem in computational
geometry previously unsolved to our knowledge, and we hope this thesis will
accommodate future studies in RNA, as well as inspire further researches on the
potential uses of polytopes' supporting functions in computational geometry.


  Spanning two decades, the Encyclopaedia of DNA Elements (ENCODE) is a
collaborative research project that aims to identify all the functional
elements in the human and mouse genomes. To best serve the scientific
community, all data generated by the consortium is shared through a web-portal
(https://www.encodeproject.org/) with no access restrictions. The fourth and
final phase of the project added a diverse set of new samples (including those
associated with human disease), and a wide range of new assays aimed at
detection, characterization and validation of functional genomic elements. The
ENCODE data portal hosts results from over 23,000 functional genomics
experiments, over 800 functional elements characterization experiments
(including in vivo transgenic enhancer assays, reporter assays and CRISPR
screens) along with over 60,000 results of computational and integrative
analyses (including imputations, predictions and genome annotations). The
ENCODE Data Coordination Center (DCC) is responsible for development and
maintenance of the data portal, along with the implementation and utilisation
of the ENCODE uniform processing pipelines to generate uniformly processed
data. Here we report recent updates to the data portal. Specifically, we have
completely redesigned the home page, improved search interface, added several
new pages to highlight collections of biologically related data (deeply
profiled cell lines, immune cells, Alzheimer's Disease, RNA-Protein
interactions, degron matrix and a matrix of experiments organised by human
donors), added single-cell experiments, and enhanced the cart interface for
visualisation and download of user-selected datasets.


  Single-cell transcriptomics techniques, such as scRNA-seq, attempt to
characterize gene expression profiles in each cell of a heterogeneous sample
individually. Due to growing amounts of data generated and the increasing
complexity of the computational protocols needed to process the resulting
datasets, the demand for dedicated training in mathematical and programming
skills may preclude the use of these powerful techniques by many teams.
  In order to help close that gap between wet-lab and dry-lab capabilities we
have developed SinglePointRNA, a shiny-based R application that provides a
graphic interface for different publicly available tools to analyze single cell
RNA-seq data.
  The aim of SinglePointRNA is to provide an accessible and transparent tool
set to researchers that allows them to perform detailed and custom analysis of
their data autonomously. SinglePointRNA is structured in a context-driven
framework that prioritizes providing the user with solid qualitative guidance
at each step of the analysis process and interpretation of the results.
Additionally, the rich user guides accompanying the software are intended to
serve as a point of entry for users to learn more about computational
techniques applied to single cell data analysis.
  The SinglePointRNA app, as well as case datasets for the different tutorials
are available at www.github.com/ScienceParkMadrid/SinglePointRNA


  Phylogenetic networks are increasingly being considered as better suited to
represent the complexity of the evolutionary relationships between species. One
class of phylogenetic networks that has received a lot of attention recently is
the class of orchard networks, which is composed of networks that can be
reduced to a single leaf using cherry reductions. Cherry reductions, also
called cherry-picking operations, remove either a leaf of a simple cherry
(sibling leaves sharing a parent) or a reticulate edge of a reticulate cherry
(two leaves whose parents are connected by a reticulate edge). In this paper,
we present a fixed-parameter tractable algorithm to solve the problem of
finding a maximum agreement cherry-reduced subnetwork (MACRS) between two
rooted binary level-1 networks. This is first exact algorithm proposed to solve
the MACRS problem. As proven in earlier work, there is a direct relationship
between finding an MACRS and calculating a distance based on cherry operations.
As a result, the proposed algorithm also provides a distance that can be used
for the comparison of level-1 networks.


  The rapidly expanding market for regenerative medicines and cell therapies
highlights the need to advance the understanding of cellular metabolisms and
improve the prediction of cultivation production process for human induced
pluripotent stem cells (iPSCs). In this paper, a metabolic kinetic model was
developed to characterize underlying mechanisms of iPSC culture process, which
can predict cell response to environmental perturbation and support process
control. This model focuses on the central carbon metabolic network, including
glycolysis, pentose phosphate pathway (PPP), tricarboxylic acid (TCA) cycle,
and amino acid metabolism, which plays a crucial role to support iPSC
proliferation. Heterogeneous measures of extracellular metabolites and multiple
isotopic tracers collected under multiple conditions were used to learn
metabolic regulatory mechanisms. Systematic cross-validation confirmed the
model's performance in terms of providing reliable predictions on cellular
metabolism and culture process dynamics under various culture conditions. Thus,
the developed mechanistic kinetic model can support process control strategies
to strategically select optimal cell culture conditions at different times,
ensure cell product functionality, and facilitate large-scale manufacturing of
regenerative medicines and cell therapies.


  Sequence alignment is common nowadays as it is used in many fields to
determine how closely two sequences are related and at times to see how little
they differ. In computational biology / Bioinformatics, there are many
algorithms developed over the course of time to not only align two sequences
quickly but also get good laboratory results from these alignments. The first
algorithms developed were based of a technique called Dynamic Programming,
which were very slow but were optimal when it comes to sensitivity. To improve
speed, more algorithms today are based of heuristic approach, by sacrificing
sensitivity. In this paper, we are going to improve on a heuristic algorithm
called MASAA (Multiple Anchor Staged Local Sequence Alignment Algorithm) and
MASAA Sensitive which we published previously. This new algorithm appropriately
called Maximum Match Subsequence Alignment Algorithm Finely Grained. The
algorithm is based on suffix tree data structure like our previous algorithms,
but to improve sensitivity, we employ adaptive seeds, and finely grained
perfect match seeds in between the already identified anchors. We tested this
algorithm on a randomly generated sequences, and Rosetta dataset where the
sequence length ranged up to 500 thousand.


  Designing protein sequences with desired biological function is crucial in
biology and chemistry. Recent machine learning methods use a surrogate
sequence-function model to replace the expensive wet-lab validation. How can we
efficiently generate diverse and novel protein sequences with high fitness? In
this paper, we propose IsEM-Pro, an approach to generate protein sequences
towards a given fitness criterion. At its core, IsEM-Pro is a latent generative
model, augmented by combinatorial structure features from a separately learned
Markov random fields (MRFs). We develop an Monte Carlo Expectation-Maximization
method (MCEM) to learn the model. During inference, sampling from its latent
space enhances diversity while its MRFs features guide the exploration in high
fitness regions. Experiments on eight protein sequence design tasks show that
our IsEM-Pro outperforms the previous best methods by at least 55% on average
fitness score and generates more diverse and novel protein sequences.


  High-throughput sequencing (HTS) technologies have revolutionized the field
of genomics, enabling rapid and cost-effective genome analysis for various
applications. However, the increasing volume of genomic data generated by HTS
technologies presents significant challenges for computational techniques to
effectively analyze genomes. To address these challenges, several
algorithm-architecture co-design works have been proposed, targeting different
steps of the genome analysis pipeline. These works explore emerging
technologies to provide fast, accurate, and low-power genome analysis.
  This paper provides a brief review of the recent advancements in accelerating
genome analysis, covering the opportunities and challenges associated with the
acceleration of the key steps of the genome analysis pipeline. Our analysis
highlights the importance of integrating multiple steps of genome analysis
using suitable architectures to unlock significant performance improvements and
reduce data movement and energy consumption. We conclude by emphasizing the
need for novel strategies and techniques to address the growing demands of
genomic data generation and analysis.


  Microtubules self-organize to structure part of the cellular cytoskeleton. As
such they give cells their shape and play a crucial role in cell division and
intracellular transport. Past studies have identified diverse spatio-temporal
patterns into which microtubules can organize when driven by motor proteins.
The question remains if there is an appropriate way to quantify these
structures and gain new knowledge about the physical principles of
self-organization in microtubule-motor mixtures. Here, we aim to approach this
problem from a complexity science perspective. We introduce an entropy-based
measure to evaluate the structural complexity of spatial patterns emerging in a
simplified agent-based computational model of a microtubule-motor interactions.
Our results demonstrate that the proposed quantifier discriminates well between
ordered, disordered, and intermediate structures. Besides, our study indicates
that the transition to steady states in such a system is likely to be
discontinuous and exhibits distinct properties of self-organized criticality.


  We investigate the spatial dynamics of two disease epidemics reaching a
three-species cyclic model. Regardless of their species, all individuals are
susceptible to being infected with two different pathogens, which spread
through person-to-person contact. The occurrence of coinfection leads to a
synergistic increase in the risk of hosts dying due to complications from
either disease. Our stochastic simulations show that departed areas inhabited
by hosts of a single pathogen arise from random initial conditions. The
single-disease spatial domains are bordered by interfaces of coinfected hosts
whose dynamics are curvature-driven. Our findings show that the coarsening
dynamics of the interface network are controlled by the fluctuations of
coinfection waves invading the single-disease territories. As the coinfection
mortality grows, the dynamics of the interface network attain the scaling
regime. We discover that organisms' infection risk is maximised if the
coinfection increases the death due to disease in $30\%$, and minimised as the
network dynamics reach the scaling regime, with species populations being
maximum. Our conclusions may help ecologists understand the dynamics of
epidemics and their impact on the stability of ecosystems.


  Small and large scale pandemics are a natural phenomenon repeatably appearing
throughout history, causing ecological and biological shifts in ecosystems and
a wide range of their habitats. These pandemics usually start with a single
strain but shortly become multi-strain due to a mutation process of the
pathogen causing the epidemic. In this study, we propose a novel
eco-epidemiological model that captures multi-species prey-predator dynamics
with a multi-strain pandemic. The proposed model extends and combines the
Lotka-Volterra prey-predator model and the Susceptible-Infectious-Recovered
(SIR) epidemiological model. We investigate the ecosystem's sensitivity and
stability during such a multi-strain pandemic through extensive simulation
relying on both synthetic cases as well as two real-world configurations. Our
results are aligned with known ecological and epidemiological findings, thus
supporting the adequacy of the proposed model in realistically capturing the
complex eco-epidemiological properties of the multi-species multi-strain
pandemic dynamics.


  Short-term forecasts of infectious disease spread are a critical component in
risk evaluation and public health decision making. While different models for
short-term forecasting have been developed, open questions about their relative
performance remain. Here, we compare short-term probabilistic forecasts of
popular mechanistic models based on the renewal equation with forecasts of
statistical time series models. Our empirical comparison is based on data of
the daily incidence of COVID-19 across six large US states over the first
pandemic year. We find that, on average, probabilistic forecasts from
statistical time series models are overall at least as accurate as forecasts
from mechanistic models. Moreover, statistical time series models better
capture volatility. Our findings suggest that domain knowledge, which is
integrated into mechanistic models by making assumptions about disease
dynamics, does not improve short-term forecasts of disease incidence. We note,
however, that forecasting is often only one of many objectives and thus
mechanistic models remain important, for example, to model the impact of
vaccines or the emergence of new variants.


  The endoplasmic reticulum (ER), a cellular organelle that forms a
cell-spanning network of tubes and sheets, is an important location of protein
synthesis and folding. When the ER experiences sustained unfolded protein
stress, IRE1 proteins embedded in the ER membrane activate and assemble into
clusters as part of the unfolded protein response (UPR). We use kinetic Monte
Carlo simulations to explore IRE1 clustering dynamics on the surface of ER
tubes. While initially growing clusters are approximately round, once a cluster
is sufficiently large a shorter interface length can be achieved by `wrapping'
around the ER tube. A wrapped cluster can grow without further interface length
increases. Relative to wide tubes, narrower tubes enable cluster wrapping at
smaller cluster sizes. Our simulations show that wrapped clusters on narrower
tubes grow more rapidly, evaporate more slowly, and require a lower protein
concentration to grow compared to equal-area round clusters on wider tubes.
These results suggest that cluster wrapping, facilitated by narrower tubes,
could be an important factor in the growth and stability of IRE1 clusters and
thus impact the persistence of the UPR, connecting geometry to signaling
behavior. This work is consistent with recent experimental observations of IRE1
clusters wrapped around narrow tubes in the ER network.


  We introduce a minimal model of multilevel selection on structured
populations, considering the interplay between game theory and population
dynamics. Through a bottleneck process, finite groups are formed with
cooperators and defectors sampled from an infinite pool. After the
fragmentation, these transient compartments grow until the carrying capacity is
attained. Eventually, all compartments are merged, well mixed and the whole
process is repeated. We show that cooperators, even if interacting only through
mean-field intra-group interactions that favor defectors, may perform well
because of the inter-group competition and the size diversity among the
compartments. These cycles of isolation and coalescence may therefore be
important in maintaining diversity among different species or strategies and
may help to understand the underlying mechanisms of the scaffolding processes
in the transition to multicellularity.


  Increased proliferation is a key driver of tumorigenesis, and quantification
of mitotic activity is a standard task for prognostication. The goal of this
systematic review is scholarly analysis of all available references on mitotic
activity in feline tumors, and to provide an overview of the measuring methods
and prognostic value. A systematic literature search in PubMed and Scopus and a
manual search in Google Scholar was conducted. All articles on feline tumors
that correlated mitotic activity with patient outcome were identified. Data
analysis revealed that of the eligible 42 articles, the mitotic count (MC,
mitotic figures per tumor area) was evaluated in 39 instances and the mitotic
index (MI, mitotic figures per tumor cells) in three instances. The risk of
bias was considered high for most studies (26/42, 62%) based on small study
populations, insufficient details of the MC/MI methods, and lack of statistical
measures for diagnostic accuracy or effect on outcome. The MC/MI methods varied
markedly between studies. A significant association of the MC with survival was
determined in 21/29 (72%) studies, while one study found an inverse effect.
There were three tumor types with at least four studies and a prognostic
association was found in 5/6 studies on mast cell tumors, 5/5 on mammary tumors
and 3/4 on soft tissue sarcomas. The MI was shown to correlate with survival by
two research groups, however a comparison to the MC was not conducted. An
updated systematic review will be needed with of new literature for different
tumor types.


  Genome data are crucial in modern medicine, offering significant potential
for diagnosis and treatment. Thanks to technological advancements, many
millions of healthy and diseased genomes have already been sequenced; however,
obtaining the most suitable data for a specific study, and specifically for
validation studies, remains challenging with respect to scale and access.
Therefore, in silico genomics sequence generators have been proposed as a
possible solution. However, the current generators produce inferior data using
mostly shallow (stochastic) connections, detected with limited computational
complexity in the training data. This means they do not take the appropriate
biological relations and constraints, that originally caused the observed
connections, into consideration. To address this issue, we propose
cancer-inspired genomics mapper model (CGMM), that combines genetic algorithm
(GA) and deep learning (DL) methods to tackle this challenge. CGMM mimics
processes that generate genetic variations and mutations to transform readily
available control genomes into genomes with the desired phenotypes. We
demonstrate that CGMM can generate synthetic genomes of selected phenotypes
such as ancestry and cancer that are indistinguishable from real genomes of
such phenotypes, based on unsupervised clustering. Our results show that CGMM
outperforms four current state-of-the-art genomics generators on two different
tasks, suggesting that CGMM will be suitable for a wide range of purposes in
genomic medicine, especially for much-needed validation studies.


  Molecular relational learning, whose goal is to learn the interaction
behavior between molecular pairs, got a surge of interest in molecular sciences
due to its wide range of applications. Recently, graph neural networks have
recently shown great success in molecular relational learning by modeling a
molecule as a graph structure, and considering atom-level interactions between
two molecules. Despite their success, existing molecular relational learning
methods tend to overlook the nature of chemistry, i.e., a chemical compound is
composed of multiple substructures such as functional groups that cause
distinctive chemical reactions. In this work, we propose a novel relational
learning framework, called CGIB, that predicts the interaction behavior between
a pair of graphs by detecting core subgraphs therein. The main idea is, given a
pair of graphs, to find a subgraph from a graph that contains the minimal
sufficient information regarding the task at hand conditioned on the paired
graph based on the principle of conditional graph information bottleneck. We
argue that our proposed method mimics the nature of chemical reactions, i.e.,
the core substructure of a molecule varies depending on which other molecule it
interacts with. Extensive experiments on various tasks with real-world datasets
demonstrate the superiority of CGIB over state-of-the-art baselines. Our code
is available at https://github.com/Namkyeong/CGIB.


  Based on the traditional VAE, a novel neural network model is presented, with
the latest molecular representation, SELFIES, to improve the effect of
generating new molecules. In this model, multi-layer convolutional network and
Fisher information are added to the original encoding layer to learn the data
characteristics and guide the encoding process, which makes the features of the
data hiding layer more aggregated, and integrates the Long Short Term Memory
neural network (LSTM) into the decoding layer for better data generation, which
effectively solves the degradation phenomenon generated by the encoding layer
and decoding layer of the original VAE model. Through experiments on zinc
molecular data sets, it is found that the similarity in the new VAE is 8.47%
higher than that of the original ones. SELFIES are better at generating a
variety of molecules than the traditional molecular representation, SELFIES.
Experiments have shown that using SELFIES and the new VAE model presented in
this paper can improve the effectiveness of generating new molecules.


  Understanding waning of vaccine-induced protection is important for both
immunology and public health. Population heterogeneities in underlying
(pre-vaccination) susceptibility and vaccine response can cause measured
vaccine effectiveness (mVE) to change over time even in the absence of pathogen
evolution and any actual waning of immune responses. We use a multi-scale
agent-based models parameterized using epidemiological and immunological data,
to investigate the effect of these heterogeneities on mVE as measured by the
hazard ratio. Based on our previous work, we consider waning of antibodies
according to a power law and link it to protection in two ways: 1) motivated by
correlates of risk data and 2) using a within-host model of stochastic viral
extinction. The effect of the heterogeneities is given by concise and
understandable formulas, one of which is essentially a generalization of
Fisher's fundamental theorem of natural selection to include higher
derivatives. Heterogeneity in underlying susceptibility accelerates apparent
waning, whereas heterogeneity in vaccine response slows down apparent waning.
Our models suggest that heterogeneity in underlying susceptibility is likely to
dominate. However, heterogeneity in vaccine response offsets <10% to >100%
(median of 29%) of this effect in our simulations. Our methodology and results
may be helpful in understanding competing heterogeneities and waning of
immunity and vaccine-induced protection. Our study suggests heterogeneity is
more likely to 'bias' mVE downwards towards faster waning of immunity but a
subtle bias in the opposite direction is also plausible.


  In this paper, we examine the computational complexity of enumeration in
certain genome rearrangement models. We first show that the Pairwise
Rearrangement problem in the Single Cut-and-Join model (Bergeron, Medvedev, &
Stoye, J. Comput. Biol. 2010) is $\#\textsf{P}$-complete under polynomial-time
Turing reductions. Next, we show that in the Single Cut or Join model (Feijao &
Meidanis, IEEE ACM Trans. Comp. Biol. Bioinf. 2011), the problem of enumerating
all medians ($\#$Median) is logspace-computable ($\textsf{FL}$), improving upon
the previous polynomial-time ($\textsf{FP}$) bound of Mikl\'os & Smith (RECOMB
2015).


  Recent advancements in specialized large-scale architectures for training
image and language have profoundly impacted the field of computer vision and
natural language processing (NLP). Language models, such as the recent ChatGPT
and GPT4 have demonstrated exceptional capabilities in processing, translating,
and generating human languages. These breakthroughs have also been reflected in
protein research, leading to the rapid development of numerous new methods in a
short time, with unprecedented performance. Language models, in particular,
have seen widespread use in protein research, as they have been utilized to
embed proteins, generate novel ones, and predict tertiary structures. In this
book chapter, we provide an overview of the use of protein generative models,
reviewing 1) language models for the design of novel artificial proteins, 2)
works that use non-Transformer architectures, and 3) applications in directed
evolution approaches.


  In this paper, we build a mathematical model for the dynamics of COVID-19 to
assess the impact of placing healthy individuals in quarantine and isolating
infected ones on the number of hospitalization and intensive care unit cases.
The proposed model is fully analyzed in order to prove the positivity of
solutions, to study the local and global stability of the disease-free
equilibria and to drive the basic and control reproduction numbers of the
model. Oman COVID-19 data is used to calibrate the model and estimate the
parameters. In particular, the published data for the year 2020 is used, when
two waves of the disease hit the country. Moreover, this period of time is
chosen when no vaccine had been introduced, but only the non-pharmaceutical
intervention (NPI) strategies were the only effective methods to control the
spread and, consequently, control the hospitalization cases to avoid pressuring
the health system. Based on the estimated parameters, the reproduction number
and contribution of different transmission routes are approximated numerically.
Sensitivity analysis is performed to identify the significant parameters in
spreading the disease. Numerical simulation is carried out to demonstrate the
effects of quarantine and isolation on the number of hospitalized cases.


  Minimal trap spaces (MTSs) capture subspaces in which the Boolean dynamics is
trapped, whatever the update mode. They correspond to the attractors of the
most permissive mode. Due to their versatility, the computation of MTSs has
recently gained traction, essentially by focusing on their enumeration. In this
paper, we address the logical reasoning on universal properties of MTSs in the
scope of two problems: the reprogramming of Boolean networks for identifying
the permanent freeze of Boolean variables that enforce a given property on all
the MTSs, and the synthesis of Boolean networks from universal properties on
their MTSs. Both problems reduce to solving the satisfiability of quantified
propositional logic formula with 3 levels of quantifiers
($\exists\forall\exists$). In this paper, we introduce a Counter-Example Guided
Refinement Abstraction (CEGAR) to efficiently solve these problems by coupling
the resolution of two simpler formulas. We provide a prototype relying on
Answer-Set Programming for each formula and show its tractability on a wide
range of Boolean models of biological networks.


  Fitness functions map large combinatorial spaces of biological sequences to
properties of interest. Inferring these multimodal functions from experimental
data is a central task in modern protein engineering. Global epistasis models
are an effective and physically-grounded class of models for estimating fitness
functions from observed data. These models assume that a sparse latent function
is transformed by a monotonic nonlinearity to emit measurable fitness. Here we
demonstrate that minimizing contrastive loss functions, such as the
Bradley-Terry loss, is a simple and flexible technique for extracting the
sparse latent function implied by global epistasis. We argue by way of a
fitness-epistasis uncertainty principle that the nonlinearities in global
epistasis models can produce observed fitness functions that do not admit
sparse representations, and thus may be inefficient to learn from observations
when using a Mean Squared Error (MSE) loss (a common practice). We show that
contrastive losses are able to accurately estimate a ranking function from
limited data even in regimes where MSE is ineffective. We validate the
practical utility of this insight by showing contrastive loss functions result
in consistently improved performance on benchmark tasks.


  Biomolecular condensates constitute a newly recognized form of spatial
organization in living cells. Although many condensates are believed to form as
a result of phase separation, the physicochemical properties that determine the
phase behavior of heterogeneous biomolecular mixtures are only beginning to be
explored. Theory and simulation provide invaluable tools for probing the
relationship between molecular determinants, such as protein and RNA sequences,
and the emergence of phase-separated condensates in such complex environments.
This review covers recent advances in the prediction and computational design
of biomolecular mixtures that phase-separate into many coexisting phases.
First, we review efforts to understand the phase behavior of mixtures with
hundreds or thousands of species using theoretical models and statistical
approaches. We then describe progress in developing analytical theories and
coarse-grained simulation models to predict multiphase condensates with the
molecular detail required to make contact with biophysical experiments. We
conclude by summarizing the challenges ahead for modeling the inhomogeneous
spatial organization of biomolecular mixtures in living cells.


  We investigate a tritrophic system in which organisms' energy depletion,
resulting from failed selection attempts, leads to a partial loss of capacity
to win the cyclic spatial game. The energy required to maintain optimal
organism fitness may be impacted by changes in environmental conditions,
increasing the death risk due to accelerated deterioration of health. We
simulate the evolutionary behavioural strategy performed by individuals of one
species, which consists in balancing efforts dedicated to reproduction and
mobility to minimise the chances of death by lack of energy. We show that the
unevenness introduced by the trade-off strategy unbalances the
rock-paper-scissors model, with the predominant species profiting from enemies'
lower birth rate of enemies. Quantifying the spatial patterns, we demonstrate
that the characteristic length scale of single-species domains decreases as
energy loss accelerates due to environmental changes. The erosion in the
spatial patterns provoked by the reproduction-mobility trade-off benefits
biodiversity, with coexistence probability rising for faster energy depletion
and higher trade-off factors. The findings have implications for ecologists
seeking to understand the impact of survival behaviour on biodiversity
promotion.


  Many major questions in the theory of evolutionary dynamics can in a
meaningful sense be mapped to analyses of stochastic trajectories in game
theoretic contexts. Often the approach is to analyze small numbers of distinct
populations and/or to assume dynamics occur within a regime of population sizes
large enough that deterministic trajectories are an excellent approximation of
reality. The addition of ecological factors, termed "eco-evolutionary
dynamics", further complicates the dynamics and results in many problems which
are intractable or impractically messy for current theoretical methods.
However, an analogous but underexplored approach is to analyze these systems
with an eye primarily towards uncertainty in the models themselves. In the
language of researchers in Reinforcement Learning and adjacent fields, a
Partially Observable Markov Process. Here we introduce a duality which maps the
complexity of accounting for both ecology and individual genotypic/phenotypic
types onto a problem of accounting solely for underlying information-theoretic
computations rather than drawing physical boundaries which do not change the
computations. Armed with this equivalence between computation and the relevant
biophysics, which we term Taak-duality, we attack the problem of "directed
evolution" in the form of a Partially Observable Markov Decision Process. This
provides a tractable case of studying eco-evolutionary trajectories of a highly
general type, and of analyzing questions of potential limits on the efficiency
of evolution in the directed case.


  The spatio-temporal organization of proteins within the cell membrane can
affect numerous biological functions, including cell signaling, communication,
and transportation. Deviations from normal spatial arrangements have been
observed in various diseases, and better understanding this process is a key
stepping-stone to advancing development of clinical interventions. However,
given the nanometer length scales involved, detecting these subtle changes has
primarily relied on complex super resolution and single molecule imaging
methods. In this work, we demonstrate an alternative fluorescent imaging
strategy for detecting protein organization based on a material that exhibits a
unique photophysical behavior known as aggregation induced emission (AIE).
Organic AIE molecules have an increase in emission signal when they are in
close proximity and the molecular motion is restricted. This property
simultaneously addresses the high background noise and low detection signal
that limit conventional widefield fluorescent imaging. To demonstrate the
potential of this approach, the fluorescent molecule sensor is conjugated to a
human epidermal growth factor receptor 2 (HER2) specific antibody and used to
investigate the spatio-temporal behavior of HER2 clustering in the membrane of
HER2-overexpressing breast cancer cells. Notably, the disruption of HER2
clusters in response to an FDA-approved monoclonal antibody therapeutic
(Trastuzumab) is successfully detected using a simple widefield fluorescent
microscope. While the sensor demonstrated here is optimized for sensing HER2
clustering, it is an easily adaptable platform. Moreover, given the
compatibility with widefield imaging, the system has the potential to be used
with high-throughput imaging techniques, accelerating investigations into
membrane protein spatio-temporal organization.


  Polygenic risk scores (PRS) summarize the combined effect of common risk
variants and are associated with breast cancer risk in patients without
identifiable monogenic risk factors. One of the most well-validated PRSs in
breast cancer to date is PRS313, which was developed from a Northern European
biobank but has shown attenuated performance in non-European ancestries. We
further investigate the generalizability of the PRS313 for American women of
European (EA), African (AFR), Asian (EAA), and Latinx (HL) ancestry within one
institution with a singular EHR system, genotyping platform, and quality
control process. We found that the PRS313 achieved overlapping Areas under the
ROC Curve (AUCs) in females of Lantix (AUC, 0.68; 95 CI, 0.65-0.71) and
European ancestry (AUC, 0.70; 95 CI, 0.69-0.71) but lower AUCs for the AFR and
EAA populations (AFR: AUC, 0.61; 95 CI, 0.56-0.65; EAA: AUC, 0.64; 95 CI,
0.60-0.680). While PRS313 is associated with Hormone Positive (HR+) disease in
European Americans (OR, 1.42; 95 CI, 1.16-1.64), for Latinx females, it may be
instead associated with Human Epidermal Growth Factor Receptor 2 (HER2+)
disease (OR, 2.52; 95 CI, 1.35-4.70) although due to small numbers, additional
studies are needed. In summary, we found that PRS313 was significantly
associated with breast cancer but with attenuated accuracy in women of African
and Asian descent within a singular health system in Los Angeles. Our work
further highlights the need for additional validation in diverse cohorts prior
to clinical implementation of polygenic risk scores.


  RNA structure and functional dynamics play fundamental roles in controlling
biological systems. Molecular dynamics simulation, which can characterize
interactions at an atomistic level, can advance the understanding on new drug
discovery, manufacturing, and delivery mechanisms. However, it is
computationally unattainable to support the development of a digital twin for
enzymatic reaction network mechanism learning, and end-to-end bioprocess design
and control. Thus, we create a hybrid ("mechanistic + machine learning") model
characterizing the interdependence of RNA structure and functional dynamics
from atomistic to macroscopic levels. To assess the proposed modeling strategy,
in this paper, we consider RNA degradation which is a critical process in
cellular biology that affects gene expression. The empirical study on RNA
lifetime prediction demonstrates the promising performance of the proposed
multi-scale bioprocess hybrid modeling strategy.


  Proteins often regulate their activities via allostery - or action at a
distance - in which the binding of a ligand at one binding site influences the
affinity for another ligand at a distal site. Although less studied than in
proteins, allosteric effects have been observed in experiments with DNA as
well. In these experiments two or more proteins bind at distinct DNA sites and
interact indirectly with each other, via a mechanism mediated by the linker DNA
molecule. We develop a mechanical model of DNA/protein interactions which
predicts three distinct mechanisms of allostery. Two of these involve an
enthalpy-mediated allostery, while a third mechanism is entropy driven. We
analyze experiments of DNA allostery and highlight the distinctive signatures
allowing one to identify which of the proposed mechanisms best fits the data.


  In this paper, we analyse a recently proposed predator-prey model with ratio
dependence and Holling type III functional response, with particular emphasis
on the dynamics close to extinction. By using Briot-Bouquet transformation we
transform the model into a system, where the extinction steady state is
represented by up to three distinct steady states, whose existence is
determined by the values of appropriate Lambert W functions. We investigate how
stability of extinction and coexistence steady states is affected by the rate
of predation, predator fecundity, and the parameter characterising the strength
of functional response. The results suggest that the extinction steady state
can be stable for sufficiently high predation rate and for sufficiently small
predator fecundity. Moreover, in certain parameter regimes, a stable extinction
steady state can coexist with a stable prey-only equilibrium or with a stable
coexistence equilibrium, and it is rather the initial conditions that determine
whether prey and predator populations will be maintained at some steady level,
or both of them will become extinct. Another possibility is for coexistence
steady state to be unstable, in which case sustained periodic oscillations
around it are observed. Numerical simulations are performed to illustrate the
behaviour for all dynamical regimes, and in each case a corresponding phase
plane of the transformed system is presented to show a correspondence with
stable and unstable extinction steady state.


  Proteins are complex biomolecules that perform a variety of crucial functions
within living organisms. Designing and generating novel proteins can pave the
way for many future synthetic biology applications, including drug discovery.
However, it remains a challenging computational task due to the large modeling
space of protein structures. In this study, we propose a latent diffusion model
that can reduce the complexity of protein modeling while flexibly capturing the
distribution of natural protein structures in a condensed latent space.
Specifically, we propose an equivariant protein autoencoder that embeds
proteins into a latent space and then uses an equivariant diffusion model to
learn the distribution of the latent protein representations. Experimental
results demonstrate that our method can effectively generate novel protein
backbone structures with high designability and efficiency. The code will be
made publicly available at
https://github.com/divelab/AIRS/tree/main/OpenProt/LatentDiff


  Given a gene-tree labeled topology $G$ and a species tree $S$, the "ancestral
configurations" at an internal node $k$ of $S$ represent the combinatorially
different sets of gene lineages that can be present at $k$ when all possible
realizations of $G$ in $S$ are considered. Ancestral configurations have been
introduced as a data structure for evaluating the conditional probability of a
gene-tree labeled topology given a species tree, and their enumeration assists
in describing the complexity of this computation. In the case that the
gene-tree labeled topology $G=t$ matches that of the species tree $S$, by
techniques of analytic combinatorics, we study distributional properties of the
"total" number of ancestral configurations measured across the different nodes
of a random labeled topology $t$ selected under the uniform and the Yule
probability models. Under both of these probabilistic scenarios, we show that
the total number $T_n$ of ancestral configurations of a random labeled topology
of $n$ taxa asymptotically follows a lognormal distribution. Over uniformly
distributed labeled topologies, the asymptotic growth of the mean and the
variance of $T_n$ are found to satisfy $\mathbb{E}_{\rm U}[T_n] \sim 2.449
\cdot 1.333^n$ and $\mathbb{V}_{\rm U}[T_n] \sim 5.050 \cdot 1.822^n$,
respectively. Under the Yule model, which assigns higher probabilities to more
balanced labeled topologies, we obtain the mean $\mathbb{E}_{\rm Y}[T_n] \sim
1.425^n$ and the variance $\mathbb{V}_{\rm Y}[T_n] \sim 2.045^n$.


  Color naming in natural languages is not arbitrary: it reflects efficient
partitions of perceptual color space modulated by the relative needs to
communicate about different colors. These psychophysical and communicative
constraints help explain why languages around the world have remarkably
similar, but not identical, mappings of colors to color terms. Languages
converge on a small set of efficient representations. But languages also
evolve, and the number of terms in a color vocabulary may change over time.
Here we show that history, i.e. the existence of an antecedent color
vocabulary, acts as a non-adaptive constraint that biases the choice of
efficient solution as a language transitions from a vocabulary of size n to n+1
terms. Moreover, as vocabularies evolve to include more terms they explore a
smaller fraction of all possible efficient vocabularies compared to
equally-sized vocabularies constructed de novo. This path dependence on the
cultural evolution of color naming presents an opportunity. Historical
constraints can be used to reconstruct ancestral color vocabularies, allowing
us to answer long-standing questions about the evolutionary sequences of color
words, and enabling us to draw inferences from phylogenetic patterns of
language change.


  In real world, mutations of genetic sequences are often accompanied by their
recombinations. Such joint phenomena are modeled by phylogenetic networks.
Nakkleh formulated the phylogenetic network reconstruction problem (PNRP) as
follows: Given a family of phylogenetic trees over a common set of taxa, is
there a unique minimal phylogenetic network whose set of spanning trees
contains the family? There are different answers to PNRP, since there are
different ways to define what a minimal network is (based on different
optimization criteria). Inspired by ideas from topological data analysis (TDA),
we devise lattice-diagram models for the visualization of phylogenetic networks
and of filtrations, called the cliquegram and the facegram, respectively, both
generalizing the dendrogram model of phylogenetic trees. Both models allow us
to solve the PNRP in a rigorous way and free of choosing optimization criteria.
The solution to the phylogenetic network and filtration reconstruction process
is obtained by taking the join operation of the dendrograms on the lattice of
cliquegrams, and of facegrams, respectively. Furthermore, we show that
computing the join-facegram from a given set of dendrograms is polynomial in
the size and number of the input trees. We propose two novel invariants of
facegrams, (i) the face-Reeb graph and (ii) the mergegram of a facegram. We
show the mergegram is 1-Lipschitz stable, while the face-Reeb graph is not. In
particular, we show that the mergegram is invariant of weak equivalences of
filtrations (a stronger form of homotopy equivalence). This new TDA-signature,
the mergegram, can be used as a computable proxy for phylogenetic networks and
also, more broadly, for filtrations of datasets, which might be of independent
interest to TDA. To illustrate the utility of those new TDA-tools to
phylogenetics, we provide experiments with artificial and benchmark biological
data.


  Machine learning (ML) approaches have been used to develop highly accurate
and efficient applications in many fields including bio-medical science.
However, even with advanced ML techniques, cancer classification using gene
expression data is still complicated because of the high dimensionality of the
datasets employed. We developed a new fuzzy gene selection technique (FGS) to
identify informative genes to facilitate cancer classification and reduce the
dimensionality of the available gene expression data. Three feature selection
methods (Mutual Information, F-ClassIf, and Chi-squared) were evaluated and
employed to obtain the score and rank for each gene. Then, using Fuzzification
and Defuzzification methods to obtain the best single score for each gene,
which aids in the identification of significant genes. Our study applied the
fuzzy measures to six gene expression datasets including four Microarray and
two RNA-seq datasets for evaluating the proposed algorithm. With our
FGS-enhanced method, the cancer classification model achieved 96.5%,96.2%,96%,
and 95.9% for accuracy, precision, recall, and f1-score respectively, which is
significantly higher than 69.2% accuracy, 57.8% precision, 66% recall, and
58.2% f1-score when the standard MLP method was used. In examining the six
datasets that were used, the proposed model demonstrates it's capacity to
classify cancer effectively.


  We report a flexible language-model based deep learning strategy, applied
here to solve complex forward and inverse problems in protein modeling, based
on an attention neural network that integrates transformer and graph
convolutional architectures in a causal multi-headed graph mechanism, to
realize a generative pretrained model. The model is applied to predict
secondary structure content (per-residue level and overall content), protein
solubility, and sequencing tasks. Further trained on inverse tasks, the model
is rendered capable of designing proteins with these properties as target
features. The model is formulated as a general framework, completely
prompt-based, and can be adapted for a variety of downstream tasks. We find
that adding additional tasks yields emergent synergies that the model exploits
in improving overall performance, beyond what would be possible by training a
model on each dataset alone. Case studies are presented to validate the method,
yielding protein designs specifically focused on structural proteins, but also
exploring the applicability in the design of soluble, antimicrobial
biomaterials. While our model is trained to ultimately perform 8 distinct
tasks, with available datasets it can be extended to solve additional problems.
In a broader sense, this work illustrates a form of multiscale modeling that
relates a set of ultimate building blocks (here, byte-level utf8 characters
that define the nature of the physical system at hand) to complex output. This
materiomic scheme captures complex emergent relationships between universal
building block and resulting properties via a synergizing learning capacity to
express a set of potentialities embedded in the knowledge used in training, via
the interplay of universality and diversity.