  Previous investigations into creative and innovation networks have suggested
that innovations often occurs at the boundary between the network's core and
periphery. In this work, we investigate the effect of global core-periphery
network structure on the speed and quality of cultural innovation. Drawing on
differing notions of core-periphery structure from [arXiv:1808.07801] and
[doi:10.1016/S0378-8733(99)00019-2], we distinguish decentralized
core-periphery, centralized core-periphery, and affinity network structure. We
generate networks of these three classes from stochastic block models (SBMs),
and use them to run an agent-based model (ABM) of collective cultural
innovation, in which agents can only directly interact with their network
neighbors. In order to discover the highest-scoring innovation, agents must
discover and combine the highest innovations from two completely parallel
technology trees. We find that decentralized core-periphery networks outperform
the others by finding the final crossover innovation more quickly on average.
We hypothesize that decentralized core-periphery network structure accelerates
collective problem-solving by shielding peripheral nodes from the local optima
known by the core community at any given time. We then build upon the "Two
Truths" hypothesis regarding community structure in spectral graph embeddings,
first articulated in [arXiv:1808.07801], which suggests that the adjacency
spectral embedding (ASE) captures core-periphery structure, while the Laplacian
spectral embedding (LSE) captures affinity. We find that, for core-periphery
networks, ASE-based resampling best recreates networks with similar performance
on the innovation SBM, compared to LSE-based resampling. Since the Two Truths
hypothesis suggests that ASE captures core-periphery structure, this result
further supports our hypothesis.


  Predicting the binding sites of the target proteins plays a fundamental role
in drug discovery. Most existing deep-learning methods consider a protein as a
3D image by spatially clustering its atoms into voxels and then feed the
voxelized protein into a 3D CNN for prediction. However, the CNN-based methods
encounter several critical issues: 1) defective in representing irregular
protein structures; 2) sensitive to rotations; 3) insufficient to characterize
the protein surface; 4) unaware of data distribution shift. To address the
above issues, this work proposes EquiPocket, an E(3)-equivariant Graph Neural
Network (GNN) for binding site prediction. In particular, EquiPocket consists
of three modules: the first one to extract local geometric information for each
surface atom, the second one to model both the chemical and spatial structure
of the protein, and the last one to capture the geometry of the surface via
equivariant message passing over the surface atoms. We further propose a dense
attention output layer to better alleviate the data distribution shift effect
incurred by the variable protein size. Extensive experiments on several
representative benchmarks demonstrate the superiority of our framework to the
state-of-the-art methods.


  Since its foundations, more than one hundred years ago, the field of
structural biology has strived to understand and analyze the properties of
molecules and their interactions by studying the structure that they take in 3D
space. However, a fundamental challenge with this approach has been the dynamic
nature of these particles, which forces us to model not a single but a whole
distribution of structures for every molecular system. This thesis proposes
Intrinsic Diffusion Modeling, a novel approach to this problem based on
combining diffusion generative models with scientific knowledge about the
flexibility of biological complexes. The knowledge of these degrees of freedom
is translated into the definition of a manifold over which the diffusion
process is defined. This manifold significantly reduces the dimensionality and
increases the smoothness of the generation space allowing for significantly
faster and more accurate generative processes. We demonstrate the effectiveness
of this approach on two fundamental tasks at the basis of computational
chemistry and biology: molecular conformer generation and molecular docking. In
both tasks, we construct the first deep learning method to outperform
traditional computational approaches achieving an unprecedented level of
accuracy for scalable programs.


  \textit{Plasmodium vivax} is the most geographically widespread
malaria-causing parasite resulting in significant associated global morbidity
and mortality. One of the factors driving this widespread phenomenon is the
ability of the parasites to remain dormant in the liver. Known as hypnozoites,
they reside in the liver following an initial exposure, before activating later
to cause further infections, referred to as relapses. As around 79-96$\%$ of
infections are attributed to relapses, we expect it will be highly impactful to
apply treatment to target the hypnozoite reservoir to eliminate \textit{P.
vivax}. Treatment with a radical cure to target the hypnozoite reservoir is a
potential tool to control or eliminate \textit{P. vivax}. We have developed a
multiscale mathematical model as a system of integro-differential equations
that captures the complex dynamics of \textit{P. vivax} hypnozoites and the
effect of hypnozoite relapse on disease transmission. Here, we use our model to
study the anticipated effect of radical cure treatment administered via a mass
drug administration (MDA) program. We implement multiple rounds of MDA with a
fixed interval between rounds, starting from different steady-state disease
prevalences. We then construct an optimisation model to obtain the optimal MDA
interval. We also incorporate mosquito seasonality in our model to study its
effect on the optimal treatment regime. We find that the effect of MDA
interventions is temporary and depends on the pre-intervention disease
prevalence (and choice of model parameters) as well as the number of MDA rounds
under consideration. We find radical cure alone may not be enough to lead to
\textit{P. vivax} elimination under our mathematical model (and choice of model
parameters) since the prevalence of infection eventually returns to pre-MDA
levels.


  PIFE was first used as an acronym for protein-induced fluorescence
enhancement, which refers to the increase in fluorescence observed upon the
interaction of a fluorophore, such as a cyanine, with a protein. This
fluorescence enhancement is due to changes in the rate of cis/trans
photoisomerisation. It is clear now that this mechanism is generally applicable
to interactions with any biomolecule and, in this review, we propose that PIFE
is thereby renamed according to its fundamental working principle as
photoisomerisation-related fluorescence enhancement, keeping the PIFE acronym
intact. We discuss the photochemistry of cyanine fluorophores, the mechanism of
PIFE, its advantages and limitations, and recent approaches to turn PIFE into a
quantitative assay. We provide an overview of its current applications to
different biomolecules and discuss potential future uses, including the study
of protein-protein interactions, protein-ligand interactions and conformational
changes in biomolecules.


  Lycopene contributes to the red-colored pigmentation of fruits and
vegetables, and it is a fat-soluble carotenoid with antioxidant properties.
Epidemiological studies have shown the significant health benefits associated
to the consumption of lycopene rich foods, because of the anti-cancer
properties. Degradative losses of lycopene during processing is a grave
concern, hence encapsulation provides a remedy. The objective of this study is
to evaluate the encapsulation efficiency of two biodegradable polymers (PLGA
and PCL) as for controlled release of lycopene in the gastrointestinal (GI)
system. The nanoparticles (NP) were synthesized by emulsion evaporation method
and physicochemical properties was determined using a Dynamic Light Scattering
spectroscopy. The results show the hydrodynamic diameter of the lycopene NP
synthesized in PCL (200 mg) and 3500 mg surfactant and sonicated for 15 min was
79.23+0.85 nm (Lowest). PLGA (500 mg) and 500 mg surfactant with 15 min
sonication was observed to have the lowest NP diameter (108.2+2.66 nm) among
the others. Significant difference result found in PDI value (0.12+0.07) when
PCL of 200mg dissolve in 3500 mg of surfactant. On the other hands the zeta
potential values were much smaller in case of PCL NP ranged between -1.3+0.046
and -4.21+0.08 mV compared to the PLGA NP -72.36+2.17 to 107.66+3.15 mV in all
experiments. Thus, NP synthesized with PCL and surfactant provide a smaller
sized nano-solution than PLGA and surfactant. As the degradation rate for PCL
is lower than PLGA so PCL can be considered as a potential biodegradable
polymer than PLGA to encapsulate lycopene.


  The identification of genomic, molecular and clinical markers predictive of
patient survival is important for developing personalized disease prevention,
diagnostic and treatment approaches. Modern omics technologies have made it
possible to investigate the prognostic impact of markers at multiple molecular
levels, including genomics, epigenomics (e.g. DNA methylation),
transcriptomics, proteomics and metabolomics, and how these potential risk
factors complement clinical characterization of patients for survival
prognosis. However, the massive sizes of the omics data sets pose challenges
for studying relationships between the molecular information and patients'
survival outcomes. We present a general workflow for survival analysis, with
emphasis on dealing with high-dimensional omics data as inputs when identifying
survival-associated omics features and validating survival models. In
particular, we focus on commonly used Cox-type penalized regressions and
hierarchical Bayesian models for feature selection in survival analysis, but
the framework and pipeline are applicable more generally. In cases where
multi-omics data are available for survival modelling, an extra caution is
needed to account for the underlying structure both within and between the
omics data sets and features. A step-by-step R tutorial using The Cancer Genome
Atlas survival and omics data for the execution and evaluation of survival
models has been made available at
\url{https://ocbe-uio.github.io/survomics/survomics.html}.


  Protein language models have excelled in a variety of tasks, ranging from
structure prediction to protein engineering. However, proteins are highly
diverse in functions and structures, and current state-of-the-art models
including the latest version of AlphaFold rely on Multiple Sequence Alignments
(MSA) to feed in the evolutionary knowledge. Despite their success, heavy
computational overheads, as well as the de novo and orphan proteins remain
great challenges in protein representation learning. In this work, we show that
MSAaugmented models inherently belong to retrievalaugmented methods. Motivated
by this finding, we introduce Retrieved Sequence Augmentation(RSA) for protein
representation learning without additional alignment or pre-processing. RSA
links query protein sequences to a set of sequences with similar structures or
properties in the database and combines these sequences for downstream
prediction. We show that protein language models benefit from the retrieval
enhancement on both structure prediction and property prediction tasks, with a
5% improvement on MSA Transformer on average while being 373 times faster. In
addition, we show that our model can transfer to new protein domains better and
outperforms MSA Transformer on de novo protein prediction. Our study fills a
much-encountered gap in protein prediction and brings us a step closer to
demystifying the domain knowledge needed to understand protein sequences. Code
is available on https://github.com/HKUNLP/RSA.


  Escherichia coli is one of many bacterial inhabitants found in human
intestines and any adaptation as a result of mutations may affect its host. A
commonly used technique employed to study these mutations is Restriction
Fragment Length Polymorphism (RFLP) and is proceeded with a suitable distance
coefficient to quantify genetic differences between 2 samples. Dice is
considered a suitable distance coefficient in RFLP analyses, while others were
left unstudied in its suitability for use. Hence, this study aims to identify
substitutes for Dice. Experimental data was obtained by subculturing E. coli
for 72 passages in 8 different adaptation media and RFLP profiles analyzed
using 20 distance coefficients. Our results suggest that Dennis, Fossum,
Matching and Russel and Rao to work as well or better than Dice. Dennis,
Matching and Fossum coefficients had highest discriminatory abilities but are
limited by the lack of upper or lower boundaries. Russel and Rao coefficient is
highly correlated with Dice coefficient (r2 = 0.998), with both higher and
lower boundaries, suggesting that Russel and Rao coefficient can be used to
substitute Dice coefficient in studying genetic distances in E. coli.


  This work evaluates the severity of COVID-19 Omicron variant in terms of the
case-fatality-rates (CFR) with respective uncertainty intervals via a
simultaneous fitting of confirmed cases and deaths in the USA, India, Brazil,
France, Germany and United Kingdom. The CFRs were calculated combining Monte
Carlo simulations and analytical methods based on Gompertz functions under the
framework of the Least Square Method and assuming that the deaths can be
described by a convolution of the confirmed cases with a common gamma function
to describe the case to death period. Linear backgrounds both for cases and
deaths were included in the fitting to account for the contributions from other
strains within the Omicron peaks. The fitting included 125 and 113
epidemiological weeks, for cases and deaths, respectively, and 64 parameters,
resulting in a chi^2 of 176.5 for 174 degrees of freedom (p = 0.434). The CFRs
with 95% confidence intervals for USA, India, Brazil, France, Germany and
United Kingdom were 0.295 (0.154-0.436)%, 0.232 (0.134-0.331)%, 0.49
(0.27-0.71)%, 0.056 (0.028-0.084)%, 0.129 (0.074-0.184)% and 0.168
(0.107-0.229)%, respectively. The case to death period was satisfactorily
described by a common gamma function with mean of 15.71 +/- 0.55 days and
coefficient of variation of 0.354 +/- 0.070. The proposed calculations provided
accurate and reliable information about the respective CFRs and the case to
death period, the latter being consistent with previous estimates for the
symptom onset to death made in the early stages of the pandemic. The CFRs thus
obtained are considerably lower than previous measurements available in the
literature, suggesting that the latter may have been overestimated, as the
probability of deaths from other strains of the virus under the generally
prominent Omicron peak, here accounted for in terms of linear backgrounds, was
not considered.


  The recent development of high-throughput sequencing creates a large
collection of multi-omics data, which enables researchers to better investigate
cancer molecular profiles and cancer taxonomy based on molecular subtypes.
Integrating multi-omics data has been proven to be effective for building more
precise classification models. Current multi-omics integrative models mainly
use early fusion by concatenation or late fusion based on deep neural networks.
Due to the nature of biological systems, graphs are a better representation of
bio-medical data. Although few graph neural network (GNN) based multi-omics
integrative methods have been proposed, they suffer from three common
disadvantages. One is most of them use only one type of connection, either
inter-omics or intra-omic connection; second, they only consider one kind of
GNN layer, either graph convolution network (GCN) or graph attention network
(GAT); and third, most of these methods lack testing on a more complex cancer
classification task. We propose a novel end-to-end multi-omics GNN framework
for accurate and robust cancer subtype classification. The proposed model
utilizes multi-omics data in the form of heterogeneous multi-layer graphs that
combines both inter-omics and intra-omic connections from established
biological knowledge. The proposed model incorporates learned graph features
and global genome features for accurate classification. We test the proposed
model on TCGA Pan-cancer dataset and TCGA breast cancer dataset for molecular
subtype and cancer subtype classification, respectively. The proposed model
outperforms four current state-of-the-art baseline models in multiple
evaluation metrics. The comparative analysis of GAT-based models and GCN-based
models reveals that GAT-based models are preferred for smaller graphs with less
information and GCN-based models are preferred for larger graphs with extra
information.


  Atlantic sea lamprey contains two corticoid receptors (CRs), CR1 and CR2,
that are identical except for a four amino acid insert (Thr-Arg-Gln-Gly) in the
CR1 DNA-binding domain (DBD). Steroids are stronger transcriptional activators
of CR2 than of CR1 suggesting that the insert reduces the transcriptional
response of lamprey CR1 to steroids. The DBD in elephant shark
mineralocorticoid receptor (MR) and glucocorticoid receptor (GR), which are
descended from a CR, lack these four amino acids, suggesting that a CR2 is
their common ancestor. To determine if, similar to lamprey CR1, the presence of
this insert in elephant shark MR and GR decreases transcriptional activation by
corticosteroids, we inserted these four CR1-specific residues into the DBD of
elephant shark MR and GR. Compared to steroid activation of wild-type elephant
shark MR and GR, cortisol, corticosterone, aldosterone, 11-deoxycorticosterone
and 11-deoxycortisol had lower transcriptional activation of these mutant MR
and GR receptors, indicating that the absence of this four-residue segment in
the DBD in wild-type elephant shark MR and GR increases transcriptional
activation by corticosteroids.


  In recent years, Reinforcement Learning (RL) has emerged as a powerful tool
for solving a wide range of problems, including decision-making and genomics.
The exponential growth of raw genomic data over the past two decades has
exceeded the capacity of manual analysis, leading to a growing interest in
automatic data analysis and processing. RL algorithms are capable of learning
from experience with minimal human supervision, making them well-suited for
genomic data analysis and interpretation. One of the key benefits of using RL
is the reduced cost associated with collecting labeled training data, which is
required for supervised learning. While there have been numerous studies
examining the applications of Machine Learning (ML) in genomics, this survey
focuses exclusively on the use of RL in various genomics research fields,
including gene regulatory networks (GRNs), genome assembly, and sequence
alignment. We present a comprehensive technical overview of existing studies on
the application of RL in genomics, highlighting the strengths and limitations
of these approaches. We then discuss potential research directions that are
worthy of future exploration, including the development of more sophisticated
reward functions as RL heavily depends on the accuracy of the reward function,
the integration of RL with other machine learning techniques, and the
application of RL to new and emerging areas in genomics research. Finally, we
present our findings and conclude by summarizing the current state of the field
and the future outlook for RL in genomics.


  Current epidemics in the biological and social domains are challenging the
standard assumptions of mathematical contagion models. Chief among them are the
complex patterns of transmission caused by heterogeneous group sizes and
infection risk varying by orders of magnitude in different settings, like
indoor versus outdoor gatherings in the COVID-19 pandemic or different
moderation practices in social media communities. However, quantifying these
heterogeneous levels of risk is difficult and most models typically ignore
them. Here, we include these novel features in an epidemic model on weighted
hypergraphs to capture group-specific transmission rates. We study analytically
the consequences of ignoring the heterogeneous transmissibility and find an
induced superlinear infection rate during the emergence of a new outbreak, even
though the underlying mechanism is a simple, linear contagion. The dynamics
produced at the individual and group levels are therefore more similar to
complex, nonlinear contagions, thus blurring the line between simple and
complex contagions in realistic settings. We support this claim by introducing
a Bayesian inference framework to quantify the nonlinearity of contagion
processes. We show that simple contagions on real weighted hypergraphs are
systematically biased toward the superlinear regime if the heterogeneity of the
weights is ignored, greatly increasing the risk of erroneous classification as
complex contagions. Our results provide an important cautionary tale for the
challenging task of inferring transmission mechanisms from incidence data. Yet,
it also paves the way for effective models that capture complex features of
epidemics through nonlinear infection rates.


  We study the spatial rock-paper-scissors model, where resource competitors'
cyclic dominance impacts organisms' energy levels. Our model assumes that
failed selection interactions can lead to energy loss, reducing the chances of
success in the spatial game and hastening decline. To prevent death by energy
insufficiency, organisms of one out of the species strategically perform a
trade-off between reproduction and mobility. When prioritising exploring more
extensive areas, organisms aim to maximise the chances of acquiring resources
to regain high energy levels. Through simulation, we examine the effect of
survival behaviour on species segregation and spatial patterns. Our outcomes
show that the trade-off between offspring generation and accelerated movement
effectively protects individuals from death due to lack of energy. Moreover,
the risk of being eliminated by an enemy in the cyclic game reduces due to the
behavioural strategy. Considering a three-state model, we quantify how the
trade-off parameter controls the organisms' energy recovery. Computing the
median organisms' survival time, we find that although individuals performing
the trade-off strategy may live longer, the organisms of other species are
negatively affected by a life expectancy reduction. Our research may elucidate
the role of adaptive survival strategies in species persistence and provide
valuable insights for ecologists.


  Starting from a stochastic individual-based description of an SIS epidemic
spreading on a random network, we study the dynamics when the size n of the
network tends to infinity. We recover in the limit an infinite-dimensional
integro-differential equation studied by Delmas, Dronnier and Zitt (2022) for
an SIS epidemic propagating on a graphon. Our work covers the case of dense and
sparse graphs, provided that the number of edges grows faster than n, but not
the case of very sparse graphs with O(n) edges. In order to establish our limit
theorem, we have to deal with both the convergence of the random graphs to the
graphon and the convergence of the stochastic process spreading on top of these
random structures: in particular, we propose a coupling between the process of
interest and an epidemic that spreads on the complete graph but with a modified
infection rate.
  Keywords: random graph, mathematical models of epidemics, measure-valued
process, large network limit, limit theorem, graphon.


  Protein-protein interaction (PPI) networks provide valuable insights into the
function of biological systems, and aligning multiple PPI networks can reveal
important functional relationships between different species. However,
assessing the quality of multiple network alignments is a challenging problem.
In this paper, we propose two new measures, the Squared GO Score (SGS) and the
Exposed G Score, to evaluate the quality of multiple network alignments while
using functional information from Gene Ontology (GO) terms. We also introduce a
$p$-value measure, the Statistical Exposed G Score, to compute the exact
significance of a multiple network alignment based on its revealed GO terms. We
also show that our measures are highly correlated with the recovered Ortholog
count, providing further evidence for their effectiveness. Our work contributes
to the development of more reliable and accurate measures for evaluating
multiple network alignments and has potential applications in predicting gene
function and identifying evolutionary relationships between different species
using multiple network alignment.


  After the recent ground-breaking advances in protein structure prediction,
one of the remaining challenges in protein machine learning is to reliably
predict distributions of structural states. Parametric models of fluctuations
are difficult to fit due to complex covariance structures between degrees of
freedom in the protein chain, often causing models to either violate local or
global structural constraints. In this paper, we present a new strategy for
modelling protein densities in internal coordinates, which uses constraints in
3D space to induce covariance structure between the internal degrees of
freedom. We illustrate the potential of the procedure by constructing a
variational autoencoder with full covariance output induced by the constraints
implied by the conditional mean in 3D, and demonstrate that our approach makes
it possible to scale density models of internal coordinates to full protein
backbones in two settings: 1) a unimodal setting for proteins exhibiting small
fluctuations and limited amounts of available data, and 2) a multimodal setting
for larger conformational changes in a high data regime.


  Antibiotics resistance has caused much complication in the treatment of
diseases, where the pathogen is no longer susceptible to specific antibiotics
and the use of such antibiotics are no longer effective for treatment. A recent
study that utilizes digital organisms suggests that complete elimination of
specific antibiotic resistance is unlikely after the disuse of antibiotics,
assuming that there are no fitness costs for maintaining resistance once
resistance are established. Fitness cost are referred to as reaction to change
in environment, where organism improves its' abilities in one area at the
expense of the other. Our goal in this study is to use digital organisms to
examine the rate of gain and loss of resistance where fitness costs have
incurred in maintaining resistance. Our results showed that GC-content based
fitness cost during de-selection by removal of antibiotic-induced selective
pressure portrayed similar trends in resistance compared to that of no fitness
cost, at all stages of initial selection, repeated de-selection and
re-introduction of selective pressure. Paired t-test suggested that prolonged
stabilization of resistance after initial loss is not statistically significant
for its difference to that of no fitness cost. This suggests that complete
elimination of specific antibiotics resistance is unlikely after the disuse of
antibiotics despite presence of fitness cost in maintaining antibiotic
resistance during the disuse of antibiotics, once a resistant pool of
micro-organism has been established.


  Owing to the analogies between the problem of wealth redistribution with
taxation in a multi-agent society, we introduce and discuss a kinetic model
describing the statistical distributions in time of the sizes of groups of
biological systems with prey-predator dynamic. While the evolution of the mean
values is shown to be driven by a classical Lotka-Volterra system of
differential equations, it is shown that the time evolution of the probability
distributions of the size of groups of the two interacting species is heavily
dependent both on a kinetic redistribution operator and the degree of
randomness present in the system. Numerical experiments are given to clarify
the time-behavior of the distributions of groups of the species.


  Researchers developed chromosome capture methods such as Hi-C to better
understand DNA's 3D folding in nuclei. The Hi-C method captures contact
frequencies between DNA segment pairs across the genome. When analyzing Hi-C
data sets, it is common to group these pairs using standard bioinformatics
methods (e.g., PCA). Other approaches handle Hi-C data as weighted networks,
where connected node represent DNA segments in 3D proximity. In this
representation, one can leverage community detection techniques developed in
complex network theory to group nodes into mesoscale communities containing
similar connection patterns. While there are several successful attempts to
analyze Hi-C data in this way, it is common to report and study the most
typical community structure. But in reality, there are often several valid
candidates. Therefore, depending on algorithm design, different community
detection methods focusing on slightly different connectivity features may have
differing views on the ideal node groupings. In fact, even the same community
detection method may yield different results if using a stochastic algorithm.
This ambiguity is fundamental to community detection and shared by most complex
networks whenever interactions span all scales in the network. This is known as
community inconsistency. This paper explores this inconsistency of 3D
communities in Hi-C data for all human chromosomes. We base our analysis on two
inconsistency metrics, one local and one global, and quantify the network
scales where the community separation is most variable. For example, we find
that TADs are less reliable than A/B compartments and that nodes with highly
variable node-community memberships are associated with open chromatin.
Overall, our study provides a helpful framework for data-driven researchers and
increases awareness of some inherent challenges when clustering Hi-C data into
3D communities.


  Collective action -- behavior that arises from the combined actions of
multiple individuals -- is observed across living beings. The question of how
and why collective action evolves has profound implications for behavioral
ecology, multicellularity, and human society. Collective action is challenging
to model mathematically, due to nonlinear fitness effects and the consequences
of spatial, group, and/or family relationships. We derive a simple condition
for collective action to be favored by natural selection. A collective's effect
on the fitness of each individual is weighted by the relatedness between them,
using a new measure of collective relatedness. If selection is weak, this
condition can be evaluated using coalescent theory. More generally, our result
applies to any synergistic social behavior, in spatial, group, and/or
family-structured populations. We use this result to obtain conditions for the
evolution of collective help among diploid siblings, subcommunities of a
network, and hyperedges of a hypergraph. We also obtain a condition for which
of two strategies is favored in a game between siblings, cousins, or other
relatives. Our work provides a rigorous basis for extending the notion of
``actor", in the study of social behavior, from individuals to collectives.


  When biological communities use signaling structures for complex
coordination, 'free-riders' emerge. The free-riding agents do not contribute to
the community resources (signals), but exploit them. Most models of such
'selfish' behavior consider free-riding as evolving through mutation and
selection. Over generations, the mutation -- which is considered to create a
stable trait -- spreads through the population. This can lead to a version of
the 'Tragedy of the Commons', where the community's coordination resource gets
fully depleted or deteriorated. In contrast to this evolutionary view, we
present a reinforcement learning model, which shows that both signaling-based
coordination and free-riding behavior can emerge within a generation, through
learning based on energy minimisation. Further, we show that there can be two
types of free-riding, and both of these are not stable traits, but dynamic
'coagulations' of agent-environment interactions. Our model thus shows how
different kinds of selfish behavior can emerge through self-organization, and
suggests that the idea of selfishness as a stable trait presumes a model based
on mutations. We conclude with a discussion of some social and policy
implications of our model.


  In this thesis we use dynamical systems to analize the consequences of
competition in the arrival of a foreign species in an ecosystem.


  The recent development of multimodal single-cell technology has made the
possibility of acquiring multiple omics data from individual cells, thereby
enabling a deeper understanding of cellular states and dynamics. Nevertheless,
the proliferation of multimodal single-cell data also introduces tremendous
challenges in modeling the complex interactions among different modalities. The
recently advanced methods focus on constructing static interaction graphs and
applying graph neural networks (GNNs) to learn from multimodal data. However,
such static graphs can be suboptimal as they do not take advantage of the
downstream task information; meanwhile GNNs also have some inherent limitations
when deeply stacking GNN layers. To tackle these issues, in this work, we
investigate how to leverage transformers for multimodal single-cell data in an
end-to-end manner while exploiting downstream task information. In particular,
we propose a scMoFormer framework which can readily incorporate external domain
knowledge and model the interactions within each modality and cross modalities.
Extensive experiments demonstrate that scMoFormer achieves superior performance
on various benchmark datasets. Remarkably, scMoFormer won a Kaggle silver medal
with the rank of 24/1221 (Top 2%) without ensemble in a NeurIPS 2022
competition. Our implementation is publicly available at Github.


  Opioid use disorder (OUD) is a chronic and relapsing condition that involves
the continued and compulsive use of opioids despite harmful consequences. The
development of medications with improved efficacy and safety profiles for OUD
treatment is urgently needed. Drug repurposing is a promising option for drug
discovery due to its reduced cost and expedited approval procedures.
Computational approaches based on machine learning enable the rapid screening
of DrugBank compounds, identifying those with the potential to be repurposed
for OUD treatment. We collected inhibitor data for four major opioid receptors
and used advanced machine learning predictors of binding affinity that fuse the
gradient boosting decision tree algorithm with two natural language processing
(NLP)-based molecular fingerprints and one traditional 2D fingerprint. Using
these predictors, we systematically analyzed the binding affinities of DrugBank
compounds on four opioid receptors. Based on our machine learning predictions,
we were able to discriminate DrugBank compounds with various binding affinity
thresholds and selectivities for different receptors. The prediction results
were further analyzed for ADMET (absorption, distribution, metabolism,
excretion, and toxicity), which provided guidance on repurposing DrugBank
compounds for the inhibition of selected opioid receptors. The pharmacological
effects of these compounds for OUD treatment need to be tested in further
experimental studies and clinical trials. Our machine learning studies provide
a valuable platform for drug discovery in the context of OUD treatment.


  With the development of computer-assisted techniques, research communities
including biochemistry and deep learning have been devoted into the drug
discovery field for over a decade. Various applications of deep learning have
drawn great attention in drug discovery, such as molecule generation, molecular
property prediction, retrosynthesis prediction, and reaction prediction. While
most existing surveys only focus on one of the applications, limiting the view
of researchers in the community. In this paper, we present a comprehensive
review on the aforementioned four aspects, and discuss the relationships among
different applications. The latest literature and classical benchmarks are
presented for better understanding the development of variety of approaches.
  We commence by summarizing the molecule representation format in these works,
followed by an introduction of recent proposed approaches for each of the four
tasks. Furthermore, we review a variety of commonly used datasets and
evaluation metrics and compare the performance of deep learning-based models.
Finally, we conclude by identifying remaining challenges and discussing the
future trend for deep learning methods in drug discovery.


  Conjugation is a central characteristic of plasmid biology and an important
mechanism of horizontal gene transfer in bacteria. However, there is little
consensus on how to accurately estimate and report plasmid conjugation rates,
in part due to the wide range of available methods. Given the similarity
between approaches, we propose general reporting guidelines for plasmid
conjugation experiments. These constitute best practices based on recent
literature about plasmid conjugation and methods to measure conjugation rates.
In addition to the general guidelines, we discuss common theoretical
assumptions underlying existing methods to estimate conjugation rates and
provide recommendations on how to avoid violating these assumptions. We hope
this will aid the implementation and evaluation of conjugation rate
measurements, and initiate a broader discussion regarding the practice of
quantifying plasmid conjugation rates.


  There has been renewed interest in understanding the mathematical structure
of ecological population models that lead to overcompensation, the process by
which a population recovers to a higher level after suffering an increase in
predation or harvesting. Here, we construct an age-structured single-species
population model that includes a Lotka-Volterra-type cannibalism interaction.
Depending on the age-dependent structure of the interaction, our model can
exhibit transient or steady-state overcompensation--as well as oscillations of
the total population--phenomena that have been observed in ecological systems.
Analytic and numerical analysis of our model reveals sufficient conditions for
overcompensation and oscillations. We also show how our structured population
PDE model can be reduced to coupled ODE models representing piecewise constant
parameter domains, providing additional mathematical insight into the emergence
of overcompensation.


  Understanding which factors informed pandemic response can help create a more
nuanced perspective on how each state of the United States handled the crisis.
To this end, we create various networks linking states together based on their
similarity in mitigation policies, politics, geographic proximity and COVID-19
case data. We use these networks to analyze the correlation between pandemic
policies and politics, location and case-load from January 2020 through March
2022. We show that the best predictors of a state's response is an aggregate
political affiliation rather than solely governor affiliation as others have
shown. Further, we illustrate that political affiliation is heavily correlated
with policy intensity from June 2020 through June 2021, but has little impact
on policy after June 2021. In contrast, geographic proximity and daily
incidence are not consistently correlated with state's having similar
mitigation policies.


  Genome-wide association study (GWAS) tests single nucleotide polymorphism
(SNP) markers across the genome to localize the underlying causal variant of a
trait. Because causal variants are seldom observed directly, a surrogate model
based on genotyped markers are widely considered. Although many methods
estimating the parameters of the surrogate model have been proposed, the
connection between the surrogate model and the true causal model is yet
investigated. In this work, we establish the connection between the surrogate
model and the true causal model. The connection shows that population structure
is accounted in GWAS by modelling the variant of interest and not the trait.
Such observation explains how environmental confounding can be partially
corrected using genetic covariates and why the previously claimed connection
between PC correction and linear mixed models is incorrect.


  Spatially annotated single-cell datasets provide unprecedented opportunities
to dissect cell-cell communication in development and disease. Heterotypic
signaling includes interactions between different cell types and is well
established in tissue development and spatial organization. Epithelial
organization requires several different programs that are tightly regulated.
Planar cell polarity is the organization of epithelial cells along the planar
axis orthogonal to the apical-basal axis. In this study, we investigate planar
cell polarity factors and explore the implications of developmental regulators
as malignant drivers. Utilizing cancer systems biology analysis, we derive gene
expression network for WNT-ligands (WNT) and their cognate frizzled (FZD)
receptors in skin cutaneous melanoma. The profiles supported by unsupervised
clustering of multiple-sequence alignments identify ligand-independent
signaling and implications for metastatic progression based on the underpinning
developmental spatial program. Omics studies and spatial biology connect
developmental programs with oncological events and explain key spatial features
of metastatic aggressiveness. Dysregulation of prominent planar cell polarity
factors such specific representative of the WNT and FZD families in malignant
melanoma recapitulates the development program of normal melanocytes but in an
uncontrolled and disorganized fashion.


  Cooperation is a crucial aspect of social life, yet understanding the nature
of cooperation and how it can be promoted is an ongoing challenge. One
mechanism for cooperation is indirect reciprocity. According to this mechanism,
individuals cooperate to maintain a good reputation. This idea is embodied in a
set of social norms called the ``leading eight''. When all information is
publicly available, these norms have two major properties. Populations that
employ these norms are fully cooperative, and they are stable against invasion
by alternative norms. In this paper, we extend the framework of the leading
eight in two directions. First, we allow social norms to be stochastic. Such
norms allow individuals to evaluate others with certain probabilities. Second,
we consider norms in which also the reputations of passive recipients can be
updated. Using this framework, we characterize all evolutionarily stable norms
that lead to full cooperation in the public information regime. When only the
donor's reputation is updated, and all updates are deterministic, we recover
the conventional model. In that case, we find two classes of stable norms: the
leading eight and the `secondary sixteen'. Stochasticity can further help to
stabilize cooperation when the benefit of cooperation is comparably small.
Moreover, updating the recipients' reputations can help populations to recover
more quickly from errors. Overall, our study highlights a remarkable trade-off
between the evolutionary stability of a norm and its robustness with respect to
errors. Norms that correct errors quickly require higher benefits of
cooperation to be stable.


  Homochirality is a signature of life on Earth yet its origins remain an
unsolved puzzle. Achieving homochirality is essential for a high-yielding
prebiotic network capable of producing functional polymers like ribonucleic
acid (RNA) and peptides. However, a prebiotically plausible and robust
mechanism to reach homochirality has not been shown to this date. The
chiral-induced spin selectivity (CISS) effect has established a strong coupling
between electron spin and molecular chirality and this coupling paves the way
for breaking the chiral molecular symmetry by spin-selective processes.
Magnetic surfaces can act as chiral agents due to the CISS effect and they can
be templates for the enantioselective crystallization of chiral molecules. Here
we studied the spin-selective crystallization of racemic ribo aminooxazoline
(RAO), an RNA precursor, on magnetite ($Fe_3O_4$) surfaces, achieving an
unprecedented enantiomeric excess of about 60$\%$. Following the initial
enrichment, we then obtained homochiral crystals of RAO after a subsequent
crystallization. Our work combines two necessary features for reaching
homochirality: chiral symmetry-breaking induced by the magnetic surface and
self-amplification by conglomerate crystallization of RAO. Our results
demonstrate a prebiotically plausible way of achieving systems level
homochirality from completely racemic starting materials.


  Coarse-graining (CG) accelerates molecular simulations of protein dynamics by
simulating sets of atoms as singular beads. Backmapping is the opposite
operation of bringing lost atomistic details back from the CG representation.
While machine learning (ML) has produced accurate and efficient CG simulations
of proteins, fast and reliable backmapping remains a challenge. Rule-based
methods produce poor all-atom geometries, needing computationally costly
refinement through additional simulations. Recently proposed ML approaches
outperform traditional baselines but are not transferable between proteins and
sometimes generate unphysical atom placements with steric clashes and
implausible torsion angles. This work addresses both issues to build a fast,
transferable, and reliable generative backmapping tool for CG protein
representations. We achieve generalization and reliability through a combined
set of innovations: representation based on internal coordinates; an
equivariant encoder/prior; a custom loss function that helps ensure local
structure, global structure, and physical constraints; and expert curation of
high-quality out-of-equilibrium protein data for training. Our results pave the
way for out-of-the-box backmapping of coarse-grained simulations for arbitrary
proteins.


  Environmental heterogeneity can drive genetic heterogeneity in expanding
populations; mutant strains may emerge that trade overall growth rate for an
improved ability to survive in patches that are hostile to the wild type. This
evolutionary dynamic is of practical importance when seeking to prevent the
emergence of damaging traits. We show that a sub-critical slow-spreading mutant
can attain dominance even when the density of patches is below their
percolation threshold and predict this transition using geometrical arguments.
This work demonstrates a phenomenon of ''assisted percolation'', where one
sub-critical process assists another to achieve super-criticality.


  We classify all additive invariants of open Petri nets: these are
$\mathbb{N}$-valued invariants which are additive with respect to sequential
and parallel composition of open Petri nets. In particular, we prove two
classification theorems: one for open Petri nets and one for monically open
Petri nets (i.e. open Petri nets whose interfaces are specified by monic maps).
Our results can be summarized as follows. The additive invariants of open Petri
nets are completely determined by their values on a particular class of
single-transition Petri nets. However, for monically open Petri nets, the
additive invariants are determined by their values on transitionless Petri nets
and all single-transition Petri nets. Our results confirm a conjecture of John
Baez (stated during the AMS' 2022 Mathematical Research Communities workshop).


  We derive an exact upper bound on the epidemic overshoot for the
Kermack-McKendrick SIR model. This maximal overshoot value of 0.2984... occurs
at $R_0^*$ = 2.151... . In considering the utility of the notion of overshoot,
a rudimentary analysis of data from the first wave of the COVID-19 pandemic in
Manaus, Brazil highlights the public health hazard posed by overshoot for
epidemics with $R_0$ near 2. Using the general analysis framework presented
within, we then consider more complex SIR models that incorporate vaccination.


  We expand the calculation of the upper bound on epidemic overshoot in SIR
models to account for nonlinear incidence. We lay out the general procedure and
restrictions to perform the calculation analytically. We show why the infected
fraction must enter linearly into the incidence term, whereas the susceptible
fraction can enter more generically, including non-monotonically. We
demonstrate the procedure by working through several examples. We also
numerically study what happens to the upper bound on overshoot when nonlinear
incidence manifests in the form of dynamics over a contact network.


  This paper conducts a numerical study of a geometrical structure called
$\epsilon$-school for predator-avoidance fish schools, based on our previous
mathematical model. Our results show that during a predator attack, the number
of $\epsilon$-school increases from one to a certain value. After the attack,
the number of $\epsilon$-school decreases in the first two predator-avoidance
patterns, but continues to increase in the third pattern. A constant value for
the number of the $\epsilon$-school is observed in the last pattern. These
suggests that when the predator is approaching, each individual in the school
focuses more on avoiding the predator, rather than on interacting with its
schoolmates. Such a trait is in agreement with real-life behavior in the
natural ecosystem.


  Autocatalytic sets are self-catalyzing and self-sustaining chemical reaction
networks that are believed to have played an important role in the origin of
life. They have been studied extensively both theoretically as well as
experimentally. This short note provides (1) a complete and formal definition
of autocatalytic sets (or RAF sets), and (2) an efficient algorithm to detect
such sets in arbitrary reaction networks. Although both have been presented in
various forms in earlier publications, this note serves as a concise and
convenient reference.


  Similarity search is one of the most fundamental computations that are
regularly performed on ever-increasing protein datasets. Scalability is of
paramount importance for uncovering novel phenomena that occur at very large
scales. We unleash the power of over 20,000 GPUs on the Summit system to
perform all-vs-all protein similarity search on one of the largest publicly
available datasets with 405 million proteins, in less than 3.5 hours, cutting
the time-to-solution for many use cases from weeks. The variability of protein
sequence lengths, as well as the sparsity of the space of pairwise comparisons,
make this a challenging problem in distributed memory. Due to the need to
construct and maintain a data structure holding indices to all other sequences,
this application has a huge memory footprint that makes it hard to scale the
problem sizes. We overcome this memory limitation by innovative matrix-based
blocking techniques, without introducing additional load imbalance.


  The Kauffman model of genetic computation highlights the importance of
criticality at the border of order and chaos. But our understanding of its
behavior is incomplete, and much of what we do know relies on heuristic
arguments. To better understand the model and obtain more rigorous insights, we
show that there are fundamental links between the critical Kauffman model and
aspects of number theory. Using these connections, we prove that the number of
attractors and the mean attractor length grow faster than previously believed.
Our work suggests that techniques from number theory, which are less familiar
to the physics community, may be the right tools for fully cracking the
Kauffman model.


  In this paper, we study the spread of a classical SIR process on a two-layer
random network, where the first layer represents the households, while the
second layer models the contacts outside the households by a random scale-free
graph. We build a three-parameter graph, called polynomial model, where the new
vertices are connected to the existing ones either uniformly, or
preferentially, or by forming random triangles. We examine the effect of the
graph's properties on the goodness of the estimation of the infection rate
$\tau$. In the classical approach, to estimate $\tau$ one needs to approximate
the number of SI edges between households, since the graph itself is supposed
to be unobservable. Our simulation study reveals that the estimation is poorer
at the beginning of the epidemic, for larger preferential attachment parameter
of the graph, and for larger $\tau$. We can improve the estimate either by (i)
adjusting it by a graph-dependent factor at a fixed time in the earlier phase
of the epidemic, or (ii) by applying a shift in the estimation of the SI edges.
Furthermore, we establish our method to be robust to changes in average
clustering of the graph model. We also extend a graph neural network (GNN)
approach for estimating contagion dynamics for our two-layered graphs. We find
that dense networks offer better training datasets. Moreover, GNN perfomance is
measured better using the $l_2$ loss function rather than cross-entropy.


  Cells of diatoms and related algae with complex plastids of red algal origin
are highly compartmentalized. These plastids are surrounded by four envelope
membranes, which also define the periplastidic compartment (PPC), the space
between the second and third membranes. The PPC corresponds to the cytosol of
the eukaryotic alga that was the ancestor of the complex plastid. Metabolic
reactions as well as cell biological processes take place in this compartment;
however, its exact function remains elusive. Automated predictions of protein
locations proved useful for genome wide explorations of metabolism in the case
of plastid proteins, but until now, no automated method for the prediction of
PPC proteins was available. Here, we present an updated version of the plastid
protein predictor ASAFind, which includes optional prediction of PPC proteins.
The new ASAFind version also accepts the output of the most recent versions of
SignalP (5.0) and TargetP (2.0) input data. Furthermore, we release a Python
script to calculate custom scoring matrices for adjustment of the ASAFind
method to other groups of algae, and included the option to run the predictions
with custom scoring matrices in a simplified score cut-off mode.


  Plastid targeted proteins of diatoms and related algae can be predicted with
high sensitivity and specificity using the ASAFind method published in 2015.
ASAFind predictions rely on SignalP predictions of endoplasmic reticulum (ER)
targeting signal peptides. Recently (in 2019), a new version of SignalP was
released, SignalP 5.0. We tested the ability of SignalP 5.0 to recognize signal
peptides of nucleus-encoded, plastid-targeted diatom pre-proteins, and to
identify the signal peptide cleavage site. The results were compared to manual
predictions of the characteristic cleavage site motif, and to previous versions
of SignalP. SignalP 5.0 is less sensitive than the previous versions of SignalP
in this specific task, and also in the detection of signal peptides of
non-plastid proteins in diatoms. However, in combination with ASAFind, the
resulting prediction performance for plastid proteins is high. In addition, we
tested the multi-location prediction tool TargetP for its suitability to
provide signal peptide information to ASAFind. The newest version, TargetP 2.0,
had the highest prediction performances for diatom signal peptides and
mitochondrial transit peptides compared to other versions of SignalP and
TargetP, thus it provides a good basis for ASAFind predictions.


  In the critical beta-splitting model of a random $n$-leaf rooted tree, clades
are recursively split into sub-clades, and a clade of $m$ leaves is split into
sub-clades containing $i$ and $m-i$ leaves with probabilities $\propto
1/(i(m-i))$. This article provides an extensive overview of structure theory
and explicit quantitative aspects. There is a canonical embedding into a
continuous-time model, that is a random tree CTCS(n) on $n$ leaves with
real-valued edge lengths, and this model turns out more convenient to study. We
show that the family $(\mbox{CTCS}(n), n \ge 2)$ is consistent under a ``delete
random leaf and prune" operation. That leads to an explicit inductive
construction of $(\mbox{CTCS}(n), n \ge 2)$ as $n$ increases. An accompanying
technical article arXiv:2302.05066 studies in detail (with perhaps surprising
precision) many distributions relating to the heights of leaves, via analytic
methods. We give alternative probabilistic proofs for some such results, which
provides an opportunity for a ``compare and contrast" discussion of the two
methodologies. We prove existence of the limit {\em fringe distribution}
relative to a random leaf, whose graphical representation is essentially the
format of the cladogram representation of biological phylogenies. We describe
informally the scaling limit process, as a process of splitting the continuous
interval $(0,1)$.
  These topics are somewhat analogous to those topics which have been well
studied in the context of the Brownian continuum random tree. Many open
problems remain.


  If the origin of life is rare and sensitive to the local conditions at the
site of its emergence, then, using the principle of mediocrity within a
multiverse framework, we may expect to find ourselves in a universe that is
better than usual at creating these necessary conditions. We use this reasoning
to investigate several origin of life scenarios to determine whether they are
compatible with the multiverse, including the prebiotic soup scenario,
hydrothermal vents, delivery of prebiotic material from impacts, and
panspermia. We find that most of these scenarios induce a preference toward
weaker-gravity universes, and that panspermia and scenarios involving solar
radiation or large impacts as a disequilibrium source are disfavored.
Additionally, we show that several hypothesized habitability criteria which are
disfavored when the origin of life is not taken into account become compatible
with the multiverse, and that the emergence of life and emergence of
intelligence cannot both be sensitive to disequilibrium production conditions.


  Folding of ribozymes into well-defined tertiary structures usually requires
divalent cations. How Mg$^{2+}$ ions direct the folding kinetics has been a
long-standing unsolved problem because experiments cannot detect the positions
and dynamics of ions. To address this problem, we used molecular simulations to
dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path
each molecule takes to reach the folded state. We quantitatively establish that
Mg$^{2+}$ binding to specific sites, coupled with counter-ion release of
monovalent cations, stimulate the formation of secondary and tertiary
structures, leading to diverse pathways that include direct rapid folding and
trapping in misfolded structures. In some molecules, key tertiary structural
elements form when Mg$^{2+}$ ions bind to specific RNA sites at the earliest
stages of the folding, leading to specific collapse and rapid folding. In
others, the formation of non-native base pairs, whose rearrangement is needed
to reach the folded state, is the rate-limiting step. Escape from energetic
traps, driven by thermal fluctuations, occurs readily. In contrast, the
transition to the native state from long-lived topologically trapped
native-like metastable states is extremely slow. Specific collapse and
formation of energetically or topologically frustrated states occur early in
the assembly process.


  We consider a population organised hierarchically with respect to size in
such a way that the growth rate of each individual depends only on the presence
of larger individuals. As a concrete example one might think of a forest, in
which the incidence of light on a tree (and hence how fast it grows) is
affected by shading of taller trees. The model is formulated as a delay
equation, more specifically a scalar renewal equation, for the population birth
rate. After discussing the well-posedness of the model, we analyse how many
stationary birth rates the equation can have in terms of the functional
parameters of the model. In particular we show that, under reasonable and
rather general assumptions, only one stationary birth rate can exist besides
the trivial one (associated to the state in which there are no individuals and
the population birth rate is zero). We give conditions for this non-trivial
stationary birth rate to exist and we analyse its stability using the principle
of linearised stability for delay equations. Finally we relate the results to
an alternative formulation of the model taking the form of a quasilinear
partial differential equation for the population size-density.


  Ecosystems are commonly organized into trophic levels -- organisms that
occupy the same level in a food chain (e.g., plants, herbivores, carnivores). A
fundamental question in theoretical ecology is how the interplay between
trophic structure, diversity, and competition shapes the properties of
ecosystems. To address this problem, we analyze a generalized Consumer Resource
Model with three trophic levels using the zero-temperature cavity method and
numerical simulations. We find that intra-trophic diversity gives rise to
``emergent competition'' between species within a trophic level due to
feedbacks mediated by other trophic levels. This emergent competition gives
rise to a crossover from a regime of top-down control (populations are limited
by predators) to a regime of bottom-up control (populations are limited by
primary producers) and is captured by a simple order parameter related to the
ratio of surviving species in different trophic levels. We show that our
theoretical results agree with empirical observations, suggesting that the
theoretical approach outlined here can be used to understand complex ecosystems
with multiple trophic levels.


  Phylogenetic trees are frequently used to model evolution. Such trees are
typically reconstructed from data like DNA, RNA, or protein alignments using
methods based on criteria like maximum parsimony (amongst others). Maximum
parsimony has been assumed to work well for data with only few state changes.
Recently, some progress has been made to formally prove this assertion. For
instance, it has been shown that each binary phylogenetic tree $T$ with $n \geq
20k$ leaves is uniquely defined by the set $A_k(T)$, which consists of all
characters with parsimony score $k$ on $T$. In the present manuscript, we show
that the statement indeed holds for all $n \geq 4k$, thus drastically lowering
the lower bound for $n$ from $20k$ to $4k$. However, it has been known that for
$n \leq 2k$ and $k \geq 3$, it is not generally true that $A_k(T)$ defines $T$.
We improve this result by showing that the latter statement can be extended
from $n \leq 2k$ to $n \leq 2k+2$. So we drastically reduce the gap of values
of $n$ for which it is unknown if trees $T$ on $n$ taxa are defined by $A_k(T)$
from the previous interval of $[2k+1,20k-1]$ to the interval $[2k+3,4k-1]$.
Moreover, we close this gap completely for the nearest neighbor interchange
(NNI) neighborhood of $T$ in the following sense: We show that as long as
$n\geq 2k+3$, no tree that is one NNI move away from $T$ (and thus very similar
to $T$) shares the same $A_k$-alignment.


  Feedback loops are essential for regulating cell proliferation and
maintaining the delicate balance between cell division and cell death. Thanks
to the exact solution of a few simple models of cell growth it is by now clear
that stochastic fluctuations play a central role in this process and that cell
growth (and in particular the robustness and stability of homeostasis) can be
properly addressed only as a stochastic process. Using epidermal homeostasis as
a prototypical example, we show that it is possible to discriminate among
different feedback strategies which turn out to be characterized by different,
experimentally testable, behaviours. In particular, we focus on the so-called
Dynamical Heterogeneity model, an epidermal homeostasis model that takes into
account two well known cellular features: the plasticity of the cells and their
adaptability to face environmental stimuli. We show that specific choices of
the parameter on which the feedback is applied may decrease the fluctuations of
the homeostatic population level and improve the recovery of the system after
an external perturbation.


  Activity and property prediction models are the central workhorses in drug
discovery and materials sciences, but currently they have to be trained or
fine-tuned for new tasks. Without training or fine-tuning, scientific language
models could be used for such low-data tasks through their announced zero- and
few-shot capabilities. However, their predictive quality at activity prediction
is lacking. In this work, we envision a novel type of activity prediction model
that is able to adapt to new prediction tasks at inference time, via
understanding textual information describing the task. To this end, we propose
a new architecture with separate modules for chemical and natural language
inputs, and a contrastive pre-training objective on data from large biochemical
databases. In extensive experiments, we show that our method CLAMP yields
improved predictive performance on few-shot learning benchmarks and zero-shot
problems in drug discovery. We attribute the advances of our method to the
modularized architecture and to our pre-training objective.


  Rich data and powerful machine learning models allow us to design drugs for a
specific protein target \textit{in silico}. Recently, the inclusion of 3D
structures during targeted drug design shows superior performance to other
target-free models as the atomic interaction in the 3D space is explicitly
modeled. However, current 3D target-aware models either rely on the voxelized
atom densities or the autoregressive sampling process, which are not
equivariant to rotation or easily violate geometric constraints resulting in
unrealistic structures. In this work, we develop a 3D equivariant diffusion
model to solve the above challenges. To achieve target-aware molecule design,
our method learns a joint generative process of both continuous atom
coordinates and categorical atom types with a SE(3)-equivariant network.
Moreover, we show that our model can serve as an unsupervised feature extractor
to estimate the binding affinity under proper parameterization, which provides
an effective way for drug screening. To evaluate our model, we propose a
comprehensive framework to evaluate the quality of sampled molecules from
different dimensions. Empirical studies show our model could generate molecules
with more realistic 3D structures and better affinities towards the protein
targets, and improve binding affinity ranking and prediction without
retraining.


  Markov state models represent a popular means to interpret molecular dynamics
trajectories in terms of memoryless transitions between metastable
conformational states. To provide a mechanistic understanding of the considered
biomolecular process, these states should reflect structurally distinct
conformations and ensure a timescale separation between fast intrastate and
slow interstate dynamics. Adopting the folding of villin headpiece (HP35) as a
well-established model problem, here we discuss the selection of suitable input
coordinates or `features', such as backbone dihedral angles and interresidue
distances. We show that dihedral angles account accurately for the structure of
the native energy basin of HP35, while the unfolded region of the free energy
landscape and the folding process are best described by tertiary contacts of
the protein. To construct a contact-based model, we consider various ways to
define and select contact distances, and introduce a low-pass filtering of the
feature trajectory as well as a correlation-based characterization of states.
Relying on input data that faithfully account for the mechanistic origin of the
studied process, the states of the resulting Markov model are clearly
discriminated by the features, describe consistently the hierarchical structure
of the free energy landscape, and$\unicode{0x2014}$as a
consequence$\unicode{0x2014}$correctly reproduce the slow timescales of the
process.


  A multistationarity region is the part of a reaction network's parameter
space that gives rise to multiple steady states. Mathematically, this region
consists of the positive parameters for which a parametrized family of
polynomial equations admits two or more positive roots. Much recent work has
focused on analyzing multistationarity regions of biologically significant
reaction networks and determining whether such regions are connected. Here we
focus on the multistationarity regions of small networks, those with few
species and few reactions. For two families of such networks -- those with one
species and up to three reactions, and those with two species and up to two
reactions -- we prove that the resulting multistationarity regions are
connected. We also give an example of a network with one species and six
reactions for which the multistationarity region is disconnected. Our proofs
rely on the formula for the discriminant of a trinomial, a classification of
small multistationary networks, and a recent result of Feliu and Telek that
partially generalizes Descartes' rule of signs.


  Along the pathogenesis of Mycobacterium Tuberculosis (MTB), hypoxia-induced
dormancy is a process involving the oxygen-depleted environment encountered
inside the lung granuloma, where bacilli enter a viable, non-replicating state
termed as latency. Affecting nearly two billion people, latent TB can linger in
the host for indefinite periods of time before resuscitating, which
significantly strains the accuracy of treatment options and patient prognosis.
Transcriptional factors thought to mediate this process have only conferred
mild growth defects, signaling that our current understanding of the MTB
genetic architecture is highly insufficient. In light of these inconsistencies,
the objective of this study was to characterize regulatory mechanisms
underlying the transition of MTB into dormancy. The project methodology
involved a three-part approach - constructing an aggregate hypoxia dataset,
inferring a gene regulatory network based on those observations, and leveraging
several downstream network analyses to make sense of it all. Results indicated
dormancy to be functionally associated with cell redox homeostasis, metal ion
cycling, and cell wall metabolism, all of which modulate essential
host-pathogen interactions. Additionally, the crosstalk between individual
regulons (Rv0821c and Rv0144; Rv1152 and Rv2359) was shown to be critical in
facilitating bacterial persistence and allowing MTB to gain control over key
micronutrients within the cell. Defense antioxidants and nutritional immunity
were also identified as future avenues to explore further. In providing some of
the first insights into the methods utilized by MTB to endure in a hypoxic
state, this research suggests a range of strategies that might aid in improved
clinical outcomes of TB treatment.


  The well known phenomenon of phase separation in synthetic polymers and
proteins has become a major topic in biophysics because it has been invoked as
a mechanism of compartment formation in cells, without the need for membranes.
Most of the coacervates (or condensates) are composed of Intrinsically
Disordered Proteins (IDPs) or regions that are structureless, often in
interaction with RNA and DNA. One of the more intriguing IDPs is the
526-residue RNA binding protein, Fused In Sarcoma (FUS), whose monomer
conformations and condensates exhibit unusual behavior that are sensitive to
solution conditions. By focussing principally on the N-terminus low complexity
domain (FUS-LC comprising residues 1-214) and other truncations, we rationalize
the findings of solid state NMR experiments, which show that FUS-LC adopts a
non-polymorphic fibril (core-1) involving residues 39-95, flanked by fuzzy
coats on both the N- and C- terminal ends. An alternate structure (core-2),
whose free energy is comparable to core-1, emerges only in the truncated
construct (residues 110-214). Both core-1 and core-2 fibrils are stabilized by
a Tyrosine ladder as well as hydrophilic interactions. The morphologies (gels,
fibrils, and glass-like behavior) adopted by FUS seem to vary greatly,
depending on the experimental conditions. The effect of phosphorylation is site
specific and affects the stability of the fibril depending on the sites that
are phosphorylated. Many of the peculiarities associated with FUS may also be
shared by other IDPs, such as TDP43 and hnRNPA2. We outline a number of
problems for which there is no clear molecular understanding.


  The rapid growth in genomic pathogen data spurs the need for efficient
inference techniques, such as Hamiltonian Monte Carlo (HMC) in a Bayesian
framework, to estimate parameters of these phylogenetic models where the
dimensions of the parameters increase with the number of sequences $N$. HMC
requires repeated calculation of the gradient of the data log-likelihood with
respect to (wrt) all branch-length-specific (BLS) parameters that traditionally
takes $\mathcal{O}(N^2)$ operations using the standard pruning algorithm. A
recent study proposes an approach to calculate this gradient in
$\mathcal{O}(N)$, enabling researchers to take advantage of gradient-based
samplers such as HMC. The CPU implementation of this approach makes the
calculation of the gradient computationally tractable for nucleotide-based
models but falls short in performance for larger state-space size models, such
as codon models. Here, we describe novel massively parallel algorithms to
calculate the gradient of the log-likelihood wrt all BLS parameters that take
advantage of graphics processing units (GPUs) and result in many fold higher
speedups over previous CPU implementations. We benchmark these GPU algorithms
on three computing systems using three evolutionary inference examples:
carnivores, dengue and yeast, and observe a greater than 128-fold speedup over
the CPU implementation for codon-based models and greater than 8-fold speedup
for nucleotide-based models. As a practical demonstration, we also estimate the
timing of the first introduction of West Nile virus into the continental Unites
States under a codon model with a relaxed molecular clock from 104 full viral
genomes, an inference task previously intractable. We provide an implementation
of our GPU algorithms in BEAGLE v4.0.0, an open source library for statistical
phylogenetics that enables parallel calculations on multi-core CPUs and GPUs.


  Catch-bonds, whereby noncovalent ligand-receptor interactions are
counterintuitively reinforced by tensile forces, play a major role in cell
adhesion under mechanical stress. A basic prerequisite for catch-bond formation
is that force-induced remodeling of ligand binding interface occurs prior to
bond rupture. However, what strategy receptor proteins utilize to meet such
specific kinetic control is still unclear, rendering the mechanistic
understanding of catch-bond an open question. Here we report a bidirectional
allostery mechanism of catch-bond for the hyaluronan (HA) receptor CD44 which
is responsible for rolling adhesion of lymphocytes and circulating tumor cells.
Binding of ligand HA allosterically reduces the threshold force for unlocking
of otherwise stably folded force-sensing element (i.e., forward allostery), so
that much smaller tensile force can trigger the conformational switching of
receptor protein to high binding-strength state via backward allosteric
coupling before bond rupture. The effect of forward allostery was further
supported by performing atomistic molecular dynamics simulations. Such
bidirectional allostery mechanism fulfills the specific kinetic control
required by catch-bond and is likely to be commonly utilized in cell adhesion.
We also revealed a slip-catch-slip triphasic pattern in force response of
CD44-HA bond arising from force-induced repartitioning of parallel dissociation
pathways. The essential thermodynamic and kinetic features of receptor proteins
for shaping the catch-bond were identified.


  Indirect reciprocity is one of the major mechanisms for the evolution of
cooperation in human societies. There are two types of indirect reciprocity:
upstream and downstream. Cooperation in downstream reciprocity follows the
pattern, 'You helped someone, and I will help you'. The direction of
cooperation is reversed in upstream reciprocity, which instead follows the
pattern, 'You helped me, and I will help someone else'. In reality, these two
types of indirect reciprocity often occur in combination. However, upstream and
downstream reciprocity have mostly been studied theoretically in isolation.
Here, we propose a new model that integrates both types. We apply the standard
giving-game framework of indirect reciprocity and analyze the model by means of
evolutionary game theory. We show that the model can result in the stable
coexistence of altruistic reciprocators and free riders in well-mixed
populations. We also found that considering inattention in the assessment rule
can strengthen the stability of this mixed equilibrium, even resulting in a
global attractor. Our results indicate that the cycles of forwarding help and
rewarding help need to be established for creating and maintaining diversity
and inclusion in a society.


  Interactions between commuting individuals can lead to large-scale spreading
of rumors, ideas, or disease, even though the commuters have no net
displacement. The emergent dynamics depend crucially on the commuting
distribution of a population, that is how the probability to travel to a
destination decays with distance from home. Applying this idea to epidemics, we
will demonstrate the qualitatively different infection dynamics emerging from
populations with different commuting distributions. If the commuting
distribution is exponentially localized, we recover a reaction-diffusion system
and observe Fisher waves traveling at a speed proportional to the
characteristic commuting distance. If the commuting distribution has a long
tail, then no finite-velocity waves can form, but we show that, in some
regimes, there is nontrivial spatial dependence that the well-mixed
approximation neglects. We discuss how, in all cases, an initial
dispersal-dominated regime can allow the disease to go undetected for a finite
amount of time before exponential growth takes over. This "offset time" is a
quantity of huge importance for epidemic surveillance and yet largely ignored
in the literature.


  We construct a reliable estimation of evolutionary parameters within the
Wright-Fisher model, which describes changes in allele frequencies due to
selection and genetic drift, from time-series data. Such data exists for
biological populations, for example via artificial evolution experiments, and
for the cultural evolution of behavior, such as linguistic corpora that
document historical usage of different words with similar meanings. Our method
of analysis builds on a Beta-with-Spikes approximation to the distribution of
allele frequencies predicted by the Wright-Fisher model. We introduce a
self-contained scheme for estimating the parameters in the approximation, and
demonstrate its robustness with synthetic data, especially in the
strong-selection and near-extinction regimes where previous approaches fail. We
further apply to allele frequency data for baker's yeast (Saccharomyces
cerevisiae), finding a significant signal of selection in cases where
independent evidence supports such a conclusion. We further demonstrate the
possibility of detecting time-points at which evolutionary parameters change in
the context of a historical spelling reform in the Spanish language.


  This paper introduces a variational formulation of natural selection, paying
special attention to the nature of "things" and the way that different "kinds"
of "things" are individuated from - and influence - each other. We use the
Bayesian mechanics of particular partitions to understand how slow phylogenetic
processes constrain - and are constrained by - fast, phenotypic processes. The
main result is a formulation of adaptive fitness as a path integral of
phenotypic fitness. Paths of least action, at the phenotypic and phylogenetic
scales, can then be read as inference and learning processes, respectively. In
this view, a phenotype actively infers the state of its econiche under a
generative model, whose parameters are learned via natural (bayesian model
selection). The ensuing variational synthesis features some unexpected aspects.
Perhaps the most notable is that it is not possible to describe or model a
population of conspecifics per se. Rather, it is necessary to consider
populations - and nested meta-populations - of different natural kinds that
influence each other. This paper is limited to a description of the
mathematical apparatus and accompanying ideas. Subsequent work will use these
methods for simulations and numerical analyses - and identify points of contact
with related mathematical formulations of evolution.


  Population ecology theory is replete with density dependent processes.
However trait-mediated or behavioral indirect interactions can both reinforce
or oppose density-dependent effects. This paper presents the first two species
competitive ODE and PDE systems where an Allee effect, which is a density
dependent process and the fear effect, which is non-consumptive and behavioral
are both present. The stability of the equilibria is discussed analytically
using the qualitative theory of ordinary differential equations. It is found
that the Allee effect and the fear effect change the extinction dynamics of the
system and the number of positive equilibrium points, but they do not affect
the stability of the positive equilibria. We also observe some special dynamics
that induce bifurcations in the system by varying the Allee or fear parameter.
Interestingly we find that the Allee effect working in conjunction with the
fear effect, can bring about several qualitative changes to the dynamical
behavior of the system with only the fear effect in place, in regimes of small
fear. That is, for small amounts of the fear parameter, it can change a
competitive exclusion type situation to a strong competition type situation. It
can also change a weak competition type situation to a bi-stability type
situation. However for large fear regimes the Allee effect reinforces the
dynamics driven by the fear effect. The analysis of the corresponding spatially
explicit model is also presented. To this end the comparison principle for
parabolic PDE is used. The conclusions of this paper have strong implications
for conservation biology, biological control as well as the preservation of
biodiversity.


  Endophytes are tiny organisms present in living tissues of distinct plants
and have been extensively studied for their endophytic microbial complement.
Roots of Rosmarinus officinalis were subjected to the isolation of endophytic
fungi and screened for antimicrobial activity against Gram-positive
(Staphylococcus aureus and Bacillus subtilis) and Gram-negative (Escherichia
coli, Pseudomonas aeruginosa, Klebsiella pneumoniae) bacteria. Genomic DNA from
active fungal strain of Trichoderma harzianum was isolated, and the internal
transcribed spacer (ITS) region was amplified using ITS4 and ITS5 primers and
sequenced for genetic inference in fungus. The crude extract of T. harzianum
isolate with Ethyl acetate was showed significant antimicrobial activity
against P. aeruginosa, S. aureus, K. pneumonia, B. subtilis and E. coli. The
antimicrobial activity was highest against P. aeruginosa at concentration of 40
microgram/ ml, followed by S. aureus and K. pneumonia at the same
concentration. The lowest antimicrobial activity was against by S. aureus at
concentration of 60 microgram/ ml. The current study is confirmed that the
antimicrobial activity is due to bioactive compounds founded in endophytic
fungi.


  Models with several levels of mixing (households, workplaces), as well as
various corresponding formulations for $R_0$, have been proposed in the
literature. However, little attention has been paid to the impact of the
distribution of the population size within social structures, effect that can
help plan effective interventions. We focus on the influence on the model
outcomes of teleworking strategies, consisting in reshaping the distribution of
workplace sizes. We consider a stochastic \emph{SIR} model with two levels of
mixing, accounting for a uniformly mixing general population, each individual
belonging also to a household and a workplace. The variance of the workplace
size distribution appears to be a good proxy for the impact of this
distribution on key outcomes of the epidemic, such as epidemic size and peak.
In particular, our findings suggest that strategies where the proportion of
individuals teleworking depends sublinearly on the size of the workplace
outperform the strategy with linear dependence. Besides, one drawback of the
model with multiple levels of mixing is its complexity, raising interest in a
reduced model. We propose an unstructured \emph{SIR} ODE-based model,
explicitly exhibiting social structure sizes. This reduced model, sharing the
same growth rate as the initial model, yields a generally satisfying
approximation of the epidemic. These results, robust to various changes in
model structure, are very promising from the perspective of implementing
effective strategies based on social distancing of specific contacts.
Furthermore, they contribute to the effort of building relevant approximations
of individual based models at intermediate scales.


  In this paper we propose an unbiased Monte Carlo maximum likelihood estimator
for discretely observed Wright-Fisher diffusions. Our approach is based on
exact simulation techniques that are of special interest for diffusion
processes defined on a bounded domain, where numerical methods typically fail
to remain within the required boundaries. We start by building unbiased maximum
likelihood estimators for scalar diffusions and later present an extension to
the multidimensional case. Consistency results of our proposed estimator are
also presented and the performance of our method is illustrated through a
numerical example.


  Several SARS-CoV-2 variants that evolved during the COVID-19 pandemic have
appeared to differ in severity, based on analyses of single-country datasets.
With decreased SARS-CoV-2 testing and sequencing, international collaborative
studies will become increasingly important for timely assessment of the
severity of newly emerged variants. The Joint WHO Regional Office for Europe
and ECDC Infection Severity Working Group was formed to produce and pilot a
standardised study protocol to estimate relative variant case-severity in
settings with individual-level SARS-CoV-2 testing and COVID-19 outcome data
during periods when two variants were co-circulating. To assess feasibility,
the study protocol and its associated statistical analysis code was applied by
local investigators in Denmark, England, Luxembourg, Norway, Portugal and
Scotland to assess the case-severity of Omicron BA.1 relative to Delta cases.
After pooling estimates using meta-analysis methods (random effects estimates),
the risk of hospital admission (adjusted hazard ratio [aHR]=0.41, 95% CI
0.31-0.54), ICU admission (aHR=0.12, 95% CI 0.05-0.27), and death (aHR=0.31,
95% CI 0.28-0.35) was lower for Omicron BA.1 compared to Delta cases. The aHRs
varied by age group and vaccination status. In conclusion, this study has
demonstrated the feasibility of conducting variant severity analyses in a
multinational collaborative framework. The results add further evidence for the
reduced severity of the Omicron BA.1 variant.


  A theory of aging based on the principles of the kinetics of chemical
reactions and the rules of natural selection of organisms is proposed. The
theory is based on a hypothesis that the biochemical processes in the organism
can be described in the terms of chemical reaction kinetics. The evolutionary
process of organisms is determined by the goal of continuing life, and the
natural selection forced organisms to develop in an optimized way for survival
and reproduction, after which any further development of the organisms did not
matter for natural selection and therefore was not regulated by it.
Accordingly, the ongoing biochemical processes in the organism after giving
birth were not stabilised and continued casual, which, in the case of complex
systems of biochemical processes, led to imbalances and biochemical processes
that did not contribute to health in the organism, resulting in aging. Based on
this view, it is necessary to look for the key biochemical processes that
regulate the vital activity of the organism. Balancing of these processes in
the period after giving birth can artificially lead to the stabilization of the
kinetics of biochemical reactions, and consequently, the continuation of life
almost indefinitely.


  The June software package is an open-source framework for the detailed
simulation of epidemics based on social interactions in a virtual population
reflecting age, gender, ethnicity, and socio-economic indicators in England. In
this paper, we present a new version of the framework specifically adapted for
Germany, which allows the simulation of the entire German population using
publicly available information on households, schools, universities,
workplaces, and mobility data for Germany. Moreover, JuneGermany incorporates
testing and vaccination strategies within the population as well as the
simultaneous handling of several different virus strains. First validation
tests of the framework have been performed for the state of Rhineland
Palatinate based on data collected between October 2020 and December 2020 and
then extrapolated to March 2021, i.e. the end of the second wave.


  As observed during the global SARS-CoV-2 (COVID-19) pandemic, high-income
countries, such as the United States, may exhibit vaccine nationalism during a
pandemic: stockpiling doses of vaccine for their own citizens and being
reluctant to distribute doses of the vaccine to lower-income countries. While
many cite moral objections to vaccine nationalism, vaccine inequity during a
pandemic could possibly worsen the global effects of the pandemic, including in
the high-income countries themselves, through the evolution of new variants of
the virus. This paper uses the COVID-19 pandemic as a case study to identify
scenarios under which it might be in a high-income nation's own interest to
donate vaccine doses to another country before its own population has been
fully vaccinated. We develop an extended SEIR
(susceptible-exposed-infectious-recovered) epidemiological model embedded in an
optimization framework and examine scenarios involving a single donor and
multiple recipient (nondonor) geographic areas. We find that policies other
than donor-first can delay the emergence of a more-contagious variant compared
to donor-first, sometimes reducing donor-country deaths in addition to total
deaths. Thus, vaccine distribution is not a zero-sum game between donor and
nondonor countries: an optimization approach can achieve dramatic reduction in
total deaths with only a small increase (and occasionally even a decrease) in
donor-country deaths. We also identify realistic scenarios under which the
optimal policy has a switching form rather than adhering to a strict priority
order. The iterative linear programming approximation approach we develop can
help confirm those instances when a priority policy is optimal and, when not
optimal, can identify superior policies.


  The spread of COVID-19 has brought challenges to health, social and economic
systems around the world. With little to no prior immunity in the global
population transmission has been driven primarily by human interaction.
However, as with common respiratory illnesses such as the flu it's suggested
that COVID-19 may become seasonal as immunity grows. Yet the effects of
meteorological conditions on the spread of COVID-19 are poorly understood with
previous studies producing contrasting results, due at least in part to limited
and inconsistent study designs. This study investigates the effect of
meteorological conditions on COVID-19 infections in England using a
spatio-temporal model applied to case counts during the initial England
lockdown. By modelling spatial and temporal effects to account for the nature
of a human transmissible virus the model isolates meteorological effects.
Inference based on 95% highest posterior density intervals shows humidity is
negatively associated with COVID-19 spread. The lack of evidence for other
weather factors affecting COVID-19 transmission shows care should be taken with
respect to seasonality when designing COVID-19 policies and public
communications.


  The foraging behavior of animals is a paradigm of target search in nature.
Understanding which foraging strategies are optimal and how animals learn them
are central challenges in modeling animal foraging. While the question of
optimality has wide-ranging implications across fields such as economy,
physics, and ecology, the question of learnability is a topic of ongoing debate
in evolutionary biology. Recognizing the interconnected nature of these
challenges, this work addresses them simultaneously by exploring optimal
foraging strategies through a reinforcement learning framework. To this end, we
model foragers as learning agents. We first prove theoretically that maximizing
rewards in our reinforcement learning model is equivalent to optimizing
foraging efficiency. We then show with numerical experiments that, in the
paradigmatic model of non-destructive search, our agents learn foraging
strategies which outperform the efficiency of some of the best known strategies
such as L\'evy walks. These findings highlight the potential of reinforcement
learning as a versatile framework not only for optimizing search strategies but
also to model the learning process, thus shedding light on the role of learning
in natural optimization processes.


  Modern high throughput sequencing technologies like metagenomic sequencing
generate millions of sequences which have to be classified based on their
taxonomic rank. Modern approaches either apply local alignment and comparison
to existing data sets like MMseqs2 or use deep neural networks as it is done in
DeepMicrobes and BERTax. Alignment-based approaches are costly in terms of
runtime, especially since databases get larger and larger. For the deep
learning-based approaches, specialized hardware is necessary for a computation,
which consumes large amounts of energy. In this paper, we propose to use
$k$-mer distributions obtained from DNA as features to classify its taxonomic
origin using machine learning approaches like the subspace $k$-nearest
neighbors algorithm, neural networks or bagged decision trees. In addition, we
propose a feature space data set balancing approach, which allows reducing the
data set for training and improves the performance of the classifiers. By
comparing performance, time, and memory consumption of our approach to those of
state-of-the-art algorithms (BERTax and MMseqs2) using several datasets, we
show that our approach improves the classification on the genus level and
achieves comparable results for the superkingdom and phylum level.
  Link:
https://es-cloud.cs.uni-tuebingen.de/d/8e2ab8c3fdd444e1a135/?p=%2FTaxonomyClassification&mode=list


  The R Package IBMPopSim (https://daphnegiorgi.github.io/IBMPopSim/) aims to
simulate the random evolution of heterogeneous populations using stochastic
Individual-Based Models (IBMs). The package enables users to simulate
population evolution, in which individuals are characterized by their age and
some characteristics, and the population is modified by different types of
events, including births/arrivals, death/exit events, or changes of
characteristics. The frequency at which an event can occur to an individual can
depend on their age and characteristics, but also on the characteristics of
other individuals (interactions). Such models have a wide range of
applications. For instance, IBMs can be used for simulating the evolution of a
heterogeneous insurance portfolio with selection or for validating mortality
forecasts. IBMPopSim overcomes the limitations of time-consuming IBMs
simulations by implementing new efficient algorithms based on thinning methods,
which are compiled using the Rcpp package while providing a user-friendly
interface.


  Pleurotus ostreatus is a common cultivated edible mushroom worldwide. The
fruiting bodies of P. ostreatus is a rich source of a $\beta$-glucans
polysaccharide. The current study aimed to investigate the effectiveness of
$\beta$-glucans as a natural polysaccharide produced by P. ostreatus as an
antioxidant, antimicrobial, and anticancer. The molecular identification of P.
ostreatus isolate was confirmed by Internal Transcribed Spacer (ITS) sequence.
The sequence alignment and phylogenetic evolutionary relationship of studied
ITS sequence were performed against some deposited sequences in GenBank. The
analysis of high-performance liquid chromatography (HPLC) as well as the result
of fourier transform infrared spectroscopy (FTIR) has confirmed the presence of
$\beta$-glucans polysaccharide in the tested samples. The percentage of
antioxidant activity of $\beta$-glucans showed a gradual increase from 8.59% to
12.36, 18.56, 23.69, 44.66 and 80.36% at the concentrations of 31.2, 64.4, 125,
250, 500, and 800 $\mu$g/ml, respectively. In addition, all concentrations of
$\beta$-glucans showed higher antioxidant activities when compared with
standard antioxidant (Vitamin C). The highest antimicrobial activity of
$\beta$-glucans polysaccharide was against P. aeruginosa with a zone of
inhibition (45 mm), while the lowest activity was against S. aureus (13 mm)
both at 100 mg/mL. The percentage of growth-inhibiting of MCF-7 a humanbreast
cancer cell line and normal WRL-68 cell line affected by $\beta$-glucans were
determined by 3-(4,5)-dimethylthiazol (-z-y1)-3,5-di-phenytetrazoliumromide
(MTT assay).