  Protoplasts represent a unique experimental system for the circulation and
formation of genetically modified plants. Here, protoplasts were isolated from
genetically modified hairy root tissues of Brassica oleracea L. induced by the
Agrobacterium rhizogenes strain (ATCC13332). The concentration of enzyme
solutions utilized for protoplast isolation was 1.5 % Cellulase YC and 0.1 %
Pectolyase Y23 in 13% mannitol solution, which resulted in high efficiency of
isolation within 8 hours, in which the protoplast yield was 2 x 104 cells ml-1
and the percentage of viability was 72%. Each protoplast has one nucleus with a
nucleation of 48%. A polymerase chain reaction (PCR) assay verified the
presence of rol B and rol C genes in hairy root tissues by detaching a single
bundle of DNA replication from these roots using a specific pair of primers.
The current study demonstrated that A. rhizogenes strain (ATCC13332) is a
vector for the incorporation of T-DNA genes into cauliflower plants, as well as
the success of the hairy roots retention of rol B and rol C genes transferred
to it.


  In this article, we will see a new approach to study the impact of a small
microscopic population of cancer cells on a macroscopic population of healthy
cells, with an example inspired by pathological hematopoiesis. Hematopoiesis is
the biological phenomenon of blood cells production by differentiation of cells
called hematopoietic stem cells (HSCs). We will study the dynamics of a
stochastic $4$-dimensional process describing the evolution over time of the
number of healthy and cancer stem cells and the number of healthy and mutant
red blood cells. The model takes into account the amplification between stem
cells and red blood cells as well as the regulation of this amplification as a
function of the number of red blood cells (healthy and mutant). A single cancer
HSC is considered while other populations are in large numbers. We assume that
the unique cancer HSC randomly switches between an active and a quiescent
state. We show the convergence in law of this process towards a piecewise
deterministic Markov process (PDMP), when the population size goes to infinity.
We then study the long time behaviour of this limit process. We show the
existence and uniqueness of an absolutely continuous invariant probability
measure with respect to the Lebesgue's measure for the limit PDMP, previously
gathered. We describe the support of the invariant probability and show that
the process converges in total variation towards it, using theory develop by M.
Benaim et al. We finally identify the invariant probability using its
infinitesimal generator. Thanks to this probabilistic approach, we obtain a
stationary system of partial differential equation describing the impact of
cancer HSC quiescent phases and regulation on the cell density of the
hematopoietic system studied.


  Discovering causal effects is at the core of scientific investigation but
remains challenging when only observational data is available. In practice,
causal networks are difficult to learn and interpret, and limited to relatively
small datasets. We report a more reliable and scalable causal discovery method
(iMIIC), based on a general mutual information supremum principle, which
greatly improves the precision of inferred causal relations while
distinguishing genuine causes from putative and latent causal effects. We
showcase iMIIC on synthetic and real-life healthcare data from 396,179 breast
cancer patients from the US Surveillance, Epidemiology, and End Results
program. More than 90\% of predicted causal effects appear correct, while the
remaining unexpected direct and indirect causal effects can be interpreted in
terms of diagnostic procedures, therapeutic timing, patient preference or
socio-economic disparity. iMIIC's unique capabilities open up new avenues to
discover reliable and interpretable causal networks across a range of research
fields.


  Cell-free protein synthesis (CFPS) systems are an attractive to complement
the usual cell-based synthesis of proteins, especially for screening
approaches. The literature describes a wide variety of CFPS systems, but their
performance is difficult to compare since the reaction components are often
used at different concentrations. Therefore, we have developed a calculation
tool based on amino acid balancing to evaluate the performance of CFPS by
determining the fractional yield as the ratio between theoretically achievable
and achieved protein molar concentration. This tool was applied to a series of
experiments from our lab and to various systems described in the literature to
identify systems that synthesize proteins very efficiently and those that still
have potential for higher yields. The well-established Escherichia coli system
showed a high efficiency in the utilization of amino acids, but interestingly,
less-attention-paid systems, such as the one based on Streptomyces lividans,
also demonstrated exceptional fractional yields of 100%, meaning complete amino
acid conversion. The methods and tools described here can quickly identify when
a system has reached its maximum or has other limitations. We believe that the
approach described here will facilitate the evaluation and optimization of
existing CFPS systems and provide the basis for the systematic development of
new CFPS systems.


  We propose a reinforcement-learning algorithm to tackle the challenge of
reconstructing phylogenetic trees. The search for the tree that best describes
the data is algorithmically challenging, thus all current algorithms for
phylogeny reconstruction use various heuristics to make it feasible. In this
study, we demonstrate that reinforcement learning can be used to learn an
optimal search strategy, thus providing a novel paradigm for predicting the
maximum-likelihood tree. Our proposed method does not require likelihood
calculation with every step, nor is it limited to greedy uphill moves in the
likelihood space. We demonstrate the use of the developed deep-Q-learning agent
on a set of unseen empirical data, namely, on unseen environments defined by
nucleotide alignments of up to 20 sequences. Our results show that the
likelihood scores of the inferred phylogenies are similar to those obtained
from widely-used software. It thus establishes a proof-of-concept that it is
beneficial to optimize a sequence of moves in the search-space, rather than
optimizing the progress made in every single move only. This suggests that a
reinforcement-learning based method provides a promising direction for
phylogenetic reconstruction.


  Contrastive learning have been widely used as pretext tasks for
self-supervised pre-trained molecular representation learning models in
AI-aided drug design and discovery. However, exiting methods that generate
molecular views by noise-adding operations for contrastive learning may face
the semantic inconsistency problem, which leads to false positive pairs and
consequently poor prediction performance. To address this problem, in this
paper we first propose a semantic-invariant view generation method by properly
breaking molecular graphs into fragment pairs. Then, we develop a
Fragment-based Semantic-Invariant Contrastive Learning (FraSICL) model based on
this view generation method for molecular property prediction. The FraSICL
model consists of two branches to generate representations of views for
contrastive learning, meanwhile a multi-view fusion and an auxiliary similarity
loss are introduced to make better use of the information contained in
different fragment-pair views. Extensive experiments on various benchmark
datasets show that with the least number of pre-training samples, FraSICL can
achieve state-of-the-art performance, compared with major existing counterpart
models.


  Generators of space-time dynamics in bioimaging have become essential to
build ground truth datasets for image processing algorithm evaluation such as
biomolecule detectors and trackers, as well as to generate training datasets
for deep learning algorithms. In this contribution, we leverage a stochastic
model, called birth-death-move (BDM) point process, in order to generate joint
dynamics of biomolecules in cells. This approach is very flexible and allows us
to model a system of particles in motion, possibly in interaction, that can
each possibly switch from a motion regime (e.g. Brownian) to another (e.g. a
directed motion), along with the appearance over time of new trajectories and
their death after some lifetime, all of these features possibly depending on
the current spatial configuration of all existing particles. We explain how to
specify all characteristics of a BDM model, with many practical examples that
are relevant for bioimaging applications. Based on real fluorescence microscopy
datasets, we finally calibrate our model to mimic the joint dynamics of
Langerin and Rab11 proteins near the plasma membrane. We show that the
resulting synthetic sequences exhibit comparable features as those observed in
real microscopy image sequences.


  Chemical reactions are the fundamental building blocks of drug design and
organic chemistry research. In recent years, there has been a growing need for
a large-scale deep-learning framework that can efficiently capture the basic
rules of chemical reactions. In this paper, we have proposed a unified
framework that addresses both the reaction representation learning and molecule
generation tasks, which allows for a more holistic approach. Inspired by the
organic chemistry mechanism, we develop a novel pretraining framework that
enables us to incorporate inductive biases into the model. Our framework
achieves state-of-the-art results on challenging downstream tasks. By
possessing chemical knowledge, our generative framework overcome the
limitations of current molecule generation models that rely on a small number
of reaction templates. In the extensive experiments, our model generates
synthesizable drug-like structures of high quality. Overall, our work presents
a significant step toward a large-scale deep-learning framework for a variety
of reaction-based applications.


  The continuing advances of omic technologies mean that it is now more
tangible to measure the numerous features collectively reflecting the molecular
properties of a sample. When multiple omic methods are used, statistical and
computational approaches can exploit these large, connected profiles.
Multi-omics is the integration of different omic data sources from the same
biological sample. In this review, we focus on correlation-based dimension
reduction approaches for single omic datasets, followed by methods for pairs of
omics datasets, before detailing further techniques for three or more omic
datasets. We also briefly detail network methods when three or more omic
datasets are available and which complement correlation-oriented tools. To aid
readers new to this area, these are all linked to relevant R packages that can
implement these procedures. Finally, we discuss scenarios of experimental
design and present road maps that simplify the selection of appropriate
analysis methods. This review will guide researchers navigate the emerging
methods for multi-omics and help them integrate diverse omic datasets
appropriately and embrace the opportunity of population multi-omics.


  A cornerstone of the classical view of tolerance is the elimination of
self-reactive T cells during negative selection in the thymus. However,
high-throughput T-cell receptor sequencing data has so far failed to detect
substantial signatures of negative selection in the observed repertoires. In
addition, quantitative estimates as well as recent experiments suggest that the
elimination of self-reactive T cells is at best incomplete. We discuss several
recent theoretical ideas that can explain tolerance while being consistent with
these observations, including collective decision making through quorum
sensing, and sensitivity to change through dynamic tuning and adaptation. We
propose that a unified quantitative theory of tolerance should combine these
elements to explain the plasticity of the immune system and its robustness to
autoimmunity.


  In freshwater ecosystems, aquatic insects that ontogenetically shift their
habitat from aquatic to terrestrial play vital roles as prey subsidies that
move nutrients and energy from aquatic to terrestrial food webs. As a result,
these subsidies negatively affect alternative terrestrial prey by enhancing
predator density. However, these aquatic insects can also transport
contamination to the terrestrial community that is primarily produced in
aquatic ecosystems. Which can reduce insectivore abundance and biomass, lower
insectivore reproductive success, and increase predation pressure on
alternative prey with consequences for aquatic and terrestrial food webs.
Motivated by this, here we consider a prey-predator model where the predator
consumes contaminated and uncontaminated prey together. We find that, at a high
level of contamination, the vulnerability of contaminated prey and predator is
determined by predation preference. More specifically, a very low predation
preference towards contaminated prey ensures predator persistence, whereas a
low, high preference excludes the predator from the system. Interestingly,
either contaminated prey or the predator exist at intermediate predation
preference due to bi-stability. Furthermore, when there is no contamination in
one of the prey, the other prey can not coexist due to apparent competition for
a specific range of predation preferences. However, when sufficient
contamination exists in one prey, alternative uncontaminated prey coexists.
With this, contamination also stabilizes and destabilizes the three species
dynamics. Our result also indicates that if the intensity of the contamination
in predator reproduction is low, then contaminated prey is more susceptible to
the contamination.


  A model of genomic sequence evolution on a species tree should include not
only a sequence substitution process, but also a coalescent process, since
different sites may evolve on different gene trees due to incomplete lineage
sorting. Chifman and Kubatko initiated the study of such models, leading to the
development of the SVDquartets methods of species tree inference. A key
observation was that symmetries in an ultrametric species tree led to
symmetries in the joint distribution of bases at the taxa. In this work, we
explore the implications of such symmetry more fully, defining new models
incorporating only the symmetries of this distribution, regardless of the
mechanism that might have produced them. The models are thus supermodels of
many standard ones with mechanistic parameterizations. We study phylogenetic
invariants for the models, and establish identifiability of species tree
topologies using them.


  Target Identification by Enzymes (TIE) problem aims to identify the set of
enzymes in a given metabolic network, such that their inhibition eliminates a
given set of target compounds associated with a disease while incurring minimum
damage to the rest of the compounds. This is an NP-complete problem, and thus
optimal solutions using classical computers fail to scale to large metabolic
networks. In this paper, we consider the TIE problem for identifying drug
targets in metabolic networks. We develop the first quantum optimization
solution, called QuTIE (Quantum optimization for Target Identification by
Enzymes), to this NP-complete problem. We do that by developing an equivalent
formulation of the TIE problem in Quadratic Unconstrained Binary Optimization
(QUBO) form, then mapping it to a logical graph, which is then embedded on a
hardware graph on a quantum computer. Our experimental results on 27 metabolic
networks from Escherichia coli, Homo sapiens, and Mus musculus show that QuTIE
yields solutions which are optimal or almost optimal. Our experiments also
demonstrate that QuTIE can successfully identify enzyme targets already
verified in wet-lab experiments for 14 major disease classes.


  We analyze plankton-nutrient food chain models composed of phytoplankton,
herbivorous zooplankton and a limiting nutrient. These models have played a key
role in understanding the dynamics of plankton in the oceanic layer. Given the
strong environmental and seasonal fluctuations that are present in the oceanic
layer, we propose a stochastic model for which we are able to fully classify
the longterm behavior of the dynamics. In order to achieve this we had to
develop new analytical techniques, as the system does not satisfy the regular
dissipativity conditions and the analysis is more subtle than in other
population dynamics models.


  We analyze the advantage of cell size control strategies in growing
populations under mortality constraints. We demonstrate a general advantage of
the adder control strategy in the presence of growth-dependent mortality, and
for different size-dependent mortality landscapes. Its advantage stems from
epigenetic heritability of cell size, which enables selection to act on the
distribution of cell sizes in a population to avoid mortality thresholds and
adapt to a mortality landscape.


  We are keenly interested in finding the limiting behaviors of total tumor
size when tumor cells are subject to the periodic repetition of therapy and
rest periods, called intermittent cyclic therapy. We hypothesize that each
tumor cell can take either therapy-sensitive or therapy-tolerant phenotype, its
phenotype transition is mainly driven by the presence or absence of
environmental stress, and such a transition is reversible. Even though those
aforementioned hypotheses make the model system simple, most of prior papers
attempted to numerically find the optimal therapeutic scheduling that minimizes
total tumor size, and there is no rigorous proof of the limiting behaviors of
total tumor size to my knowledge. Here we present such long-waited
mathematically rigorous results. In the first part of the paper, we present the
derivation of total tumor size reduction criterion and prove two theorems of
two different limiting behaviors of total tumor size under two different
therapy strategies, one leading to an asymptotic finite tumor size according to
an iterated map method and anther leading to asymptotically diminishing of
total tumor size. In the second part of the paper, we discuss the effects of
the intratumoral competition between sensitive and tolerant phenotypes on the
total tumor size reduction criterion.


  Social network structures play an important role in the lives of humans and
non-human animals by affecting wellbeing, the spread of disease and
information, and evolutionary processes. Nevertheless, we still lack a good
understanding of how these structures emerge from individual behaviour. Here we
present a general model for the emergence of social structures, which is based
on a key aspect of real social systems observed across species, namely social
preferences for traits (individual characteristics such as age, sex, etc.). We
first show that the model can generate diverse artificial social structures,
and consider its potential for being combined with real network data. We then
use the model to gain fundamental insights into how two main categories of
social preferences (similarity and popularity) affect social structure and
function. The results show that the types of social preference, in combination
with the types of trait they are used with, can have important consequences for
the spread of information and disease, and the robustness of social structures
against fragmentation. The results also suggest that symmetric degree
distributions could be expected to be common in social networks. More
generally, the study implies that trait-based social preferences can have
consequences for social systems that go far beyond their effect on direct
benefits from social partners. We discuss the implications of the results for
social evolution.


  Predator prey interactions are one of ecology's central research themes, but
with many interdisciplinary implications across the social and natural
sciences. Here we consider an often-overlooked species in these interactions,
namely parasites. We first show that a simple predator prey parasite model,
inspired by the classical Lotka Volterra equations, fails to produce a stable
coexistence of all three species, thus failing to provide a biologically
realistic outcome. To improve this, we introduce free space as a relevant
eco-evolutionary component in a new mathematical model that uses a
game-theoretical payoff matrix to describe a more realistic setup. We then show
that the consideration of free space stabilizes the dynamics by means of cyclic
dominance that emerges between the three species. We determine the parameter
regions of coexistence as well as the types of bifurcations leading to it by
means of analytical derivations as well as by means of numerical simulations.
We conclude that the consideration of free space as a finite resource reveals
the limits of biodiversity in predator prey parasite interactions, and it may
also help us in the determination of factors that promote a healthy biota.


  We present a microscopic model of replicator dynamics with strategy-dependent
time delays. In such a model, new players are born from parents who interacted
and received payoffs in the past. In the case of small delays, we use Taylor
expansion to get ordinary differential equations for frequencies of strategies
with time delays as parameters. We apply our technique to get analytic
expressions for interior stationary states in two games: Snowdrift and
Stag-hunt. We show that interior stationary states depend continuously upon
time delays. Our analytic formulas for stationary states approximate well exact
numerical results for small time delays.


  Tumor heterogeneity is a complex and widely recognized trait that poses
significant challenges in developing effective cancer therapies. In particular,
many tumors harbor a variety of subpopulations with distinct therapeutic
response characteristics. Characterizing this heterogeneity by determining the
subpopulation structure within a tumor enables more precise and successful
treatment strategies. In our prior work, we developed PhenoPop, a computational
framework for unravelling the drug-response subpopulation structure within a
tumor from bulk high-throughput drug screening data. However, the deterministic
nature of the underlying models driving PhenoPop restricts the model fit and
the information it can extract from the data. As an advancement, we propose a
stochastic model based on the linear birth-death process to address this
limitation. Our model can formulate a dynamic variance along the horizon of the
experiment so that the model uses more information from the data to provide a
more robust estimation. In addition, the newly proposed model can be readily
adapted to situations where the experimental data exhibits a positive time
correlation. We test our model on simulated data (in silico) and experimental
data (in vitro), which supports our argument about its advantages.


  Mutually beneficial interactions between plant and pollinators play an
essential role in the biodiversity, stability of the ecosystem and crop
production. Despite their immense importance, rapid decline events of
pollinators are common worldwide in past decades. Excessive use of chemical
pesticides is one of the most important threat to pollination in the current
era of anthropogenic changes. Pesticides are applied to the plants to increase
their growth by killing harmful pests and pollinators accumulates toxic
pesticides from the interacting plants directly from the nectar and pollen.
This has a significant adverse effect on the pollinator growth and the
mutualism which in turn can cause an abrupt collapse of the community however
predicting the fate of such community dynamics remains a blur under the
alarming rise in the dependency of chemical pesticides. We mathematically
modeled the influence of pesticides in a multispecies mutualistic community and
used 105 real plant-pollinator networks sampled worldwide as well as simulated
networks, to assess its detrimental effect on the plant-pollinator mutualistic
networks. Our results indicate that the persistence of the community is
strongly influenced by the level of pesticide and catastrophic and irreversible
community collapse may occur due to pesticide. Furthermore, a species rich,
highly nested community with low connectance and modularity has greater
potential to function under the influence of pesticide. We finally proposed a
realistic intervention strategy which involves the management of the pesticide
level of one targeted plant from the community. We show that our intervention
strategy can significantly delay the collapse of the community. Overall our
study can be considered as the first attempt to understand the consequences of
the chemical pesticide on a plant-pollinator mutualistic community.


  We apply a perturbative Doi--Peliti field-theoretical analysis to the
stochastic spatially extended symmetric Rock-Paper-Scissors (RPS) and
May--Leonard (ML) models, in which three species compete cyclically. Compared
to the two-species Lotka--Volterra predator-prey (LV) model, according to
numerical simulations, these cyclical models appear to be less affected by
intrinsic stochastic fluctuations. Indeed, we demonstrate that the qualitative
features of the ML model are insensitive to intrinsic reaction noise. In
contrast, and although not yet observed in numerical simulations, we find that
the RPS model acquires significant fluctuation-induced renormalizations in the
perturbative regime, similar to the LV model. We also study the formation of
spatio-temporal structures in the framework of stability analysis and provide a
clearcut explanation for the absence of spatial patterns in the RPS system,
whereas the spontaneous emergence of spatio-temporal structures features
prominently in the LV and the ML models.


  The number of available genomes of prokaryotic organisms is rapidly growing
enabling comparative genomics studies. The comparison of genomes of organisms
with a common phenotype, habitat or phylogeny often shows that these genomes
share some common contents. Collecting rules expressing common genome traits
depending on given factors is useful, as such rules could be used for quality
control or for identifying interesting exceptions and formulating hypothesis.
Automatizing the rules verification using computation tools requires the
definition of a representation schema. In this study, we present EGC (Expected
Genome Contents), a flat-text file format for the representation of expectation
rules about the content of prokaryotic genomes. A parser for the EGC format has
been implemented using the TextFormats software library, accompanied by a set
of related Python packages.


  Figuring out small molecule binding sites in target proteins, in the
resolution of either pocket or residue, is critical in many virtual and real
drug-discovery scenarios. Since it is not always easy to find such binding
sites based on domain knowledge or traditional methods, different deep learning
methods that predict binding sites out of protein structures have been
developed in recent years. Here we present a new such deep learning algorithm,
that significantly outperformed all state-of-the-art baselines in terms of the
both resolutions$\unicode{x2013}$pocket and residue. This good performance was
also demonstrated in a case study involving the protein human serum albumin and
its binding sites. Our algorithm included new ideas both in the model
architecture and in the training method. For the model architecture, it
incorporated SE(3)-invariant geometric self-attention layers that operate on
top of residue-level CNN outputs. This residue-level processing of the model
allowed a transfer learning between the two resolutions, which turned out to
significantly improve the binding pocket prediction. Moreover, we developed
novel augmentation method based on protein homology, which prevented our model
from over-fitting. Overall, we believe that our contribution to the literature
is twofold. First, we provided a new computational method for binding site
prediction that is relevant to real-world applications, as shown by the good
performance on different benchmarks and case study. Second, the novel ideas in
our method$\unicode{x2013}$the model architecture, transfer learning and the
homology augmentation$\unicode{x2013}$would serve as useful components in
future works.


  Living organisms continuously harness energy to perform complex functions for
their adaptation and survival while part of that energy is dissipated in the
form of heat or chemical waste. Determining the energetic cost and the
efficiency of specific cellular processes remains a largely open problem. Here,
we analyze the efficiency of mitochondrial adenosine triphosphate (ATP)
production through the tricarboxylic acid (TCA) cycle and oxidative
phosphorylation that generates most of the cellular chemical energy in
eukaryotes. The regulation of this pathway by calcium signaling represents a
well-characterized example of a regulatory cross-talk that can affect the
energetic output of a metabolic pathway, but its concrete energetic impact
remains elusive. On the one hand, calcium enhances ATP production by activating
key enzymes of the TCA cycle, but on the other hand calcium homeostasis depends
on ATP availability. To evaluate how calcium signaling impacts the efficiency
of mitochondrial metabolism, we propose a detailed kinetic model describing the
calcium-mitochondria cross-talk and we analyze it using a nonequilibrium
thermodynamic approach: after identifying the effective reactions driving
mitochondrial metabolism out of equilibrium, we quantify the thermodynamic
efficiency of the metabolic machinery for different physiological conditions.
We find that calcium oscillations increase the efficiency with a maximum close
to substrate-limited conditions, suggesting a compensatory effect of calcium
signaling on the energetics of mitochondrial metabolism.


  We study the spread of multi-competitive viruses over a (possibly)
time-varying network of individuals accounting for the presence of shared
infrastructure networks that further enables transmission of the virus. We
establish a sufficient condition for exponentially fast eradication of a virus
for: 1) time-invariant graphs, 2) time-varying graphs with symmetric
interactions between individuals and homogeneous virus spread across the
network (same healing and infection rate for all individuals), and 3) directed
and slowly varying graphs with heterogeneous virus spread (not necessarily same
healing and infection rates for all individuals) across the network. Numerical
examples illustrate our theoretical results and indicate that, for the
time-varying case, violation of the aforementioned sufficient conditions could
lead to the persistence of a virus.


  The first important step in extracting DNA characters is using the output
data of MinION devices in the form of electrical current signals. Various
cutting-edge base callers use this data to detect the DNA characters based on
the input. In this paper, we discuss several shortcomings of prior base callers
in the case of time-critical applications, privacy-aware design, and the
problem of catastrophic forgetting. Next, we propose the LRDB model, a
lightweight open-source model for private developments with a better
read-identity (0.35% increase) for the target bacterial samples in the paper.
We have limited the extent of training data and benefited from the transfer
learning algorithm to make the active usage of the LRDB viable in critical
applications. Henceforth, less training time for adapting to new DNA samples
(in our case, Bacterial samples) is needed. Furthermore, LRDB can be modified
concerning the user constraints as the results show a negligible accuracy loss
in case of using fewer parameters. We have also assessed the noise-tolerance
property, which offers about a 1.439% decline in accuracy for a 15dB noise
injection, and the performance metrics show that the model executes in a medium
speed range compared with current cutting-edge models.


  Simulations of biomolecules have enormous potential to inform our
understanding of biology but require extremely demanding calculations. For over
twenty years, the Folding@home distributed computing project has pioneered a
massively parallel approach to biomolecular simulation, harnessing the
resources of citizen scientists across the globe. Here, we summarize the
scientific and technical advances this perspective has enabled. As the
project's name implies, the early years of Folding@home focused on driving
advances in our understanding of protein folding by developing statistical
methods for capturing long-timescale processes and facilitating insight into
complex dynamical processes. Success laid a foundation for broadening the scope
of Folding@home to address other functionally relevant conformational changes,
such as receptor signaling, enzyme dynamics, and ligand binding. Continued
algorithmic advances, hardware developments such as GPU-based computing, and
the growing scale of Folding@home have enabled the project to focus on new
areas where massively parallel sampling can be impactful. While previous work
sought to expand toward larger proteins with slower conformational changes, new
work focuses on large-scale comparative studies of different protein sequences
and chemical compounds to better understand biology and inform the development
of small molecule drugs. Progress on these fronts enabled the community to
pivot quickly in response to the COVID-19 pandemic, expanding to become the
world's first exascale computer and deploying this massive resource to provide
insight into the inner workings of the SARS-CoV-2 virus and aid the development
of new antivirals. This success provides a glimpse of what's to come as
exascale supercomputers come online, and Folding@home continues its work.


  We have developed a mathematical model that captures stress-induced
mutagenesis, a fundamental aspect of pathogenic and neoplastic evolutionary
dynamics, on the fitness landscape with multiple relevant genetic traits as a
high-dimensional Euclidean space. In this framework, stress-induced mutagenesis
manifests as a heterogeneous diffusion process. We show how increasing
mutations, and thus reducing exploitation, in a static ecology with fixed
carrying capacity and maximum growth rates, can paradoxically boost population
size. Remarkably, this unexpected biophysical phenomenon applies universally to
any number of traits.


  This work studies fundamental questions regarding the optimal design of
antimicrobial treatment protocols, using standard pharmacodynamic and
pharmacokinetic mathematical models. We consider the problem of designing an
antimicrobial treatment schedule to achieve eradication of a microbial
infection, while minimizing the area under the time-concentration curve (AUC).
We first solve this problem under the assumption that an arbitrary
antimicrobial concentration profile may be chosen, and prove that the 'ideal'
concentration profile consists of a constant concentration over a finite time
duration, where explicit expressions for the optimal concentration and the time
duration are given in terms of the pharmacodynamic parameters. Since
antimicrobial concentration profiles are induced by a dosing schedule and the
antimicrobial pharmacokinetics, the ideal concentration profile is not strictly
feasible. We therefore also investigate the possibility of achieving outcomes
which are close to those provided by the ideal concentration profile,using a
bolus+continuous dosing schedule, which consists of a loading dose followed by
infusion of the antimicrobial at a constant rate. We explicitly find the
optimal bolus+continuous dosing schedule, and show that, for realistic
parameter ranges, this schedule achieves results which are nearly as efficient
as those attained by the ideal concentration profile. The optimality results
obtained here provide a baseline and reference point for comparison and
evaluation of antimicrobial treatment plans.


  A novel method combining maximum entropy principle, the Bayesian-inference of
ensembles approach, and the optimization of empirical forward models is
presented. Here we focus on the Karplus parameters for RNA systems, which
relate the dihedral angles of $\gamma$, $\beta$, and the dihedrals in the sugar
ring to the corresponding $^3J$-coupling signal between coupling protons.
Extensive molecular simulations are performed on a set of RNA tetramers and
hexamers and combined with available nucleic-magnetic-resonance data. Within
the new framework, the sampled structural dynamics can be reweighted to match
experimental data while the error arising from inaccuracies in the forward
models can be corrected simultaneously and consequently does not leak into the
reweighted ensemble. Carefully crafted cross-validation procedure and
regularization terms enable obtaining transferable Karplus parameters. Our
approach identifies the optimal regularization strength and new sets of Karplus
parameters balancing good agreement between simulations and experiments with
minimal changes to the original ensemble.


  Non-equilibrium and instability features of prey-predator-like systems
associated to topological quantum domains emerging from a quantum phase-space
description are investigated in the framework of the Weyl-Wigner quantum
mechanics. Reporting about the generalized Wigner flow for one-dimensional
Hamiltonian systems, $\mathcal{H}(x,\,k)$, constrained by $\partial^2
\mathcal{H} / \partial x \, \partial k = 0$, the prey-predator dynamics driven
by Lotka-Volterra (LV) equations is mapped onto the Heisenberg-Weyl
non-commutative algebra, $[x,\,k] = i$, where the canonical variables $x$ and
$k$ are related to the two-dimensional LV parameters, $y = e^{-x}$ and $z =
e^{-k}$. From the non-Liouvillian pattern driven by the associated Wigner
currents, hyperbolic equilibrium and stability parameters for the
prey-predator-like dynamics are then shown to be affected by quantum
distortions over the classical background, in correspondence with
non-stationarity and non-Liouvillianity properties quantified in terms of
Wigner currents and Gaussian ensemble parameters. As an extension, considering
the hypothesis of discretizing the time parameter, non-hyperbolic bifurcation
regimes are identified and quantified in terms of $z-y$ anisotropy and Gaussian
parameters. The bifurcation diagrams exhibit, for quantum regimes, chaotic
patterns highly dependent on Gaussian localization. Besides exemplifying a
broad range of applications of the generalized Wigner information flow
framework, our results extend, from the continuous (hyperbolic regime) to
discrete (chaotic regime) domains, the procedure for quantifying the influence
of quantum fluctuations over equilibrium and stability scenarios of LV driven
systems.


  A new three dimensional approach to the chaos game representation of protein
sequences is explored in this thesis. The basics of DNA, the synthesis of
proteins from DNA, protein structure and functionality and sequence alignment
techniques are presented. The mathematical background needed for understanding
the chaos game representation and fractal analysis are briefly discussed.
  An account of the existing literature on the chaos game representation of DNA
sequences and a detailed account of the chaos game representation of protein
sequences in two dimensions with its advantages and limitations are presented.
We explore a new three dimensional approach to the chaos game representation of
protein sequences (3D-CGR) and study its ability a) to determine protein
sequence similarity and differences, b) to study the effect of dinucleotide
biases at amino acid level on the 3D-CGR derived protein homology, and c) to
identify sequence similarity based on shuffled motifs that could be used for
studying protein evolution due to exon shuffling.


  We study analytically and numerically a model metabolic cycle composed of an
arbitrary number of species of catalytically active particles. Each species
converts a substrate into a product, the latter being used as the substrate by
the next species in the cycle. Through a combination of catalytic activity and
chemotactic mobility, the active particles develop effective non-reciprocal
interactions with particles belonging to neighbouring species in the cycle. We
find that such model metabolic cycles are able to self-organize through a
macroscopic instability, with a strong dependence on the number of species they
incorporate. The parity of that number has a key influence: cycles containing
an even number of species are able to minimize repulsion between their
component particles by aggregating all even-numbered species in one cluster,
and all odd-numbered species in another. Such a grouping is not possible if the
cycle contains an odd number of species, which can lead to oscillatory steady
states in the case of chasing interactions.


  Tumour heterogeneity in breast cancer poses challenges in predicting outcome
and response to therapy. Spatial transcriptomics technologies may address these
challenges, as they provide a wealth of information about gene expression at
the cell level, but they are expensive, hindering their use in large-scale
clinical oncology studies. Predicting gene expression from hematoxylin and
eosin stained histology images provides a more affordable alternative for such
studies. Here we present BrST-Net, a deep learning framework for predicting
gene expression from histopathology images using spatial transcriptomics data.
Using this framework, we trained and evaluated 10 state-of-the-art deep
learning models without utilizing pretrained weights for the prediction of 250
genes. To enhance the generalisation performance of the main network, we
introduce an auxiliary network into the framework. Our methodology outperforms
previous studies, with 237 genes identified with positive correlation,
including 24 genes with a median correlation coefficient greater than 0.50.
This is a notable improvement over previous studies, which could predict only
102 genes with positive correlation, with the highest correlation values
ranging from 0.29 to 0.34.


  The search for biosignatures necessitates developing our understanding of
life under different conditions. If life can influence the climate evolution of
its planet then understanding the behaviour of life-climate feedbacks under
extreme conditions is key to determine the 'edges' of the habitable zone.
Additionally understanding the behaviour of a temperature limited biosphere
will help towards formulating biosignature predictions for alien life living
under conditions very different to those on Earth. Towards this aim, we extend
the 'ExoGaia Model' - an abstract model of microbial life living on a highly
simplified 0-dimensional planet. Via their metabolisms, microbes influence the
atmospheric composition and therefore the temperature of the planet and
emergent feedback loops allow microbes to regulate their climate and maintain
long term habitability. Here, we adapt the ExoGaia model to include temperature
adaptation of the microbes by allowing different species to have different
temperature 'preferences'. We find that rather than adapting towards the
planet's abiotic conditions the biosphere tends to more strongly influence the
climate of its planet, suggesting that the surface temperature of an inhabited
planet might be significantly different from that predicted using abiotic
models. We find that the success rate for microbial establishment on planets is
improved when adaptation is allowed. However, planetary abiotic context is
important for determining whether overall survival prospects for life will be
improved or degraded. These results indicate the necessity to develop an
understanding of life living under different limiting regimes to form
predictions for the boundaries of the habitable zone.


  Although promising for biomedicine, the clinical translation of inorganic
nanoparticles (NPs) is limited by low biocompatibility and stability in
biological fluids. A common strategy to circumvent this drawback consists in
disguising the active inorganic core with a lipid bilayer coating, reminiscent
of the structure of the cell membrane to redefine the chemical and biological
identity of NPs. While recent reports introduced membrane coating procedures
for NPs, a robust and accessible method to quantify the integrity of the
bilayer coverage is not yet available. To fill this gap, we prepared SiO2
nanoparticles (SiO2NPs) with different membrane coverage degrees and monitored
their interaction with AuNPs by combining microscopic, scattering, and optical
techniques. The membrane-coating on SiO2NPs induces spontaneous clustering of
AuNPs, whose extent depends on the coating integrity. Remarkably, we discovered
a linear correlation between the membrane coverage and a spectral descriptor
for the AuNPs plasmonic resonance, spanning a wide range of coating yields.
These results provide a fast and cost-effective assay to monitor the
compatibilization of NPs with biological environments, essential for bench
tests and scale-up. In addition, we introduce a robust and scalable method to
prepare SiO2NPsAuNPs hybrids through spontaneous self assembly, with a high
fidelity structural control mediated by a lipid bilayer.


  Molecular fingerprints are widely used for predicting chemical properties,
and selecting appropriate fingerprints is important. We generate new
fingerprints based on the assumption that a performance of prediction using a
more effective fingerprint is better. We generate effective interaction
fingerprints that are the product of multiple base fingerprints. It is
difficult to evaluate all combinations of interaction fingerprints because of
computational limitations. Against this problem, we transform a problem of
searching more effective interaction fingerprints into a quadratic
unconstrained binary optimization problem. In this study, we found effective
interaction fingerprints using QM9 dataset.


  China adjusted the zero-COVID strategy in late 2022, triggering an
unprecedented Omicron wave. We estimated the time-varying reproduction numbers
of 32 provincial-level administrative divisions from December 2022 to January
2023. We found that the pooled estimate of initial reproduction numbers is 4.74
(95% CI: 4.41, 5.07).


  We study a three-species cyclic model whose organisms are vulnerable to
contamination with an infectious disease which propagates person-to-person. We
consider that individuals of one species perform an evolutionary
self-preservation strategy by reducing the mobility rate to minimise infection
risk whenever an epidemic outbreak reaches the neighbourhood. Running
stochastic simulations, we quantify the changes in spatial patterns induced by
unevenness in the cyclic game introduced by the mobility restriction strategy
of organisms of one out of the species. Our findings show that variations in
disease virulence impact the benefits of dispersal limitation reaction, with
the relative reduction of the organisms' infection risk accentuating in surges
of less contagious or deadlier diseases. The effectiveness of the mobility
restriction tactic depends on the tolerable fraction of infected neighbours
used as a trigger of the defensive strategy and the deceleration level. If each
organism promptly reacts to the arrival of the first viral vectors in its
surroundings with strict mobility reduction, contamination risk decreases
significantly. Our conclusions may help biologists understand the impact of
evolutionary defensive strategies in ecosystems during an epidemic.


  DNA supercoiling is central to many fundamental processes of living
organisms. Its average level along the chromosome and over time reflects the
dynamic equilibrium of opposite activities of topoisomerases, which are
required to relax mechanical stresses that are inevitably produced during DNA
replication and gene transcription. Supercoiling affects all scales of the
spatio-temporal organization of bacterial DNA, from the base pair to the large
scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s
and 1970s, respectively, the first physical models were proposed concomitantly
in order to predict the deformation properties of the double helix. About
fifteen years later, polymer physics models demonstrated on larger scales the
plectonemic nature and the tree-like organization of supercoiled DNA. Since
then, many works have tried to establish a better understanding of the multiple
structuring and physiological properties of bacterial DNA in thermodynamic
equilibrium and far from equilibrium. The purpose of this essay is to address
upcoming challenges by thoroughly exploring the relevance, predictive capacity,
and limitations of current physical models, with a specific focus on structural
properties beyond the scale of the double helix. We discuss more particularly
the problem of DNA conformations, the interplay between DNA supercoiling with
gene transcription and DNA replication, its role on nucleoid formation and,
finally, the problem of scaling up models. Our primary objective is to foster
increased collaboration between physicists and biologists. To achieve this, we
have reduced the respective jargon to a minimum and we provide some explanatory
background material for the two communities.


  This review article provides an overview of structurally oriented,
experimental datasets that can be used to benchmark protein force fields,
focusing on data generated by nuclear magnetic resonance (NMR) spectroscopy and
room temperature (RT) protein crystallography. We discuss why these observables
are useful for assessing force field accuracy, how they can be calculated from
simulation trajectories, and statistical issues that arise when comparing
simulations with experiment. The target audience for this article is
computational researchers and trainees who develop, benchmark, or use protein
force fields for molecular simulations.


  Understanding the molecular processes that drive cellular life is a
fundamental question in biological research. Ambitious programs have gathered a
number of molecular datasets on large populations. To decipher the complex
cellular interactions, recent work has turned to supervised machine learning
methods. The scientific questions are formulated as classical learning problems
on tabular data or on graphs, e.g. phenotype prediction from gene expression
data. In these works, the input features on which the individual predictions
are predominantly based are often interpreted as indicative of the cause of the
phenotype, such as cancer identification. Here, we propose to explore the
relevance of the biomarkers identified by Integrated Gradients, an
explainability method for feature attribution in machine learning. Through a
motivating example on The Cancer Genome Atlas, we show that ranking features by
importance is not enough to robustly identify biomarkers. As it is difficult to
evaluate whether biomarkers reflect relevant causes without known ground truth,
we simulate gene expression data by proposing a hierarchical model based on
Latent Dirichlet Allocation models. We also highlight good practices for
evaluating explanations for genomics data and propose a direction to derive
more insights from these explanations.


  In this paper we develop explicit and semi-implicit second-order
high-resolution finite difference schemes for a structured
coagulation-fragmentation model formulated on the space of Radon measures. We
prove the convergence of each of the two schemes to the unique weak solution of
the model. We perform numerical simulations to demonstrate that the second
order accuracy is achieved by both schemes.


  Understanding how cooperative behaviours can emerge from competitive
interactions is an open problem in biology and social sciences. While
interactions are usually modelled as pairwise networks, the units of many
real-world systems can also interact in groups of three or more. Here, we
introduce a general framework to extend pairwise games to higher-order
networks. By studying social dilemmas on hypergraphs with a tunable structure,
we find an explosive transition to cooperation triggered by a critical number
of higher-order games. The associated bistable regime implies that an initial
critical mass of cooperators is also required for the emergence of prosocial
behavior. Our results show that higher-order interactions provide a novel
explanation for the survival of cooperation.


  Most widely used ligand docking methods assume a rigid protein structure.
This leads to problems when the structure of the target protein deforms upon
ligand binding. In particular, the ligand's true binding pose is often scored
very unfavorably due to apparent clashes between ligand and protein atoms,
which lead to extremely high values of the calculated van der Waals energy
term. Traditionally, this problem has been addressed by explicitly searching
for receptor conformations to account for the flexibility of the receptor in
ligand binding. Here we present a deep learning model trained to take receptor
flexibility into account implicitly when predicting van der Waals energy. We
show that incorporating this machine-learned energy term into a
state-of-the-art physics-based scoring function improves small molecule ligand
pose prediction results in cases with substantial protein deformation, without
degrading performance in cases with minimal protein deformation. This work
demonstrates the feasibility of learning effects of protein flexibility on
ligand binding without explicitly modeling changes in protein structure.


  Recent years have seen rapid development of descriptor generation based on
representation learning of extremely diverse molecules, especially those that
apply natural language processing (NLP) models to SMILES, a literal
representation of molecular structure. However, little research has been done
on how these models understand chemical structure. To address this black box,
we investigated the relationship between the learning progress of SMILES and
chemical structure using a representative NLP model, the Transformer. We show
that while the Transformer learns partial structures of molecules quickly, it
requires extended training to understand overall structures. Consistently, the
accuracy of molecular property predictions using descriptors generated from
models at different learning steps was similar from the beginning to the end of
training. Furthermore, we found that the Transformer requires particularly long
training to learn chirality and sometimes stagnates with low performance due to
misunderstanding of enantiomers. These findings are expected to deepen the
understanding of NLP models in chemistry.


  In this article, we develop a simple mathematical GNU Octave/MATLAB code that
is easy to modify for the simulation of mathematical models governed by
fractional-order differential equations, and for the resolution of
fractional-order optimal control problems through Pontryagin's maximum
principle (indirect approach to optimal control). For this purpose, a
fractional-order model for the respiratory syncytial virus (RSV) infection is
considered. The model is an improvement of one first proposed by the authors in
[Chaos Solitons Fractals 117 (2018), 142--149]. The initial value problem
associated with the RSV infection fractional model is numerically solved using
Garrapa's fde12 solver and two simple methods coded here in Octave/MATLAB: the
fractional forward {Euler's} method and the predict-evaluate-correct-evaluate
(PECE) method of Adams--Bashforth--Moulton. A fractional optimal control
problem is then formulated having treatment as the control. The fractional
Pontryagin maximum principle is used to characterize the fractional optimal
control and the extremals of the problem are determined numerically through the
implementation of the forward-backward PECE method. The implemented algorithms
are available on GitHub and, at the end of the paper, in appendixes, both for
the uncontrolled initial value problem as well as for the fractional optimal
control problem, using the free GNU Octave computing software and assuring
compatibility with~MATLAB.


  In this work we present a deep learning approach to conduct hypothesis-free,
transcriptomics-based matching of drugs for diseases. Our proposed neural
network architecture is trained on approved drug-disease indications, taking as
input the relevant disease and drug differential gene expression profiles, and
learns to identify novel indications. We assemble an evaluation dataset of
disease-drug indications spanning 68 diseases and evaluate in silico our
approach against the most widely used transcriptomics-based matching baselines,
CMap and the Characteristic Direction. Our results show a more than 200%
improvement over both baselines in terms of standard retrieval metrics. We
further showcase our model's ability to capture different genes' expressions
interactions among drugs and diseases. We provide our trained models, data and
code to predict with them at https://github.com/healx/dgem-nn-public.


  The interaction of phase-separating systems with chemical reactions is of
great interest in various contexts, from biology to material science. In
biology, phase separation is thought to be the driving force behind the
formation of biomolecular condensates, i.e. organelles without a membrane that
are associated with cellular metabolism, stress response, and development. RNA,
proteins, and small molecules participating in the formation of condensates are
also involved in a variety of biochemical reactions: how do the chemical
reaction dynamics influence the process of phase separation? Here we are
interested in finding chemical reactions that can arrest the growth of
condensates, generating stable spatial patterns of finite size (microphase
separation), in contrast with the otherwise spontaneous (unstable) growth of
condensates. We consider a classical continuum model for phase separation
coupled to a chemical reaction network (CRN), and we seek conditions for the
emergence of stable oscillations of the solution in space. Given reaction
dynamics with uncertain rate constants, but known structure, we derive easily
computable conditions to assess whether microphase separation is impossible,
possible for some parameter values, or robustly guaranteed for all parameter
values within given bounds. Our results establish a framework to evaluate which
classes of CRNs favor the emergence of condensates with finite size, a question
that is broadly relevant to understanding and engineering life.


  Pretrained protein structure models without labels are crucial foundations
for the majority of protein downstream applications. The conventional structure
pretraining methods follow the mature natural language pretraining methods such
as denoised reconstruction and masked language modeling but usually destroy the
real representation of spatial structures. The other common pretraining methods
might predict a fixed set of predetermined object categories, where a
restricted supervised manner limits their generality and usability as
additional labeled data is required to specify any other protein concepts. In
this work, we introduce a novel unsupervised protein structure representation
pretraining with a robust protein language model. In particular, we first
propose to leverage an existing pretrained language model to guide structure
model learning through an unsupervised contrastive alignment. In addition, a
self-supervised structure constraint is proposed to further learn the intrinsic
information about the structures. With only light training data, the pretrained
structure model can obtain better generalization ability. To quantitatively
evaluate the proposed structure models, we design a series of rational
evaluation methods, including internal tasks (e.g., contact map prediction,
distribution alignment quality) and external/downstream tasks (e.g., protein
design). The extensive experimental results conducted on multiple tasks and
specific datasets demonstrate the superiority of the proposed
sequence-structure transformation framework.


  The goal of most materials discovery is to discover materials that are
superior to those currently known. Fundamentally, this is close to
extrapolation, which is a weak point for most machine learning models that
learn the probability distribution of data. Herein, we develop AI-driven
combinatorial chemistry, which is a rule-based inverse molecular designer that
does not rely on data. Since our model has the potential to generate all
possible molecular structures that can be obtained from combinations of
molecular fragments, unknown materials with superior properties can be
discovered. We theoretically and empirically demonstrate that our model is more
suitable for discovering better materials than probability
distribution-learning models. In an experiment aimed at discovering molecules
that hit seven target properties, our model discovered 1,315 of all
target-hitting molecules and 7,629 of five target-hitting molecules out of
100,000 trials, whereas the probability distribution-learning models failed. To
illustrate the performance in actual problems, we also demonstrate that our
models work well on two practical applications: discovering protein docking
materials and HIV inhibitors.


  Microbes are often discussed in terms of dichotomies such as
copiotrophic/oligotrophic and fast/slow-growing microbes, defined using the
characterisation of microbial growth in isolated cultures. The dichotomies are
usually qualitative and/or study-specific, sometimes precluding clear-cut
results interpretation. We are able to interpret microbial dichotomies as life
history strategies by combining ecology theory with Monod curves, a classical
laboratory tool of bacterial physiology. Monod curves relate the specific
growth rate of a microbe with the concentration of a limiting nutrient, and
provide quantities that directly correspond to key ecological parameters in
McArthur and Wilsons r/K selection theory, Tilmans resource competition and
community structure theory and Grimes triangle of life strategies. The
resulting model allows us to reconcile the copiotrophic/oligotrophic and
fast/slow-growing dichotomies as different subsamples of a life history
strategy triangle that also includes r/K strategists. We analyzed some
ecological context by considering the known viable carbon sources for
heterotrophic microbes in the framework of community structure theory. This
partly explains the microbial diversity observed using metagenomics. In sum,
ecology theory in combination with Monod curves can be a unifying quantitative
framework for the study of natural microbial communities, calling for the
integration of modern laboratory and field experiments.


  The human gut microbiota is a complex ecosystem that affects a range of human
physiology. In order to explore the dynamics of the human gut microbiota, we
used a system of ordinary differential equations to model mathematically the
biomass of three microorganism populations: $\textit{Bacteroides
thetaiotaomicron}$, $\textit{Eubacterium rectale}$, and
$\textit{Methanobrevibacter smithii}$. Additionally, we modeled the
concentrations of relevant nutrients necessary to sustain these populations
over time. Our model highlights the interactions and the competition among
these three species. These three microorganisms were specifically chosen due to
the system's end product, butyrate, which is a short chain fatty acid that aids
in developing and maintaining the intestinal barrier in the human gut. The
basis of our mathematical model assumes the gut is structured such that
bacteria and nutrients exit the gut at a rate proportional to its volume, the
rate of volumetric flow, and the biomass or concentration of the particular
population or nutrient. We performed global sensitivity analyses using Sobol'
sensitivities to estimate the relative importance of model parameters on
simulation results.


  Identifiability of the discrete tree parameter is a key property for
phylogenetic models since it is necessary for statistically consistent
estimation of the tree from sequence data. Algebraic methods have proven to be
very effective at showing that tree and network parameters of phylogenetic
models are identifiable, especially when the underlying models are group-based.
However, since group-based models are time-reversible, only the unrooted tree
topology is identifiable and the location of the root is not. In this note we
show that the rooted tree parameter of the Cavender-Farris-Neyman Model with a
Molecular Clock is generically identifiable by using the invariants of the
model which were characterized by Coons and Sullivant.


  The branching of an RNA molecule is an important structural characteristic
yet difficult to predict correctly, especially for longer sequences. Using
plane trees as a combinatorial model for RNA folding, we consider the
thermodynamic cost, known as the barrier height, of transitioning between
branching configurations. Using branching skew as a coarse energy
approximation, we characterize various types of paths in the discrete
configuration landscape. In particular, we give sufficient conditions for a
path to have both minimal length and minimal branching skew. The proofs offer
some biological insights, notably the potential importance of both hairpin
stability and domain architecture to higher resolution RNA barrier height
analyses.


  One central goal of systems biology is to infer biochemical regulations from
large-scale OMICS data. Many aspects of cellular physiology and organism
phenotypes could be understood as a result of the metabolic interaction network
dynamics. Previously, we have derived a mathematical method addressing this
problem using metabolomics data for the inverse calculation of a biochemical
Jacobian network. However, these algorithms for this inference are limited by
two issues: they rely on structural network information that needs to be
assembled manually, and they are numerically unstable due to ill-conditioned
regression problems, which makes them inadequate for dealing with large-scale
metabolic networks. In this work, we present a novel regression-loss based
inverse Jacobian algorithm and related workflow COVRECON. It consists of two
parts: a, Sim-Network and b, Inverse differential Jacobian evaluation.
Sim-Network automatically generates an organism-specific enzyme and reaction
dataset from Bigg and KEGG databases, which is then used to reconstruct the
Jacobian's structure for a specific metabolomics dataset. Instead of directly
solving a regression problem, the new inverse differential Jacobian part is
based on a more robust approach and rates the biochemical interactions
according to their relevance from large-scale metabolomics data. This approach
is illustrated by in silico stochastic analysis with different-sized metabolic
networks from the BioModels database. The advantages of COVRECON are that 1) it
automatically reconstructs a data-driven superpathway metabolic interaction
model; 2) more general network structures can be considered; 3) the new inverse
algorithms improve stability, decrease computation time, and extend to
large-scale models


  13C-Metabolic Flux Analysis (13C-MFA) and Flux Balance Analysis (FBA) are
widely used to investigate the operation of biochemical networks in both
biological and biotechnological research. Both of these methods use metabolic
reaction network models of metabolism operating at steady state, so that
reaction rates (fluxes) and the levels of metabolic intermediates are
constrained to be invariant. They provide estimated (MFA) or predicted (FBA)
values of the fluxes through the network in vivo, which cannot be measured
directly. A number of approaches have been taken to test the reliability of
estimates and predictions from constraint-based methods and to decide on and/or
discriminate between alternative model architectures. Despite advances in other
areas of the statistical evaluation of metabolic models, validation and model
selection methods have been underappreciated and underexplored. We review the
history and state-of-the-art in constraint-based metabolic model validation and
model selection. Applications and limitations of the X2-test of
goodness-of-fit, the most widely used quantitative validation and selection
approach in 13C-MFA, are discussed, and complementary and alternative forms of
validation and selection are proposed. A combined model validation and
selection framework for 13C-MFA incorporating metabolite pool size information
that leverages new developments in the field is presented and advocated for.
Finally, we discuss how the adoption of robust validation and selection
procedures can enhance confidence in constraint-based modeling as a whole and
ultimately facilitate more widespread use of FBA in biotechnology in
particular.


  We discuss the inadequacy of covariances/correlations and other measures in
L-2 as relative distance metrics. We propose a computationally simple heuristic
to transform a map based on standard principal component analysis (PCA) (when
the variables are asymptotically Gaussian) into an entropy-based map where
distances are based on mutual information (MI).


  The accuracy of classical force fields (FFs) has been shown to be limited for
the simulation of cation-protein systems despite their importance in
understanding the processes of life. Improvements can result from optimizing
the parameters of classical FFs or by extending the FF formulation by terms
describing charge transfer and polarization effects. In this work, we introduce
our implementation of the CTPOL model in OpenMM, which extends the classical
additive FF formula by adding charge transfer (CT) and polarization (POL).
Furthermore, we present an open-source parameterization tool, called FFAFFURR
that enables the (system specific) parameterization of OPLS-AA and CTPOL
models. The performance of our workflow was evaluated by its ability to
reproduce quantum chemistry energies and by molecular dynamics simulations of a
Zinc finger protein.


  Power law distributions are widely observed in chemical physics, geophysics,
biology, and beyond. The independent variable x of these distributions has an
obligatory lower bound and in many cases also an upper bound. Estimating these
bounds from sample data is notoriously difficult, with a recent method
involving O(N^3) operations, where N denotes sample size. Here I develop an
approach for estimating the lower and upper bounds that involves O(N)
operations. The approach centers on calculating the mean values, x_min and
x_max, of the smallest x and the largest x in N-point samples. A fit of x_min
or x_max as a function of N yields the estimate for the lower or upper bound.
Application to synthetic data demonstrates the accuracy and reliability of this
approach.


  Genetic information is encoded in a linear sequence of nucleotides,
represented by letters ranging from thousands to billions. Mutations refer to
changes in the DNA or RNA nucleotide sequence. Thus, mutation detection is
vital in all areas of biology and medicine. Careful monitoring of
virulence-enhancing mutations is essential. However, an enormous amount of
classical computing power is required to analyze genetic sequences of this
size. Inspired by human perception of vision and pixel representation of images
on quantum computers, we leverage these techniques to implement a pairwise
sequence analysis. The methodology has a potential advantage over classical
approaches and can be further applied to identify mutations and other
modifications in genetic sequences. We present a method to display and analyze
the similarity between two genome sequences on a quantum computer where a
similarity score is calculated to determine the similarity between nucleotides.


  Phylogenetic and discrete-trait evolutionary inference depend heavily on an
appropriate characterization of the underlying character substitution process.
In this paper, we present random-effects substitution models that extend common
continuous-time Markov chain models into a richer class of processes capable of
capturing a wider variety of substitution dynamics. As these random-effects
substitution models often require many more parameters than their usual
counterparts, inference can be both statistically and computationally
challenging. Thus, we also propose an efficient approach to compute an
approximation to the gradient of the data likelihood with respect to all
unknown substitution model parameters. We demonstrate that this approximate
gradient enables scaling of sampling-based inference, namely Bayesian inference
via Hamiltonian Monte Carlo, under random-effects substitution models across
large trees and state-spaces. Applied to a dataset of 583 SARS-CoV-2 sequences,
an HKY model with random-effects shows strong signals of nonreversibility in
the substitution process, and posterior predictive model checks clearly show
that it is a more adequate model than a reversible model. When analyzing the
pattern of phylogeographic spread of 1441 influenza A virus (H3N2) sequences
between 14 regions, a random-effects phylogeographic substitution model infers
that air travel volume adequately predicts almost all dispersal rates. A
random-effects state-dependent substitution model reveals no evidence for an
effect of arboreality on the swimming mode in the tree frog subfamily Hylinae.
Simulations reveal that random-effects substitution models can accommodate both
negligible and radical departures from the underlying base substitution model.
We show that our gradient-based inference approach is over an order of
magnitude more time efficient than conventional approaches.


  Scientists have been trying to identify all of the genes in the human genome
since the initial draft of the genome was published in 2001. Over the
intervening years, much progress has been made in identifying protein-coding
genes, and the estimated number has shrunk to fewer than 20,000, although the
number of distinct protein-coding isoforms has expanded dramatically. The
invention of high-throughput RNA sequencing and other technological
breakthroughs have led to an explosion in the number of reported non-coding RNA
genes, although most of them do not yet have any known function. A combination
of recent advances offers a path forward to identifying these functions and
towards eventually completing the human gene catalogue. However, much work
remains to be done before we have a universal annotation standard that includes
all medically significant genes, maintains their relationships with different
reference genomes, and describes clinically relevant genetic variants.


  Gene regulatory relationships can be abstracted as a gene regulatory network
(GRN), which plays a key role in characterizing complex cellular processes and
pathways. Recently, graph neural networks (GNNs), as a class of deep learning
models, have emerged as a useful tool to infer gene regulatory relationships
from gene expression data. However, deep learning models have been found to be
vulnerable to noise, which greatly hinders the adoption of deep learning in
constructing GRNs, because high noise is often unavoidable in the process of
gene expression measurement. Can we preferably prototype a robust GNN for
constructing GRNs? In this paper, we give a positive answer by proposing a
Quadratic Graph Attention Network (Q-GAT) with a dual attention mechanism. We
study the changes in the predictive accuracy of Q-GAT and 9 state-of-the-art
baselines by introducing different levels of adversarial perturbations.
Experiments in the E. coli and S. cerevisiae datasets suggest that Q-GAT
outperforms the state-of-the-art models in robustness. Lastly, we dissect why
Q-GAT is robust through the signal-to-noise ratio (SNR) and interpretability
analyses. The former informs that nonlinear aggregation of quadratic neurons
can amplify useful signals and suppress unwanted noise, thereby facilitating
robustness, while the latter reveals that Q-GAT can leverage more features in
prediction thanks to the dual attention mechanism, which endows Q-GAT with the
ability to confront adversarial perturbation. We have shared our code in
https://github.com/Minorway/Q-GAT_for_Robust_Construction_of_GRN for readers'
evaluation.


  Developing a mathematical understanding of autocatalysis in reaction networks
has both theoretical and practical implications. We review definitions of
autocatalytic networks and prove some properties for minimal autocatalytic
subnetworks (MASs). We show that it is possible to classify MASs in equivalence
classes, and develop mathematical results about their behavior. We also provide
linear-programming algorithms to exhaustively enumerate them and a scheme to
visualize their polyhedral geometry and combinatorics. We then define cluster
chemical reaction networks, a framework for coarse-graining real chemical
reactions with positive integer conservation laws. We find that the size of the
list of minimal autocatalytic subnetworks in a maximally connected cluster
chemical reaction network with one conservation law grows exponentially in the
number of species. We end our discussion with open questions concerning an
ecosystem of autocatalytic subnetworks and multidisciplinary opportunities for
future investigation.


  In contemporary society, social networks accelerate decision dynamics causing
a rapid switch of opinions in a number of fields, including the prevention of
infectious diseases by means of vaccines. This means that opinion dynamics can
nowadays be much faster than the spread of epidemics. Hence, we propose a
Susceptible-Infectious-Removed epidemic model coupled with an evolutionary
vaccination game embedding the public health system efforts to increase vaccine
uptake. This results in a global system ``epidemic model + evolutionary game''.
The epidemiological novelty of this work is that we assume that the switching
to the strategy ``pro vaccine'' depends on the incidence of the disease. As a
consequence of the above-mentioned accelerated decisions, the dynamics of the
system acts on two different scales: a fast scale for the vaccine decisions and
a slower scale for the spread of the disease. Another, and more methodological,
element of novelty is that we apply Geometrical Singular Perturbation Theory
(GSPT) to such a two-scale model and we then compare the geometric analysis
with the Quasi-Steady-State Approximation (QSSA) approach, showing a
criticality in the latter. Later, we apply the GSPT approach to the disease
prevalence-based model already studied in (Della Marca and d'Onofrio, Comm Nonl
Sci Num Sim, 2021) via the QSSA approach by considering medium-large values of
the strategy switching parameter.


  Background: The aim of this study is to use bioinformatics to discover the
biomarkers associated with patients with Pancreatic Ductal
Adenocarcinoma(PDAC).
  Material and Methods: GSE28735, GSE15471, and GSE62452 are gene microarray
datasets drived from the GEO database, included 153 PDAC samples and 145 normal
samples. By analyzing both Gene Ontology (GO) and the Kyoto Encyclopedia of
Genes (KEGG), for screening DEGs has provided information about their
biological function. Protein-protein interactions (PPI) of DEGs were analyzed
using the Search Tool for the Retrieval of Interacting Genes database (STRING)
and visualized by Cytoscape. UALCAN was also used to perform prognostic
analyses.
  Results: In Pancreatic ductal adenocarcinoma, we discovered 2264 upregulated
DEGs (uDEGs) and 723 downregulated DEGs (dDEGs). The Gene Ontology (GO)
indicated that extracellular matrix organization, extracellular structure
organization, collagen catabolic process and also the uDEGs are enriched in
focal adhesion, PI3K-Akt signaling pathway, ECM-receptor interaction by KEGG
analysis. The top ten hub genes were identified from the PPI network: COL1A1,
COL3A1, COL1A2, FN1, COL5A2, ITGA3, FBN1, MET, COL6A3, and BGN. These 10 hub
genes are highly upregulated in pancreatic ductal adenocarcinoma, according to
GEPIA research. The UALCAN prognostic analysis of the 10 hub genes revealed
that two of the elevated genes, ITGA3 and MET, significantly shortened PDAC
patient survival time. Pancreatic ductal adenocarcinoma is linked to Focal
adhesion, the PI3K-Akt signaling pathway, and ECM-receptor interaction,
according to Module analysis from the PPI network.
  Conclusion: This investigation identified hub genes and important signal
pathways, which adds to our knowledge of molecular mechanisms and could be
exploited as diagnostic and therapeutic biomarkers for Pancreatic ductal
adenocarcinoma.


  Noise is a ubiquitous feature of the physical world. As a result, the first
prerequisite of life is fault tolerance: maintaining integrity of state despite
external bombardment. Recent experimental advances have revealed that
biological systems achieve fault tolerance by implementing mathematically
intricate error-correcting codes and by organizing in a modular fashion that
physically separates functionally distinct subsystems. These elaborate
structures represent a vanishing volume in the massive genetic configuration
space. How is it possible that the primitive process of evolution, by which all
biological systems evolved, achieved such unusual results? In this work,
through experiments in Boolean networks, we show that the simultaneous presence
of error correction and modularity in biological systems is no coincidence.
Rather, it is a typical co-occurrence in noisy dynamic systems undergoing
evolution. From this, we deduce the principle of error correction enhanced
evolvability: systems possessing error-correcting codes are more effectively
improved by evolution than those without.


  Biomolecular communication demands that interactions between parts of a
molecular system act as scaffolds for message transmission. It also requires an
evolving and organized system of signs - a communicative agency - for creating
and transmitting meaning. Here I explore the need to dissect biomolecular
communication with retrodiction approaches that make claims about the past
given information that is available in the present. While the passage of time
restricts the explanatory power of retrodiction, the use of molecular structure
in biology offsets information erosion. This allows description of the gradual
evolutionary rise of structural and functional innovations in RNA and proteins.
The resulting chronologies can also describe the gradual rise of molecular
machines of increasing complexity and computation capabilities. For example,
the accretion of rRNA substructures and ribosomal proteins can be traced in
time and placed within a geological timescale. Phylogenetic, algorithmic and
theoretical-inspired accretion models can be reconciled into a congruent
evolutionary model. Remarkably, the time of origin of enzymes, functional RNA,
non-ribosomal peptide synthetase (NRPS) complexes, and ribosomes suggest they
gradually climbed Chomsky's hierarchy of formal grammars, supporting the
gradual complexification of machines and communication in molecular biology.
Future retrodiction approaches and in-depth exploration of theoretical models
of computation will need to confirm such evolutionary progression.


  Explainable and interpretable unsupervised machine learning helps understand
the underlying structure of data. We introduce an ensemble analysis of machine
learning models to consolidate their interpretation. Its application shows that
restricted Boltzmann machines compress consistently into a few bits the
information stored in a sequence of five amino acids at the start or end of
$\alpha$-helices or $\beta$-sheets. The weights learned by the machines reveal
unexpected properties of the amino acids and the secondary structure of
proteins: (i) His and Thr have a negligible contribution to the amphiphilic
pattern of $\alpha$-helices; (ii) there is a class of $\alpha$-helices
particularly rich in Ala at their end; (iii) Pro occupies most often slots
otherwise occupied by polar or charged amino acids, and its presence at the
start of helices is relevant; (iv) Glu and especially Asp on one side, and Val,
Leu, Iso, and Phe on the other, display the strongest tendency to mark
amphiphilic patterns, i.e., extreme values of an "effective hydrophobicity",
though they are not the most powerful (non) hydrophobic amino acids.


  The Markovian approach, which assumes exponentially distributed
interinfection times, is dominant in epidemic modeling. However, this
assumption is unrealistic as an individual's infectiousness depends on its
viral load and varies over time. In this paper, we present a
Susceptible-Infected-Recovered-Vaccinated-Susceptible epidemic model
incorporating non-Markovian infection processes. The model can be easily
adapted to accurately capture the generation time distributions of emerging
infectious diseases, which is essential for accurate epidemic prediction. We
observe noticeable variations in the transient behavior under different
infectiousness profiles and the same basic reproduction number R0. The
theoretical analyses show that only R0 and the mean immunity period of the
vaccinated individuals have an impact on the critical vaccination rate needed
to achieve herd immunity. A vaccination level at the critical vaccination rate
can ensure a very low incidence among the population in case of future
epidemics, regardless of the infectiousness profiles.


  Understanding the base pairing of an RNA sequence provides insight into its
molecular structure.By mining suboptimal sampling data, RNAprofiling 1.0
identifies the dominant helices in low-energy secondary structures as features,
organizes them into profiles which partition the Boltzmann sample, and
highlights key similarities/differences among the most informative, i.e.
selected, profiles in a graphical format. Version 2.0 enhances every step of
this approach. First, the featured substructures are expanded from helices to
stems. Second, profile selection includes low-frequency pairings similar to
featured ones. In conjunction, these updates extend the utility of the method
to sequences up to length 600, as evaluated over a sizable dataset. Third,
relationships are visualized in a decision tree which highlights the most
important structural differences. Finally, this cluster analysis is made
accessible to experimental researchers in a portable format as an interactive
webpage, permitting a much greater understanding of trade-offs among different
possible base pairing combinations.


  The molecular origins of proteins' functions are a combinatorial search
problem in the proteins' sequence space, which requires enormous resources to
solve. However, evolution has already solved this optimization problem for us,
leaving behind suboptimal solutions along the way. Comparing suboptimal
proteins along the evolutionary pathway, or ancestors, with more optimal modern
proteins can lead us to the exact molecular origins of a particular function.
In this paper, we study the long-standing question of the selectivity of
Imatinib, an anti-cancer kinase inhibitor drug. We study two related kinases,
Src and Abl, and four of their common ancestors, to which Imatinib has
significantly different affinities. Our results show that the orientation of
the N-lobe with respect to the C-lobe varies between the kinases along their
evolutionary pathway and is consistent with Imatinib's inhibition constants as
measured experimentally. The conformation of the DFG-motif (Asp-Phe-Gly) and
the structure of the P-loop also seem to have different stable conformations
along the evolutionary pathway, which is aligned with Imatinib's affinity.


  The high-throughput short-reads RNA-seq protocols often produce paired-end
reads, with the middle portion of the fragments being unsequenced. We explore
if the full-length fragments can be computationally reconstructed from the
sequenced two ends in the absence of the reference genome - a problem here we
refer to as de novo bridging. Solving this problem provides longer, more
informative RNA-seq reads, and benefits downstream RNA-seq analysis such as
transcript assembly, expression quantification, and splicing differential
analysis. However, de novo bridging is a challenging and complicated task owing
to alternative splicing, transcript noises, and sequencing errors. It remains
unclear if the data provides sufficient information for accurate bridging, let
alone efficient algorithms that determine the true bridges. Methods have been
proposed to bridge paired-end reads in the presence of reference genome (called
reference-based bridging), but the algorithms are far away from scaling for de
novo bridging as the underlying compacted de Bruijn graph(cdBG) used in the
latter task often contains millions of vertices and edges. We designed a new
truncated Dijkstra's algorithm for this problem, and proposed a novel algorithm
that reuses the shortest path tree to avoid running the truncated Dijkstra's
algorithm from scratch for all vertices for further speeding up. These
innovative techniques result in scalable algorithms that can bridge all
paired-end reads in a cdBG with millions of vertices. Our experiments showed
that paired-end RNA-seq reads can be accurately bridged to a large extent. The
resulting tool is freely available at
https://github.com/Shao-Group/rnabridge-denovo.


  Publicly available collections of drug-like molecules have grown to comprise
10s of billions of possibilities in recent history due to advances in chemical
synthesis. Traditional methods for identifying ``hit'' molecules from a large
collection of potential drug-like candidates have relied on biophysical theory
to compute approximations to the Gibbs free energy of the binding interaction
between the drug to its protein target. A major drawback of the approaches is
that they require exceptional computing capabilities to consider for even
relatively small collections of molecules.
  Hyperdimensional Computing (HDC) is a recently proposed learning paradigm
that is able to leverage low-precision binary vector arithmetic to build
efficient representations of the data that can be obtained without the need for
gradient-based optimization approaches that are required in many conventional
machine learning and deep learning approaches. This algorithmic simplicity
allows for acceleration in hardware that has been previously demonstrated for a
range of application areas. We consider existing HDC approaches for molecular
property classification and introduce two novel encoding algorithms that
leverage the extended connectivity fingerprint (ECFP) algorithm.
  We show that HDC-based inference methods are as much as 90 times more
efficient than more complex representative machine learning methods and achieve
an acceleration of nearly 9 orders of magnitude as compared to inference with
molecular docking. We demonstrate multiple approaches for the encoding of
molecular data for HDC and examine their relative performance on a range of
challenging molecular property prediction and drug-protein binding
classification tasks. Our work thus motivates further investigation into
molecular representation learning to develop ultra-efficient pre-screening
tools.


  Motivation: There is accumulating evidence showing the important roles of
bacteriophages (phages) in regulating the structure and functions of the
microbiome. However, lacking an easy-to-use and integrated phage analysis
software hampers microbiome-related research from incorporating phages in the
analysis. Results: In this work, we developed a web server, PhaBOX, which can
comprehensively identify and analyze phage contigs in metagenomic data. It
supports integrated phage analysis, including phage contig identification from
the metagenomic assembly, lifestyle prediction, taxonomic classification, and
host prediction. Instead of treating the algorithms as a black box, PhaBOX also
supports visualization of the essential features for making predictions. The
web server is designed with a user-friendly graphical interface that enables
both informatics-trained and non-specialist users to analyze phages in
microbiome data with ease. Availability: The web server of PhaBOX is available
via: https://phage.ee.cityu.edu.hk. The source code of PhaBOX is available at:
https://github.com/KennthShang/PhaBOX Contact: yannisun@cityu.edu.hk


  Climate change and the resulting increased drought frequencies pose
considerable threats to forest herb populations, especially when compounded by
additional environmental challenges. Specifically, habitat fragmentation may
disrupt climate adaptation and cause shifts in mating systems. To examine this,
we conducted a garden experiment with Primula elatior offspring from 24
populations across a climate and landscape fragmentation gradient. We evaluated
vegetative, regulatory, and reproductive traits under different soil moisture
regimes, assessing local adaptation and phenotypic plasticity. We also
conducted a field study in 60 populations along the same gradient to examine
potential breakdown of reciprocal herkogamy. Our results showed an evolutionary
shift from drought avoidance in southern populations to drought tolerance in
northern populations for large, connected populations. However, fragmentation
disrupted climate clines and adaptive responses to drought in key traits
related to growth, biomass allocation and water regulation. Our findings also
indicate the beginning of an evolutionary breakdown in reciprocal herkogamy.
These disruptions resulted in significantly reduced flowering investment,
especially in southern fragmented populations. These findings provide new
evidence of how habitat fragmentation disrupts climate adaptation and drought
tolerance in Primula elatior, emphasizing the need to account for habitat
fragmentation in conservation strategies to preserve resilient forest herb
populations amidst global changes.


  This work presents a model combining the simplest communicable and
non-communicable disease models. The latter is, by far, the leading cause of
sickness and death in the World, and introduces basal heterogeneity in
populations where communicable diseases evolve. The model can be interpreted as
a risk-structured model, another way of accounting for population
heterogeneity.
  Our results show that considering the non-communicable disease (in the end,
heterogeneous populations) allows the communicable disease to become endemic
even if the basic reproduction number is less than $1$. This feature is known
as subcritical bifurcation. Furthermore, ignoring the non-communicable disease
dynamics results in overestimating the reproduction number and, thus, giving
wrong information about the actual number of infected individuals. We calculate
sensitivity indices and derive interesting epidemic-control information.


  Drug overdose deaths continue to increase in the United States for all major
drug categories. Over the past two decades the total number of overdose
fatalities has increased more than five-fold; since 2013 the surge in overdose
rates is primarily driven by fentanyl and methamphetamines. Different drug
categories and factors such as age, gender, and ethnicity are associated with
different overdose mortality characteristics that may also change in time. For
example, the average age at death from a drug overdose has decreased from 1940
to 1990 while the overall mortality rate has steadily increased. To provide
insight into the population-level dynamics of drug-overdose mortality, we
develop an age-structured model for drug addiction. Using an augmented ensemble
Kalman filter (EnKF), we show through a simple example how our model can be
combined with synthetic observation data to estimate mortality rate and an
age-distribution parameter. Finally, we use an EnKF to combine our model with
observation data on overdose fatalities in the United States from 1999 to 2020
to forecast the evolution of overdose trends and estimate model parameters.


  Given the large size and complexity of most biochemical regulation and
signaling networks, there is a non-trivial relationship between the micro-level
logic of component interactions and the observed macro-dynamics. Here we
address this issue by formalizing the existing concept of pathway modules,
which are sequences of state updates that are guaranteed to occur (barring
outside interference) in the dynamics of automata networks after the
perturbation of a subset of driver nodes. We present a novel algorithm to
automatically extract pathway modules from networks and we characterize the
interactions that may take place between modules. This methodology uses only
the causal logic of individual node variables (micro-dynamics) without the need
to compute the dynamical landscape of the networks (macro-dynamics).
Specifically, we identify complex modules, which maximize pathway length and
require synergy between their components. This allows us to propose a new take
on dynamical modularity that partitions complex networks into causal pathways
of variables that are guaranteed to transition to specific states given a
perturbation to a set of driver nodes. Thus, the same node variable can take
part in distinct modules depending on the state it takes. Our measure of
dynamical modularity of a network is then inversely proportional to the overlap
among complex modules and maximal when complex modules are completely
decouplable from one another in the network dynamics. We estimate dynamical
modularity for several genetic regulatory networks, including the Drosophila
melanogaster segment-polarity network. We discuss how identifying complex
modules and the dynamical modularity portrait of networks explains the
macro-dynamics of biological networks, such as uncovering the (more or less)
decouplable building blocks of emergent computation (or collective behavior) in
biochemical regulation and signaling.


  Protein language models (pLMs), pre-trained via causal language modeling on
protein sequences, have been a promising tool for protein sequence design. In
real-world protein engineering, there are many cases where the amino acids in
the middle of a protein sequence are optimized while maintaining other
residues. Unfortunately, because of the left-to-right nature of pLMs, existing
pLMs modify suffix residues by prompting prefix residues, which are
insufficient for the infilling task that considers the whole surrounding
context. To find the more effective pLMs for protein engineering, we design a
new benchmark, Secondary structureE InFilling rEcoveRy, SEIFER, which
approximates infilling sequence design scenarios. With the evaluation of
existing models on the benchmark, we reveal the weakness of existing language
models and show that language models trained via fill-in-middle transformation,
called ProtFIM, are more appropriate for protein engineering. Also, we prove
that ProtFIM generates protein sequences with decent protein representations
through exhaustive experiments and visualizations.