  We analysed publicly available data on place of occurrence of COVID-19 deaths
from national statistical agencies in the UK between March 9 2020 and February
28 2021. We introduce a modified Weibull model that describes the deaths due to
COVID-19 at a national and place of occurrence level. We observe similar trends
in the UK where deaths due to COVID-19 first peak in Homes, followed by
Hospitals and Care Homes 1-2 weeks later in the first and second waves. This is
in line with the infectious period of the disease, indicating a possible
transmission vehicle between the settings. Our results show that the first wave
is characterised by fast growth and a slow reduction after the peak in deaths
due to COVID-19. The second and third waves have the converse property, with
slow growth and a rapid decrease from the peak. This difference may result from
behavioural changes in the population (social distancing, masks, etc). Finally,
we introduce a double logistic model to describe the dynamic proportion of
COVID-19 deaths occurring in each setting. This analysis reveals that the
proportion of COVID-19 deaths occurring in Care Homes increases from the start
of the pandemic and past the peak in total number of COVID-19 deaths in the
first wave. After the catastrophic impact in the first wave, the proportion of
COVID-19 deaths occurring in Care Homes gradually decreased from is maximum
after the first wave indicating residence were better protected in the second
and third waves compared to the first.


  We propose a qualitative analysis of a recent fractional-order COVID-19
model. We start by showing that the model is mathematically and biologically
well posed. Then, we give a proof on the global stability of the disease free
equilibrium point. Finally, some numerical simulations are performed to ensure
stability and convergence of the disease free equilibrium point.


  Motivation: Alignment-free (AF) distance/similarity functions are a key tool
for sequence analysis. Experimental studies on real datasets abound and, to
some extent, there are also studies regarding their control of false positive
rate (Type I error). However, assessment of their power, i.e., their ability to
identify true similarity, has been limited to some members of the D2 family by
experimental studies on short sequences, not adequate for current applications,
where sequence lengths may vary considerably. Such a State of the Art is
methodologically problematic, since information regarding a key feature such as
power is either missing or limited. Results: By concentrating on a
representative set of word-frequency based AF functions, we perform the first
coherent and uniform evaluation of the power, involving also Type I error for
completeness. Two Alternative models of important genomic features (CIS
Regulatory Modules and Horizontal Gene Transfer), a wide range of sequence
lengths from a few thousand to millions, and different values of k have been
used. As a result, we provide a characterization of those AF functions that is
novel and informative. Indeed, we identify weak and strong points of each
function considered, which may be used as a guide to choose one for analysis
tasks. Remarkably, of the fifteen functions that we have considered, only four
stand out, with small differences between small and short sequence length
scenarios. Finally, in order to encourage the use of our methodology for
validation of future AF functions, the Big Data platform supporting it is
public.


  Does regulation in the genome use collective behavior, similar to the way the
brain or deep neural networks operate? Here I make the case for why having a
genomic network capable of a high level of computation would be strongly
selected for, and suggest how it might arise from biochemical processes that
succeed in regulating in a collective manner, very different than the usual way
we think about genetic regulation.


  This paper introduces a new multi-strain epidemic model with reinfection and
cross-immunity to provide insights into the resurgence of the COVID-19 epidemic
in an area with reportedly high seroprevalence due to a largely unmitigated
outbreak: the state of Amazonas, Brazil. Although high seroprevalence could
have been expected to trigger herd immunity and prevent further waves in the
state, we have observed persistent levels of infection after the first wave and
eventually the emergence of a second viral strain just before an augmented
second wave. Our experiments suggest that the persistent levels of infection
after the first wave may be due to reinfection, whereas the higher peak at the
second wave can be explained by the emergence of the second variant and a low
level of cross-immunity between the original and the second variant. Finally,
the proposed model provides insights into the effect of reinfection and
cross-immunity on the long-term spread of an unmitigated epidemic.


  Testing is recommended for all close contacts of confirmed COVID-19 patients.
However, existing group testing methods are oblivious to the circumstances of
contagion provided by contact tracing. Here, we build upon a well-known
semi-adaptive pool testing method, Dorfman's method with imperfect tests, and
derive a simple group testing method based on dynamic programming that is
specifically designed to use the information provided by contact tracing.
Experiments using a variety of reproduction numbers and dispersion levels,
including those estimated in the context of the COVID-19 pandemic, show that
the pools found using our method result in a significantly lower number of
tests than those found using standard Dorfman's method, especially when the
number of contacts of an infected individual is small. Moreover, our results
show that our method can be more beneficial when the secondary infections are
highly overdispersed.


  Background. The supertree problem, i.e., the task of finding a common
refinement of a set of rooted trees is an important topic in mathematical
phylogenetics. The special case of a common leaf set $L$ is known to be
solvable in linear time. Existing approaches refine one input tree using
information of the others and then test whether the results are isomorphic.
  Results. A linear-time algorithm, LinCR, for constructing the common
refinement $T$ of $k$ input trees with a common leaf set is proposed that
explicitly computes the parent function of $T$ in a bottom-up approach.
  Conclusion. LinCR is simpler to implement than other asymptotically optimal
algorithms for the problem and outperforms the alternatives in empirical
comparisons.
  Availability. An implementation of LinCR in Python is freely available at
https://github.com/david-schaller/tralda.


  Using brain imaging quantitative traits (QTs) to identify the genetic risk
factors is an important research topic in imaging genetics. Many efforts have
been made via building linear models, e.g. linear regression (LR), to extract
the association between imaging QTs and genetic factors such as single
nucleotide polymorphisms (SNPs). However, to the best of our knowledge, these
linear models could not fully uncover the complicated relationship due to the
loci's elusive and diverse impacts on imaging QTs. Though deep learning models
can extract the nonlinear relationship, they could not select relevant genetic
factors. In this paper, we proposed a novel multi-task deep feature selection
(MTDFS) method for brain imaging genetics. MTDFS first adds a multi-task
one-to-one layer and imposes a hybrid sparsity-inducing penalty to select
relevant SNPs making significant contributions to abnormal imaging QTs. It then
builds a multi-task deep neural network to model the complicated associations
between imaging QTs and SNPs. MTDFS can not only extract the nonlinear
relationship but also arms the deep neural network with the feature selection
capability. We compared MTDFS to both LR and single-task DFS (DFS) methods on
the real neuroimaging genetic data. The experimental results showed that MTDFS
performed better than both LR and DFS in terms of the QT-SNP relationship
identification and feature selection. In a word, MTDFS is powerful for
identifying risk loci and could be a great supplement to the method library for
brain imaging genetics.


  The Beyond-Linear-Use-of-Equation-Superposition (BLUES) function method is
extended to coupled nonlinear ordinary differential equations and applied to
the epidemiological SIRS model with vaccination. Accurate analytic
approximations are obtained for the time evolution of the susceptible and
infected population fractions. The results are compared with those obtained
with alternative methods, notably Adomian decomposition, variational iteration
and homotopy perturbation. In contrast with these methods, the BLUES iteration
converges rapidly, globally, and captures the exact asymptotic behavior for
long times. The time of the infection peak is calculated using the BLUES
approximants and the results are compared with numerical solutions, which
indicate that the method is able to generate useful analytic expressions that
coincide with the (numerically) exact ones already for a small number of
iterations.


  Clinical decision-making in oncology involves multimodal data such as
radiology scans, molecular profiling, histopathology slides, and clinical
factors. Despite the importance of these modalities individually, no deep
learning framework to date has combined them all to predict patient prognosis.
Here, we predict the overall survival (OS) of glioma patients from diverse
multimodal data with a Deep Orthogonal Fusion (DOF) model. The model learns to
combine information from multiparametric MRI exams, biopsy-based modalities
(such as H&E slide images and/or DNA sequencing), and clinical variables into a
comprehensive multimodal risk score. Prognostic embeddings from each modality
are learned and combined via attention-gated tensor fusion. To maximize the
information gleaned from each modality, we introduce a multimodal
orthogonalization (MMO) loss term that increases model performance by
incentivizing constituent embeddings to be more complementary. DOF predicts OS
in glioma patients with a median C-index of 0.788 +/- 0.067, significantly
outperforming (p=0.023) the best performing unimodal model with a median
C-index of 0.718 +/- 0.064. The prognostic model significantly stratifies
glioma patients by OS within clinical subsets, adding further granularity to
prognostic clinical grading and molecular subtyping.


  Drug-target interaction (DTI) prediction plays a crucial role in drug
discovery, and deep learning approaches have achieved state-of-the-art
performance in this field. We introduce an ensemble of deep learning models
(EnsembleDLM) for DTI prediction. EnsembleDLM only uses the sequence
information of chemical compounds and proteins, and it aggregates the
predictions from multiple deep neural networks. This approach not only achieves
state-of-the-art performance in Davis and KIBA datasets but also reaches
cutting-edge performance in the cross-domain applications across different
bio-activity types and different protein classes. We also demonstrate that
EnsembleDLM achieves a good performance (Pearson correlation coefficient and
concordance index > 0.8) in the new domain with approximately 50% transfer
learning data, i.e., the training set has twice as much data as the test set.


  Planning for the protection of species often involves difficult choices about
which species to prioritize, given constrained resources. One way of
prioritizing species is to consider their "evolutionary distinctiveness", i.e.
their relative evolutionary isolation on a phylogenetic tree. Several
evolutionary isolation metrics or phylogenetic diversity indices have been
introduced in the literature, among them the so-called Fair Proportion index
(also known as the "evolutionary distinctiveness" score). This index apportions
the total diversity of a tree among all leaves, thereby providing a simple
prioritization criterion for conservation.
  Here, we focus on the prioritization order obtained from the Fair Proportion
index and analyze the effects of species extinction on this ranking. More
precisely, we analyze the extent to which the ranking order may change when
some species go extinct and the Fair Proportion index is re-computed for the
remaining taxa. We show that for each phylogenetic tree, there are edge lengths
such that the extinction of one leaf per cherry completely reverses the
ranking. Moreover, we show that even if only the lowest ranked species goes
extinct, the ranking order may drastically change. We end by analyzing the
effects of these two extinction scenarios (extinction of the lowest ranked
species and extinction of one leaf per cherry) for a collection of empirical
and simulated trees. In both cases, we can observe significant changes in the
prioritization orders, highlighting the empirical relevance of our theoretical
findings.


  Malaria is a mosquito-borne disease with a devastating global impact.
Plasmodium vivax is a major cause of human malaria beyond sub-Saharan Africa.
Relapsing infections, driven by a reservoir of liver-stage parasites known as
hypnozoites, present unique challenges for the control of P. vivax malaria.
Following indeterminate dormancy periods, hypnozoites may activate to trigger
relapses. Clearance of the hypnozoite reservoir through drug treatment (radical
cure) has been proposed as a potential tool for the elimination of P. vivax
malaria. Here, we introduce a stochastic, within-host model to jointly
characterise hypnozoite and infection dynamics for an individual in a general
transmission setting, allowing for radical cure. We begin by extending an
existing activation-clearance model for a single hypnozoite, adapted to both
short- and long-latency strains, to include drug treatment. We then embed this
activation-clearance model in an epidemiological framework accounting for
repeated mosquito inoculation and the administration of radical cure. By
constructing an open network of infinite server queues, we derive analytic
expressions for several quantities of epidemiological significance, including
the size of the hypnozoite reservoir; the relative contribution of relapses to
the infection burden; the distribution of multiple infections; the cumulative
number of recurrences over time, and the time to first recurrence following
drug treatment. By deriving, rather than assuming parameteric forms, we
characterise the transient dynamics of the hypnozoite reservoir following
radical cure more accurately than previous approaches. To yield
population-level insights, our analytic within-host distributions can be
embedded in multiscale models. Our work thus contributes to the epidemiological
understanding of the effects of radical cure on P. vivax malaria.


  We present a methodology for systematically extending epidemic models to
multilevel and multiscale spatio-temporal pandemic ones. Our approach builds on
the use of coloured stochastic and continuous Petri nets facilitating the sound
component-based extension of basic SIR models to include population
stratification and also spatio-geographic information and travel connections,
represented as graphs, resulting in robust stratified pandemic metapopulation
models. This method is inherently easy to use, producing scalable and reusable
models with a high degree of clarity and accessibility which can be read either
in a deterministic or stochastic paradigm. Our method is supported by a
publicly available platform PetriNuts; it enables the visual construction and
editing of models; deterministic, stochastic and hybrid simulation as well as
structural and behavioural analysis. All the models are available as
supplementary material, ensuring reproducibility.


  Structural biology aims at characterizing the structural and dynamic
properties of biological macromolecules at atomic details. Gaining insight into
three dimensional structures of biomolecules and their interactions is critical
for understanding the vast majority of cellular processes, with direct
applications in health and food sciences. Since 2010, the WeNMR project
(www.wenmr.eu) has implemented numerous web-based services to facilitate the
use of advanced computational tools by researchers in the field, using the high
throughput computing infrastructure provided by EGI. These services have been
further developed in subsequent initiatives under H2020 projects and are now
operating as Thematic Services in the European Open Science Cloud (EOSC) portal
(www.eosc-portal.eu), sending >12 millions of jobs and using around 4000
CPU-years per year. Here we review 10 years of successful e-infrastructure
solutions serving a large worldwide community of over 23,000 users to date,
providing them with user-friendly, web-based solutions that run complex
workflows in structural biology. The current set of active WeNMR portals are
described, together with the complex backend machinery that allows distributed
computing resources to be harvested efficiently.


  SARS-CoV-2 is the third betacoronavirus to enter the human population in the
past 20 years, revealing a concerning pattern. Clearly, preventing a future
pandemic from such viruses is a critical priority. Previous studies have shown
that shRNAs can be powerful suppressors of RNA viruses in transgenic animals
and substantially reduce transmission. Thus, we propose the introduction of
anti-betacoronavirus shRNAs using CRISPR/CAS9 gene drive into the horseshoe bat
population, the natural reservoir of those viruses, to combat this pandemic
threat at its source. Importantly, our approach is not expected to create any
harm to bats and can benefit other animals in the ecosystem that contract
betacoronaviruses from bats. We map the ethical and the technical aspects and
suggest guidelines for moving forward with this proposal.


  Biological phenotypes are products of complex evolutionary processes in which
selective forces influence multiple biological trait measurements in unknown
ways. Phylogenetic factor analysis disentangles these relationships across the
evolutionary history of a group of organisms. Scientists seeking to employ this
modeling framework confront numerous modeling and implementation decisions, the
details of which pose computational and replicability challenges. General and
impactful community employment requires a data scientific analysis plan that
balances flexibility, speed and ease of use, while minimizing model and
algorithm tuning. Even in the presence of non-trivial phylogenetic model
constraints, we show that one may analytically address latent factor
uncertainty in a way that (a) aids model flexibility, (b) accelerates
computation (by as much as 500-fold) and (c) decreases required tuning. We
further present practical guidance on inference and modeling decisions as well
as diagnosing and solving common problems in these analyses. We codify this
analysis plan in an automated pipeline that distills the potentially
overwhelming array of modeling decisions into a small handful of (typically
binary) choices. We demonstrate the utility of these methods and analysis plan
in four real-world problems of varying scales.


  Cryo-EM reconstruction algorithms seek to determine a molecule's 3D density
map from a series of noisy, unlabeled 2D projection images captured with an
electron microscope. Although reconstruction algorithms typically model the 3D
volume as a generic function parameterized as a voxel array or neural network,
the underlying atomic structure of the protein of interest places well-defined
physical constraints on the reconstructed structure. In this work, we exploit
prior information provided by an atomic model to reconstruct distributions of
3D structures from a cryo-EM dataset. We propose Cryofold, a generative model
for a continuous distribution of 3D volumes based on a coarse-grained model of
the protein's atomic structure, with radial basis functions used to model atom
locations and their physics-based constraints. Although the reconstruction
objective is highly non-convex when formulated in terms of atomic coordinates
(similar to the protein folding problem), we show that gradient descent-based
methods can reconstruct a continuous distribution of atomic structures when
initialized from a structure within the underlying distribution. This approach
is a promising direction for integrating biophysical simulation, learned neural
models, and experimental data for 3D protein structure determination.


  We consider the problem of binomiality of the steady state ideals of
biochemical reaction networks. We are interested in finding polynomial
conditions on the parameters such that the steady state ideal of a chemical
reaction network is binomial under every specialisation of the parameters if
the conditions on the parameters hold. We approach the binomiality problem
using Comprehensive Gr\"obner systems. Considering rate constants as
parameters, we compute comprehensive Gr\"obner systems for various reactions.
In particular, we make automatic computations on n-site phosphorylations and
biomodels from the Biomodels repository using the grobcov library of the
computer algebra system Singular.


  During the COVID-19 pandemic, conflicting opinions on physical distancing
swept across social media, affecting both human behavior and the spread of
COVID-19. Inspired by such phenomena, we construct a two-layer multiplex
network for the coupled spread of a disease and conflicting opinions. We model
each process as a contagion. On one layer, we consider the concurrent evolution
of two opinions -- pro-physical-distancing and anti-physical-distancing -- that
compete with each other and have mutual immunity to each other. The disease
evolves on the other layer, and individuals are less likely (respectively, more
likely) to become infected when they adopt the pro-physical-distancing
(respectively, anti-physical-distancing) opinion. We develop approximations of
mean-field type by generalizing monolayer pair approximations to multilayer
networks; these approximations agree well with Monte Carlo simulations for a
broad range of parameters and several network structures. Through numerical
simulations, we illustrate the influence of opinion dynamics on the spread of
the disease from complex interactions both between the two conflicting opinions
and between the opinions and the disease. We find that lengthening the duration
that individuals hold an opinion may help suppress disease transmission, and we
demonstrate that increasing the cross-layer correlations or intra-layer
correlations of node degrees may lead to fewer individuals becoming infected
with the disease.


  A rooted tree $T$ with vertex labels $t(v)$ and set-valued edge labels
$\lambda(e)$ defines maps $\delta$ and $\varepsilon$ on the pairs of leaves of
$T$ by setting $\delta(x,y)=q$ if the last common ancestor $\text{lca}(x,y)$ of
$x$ and $y$ is labeled $q$, and $m\in \varepsilon(x,y)$ if $m\in\lambda(e)$ for
at least one edge $e$ along the path from $\text{lca}(x,y)$ to $y$. We show
that a pair of maps $(\delta,\varepsilon)$ derives from a tree $(T,t,\lambda)$
if and only if there exists a common refinement of the (unique) least-resolved
vertex labeled tree $(T_{\delta},t_{\delta})$ that explains $\delta$ and the
(unique) least resolved edge labeled tree
$(T_{\varepsilon},\lambda_{\varepsilon})$ that explains $\varepsilon$ (provided
both trees exist). This result remains true if certain combinations of labels
at incident vertices and edges are forbidden.


  According to the fundamental principle of evolutionary game theory, the more
successful strategy in a population should spread. Hence, during a strategy
imitation process a player compares its payoff value to the payoff value held
by a competing strategy. But this information is not always accurate. To avoid
ambiguity a learner may therefore decide to collect a more reliable statistics
by averaging the payoff values of its opponents in the neighborhood, and makes
a decision afterwards. This simple alteration of the standard microscopic
protocol significantly improves the cooperation level in a population.
Furthermore, the positive impact can be strengthened by increasing the role of
the environment and the size of the evaluation circle. The mechanism that
explains this improvement is based on a self-organizing process which reveals
the detrimental consequence of defector aggregation that remains partly hidden
during face-to-face comparisons. Notably, the reported phenomenon is not
limited to lattice populations but remains valid also for systems described by
irregular interaction networks.


  Recently a novel framework has been proposed for designing the molecular
structure of chemical compounds using both artificial neural networks (ANNs)
and mixed integer linear programming (MILP). In the framework, we first define
a feature vector $f(C)$ of a chemical graph $C$ and construct an ANN that maps
$x=f(C)$ to a predicted value $\eta(x)$ of a chemical property $\pi$ to $C$.
After this, we formulate an MILP that simulates the computation process of
$f(C)$ from $C$ and that of $\eta(x)$ from $x$. Given a target value $y^*$ of
the chemical property $\pi$, we infer a chemical graph $C^\dagger$ such that
$\eta(f(C^\dagger))=y^*$ by solving the MILP. In this paper, we use linear
regression to construct a prediction function $\eta$ instead of ANNs. For this,
we derive an MILP formulation that simulates the computation process of a
prediction function by linear regression. The results of computational
experiments suggest our method can infer chemical graphs with around up to 50
non-hydrogen atoms.


  We formulate a compartmental model for the propagation of a respiratory
disease in a patchy environment. The patches are connected through the mobility
of individuals, and we assume that disease transmission and recovery are
possible during travel. Moreover, the migration terms are assumed to depend on
the distance between patches and the perceived severity of the disease. The
positivity and boundedness of the model solutions are discussed. We
analytically show the existence and global asymptotic stability of the
disease-free equilibrium. We study three different network topologies
numerically and find that underlying network structure is crucial for disease
transmission. Further numerical simulations reveal that infection during travel
has the potential to change the stability of disease-free equilibrium from
stable to unstable. The coupling strength and transmission coefficients are
also very crucial in disease propagation. Different exit screening scenarios
indicate that the patch with the highest prevalence may have adverse effects
but other patches will be benefited from exit screening. Furthermore, while
studying the multi-strain dynamics, it is observed that two co-circulating
strains will not persist simultaneously in the community but only one of the
strains may persist in the long run. Transmission coefficients corresponding to
the second strain are very crucial and show threshold like behavior with
respect to the equilibrium density of the second strain.


  We report on new stability conditions for evolutionary dynamics in the
context of population games. We adhere to the prevailing framework consisting
of many agents, grouped into populations, that interact noncooperatively by
selecting strategies with a favorable payoff. Each agent is repeatedly allowed
to revise its strategy at a rate referred to as revision rate. Previous
stability results considered either that the payoff mechanism was a memoryless
potential game, or allowed for dynamics (in the payoff mechanism) at the
expense of precluding any explicit dependence of the agents' revision rates on
their current strategies. Allowing the dependence of revision rates on
strategies is relevant because the agents' strategies at any point in time are
generally unequal. To allow for strategy-dependent revision rates and payoff
mechanisms that are dynamic (or memoryless games that are not potential), we
focus on an evolutionary dynamics class obtained from a straightforward
modification of one that stems from the so-called impartial pairwise comparison
strategy revision protocol. Revision protocols consistent with the modified
class retain from those in the original one the advantage that the agents
operate in a fully decentralized manner and with minimal information
requirements - they need to access only the payoff values (not the mechanism)
of the available strategies. Our main results determine conditions under which
system-theoretic passivity properties are assured, which we leverage for
stability analysis.


  Drug repurposing provides an opportunity to redeploy drugs, which ideally are
already approved for use in humans, for the treatment of other diseases. For
example, the repurposing of dexamethasone and baricitinib has played a crucial
role in saving patient lives during the ongoing SARS-CoV-2 pandemic. There
remains a need to expand therapeutic approaches to prevent life-threatening
complications in patients with COVID-19. Using an in silico approach based on
structural similarity to drugs already in clinical trials for COVID-19,
potential drugs were predicted for repurposing. For a subset of identified
drugs with different targets to their corresponding COVID-19 clinical trial
drug, a mechanism of action analysis was applied to establish whether they
might have a role in inhibiting the replication of SARS-CoV-2. Of sixty drugs
predicted in this study, two with the potential to inhibit SARS-CoV-2
replication were identified using mechanism of action analysis. Triamcinolone
is a corticosteroid that is structurally similar to dexamethasone; gallopamil
is a calcium channel blocker that is structurally similar to verapamil. In
silico approaches indicate possible mechanisms of action for both drugs in
inhibiting SARS-CoV-2 replication. The identification of these drugs as
potentially useful for patients with COVID-19 who are at a higher risk of
developing severe disease supports the use of in silico approaches to
facilitate quick and cost-effective drug repurposing. Such drugs could expand
the number of treatments available to patients who are not protected by
vaccination.


  Short Read Alignment Mapping Metrics (SRAMM): is an efficient and versatile
command line tool providing additional short read mapping metrics, filtering,
and graphs. Short read aligners report MAPing Quality (MAPQ), but these methods
generally are neither standardized nor well described in literature or software
manuals. Additionally, third party mapping quality programs are typically
computationally intensive or designed for specific applications. SRAMM
efficiently generates multiple different concept-based mapping scores to
provide for an informative post alignment examination and filtering process of
aligned short reads for various downstream applications. SRAMM is compatible
with Python 2.6+ and Python 3.6+ on all operating systems. It works with any
short read aligner that generates SAM/BAM/CRAM file outputs and reports 'AS'
tags. It is freely available under the MIT license at
http://github.com/achon/sramm.


  Non-pharmaceutical interventions(NPIs) play an important role in the early
stage control of COVID-19 pandemic. Vaccination is considered to be the
inevitable course to stop the spread of SARS-CoV-2. Based on the mechanism, a
SVEIR COVID-19 model with vaccination and NPIs is proposed. By means of the
basic reproduction number $R_{0}$, it is shown that the disease-free
equilibrium is globally attractive if $\mathscr{R}_{0}<1$, and COVID-19 is
uniform persistence if $\mathscr{R}_{0}>1$. Taking Indian dates for example in
the numerical simulation, we find that our dynamical results fits well with the
statistical dates. Consequently, we forecast the spreading trend of COVID-19
pandemic in India. Furthermore, our results imply that improving the intensity
of NPIs will greatly reduce the number of confirmed cases. Especially, NPIs are
indispensable even if all the people were vaccinated when the efficiency of
vaccine is relatively low. By simulating the relation ships of the basic
reproduction number $\mathscr{R}_{0}$, the vaccination rate and the efficiency
of vaccine, we find that it is impossible to achieve the herd immunity without
NPIs when the efficiency of vaccine is lower than $76.9\%$. Therefore, the herd
immunity area is defined by the evolution of relationships between the
vaccination rate and the efficiency of vaccine. In the study of two patchy, we
give the conditions for India and China to be open to navigation. Furthermore,
an appropriate dispersal of population between India and China is obtained. A
discussion completes the paper.


  Autocatalysis is a deceptively simple concept, referring to the situation
that a chemical species $X$ catalyzes its own formation. From the perspective
of chemical kinetics, autocatalysts show a regime of super-linear growth. Given
a chemical reaction network, however, it is not at all straightforward to
identify species that are autocatalytic in the sense that there is a
sub-network that takes $X$ as input and produces more than one copy of $X$ as
output. The difficulty arises from the need to distinguish autocatalysis e.g.
from the superposition of a cycle that consumes and produces equal amounts of
$X$ and a pathway that produces $X$. To deal with this issue, a number of
competing notions, such as exclusive autocatalysis and autocatalytic cycles,
have been introduced. A closer inspection of concepts and their usage by
different authors shows, however, that subtle differences in the definitions
often makes conceptually matching ideas difficult to bring together formally.
In this contribution we make some of the available approaches comparable by
translating them into a common formal framework that uses integer hyperflows as
a basis to study autocatalysis in large chemical reaction networks. As an
application we investigate the prevalence of autocatalysis in metabolic
networks.


  Ecosystems around the world are at risk of critical transitions due to
increasing anthropogenic pressures and climate change. Yet it is unclear where
the risks are higher or where in the world ecosystems are more vulnerable. Here
I measure resilience of primary productivity proxies for marine and terrestrial
ecosystems globally. Up to 29% of global terrestrial ecosystem, and 24% marine
ones, show symptoms of resilience loss. These symptoms are shown in all biomes,
but Arctic tundra and boreal forest are the most affected, as well as the
Indian Ocean and Eastern Pacific. Although the results are likely an
underestimation, they enable the identification of risk areas as well as the
potential synchrony of some transitions, helping prioritize areas for
management interventions and conservation.


  The interplay of biological, social, structural and random factors makes
disease forecasting extraordinarily complex. The course of an epidemic exhibits
average growth dynamics determined by features of the pathogen and the
population, yet also features significant variability reflecting the stochastic
nature of disease spread. In this work, we reframe a stochastic branching
process analysis in terms of probability generating functions and compare it to
continuous time epidemic simulations on networks. In doing so, we predict the
diversity of emerging epidemic courses on both homogeneous and heterogeneous
networks. We show how the challenge of inferring the early course of an
epidemic falls on the randomness of disease spread more so than on the
heterogeneity of contact patterns. We provide an analysis which helps quantify,
in real time, the probability that an epidemic goes supercritical or
conversely, dies stochastically. These probabilities are often assumed to be
one and zero, respectively, if the basic reproduction number, or R0, is greater
than 1, ignoring the heterogeneity and randomness inherent to disease spread.
This framework can give more insight into early epidemic spread by weighting
standard deterministic models with likelihood to inform pandemic preparedness
with probabilistic forecasts.


  Whole genome sequencing (WGS) is quickly becoming the customary means for
identification of antimicrobial resistance (AMR) due to its ability to obtain
high resolution information about the genes and mechanisms that are causing
resistance and driving pathogen mobility. By contrast, traditional phenotypic
(antibiogram) testing cannot easily elucidate such information. Yet development
of AMR prediction tools from genotype-phenotype data can be biased, since
sampling is non-randomized. Sample provenience, period of collection, and
species representation can confound the association of genetic traits with AMR.
Thus, prediction models can perform poorly on new data with sampling
distribution shifts. In this work -- under an explicit set of causal
assumptions -- we evaluate the effectiveness of propensity-based rebalancing
and confounding adjustment on AMR prediction using genotype-phenotype AMR data
from the Pathosystems Resource Integration Center (PATRIC). We select bacterial
genotypes (encoded as k-mer signatures, i.e. DNA fragments of length k),
country, year, species, and AMR phenotypes for the tetracycline drug class,
preparing test data with recent genomes coming from a single country. We test
boosted logistic regression (BLR) and random forests (RF) with/without
bias-handling. On 10,936 instances, we find evidence of species, location and
year imbalance with respect to the AMR phenotype. The crude versus
bias-adjusted change in effect of genetic signatures on AMR varies but only
moderately (selecting the top 20,000 out of 40+ million k-mers). The area under
the receiver operating characteristic (AUROC) of the RF (0.95) is comparable to
that of BLR (0.94) on both out-of-bag samples from bootstrap and the external
test (n=1,085), where AUROCs do not decrease. We observe a 1%-5% gain in AUROC
with bias-handling compared to the sole use of genetic signatures. ...


  How information is encoded in bio-molecular sequences is difficult to
quantify since such an analysis usually requires sampling an exponentially
large genetic space. Here we show how information theory reveals both robust
and compressed encodings in the largest complete genotype-phenotype map (over 5
trillion sequences) obtained to date.


  Phenotype transition takes place in many biological processes such as
differentiation, and understanding how a cell reprograms its global gene
expression profile is a problem of rate theories. A cell phenotype transition
accompanies with switching of expression rates of clusters of genes, analogous
to domain flipping in an Ising system. Here through analyzing single cell RNA
sequencing data in the framework of transition path theory, we set to study how
such a genome-wide expression program switching proceeds in three different
cell transition processes. For each process after reconstructing a Markov
transition model in the cell state space, we formed an ensemble of shortest
paths connecting the initial and final cell states, reconstructed a reaction
coordinate describing the transition progression, and inferred the gene
regulation network (GRN) along the reaction coordinate. In all three processes
we observed common pattern that the frustration of gene regulatory network
(GRN), defined as overall confliction between the regulation received by genes
and their expression states, first increases then decreases when approaching a
new phenotype. The results support a mechanism of concerted silencing of genes
that are active in the initial phenotype and activation of genes that are
active in the final phenotype.


  In the machine learning culture, an independent test set is required for
proper model verification. Failures in model verification, including test set
omission and test set leakage, make it impossible to know whether or not a
trained model is fit for purpose. In this article, we present a systematic
review and quantitative analysis of human gut microbiome classification
studies, conducted to measure the frequency and impact of test set omission and
test set leakage on area under the receiver operating curve (AUC) reporting.
Among 102 articles included for analysis, we find that only 12% of studies
report a bona fide test set AUC, meaning that the published AUCs for 88% of
studies cannot be trusted at face value. Our findings cast serious doubt on the
general validity of research claiming that the gut microbiome has high
diagnostic or prognostic potential in human disease.


  Background:Typically, proteins perform key biological functions by
interacting with each other. As a consequence, predicting which protein pairs
interact is a fundamental problem. Experimental methods are slow, expensive,
and may be error prone.Many computational methods have been proposed to
identify candidate interacting pairs. When accurate, they can serve as an
inexpensive, preliminary filtering stage, to be followed by downstream
experimental validation. Among such methods, sequence-based ones are very
promising.Results:We present, a new algorithm that leverages both topological
and biological information to predict protein-protein interactions. We
comprehensively compare our Framework with state-of-the-art approaches on
reliable PPIs datasets, showing that they have competitive or higher accuracy
on biologically validated test sets.Conclusion:We shown that topological plus
sequence-based computational methods can effectively predict the entire human
interactome compared with methods that leverage only one source of biological
information.


  Mechanical loading generally weakens adhesive structures and eventually leads
to their rupture. However, biological systems can adapt to loads by
strengthening adhesions, which is essential for maintaining the integrity of
tissue and whole organisms. Inspired by cellular focal adhesions, we suggest
here a generic, molecular mechanism that allows adhesion systems to harness
applied loads for self-stabilization under non-equilibrium conditions --
without any active feedback involved. The mechanism is based on conformation
changes of adhesion molecules that are dynamically exchanged with a reservoir.
Tangential loading drives the occupation of some stretched conformation states
out of equilibrium, which, for thermodynamic reasons, leads to association of
further molecules with the adhesion cluster. Self-stabilization robustly
increases adhesion lifetimes in broad parameter ranges. Unlike for catch-bonds,
bond dissociation rates do not decrease with force. The self-stabilization
principle can be realized in many ways in complex adhesion-state networks; we
show how it naturally occurs in cellular adhesions involving the adaptor
proteins talin and vinculin.


  Allostery is a form of protein regulation, where ligands that bind sites
located apart from the active site can modify the activity of the protein. The
molecular mechanisms of allostery have been extensively studied, because
allosteric sites are less conserved than active sites, and drugs targeting them
are more specific than drugs binding the active sites. Here we quantify the
importance of allostery in genetic disease. We show that 1) known allosteric
proteins are central in disease networks, and contribute to genetic disease and
comorbidities much more than non-allosteric proteins, in many major disease
types like hematopoietic diseases, cardiovascular diseases, cancers, diabetes,
or diseases of the central nervous system. 2) variants from cancer genome-wide
association studies are enriched near allosteric proteins, indicating their
importance to polygenic traits; and 3) the importance of allosteric proteins in
disease is due, at least partly, to their central positions in protein-protein
interaction networks, and probably not due to their dynamical properties.


  There is great interest to develop artificial intelligence-based
protein-ligand affinity models due to their immense applications in drug
discovery. In this paper, PointNet and PointTransformer, two pointwise
multi-layer perceptrons have been applied for protein-ligand affinity
prediction for the first time. Three-dimensional point clouds could be rapidly
generated from the data sets in PDBbind-2016, which contain 3 772 and 11 327
individual point clouds derived from the refined or/and general sets,
respectively. These point clouds were used to train PointNet or
PointTransformer, resulting in protein-ligand affinity prediction models with
Pearson correlation coefficients R = 0.831 or 0.859 from the larger point
clouds respectively, based on the CASF-2016 benchmark test. The analysis of the
parameters suggests that the two deep learning models were capable to learn
many interactions between proteins and their ligands, and these key atoms for
the interaction could be visualized in point clouds. The protein-ligand
interaction features learned by PointTransformer could be further adapted for
the XGBoost-based machine learning algorithm, resulting in prediction models
with an average Rp of 0.831, which is on par with the state-of-the-art machine
learning models based on PDBbind database. These results suggest that point
clouds derived from the PDBbind datasets are useful to evaluate the performance
of 3D point clouds-centered deep learning algorithms, which could learn
critical protein-ligand interactions from natural evolution or medicinal
chemistry and have wide applications in studying protein-ligand interactions.


  We investigate the impact of parity on the abundance of weak species in the
context of the simplest generalization of the rock-paper-scissors model to an
arbitrary number of species -- we consider models with a total number of
species ($N_S$) between 3 and 12, having one or more (weak) species
characterized by a reduced predation probability (by a factor of ${\mathcal
P}_w$ with respect to the other species). We show, using lattice based spatial
stochastic simulations with random initial conditions, large enough for
coexistence to prevail, that parity effects are significant. We find that the
performance of weak species is dependent on whether the total number of species
is even or odd, especially for $N_S \le 8$, with odd numbers of species being
on average more favourable to weak species than even ones. We further show
that, despite the significant dispersion observed among individual models, a
weak species has on average a higher abundance than a strong one if ${\mathcal
P}_w$ is sufficiently smaller than unity -- the notable exception being the
four species case.


  The amount and variety of data is increasing drastically for several years.
These data are often represented as networks, which are then explored with
approaches arising from network theory. Recent years have witnessed the
extension of network exploration methods to leverage more complex and richer
network frameworks. Random walks, for instance, have been extended to explore
multilayer networks. However, current random walk approaches are limited in the
combination and heterogeneity of network layers they can handle. New analytical
and numerical random walk methods are needed to cope with the increasing
diversity and complexity of multilayer networks. We propose here MultiXrank, a
Python package that enables Random Walk with Restart (RWR) on any kind of
multilayer network with an optimized implementation. This package is supported
by a universal mathematical formulation of the RWR. We evaluated MultiXrank
with leave-one-out cross-validation and link prediction, and introduced
protocols to measure the impact of the addition or removal of multilayer
network data on prediction performances. We further measured the sensitivity of
MultiXrank to input parameters by in-depth exploration of the parameter space.
Finally, we illustrate the versatility of MultiXrank with different use-cases
of unsupervised node prioritization and supervised classification in the
context of human genetic diseases.


  The Fisher-KPP model, and generalisations thereof, is a simple
reaction-diffusion models of biological invasion that assumes individuals in
the population undergo linear diffusion with diffusivity $D$, and logistic
proliferation with rate $\lambda$. Biologically-relevant initial conditions
lead to long-time travelling wave solutions that move with speed
$c=2\sqrt{\lambda D}$. Despite these attractive features, there are several
biological limitations of travelling wave solutions of the Fisher-KPP model.
First, these travelling wave solutions do not predict a well-defined invasion
front. Second, biologically-relevant initial conditions lead to travelling
waves that move with speed $c=2\sqrt{\lambda D} > 0$. This means that, for
biologically-relevant initial data, the Fisher-KPP model can not be used to
study invasion with $c \ne 2\sqrt{\lambda D}$, or retreating travelling waves
with $c < 0$. Here, we reformulate the Fisher-KPP model as a moving boundary
problem on $x < s(t)$, and we show that this reformulated model alleviates the
key limitations of the Fisher-KPP model. Travelling wave solutions of the
moving boundary problem predict a well-defined front, and can propagate with
any wave speed, $-\infty < c < \infty$. Here, we establish these results using
a combination of high-accuracy numerical simulations of the time-dependent
partial differential equation, phase plane analysis and perturbation methods.
All software required to replicate this work is available on GitHub.


  Axonal growth and guidance at the ventral floor plate is here followed
$\textit{in vivo}$ in real time at high resolution by light-sheet microscopy
along several hundred micrometers of the zebrafish spinal cord. The recordings
show the strikingly stereotyped spatio-temporal control that governs midline
crossing. Commissural axons are observed crossing the ventral floor plate
midline perpendicularly at about 20 microns/h, in a manner dependent on the
Robo3 receptor and with a growth rate minimum around the midline, confirming
previous observations. At guidance points, commissural axons are seen to
decrease their growth rate and growth cones increase in size. Commissural
filopodia appear to interact with the nascent neural network, and thereby
trigger immediate plastic and reversible sinusoidal-shaped bending movements of
neighboring commissural shafts. Ipsilateral axons extend concurrently, but
straight and without bends, at three to six times higher growth rates than
commissurals, indicating they project their path on a substrate-bound surface
rather than relying on diffusible guidance cues. Growing axons appeared to be
under stretch, an observation that is of relevance for tension-based models of
cortical morphogenesis. The \textit{in vivo} observations provide for a
discussion of the current distinction between substrate-bound and diffusible
guidance cues. The study applies the transparent zebrafish model that provides
an experimental model system to explore further the cellular, molecular and
physical mechanisms involved during axonal growth, guidance and midline
crossing through a combination of $\textit{in vitro}$ and $\textit{in vivo}$
approaches.


  Recent studies have shown that the majority of published computational models
in systems biology and physiology are not repeatable or reproducible. There are
a variety of reasons for this. One of the most likely reasons is that given how
busy modern researchers are and the fact that no credit is given to authors for
publishing repeatable work, it is inevitable that this will be the case. The
situation can only be rectified when government agencies, universities and
other research institutions change policies and that journals begin to insist
that published work is in fact at least repeatable if not reproducible. In this
chapter guidelines are described that can be used by researchers to help make
sure their work is repeatable. A scoring system is suggested that authors can
use to determine how well they are doing.


  Microbiome studies have recently transitioned from experimental designs with
a few hundred samples to designs spanning tens of thousands of samples. Modern
studies such as the Earth Microbiome Project (EMP) afford the statistics
crucial for untangling the many factors that influence microbial community
composition. Analyzing those data used to require access to a compute cluster,
making it both expensive and inconvenient. We show that recent improvements in
both hardware and software now allow to compute key bioinformatics tasks on
EMP-sized data in minutes using a gaming-class laptop, enabling much faster and
broader microbiome science insights.


  Numerous real-world systems can be naturally modeled as multilayer networks,
enabling an efficient way to characterize those complex systems. Much evidence
in the context of system biology indicated that the collections between
different molecular networks can dramatically impact the global network
functions. Here, we focus on the molecular multiplex networks coupled by the
transcriptional regulatory network (TRN) and protein-protein interaction (PPI)
network, exploring the controllability and energy requiring in these types of
molecular multiplex networks. We find that the driver nodes tend to avoid
essential or pathogen-related genes. Yet, imposing the external inputs to these
essential or pathogen-related genes can remarkably reduce the energy cost,
implying their crucial role in network control. Moreover, we find that lower
minimal driver nodes as well as energy requiring are associated with
disassortative coupling between TRN and PPI networks. Our findings in several
species provide comprehensive understanding of genes' roles in biology and
network control.


  Renewal equations are a popular approach used in modelling the number of new
infections, i.e., incidence, in an outbreak. We develop a stochastic model of
an outbreak based on a time-varying variant of the Crump-Mode-Jagers branching
process. This model accommodates a time-varying reproduction number and a
time-varying distribution for the generation interval. We then derive
renewal-like integral equations for incidence, cumulative incidence and
prevalence under this model. We show that the equations for incidence and
prevalence are consistent with the so-called back-calculation relationship. We
analyse two particular cases of these integral equations, one that arises from
a Bellman-Harris process and one that arises from an inhomogeneous Poisson
process model of transmission. We also show that the incidence integral
equations that arise from both of these specific models agree with the renewal
equation used ubiquitously in infectious disease modelling. We present a
numerical discretisation scheme to solve these equations, and use this scheme
to estimate rates of transmission from serological prevalence of SARS-CoV-2 in
the UK and historical incidence data on Influenza, Measles, SARS and Smallpox.


  Suppose we have $n$ different types of self-replicating entity, with the
population $P_i$ of the $i$th type changing at a rate equal to $P_i$ times the
fitness $f_i$ of that type. Suppose the fitness $f_i$ is any continuous
function of all the populations $P_1, \dots, P_n$. Let $p_i$ be the fraction of
replicators that are of the $i$th type. Then $p = (p_1, \dots, p_n)$ is a
time-dependent probability distribution, and we prove that its speed as
measured by the Fisher information metric equals the variance in fitness. In
rough terms, this says that the speed at which information is updated through
natural selection equals the variance in fitness. This result can be seen as a
modified version of Fisher's fundamental theorem of natural selection. We
compare it to Fisher's original result as interpreted by Price, Ewens and
Edwards.


  Dengue represents a serious public health problem in tropical and subtropical
regions worldwide. The number of dengue cases and its geographical expansion
has increased in recent decades, driven mostly after by social and
environmental factors. In Costa Rica, it has been endemic since it was first
introduced in 1993. In this article, wavelet analyzes (wavelet power spectrum
and wavelet coherence) were performed to detect and quantify dengue periodicity
and describe patterns of synchrony between dengue incidence and climatic and
environmental factors: Normalized Difference Water Index, Enhanced Vegetation
Index, Normalized Difference Vegetation Index, Tropical North Atlantic indices,
Land Surface Temperature, and El Ni\~no Southern Oscillation indices in 32
different cantons, using dengue surveillance from 2000 to 2019. Results showed
that the dengue dominant cycles are in periods of 1, 2, and 3 years. The
wavelet coherence analysis showed that the vegetation indices are correlated
with dengue incidence in places located in the central and Northern Pacific of
the country in the period of 1 year. Climatic variables such as El Ni\~no 3,
3.4, 4, showed a strong correlation with dengue incidence in the period of 3
years and the Tropical North Atlantic is correlated with dengue incidence in
the period of 1 year. Land Surface Temperature showed a strong correlation with
dengue time series in the 32 cantons.


  In this paper, we study analytically the statistics of the number of
equilibria in pairwise social dilemma evolutionary games with mutation where a
game's payoff entries are random variables. Using the replicator-mutator
equations, we provide explicit formulas for the probability distributions of
the number of equilibria as well as other statistical quantities. This analysis
is highly relevant assuming that one might know the nature of a social dilemma
game at hand (e.g., cooperation vs coordination vs anti-coordination), but
measuring the exact values of its payoff entries is difficult. Our delicate
analysis shows clearly the influence of the mutation probability on these
probability distributions, providing insights into how varying this important
factor impacts the overall behavioural or biological diversity of the
underlying evolutionary systems.


  We use mobile device data to construct empirical interpersonal physical
contact networks in the city of Portland, Oregon, both before and after social
distancing measures were enacted during the COVID-19 pandemic. These networks
reveal how social distancing measures and the public's reaction to the
incipient pandemic affected the connectivity patterns within the city. We find
that as the pandemic developed there was a substantial decrease in the number
of individuals with many contacts. We further study the impact of these
different network topologies on the spread of COVID-19 by simulating an SEIR
epidemic model over these networks, and find that the reduced connectivity
greatly suppressed the epidemic. We then investigate how the epidemic responds
when part of the population is vaccinated, and we compare two vaccination
distribution strategies, both with and without social distancing. Our main
result is that the heavy-tailed degree distribution of the contact networks
causes a targeted vaccination strategy that prioritizes high-contact
individuals to reduce the number of cases far more effectively than a strategy
that vaccinates individuals at random. Combining both targeted vaccination and
social distancing leads to the greatest reduction in cases, and we also find
that the marginal benefit of a targeted strategy as compared to a random
strategy exceeds the marginal benefit of social distancing for reducing the
number of cases. These results have important implications for ongoing vaccine
distribution efforts worldwide.


  Many proteins undergo conformational changes during their activity. A full
understanding of the function of these proteins can only be obtained if
different conformations and transitions between them can be monitored in
aqueous solution, with adequate temporal resolution and, ideally, on a
single-molecule level. Interrogating conformational dynamics of single proteins
remains, however, exquisitely challenging and typically requires site-directed
chemical modification combined with rigorous minimization of possible
artifacts. These obstacles limit the number of single-protein investigations.
The work presented here introduces an approach that traps single unmodified
proteins from solution in a plasmonic hotspot and makes it possible to assign
changes in refractive index to changes in protein conformation while monitoring
these changes for minutes to hours with a temporal resolution at least as fast
as 40 microseconds. The resulting single molecule data reveals that adenylate
kinase employs a hidden enzymatic sub-cycle during catalysis, that citrate
synthase populates a previously unknown intermediate conformation, which is
more important for its enzymatic activity than its well-known open
conformation, that hemoglobin transitions in several steps from its
deoxygenated and rigid T state to its oxygenated and flexible R state, and that
apo-calmodulin thermally unfolds and refolds in steps that correspond to
conformational changes of individual protein domains.


  We investigate the distribution of numbers of secondary cases in households
in the Office for National Statistics COVID-19 Infection Survey (ONS CIS),
stratified by timing of vaccination and infection in the households. This shows
a total effect of a statistically significant approximate halving of the
secondary attack rate in households following vaccination.


  An outbreak of the Delta (B.1.617.2) variant of SARS-CoV-2 that began around
mid-June 2021 in Sydney, Australia, quickly developed into a nation-wide
epidemic. The ongoing epidemic is of major concern as the Delta variant is more
infectious than previous variants that circulated in Australia in 2020. Using a
re-calibrated agent-based model, we explored a feasible range of
non-pharmaceutical interventions, including case isolation, home quarantine,
school closures, and stay-at-home restrictions (i.e., "social distancing"). Our
modelling indicated that the levels of reduced interactions in workplaces and
across communities attained in Sydney and other parts of the nation were
inadequate for controlling the outbreak. A counter-factual analysis suggested
that if 70% of the population followed tight stay-at-home restrictions, then at
least 45 days would have been needed for new daily cases to fall from their
peak to below ten per day. Our model predicted that, under a progressive
vaccination rollout, if 40-50% of the Australian population follow stay-at-home
restrictions, the incidence will peak by mid-October 2021: the peak in
incidence across the nation was indeed observed in mid-October. We also
quantified an expected burden on the healthcare system and potential fatalities
across Australia.


  Non-pharmacologic interventions (NPIs) are one method to mitigate the spread
and effects of the COVID-19 pandemic in the United States. NPIs promote
protective actions to reduce exposure risk and can reduce mobility patterns
within communities. Growing research literature suggests that socially
vulnerable populations are disproportionately impacted with higher infection
and higher fatality rates of COVID-19, though there is limited understanding of
the underlying mechanisms to this health disparity. Thus, the research examines
two distinct and complimentary datasets at a granular scale for five urban
locations. Through statistical and spatial analyses, the research extensively
investigates the exposure risk reduction of socially vulnerable populations due
to NPIs. The mobility dataset tracks population movement across ZIP codes; it
is used for an origin-destination network analysis. The population activity
dataset is based on the number of visits from census block groups (CBG) to
points of interest (POIs), such as grocery stores, restaurants, education
centers, and medical facilities; it is used for network analysis of
population-facilities interactions. The mobility dataset showed that, after the
implementation of NPIs, socially vulnerable populations engaged in increased
mobility in the form of inflow between ZIP code areas. Similarly, population
activity analysis showed an increased exposure risk for socially vulnerable
populations based on a greater number of inflow visits of CBGs to POIs, which
increases the risk of contact at POIs, and a greater number of outflow visits
from POIs to home CBGs, which increases risk of transmission within CBGs. These
findings can assist emergency planners and public health officials in
comprehending how different groups are able to implement protective actions and
can inform more equitable and data-driven NPI policies for future epidemics.


  Cryo-electron tomography (cryo-ET) is uniquely suited to precisely localize
macromolecular complexes in situ, that is in a close-to-native state within
their cellular compartments, in three-dimensions at high resolution. Point
pattern analysis (PPA) allows quantitative characterization of the spatial
organization of particles. However, current implementations of PPA functions
are not suitable for applications to cryo-ET data because they do not consider
the real, typically irregular 3D shape of cellular compartments and molecular
complexes. Here, we designed and implemented first and the second-order, uni-
and bivariate PPA functions in a Python package for statistical spatial
analysis of particles located in three dimensional regions of arbitrary shape,
such as those encountered in cellular cryo-ET imaging (PyOrg).
  To validate the implemented functions, we applied them to specially designed
synthetic datasets. This allowed us to find the algorithmic solutions that
provide the best accuracy and computational performance, and to evaluate the
precision of the implemented functions. Applications to experimental data
showed that despite the higher computational demand, the use of the
second-order functions is advantageous to the first-order ones, because they
allow characterization of the particle organization and statistical inference
over a range of distance scales, as well as the comparative analysis between
experimental groups comprising multiple tomograms.
  Altogether, PyOrg is a versatile, precise, and efficient open-source software
for reliable quantitative characterization of macromolecular organization
within cellular compartments imaged in situ by cryo-ET, as well as to other 3D
imaging systems where real-size particles are located within regions possessing
complex geometry.


  The Leloir galactose utilization or GAL pathway of budding yeasts, including
that of the baker's yeast Saccharomyces cerevisiae and the opportunistic human
pathogen Candida albicans, breaks down the sugar galactose for energy and
biomass production. The GAL pathway has long served as a model system for
understanding how eukaryotic metabolic pathways, including their modes of
regulation, evolve. More recently, the physical linkage of the structural genes
GAL1, GAL7, and GAL10 in diverse budding yeast genomes has been used as a model
for understanding the evolution of gene clustering. In this opinion, we
summarize exciting recent work on three different aspects of this iconic
pathway's evolution: gene cluster organization, GAL gene regulation, and the
population genetics of the GAL pathway.


  Spike proteins, 1200 amino acids, are divided into two nearly equal parts, S1
and S2. We review here phase transition theory, implemented quantitatively by
thermodynamic scaling. The theory explains the evolution of Coronavirus
extremely high contagiousness caused by a few mutations from CoV2003 to CoV2019
identified among hundreds in S1. The theory previously predicted the
unprecedented success of spike-based vaccines. Here we analyze impressive
successes by McClellan et al., 2020, in stabilizing their original S2P vaccine
to Hexapro. Hexapro has expanded the two proline mutations of S2P, 2017, to six
combined proline mutations in S2. Their four new mutations are the result of
surveying 100 possibilities in their detailed structure-based context Our
analysis, based on only sparse publicly available data, suggests new proline
mutations could improve the Hexapro combination to Octapro or beyond.


  Many experiments show that protein condensates formed by liquid-liquid phase
separation exhibit aging rheological properties. Quantitatively, recent
experiments by Jawerth et al. (Science 370, 1317, 2020) show that protein
condensates behave as aging Maxwell fluids with an increasing relaxation time
as the condensates age. Despite the universality of this aging phenomenon, a
theoretical understanding of this aging behavior is lacking. In this work, we
propose a mesoscopic model of protein condensates in which a phase transition
from aging phase to non-aging phase occurs as the control parameter changes,
such as temperature. The model predicts that protein condensates behave as
viscoelastic Maxwell fluids at all ages, with the macroscopic viscosity
increasing over time. The model also predicts that protein condensates are
non-Newtonian fluids under a constant shear rate with the shear stress
increasing over time. Our model successfully explains multiple existing
experimental observations and also makes general predictions that are
experimentally testable.


  Rapid advance of experimental techniques provides an unprecedented in-depth
view into complex developmental processes. Still, little is known on how the
complexity of multicellular organisms evolved by elaborating developmental
programs and inventing new cell types. A hurdle to understanding developmental
evolution is the difficulty of even describing the intertwined network of
spatiotemporal processes underlying the development of complex multicellular
organisms. Nonetheless, an overview of developmental trajectories can be
obtained from cell type lineage maps. Here, we propose that these lineage maps
can also reveal how developmental programs evolve: the modes of evolving new
cell types in an organism should be visible in its developmental trajectories,
and therefore in the geometry of its cell type lineage map. This idea is
demonstrated using a parsimonious generative model of developmental programs,
which allows us to reliably survey the universe of all possible programs and
examine their topological features. We find that, contrary to belief, tree-like
lineage maps are rare and lineage maps of complex multicellular organisms are
likely to be directed acyclic graphs where multiple developmental routes can
converge on the same cell type. While cell type evolution prescribes what
developmental programs come into existence, natural selection prunes those
programs which produce low-functioning organisms. Our model indicates that
additionally, lineage map topologies are correlated with such a functional
property: the ability of organisms to regenerate.


  Many factors influence biomolecules binding, and its assessment constitutes
an elusive challenge in computational structural biology. In this respect, the
evaluation of shape complementarity at molecular interfaces is one of the main
factors to be considered. We focus on the particular case of antibody-antigen
complexes to quantify the complementarities occurring at molecular interfaces.
We relied on a method we recently developed, which employs the 2D Zernike
descriptors, to characterize investigated regions with an ordered set of
numbers summarizing the local shape properties. Collected a structural dataset
of antibody-antigen complexes, we applied this method and we statistically
distinguished, in terms of shape complementarity, pairs of interacting regions
from non-interacting ones. Thus, we set up a novel computational strategy based
on \textit{in-silico} mutagenesis of antibody binding site residues. We
developed a Monte Carlo procedure to increase the shape complementarity between
the antibody paratope and a given epitope on a target protein surface. We
applied our protocol against several molecular targets in SARS-CoV-2 spike
protein, known to be indispensable for viral cell invasion. We, therefore,
optimized the shape of template antibodies for the interaction with such
regions. As the last step of our procedure, we performed an independent
molecular docking validation of the results of our Monte Carlo simulations.


  Ecology and evolution are inherently linked, and studying a mathematical
model that considers both holds promise of insightful discoveries related to
the dynamics of cooperation. In the present article, we use the prisoner's
dilemma (PD) game as a basis for long-term apprehension of the essential social
dilemma related to cooperation among unrelated individuals. We upgrade the
contemporary PD game with an inclusion of evolution-induced act of punishment
as a third competing strategy in addition to the traditional cooperators and
defectors. In a population structure, the abundance of ecologically-viable free
space often regulates the reproductive opportunities of the constituents.
Hence, additionally, we consider the availability of free space as an
ecological footprint, thus arriving at a simple eco-evolutionary model, which
displays fascinating complex dynamics. As possible outcomes, we report the
individual dominance of cooperators and defectors as well as a plethora of
mixed states, where different strategies coexist followed by maintaining the
diversity in a socio-ecological framework. These states can either be steady or
oscillating, whereby oscillations are sustained by cyclic dominance among
different combinations of cooperators, defectors, and punishers. We also
observe a novel route to cyclic dominance where cooperators, punishers, and
defectors enter a coexistence via an inverse Hopf bifurcation that is followed
by an inverse period doubling route.


  We study the dynamics of membrane vesicle motor transport into dendritic
spines, which are bulbous intracellular compartments in neurons that play a key
role in transmitting signals between neurons. We consider the stochastic analog
of the vesicle transport model in [Park and Fai, The Dynamics of Vesicles
Driven Into Closed Constrictions by Molecular Motors. Bull. Math. Biol. 82, 141
(2020)]. The stochastic version, which may be considered as an agent-based
model, relies mostly on the action of individual myosin motors to produce
vesicle motion. To aid in our analysis, we coarse-grain this agent-based model
using a master equation combined with a partial differential equation
describing the probability of local motor positions. We confirm through
convergence studies that the coarse-graining captures the essential features of
bistability in velocity (observed in experiments) and waiting-time
distributions to switch between steady-state velocities. Interestingly, these
results allow us to reformulate the translocation problem in terms of the mean
first passage time for a run-and-tumble particle moving on a finite domain with
absorbing boundaries at the two ends. We conclude by presenting numerical and
analytical calculations of vesicle translocation.


  In the study of rooted phylogenetic networks, analyzing the set of rooted
phylogenetic trees that are embedded in such a network is a recurring task.
From an algorithmic viewpoint, this analysis almost always requires an
exhaustive search of a particular multiset $S$ of rooted phylogenetic trees
that are embedded in a rooted phylogenetic network $\mathcal{N}$. Since the
size of $S$ is exponential in the number of reticulations of $\mathcal{N}$, it
is consequently of interest to keep this number as small as possible but
without loosing any element of $S$. In this paper, we take a first step towards
this goal by introducing the notion of a non-essential arc of $\mathcal{N}$,
which is an arc whose deletion from $\mathcal{N}$ results in a rooted
phylogenetic network $\mathcal{N}'$ such that the sets of rooted phylogenetic
trees that are embedded in $\mathcal{N}$ and $\mathcal{N}'$ are the same. We
investigate the popular class of tree-child networks and characterize which
arcs are non-essential. This characterization is based on a family of directed
graphs. Using this novel characterization, we show that identifying and
deleting all non-essential arcs in a tree-child network takes time that is
cubic in the number of leaves of the network. Moreover, we show that deciding
if a given arc of an arbitrary phylogenetic network is non-essential is
$\Pi_2^P$-complete.


  A major open problem in biophysics is to understand the highly heterogeneous
transport of many structures inside living cells, such as endosomes. We find
that mathematically it is described by spatio-temporal heterogeneous fractional
Brownian motion (hFBM) which is defined as FBM with a randomly switching
anomalous exponent and random generalized diffusion coefficient. Using a
comprehensive local analysis of a large ensemble of experimental endosome
trajectories (> 10^5), we show that their motion is characterized by power-law
probability distributions of displacements and displacement increments,
exponential probability distributions of local anomalous exponents and
power-law probability distributions of local generalized diffusion coefficients
of endosomes which are crucial ingredients of spatio-temporal hFBM. The
increased sensitivity of deep learning neural networks for FBM characterisation
corroborates the development of this multi-fractal analysis. Our findings are
an important step in understanding endosome transport. We also provide a
powerful tool for studying other heterogeneous cellular processes.


  Phylogenetic Diversity (PD) is a prominent quantitative measure of the
biodiversity of a collection of present-day species (taxa). This measure is
based on the evolutionary distance among the species in the collection. Loosely
speaking, if $\mathcal{T}$ is a rooted phylogenetic tree whose leaf set $X$
represents a set of species and whose edges have real-valued lengths (weights),
then the PD score of a subset $S$ of $X$ is the sum of the weights of the edges
of the minimal subtree of $\mathcal{T}$ connecting the species in $S$. In this
paper, we define several natural variants of the PD score for a subset of taxa
which are related by a known rooted phylogenetic network. Under these variants,
we explore, for a positive integer $k$, the computational complexity of
determining the maximum PD score over all subsets of taxa of size $k$ when the
input is restricted to different classes of rooted phylogenetic networks


  Alterations in the human gut bacteriome can be associated with human health
issues, such as type-2 diabetes and cardiovascular disease. Both external and
internal factors can drive changes in the composition and in the interactions
of the human gut bacteriome, impacting negatively on the host cells. In this
paper, we focus on the human gut bacteriome metabolism and we propose a
two-layer network system to investigate its dynamics. Furthermore, we develop
an in-silico simulation model (virtual GB), allowing us to study the impact of
the metabolite exchange through molecular communications in the human gut
bacteriome network system. Our results show that the regulation of molecular
inputs can strongly affect bacterial population growth and create an unbalanced
network, as shown by the shift in the node weights based on the molecular
signals that are produced. Additionally, we show that the metabolite molecular
communication production is greatly affected when directly manipulating the
composition of the human gut bacteriome network in the virtual GB. These
results indicate that our human GB interaction model can help to identify
hidden behaviors of the human gut bacteriome depending on the molecular signal
interactions. Moreover, the virtual GB can support the research and development
of novel medical treatments based on the accurate control of bacterial growth
and exchange of metabolites.


  The far-reaching consequences of ecological interactions in the dynamics of
biological communities remain an intriguing subject. For decades, competition
has been a cornerstone in ecological processes, but mounting evidence shows
that cooperation does also contribute to the structure of biological
communities. Here, we propose a simple deterministic model for the study of the
effects of facilitation and competition in the dynamics of such systems. The
simultaneous inclusion of both effects produces rich dynamics and captures the
context-dependence observed in the formation of ecological communities. The
approach reproduces relevant aspects of primary and secondary plant succession,
the effect invasive species, and the survival of rare species. The model also
takes into account the role of the ecological priority effect and stress the
crucial role of facilitation in conservation efforts and species coexistence.


  The deviation between chronological age and age predicted from neuroimaging
data has been identified as a sensitive risk-marker of cross-disorder brain
changes, growing into a cornerstone of biological age-research. However,
Machine Learning models underlying the field do not consider uncertainty,
thereby confounding results with training data density and variability. Also,
existing models are commonly based on homogeneous training sets, often not
independently validated, and cannot be shared due to data protection issues.
Here, we introduce an uncertainty-aware, shareable, and transparent Monte-Carlo
Dropout Composite-Quantile-Regression (MCCQR) Neural Network trained on
N=10,691 datasets from the German National Cohort. The MCCQR model provides
robust, distribution-free uncertainty quantification in high-dimensional
neuroimaging data, achieving lower error rates compared to existing models
across ten recruitment centers and in three independent validation samples
(N=4,004). In two examples, we demonstrate that it prevents spurious
associations and increases power to detect accelerated brain-aging. We make the
pre-trained model publicly available.


  Freshwater Unionid bivalves have recently faced ecological upheaval through
pollution, barriers to dispersal, human harvesting, and changes in fish-host
prevalence. Currently, over 70% of species are threatened, endangered or
extinct. To characterize the genetic response to these recent selective
pressures, we collected population genetic data for one successful bivalve
species, Megalonaias nervosa. We identify megabase sized regions that are
nearly monomorphic across the population, a signal of strong, recent selection
reshaping genetic diversity. These signatures of selection encompass a total of
73Mb, greater response to selection than is commonly seen in population genetic
models. We observe 102 duplicate genes with high dN/dS on terminal branches
among regions with sweeps, suggesting that gene duplication is a causative
mechanism of recent adaptation in M. nervosa. Genes in sweeps reflect
functional classes known to be important for Unionid survival, including
anticoagulation genes important for fish host parasitization, detox genes,
mitochondria management, and shell formation. We identify selective sweeps in
regions with no known functional impacts, suggesting mechanisms of adaptation
that deserve greater attention in future work on species survival. In contrast,
polymorphic transposable element insertions appear to be detrimental and
underrepresented among regions with sweeps. TE site frequency spectra are
skewed toward singleton variants, and TEs among regions with sweeps are present
only at low frequency. Our work suggests that duplicate genes are an essential
source of genetic novelty that has helped this successful species succeed in
environments where others have struggled. These results suggest that gene
duplications deserve greater attention in non-model population genomics,
especially in species that have recently faced sudden environmental challenges.


  We develop a novel technique to exploit the extensive data sets provided by
underwater neutrino telescopes to gain information on bioluminescence in the
deep sea. The passive nature of the telescopes gives us the unique opportunity
to infer information on bioluminescent organisms without actively interfering
with them. We propose a statistical method that allows us to reconstruct the
light emission of individual organisms, as well as their location and movement.
A mathematical model is built to describe the measurement process of underwater
neutrino telescopes and the signal generation of the biological organisms. The
Metric Gaussian Variational Inference algorithm is used to reconstruct the
model parameters using photon counts recorded by the neutrino detectors. We
apply this method to synthetic data sets and data collected by the ANTARES
neutrino telescope. The telescope is located 40 km off the French coast and
fixed to the sea floor at a depth of 2475 m. The runs with synthetic data
reveal that we can reliably model the emitted bioluminescent flashes of the
organisms. Furthermore, we find that the spatial resolution of the localization
of light sources highly depends on the configuration of the telescope. Precise
measurements of the efficiencies of the detectors and the attenuation length of
the water are crucial to reconstruct the light emission. Finally, the
application to ANTARES data reveals the first precise localizations of
bioluminescent organisms using neutrino telescope data.


  During organismal development, biochemical reaction networks sense and
respond to mechanical forces to coordinate embryonic patterning with embryo
morphogenesis. Factors such as cortical tension, cell density, and matrix
mechanical properties influence differentiation and cell fate decisions by
modulating gene regulatory signaling networks. A major goal in synthetic
development is to construct gene regulatory circuits that program the
patterning and morphogenesis of synthetic multicellular structures. However, in
the synthetic context, little is known regarding how the physical properties of
the growth environment impact the behavior of synthetic gene circuits. Here, we
exploit physical-chemical coupling observed in a synthetic patterning circuit
in order to control the size and spatial distribution of patterned synthetic
cell sheets. We show that cell density attenuates the propagation of signal
between neighboring cells in a multicellular sheet containing a
contact-dependent patterning circuit based on the synNotch signaling system.
Density-dependent attenuation leads to a signal propagation wave that exhibits
distinct qualitative phases of persistent propagation, transient propagation,
and no propagation. Through computational modeling, we demonstrate that cell
growth parameters determine the phase of propagation observed within a growing
cell sheet. Using growth-modulating drugs and spatial density gradients, we
control the size of synNotch-activated cell populations and generate
tissue-scale activation gradients and kinematic waves. Our study reveals that
density-dependent synNotch activity can be exploited to control a synthetic
multicellular patterning circuit. More broadly, we show that synthetic gene
circuits can be critically impacted by their physical context, providing an
alternate means for programming circuit behavior.


  Testing individuals for pathogens can affect the spread of epidemics.
Understanding how individual-level processes of sampling and reporting test
results can affect community- or population-level spread is a dynamical
modeling question. The effect of testing processes on epidemic dynamics depends
on factors underlying implementation, particularly testing intensity and on
whom testing is focused. Here, we use a simple model to explore how the
individual-level effects of testing might directly impact population-level
spread. Our model development was motivated by the COVID-19 epidemic, but has
generic epidemiological and testing structures. To the classic SIR framework we
have added a per capita testing intensity, and compartment-specific testing
weights, which can be adjusted to reflect different testing emphases --
surveillance, diagnosis, or control. We derive an analytic expression for the
relative reduction in the basic reproductive number due to testing,
test-reporting and related isolation behaviours. Intensive testing and fast
test reporting are expected to be beneficial at the community level because
they can provide a rapid assessment of the situation, identify hot spots, and
may enable rapid contact-tracing. Direct effects of fast testing at the
individual level are less clear, and may depend on how individuals' behaviour
is affected by testing information. Our simple model shows that under some
circumstances both increased testing intensity and faster test reporting can
reduce the effectiveness of control, and allows us to explore the conditions
under which this occurs. Conversely, we find that focusing testing on infected
individuals always acts to increase effectiveness of control.


  Maintenance of epidermal thickness is critical to the barrier function of the
skin. Decreased tissue thickness, specifically in the stratum corneum (the
outermost layer of the tissue), causes discomfort and inflammation, and is
related to several severe diseases of the tissue. In order to maintain both
stratum corneum thickness and overall tissue thickness it is necessary for the
system to balance cell proliferation and cell loss. Cell proliferation in the
epidermis occurs in the basal layer and causes constant upwards movement in the
tissue. Cell loss occurs when dead cells at the top of the tissue are lost to
the environment through a process called desquamation. Desquamation is thought
to occur through a gradual reduction in adhesion between cells, due to the
cleaving of adhesion proteins by enzymes, in the stratum corneum.
  In this paper we will investigate combining a (mass action) subcellular model
of desquamation with a three dimensional (cell centre based) multicellular
model of the interfollicular epidermis to better understand maintenance of
epidermal thickness. These investigations show that hypothesised biological
models for the degradation of cell-cell adhesion from the literature are able
to provide a consistent rate of cell loss in the multicellular model. This loss
balances proliferation, and hence maintains a homeostatic tissue thickness.
Moreover, we find that multiple proliferative cell populations in the basal
layer can be represented by a single proliferative cell population, simplifying
investigations with this model.
  The model is used to investigate a disorder (Netherton Syndrome) which
disrupts desquamation. The model shows how biochemical changes can cause
disruptions to the tissue, resulting in a reduced tissue thickness and
consequently diminishing the protective role of the tissue. A hypothetical
treatment result is also investigated. [ABR]


  We present a data set from a first-principles study of amino-methylated and
acetylated (capped) dipeptides of the 20 proteinogenic amino acids - including
alternative possible side chain protonation states and their interactions with
selected divalent cations (Ca$^{2+}$, Mg$^{2+}$ and Ba$^{2+}$). The data covers
21,909 stationary points on the respective potential-energy surfaces in a wide
relative energy range of up to 4 eV (390 kJ/mol). Relevant properties of
interest, like partial charges, were derived for the conformers. The motivation
was to provide a solid data basis for force field parameterization and further
applications like machine learning or benchmarking. In particular the process
of creating all this data on the same first-principles footing, i.e.
density-functional theory calculations employing the generalized gradient
approximation with a van der Waals correction, makes this data suitable for
data-driven force field development. To make the data accessible across domain
borders and to machines, we formalized the metadata in an ontology.


  Cyclic dominance is frequently believed to be a mechanism that maintains
diversity of competing species. But this delicate balance could also be fragile
if some of the members is weakened because an extinction of a species will
involve the annihilation of its predator hence leaving only a single species
alive. To check this expectation we here introduce a fourth species which
chases exclusively a single member of the basic model composed by three
cyclically dominant species. Interestingly, the coexistence is not necessarily
broken and we have detected three consecutive phase transitions as we vary only
the invasion strength of the fourth pestilent species. The resulting phases are
analyzed by different techniques including the study of the Hamming distance
density profiles. Some of our observations strengthen previous findings about
cyclically dominant system, but they also offer new revelations and
counter-intuitive phenomenon, like supporting pestilent species may result in
its extinction, hence enriching our understanding about these very simple but
still surprisingly complex systems.


  Many of the molecular mechanisms underlying the pathological aggregation of
proteins observed in neurodegenerative diseases are still not fully understood.
Among the diseases associated with protein aggregates, for example, Amyotrophic
Lateral Sclerosis (ALS) is of relevant importance. Although understanding the
processes that cause the disease is still an open challenge, its relationship
with protein aggregation is widely known. In particular, human TDP-43, an
RNA/DNA binding protein, is a major component of pathological cytoplasmic
inclusions described in ALS patients. The deposition of the phosphorylated
full-length TDP-43 in spinal cord cells has been widely studied, and it has
been shown that the brain cortex presents an accumulation of phosphorylated
C-terminal fragments (CTFs). Even if it is debated whether CTFs represent a
primary cause of ALS, they are a hallmark of TDP-43 related neurodegeneration
in the brain. Here, we investigate the CTFs aggregation process, providing a
possible computational model of interaction based on the evaluation of shape
complementarity at the interfaces. To this end, extensive Molecular Dynamics
(MD) simulations were conducted for different types of fragments with the aim
of exploring the equilibrium configurations. Adopting a newly developed
approach based on Zernike polynomials, for finding complementary regions of the
molecular surface, we sampled a large set of exposed portions of the molecular
surface of CTFs structures as obtained from MD simulations. The analysis
proposes a set of possible associations between the CTFs, which could drive the
aggregation process of the CTFs.


  The application of deep learning to generative molecule design has shown
early promise for accelerating lead series development. However, questions
remain concerning how factors like training, dataset, and seed bias impact the
technology's utility to medicine and computational chemists. In this work, we
analyze the impact of seed and training bias on the output of an
activity-conditioned graph-based variational autoencoder (VAE). Leveraging a
massive, labeled dataset corresponding to the dopamine D2 receptor, our
graph-based generative model is shown to excel in producing desired conditioned
activities and favorable unconditioned physical properties in generated
molecules. We implement an activity swapping method that allows for the
activation, deactivation, or retention of activity of molecular seeds, and we
apply independent deep learning classifiers to verify the generative results.
Overall, we uncover relationships between noise, molecular seeds, and training
set selection across a range of latent-space sampling procedures, providing
important insights for practical AI-driven molecule generation.


  We develop a method for analyzing spatial and spatiotemporal anomalies in
geospatial data using topological data analysis (TDA). To do this, we use
persistent homology (PH), which allows one to algorithmically detect geometric
voids in a data set and quantify the persistence of such voids. We construct an
efficient filtered simplicial complex (FSC) such that the voids in our FSC are
in one-to-one correspondence with the anomalies. Our approach goes beyond
simply identifying anomalies; it also encodes information about the
relationships between anomalies. We use vineyards, which one can interpret as
time-varying persistence diagrams (which are an approach for visualizing PH),
to track how the locations of the anomalies change with time. We conduct two
case studies using spatially heterogeneous COVID-19 data. First, we examine
vaccination rates in New York City by zip code at a single point in time.
Second, we study a year-long data set of COVID-19 case rates in neighborhoods
of the city of Los Angeles.


  The usual problem for group testing is this: For a given number of
individuals and a given prevalence, how many tests T* are required to find
every infected individual? In real life, however, the problem is usually
different: For a given number of individuals, a given prevalence, and a limited
number of tests T much smaller than T*, how can these tests best be used? In
this conference paper, we outline some recent results on this problem for two
models. First, the "practical" model, which is relevant for screening for
COVID-19 and has tests that are highly specific but imperfectly sensitive,
shows that simple algorithms can be outperformed at low prevalence and high
sensitivity. Second, the "theoretical" model of very low prevalence with
perfect tests gives interesting new mathematical results.


  Among the diverse and important applications that networks currently have is
the modeling of infectious diseases. Immunization, or the process of protecting
nodes in the network, plays a key role in stopping diseases from spreading.
Hence the importance of having tools or strategies that allow the solving of
this challenge. In this work, we evaluate the effectiveness of the
DIL-W^{\alpha} ranking in immunizing nodes in an edge-weighted network. The
network is obtained from a real database and the spread of COVID-19 was modeled
with the classic SIR model. We apply the protection to the network, according
to the importance ranking list produced by DIL-W^{\alpha}.