  Cryogenic electron microscopy (cryo-EM) has transformed structural biology by
allowing to reconstruct 3D biomolecular structures up to near-atomic
resolution. However, the 3D reconstruction process remains challenging, as the
3D structures may exhibit substantial shape variations, while the 2D image
acquisition suffers from a low signal-to-noise ratio, requiring to acquire very
large datasets that are time-consuming to process. Current reconstruction
methods are precise but computationally expensive, or faster but lack a
physically-plausible model of large molecular shape variations. To fill this
gap, we propose CryoChains that encodes large deformations of biomolecules via
rigid body transformation of their chains, while representing their finer shape
variations with the normal mode analysis framework of biophysics. Our synthetic
data experiments on the human GABA\textsubscript{B} and heat shock protein show
that CryoChains gives a biophysically-grounded quantification of the
heterogeneous conformations of biomolecules, while reconstructing their 3D
molecular structures at an improved resolution compared to the current fastest,
interpretable deep learning method.


  Recent advances in microscopy and 3D reconstruction methods have allowed for
characterization of cellular morphology in unprecedented detail, including the
irregular geometries of intracellular subcompartments such as membrane-bound
organelles. These geometries are now compatible with predictive modeling of
cellular function. Biological cells respond to stimuli through sequences of
chemical reactions generally referred to as cell signaling pathways. The
propagation and reaction of chemical substances in cell signaling pathways can
be represented by coupled nonlinear systems of reaction-transport equations.
These reaction pathways include numerous chemical species that react across
boundaries or interfaces (e.g., the cell membrane and membranes of organelles
within the cell) and domains (e.g., the bulk cell volume and the interior of
organelles). Such systems of multi-dimensional partial differential equations
(PDEs) are notoriously difficult to solve because of their high dimensionality,
non-linearities, strong coupling, stiffness, and potential instabilities. In
this work, we describe Spatial Modeling Algorithms for Reactions and Transport
(SMART), a high-performance finite-element-based simulation package for model
specification and numerical simulation of spatially-varying reaction-transport
processes. SMART is based on the FEniCS finite element library, provides a
symbolic representation framework for specifying reaction pathways, and
supports geometries in 2D and 3D including large and irregular cell geometries
obtained from modern ultrastructural characterization methods.


  From 08-12 August, 2022, 32 individuals participated in a workshop, Stability
and Fluctuations in Complex Ecological Systems, at the Lorentz Center, located
in Leiden, The Netherlands. An interdisciplinary dialogue between ecologists,
mathematicians, and physicists provided a foundation of important problems to
consider over the next 5-10 years. This paper outlines eight areas including
(1) improving our understanding of the effect of scale, both temporal and
spatial, for both deterministic and stochastic problems; (2) clarifying the
different terminologies and definitions used in different scientific fields;
(3) developing a comprehensive set of data analysis techniques arising from
different fields but which can be used together to improve our understanding of
existing data sets; (4) having theoreticians/computational scientists
collaborate closely with empirical ecologists to determine what new data should
be collected; (5) improving our knowledge of how to protect and/or restore
ecosystems; (6) incorporating socio-economic effects into models of ecosystems;
(7) improving our understanding of the role of deterministic and stochastic
fluctuations; (8) studying the current state of biodiversity at the functional
level, taxa level and genome level.


  Opioid Use Disorder (OUD) has emerged as a significant global public health
issue, with complex multifaceted conditions. Due to the lack of effective
treatment options for various conditions, there is a pressing need for the
discovery of new medications. In this study, we propose a deep generative model
that combines a stochastic differential equation (SDE)-based diffusion modeling
with the latent space of a pretrained autoencoder model. The molecular
generator enables efficient generation of molecules that are effective on
multiple targets, specifically the mu, kappa, and delta opioid receptors.
Furthermore, we assess the ADMET (absorption, distribution, metabolism,
excretion, and toxicity) properties of the generated molecules to identify
drug-like compounds. To enhance the pharmacokinetic properties of some lead
compounds, we employ a molecular optimization approach. We obtain a diverse set
of drug-like molecules. We construct binding affinity predictors by integrating
molecular fingerprints derived from autoencoder embeddings, transformer
embeddings, and topological Laplacians with advanced machine learning
algorithms. Further experimental studies are needed to evaluate the
pharmacological effects of these drug-like compounds for OUD treatment. Our
machine learning platform serves as a valuable tool in designing and optimizing
effective molecules for addressing OUD.


  Models of diffusive processes that occur on evolving domains are frequently
employed to describe biological and physical phenomena, such as diffusion
within expanding tissues or substrates. Previous investigations into these
models either report numerical solutions or require an assumption of linear
diffusion to determine exact solutions. Unfortunately, numerical solutions do
not reveal the relationship between the model parameters and the solution
features. Additionally, experimental observations typically report the presence
of sharp fronts, which are not captured by linear diffusion. Here we address
both limitations by presenting exact sharp-fronted solutions to a model of
degenerate nonlinear diffusion on a growing domain. We obtain the solution by
identifying a series of transformations that converts the model of a nonlinear
diffusive process on an evolving domain to a nonlinear diffusion equation on a
fixed domain, which admits known exact solutions for certain choices of
diffusivity functions. We determine expressions for critical time scales and
domain growth rates such that the diffusive population never reaches the domain
boundaries and hence the solution remains valid.


  In this work we propose a stochastic differential equation (SDE) for
modelling health related quality of life (HRQoL) over a lifespan. HRQoL is
assumed to be bounded between 0 and 1, equivalent to death and perfect health,
respectively. Drift and diffusion parameters of the SDE are chosen to mimic
decreasing HRQoL over life and ensuring epidemiological meaningfulness. The
Euler-Maruyama method is used to simulate trajectories of individuals in a
population of n = 1000 people. Age of death of an individual is simulated as a
stopping time with Weibull distribution conditioning the current value of HRQoL
as time-varying covariate. The life expectancy and health adjusted life years
are compared to the corresponding values for German women.


  Recent evidence suggests that nongenetic (epigenetic) mechanisms play an
important role at all stages of cancer evolution. In many cancers, these
mechanisms have been observed to induce dynamic switching between two or more
cell states, which commonly show differential responses to drug treatments. To
understand how these cancers evolve over time, and how they respond to
treatment, we need to understand the state-dependent rates of cell
proliferation and phenotypic switching. In this work, we propose a rigorous
statistical framework for estimating these parameters, using data from commonly
performed cell line experiments, where phenotypes are sorted and expanded in
culture. The framework explicitly models the stochastic dynamics of cell
division, cell death and phenotypic switching, and it provides likelihood-based
confidence intervals for the model parameters. The input data can be either the
fraction of cells or the number of cells in each state at one or more time
points. Through a combination of theoretical analysis and numerical
simulations, we show that when cell fraction data is used, the rates of
switching may be the only parameters that can be estimated accurately. On the
other hand, using cell number data enables accurate estimation of the net
division rate for each phenotype, and it can even enable estimation of the
state-dependent rates of cell division and cell death. We conclude by applying
our framework to a publicly available dataset.


  We consider a single genetic locus with two alleles $A_1$ and $A_2$ in a
large haploid population. The locus is subject to selection and two-way, or
recurrent, mutation. Assuming the allele frequencies follow a Wright-Fisher
diffusion and have reached stationarity, we describe the asymptotic behaviors
of the conditional gene genealogy and the latent mutations of a sample with
known allele counts, when the count $n_1$ of allele $A_1$ is fixed, and when
either or both the sample size $n$ and the selection strength
$\lvert\beta\rvert$ tend to infinity. Our study extends previous work under
neutrality to the case of non-neutral rare alleles, asserting that when
selection is not too strong relative to the sample size, even if it is strongly
positive or strongly negative in the usual sense ($\beta\to -\infty$ or
$\beta\to +\infty$), the number of latent mutations of the $n_1$ copies of
allele $A_1$ follows the same distribution as the number of alleles in the
Ewens sampling formula. On the other hand, very strong positive selection
relative to the sample size leads to neutral gene genealogies with a single
ancient latent mutation. We also demonstrate robustness of our asymptotic
results against changing population sizes, when one of $\lvert\beta\rvert$ or
$n$ is large.


  Proteolysis-Targeting Chimeras (PROTACs) represent a novel class of small
molecules which are designed to act as a bridge between an E3 ligase and a
disease-relevant protein, thereby promoting its subsequent degradation. PROTACs
are composed of two protein binding "active" domains, linked by a "linker"
domain. The design of the linker domain is challenging due to geometric and
chemical constraints given by its interactions, and the need to maximize
drug-likeness. To tackle these challenges, we introduce ShapeLinker, a method
for de novo design of linkers. It performs fragment-linking using reinforcement
learning on an autoregressive SMILES generator. The method optimizes for a
composite score combining relevant physicochemical properties and a novel,
attention-based point cloud alignment score. This new method successfully
generates linkers that satisfy both relevant 2D and 3D requirements, and
achieves state-of-the-art results in producing novel linkers assuming a target
linker conformation. This allows for more rational and efficient PROTAC design
and optimization. Code and data are available at
https://github.com/aivant/ShapeLinker.


  Logical modeling is a powerful tool in biology, offering a system-level
understanding of the complex interactions that govern biological processes. A
gap that hinders the scalability of logical models is the need to specify the
update function of every vertex in the network depending on the status of its
predecessors. To address this, we introduce in this paper the concept of strong
regulation, where a vertex is only updated to active/inactive if all its
predecessors agree in their influences; otherwise, it is set to ambiguous. We
explore the interplay between active, inactive, and ambiguous influences in a
network. We discuss the existence of phenotype attractors in such networks,
where the status of some of the variables is fixed to active/inactive, while
the others can have an arbitrary status, including ambiguous.


  Modern high-throughput sequencing assays efficiently capture not only gene
expression and different levels of gene regulation but also a multitude of
genome variants. Focused analysis of alternative alleles of variable sites at
homologous chromosomes of the human genome reveals allele-specific gene
expression and allele-specific gene regulation by assessing allelic imbalance
of read counts at individual sites. Here we formally describe an advanced
statistical framework for detecting the allelic imbalance in allelic read
counts at single-nucleotide variants detected in diverse omics studies
(ChIP-Seq, ATAC-Seq, DNase-Seq, CAGE-Seq, and others). MIXALIME accounts for
copy-number variants and aneuploidy, reference read mapping bias, and provides
several scoring models to balance between sensitivity and specificity when
scoring data with varying levels of experimental noise-caused overdispersion.


  Loop extrusion by motor proteins mediates the attractive interactions in
chromatin on the length scale of megabases, providing the polymer with a
well-defined structure and at the same time determining its dynamics. The mean
square displacement of chromatin loci varies from a Rouse-like scaling to a
more constrained subdiffusion, depending on cell type, genomic region and time
scale. With a simple polymeric model, we show that such a Rouse-like dynamics
occurs when the parameters of the model are chosen so that contacts are local
along the chain, while in presence of non-local contacts, we observe
subdiffusion at short time scales with exponents smaller than 0.5. Such
exponents are independent of the detailed choice of the parameters and build a
master curve that depends only on the mean locality of the resulting contacts.
We compare the loop-extrusion model with a polymeric model with static links,
showing that also in this case only the presence of nonlocal contacts can
produce low-exponent subdiffusion. We interpret these results in terms of a
simple analytical model.


  Shaped by natural selection and other evolutionary forces, an organism's
evolutionary history is reflected through its genome sequence, content of
functional elements and organization. Consequently, organisms connected through
phylogeny, metabolic or morphological traits, geographical proximity, or
habitat features are likely to exhibit similarities in their genomes. These
similarities give rise to expectations about the content of genomes within
these organism groups.
  Such expectations are often informally expressed in scientific literature,
focusing on the analysis of individual genomes or comparisons among related
groups of organisms. Our objective is to develop a system for formalized
expectations as rules, facilitating automated verification, and evaluation of
newly sequenced genomes.
  In this study, we present a database comprising rules manually extracted from
scientific literature. Furthermore, we explore the feasibility of automatizing
the extraction and analysis process using large language models, such as GPT3.5
and GPT4.
  We have developed a web application, EGCWebApp, which enables users to
visualize and edit the rules. Additionally, we provided a Python library and
command-line tools collection, egctools, to further extend the functionality
for processing and managing these rules.


  This paper introduces an innovative model for infectious diseases in
predator-prey populations. We not only prove the existence of global
non-negative solutions but also establish essential criteria for the system's
decline and sustainability. Furthermore, we demonstrate the presence of a Borel
invariant measure, adding a new dimension to our understanding of the system.
To illustrate the practical implications of our findings, we present numerical
results. With our model's comprehensive approach, we aim to provide valuable
insights into the dynamics of infectious diseases and their impact on
predator-prey populations.


  Motivation:Classification of gene trees is an important task both in analysis
of multi-locus phylogenetic data, and assessment of the convergence of Markov
Chain Monte Carlo (MCMC) analyses used in Bayesian phylogenetic tree
reconstruction. The logistic regression model is one of the most popular
classification models in statistical learning, thanks to its computational
speed and interpretability. However, it is not appropriate to directly apply
the logistic regression model to a set of phylogenetic trees with the same set
of leaf labels, as the space of phylogenetic trees is not Euclidean. Results:
It is well-known in tropical geometry and phylogenetics that the space of
phylogenetic trees is a tropical linear space in terms of the max-plus algebra.
Therefore, in this paper, we propose an analogue approach of the logistic
regression model in the setting of tropical geometry. In our proposed method,
we consider two cases: where the numbers of the species trees are fixed as one
and two, and we estimate the species tree(s) from a sample of gene trees
distributed over the space of ultrametrics, which is a tropical linear space.
we show that both models are statistically consistent and bounds on the
generalization error of both models are derived. Finally, we conduct
computational experiments on simulated data generated by the multi-species
coalescent model and apply our model to African coelacanth genomes to infer the
species tree.


  The identification of compound-protein interactions (CPI) plays a critical
role in drug screening, drug repurposing, and combination therapy studies. The
effectiveness of CPI prediction relies heavily on the features extracted from
both compounds and target proteins. While various prediction methods employ
different feature combinations, both molecular-based and network-based models
encounter the common obstacle of incomplete feature representations. Thus, a
promising solution to this issue is to fully integrate all relevant CPI
features. This study proposed a novel model named MCPI, which is designed to
improve the prediction performance of CPI by integrating multiple sources of
information, including the PPI network, CCI network, and structural features of
CPI. The results of the study indicate that the MCPI model outperformed other
existing methods for predicting CPI on public datasets. Furthermore, the study
has practical implications for drug development, as the model was applied to
search for potential inhibitors among FDA-approved drugs in response to the
SARS-CoV-2 pandemic. The prediction results were then validated through the
literature, suggesting that the MCPI model could be a useful tool for
identifying potential drug candidates. Overall, this study has the potential to
advance our understanding of CPI and guide drug development efforts.


  We study the SIR model of epidemics on positively correlated heterogeneous
networks with population variability, and explore the dependence of the final
outbreak size on the network heterogeneity strength and basic reproduction
number $R_0$ -- the ratio between the infection and recovery rates per
individual. We reveal a critical value $R_0^c$, above which the maximal
outbreak size is obtained at zero heterogeneity, but below which, the maximum
is obtained at finite heterogeneity strength. This second-order phase
transition, universal for all network distributions with finite standardized
moments indicates that, network heterogeneity can greatly increase the final
outbreak size. We also show that this effect can be enhanced by adding
population heterogeneity, in the form of varying inter-individual
susceptibility and infectiousness. Notably, our results provide key insight as
to the predictability of the well-mixed SIR model for the final outbreak size,
in realistic scenarios.


  Modelling gene-gene epistatic interactions when computing genetic risk scores
is not a well-explored subfield of genetics and could have potential to improve
risk stratification in practice. Though applications of machine learning (ML)
show promise as an avenue of improvement for current genetic risk assesments,
they frequently suffer from the problem of two many features and to little
data. We propose a method that when combined with ML allows information from
individual genetic contributors to be preserved while incorporating information
on their interactions in a single feature. This allows second-order analysis,
while simultaneously increasing the number of input features to ML models as
little as possible. We presented three methods that can be utilized to account
for genetic interactions. We found that interaction methods that preserved
information from the constituent SNPs performed significantly better than the
simplest interaction method. Since the currently available ML methods are able
to account for complex interactions, utilizing raw SNP genotypes alone is
sufficient because the simplest model outperforms all the interaction methods
Given that understanding and accounting for epistatic interactions is one of
the most promising avenues for increasing explained variability in heritable
disease, this work represents a first step toward an algorithmic interaction
method that preserves the information in each component. This is relevant not
only because of potential improvements in model quality, but also because
explicit interaction terms allow a human readable interpretation of potential
interaction pathways within the disease.


  Compartmentalised biochemical reactions are a ubiquitous building block of
biological systems. The interplay between chemical and compartmental dynamics
can drive rich and complex dynamical behaviors that are difficult to analyse
mathematically -- especially in the presence of stochasticity. We have recently
proposed an effective moment equation approach to study the statistical
properties of compartmentalised biochemical systems. So far, however, this
approach is limited to polynomial rate laws and moreover, it relies on suitable
moment closure approximations, which can be difficult to find in practice. In
this work we propose a systematic method to derive closed moment dynamics for
compartmentalised biochemical systems. We show that for the considered class of
systems, the moment equations involve expectations over functions that
factorize into two parts, one depending on the molecular content of the
compartments and one depending on the compartment number distribution. Our
method exploits this structure and approximates each function with suitable
polynomial expansions, leading to a closed system of moment equations. We
demonstrate the method using three systems inspired by cell populations and
organelle networks and study its accuracy across different dynamical regimes.


  Molecular representation learning (MRL) is a fundamental task for drug
discovery. However, previous deep-learning (DL) methods focus excessively on
learning robust inner-molecular representations by mask-dominated pretraining
framework, neglecting abundant chemical reactivity molecular relationships that
have been demonstrated as the determining factor for various molecular property
prediction tasks. Here, we present MolCAP to promote MRL, a graph pretraining
Transformer based on chemical reactivity (IMR) knowledge with prompted
finetuning. Results show that MolCAP outperforms comparative methods based on
traditional molecular pretraining framework, in 13 publicly available molecular
datasets across a diversity of biomedical tasks. Prompted by MolCAP, even basic
graph neural networks are capable of achieving surprising performance that
outperforms previous models, indicating the promising prospect of applying
reactivity information for MRL. In addition, manual designed molecular templets
are potential to uncover the dataset bias. All in all, we expect our MolCAP to
gain more chemical meaningful insights for the entire process of drug
discovery.


  Traditionally, population models distinguish individuals on the basis of
their current state. Given a distribution, a discrete time model then specifies
(precisely in deterministic models, probabilistically in stochastic models) the
population distribution at the next time point. The renewal equation
alternative concentrates on newborn individuals and the model specifies the
production of offspring as a function of age. This has two advantages: (i) as a
rule, there are far fewer birth states than individual states in general, so
the dimension is often low; (ii) it relates seamlessly to the next-generation
matrix and the basic reproduction number. Here we start from the renewal
equation for the births and use results of Feller and Thieme to characterise
the asymptotic large time behaviour. Next we explicitly elaborate the
relationship between the two bookkeeping schemes. This allows us to transfer
the characterisation of the large time behaviour to traditional
structured-population models.


  Viruses present an amazing genetic variability. An ensemble of infecting
viruses, also called a viral quasispecies, is a cloud of mutants centered
around a specific genotype. The simplest model of evolution, whose equilibrium
state is described by the quasispecies equation, is the Moran--Kingman model.
For the sharp peak landscape, we perform several exact computations and we
derive several exact formulas. We obtain also an exact formula for the
quasispecies distribution, involving a series and the mean fitness. A very
simple formula for the mean Hamming distance is derived, which is exact and
which do not require a specific asymptotic expansion (like sending the length
of the macromolecules to $\infty$ or the mutation probability to $0$). We try
also to extend these formulas to a general fitness landscape. We obtain an
equation involving the covariance of the fitness and the Hamming class number
in the quasispecies distribution. With the help of these formulas, we discuss
the phenomenon of the error threshold and the notion of quasispecies. We
recover the limiting quasipecies distribution in the long chain regime. We go
beyond the sharp peak landscape and we consider fitness landscapes having
finitely many peaks and a plateau--type landscape. We finally prove rigorously
within this framework the possible occurrence of the survival of the flattest,
a phenomenon which has been previously discovered by Wilke, Wang, Ofria, Lenski
and Adami and which has been investigated in several works.


  Single-cell omics enable the profiles of cells, which contain large numbers
of biological features, to be quantified. Cluster analysis, a dimensionality
reduction process, is used to reduce the dimensions of the data to make it
computationally tractable. In these analyses, cells are represented as vectors
in $n$-Dimensional space, where each dimension corresponds to a certain cell
feature. The distance between cells is used as a surrogate measure of
similarity, providing insight into the cell's state, function, and genetic
mechanisms. However, as cell profiles are clustered in 3D or higher-dimensional
space, it remains unknown which distance metric provides the most accurate
spatiotemporal representation of similarity, limiting the interpretability of
the data. I propose and prove a generalized proposition and set of corollaries
that serve as a criterion to determine which of the standard distance measures
is most accurate for conveying cell profile heterogeneity. Each distance method
is evaluated via statistical, geometric, and topological proofs, which are
formalized into a set of criteria. In this paper, I present the putative,
first-ever method to elect the most accurate and precise distance metrics with
any profiling modality, which are determined to be the Wasserstein distance and
cosine similarity metrics, respectively, in general cases. I also identify
special cases in which the criterion may select non-standard metrics. Combining
the metric properties selected by the criterion, I develop a novel, custom,
optimal distance metric that demonstrates superior computational efficiency,
peak annotation, motif identification, and footprinting for transcription
factor binding sites when compared with leading methods.


  High-content cellular imaging, transcriptomics, and proteomics data provide
rich and complementary views on the molecular layers of biology that influence
cellular states and function. However, the biological determinants through
which changes in multi-omics measurements influence cellular morphology have
not yet been systematically explored, and the degree to which cell imaging
could potentially enable the prediction of multi-omics directly from cell
imaging data is therefore currently unclear. Here, we address the question of
whether it is possible to predict bulk multi-omics measurements directly from
cell images using Image2Omics -- a deep learning approach that predicts
multi-omics in a cell population directly from high-content images stained with
multiplexed fluorescent dyes. We perform an experimental evaluation in
gene-edited macrophages derived from human induced pluripotent stem cell
(hiPSC) under multiple stimulation conditions and demonstrate that Image2Omics
achieves significantly better performance in predicting transcriptomics and
proteomics measurements directly from cell images than predictors based on the
mean observed training set abundance. We observed significant predictability of
abundances for 5903 (22.43%; 95% CI: 8.77%, 38.88%) and 5819 (22.11%; 95% CI:
10.40%, 38.08%) transcripts out of 26137 in M1 and M2-stimulated macrophages
respectively and for 1933 (38.77%; 95% CI: 36.94%, 39.85%) and 2055 (41.22%;
95% CI: 39.31%, 42.42%) proteins out of 4986 in M1 and M2-stimulated
macrophages respectively. Our results show that some transcript and protein
abundances are predictable from cell imaging and that cell imaging may
potentially, in some settings and depending on the mechanisms of interest and
desired performance threshold, even be a scalable and resource-efficient
substitute for multi-omics measurements.


  External flows of energy, entropy, and matter can cause sudden transitions in
the stability of biological and industrial systems, fundamentally altering
their dynamical function. How might we control and design these transitions in
chemical reaction networks? Here, we analyze transitions giving rise to complex
behavior in random reaction networks subject to external driving forces. In the
absence of driving, we characterize the uniqueness of the steady state and
identify the percolation of a giant connected component in these networks as
the number of reactions increases. When subject to chemical driving (influx and
outflux of chemical species), the steady state can undergo bifurcations,
leading to multistability or oscillatory dynamics. By quantifying the
prevalence of these bifurcations, we show how chemical driving and network
sparsity tend to promote the emergence of these complex dynamics and increased
rates of entropy production. We show that catalysis also plays an important
role in the emergence of complexity, strongly correlating with the prevalence
of bifurcations. Our results suggest that coupling a minimal number of chemical
signatures with external driving can lead to features present in biochemical
processes and abiogenesis.


  The discovery of general principles underlying the complexity and diversity
of cellular and developmental systems is a central and long-standing aim of
biology. Whilst new technologies collect data at an ever-accelerating rate,
there is growing concern that conceptual progress is not keeping pace. We
contend that this is due to a paucity of appropriate conceptual frameworks to
serve as a basis for general theories of mesoscale biological phenomena. In
exploring this issue, we have developed a foundation for one such framework,
termed the Core and Periphery (C&P) hypothesis, which reveals hidden generality
across the diverse and complex behaviors exhibited by cells and tissues. Here,
we present the C&P concept, provide examples of its applicability across
multiple scales, argue its consistency with evolution, and discuss key
implications and open questions. We propose that the C&P hypothesis could
unlock new avenues of conceptual progress in cell and developmental biology.


  A number of theoretical models have been developed in recent years modelling
epidemic spread in educational settings such as universities to help inform
re-opening strategies during the Covid-19 pandemic. However, these studies have
had differing conclusions as to the most effective non-pharmaceutical
interventions. They also largely assumed permanent acquired immunity, meaning
we have less understanding of how disease dynamics will play out when immunity
wanes. Here we complement these studies by developing and analysing a
stochastic simulation model of disease spread on a university campus where we
allow immunity to wane, expoloring the effectiveness of different
interventions. We find that the two most effective interventions to limit the
severity of a disease outbreak are reducing extra-household mixing and
surveillance testing backed-up by a moderate isolation period. We find that
contact tracing only has a limited effect, while reducing class sizes only has
much effect if extra-household mixing is already low. We identify a range of
measures that can not only limit an outbreak but prevent it entirely, and also
comment on the variation in measures of severity that emerge from our
stochastic simulations. We hope that our model may help in designing effective
strategies for universities in future disease outbreaks.


  Whenever applicable, the Stochastic Gradient Descent (SGD) has shown itself
to be unreasonably effective. Instead of underperforming and getting trapped in
local minima due to the batch noise, SGD leverages it to learn to generalize
better and find minima that are good enough for the entire dataset. This led to
numerous theoretical and experimental investigations, especially in the context
of Artificial Neural Networks (ANNs), leading to better machine learning
algorithms. However, SGD is not applicable in a non-differentiable setting,
leaving all that prior research off the table.
  In this paper, we show that a class of evolutionary algorithms (EAs) inspired
by the Gillespie-Orr Mutational Landscapes model for natural evolution is
formally equivalent to SGD in certain settings and, in practice, is well
adapted to large ANNs. We refer to such EAs as Gillespie-Orr EA class (GO-EAs)
and empirically show how an insight transfer from SGD can work for them. We
then show that for ANNs trained to near-optimality or in the transfer learning
setting, the equivalence also allows transferring the insights from the
Mutational Landscapes model to SGD.
  We then leverage this equivalence to experimentally show how SGD and GO-EAs
can provide mutual insight through examples of minima flatness, transfer
learning, and mixing of individuals in EAs applied to large models.


  In the study of therapeutic strategies for the treatment of cancer,
eco-evolutionary dynamics are of particular interest, since characteristics of
the tumour population, interaction with the environment and effects of the
treatment, influence the geometric and epigenetic characterization of the
tumour with direct consequences on the efficacy of the therapy and possible
relapses. In particular, when considering radiotherapy, oxygen concentration
plays a central role both in determining the effectiveness of the treatment and
the selective pressure due to hypoxia. We propose a mathematical model, settled
in the framework of epigenetically-structured population dynamics and
formulated in terms of systems of coupled non-linear integro-differential
equations, that aims to catch these phenomena and to provide a predictive tool
for the tumour mass evolution and therapeutic effects. The outcomes of the
simulations show how the model is able to explain the impact of environmental
selection and therapies on the evolution of the mass, motivating observed
dynamics such as relapses and therapeutic failures. Furthermore it offers a
first hint for the development of therapies which can be adapted to overcome
problems of resistance and relapses.


  Nucleic acids are highly deformable helical molecules constantly stretched,
twisted and bent in their biological functioning. Single molecule experiments
have shown that double stranded (ds)-RNA and standard ds-DNA have opposite
twist-stretch patterns and stretching properties when overwound under a
constant applied load. The key structural features of the A-form RNA and B-form
DNA helices are here incorporated in a three-dimensional mesoscopic Hamiltonian
model which accounts for the radial, bending and twisting fluctuations of the
base pairs. Using path integral techniques which sum over the ensemble of the
base pair fluctuations, I compute the average helical repeat of the molecules
as a function of the load. The obtained twist-stretch relations and stretching
properties, for short A- and B-helical fragments, are consistent with the
opposite behaviors observed in kilo-base long molecules.


  Urinary Tract Infection (UTIs) is referred as one of the most common
infection in medical sectors worldwide and antimicrobial resistance (AMR) is
also a global threat to human that is related with many diseases. As
antibiotics used for the treatment of infectious diseases, the rate of
resistance is increasing day by day. Gram positive pathogens are commonly found
in urine sample collected from different age groups of people, associated with
UTI. The study was conducted in a diagnostic center in Dhaka, Bangladesh with
total 1308 urine samples from November 2021 to April 2022. Gram positive
pathogens were isolated and antimicrobial susceptibility tests were done. From
total 121 samples of gram positive bacteria the highest prevalence rate of UTIs
was found in age group of 21-30 year. Mostly Enterococcus spp. (33.05%)
Staphylococcus aureus (27.27%), Streptococcus spp. (20.66%), Beta-hemolytic
streptococci (19.00%) were found as causative agents of UTI compared to others.
The majority of isolates have been detected as multi-drug resistant (MDR). The
higher percentage of antibiotic resistance were found against Azithromycin
(75%), and cefixime (64.46%). This research focused on the regular basis of
surveillance for the Gram-positive bacteria antibiotic susceptibility to
increase awareness about the use of proper antibiotic thus minimize the drug
resistance.


  We present a model, based on symmetry and geometry, for proteins. Using
elementary ideas from mathematics and physics, we derive the geometries of
discrete helices and sheets. We postulate a compatible solvent-mediated
emergent pairwise attraction that assembles these building blocks, while
respecting their individual symmetries. Instead of seeking to mimic the
complexity of proteins, we look for a simple abstraction of reality that yet
captures the essence of proteins. We employ analytic calculations and detailed
Monte Carlo simulations to explore some consequences of our theory. The
predictions of our approach are in accord with experimental data. Our framework
provides a rationalization for understanding the common characteristics of
proteins. Our results show that the free energy landscape of a globular protein
is pre-sculpted at the backbone level, sequences and functionalities evolve in
the fixed backdrop of the folds determined by geometry and symmetry, and that
protein structures are unique in being simultaneously characterized by
stability, diversity, and sensitivity.


  Recently, we presented a framework for understanding protein structure based
on the idea that simple constructs of holding hands or touching of objects can
be used to rationalize the common characteristics of globular proteins. We
developed a consistent approach for understanding the formation of the two key
common building blocks of helices and sheets as well as the compatible assembly
of secondary structures into the tertiary structure through the notion of
poking pairwise interactions. Here we benchmark our predictions with a detailed
analysis of structural data of over 4000 proteins from the Protein Data Bank.
We also present the results of detailed computer simulations of a simplified
model demonstrating a pre-sculpted free energy landscape, determined by
geometry and symmetry, comprising numerous minima corresponding to putative
native state structures. We explore the consequences of our model. Our results
suggest that symmetry and geometry are a powerful guide to capture the
simplicity underlying protein complexity.


  Natural microbial populations often have complex spatial structures. This can
impact their evolution, in particular the ability of mutants to take over.
While mutant fixation probabilities are known to be unaffected by sufficiently
symmetric structures, evolutionary graph theory has shown that some graphs can
amplify or suppress natural selection, in a way that depends on microscopic
update rules. We propose a model of spatially structured populations on graphs
directly inspired by batch culture experiments, alternating within-deme growth
on nodes and migration-dilution steps, and yielding successive bottlenecks.
This setting bridges models from evolutionary graph theory with Wright-Fisher
models. Using a branching process approach, we show that spatial structure with
frequent migrations can only yield suppression of natural selection. More
precisely, in this regime, circulation graphs, where the total incoming
migration flow equals the total outgoing one in each deme, do not impact
fixation probability, while all other graphs strictly suppress selection.
Suppression becomes stronger as the asymmetry between incoming and outgoing
migrations grows. Amplification of natural selection can nevertheless exist in
a restricted regime of rare migrations and very small fitness advantages, where
we recover the predictions of evolutionary graph theory for the star graph.


  Myosin II plays a pivotal role in muscle contraction by generating force
through the cooperative action of multiple motors on actin filaments. In this
study, we integrate the nonlinear elasticity of the neck linker in individual
myosin II and comprehensively investigate the evolution of cooperativity and
dynamics at {\it microstate} and {\it mesostate} levels using a combined model
of single and multiple motors. We find that a substantial proportion of
actin-bound motors reside in the {\it mid-} and {\it post-power stroke} states,
and our nonlinear model reveals their increased capacity for load sharing.
Additionally, we systematically explore the impact of mechanical load and ATP
concentration on myosin II motors. Notably, we observe that the average net
distance of actin undergoes a transition from a weak load-sensitive regime at
low ATP concentrations to a load-sensitive regime at higher ATP concentrations.
Furthermore, increasing the load or raising the ATP concentration to saturation
can enhance the efficiency and output power of myosin filament. Moreover, the
efficiency of the myosin filament increases with the power stroke strength,
reaching a maximum at a specific range, and subsequently declining beyond that
threshold. Finally, we explore the mean run time/length and mean existence
probability of myosin filament, shedding light on its overall behavior.


  Ocean acidification, a direct consequence of increased carbon dioxide (CO2)
emissions, has emerged as a critical area of concern within the scientific
community. The world's oceans absorb approximately one-third of human-caused
CO2 emissions, leading to chemical reactions that reduce seawater pH, carbonate
ion concentration, and saturation states of biologically important calcium
carbonate minerals. This process, known as ocean acidification, has
far-reaching implications for marine ecosystems, particularly for marine
organisms such as fish, whose migratory patterns are integral to the health and
function of these ecosystems.


  The BReast CAncer type 2 susceptibility protein (BRCA2) responds to DNA
damage by participating in homology-directed repair. BRCA2 deficiency
culminates in defective DNA damage repair (DDR) that when prolonged leads to
the accumulation of DNA damage causing cancer or apoptosis. Oxidative stress
promotes DNA damage and apoptosis and is a common mechanism through which
cardiovascular risk factors lead to endothelial dysfunction and
atherosclerosis. Herein, we show that endothelial BRCA2 plays a protective role
against atherosclerosis under hypercholesterolemic stress. We successfully
generated and characterized endothelial cell (EC)-specific BRCA2 knockout
(BRCA2endo) mice. To study the effect of EC-specific BRCA2-loss in
atherosclerosis, we generated and characterized BRCA2endo mice on
apolipoprotein E null background (ApoE-/-), fed them with high-fat diet (HFD)
and evaluated atherosclerosis. Baseline phenotyping of BRCA2endo mice did not
show any adverse effects in terms of DNA damage and apoptosis as well as
cardiac and metabolic function. However, using HFD-fed apolipoprotein E
knockout (ApoE-/-) background, we demonstrated that EC-specific loss of BRCA2
resulted in aortic plaque deposition and splenomegaly. Comparison of RNA
sequencing data from aortas of EC-specific BRCA2-deficient ApoE-/- and
BRCA2-intact ApoE-/- mice revealed a total of 530 significantly differentially
expressed genes with Protein Folding Response and Lipid Metabolism as the most
affected pathways. This study provides foundational knowledge regarding BRCA2
status and function in the cardiovascular system, and highlights the potential
of BRCA2 as a novel therapeutic target in prevention and treatment of
atherosclerosis. Our data indicate that BRCA2 mutation carriers may be at a
previously unrecognized risk of atherosclerosis in addition to breast and
ovarian cancer.


  The onset of the COVID-19 pandemic drove a widespread, often uncoordinated
effort by research groups to develop mathematical models of SARS-CoV-2 to study
its spread and inform control efforts. The urgent demand for insight at the
outset of the pandemic meant early models were typically either simple or
repurposed from existing research agendas. Our group predominantly uses
agent-based models (ABMs) to study fine-scale intervention scenarios. These
high-resolution models are large, complex, require extensive empirical data,
and are often more detailed than strictly necessary for answering qualitative
questions like "Should we lockdown?" During the early stages of an
extraordinary infectious disease crisis, particularly before clear empirical
evidence is available, simpler models are more appropriate. As more detailed
empirical evidence becomes available, however, and policy decisions become more
nuanced and complex, fine-scale approaches like ours become more useful.
  In this manuscript, we discuss how our group navigated this transition as we
modeled the pandemic. The role of modelers often included nearly real-time
analysis, and the massive undertaking of adapting our tools quickly. We were
often playing catch up with a firehose of evidence, while simultaneously
struggling to do both academic research and real-time decision support, under
conditions conducive to neither. By reflecting on our experiences of responding
to the pandemic and what we learned from these challenges, we can better
prepare for future demands.


  Advances in metagenomic methods have revealed an astonishing number and
diversity of microbial lifeforms, most of which are rare relative to the
most-abundant, dominant taxa. The ecological and evolutionary mechanisms that
generate and sustain broadly observed microbial diversity patterns remain
debated. One possibility is that complex interactions between numerous taxa are
a main driver of the composition of a microbial community. Lotka-Volterra
equations with disordered interactions between species offer a minimal yet rich
modelling framework to investigate this hypothesis. We consider communities
with strong, mostly competitive interactions, where species-rich coexistence
equilibria are typically unstable. When species extinction is prevented by a
small rate of immigration, one generically finds a sustained chaotic phase,
where all species participate in a continuous turnover of who is rare and who
is dominant. The distribution of rare species' abundances -- in a snapshot of
the whole community, and for each species individually in time -- follows a
distribution with a prominent power-law trend with exponent $\nu>1$. We
formulate a focal-species model in terms of a logistic growth equation with
coloured noise that reproduces dynamical features of the disordered
Lotka-Volterra model. With its use, we discover that $\nu$ is mainly determined
by three effective parameters of the dominant community, such as its timescale
of turnover. Approximate proportionalities between the effective parameters
constrain the variation of $\nu$ across the range of interaction statistics
resulting in chaotic turnover. We discuss our findings in the context of marine
plankton communities, where chaos, boom-bust dynamics, and a power-law
abundance distribution have been observed.


  Developing accurate and efficient coarse-grained representations of proteins
is crucial for understanding their folding, function, and interactions over
extended timescales. Our methodology involves simulating proteins with
molecular dynamics and utilizing the resulting trajectories to train a neural
network potential through differentiable trajectory reweighting. Remarkably,
this method requires only the native conformation of proteins, eliminating the
need for labeled data derived from extensive simulations or memory-intensive
end-to-end differentiable simulations. Once trained, the model can be employed
to run parallel molecular dynamics simulations and sample folding events for
proteins both within and beyond the training distribution, showcasing its
extrapolation capabilities. By applying Markov State Models, native-like
conformations of the simulated proteins can be predicted from the
coarse-grained simulations. Owing to its theoretical transferability and
ability to use solely experimental static structures as training data, we
anticipate that this approach will prove advantageous for developing new
protein force fields and further advancing the study of protein dynamics,
folding, and interactions.


  Critical birth-death processes with $n$ distinct types are investigated. Each
type $i$ cell divides independently $(i)\to(i)+(i)$ or mutates $(i)\to(i+1)$ at
the same rate. The total number of cells grows exponentially as a Yule process
until the maximal type $n$ cells appear, which cannot mutate but die at rate
one. The last type makes the process critical and hence after the exponentially
growing phase eventually all cells die with probability one. The process mimics
the accumulation of mutations in a growing population where too many mutations
are lethal. This has applications for understanding the mutational burden and
so-called error catastrophe in cancer, bacteria, or virus. We present
large-time asymptotic results for the general $n$-type critical birth-death
process. We find that the mass function of the number of cells of type $k$ has
algebraic and stationary tail $(\text{size})^{-1-\chi_k}$, with
$\chi_k=2^{1-k}$, for $k=2,\dots,n$, in sharp contrast to the exponential tail
of the first type. The same exponents describe the tail of the asymptotic
survival probability $(\text{time})^{-\chi_n}$. We discuss the consequences and
applications of the results for studying extinction due to mutational burden in
biological populations.


  Pool sequencing is an efficient method for capturing genome-wide allele
frequencies from multiple individuals, with broad applications such as studying
adaptation in Evolve-and-Resequence experiments, monitoring of genetic
diversity in wild populations, and genotype-to-phenotype mapping. Here, we
present grenedalf, a command line tool written in C++ that implements common
population genetic statistics such as $\theta$, Tajima's D, and FST for Pool
sequencing. It is orders of magnitude faster than current tools, and is focused
on providing usability and scalability, while also offering a plethora of input
file formats and convenience options.


  In the last decades, the capacity to generate large amounts of data in
science and engineering applications has been growing steadily. Meanwhile, the
progress in machine learning has turned it into a suitable tool to process and
utilise the available data. Nonetheless, many relevant scientific and
engineering problems present challenges where current machine learning methods
cannot yet efficiently leverage the available data and resources. For example,
in scientific discovery, we are often faced with the problem of exploring very
large, high-dimensional spaces, where querying a high fidelity, black-box
objective function is very expensive. Progress in machine learning methods that
can efficiently tackle such problems would help accelerate currently crucial
areas such as drug and materials discovery. In this paper, we propose the use
of GFlowNets for multi-fidelity active learning, where multiple approximations
of the black-box function are available at lower fidelity and cost. GFlowNets
are recently proposed methods for amortised probabilistic inference that have
proven efficient for exploring large, high-dimensional spaces and can hence be
practical in the multi-fidelity setting too. Here, we describe our algorithm
for multi-fidelity active learning with GFlowNets and evaluate its performance
in both well-studied synthetic tasks and practically relevant applications of
molecular discovery. Our results show that multi-fidelity active learning with
GFlowNets can efficiently leverage the availability of multiple oracles with
different costs and fidelities to accelerate scientific discovery and
engineering design.


  With the development of network science, Turing pattern has been proven to be
formed in discrete media such as complex networks, opening up the possibility
of exploring it as a generation mechanism in the context of biology, chemistry,
and physics. Turing instability in the predator-prey system has been widely
studied in recent years. We hope to use the predator-prey interaction
relationship in biological populations to explain the influence of network
topology on pattern formation. In this paper, we establish a predator-prey
model with weak Allee effect, analyze and verify the Turing instability
conditions on the large ER (Erd\"{o}s-R\'{e}nyi) random network with the help
of Turing stability theory and numerical experiments, and obtain the Turing
instability region. The results indicate that diffusion plays a decisive role
in the generation of spatial patterns, whether in continuous or discrete media.
For spatiotemporal patterns, different initial values can also bring about
changes in the pattern. When we analyze the model based on the network
framework, we find that the average degree of the network has an important
impact on the model, and different average degrees will lead to changes in the
distribution pattern of the population.


  Human T-cell leukemia virus type 1 (HTLV-1) causes adult T-cell leukemia
(ATL) and HTLV-1-associated myelopathy (HAM) after a long latent period in a
fraction of infected individuals. These HTLV-1-infected cells typically have
phenotypes similar to that of CD4${^+}$ T cells, but the cell status is not
well understood. To extract the inherent information of HTLV-1-infected CD4$^+$
cells, we integratively analyzed the ATAC-seq and RNA-seq data of infected
cells. Compared to CD4${^+}$ T cells from healthy donors, we found anomalous
chromatin accessibility in HTLV-1-infected CD4${^+}$ cells derived from ATL
cases in terms of location and sample-to-sample fluctuations in open chromatin
regions. Further, by focusing on systematically selected genes near the open
chromatin regions, all the gene expressions in ATL cases were found to be
distinct from those of healthy CD4$^+$ T cells. Based on a further analysis of
chromatin accessibility, we detected TLL1 (Tolloid Like 1) as one of the key
genes that exhibit unique gene expressions in ATL cases. A luciferase assay
indicated that TLL1 has a strong regulatory effect on TGF-$\beta$. Overall,
this study provides results about the status of HTLV-1 infected cells, which
are qualitatively consistent across the different scales of chromatin
accessibility, transcription, and immunophenotype.


  Natural language is expected to be a key medium for various human-machine
interactions in the era of large language models. When it comes to the
biochemistry field, a series of tasks around molecules (e.g., property
prediction, molecule mining, etc.) are of great significance while having a
high technical threshold. Bridging the molecule expressions in natural language
and chemical language can not only hugely improve the interpretability and
reduce the operation difficulty of these tasks, but also fuse the chemical
knowledge scattered in complementary materials for a deeper comprehension of
molecules. Based on these benefits, we propose the conversational molecular
design, a novel task adopting natural language for describing and editing
target molecules. To better accomplish this task, we design ChatMol, a
knowledgeable and versatile generative pre-trained model, enhanced by injecting
experimental property information, molecular spatial knowledge, and the
associations between natural and chemical languages into it. Several typical
solutions including large language models (e.g., ChatGPT) are evaluated,
proving the challenge of conversational molecular design and the effectiveness
of our knowledge enhancement method. Case observations and analysis are
conducted to provide directions for further exploration of natural-language
interaction in molecular discovery.


  Pre-trained models have been successful in many protein engineering tasks.
Most notably, sequence-based models have achieved state-of-the-art performance
on protein fitness prediction while structure-based models have been used
experimentally to develop proteins with enhanced functions. However, there is a
research gap in comparing structure- and sequence-based methods for predicting
protein variants that are better than the wildtype protein. This paper aims to
address this gap by conducting a comparative study between the abilities of
equivariant graph neural networks (EGNNs) and sequence-based approaches to
identify promising amino-acid mutations. The results show that our proposed
structural approach achieves a competitive performance to sequence-based
methods while being trained on significantly fewer molecules. Additionally, we
find that combining assay labelled data with structure pre-trained models
yields similar trends as with sequence pre-trained models.
  Our code and trained models can be found at:
https://github.com/semiluna/partIII-amino-acid-prediction.


  Possessing the five-year durability rate of nearly 5%, currently, the fourth
leading cause for cancer-related deaths is pancreatic cancer. Previously,
several works have resolved that early diagnosis performs a meaningful function
in enhancing the durability rate and diverse online tools have been utilized to
distinguish prognostic biomarker which is a lengthy process. We believe that
the statistical feature selection method can produce a better and faster result
here. To authenticate our statement, we picked three different mRNA microarray
(GSE15471, GSE28735, and GSE16515) and a microRNA (GSE41372) dataset for
identification of differentially expressed genes (DEGs) and differentially
expressed microRNAs (DEMs). By adopting some feature selecting methods, 178
DEGs and 16 DEMs were elected. After identifying target genes of DEMs, we
selected two DEGs (ECT2 and NRP2) which were also identified among DEMs target
genes. Moreover, overall durability report established that ECT2 and NRP2 were
associated with poor overall survival. Hence, we concluded that for pancreatic
cancer, statistical feature selection approaches certainly perform better for
biomarker identification than pre-defined online programs, and here, ECT2 and
NRP2 can act as possible prognostic biomarkers. All the resources, programs and
snippets of our literature can be discovered at
https://github.com/Srizon143005/PancreaticCancerBiomarkers.


  Voltage-gated sodium (Nav) channels constitute a prime target for drug design
and discovery, given their implication in various diseases such as epilepsy,
migraine and ataxia to name a few. In this regard, performing morphological
analysis is a crucial step in comprehensively understanding their biological
function and mechanism, as well as in uncovering subtle details of their
mechanism that may be elusive to experimental observations. Despite their
tremendous therapeutic potential, drug design resources are deficient,
particularly in terms of accurate and comprehensive geometric information. This
paper presents a geometric dataset of molecular surfaces that are
representative of Nav channels in mammals. For each structure we provide three
representations and a number of geometric measures, including length, volume
and straightness of the recognized channels. To demonstrate the effective use
of GEO-Nav, we have tested it on two methods belonging to two different
categories of approaches: a sphere-based and a tessellation-based method.


  We resolve difficulties in training and sampling from a discrete generative
model by learning a smoothed energy function, sampling from the smoothed data
manifold with Langevin Markov chain Monte Carlo (MCMC), and projecting back to
the true data manifold with one-step denoising. Our Discrete Walk-Jump Sampling
formalism combines the maximum likelihood training of an energy-based model and
improved sample quality of a score-based model, while simplifying training and
sampling by requiring only a single noise level. We evaluate the robustness of
our approach on generative modeling of antibody proteins and introduce the
distributional conformity score to benchmark protein generative models. By
optimizing and sampling from our models for the proposed distributional
conformity score, 97-100% of generated samples are successfully expressed and
purified and 35% of functional designs show equal or improved binding affinity
compared to known functional antibodies on the first attempt in a single round
of laboratory experiments. We also report the first demonstration of long-run
fast-mixing MCMC chains where diverse antibody protein classes are visited in a
single MCMC chain.


  Objective: COVID-19 has spread worldwide and made a huge influence across the
world. Modeling the infectious spread situation of COVID-19 is essential to
understand the current condition and to formulate intervention measurements.
Epidemiological equations based on the SEIR model simulate disease development.
The traditional parameter estimation method to solve SEIR equations could not
precisely fit real-world data due to different situations, such as social
distancing policies and intervention strategies. Additionally, learning-based
models achieve outstanding fitting performance, but cannot visualize
mechanisms. Methods: Thus, we propose a deep dynamic epidemiological (DDE)
method that combines epidemiological equations and deep-learning advantages to
obtain high accuracy and visualization. The DDE contains deep networks to fit
the effect function to simulate the ever-changing situations based on the
neural ODE method in solving variants' equations, ensuring the fitting
performance of multi-level areas. Results: We introduce four SEIR variants to
fit different situations in different countries and regions. We compare our DDE
method with traditional parameter estimation methods (Nelder-Mead, BFGS,
Powell, Truncated Newton Conjugate-Gradient, Neural ODE) in fitting the
real-world data in the cases of countries (the USA, Columbia, South Africa) and
regions (Wuhan in China, Piedmont in Italy). Our DDE method achieves the best
Mean Square Error and Pearson coefficient in all five areas. Further, compared
with the state-of-art learning-based approaches, the DDE outperforms all
techniques, including LSTM, RNN, GRU, Random Forest, Extremely Random Trees,
and Decision Tree. Conclusion: DDE presents outstanding predictive ability and
visualized display of the changes in infection rates in different regions and
countries.


  Over the past years, the ever-growing trend on data storage demand, more
specifically for "cold" data (i.e. rarely accessed), has motivated research for
alternative systems of data storage. Because of its biochemical
characteristics, synthetic DNA molecules are now considered as serious
candidates for this new kind of storage. This paper introduces a novel
arithmetic coder for DNA data storage, and presents some results on a lossy
JPEG 2000 based image compression method adapted for DNA data storage that uses
this novel coder.
  The DNA coding algorithms presented here have been designed to efficiently
compress images, encode them into a quaternary code, and finally store them
into synthetic DNA molecules. This work also aims at making the compression
models better fit the problematic that we encounter when storing data into DNA,
namely the fact that the DNA writing, storing and reading methods are error
prone processes.
  The main take away of this work is our arithmetic coder and it's integration
into a performant image codec.


  Protein loop modeling is the most challenging yet highly non-trivial task in
protein structure prediction. Despite recent progress, existing methods
including knowledge-based, ab initio, hybrid and deep learning (DL) methods
fall significantly short of either atomic accuracy or computational efficiency.
Moreover, an overarching focus on backbone atoms has resulted in a dearth of
attention given to side-chain conformation, a critical aspect in a host of
downstream applications including ligand docking, molecular dynamics simulation
and drug design. To overcome these limitations, we present KarmaLoop, a novel
paradigm that distinguishes itself as the first DL method centered on full-atom
(encompassing both backbone and side-chain heavy atoms) protein loop modeling.
Our results demonstrate that KarmaLoop considerably outperforms conventional
and DL-based methods of loop modeling in terms of both accuracy and efficiency,
with the average RMSD improved by over two-fold compared to the second-best
baseline method across different tasks, and manifests at least two orders of
magnitude speedup in general. Consequently, our comprehensive evaluations
indicate that KarmaLoop provides a state-of-the-art DL solution for protein
loop modeling, with the potential to hasten the advancement of protein
engineering, antibody-antigen recognition, and drug design.


  Recent research on predicting the binding affinity between drug molecules and
proteins use representations learned, through unsupervised learning techniques,
from large databases of molecule SMILES and protein sequences. While these
representations have significantly enhanced the predictions, they are usually
based on a limited set of modalities, and they do not exploit available
knowledge about existing relations among molecules and proteins. In this study,
we demonstrate that by incorporating knowledge graphs from diverse sources and
modalities into the sequences or SMILES representation, we can further enrich
the representation and achieve state-of-the-art results for drug-target binding
affinity prediction in the established Therapeutic Data Commons (TDC)
benchmarks. We release a set of multimodal knowledge graphs, integrating data
from seven public data sources, and containing over 30 million triples. Our
intention is to foster additional research to explore how multimodal knowledge
enhanced protein/molecule embeddings can improve prediction tasks, including
prediction of binding affinity. We also release some pretrained models learned
from our multimodal knowledge graphs, along with source code for running
standard benchmark tasks for prediction of biding affinity.


  Molecule property prediction has gained significant attention in recent
years. The main bottleneck is the label insufficiency caused by expensive lab
experiments. In order to alleviate this issue and to better leverage textual
knowledge for tasks, this study investigates the feasibility of employing
natural language instructions to accomplish molecule-related tasks in a
zero-shot setting. We discover that existing molecule-text models perform
poorly in this setting due to inadequate treatment of instructions and limited
capacity for graphs. To overcome these issues, we propose GIMLET, which unifies
language models for both graph and text data. By adopting generalized position
embedding, our model is extended to encode both graph structures and
instruction text without additional graph encoding modules. GIMLET also
decouples encoding of the graph from tasks instructions in the attention
mechanism, enhancing the generalization of graph features across novel tasks.
We construct a dataset consisting of more than two thousand molecule tasks with
corresponding instructions derived from task descriptions. We pretrain GIMLET
on the molecule tasks along with instructions, enabling the model to transfer
effectively to a broad range of tasks. Experimental results demonstrate that
GIMLET significantly outperforms molecule-text baselines in instruction-based
zero-shot learning, even achieving closed results to supervised GNN models on
tasks such as toxcast and muv.


  The honeybee plays an extremely important role in ecosystem stability and
diversity and in the production of bee pollinated crops. Honey bees and other
pollinators are under threat from the combined effects of nutritional stress,
parasitism, pesticides, and climate change that impact the timing, duration,
and variability of seasonal events. To understand how parasitism and
seasonality influence honey bee colonies separately and interactively, we
developed a non-autonomous nonlinear honeybee-parasite interaction differential
equation model that incorporates seasonality into the egg-laying rate of the
queen. Our theoretical results show that parasitism negatively impacts the
honey bee population either by decreasing colony size or destabilizing
population dynamics through supercritical or subcritical Hopf-bifurcations
depending on conditions. Our bifurcation analysis and simulations suggest that
seasonality alone may have positive or negative impacts on the survival of
honey bee colonies. More specifically, our study indicates that (1) the timing
of the maximum egg-laying rate seems to determine when seasonality has positive
or negative impacts; and (2) when the period of seasonality is large it can
lead to the colony collapsing. Our study further suggests that the synergistic
influences of parasitism and seasonality can lead to complicated dynamics that
may positively and negatively impact the honey bee colony's survival. Our work
partially uncovers the intrinsic effects of climate change and parasites, which
potentially provide essential insights into how best to maintain or improve a
honey bee colony's health.


  Spin-glasses constitute a well-grounded framework for evolutionary models. Of
particular interest for (some of) these models is the lack of self-averaging of
their order parameters (e.g. the Hamming distance between the genomes of two
individuals), even in asymptotic limits, much as like the behavior of the
overlap between the configurations of two replica in mean-field spin-glasses.
In the latter, this lack of self-averaging is related to peculiar fluctuations
of the overlap, known as Ghirlanda-Guerra identities and Aizenman-Contucci
polynomials, that cover a pivotal role in describing the ultrametric structure
of the spin-glass landscape. As for evolutionary models, such identities may
therefore be related to a taxonomic classification of individuals, yet a full
investigation on their validity is missing. In this paper, we study ultrametric
identities in simple cases where solely random mutations take place, while
selective pressure is absent, namely in {\em flat landscape} models. In
particular, we study three paradigmatic models in this setting: the {\em one
parent model} (which, by construction, is ultrametric at the level of single
individuals), the {\em homogeneous population model} (which is replica
symmetric), and the {\em species formation model} (where a broken-replica
scenario emerges at the level of species). We find analytical and numerical
evidence that in the first and in the third model nor the Ghirlanda-Guerra
neither the Aizenman-Contucci constraints hold, rather a new class of
ultrametric identities is satisfied; in the second model all these constraints
hold trivially. Very preliminary results on a real biological human genome
derived by {\em The 1000 Genome Project Consortium} and on two artificial human
genomes (generated by two different types neural networks) seem in better
agreement with these new identities rather than the classic ones.


  In this article, we have considered a planar slow-fast modified Leslie-Gower
predator-prey model with a weak Allee effect in the predator, based on the
natural assumption that the prey reproduces far more quickly than the predator.
We present a thorough mathematical analysis demonstrating the existence of
homoclinic orbits, heteroclinic orbits, singular Hopf bifurcation, canard limit
cycles, relaxation oscillations, the birth of canard explosion by combining the
normal form theory of slow-fast systems, Fenichel's theorem and blow-up
technique near non-hyperbolic point. We have obtained very rich dynamical
phenomena of the model, including the saddle-node, Hopf, transcritical
bifurcation, generalized Hopf, cusp point, homoclinic orbit, heteroclinic
orbit, and Bogdanov-Takens bifurcations. Moreover, we have investigated the
global stability of the unique positive equilibrium, as well as bistability,
which shows that the system can display either 'prey extinction', 'stable
coexistence', or 'oscillating coexistence' depending on the initial population
size and values of the system parameters. The theoretical findings are verified
by numerical simulations.


  It is well known that the confirmed COVID-19 infection is only a fraction of
the true fraction. In this paper we use an artificial neural network to learn
the connection between the confirmed infection count, the testing data, and the
true infection count. The true infection count in the training set is obtained
by backcasting from the death count and the infection fatality ratio (IFR).
Multiple factors are taken into consideration in the estimation of IFR. We also
calibrate the recovered true COVID-19 case count with an SEIR model.


  We investigate the Generalized Lotka-Volterra (GLV) equations, a central
model in theoretical ecology, where species interactions are assumed to be
fixed over time and heterogeneous (quenched noise). Recent studies have
suggested that the stability properties and abundance distributions of large
disordered GLV systems depend, in the simplest scenario, solely on the mean and
variance of the distribution of species interactions. However, empirical
communities deviate from this level of universality. In this article, we
present a generalized version of the dynamical mean field theory for
non-Gaussian interactions that can be applied to various models, including the
GLV equations. Our results show that the generalized mean field equations have
solutions which depend on all cumulants of the distribution of species
interactions, leading to a breakdown of universality. We leverage on this
informative breakdown to extract microscopic interaction details from the
macroscopic distribution of densities which are in agreement with empirical
data. Specifically, in the case of sparse interactions, which we analytically
investigate, we establish a simple relationship between the distribution of
interactions and the distribution of species population densities.


  The theory of alternative stable states and tipping points has garnered a lot
of attention in the last decades. It predicts potential critical transitions
from one ecosystem state to a completely different state under increasing
environmental stress. However, typically ecosystem models that predict tipping
do not resolve space explicitly. As ecosystems are inherently spatial, it is
important to understand the effects of incorporating spatial processes in
models, and how those insights translate to the real world. Moreover, spatial
ecosystem structures, such as vegetation patterns, are important in the
prediction of ecosystem response in the face of environmental change. Models
and observations from real savanna ecosystems and drylands have suggested that
they may exhibit both tipping behavior as well as spatial pattern formation.
Hence, in this paper, we use mathematical models of humid savannas and drylands
to illustrate several pattern formation phenomena that may arise when
incorporating spatial dynamics in models that exhibit tipping without resolving
space. We argue that such mechanisms challenge the notion of large scale
critical transitions in response to global change and reveal a more resilient
nature of spatial ecosystems.


  Numerous epidemiological studies suggest febrile infections could confer
long-term immunity to certain types of cancers, though the precise mechanisms
for this phenomenon remain unclear. Systemic heat-shock responses to fever may
be key to understanding the overlapping outcomes of immune responses to
infection and cancer. To investigate this hypothesis, we performed epitope
discovery between heat-shock proteins (HSP) and cancer-associated antigens
(CAA) and annotated the results with experimentally validated epitopes in the
Immune Epitope Database (IEDB) (Vita et al., 2019). Further, epitopes were
matched with their homologs in human pathogens. Results identified 94 epitopes
shared between HSPs and CAAs, with experimental evidence of presentation at MHC
molecules and with high homology to several epitopes of human pathogens. The
identified epitopes can be used as candidates for designing cancer vaccines.
They may also be used to identify autoreactive antibodies or TCR specificities
that, as antibody drugs and cell therapies, would reproduce the effect of
febrile infection in conferring cancer immunity. Our results support the
hypothesis that the loss of self-tolerance to HSPs during febrile infection
confers tumor immunity through molecular mimicry.


  We examine here the effects of recurrent vaccination and waning immunity on
the establishment of an endemic equilibrium in a population. An
individual-based model that incorporates memory effects for transmission rate
during infection and subsequent immunity is introduced, considering
stochasticity at the individual level. By letting the population size going to
infinity, we derive a set of equations describing the large scale behavior of
the epidemic. The analysis of the model's equilibria reveals a criterion for
the existence of an endemic equilibrium, which depends on the rate of immunity
loss and the distribution of time between booster doses. The outcome of a
vaccination policy in this context is influenced by the efficiency of the
vaccine in blocking transmissions and the distribution pattern of booster doses
within the population. Strategies with evenly spaced booster shots at the
individual level prove to be more effective in preventing disease spread
compared to irregularly spaced boosters, as longer intervals without
vaccination increase susceptibility and facilitate more efficient disease
transmission. We provide an expression for the critical fraction of the
population required to adhere to the vaccination policy in order to eradicate
the disease, that resembles a well-known threshold for preventing an outbreak
with an imperfect vaccine. We also investigate the consequences of unequal
vaccine access in a population and prove that, under reasonable assumptions,
fair vaccine allocation is the optimal strategy to prevent endemicity.


  Highly-diverse ecosystems exhibit a broad distribution of population sizes
and species turnover, where species at high and low abundances are exchanged
over time. We show that these two features generically emerge in the
fluctuating phase of many-variable model ecosystems with disordered species
interactions, when species are supported by migration from outside the system
at a small rate. We show that these and other phenomena can be understood
through the existence of a scaling regime in the limit of small migration, in
which large fluctuations and long timescales emerge. We construct an exact
analytical theory for this asymptotic regime, that provides scaling predictions
on timescales and abundance distributions that are verified exactly in
simulations. In this regime, a clear separation emerges between rare and
abundant species at any given time, despite species moving back and forth
between the rare and abundant subsets. The number of abundant species is found
to lie strictly below a well-known stability bound, maintaining the system away
from marginality. At the same time, other measures of diversity, which also
include some of the rare species, go above this bound.
  In the asymptotic limit where the migration rate goes to zero, trajectories
of individual species abundances are described by non-Markovian jump-diffusion
processes, which proceeds as follows: A rare species remains so for some time,
then experiences a jump in population sizes after which it becomes abundant (a
species turnover event) and later sees its population size gradually decreasing
again until rare, due to the competition with other species. The asymmetry of
abundance trajectories under time-reversal is maintained at small but finite
migration rate. These features may serve as fingerprints of endogenous
fluctuations in highly-diverse ecosystems.


  In the current manuscript, we consider a predator-prey model where the
predator is modeled as a generalist using a modified Leslie-Gower scheme, and
the prey exhibits group defence via a generalised response. We show that the
model could exhibit finite time blow-up, contrary to the current literature
(Eur. Phys. J. Plus 137, 28). We also propose a new concept via which the
predator population blows up in finite time while the prey population quenches
in finite time. The blow-up and quenching times are proved to be one and the
same. Our analysis is complemented by numerical findings. This includes a
numerical description of the basin of attraction for large data blow-up
solutions, as well as several rich bifurcations leading to multiple limit
cycles, both in co-dimension one and two. Lastly, we posit a delayed version of
the model with globally existing solutions for any initial data.


  Geometric deep learning has demonstrated a great potential in non-Euclidean
data analysis. The incorporation of geometric insights into learning
architecture is vital to its success. Here we propose a curvature-enhanced
graph convolutional network (CGCN) for biomolecular interaction prediction, for
the first time. Our CGCN employs Ollivier-Ricci curvature (ORC) to characterize
network local structures and to enhance the learning capability of GCNs. More
specifically, ORCs are evaluated based on the local topology from node
neighborhoods, and further used as weights for the feature aggregation in
message-passing procedure. Our CGCN model is extensively validated on fourteen
real-world bimolecular interaction networks and a series of simulated data. It
has been found that our CGCN can achieve the state-of-the-art results. It
outperforms all existing models, as far as we know, in thirteen out of the
fourteen real-world datasets and ranks as the second in the rest one. The
results from the simulated data show that our CGCN model is superior to the
traditional GCN models regardless of the positive-to-negativecurvature ratios,
network densities, and network sizes (when larger than 500).


  In recent years, AI-assisted drug design methods have been proposed to
generate molecules given the pockets' structures of target proteins. Most of
them are atom-level-based methods, which consider atoms as basic components and
generate atom positions and types. In this way, however, it is hard to generate
realistic fragments with complicated structures. To solve this, we propose
D3FG, a functional-group-based diffusion model for pocket-specific molecule
generation and elaboration. D3FG decomposes molecules into two categories of
components: functional groups defined as rigid bodies and linkers as mass
points. And the two kinds of components can together form complicated fragments
that enhance ligand-protein interactions.
  To be specific, in the diffusion process, D3FG diffuses the data distribution
of the positions, orientations, and types of the components into a prior
distribution; In the generative process, the noise is gradually removed from
the three variables by denoisers parameterized with designed equivariant graph
neural networks. In the experiments, our method can generate molecules with
more realistic 3D structures, competitive affinities toward the protein
targets, and better drug properties. Besides, D3FG as a solution to a new task
of molecule elaboration, could generate molecules with high affinities based on
existing ligands and the hotspots of target proteins.


  Drug-target interaction (DTI) prediction, which aims at predicting whether a
drug will be bounded to a target, have received wide attention recently, with
the goal to automate and accelerate the costly process of drug design. Most of
the recently proposed methods use single drug-drug similarity and target-target
similarity information for DTI prediction, which are unable to take advantage
of the abundant information regarding various types of similarities between
them. Very recently, some methods are proposed to leverage multi-similarity
information, however, they still lack the ability to take into consideration
the rich topological information of all sorts of knowledge bases where the
drugs and targets reside in. More importantly, the time consumption of these
approaches is very high, which prevents the usage of large-scale network
information. We thus propose a network-based drug-target interaction prediction
approach, which applies probabilistic soft logic (PSL) to meta-paths on a
heterogeneous network that contains multiple sources of information, including
drug-drug similarities, target-target similarities, drug-target interactions,
and other potential information. Our approach is based on the PSL graphical
model and uses meta-path counts instead of path instances to reduce the number
of rule instances of PSL. We compare our model against five methods, on three
open-source datasets. The experimental results show that our approach
outperforms all the five baselines in terms of AUPR score and AUC score.


  We consider parametrized systems of generalized polynomial inequalities (with
real exponents) in $n$ positive variables $x$ and involving $m$ monomial terms
$c \circ x^B$ (with exponent matrix $B$, positive parameters $c$, and
component-wise product $\circ$). More generally, we study solutions $x \in
\mathbb{R}^n_>$ to $(c \circ x^B) \in C$ for a "coefficient cone" $C \subseteq
\mathbb{R}^m_>$. Our framework encompasses systems of generalized polynomial
equations, as studied in real fewnomial and reaction network theory.
  We identify the relevant geometric objects of the problem, namely a bounded
set $P$ arising from the coefficient cone and two subspaces representing
monomial differences and dependencies. The dimension of the latter subspace
(the monomial dependency $d$) is crucial.
  As our main result, we rewrite the problem in terms of $d$ binomial
conditions on the "coefficient set" $P$, involving $d$ monomials in the
parameters. In particular, we obtain a classification of polynomial inequality
systems. If $d=0$, solutions exist (for all $c$) and can be parametrized
explicitly, thereby generalizing monomial parametrizations (of the solutions).
Even if $d>0$, solutions on the coefficient set can often be determined more
easily than solutions of the original system.
  Our framework allows a unified treatment of multivariate polynomial
inequalities, based on linear algebra and convex geometry. We demonstrate its
novelty and relevance in three examples from real fewnomial and reaction
network theory. For two mass-action systems, we parametrize the set of
equilibria and the region for multistationarity, respectively, and for
univariate trinomials, we provide a "solution formula" involving discriminants
and "roots".


  Advances in deep generative models shed light on de novo molecule generation
with desired properties. However, molecule generation targeted for dual protein
targets still faces formidable challenges including protein 3D structure data
requisition for model training, auto-regressive sampling, and model
generalization for unseen targets. Here, we proposed DiffDTM, a novel
conditional structure-free deep generative model based on a diffusion model for
dual targets based molecule generation to address the above issues.
Specifically, DiffDTM receives protein sequences and molecular graphs as inputs
instead of protein and molecular conformations and incorporates an information
fusion module to achieve conditional generation in a one-shot manner. We have
conducted comprehensive multi-view experiments to demonstrate that DiffDTM can
generate drug-like, synthesis-accessible, novel, and high-binding affinity
molecules targeting specific dual proteins, outperforming the state-of-the-art
(SOTA) models in terms of multiple evaluation metrics. Furthermore, we utilized
DiffDTM to generate molecules towards dopamine receptor D2 and
5-hydroxytryptamine receptor 1A as new antipsychotics. The experimental results
indicate that DiffDTM can be easily plugged into unseen dual targets to
generate bioactive molecules, addressing the issues of requiring insufficient
active molecule data for training as well as the need to retrain when
encountering new targets.


  For population systems modeled by age-structured hyperbolic partial
differential equations (PDEs), we redesign the existing feedback laws, designed
under the assumption that the dilution input is directly actuated, to the more
realistic case where dilution is governed by actuation dynamics (modeled simply
by an integrator). In addition to the standard constraint that the population
density must remain positive, the dilution dynamics introduce constraints of
not only positivity of dilution, but possibly of given positive lower and upper
bounds on dilution. We present several designs, of varying complexity, and with
various measurement requirements, which not only ensure global asymptotic (and
local exponential) stabilization of a desired positive population density
profile from all positive initial conditions, but do so without violating the
constraints on the dilution state. To develop the results, we exploit the
relation between first-order hyperbolic PDEs and an equivalent representation
in which a scalar input-driven mode is decoupled from input-free
infinite-dimensional internal dynamics represented by an integral delay system.


  Genome-wide association studies (GWAS) are commonly employed to study the
genetic basis of complex traits and diseases, and a key question is how much
heritability could be explained by all variants in GWAS. One widely used
approach that relies on summary statistics only is LD score regression (LDSC),
however the approach requires certain assumptions on the SNP effects (all SNPs
contribute to heritability and each SNP contributes equal variance). More
flexible modeling methods may be useful. We previously developed an approach
recovering the true z-statistics from a set of observed z-statistics with an
empirical Bayes approach, using only summary statistics. However, methods for
standard error (SE) estimation are not available yet, limiting the
interpretation of results and applicability of the approach. In this study we
developed several resampling-based approaches to estimate the SE of SNP-based
heritability, including two jackknife and three parametric bootstrap methods.
Simulations showed that delete-d-jackknife and parametric bootstrap approaches
provide good estimates of the SE. Particularly, the parametric bootstrap
approaches yield the lowest root-mean-squared-error (RMSE) of the true SE. In
addition, we applied our method to estimate SNP-based heritability of 12
immune-related traits (levels of cytokines and growth factors) to shed light on
their genetic architecture. We also implemented the methods to compute the sum
of heritability explained and the corresponding SE in an R package SumVg,
available at https://github.com/lab-hcso/Estimating-SE-of-total-heritability/ .
In conclusion, SumVg may provide a useful alternative tool for SNP heritability
and SE estimates, which does not rely on distributional assumptions of SNP
effects.


  Malaria is a parasitic disease, transmitted by mosquito vectors. Plasmodium
vivax presents particular challenges for disease control, in light of an
undetectable reservoir of latent parasites (hypnozoites) within the host liver.
Superinfection, which is driven by temporally proximate mosquito inoculation
and/or hypnozoite activation events, is an important feature of P. vivax. Here,
we present a model of hypnozoite accrual and superinfection for P. vivax. To
couple host and vector dynamics, we construct a density-dependent Markov
population process with countably many types, for which disease extinction is
shown to occur almost surely. We also establish a functional law of large
numbers, taking the form of an infinite-dimensional system of ordinary
differential equations that can also be recovered under the hybrid
approximation or a standard compartment modelling approach. Recognising that
the subset of these equations that models the infection status of human hosts
has precisely the same form as the Kolmogorov forward equations for a Markovian
network of infinite server queues with an inhomogeneous batch arrival process,
we use physical insight into the evolution of the latter to write down a
time-dependent multivariate generating function for the solution. We use this
characterisation to collapse the infinite-compartment model into a single
integrodifferential equation (IDE) governing the intensity of mosquito-to-human
transmission. Through a steady state analysis, we recover a threshold
phenomenon for this IDE in terms of a bifurcation parameter $R_0$, with the
disease-free equilibrium shown to be uniformly asymptotically stable if $R_0<1$
and an endemic equilibrium solution emerging if $R_0>1$. Our work provides a
theoretical basis to explore the epidemiology of P. vivax, and introduces a
general strategy for constructing tractable population-level models of malarial
superinfection.


  Diverse and complex intervention policies deployed over the last years have
shown varied effectiveness in controlling the COVID-19 pandemic. However, a
systematic analysis and modelling of the combined effects of different viral
lineages and complex intervention policies remains a challenge. Using
large-scale agent-based modelling and a high-resolution computational
simulation matching census-based demographics of Australia, we carried out a
systematic comparative analysis of several COVID-19 pandemic scenarios. The
scenarios covered two most recent Australian census years (2016 and 2021),
three variants of concern (ancestral, Delta and Omicron), and five
representative intervention policies. In addition, we introduced pandemic
Lorenz curves measuring an unequal distribution of the pandemic severity across
local areas. We quantified nonlinear effects of population heterogeneity on the
pandemic severity, highlighting that (i) the population growth amplifies
pandemic peaks, (ii) the changes in population size amplify the peak incidence
more than the changes in density, and (iii) the pandemic severity is
distributed unequally across local areas. We also examined and delineated the
effects of urbanisation on the incidence bimodality, distinguishing between
urban and regional pandemic waves. Finally, we quantified and examined the
impact of school closures, complemented by partial interventions, and
identified the conditions when inclusion of school closures may decisively
control the transmission. Our results suggest that (a) public health response
to long-lasting pandemics must be frequently reviewed and adapted to
demographic changes, (b) in order to control recurrent waves, mass-vaccination
rollouts need to be complemented by partial NPIs, and (c) healthcare and
vaccination resources need to be prioritised towards the localities and regions
with high population growth and/or high density.


  Molecular conformer generation (MCG) is an important task in cheminformatics
and drug discovery. The ability to efficiently generate low-energy 3D
structures can avoid expensive quantum mechanical simulations, leading to
accelerated virtual screenings and enhanced structural exploration. Several
generative models have been developed for MCG, but many struggle to
consistently produce high-quality conformers. To address these issues, we
introduce CoarsenConf, which coarse-grains molecular graphs based on torsional
angles and integrates them into an SE(3)-equivariant hierarchical variational
autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained
atomic coordinates of subgraphs connected via rotatable bonds, creating a
variable-length coarse-grained latent representation. Our model uses a novel
aggregated attention mechanism to restore fine-grained coordinates from the
coarse-grained latent representation, enabling efficient generation of accurate
conformers. Furthermore, we evaluate the chemical and biochemical quality of
our generated conformers on multiple downstream applications, including
property prediction and oracle-based protein docking. Overall, CoarsenConf
generates more accurate conformer ensembles compared to prior generative
models.


  Generation of molecules with desired chemical and biological properties such
as high drug-likeness, high binding affinity to target proteins, is critical
for drug discovery. In this paper, we propose a probabilistic generative model
to capture the joint distribution of molecules and their properties. Our model
assumes an energy-based model (EBM) in the latent space. Conditional on the
latent vector, the molecule and its properties are modeled by a molecule
generation model and a property regression model respectively. To search for
molecules with desired properties, we propose a sampling with gradual
distribution shifting (SGDS) algorithm, so that after learning the model
initially on the training data of existing molecules and their properties, the
proposed algorithm gradually shifts the model distribution towards the region
supported by molecules with desired values of properties. Our experiments show
that our method achieves very strong performances on various molecule design
tasks.