  Mumbai, amongst the most densely populated cities in the world, has witnessed
the fourth largest number of cases and the largest number of deaths among all
the cities in India (as of 28th October 2020). Along with the rest of India,
lockdowns (of varying degrees) have been in effect in Mumbai since March 25,
2020. Given the large economic toll on the country from the lockdown and the
related restrictions on mobility of people and goods, swift opening of the
economy especially in a financial hub such as Mumbai becomes critical. In this
report, we use the IISc-TIFR agent based simulator to develop long term
projections for Mumbai under realistic scenarios related to Mumbai's opening of
the workplaces, or equivalently, the economy, and the associated public
transportation through local trains and buses. These projections were developed
taking into account a possible second wave if the economy and the local trains
are fully opened either on November 1, 2020 or on January 1, 2021. The impact
on infection spread in Mumbai if the schools and colleges open on January first
week 2021 is also considered. We also try to account for the increased
intermingling amongst the population during the Ganeshotsav festival as well as
around the Navratri/Dussehra and Diwali festival. Our conclusion, based on our
simulations, is that the impact of fully opening up the economy on November 1
is manageable provided reasonable medical infrastructure is in place. Further,
schools and colleges opening in January do not lead to excessive increase in
infections.
  The report also explores the relative effectiveness of contact tracing vs
containment zones, and also includes very rudimentary results of the effect of
vaccinating the elderly population in February 2021.


  Local ancestry inference (LAI) identifies the ancestry of each segment of an
individual's genome and is an important step in medical and population genetic
studies of diverse cohorts. Several techniques have been used for LAI,
including Hidden Markov Models and Random Forests. Here, we formulate the LAI
task as an image segmentation problem and develop a new LAI tool using a deep
convolutional neural network with an encoder-decoder architecture. We train our
model using complete genome sequences from 982 unadmixed individuals from each
of five continental ancestry groups, and we evaluate it using simulated admixed
data derived from an additional 279 individuals selected from the same
populations. We show that our model is able to learn admixture as a zero-shot
task, yielding ancestry assignments that are nearly as accurate as those from
the existing gold standard tool, RFMix.


  Identifying potential drug targets using metabolic modeling requires
integrating multiple modeling methods and heterogenous biological datasets,
which can be challenging without sophisticated tools. We developed COMO, a
user-friendly pipeline that integrates multi-omics data processing,
context-specific metabolic model development, simulations, drug databases, and
disease data to aid drug discovery. COMO can be installed as a Docker image and
includes intuitive instructions within a Jupyter Lab environment. It provides a
comprehensive solution for multi-omics integration of bulk and single-cell
RNA-seq, microarrays, and proteomics to develop context-specific metabolic
models. Using public databases, open-source solutions for model construction,
and a streamlined approach for predicting repurposable drugs, COMO empowers
researchers to investigate low-cost alternatives and novel disease treatments.
As a case study, we used the pipeline to construct metabolic models of B cells,
which simulate and analyze them to predict 25 and 23 metabolic drug targets for
rheumatoid arthritis and systemic lupus erythematosus, respectively. COMO can
be used to construct models for any cell or tissue type and identify drugs for
any human disease. The pipeline has the potential to improve the health of the
global community cost-effectively by providing high-confidence targets to
pursue in preclinical and clinical studies.


  In this study, we analyze the effectiveness of measures aimed at finding and
isolating infected individuals to contain epidemics like COVID-19, as the
suppression induced over the effective reproduction number. We develop a
mathematical model to compute the relative suppression of the effective
reproduction number of an epidemic that such measures produce. This outcome is
expressed as a function of a small set of parameters that describe the main
features of the epidemic and summarize the effectiveness of the isolation
measures. In particular, we focus on the impact when a fraction of the
population uses a mobile application for epidemic control. Finally, we apply
the model to COVID-19, providing several computations as examples, and a link
to a public repository to run custom calculations. These computations display
in a quantitative manner the importance of recognizing infected individuals
from symptoms and contact-tracing information, and isolating them as early as
possible. The computations also assess the impact of each variable on the
suppression of the epidemic.


  The infection fatality rate (IFR) of the Coronavirus Disease 2019 (COVID-19)
is one of the most discussed figures in the context of this pandemic. Using
German COVID-19 surveillance data and age-group specific IFR estimates from
multiple international studies, this work investigates time-dependent
variations in effective IFR over the course of the pandemic. Three different
methods for estimating (effective) IFRs are presented: (a) population-averaged
IFRs based on the assumption that the infection risk is independent of age and
time, (b) effective IFRs based on the assumption that the age distribution of
confirmed cases approximately reflects the age distribution of infected
individuals, and (c) effective IFRs accounting for age- and time-dependent dark
figures of infections. Results show that effective IFRs in Germany are
estimated to vary over time, as the age distributions of confirmed cases and
estimated infections are changing during the course of the pandemic. In
particular during the first and second waves of infections in spring and
autumn/winter 2020, there has been a pronounced shift in the age distribution
of confirmed cases towards older age groups, resulting in larger effective IFR
estimates. The temporary increase in effective IFR during the first wave is
estimated to be smaller but still remains when adjusting for age- and
time-dependent dark figures. A comparison of effective IFRs with observed CFRs
indicates that a substantial fraction of the time-dependent variability in
observed mortality can be explained by changes in the age distribution of
infections. Furthermore, a vanishing gap between effective IFRs and observed
CFRs is apparent after the first infection wave, while a moderately increasing
gap can be observed during the second wave. Further research is warranted to
obtain timely age-stratified IFR estimates.


  We study Brownian motion in a space with a high density of moving obstacles
in 1, 2 and 3 dimensions. Our tracers diffuse anomalously over many decades in
time, before reaching a diffusive steady state with an effective diffusion
constant $D_\mathrm{eff}$ that depends on the obstacle density and diffusivity.
The scaling of $D_\mathrm{eff}$, above and below a critical regime at the
percolation point for void space, is characterized by two critical exponents:
the conductivity $\mu$, also found in models with frozen obstacles, and $\psi$,
which quantifies the effect of obstacle diffusivity.


  Cultural phylogenies, or "trees" of culture, are typically built using
methods from biology that use similarities and differences in artifacts to
infer the historical relationships between the populations that produced them.
While these methods have yielded important insights, particularly in
linguistics, researchers continue to debate the extent to which cultural
phylogenies are tree-like or reticulated due to high levels of horizontal
transmission. In this study, we propose a novel method for phylogenetic
reconstruction using dynamic community detection that explicitly accounts for
transmission between lineages. We used data from 1,498,483 collaborative
relationships between electronic music artists to construct a cultural
phylogeny based on observed population structure. The results suggest that,
although the phylogeny is fundamentally tree-like, horizontal transmission is
common and populations never become fully isolated from one another. In
addition, we found evidence that electronic music diversity has increased
between 1975 and 1999. The method used in this study is available as a new R
package called DynCommPhylo. Future studies should apply this method to other
cultural systems such as academic publishing and film, as well as biological
systems where high resolution reproductive data is available, to assess how
levels of reticulation in evolution vary across domains.


  In this work, we review the figures used to characterize an epidemic outbreak
most. Particular attention is drawn to epidemic spreading at time-varying
transition rates. A time-varying SIR-like model is used to describe the
epidemic dynamics, from which closed equations relating parameters and key
quantities like reproduction ratio and doubling time are derived. The
definition and computation of these metrics are revisited in the context of the
general solution to the time-varying model dynamics, focusing on the
similarities and differences with the time-invariant case. Further, the
prediction of these metrics, that is of the disease evolution, as response to
different scenarios is also investigated.


  Processing time-dependent information requires cells to quantify the duration
of past regulatory events and program the time span of future signals. At the
single-cell level, timer mechanisms can be implemented with genetic circuits:
sets of genes connected to achieve specific tasks. However, such systems are
difficult to implement in single cells due to saturation in molecular
components and stochasticity in the limited intracellular space. Multicellular
implementations, on the other hand, outsource some of the components of
information-processing circuits to the extracellular space, and thereby might
escape those constraints. Here we develop a theoretical framework, based on a
trilinear coordinate representation, to study the collective behavior of a
cellular population composed of three cell types under stationary conditions.
This framework reveals that distributing different processes (in our case the
production, detection and degradation of a time-encoding signal) across
distinct cell types enables the robust implementation of a multicellular timer.
Our analysis also shows the circuit to be easily tunable by varying the
cellular composition of the consortium.


  Motivation: Estimating model parameters from experimental observations is one
of the key challenges in systems biology and can be computationally very
expensive. While the Julia programming language was recently developed as a
high-level and high-performance language for scientific computing, systems
biologists have only started to realise its potential. For instance, we have
recently used Julia to cut down the optimization time of a microbial community
model by a factor of 140. To facilitate access of the systems biology community
to the efficient nonlinear solvers used for this optimisation, we developed
SBML2Julia. SBML2Julia translates optimisation problems specified in SBML and
TSV files (PEtab format) into Julia for Mathematical Programming (JuMP),
executes the optimization and returns the results in tabular format.
  Availability and implementation: SBML2Julia is freely available under the MIT
license. It comes with a command line interface and Python API. Internally,
SBML2Julia calls the Julia LTS release v1.0.5 for optimisation. All necessary
dependencies can be pulled from Docker Hub
(https://hub.docker.com/repository/docker/paulflang/sbml2julia). Source code
and documentation are available at https://github.com/paulflang/SBML2Julia.


  In this article, we analyze the spread of cases resulting from community
transmission of COVID-19 in Senegal in order to identify statistical
associations. The identification and knowledge of the factors associated with
this community transmission can be a decision support tool to limit the spread
of the disease. We estimate parameters and evaluate the growth factor,
community rate, weekly increase and daily difference, and make forecasting to
help on how to find concrete actions to control the situation.


  Fast and effective drug discovery processes rely on rational drug design to
circumvent the tedious and expensive trial and error approach. However,
accurate predictions of new remedies, which are often enzyme inhibitors,
require a clear understanding of the nature and function of the key players
governing the interaction between the drug candidate and its target. Here, we
propose to calculate electric fields to explicitly link structure to function
in molecular dynamics simulations, a method that can easily be integrated
within the rational drug discovery workflow. By projecting the electric fields
onto specific bonds, we can identify the system components that are at the
origin of stabilizing intermolecular interactions (covalent and non-covalent)
in the active site. This helps to significantly narrow the exploration space
when predicting new inhibitors. To illustrate this method, we characterize the
binding of the non-covalent inhibitor X77 to the main protease of SARS-CoV-2, a
particularly time-sensitive drug discovery problem. With electric field
calculations, we were able to identify 3 key residues (Asn-142, Met-165 and
Glu-166), that have functional consequences on X77. This contrasts with the
nearly 20 residues reported in previous studies as being in close contact with
inhibitors in the active site of the protease. As a result, the search for new
non-covalent inhibitors can now be accelerated by techniques that look to
optimize the interaction between candidate molecules and these residues.


  Serial electron diffraction (SerialED) is an emerging technique, which
applies the snapshot data-collection mode of serial X-ray crystallography to
three-dimensional electron diffraction (3D ED), forgoing the conventional
rotation method. Similarly to serial X-ray crystallography, this approach leads
to almost complete absence of radiation damage effects even for the most
sensitive samples, and allows for a high level of automation. However, SerialED
also necessitates new techniques of data processing, which combine existing
pipelines for rotation electron diffraction and serial X-ray crystallography
with some more particular solutions for challenges arising in SerialED
specifically. Here, we introduce our analysis pipeline for SerialED data, and
its implementation using the CrystFEL and diffractem program packages. Detailed
examples are provided in extensive supplementary code.


  The study of biological processes can greatly benefit from tools that
automatically predict gene functions or directly cluster genes based on shared
functionality. Existing data mining methods predict protein functionality by
exploiting data obtained from high-throughput experiments or meta-scale
information from public databases. Most existing prediction tools are targeted
at predicting protein functions that are described in the gene ontology (GO).
However, in many cases biologists wish to discover functionally related genes
for which GO terms are inadequate. In this paper, we introduce a rich set of
features and use them in conjunction with semisupervised learning approaches in
order to expand an initial set of seed genes to a larger cluster of
functionally related genes. Among all the semi-supervised methods that were
evaluated, the framework of learning with positive and unlabeled examples (LPU)
is shown to be especially appropriate for mining functionally related genes.
When evaluated on experimentally validated benchmark data, the LPU approaches1
significantly outperform a standard supervised learning algorithm as well as an
established state-of-the-art method. Given an initial set of seed genes, our
best performing approach could be used to mine functionally related genes in a
wide range of organisms.


  Subclasses of lymphocytes carry different functional roles to work together
to produce an immune response and lasting immunity. Additionally to these
functional roles, T and B-cell lymphocytes rely on the diversity of their
receptor chains to recognize different pathogens. The lymphocyte subclasses
emerge from common ancestors generated with the same diversity of receptors
during selection processes. Here we leverage biophysical models of receptor
generation with machine learning models of selection to identify specific
sequence features characteristic of functional lymphocyte repertoires and
subrepertoires. Specifically using only repertoire level sequence information,
we classify CD4$^+$ and CD8$^+$ T-cells, find correlations between receptor
chains arising during selection and identify T-cells subsets that are targets
of pathogenic epitopes. We also show examples of when simple linear classifiers
do as well as more complex machine learning methods.


  Mathematical modelling of the spread of epidemics has been an interesting
challenge in the field of epidemiology. The SIR Model proposed by Kermack and
McKendrick in 1927 is a prototypical model of epidemiology. However, it has its
limitations. In this paper, we show two independent ways of generalizing this
model, the first one if the vaccine isn't discovered or ready to use and the
next one, if the vaccine is discovered and ready to use. In the first part, we
have pointed out a major over-simplification, i.e., assumption of variation of
the time derivatives of the variables with the linear or quadratic powers of
the individual variables and introduce two new parameters to incorporate
further nonlinearity in the number of infected people in the model. As a result
of this, we show how this additional nonlinearity, in the newly introduced
parameters, can bring a significant shift in the peak time of infection, i.e.,
the time at which the infected population reaches maximum. We show that in
special cases, even we can get a transition from epidemic to a non-epidemic
stage of a particular infectious disease. We further study one such special
case and treat it as a problem of phase transition. Then, we investigate all
the necessary parameters of this phase transition, like the order parameter and
critical exponent. We observe that $O_p \sim (q_c-q)^{\beta}$. {\it As far as
we know the phase transition and its quantification in terms of the scaling
behaviour is not yet know in the context of pandemic}. In the second part, we
incorporate in the model, a consideration of artificial herd immunity and show
how we can decrease the peak time of infection with a subsequent decrease in
the maximum number of infected people.


  Here, we report analysis and summary of research in the field of localization
microscopy for optical imaging. We introduce the basic elements of
super-resolved localization microscopy methods for PALM and STORM, commonly
used both in vivo and in vitro, discussing the core essentials of background
theory, instrumentation and computational algorithms. We discuss the resolution
limit of light microscopy and the mathematical framework for localizing
fluorescent dyes in space beyond this limit, including the precision obtainable
as a function of the amount of light emitted from a dye, and how it leads to a
fundamental compromise between spatial and temporal precision. The properties
of a "good dye" are outlined, as are the features of PALM and STORM
super-resolution microscopy and adaptations that may need to be made to
experimental protocols to perform localization determination. We analyse
briefly some of the methods of modern super-resolved optical imaging that work
through reshaping point spread functions and how they utilize aspects of
localization microscopy, such as stimulated depletion (STED) methods and
MINFLUX, and summarize modern methods that push localization into 3D using
non-Gaussian point spread functions. We report on current methods for analyzing
localization data including determination of 2D and 3D diffusion constants,
molecular stoichiometries, and performing cluster analysis with cutting-edge
techniques, and finally discuss how these techniques may be used to enable
important insight into a range of biological processes.


  A central feature of an emerging infectious disease in a pandemic scenario is
the spread through geographical scales and the impacts on different locations
according to the adopted mitigation protocols. We investigated a stochastic
epidemic model with the metapopulation approach in which patches represent
municipalities. Contagion follows a stochastic compartmental model for
municipalities; the latter, in turn, interact with each other through recurrent
mobility. As a case of study, we consider the epidemic of COVID-19 in Brazil
performing data-driven simulations. Properties of the simulated epidemic curves
have very broad distributions across different geographical locations and
scales, from states, passing through intermediate and immediate regions down to
municipality levels. Correlations between delay of the epidemic outbreak and
distance from the respective capital cities were predicted to be strong in
several states and weak in others, signaling influences of multiple epidemic
foci propagating towards the inland cities. Responses of different regions to a
same mitigation protocol can vary enormously implying that the policies of
combating the epidemics must be engineered according to the region' specificity
but integrated with the overall situation. Real series of reported cases
confirm the qualitative scenarios predicted in simulations. Even though we
restricted our study to Brazil, the prospects and model can be extended to
other geographical organizations with heterogeneous demographic distributions.


  Human Papillomaviruses (HPVs) are involved in the etiology of anogenital and
head and neck cancers. The HPV DNA prevalence greatly differs by anatomical
site. Indeed, the high rates of viral DNA prevalence in anal and cer-vical
carcinomas contrast with the lower fraction of cancer cases attributable to
HPVs in other anatomical sites, chiefly the vulva, the penis and head and neck.
Here we analyzed 2635 Formalin Fixed Paraffin Embedded surgical samples that
had previously tested negative for the presence of HPVs DNA using the
SPF10/DEIA procedure, in order to identify the presence of other PVs not
explicitly targeted by standard molecular epidemiologic approaches. All samples
were reanalyzed using five broad-PV PCR primer sets (CP1/2, FAP6064/FAP64,
SKF/SKR, MY9/MY11, MFI/MFII) targeting the main PV main clades. In head and
neck carcinoma samples (n = 1141), we recovered DNA from two BetaHPVs, namely
HPV20 and HPV21, and from three cutaneous AlphaPVs, namely HPV2, HPV57 and
HPV61. In vulvar squamous cell carcinoma samples (n = 902), we found one of the
samples containing DNA of one cutaneous HPV, namely HPV2, and 29 samples
contained DNA from essentially mucosal HPVs. In penile squamous cell carcinoma
samples (n = 592), we retrieved the DNA of HPV16 in 16 samples. Our results
show first that the SPF10/DEIA is very sensitive, as we recovered only 2.1%
(55/2635) false negative results; second, that although the DNA of cutaneous
HPVs can be detected in cancer samples, their relative contribution remains
anyway minor (0.23%; 6/2635) and may be neglected for screening and vaccination
purposes; and third, their contribution to malignancy is not necessarily
warranted and needs to be elucidated.


  The protein-protein interactions (PPIs) of 14-3-3 proteins are a model system
for studying PPI stabilization. The complex natural product Fusicoccin A
stabilizes many 14-3-3 PPIs but is not amenable for use in SAR studies,
motivating the search for more drug-like chemical matter. However, drug-like
14-3-3 PPI stabilizers enabling such study have remained elusive. An X-ray
crystal structure of a PPI in complex with an extremely low potency stabilizer
uncovered an unexpected non-protein interacting, ligand-chelated Mg 2+ leading
to the discovery of metal ion-dependent 14-3-3 PPI stabilization potency. This
originates from a novel chelation-controlled bioactive conformation
stabilization effect. Metal chelation has been associated with pan-assay
interference compounds (PAINS) and frequent hitter behavior, but chelation can
evidently also lead to true potency gains and find use as a medicinal chemistry
strategy to guide compound optimization. To demonstrate this, we exploited the
effect to design the first potent, selective and drug-like 14-3-3 PPI
stabilizers.


  Despite considerable progress, ab initio protein structure prediction remains
suboptimal. A crowdsourcing approach is the online puzzle video game Foldit,
that provided several useful results that matched or even outperformed
algorithmically computed solutions. Using Foldit, the WeFold crowd had several
successful participations in the Critical Assessment of Techniques for Protein
Structure Prediction. Based on the recent Foldit standalone version, we trained
a deep reinforcement neural network called DeepFoldit to improve the score
assigned to an unfolded protein, using the Q-learning method with experience
replay. This paper is focused on model improvement through hyperparameter
tuning. We examined various implementations by examining different model
architectures and changing hyperparameter values to improve the accuracy of the
model. The new model hyper-parameters also improved its ability to generalize.
Initial results, from the latest implementation, show that given a set of small
unfolded training proteins, DeepFoldit learns action sequences that improve the
score both on the training set and on novel test proteins. Our approach
combines the intuitive user interface of Foldit with the efficiency of deep
reinforcement learning.


  Recently, there has been great interest in learning how to best represent
proteins, specifically with fixed-length embeddings. Deep learning has become a
popular tool for protein representation learning as a model's hidden layers
produce potentially useful vector embeddings. TAPE introduced a number of
benchmark tasks and showed that semi-supervised learning, via pretraining
language models on a large protein corpus, improved performance on downstream
tasks. Two of the tasks (fluorescence prediction and stability prediction)
involve learning fitness landscapes. In this paper, we show that CNN models
trained solely using supervised learning both compete with and sometimes
outperform the best models from TAPE that leverage expensive pretraining on
large protein datasets. These CNN models are sufficiently simple and small that
they can be trained using a Google Colab notebook. We also find for the
fluorescence task that linear regression outperforms our models and the TAPE
models. The benchmarking tasks proposed by TAPE are excellent measures of a
model's ability to predict protein function and should be used going forward.
However, we believe it is important to add baselines from simple models to put
the performance of the semi-supervised models that have been reported so far
into perspective.


  This paper proposes a workflow to identify genes that respond to specific
treatments in plants. The workflow takes as input the RNA sequencing read
counts and phenotypical data of different genotypes, measured under control and
treatment conditions. It outputs a reduced group of genes marked as relevant
for treatment response. Technically, the proposed approach is both a
generalization and an extension of WGCNA. It aims to identify specific modules
of overlapping communities underlying the co-expression network of genes.
Module detection is achieved by using Hierarchical Link Clustering. The
overlapping nature of the systems' regulatory domains that generate
co-expression can be identified by such modules. LASSO regression is employed
to analyze phenotypic responses of modules to treatment.
  Results. The workflow is applied to rice (Oryza sativa), a major food source
known to be highly sensitive to salt stress. The workflow identifies 19 rice
genes that seem relevant in the response to salt stress. They are distributed
across 6 modules: 3 modules, each grouping together 3 genes, are associated to
shoot K content; 2 modules of 3 genes are associated to shoot biomass; and 1
module of 4 genes is associated to root biomass. These genes represent target
genes for the improvement of salinity tolerance in rice.
  Conclusion. A more effective framework to reduce the search-space for target
genes that respond to a specific treatment is introduced. It facilitates
experimental validation by restraining efforts to a smaller subset of genes of
high potential relevance.


  Based on the BioBricks standard, restriction synthesis is a novel catabolic
iterative DNA synthesis method that utilizes endonucleases to synthesize a
query sequence from a reference sequence. In this work, the reference sequence
is built from shorter subsequences by classifying them as applicable or
inapplicable for the synthesis method using three different machine learning
methods: Support Vector Machines (SVMs), random forest, and Convolution Neural
Networks (CNNs). Before applying these methods to the data, a series of feature
selection, curation, and reduction steps are applied to create an accurate and
representative feature space. Following these preprocessing steps, three
different pipelines are proposed to classify subsequences based on their
nucleotide sequence and other relevant features corresponding to the
restriction sites of over 200 endonucleases. The sensitivity using SVMs, random
forest, and CNNs are 94.9%, 92.7%, 91.4%, respectively. Moreover, each method
scores lower in specificity with SVMs, random forest, and CNNs resulting in
77.4%, 85.7%, and 82.4%, respectively. In addition to analyzing these results,
the misclassifications in SVMs and CNNs are investigated. Across these two
models, different features with a derived nucleotide specificity visually
contribute more to classification compared to other features. This observation
is an important factor when considering new nucleotide sensitivity features for
future studies.


  The amyloid state of proteins is widely studied with relevancy in neurology,
biochemistry, and biotechnology. In contrast with amorphous aggregation, the
amyloid state has a well-defined structure, consisting of parallel and
anti-parallel $\beta$-sheets in a periodically repeated formation. The
understanding of the amyloid state is growing with the development of novel
molecular imaging tools, like cryogenic electron microscopy. Sequence-based
amyloid predictors were developed by using mostly artificial neural networks
(ANNs) as the underlying computational techniques. From a good neural
network-based predictor, it is a very difficult task to identify those
attributes of the input amino acid sequence, which implied the decision of the
network. Here we present a Support Vector Machine (SVM)-based predictor for
hexapeptides with correctness higher than 84\%, i.e., it is at least as good as
the published ANN-based tools. Unlike the artificial neural networks, the
decision of the SVMs are much easier to analyze, and from a good predictor, we
can infer rich biochemical knowledge.
  Availability and Implementation: The Budapest Amyloid Predictor webserver is
freely available at https://pitgroup.org/bap.


  In the current scenario of climate change, cold-adapted insects are among the
most threatened organisms in high-altitude habitats of the Alps. Upslope shifts
and changes in phenology are two of the most investigated responses to climate
change, but there is an increasing interest in evaluating the presence of
high-altitude landforms acting as refugia. Nebria germari Heer, 1837
(Coleoptera: Carabidae) is a hygrophilic and cold-adapted species that still
exhibits large populations on supraglacial debris of the Eastern Alps. This
work aims at describing the ecology and phenology of the populations living on
supraglacial debris. To this end, we analysed the populations from three
Dolomitic glaciers whose surfaces are partially covered by stony debris. We
found that supraglacial debris is characterised by more stable colder and
wetter conditions than the surrounding debris slopes and by a shorter snow-free
period. The populations found on supraglacial debris were spring breeders,
differently from those documented in the 1980s on Dolomitic high alpine
grasslands, which were reported as autumn breeders. Currently Nebria germari
seems therefore to find a suitable habitat on supraglacial debris, where
micrometeorological conditions are appropriate for its life-cycle and
competition and predation are reduced.


  The question of why we age is a fundamental one. It is about who we are, and
it also might have critical practical aspects as we try to find ways to age
slower. Or to not age at all. Different reasons point at distinct strategies
for the research of anti-ageing drugs. While the main reason why biological
systems work as they do is evolution, for quite a while, it was believed that
aging required another explanation. Aging seems to harm individuals so much
that even if it has group benefits, those benefits were unlikely to be enough.
That has led many scientists to propose non-evolutionary explanations as to why
we age. But those theories seem to fail at explaining all the data on how
species age. Here, I will show that the insistence of finding the one idea that
explains it all might be at the root of the difficulty of getting a full
picture. By exploring an evolutionary model of aging where locality and
temporal changes are fundamental aspects of the problem, I will show that
environmental change causes the barrier for group advantages to become much
weaker. That weakening might help small group advantages to add up to the point
they could make an adaptive difference. To answer why we age, we might have to
abandon asking which models are correct. The full answer might come from
considering how much each hypothesis behind each existing model, evolutionary
and non-evolutionary ones, contributes to the real world's solution.


  We developed a new base caller DeepNano-coral for nanopore sequencing, which
is optimized to run on the Coral Edge Tensor Processing Unit, a small
USB-attached hardware accelerator. To achieve this goal, we have designed new
versions of two key components used in convolutional neural networks for speech
recognition and base calling. In our components, we propose a new way of
factorization of a full convolution into smaller operations, which decreases
memory access operations, memory access being a bottleneck on this device.
DeepNano-coral achieves real-time base calling during sequencing with the
accuracy slightly better than the fast mode of the Guppy base caller and is
extremely energy efficient, using only 10W of power. Availability:
https://github.com/fmfi-compbio/coral-basecaller


  Transcription factors are proteins that regulate gene activity by activating
or repressing gene transcription. A special class of transcriptional repressors
operates via a short-range mechanism, making local DNA regions inaccessible to
binding by activators, and thus providing an indirect repressive action on the
target gene. This mechanism is commonly modeled assuming that repressors
interact with DNA under thermodynamic equilibrium and neglecting some
configurations of the gene regulatory region. We elaborate on a more general
nonequilibrium model of short-range repression using the graph formalism for
transitions between gene states, and we apply analytical calculations to
compare it with the equilibrium model in terms of the repression strength and
expression noise. In contrast to the equilibrium approach, the new model allows
us to separate two basic mechanisms of short-range repression. The first
mechanism is associated with the recruiting of factors that mediate chromatin
condensation, and the second one concerns the blocking of factors that mediate
chromatin loosening. The nonequilibrium model demonstrates better performance
on previously published gene expression data obtained for transcription factors
controlling Drosophila development, and furthermore it predicts that the first
repression mechanism is the most favorable in this system. The presented
approach can be scaled to larger gene networks and can be used to infer
specific modes and parameters of transcriptional regulation from gene
expression data.


  A biological system achieve homeostasis when there is a regulated quantity
that is maintained within a narrow range of values. Here we consider
homeostasis as a phenomenon of network dynamics. In this context, we improve a
general theory for the analysis of homeostasis in network dynamical systems
with distinguished input and output nodes, called `input-output networks'. The
theory allows one to define `homeostasis types' of a given network in a `model
independent' fashion, in the sense that the classification depends on the
network topology rather than on the specific model equations. Each `homeostasis
type' represents a possible mechanism for generating homeostasis and is
associated with a suitable `subnetwork motif' of the original network. Our
contribution is an extension of the theory to the case of networks with
multiple input nodes. To showcase our theory, we apply it to bacterial
chemotaxis, a paradigm for homeostasis in biochemical systems. By considering a
representative model of Escherichia coli chemotaxis, we verify that the
corresponding abstract network has multiple input nodes. Thus showing that our
extension of the theory allows for the inclusion of an important class of
models that were previously out of reach. Moreover, from our abstract point of
view, the occurrence of homeostasis in the studied model is caused by a new
mechanism, called input counterweight homeostasis. This new homeostasis
mechanism was discovered in the course of our investigation and is generated by
a balancing between the several input nodes of the network -- therefore, it
requires the existence of at least two input nodes to occur. Finally, the
framework developed here allows one to formalize a notion of `robustness' of
homeostasis based on the concept of `genericity' from the theory dynamical
systems. We discuss how this kind of robustness of homeostasis appears in the
chemotaxis model.


  Turing instabilities of reaction-diffusion systems can only arise if the
diffusivities of the chemical species are sufficiently different. This
threshold is unphysical in most systems with $N=2$ diffusing species, forcing
experimental realizations of the instability to rely on fluctuations or
additional nondiffusing species. Here we ask whether this diffusive threshold
lowers for $N>2$ to allow "true" Turing instabilities. Inspired by May's
analysis of the stability of random ecological communities, we analyze the
probability distribution of the diffusive threshold in reaction-diffusion
systems defined by random matrices describing linearized dynamics near a
homogeneous fixed point. In the numerically tractable cases $N\leqslant 6$, we
find that the diffusive threshold becomes more likely to be smaller and
physical as $N$ increases and that most of these many-species instabilities
cannot be described by reduced models with fewer species.


  Drug combinations play an important role in therapeutics due to its better
efficacy and reduced toxicity. Recent approaches have applied machine learning
to identify synergistic combinations for cancer, but they are not applicable to
new diseases with limited combination data. Given that drug synergy is closely
tied to biological targets, we propose a \emph{biological bottleneck} model
that jointly learns drug-target interaction and synergy. The model consists of
two parts: a drug-target interaction and target-disease association module.
This design enables the model to \emph{explain} how a biological target affects
drug synergy. By utilizing additional biological information, our model
achieves 0.78 test AUC in drug synergy prediction using only 90 COVID drug
combinations for training. We experimentally tested the model predictions in
the U.S. National Center for Advancing Translational Sciences (NCATS)
facilities and discovered two novel drug combinations (Remdesivir + Reserpine
and Remdesivir + IQ-1S) with strong synergy in vitro.


  In this paper, we firstly propose SQIARD and SIARD models to investigate the
transmission of COVID-19 with quarantine, infected and asymptomatic infected,
and discuss the relation between the respective basic reproduction number $R_0,
R_Q$ and the stability of the equilibrium points of model. Secondly, after
training the related data parameters, in our numerical simulations, we
respectively conduct the forecast of the data of US, South Korea, Brazil,
India, Russia and Italy, and the effect of prediction of the epidemic situation
in each country. Furthermore, we apply US data to compare SQIARD with SIARD,
and display the effects of predictions.


  Chronic kidney disease (CKD) is a gradual loss of renal function over time,
and it increases the risk of mortality, decreased quality of life, as well as
serious complications. The prevalence of CKD has been increasing in the last
couple of decades, which is partly due to the increased prevalence of diabetes
and hypertension. To accurately detect CKD in diabetic patients, we propose a
novel framework to learn sparse longitudinal representations of patients'
medical records. The proposed method is also compared with widely used
baselines such as Aggregated Frequency Vector and Bag-of-Pattern in Sequences
on real EHR data, and the experimental results indicate that the proposed model
achieves higher predictive performance. Additionally, the learned
representations are interpreted and visualized to bring clinical insights.


  Pairwise alignment of DNA sequencing data is a ubiquitous task in
bioinformatics and typically represents a heavy computational burden.
State-of-the-art approaches to speed up this task use hashing to identify short
segments (k-mers) that are shared by pairs of reads, which can then be used to
estimate alignment scores. However, when the number of reads is large,
accurately estimating alignment scores for all pairs is still very costly.
Moreover, in practice, one is only interested in identifying pairs of reads
with large alignment scores. In this work, we propose a new approach to
pairwise alignment estimation based on two key new ingredients. The first
ingredient is to cast the problem of pairwise alignment estimation under a
general framework of rank-one crowdsourcing models, where the workers'
responses correspond to k-mer hash collisions. These models can be accurately
solved via a spectral decomposition of the response matrix. The second
ingredient is to utilise a multi-armed bandit algorithm to adaptively refine
this spectral estimator only for read pairs that are likely to have large
alignments. The resulting algorithm iteratively performs a spectral
decomposition of the response matrix for adaptively chosen subsets of the read
pairs.


  A timely understanding of the spatiotemporal pattern and development trend of
COVID-19 is critical for timely prevention and control. However, the
under-reporting of cases is widespread in fields associated with public health.
It is also possible to draw biased inferences and formulate inappropriate
prevention and control policies if the phenomenon of under-reporting is not
taken into account. Therefore, in this paper, a novel framework was proposed to
explore the impact of under-reporting on COVID-19 spatiotemporal distributions,
and empirical analysis was carried out using infection data of healthcare
workers in Wuhan and Hubei (excluding Wuhan). The results show that (1) the
lognormal distribution was the most suitable to describe the evolution of
epidemic with time; (2) the estimated peak infection time of the reported cases
lagged the peak infection time of the healthcare worker cases, and the
estimated infection time interval of the reported cases was smaller than that
of the healthcare worker cases. (3) The impact of under-reporting cases on the
early stages of the pandemic was greater than that on its later stages, and the
impact on the early onset area was greater than that on the late onset area.
(4) Although the number of reported cases was lower than the actual number of
cases, a high spatial correlation existed between the cumulatively reported
cases and healthcare worker cases. The proposed framework of this study is
highly extensible, and relevant researchers can use data sources from other
counties to carry out similar research.


  De novo genome assembly focuses on finding connections between a vast amount
of short sequences in order to reconstruct the original genome. The central
problem of genome assembly could be described as finding a Hamiltonian path
through a large directed graph with a constraint that an unknown number of
nodes and edges should be avoided. However, due to local structures in the
graph and biological features, the problem can be reduced to graph
simplification, which includes removal of redundant information. Motivated by
recent advancements in graph representation learning and neural execution of
algorithms, in this work we train the MPNN model with max-aggregator to execute
several algorithms for graph simplification. We show that the algorithms were
learned successfully and can be scaled to graphs of sizes up to 20 times larger
than the ones used in training. We also test on graphs obtained from real-world
genomic data---that of a lambda phage and E. coli.


  Understanding the relationship between complexity and stability in large
dynamical systems -- such as ecosystems -- remains a key open question in
complexity theory which has inspired a rich body of work developed over more
than fifty years. The vast majority of this theory addresses asymptotic linear
stability around equilibrium points, but the idea of `stability' in fact has
other uses in the empirical ecological literature. The important notion of
`temporal stability' describes the character of fluctuations in population
dynamics, driven by intrinsic or extrinsic noise. Here we apply tools from
random matrix theory to the problem of temporal stability, deriving analytical
predictions for the fluctuation spectra of complex ecological networks. We show
that different network structures leave distinct signatures in the spectrum of
fluctuations, and demonstrate the application of our theory to the analysis
ecological timeseries data of plankton abundances.


  Systemic lupus erythematosus (SLE) is the tenth leading cause of death in
females 15-24 years old in the US. The diversity of symptoms and immune
pathways expressed in SLE patients causes difficulties in treating SLE as well
as in new clinical trials. This study used unsupervised learning on gene
expression data from adult SLE patients to separate patients into clusters. The
dimensionality of the gene expression data was reduced by three separate
methods (PCA, UMAP, and a simple linear autoencoder) and the results from each
of these methods were used to separate patients into six clusters with k-means
clustering.
  The clusters revealed three separate immune pathways in the SLE patients that
caused SLE. These pathways were: (1) high interferon levels, (2) high
autoantibody levels, and (3) dysregulation of the mitochondrial apoptosis
pathway. The first two pathways have been extensively studied in SLE. However,
mitochondrial apoptosis has not been investigated before to the best of our
knowledge as a standalone cause of SLE, independent of autoantibody production,
indicating that mitochondrial proteins could lead to a new set of therapeutic
targets for SLE in future research.


  Antagonistic interactions in biological systems, which occur when one
perturbation blunts the effect of another, are typically interpreted as
evidence that the two perturbations impact the same cellular pathway or
function. Yet, this interpretation ignores extreme antagonistic interactions
wherein an otherwise deleterious perturbation compensates for the function lost
due to a prior perturbation. Here, we report on gene-environment interactions
involving genetic mutations that are deleterious in a permissive environment
but beneficial in a specific environment that restricts growth. These extreme
antagonistic interactions constitute gene-environment analogs of synthetic
rescues previously observed for gene-gene interactions. Our approach uses two
independent adaptive evolution steps to address the lack of experimental
methods to systematically identify such extreme interactions. We apply the
approach to Escherichia coli by successively adapting it to defined glucose
media without and with the antibiotic rifampicin. The approach identified
multiple mutations that are beneficial in the presence of rifampicin and
deleterious in its absence. The analysis of transcription shows that the
antagonistic adaptive mutations repress a stringent response-like
transcriptional program, while non-antagonistic mutations have an opposite
transcriptional profile. Our approach represents a step toward the systematic
characterization of extreme antagonistic gene-drug interactions, which can be
used to identify targets to select against antibiotic resistance.


  By fitting a compartment ODE model for Covid-19 propagation to cumulative
case and death data for US states and European countries, we find that the case
mortality rate seems to have decreased by at least 80% in most of the US and at
least 90% in most of Europe. These are much larger and faster changes than
reported in empirical studies, such as the 18% decrease in mortality found for
the New York City hospital system from March to August 2020 (Horwitz et al,
Trends in Covid-19 risk-adjusted mortality rates, J. Hosp. Med. 2020). Our
reported decreases surprisingly do not have strong correlations to other model
parameters (such as contact rate) or other standard state/national metrics such
as population density, GDP, and median age. Almost all the decreases occurred
between mid-April and mid-June, which unexpectedly corresponds to the time when
many state and national lockdowns were released resulting in surges of new
cases. Several plausible causes for this drop are examined, such as
improvements in treatment, face mask wearing, a new virus strain, and
potentially changing demographics of infected patients, but none are
overwhelmingly convincing given the currently available evidence.


  The native state structures of globular proteins are stable and well-packed
indicating that self-interactions are favored over protein-solvent interactions
under folding conditions. We use this as a guiding principle to derive the
geometry of the building blocks of protein structures, alpha-helices and
strands assembled into beta-sheets, with no adjustable parameters, no amino
acid sequence information, and no chemistry. There is an almost perfect fit
between the dictates of mathematics and physics and the rules of quantum
chemistry. Our theory establishes an energy landscape that channels protein
evolution by providing sequence-independent platforms for elaborating
sequence-dependent functional diversity. Our work highlights the vital role of
discreteness in life and has implications for the creation of artificial life
and on the nature of life elsewhere in the cosmos.


  Recent studies have revealed effects of various tree species on soil physical
and chemical properties. However, effects of various tree species on
composition and activity of soil microbiota and the relevant controls remain
poorly understood. We evaluated the influence of tree species associated with
two different mycorrhizal types, ectomycorrhiza (EcM) and arbuscular mycorrhiza
(AM), on growth, biomass and metabolic activity of soil fungal and bacterial
communities using common garden tree species experiments throughout Denmark.
The soil microbial communities differed between six European tree species as
well as between EcM (beech, lime, oak and spruce) and AM (ash and maple) tree
species. The EcM tree species had higher fungal biomass, fungal growth and
bacterial biomass, while AM species showed higher bacterial growth. The results
indicated that microbial community composition and functioning differed between
groups of tree species with distinct litter qualities that generate soil C/N
ratio and soil pH differences. The mycorrhizal association only partly
explained litter quality and soil microbial species differences since lime was
more similar to AM tree species. In addition, our results indicated that tree
species-mediated soil pH and C/N ratio were the most important variables
shaping microbial communities with a positive effect on bacterial and a
negative effect on fungal growth rates. The results suggest that tree
species-mediated microbial community composition and activity may be important
drivers of the different vertical soil C distribution previously observed in AM
and EcM tree species.


  In this paper, a mathematical analysis of the global dynamics of a viral
infection model in vivo is carried out. We study the dynamics of a hepatitis C
virus (HCV) model, under therapy, that considers both extracellular and
intracellular levels of infection. At present most mathematical modeling of
viral kinetics after treatment only addresses the process of infection of a
cell by the virus and the release of virions by the cell, while the processes
taking place inside the cell are not included. We prove that the solutions of
the new model with positive initial values are positive, exist globally in time
and are bounded. The model has two virus-free steady states. They are
distinguished by the fact that viral RNA is absent inside the cells in the
first state and present inside the cells in the second. There are basic
reproduction numbers associated to each of these steady states. If the basic
reproduction number of the first steady state is less than one then that state
is asymptotically stable. If the basic reproduction number of the first steady
state is greater than one and that of the second less than one then the second
steady state is asymptotically stable. If both basic reproduction numbers are
greater than one then we obtain various conclusions which depend on different
restrictions on the parameters of the model. Under increasingly strong
assumptions we prove that there is at least one positive steady state (infected
equilibrium), that there is a unique positive steady state and that the
positive steady state is stable. We also give a condition under which every
positive solution converges to a positive steady state. This is proved by
methods of Li and Muldowney. Finally, we illustrate the theoretical results by
numerical simulations.


  We introduce a Python package for modeling and studying the spread of
infectious diseases using an agent-based SEIR style epidemiological model with
a focus on university campuses. This document explains the epidemiological
model used in the package and gives examples highlighting the ways that the
package can be used.


  Elongation is a fundament process in amyloid fiber growth, which is normally
characterized by a linear relationship between the fiber elongation rate and
the monomer concentration. However, in high concentration regions, a sub-linear
dependence was often observed, which could be explained by a universal
saturation mechanism. In this paper, we modeled the saturated elongation
process through a Michaelis-Menten like mechanism, which is constituted by two
sub-steps -- unspecific association and dissociation of a monomer with the
fibril end, and subsequent conformational change of the associated monomer to
fit itself to the fibrillar structure. Typical saturation concentrations were
found to be $7-70\mu M$ for A$\beta$40, $\alpha$-synuclein and etc.
Furthermore, by using a novel Hamiltonian formulation, analytical solutions
valid for both weak and strong saturated conditions were constructed and
applied to the fibrillation kinetics of $\alpha$-synuclein and silk fibroin.


  We consider a compartmental model from which we incorporate a time-dependent
health care capacity having a logistic growth. This allows us to take into
account the Senegalese authorities response in anticipating the growing number
of infected cases. We highlight the importance of anticipation and timing to
avoid overwhelming that could impact considerably the treatment of patients and
the well-being of health care workers. A condition, depending on the health
care capacity and the flux of new hospitalized individuals, to avoid possible
overwhelming is provided. We also use machine learning approach to project
forward the cumulative number of cases from March 02, 2020, until 1st December,
2020.


  For a stochastic COVID-19 model with jump-diffusion, we prove the existence
and uniqueness of the global positive solution. We also investigate some
conditions for the extinction and persistence of the disease. We calculate the
threshold of the stochastic epidemic system which determines the extinction or
permanence of the disease at different intensities of the stochastic noises.
This threshold is denoted by $\xi$ which depends on the white and jump noises.
The effects of these noises on the dynamics of the model are studied. The
numerical experiments show that the random perturbation introduced in the
stochastic model suppresses disease outbreaks as compared to its deterministic
counterpart. In other words, the impact of the noises on the extinction and
persistence is high. When the noise is large or small, our numerical findings
show that the COVID-19 vanishes from the population if $\xi <1;$ whereas the
epidemic can't go out of control if $\xi >1.$ From this, we observe that white
noise and jump noise have a significant effect on the spread of COVID-19
infection, i.e., we can conclude that the stochastic model is more realistic
than the deterministic one. Finally, to illustrate this phenomenon, we put some
numerical simulations.


  Since the recent introduction of several viable vaccines for SARS-CoV-2,
vaccination uptake has become the key factor that will determine our success in
containing the COVID-19 pandemic. We argue that game theory and social network
models should be used to guide decisions pertaining to vaccination programmes
for the best possible results. In the months following the introduction of
vaccines, their availability and the human resources needed to run the
vaccination programmes have been scarce in many countries. Vaccine hesitancy is
also being encountered from some sections of the general public. We emphasize
that decision-making under uncertainty and imperfect information, and with only
conditionally optimal outcomes, is a unique forte of established game-theoretic
modelling. Therefore, we can use this approach to obtain the best framework for
modelling and simulating vaccination prioritization and uptake that will be
readily available to inform important policy decisions for the optimal control
of the COVID-19 pandemic.


  Cell type-specific gene expression patterns are represented as memory states
of a Hopfield neural network model. It is shown that order parameters of this
model can be interpreted as concentrations of master transcription regulators
that form concurrent positive feedback loops with a large number of downstream
regulated genes. The order parameter free energy then defines an epigenetic
landscape in which local minima correspond to stable cell states. The model is
applied to gene expression data in the context of hematopoiesis.


  In this work, we adapt the epidemiological SIR model to study the evolution
of the dissemination of COVID-19 in Germany and Brazil (nationally, in the
State of Paraiba, and in the City of Campina Grande). We prove the well
posedness and the continuous dependence of the model dynamics on its
parameters. We also propose a simple probabilistic method for the evolution of
the active cases that is instrumental for the automatic estimation of
parameters of the epidemiological model. We obtained statistical estimates of
the active cases based the probabilistic method and on the confirmed cases
data. From this estimated time series we obtained a time-dependent contagion
rate, which reflects a lower or higher adherence to social distancing by the
involved populations. By also analysing the data on daily deaths, we obtained
the daily lethality and recovery rates. We then integrate the equations of
motion of the model using these time-dependent parameters. We validate our
epidemiological model by fitting the official data of confirmed, recovered,
death, and active cases due to the pandemic with the theoretical predictions.
We obtained very good fits of the data with this method. The automated
procedure developed here could be used for basically any population with a
minimum of extra work. Finally, we also propose and validate a forecasting
method based on Markov chains for the evolution of the epidemiological data for
up to two weeks.


  In this work, we aim to study that the dynamics behavior for cumulative
number of SARS-CoV-2 pandemic can provide information on the overall behavior
of the spread over daily time.The cumulative data can be synthesized in an
empirical form obtained from a Susceptible-Infected-Recovered (SIR) model
previously studied on a Euclidean network. From the study we carried out, we
can conclude that the SIR model on the Euclidean network can reproduce data
from several countries with a deviation of precision for given parameter
values. This gives an insight into the different agents that influence the
behavior of SARS-CoV-2 especially during the virus mutation period. We are thus
trying to analyze the effect of genetic mutations in different countries, and
how a specific mutation can make the virus more contagious.


  Aim of this paper is to provide the first comprehensive synthesis about
ground beetle (Coleoptera: Carabidae) distribution in high altitude habitats.
Specifically, the attention is focused on the species assemblages living on the
most common ice-related mountain landforms (glaciers, debris-covered glaciers,
glacier forelands and rock glaciers) and on the challenges, threats and
opportunities carabids living in these habitats have to face in relation to the
ongoing climate warming. The suggested role of the ice-related alpine
landforms, as present climatic refugia for cold-adapted ground beetles, is
discussed. Finally, the needs to develop a large-scale High-alpine Biodiversity
Monitoring Program to describe how the current climate change is shaping the
distribution of high altitude specialists, is highlighted.


  Genetic information and environmental factors determine the path of an
individuals life and therefore, the evolution of its entire species. We have
succeeded in proposing and studying a model that captures this idea. In our
model, a renewable resource extended throughout the environment provides the
energy necessary to sustain life, including movement and reproduction. Since
the resource doesn't regrow immediately, it generates competition between
individuals and therefore provides a natural selection pressure from which
evolution of the genetic traits is observed. As a result of this, several
phenomena characteristic of living systems emerge from this model without
having to introduce them explicitly. These include speciation and punctuated
equilibrium, competitive exclusion, and altruistic behaviour from selfish
rules.


  A major challenge in the pharmaceutical industry is to design novel molecules
with specific desired properties, especially when the property evaluation is
costly. Here, we propose MNCE-RL, a graph convolutional policy network for
molecular optimization with molecular neighborhood-controlled embedding
grammars through reinforcement learning. We extend the original
neighborhood-controlled embedding grammars to make them applicable to molecular
graph generation and design an efficient algorithm to infer grammatical
production rules from given molecules. The use of grammars guarantees the
validity of the generated molecular structures. By transforming molecular
graphs to parse trees with the inferred grammars, the molecular structure
generation task is modeled as a Markov decision process where a policy gradient
strategy is utilized. In a series of experiments, we demonstrate that our
approach achieves state-of-the-art performance in a diverse range of molecular
optimization tasks and exhibits significant superiority in optimizing molecular
properties with a limited number of property evaluations.


  The SIR(D) epidemiological model is defined through a system of
transcendental equations, not solvable by elementary functions. In the present
paper those equations are successfully replaced by approximate ones, whose
solutions are given explicitly in terms of elementary functions, originating,
piece-wisely, from generalized logistic functions: they ensure {\em exact} (in
the numerical sense) asymptotic values, besides to be quite practical to use,
for example with fit to data algorithms; moreover they unveil a useful feature,
that in fact, at least with very strict approximation, is also owned by the
(numerical) solutions of the {\em exact} equations. The novelties in the work
are: the way the approximate equations are obtained, using simple, analytic
geometry considerations; the easy and practical formulation of the final
approximate solutions; the mentioned useful feature, never disclosed before.
The work's method and result prove to be robust over a range of values of the
well known non-dimensional parameter called {\em basic reproduction ratio},
that covers at least all the known epidemic cases, from influenza to measles:
this is a point which doesn't appear much discussed in analogous works.


  It is challenging to parameterize the force field for calcium ions (Ca2+) in
calcium-binding proteins because of their unique coordination chemistry that
involves the surrounding atoms required for stability. In this work, we
observed wide variation in Ca2+ binding loop conformations of the Ca2+-binding
protein calmodulin (CaM), which adopts the most populated ternary structures
determined from the MD simulations, followed by ab initio quantum mechanical
(QM) calculations on all twelve amino acids in the loop that coordinate Ca2+ in
aqueous solution. Ca2+ charges were derived by fitting to the electrostatic
potential (ESP) in the context of a classical or polarizable force field (PFF).
We discovered that the atomic radius of Ca2+ in conventional force fields is
too large for the QM calculation to capture the variation in the coordination
geometry of Ca2+ in its ionic form, leading to unphysical charges.
Specifically, we found that the fitted atomic charges of Ca2+ in the context of
PFF depend on the coordinating geometry of electronegative atoms from the amino
acids in the loop. Although nearby water molecules do not influence the atomic
charge of Ca2+, they are crucial for compensating for the coordination of Ca2+
due to the conformational flexibility in the EF-hand loop. Our method advances
the development of force fields for metal ions and protein binding sites in
dynamic environments.


  Mosquitoes are vectors of viral diseases with epidemic potential in many
regions of the world, and in absence of vaccines or therapies, their control is
the main alternative. Chemical control through insecticides has been one of the
conventional strategies, but induces insecticide resistance, which may affect
other insects and cause ecological damage. Biological control, through the
release of mosquitoes infected by the maternally inherited bacterium Wolbachia,
which inhibits their vector competence, has been proposed as an alternative.
The effects of both techniques may be intermingled in practice: prior
insecticide spraying may debilitate wild population, so facilitating subsequent
invasion by the bacterium; but the latter may also be hindered by the release
of susceptible mosquitoes in an environment where the wild population became
resistant, as a result of preexisting undesired exposition to insecticide. To
tackle such situations, we propose here a unifying model allowing to account
for the cross effects of both control techniques, and based on the latter,
design release strategies able to infect a wild population. The latter are
feedback laws, whose stabilizing properties are studied.


  The Susceptible-Infected-Recovered (SIR) model is the cornerstone of
epidemiological models. However, this specification depends on two parameters
only, which implies a lack of flexibility and the difficulty to replicate the
volatile reproduction numbers observed in practice. We extend the classic SIR
model by introducing nonlinear stochastic transmission, to get a stochastic SIR
model. We derive its exact solution and discuss the condition for herd
immunity. The stochastic SIR model corresponds to a population of infinite
size. When the population size is finite, there is also sampling uncertainty.
We propose a state-space framework under which we analyze the relative
magnitudes of the observational and stochastic epidemiological uncertainties
during the evolution of the epidemic. We also emphasize the lack of robustness
of the notion of herd immunity when the SIR model is time discretized.


  We introduce an asymmetric noisy voter model to study the joint effect of
immigration and a competition-dispersal tradeoff in the dynamics of two species
competing for space in regular lattices. Individuals of one species can invade
a nearest-neighbor site in the lattice, while individuals of the other species
are able to invade sites at any distance but are less competitive locally,
i.e., they establish with a probability $g \le 1$. The model also accounts for
immigration, modeled as an external noise that may spontaneously replace an
individual at a lattice site by another individual of the other species. This
combination of mechanisms gives rise to a rich variety of outcomes for species
competition, including exclusion of either species, mono-stable coexistence of
both species at different population proportions, and bi-stable coexistence
with proportions of populations that depend on the initial condition.
Remarkably, in the bi-stable phase, the system undergoes a discontinuous
transition as the intensity of immigration overcomes a threshold, leading to a
half loop dynamics associated to a cusp catastrophe, which causes the
irreversible loss of the species with the shortest dispersal range.


  How cells sense and respond to mechanical stimuli remains an open question.
Recent advances have identified the translocation of Yes-associated protein
(YAP) between nucleus and cytoplasm as a central mechanism for sensing
mechanical forces and regulating mechanotransduction. We formulate a
spatiotemporal model of the mechanotransduction signalling pathway that
includes coupling of YAP with the cell force-generation machinery through the
Rho family of GTPases. Considering the active and inactive forms of a single
Rho protein (GTP/GDP-bound) and of YAP (non-phosphorylated/phosphorylated), we
study the cross-talk between cell polarization due to active Rho and YAP
activation through its nuclear localization.
  For fixed mechanical stimuli, our model predicts stationary
nuclear-to-cytoplasmic YAP ratios consistent with experimental data at varying
adhesive cell area. We further predict damped and even sustained oscillations
in the YAP nuclear-to-cytoplasmic ratio by accounting for recently reported
positive and negative YAP-Rho feedback. Extending the framework to time-varying
mechanical stimuli that simulate cyclic stretching and compression, we show
that the YAP nuclear-to-cytoplasmic ratio's time dependence follows that of the
cyclic mechanical stimulus. The model presents one of the first frameworks for
understanding spatiotemporal YAP mechanotransduction, providing several
predictions of possible YAP localization dynamics, and suggesting new
directions for experimental and theoretical studies.


  We present Silmarillion, a novel, inclusive system for digital contact
tracing and epidemic risk notification, which simultaneously provides utility
as well as security. Silmarillion relies on a low-cost infrastructure of
strategically placed beacons, inexpensive and low-maintenance user devices like
dongles (if smartphones are inaccessible), and a backend that assists in
epidemiological analysis and risk dissemination. Unlike today's
smartphone-based contact tracing systems, Silmarillion records encounters
between users' devices and beacons installed in well-known and strategic
locations, which enables capturing contextual information relevant for
epidemiological analysis, individual risk prediction, and prevention of spread
of false risk information. Silmarillion keeps a user's encounter history local
to their device and allows the user to control what information they share with
the backend. Furthermore, it provides differential privacy for patients in risk
dissemination and information-theoretic privacy for users receiving the risk
information. We have evaluated a prototype of Silmarillion using small IoT
boards and show that the battery and speed of risk dissemination is adequate
for a practical deployment. Furthermore, we ran a small-scale deployment within
a university building, demonstrating Silmarillion's practicality.


  The spread of the coronavirus (COVID-19), starting in late 2019, has
determined in Italy several interventions aimed to prevent saturation of the
health system. We have examined the effects of such measures by proposing a
mean-field model describing the spread of the infection based on a simple
diffusion process where all the observable variables (number of people still
positive to the infection, hospitalized and fatalities cases, healed people,
and total number of people that has contracted the infection) depend on average
parameters, namely diffusion coefficient, infection cross-section, and
population density. Although this model is less sophisticated than other models
in the literature, it allows us to directly relate the trend of the epidemic
statistical information (hospitalized cases, number of fatalities, number of
infected people, etc.) to a well defined observable physical quantity: the
average number of people that any individual meets per day. The model fits very
well the epidemic data, and allows us to strictly relate the time evolution of
the number of hospitalized case and fatalities of the outbreak to the change of
people mobility, consequent to the implementation of progressive restrictions
in Italy, running until the present days (November the 15th, 2020).


  Longest Run Subsequence is a problem introduced recently in the context of
the scaffolding phase of genome assembly (Schrinner et al., WABI 2020). The
problem asks for a maximum length subsequence of a given string that contains
at most one run for each symbol (a run is a maximum substring of consecutive
identical symbols). The problem has been shown to be NP-hard and to be
fixed-parameter tractable when the parameter is the size of the alphabet on
which the input string is defined. In this paper we further investigate the
complexity of the problem and we show that it is fixed-parameter tractable when
it is parameterized by the number of runs in a solution, a smaller parameter.
Moreover, we investigate the kernelization complexity of Longest Run
Subsequence and we prove that it does not admit a polynomial kernel when
parameterized by the size of the alphabet or by the number of runs. Finally, we
consider the restriction of Longest Run Subsequence when each symbol has at
most two occurrences in the input string and we show that it is APX-hard.


  The blue crab (Callinectes sapidus) is a dominant ecological species of high
commercial value. Spawning stock and recruitment of the Chesapeake Bay
population declined by 80% in the 1990s. After severe management actions were
implemented in 2008, female abundance rebounded to pre-1994 levels and
stabilized. The stepwise decline in the early 1990s, followed by a consistently
low level of abundance for 15 y and a jump to high abundance after 2008,
suggested the existence of alternative stable states. Alternatively, high
fishing pressure combined with low recruitment in 1992 could have triggered a
proportional decline in the population, followed by a population increase in
2008 due to rigorous management actions that reduced fishing. We evaluated
these alternatives with a stage-structured dynamic population model using
ordinary differential equations. In addition, stock assessment models assume
that fishing and mortality are independent of density. Hence, we also
investigated the role of density-dependent predation, cannibalism and fishing
in blue crab population dynamics. We conclude that for the blue crab population
in Chesapeake Bay: (1) bistable positive states are not likely with
biologically realistic parameter values; (2) hyperbolic (depensatory) fishing
will not produce extinction at the range of population densities observed in
the bay; and (3) crabs can survive a higher fishing rate under the more
realistic assumption of sigmoidal (density-dependent) predation and cannibalism
than under constant (density-independent) predation and cannibalism. These
collectively indicate that the blue crab population in Chesapeake Bay is
resilient to a range of biotic and abiotic disturbances.


  Families are places of affection and cooperation, but also of conflict. In
his famous paper Parent-Offspring Conflict, Robert L. Trivers builds upon W. D.
Hamilton's concept of inclusive fitness to argue for genetic conflict in
parent-offspring relationships, and to derive numerical predictions on the
intensity of the conflict. We propose a mathematical model of game theory that
depicts how each member of a family allocates her resource budget to maximize
her inclusive fitness; this latter is made of the sum of personal fitness plus
the sum of relatives fitnesses weighted by Wright's coefficients of
relationship. We define an optimal allocation profile as a Nash equilibrium,
and we characterize the solutions in function of resource budgets, coefficients
of relationship and derivatives of personal fitnesses.


  Since the beginning of the corona pandemic in March 2020, various parameters
for describing the spread of the disease have been specified for Germany in
addition to the daily infection figures (new infections and total infections),
which are also used for political decisions. In addition to excess mortality
and the weekly incidence, these include the doubling time $T_2$ and the
reproduction number $R_t$. For the latter, various estimates can be found on
the website of the Robert-Koch-Institute, see \cite{EstR:RKI}, which are
calculated from the case numbers for all of Germany; local differences are not
taken into account here. In the present article, the calculations of the RKI on
a local level are examined using the example of northern Rhineland-Palatinate
and its districts. Here, not the reporting date but the onset of illness is
used as a reference for the calculation of $R_t$. For cases where the onset of
illness is not known, an adjusted generalized extreme value distribution (GEV)
is first fitted to the data for which the reporting delay (difference between
the onset of illness and the reporting date) is available and examined for
further characteristics such as local as well as demographic differences. This
GEV distribution is then used to calculate the reporting delays of incomplete
data points. The calculation of the daily value of $R_t$ between the end of
February and the end of October showed a similar course of the reproductive
rate compared to the nationwide figures. Expectably larger statistical
fluctuations were observed in the summer, mainly due to lower case numbers. The
values for northern Rhineland-Palatinate have been consistently above $1$ since
about mid-September. The calculations can also be transferred to other regions
and administrative districts.


  We propose an SEIR-type meta-population model to simulate and monitor the
Covid-19 epidemic evolution. The basic model consists of seven compartments,
namely susceptible (S), exposed (E), three infective classes, recovered (R),
and deceased (D). We define these compartments for n age and gender groups in m
different spatial locations. So, the resulting model has, for each age group,
gender, and place, all epidemiological classes. The mixing between them is
accomplished by means of time-dependent infection rate matrices. The model is
calibrated with the curve of daily new infections in New York City and its
boroughs, including census data, and the proportions of infections,
hospitalizations, and deaths for each age range. We end up with a model that
matches the reported curves and predicts accurately infection information for
different places and age classes.


  Intuitively, sampling is likely to be more efficient for prevalence
estimation, if the cases (or positives) have a relatively higher representation
in the sample than in the population. In case the virus is transmitted via
personal contacts, contact tracing of the observed cases (but not noncases), to
be referred to as \emph{adaptive network tracing}, can generate a higher yield
of cases than random sampling from the population. The efficacy of relevant
designs for cross-sectional and change estimation is investigated. The
availability of these designs allows one unite tracing for combating the
epidemic and sampling for estimating the prevalence in a single endeavour.


  We propose a compartmental mathematical model for the spread of the COVID-19
disease, showing its usefulness with respect to the pandemic in Portugal, from
the first recorded case in the country till the end of the three states of
emergency. New results include the compartmental model, described by a system
of seven ordinary differential equations; proof of positivity and boundedness
of solutions; investigation of equilibrium points and their stability analysis;
computation of the basic reproduction number; and numerical simulations with
official real data from the Portuguese health authorities. Besides completely
new, the proposed model allows to describe quite well the spread of COVID-19 in
Portugal, fitting simultaneously not only the number of active infected
individuals but also the number of hospitalized individuals, respectively with
a $L^2$ error of $9.2152e-04$ and $1.6136e-04$ with respect to the initial
population. Such results are very important, from a practical point of view,
and far from trivial from a mathematical perspective. Moreover, the obtained
value for the basic reproduction number is in agreement with the one given by
the Portuguese authorities at the end of the emergency states.


  We consider the qualitative behavior of a mathematical model for transmission
dynamics with two nonlinear stages of contagion. The proposed model is inspired
by phenomena occurring in epidemiology (spread of infectious diseases) or
social dynamics (spread of opinions, behaviors, ideas), and described by a
compartmental approach. Upon contact with a promoter (contagious individual), a
naive (susceptible) person can either become promoter himself or become
$\textit{weakened}$, hence more vulnerable. Weakened individuals become
contagious when they experience a second contact with members of the promoter
group. After a certain time in the contagious compartment, individuals become
inactive (are insusceptible and cannot spread) and are removed from the chain
of transmission. We combine this two-stage contagion process with renewal of
the naive population, modeled by means of transitions from the weakened or the
inactive status to the susceptible compartment. This leads to rich dynamics,
showing for instance coexistence and bistability of equilibria and periodic
orbits. Properties of (nontrivial) equilibria are studied analytically. In
addition, a numerical investigation of the parameter space reveals numerous
bifurcations, showing that the dynamics of such a system can be more complex
than those of classical epidemiological ODE models.


  The draft whole-genome sequence of the Japanese rhinoceros beetle, Trypoxylus
dichotomus was obtained using long-read PacBio sequence technology. The final
assembled genome consisted of 739 Mbp in 2,347 contigs, with 24.5x mean
coverage and a G+C content of 35.99%.


  COVID-19 is an infection caused by SARS-CoV-2 (Severe Acute Respiratory
Syndrome coronavirus 2), which has caused a global outbreak. Current research
efforts are focused on the understanding of the molecular mechanisms involved
in SARS-CoV-2 infection in order to propose drug-based therapeutic options.
Transcriptional changes due to epigenetic regulation are key host cell
responses to viral infection and have been studied in SARS-CoV and MERS-CoV;
however, such changes are not fully described for SARS-CoV-2. In this study, we
analyzed multiple transcriptomes obtained from cell lines infected with
MERS-CoV, SARS-CoV and SARS-CoV-2, and from COVID-19 patient-derived samples.
Using integrative analyses of gene co-expression networks and de-novo pathway
enrichment, we characterize different gene modules and protein pathways
enriched with Transcription Factors or Epifactors relevant for SARS-CoV-2
infection. We identified EP300, MOV10, RELA and TRIM25 as top candidates, and
more than 60 additional proteins involved in the epigenetic response during
viral infection that have therapeutic potential. Our results show that
targeting the epigenetic machinery could be a feasible alternative to treat
COVID-19.


  Emerging epidemics and local infection clusters are initially prone to
stochastic effects that can substantially impact the epidemic trajectory. While
numerous studies are devoted to the deterministic regime of an established
epidemic, mathematical descriptions of the initial phase of epidemic growth are
comparatively rarer. Here, we review existing mathematical results on the
epidemic size over time, and derive new results to elucidate the early dynamics
of an infection cluster started by a single infected individual. We show that
the initial growth of epidemics that eventually take off is accelerated by
stochasticity. These results are critical to improve early cluster detection
and control. As an application, we compute the distribution of the first
detection time of an infected individual in an infection cluster depending on
the testing effort, and estimate that the SARS-CoV-2 variant of concern Alpha
detected in September 2020 first appeared in the United Kingdom early August
2020. We also compute a minimal testing frequency to detect clusters before
they exceed a given threshold size. These results improve our theoretical
understanding of early epidemics and will be useful for the study and control
of local infectious disease clusters.


  We investigate the phenomenology emerging from a 2-species dynamics under the
scenario of a quasi-neutral competition within a metapopulation framework. We
employ stochastic and deterministic approaches, namely spatially-constrained
individual-based Monte Carlo simulations and coupled mean-field ODEs. Our
results show the multifold interplay between competition, birth-death dynamics
and spatial constraints induces a nonmonotonic relation between the ecological
majority-minority switching and the diffusion between patches. This means that
diffusion can set off birth-death ratios and enhance the preservation of a
species.


  Tipping points have been shown to be ubiquitous, both in models and
empirically in a range of physical and biological systems. The question of how
tipping points cascade through systems has been less well studied and is an
important one. A study of noise-induced tipping, in particular, could provide
key insights into tipping cascades. Here, we consider a specific example of a
simple model system that could have cascading tipping points. This model
consists of two interacting populations with underlying Allee effects and
stochastic dynamics, in separate patches connected by dispersal, which can
generate bistability. From an ecological standpoint, we look for rescue effects
whereby one population can prevent the collapse of a second population. As a
way to investigate the stochastic dynamics, we use an individual-based modeling
approach rooted in chemical reaction network theory. Then, using
continuous-time Markov chains and the theory of first passage times, we
essentially approximate, or emulate, the original high-dimensional model by a
Markov chain with just four states, where each state corresponds to a
combination of population thresholds. Analysis of this reduced model shows when
the system is likely to recover, as well as when tipping cascades through the
whole system.


  The COVID-19 pandemic has caused unprecedented disruption, particularly in
retail. Where essential demand cannot be fulfilled online, or where more
stringent measures have been relaxed, customers must visit shop premises in
person. This naturally gives rise to some risk of susceptible individuals
(customers or staff) becoming infected. It is essential to minimise this risk
as far as possible while retaining economic viability of the shop. We therefore
explore and compare the spread of COVID-19 in different shopping situations
involving person-to-person interactions: (i) free-flowing, unstructured
shopping; (ii) structured shopping (e.g. a queue). We examine which of (i) or
(ii) may be preferable for minimising the spread of COVID-19 in a given shop,
subject to constraints such as the geometry of the shop; compliance of the
population to local guidelines; and additional safety measures which may be
available to the organisers of the shop. We derive a series of conclusions,
such as unidirectional free movement being preferable to bidirectional
shopping, and that the number of servers should be maximised as long as they
can be well protected from infection.


  By the end of May 2020, all states in the US have eased their COVID-19
mitigation measures. Different states adopted markedly different policies and
timing for reopening. An important question remains in how the relaxation of
mitigation measures is related to the number of casualties. To address this
question, we compare the actual data to a hypothetical case in which the
mitigation measures are left intact using a projection of the data from before
mitigation measures were eased. We find that different states have shown
significant differences in the number of deaths, possibly due to their
different policies and reopening schedules. Our study provides a gauge for the
effectiveness of the approaches by different state governments and can serve as
a guide for implementing best policies in the future. It also indicates that
the face mask mandate seems to correlate with the change in the death count
more dramatically than other measures.


  In this paper, we develop an extension of standard epidemiological models,
suitable for COVID-19. This extension incorporates the transmission due to
pre-symptomatic or asymptomatic carriers of the virus. Furthermore, this model
also captures the spread of the disease due to the movement of people to/from
different administrative boundaries within a country. The model describes the
probabilistic rise in the number of confirmed cases due to the concomitant
effects of (incipient) human transmission and multiple compartments. The
associated parameters in the model can help architect the public health policy
and operational management of the pandemic. For instance, this model
demonstrates that increasing the testing for symptomatic patients does not have
any major effect on the progression of the pandemic, but testing rate of the
asymptomatic population has an extremely crucial role to play. The model is
executed using the data obtained for the state of Chhattisgarh in the Republic
of India. The model is shown to have significantly better predictive capability
than the other epidemiological models. This model can be readily applied to any
administrative boundary (state or country). Moreover, this model can be applied
for any other epidemic as well.


  Generally, carbohydrate-active enzymes are studied using chromogenic
substrates that provide quick and easy color-based detection of enzyme-mediated
hydrolysis. In the case of feruloyl esterases, commercially available
chromogenic ferulate derivatives are both costly and limited in terms of their
experimental application. In this study, we describe solutions for these two
issues, using a chemoenzymatic approach to synthesize different ferulate
compounds. The overall synthetic routes towards commercially available
5-bromo-4-chloro-3-indolyl and 4-nitrophenyl
O-5-feruloyl-$\alpha$-l-arabinofuranosides 1a and 1b were significantly
shortened (7-8 steps reduced to 4-6) and transesterification yields enhanced
(from 46 to 73% for 1a and 47 to 86 % for 1b). This was achieved using
enzymatic (immobilized Lipolase 100T from Thermomyces lanuginosus)
transesterification of unprotected vinyl ferulate to the primary hydroxyl group
of $\alpha$-l-arabinofuranosides. Moreover, a novel feruloylated-butanetriol
4-nitrocatechol-1-yl analog 12, containing a cleavable hydroxylated linker was
also synthesized in 29% overall yield in 3 steps (convergent synthesis). The
latter route combined regioselective functionalization of 4-nitrocatechol and
enzymatic transferuloylation. The use of 12 as a substrate to characterize type
A feruloyl esterase from Aspergillus niger reveals the advantages of this
substrate for the characterizations of feruloyl esterases.


  The spread of African swine fever (ASF) poses a grave threat to the global
swine industry. Understanding transmission dynamics, such as through
mechanistic modeling, is critical for designing effective control strategies.
Articles were examined across multiple epidemiological and model
characteristics. Model filiation was determined through creation of a
neighbor-joined tree using phylogenetic software. Of 34 four articles
qualifying for inclusion, four main modelling objectives were identified:
estimating transmission parameters (11 studies), assessing determinants of
transmission (7), examining consequences of hypothetical outbreaks (5), and
assessing alternative control strategies (11). Estimated transmission
parameters varied widely as did parameter assumptions between models.
Uncertainties on epidemiological and ecological parameters were usually
accounted for to assess the impact on the modeled infection trajectory. Almost
all models are host specific, being developed for either domestic pigs or wild
boar despite the fact that spillover events between domestic pigs and wild boar
are evidenced to play an important role in ASF outbreaks. The development of
models incorporating such transmission routes is crucial. All compared control
strategies were defined a priori, and future models should be built to identify
optimal contributions across many control methods. Further, control strategies
were examined in competition, opposed to how they would be synergistically
implemented. While comparing strategies is beneficial for identifying a
rank-order efficacy of control methods, this structure does not necessarily
determine the most effective combination of all available strategies. In order
for ASFV models to effectively support decision-making in controlling ASFV
globally, these modelling limitations need to be addressed.


  Efficient molecular featurization is one of the major issues for machine
learning models in drug design. Here we propose persistent Ricci curvature
(PRC), in particular Ollivier persistent Ricci curvature (OPRC), for the
molecular featurization and feature engineering, for the first time. Filtration
process proposed in persistent homology is employed to generate a series of
nested molecular graphs. Persistence and variation of Ollivier Ricci curvatures
on these nested graphs are defined as Ollivier persistent Ricci curvature.
Moreover, persistent attributes, which are statistical and combinatorial
properties of OPRCs during the filtration process, are used as molecular
descriptors, and further combined with machine learning models, in particular,
gradient boosting tree (GBT). Our OPRC-GBT model is used in the prediction of
protein-ligand binding affinity, which is one of key steps in drug design.
Based on three most-commonly used datasets from the well-established
protein-ligand binding databank, i.e., PDBbind, we intensively test our model
and compare with existing models. It has been found that our model are better
than all machine learning models with traditional molecular descriptors.


  Correlations between high life expectancy and low lifespan inequality are
frequently observed. A recent article seeks to explain this phenomenon by
proposing that a mortality improvement maps to life expectancy and relative
lifespan equality through the same weights w(x), with an extra weight W(x)
applied for relative lifespan equality, specific for the equality indicator
used. From statistical theory, the claim that changes of life expectancy and
lifespan (in)equality map through the same weights is unexpected.
  The current note explains this phenomenon by showing that W(x) undoes (part
of) w(x). Thus, while the change in relative lifespan equality is proportional
to the product of both weights, w(x)W(x), it is proportional to neither w(x)
nor W(x). As a result, some simplification is possible that gives way to an
intuitive understanding and a more direct interpretation of the perturbation
analysis proposed.