  Single cell combinatorial indexing RNA sequencing (sci-RNA-seq) is a powerful
method for recovering gene expression data from an exponentially scalable
number of individual cells or nuclei. However, sci-RNA-seq is a complex
protocol that has historically exhibited variable performance on different
tissues, as well as lower sensitivity than alternative methods. Here we report
a simplified, optimized version of the three-level sci-RNA-seq protocol that is
faster, higher yield, more robust, and more sensitive, than the original
sci-RNA-seq3 protocol, with reagent costs on the order of 1 cent per cell or
less. We showcase the optimized protocol via whole organism analysis of an
E16.5 mouse embryo, profiling ~380,000 nuclei in a single experiment. Finally,
we introduce a "tiny sci-*" protocol for experiments where input is extremely
limited.


  The field of machine learning for drug discovery is witnessing an explosion
of novel methods. These methods are often benchmarked on simple physicochemical
properties such as solubility or general druglikeness, which can be readily
computed. However, these properties are poor representatives of objective
functions in drug design, mainly because they do not depend on the candidate's
interaction with the target. By contrast, molecular docking is a widely
successful method in drug discovery to estimate binding affinities. However,
docking simulations require a significant amount of domain knowledge to set up
correctly which hampers adoption. To this end, we present DOCKSTRING, a bundle
for meaningful and robust comparison of ML models consisting of three
components: (1) an open-source Python package for straightforward computation
of docking scores; (2) an extensive dataset of docking scores and poses of more
than 260K ligands for 58 medically-relevant targets; and (3) a set of
pharmaceutically-relevant benchmark tasks including regression, virtual
screening, and de novo design. The Python package implements a robust ligand
and target preparation protocol that allows non-experts to obtain meaningful
docking scores. Our dataset is the first to include docking poses, as well as
the first of its size that is a full matrix, thus facilitating experiments in
multiobjective optimization and transfer learning. Overall, our results
indicate that docking scores are a more appropriate evaluation objective than
simple physicochemical properties, yielding more realistic benchmark tasks and
molecular candidates.


  Multi-trait genome-wide association studies (GWAS) use multi-variate
statistical methods to identify associations between genetic variants and
multiple correlated traits simultaneously, and have higher statistical power
than independent univariate analyses of traits. Reverse regression, where
genotypes of genetic variants are regressed on multiple traits simultaneously,
has emerged as a promising approach to perform multi-trait GWAS in
high-dimensional settings where the number of traits exceeds the number of
samples. We analyzed different machine learning methods (ridge regression,
naive Bayes/independent univariate, random forests and support vector machines)
for reverse regression in multi-trait GWAS, using genotypes, gene expression
data and ground-truth transcriptional regulatory networks from the DREAM5
SysGen Challenge and from a cross between two yeast strains to evaluate
methods. We found that genotype prediction performance, in terms of root mean
squared error (RMSE), allowed to distinguish between genomic regions with high
and low transcriptional activity. Moreover, model feature coefficients
correlated with the strength of association between variants and individual
traits, and were predictive of true trans-eQTL target genes, with complementary
findings across methods. Code to reproduce the analysis is available at
https://github.com/michoel-lab/Reverse-Pred-GWAS


  We explore the relationship between Eulerian and Lagrangian approaches for
modeling movement in vector-borne diseases for discrete space. In the Eulerian
approach we account for the movement of hosts explicitly through movement rates
captured by a graph Laplacian matrix $L$. In the Lagrangian approach we only
account for the proportion of time that individuals spend in foreign patches
through a mixing matrix $P$. We establish a relationship between an Eulerian
model and a Lagrangian model for the hosts in terms of the matrices $L$ and
$P$. We say that the two modeling frameworks are consistent if for a given
matrix $P$, the matrix $L$ can be chosen so that the residence times of the
matrix $P$ and the matrix $L$ match. We find a sufficient condition for
consistency, and examine disease quantities such as the final outbreak size and
basic reproduction number in both the consistent and inconsistent cases. In the
special case of a two-patch model, we observe how similar values for the basic
reproduction number and final outbreak size can occur even in the inconsistent
case. However, there are scenarios where the final sizes in both approaches can
significantly differ by means of the relationship we propose.


  Gene expression via transcription-translation is the most fundamental
reaction to sustain biological systems, and complex reactions such as this one
occur in a small compartment of living cells. There is increasing evidence that
t physical effects, such as molecular crowding or excluded volume effects of
transcriptional-translational machinery, affect the yield of reaction products.
On the other hand, transcriptional feedback that controls gene expression
during mRNA synthesis is also a vital mechanism that regulates protein
synthesis in cells. However, the excluded volume effect of spatial constraints
on feedback regulation is not well understood. Here, we study the confinement
effect on transcriptional autoregulatory feedbacks of gene expression reactions
using a theoretical model. The excluded volume effects between molecules and
the membrane interface suppress the gene expression in a small cell-sized
compartment. We find that negative feedback regulation at the transcription
step mitigates this size-induced gene repression and alters the scaling
relation of gene expression level on compartment volume, approaching the
regular scaling relation without the steric effect. This recovery of regular
size-scaling of gene expression does not appear in positive feedback
regulation, suggesting that negative autoregulatory feedback is crucial for
maintaining reaction products constant regardless of compartment size in
heterogeneous cell populations.


  Background:Gene regulations often change over time rather than being
constant. But many of gene regulatory networks extracted from databases are
static. The tumor suppressor gene $P53$ is involved in the pathogenesis of many
tumors, and its inhibition effects occur after a certain period. Therefore, it
is of great significance to elucidate the regulation mechanism over time
points. Result:A qualitative method for representing dynamic gene regulatory
network is developed, called CdtGRN. It adopts the combination of convolutional
neural networks(CNN) and fully connected networks(DNN) as the core mechanism of
prediction. The ionizing radiation Affymetrix dataset (E-MEXP-549) was obtained
at ArrayExpress, by microarray gene expression levels predicting relations
between regulation. CdtGRN is tested against a time-delayed gene regulatory
network with $22,284$ genes related to $P53$. The accuracy of CdtGRN reaches
92.07$\%$ on the classification of conservative verification set, and a kappa
coefficient reaches $0.84$ and an average AUC accuracy is 94.25$\%$. This
resulted in the construction of. Conclusion:The algorithm and program we
developed in our study would be useful for identifying dynamic gene regulatory
networks, and objectively analyze the delay of the regulatory relationship by
analyzing the gene expression levels at different time points. The time-delayed
gene regulatory network of $P53$ is also inferred and represented
qualitatively, which is helpful to understand the pathological mechanism of
tumors.


  This study included isolate petroleum hydrocarbons degradable bacteria and
develops a consortium or a mixture of bacteria with high biodegradation
capabilities which can be used in biological treatment units of the
contaminated soils before release. In this studyt ten bacterial strains were
isolated from soils contaminated with crude oil, by primary and secondary
screening, and by using the sterile saline solution with 1% of crude oil, five
bacterial isolates that can degrade oil were identified. The bacterial isolates
were isolated from polluted soils with crude oil, the samples were collected
from different areas of the Baiji refinery, and samples of heavy crude oil
extracted from Qayyarah fields were used in designing the biodegradation
experiments, the five isolates were diagnosed based on phenotypic, culture and
biochemical characterizes, and detection of the gene sequence of bacterial
isolates (16SrRNA). The (16SrRNA) gene sequence of the bacterial isolates was
recorded under the accession numbers LC596402, LC596403, LC596406, LC596404,
LC596405, for the genus (AM-I-1, AM-I-2, AM-I-5) and the species (AM-I-3,
AM-I-4) that following Bacillus Sp respectively, in the NCBI's GenBank, the
efficiency of the five isolates were examined for utilizing petroleum
hydrocarbons using a sterile mineral salt medium MS supplemented with crude oil
as the sole source of carbon with different concentrations, the results showed
that the decomposition of petroleum hydrocarbons at the concentrations (0.5, 1,
1.5, 2, 3) % reached (57.32, 69.36,63.71, 75.20, 68.60) %, respectively, for
mixed isolates, depending on the results of the gas chromatography (GC)
analysis.


  The articulatory-acoustic relationship is many-to-one and non linear and this
is a great limitation for studying speech production. A simplification is
proposed to set a bijection between the vowel space (f1, f2) and the parametric
space of different vocal tract models. The generic area function model is based
on mixtures of cosines allowing the generation of main vowels with two
formulas. Then the mixture function is transformed into a coordination function
able to deal with articulatory parameters. This is shown that the coordination
function acts similarly with the Fant's model and with the 4-Tube DRM derived
from the generic model.


  Searching for novel molecular compounds with desired properties is an
important problem in drug discovery. Many existing frameworks generate
molecules one atom at a time. We instead propose a flexible editing paradigm
that generates molecules using learned molecular fragments--meaningful
substructures of molecules. To do so, we train a variational autoencoder (VAE)
to encode molecular fragments in a coherent latent space, which we then utilize
as a vocabulary for editing molecules to explore the complex chemical property
space. Equipped with the learned fragment vocabulary, we propose Fragment-based
Sequential Translation (FaST), which learns a reinforcement learning (RL)
policy to iteratively translate model-discovered molecules into increasingly
novel molecules while satisfying desired properties. Empirical evaluation shows
that FaST significantly improves over state-of-the-art methods on benchmark
single/multi-objective molecular optimization tasks.


  In this research, we provide a mathematical analysis for the novel
coronavirus responsible for COVID-19, which continues to be a big source of
threat for humanity. Our fractional-order analysis is carried out using a
non-singular kernel type operator known as the Atangana--Baleanu--Caputo (ABC)
derivative. We parametrize the model adopting available information of the
disease from Pakistan in the period 9th April to 2nd June 2020. We obtain the
required solution with the help of a hybrid method, which is a combination of
the decomposition method and the Laplace transform. Furthermore, a sensitivity
analysis is carried out to evaluate the parameters that are more sensitive to
the basic reproduction number of the model. Our results are compared with the
real data of Pakistan and numerical plots are presented at various fractional
orders.


  In this article we present a new p-adic generalization of the Eigen-Schuster
model where the genomes (sequences) are represented by words written in the
alphabet 0,1,,p-1, where p is a prime number, with a time variable length. The
time evolution of the concentration of a sequence is controlled by a p-adic
evolution equation. The long term behavior of the concentration of a sequence
depends on a fitness function f, a mutation measure Q, and an initial
concentration distribution. The new model provides essentially two types of
asymptotic scenarios for evolution. If the complexity of sequences grows at the
right pace, then in the long term the survival is assured. This agrees with the
fact that larger genome size improves the replication fidelity. In other case,
the sequences cannot copy themselves with sufficiently fidelity, and in the
long term they will not survive. Eigen's paradox is one, among the infinitely
many, possible scenarios of the evolution in the long term. The mathematical
formulation of this fact requires solving the Cauchy problem for the p-adic
Eigen-Schuster model in a rigorous mathematical way. This requires imposing
restrictions on the fitness function and on the mutation measure, among other
conditions. The study of the mentioned initial value problem requires
techniques of p-adic wavelets and heat kernels developed in the last 35 years.


  Genome sequence analysis plays a pivotal role in enabling many medical and
scientific advancements in personalized medicine, outbreak tracing, and
forensics. However, the analysis of genome sequencing data is currently
bottlenecked by the computational power and memory bandwidth limitations of
existing systems. In this dissertation, we propose four major works, where we
characterize the real-system behavior of the genome sequence analysis pipeline
and its associated tools, expose the bottlenecks and tradeoffs, and co-design
fast and efficient algorithms along with scalable and energy-efficient
customized hardware accelerators for the key bottlenecks to enable faster
genome sequence analysis.
  First, we comprehensively analyze the tools in the genome assembly pipeline
for long reads in multiple dimensions, uncovering bottlenecks and tradeoffs
that different combinations of tools and different underlying systems lead to.
Second, we propose GenASM, an acceleration framework that builds upon
bitvector-based approximate string matching to accelerate multiple steps of the
genome sequence analysis pipeline. We co-design our highly-parallel, scalable
and memory-efficient algorithms with low-power and area-efficient hardware
accelerators. Third, we implement an FPGA-based prototype for GenASM, where
state-of-the-art 3D-stacked memory offers high memory bandwidth and FPGA
resources offer high parallelism. Fourth, we propose SeGraM, the first hardware
acceleration framework for sequence-to-graph mapping and alignment. We
co-design algorithms and accelerators for memory-efficient minimizer-based
seeding and bitvector-based, highly-parallel sequence-to-graph alignment.
  Overall, we demonstrate that genome sequence analysis can be accelerated by
co-designing scalable and energy-efficient customized accelerators along with
efficient algorithms for the key steps of genome sequence analysis.


  In this work we present a system based on a Deep Learning approach, by using
a Convolutional Neural Network, capable of classifying protein chains of amino
acids based on the protein description contained in the Protein Data Bank. Each
protein is fully described in its chemical-physical-geometric properties in a
file in XML format. The aim of the work is to design a prototypical Deep
Learning machinery for the collection and management of a huge amount of data
and to validate it through its application to the classification of a sequences
of amino acids. We envisage applying the described approach to more general
classification problems in biomolecules, related to structural properties and
similarities.


  In this work we set out to find a method to classify protein structures using
a Deep Learning methodology. Our Artificial Intelligence has been trained to
recognize complex biomolecule structures extrapolated from the Protein Data
Bank (PDB) database and reprocessed as images; for this purpose various tests
have been conducted with pre-trained Convolutional Neural Networks, such as
InceptionResNetV2 or InceptionV3, in order to extract significant features from
these images and correctly classify the molecule. A comparative analysis of the
performances of the various networks will therefore be produced.


  We develop an analysis of the cryptocurrency market borrowing methods and
concepts from ecology. This approach makes it possible to identify specific
diversity patterns and their variation, in close analogy with ecological
systems, and to characterize the cryptocurrency market in an effective way. At
the same time, it shows how non-biological systems can have an important role
in contrasting different ecological theories and in testing the use of neutral
models. The study of the cryptocurrencies abundance distribution and the
evolution of the community structure strongly indicates that these statistical
patterns are not consistent with neutrality. In particular, the necessity to
increase the temporal change in community composition when the number of
cryptocurrencies grows, suggests that their interactions are not necessarily
weak. The analysis of the intraspecific and interspecific interdependency
supports this fact and demonstrates the presence of a market sector influenced
by mutualistic relations. These latest findings challenge the hypothesis of
weakly interacting symmetric species, the postulate at the heart of neutral
models.


  We formulate a general age-of-infection epidemic model with two pathways: the
symptomatic infections and the asymptomatic infections. We then calculate the
basic reproduction number $\mathcal{R}_0$ and establish the final size
relation. It is shown that the ratio of accumulated counts of symptomatic
patients and asymptomatic patients is determined by the symptomatic ratio $f$
which is defined as the probability of eventually becoming symptomatic after
being infected. We also formulate and study a general age-of-infection model
with disease deaths and with two infection pathways. The final size relation is
investigated, and the upper and lower bounds for final epidemic size are given.
Several numerical simulations are performed to verify the analytical results.


  Network science has increasingly become central to the field of epidemiology
and our ability to respond to infectious disease threats. However, many
networks derived from modern datasets are not just large, but dense, with a
high ratio of edges to nodes. This includes human mobility networks where most
locations have a large number of links to many other locations. Simulating
large-scale epidemics requires substantial computational resources and in many
cases is practically infeasible. One way to reduce the computational cost of
simulating epidemics on these networks is sparsification, where a
representative subset of edges is selected based on some measure of their
importance. We test several sparsification strategies, ranging from naive
thresholding to random sampling of edges, on mobility data from the U.S.
Following recent work in computer science, we find that the most accurate
approach uses the effective resistances of edges, which prioritizes edges that
are the only efficient way to travel between their endpoints. The resulting
sparse network preserves many aspects of the behavior of an SIR model,
including both global quantities, like the epidemic size, and local details of
stochastic events, including the probability each node becomes infected and its
distribution of arrival times. This holds even when the sparse network
preserves fewer than $10\%$ of the edges of the original network. In addition
to its practical utility, this method helps illuminate which links of a
weighted, undirected network are most important to disease spread.


  Short peptides with antimicrobial activity have therapeutic potential for
treating bacterial infections. Mechanisms of actions for antimicrobial peptides
require binding the biological membrane of their target, which often represents
a key mechanistic step. A multitude of data-driven approaches have been
developed to predict potential antimicrobial peptide sequences; however, these
methods are usually agnostic to the physical interactions between the peptide
and the membrane. Towards developing higher throughput screening methodologies,
here we use Markov State Modeling and all-atom molecular dynamics simulations
to quantify the membrane binding and insertion kinetics of three prototypical
and antimicrobial peptides (alpha-helical magainin 2 and PGLa and beta-hairpin
tachyplesin 1). By leveraging a set of collective variables that capture the
essential physics of the amphiphilic and cationic peptide-membrane interactions
we reveal how the slowest kinetic process of membrane insertion is the dynamic
rolling of the peptide from a prebound to fully inserted state. These results
add critical details to how antimicrobial peptides insert into bacterial
membranes.


  Transmission models for infectious diseases are typically formulated in terms
of dynamics between individuals or groups with processes such as disease
progression or recovery for each individual captured phenomenologically,
without reference to underlying biological processes. Furthermore, the
construction of these models is often monolithic: they don't allow one to
readily modify the processes involved or include the new ones, or to combine
models at different scales. We show how to construct a simple model of immune
response to a respiratory virus and a model of transmission using an easily
modifiable set of rules allowing further refining and merging the two models
together. The immune response model reproduces the expected response curve of
PCR testing for COVID-19 and implies a long-tailed distribution of
infectiousness reflective of individual heterogeneity. This immune response
model, when combined with a transmission model, reproduces the previously
reported shift in the population distribution of viral loads along an epidemic
trajectory.


  Investigation of age-related genes is of great importance for multiple
purposes, for instance, improving our understanding of the mechanism of ageing,
increasing life expectancy, age prediction, and other healthcare applications.
In his work, starting with a set of 27,142 genes, we develop an
information-theoretic framework for identifying genes that are associated with
aging by applying unsupervised and semi-supervised learning techniques on human
dermal fibroblast gene expression data. First, we use unsupervised learning and
apply information-theoretic measures to identify key features for effective
representation of gene expression values in the transcriptome data. Using the
identified features, we perform clustering on the data. Finally, we apply
semi-supervised learning on the clusters using different distance measures to
identify novel genes that are potentially associated with aging. Performance
assessment for both unsupervised and semi-supervised methods show the
effectiveness of the framework.


  The study of signatures of aging in terms of genomic biomarkers can be
uniquely helpful in understanding the mechanisms of aging and developing models
to accurately predict the age. Prior studies have employed gene expression and
DNA methylation data aiming at accurate prediction of age. In this line, we
propose a new framework for human age estimation using information from human
dermal fibroblast gene expression data. First, we propose a new spatial
representation as well as a data augmentation approach for gene expression
data. Next in order to predict the age, we design an architecture of neural
network and apply it to this new representation of the original and augmented
data, as an ensemble classification approach. Our experimental results suggest
the superiority of the proposed framework over state-of-the-art age estimation
methods using DNA methylation and gene expression data.


  The Black Death is regarded as a turning point in late medieval European
history. Recent studies have shown that even regions that have so far been
perceived in the literature as not or only marginally affected by the epidemic,
suffered from its profound demographic and economic consequences. The scale and
geographical range of the plague in Central Europe, the Kingdom of Poland
included remains, however, a matter of dispute and from the beginning scholar's
views on this matter have been divided. What is particularly important, the
outbreak of the plague in Western Europe coincided with the reign of Casimir of
the Piast dynasty, the only ruler of Poland to receive the nickname the Great,
who is associated with the modernization and extraordinary development of his
kingdom which was entering the golden age of its history.


  Recent developments in machine-learning based molecular fragment linking have
demonstrated the importance of informing the generation process with structural
information specifying the relative orientation of the fragments to be linked.
However, such structural information has not yet been provided in the form of a
complete relative coordinate system. Mathematical details for a decoupled set
of bond lengths, bond angles and torsion angles are elaborated and the
coordinate system is demonstrated to be complete. Significant impact on the
quality of the generated linkers is demonstrated numerically. The amount of
reliable information within the different types of degrees of freedom is
investigated. Ablation studies and an information-theoretical analysis are
performed. The presented benefits suggest the application of a complete and
decoupled relative coordinate system as a standard good practice in linker
design.


  During epidemics, estimation of the effective reproduction number (ERN)
associated with infectious disease is a challenging topic for policy
development and medical resource management. There is still an open question
about the dominant factors to characterize in corona virus disease 2019
(COVID-19), although recent studies based on nonlinear regression with machine
learning suggested mobility. The emergence of new viral variants is common in
widespread pandemics. However, understanding the potential ERN of new variants
is required for policy revision, including lockdown constraints. In this study,
we proposed time-averaged mobility at transit stations as a surrogate to
correlate with ERN using data from three prefectures in Japan. The latency and
duration to average over the mobility were 6-8 days and 6-7, respectively (R2
was 0.109-0.512 in Tokyo, 0.365-0.607 in Osaka, and 0.317-0.631 in Aichi). The
same linear correlation was confirmed in Singapore and London. The
mobility-adjusted ERN of the alpha variant was 15%-30%, and was 20%-40% higher
than the standard type in Osaka, Aichi, and London. Similarly, the ERN of the
delta variant was 20%-40% higher than that of the standard type in Osaka and
Aichi. The proposed metric would be useful for proper evaluation of the
infectivity of different variants in terms of ERN.


  In individual SARS-CoV-2 outbreaks, the count of confirmed cases and deaths
follow a Gompertz growth function for locations of very different sizes. This
lack of dependence on region size leads us to hypothesize that virus spread
depends on universal properties of the network of social interactions. We test
this hypothesis by simulating the propagation of a virus on networks of
different topologies. Our main finding is that Gompertz growth observed for
early outbreaks occurs only for a scale-free network, in which nodes with many
more neighbors than average are common. These nodes that have very many
neighbors are infected early in the outbreak and then spread the infection very
rapidly. When these nodes are no longer infectious, the remaining nodes that
have most neighbors take over and continue to spread the infection. In this
way, the rate of spread is fastest at the very start and slows down
immediately. Geometrically it is seen that the "surface" of the epidemic, the
number of susceptible nodes in contact with the infected nodes, starts to
rapidly decrease very early in the epidemic and as soon as the larger nodes
have been infected. In our simulation, the speed and impact of an outbreak
depend on three parameters: the average number of contacts each node makes, the
probability of being infected by a neighbor, and the probability of recovery.
Intelligent interventions to reduce the impact of future outbreaks need to
focus on these critical parameters in order to minimize economic and social
collateral damage.


  Proteins are the common constituents of all living cells. They are molecular
machines that interact with each other as well as with other cell products and
carry out a dizzying array of functions with distinction. These interactions
follow from their native state structures and therefore understanding
sequence-structure relationships is of fundamental importance. What is quite
remarkable about proteins is that their understanding necessarily straddles
several disciplines. The importance of geometry in defining protein native
state structure, the constraints placed on protein behavior by mathematics and
physics, the need for proteins to obey the laws of quantum chemistry, and the
rich role of evolution and biology all come together in defining protein
science. Here we review ideas from the literature and present an
interdisciplinary framework that aims to marry ideas from Plato and Darwin and
demonstrates an astonishing consilience between disciplines in describing
proteins. We discuss the consequences of this framework on protein behavior.


  Controlling an epidemiological model is often performed using optimal control
theory techniques for which the solution depends on the equations of the
control system, objective functional and possible state and/or control
constraints. In this paper, we propose a model-free control approach based on
an algorithm that operates in 'real-time' and drives the state solution
according to a direct feedback on the state solution that is aimed to be
minimized, and without knowing explicitly the equations of the control system.
We consider a concrete epidemic problem of minimizing the number of HIV
infected individuals, through the preventive measure pre-exposure prophylaxis
(PrEP) given to susceptible individuals. The solutions must satisfy control and
mixed state-control constraints that represent the limitations on PrEP
implementation. Our model-free based control algorithm allows to close the loop
between the number of infected individuals with HIV and the supply of PrEP
medication 'in real time', in such manner that the number of infected
individuals is asymptotically reduced and the number of individuals under PrEP
medication is below a fixed constant value. We prove the efficiency of our
approach and compare the model-free control solutions with the ones obtained
using a classical optimal control approach via Pontryagin maximum principle.
The performed numerical simulations allow us to conclude that the model-free
based control strategy highlights new and interesting performances compared
with the classical optimal control approach.


  The physical-chemical property of molecular collision cross section (CCS) is
increasingly used to assist in small molecule identification; however,
questions remain regarding the extent of its true utility in contributing to
such identifications, especially given its correlation with mass. To
investigate the contribution of CCS to uniqueness within a given library, we
measured its discriminatory capacity as a function of error in CCS values (from
measurement or prediction), CCS variance, parent mass, mass error, and/or
reference database size using a multi-directional grid search. While
experimental CCS databases exist, they are currently small; thus, we used a CCS
prediction tool, DarkChem, to provide theoretical CCS values for use in this
study. These predicted CCS values were then modified to mirror experimental
variance. By augmenting our search within a library based on mass alone with
CCS at a variety of accuracies, we found that, (i) the use of multiple adducts
(i.e. alternative ionized forms of the same parent compound) for the same
molecule, compared to using a single adduct, greatly improves specificity and
(ii) even a single CCS leads to a significant specificity boost when low CCS
error (e.g. 1% composite error) can be achieved. Based on these results, we
recommend using multiple adducts to build up evidence of presence, as each
adduct supplies additional information per dimension. Additionally, the utility
of ion mobility spectrometry when coupled with mass spectrometry should still
be considered, regardless of whether CCS is considered as an identification
metric, due to advantages such as increased peak resolution, sensitivity (e.g.
from reducing load on the detector at any given time), improvements in
data-independent MS/MS spectra acquisition, and cleaner tandem mass spectral
fragmentation patterns.


  Accurate molecular identification of metabolites can unlock new areas of the
molecular universe and allow greater insight into complex biological and
environmental systems than currently possible. Analytical approaches for
measuring the metabolome, such as NMR spectroscopy, and separation techniques
coupled with mass spectrometry, such as LC-IMS-MS, have risen to this challenge
by yielding rich experimental data that can be queried by cross-reference with
similar information for known standards in reference libraries. Confident
identification of molecules in metabolomics studies, though, is often limited
by the diversity of available data across chemical space, the unavailability of
authentic reference standards, and the corresponding lack of comprehensiveness
of standard reference libraries. The In Silico Chemical Library Engine (ISiCLE)
addresses theses hindrances by providing a first-principles, cheminformatics
pipeline that yields collisional cross section (CCS) values for any given
molecule and without the need for training data. In this program, chemical
identifiers undergo MD simulations, quantum chemical transformations, and ion
mobility calculations for the generation of predicted CCS values. Here, we
present a new module for ISiCLE that addresses the sensitivity of CCS
predictions to ionization site location. An update to adduct creation methods
is proposed concerning a transition from pKa and pKb led predictions to a
Gibb's free energy (GFE) based determinacy of true ionization site location. A
validation set of experimentally confirmed molecular protonation sites was
assembled from literature and cross-referenced with the respective pKb
predicted locations and GFE values for all potential ionization site
placements. Upon evaluation of the two methods, the lowest GFE value was found
to predict the true ionization site location with 100% accuracy while pKb had
less accuracy.


  In this study, we propose a time-dependent
Susceptible-Exposed-Infected-Recovered (SEIR) model for the analysis of the
SARS-CoV-2 epidemic outbreak in three different countries, the United States of
America, Italy and Iceland using public data inherent the numbers of the
epidemic wave. Since several types and grades of actions were adopted by the
governments, including travel restrictions, social distancing, or limitation of
movement, we want to investigate how these measures can affect the epidemic
curve of the infectious population. The parameters of interest for the SEIR
model were estimated employing a composite likelihood approach. Moreover,
standard errors have been corrected for temporal dependence. The adoption of
restrictive measures results in flatten epidemic curves, and the future
evolution indicated a decrease in the number of cases.


  Functional groups and moieties are chemical descriptors of biomolecules that
can be used to interpret their properties and functions, leading to the
understanding of chemical or biological mechanisms. These chemical building
blocks, or sub-structures, enable the identification of common molecular
subgroups, assessing the structural similarities and critical interactions
among a set of biological molecules with known activities, and designing novel
compounds with similar chemical properties. Here, we introduce a Python-based
tool, SPECTRe (Substructure Processing, Enumeration, and Comparison Tool
Resource), designed to provide all substructures in a given molecular
structure, regardless of the molecule size, employing efficient enumeration and
generation of substructures represented in a human-readable SMILES format
through the use of classical graph traversal (breadth-first and depth-first
search) algorithms. We demonstrate the application of SPECTRe for a set of
10,375 molecules in the molecular weight range 27 to 350 Da (<=26 non-hydrogen
atoms), spanning a wide array of structure-based chemical functionalities and
chemical classes. We found that the substructure count as a measure of
molecular complexity depends strongly on the number of unique atom and bond
types present, degree of branching, and presence of rings. The substructure
counts are found to be similar for a set of molecules belonging to particular
chemical classes and classified based on the characteristic features of certain
topologies. We demonstrate that SPECTRe shows promise to be useful in many
applications of cheminformatics such as virtual screening for drug discovery,
property prediction, fingerprint-based molecular similarity searching, and data
mining for identifying frequent substructures.


  In this study, we apply a novel combination of close proximity and room-scale
risk assessment approaches for people sharing public transport environments to
predict their contagion risk due to SARS-CoV-2 respiratory infection. In
particular, the individual infection risk of susceptible subjects and the
transmissibility of SARS-CoV-2 (expressed through the reproduction number) are
evaluated for two types of buses, differing in terms of exposure time and
crowding index: urban and long-distance buses. Infection risk and reproduction
number are calculated for different scenarios as a function of the ventilation
rates (both measured and estimated according to standards), crowding indexes,
and travel times. The results show that for urban buses, the close proximity
contribution significantly affects the maximum occupancy to maintain a
reproductive number of < 1. In particular, full occupancy of the bus would be
permitted only for an infected subject breathing, whereas for an infected
subject speaking, masking would be required. For long-distance buses, full
occupancy of the bus can be maintained only if specific mitigation solutions
are simultaneously applied. For example, for an infected person speaking for 1
h, appropriate filtration of the recirculated air and simultaneous use of FFP2
masks would permit full occupancy of the bus for a period of almost 8 h.
Otherwise, a high percentage of immunized persons (> 80%) would be needed.


  This article addresses the emergence of collective agent behaviors from
simple and realistic decisions made at the level of the local Voronoi diagram,
i.e., the Voronoi region (personal space) of agents and their neighboring
agents in the sense of the Voronoi diagram. We present a discrete time model in
2D in which individual agents are aware of their local Voronoi environment and
may seek static target locations. Our deterministic model incorporates features
common to many other "active matter" models like collision avoidance, alignment
among agents, and homing toward targets. However, it is novel in key respects:
the model combines topological and metrical features in a natural manner based
upon the agents local perception of their Voronoi cell. With at most two
control parameters, it is shown numerically to capture a wide range of
collective behaviors. It is also shown that the Voronoi topology facilitates
the computation of several novel observables for quantifying discrete
collective behaviors.


  A general problem of current interest is the analysis of diffusion problems
in singularly perturbed domains, within which small subdomains are removed from
the domain interior and boundary conditions imposed on the resulting holes. One
major application is to intracellular diffusion, where the holes could
represent organelles or biochemical substrates. In this paper we use a
combination of matched asymptotic analysis and Green's function methods to
calculate the so-called accumulation time for relaxation to steady state in 2D
domains. The standard measure of the relaxation rate is in terms of the
principal nonzero eigenvalue of the negative Laplacian. However, this global
measure does not account for possible differences in the relaxation rate at
different spatial locations, is independent of the initial conditions, and
relies on the assumption that the eigenvalues have sufficiently large spectral
gaps. As previously established for diffusion-based morphogen gradient
formation, the accumulation time provides a better measure of the relaxation
process.


  Motivated by empirical arguments that are well-known from the genome-wide
association studies (GWAS) literature, we study the statistical properties of
linear mixed models (LMMs) applied to GWAS. First, we study the sensitivity of
LMMs to the inclusion of a candidate SNP in the kinship matrix, which is often
done in practice to speed up computations. Our results shed light on the size
of the error incurred by including a candidate SNP, providing a justification
to this technique in order to trade-off velocity against veracity. Second, we
investigate how mixed models can correct confounders in GWAS, which is widely
accepted as an advantage of LMMs over traditional methods. We consider two
sources of confounding factors, population stratification and environmental
confounding factors, and study how different methods that are commonly used in
practice trade-off these two confounding factors differently.


  Species survival in the $(3, 1)$ May-Leonard system is determined by the
mobility, with a critical mobility threshold between long-term coexistence and
extinction. We show experimentally that the critical mobility threshold is
determined by the topology of the graph, with the critical mobility threshold
of the periodic lattice twice that of a topological 2-sphere. We use
topological techniques to monitor the evolution of patterns in various graphs
and estimate their mean persistence time, and show that the difference in
critical mobility threshold is due to a specific pattern that can form on the
lattice but not on the sphere. Finally, we release the software toolkit we
developed to perform these experiments to the community.


  Difficulties identifying appropriate biodiversity impact metrics remain a
major barrier to inclusion of biodiversity considerations in environmentally
responsible investment. We propose and analyse a simple science-based local
metric: the sum of proportional changes in local species abundances relative to
their global species abundances, with a correction for species close to
extinction. As we show, this metric quantifies changes in the mean long-term
global survival probability of species. It links mathematically to a widely
cited global biodiversity indicator, the Living Planet Index, for which we
propose an improved formula that directly addresses the known problem of
singularities caused by extinctions. We show that, in an ideal market, trade in
our metric would lead to near-optimal allocation of resources to species
conservation. We further show that the metric is closely related to several
other metrics and indices already in use. Barriers to adoption are therefore
low. Used in conjunction with metrics addressing ecosystem functioning and
services, potential areas of application include biodiversity related financial
disclosures and voluntary or legislated no net biodiversity loss policies.


  Unfolded protein aggregation in cellular system is a problem causing various
types of diseases depending on which type unfolded proteins aggregate. This
phenomenon of aggregation may take place during production, storage, shipment
or delivery in the cellular medium. In the present work, we studied a
simplified and extended version of unfolded protein aggregation model by Lumry
and Eyring using stochastic approach. We solved analytically the Master
equation of the model for the probability distribution $P(x,t)$ of the unfolded
protein population and the solution was found to be time dependent complex
binomial distribution. In the large population limit $P(x,t)\sim
\Lambda(x,t)\times Pois(x,t)$. Further, the distribution became Normal
distribution at large population and mean of the distribution limit:
$P(x,t)\sim\Lambda(x,t)\times N(\langle qx\rangle,\langle qx\rangle)$. The
fluctuations inherent in the dynamics measured by Fano factor can have
sub-Poisson, Poisson and super-Poisson at different situations.


  We consider a general class of Markovian models describing the growth in a
randomly fluctuating environment of a clonal biological population having
several phenotypes related by stochastic switching. Phenotypes differ e.g. by
the level of gene expression for a population of bacteria. The time-averaged
growth rate of the population, $\Lambda$, is self-averaging in the limit of
infinite times; it may be understood as the fitness of the population in a
context of Darwinian evolution.
  The observation time $T$ being however typically finite, the growth rate
fluctuates. For $T$ finite but large, we obtain the variance of the
time-averaged growth rate as the maximum of a functional based on the
stationary probability distribution for the phenotypes. This formula is
general. In the case of two states, the stationary probability was computed by
Hufton, Lin and Galla \cite{HufLin2}, allowing for an explicit expression which
can be checked numerically.


  Stability is a key ingredient of protein fitness and its modification through
targeted mutations has applications in various fields such as protein
engineering, drug design and deleterious variant interpretation. Many studies
have been devoted over the past decades to building new, more effective methods
for predicting the impact of mutations on protein stability, based on the
latest developments in artificial intelligence (AI). We discuss their features,
algorithms, computational efficiency, and accuracy estimated on an independent
test set. We focus on a critical analysis of their limitations, the recurrent
biases towards the training set, their generalizability and interpretability.
We found that the accuracy of the predictors has stagnated at around 1 kcal/mol
for over 15 years. We conclude by discussing the challenges that need to be
addressed to reach improved performance.


  A completely non-generous and reciprocal population of players can create a
robust cooperating state that cannot be invaded by always defecting free riders
if the interactions among players are repeated for long enough. However, strict
non-generosity and strict reciprocity are ideal concepts, and may not even be
desirable sometimes. Therefore, to what extent generosity or non-reciprocity
can be allowed while still not be swamped by the mutants, is a natural
question. In this paper, we not only ask this question but furthermore ask how
generosity comparatively fares against non-reciprocity in this context. For
mathematical concreteness, we work within the framework of multiplayer repeated
prisoner's dilemma game with reactive strategies in a finite and an infinite
population; and explore the aforementioned questions through the effects of the
benefit to cost ratio, the interaction group size, and the population size.


  Phylogenetic (i.e. leaf-labeled) trees play a fundamental role in
evolutionary research. A typical problem is to reconstruct such trees from data
like DNA alignments (whose columns are often referred to as characters), and a
simple optimization criterion for such reconstructions is maximum parsimony. It
is generally assumed that this criterion works well for data in which state
changes are rare. In the present manuscript, we prove that each phylogenetic
tree $T$ with $n\geq 20 k$ leaves is uniquely defined by the set $A_k(T)$,
which consists of all characters with parsimony score $k$ on $T$. This can be
considered as a promising first step towards showing that maximum parsimony as
a tree reconstruction criterion is justified when the number of changes in the
data is relatively small.


  According to the World Health Organization the global incidence rate of
dengue infections have risen drastically in recent years. It is estimated that
globally the number of new infections is in the range of $100$ to $400$ million
per annum. At the same time a number of recent studies reported the existence
of Wolbachia strains, which inhibit dengue virus replication in mosquito
species that are primary vectors for dengue transmission. In this study we
focus on the impact of Wolbachia on dengue transmission dynamics in an
age-structured human population. We introduce a mathematical model, which takes
into account age-related effects, such as age-dependent human recovery and
mortality rates, as well as age-dependent vector to human dengue transmission
efficacy. We deduce an explicit formula for the basic reproduction number
$\mathcal{R}_0$, which allows us to assess the impact of the above mentioned
age-related effects on the local asymptotic stability of the dengue free
equilibrium. The formula we deduce for $\mathcal{R}_0$ also shows the intricate
relationship between human demography and the presence of a dengue inhibiting
Wolbachia strain.


  Genetic circuit design is a well-studied problem in synthetic biology. Ever
since the first genetic circuits -- the repressilator and the toggle switch --
were designed and implemented, many advances have been made in this area of
research. The current review systematically organizes a number of key works in
this domain by employing the versatile framework of generalized morphological
analysis. Literature in the area has been mapped based on (a) the design
methodologies used, ranging from brute-force searches to control-theoretic
approaches, (b) the modelling techniques employed, (c) various circuit
functionalities implemented, (d) key design characteristics, and (e) the
strategies used for the robust design of genetic circuits. We conclude our
review with an outlook on multiple exciting areas for future research, based on
the systematic assessment of key research gaps that have been readily
unravelled by our analysis framework.


  Tandem mass spectra capture fragmentation patterns that provide key
structural information about a molecule. Although mass spectrometry is applied
in many areas, the vast majority of small molecules lack experimental reference
spectra. For over seventy years, spectrum prediction has remained a key
challenge in the field. Existing deep learning methods do not leverage global
structure in the molecule, potentially resulting in difficulties when
generalizing to new data. In this work we propose a new model, MassFormer, for
accurately predicting tandem mass spectra. MassFormer uses a graph transformer
architecture to model long-distance relationships between atoms in the
molecule. The transformer module is initialized with parameters obtained
through a chemical pre-training task, then fine-tuned on spectral data.
MassFormer outperforms competing approaches for spectrum prediction on multiple
datasets, and is able to recover prior knowledge about the effect of collision
energy on the spectrum. By employing gradient-based attribution methods, we
demonstrate that the model can identify relationships between fragment peaks.
To further highlight MassFormer's utility, we show that it can match or exceed
existing prediction-based methods on two spectrum identification tasks. We
provide open-source implementations of our model and baseline approaches, with
the goal of encouraging future research in this area.


  Mitochondrial and nuclear genomes must be co-adapted to ensure proper
cellular respiration and energy production. Mito-nuclear incompatibility
reduces individual fitness and induces hybrid infertility, suggesting a
possible role in reproductive barriers and speciation. Here we develop a
birth-death model for evolution in spatially extended populations under
selection for mito-nuclear co-adaptation. Mating is constrained by physical and
genetic proximity, and offspring inherit nuclear genomes from both parents,
with recombination. The model predicts macroscopic patterns including a
community's long-term species diversity, its species abundance distribution,
speciation and extinction rates, as well as intra- and inter-specific genetic
variation. We explore how these long-term outcomes depend upon the microscopic
parameters of reproduction: individual fitness governed by mito-nuclear
compatibility, constraints on mating compatibility, and ecological carrying
capacity. We find that strong selection for mito-nuclear compatibility reduces
the equilibrium number of species after a radiation, increases the species'
abundances, while simultaneously increasing both speciation and extinction
rates. The negative correlation between species diversity and diversification
rates in our model agrees with the broad empirical pattern of lower species
diversity and higher speciation/extinction rates in temperate regions, compared
to the tropics. We therefore suggest that these empirical patterns may be
caused in part by latitudinal variation in metabolic demands, and corresponding
variation in selection on mito-nuclear function.


  Antipredator behavior is present in many biological systems where individuals
collectively react to an imminent attack. The antipredator response may
influence spatial pattern formation and ecosystem stability but requires an
organism's cost to contribute to the collective effort. We investigate a
nonhierarchical tritrophic system, whose predator-prey interactions are
described by the rock-paper-scissors game rules. In our spatial stochastic
simulations, the radius of antipredator response defines the maximum prey group
size that disturbs the predator's action, determining the individual cost to
participate in antipredator strategies. We consider that each organism
contributes equally to the collective effort, having its mobility limited by
the proportion of energy devoted to the antipredator reaction. Our outcomes
show that the antipredator response leads to spiral patterns, with the
segregation of organisms of the same species occupying departed spatial
domains. We found that a less localized antipredator response increases the
average size of the single-species patches, improving the protection of
individuals against predation. Finally, our findings show that although the
increase of the predation risk for a more localized antipredator response, the
high mobility constraining benefits species coexistence. Our results may help
ecologists understand the mechanisms leading to the stability of biological
systems where locality is crucial to behavioral interactions among species.


  Data on high-dimensional spheres arise frequently in many disciplines either
naturally or as a consequence of preliminary processing and can have intricate
dependence structure that needs to be understood. We develop exploratory factor
analysis of the projected normal distribution to explain the variability in
such data using a few easily interpreted latent factors. Our methodology
provides maximum likelihood estimates through a novel fast alternating
expectation profile conditional maximization algorithm. Results on simulation
experiments on a wide range of settings are uniformly excellent. Our
methodology provides interpretable and insightful results when applied to
tweets with the $\#MeToo$ hashtag in early December 2018, to time-course
functional Magnetic Resonance Images of the average pre-teen brain at rest, to
characterize handwritten digits, and to gene expression data from cancerous
cells in the Cancer Genome Atlas.


  The composition of cellular metabolism is different across species. Empirical
data reveal that bacterial species contain similar numbers of metabolic
reactions but that the cross-species popularity of reactions is so heterogenous
that some reactions are found in all the species while others are in just few
species, characterized by a power-law distribution with the exponent one.
Introducing an evolutionary model concretizing the stochastic recruitment of
chemical reactions into the metabolism of different species at different times
and their inheritance to descendants, we demonstrate that the exponential
growth of the number of species containing a reaction and the saturated
recruitment rate of brand-new reactions lead to the empirically identified
power-law popularity distribution. Furthermore, the structural characteristics
of metabolic networks and the species' phylogeny in our simulations agree well
with empirical observations.


  We introduce a model of parasite infection in a cell population, where cells
can be infected, either at birth through maternal transmission, from a contact
with the parasites reservoir, or because of the parasites released in the cell
medium after the lyses of infected cells. Inside the cells and between
infection events, the quantity of parasites evolves as a general non linear
branching process. We study the long time behaviour of the infection.


  We predict the future course of ongoing susceptible-infected-susceptible
(SIS) epidemics on regular, Erd\H{o}s-R\'{e}nyi and Barab\'asi-Albert networks.
It is known that the contact network influences the spread of an epidemic
within a population. Therefore, observations of an epidemic, in this case at
the population-level, contain information about the underlying network. This
information, in turn, is useful for predicting the future course of an ongoing
epidemic. To exploit this in a prediction framework, the exact high-dimensional
stochastic model of an SIS epidemic on a network is approximated by a
lower-dimensional surrogate model. The surrogate model is based on a
birth-and-death process; the effect of the underlying network is described by a
parametric model for the birth rates. We demonstrate empirically that the
surrogate model captures the intrinsic stochasticity of the epidemic once it
reaches a point from which it will not die out. Bayesian parameter inference
allows for uncertainty about the model parameters and the class of the
underlying network to be incorporated directly into probabilistic predictions.
An evaluation of a number of scenarios shows that in most cases the resulting
prediction intervals adequately quantify the prediction uncertainty. As long as
the population-level data is available over a long-enough period, even if not
sampled frequently, the model leads to excellent predictions where the
underlying network is correctly identified and prediction uncertainty mainly
reflects the intrinsic stochasticity of the spreading epidemic. For predictions
inferred from shorter observational periods, uncertainty about parameters and
network class dominate prediction uncertainty. The proposed method relies on
minimal data and is numerically efficient, which makes it attractive either as
a standalone inference and prediction scheme or in conjunction with other
methods.


  In the recent paper by Mu et al (J. Phys. Chem. Lett. 2021, 12, 2900) the
maximal cooperativity of human glucokinase is explained by the "kinetic
resonance" effect derived from a minimal three-state model. However, a closer
inspection of the latter shows that this effect seems to be rather a particular
quantitative coincidence than a general phenomenon that reflects a physical
mechanism of conformational regulation.


  We study the impact of contact heterogeneity on epidemic dynamics. A system
characterized by multiple susceptible populations is considered. The
description of the spread of an infectious disease is obtained through the
study of a system of Boltzmann-type equations for the number densities of
social contacts of the introduced compartments. A macroscopic system of
equations characterizing observable effects of the epidemic is then derived to
assess the impact of contact heterogeneity.


  Motivated by the increasing number of COVID-19 cases that have been observed
in many countries after the vaccination and relaxation of non-pharmaceutical
interventions, we propose a mathematical model on time-varying networks for the
spread of recurrent epidemic diseases in a partially vaccinated population. The
model encapsulates several realistic features, such as the different
effectiveness of the vaccine against transmission and development of severe
symptoms, testing practices, the possible implementation of non-pharmaceutical
interventions to reduce the transmission, isolation of detected individuals,
and human behaviour. Using a mean-field approach, we analytically derive the
epidemic threshold of the model and, if the system is above such a threshold,
we compute the epidemic prevalence at the endemic equilibrium. These
theoretical results show that precautious human behaviour and effective testing
practices are key toward avoiding epidemic outbreaks. Interestingly, we found
that, in many realistic scenarios, vaccination is successful in mitigating the
outbreak by reducing the prevalence of seriously ill patients, but it could be
a double-edged sword, whereby in some cases it might favour resurgent
outbreaks, calling for higher testing rates, more cautiousness and
responsibility among the population, or the reintroduction of
non-pharmaceutical interventions to achieve complete eradication.


  Mathematical models are instrumental to forecast the spread of pathogens and
to evaluate the effectiveness of non-pharmaceutical measures. A plethora of
optimal strategies has been recently developed to minimize either the infected
peak prevalence (IPP) or the epidemic final size (EFS). While most of the
control strategies optimize a simple cost function along a fixed finite-time
horizon, no consensus has been reached about how to simultaneously handle the
IPP, the EFS, and the avoiding of new cycles of infections rebounding. In this
work, based on the characterization of the dynamical behaviour of SIR-type
models under control actions (including the stability of equilibrium sets, in
terms of the herd immunity), it is studied how to minimize the EFS while
keeping - at any time - the IPP controlled. A procedure is proposed to tailor
non-pharmaceutical interventions by separating transient from stationary
control objectives and the potential benefits of the strategy are illustrated
by a detailed analysis and simulation results related to the COVID-19 pandemic.


  This article explores mathematical models for understanding the evolution of
contagious diseases. The most widely known set of models are the compartmental
ones, which are based on a set of differential equations. But these are not the
only models. This review visits many different families of models.
Additionally, we show these families, not as unrelated entities, but following
a common thread in which the problems or assumptions of a model are solved or
generalized by another model. In this way, we can understand their
relationships, assumptions, simplifications, and, ultimately, limitations.
Prompted by the current Covid19 pandemic, we have a special focus on spread
forecasting. We illustrate the difficulties encountered to do realistic
predictions. In general, they are only approximations to a reality whose
biological and societal complexity is much larger. Particularly troublesome are
the large underlying variability, the problem's time-varying nature, and the
difficulty to estimate the required parameters for a faithful model.
Additionally, we will also see that these models have a multiplicative nature
implying that small errors in the system parameters cause a huge uncertainty in
the prediction. Stochastic or agent-based models can overcome some of the
modeling problems of systems based on differential or stochastic equations.
Their main difficulty is that they are as accurate and realistic as the data
available to estimate their detailed parametrization, and very often this
detailed data is not at the modeller's disposal. Although the predictive power
of mathematical models to forecast the evolution of a contagious disease is
very limited, these models are still very useful to plan interventions as they
can calculate their impact if all other parameters stay fixed. They are also
very useful to understand the properties of disease propagation in complex
systems.


  Through the use of cutting-edge unsupervised classification techniques from
statistics and machine learning, we characterise symptom phenotypes among
symptomatic SARS-CoV-2 PCR-positive community cases. We first analyse each
dataset in isolation and across age bands, before using methods that allow us
to compare multiple datasets. While we observe separation due to the total
number of symptoms experienced by cases, we also see a separation of symptoms
into gastrointestinal, respiratory and other types, and different symptom
co-occurrence patterns at the extremes of age. In this way, we are able to
demonstrate the deep structure of symptoms of COVID-19 without usual biases due
to study design. This is expected to have implications for the identification
and management of community SARS-CoV-2 cases and could be further applied to
symptom-based management of other diseases and syndromes.


  May and Leonard (SIAM J. Appl. Math 1975) introduced a three-species
Lotka-Volterra type population model that exhibits heteroclinic cycling. Rather
than producing a periodic limit cycle, the trajectory takes longer and longer
to complete each "cycle", passing closer and closer to unstable fixed points in
which one population dominates and the others approach zero. Aperiodic
heteroclinic dynamics have subsequently been studied in ecological systems
(side-blotched lizards; colicinogenic E. coli), in the immune system, in neural
information processing models ("winnerless competition"), and in models of
neural central pattern generators. Yet as May and Leonard observed
"Biologically, the behavior (produced by the model) is nonsense. Once it is
conceded that the variables represent animals, and therefore cannot fall below
unity, it is clear that the system will, after a few cycles, converge on some
single population, extinguishing the other two." Here, we explore different
ways of introducing discrete stochastic dynamics based on May and Leonard's ODE
model, with application to ecological population dynamics, and to a neuromotor
central pattern generator system. We study examples of several quantitatively
distinct asymptotic behaviors, including total extinction of all species,
extinction to a single species, and persistent cyclic dominance with finite
mean cycle length.


  Widespread application of insecticide remains the primary form of control for
Chagas disease in Central America, despite only temporarily reducing domestic
levels of the endemic vector Triatoma dimidiata and having little long-term
impact. Recently, an approach emphasizing community feedback and housing
improvements has been shown to yield lasting results. However, the additional
resources and personnel required by such an intervention likely hinders its
widespread adoption. One solution to this problem would be to target only a
subset of houses in a community while still eliminating enough infestations to
interrupt disease transfer. Here we develop a sequential sampling framework
that adapts to information specific to a community as more houses are visited,
thereby allowing us to efficiently find homes with domiciliary vectors while
minimizing sampling bias. The method fits Bayesian geostatistical models to
make spatially informed predictions, while gradually transitioning from
prioritizing houses based on prediction uncertainty to targeting houses with a
high risk of infestation. A key feature of the method is the use of a single
exploration parameter, $\alpha$, to control the rate of transition between
these two design targets. In a simulation study using empirical data from five
villages in southeastern Guatemala, we test our method using a range of values
for $\alpha$, and find it can consistently select fewer homes than random
sampling, while still bringing the village infestation rate below a given
threshold. We further find that when additional socioeconomic information is
available, much larger savings are possible, but that meeting the target
infestation rate is less consistent, particularly among the less exploratory
strategies. Our results suggest new options for implementing long-term T.
dimidiata control.


  We introduce and investigate the space of \emph{equidistant}
$X$-\emph{cactuses}. These are rooted, arc weighted, phylogenetic networks with
leaf set $X$, where $X$ is a finite set of species, and all leaves have the
same distance from the root. The space contains as a subset the space of
ultrametric trees on $X$ that was introduced by Gavryushkin and Drummond. We
show that equidistant-cactus space is a CAT(0)-metric space which implies, for
example, that there are unique geodesic paths between points. As a key step to
proving this, we present a combinatorial result concerning \emph{ranked} rooted
$X$-cactuses. In particular, we show that such networks can be encoded in terms
of a pairwise compatibility condition arising from a poset of collections of
pairs of subsets of $X$ that satisfy certain set-theoretic properties. As a
corollary, we also obtain an encoding of ranked, rooted $X$-trees in terms of
partitions of $X$, which provides an alternative proof that the space of
ultrametric trees on $X$ is CAT(0). As with spaces of phylogenetic trees, we
expect that our results should provide the basis for and new directions in
performing statistical analyses for collections of phylogenetic networks with
arc lengths.


  In high-throughput functional genomic screens, each gene product is commonly
assumed to exhibit a singular biological function within a defined protein
complex or pathway. In practice, a single gene perturbation may induce multiple
cascading functional outcomes, a genetic principle known as pleiotropy. Here,
we model pleiotropy in fitness screen collections by representing each gene
perturbation as the sum of multiple perturbations of biological functions, each
harboring independent fitness effects inferred empirically from the data. Our
approach ('Webster') recovered pleiotropic functions for DNA damage proteins
from genotoxic fitness screens, untangled distinct signaling pathways upstream
of shared effector proteins from cancer cell fitness screens, and learned
aspects of the cellular hierarchy in an unsupervised manner. Modeling compound
sensitivity profiles in terms of genetically defined functions recovered
compound mechanisms of action. Our approach establishes a sparse approximation
mechanism for unraveling complex genetic architectures underlying
high-dimensional gene perturbation readouts.


  The metabolic activity of microbes has played an essential role in the
evolution and persistence of life on Earth. Microbial metabolism plays a
primary role in the flow of carbon, nitrogen and other elements through the
biosphere on a global scale. Microbes perform these metabolic activities in the
context of complex communities comprised of many species that interact in
dynamic and spatially-structured environments. Molecular genetics has revealed
many of the metabolic pathways microbes utilize to generate energy and biomass.
However, most of this knowledge is derived from model organisms, so we have a
limited view of role of the massive genomic diversity in the wild on metabolic
phenotypes. Further, we are only beginning to glimpse the principles governing
how the metabolism of a community emerges from the collective action of its
constituent members. As a result, one of the biggest challenges in the field is
to understand how the metabolic activity of a community emerges from the
genomic structure of the constituents. Here we propose an approach to this
problem that rests on the quantitative analysis of metabolic activity in
ensembles of microbial communities. We propose quantifying metabolic fluxes in
diverse communities, either in the laboratory or the wild. We suggest that
using sequencing data to quantify the genomic, taxonomic or transcriptional
variation across an ensemble of communities can reveal low-dimensional
descriptions of community structure that can explain or predict their emergent
metabolic activity. We survey the types of communities for which this approach
might be best suited and then review the analytical techniques available for
quantifying metabolite dynamics in communities. Finally, we discuss what types
of data analysis approaches might be lucrative for learning the
structure-function mapping in communities from these data.


  AlphaFold predicts protein structures from the amino acid sequence at or near
experimental resolution, solving the 50-year-old protein folding challenge,
leading to progress by transforming large-scale genomics data into protein
structures. AlphaFold will also greatly change the scientific research model
from low-throughput to high-throughput manner. The AlphaFold framework is a
mixture of two types of workloads: MSA construction based on CPUs and model
inference on GPUs. The first CPU stage dominates the overall runtime, taking
hours for a single protein due to the large database sizes and I/O bottlenecks.
However, GPUs in this CPU stage remain idle, resulting in low GPU utilization
and restricting the capacity of large-scale structure predictions. Therefore,
we proposed ParaFold, an open-source parallel version of AlphaFold for high
throughput protein structure predictions. ParaFold separates the CPU and GPU
parts to enable large-scale structure predictions. ParaFold also effectively
reduces the CPU and GPU runtime with two optimizations without compromising the
quality of prediction results: using multi-threaded parallelism on CPUs and
using optimized JAX compilation on GPUs. We evaluated ParaFold with three
datasets of different size and protein lengths. We evaluated the accuracy and
efficiency of optimizations on CPUs and GPUs, and showed the large-scale
prediction capability by running ParaFold inferences of 19,704 small proteins
in five hours on one NVIDIA DGX-2. Using the JAX compile optimization, ParaFold
attained a 13.8X average speedup over AlphaFold. ParaFold offers a rapid and
effective approach for high-throughput structure predictions, leveraging the
predictive power by running on supercomputers, with shorter time, and at a
lower cost. The development of ParaFold will greatly speed up high-throughput
studies and render the protein "structure-omics" feasible.


  The bacterium Group A Streptococcus (Streptococcus pyogenes, GAS) is a
human-specific pathogen and a major cause of global morbidity and mortality.
Despite decades of research our knowledge of GAS infection and immunity is
incomplete, hampering vaccine design and other efforts to reduce disease
prevalence. Epidemiological studies indicate positive associations between the
prevalence of GAS-related disease, the diversity of circulating strains and the
degree of poverty in host populations. However, the infection and immune
mechanisms underlying these associations are not clear. In this work, we use an
agent-based model to demonstrate that observed diversity and prevalence are
best accounted for by the hypothesis that GAS infection confers enduring
strain-specific immunity, with reduced or absent cross-protection against
infection by other strains. Our results suggest that the success of GAS
vaccines will depend on their ability to elicit long-lasting cross-protective
immunity over multiple strain types.


  Computational protein design, i.e. inferring novel and diverse protein
sequences consistent with a given structure, remains a major unsolved
challenge. Recently, deep generative models that learn from sequences alone or
from sequences and structures jointly have shown impressive performance on this
task. However, those models appear limited in terms of modeling structural
constraints, capturing enough sequence diversity, or both. Here we consider
three recently proposed deep generative frameworks for protein design: (AR) the
sequence-based autoregressive generative model, (GVP) the precise
structure-based graph neural network, and Fold2Seq that leverages a fuzzy and
scale-free representation of a three-dimensional fold, while enforcing
structure-to-sequence (and vice versa) consistency. We benchmark these models
on the task of computational design of antibody sequences, which demand
designing sequences with high diversity for functional implication. The
Fold2Seq framework outperforms the two other baselines in terms of diversity of
the designed sequences, while maintaining the typical fold.


  In the context of population genetics, active information can be extended to
measure the change of information of a given event (e.g., fixation of an
allele) from a neutral model in which only genetic drift is taken into account
to a non-neutral model that includes other sources of frequency variation
(e.g., selection and mutation). In this paper we illustrate active information
in population genetics through the Wright-Fisher model.


  Due to the relevance for conservation biology, there is an increasing
interest to extend evolutionary genomics models to plant, animal or microbial
species. However, this requires to understand the effect of life-history traits
absent in humans on genomic evolution. In this context, it is fundamentally of
interest to generalize the replicator equation, which is at the heart of most
population genomics models. However, as the inclusion of life-history traits
generates models with a large state space, the analysis becomes involving. We
focus, here, on quiescence and seed banks, two features common to many plant,
invertebrate and microbial species. We develop a method to obtain a
low-dimensional replicator equation in the context of evolutionary game theory,
based on two assumptions: (1) the life-history traits are {\it per se} neutral,
and (2) frequency-dependent selection is weak. We use the results to
investigate the evolution and maintenance of cooperation based on the
Prisoner's dilemma. We first consider the generalized replicator equation, and
then refine the investigation using adaptive dynamics. It turns out that,
depending on the structure and timing of the quiescence/dormancy life-history
trait, cooperation in a homogeneous population can be stabilized. We finally
discuss and highlight the relevance of these results for plant, invertebrate
and microbial communities.


  Multivalent cargo that can interact with substrates via multiple interaction
sites exhibit shared characteristics despite being found in different systems
at different length-scales. Here, a general analytical model has been developed
that can describe the motion of multivalent cargo as a response to position
dependence in the binding and unbinding rates of their interaction sites. Cargo
exhibit both an effective diffusivity and velocity, which acts in the direction
of increasing cargo-substrate binding rate and decreasing cargo-substrate
unbinding rate. This model can reproduce previously published experimental
findings using only the binding and unbinding rate distributions of cargo
interaction sites, and without any further parameter fitting. Extension of the
cargo binding model to two dimensions reveals an effective velocity with the
same properties as that derived for the $1$D case.


  Motivated by chemical reaction rules, we introduce a rule-based
epidemiological framework for the mathematical modelling of the COVID-19
pandemic which lies the foundations for further study of this class of models.
Here we stress that we do not have a specific model in mind, but a whole
collection of models which can be transformed into each other, or represent
different aspects of the pandemic. Each single model is based on a mathematical
formulation that we call a rule-based system. They have several advantages, for
example that they can be both interpreted both stochastically and
deterministically without changing the model structure. Rule-based systems
should be easier to communicate to non-specialists, when compared to
differential equations. Due to their combinatorial nature, the rule-based model
framework we propose is ideal for systematic mathematical modelling, systematic
links to statistics, data analysis in general and machine learning. Most
importantly, the framework can be used inside a Bayesian model selection
algorithm, which can extend unsupervised machine learning technology to gain
mechanistic insight through mathematical models. We first introduce the
Bayesian model selection algorithm and give an introduction to rule-based
systems. We propose precursor models, all consistent with our rule-based
framework and relevant for the COVID-19 pandemic. We complement the models with
numerical results and, for some of them, we include analytical results such as
the calculation of the number of new infections generated by an infected person
during their infectious period. We conclude the article with an in depth
discussion of the role of mathematical modelling through the COVID-19 pandemic.


  Recent advances in molecular machine learning, especially deep neural
networks such as Graph Neural Networks (GNNs) for predicting structure activity
relationships (SAR) have shown tremendous potential in computer-aided drug
discovery. However, the applicability of such deep neural networks are limited
by the requirement of large amounts of training data. In order to cope with
limited training data for a target task, transfer learning for SAR modeling has
been recently adopted to leverage information from data of related tasks. In
this work, in contrast to the popular parameter-based transfer learning such as
pretraining, we develop novel deep transfer learning methods TAc and TAc-fc to
leverage source domain data and transfer useful information to the target
domain. TAc learns to generate effective molecular features that can generalize
well from one domain to another, and increase the classification performance in
the target domain. Additionally, TAc-fc extends TAc by incorporating novel
components to selectively learn feature-wise and compound-wise transferability.
We used the bioassay screening data from PubChem, and identified 120 pairs of
bioassays such that the active compounds in each pair are more similar to each
other compared to its inactive compounds. Our experiments clearly demonstrate
that TAc achieves significant improvement over all baselines across a large
number of target tasks. Furthermore, although TAc-fc achieves slightly worse
ROC-AUC on average compared to TAc, TAc-fc still achieves the best performance
on more tasks in terms of PR-AUC and F1 compared to other methods. In summary,
TAc-fc is also found to be a strong model with competitive or even better
performance than TAc on a notable number of target tasks.


  Ancestral state reconstruction is one of the most important tasks in
evolutionary biology. Conditions under which we can reliably reconstruct the
ancestral state have been studied for both discrete and continuous traits.
However, the connection between these results is unclear, and it seems that
each model needs different conditions. In this work, we provide a unifying
theory on the consistency of ancestral state reconstruction for various types
of trait evolution models. Notably, we show that for a sequence of nested trees
with bounded heights, the necessary and sufficient conditions for the existence
of a consistent ancestral state reconstruction method under discrete models,
the Brownian motion model, and the threshold model are equivalent. When tree
heights are unbounded, we provide a simple counter-example to show that this
equivalence is no longer valid.


  We introduce a model to study the impact of catastrophes on evolutionary
paths. If we do not allow catastrophes the number of changes in the maximum
fitness of a population grows logarithmically with respect to time. Allowing
catastrophes (no matter how rare) yields a drastically different behavior. When
catastrophes are possible the number of changes in the maximum fitness of the
population grows linearly with time. Moreover, the evolutionary paths are a lot
less predictable when catastrophes are possible. Our results can be seen as
supporting the hypothesis that catastrophes speed up evolution by disrupting
dominant species and creating space for new species to emerge and evolve.


  Multi-strain pathogens such as Group A Streptococcus, Streptococcus
pneumoniae, and Staphylococcus aureus cause millions of infections each year
with a substantial health burden. Control of multi-strain pathogens can be
complicated by the high strain diversity often observed in endemic settings. It
is not well understood how high strain diversity is maintained in populations,
given that they compete with each other both directly (within an individual
host) and indirectly (via host immunity). Previous modelling studies have
investigated how indirect competition affects the prevalence and diversity of
strains. However, these studies often make simplifying assumptions about the
direct competition that occurs within hosts. Currently, little data is
available to validate these assumptions, hence there is a need to clarify how
sensitive model outputs are to these assumptions. In this study, we compare the
dynamics of multi-strain pathogens under different assumptions about direct
competition between strains using an agent-based model. We find that the
assumptions made about direct competition can affect the epidemiological
dynamics, particularly when there is no long-term immunity following infections
and a low rate of importation of non-circulating strains. Our results suggest
that while direct and indirect competition can each decrease strain diversity
when they act in isolation, they may increase strain diversity when they act
together. This finding highlights the importance of examining sensitivity to
assumptions about strain competition. In particular, omitting consideration of
direct competition can lead to inaccurate estimates of the likely effectiveness
of control strategies as changes in strain diversity shift the level of direct
strain competition.


  A selfish learner seeks to maximize their own success, disregarding others.
When success is measured as payoff in a game played against another learner,
mutual selfishness typically fails to produce the optimal outcome for a pair of
individuals. However, learners often operate in populations, and each learner
may have a limited duration of interaction with any other individual. Here, we
compare selfish learning in stable pairs to selfish learning with stochastic
encounters in a population. We study gradient-based optimization in repeated
games like the prisoner's dilemma, which feature multiple Nash equilibria, many
of which are suboptimal. We find that myopic, selfish learning, when
distributed in a population via ephemeral encounters, can reverse the dynamics
that occur in stable pairs. In particular, when there is flexibility in partner
choice, selfish learning in large populations can produce optimal payoffs in
repeated social dilemmas. This result holds for the entire population, not just
for a small subset of individuals. Furthermore, as the population size grows,
the timescale to reach the optimal population payoff remains finite in the
number of learning steps per individual. While it is not universally true that
interacting with many partners in a population improves outcomes, this form of
collective learning achieves optimality for several important classes of social
dilemmas. We conclude that na\"{i}ve learning can be surprisingly effective in
populations of individuals navigating conflicts of interest.


  Advances in next-generation metagenome sequencing have the potential to
revolutionize the point-of-care diagnosis of novel pathogen infections, which
could help prevent potential widespread transmission of diseases. Given the
high volume of metagenome sequences, there is a need for scalable frameworks to
analyze and segment metagenome sequences from clinical samples, which can be
highly imbalanced. There is an increased need for learning robust
representations from metagenome reads since pathogens within a family can have
highly similar genome structures (some more than 90%) and hence enable the
segmentation and identification of novel pathogen sequences with limited
labeled data. In this work, we propose Metagenome2Vec - a contextualized
representation that captures the global structural properties inherent in
metagenome data and local contextualized properties through self-supervised
representation learning. We show that the learned representations can help
detect six (6) related pathogens from clinical samples with less than 100
labeled sequences. Extensive experiments on simulated and clinical metagenome
data show that the proposed representation encodes compositional properties
that can generalize beyond annotations to segment novel pathogens in an
unsupervised setting.


  Multitype branching processes are ideal for studying the population dynamics
of stem cell populations undergoing mutation accumulation over the years
following transplant. In such stochastic models, several quantities are of
clinical interest as insertional mutagenesis carries the potential threat of
leukemogenesis following gene therapy with autologous stem cell
transplantation. In this paper, we develop a three-type branching process model
describing accumulations of mutations in a population of stem cells
distinguished by their ability for long-term self-renewal. Our outcome of
interest is the appearance of a double-mutant cell, which carries a high
potential for leukemic transformation. In our model, a single-hit mutation
carries a slight proliferative advantage over a wild-type stem cells. We
compute marginalized transition probabilities that allow us to capture
important quantitative aspects of our model, including the probability of
observing a double-hit mutant and relevant moments of a single-hit mutation
population over time. We thoroughly explore the model behavior numerically,
varying birth rates across the initial sizes and populations of wild type stem
cells and single-hit mutants, and compare the probability of observing a
double-hit mutant under these conditions. We find that increasing the number of
single-mutants over wild-type particles initially present has a large effect on
the occurrence of a double-mutant, and that it is relatively safe for
single-mutants to be quite proliferative, provided the lentiviral gene addition
avoids creating single mutants in the original insertion process. Our approach
is broadly applicable to an important set of questions in cancer modeling and
other population processes involving multiple stages, compartments, or types.


  Viruses are right at the interface of inanimate matter and life. However,
recent experiments [T. Sakai, et al., J.~Virol.~{\bf 92}, e01522-17 (2018)]
have shown that some influenza strains can actively roll on glycan-covered
surfaces. In a previous letter [F. Ziebert and I. M. Kuli\'{c}, Phys. Rev.
Lett. {\bf 126}, 218101 (2021)] we suggested this to be a form of viral surface
metabolism: a collection of spike proteins that attach to and cut the glycans
act as a self-organized mechano-chemical motor. Here we study in more depth the
physics of the emergent self-rolling states. We give scaling arguments how the
motion arises, substantiated by a detailed analytical theory that yields the
full torque-angular velocity relation of the self-organized motor. Stochastic
Gillespie simulations are used to validate the theory and to quantify
stochastic effects like virus detachment and reversals of its direction.
Finally, we also cross-check several approximations made previously and show
that the proposed mechanism is very robust. All these results point together to
the statistical inevitability of viral rolling in presence of enzymatic
activity.


  The aim of this paper is to study the long-term consequence on longevity of a
mortality shock. We adopt an historical and modeling approach to study how the
population evolution following a mortality shock such as the COVID-19 pandemic
could impact future mortality rates. In the first of part the paper, we study
the several cholera epidemics in France and in England starting from the 1830s,
and their impact on the major development of public health at the end of the
nineteenth century. In the second part, we present the mathematical modeling of
stochastic Individual-Based models. Using the R package IBMPopSim, this
flexible framework is then applied to simulate the long-term impact of a
mortality shock, using a toy model where nonlinear population compositional
changes affect future mortality rates.


  Kuhn (1962) proposed an evolutionary model to explain how scientific
knowledge is built, based on the concept of paradigm. Even though Kuhn's model
is general, it has been applied only to a few topics in evolutionary biology,
especially broad-based paradigms. Our goal here is to analyze a small-scale
paradigm change that occurred about the mating system of a Mediterranean shrub:
P. angustifolia (Oleaceae) through the lens of Kuhn's model. We first summarize
the different steps of the paradigm change and replace them in the more general
context of the sex ratio theory. Second, we show how the different steps of the
paradigm changes can be interpreted by Kuhnian concepts and tools. Finally, we
discuss the actual and future status of the new paradigm.


  Despite its high and direct impact on nearly all biological processes, the
underlying structure of gene-gene interaction networks is investigated so far
according to pair connections. To address this, we explore the gene interaction
networks of the yeast Saccharomyces cerevisiae beyond pairwise interaction
using the structural balance theory (SBT). Specifically, we ask whether
essential and nonessential gene interaction networks are structurally balanced.
We study triadic interactions in the weighted signed undirected gene networks
and observe that balanced and unbalanced triads are over and underrepresented
in both networks, thus beautifully in line with the strong notion of balance.
Moreover, we note that the energy distribution of triads is significantly
different in both essential and nonessential networks compared with the
shuffled networks. Yet, this difference is greater in the essential network
regarding the frequency as well as the energy of triads. Additionally, results
demonstrate that triads in the essential gene network are more interconnected
through sharing common links, while in the nonessential network they tend to be
isolated. Last but not least, we investigate the contribution of all-length
signed walks and its impact on the degree of balance. Our findings reveal that
interestingly when considering longer cycles the nonessential gene network is
more balanced compared to the essential network.


  Minimizers and convolutional neural networks (CNNs) are two quite distinct
popular techniques that have both been employed to analyze categorical
biological sequences. At face value, the methods seem entirely dissimilar.
Minimizers use min-wise hashing on a rolling window to extract a single
important k-mer feature per window. CNNs start with a wide array of randomly
initialized convolutional filters, paired with a pooling operation, and then
multiple additional neural layers to learn both the filters themselves and how
they can be used to classify the sequence.
  Here, our main result is a careful mathematical analysis of hash function
properties showing that for sequences over a categorical alphabet, random
Gaussian initialization of convolutional filters with max-pooling is equivalent
to choosing a minimizer ordering such that selected k-mers are (in Hamming
distance) far from the k-mers within the sequence but close to other
minimizers. In empirical experiments, we find that this property manifests as
decreased density in repetitive regions, both in simulation and on real human
telomeres. We additionally train from scratch a CNN embedding of synthetic
short-reads from the SARS-CoV-2 genome into 3D Euclidean space that locally
recapitulates the linear sequence distance of the read origins, a modest step
towards building a deep learning assembler, though it is at present too slow to
be practical. In total, this manuscript provides a partial explanation for the
effectiveness of CNNs in categorical sequence analysis.


  Some dividing cells sense their shape by becoming polarized along their long
axis. Cell polarity is controlled in part by polarity proteins like Rho GTPases
cycling between active membrane-bound forms and inactive cytosolic forms,
modeled as a "wave-pinning" reaction-diffusion process. Does shape sensing
emerge from wave-pinning? We show that wave pinning senses the cell's long
axis. Simulating wave-pinning on a curved surface, we find that high-activity
domains migrate to peaks and troughs of the surface. For smooth surfaces, a
simple rule of minimizing the domain perimeter while keeping its area fixed
predicts the final position of the domain and its shape. However, when we
introduce roughness to our surfaces, shape sensing can be disrupted, and
high-activity domains can become localized to locations other than the global
peaks and valleys of the surface. On rough surfaces, the domains of the
wave-pinning model are more robust in finding the peaks and troughs than the
minimization rule, though both can become trapped in steady states away from
the peaks and valleys. We can control the robustness of shape sensing by
altering the Rho GTPase diffusivity and the domain size. We also find that the
shape sensing properties of cell polarity models can explain how domains
localize to curved regions of deformed cells. Our results help to understand
the factors that allow cells to sense their shape - and the limits that
membrane roughness can place on this process.