  In this paper, we address the modeling issues of cell movement and division
with a special focus on the phenomenon of volume exclusion in a lattice-based,
exact stochastic simulation framework. We propose a new exact method, called
Reduced Rate Method -- RRM, that is substantially quicker than the previously
used exclusion method, for large number of cells. In addition, we introduce
three novel reaction types: the contact-inhibited, the contact-promoted, and
the spontaneous reactions. To the best of our knowledge, these reaction types
have not been taken into account in lattice-based stochastic simulations of
cell cultures. These new types of events may be easily applied to complicated
systems, enabling the generation of biologically feasible stochastic cell
culture simulations. Furthermore, we show that the exclusion algorithm and our
RRM algorithm are mathematically equivalent in the sense that the next reaction
to be realized and the corresponding sojourn time both belong to the same
reaction and time distributions in the two approaches -- even with the newly
introduced reaction types.
  Exact, agent-based, stochastic methods of cell culture simulations seem to be
undervalued and are mostly used as benchmarking tools to validate deterministic
approximations of the corresponding stochastic models. Our proposed methods are
exact, they are easy to implement, have a high predictive value, and can be
conveniently extended with new features. Therefore, these approaches promise a
great potential.


  Quantification of heritability is a fundamental desideratum in genetics,
which allows an assessment of the contribution of additive genetic variation to
the variability of a trait of interest. The traditional computational
approaches for assessing the heritability of a trait have been developed in the
field of quantitative genetics. However, the rise of modern population genomics
with large sample sizes has led to the development of several new machine
learning based approaches to inferring heritability. In this paper, we
systematically summarize recent advances in machine learning which can be used
to infer heritability. We focus on an application of these methods to bacterial
genomes, where heritability plays a key role in understanding phenotypes such
as antibiotic resistance and virulence, which are particularly important due to
the rising frequency of antimicrobial resistance. By designing a heritability
model incorporating realistic patterns of genome-wide linkage disequilibrium
for a frequently recombining bacterial pathogen, we test the performance of a
wide spectrum of different inference methods, including also GCTA. In addition
to the synthetic data benchmark, we present a comparison of the methods for
antibiotic resistance traits for multiple bacterial pathogens. Insights from
the benchmarking and real data analyses indicate a highly variable performance
of the different methods and suggest that heritability inference would likely
benefit from tailoring of the methods to the specific genetic architecture of
the target organism.


  The Yule branching process is a classical model for the random generation of
gene tree topologies in population genetics. It generates binary ranked trees
-- also called "histories" -- with a finite number $n$ of leaves. We study the
lengths $\ell_1 > \ell_2 > ... > \ell_k > ...$ of the external branches of a
Yule generated random history of size $n$, where the length of an external
branch is defined as the rank of its parent node. When $n \rightarrow \infty$,
we show that the random variable $\ell_k$, once rescaled as
$\frac{n-\ell_k}{\sqrt{n/2}}$, follows a $\chi$-distribution with $2k$ degrees
of freedom, with mean $\mathbb E(\ell_k) \sim n$ and variance $\mathbb
V(\ell_k) \sim n \big(k-\frac{\pi k^2}{16^k} \binom{2k}{k}^2\big)$. Our results
contribute to the study of the combinatorial features of Yule generated gene
trees, in which external branches are associated with singleton mutations
affecting individual gene copies.


  Molecular property prediction is of crucial importance in many disciplines
such as drug discovery, molecular biology, or material and process design. The
frequently employed quantitative structure-property/activity relationships
(QSPRs/QSARs) characterize molecules by descriptors which are then mapped to
the properties of interest via a linear or nonlinear model. In contrast, graph
neural networks, a novel machine learning method, directly work on the
molecular graph, i.e., a graph representation where atoms correspond to nodes
and bonds correspond to edges. GNNs allow to learn properties in an end-to-end
fashion, thereby avoiding the need for informative descriptors as in
QSPRs/QSARs. GNNs have been shown to achieve state-of-the-art prediction
performance on various property predictions tasks and represent an active field
of research. We describe the fundamentals of GNNs and demonstrate the
application of GNNs via two examples for molecular property prediction.


  Membrane proteins typically deform the surrounding lipid bilayer membrane,
which can play an important role in the function, regulation, and organization
of membrane proteins. Membrane elasticity theory provides a beautiful
description of protein-induced lipid bilayer deformations, in which all
physical parameters can be directly determined from experiments. Analytic
treatments of the membrane elasticity theory of protein-induced lipid bilayer
deformations have largely focused on idealized protein shapes with circular
cross section, and on perturbative solutions for proteins with non-circular
cross section. We develop here a boundary value method (BVM) that permits the
construction of non-perturbative analytic solutions of protein-induced lipid
bilayer deformations for non-circular protein cross sections, for constant as
well as variable boundary conditions along the bilayer-protein interface. We
apply this BVM to protein-induced lipid bilayer thickness deformations. Our BVM
reproduces available analytic solutions for proteins with circular cross
section and yields, for proteins with non-circular cross section, excellent
agreement with numerical, finite element solutions. On this basis, we formulate
a simple analytic approximation of the bilayer thickness deformation energy
associated with general protein shapes and show that, for modest deviations
from rotational symmetry, this analytic approximation is in good agreement with
BVM solutions. Using the BVM, we survey the dependence of protein-induced lipid
bilayer thickness deformations on protein shape, and thus explore how the
coupling of protein shape and bilayer thickness deformations affects protein
oligomerization and transitions in protein conformational state.


  Rapid urbanization is a major cause of habitat and biodiversity loss and
human-animal conflict. While urbanization is inevitable, we need to develop a
good understanding of the urban ecosystem and the urban-adapted species in
order to ensure sustainable cities for our future. Scavengers play a major role
in urban ecosystems, and often, urban adaptation involves a shift towards
scavenging behaviour in wild animals. We carried out an experiment at different
sites in the state of West Bengal, India, to identify the scavenging guild
within urban habitats, in response to human provided food. Our study revealed a
total of 17 different vertebrate species were identified across sites over 498
sessions of observations. We carried out network analysis to understand the
dynamics of the system, and found that the free-ranging dog and common mynah
were key species within the scavenging networks. This study revealed the
complexity of scavenging networks within human-dominated habitats.


  Progress in molecular methods has enabled the monitoring of bacterial
populations in time. Nevertheless, understanding community dynamics and its
links with ecosystem functioning remains challenging due to the tremendous
diversity of microorganisms. Conceptual frameworks that make sense of
time-series of taxonomically-rich bacterial communities, regarding their
potential ecological function, are needed. A key concept for organizing
ecological functions is the niche, the set of strategies that enable a
population to persist and define its impacts on the surroundings. Here we
present a framework based on manifold learning, to organize genomic information
into potentially occupied bacterial metabolic niches over time. We apply the
method to re-construct the dynamics of putatively occupied metabolic niches
using a long-term bacterial time-series from the Baltic Sea, the Linnaeus
Microbial Observatory (LMO). The results reveal a relatively low-dimensional
space of occupied metabolic niches comprising groups of taxa with similar
functional capabilities. Time patterns of occupied niches were strongly driven
by seasonality. Some metabolic niches were dominated by one bacterial taxon
whereas others were occupied by multiple taxa, and this depended on season.
These results illustrate the power of manifold learning approaches to advance
our understanding of the links between community composition and functioning in
microbial systems.


  Objectives The review is dedicated to evaluate the current monkeypox virus
(MPXV) detection methods, discuss their pros and cons, and provide recommended
solutions to the problems.
  Methods The literature for this review is identified through searches in
PubMed, Web of Science, Google Scholar, ResearchGate, and Science Direct
advanced search for articles published in English without any start date until
June, 2022, by use of the terms "monkeypox virus" or "poxvirus" along with
"diagnosis"; "PCR"; "real-time PCR"; "LAMP"; "RPA"; "immunoassay";
"reemergence"; "biothreat"; "endemic", and "multi-country outbreak" and also,
by tracking citations of the relevant papers. The most relevant articles are
included in the review.
  Results Our literature review shows that PCR is the gold standard method for
MPXV detection. In addition, loop-mediated isothermal amplification (LAMP) and
recombinase polymerase amplification (RPA) have been reported as alternatives
to PCR. Immunodiagnostics, whole particle detection, and image-based detection
are the non-nucleic acid-based MPXV detection modalities.
  Conclusions PCR is easy to leverage and adapt for a quick response to an
outbreak, but the PCR-based MPXV detection approaches may not be suitable for
marginalized settings. Limited progress has been made towards innovations in
MPXV diagnostics, providing room for the development of novel detection
techniques for this virus.


  The cell cycle is one of the most fundamental biological processes important
for understanding normal physiology and various pathologies such as cancer.
Single cell RNA sequencing technologies give an opportunity to analyse the cell
cycle transcriptome dynamics in an unprecedented range of conditions (cell
types and perturbations), with thousands of publicly available datasets. Here
we review the main computational tasks in such analysis: 1) identification of
cell cycle phases, 2) pseudotime inference, 3) identification and profiling of
cell cycle-related genes, 4) removing cell cycle effect, 5) identification and
analysis of the G0 (quiescent) cells. We review seventeen software packages
that are available today for the cell cycle analysis using scRNA-seq data.
Despite huge progress achieved, none of the packages can produce complete and
reliable results with respect to all aforementioned tasks. One of the major
difficulties for existing packages is distinguishing between two patterns of
cell cycle transcriptomic dynamics: normal and characteristic for embryonic
stem cells (ESC), with the latter one shared by many cancer cell lines.
Moreover, some cell lines are characterized by a mixture of two subpopulations,
one following the standard and one ESC-like cell cycle, which makes the
analysis even more challenging. In conclusion, we discuss the difficulties of
the analysis of cell cycle-related single cell transcriptome and provide
certain guidelines for the use of the existing methods.


  Transition metal chromophores with earth-abundant transition metals are an
important design target for their applications in lighting and non-toxic
bioimaging, but their design is challenged by the scarcity of complexes that
simultaneously have optimal target absorption energies in the visible region as
well as well-defined ground states. Machine learning (ML) accelerated discovery
could overcome such challenges by enabling screening of a larger space, but is
limited by the fidelity of the data used in ML model training, which is
typically from a single approximate density functional. To address this
limitation, we search for consensus in predictions among 23 density functional
approximations across multiple rungs of Jacobs ladder. To accelerate the
discovery of complexes with absorption energies in the visible region while
minimizing MR character, we use 2D efficient global optimization to sample
candidate low-spin chromophores from multi-million complex spaces. Despite the
scarcity (i.e., approx. 0.01\%) of potential chromophores in this large
chemical space, we identify candidates with high likelihood (i.e., > 10\%) of
computational validation as the ML models improve during active learning,
representing a 1,000-fold acceleration in discovery. Absorption spectra of
promising chromophores from time-dependent density functional theory verify
that 2/3 of candidates have the desired excited state properties. The
observation that constituent ligands from our leads have demonstrated
interesting optical properties in the literature exemplifies the effectiveness
of our construction of a realistic design space and active learning approach.


  Single-step retrosynthesis is the cornerstone of retrosynthesis planning,
which is a crucial task for computer-aided drug discovery. The goal of
single-step retrosynthesis is to identify the possible reactants that lead to
the synthesis of the target product in one reaction. By representing organic
molecules as canonical strings, existing sequence-based retrosynthetic methods
treat the product-to-reactant retrosynthesis as a sequence-to-sequence
translation problem. However, most of them struggle to identify diverse
chemical reactions for a desired product due to the deterministic inference,
which contradicts the fact that many compounds can be synthesized through
various reaction types with different sets of reactants. In this work, we aim
to increase reaction diversity and generate various reactants using discrete
latent variables. We propose a novel sequence-based approach, namely
RetroDVCAE, which incorporates conditional variational autoencoders into
single-step retrosynthesis and associates discrete latent variables with the
generation process. Specifically, RetroDVCAE uses the Gumbel-Softmax
distribution to approximate the categorical distribution over potential
reactions and generates multiple sets of reactants with the variational
decoder. Experiments demonstrate that RetroDVCAE outperforms state-of-the-art
baselines on both benchmark dataset and homemade dataset. Both quantitative and
qualitative results show that RetroDVCAE can model the multi-modal distribution
over reaction types and produce diverse reactant candidates.


  The ability to characterize proteins at sequence-level resolution is vital to
biological research. Currently, the leading method for protein sequencing is by
liquid chromatography mass spectrometry (LC-MS) whereas proteins are reduced to
their constituent peptides by enzymatic digest and subsequently analyzed on an
LC-MS instrument. The short peptide sequences that result from this analysis
are used to characterize the original protein content of the sample. Here we
present PASS, a de novo assembler for short peptide sequences that can be used
to reconstruct large portions of protein targets, a step that can facilitate
downstream sample characterization efforts. We show how, with adequate peptide
sequence coverage and little-to-no additional sequence processing, PASS
reconstructs protein sequences into relatively large (100 amino acid or longer)
contigs having high (93.1 - 99.1%) sequence identity to reference antibody
light and heavy chain proteins. Availability: PASS is released under the GNU
General Public License Version 3 (GPLv3) and is publicly available from
https://github.com/warrenlr/PASS


  Understanding the inherent timescales of large bubbles in DNA is critical to
a thorough comprehension of its physicochemical characteristics, as well as
their potential role on helix opening and biological function. In this work we
employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study
relaxation dynamics of large bubbles in homopolymer DNA, using simulations up
to the microsecond time scale. By studying energy autocorrelation functions of
relatively large bubbles inserted into thermalised DNA molecules, we extract
characteristic relaxation times from the equilibration process for both
adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of
different amplitudes and widths are investigated through extensive statistics
and appropriate fittings of their relaxation. Characteristic relaxation times
increase with bubble height and width. We show that, within the model,
relaxation times are two orders of magnitude longer in GC sequences than in AT
sequences. Overall, our results confirm that large bubbles leave a lasting
impact on the molecule's dynamics, for times between 0.4-200ns depending on the
homopolymer type and bubble shape, thus clearly affecting long-time evolutions
of the molecule.


  Molecular property prediction is a fundamental task in the drug and material
industries. Physically, the properties of a molecule are determined by its own
electronic structure, which is a quantum many-body system and can be exactly
described by the Schr"odinger equation. Full-range many-body interactions
between electrons have been proven effective in obtaining an accurate solution
of the Schr"odinger equation by classical computational chemistry methods,
although modeling such interactions consumes an expensive computational cost.
Meanwhile, deep learning methods have also demonstrated their competence in
molecular property prediction tasks. Inspired by the classical computational
chemistry methods, we design a novel method, namely GEM-2, which
comprehensively considers full-range many-body interactions in molecules.
Multiple tracks are utilized to model the full-range interactions between the
many-bodies with different orders, and a novel axial attention mechanism is
designed to approximate the full-range interaction modeling with much lower
computational cost. Extensive experiments demonstrate the overwhelming
superiority of GEM-2 over multiple baseline methods in quantum chemistry and
drug discovery tasks. The ablation studies also verify the effectiveness of the
full-range many-body interactions.


  Since the seminal work of Powell, the relationships between the population
growth rate, the probability distributions of generation time, and the
distribution of cell age have been known for the bacterial population in a
steady state of exponential growth. Here, we generalize these relationships to
include an unsteady (transient) state for both the batch culture and the mother
machine experiment. In particular, we derive a time-dependent Euler-Lotka
equation (relating the generation time distributions to the population growth
rate) and a generalization of the inequality between the mean generation time
and the population doubling time. To do this, we use a model proposed by
Lebowitz and Rubinow, in which each cell is described by its age and generation
time. We show that our results remain valid for a class of more complex models
that use other state variables in addition to cell age and generation time, as
long as the integration of these additional variables reduces the model to
Lebowitz-Rubinow form. As an application of this formalism, we calculate the
fitness landscapes for phenotypic traits (cell age, generation time) in a
population that is not growing exponentially. We clarify that the known fitness
landscape formula for the cell age as a phenotypic trait is an approximation to
the exact time-dependent formula.


  Deterministic approximations to stochastic Susceptible-Infectious-Susceptible
models typically predict a stable endemic steady-state when above threshold.
This can be hard to relate to the underlying stochastic dynamics, which has no
endemic steady-state but can exhibit approximately stable behaviour. Here we
relate the approximate models to the stochastic dynamics via the definition of
the quasi-stationary distribution (QSD), which captures this approximately
stable behaviour. We develop a system of ordinary differential equations that
approximate the number of infected individuals in the QSD for arbitrary contact
networks and parameter values. When the epidemic level is high, these QSD
approximations coincide with the existing approximation methods. However, as we
approach the epidemic threshold, the models deviate, with these models
following the QSD and the existing methods approaching the all susceptible
state. Through consistently approximating the QSD, the proposed methods provide
a more robust link to the stochastic models.


  Due to stochastic fluctuations arising from finite population size, known as
genetic drift, the ability of a population to explore a rugged fitness
landscape depends on its size. In the weak mutation regime, while the mean
steady-state fitness increases with population size, we find that the height of
the first fitness peak encountered when starting from a random genotype
displays various behaviors versus population size, even among small and simple
rugged landscapes. We show that the accessibility of the different fitness
peaks is key to determining whether this height overall increases or decreases
with population size. Furthermore, there is often a finite population size that
maximizes the height of the first fitness peak encountered when starting from a
random genotype. This holds across various classes of model rugged landscapes
with sparse peaks, and in some experimental and experimentally-inspired ones.
Thus, early adaptation in rugged fitness landscapes can be more efficient and
predictable for relatively small population sizes than in the large-size limit.


  Antibody design is valuable for therapeutic usage and biological research.
Existing deep-learning-based methods encounter several key issues: 1)
incomplete context for Complementarity-Determining Regions (CDRs) generation;
2) incapability of capturing the entire 3D geometry of the input structure; 3)
inefficient prediction of the CDR sequences in an autoregressive manner. In
this paper, we propose Multi-channel Equivariant Attention Network (MEAN) to
co-design 1D sequences and 3D structures of CDRs. To be specific, MEAN
formulates antibody design as a conditional graph translation problem by
importing extra components including the target antigen and the light chain of
the antibody. Then, MEAN resorts to E(3)-equivariant message passing along with
a proposed attention mechanism to better capture the geometrical correlation
between different components. Finally, it outputs both the 1D sequences and 3D
structure via a multi-round progressive full-shot scheme, which enjoys more
efficiency and precision against previous autoregressive approaches. Our method
significantly surpasses state-of-the-art models in sequence and structure
modeling, antigen-binding CDR design, and binding affinity optimization.
Specifically, the relative improvement to baselines is about 23% in
antigen-binding CDR design and 34% for affinity optimization.


  We characterize completely the capacity for (nondegenerate) multistationarity
of mass action reaction networks with one-dimensional stoichiometric subspace
in terms of reaction structure. Specifically, we show that networks with two or
more source complexes have the capacity for multistationarity if and only if
they have both patterns $(\to, \gets)$ and $(\gets, \to)$ in some 1D
projections. Moreover, we specify the classes of networks for which only
degenerate multiple steady states may occur. In particular, we characterize the
capacity for nondegenerate multistationarity of small networks composed of one
irreversible and one reversible reaction, or two reversible reactions


  Interactions between SARS-CoV-2 and human proteins (SARS-CoV-2 PPIs) cause
information transfer through biochemical pathways that contribute to the
immunopathology of COVID-19. Here, we present a communication network model of
the immune system to compute the information transferred by the viral proteins
using the available SARS-CoV-2 PPIs data. The amount of transferred information
depends on the reference state of the immune system, or the state without
SARS-CoV-2 PPIs, and can quantify how many variables of the immune system are
controlled by the viral proteins. The information received by the immune system
proteins from the viral proteins is useful to identify the biological processes
(BPs) susceptible to dysregulation, and also to estimate the duration of viral
PPIs necessary for the dysregulation to occur. We found that computing the drop
in information from viral PPIs due to drugs provides a direct measure for the
efficacy of therapies.


  Self-supervised learning (SSL) is a method that learns the data
representation by utilizing supervision inherent in the data. This learning
method is in the spotlight in the drug field, lacking annotated data due to
time-consuming and expensive experiments. SSL using enormous unlabeled data has
shown excellent performance for molecular property prediction, but a few issues
exist. (1) Existing SSL models are large-scale; there is a limitation to
implementing SSL where the computing resource is insufficient. (2) In most
cases, they do not utilize 3D structural information for molecular
representation learning. The activity of a drug is closely related to the
structure of the drug molecule. Nevertheless, most current models do not use 3D
information or use it partially. (3) Previous models that apply contrastive
learning to molecules use the augmentation of permuting atoms and bonds.
Therefore, molecules having different characteristics can be in the same
positive samples. We propose a novel contrastive learning framework,
small-scale 3D Graph Contrastive Learning (3DGCL) for molecular property
prediction, to solve the above problems. 3DGCL learns the molecular
representation by reflecting the molecule's structure through the pre-training
process that does not change the semantics of the drug. Using only 1,128
samples for pre-train data and 1 million model parameters, we achieved the
state-of-the-art or comparable performance in four regression benchmark
datasets. Extensive experiments demonstrate that 3D structural information
based on chemical knowledge is essential to molecular representation learning
for property prediction.


  Learning accurate drug representations is essential for task such as
computational drug repositioning. A drug hierarchy is a valuable source that
encodes knowledge of relations among drugs in a tree-like structure where drugs
that act on the same organs, treat the same disease, or bind to the same
biological target are grouped together. However, its utility in learning drug
representations has not yet been explored, and currently described drug
representations cannot place novel molecules in a drug hierarchy. Here, we
develop a semi-supervised drug embedding that incorporates two sources of
information: (1) underlying chemical grammar that is inferred from chemical
structures of drugs and drug-like molecules (unsupervised), and (2)
hierarchical relations that are encoded in an expert-crafted hierarchy of
approved drugs (supervised). We use the Variational Auto-Encoder (VAE)
framework to encode the chemical structures of molecules and use the drug-drug
similarity information obtained from the hierarchy to induce the clustering of
drugs in hyperbolic space. The hyperbolic space is amenable for encoding
hierarchical relations. Our qualitative results support that the learned drug
embedding can induce the hierarchical relations among drugs. We demonstrate
that the learned drug embedding can be used for drug repositioning.


  SARS-COV-2 is a positive single-strand RNA-based macromolecule that has
caused the death of more than 6.3 million people since June 2022. Moreover, by
disturbing global supply chains through lockdown, the virus has indirectly
caused devastating damage to the global economy. It is vital to design and
develop drugs for this virus and its various variants. In this paper, we
developed an in-silico study-based hybrid framework to repurpose existing
therapeutic agents in finding drug-like bioactive molecules that would cure
Covid-19. We employed the Lipinski rules on the retrieved molecules from the
ChEMBL database and found 133 drug-likeness bioactive molecules against SARS
coronavirus 3CL Protease. Based on standard IC50, the dataset was divided into
three classes active, inactive, and intermediate. Our comparative analysis
demonstrated that the proposed Extra Tree Regressor (ETR) based QSAR model has
improved prediction results related to the bioactivity of chemical compounds as
compared to Gradient Boosting, XGBoost, Support Vector, Decision Tree, and
Random Forest based regressor models. ADMET analysis is carried out to identify
thirteen bioactive molecules with ChEMBL IDs 187460, 190743, 222234, 222628,
222735, 222769, 222840, 222893, 225515, 358279, 363535, 365134 and 426898.
These molecules are highly suitable drug candidates for SARS-COV-2 3CL
Protease. In the next step, the efficacy of bioactive molecules is computed in
terms of binding affinity using molecular docking and then shortlisted six
bioactive molecules with ChEMBL IDs 187460, 222769, 225515, 358279, 363535, and
365134. These molecules can be suitable drug candidates for SARS-COV-2. It is
anticipated that the pharmacologist/drug manufacturer would further investigate
these six molecules to find suitable drug candidates for SARS-COV-2. They can
adopt these promising compounds for their downstream drug development stages.


  This paper conducts a systematic statistical analysis of the characteristics
of the geographical empirical distributions for the numbers of both cumulative
and daily confirmed COVID-19 cases and deaths at county, city, and state levels
over a time span from January 2020 to June 2022. The mathematical heavy-tailed
distributions can be used for fitting the empirical distributions observed in
different temporal stages and geographical scales. The estimations of the shape
parameter of the tail distributions using the Generalized Pareto Distribution
also support the observations of the heavy-tailed distributions. According to
the characteristics of the heavy-tailed distributions, the evolution course of
the geographical empirical distributions can be divided into three distinct
phases, namely the power-law phase, the lognormal phase I, and the lognormal
phase II. These three phases could serve as an indicator of the severity degree
of the COVID-19 pandemic within an area. The empirical results suggest
important intrinsic dynamics of a human infectious virus spread in the human
interconnected physical complex network. The findings extend previous empirical
studies and could provide more strict constraints for current mathematical and
physical modeling studies, such as the SIR model and its variants based on the
theory of complex networks.


  Function diversity, namely, the range of tasks individuals can perform, is
essential to productive organizations. This concept has been studied in
disparate discipline contexts, while general patterns and mechanisms remain
unclear. Here, we first analyze over five thousand observations of top-down
organizations -- US federal agencies, Norwegian companies, and US universities,
and find that the number of distinct functions scales with organizational size,
approximately as a power law with an exponent of 1/2. Further, we find common
patterns in the distribution of function abundance within organizations. This
universality suggests that human organizations, despite differences in their
purpose, structure, and culture, may share common mechanisms for creating
specializations. Additionally, we find that cities -- bottom-up organizations
-- differ from top-down organizations and exhibit logarithmic scaling. We
discuss potential avenues for modeling the mechanisms for these observations
using history-dependent random processes, and offer several criteria for model
selection.


  We analyze mathematical models for COVID-19 with discrete time delays and
vaccination. Sufficient conditions for the local stability of the endemic and
disease-free equilibrium points are proved for any positive time delay. The
stability results are illustrated through numerical simulations performed in
MATLAB.


  Understanding the spread of COVID-19 has been the subject of numerous
studies, highlighting the significance of reliable epidemic models. Here, we
introduce a novel epidemic model using a latent Hawkes process with temporal
covariates for modelling the infections. Unlike other models, we model the
reported cases via a probability distribution driven by the underlying Hawkes
process. Modelling the infections via a Hawkes process allows us to estimate by
whom an infected individual was infected. We propose a Kernel Density Particle
Filter (KDPF) for inference of both latent cases and reproduction number and
for predicting the new cases in the near future. The computational effort is
proportional to the number of infections making it possible to use particle
filter type algorithms, such as the KDPF. We demonstrate the performance of the
proposed algorithm on synthetic data sets and COVID-19 reported cases in
various local authorities in the UK, and benchmark our model to alternative
approaches.


  This survey for mathematicians summarizes several works by the author on
protein geometry and protein function with applications to viral glycoproteins
in general and the spike glycoprotein of the SARS-CoV-2 virus in particular.
Background biology and biophysics are sketched. This body of work culminates in
a postulate that protein secondary structure regulates mutation, with backbone
hydrogen bonds materializing in critical regions to avoid mutation, and
disappearing from other regions to enable it.


  The primary goal of this research is to investigate the impact of delay on
the dynamics of the Susceptible-Exposed-Infected-Recovered-Death and
Susceptible (SEIRDS) model, to which we add a stochastic term to account for
uncertainty in COVID-19 parameter estimations. We run two models, one
deterministic and one stochastic, and show that their solutions exist and are
unique. We also numerically investigate the impact of immunity loss on the
emerging time of a new wave, as well as the necessary condition for the
extinction and persistence of the disease.


  The COVID-19 pandemic has rapidly spread around the world and burdened public
health in almost all countries involving France. After the spread of
SARS-CoV-2, France harvested many deaths in total. In this paper, we develop
models with integer and fractional orders to investigate the dynamics of
COVID-19 transmission in French hospitals and intensive care units (ICUs).
Moreover, this paper aims to explore the impact of precautionary measures on
the total infected cases in hospitals and ICUs of COVID-19 for the entire
France by using available actual data.


  We prove that the exact closure of SIR pairwise epidemic equations on a
configuration model network is possible if and only if the degree distribution
is Poisson, Binomial, or Negative Binomial. The proof relies on establishing,
for these specific degree distributions, the equivalence of the closed pairwise
model and the so-called dynamical survival analysis (DSA) edge-based model
which was previously shown to be exact. Indeed, as we show here, the DSA model
is equivalent to the well-known edge-based Volz model. We use this result to
provide reductions of the closed pairwise and Volz models to the same single
equation involving only susceptibles, which has a useful statistical
interpretation in terms of the times to infection. We illustrate our findings
with some numerical examples.


  The pharmaceutical success of atorvastatin (ATV), a widely employed drug
against the "bad" cholesterol (LDL) and cardiovascular diseases, traces back to
its ability to scavenge free radicals. Unfortunately, information on its
antioxidant properties is missing or unreliable. Here, we report detailed
quantum chemical results for ATV and its ortho- and para-hydroxy metabolites
(o-ATV, p-ATV) in the methanolic phase. They comprise global reactivity
indices, bond order indices, and spin densities as well as all relevant
enthalpies of reaction (bond dissociation BDE, ionization IP and electron
attachment EA, proton detachment PDE and proton affinity PA, and electron
transfer ETE). With these properties in hand, we can provide the first
theoretical explanation of the experimental finding that, due to their free
radical scavenging activity, ATV hydroxy metabolites rather than the parent
ATV, have substantial inhibitory effect on LDL and the like. Surprisingly
(because it is contrary to the most cases currently known), we unambiguously
found that HAT (direct hydrogen atom transfer) rather than SPLET (sequential
proton loss electron transfer) or SET-PT (stepwise electron transfer proton
transfer) is the thermodynamically preferred pathway by which o-ATV and p-ATV
in methanolic phase can scavenge DPPH$^\bullet$ (1,1-diphenyl-2-picrylhydrazyl)
radicals. From a quantum chemical perspective, the ATV's species investigated
are surprising because of the nontrivial correlations between bond dissociation
energies, bond lengths, bond order indices and pertaining stretching
frequencies, which do not fit the framework of naive chemical intuition.


  Genome-wide Association Studies (GWASes) identify genomic variations that are
statistically associated with a trait, such as a disease, in a group of
individuals. Unfortunately, careless sharing of GWAS statistics might give rise
to privacy attacks. Several works attempted to reconcile secure processing with
privacy-preserving releases of GWASes. However, we highlight that these
approaches remain vulnerable if GWASes utilize overlapping sets of individuals
and genomic variations. In such conditions, we show that even when relying on
state-of-the-art techniques for protecting releases, an adversary could
reconstruct the genomic variations of up to 28.6% of participants, and that the
released statistics of up to 92.3% of the genomic variations would enable
membership inference attacks. We introduce I-GWAS, a novel framework that
securely computes and releases the results of multiple possibly interdependent
GWASes. I-GWAS continuously releases privacy-preserving and noise-free GWAS
results as new genomes become available.


  Biological living systems in general exhibit complex and diverse dynamics.
The latter, in particular, is essential, since diversification increases the
odds of survival of an organism while reducing the risk of extinction of the
population. Primarily, diversification is a consequence of the randomness in
the replication process of a biological cell, which eventually manifests into a
motley set of macroscopic features of an individual. These heterogeneous
features of individuals constitutes for diversity in population. Cancer is a
prime example of such a complex system where the transformed cells exhibit
plethora of disparate features, which in turn makes modeling their dynamics
quite challenging. In this paper we consider cancer as a prototype of a complex
living system and provide two contrasting perspective for studying and modeling
it. Based on this we illicit a deeper role of diversification in the evolution
of cancer. Following this, we ask ourselves how can we model these diverse
dynamics in a multiscale setting. We address the shortcoming of the existing
multiscale modeling technique by providing an abstract but mathematically
rigorous framework for deducing stochastic evolution equations at the
macroscopic level starting from a microscopic description of the involved
dynamics. We achieve this by making use of the connection between stochastic
process and the semigroup operator generated by them. In particular, we look at
the semigroups generated by Levy processes and their connection with the
characteristic functions and Levy symbols. The latter turns out to represent
pseudo-differential operators using which we eventually provide a mechanism for
constructing stochastic evolution equations. Altogether, this provides the
framework for modeling diverse dynamics at the macroscale starting from the
microscale.


  We address the task of controlled generation of small molecules, which
entails finding novel molecules with desired properties under certain
constraints (e.g., similarity to a reference molecule). Here we introduce
MolMIM, a probabilistic auto-encoder for small molecule drug discovery that
learns an informative and clustered latent space. MolMIM is trained with Mutual
Information Machine (MIM) learning, and provides a fixed length representation
of variable length SMILES strings. Since encoder-decoder models can learn
representations with ``holes'' of invalid samples, here we propose a novel
extension to the training procedure which promotes a dense latent space, and
allows the model to sample valid molecules from random perturbations of latent
codes. We provide a thorough comparison of MolMIM to several variable-size and
fixed-size encoder-decoder models, demonstrating MolMIM's superior generation
as measured in terms of validity, uniqueness, and novelty. We then utilize
CMA-ES, a naive black-box and gradient free search algorithm, over MolMIM's
latent space for the task of property guided molecule optimization. We achieve
state-of-the-art results in several constrained single property optimization
tasks as well as in the challenging task of multi-objective optimization,
improving over previous success rate SOTA by more than 5\% . We attribute the
strong results to MolMIM's latent representation which clusters similar
molecules in the latent space, whereas CMA-ES is often used as a baseline
optimization method. We also demonstrate MolMIM to be favourable in a compute
limited regime, making it an attractive model for such cases.


  Due to the aging of the world population and westernization of lifestyles,
the prevalence of neurodegenerative diseases such as Alzheimer's disease (AD)
and Parkinson's disease (PD) is rapidly rising and is expected to put a strong
socioeconomic burden on health systems worldwide. Due to the limited success of
clinical trials of therapies against neurodegenerative diseases, research has
extended its scope to a systems medicine point of view, with a particular focus
on the gastrointestinal-brain axis as a potential main actor in disease
development and progression. Microbiome as well as metabolome studies along the
gastrointestinal-brain axis have already revealed important insights into
disease pathomechanisms. Both the microbiome and metabolome can be easily
manipulated by dietary and lifestyle interventions, and might thus offer novel,
readily available therapeutic options to prevent the onset as well as the
progression of PD and AD. This review summarizes our current knowledge on the
association between microbiota, metabolites, and neurodegeneration in light of
the gastrointestinal-brain axis. In this context, we also illustrate
state-of-the art methods of microbiome and metabolome research as well as
metabolic modeling that facilitate the identification of disease
pathomechanisms. We conclude our review with therapeutic options to modulate
microbiome composition to prevent or delay neurodegeneration and illustrate
potential future research directions to fight PD and AD.


  In this article, we develop a deterministic compartmental model including the
spatio-temporal aspect and the wearing of masks by the general public for the
spread of the Covid-19 epidemic. This model is based on the SIR model of
Kermack and McKendrick and includes mobility terms that model the movement of
individuals wearing masks between Guadeloupe andMartinique. The objective of
the model is to analyze the effects of mobility and the wearing of masks by the
general public on the spread of the Covid-19 epidemic in these two
regions.Numerical simulations of this model show the importance of wearing
masks by the general public, even if they are not very inefficient, and
indicate that the model is capable of qualitatively simulating the propagation
trends of the Covid-19 epidemic. This work is part ofmy doctoral thesis, the
subject of which is entitled {\guillemotleft} stochastic modeling of Covid-19
in an island environment taking into account the spatio-temporal aspect
{\guillemotright}.


  While generative models have recently become ubiquitous in many scientific
areas, less attention has been paid to their evaluation. For molecular
generative models, the state-of-the-art examines their output in isolation or
in relation to its input. However, their biological and functional properties,
such as ligand-target interaction is not being addressed. In this study, a
novel biologically-inspired benchmark for the evaluation of molecular
generative models is proposed. Specifically, three diverse reference datasets
are designed and a set of metrics are introduced which are directly relevant to
the drug discovery process. In particular we propose a recreation metric, apply
drug-target affinity prediction and molecular docking as complementary
techniques for the evaluation of generative outputs. While all three metrics
show consistent results across the tested generative models, a more detailed
comparison of drug-target affinity binding and molecular docking scores
revealed that unimodal predictiors can lead to erroneous conclusions about
target binding on a molecular level and a multi-modal approach is thus
preferrable. The key advantage of this framework is that it incorporates prior
physico-chemical domain knowledge into the benchmarking process by focusing
explicitly on ligand-target interactions and thus creating a highly efficient
tool not only for evaluating molecular generative outputs in particular, but
also for enriching the drug discovery process in general.


  Under constant selection, each trait has a fixed fitness, and small mutation
rates allow populations to efficiently exploit the optimal trait. Therefore it
is reasonable to expect mutation rates will evolve downwards. However, we find
this need not be the case, examining several models of mutation. While upwards
evolution of mutation rate has been found with frequency or time dependent
fitness, we demonstrate its possibility in a much simpler context. This work
uses adaptive dynamics to study the evolution of mutation rate, and the
replicator-mutator equation to model trait evolution. Our approach differs from
previous studies by considering a wide variety of methods to represent
mutation. We use a finite string approach inspired by genetics, as well as a
model of local mutation on a discretization of the unit intervals, handling
mutation beyond the endpoints in three ways. The main contribution of this work
is a demonstration that the evolution of mutation rate can be significantly
more complicated than what is usually expected in relatively simple models.


  Single-strand breaks (SSBs) are the major DNA damage in the genome arising
spontaneously as the outcome of genotoxins and intermediates of DNA
transactions. SSBs play a crucial role in various biological processes and show
a non-random distribution in the genome. Several SSB detection approaches such
as S1 END-seq and SSiNGLe-ILM emerged to characterize the genomic landscape of
SSB with nucleotide resolution. However, these sequencing-based methods are
costly and unfeasible for large-scale analysis of diverse species. Thus, we
proposed the first computational approach, SSBlazer, which is an explainable
and scalable deep learning framework for genome-wide nucleotide-resolution SSB
site prediction. We demonstrated that SSBlazer can accurately predict SSB sites
and sufficiently alleviate false positives by constructing an imbalanced
dataset to simulate the realistic SSB distribution. The model interpretation
analysis reveals that SSBlazer captures the pattern of individual CpG in
genomic context and the motif of TGCC in the center region as critical
features. Besides, SSBlazer is a lightweight model with robust cross-species
generalization ability in the cross-species evaluation, which enables the
large-scale genome-wide application in diverse species. Strikingly, the
putative SSB genomic landscapes of 216 vertebrates reveal a negative
correlation between SSB frequency and evolutionary hierarchy, suggesting that
the genome tends to be integrity during evolution.


  Understanding the molecular mechanisms driving the binding between
bio-molecules is a crucial challenge in molecular biology. In this respect,
characteristics like the preferentially hydrophobic composition of the binding
interfaces, the role of van der Waals interactions (short range forces), and
the consequent shape complementarity between the interacting molecular surfaces
are well established. However, no consensus has yet been reached on how and how
much electrostatic participates in the various stages of protein-protein
interactions. Here, we perform extensive analyses on a large dataset of protein
complexes for which both experimental binding affinity and pH data were
available. We found that (i) although different classes of dimers do not
present marked differences in the amino acid composition and charges
disposition in the binding region, (ii) homodimers with identical binding
region show higher electrostatic compatibility with respect to both homodimers
with non-identical binding region and heterodimers. The level of electrostatic
compatibility also varies with the pH of the complex, reaching the lowest
values for low pH. Interestingly, (iii) shape and electrostatic complementarity
behave oppositely when one stratifies the complexes by their binding affinity.
Conversely, complexes with low values of binding affinity exploit Coulombic
complementarity to acquire specificity, suggesting that electrostatic
complementarity may play a greater role in transient (or less stable)
complexes. In light of these results, (iv) we provide a fast and efficient
method to measure electrostatic complementarity without the need of knowing the
complex structure. Expanding the electrostatic potential on a basis of 2D
orthogonal polynomials, we can discriminate between transient and permanent
protein complexes with an AUC of the ROC of 0.8.


  Ultrametric matrices have a rich structure that is not apparent from their
definition. Notably, the subclass of strictly ultrametric matrices are
covariance matrices of certain weighted rooted binary trees. In applications,
these matrices can be large and dense, making them difficult to store and
handle. In this manuscript, we exploit the underlying tree structure of these
matrices to sparsify them via a similarity transformation based on Haar-like
wavelets. We show that, with overwhelmingly high probability, only an
asymptotically negligible fraction of the off-diagonal entries in random but
large strictly ultrametric matrices remain non-zero after the transformation;
and develop a fast algorithm to compress such matrices directly from their tree
representation. We also identify the subclass of matrices diagonalized by the
wavelets and supply a sufficient condition to approximate the spectrum of
strictly ultrametric matrices outside this subclass. Our methods give
computational access to a covariance model of the microbiologists' Tree of
Life, which was previously inaccessible due to its size, and motivate defining
a new but wavelet-based phylogenetic $\beta$-diversity metric. Applying this
metric to a metagenomic dataset demonstrates that it can provide novel insight
into noisy high-dimensional samples and localize speciation events that may be
most important in determining relationships between environmental factors and
microbial composition.


  Pairwise sequence alignment is a very time-consuming step in common
bioinformatics pipelines. Speeding up this step requires heuristics, efficient
implementations, and/or hardware acceleration. A promising candidate for all of
the above is the recently proposed GenASM algorithm. We identify and address
three inefficiencies in the GenASM algorithm: it has a high amount of data
movement, a large memory footprint, and does some unnecessary work. We propose
Scrooge, a fast and memory-frugal genomic sequence aligner. Scrooge includes
three novel algorithmic improvements which reduce the data movement, memory
footprint, and the number of operations in the GenASM algorithm. We provide
efficient open-source implementations of the Scrooge algorithm for CPUs and
GPUs, which demonstrate the significant benefits of our algorithmic
improvements. For long reads, the CPU version of Scrooge achieves a 20.1x,
1.7x, and 2.1x speedup over KSW2, Edlib, and a CPU implementation of GenASM,
respectively. The GPU version of Scrooge achieves a 4.0x 80.4x, 6.8x, 12.6x and
5.9x speedup over the CPU version of Scrooge, KSW2, Edlib, Darwin-GPU, and a
GPU implementation of GenASM, respectively. We estimate an ASIC implementation
of Scrooge to use 3.6x less chip area and 2.1x less power than a GenASM ASIC
while maintaining the same throughput. Further, we systematically analyze the
throughput and accuracy behavior of GenASM and Scrooge under various
configurations. As the best configuration of Scrooge depends on the computing
platform, we make several observations that can help guide future
implementations of Scrooge. Availability: https://github.com/CMU-SAFARI/Scrooge


  The persistence of biodiversity of species is a challenging proposition in
ecological communities in the face of Darwinian selection. The present article
investigates beyond the pairwise competitive interactions and provides a novel
perspective for understanding the influence of higher-order interactions on the
evolution of social phenotypes. Our simple model yields a prosperous outlook to
demonstrate the impact of perturbations on intransitive competitive
higher-order interactions. Using a mathematical technique, we show how alone
the perturbed interaction network can quickly determine the coexistence
equilibrium of competing species instead of solving a large system of ordinary
differential equations. It is possible to split the system into multiple
feasible cluster states depending on the number of perturbations. Our analysis
also reveals the ratio between the unperturbed and perturbed species is
inversely proportional to the amount of employed perturbation. Our results
suggest that nonlinear dynamical systems and interaction topologies can be
interplayed to comprehend species' coexistence under adverse conditions.
Particularly our findings signify that less competition between two species
increases their abundance and outperforms others.


  The conventional approach for analyzing gene expression data involves
clustering algorithms. Cluster analyses provide partitioning of the set of
genes that can predict biological classification based on its similarity in
n-dimensional space. In this study, we investigate whether network analysis
will provide an advantage over the traditional approach. We identify the
advantages and disadvantages of using the value-based and the rank-based
construction in creating a graph representation of the original gene-expression
data in a time-series format. We tested four community detection algorithms,
namely, the Clauset-Newman-Moore (greedy), Louvain, Leiden, and Girvan-Newman
algorithms in predicting the 5 functional groups of genes. We used the Adjusted
Rand Index to assess the quality of the predicted communities with respect to
the biological classifications. We showed that Girvan-Newman outperforms the 3
modularity-based algorithms in both value-based and ranked-based constructed
graphs. Moreover, we also show that when compared to the conventional
clustering algorithms such as K-means, Spectral, Birch, and Agglomerative
algorithms, we obtained a higher ARI with Girvan-Newman. This study also
provides a tool for graph construction, visualization, and community detection
for further analysis of gene expression data.


  Natural language processing (NLP) is an area of artificial intelligence that
applies information technologies to process the human language, understand it
to a certain degree, and use it in various applications. This area has rapidly
developed in the last few years and now employs modern variants of deep neural
networks to extract relevant patterns from large text corpora. The main
objective of this work is to survey the recent use of NLP in the field of
pharmacology. As our work shows, NLP is a highly relevant information
extraction and processing approach for pharmacology. It has been used
extensively, from intelligent searches through thousands of medical documents
to finding traces of adversarial drug interactions in social media. We split
our coverage into five categories to survey modern NLP methodology, commonly
addressed tasks, relevant textual data, knowledge bases, and useful programming
libraries. We split each of the five categories into appropriate subcategories,
describe their main properties and ideas, and summarize them in a tabular form.
The resulting survey presents a comprehensive overview of the area, useful to
practitioners and interested observers.


  Accurate prediction of the protein-ligand binding affinities is an essential
challenge in the structure-based drug design. Despite recent advance in
data-driven methods in affinity prediction, their accuracy is still limited,
partially because they only take advantage of static crystal structures while
the actual binding affinities are generally depicted by the thermodynamic
ensembles between proteins and ligands. One effective way to approximate such a
thermodynamic ensemble is to use molecular dynamics (MD) simulation. Here, we
curated an MD dataset containing 3,218 different protein-ligand complexes, and
further developed Dynaformer, which is a graph-based deep learning model.
Dynaformer was able to accurately predict the binding affinities by learning
the geometric characteristics of the protein-ligand interactions from the MD
trajectories. In silico experiments demonstrated that our model exhibits
state-of-the-art scoring and ranking power on the CASF-2016 benchmark dataset,
outperforming the methods hitherto reported. Moreover, we performed a virtual
screening on the heat shock protein 90 (HSP90) using Dynaformer that identified
20 candidates and further experimentally validated their binding affinities. We
demonstrated that our approach is more efficient, which can identify 12 hit
compounds (two were in the submicromolar range), including several newly
discovered scaffolds. We anticipate this new synergy between large-scale MD
datasets and deep learning models will provide a new route toward accelerating
the early drug discovery process.


  All-atom simulations have become increasingly popular to study conformational
and dynamical properties of nucleic acids as they are accurate and provide high
spatial and time resolutions. This high resolution however comes at a heavy
computational cost and within the time scales of simulations nucleic acids
weakly fluctuate around their ideal structure exploring a limited set of
conformations. We introduce the RBB-NA algorithm which is capable of
controlling rigid base parameters in all-atom simulations of Nucleic Acids.
With suitable biasing potentials this algorithm can "force" a DNA or RNA
molecule to assume specific values of the six rotational (tilt, roll, twist,
buckle, propeller, opening) and/or the six translational parameters (shift,
slide, rise, shear, stretch, stagger). The algorithm enables the use of
advanced sampling techniques to probe the structure and dynamics of locally
strongly deformed Nucleic Acids. We illustrate its performance showing some
examples in which DNA is strongly twisted, bent or locally buckled. In these
examples RBB-NA reproduces well the unconstrained simulations data and other
known features of DNA mechanics, but it also allows one to explore the
anharmonic behavior characterizing the mechanics of nucleic acids in the high
deformation regime.


  In this paper we present a new approach to solve semi-supervised
classification tasks for biomedical applications, involving a supervised
autoencoder network. We create a network architecture that encodes labels into
the latent space of an autoencoder, and define a global criterion combining
classification and reconstruction losses. We train the Semi-Supervised
AutoEncoder (SSAE) on labelled data using a double descent algorithm. Then, we
classify unlabelled samples using the learned network thanks to a softmax
classifier applied to the latent space which provides a classification
confidence score for each class.
  We implemented our SSAE method using the PyTorch framework for the model,
optimizer, schedulers, and loss functions. We compare our semi-supervised
autoencoder method (SSAE) with classical semi-supervised methods such as Label
Propagation and Label Spreading, and with a Fully Connected Neural Network
(FCNN). Experiments show that the SSAE outperforms Label Propagation and
Spreading and the Fully Connected Neural Network both on a synthetic dataset
and on two real-world biological datasets.


  Malaria is a vector-borne disease that exacts a grave toll in the Global
South. The epidemiology of Plasmodium vivax, the most geographically expansive
agent of human malaria, is characterised by the accrual of a reservoir of
dormant parasites known as hypnozoites. Relapses, arising from hypnozoite
activation events, comprise the majority of the blood-stage infection burden,
with implications for the acquisition of immunity and the distribution of
superinfection. Here, we construct a hybrid transmission model for P. vivax
that concurrently accounts for the accrual of the hypnozoite reservoir,
(blood-stage) superinfection and the acquisition of immunity. We begin by
analytically characterising within-host dynamics as a function of
mosquito-to-human transmission intensity, extending our previous model
(comprising an open network of infinite server queues) to capture a discretised
immunity level. To model transmission-blocking and antidisease immunity, we
allow for geometric decay in the respective probabilities of successful
human-to-mosquito transmission and symptomatic blood-stage infection as a
function of this immunity level. Under a hybrid approximation -- whereby
probabilistic within-host distributions are cast as expected population-level
proportions -- we couple host and vector dynamics to recover a deterministic
compartmental model in line with Ross-Macdonald theory. We then perform a
steady-state analysis for this compartmental model, informed by the (analytic)
distributions derived at the within-host level. To characterise transient
dynamics, we derive a reduced system of integrodifferential equations (IDEs),
likewise informed by our within-host queueing network, allowing us to recover
population-level distributions for various quantities of epidemiological
interest. Our model provides insights into important, but poorly understood,
epidemiological features of P. vivax.


  Therapeutic modulation of immune states is central to the treatment of human
disease. However, how drugs and drug combinations impact the diverse cell types
in the human immune system remains poorly understood at the transcriptome
scale. Here, we apply single-cell mRNA-seq to profile the response of human
immune cells to 502 immunomodulatory drugs alone and in combination. We develop
a unified mathematical model that quantitatively describes the transcriptome
scale response of myeloid and lymphoid cell types to individual drugs and drug
combinations through a single inferred regulatory network. The mathematical
model reveals how drug combinations generate novel, macrophage and T-cell
states by recruiting combinations of gene expression programs through both
additive and non-additive drug interactions. A simplified drug response algebra
allows us to predict the continuous modulation of immune cell populations
between activated, resting and hyper-inhibited states through combinatorial
drug dose titrations. Our results suggest that transcriptome-scale mathematical
models could enable the design of therapeutic strategies for programming the
human immune system using combinations of therapeutics.


  The cambium is a meristematic tissue in plant stems. Here, cell divisions
occur that are required for radial growth of plant stems. Daughters of cell
divisions within the cambium differentiate into woody xylem cells towards the
inside of the stem, or phloem towards the outside. As such, a pattern of
xylem-cambium-phloem is present along the radial axis of the stem. A
ligand-receptor pair, TDIF-PXY promotes cell division in the cambium, as do the
phytohormones, cytokinin and auxin. An auxin response factor, MP, has been
proposed to initiate cambial cell divisions by promoting PXY expression,
however, MP has also been reported to repress cambial cell divisions later in
development where TDIF-PXY complexes are also reported to suppress MP activity.
Here, we used a mathematical modelling approach to investigate how MP cell
division-promoting activity and cell division-repressing activity might be
integrated into the same network as a negative feedback loop. In our model,
this feedback loop improved the ability of the cambium to pattern correctly and
was found to be required for normal patterning when MP was stable. The
implications of this model in early and late cambium development are discussed.


  Protein-protein interaction (PPI) networks consist of the physical and/or
functional interactions between the proteins of an organism. Since the
biophysical and high-throughput methods used to form PPI networks are
expensive, time-consuming, and often contain inaccuracies, the resulting
networks are usually incomplete. In order to infer missing interactions in
these networks, we propose a novel class of link prediction methods based on
continuous-time classical and quantum random walks. In the case of quantum
walks, we examine the usage of both the network adjacency and Laplacian
matrices for controlling the walk dynamics. We define a score function based on
the corresponding transition probabilities and perform tests on four real-world
PPI datasets. Our results show that continuous-time classical random walks and
quantum walks using the network adjacency matrix can successfully predict
missing protein-protein interactions, with performance rivalling the state of
the art.


  The Coefficient of Parentage (COP) between two individuals is the expected
inbreeding of their offspring. Originally exploited by animal breeders, is now
a routine calculation among plant breeders as part of crop improvement
programs. Here we show that the COP between strains requires a different
calculation than the used to calculate the COP between individuals. Failure to
do so may result in an overestimation of the amount of inbreeding. Here we
provide a simple methodology to calculate the correct coefficient of parentage
between strains.


  The Goldman-Hodgkin-Katz (GHK) equations have been widely applied to ion
channel studies, simulations, and model developments. However, they are
constructed under a constant electric field, causing them to have a low degree
of approximation in the prediction of ionic fluxes, electric currents, and
membrane potentials. In this paper, the equations are extended from the
constant electric field to the nonlinear electric field induced by charges from
an ionic solution and an ion channel protein. Furthermore, a novel numerical
quadrature scheme is developed to estimate one major parameter, called the
extension parameter, of the extended GHK equations in terms of a set of
electrostatic potential values. To this end, the extended GHK equations become
a bridge between the "macroscopic" ion channel kinetics and the "microscopic"
electrostatic potential values across a cell membrane. To generate a set of
required electrostatic potential values, a nonlinear finite element iterative
scheme for solving a one-dimensional Poisson-Nernst-Planck ion channel model is
developed and implemented as a Python software package. This package is then
used to do numerical studies on the extended GHK equations, the numerical
quadrature scheme, and the nonlinear iterative scheme. Numerical results
confirm the importance of considering charge effects in the calculation of
ionic fluxes. They also validate the high numerical accuracy of the numerical
quadrature scheme, the fast convergence rate of the nonlinear iterative scheme,
and the high performance of the software package.


  Epidemic models are a valuable tool in the decision making process. Once a
mathematical model for an epidemics has been established, the very next step is
calculating a mathematical expression for the basic reproductive number, $R_0$,
which is the average number of infections caused by an individual that is
introduced in a population of susceptibles. Finding a mathematical expression
for $R_0$ is important because it allows to analyze the effect of the different
parameters in the model on $R_0$ so that we can act on them to keep $R_0 < 1$,
so that the epidemic fades out. In this work we show how to calculate $R_0$ in
complicated epidemic models by using only basic concepts of Markov chains.


  Epidemic forecasting is the key to effective control of epidemic transmission
and helps the world mitigate the crisis that threatens public health. To better
understand the transmission and evolution of epidemics, we propose EpiGNN, a
graph neural network-based model for epidemic forecasting. Specifically, we
design a transmission risk encoding module to characterize local and global
spatial effects of regions in epidemic processes and incorporate them into the
model. Meanwhile, we develop a Region-Aware Graph Learner (RAGL) that takes
transmission risk, geographical dependencies, and temporal information into
account to better explore spatial-temporal dependencies and makes regions aware
of related regions' epidemic situations. The RAGL can also combine with
external resources, such as human mobility, to further improve prediction
performance. Comprehensive experiments on five real-world epidemic-related
datasets (including influenza and COVID-19) demonstrate the effectiveness of
our proposed method and show that EpiGNN outperforms state-of-the-art baselines
by 9.48% in RMSE.


  Predicting protein-protein interactions from sequences is an important goal
of computational biology. Various sources of information can be used to this
end. Starting from the sequences of two interacting protein families, one can
use phylogeny or residue coevolution to infer which paralogs are specific
interaction partners within each species. We show that these two signals can be
combined to improve the performance of the inference of interaction partners
among paralogs. For this, we first align the sequence-similarity graphs of the
two families through simulated annealing, yielding a robust partial pairing. We
next use this partial pairing to seed a coevolution-based iterative pairing
algorithm. This combined method improves performance over either separate
method. The improvement obtained is striking in the difficult cases where the
average number of paralogs per species is large or where the total number of
sequences is modest.


  COVID-19 will be a continuous threat to human population despite having a few
vaccines at hand until we reach the endemic state through natural herd immunity
and total immunization through universal vaccination. However, the vaccine acts
as a practical tool for reducing the massive public health problem and the
emerging economic consequences that the continuing COVID -19 epidemic is
causing worldwide, while the vaccine efficacy wanes. In this work, we propose
and analyze an epidemic model of
Susceptible-Exposed-Infected-Recovered-Vaccinated (SEIRV) population taking
into account the rate of vaccination and vaccine waning. The dynamics of the
model has been investigated, and the condition for a disease-free endemic
equilibrium state is obtained. Further, the analysis is extended to study the
COVID-19 spread in India by considering the availability of vaccines and the
related critical parameters such as vaccination rate, vaccine efficacy and
waning of vaccine's impact on deciding the emerging fate of this epidemic. We
have also discussed the conditions for herd immunity due to vaccinated
individuals among the people. Our results highlight the importance of vaccines,
the effectiveness of booster vaccination in protecting people from infection,
and their importance in epidemic and pandemic modelling.


  The viral load is known to be a chief predictor of the risk of transmission
of infectious diseases. In this work, we investigate the role of the
individuals' viral load in the disease transmission by proposing a new
susceptible-infectious-recovered epidemic model for the densities and mean
viral loads of each compartment. To this aim, we formally derive the
compartmental model from an appropriate microscopic one. Firstly, we consider a
multi-agent system in which individuals are identified by the epidemiological
compartment to which they belong and by their viral load. Microscopic rules
describe both the switch of compartment and the evolution of the viral load. In
particular, in the binary interactions between susceptible and infectious
individuals, the probability for the susceptible individual to get infected
depends on the viral load of the infectious individual. Then, we implement the
prescribed microscopic dynamics in appropriate kinetic equations, from which
the macroscopic equations for the densities and viral load momentum of the
compartments are eventually derived. In the macroscopic model, the rate of
disease transmission turns out to be a function of the mean viral load of the
infectious population. We analytically and numerically investigate the case
that the transmission rate linearly depends on the viral load, which is
compared to the classical case of constant transmission rate. A qualitative
analysis is performed based on stability and bifurcation theory. Finally,
numerical investigations concerning the model reproduction number and the
epidemic dynamics are presented.


  Foraging and acquiring of food is a delicate balance between managing the
costs, both energy and social, and individual preferences. Previous research on
the solitary foraging of free ranging dogs showed that they prioritized the
nutritionally highest valued food patch first but do not ignore other less
valuable food either, displaying typical scavenger behaviour. The current
experiment was carried out on groups of dogs with the same set up to see the
change in foraging strategies, if any, under the influence of social cost like
intra-group competition. We found multiple differences between the strategies
of dogs foraging alone versus in groups with competition playing an implicit
role in the decision making of dogs when foraging in groups. Dogs were able to
continually assess and evaluate the available resources in a patch and adjust
their behaviour accordingly. Foraging in groups also provided benefits of
reduced individual vigilance. The various decisions and choices made seemed to
have a basis in the optimal foraging theory wherein the dogs harvested the
nutritionally richest patch possible with the least risk and cost involved but
was willing to compromise if that was not possible. This underscores the
cognitive, quick decision-making abilities and adaptable behaviour of these
dogs.


  The disruption of circadian rhythm is a cardinal symptom for Alzheimer's
disease (AD) patients. The full circadian rhythm orchestration of gene
expression in the human brain and its inherent associations with AD remain
largely unknown. We present a novel comprehensive approach, PRIME, to detect
and analyze rhythmic oscillation patterns in untimed high-dimensional gene
expression data across multiple datasets. To demonstrate the utility of PRIME,
firstly, we validate it by a time course expression dataset from mouse liver as
a cross-species and cross-organ validation. Then, we apply it to study
oscillation patterns in untimed genome-wide gene expression from 19 human brain
regions of controls and AD patients. Our findings reveal clear, synchronized
oscillation patterns in 15 pairs of brain regions of control, while these
oscillation patterns either disappear or dim for AD. It is worth noting that
PRIME discovers the circadian rhythmic patterns without requiring the sample's
timestamps. The codes for PRIME, along with codes to reproduce the figures in
this paper, are available at https://github.com/xinxingwu-uk/PRIME.


  The global rise in anthropogenic reactive nitrogen (N) and the negative
impacts of N deposition on terrestrial plant diversity are well-documented. The
R* theory of resource competition predicts reversible decreases in plant
diversity in response to N loading. However, empirical evidence for the
reversibility of N-induced biodiversity loss is mixed. In a long-term
N-enrichment experiment in Minnesota, a low-diversity state that emerged during
N addition has persisted for decades after additions ceased. Hypothesized
mechanisms preventing recovery of biodiversity include nutrient recycling,
insufficient external seed supply, and litter inhibition of plant growth. Here
we present an ODE model that unifies these mechanisms, produces bistability at
intermediate N inputs, and qualitatively matches the observed hysteresis at
Cedar Creek. Key features of the model, including native species' growth
advantage in low-N conditions and limitation by litter accumulation, generalize
from Cedar Creek to North American grasslands. Our results suggest that
effective biodiversity restoration in these systems may require management
beyond reducing N inputs, such as burning, grazing, haying, and seed additions.
By coupling resource competition with an additional inter-specific inhibitory
process, the model also illustrates a general mechanism for bistability and
hysteresis that may occur in multiple ecosystem types.


  Summary: Mapping ancient DNA to a reference genome is challenging as it
involves numerous steps, is time-consuming and has to be repeated within a
study to assess the quality of extracts and libraries; as a result, the mapping
needs to be automatized to handle large amounts of data in a reproducible way.
We present mapache, a flexible, robust, and scalable pipeline to map, quantify
and impute ancient and present-day DNA in a reproducible way. Mapache is
implemented in the workflow manager Snakemake and is optimized for low-space
consumption, allowing to efficiently (re)map large data sets such as reference
panels and multiple extracts and libraries.


  Saffron (Crocus sativus L.) is a triploid, sterile, monocot plant belonging
to the family Iridaceae, sub-family Crocoideae. C.sativus only blooms once a
year and should be collected within a very short duration, the stigmas of
Saffron flowers are harvested manually and subjected to desiccation then have
been used as a spice. It has been also used as a drug to treat tumors, cancer,
chronic uterine hemorrhage, insomnia, scarlet fever, smallpox, colds, and
cardiovascular disorders. It has been shown that saffron is a protective agent
against chromosomal damage. Saffron has been vegetatively propagated by corm,
each mother corm produces 7-8 cormlets each year. The main colors of saffron,
crocetin, and crocetin glycosides, and the main flavors, picrocrocin, Safranal
is the main component of aroma and Its bitter taste is related to
Glycosidepicrocrocin that is derived from the oxidative cleavage of the
carotenoid, zeaxanthin.by zeaxanthin cleavage dioxygenase (ZCD). We
investigated gene expression of ZCD in vitro by using tissue culture of
perianth obtained from immature flora buds of Ghaen of Khorasan province, Iran,
on MS medium supplemented by 10 mg/L NAA and BAP. RNA of each sample was
extracted using RNx method; followed by RT-PCR techniques. The results
indicated that ZCD was present in perianth of all cultured samples of mentioned
areas. Investigation of this pathway which controls the saffron apocarotenoid
pigments in perianth is important to produce saffron with high quality and
quantity. The results showed that the expression of ZCD in perianth whether
through mevalonic acid or non-mevalonic acid needs more investigation


  Since COVID-19 vaccines became available, no studies have quantified how
different disaster evacuation strategies can mitigate pandemic risks in
shelters. Therefore, we applied an age-structured epidemiological model, known
as the Susceptible-Exposed-Infectious-Recovered (SEIR) model, to investigate to
what extent different vaccine uptake levels and the Diversion protocol
implemented in Taiwan decrease infections and delay pandemic peak occurrences.
Taiwan's Diversion protocol involves diverting those in self-quarantine due to
exposure, thus preventing them from mingling with the general public at a
congregate shelter. The Diversion protocol, combined with sufficient vaccine
uptake, can decrease the maximum number of infections and delay outbreaks
relative to scenarios without such strategies. When the diversion of all
exposed people is not possible, or vaccine uptake is insufficient, the
Diversion protocol is still valuable. Furthermore, a group of evacuees that
consists primarily of a young adult population tends to experience pandemic
peak occurrences sooner and have up to 180% more infections than does a
majority elderly group when the Diversion protocol is implemented. However,
when the Diversion protocol is not enforced, the majority elderly group suffers
from up to 20% more severe cases than the majority young adult group.


  Dynein motors exhibit catch bonding, where the unbinding rate of the motors
from microtubule filaments decreases with increasing opposing load. The
implications of this catch bond on the transport properties of dynein-driven
cargo are yet to be fully understood. In this context, optical trapping assays
constitute an important means of accurately measuring the forces generated by
molecular motor proteins. We investigate, using theory and stochastic
simulations, the transport properties of cargo transported by catch bonded
dynein molecular motors - both singly and in teams - in a harmonic potential,
which mimics the variable force experienced by cargo in an optical trap. We
estimate the biologically relevant measures of first passage time - the time
during which the cargo remains bound to the microtubule and detachment force
-the force at which the cargo unbinds from the microtubule, using both
two-dimensional and one-dimensional force balance frameworks. Our results
suggest that even for cargo transported by a single motor, catch bonding may
play a role depending on the force scale which marks the onset of the catch
bond. By comparing with experimental measurements on single dynein-driven
transport, we estimate realistic bounds of this catch bond force scale.
Generically, catch bonding results in increased persistent motion, and can also
generate non-monotonic behaviour of first passage times. For cargo transported
by multiple motors, emergent collective effects due to catch bonding can result
in non-trivial re-entrant phenomena wherein average first passage times and
detachment forces exhibit non-monotonic behaviour as a function of the stall
force and the motor velocity.


  Protein structures and functions are determined by a contiguous arrangement
of amino acid sequences. Designing novel protein sequences and structures with
desired geometry and functions is a complex task with large state spaces. Here
we develop a novel protein design pipeline consisting of two deep learning
algorithms, ProteinSolver and AlphaFold2. ProteinSolver is a deep graph neural
network that generates amino acid sequences such that the forces between
interacting amino acids are favorable and compatible with the fold while
AlphaFold2 is a deep learning algorithm that predicts the protein structures
from protein sequences. We present forty de novo designed binding sites of the
PTP1B and P53 proteins with high precision, out of which thirty proteins are
novel. Using ProteinSolver and AlphaFold2 in conjunction, we can trim the
exploration of the large protein conformation space, thus expanding the ability
to find novel and diverse de novo protein designs.


  Modeling and control of epidemics such as the novel Corona virus have assumed
paramount importance at a global level. A natural and powerful dynamical
modeling framework to use in this context is a continuous time Markov decision
process (CTMDP) that encompasses classical compartmental paradigms such as the
Susceptible-Infected-Recovered (SIR) model. The challenges with CTMDP based
models motivate the need for a more efficient approach and the mean field
approach offers an effective alternative. The mean field approach computes the
collective behavior of a dynamical system comprising numerous interacting nodes
(where nodes represent individuals in the population). This paper (a) presents
an overview of the mean field approach to epidemic modeling and control and (b)
provides a state-of-the-art update on recent advances on this topic. Our
discussion in this paper proceeds along two specific threads. The first thread
assumes that the individual nodes faithfully follow a socially optimal control
policy prescribed by a regulatory authority. The second thread allows the
individual nodes to exhibit independent, strategic behavior. In this case, the
strategic interaction is modeled as a mean field game and the control is based
on the associated mean field Nash equilibria. In this paper, we start with a
discussion of modeling of epidemics using an extended compartmental model -
SIVR and provide an illustrative example. We next provide a review of relevant
literature, using a mean field approach, on optimal control of epidemics,
dealing with how a regulatory authority may optimally contain epidemic spread
in a population. Following this, we provide an update on the literature on the
use of the mean field game based approach in the study of epidemic spread and
control. We conclude the paper with relevant future research directions.


  We investigate a cyclic game system where organisms face an epidemic beyond
being threatened by natural enemies. As a survival strategy, individuals of one
out of the species usually safeguard themselves by approaching the enemies of
their enemies and performing social distancing to escape contamination when an
outbreak affects the neighbourhood. We simulate how the survival movement
strategy to local epidemic surges must adapt if a pathogen mutation makes the
disease deadlier. We study the spatial distribution of local outbreaks and
observe the influence of disease mortality on individuals' spatial
organisation. We show that adapting the survival movement strategy for a high
mortality disease demands an altruistic behaviour of the organisms since their
death risk increases. Despite weakening the disease transmission chain, which
benefits the species, abandoning refuges provided by safeguarding social
interaction increases the vulnerability to being eliminated in the cyclic game.
Considering that not all individuals exhibit altruism, we find the relative
growth in the species density as a function of the proportion of individuals
behaving altruistically. Our results may be helpful for biologists and data
scientists to understand how adaptive altruistic processes can affect
population dynamics in complex systems.


  Aging is thought to be a consequence of intrinsic breakdowns in how genetic
information is processed. But mounting experimental evidence suggests that
aging can be slowed. To help resolve this mystery, I derive a mortality
equation which characterizes the dynamics of an evolving population with a
given maximum age. Remarkably, while the spectrum of eigenvalues that govern
the evolution depends on the fitness, how they change with the maximum age is
independent of fitness. This makes it possible to establish the conditions
under which programmed aging can provide an evolutionary benefit.


  Developing and maintaining life requires a lot of computation. This is done
by gene regulatory networks. But we have little understanding of how this
computation is organized. I show that there is a direct correspondence between
the structural and functional building blocks of regulatory networks, which I
call regulatory motifs. I derive a simple bound on the range of function that
these motifs can perform, in terms of the local network structure. I prove that
this range is a small fraction of all possible functions, which severely
constrains global network behavior. Part of this restriction is due to
redundancy in the function that regulatory motifs can achieve - there are many
ways to perform the same task. Regulatory motifs help us understanding how
genetic computation is organized and what it can achieve.


  The decisive role of Embryology in understanding the evolution of animal
forms is founded and deeply rooted in the history of science. It is recognized
that the emergence of multicellularity would not have been possible without the
formation of the first embryo. We speculate that biophysical phenomena and the
surrounding environment of the Ediacaran ocean were instrumental in co-opting a
neoplastic functional module (NFM) within the nucleus of the first zygote.
Thus, the neoplastic process, understood here as a biological phenomenon with
profound embryologic implications, served as the evolutionary engine that
favored the formation of the first embryo and cancerous diseases and allowed to
coherently create and recreate body shapes in different animal groups during
evolution. In this article, we provide a deep reflection on the Physics of the
first embryogenesis and its contribution to the exaptation of additional NFM
components, such as the extracellular matrix. Knowledge of NFM components,
structure, dynamics, and origin advances our understanding of the numerous
possibilities and different innovations that embryos have undergone to create
animal forms via Neoplasia during evolutionary radiation. The developmental
pathways of Neoplasia have their origins in ctenophores and were consolidated
in mammals and other apical groups.


  Minimum flow decomposition (MFD) -- the problem of finding a minimum set of
weighted source-to-sink paths that perfectly decomposes a flow -- is a
classical problem in Computer Science, and variants of it are powerful models
in different fields such as Bioinformatics and Transportation. Even on acyclic
graphs, the problem is NP-hard, and most practical solutions have been via
heuristics or approximations. While there is an extensive body of research on
acyclic graphs, currently, there is no \emph{exact} solution on graphs with
cycles. In this paper, we present the first ILP formulation for three natural
variants of the MFD problem in graphs with cycles, asking for a decomposition
consisting only of weighted source-to-sink paths or cycles, trails, and walks,
respectively. On three datasets of increasing levels of complexity from both
Bioinformatics and Transportation, our approaches solve any instance in under
10 minutes. Our implementations are freely available at
github.com/algbio/MFD-ILP.


  Adaptation to changing environments is a universal feature of life and can
involve the organism modifying itself in response to the environment as well as
actively modifying the environment to control selection pressures. The latter
case couples the organism to environment. Then, how quickly should the organism
change in response to the environment? We formulate this question in terms of
how memory duration scales with environmental rate of change when there are
trade-offs in remembering vs. forgetting. We derive a universal scaling law for
optimal memory duration, taking into account memory precision as well as two
components of environmental volatility, bias and stability. We find sublinear
scaling with any amount of environmental volatility. We use a memory complexity
measure to explore the strategic conditions (game dynamics) favoring actively
reducing environmental volatility -- outsourcing memory through niche
construction -- over investing in neural tissue. We predict stabilizing niche
construction will evolve when neural tissue is costly, the environment is
variable, and it is beneficial to be able to encode a rich repertoire of
environmental states.


  Transcription Factors (TFs) are proteins that regulate gene expression. The
regulation mechanism is via the binding of a TF to a specific part of the gene
associated with it, the TF's target. The target of a specific TF corresponds to
a vanishingly small part of the entire DNA, where at the same time the search
must end in a matter of tens of seconds at most for its biological purpose to
be fulfilled - this makes the search a problem of high interest. Facilitated
Diffusion is a mechanism used in nature for a robust and efficient search
process. This mechanism combines 1D diffusion along the DNA and "excursions" of
diffusion in 3D that help the TF to quickly arrive at distant parts of the DNA.
In this paper we provide a derivation concerning this mechanism that links this
search process to fundamental concepts in probability theory (conditional
probability).


  We study a three-species cyclic game system where organisms face a contagious
disease whose virulence may change by a pathogen mutation. As a responsive
defence strategy, organisms' mobility is restricted to reduce disease
dissemination in the system. The impact of the collective self-preservation
strategy on the disease infection risk is investigated by performing stochastic
simulations of the spatial version of the rock-paper-scissors game. Our
outcomes show that the mobility control strategy induces plasticity in the
spatial patterns with groups of organisms of the same species inhabiting
spatial domains whose characteristic length scales depend on the level of
dispersal restrictions. The spatial organisation plasticity allows the
ecosystems to adapt to minimise the individuals' disease contamination risk if
an eventual pathogen alters the disease virulence. We discover that if a
pathogen mutation makes the disease more transmissible or less lethal, the
organisms benefit more if the mobility is not strongly restricted, thus forming
large spatial domains. Conversely, the benefits of protecting against a
pathogen causing a less contagious or deadlier disease are maximised if the
average size of groups of individuals of the same species is significantly
limited, reducing the dimensions of groups of organisms significantly. Our
findings may help biologists understand the effects of dispersal control as a
conservation strategy in ecosystems affected by epidemic outbreaks.


  Analyzing a functional genomics experiment, such as ATAC-, ChIP- or
RNA-sequencing, requires reference data including a genome assembly and gene
annotation. These resources can generally be retrieved from different
organizations and in different versions. Most bioinformatic workflows require
the user to supply this genomic data manually, which can be a tedious and
error-prone process.
  Here we present genomepy, which can search, download, and preprocess the
right genomic data for your analysis. Genomepy can search genomic data on NCBI,
Ensembl, UCSC and GENCODE, and compare available gene annotations to enable an
informed decision. The selected genome and gene annotation can be downloaded
and preprocessed with sensible, yet controllable, defaults. Additional
supporting data can be automatically generated or downloaded, such as aligner
indexes, genome metadata and blacklists.
  Genomepy is freely available at https://github.com/vanheeringen-lab/genomepy
under the MIT license and can be installed through pip or bioconda.


  Motivation: The development of novel compounds targeting proteins of interest
is one of the most important tasks in the pharmaceutical industry. Deep
generative models have been applied to targeted molecular design and have shown
promising results. Recently, target-specific molecule generation has been
viewed as a translation between the protein language and the chemical language.
However, such a model is limited by the availability of interacting
protein-ligand pairs. On the other hand, large amounts of unlabeled protein
sequences and chemical compounds are available and have been used to train
language models that learn useful representations. In this study, we propose
exploiting pretrained biochemical language models to initialize (i.e. warm
start) targeted molecule generation models. We investigate two warm start
strategies: (i) a one-stage strategy where the initialized model is trained on
targeted molecule generation (ii) a two-stage strategy containing a
pre-finetuning on molecular generation followed by target specific training. We
also compare two decoding strategies to generate compounds: beam search and
sampling.
  Results: The results show that the warm-started models perform better than a
baseline model trained from scratch. The two proposed warm-start strategies
achieve similar results to each other with respect to widely used metrics from
benchmarks. However, docking evaluation of the generated compounds for a number
of novel proteins suggests that the one-stage strategy generalizes better than
the two-stage strategy. Additionally, we observe that beam search outperforms
sampling in both docking evaluation and benchmark metrics for assessing
compound quality.
  Availability and implementation: The source code is available at
https://github.com/boun-tabi/biochemical-lms-for-drug-design and the materials
are archived in Zenodo at https://doi.org/10.5281/zenodo.6832145