  Antibodies are versatile proteins that can bind to pathogens and provide
effective protection for human body. Recently, deep learning-based
computational antibody design has attracted popular attention since it
automatically mines the antibody patterns from data that could be complementary
to human experiences. However, the computational methods heavily rely on
high-quality antibody structure data, which is quite limited. Besides, the
complementarity-determining region (CDR), which is the key component of an
antibody that determines the specificity and binding affinity, is highly
variable and hard to predict. Therefore, the data limitation issue further
raises the difficulty of CDR generation for antibodies. Fortunately, there
exists a large amount of sequence data of antibodies that can help model the
CDR and alleviate the reliance on structure data. By witnessing the success of
pre-training models for protein modeling, in this paper, we develop the
antibody pre-training language model and incorporate it into the
(antigen-specific) antibody design model in a systemic way. Specifically, we
first pre-train an antibody language model based on the sequence data, then
propose a one-shot way for sequence and structure generation of CDR to avoid
the heavy cost and error propagation from an autoregressive manner, and finally
leverage the pre-trained antibody model for the antigen-specific antibody
generation model with some carefully designed modules. Through various
experiments, we show that our method achieves superior performances over
previous baselines on different tasks, such as sequence and structure
generation and antigen-binding CDR-H3 design.


  Animals use various processes to inform themselves about their environment
and make decisions about how to move and form their territory. In some cases,
populations inform themselves of competing groups through observations at
distances, scent markings, or memories of locations where an individual has
encountered competing populations. As the process of gathering this information
is inherently nonlocal, mechanistic models that include nonlocal terms have
been proposed to investigate the movement of species. Naturally, these models
present analytical and computational challenges. In this work we study a
multi-species model with nonlocal advection. We introduce an efficient
numerical scheme using spectral methods to compute solutions of a nonlocal
reaction-advection-diffusion system for a large number of interacting species.
Moreover, we investigate the effects that the parameters and interaction
potentials have on the population densities. Finally, we propose a method using
maximum likelihood estimation to determine the most important factors driving
species' movements and test this method using synthetic data.


  Genotype-to-phenotype mappings translate genotypic variations such as
mutations into phenotypic changes. Neutrality is the observation that some
mutations do not lead to phenotypic changes. Studying the search trajectories
in genotypic and phenotypic spaces, especially through neutral mutations, helps
us to better understand the progression of evolution and its algorithmic
behaviour. In this study, we visualise the search trajectories of a genetic
programming system as graph-based models, where nodes are genotypes/phenotypes
and edges represent their mutational transitions. We also quantitatively
measure the characteristics of phenotypes including their genotypic abundance
(the requirement for neutrality) and Kolmogorov complexity. We connect these
quantified metrics with search trajectory visualisations, and find that more
complex phenotypes are under-represented by fewer genotypes and are harder for
evolution to discover. Less complex phenotypes, on the other hand, are
over-represented by genotypes, are easier to find, and frequently serve as
stepping-stones for evolution.


  All cognitive agents are composite beings. Specifically, complex living
agents consist of cells, which are themselves competent sub-agents navigating
physiological and metabolic spaces. Behavior science, evolutionary
developmental biology, and the field of machine intelligence all seek an answer
to the scaling of biological cognition: what evolutionary dynamics enable
individual cells to integrate their activities to result in the emergence of a
novel, higher-level intelligence that has goals and competencies that belong to
it and not to its parts? Here, we report the results of simulations based on
the TAME framework, which proposes that evolution pivoted the collective
intelligence of cells during morphogenesis of the body into traditional
behavioral intelligence by scaling up the goal states at the center of
homeostatic processes. We tested the hypothesis that a minimal evolutionary
framework is sufficient for small, low-level setpoints of metabolic homeostasis
in cells to scale up into collectives (tissues) which solve a problem in
morphospace: the organization of a body-wide positional information axis (the
classic French Flag problem). We found that these emergent morphogenetic agents
exhibit a number of predicted features, including the use of stress propagation
dynamics to achieve its target morphology as well as the ability to recover
from perturbation (robustness) and long-term stability (even though neither of
these was directly selected for). Moreover we observed unexpected behavior of
sudden remodeling long after the system stabilizes. We tested this prediction
in a biological system - regenerating planaria - and observed a very similar
phenomenon. We propose that this system is a first step toward a quantitative
understanding of how evolution scales minimal goal-directed behavior
(homeostatic loops) into higher-level problem-solving agents in morphogenetic
and other spaces.


  The study of infectious disease has been of interest to ecologists since
long. The initiation of epidemic and the long term disease dynamics are largely
influenced by the nature of the underlying consumer (host)-resource dynamics.
Ecological traits of such systems may be often modulated by toxins released in
the environment due to ongoing anthropogenic activities. This, in addition to
toxin-mediated alteration of epidemiological traits, has a significant impact
on disease progression in ecosystems which is quite less studied. In order to
address this, we consider a mathematical model of disease transmission in
consumer population where multiple traits are affected by environmental toxins.
Long term dynamics show that the level of environmental toxin determines
disease persistence, and increasing toxin may even eradicate the disease in
certain circumstances. Furthermore, our results demonstrate bistability between
different ecosystem states and the possibility of an abrupt transition from
disease-free coexistence to disease-induced extinction of consumers. Overall
the results from this study will help us gain fundamental insights into disease
propagation in natural ecosystems in the face of present anthropogenic changes.


  Chaperonins are biological nanomachines that help newly translated proteins
to fold by rescuing them from kinetically trapped misfolded states. Protein
folding assistance by the chaperonin machinery is obligatory in vivo for a
subset of proteins in the bacterial proteome. Chaperonins are large oligomeric
complexes, with unusual seven fold symmetry (group I) or eight/nine fold
symmetry (group II), that form double-ring constructs, enclosing a central
folding chamber. Dramatic large-scale conformational changes, that take place
during ATP-driven cycles, allow chaperonins to bind misfolded proteins,
encapsulate them into the expanded cavity and release them back into the
cellular environment, regardless of whether they are folded or not. The theory
associated with the iterative annealing mechanism, which incorporated the
conformational free energy landscape description of protein folding,
\textit{quantitatively} explains most, if not all, the available data.
Misfolded conformations are associated with low energy minima in a rugged
energy landscape. Random disruptions of these low energy conformations result
in higher free energy, less folded, conformations that can stochastically
partition into the native state. Group I chaperonins (GroEL homologues in
eubacteria and endosymbiotic organelles), recognize a large number of misfolded
proteins non-specifically and operate through highly coordinated cooperative
motions. By contrast, the less well understood group II chaperonins (CCT in
Eukarya and thermosome/TF55 in Archaea), assist a selected set of substrate
proteins. Chaperonins are implicated in bacterial infection, autoimmune
disease, as well as protein aggregation and degradation diseases. Understanding
the chaperonin mechanism and their substrates is important not only for the
fundamental aspect of cellular protein folding, but also for effective
therapeutic strategies.


  Early estimates of the transmission properties of a newly emerged pathogen
are critical to an effective public health response, and are often based on
limited outbreak data. Here, we use simulations to investigate a potential
source of bias in such estimates, arising from correlations between the viral
load of cases in transmission chains. We show that this mechanism can affect
estimates of fundamental transmission properties characterising the spread of a
virus. Our computational model simulates a disease transmission mechanism in
which the viral load of the infector at the time of transmission influences the
infectiousness of the infectee. These correlations in transmission pairs
produce a population-level decoherence process during which the distributions
of initial viral loads in each subsequent generation converge to a steady
state. We find that outbreaks arising from index cases with low initial viral
loads give rise to early estimates of transmission properties that are subject
to large biases. These findings demonstrate the potential for bias arising from
transmission mechanics to affect estimates of the transmission properties of
newly emerged viruses.


  A deterministic model with testing of infected individuals has been proposed
to investigate the potential consequences of the impact of testing strategy.
The model exhibits global dynamics concerning the disease-free and a unique
endemic equilibrium depending on the basic reproduction number when the
recruitment of infected individuals is zero; otherwise, the model does not have
a disease-free equilibrium, and disease never dies out in the community. Model
parameters have been estimated using the maximum likelihood method with respect
to the data of early COVID-19 outbreak in India. The practical identifiability
analysis shows that the model parameters are estimated uniquely. The
consequences of the testing rate for the weekly new cases of early COVID-19
data in India tell that if the testing rate is increased by 20% and 30% from
its baseline value, the weekly new cases at the peak are decreased by 37.63%
and 52.90%; and it also delayed the peak time by four and fourteen weeks,
respectively. Similar findings are obtained for the testing efficacy that if it
is increased by 12.67% from its baseline value, the weekly new cases at the
peak are decreased by 59.05% and delayed the peak by 15 weeks. Therefore, a
higher testing rate and efficacy reduce the disease burden by tumbling the new
cases, representing a real scenario. It is also obtained that the testing rate
and efficacy reduce the epidemic's severity by increasing the final size of the
susceptible population. The testing rate is found more significant if testing
efficacy is high. Global sensitivity analysis using partial rank correlation
coefficients (PRCCs) and Latin hypercube sampling (LHS) determine the key
parameters that must be targeted to worsen/contain the epidemic.


  The SIRS model with constant vaccination and immunity waning rates is well
known to show a transition from a disease-free to an endemic equilibrium as the
basic reproduction number $r_0$ is raised above threshold. It is shown that
this model maps to Hethcote's classic endemic model originally published in
1973. In this way one obtains unifying formulas for a whole class of models
showing endemic bifurcation. In particular, if the vaccination rate is smaller
than the recovery rate and $r_- < r_0 < r_+$ for certain upper and lower bounds
$r_\pm$, then trajectories spiral into the endemic equilibrium via damped
infection waves. Latest data of the SARS-CoV-2 Omicron variant suggest that
according to this simplified model continuous vaccination programs will not be
capable to escape the oscillating endemic phase. However, in view of the strong
damping factors predicted by the model, in reality these oscillations will
certainly be overruled by time-dependent contact behaviors.


  Comparing how different populations have suffered under COVID-19 is a core
part of ongoing investigations into how public policy and social inequalities
influence the number of and severity of COVID-19 cases. But COVID-19 incidence
can vary multifold from one subpopulation to another, including between
neighborhoods of the same city, making comparisons of case rates deceptive. At
the same time, although epidemiological heterogeneities are increasingly
well-represented in mathematical models of disease spread, fitting these models
to real data on case numbers presents a tremendous challenge, as does
interpreting the models to answer questions such as: Which public health
policies achieve the best outcomes? Which social sacrifices are most worth
making? Here we compare COVID-19 case-curves between different US states, by
clustering case surges between March 2020 and March 2021 into groups with
similar dynamics. We advance the hypothesis that each surge is driven by a
subpopulation of COVID-19 contacting individuals, and make detecting the size
of that population a step within our clustering algorithm. Clustering reveals
that case trajectories in each state conform to one of a small number (4-6) of
archetypal dynamics. Our results suggest that while the spread of COVID-19 in
different states is heterogeneous, there are underlying universalities in the
spread of the disease that may yet be predictable by models with reduced
mathematical complexity. These universalities also prove to be surprisingly
robust to school closures, which we choose as a common, but high social cost,
public health measure.


  SIRS epidemic models assume that individual immunity (from infection and
vaccination) wanes in one big leap, from complete immunity to complete
susceptibility. For many diseases immunity on the contrary wanes gradually,
something that's become even more evident during COVID-19 pandemic where also
recently infected have a reinfection risk, and where booster vaccines are given
to increase immunity. This paper considers an epidemic model allowing for such
gradual waning of immunity (either linear or exponential waning) thereby
extending SIRS epidemics, and also incorporates vaccination. The two versions
for gradual waning of immunity are compared with the classic SIRS epidemic,
where the three models are calibrated by having the same \emph{average
cumulative immunity}. All models are shown to have identical basic reproduction
number $R_0$. However, if no prevention is put in place, the exponential waning
model has highest prevalence and the classic SIRS model has lowest. Similarly,
the amount of vaccine supply needed to reach and maintain herd immunity is
highest for the model with exponential decay of immunity and lowest for the
classic SIRS model. consequently, if truth lies close to exponential (or
linear) decay of immunity, expressions based on the SIRS epidemic will
underestimate the endemic level and the critical vaccine supply will not be
sufficient to reach and maintain herd immunity. For parameter choices fitting
to COVID-19, the critical amount of vaccine supply is about 50% higher if
immunity wanes linearly, and more than 150% higher when immunity wanes
exponentially, as compared to the classic SIRS epidemic model.


  We study the stochastic spatial Lotka-Volterra (LV) model for predator-prey
interaction subject to a periodically varying carrying capacity. The LV model
with on-site lattice occupation restrictions that represent finite food
resources for the prey exhibits a continuous active-to-absorbing phase
transition. The active phase is sustained by spatio-temporal patterns in the
form of pursuit and evasion waves. Monte Carlo simulations on a two-dimensional
lattice are utilized to investigate the effect of seasonal variations of the
environment on species coexistence. The results of our simulations are also
compared to a mean-field analysis. We find that the parameter region of
predator and prey coexistence is enlarged relative to the stationary situation
when the carrying capacity varies periodically. The stationary regime of our
periodically varying LV system shows qualitative agreement between the
stochastic model and the mean-field approximation. However, under periodic
carrying capacity switching environments, the mean-field rate equations predict
period-doubling scenarios that are washed out by internal reaction noise in the
stochastic lattice model. Utilizing visual representations of the lattice
simulations and dynamical correlation functions, we study how the pursuit and
evasion waves are affected by ensuing resonance effects. Correlation function
measurements indicate a time delay in the response of the system to sudden
changes in the environment. Resonance features are observed in our simulations
that cause prolonged persistent spatial correlations. Different effective
static environments are explored in the extreme limits of fast- and slow
periodic switching. The analysis of the mean-field equations in the
fast-switching regime enables a semi-quantitative description of the stationary
state.


  Strongly opinionated minorities can have a dramatic impact on the opinion
dynamics of a large population. Two factions of inflexible minorities,
polarised into two competing opinions, could lead the entire population to
persistent indecision. Equivalently, populations can remain undecided when
individuals sporadically change their opinion based on individual information
rather than social information. Our analysis compares the cross-inhibition
model with the voter model for decisions between equally good alternatives, and
with the weighted voter model for decisions among alternatives characterised by
different qualities. Here we show that cross-inhibition, differently from the
other two models, is a simple mechanism, ubiquitous in collective biological
systems, that allows the population to reach a stable majority for one
alternative even in the presence of asocial behaviour. The results predicted by
the mean-field models are confirmed by experiments with swarms of 100 locally
interacting robots. This work suggests an answer to the longstanding question
of why inhibitory signals are widespread in natural systems of collective
decision making, and, at the same time, it proposes an efficient mechanism for
designing resilient swarms of minimalistic robots.


  Signed networks, i.e., networks with positive and negative edges, commonly
arise in various domains from social media to epidemiology. Modeling signed
networks has many practical applications, including the creation of synthetic
data sets for experiments where obtaining real data is difficult. Influential
prior works proposed and studied various graph topology models, as well as the
problem of selecting the most fitting model for different application domains.
However, these topology models are typically unsigned.
  In this work, we pose a novel Maximum-Likelihood-based optimization problem
for modeling signed networks given their topology and showcase it in the
context of gene regulation. Regulatory interactions of genes play a key role in
organism development, and when broken can lead to serious organism
abnormalities and diseases. Our contributions are threefold: First, we design a
new class of signage models for a given topology. Based on the parameter
setting, we discuss its biological interpretations for gene regulatory networks
(GRNs). Second, we design algorithms computing the Maximum Likelihood --
depending on the parameter setting, our algorithms range from closed-form
expressions to MCMC sampling. Third, we evaluated the results of our algorithms
on synthetic datasets and real-world large GRNs. Our work can lead to the
prediction of unknown gene regulations, the generation of biological
hypotheses, and realistic GRN benchmark datasets.


  Deep Learning and big data have shown tremendous success in bioinformatics
and computational biology in recent years; artificial intelligence methods have
also significantly contributed in the task of protein function classification.
This review paper analyzes the recent developments in approaches for the task
of predicting protein function using deep learning. We explain the importance
of determining the protein function and why automating the following task is
crucial. Then, after reviewing the widely used deep learning techniques for
this task, we continue our review and highlight the emergence of the modern
State of The Art (SOTA) deep learning models which have achieved groundbreaking
results in the field of computer vision, natural language processing and
multi-modal learning in the last few years. We hope that this review will
provide a broad view of the current role and advances of deep learning in
biological sciences, especially in predicting protein function tasks and
encourage new researchers to contribute to this area.


  Accurate genome sequencing can improve our understanding of biology and the
genetic basis of disease. The standard approach for generating DNA sequences
from PacBio instruments relies on HMM-based models. Here, we introduce
Distilled DeepConsensus - a distilled transformer-encoder model for sequence
correction, which improves upon the HMM-based methods with runtime constraints
in mind. Distilled DeepConsensus is 1.3x faster and 1.5x smaller than its
larger counterpart while improving the yield of high quality reads (Q30) over
the HMM-based method by 1.69x (vs. 1.73x for larger model). With improved
accuracy of genomic sequences, Distilled DeepConsensus improves downstream
applications of genomic sequence analysis such as reducing variant calling
errors by 39% (34% for larger model) and improving genome assembly quality by
3.8% (4.2% for larger model). We show that the representations learned by
Distilled DeepConsensus are similar between faster and slower models.


  The Stag-hunt game is a prototype for social contracts. Adopting a new and
better social contract is usually challenging because the current one is
already widely adopted and stable due to deviants' sanctions. Thus, how does a
population shift from the current social contract to a better one? In other
words, how can a social system leave a local social optimum configuration to
achieve an optimum global state? Here, we investigate the effect of promoting
diversity on the evolution of social contracts. We considered group-structured
populations where individuals play the Stag-hunt game in all groups. We model
the diversity incentive as a Snow-drift game played in a single focus group
where the individual is more prone to adopt a deviant norm. We show that
moderate diversity incentives can change the system dynamics, leading the whole
population to move from the locally optimal social normal to the globally
optimal one. Thus, an initial fraction of adopters of the new norm can drive
the system toward the new social optimum norm. After the new social contract
becomes the new equilibrium, it remains stable even without the incentive. The
results are obtained using Monte Carlo simulations and analytical approximation
methods.


  The participation of astrocytes in brain computation was formally
hypothesized in 1992, coinciding with the discovery that these glial cells
display a complex form of Ca2+ excitability. This fostered conceptual advances
centered on the notion of reciprocal interactions between neurons and
astrocytes, which permitted a critical leap forward in uncovering many roles of
astrocytes in brain circuits, and signaled the rise of a major new force in
neuroscience: that of glial biology. In the past decade, a multitude of
unconventional and disparate functions of astrocytes have been documented that
are not predicted by these canonical models and that are challenging to piece
together into a holistic and parsimonious picture. This highlights a disconnect
between the rapidly evolving field of astrocyte biology and the conceptual
frameworks guiding it, and emphasizes the need for a careful reconsideration of
how we theorize the functional position of astrocytes in brain circuitry. Here,
we propose a unifying, highly transferable, data-driven, and
computationally-relevant conceptual framework for astrocyte biology, which we
coin contextual guidance. It describes astrocytes as contextual gates that
decode multiple environmental factors to shape neural circuitry in an adaptive,
state-dependent fashion. This paradigm is organically inclusive of all
fundamental features of astrocytes, many of which have remained unaccounted for
in previous theories. We find that this new concept provides an intuitive and
powerful theoretical space to improve our understanding of brain function and
computational models thereof across scales because it depicts astrocytes as a
hub for circumstantial inputs into relevant specialized circuits that permits
adaptive behaviors at the network and organism level.


  Despite being self-supervised, protein language models have shown remarkable
performance in fundamental biological tasks such as predicting impact of
genetic variation on protein structure and function. The effectiveness of these
models on diverse set of tasks suggests that they learn meaningful
representations of fitness landscape that can be useful for downstream clinical
applications. Here, we interrogate the use of these language models in
characterizing known pathogenic mutations in curated, medically actionable
genes through an exhaustive search of putative compensatory mutations on each
variant's genetic background. Systematic analysis of the predicted effects of
these compensatory mutations reveal unappreciated structural features of
proteins that are missed by other structure predictors like AlphaFold. While
deep mutational scan experiments provide an unbiased estimate of the mutational
landscape, we encourage the community to generate and curate rescue mutation
experiments to inform the design of more sophisticated co-masking strategies
and leverage large language models more effectively for downstream clinical
prediction tasks.


  The ability to design and optimize biological sequences with specific
functionalities would unlock enormous value in technology and healthcare. In
recent years, machine learning-guided sequence design has progressed this goal
significantly, though validating designed sequences in the lab or clinic takes
many months and substantial labor. It is therefore valuable to assess the
likelihood that a designed set contains sequences of the desired quality (which
often lies outside the label distribution in our training data) before
committing resources to an experiment. Forecasting, a prominent concept in many
domains where feedback can be delayed (e.g. elections), has not been used or
studied in the context of sequence design. Here we propose a method to guide
decision-making that forecasts the performance of high-throughput libraries
(e.g. containing $10^5$ unique variants) based on estimates provided by models,
providing a posterior for the distribution of labels in the library. We show
that our method outperforms baselines that naively use model scores to estimate
library performance, which are the only tool available today for this purpose.


  Switch-like motifs are among the basic building blocks of biochemical
networks. A common motif that can serve as an ultrasensitive switch consists of
two enzymes acting antagonistically on a substrate, one making and the other
removing a covalent modification. To work as a switch, such covalent
modification cycles must be held out of thermodynamic equilibrium by continuous
expenditure of energy. Here, we exploit the linear framework for timescale
separation to establish tight bounds on the performance of any
covalent-modification switch, in terms of the chemical potential difference
driving the cycle. The bounds apply to arbitrary enzyme mechanisms, not just
Michaelis-Menten, with arbitrary rate constants, and thereby reflect
fundamental physical constraints on covalent switching.


  Resolving the structural variability of proteins is often key to
understanding the structure-function relationship of those macromolecular
machines. Single particle analysis using Cryogenic electron microscopy
(CryoEM), combined with machine learning algorithms, provides a way to reveal
the dynamics within the protein system from noisy micrographs. Here, we
introduce an improved computational method that uses Gaussian mixture models
for protein structure representation and deep neural networks for conformation
space embedding. By integrating information from molecular models into the
heterogeneity analysis, we can resolve complex protein conformational changes
at near atomic resolution and present the results in a more interpretable form.


  The SARS-CoV-2 coronavirus is the cause of the COVID-19 disease in humans.
Like many coronaviruses, it can adapt to different hosts and evolve into
different lineages. It is well-known that the major SARS-CoV-2 lineages are
characterized by mutations that happen predominantly in the spike protein.
Understanding the spike protein structure and how it can be perturbed is vital
for understanding and determining if a lineage is of concern. These are crucial
to identifying and controlling current outbreaks and preventing future
pandemics. Machine learning (ML) methods are a viable solution to this effort,
given the volume of available sequencing data, much of which is unaligned or
even unassembled. However, such ML methods require fixed-length numerical
feature vectors in Euclidean space to be applicable. Similarly, euclidean space
is not considered the best choice when working with the classification and
clustering tasks for biological sequences. For this purpose, we design a method
that converts the protein (spike) sequences into the sequence similarity
network (SSN). We can then use SSN as an input for the classical algorithms
from the graph mining domain for the typical tasks such as classification and
clustering to understand the data. We show that the proposed alignment-free
method is able to outperform the current SOTA method in terms of clustering
results. Similarly, we are able to achieve higher classification accuracy using
well-known Node2Vec-based embedding compared to other baseline embedding
approaches.


  The DNA duplex may be locally strongly bent in complexes with proteins, for
example, with polymerases or in a nucleosome. At such bends, the DNA helix is
locally in the non-canonical forms A (with a narrow major groove and a large
amount of north sugars) or C (with a narrow minor groove and a large share of
BII phosphates). To model the formation of such complexes by molecular dynamics
methods, the force field is required to reproduce these conformational
transitions for a naked DNA. We analyzed the available experimental data on the
B-C and B-A transitions under the conditions easily implemented in modeling: in
an aqueous NaCl solution. We selected six DNA duplexes which conformations at
different salt concentrations are known reliably enough. At low salt
concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly
shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a
strong and salt concentration dependent bias toward the A-form. The polymers
poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt
concentrations. The reproduction of the behavior of these oligomers can serve
as a test for the balance of interactions between the base stacking and the
conformational flexibility of the sugar-phosphate backbone in a DNA force
field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids,
and we failed to reproduce the experiment. In all the force fields, the salt
concentration dependence is very weak. The known B-philicity of the AMBER force
field proved to result from the B-philicity of its excessively strong base
stacking. In the CHARMM force field, the B-form is a result of a fragile
balance between the A-philic base stacking (especially for G:C pairs) and the
C-philic backbone. Finally, we analyzed some recent simulations of the LacI-,
SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force
field.


  Macromolecules change their shape (conformation) in the process of carrying
out their functions. The imaging by cryo-electron microscopy of rapidly-frozen,
individual copies of macromolecules (single particles) is a powerful and
general approach to understanding the motions and energy landscapes of
macromolecules. Widely-used computational methods already allow the recovery of
a few distinct conformations from heterogeneous single-particle samples, but
the treatment of complex forms of heterogeneity such as the continuum of
possible transitory states and flexible regions remains largely an open
problem. In recent years there has been a surge of new approaches for treating
the more general problem of continuous heterogeneity. This paper surveys the
current state of the art in this area.


  We consider reaction networks that admit a singular perturbation reduction in
a certain parameter range. The focus of this paper is on deriving "small
parameters" (briefly for small perturbation parameters), to gauge the accuracy
of the reduction, in a manner that is consistent, amenable to computation and
permits an interpretation in chemical or biochemical terms. Our work is based
on local timescale estimates via ratios of the real parts of eigenvalues of the
Jacobian near critical manifolds. This approach modifies the one introduced by
Segel and Slemrod and is familiar from computational singular perturbation
theory. While parameters derived by this method cannot provide universal
quantitative estimates for the accuracy of a reduction, they represent a
critical first step toward this end. Working directly with eigenvalues is
generally unfeasible, and at best cumbersome. Therefore we focus on the
coefficients of the characteristic polynomial to derive parameters, and relate
them to timescales. Thus, we obtain distinguished parameters for systems of
arbitrary dimension, with particular emphasis on reduction to dimension one. As
a first application, we discuss the Michaelis--Menten reaction mechanism system
in various settings, with new and perhaps surprising results. We proceed to
investigate more complex enzyme catalyzed reaction mechanisms (uncompetitive,
competitive inhibition and cooperativity) of dimension three, with reductions
to dimension one and two. The distinguished parameters we derive for these
three-dimensional systems are new. In fact, no rigorous derivation of small
parameters seems to exist in the literature so far. Numerical simulations are
included to illustrate the efficacy of the parameters obtained, but also to
show that certain limitations must be observed.


  We study Covid-19 spreading dynamics underlying 84 curves of daily Covid-19
infection rates pertaining to 84 districts belonging to the largest seven
cities in Taiwan during her pristine surge period. Our computational
developments begin with selecting and extracting 18 features from each smoothed
district-specific curve. This step of computing effort allows unstructured data
to be converted into structured data, with which we then demonstrate asymmetric
growth and decline dynamics among all involved curves. Specifically, based on
Theoretical Information measurements of conditional entropy and mutual
information, we compute major factors of order-1 and order-2 that reveal
significant effects on affecting the curves' peak value and curvature at peak,
which are two essential features characterizing all the curves. Further, we
investigate and demonstrate major factors determining the geographic and
social-economic induced behavioral effects by encoding each of these 84
districts with two binary characteristics: North-vs-South and
Unban-vs-suburban. Furthermore, based on this data-driven knowledge on the
district scale, we go on to study fine-scale behavioral effects on infectious
disease spreading through similarity among 96 age-group-specific curves of
daily infection rate within 12 urban districts of Taipei and 12 suburban
districts of New Taipei City, which counts for almost one-quarter of the island
nation's total population. We conclude that human living, traveling, and
working behaviors do implicitly affect the spreading dynamics of Covid-19
across Taiwan profoundly.


  We modelled emergence and spread of the Omicron variant of SARS-CoV-2 in
Australia between December 2021 and June 2022. This pandemic stage exhibited a
diverse epidemiological profile with emergence of co-circulating sub-lineages
of Omicron, further complicated by differences in social distancing behaviour
which varied over time. Our study delineated distinct phases of the
Omicron-associated pandemic stage, and retrospectively quantified the adoption
of social distancing measures, fluctuating over different time periods in
response to the observable incidence dynamics. We also modelled the
corresponding disease burden, in terms of hospitalisations, intensive care unit
occupancy, and mortality. Supported by good agreement between simulated and
actual health data, our study revealed that the nonlinear dynamics observed in
the daily incidence and disease burden were determined not only by introduction
of sub-lineages of Omicron, but also by the fluctuating adoption of social
distancing measures. Our high-resolution model can be used in design and
evaluation of public health interventions during future crises.


  Generating molecules that bind to specific proteins is an important but
challenging task in drug discovery. Previous works usually generate atoms in an
auto-regressive way, where element types and 3D coordinates of atoms are
generated one by one. However, in real-world molecular systems, the
interactions among atoms in an entire molecule are global, leading to the
energy function pair-coupled among atoms. With such energy-based consideration,
the modeling of probability should be based on joint distributions, rather than
sequentially conditional ones. Thus, the unnatural sequentially auto-regressive
modeling of molecule generation is likely to violate the physical rules, thus
resulting in poor properties of the generated molecules. In this work, a
generative diffusion model for molecular 3D structures based on target proteins
as contextual constraints is established, at a full-atom level in a
non-autoregressive way. Given a designated 3D protein binding site, our model
learns the generative process that denoises both element types and 3D
coordinates of an entire molecule, with an equivariant network. Experimentally,
the proposed method shows competitive performance compared with prevailing
works in terms of high affinity with proteins and appropriate molecule sizes as
well as other drug properties such as drug-likeness of the generated molecules.


  Motivated by the question of optimal vaccine allocation strategies in
heterogeneous population for epidemic models, we study various properties of
the \emph{effective reproduction number}. In the simplest case, given a fixed,
non-negative matrix $K$, this corresponds mathematically to the study of the
spectral radius $R_e(\eta)$ of the matrix product $\mathrm{Diag}(\eta)K$, as a
function of $\eta\in\mathbb{R}_+^n$. The matrix $K$ and the vector $\eta$ can
be interpreted as a next-generation operator and a vaccination strategy. This
can be generalized in an infinite dimensional case where the matrix $K$ is
replaced by a positive integral compact operator, which is composed with a
multiplication by a non-negative function $\eta$.
  We give sufficient conditions for the function $R_e$ to be convex or a
concave. Eventually, we provide equivalence properties on models which ensure
that the function $R_e$ is unchanged.


  Phylogenetic networks provide a means of describing the evolutionary history
of sets of species believed to have undergone hybridization or gene flow during
their evolution. The mutation process for a set of such species can be modeled
as a Markov process on a phylogenetic network. Previous work has shown that a
site-pattern probability distributions from a Jukes-Cantor phylogenetic network
model must satisfy certain algebraic invariants. As a corollary, aspects of the
phylogenetic network are theoretically identifiable from site-pattern
frequencies. In practice, because of the probabilistic nature of sequence
evolution, the phylogenetic network invariants will rarely be satisfied, even
for data generated under the model. Thus, using network invariants for
inferring phylogenetic networks requires some means of interpreting the
residuals, or deviations from zero, when observed site-pattern frequencies are
substituted into the invariants. In this work, we propose a method of utilizing
invariant residuals and support vector machines to infer 4-leaf level-one
phylogenetic networks, from which larger networks can be reconstructed. Given
data for a set of species, the support vector machine is first trained on model
data to learn the patterns of residuals corresponding to different network
structures to classify the network that produced the data. We demonstrate the
performance of our method on simulated data from the specified model and
primate data.


  The 6th wave of COVID-19 in Tokyo continued for the longest period of
infection (about 190 days from late Nov. 2021), and the 7th wave, which
occurred in mid-May 2022, was the largest wave ever (cumulative 1.7 million
people). In order to elucidate their factors, the infection wave was analyzed
by using the Avrami equation. The main component of the 6th wave was formed by
the coupling of increased human interaction due to the New Year holidays and
the invasion of the new virus variant Omicron BA.1. After that, side waves were
formed by the coupling of the invasion of the new virus variant Omicron BA.2
and the human interaction in the consecutive holidays in February, March, and
May. These side waves caused the 6th wave not to converge for a long time. The
outbreak of the main component of the 7th wave occurred by the coupling of the
invasion of the new virus variant Omicron BA.5 and multiple social factors,
followed by human interaction during the July holidays. Based on the results
that the domain growth rate $K$ and the infection rise time $t_\mathrm{on}$
were almost independent of the initial susceptible $D_\mathrm{s}$, dense
nucleation followed by a near growth model was deduced. The quantity $K \cdot
t_\mathrm{on}$ was considered to represent the infectivity of the virus.


  All-cause mortality is a very coarse grain, albeit very reliable, index to
check the health implications of lifestyle determinants, systemic threats and
socio-demographic factors. In this work we adopt a statistical-mechanics
approach to the analysis of temporal fluctuations of all-cause mortality,
focusing on the correlation structure of this index across different regions of
Italy. The correlation network among the 20 Italian regions was reconstructed
using temperature oscillations and travellers' flux (as a function of distance
and region's attractiveness, based on GDP), allowing for a separation between
infective and non-infective death causes. The proposed approach allows
monitoring of emerging systemic threats in terms of anomalies of correlation
network structure.


  The popularity of molecular computation has given rise to several models of
abstraction, one of the more recent ones being Chemical Reaction Networks
(CRNs). These are equivalent to other popular computational models, such as
Vector Addition Systems and Petri-Nets, and restricted versions are equivalent
to Population Protocols. This paper continues the work on core reachability
questions related to Chemical Reaction Networks; given two configurations, can
one reach the other according to the system's rules? With no restrictions,
reachability was recently shown to be Ackermann-complete, this resolving a
decades-old problem.
  Here, we fully characterize monotone reachability problems based on various
restrictions such as the rule size, the number of rules that may create a
species (k-source) or consume a species (k-consuming), the volume, and whether
the rules have an acyclic production order (feed-forward). We show
PSPACE-completeness of reachability with only bimolecular reactions with
two-source and two-consuming rules. This proves hardness of reachability in
Population Protocols, which was unknown. Further, this shows reachability in
CRNs is PSPACE-complete with size-2 rules, which was previously only known with
size-5 rules. This is achieved using techniques within the motion planning
framework.
  We provide many important results for feed-forward CRNs where rules are
single-source or single-consuming. We show that reachability is solvable in
polynomial time if the system does not contain special void or autogenesis
rules. We then fully characterize all systems of this type and show that if you
allow void/autogenesis rules, or have more than one source and one consuming,
the problems become NP-complete. Finally, we show several interesting special
cases of CRNs based on these restrictions or slight relaxations and note future
significant open questions related to this taxonomy.


  Machine learning has become a promising approach for molecular modeling.
Positional quantities, such as interatomic distances and bond angles, play a
crucial role in molecule physics. The existing works rely on careful manual
design of their representation. To model the complex nonlinearity in predicting
molecular properties in an more end-to-end approach, we propose to encode the
positional quantities with a learnable embedding that is continuous and
differentiable. A regularization technique is employed to encourage embedding
smoothness along the physical dimension. We experiment with a variety of
molecular property and force field prediction tasks. Improved performance is
observed for three different model architectures after plugging in the proposed
positional encoding method. In addition, the learned positional encoding allows
easier physics-based interpretation. We observe that tasks of similar physics
have the similar learned positional encoding.


  Mendel performed his experiments from 1856 and 1863, presented his results in
two meetings of the Natural Science Society in Brunn in February and March of
1865, and finally published them in the iconical paper of 1866, Versuche uber
Plflanzenhybriden. Two main translations to Spanish are available: the one done
by Prevosti in 1977, and a more recent one by Gomez-Graciani in 2014. Both were
translated from the English version of Druery and Bateson. Here, we present the
first direct translation from German to Spanish of this paper. This new version
corrects some errors detected in the previous versions, uses a more agile style
and includes the Darwinized terms pointed out by Abbot and Faribanks in a
recent paper (Genetics, 204(2):401-405, 2016; doi:
10.1534/genetics.116.194613).


  We consider the Markovian SIR epidemic model, but where individuals either
recover naturally or are diagnosed, the latter implying isolation and subject
to being contact traced. We focus on the novel digital contact tracing based on
individuals using a tracing app and consider both the situation when this is
the only form of contact tracing, and when combined with manual contact
tracing. It is proven that the epidemic process converges to a limiting process
as the population size $n$ grows large. The limiting process is not a branching
process (as is common for epidemic processes) because infectees may be contact
traced "from" their infectors, thus creating dependencies. But by instead
treating to-be-traced components as macro-individuals allows a multi-type
branching process interpretation, thus giving an expression for the
reproduction number $R$. This reproduction number (for the digital-only case)
is then converted to an individual reproduction number $R^{(ind)}$ being easier
to interpret. It is proven that, contrary to $R$, $R^{(ind)}$ decays
monotonically with $\pi$ (the fraction of app-users), and the two reproduction
numbers are proven to have the same threshold at 1. The digital (only) contact
tracing is then compared with earlier results for manual (only) contact tracing
where a fraction $p$ of all contacts are traced. It is proven that the critical
fraction app-users $\pi_c$ for which digital contact tracing is reduced to
$R=1$, and the critical fraction manually contact traced $p_c$ to reach $R=1$,
satisfy $\pi_c>p_c$.


  Complex system stability can be studied via linear stability analysis using
Random Matrix Theory (RMT) or via feasibility (requiring positive equilibrium
abundances). Both approaches highlight the importance of interaction structure.
Here we show, analytically and numerically, how RMT and feasibility approaches
can be complementary. In generalised Lotka-Volterra (GLV) models with random
interaction matrices, feasibility increases when predator-prey interactions
increase; increasing competition/mutualism has the opposite effect. These
changes have crucial impact on the stability of the GLV model.


  Post-transcriptional modifications are crucial for RNA function, with roles
ranging from the stabilization of functional RNA structures to modulation of
RNA--protein interactions. Additionally, artificially modified RNAs have been
suggested as optimal oligonucleotides for therapeutic purposes. The impact of
chemical modifications on secondary structure has been rationalized for some of
the most common modifications. However, the characterization of how the
modifications affect the three-dimensional RNA structure and dynamics and its
capability to bind proteins is still highly challenging. Molecular dynamics
simulations, coupled with enhanced sampling methods and integration of
experimental data, provide a direct access to RNA structural dynamics. In the
context of RNA chemical modifications, alchemical simulations where a wild type
nucleotide is converted to a modified one are particularly common. In this
Chapter, we review recent molecular dynamics studies of modified
ribonucleotides. We discuss the technical aspects of the reviewed works,
including the employed force fields, enhanced sampling methods, and alchemical
methods, in a way that is accessible to experimentalists. Finally, we provide
our perspective on this quickly growing field of research. The goal of this
Chapter is to provide a guide for experimentalists to understand molecular
dynamics works and, at the same time, give to molecular dynamics experts a
solid review of published articles that will be a useful starting point for new
research.


  Climate change has been shown to affect fitness-related traits in a wide
range of taxa; for instance, warming leads to phenological advancements in many
plant and animal species. The influence of climate change on social and
secondary sexual traits, that are associated with fitness due to their role as
quality signals, is however unknown. Here, we use more than 5800 observations
collected on two Mediterranean blue tit subspecies (Cyanistes caeruleus
caeruleus and C.c. ogliastrae) to explore whether blue crown and yellow breast
patch colourations have changed over the past 15 years. Our data suggests that
colouration has become duller and less chromatic in both sexes. In addition, in
the Corsican C.c. ogliastrae, but not in the mainland C.c. caeruleus, the
decrease is associated with an increase in temperature at moult. Quantitative
genetic analyses do not reveal any microevolutionary change in the colour
traits along the study period, strongly suggesting that the observed change
over time was caused by a plastic response to the environmental conditions.
Overall, this study suggests that ornamental colourations could become less
conspicuous due to warming, revealing climate change effects on sexual and
social ornaments and calling for further research on the proximate mechanisms
behind these effects.


  In this paper we used an adapted version of an existing simulation model of
SARS-CoV-2 transmission in Scotland to investigate the rise of the Omicron
variant of concern, in order to evaluate plausible scenarios for transmission
advantage and vaccine immune escape relative to the Delta variant. We also
explored possible outcomes of different levels of imposed non-pharmaceutical
intervention. The initial results of these scenarios were used to inform the
Scottish Government in the early outbreak stages of the Omicron variant.
  We use an explicitly spatial agent-based simulation model combined with
spatially fine-grained COVID-19 observation data from Public Health Scotland.
Using the model with parameters fit over the Delta variant epidemic, some
initial assumptions about Omicron transmission advantage and vaccine escape,
and a simple growth rate fitting procedure, we were able to capture the initial
outbreak dynamics for Omicron. We also find the modelled dynamics hold up to
retrospective scrutiny.
  We found that the modelled imposition of extra non-pharmaceutical
interventions planned by the Scottish Government at the time would likely have
little effect in light of the transmission advantage held by the Omicron
variant and the fact that the planned interventions would have occurred too
late in the outbreak's trajectory. Finally, we found that any assumptions made
about the projected distribution of vaccines in the model population had little
bearing on the outcome, in terms of outbreak size and timing, rather that the
detailed landscape of immunity prior to the outbreak was of far greater
importance.


  Costly punishment has been suggested as a key mechanism for stabilizing
cooperation in one-shot games. However, recent studies have revealed that the
effectiveness of costly punishment can be diminished by second-order free
riders (i.e., cooperators who never punish defectors) and antisocial punishers
(i.e., defectors who punish cooperators). In a two-stage prisoner's dilemma
game, players not only need to choose between cooperation and defection in the
first stage, but also need to decide whether to punish their opponent in the
second stage. Here, we extend the theory of punishment in one-shot games by
introducing simple bots, who consistently choose prosocial punishment and do
not change their actions over time. We find that this simple extension of the
game allows prosocial punishment to dominate in well-mixed and networked
populations, and that the minimum fraction of bots required for the dominance
of prosocial punishment monotonically increases with increasing dilemma
strength. Furthermore, if humans possess a learning bias toward a "copy the
majority" rule or if bots are present at higher degree nodes in scale-free
networks, the fully dominance of prosocial punishment is still possible at a
high dilemma strength. These results indicate that introducing bots can be a
significant factor in establishing prosocial punishment. We therefore, provide
a novel explanation for the evolution of prosocial punishment, and note that
the contrasting results that emerge from the introduction of different types of
bots also imply that the design of the bots matters.


  Although substantial efforts have been made using graph neural networks
(GNNs) for AI-driven drug discovery (AIDD), effective molecular representation
learning remains an open challenge, especially in the case of insufficient
labeled molecules. Recent studies suggest that big GNN models pre-trained by
self-supervised learning on unlabeled datasets enable better transfer
performance in downstream molecular property prediction tasks. However, the
approaches in these studies require multiple complex self-supervised tasks and
large-scale datasets, which are time-consuming, computationally expensive, and
difficult to pre-train end-to-end. Here, we design a simple yet effective
self-supervised strategy to simultaneously learn local and global information
about molecules, and further propose a novel bi-branch masked graph transformer
autoencoder (BatmanNet) to learn molecular representations. BatmanNet features
two tailored complementary and asymmetric graph autoencoders to reconstruct the
missing nodes and edges, respectively, from a masked molecular graph. With this
design, BatmanNet can effectively capture the underlying structure and semantic
information of molecules, thus improving the performance of molecular
representation. BatmanNet achieves state-of-the-art results for multiple drug
discovery tasks, including molecular properties prediction, drug-drug
interaction, and drug-target interaction, on 13 benchmark datasets,
demonstrating its great potential and superiority in molecular representation
learning.


  In recent years, artificial intelligence has played an important role on
accelerating the whole process of drug discovery. Various of molecular
representation schemes of different modals (e.g. textual sequence or graph) are
developed. By digitally encoding them, different chemical information can be
learned through corresponding network structures. Molecular graphs and
Simplified Molecular Input Line Entry System (SMILES) are popular means for
molecular representation learning in current. Previous works have done attempts
by combining both of them to solve the problem of specific information loss in
single-modal representation on various tasks. To further fusing such
multi-modal imformation, the correspondence between learned chemical feature
from different representation should be considered. To realize this, we propose
a novel framework of molecular joint representation learning via Multi-Modal
information of SMILES and molecular Graphs, called MMSG. We improve the
self-attention mechanism by introducing bond level graph representation as
attention bias in Transformer to reinforce feature correspondence between
multi-modal information. We further propose a Bidirectional Message
Communication Graph Neural Network (BMC GNN) to strengthen the information flow
aggregated from graphs for further combination. Numerous experiments on public
property prediction datasets have demonstrated the effectiveness of our model.


  The study of compounds with antioxidant capabilities is of great interest to
the scientific community, as it has implications in several areas, from
Agricultural Sciences to Biological Sciences, including Food Engineering,
Medicine and Pharmacy. In applications related to human health, it is known
that antioxidant activity can delay or inhibit oxidative damage to cells,
reducing damage caused by free radicals, helping in the treatment, or even
preventing or postponing the onset of various diseases. Among the compounds
that have antioxidant properties, there are several classes of Phenolic
Compounds, which include several compounds with different chemical structures.
In this work, based on the molecular branching of compounds and their
intramolecular charge distributions, and using Molecular Topology, we propose a
significant topological-mathematical model to evaluate the potential of
candidate compounds to have an antioxidant function.


  Machine learning approaches have become popular for molecular modeling tasks,
including molecular force fields and properties prediction. Traditional
supervised learning methods suffer from scarcity of labeled data for particular
tasks, motivating the use of large-scale dataset for other relevant tasks. We
propose to pretrain neural networks on a dataset of 86 millions of molecules
with atom charges and 3D geometries as inputs and molecular energies as labels.
Experiments show that, compared to training from scratch, fine-tuning the
pretrained model can significantly improve the performance for seven molecular
property prediction tasks and two force field tasks. We also demonstrate that
the learned representations from the pretrained model contain adequate
information about molecular structures, by showing that linear probing of the
representations can predict many molecular information including atom types,
interatomic distances, class of molecular scaffolds, and existence of molecular
fragments. Our results show that supervised pretraining is a promising research
direction in molecular modeling


  The Voter model is a well-studied stochastic process that models the invasion
of a novel trait $A$ (e.g., a new opinion, social meme, genetic mutation,
magnetic spin) in a network of individuals (agents, people, genes, particles)
carrying an existing resident trait $B$. Individuals change traits by
occasionally sampling the trait of a neighbor, while an invasion bias
$\delta\geq 0$ expresses the stochastic preference to adopt the novel trait $A$
over the resident trait $B$. The strength of an invasion is measured by the
probability that eventually the whole population adopts trait $A$, i.e., the
fixation probability. In more realistic settings, however, the invasion bias is
not ubiquitous, but rather manifested only in parts of the network. For
instance, when modeling the spread of a social trait, the invasion bias
represents localized incentives. In this paper, we generalize the standard
biased Voter model to the positional Voter model, in which the invasion bias is
effectuated only on an arbitrary subset of the network nodes, called biased
nodes. We study the ensuing optimization problem, which is, given a budget $k$,
to choose $k$ biased nodes so as to maximize the fixation probability of a
randomly occurring invasion. We show that the problem is NP-hard both for
finite $\delta$ and when $\delta \rightarrow \infty$ (strong bias), while the
objective function is not submodular in either setting, indicating strong
computational hardness. On the other hand, we show that, when
$\delta\rightarrow 0$ (weak bias), we can obtain a tight approximation in
$O(n^{2\omega})$ time, where $\omega$ is the matrix-multiplication exponent. We
complement our theoretical results with an experimental evaluation of some
proposed heuristics.


  Network Medicine has improved the mechanistic understanding of disease,
offering quantitative insights into disease mechanisms, comorbidities, and
novel diagnostic tools and therapeutic treatments. Yet, most network-based
approaches rely on a comprehensive map of protein-protein interactions,
ignoring interactions mediated by non-coding RNAs (ncRNAs). Here, we
systematically combine experimentally confirmed binding interactions mediated
by ncRNA with protein-protein interactions, constructing the first
comprehensive network of all physical interactions in the human cell. We find
that the inclusion of ncRNA, expands the number of genes in the interactome by
46% and the number of interactions by 107%, significantly enhancing our ability
to identify disease modules. Indeed, we find that 132 diseases, lacked a
statistically significant disease module in the protein-based interactome, but
have a statistically significant disease module after inclusion of
ncRNA-mediated interactions, making these diseases accessible to the tools of
network medicine. We show that the inclusion of ncRNAs helps unveil
disease-disease relationships that were not detectable before and expands our
ability to predict comorbidity patterns between diseases. Taken together, we
find that including non-coding interactions improves both the breath and the
predictive accuracy of network medicine.


  A central problem in computational biophysics is protein structure
prediction, i.e., finding the optimal folding of a given amino acid sequence.
This problem has been studied in a classical abstract model, the HP model,
where the protein is modeled as a sequence of H (hydrophobic) and P (polar)
amino acids on a lattice. The objective is to find conformations maximizing H-H
contacts. It is known that even in this reduced setting, the problem is
intractable (NP-hard). In this work, we apply deep reinforcement learning (DRL)
to the two-dimensional HP model. We can obtain the conformations of best known
energies for benchmark HP sequences with lengths from 20 to 50. Our DRL is
based on a deep Q-network (DQN). We find that a DQN based on long short-term
memory (LSTM) architecture greatly enhances the RL learning ability and
significantly improves the search process. DRL can sample the state space
efficiently, without the need of manual heuristics. Experimentally we show that
it can find multiple distinct best-known solutions per trial. This study
demonstrates the effectiveness of deep reinforcement learning in the HP model
for protein folding.


  Heterogeneity is a hallmark of complex diseases. Regression-based
heterogeneity analysis, which is directly concerned with outcome-feature
relationships, has led to a deeper understanding of disease biology. Such an
analysis identifies the underlying subgroup structure and estimates the
subgroup-specific regression coefficients. However, most of the existing
regression-based heterogeneity analyses can only address disjoint subgroups;
that is, each sample is assigned to only one subgroup. In reality, some samples
have multiple labels, for example, many genes have several biological
functions, and some cells of pure cell types transition into other types over
time, which suggest that their outcome-feature relationships (regression
coefficients) can be a mixture of relationships in more than one subgroups, and
as a result, the disjoint subgrouping results can be unsatisfactory. To this
end, we develop a novel approach to regression-based heterogeneity analysis,
which takes into account possible overlaps between subgroups and high data
dimensions. A subgroup membership vector is introduced for each sample, which
is combined with a loss function. Considering the lack of information arising
from small sample sizes, an $l_2$ norm penalty is developed for each membership
vector to encourage similarity in its elements. A sparse penalization is also
applied for regularized estimation and feature selection. Extensive simulations
demonstrate its superiority over direct competitors. The analysis of Cancer
Cell Line Encyclopedia data and lung cancer data from The Cancer Genome Atlas
shows that the proposed approach can identify an overlapping subgroup structure
with favorable performance in prediction and stability.


  We consider a flexible Bayesian evidence synthesis approach to model the
age-specific transmission dynamics of COVID-19 based on daily mortality counts.
The temporal evolution of transmission rates in populations containing multiple
types of individual is reconstructed via an appropriate dimension-reduction
formulation driven by independent diffusion processes. A suitably tailored
compartmental model is used to learn the latent counts of infection, accounting
for fluctuations in transmission influenced by public health interventions and
changes in human behaviour. The model is fitted to freely available COVID-19
data sources from the UK, Greece and Austria and validated using a large-scale
seroprevalence survey in England. In particular, we demonstrate how model
expansion can facilitate evidence reconciliation at a latent level. The code
implementing this work is made freely available via the Bernadette R package.


  The relationship between genotype and phenotype plays a crucial role in
determining the function and robustness of biological systems. Here the
evolution progresses through the change in genotype, whereas the selection is
based on the phenotype, and genotype-phenotype relation also evolves. Theory
for such phenotypic evolution remains poorly-developed, in contrast to
evolution under the fitness landscape determined by genotypes. Here we provide
statistical-physics formulation of this problem by introducing replicas for
genotype and phenotype. We apply it to an evolution model, in which phenotypes
are given by spin configurations; genotypes are interaction matrix for spins to
give the Hamiltonian, and the fitness depends only on the configuration of a
subset of spins called target. We describe the interplay between the genetic
variations and phenotypic variances by noise in this model by our new approach
that extends the replica theory for spin-glasses to include spin-replica for
phenotypes and coupling-replica for genotypes. Within this framework we obtain
a phase diagram of the evolved phenotypes against the noise and selection
pressure, where each phase is distinguished by the fitness and overlaps for
genotypes and phenotypes. Among the phases, robust fitted phase, relevant to
biological evolution, is achieved under the intermediate level of noise
(temperature), where robustness to noise and to genetic mutation are
correlated, as a result of replica symmetry. We also find a trade-off between
maintaining a high fitness level of phenotype and acquiring a robust pattern of
genes as well as the dependence of this trade-off on the ratio between the size
of the functional (target) part to that of the remaining non-functional
(non-target) one. The selection pressure needed to achieve high fitness
increases with the fraction of target spins.


  We consider the simple epidemiological SIS model for a general heterogeneous
population introduced by Lajmanovich and Yorke (1976) in finite dimension, and
its infinite dimensional generalization we introduced in previous works. In
this model the basic reproducing number $R_0$ is given by the spectral radius
of an integral operator. If $R_0>1$, then there exists a maximal endemic
equilibrium. In this very general heterogeneous SIS model, we prove that
vaccinating according to the profile of this maximal endemic equilibrium
ensures herd immunity. Moreover, this vaccination strategy is critical: the
resulting effective reproduction number is exactly equal to one. As an
application, we estimate that if $R_0 = 2$ in an age-structured community with
mixing rates fitted to social activity, applying this strategy would require
approximately 29% less vaccine doses than the strategy which consists in
vaccinating uniformly a proportion $1 - 1/R_0$ of the population.
  From a dynamical systems point of view, we prove that the non-maximality of
an equilibrium $g$ is equivalent to its linear instability in the original
dynamics, and to the linear instability of the disease-free state in the
modified dynamics where we vaccinate according to $g$.


  International initiatives such as METABRIC (Molecular Taxonomy of Breast
Cancer International Consortium) have collected several multigenomic and
clinical data sets to identify the undergoing molecular processes taking place
throughout the evolution of various cancers. Numerous Machine Learning and
statistical models have been designed and trained to analyze these types of
data independently, however, the integration of such differently shaped and
sourced information streams has not been extensively studied. To better
integrate these data sets and generate meaningful representations that can
ultimately be leveraged for cancer detection tasks could lead to giving
well-suited treatments to patients. Hence, we propose a novel learning pipeline
comprising three steps - the integration of cancer data modalities as graphs,
followed by the application of Graph Neural Networks in an unsupervised setting
to generate lower-dimensional embeddings from the combined data, and finally
feeding the new representations on a cancer sub-type classification model for
evaluation. The graph construction algorithms are described in-depth as
METABRIC does not store relationships between the patient modalities, with a
discussion of their influence over the quality of the generated embeddings. We
also present the models used to generate the lower-latent space
representations: Graph Neural Networks, Variational Graph Autoencoders and Deep
Graph Infomax. In parallel, the pipeline is tested on a synthetic dataset to
demonstrate that the characteristics of the underlying data, such as homophily
levels, greatly influence the performance of the pipeline, which ranges between
51\% to 98\% accuracy on artificial data, and 13\% and 80\% on METABRIC. This
project has the potential to improve cancer data understanding and encourages
the transition of regular data sets to graph-shaped data.


  The diverse metabolic pathways are fundamental to all living organisms, as
they harvest energy, synthesize biomass components, produce molecules to
interact with the microenvironment, and neutralize toxins. While discovery of
new metabolites and pathways continues, the prediction of pathways for new
metabolites can be challenging. It can take vast amounts of time to elucidate
pathways for new metabolites; thus, according to HMDB only 60% of metabolites
get assigned to pathways. Here, we present an approach to identify pathways
based on metabolite structure. We extracted 201 features from SMILES
annotations, and identified new metabolites from PubMed abstracts and HMDB.
After applying clustering algorithms to both groups of features, we quantified
correlations between metabolites, and found the clusters accurately linked 92%
of known metabolites to their respective pathways. Thus, this approach could be
valuable for predicting metabolic pathways for new metabolites.


  We discover a robust self-supervised strategy tailored towards molecular
representations for generative masked language models through a series of
tailored, in-depth ablations. Using this pre-training strategy, we train
BARTSmiles, a BART-like model with an order of magnitude more compute than
previous self-supervised molecular representations. In-depth evaluations show
that BARTSmiles consistently outperforms other self-supervised representations
across classification, regression, and generation tasks setting a new
state-of-the-art on 11 tasks. We then quantitatively show that when applied to
the molecular domain, the BART objective learns representations that implicitly
encode our downstream tasks of interest. For example, by selecting seven
neurons from a frozen BARTSmiles, we can obtain a model having performance
within two percentage points of the full fine-tuned model on task Clintox.
Lastly, we show that standard attribution interpretability methods, when
applied to BARTSmiles, highlight certain substructures that chemists use to
explain specific properties of molecules. The code and the pretrained model are
publicly available.


  We applied coherence analysis to molecular dynamics simulations of the plant
protein crambin, a thionin storage protein found in Abyssinian cabbage.
Coherence analysis was developed by engineers to identify linear interactions,
without statistical assumptions. Coherence is greater than 0.9 between the
displacement of oxygen and nitrogen atoms of H bonds in alpha helices for
frequencies between 0.391 GHz and 5.08 GHz (corresponding reciprocally to times
of 2.56 ns and 0.197 ns). These H bonds act much like a linear system.
Unrelated atoms have uncorrelated motions and much smaller coherence, say 0.02.
Groups of atoms (that form a layer of an alpha helix) were averaged and the
coherence function of two groups was evaluated. Layers of the alpha helix form
a linear system, suggesting that the harmonic analysis of classical molecular
dynamics can successfully describe the allosteric interactions of the layers of
an alpha helix.


  In the genomic analysis, it is significant while challenging to identify
markers associated with cancer outcomes or phenotypes. Based on the biological
mechanisms of cancers and the characteristics of datasets as well, this paper
proposes a novel integrative interaction approach under the semiparametric
model, in which the genetic factors and environmental factors are included as
the parametric and nonparametric components, respectively. The goal of this
approach is to identify the genetic factors and gene-gene interactions
associated with cancer outcomes, and meanwhile, estimate the nonlinear effects
of environmental factors. The proposed approach is based on the threshold
gradient directed regularization (TGDR) technique. Simulation studies indicate
that the proposed approach outperforms in the identification of main effects
and interactions, and has favorable estimation and prediction accuracy compared
with the alternative methods. The analysis of non-small-cell lung carcinomas
(NSCLC) datasets from The Cancer Genome Atlas (TCGA) are conducted, showing
that the proposed approach can identify markers with important implications and
have favorable performance in prediction accuracy, identification stability,
and computation cost.


  With advanced imaging, sequencing, and profiling technologies, multiple omics
data become increasingly available and hold promises for many healthcare
applications such as cancer diagnosis and treatment. Multimodal learning for
integrative multi-omics analysis can help researchers and practitioners gain
deep insights into human diseases and improve clinical decisions. However,
several challenges are hindering the development in this area, including the
availability of easily accessible open-source tools. This survey aims to
provide an up-to-date overview of the data challenges, fusion approaches,
datasets, and software tools from several new perspectives. We identify and
investigate various omics data challenges that can help us understand the field
better. We categorize fusion approaches comprehensively to cover existing
methods in this area. We collect existing open-source tools to facilitate their
broader utilization and development. We explore a broad range of omics data
modalities and a list of accessible datasets. Finally, we summarize future
directions that can potentially address existing gaps and answer the pressing
need to advance multimodal learning for multi-omics data analysis.


  Energy flow through consumer-resource interactions is largely determined by
body size. Allometric relationships govern the dynamics of populations by
impacting rates of reproduction, as well as alternative sources of mortality,
which have differential impacts on smaller to larger organisms. Here we derive
and investigate the timescales associated with four alternative sources of
mortality for terrestrial mammals: mortality from starvation, mortality
associated with aging, mortality from consumption by predators, and mortality
introduced by anthropogenic subsidized harvest. The incorporation of these
allometric relationships into a minimal consumer-resource model illuminates
central constraints that may contribute to the structure of mammalian
communities. Our framework reveals that while starvation largely impacts
smaller-bodied species, the allometry of senescence is expected to be more
difficult to observe. In contrast, external predation and subsidized harvest
have greater impacts on the populations of larger-bodied species. Moreover, the
inclusion of predation mortality reveals mass thresholds for mammalian
herbivores, where dynamic instabilities may limit the feasibility of megafaunal
populations. We show how these thresholds vary with alternative predator-prey
mass relationships, which are not well understood within terrestrial systems.
Finally, we use our framework to predict the harvest pressure required to
induce mass-specific extinctions, which closely align with previous estimates
of anthropogenic megafaunal exploitation in both paleontological and historical
contexts. Together our results underscore the tenuous nature of megafaunal
populations, and how different sources of mortality may contribute to their
ephemeral nature over evolutionary time.


  Motivation: The abundance of gene flow in the Tree of Life challenges the
notion that evolution can be represented with a fully bifurcating process, as
this process cannot capture important biological realities like hybridization,
introgression, or horizontal gene transfer. Coalescent-based network methods
are increasingly popular, yet not scalable for big data, because they need to
perform a heuristic search in the space of networks as well as numerical
optimization that can be NP-hard.
  Results: Here, we introduce a novel method to reconstruct phylogenetic
networks based on algebraic invariants. While there is a long tradition of
using algebraic invariants in phylogenetics, our work is the first to define
phylogenetic invariants on concordance factors (frequencies of 4-taxon splits
in the input gene trees) to identify level-1 phylogenetic networks under the
multispecies coalescent model. Our novel inference methodology is
optimization-free as it only requires the evaluation of polynomial equations,
and as such, it bypasses the traversal of network space, yielding a
computational speed at least 10 times faster than the fastest-to-date network
methods. We illustrate the accuracy and speed of our new method on a variety of
simulated scenarios as well as in the estimation of a phylogenetic network for
the genus Canis.
  Availability and Implementation: We implement our novel theory on an
open-source publicly available Julia package PhyloDiamond.jl available at
https://github.com/solislemuslab/PhyloDiamond.jl with broad applicability
within the evolutionary biology community.
  Contact: solislemus@wisc.edu


  Heterogeneity is a hallmark of many complex diseases. There are multiple ways
of defining heterogeneity, among which the heterogeneity in genetic
regulations, for example GEs (gene expressions) by CNVs (copy number
variations) and methylation, has been suggested but little investigated.
Heterogeneity in genetic regulations can be linked with disease severity,
progression, and other traits and is biologically important. However, the
analysis can be very challenging with the high dimensionality of both sides of
regulation as well as sparse and weak signals. In this article, we consider the
scenario where subjects form unknown subgroups, and each subgroup has unique
genetic regulation relationships. Further, such heterogeneity is "guided" by a
known biomarker. We develop an MSF (Multivariate Sparse Fusion) approach, which
innovatively applies the penalized fusion technique to simultaneously determine
the number and structure of subgroups and regulation relationships within each
subgroup. An effective computational algorithm is developed, and extensive
simulations are conducted. The analysis of heterogeneity in the GE-CNV
regulations in melanoma and GE-methylation regulations in stomach cancer using
the TCGA data leads to interesting findings.


  The analysis of cancer omics data is a "classic" problem, however, still
remains challenging. Advancing from early studies that are mostly focused on a
single type of cancer, some recent studies have analyzed data on multiple
"related" cancer types/subtypes, examined their commonality and difference, and
led to insightful findings. In this article, we consider the analysis of
multiple omics datasets, with each dataset on one type/subtype of "related"
cancers. A Community Fusion (CoFu) approach is developed, which conducts marker
selection and model building using a novel penalization technique,
informatively accommodates the network community structure of omics
measurements, and automatically identifies the commonality and difference of
cancer omics markers. Simulation demonstrates its superiority over direct
competitors. The analysis of TCGA lung cancer and melanoma data leads to
interesting findings


  In cancer research, high-throughput profiling has been extensively conducted.
In recent studies, the integrative analysis of data on multiple cancer patient
groups/subgroups has been conducted. Such analysis has the potential to reveal
the genomic commonality as well as difference across groups/subgroups. However,
in the existing literature, methods with a special attention to the genomic
commonality and difference are very limited. In this study, a novel estimation
and marker selection method based on the sparse boosting technique is developed
to address the commonality/difference problem. In terms of technical
innovation, a new penalty and computation of increments are introduced. The
proposed method can also effectively accommodate the grouping structure of
covariates. Simulation shows that it can outperform direct competitors under a
wide spectrum of settings. The analysis of two TCGA (The Cancer Genome Atlas)
datasets is conducted, showing that the proposed analysis can identify markers
with important biological implications and have satisfactory prediction and
stability.


  Next-generation sequencing technologies have enhanced the scope of
Internet-of-Things (IoT) to include genomics for personalized medicine through
the increased availability of an abundance of genome data collected from
heterogeneous sources at a reduced cost. Given the sheer magnitude of the
collected data and the significant challenges offered by the presence of highly
similar genomic structure across species, there is a need for robust, scalable
analysis platforms to extract actionable knowledge such as the presence of
potentially zoonotic pathogens. The emergence of zoonotic diseases from novel
pathogens, such as the influenza virus in 1918 and SARS-CoV-2 in 2019 that can
jump species barriers and lead to pandemic underscores the need for scalable
metagenome analysis. In this work, we propose MG2Vec, a deep learning-based
solution that uses the transformer network as its backbone, to learn robust
features from raw metagenome sequences for downstream biomedical tasks such as
targeted and generalized pathogen detection. Extensive experiments on four
increasingly challenging, yet realistic diagnostic settings, show that the
proposed approach can help detect pathogens from uncurated, real-world clinical
samples with minimal human supervision in the form of labels. Further, we
demonstrate that the learned representations can generalize to completely
unrelated pathogens across diseases and species for large-scale metagenome
analysis. We provide a comprehensive evaluation of a novel representation
learning framework for metagenome-based disease diagnostics with deep learning
and provide a way forward for extracting and using robust vector
representations from low-cost next generation sequencing to develop
generalizable diagnostic tools.


  Antimicrobial resistance is one of the biggest health problem, especially in
the current period of COVID-19 pandemic. Due to the unique membrane-destruction
bactericidal mechanism, antimicrobial peptide-mimetic copolymers are paid more
attention and it is urgent to find more potential candidates with
broad-spectrum antibacterial efficacy and low toxicity. Artificial intelligence
has shown significant performance on small molecule or biotech drugs, however,
the higher-dimension of polymer space and the limited experimental data
restrict the application of existing methods on copolymer design. Herein, we
develop a universal random copolymer inverse design system via multi-model
copolymer representation learning, knowledge distillation and reinforcement
learning. Our system realize a high-precision antimicrobial activity prediction
with few-shot data by extracting various chemical information from multi-modal
copolymer representations. By pre-training a scaffold-decorator generative
model via knowledge distillation, copolymer space are greatly contracted to the
near space of existing data for exploration. Thus, our reinforcement learning
algorithm can be adaptive for customized generation on specific scaffolds and
requirements on property or structures. We apply our system on collected
antimicrobial peptide-mimetic copolymers data, and we discover candidate
copolymers with desired properties.


  The Luria--Delbr\"uck mutation model is a cornerstone of evolution theory and
has been mathematically formulated in a number of ways. In this paper we
illustrate how this model of mutation rates can be derived by means of
classical statistical mechanics tools, in particular by modeling the phenomenon
resorting to methodologies borrowed from classical kinetic theory of rarefied
gases. The aim is to construct a linear kinetic model that can reproduce the
Luria--Delbr\"uck distribution starting from the elementary interactions that
qualitatively and quantitatively describe the variation of mutated cells. The
kinetic description is easily adaptable to different situations and makes it
possible to clearly identify the differences between the elementary variations
leading to the formulations of Luria--Delbr\"uck, Lea--Coulson, and Kendall,
respectively. The kinetic approach additionally emphasizes basic principles
which not only help to unify existing results but also allow for useful
extensions.


  Carica papaya (CP) leaf is long known for its traditional pharmacological
effects against dengue virus and malaria. Therefore, CP could also be a
potential solution for the treatment of other infectious diseases, such as
coronavirus. In this study, we evaluate the potential effect of the ethanolic
CP leaf extract in inhibiting the enzymatic activity of three protein targets
in SARS-CoV-2 life cycle, which include the 3-chymotrypsin-like protease
(3CLpro), papain-like protease (PLpro) and the human transmembrane protein
serine 2 (TMPRSS2). Results demonstrate that CP leaf extract inhibits 3CLpro
and PLpro significantly with IC50 of 0.02 microgram/mL and 0.06 microgram/mL,
respectively, but it is inactive towards TMPRSS2. Phenol,
2-methyl-5-(1,2,2-trimethylcyclopentyl)-(S)- (17a) and beta-mannofuranoside,
farnesyl- (21a) were identified in the extract using GC-MS. These two compounds
demonstrated a stronger binding affinity towards the main proteases than
TMPRSS2 during the docking simulation, which agrees with the in-vitro study.
Further pharmacophore mapping suggests that 17a has a fit score higher than 21a
to the SARS-CoV-2 3CLpro pharmacophore model concluding that CP leaf extract
has the potential to be developed as a herbal SARS-CoV-2 antiviral agent.


  Exploiting the mathematical curiosity of intransitive dice, we present a
simple theoretical model for co-evolution that captures scales ranging from the
genome of the individual to the system-wide emergence of species diversity. We
study a set of evolving agents that interact competitively in a closed system,
in which both the dynamics of mutations and competitive advantage emerge
directly from interpreting a genome as the sides of a die. The model
demonstrates sympatric speciation where new species evolve from existing ones
while in contact with the entire ecosystem. Allowing free mutations both in the
genomes and the mutation rates, we find, in contrast to hierarchical models of
fitness, the emergence of a metastable state of finite mutation rate and
diversity.


  In 1995 I wrote a paper: "A Speed Limit for Evolution" whose main result was
that evolution must proceed rather slowly, in accordance with the earlier views
and intuitions of many authors. The abstract of the paper said: "The genetic
information expressed in some part of the phenotype of a species cannot
increase faster than a given rate, determined by the selection pressure on that
part. This rate is typically a small fraction of a bit per generation". This
result was derived in the presence of sexual reproduction and other effects
such as temporarily isolated sub-populations. In 1999 David Mackay published a
paper which apparently contradicted this result. In the abstract, he wrote "We
find striking differences between populations that have recombination and
populations that do not. If variation is produced by mutation alone, then the
entire population gains up to roughly 1 bit per generation. If variation is
created by recombination, the population can gain of the order of sqrt(G) bits
per generation." Mackay proposed that there were outstanding evolutionary
benefits to sexual reproduction, and that my result was too low by a very large
factor. He later repeated this result in a textbook he wrote in 2003. The
purpose of this note is to show that the key assumption of Mackays model, that
"fitness is a strictly additive trait" is so unrealistic as to render his
results irrelevant to any actual life form. In consequence, the speed limit I
derived is still valid, and has important consequences for human cognitive
evolution.


  Myriad viruses use positive-strand RNA molecules as their genomes. Far from
being only a repository of genetic material, viral RNA performs numerous other
functions mediated by its physical structure and chemical properties. In this
chapter, we focus on its structure and discuss how long RNA molecules can be
treated as branched polymers through planar graphs. We describe the major
results that can be obtained by this approach, in particular the observation
that viral RNA genomes have a characteristic compactness that sets them aside
from similar random RNAs. We also discuss how different parameters used in the
current RNA folding software influence the resulting structures and how they
can be related to experimentally observable quantities. Finally, we show how
the connection to branched polymers can be extended to take advantage of known
results from polymer physics and can be further moulded to include additional
interactions, such as excluded volume or electrostatics.


  In models for the evolution of predation from initially purely competitive
species interactions, the propensity of predation is most often assumed to be a
direct consequence of the relative morphological and physiological traits of
interacting species. Here we explore a model in which predation ability is an
independently evolving phenotypic feature, so that even when the relative
morphological or physiological traits allow for predation, predation only
occurs if the predation ability of individuals has independently evolved to
high enough values. In addition to delineating the conditions for the
evolutionary emergence of predation, the model reproduces stationary and
non-stationary multilevel food webs with the top predators not necessarily
having size superiority.


  Neglected tropical diseases (NTDs) continue to affect the livelihood of
individuals in countries in the Southeast Asia and Western Pacific region.
These diseases have been long existing and have caused devastating health
problems and economic decline to people in low- and middle-income (developing)
countries. An estimated 1.7 billion of the world's population suffer one or
more NTDs annually, this puts approximately one in five individuals at risk for
NTDs. In addition to health and social impact, NTDs inflict significant
financial burden to patients, close relatives, and are responsible for billions
of dollars lost in revenue from reduced labor productivity in developing
countries alone. There is an urgent need to better improve the control and
eradication or elimination efforts towards NTDs. This can be achieved by
utilizing machine learning tools to better the surveillance, prediction and
detection program, and combat NTDs through the discovery of new therapeutics
against these pathogens. This review surveys the current application of machine
learning tools for NTDs and the challenges to elevate the state-of-the-art of
NTDs surveillance, management, and treatment.


  Nuclear magnetic resonance (NMR) spectroscopy has become a formidable tool
for biochemistry and medicine. Although J-coupling carries essential structural
information it may also limit the spectral resolution. Homonuclear decoupling
remains a challenging problem. In this work, we introduce a new approach that
uses a specific coupling value as prior knowledge, and Hankel property of
exponential NMR signal to achieve the broadband heteronuclear decoupling using
the low-rank method. Our results on synthetic and realistic HMQC spectra
demonstrate that the proposed method not only effectively enhances resolution
by decoupling, but also maintains sensitivity and suppresses spectral
artefacts. The approach can be combined with the non-uniform sampling, which
means that the resolution can be further improved without any extra acquisition
time


  Long-term behaviors of biochemical systems are described by steady states in
deterministic models and stationary distributions in stochastic models.
Obtaining their analytic solutions can be done for limited cases, such as
linear or finite-state systems, as it generally requires solving many coupled
equations. Interestingly, analytic solutions can be easily obtained when
underlying networks have special topologies, called weak reversibility (WR) and
zero deficiency (ZD), and the kinetic law follows a generalized form of
mass-action kinetics. However, such desired topological conditions do not hold
for the majority of cases. Thus, translating networks to have WR and ZD while
preserving the original dynamics was proposed. Yet, this approach is limited
because manually obtaining the desired network translation among the large
number of candidates is challenging. Here, we prove necessary conditions for
having WR and ZD after translation, and based on these conditions, we develop a
user-friendly computational package, TOWARDZ, that automatically and
efficiently identifies translated networks with WR and ZD. This allows us to
quantitatively examine how likely it is to obtain WR and ZD after translation
depending on the number of species and reactions. Importantly, we also describe
how our package can be used to analytically derive steady states of
deterministic models and stationary distributions of stochastic models. TOWARDZ
provides an effective tool to analyze biochemical systems.