  COVID-19 infections have well described systemic manifestations, especially
respiratory problems. There are currently no specific treatments or vaccines
against the current strain. With higher case numbers, a range of neurological
symptoms are becoming apparent. The mechanisms responsible for these are not
well defined, other than those related to hypoxia and microthrombi. We
speculate that sustained systemic immune activation seen with SARS-CoV-2 may
also cause secondary autoimmune activation in the CNS. Patients with chronic
neurological diseases may be at higher risk because of chronic secondary
respiratory disease and potentially poor nutritional status. Here, we review
the impact of COVID-19 on people with chronic neurological diseases and
potential mechanisms. We believe special attention to protecting people with
neurodegenerative disease is warranted. We are concerned about a possible
delayed pandemic in the form of an increased burden of neurodegenerative
disease after acceleration of pathology by systemic COVID-19 infections.


  Since the epidemic began in November 2019, no viable medicine against
SARS-CoV-2 has been discovered. The typical medication discovery strategy
requires several years of rigorous research and development as well as a
significant financial commitment, which is not feasible in the face of the
current epidemic. Through molecular docking and dynamic simulation studies, we
used the FDA-approved drug mezonavir against the most important viral targets,
including spike (S) glycoprotein, Transmembrane serine protease 2 (TMPRSS2),
RNA-dependent RNA polymerase (RdRp), Main protease (Mpro), human
angiotensin-converting enzyme 2 (ACE-2), and furin. These targets are critical
for viral replication and infection propagation because they play a key role in
replication/transcription and host cell recognition. Molecular docking revealed
that the antiviral medication mozenavir showed a stronger affinity for
SARS-CoV-2 target proteins than reference medicines in this investigation. We
discovered that mozenavir increases the complex's stability and validates the
molecular docking findings using molecular dynamics modelling. Furin, a target
protein of COVID-19, has a greater binding affinity (-12.04 kcal/mol) than
other COVID-19 target proteins, forming different hydrogen bonds and polar and
hydrophobic interactions, suggesting that it might be used as an antiviral
treatment against SARS-CoV-2. Overall, the present in silico results will be
valuable in identifying crucial targets for subsequent experimental
investigations that might help combat COVID-19 by blocking the protease furin's
proteolytic activity.


  Machine-learning models that learn from data to predict how protein sequence
encodes function are emerging as a useful protein engineering tool. However,
when using these models to suggest new protein designs, one must deal with the
vast combinatorial complexity of protein sequences. Here, we review how to use
a sequence-to-function machine-learning surrogate model to select sequences for
experimental measurement. First, we discuss how to select sequences through a
single round of machine-learning optimization. Then, we discuss sequential
optimization, where the goal is to discover optimized sequences and improve the
model across multiple rounds of training, optimization, and experimental
measurement.


  The search for more realistic models for interacting species has produced
many adaptations of the original Lotka-Volterra equations, such as the
inclusion of the Allee effect and the different Holling's types of functional
response. In the present work we show that a correct implementation of both
ideas together requires a careful formulation. We focus our work in the fact
that a density dependent carrying capacity, combined with the Allee effect, can
lead to meaningless effects. We illustrate the difficulties in predator-prey
and two-species competition models, together with our proposed solution of the
careful inclusion of the corresponding cubic terms.


  Effective molecular representation learning is of great importance to
facilitate molecular property prediction, which is a fundamental task for the
drug and material industry. Recent advances in graph neural networks (GNNs)
have shown great promise in applying GNNs for molecular representation
learning. Moreover, a few recent studies have also demonstrated successful
applications of self-supervised learning methods to pre-train the GNNs to
overcome the problem of insufficient labeled molecules. However, existing GNNs
and pre-training strategies usually treat molecules as topological graph data
without fully utilizing the molecular geometry information. Whereas, the
three-dimensional (3D) spatial structure of a molecule, a.k.a molecular
geometry, is one of the most critical factors for determining molecular
physical, chemical, and biological properties. To this end, we propose a novel
Geometry Enhanced Molecular representation learning method (GEM) for Chemical
Representation Learning (ChemRL). At first, we design a geometry-based GNN
architecture that simultaneously models atoms, bonds, and bond angles in a
molecule. To be specific, we devised double graphs for a molecule: The first
one encodes the atom-bond relations; The second one encodes bond-angle
relations. Moreover, on top of the devised GNN architecture, we propose several
novel geometry-level self-supervised learning strategies to learn spatial
knowledge by utilizing the local and global molecular 3D structures. We compare
ChemRL-GEM with various state-of-the-art (SOTA) baselines on different
molecular benchmarks and exhibit that ChemRL-GEM can significantly outperform
all baselines in both regression and classification tasks. For example, the
experimental results show an overall improvement of 8.8% on average compared to
SOTA baselines on the regression tasks, demonstrating the superiority of the
proposed method.


  The robustness of an ecological network quantifies the resilience of the
ecosystem it represents to species loss. It corresponds to the proportion of
species that are disconnected from the rest of the network when extinctions
occur sequentially. Classically, the robustness is calculated for a given
network, from the simulation of a large number of extinction sequences. The
link between network structure and robustness remains an open question. Setting
a joint probabilistic model on the network and the extinction sequences allows
analysis of this relation.
  Bipartite stochastic block models have proven their ability to model
bipartite networks e.g. plant-pollinator networks: species are divided into
blocks and interaction probabilities are determined by the blocks of
membership. Analytical expressions of the expectation and variance of
robustness are obtained under this model, for different distributions of
primary extinction sequences. The impact of the network structure on the
robustness is analyzed through a set of properties and numerical illustrations.
The analysis of a collection of bipartite ecological networks allows us to
compare the empirical approach to our probabilistic approach, and illustrates
the relevance of the latter when it comes to computing the robustness of a
partially observed or incompletely sampled network.


  Strategies adopted globally to mitigate the threat of COVID-19 have primarily
involved lockdown measures with substantial economic and social costs with
varying degrees of success. Morbidity patterns of COVID-19 variants have a
strong association with age, while restrictive lockdown measures have
association with negative mental health outcomes in some age groups. Reduced
economic prospects may also afflict some age cohorts more than others.
Motivated by this, we propose a model to describe COVID-19 community spread
incorporating the role of age-specific social interactions. Through a flexible
parameterisation of an age-structured deterministic Susceptible Exposed
Infectious Removed (SEIR) model, we provide a means for characterising
different forms of lockdown which may impact specific age groups differently.
Social interactions are represented through age group to age group contact
matrices, which can be trained using available data and are thus locally
adapted. This framework is easy to interpret and suitable for describing
counterfactual scenarios, which could assist policy makers with regard to
minimising morbidity balanced with the costs of prospective suppression
strategies. Our work originates from an Irish context and we use disease
monitoring data from February 29th 2020 to January 31st 2021 gathered by Irish
governmental agencies. We demonstrate how Irish lockdown scenarios can be
constructed using the proposed model formulation and show results of
retrospective fitting to incidence rates and forward planning with relevant
``what if/instead of'' lockdown counterfactuals with uncertainty
quantification. Our formulation is agnostic to a specific locale, in that
lockdown strategies in other regions can be straightforwardly encoded using
this model. The methods we describe are made publicly available online through
an accessible and easy to use web interface.


  Transcription and translation retrieve and operationalize gene encoded
information in cells. These processes are not instantaneous and incur
significant delays. In this paper we study Goodwin models of both inducible and
repressible operons with state-dependent delays. The paper provides
justification and derivation of the model, detailed analysis of the appropriate
setting of the corresponding dynamical system, and extensive numerical analysis
of its dynamics. Comparison with constant delay models shows significant
differences in dynamics that include existence of stable periodic orbits in
inducible systems and multistability in repressible systems. A combination of
parameter space exploration, numerics, analysis of steady state linearization
and bifurcation theory indicates the likely presence of Shilnikov-type
homoclinic bifurcations in the repressible operon model.


  Motivation: Manual curation of genome-scale reconstructions is laborious, yet
existing automated curation tools typically do not take species-specific
experimental data and manually refined genome annotations into account.
Results: We developed DEMETER, a COBRA Toolbox extension that enables the
efficient simultaneous refinement of thousands of draft genome-scale
reconstructions while ensuring adherence to the quality standards in the field,
agreement with available experimental data, and refinement of pathways based on
manually refined genome annotations. Availability: DEMETER and tutorials are
available at https://github.com/opencobra/cobratoolbox.


  We present a new non-negative matrix factorization model for $(0,1)$
bounded-support data based on the doubly non-central beta (DNCB) distribution,
a generalization of the beta distribution. The expressiveness of the DNCB
distribution is particularly useful for modeling DNA methylation datasets,
which are typically highly dispersed and multi-modal; however, the model
structure is sufficiently general that it can be adapted to many other domains
where latent representations of $(0,1)$ bounded-support data are of interest.
Although the DNCB distribution lacks a closed-form conjugate prior, several
augmentations let us derive an efficient posterior inference algorithm composed
entirely of analytic updates. Our model improves out-of-sample predictive
performance on both real and synthetic DNA methylation datasets over
state-of-the-art methods in bioinformatics. In addition, our model yields
meaningful latent representations that accord with existing biological
knowledge.


  Phylogenetic networks represent evolutionary history of species and can
record natural reticulate evolutionary processes such as horizontal gene
transfer and gene recombination. This makes phylogenetic networks a more
comprehensive representation of evolutionary history compared to phylogenetic
trees. Stochastic processes for generating random trees or networks are
important tools in evolutionary analysis, especially in phylogeny
reconstruction where they can be utilized for validation or serve as priors for
Bayesian methods. However, as more network generators are developed, there is a
lack of discussion or comparison for different generators. To bridge this gap,
we compare a set of phylogenetic network generators by profiling topological
summary statistics of the generated networks over the number of reticulations
and comparing the topological profiles.


  Motivation: In spite of great success of genome-wide association studies
(GWAS), multiple challenges still remain. First, complex traits are often
associated with many single nucleotide polymorphisms (SNPs), each with small or
moderate effect sizes. Second, our understanding of the functional mechanisms
through which genetic variants are associated with complex traits is still
limited. To address these challenges, we propose GPA-Tree and it simultaneously
implements association mapping and identifies key combinations of functional
annotations related to risk-associated SNPs by combining a decision tree
algorithm with a hierarchical modeling framework. Results: First, we
implemented simulation studies to evaluate the proposed GPA-Tree method and
compared its performance with existing statistical approaches. The results
indicate that GPA-Tree outperforms existing statistical approaches in detecting
risk-associated SNPs and identifying the true combinations of functional
annotations with high accuracy. Second, we applied GPA-Tree to a systemic lupus
erythematosus (SLE) GWAS and functional annotation data including GenoSkyline
and GenoSkylinePlus. The results from GPA-Tree highlight the dysregulation of
blood immune cells, including but not limited to primary B, memory helper T,
regulatory T, neutrophils and CD8+ memory T cells in SLE. These results
demonstrate that GPA-Tree can be a powerful tool that improves association
mapping while facilitating understanding of the underlying genetic architecture
of complex traits and potential mechanisms linking risk-associated SNPs with
complex traits.


  One regards spaces of trees as stratified spaces, to study distributions of
phylogenetic trees. Stratified spaces with may have cycles, however spaces of
trees with a fixed number of leafs are contractible. Spaces of trees with three
leafs, in particular, are spiders with three legs. One gives an elementary
proof of the stickiness of intrinsic sample means on spiders. One also
represents four leafs tree data in terms of an associated Petersen graph. One
applies such ideas to analyze RNA sequences of SARS-CoV-2 from multiple
sources, by building samples of trees and running nonparametric statistics for
intrinsic means on tree spaces with three and four leafs. SARS-CoV-2 are also
used to built trees with leaves consisting in addition to other related
coronaviruses.


  Signaling pathways serve to communicate information about extracellular
conditions into the cell, to both the nucleus and cytoplasmic processes to
control cell responses. Genetic mutations in signaling network components are
frequently associated with cancer and can result in cells acquiring an ability
to divide and grow uncontrollably. Because signaling pathways play such a
significant role in cancer initiation and advancement, their constituent
proteins are attractive therapeutic targets. In this review, we discuss how
signaling pathway modeling can assist with identifying effective drugs for
treating diseases, such as cancer. An achievement that would facilitate the use
of such models is their ability to identify controlling biochemical parameters
in signaling pathways, such as molecular abundances and chemical reaction
rates, because this would help determine effective points of attack by
therapeutics.


  We propose a benchmark to study surrogate model accuracy for protein-ligand
docking. We share a dataset consisting of 200 million 3D complex structures and
2D structure scores across a consistent set of 13 million "in-stock" molecules
over 15 receptors, or binding sites, across the SARS-CoV-2 proteome. Our work
shows surrogate docking models have six orders of magnitude more throughput
than standard docking protocols on the same supercomputer node types. We
demonstrate the power of high-speed surrogate models by running each target
against 1 billion molecules in under a day (50k predictions per GPU seconds).
We showcase a workflow for docking utilizing surrogate ML models as a
pre-filter. Our workflow is ten times faster at screening a library of
compounds than the standard technique, with an error rate less than 0.01\% of
detecting the underlying best scoring 0.1\% of compounds. Our analysis of the
speedup explains that to screen more molecules under a docking paradigm,
another order of magnitude speedup must come from model accuracy rather than
computing speed (which, if increased, will not anymore alter our throughput to
screen molecules). We believe this is strong evidence for the community to
begin focusing on improving the accuracy of surrogate models to improve the
ability to screen massive compound libraries 100x or even 1000x faster than
current techniques.


  The biochemical reaction networks that regulate living systems are all
stochastic to varying degrees. The resulting randomness affects biological
outcomes at multiple scales, from the functional states of single proteins in a
cell to the evolutionary trajectory of whole populations. Controlling how the
distribution of these outcomes changes over time -- via external interventions
like time-varying concentrations of chemical species -- is a complex challenge.
In this work, we show how counterdiabatic (CD) driving, first developed to
control quantum systems, provides a versatile tool for steering biological
processes. We develop a practical graph-theoretic framework for CD driving in
discrete-state continuous-time Markov networks. Though CD driving is limited to
target trajectories that are instantaneous stationary states, we show how to
generalize the approach to allow for non-stationary targets and local control
-- where only a subset of system states are targeted. The latter is
particularly useful for biological implementations where there may be only a
small number of available external control knobs, insufficient for global
control. We derive simple graphical criteria for when local versus global
control is possible. Finally, we illustrate the formalism with global control
of a genetic regulatory switch and local control in chaperone-assisted protein
folding. The derived control protocols in the chaperone system closely resemble
natural control strategies seen in experimental measurements of heat shock
response in yeast and E. coli.


  We study species abundance in the empirical plant-pollinator mutualistic
networks exhibiting broad degree distributions, with uniform intra-group
competition assumed, by the Lotka-Volterra equation. The stability of a fixed
point is found to be identified by the signs of its non-zero components and
those of its neighboring fixed points. Taking the annealed approximation, we
derive the non-zero components to be formulated in terms of degrees and the
rescaled interaction strengths, which lead us to find different stable fixed
points depending on parameters, and we obtain the phase diagram. The selective
extinction phase finds small-degree species extinct and effective interaction
reduced, maintaining stability and hindering the onset of instability. The
non-zero minimum species abundances from different empirical networks show data
collapse when rescaled as predicted theoretically.


  In this paper we introduce a space-dependent multiscale model to describe the
spatial spread of an infectious disease under uncertain data with particular
interest in simulating the onset of the COVID-19 epidemic in Italy. While virus
transmission is ruled by a SEIAR type compartmental model, within our approach
the population is given by a sum of commuters moving on a extra-urban scale and
non commuters interacting only on the smaller urban scale. A transport dynamic
of the commuter population at large spatial scales, based on kinetic equations,
is coupled with a diffusion model for non commuters at the urban scale. Thanks
to a suitable scaling limit, the kinetic transport model used to describe the
dynamics of commuters, within a given urban area coincides with the diffusion
equations that characterize the movement of non-commuting individuals. Because
of the high uncertainty in the data reported in the early phase of the
epidemic, the presence of random inputs in both the initial data and the
epidemic parameters is included in the model. A robust numerical method is
designed to deal with the presence of multiple scales and the uncertainty
quantification process. In our simulations, we considered a realistic
geographical domain, describing the Lombardy region, in which the size of the
cities, the number of infected individuals, the average number of daily
commuters moving from one city to another, and the epidemic aspects are taken
into account through a calibration of the model parameters based on the actual
available data. The results show that the model is able to describe correctly
the main features of the spatial expansion of the first wave of COVID-19 in
northern Italy.


  Signal transduction is an important and basic mechanism to cell life
activities. The stochastic state transition of receptor induces the release of
signaling molecular, which triggers the state transition of other receptors. It
constructs a nonlinear sigaling network, and leads to robust switchlike
properties which are critical to biological function. Network architectures and
state transitions of receptor affect the performance of this biological
network. In this work, we perform a study of nonlinear signaling on biological
polymorphic network by analyzing network dynamics of the Ca$^{2+}$ induced
Ca$^{2+}$ release mechanism, where fast and slow processes are involved and the
receptor has four conformational states. Three types of networks,
Erd\"os-R\'enyi network, Watts-Strogatz network and BaraB\'asi-Albert network,
are considered with different parameters. The dynamics of the biological
networks exhibit different patterns at different time scales. At short time
scale, the second open state is essential to reproduce the quasi-bistable
regime, which emerges at a critical strength of connection for all three states
involved in the fast processes and disappears at another critical point. The
pattern at short time scale is not sensitive to the network architecture. At
long time scale, only monostable regime is observed, and difference of network
architectures affects the results more seriously. Our finding identifies
features of nonlinear signaling networks with multistate that may underlie
their biological function.


  The COVID-19 pandemic has initiated an unprecedented worldwide effort to
characterize its evolution through the mapping of mutations of the coronavirus
SARS-CoV-2. The early identification of mutations that could confer adaptive
advantages to the virus, such as higher infectivity or immune evasion, is of
paramount importance. However, the large number of currently available genomes
precludes the efficient use of phylogeny-based methods. Here we present
CoVtRec, a fast and scalable Topological Data Analysis approach for the
surveillance of emerging adaptive mutations in large genomic datasets. Our
method overcomes limitations of state-of-the-art phylogeny-based approaches by
quantifying the potential adaptiveness of mutations merely by their topological
footprint in the genome alignment, without resorting to the reconstruction of a
single optimal phylogenetic tree. Analyzing millions of SARS-CoV-2 genomes from
GISAID, we find a correlation between topological signals and adaptation to the
human host. By leveraging the stratification by time in sequence data, our
method enables the high-resolution longitudinal analysis of topological signals
of adaptation. We characterize the convergent evolution of the coronavirus
throughout the whole pandemic to date, report on emerging potentially adaptive
mutations, and pinpoint mutations in Variants of Concern that are likely
associated with positive selection. Our approach can improve the surveillance
of mutations of concern and guide experimental studies.


  Accurate modeling of the solvent environment for biological molecules is
crucial for computational biology and drug design. A popular approach to
achieve long simulation time scales for large system sizes is to incorporate
the effect of the solvent in a mean-field fashion with implicit solvent models.
However, a challenge with existing implicit solvent models is that they often
lack accuracy or certain physical properties compared to explicit solvent
models, as the many-body effects of the neglected solvent molecules is
difficult to model as a mean field. Here, we leverage machine learning (ML) and
multi-scale coarse graining (CG) in order to learn implicit solvent models that
can approximate the energetic and thermodynamic properties of a given explicit
solvent model with arbitrary accuracy, given enough training data. Following
the previous ML--CG models CGnet and CGSchnet, we introduce ISSNet, a graph
neural network, to model the implicit solvent potential of mean force. ISSNet
can learn from explicit solvent simulation data and be readily applied to MD
simulations. We compare the solute conformational distributions under different
solvation treatments for two peptide systems. The results indicate that ISSNet
models can outperform widely-used generalized Born and surface area models in
reproducing the thermodynamics of small protein systems with respect to
explicit solvent. The success of this novel method demonstrates the potential
benefit of applying machine learning methods in accurate modeling of solvent
effects for in silico research and biomedical applications.


  The structure of an RNA molecule plays a significant role in its biological
function. Predicting structure given a one dimensional sequence of RNA
nucleotide bases is a difficult and important problem. Many computer programs
(known as in silico) are available for predicting 2-dimensional (secondary)
structures however 3-dimensional (tertiary) structure prediction is much more
difficult mainly due to the far greater number of feasible solutions and fewer
experimental data on the thermodynamic energies of 3D structures. It is also
challenging to verify the most likely three dimensional structure even with the
availability of sophisticated x-ray crystallography and nuclear magnetic
resonance imaging technologies. In this paper we develop three dimensional RNA
folding predictions by adding penalty and reward parameters to a previous two
dimensional approach based on Quadratic Unconstrained Binary Optimization
(QUBO) models. These parameters provide flexibility in the amount of three
dimensional folding allowed. We address the problem of multiple near-optimal
structures via a new weighted similarity structure measure and illustrate
folding pathways via progressively improving local optimal solutions. The
problems are solved via a new commercial QUBO solver AlphaQUBO (Meta-Analytics,
2020) that solves problems having hundreds of thousands of binary variables.


  Mathematical models are widely recognized as an important tool for analyzing
and understanding the dynamics of infectious disease outbreaks, predict their
future trends, and evaluate public health intervention measures for disease
control and elimination. We propose a novel stochastic metapopulation
state-space model for COVID-19 transmission, based on a discrete-time
spatio-temporal susceptible/exposed/infected/recovered/deceased (SEIRD) model.
The proposed framework allows the hidden SEIRD states and unknown transmission
parameters to be estimated from noisy, incomplete time series of reported
epidemiological data, by application of unscented Kalman filtering (UKF),
maximum-likelihood adaptive filtering, and metaheuristic optimization.
Experiments using both synthetic data and real data from the Fall 2020 Covid-19
wave in the state of Texas demonstrate the effectiveness of the proposed model.


  During a long-running pandemic a pathogen can mutate, producing new strains
with different epidemiological parameters. Existing approaches to epidemic
modelling only consider one virus strain. We have developed a modified SEIR
model to simulate multiple virus strains within the same population. As a case
study, we investigate the potential effects of SARS-CoV-2 strain B.1.1.7 on the
city of Moscow. Our analysis indicates a high risk of a new wave of infections
in September-October 2021 with up to 35 000 daily infections at peak. We
open-source our code and data.


  Gene regulatory networks are comprised of biochemical reactions, which are
inherently stochastic. Each reaction channel contributes to this stochasticity
in different measure. In this paper we study the stochastic dynamics of protein
interaction networks (PIN) that are made up of monomers and dimers. The network
is defined by the dimers, which are formed by hybridizing two monomers. The
size of a PIN was defined as the number of monomers that interacts with at
least one other monomer (including itself). We generated 4200 random PIN of
sizes between 2 and 8 (600 per size) and simulated via the Gillespie algorithm
the stochastic evolution of copy numbers of all monomers and dimers until they
reached a steady state. The simulations revealed that the Fano factors of both
monomers and dimers in all networks and for all time points were close to one,
either from below or above. Only 10% of Fano factors for monomers were above
1.3 and 10% of Fano factors for dimers were above 1.17, with 5.54 and 3.47 as
the maximum value recorded for monomers and dimers, respectively. These
findings suggests that PIN in real biological setting contribute to the overall
stochastic noise that is close to Poisson. Our results also show a correlation
between stochastic noise, network size and network connectivity: for monomers,
the Fano factors tend towards 1 from above, while the Fano factors for dimers
tend towards 1 from below. For monomers, this tendency is amplified with
increased network connectivity.


  Summary: The CoV-Spectrum website supports the identification of new
SARS-CoV-2 variants of concern and the tracking of known variants. Its flexible
amino acid and nucleotide mutation search allows querying of variants before
they are designated by a lineage nomenclature system. The platform brings
together SARS-CoV-2 data from different sources and applies analyses. Results
include the proportion of different variants over time, their demographic and
geographic distributions, common mutations, hospitalization and mortality
probabilities, estimates for transmission fitness advantage and insights
obtained from wastewater samples.
  Availability and Implementation: CoV-Spectrum is available at
https://cov-spectrum.ethz.ch. The code is released under the GPL-3.0 license at
https://github.com/cevo-public/cov-spectrum-website


  In Southern Italy, since 2013, there has been an ongoing Olive Quick Decline
Syndrome (OQDS) outbreak, due to the bacterium Xylella fastidiosa. In a couple
of previous papers, the authors have proposed a mathematical approach for
identifying possible control strategies for eliminating or at least reduce the
economic impact of such event. The main players involved in OQDS are
represented by the insect vector, Philaenus spumarius, its host plants (olive
trees and weeds) and the bacterium, X. fastidiosa. A basic mathematical model
has been expressed in terms of a system of ordinary differential equations; a
preliminary analysis already provided interesting results about possible
control strategies within an integrated pest management framework, not
requiring the removal of the productive resource represented by the olive
trees. The same conjectures have been later confirmed by analyzing the impact
of possible spatial heterogeneities on controlling a X. fastidiosa epidemic.
These encouraging facts have stimulated a more detailed and rigorous
mathematical analysis of the same system, as presented in this paper. A clear
picture of the possible steady states (equilibria) and their stability
properties has been outlined, within a variety of different parameter
scenarios, for the original spatially homogeneous ecosystem. The results
obtained here confirm, in a mathematically rigorous way, what had been
conjectured in the previous papers, i.e. that the removal of a suitable amount
of weed biomass (reservoir of the juvenile stages of the insect vector of X.
fastidiosa from olive orchards and surrounding areas is the most acceptable
strategy to control the spread of the OQDS. In addition, as expected, the
adoption of more resistant olive tree cultivars has been shown to be a good
strategy, though less cost-effective, in controlling the pathogen.


  The cost of DNA sequencing has dropped exponentially over the past decade,
making genomic data accessible to a growing number of scientists. In
bioinformatics, localization of short DNA sequences (reads) within large
genomic sequences is commonly facilitated by constructing index data structures
which allow for efficient querying of substrings. Recent metagenomic
classification pipelines annotate reads with taxonomic labels by analyzing
their $k$-mer histograms with respect to a reference genome database. CPU-based
index construction is often performed in a preprocessing phase due to the
relatively high cost of building irregular data structures such as hash maps.
However, the rapidly growing amount of available reference genomes establishes
the need for index construction and querying at interactive speeds. In this
paper, we introduce MetaCache-GPU -- an ultra-fast metagenomic short read
classifier specifically tailored to fit the characteristics of CUDA-enabled
accelerators. Our approach employs a novel hash table variant featuring
efficient minhash fingerprinting of reads for locality-sensitive hashing and
their rapid insertion using warp-aggregated operations. Our performance
evaluation shows that MetaCache-GPU is able to build large reference databases
in a matter of seconds, enabling instantaneous operability, while popular
CPU-based tools such as Kraken2 require over an hour for index construction on
the same data. In the context of an ever-growing number of reference genomes,
MetaCache-GPU is the first metagenomic classifier that makes analysis pipelines
with on-demand composition of large-scale reference genome sets practical. The
source code is publicly available at https://github.com/muellan/metacache .


  Learning meaningful representations of data that can address challenges such
as batch effect correction and counterfactual inference is a central problem in
many domains including computational biology. Adopting a Conditional VAE
framework, we show that marginal independence between the representation and a
condition variable plays a key role in both of these challenges. We propose the
Contrastive Mixture of Posteriors (CoMP) method that uses a novel misalignment
penalty defined in terms of mixtures of the variational posteriors to enforce
this independence in latent space. We show that CoMP has attractive theoretical
properties compared to previous approaches, and we prove counterfactual
identifiability of CoMP under additional assumptions. We demonstrate
state-of-the-art performance on a set of challenging tasks including aligning
human tumour samples with cancer cell-lines, predicting transcriptome-level
perturbation responses, and batch correction on single-cell RNA sequencing
data. We also find parallels to fair representation learning and demonstrate
that CoMP is competitive on a common task in the field.


  We derive an alternative expression for a delayed logistic equation in which
the rate of change in the population involves a growth rate that depends on the
population density during an earlier time period. In our formulation, the delay
in the growth term is consistent with the rate of instantaneous decline in the
population given by the model. Our formulation is a modification of [Arino et
al., J.~Theoret.~Biol.~241(1):109--119, 2006] by taking the intraspecific
competition between the adults and juveniles into account. We provide a
complete global analysis showing that no sustained oscillations are possible. A
threshold giving the interface between extinction and survival is determined in
terms of the parameters in the model. The theory of chain transitive sets and
the comparison theorem for cooperative delay differential equations are used to
determine the global dynamics of the model.
  We extend our delayed logistic equation to a system modeling the competition
of two species. For the competition model, we provide results on local
stability, bifurcation diagrams, and adaptive dynamics. Assuming that the
species with shorter delay produces fewer offspring at a time than the species
with longer delay, we show that there is a critical value, $\tau^*$, such that
the evolutionary trend is for the delay to approach $\tau^*$.


  Sequencing costs currently prohibit the application of single-cell mRNA-seq
to many biological and clinical analyses. Targeted single-cell mRNA-sequencing
reduces sequencing costs by profiling reduced gene sets that capture biological
information with a minimal number of genes. Here, we introduce an active
learning method (ActiveSVM) that identifies minimal but highly-informative gene
sets that enable the identification of cell-types, physiological states, and
genetic perturbations in single-cell data using a small number of genes. Our
active feature selection procedure generates minimal gene sets from single-cell
data through an iterative cell-type classification task where misclassified
cells are examined at each round of analysis to identify maximally informative
genes through an `active' support vector machine (ActiveSVM) classifier. By
focusing computational resources on misclassified cells, ActiveSVM scales to
analyze data sets with over a million single cells. We demonstrate that
ActiveSVM feature selection identifies gene sets that enable ~90% cell-type
classification accuracy across a variety of data sets including cell atlas and
disease characterization data sets. The method generalizes to reveal genes that
respond to genetic perturbations and to identify region specific gene
expression patterns in spatial transcriptomics data. The discovery of small but
highly informative gene sets should enable substantial reductions in the number
of measurements necessary for application of single-cell mRNA-seq to clinical
tests, therapeutic discovery, and genetic screens.


  The world is still fighting against COVID-19, which has been lasting for more
than a year. Till date, it has been a greatest challenge to human beings in
fighting against COVID-19 since, the pathogen SARS-COV-2 that causes COVID-19
has significant biological and transmission characteristics when compared to
SARS-COV and MERS-COV pathogens. In spite of many control strategies that are
implemented to reduce the disease spread, there is a rise in the number of
infected cases around the world. Hence, a mathematical model which can describe
the real nature and impact of COVID-19 is necessary for the better
understanding of disease transmission dynamics of COVID-19. This article
proposes a new compartmental SEIRUC mathematical model, which includes the new
state called convalesce (C). The basic reproduction number $\mathcal{R}_0$ is
identified for the proposed model. The stability analysis are performed for the
disease free equilibrium ($\mathcal{E}_0$) as well for the endemic equilibrium
($\mathcal{E}_*$) by using the Routh-Hurwitz criterion. The graphical
illustrations of the proposed mathematical results are provided to validate the
theoretical results.


  The effect of polyvalent molecular cations, such as spermine, on the
condensation of DNA into very well-defined toroidal shapes have been well
studied and understood. However, a great effort has been made trying to obtain
similar condensed structures from either ssRNA or dsRNA, which the latter
carries similar negative charge density as dsDNA, although it adopts a
different helical form. But the analogous condensation of RNA molecules into
well-defined structures has so far been elusive. In this work, we show that
ssRNA molecules can easily be condensed into nanoring structures on a mice
surface, where each nanoring structure is formed mostly by a single RNA
molecule. The condensation occurs in a concentration range of different atomic
cations, from monovalent to trivalent. The structures of the RNA nanorings on
mica surfaces were oberved by atomic force microscopy (AFM). The samples were
observed in tapping mode and were prepared by drop evaporation of a solution of
RNA in the presence of one type of the different cations used. As far as we
know, this is the first time that nanorings or any other well-defined condensed
RNA structures have been reported. The RNA nanorings formation can be
understood by an energy competition between the hydrogen bonding forming
hairpin stems, weakened by the salts, and hairpin loops. This results may have
an important biological relevance, since it has been proposed that RNA is the
oldest genome coding molecule and the formation of these structures could have
given it stability against degradation in primeval times. Even more, the
nanoring structures could have the potential to be used as biosensors and
functionalized nanodevices.


  A highly organized and densely packed lattice of molecular machinery within
the sarcomeres of muscle cells powers contraction. Although many of the
proteins that drive contraction have been studied extensively, the mechanical
impact of fluid shearing within the lattice of molecular machinery has received
minimal attention. It was recently proposed that fluid flow augments substrate
transport in the sarcomere, however, this analysis used analytical models of
fluid flow in the molecular machinery that could not capture its full
complexity. By building a finite element model of the sarcomere, we estimate
the explicit flow field, and contrast it with analytical models. Our results
demonstrate that viscous drag forces on sliding filaments are surprisingly
small in contrast to the forces generated by single myosin molecular motors.
This model also indicates that the energetic cost of fluid flow through viscous
shearing with lattice proteins is likely minimal. The model also highlights a
steep velocity gradient between sliding filaments and demonstrates that the
maximal radial fluid velocity occurs near the tips of the filaments. To our
knowledge, this is the first computational analysis of fluid flow within the
highly structured sarcomere.


  There are multiple testing methods to ascertain an infection in an individual
and they vary in their performances, cost and delay. Unfortunately, better
performing tests are sometimes costlier and time consuming and can only be done
for a small fraction of the population. On the other hand, greater number of
individuals can be tested using a cheaper, rapid test, but may only provide
less reliable results. In this work, we studied the interplay between cost and
delay of the tests as well the additional advantages offered by partial and
complete lockdowns. To understand the influence of different test strategies,
we implemented them on realistic random social networks with a COVID-19-like
epidemic in progression. Specifically, we compared the performance of two tests
mimicking the characteristics of popular tests implemented for COVID-19
detection. We present procedures and intuitive understanding to ascertain the
optimum combination of the tests to minimize the peak infection as well as
total quarantine days when the number of tests is constrained by a fixed total
budget.


  The COVID-19 pandemic has influenced the lives of people globally. In the
past year many researchers have proposed different models and approaches to
explore in what ways the spread of the disease could be mitigated. One of the
models that have been used a great deal is the
Susceptible-Exposed-Infectious-Recovered (SEIR) model. Some researchers have
modified the traditional SEIR model, and proposed new versions of it. However,
to the best of our knowledge, the state-of-the-art papers have not considered
the effect of different vaccine types, meaning single shot and double shot
vaccines, in their SEIR model. In this paper, we propose a modified version of
the SEIR model which takes into account the effect of different vaccine types.
We compare how different policies for the administration of the vaccine can
influence the rate at which people are exposed to the disease, get infected,
recover, and pass away. Our results suggest that taking the double shot vaccine
such as Pfizer-BioNTech and Moderna does a better job at mitigating the spread
and fatality rate of the disease compared to the single shot vaccine, due to
its higher efficacy.


  The dynamics of epidemics depend on how people's behavior changes during an
outbreak. At the beginning of the epidemic, people do not know about the virus,
then, after the outbreak of epidemics and alarm, they begin to comply with the
restrictions and the spreading of epidemics may decline. Over time, some people
get tired/frustrated by the restrictions and stop following them (exhaustion),
especially if the number of new cases drops down. After resting for a while,
they can follow the restrictions again. But during this pause the second wave
can come and become even stronger then the first one. Studies based on SIR
models do not predict the observed quick exit from the first wave of epidemics.
Social dynamics should be considered. The appearance of the second wave also
depends on social factors. Many generalizations of the SIR model have been
developed that take into account the weakening of immunity over time, the
evolution of the virus, vaccination and other medical and biological details.
However, these more sophisticated models do not explain the apparent
differences in outbreak profiles between countries with different intrinsic
socio-cultural features. In our work, a system of models of the COVID-19
pandemic is proposed, combining the dynamics of social stress with classical
epidemic models. Social stress is described by the tools of sociophysics. The
combination of a dynamic SIR-type model with the classical triad of stages of
the general adaptation syndrome, alarm-resistance-exhaustion, makes it possible
to describe with high accuracy the available statistical data for 13 countries.
The sets of kinetic constants corresponding to optimal fit of model to data
were found. They characterize the ability of society to mobilize efforts
against epidemics and maintain this concentration over time, and can further
help in the development of strategies specific to a particular society.


  Population dynamics and evolutionary genetics underly the structure of
ecosystems, changing on the same timescale for interacting species with rapid
turnover, such as virus (e.g. HIV) and immune response. Thus, an important
problem in mathematical modeling is to connect ecology, evolution and genetics,
which often have been treated separately. Here, extending analysis of multiple
virus and immune response populations in a resource - prey (consumer) -
predator model from Browne and Smith \cite{browne2018dynamics}, we show that
long term dynamics of viral mutants evolving resistance at distinct epitopes
(viral proteins targeted by immune responses) are governed by epistasis in the
virus fitness landscape. In particular, the stability of persistent equilibrium
virus-immune (prey-predator) network structures, such as nested and one-to-one,
and bifurcations are determined by a collection of circuits defined by
combinations of viral fitnesses that are minimally additive within a hypercube
of binary sequences representing all possible viral epitope sequences ordered
according to immunodominance hierarchy. Numerical solutions of our ordinary
differential equation system, along with an extended stochastic version
including random mutation, demonstrate how pairwise or multiplicative epistatic
interactions shape viral evolution against concurrent immune responses and
convergence to the multi-variant steady state predicted by theoretical results.
Furthermore, simulations illustrate how periodic infusions of subdominant
immune responses can induce a bifurcation in the persistent viral strains,
offering superior host outcome over an alternative strategy of immunotherapy
with strongest immune response.


  Several mathematical models in SARS-CoV-2 have shown how target-cell model
can help to understand the spread of the virus in the host and how potential
candidates of antiviral treatments can help to control the virus. Concepts as
equilibrium and stability show to be crucial to qualitatively determine the
best alternatives to schedule drugs, according to effectivity in inhibiting the
virus infection and replication rates. Important biological events such as
rebounds of the infections (when antivirals are incorrectly interrupted) can
also be explained by means of a dynamic study of the target-cell model. In this
work, a full characterization of the dynamical behavior of the target-cell
models under control actions is made and, based on this characterization, the
optimal fixed-dose antiviral schedule that produces the smallest amount of dead
cells (without viral load rebounds) is computed. Several simulation results -
performed by considering real patient data - show the potential benefits of
both, the model characterization and the control strategy.


  Models based on machine learning can enable accurate and fast molecular
property predictions, which is of interest in drug discovery and material
design. Various supervised machine learning models have demonstrated promising
performance, but the vast chemical space and the limited availability of
property labels make supervised learning challenging. Recently, unsupervised
transformer-based language models pretrained on a large unlabelled corpus have
produced state-of-the-art results in many downstream natural language
processing tasks. Inspired by this development, we present molecular embeddings
obtained by training an efficient transformer encoder model, MoLFormer, which
uses rotary positional embeddings. This model employs a linear attention
mechanism, coupled with highly distributed training, on SMILES sequences of 1.1
billion unlabelled molecules from the PubChem and ZINC datasets. We show that
the learned molecular representation outperforms existing baselines, including
supervised and self-supervised graph neural networks and language models, on
several downstream tasks from ten benchmark datasets. They perform
competitively on two others. Further analyses, specifically through the lens of
attention, demonstrate that MoLFormer trained on chemical SMILES indeed learns
the spatial relationships between atoms within a molecule. These results
provide encouraging evidence that large-scale molecular language models can
capture sufficient chemical and structural information to predict various
distinct molecular properties, including quantum-chemical properties.


  Of the many modern approaches to calculating evolutionary distance via models
of genome rearrangement, most are tied to a particular set of genomic modelling
assumptions and to a restricted class of allowed rearrangements. The "position
paradigm", in which genomes are represented as permutations signifying the
position (and orientation) of each region, enables a refined model-based
approach, where one can select biologically plausible rearrangements and assign
to them relative probabilities/costs. Here, one must further incorporate any
underlying structural symmetry of the genomes into the calculations and ensure
that this symmetry is reflected in the model. In our recently-introduced
framework of {\em genome algebras}, each genome corresponds to an element that
simultaneously incorporates all of its inherent physical symmetries. The
representation theory of these algebras then provides a natural model of
evolution via rearrangement as a Markov chain. Whilst the implementation of
this framework to calculate distances for genomes with `practical' numbers of
regions is currently computationally infeasible, we consider it to be a
significant theoretical advance: one can incorporate different genomic
modelling assumptions, calculate various genomic distances, and compare the
results under different rearrangement models. The aim of this paper is to
demonstrate some of these features.


  The language commonly used in human genetics can inadvertently pose problems
for multiple reasons. Terms like "ancestry", "ethnicity", and other ways of
grouping people can have complex, often poorly understood, or multiple meanings
within the various fields of genetics, between different domains of biological
sciences and medicine, and between scientists and the general public.
Furthermore, some categories in frequently used datasets carry scientifically
misleading, outmoded or even racist perspectives derived from the history of
science. Here, we discuss examples of problematic lexicon in genetics, and how
commonly used statistical practices to control for the non-genetic environment
may exacerbate difficulties in our terminology, and therefore understanding.
Our intention is to stimulate a much-needed discussion about the language of
genetics, to begin a process to clarify existing terminology, and in some cases
adopt a new lexicon that both serves scientific insight, and cuts us loose from
various aspects of a pernicious past.


  Here, we discuss a collection of cutting-edge techniques and applications in
use today by some of the leading experts in the field of correlative approaches
in single-molecule biophysics. A key difference in emphasis, compared with
traditional single-molecule biophysics approaches detailed previously, is on
the emphasis of the development and use of complex methods which explicitly
combine multiple approaches to increase biological insights at the
single-molecule level. These so-called correlative single-molecule biophysics
methods rely on multiple, orthogonal tools and analysis, as opposed to any one
single driving technique. Importantly, they span both in vivo and in vitro
biological systems as well as the interfaces between theory and experiment in
often highly integrated ways, very different to earlier traditional
non-integrative approaches. The first applications of correlative
single-molecule methods involved adaption of a range of different experimental
technologies to the same biological sample whose measurements were
synchronised. However, now we find a greater flora of integrated methods
emerging that include approaches applied to different samples at different
times and yet still permit useful molecular-scale correlations to be performed.
The resultant findings often enable far greater precision of length and time
scales of measurements, and a more understanding of the interplay between
different processes in the same cell. Many new correlative single-molecule
biophysics techniques also include more complex, physiologically relevant
approaches as well as increasing number that combine advanced computational
methods and mathematical analysis with experimental tools. Here we review the
motivation behind the development of correlative single-molecule microscopy
methods, its history and recent progress in the field.


  Despite commercial airlines mandating masks, there have been multiple
documented events of COVID-19 superspreading on flights. We used available data
from three flights, including cabin layout and seat locations of infected and
uninfected passengers, to suggest interventions to mitigate COVID-19
superspreading events during air travel. Specifically, we studied: 1) London to
Hanoi with 201 passengers, including 13 secondary infections among passengers;
2) Singapore to Hangzhou with 321 passengers, including 12 to 14 secondary
infections; 3) a non-superspreading event on a private jet in Japan with 9
passengers and no secondary infections. We show that inclusion of passenger
movement better explains the infection spread patterns than conventional models
do. We also found that FFP2/N95 mask usage would have reduced infection by
95-100%, while cloth masks would have reduced it by only 40-80%. This suggests
the need for more stringent guidelines to reduce aviation-related
superspreading events of COVID-19.


  In contrast to the common assumption in epidemic models that the rate of
infection between individuals is constant, in reality, an individual's viral
load determines their infectiousness. We compare the average and individual
reproductive numbers and epidemic dynamics for a model incorporating
time-dependent infectiousness and a standard SIR model for both fully-mixed and
category-mixed populations. We find that the reproductive number only depends
on the total infectious exposure and the largest eigenvalue of the mixing
matrix and that these two effects are independent of each other. When we
compare our time-dependent mean-field model to the SIR model with equivalent
rates, the epidemic peak is advanced and modifying the infection rate function
has a strong effect on the time dynamics of the epidemic. We also observe
behavior akin to a traveling wave as individuals transition through infectious
states.


  Within epidemiological modeling, the majority of analyses assume a single
epidemic process for generating ground-truth data. However, this assumed data
generation process can be unrealistic, since data sources for epidemics are
often aggregated across geographic regions and communities. As a result,
state-of-the-art models for estimating epidemiological parameters,
e.g.~transmission rates, can be inappropriate when faced with complex systems.
Our work empirically demonstrates some limitations of applying epidemiological
models to aggregated datasets. We generate three complex outbreak scenarios by
combining incidence curves from multiple epidemics that are independently
simulated via SEIR models with different sets of parameters. Using these
scenarios, we assess the robustness of a state-of-the-art Bayesian inference
method that estimates the epidemic trajectory from viral load surveillance
data. We evaluate two data-generating models within this Bayesian inference
framework: a simple exponential growth model and a highly flexible Gaussian
process prior model. Our results show that both models generate accurate
transmission rate estimates for the combined incidence curve at the cost of
generating biased estimates for each underlying epidemic, reflecting highly
heterogeneous underlying population dynamics. The exponential growth model,
while interpretable, is unable to capture the complexity of the underlying
epidemics. With sufficient surveillance data, the Gaussian process prior model
captures the shape of complex trajectories, but is imprecise for periods of low
data coverage. Thus, our results highlight the potential pitfalls of neglecting
complexity and heterogeneity in the data generation process, which can mask
underlying location- and population-specific epidemic dynamics.


  In the context of a pandemic like COVID-19, and until most people are
vaccinated, proactive testing and interventions have been proved to be the only
means to contain the disease spread. Recent academic work has offered
significant evidence in this regard, but a critical question is still open: Can
we accurately identify all new infections that happen every day, without this
being forbiddingly expensive, i.e., using only a fraction of the tests needed
to test everyone everyday (complete testing)? Group testing offers a powerful
toolset for minimizing the number of tests, but it does not account for the
time dynamics behind the infections. Moreover, it typically assumes that people
are infected independently, while infections are governed by community spread.
Epidemiology, on the other hand, does explore time dynamics and community
correlations through the well-established continuous-time SIR stochastic
network model, but the standard model does not incorporate discrete-time
testing and interventions. In this paper, we introduce a "discrete-time SIR
stochastic block model" that also allows for group testing and interventions on
a daily basis. Our model can be regarded as a discrete version of the
continuous-time SIR stochastic network model over a specific type of weighted
graph that captures the underlying community structure. We analyze that model
w.r.t. the minimum number of group tests needed everyday to identify all
infections with vanishing error probability. We find that one can leverage the
knowledge of the community and the model to inform nonadaptive group testing
algorithms that are order-optimal, and therefore achieve the same performance
as complete testing using a much smaller number of tests.


  This paper proposes a model to explain the potential role of inter-group
conflicts in determining the rise and fall of signaling norms. Individuals in a
population are characterized by high and low productivity types and they are
matched in pairs to form social relationships such as mating or foraging
relationships. In each relationship, an individual's payoff is increasing in
its own type and its partner's type. Hence, the payoff structure of a
relationship does not resemble a dilemma situation. Assume that types are not
observable. In one population, assortative matching according to types is
sustained by signaling. In the other population, individuals do not signal and
they are randomly matched. Types evolve within each population. At the same
time, the two populations may engage in conflicts. Due to assortative matching,
high types grow faster in the population with signaling, yet they bear the cost
of signaling, which lowers their population's fitness in the long run. Through
simulations, we show that the survival of the signaling population depends
crucially on the timing and the efficiency of weapon used in inter-group
conflicts.


  Understanding dynamics of an outbreak like that of COVID-19 is important in
designing effective control measures. This study aims to develop an agent based
model that compares changes in infection progression by manipulating different
parameters in a synthetic population. Model input includes population
characteristics like age, sex, working status etc. of each individual and other
factors influencing disease dynamics. Depending on number of epicentres of
infection, location of primary cases, sensitivity, proportion of asymptomatic
and frequency or duration of lockdown, our simulator tracks every individual
and hence infection progression through community over time.
  In a closed community of 10000 people, it is seen that without any lockdown,
number of cases peak around 6th week and wanes off around 15th week. If primary
case is located inside dense population cluster like slums, cases peak early
and wane off slowly. With introduction of lockdown, cases peak at slower rate.
If sensitivity of identifying infection decreases, cases and deaths increase.
Number of cases declines with increase in proportion of asymptomatic cases. The
model is robust and provides reproducible estimates with realistic parameter
values. It also guides in identifying measures to control outbreak in a
community. It is flexible in accommodating different parameters like
infectivity period, yield of testing, socio-economic strata, daily travel,
awareness level, population density, social distancing, lockdown etc. and can
be tailored to study other infections with similar transmission pattern.


  In order to slow the spread of the CoViD-19 pandemic, governments around the
world have enacted a wide set of policies limiting the transmission of the
disease. Initially, these focused on non-pharmaceutical interventions; more
recently, vaccinations and large-scale rapid testing have started to play a
major role. The objective of this study is to explain the quantitative effects
of these policies on determining the course of the pandemic, allowing for
factors like seasonality or virus strains with different transmission profiles.
To do so, the study develops an agent-based simulation model, which is
estimated using data for the second and the third wave of the CoViD-19 pandemic
in Germany. The paper finds that during a period where vaccination rates rose
from 5% to 40%, large-scale rapid testing had the largest effect on reducing
infection numbers. Frequent large-scale rapid testing should remain part of
strategies to contain CoViD-19; it can substitute for many non-pharmaceutical
interventions that come at a much larger cost to individuals, society, and the
economy.


  To protect the location of camera trap data containing sensitive, high-target
species, many ecologists randomly obfuscate the latitude and longitude of the
camera when publishing their data. For example, they may publish a random
location within a 1km radius of the true camera location for each camera in
their network. In this paper, we investigate the robustness of geo-obfuscation
for maintaining camera trap location privacy, and show via a case study that a
few simple, intuitive heuristics and publicly available satellite rasters can
be used to reduce the area likely to contain the camera by 87% (assuming random
obfuscation within 1km), demonstrating that geo-obfuscation may be less
effective than previously believed.


  Closed-form, analytical approximations for electrostatic properties of
molecules are of unique value, as these can provide computational speed,
versatility, and physical insight. Here, we derive a simple, closed-form
formula for the apparent surface charge (ASC), as well as for the electric
field, generated by a molecular charge distribution in aqueous solution. The
approximation, with no fitted parameters, is tested against numerical solutions
of the Poisson equation, where it yields a significant speed-up. For neutral
small molecules, the hydration free energies estimated from the closed-form ASC
formula are within 0.8 kcal/mol RMSD from the numerical Poisson reference; the
electric field at the surface is in quantitative agreement with the reference.
Performance of the approximation is also tested on larger structures, including
a protein, a DNA fragment, and a viral receptor-target complex. For all
structures tested, a near quantitative agreement with the numerical Poisson
reference is achieved, except in regions of high negative curvature, where the
new approximation is still qualitatively correct. A unique efficiency feature
of the proposed "source-based" closed-form approximation is that the ASC and
electric field can be estimated individually at any point or surface patch,
without the need to obtain the full global solution. An open source software
implementation of the method is available:
http://people.cs.vt.edu/~onufriev/CODES/aasc.zip


  Reversible covalent kinase inhibitors (RCKIs) are a class of novel kinase
inhibitors attracting increasing attention because they simultaneously show the
selectivity of covalent kinase inhibitors, yet avoid permanent
protein-modification-induced adverse effects. Over the last decade, RCKIs have
been reported to target different kinases, including atypical kinases.
Currently, three RCKIs are undergoing clinical trials to treat specific
diseases, for example, Pemphigus, an autoimmune disorder. In this perspective,
first, RCKIs are systematically summarized, including characteristics of
electrophilic groups, chemical scaffolds, nucleophilic residues, and binding
modes. Second, we provide insights into privileged electrophiles, the
distribution of nucleophiles and hence effective design strategies for RCKIs.
Finally, we provide a brief perspective on future design strategies for RCKIs,
including those that target proteins other than kinases.


  In this paper we introduce a compartmental epidemic model describing the
transmission of the COVID-19 disease in presence of non-mandatory vaccination.
The model takes into account the hesitancy and refusal of vaccination. To this
aim, we employ the information index, which mimics the idea that individuals
take their decision on vaccination based not only on the present but also on
the past information about the spread of the disease. Theoretical analysis and
simulations show clearly as a voluntary vaccination can certainly reduce the
impact of the disease but it is unable to eliminate it. We also show how the
information-related parameters affect the dynamics of the disease. In
particular, the hesitancy and refusal of vaccination is better contained in
case of large information coverage and small memory characteristic time.
Finally, the possible influence of seasonality is also investigated.


  Methods for evaluating the quality of genomic and metagenomic data are
essential to aid genome assembly and to correctly interpret the results of
subsequent analyses. BUSCO estimates the completeness and redundancy of
processed genomic data based on universal single-copy orthologs. Here we
present new functionalities and major improvements of the BUSCO software, as
well as the renewal and expansion of the underlying datasets in sync with the
OrthoDB v10 release. Among the major novelties, BUSCO now enables phylogenetic
placement of the input sequence to automatically select the most appropriate
dataset for the assessment, allowing the analysis of metagenome-assembled
genomes of unknown origin. A newly-introduced genome workflow increases the
efficiency and runtimes especially on large eukaryotic genomes. BUSCO is the
only tool capable of assessing both eukaryotic and prokaryotic species, and can
be applied to various data types, from genome assemblies and metagenomic bins,
to transcriptomes and gene sets.


  We propose a novel approach for modelling bat motion dynamics and use it to
predict roost locations using data from static acoustic detectors.
Specifically, radio tracking studies of Greater Horseshoe bats demonstrate that
bat movement can be split into two phases: dispersion and return. Dispersion is
easily understood and can be modelled as simple random motion. The return phase
is much more complex, as it requires intelligent directed motion and results in
all agents returning home in a stereotypical manner. Critically, combining
reaction-diffusion theory and domain shrinking we deterministically and
stochastically model a ``leap-frogging'' motion, which fits favourably with the
observed tracking data.


  We introduce and analyze an extension of the mixed feedback loop model of
Fran\c{c}ois and Hakim. Our extension includes an additional promoter state and
allows for reversible protein sequestration, which was absent from the original
studies of the mixed feedback loop model. Motivated by experimental
observations that link DNA methylation with circadian gene expression, we use
our extended model to investigate the role of DNA methylation(s) in the
mammalian circadian clock. We extend the perturbation analysis of Fran\c{c}ois
and Hakim to determine how methylation affects the presence and the periodicity
of oscillations. We derive a modified Goodwin oscillator model as an
approximation to show that although methylation contributes to period control,
excessive methylation can abolish rhythmicity.


  We study the spread of information on multi-type directed random graphs. In
such graphs the vertices are partitioned into distinct types (communities) that
have different transmission rates between themselves and with other types. We
construct multivariate generating functions and use multi-type branching
processes to derive an equation for the size of the large out-components in
multi-type random graphs with a general class of degree distributions. We use
our methods to analyse the spread of epidemics and verify the results with
population based simulations


  The expected number of secondary infections arising from each index case,
referred to as the reproduction or $R$ number, is a vital summary statistic for
understanding and managing epidemic diseases. There are many methods for
estimating $R$; however, few explicitly model heterogeneous disease
reproduction, which gives rise to superspreading within the population. We
propose a parsimonious discrete-time branching process model for epidemic
curves that incorporates heterogeneous individual reproduction numbers. Our
Bayesian approach to inference illustrates that this heterogeneity results in
less certainty on estimates of the time-varying cohort reproduction number
$R_t$. We apply these methods to a COVID-19 epidemic curve for the Republic of
Ireland and find support for heterogeneous disease reproduction. Our analysis
allows us to estimate the expected proportion of secondary infections
attributable to the most infectious proportion of the population. For example,
we estimate that the 20% most infectious index cases account for approximately
75-98% of the expected secondary infections with 95% posterior probability. In
addition, we highlight that heterogeneity is a vital consideration when
estimating $R_t$.


  This paper uses Covasim, an agent-based model (ABM) of COVID-19, to evaluate
and scenarios of epidemic spread in New York State (USA), the UK, and the
Novosibirsk region (Russia). Epidemiological parameters such as contagiousness
(virus transmission rate), initial number of infected people, and probability
of being tested depend on the region's demographic and geographical features,
the containment measures introduced; they are calibrated to data about COVID-19
spread in the region of interest. At the first stage of our study,
epidemiological data (numbers of people tested, diagnoses, critical cases,
hospitalizations, and deaths) for each of the mentioned regions were analyzed.
The data were characterized in terms of seasonality, stationarity, and
dependency spaces, and were extrapolated using machine learning techniques to
specify unknown epidemiological parameters of the model. At the second stage,
the Optuna optimizer based on the tree Parzen estimation method for objective
function minimization was applied to determine the model's unknown parameters.
The model was validated with the historical data of 2020. The modeled results
of COVID-19 spread in New York State, the UK and the Novosibirsk region have
demonstrated that if the level of testing and containment measures is
preserved, the number of positive cases in New York State and the Novosibirsk
region will remain the same during March of 2021, while in the UK it will
reduce. Due to the features of the data for the Novosibirsk region (two
datasets are characterized as stationary series with probability of 1), the
forecast precision is relatively high for the number of hospitalizations, but
is lower for new cases of COVID-19.


  To better understand the potential of drug repurposing in COVID-19, we
analyzed control strategies over essential host factors for SARS-CoV-2
infection. We constructed comprehensive directed protein-protein interaction
networks integrating the top ranked host factors, drug target proteins, and
directed protein-protein interaction data. We analyzed the networks to identify
drug targets and combinations thereof that offer efficient control over the
host factors. We validated our findings against clinical studies data and
bioinformatics studies. Our method offers a new insight into the molecular
details of the disease and into potentially new therapy targets for it. Our
approach for drug repurposing is significant beyond COVID-19 and may be applied
also to other diseases.


  Malaria is an infectious disease with an immense global health burden.
Plasmodium vivax is the most geographically widespread species of malaria.
Relapsing infections, caused by the activation of liver-stage parasites known
as hypnozoites, are a critical feature of the epidemiology of Plasmodium vivax.
Hypnozoites remain dormant in the liver for weeks or months after inoculation,
but cause relapsing infections upon activation. Here, we introduce a dynamic
probability model of the activation-clearance process governing both potential
relapses and the size of the hypnozoite reservoir. We begin by modelling
activation-clearance dynamics for a single hypnozoite using a continuous-time
Markov chain. We then extend our analysis to consider activation-clearance
dynamics for a single mosquito bite, which can simultaneously establish
multiple hypnozoites, under the assumption of independent hypnozoite behaviour.
We derive analytic expressions for the time to first relapse and the time to
hypnozoite clearance for mosquito bites establishing variable numbers of
hypnozoites, both of which are quantities of epidemiological significance. Our
results extend those in the literature, which were limited due to an assumption
of non-independence. Our within-host model can be embedded readily in
multi-scale models and epidemiological frameworks, with analytic solutions
increasing the tractability of statistical inference and analysis. Our work
therefore provides a foundation for further work on immune development and
epidemiological-scale analysis, both of which are important for achieving the
goal of malaria elimination.


  Designing novel protein sequences for a desired 3D topological fold is a
fundamental yet non-trivial task in protein engineering. Challenges exist due
to the complex sequence--fold relationship, as well as the difficulties to
capture the diversity of the sequences (therefore structures and functions)
within a fold. To overcome these challenges, we propose Fold2Seq, a novel
transformer-based generative framework for designing protein sequences
conditioned on a specific target fold. To model the complex sequence--structure
relationship, Fold2Seq jointly learns a sequence embedding using a transformer
and a fold embedding from the density of secondary structural elements in 3D
voxels. On test sets with single, high-resolution and complete structure inputs
for individual folds, our experiments demonstrate improved or comparable
performance of Fold2Seq in terms of speed, coverage, and reliability for
sequence design, when compared to existing state-of-the-art methods that
include data-driven deep generative models and physics-based RosettaDesign. The
unique advantages of fold-based Fold2Seq, in comparison to a structure-based
deep model and RosettaDesign, become more evident on three additional
real-world challenges originating from low-quality, incomplete, or ambiguous
input structures. Source code and data are available at
https://github.com/IBM/fold2seq.


  A precise estimate of allele and haplotype polymorphism is of great interest
in theoretical population genetics, but also has practical applications, such
as bone marrow registries management. Allele polymorphism is driven mainly by
point mutations, while haplotype polymorphism is also affected by
recombination. Current estimates treat recombination as mutations in an
infinite site model. We here show that even in the simple case of two loci in a
haploid individual, for a finite population, most recombination events produce
existing haplotypes, and as such are silent. Silent recombination considerably
reduces the total number of haplotypes expected from the infinite site model
for populations that are not much larger than one over the mutation rate.
Moreover, in contrast with mutations, the number of haplotypes does not grow
linearly with the population size. We hence propose a more accurate estimate of
the total number of haplotypes that takes into account silent recombination. We
study large-scale Human Leukocyte Antigen (HLA) haplotype frequencies from
human populations to show that the current estimated recombination rate in the
HLA region is underestimated.


  Liquid-liquid phase separation is emerging as a crucial phenomenon in several
fundamental cell processes. A range of eukaryotic systems exhibit liquid
condensates. However, their function in bacteria, which in general lack
membrane-bound compartments, remains less clear. Here, we used high-resolution
optical microscopy to observe single bacterial aggresomes, nanostructured
intracellular assemblies of proteins, to undercover their role in cell stress.
We find that proteins inside aggresomes are mobile and undergo dynamic
turnover, consistent with a liquid state. Our observations are in quantitative
agreement with phase-separated liquid droplet formation driven by interacting
proteins under thermal equilibrium that nucleate following diffusive collisions
in the cytoplasm. We have discovered aggresomes in multiple species of
bacteria, and show that these emergent, metastable liquid-structured protein
assemblies increase bacterial fitness by enabling cells to tolerate
environmental stresses.


  To forecast the time dynamics of an epidemic, we propose a discrete
stochastic model that unifies and generalizes previous approaches to the
subject. Viewing a given population of individuals or groups of individuals
with given health state attributes as living in and moving between the nodes of
a graph, we use Monte-Carlo Markov Chain techniques to simulate the movements
and health state changes of the individuals according to given probabilities of
stay that have been preassigned to each of the nodes. We utilize this model to
either capture and predict the future geographic evolution of an epidemic in
time, or the evolution of an epidemic inside a heterogeneous population which
is divided into homogeneous sub-populations, or, more generally, its evolution
in a combination or superposition of the previous two contexts. We also prove
that when the size of the population increases and a natural hypothesis is
satisfied, the stochastic process associated to our model converges to a
deterministic process. Indeed, when the length of the time step used in the
discrete model converges to zero, in the limit this deterministic process is
driven by a differential equation yielding the evolution of the expectation
value of the number of infected as a function of time. In the second part of
the paper, we apply our model to study the evolution of the Covid-19 epidemic.
We deduce a decomposition of the function yielding the number of infectious
individuals into "wavelets", which allows to trace in time the expectation
value for the number of infections inside each sub-population. Within this
framework, we also discuss possible causes for the occurrence of multiple
epidemiological waves.


  Dormancy is a costly adaptive strategy that is widespread among living
organisms inhabiting diverse environments. We explore mathematical models of
predator-prey systems, in order to assess the impact of prey dormancy on the
competition between two types of prey, a perennially active (PA) and capable of
entering dormancy (dormitive). Both the active form and the dormant form of the
dormitive prey are individually at a disadvantage compared to the PA prey and
would go extinct due to their low growth rate, energy waste on the production
of dormant prey, and inability of the latter to grow autonomously. However, the
dormitive prey can paradoxically outcompete the PA prey with superior traits
and even cause its extinction by alternating between the two losing strategies.
This outcome recapitulates the game-theoretic Parrondo's paradox, where two
losing strategies combine to achieve a winning outcome. We observed higher
fitness of the dormitive prey in rich environments because a large predator
population in a rich environment cannot be supported by the prey without
adopting an evasive strategy, that is, dormancy. In such environments,
populations experience large-scale fluctuations, which can be survived by
dormitive but not by PA prey. Dormancy of the prey appears to be a natural
evolutionary response to self-destructive over-predation that stabilizes
evolving predator-prey systems through Parrondo's paradox.


  Background: With the rapid growth of massively parallel sequencing
technologies, still more laboratories are utilizing sequenced DNA fragments for
genomic analyses. Interpretation of sequencing data is, however, strongly
dependent on bioinformatics processing, which is often too demanding for
clinicians and researchers without a computational background. Another problem
represents the reproducibility of computational analyses across separated
computational centers with inconsistent versions of installed libraries and
bioinformatics tools.
  Results: We propose an easily extensible set of computational pipelines,
called SnakeLines, for processing sequencing reads; including mapping,
assembly, variant calling, viral identification, transcriptomics, metagenomics,
and methylation analysis. Individual steps of an analysis, along with methods
and their parameters can be readily modified in a single configuration file.
Provided pipelines are embedded in virtual environments that ensure isolation
of required resources from the host operating system, rapid deployment, and
reproducibility of analysis across different Unix-based platforms.
  Conclusion: SnakeLines is a powerful framework for the automation of
bioinformatics analyses, with emphasis on a simple set-up, modifications,
extensibility, and reproducibility.
  Keywords: Computational pipeline, framework, massively parallel sequencing,
reproducibility, virtual environment


  Aldosterone, the main physiological mineralocorticoid in humans and other
terrestrial vertebrates, first appears in lungfish, which are lobe-finned fish
that are forerunners of terrestrial vertebrates. Aldosterone activation of the
MR regulates internal homeostasis of water, sodium and potassium, which was
critical in the conquest of land by vertebrates. We studied transcriptional
activation of the slender African lungfish MR by aldosterone, other
corticosteroids and progesterone and find that aldosterone,
11-deoxycorticosterone, 11-deoxycortisol and progesterone have half-maximal
responses (EC50s) below 1 nM and are potential physiological
mineralocorticoids. In contrast, EC50s for corticosterone and cortisol were 23
nM and 66 nM, respectively. Unexpectedly, truncated lungfish MR, consisting of
the DNA-binding, hinge and steroid-binding domains, had a stronger response to
corticosteroids and progesterone than full-length lungfish MR, indicating that
the N-terminal domain represses steroid activation of lungfish MR, unlike human
MR in which the N-terminal domain contains an activation function. BLAST
searches of GenBank did not retrieve a GR ortholog, leading us to test
dexamethasone and triamcinolone for activation of lungfish MR. At 10 nM, both
synthetic glucocorticoids are about 4-fold stronger than 10 nM aldosterone in
activating full-length lungfish MR, leading us to propose that lungfish MR also
functions as a GR.


  Cryo-electron tomography (Cryo-ET) is a 3D imaging technique that enables the
systemic study of shape, abundance, and distribution of macromolecular
structures in single cells in near-atomic resolution. However, the systematic
and efficient $\textit{de novo}$ recognition and recovery of macromolecular
structures captured by Cryo-ET are very challenging due to the structural
complexity and imaging limits. Even macromolecules with identical structures
have various appearances due to different orientations and imaging limits, such
as noise and the missing wedge effect. Explicitly disentangling the semantic
features of macromolecules is crucial for performing several downstream
analyses on the macromolecules. This paper has addressed the problem by
proposing a 3D Spatial Variational Autoencoder that explicitly disentangle the
structure, orientation, and shift of macromolecules. Extensive experiments on
both synthesized and real cryo-ET datasets and cross-domain evaluations
demonstrate the efficacy of our method.


  We consider an exponentially growing population of cells undergoing mutations
and ask about the effect of reproductive fluctuations (genetic drift) on its
long-term evolution. We combine first step analysis with the stochastic
dynamics of a birth-death process to analytically calculate the probability
that the parent of a given genotype will go extinct. We compare the results
with numerical simulations and show how this turnover of genetic clones can be
used to infer the rates underlying the population dynamics. Our work is
motivated by growing populations of tumour cells, the epidemic spread of
viruses, and bacterial growth.


  Mathematical models are formal and simplified representations of the
knowledge related to a phenomenon. In classical epidemic models, a neglected
aspect is the heterogeneity of disease transmission and progression linked to
the viral load of each infectious individual. Here, we attempt to investigate
the interplay between the evolution of individuals' viral load and the epidemic
dynamics from a theoretical point of view. In the framework of multi-agent
systems, we propose a particle stochastic model describing the infection
transmission through interactions among agents and the individual physiological
course of the disease. Agents have a double microscopic state: a discrete
label, that denotes the epidemiological compartment to which they belong and
switches in consequence of a Markovian process, and a microscopic trait,
representing a normalized measure of their viral load, that changes in
consequence of binary interactions or interactions with a background.
Specifically, we consider Susceptible--Infected--Removed--like dynamics where
infectious individuals may be isolated from the general population and the
isolation rate may depend on the viral load sensitivity and frequency of tests.
We derive kinetic evolution equations for the distribution functions of the
viral load of the individuals in each compartment, whence, via suitable
upscaling procedures, we obtain a macroscopic model for the densities and viral
load momentum. We perform then a qualitative analysis of the ensuing
macroscopic model, and we present numerical tests in the case of both constant
and viral load-dependent isolation control. Also, the matching between the
aggregate trends obtained from the macroscopic descriptions and the original
particle dynamics simulated by a Monte Carlo approach is investigated.


  In this paper we consider a SEIRD epidemic model for a population composed by
two groups with asymmetric interaction. Given two sets of parameters for the
model, we estimate the distance of the solutions for the second group based on
the distance of solutions for the first one. We also study the final size of
the epidemic for each group. We illustrate our results with the spread of
coronavirus disease 2019 (COVID-19) pandemic in the New York County (USA) for
the initial stage of the contamination, and in the the cities of Petrolina and
Juazeiro (Brazil).


  The dynamics of a mathematical model of the Calvin cycle, which is part of
photosynthesis, is analysed. Since diffusion of ATP is included in the model a
system of reaction-diffusion equations is obtained. It is proved that for a
suitable choice of parameters there exist spatially inhomogeneous positive
steady states, in fact infinitely many of them. It is also shown that all
positive steady states, homogeneous and inhomogeneous, are nonlinearly
unstable. The only smooth steady state which could be stable is a trivial one,
where all concentrations except that of ATP are zero. It is found that in the
spatially homogeneous case there are steady states with the property that the
linearization about that state has eigenvalues which are not real, indicating
the presence of oscillations. Numerical simulations exhibit solutions for which
the concentrations are not monotone functions of time.


  We propose a physical theory underlying the temporal evolution of competing
virus variants that relies on the existence of (quasi) fixed points capturing
the large time scale invariance of the dynamics. To motivate our result we
first modify the time-honoured compartmental models of the SIR type to account
for the existence of competing variants and then show how their evolution can
be naturally re-phrased in terms of flow equations ending at quasi fixed
points. As the natural next step we employ (near) scale invariance to organise
the time evolution of the competing variants within the effective description
of the epidemic Renormalization Group framework. We test the resulting theory
against the time evolution of COVID-19 virus variants that validate the theory
empirically.


  In temperate climates, infection rates of enveloped viruses peak during the
winter. While these seasonal trends are established in influenza and human
coronaviruses, the mechanisms driving the variation remain poorly understood
and thus difficult to extend to similar viruses like SARS-CoV-2. In this study,
we use the Quadrature-based model of Respiratory Aerosol and Droplets (QuaRAD)
to explore the sensitivity of airborne transmission to the seasonal variation
in indoor relative humidity across the wide range of relevant conditions, using
SARS-CoV-2 as an example. Relative humidity impacts the evaporation rate and
equilibrium size of airborne particles, which in turn may impact particle
removal rates and virion viability. Across a large ensemble of scenarios, we
found that the dry indoor conditions typical of the winter season lead to
slower inactivation than in the more humid summer season; in poorly ventilated
spaces, this reduction in inactivation rates increases the concentration of
active virions, but this effect was important when the susceptible person was
farther than 2 m downwind of the infectious person. On the other hand, changes
in particle settling velocity with relative humidity did not significantly
affect the removal or travel distance of virus-laden scenarios.


  MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play an
important role in post-transcriptional gene regulation. However, the
experimental determination of miRNA sequence and structure is both expensive
and time-consuming. Therefore, computational and machine learning-based
approaches have been adopted to predict novel microRNAs. With the involvement
of data science and machine learning in biology, multiple research studies have
been conducted to find microRNAs with different computational methods and
different miRNA features. Multiple approaches are discussed in detail
considering the learning algorithm/s used, features considered, dataset/s used
and the criteria used in evaluations. This systematic review focuses on the
machine learning methods developed for miRNA identification in plants. This
will help researchers to gain a detailed idea about past studies and identify
novel paths that solve drawbacks occurred in past studies. Our findings
highlight the need for plant-specific computational methods for miRNA
identification.


  Relative lifetimes of inherent double stranded DNA openings with lengths up
to ten base pairs are presented for different gene promoters and corresponding
mutants that either increase or decrease transcriptional activity, in the
framework of the Peyrard-Bishop-Dauxois model. Extensive microcanonical
simulations are used, with energies corresponding to physiological temperature.
The bubble lifetime profiles along the DNA sequences demonstrate a significant
reduction of the average lifetime at the mutation sites when the mutated
promoter decreases transcription, while a corresponding enhancement of the
bubble lifetime is observed in the case of mutations leading to increased
transcription. The relative difference of bubble lifetimes between the mutated
and the wild type promoters at the position of mutation varies from 20% to more
than 30% as the bubble length is decreasing


  The concept of k-spectrum for genomes is here investigated as a basic tool to
analyze genomes. Related spectral notions based on k-mers are introduced with
some related mathematical properties which are relevant for informational
analysis of genomes. Procedures to generate spectral segmentations of genomes
are provided and are tested (under several values of length k for k-mers) on
cases of real genomes, such as some human chromosomes and Saccharomyces
cerevisiae.


  Gene expression datasets consist of thousand of genes with relatively small
samplesizes (i.e. are large-$p$-small-$n$). Moreover, dependencies of various
orders co-exist in the datasets. In the Undirected probabilistic Graphical
Model (UGM) framework the Glasso algorithm has been proposed to deal with high
dimensional micro-array datasets forcing sparsity. Also, modifications of the
default Glasso algorithm are developed to overcome the problem of complex
interaction structure. In this work we advocate the use of a simple score-based
Hill Climbing algorithm (HC) that learns Gaussian Bayesian Networks (BNs)
leaning on Directed Acyclic Graphs (DAGs). We compare HC with Glasso and its
modifications in the UGM framework on their capability to reconstruct GRNs from
micro-array data belonging to the Escherichia Coli genome. We benefit from the
analytical properties of the Joint Probability Density (JPD) function on which
both directed and undirected PGMs build to convert DAGs to UGMs.
  We conclude that dependencies in complex data are learned best by the HC
algorithm, presenting them most accurately and efficiently, simultaneously
modelling strong local and weaker but significant global connections coexisting
in the gene expression dataset. The HC algorithm adapts intrinsically to the
complex dependency structure of the dataset, without forcing a specific
structure in advance. On the contrary, Glasso and modifications model
unnecessary dependencies at the expense of the probabilistic information in the
network and of a structural bias in the JPD function that can only be relieved
including many parameters.