  The expansion of deleted mitochondrial DNA (mtDNA) molecules has been linked
to ageing, particularly in skeletal muscle fibres; its mechanism has remained
unclear for three decades. Previous accounts assigned a replicative advantage
to the deletions, but there is evidence that cells can, instead, selectively
remove defective mtDNA. We present a spatial model that, without a replicative
advantage, but instead through a combination of enhanced density for mutants
and noise, produces a wave of expanding mutations with wave speed consistent
with experimental data, unlike a standard model based on replicative advantage.
We provide a formula that predicts that the wave speed drops with copy number,
in agreement with experimental data. Crucially, our model yields travelling
waves of mutants even if mutants are preferentially eliminated. Justified by
this exemplar of how noise, density and spatial structure affect muscle ageing,
we introduce the mechanism of stochastic survival of the densest, an
alternative to replicative advantage, that may underpin other phenomena, like
the evolution of altruism.


  The Susceptible-Infected-Recovered (SIR) epidemic model as well as its
generalizations are extensively used for the study of the spread of infectious
diseases, and for the understanding of the dynamical evolution of epidemics.
From SIR type models only the model without vital dynamics has an exact
analytic solution, which can be obtained in an exact parametric form. The SIR
model with vital dynamics, the simplest extension of the basic SIR model, does
not admit a closed form representation of the solution. However, in order to
perform the comparison with the epidemiological data accurate representations
of the time evolution of the SIR model with vital dynamics would be very
useful. In the present paper, we obtain first the basic evolution equation of
the SIR model with vital dynamics, which is given by a strongly nonlinear
second order differential equation. Then we obtain a series representation of
the solution of the model, by using the Adomian and Laplace-Adomian
Decomposition Methods to solve the dynamical evolution equation of the model.
The solutions are expressed in the form of infinite series. The series
representations of the time evolution of the SIR model with vital dynamics are
compared with the exact numerical solutions of the model, and we find that, at
least for a specific range of parameters, there is a good agreement between the
Adomian and Laplace-Adomian semianalytical solutions, containing only a small
number of terms, and the numerical results.


  The decision of whether or not to vaccinate is a complex one. It involves the
contribution both to a social good -- herd immunity -- and to one's own well
being. It is informed by social influence, personal experience, education, and
mass media. In our work, we investigate a situation in which individuals make
their choice based on how social neighbourhood responded to previous epidemics.
We do this by proposing a minimalistic model using components from game theory,
network theory and the modelling of epidemic spreading, and opinion dynamics.
Individuals can use the information about the neighbourhood in two ways --
either they follow the majority or the best-performing neighbour. Furthermore,
we let individuals learn which of these two decision-making strategies to
follow from their experience. Our results show that the flexibility of
individuals to chose how to integrate information from the neighbourhood
increases the vaccine uptake and decreases the epidemic severity if the
following conditions are fulfilled. First, the initial fraction of individuals
who imitate the neighbourhood majority should be limited, and second, the
memory of previous outbreaks should be sufficiently long. These results have
implications for the acceptance of novel vaccines and raising awareness about
vaccination, while also pointing to promising future research directions.


  In this work, we provide a general mathematical formalism to study the
optimal control of an epidemic, such as the COVID-19 pandemic, via incentives
to lockdown and testing. In particular, we model the interplay between the
government and the population as a principal-agent problem with moral hazard,
\`a la Cvitani\'c, Possama\"i, and Touzi [27], while an epidemic is spreading
according to dynamics given by compartmental stochastic SIS or SIR models, as
proposed respectively by Gray, Greenhalgh, Hu, Mao, and Pan [45] and Tornatore,
Buccellato, and Vetro [88]. More precisely, to limit the spread of a virus, the
population can decrease the transmission rate of the disease by reducing
interactions between individuals. However, this effort, which cannot be
perfectly monitored by the government, comes at social and monetary cost for
the population. To mitigate this cost, and thus encourage the lockdown of the
population, the government can put in place an incentive policy, in the form of
a tax or subsidy. In addition, the government may also implement a testing
policy in order to know more precisely the spread of the epidemic within the
country, and to isolate infected individuals. In terms of technical results, we
demonstrate the optimal form of the tax, indexed on the proportion of infected
individuals, as well as the optimal effort of the population, namely the
transmission rate chosen in response to this tax. The government's optimisation
problem then boils down to solving an Hamilton-Jacobi-Bellman equation.
Numerical results confirm that if a tax policy is implemented, the population
is encouraged to significantly reduce its interactions. If the government also
adjusts its testing policy, less effort is required on the population side,
individuals can interact almost as usual, and the epidemic is largely contained
by the targeted isolation of positively-tested individuals.


  Synthetic biology constructs often rely upon the introduction of "circuit"
genes into host cells, in order to express novel proteins and thus endow the
host with a desired behavior. The expression of these new genes "consumes"
existing resources in the cell, such as ATP, RNA polymerase, amino acids, and
ribosomes. Ribosomal competition among strands of mRNA may be described by a
system of nonlinear ODEs called the Ribosomal Flow Model (RFM). The competition
for resources between host and circuit genes can be ameliorated by splitting
the ribosome pool by use of orthogonal ribosomes, where the circuit genes are
exclusively translated by mutated ribosomes. In this work, the RFM system is
extended to include orthogonal ribosome competition. This Orthogonal Ribosomal
Flow Model (ORFM) is proven to be stable through the use of Robust Lyapunov
Functions. The optimization problem of maximizing the weighted protein
translation rate by adjusting allocation of ribosomal species is formulated and
implemented.


  The COVID-19 pandemic has forced policy makers to decree urgent confinements
to stop a rapid and massive contagion. However, after that stage, societies are
being forced to find an equilibrium between the need to reduce contagion rates
and the need to reopen their economies. The experience hitherto lived has
provided data on the evolution of the pandemic, in particular the population
dynamics as a result of the public health measures enacted. This allows the
formulation of forecasting mathematical models to anticipate the consequences
of political decisions. Here we propose a model to do so and apply it to the
case of Portugal. With a mathematical deterministic model, described by a
system of ordinary differential equations, we fit the real evolution of
COVID-19 in this country. After identification of the population readiness to
follow social restrictions, by analyzing the social media, we incorporate this
effect in a version of the model that allow us to check different scenarios.
This is realized by considering a Monte Carlo discrete version of the previous
model coupled via a complex network. Then, we apply optimal control theory to
maximize the number of people returning to "normal life" and minimizing the
number of active infected individuals with minimal economical costs while
warranting a low level of hospitalizations. This work allows testing various
scenarios of pandemic management (closure of sectors of the economy,
partial/total compliance with protection measures by citizens, number of beds
in intensive care units, etc.), ensuring the responsiveness of the health
system, thus being a public health decision support tool.


  RNA 3D motifs are recurrent substructures, modelled as networks of base pair
interactions, which are crucial for understanding structure-function
relationships. The task of automatically identifying such motifs is
computationally hard, and remains a key challenge in the field of RNA
structural biology and network analysis. State of the art methods solve special
cases of the motif problem by constraining the structural variability in
occurrences of a motif, and narrowing the substructure search space. Here, we
relax these constraints by posing the motif finding problem as a graph
representation learning and clustering task. This framing takes advantage of
the continuous nature of graph representations to model the flexibility and
variability of RNA motifs in an efficient manner. We propose a set of node
similarity functions, clustering methods, and motif construction algorithms to
recover flexible RNA motifs. Our tool, VeRNAl can be easily customized by users
to desired levels of motif flexibility, abundance and size. We show that VeRNAl
is able to retrieve and expand known classes of motifs, as well as to propose
novel motifs.


  There remains an urgent need to identify existing drugs that might be
suitable for treating patients suffering from COVID-19 infection. Drugs rarely
act at a single molecular target, with off target effects often being
responsible for undesirable side effects and sometimes, beneficial synergy
between targets for a specific illness. Off target activities have also led to
blockbuster drugs in some cases, e.g. Viagra for erectile dysfunction and
Minoxidil for male pattern hair loss. Drugs already in use or in clinical
trials plus approved natural products constitute a rich resource for discovery
of therapeutic agents that can be repurposed for existing and new conditions,
based on the rationale that they have already been assessed for safety in man.
A key question then is how to rapidly and efficiently screen such compounds for
activity against new pandemic pathogens such as COVID-19. Here we show how a
fast and robust computational process can be used to screen large libraries of
drugs and natural compounds to identify those that may inhibit the main
protease of SARS-Cov-2 (3CL pro, Mpro). We show how the resulting shortlist of
candidates with strongest binding affinities is highly enriched in compounds
that have been independently identified as potential antivirals against
COVID-19. The top candidates also include a substantial number of drugs and
natural products not previously identified as having potential COVID-19
activity, thereby providing additional targets for experimental validation.
This in silico screening pipeline may also be useful for repurposing of
existing drugs and discovery of new drug candidates against other medically
important pathogens and for use in future pandemics.


  Social hierarchy is important that can not be ignored in human socioeconomic
activities and in the animal world. Here we incorporate this factor into the
evolutionary game to see what impact it could have on the cooperation outcome.
The probabilistic strategy adoption between two players is then not only
determined by their payoffs, but also by their hierarchy difference -- players
in the high rank are more likely to reproduce their strategies than the peers
in the low rank. Through simulating the evolution of Prisoners' dilemma game
with three hierarchical distributions, we find that the levels of cooperation
are enhanced in all cases, and the enhancement is optimal in the uniform case.
The enhancement is due to the fact that the presence of hierarchy facilitates
the formation of cooperation clusters with high-rank players acting as the
nucleation cores. This mechanism remains valid on Barab\'asi-Albert scale-free
networks, in particular the cooperation enhancement is maximal when the hubs
are of higher social ranks. We also study a two-hierarchy model, where similar
cooperation promotion is revealed and some theoretical analyses are provided.
Our finding may partially explain why the social hierarchy is so ubiquitous on
this planet.


  In this work, we present a paradigm bridging electromagnetic (EM) and
molecular communication through a stimuli-responsive intra-body model. It has
been established that protein molecules, which play a key role in governing
cell behavior, can be selectively stimulated using Terahertz (THz) band
frequencies. By triggering protein vibrational modes using THz waves, we induce
changes in protein conformation, resulting in the activation of a controlled
cascade of biochemical and biomechanical events. To analyze such an
interaction, we formulate a communication system composed of a nanoantenna
transmitter and a protein receiver. We adopt a Markov chain model to account
for protein stochasticity with transition rates governed by the nanoantenna
force. Both two-state and multi-state protein models are presented to depict
different biological configurations. Closed form expressions for the mutual
information of each scenario is derived and maximized to find the capacity
between the input nanoantenna force and the protein state. The results we
obtain indicate that controlled protein signaling provides a communication
platform for information transmission between the nanoantenna and the protein
with a clear physical significance. The analysis reported in this work should
further research into the EM-based control of protein networks.


  With the emergence of a new pandemic worldwide, a novel strategy to approach
it has emerged. Several initiatives under the umbrella of "open science" are
contributing to tackle this unprecedented situation. In particular, the "R
Language and Environment for Statistical Computing" offers an excellent tool
and ecosystem for approaches focusing on open science and reproducible results.
Hence it is not surprising that with the onset of the pandemic, a large number
of R packages and resources were made available for researches working in the
pandemic. In this paper, we present an R package that allows users to access
and analyze worldwide data from resources publicly available. We will introduce
the covid19.analytics package, focusing in its capabilities and presenting a
particular study case where we describe how to deploy the "COVID19.ANALYTICS
Dashboard Explorer".


  The main reasons for the ongoing COVID-19 (coronavirus disease 2019) pandemic
are the unavailability of recommended efficacious drugs or vaccines along with
the human to human transmission nature of SARS-CoV-2 virus. So, there is urgent
need to search appropriate therapeutic approach by repurposing approved drugs.
In this communication, molecular docking analyses of two influenza antiviral
drugs baloxavir acid (BXA) and baloxavir marboxil (BXM) were performed with the
three therapeutic target proteins of severe acute respiratory syndrome
coronavirus 2 (SARS-CoV-2), i.e., main protease (Mpro), papain-like protease
(PLpro) and RNA-dependent RNA polymerase (RdRp). The molecular docking results
of both the drugs BXA and BXM were analysed and compared. The investigational
drug BXA binds at the active site of Mpro and RdRp, whereas the approved drug
BXM binds only at the active site of RdRp. Also, comparison of dock score
revealed that BXA is binding more effectively at the active site of RdRp than
BXM. The computational molecular docking revealed that the drug BXA may be more
effective against COVID-19 as compared to BXM.


  We introduce a mathematical description of the impact of sociality in the
spread of infectious diseases by integrating an epidemiological dynamics with a
kinetic modeling of population-based contacts. The kinetic description leads to
study the evolution over time of Boltzmann-type equations describing the number
densities of social contacts of susceptible, infected and recovered
individuals, whose proportions are driven by a classical SIR-type compartmental
model in epidemiology. Explicit calculations show that the spread of the
disease is closely related to moments of the contact distribution. Furthermore,
the kinetic model allows to clarify how a selective control can be assumed to
achieve a minimal lockdown strategy by only reducing individuals undergoing a
very large number of daily contacts. We conduct numerical simulations which
confirm the ability of the model to describe different phenomena characteristic
of the rapid spread of an epidemic. Motivated by the COVID-19 pandemic, a last
part is dedicated to fit numerical solutions of the proposed model with
infection data coming from different European countries.


  Gene regulatory networks (GRNs) play a central role in cellular
decision-making. Understanding their structure and how it impacts their
dynamics constitutes thus a fundamental biological question. GRNs are
frequently modeled as Boolean networks, which are intuitive, simple to
describe, and can yield qualitative results even when data is sparse. We
assembled the largest repository of expert-curated Boolean GRN models. A
meta-analysis of this diverse set of models reveals several design principles.
GRNs exhibit more canalization, redundancy and stable dynamics than expected.
Moreover, they are enriched for certain recurring network motifs. This raises
the important question why evolution favors these design mechanisms.


  Mathematical models are increasing adopted for setting targets for disease
prevention and control. As model-informed policies are implemented, however,
the inaccuracies of some forecasts become apparent, for example overprediction
of infection burdens and overestimation of intervention impacts. Here, we
attribute these discrepancies to methodological limitations in capturing the
heterogeneities of real-world systems. The mechanisms underpinning single
factors for infection and their interactions determine individual propensities
to acquire disease. These are potentially so numerous that to attain a full
mechanistic description may be unfeasible. To contribute constructively to the
development of health policies, model developers either leave factors out
(reductionism) or adopt a broader but coarse description (holism). In our view,
predictive capacity requires holistic descriptions of heterogeneity which are
currently underutilised in infectious disease epidemiology but common in other
disciplines.


  Learning on 3D structures of large biomolecules is emerging as a distinct
area in machine learning, but there has yet to emerge a unifying network
architecture that simultaneously leverages the graph-structured and geometric
aspects of the problem domain. To address this gap, we introduce geometric
vector perceptrons, which extend standard dense layers to operate on
collections of Euclidean vectors. Graph neural networks equipped with such
layers are able to perform both geometric and relational reasoning on efficient
and natural representations of macromolecular structure. We demonstrate our
approach on two important problems in learning from protein structure: model
quality assessment and computational protein design. Our approach improves over
existing classes of architectures, including state-of-the-art graph-based and
voxel-based methods. We release our code at https://github.com/drorlab/gvp.


  During a pandemic, there are conflicting demands arising from public health
and economic cost. Lockdowns are a common way of containing infections, but
they adversely affect the economy. We study the question of how to minimise the
economic damage of a lockdown while still containing infections. Our analysis
is based on the SIR model, which we analyse using a clock set by the virus.
This use of the "virus time" permits a clean mathematical formulation of our
problem. We optimise the economic cost for a fixed health cost and arrive at a
strategy for navigating the pandemic. This involves adjusting the level of
lockdowns in a controlled manner so as to minimise the economic cost.


  Cyclic dominance offers not just a way to maintain biodiversity, but also
serves as a sort of defensive alliance against an external invader.
Interestingly, a new level of competition can be observed when two cyclic loops
are present. Here the inner invasion speed plays a decisive role on the
evolutionary outcome, because faster invasion rate provides an evolutionary
advantage to an alliance. In this Letter we demonstrate that the heterogeneity
of inner invasion rates makes an alliance vulnerable against a loop where group
members are equal. Quite surprisingly, a loop where invasion rates are uniform
can still dominate an alliance formed by heteregeneous rates even if the
average speed of invasion is significantly higher in the latter group. At a
specific range of parameter space, when intergroup invasion or the average
inner invasion is moderate, the heterogeneous alliance with higher internal
invasion speed may prevail, or the system terminates onto a novel 4- or
5-species solution.


  A hallmark of topological phases is the occurrence of topologically protected
modes at the system`s boundary. Here we find topological phases in the
antisymmetric Lotka-Volterra equation (ALVE). The ALVE is a nonlinear dynamical
system and describes, e.g., the evolutionary dynamics of a rock-paper-scissors
cycle. On a one-dimensional chain of rock-paper-scissor cycles, topological
phases become manifest as robust polarization states. At the transition point
between left and right polarization, solitonic waves are observed. This
topological phase transition lies in symmetry class $D$ within the "ten-fold
way" classification as also realized by 1D topological superconductors.


  Epidemic spread in a population is traditionally modeled via
compartmentalized models which represent the free evolution of disease in
absence of any intervention policies. In addition, these models assume full
observability of disease cases and do not account for under-reporting. We
present a mathematical model, namely PolSIRD, which accounts for the
under-reporting by introducing an observation mechanism. It also captures the
effects of intervention policies on the disease spread parameters by leveraging
intervention policy data along with the reported disease cases. Furthermore, we
allow our recurrent model to learn the initial hidden state of all compartments
end-to-end along with other parameters via gradient-based training. We apply
our model to the spread of the recent global outbreak of COVID-19 in the United
States, where our model outperforms the methods employed by the CDC in
predicting the spread. We also provide counterfactual simulations from our
model to analyze the effect of lifting the intervention policies prematurely
and our model correctly predicts the second wave of the epidemic.


  We develop a trait-based model founded on the hypothesis that biological
systems evolve and organize to maximize entropy production by dissipating
chemical and electromagnetic potentials over longer time scales than abiotic
processes by implementing temporal strategies. A marine food web consisting of
phytoplankton, bacteria and consumer functional groups is used to explore how
temporal strategies, or the lack there of, change entropy production in a
shallow pond that receives a continuous flow of reduced organic carbon plus
inorganic nitrogen and illumination from solar radiation with diel and seasonal
dynamics. Results show that a temporal strategy that employs an explicit
circadian clock produces more entropy than a passive strategy that uses
internal carbon storage or a balanced growth strategy that requires
phytoplankton to grow with fixed stoichiometry. When the community is forced to
operate at high specific growth rates near 2 d-1, the optimization-guided model
selects for phytoplankton ecotypes that exhibit complementary for winter versus
summer environmental conditions to increase entropy production. We also present
a new type of trait-based modeling where trait values are determined by
maximizing entropy production rather than by random selection.


  We model and calculate the fraction of infected population necessary to reach
herd immunity, taking into account the heterogeneity in infectiousness and
susceptibility, as well as the correlation between those two parameters. We
show that these cause the effective reproduction number to decrease more
rapidly, and consequently have a drastic effect on the estimate of the
necessary percentage of the population that has to contract the disease for
herd immunity to be reached. We quantify the difference between the size of the
infected population when the effective reproduction number decreases below 1
vs. the ultimate fraction of population that had contracted the disease. This
sheds light on an important distinction between herd immunity and the end of
the disease and highlights the importance of limiting the spread of the disease
even if we plan to naturally reach herd immunity. We analyze the effect of
various lock-down scenarios on the resulting final fraction of infected
population. We discuss implications to COVID-19 and other pandemics and compare
our theoretical results to population-based simulations. We consider the
dependence of the disease spread on the architecture of the infectiousness
graph and analyze different graph architectures and the limitations of the
graph models.


  In this technical report, we describe a new version of SimpleSBML which
provides an easier to use interface to python-libSBML allowing users of Python
to more easily construct, edit, and inspect SBML based models. The most
commonly used package for constructing SBML models in Python is python-libSBML
based on the C/C++ library libSBML. python-libSBML is a comprehensive library
with a large range of options but can be difficult for new users to learn and
requires long scripts to create even the simplest models. Inspecting existing
SBML models can also be difficult due to the complexity of the underlying
object model. Instead, we present SimpleSBML, a package that allows users to
add and inspect species, parameters, reactions, events, and rules to a libSBML
model with only one command for each. Models can be exported to SBML format,
and SBML files can be imported and converted to SimpleSBML commands making it
very easy to edit the original SBML model. In the new version, a range of `get'
methods is provided that allows users to inspect existing SBML models without
having to understand the underlying object model used by libSBML.


  According to the criticality hypothesis, collective biological systems should
operate in a special parameter region, close to so-called critical points,
where the collective behavior undergoes a qualitative change between different
dynamical regimes. Critical systems exhibit unique properties, which may
benefit collective information processing such as maximal responsiveness to
external stimuli. Besides neuronal and gene-regulatory networks, recent
empirical data suggests that also animal collectives may be examples of
self-organized critical systems. However, open questions about
self-organization mechanisms in animal groups remain: Evolutionary adaptation
towards a group-level optimum (group-level selection), implicitly assumed in
the "criticality hypothesis", appears in general not reasonable for
fission-fusion groups composed of non-related individuals. Furthermore,
previous theoretical work relies on non-spatial models, which ignore
potentially important self-organization and spatial sorting effects. Using a
generic, spatially-explicit model of schooling prey being attacked by a
predator, we show first that schools operating at criticality perform best.
However, this is not due to optimal response of the prey to the predator, as
suggested by the "criticality hypothesis", but rather due to the spatial
structure of the prey school at criticality. Secondly, by investigating
individual-level evolution, we show that strong spatial self-sorting effects at
the critical point lead to strong selection gradients, and make it an
evolutionary unstable state. Our results demonstrate the decisive role of
spatio-temporal phenomena in collective behavior, and that individual-level
selection is in general not a viable mechanism for self-tuning of unrelated
animal groups towards criticality.


  Protein folding is a phenomenon that has been studied for about 50 years and
still remains as an unsolved problem. The main feature of this process is that
it occurs as an all or none process, so a protein, can jump directly between
folded and unfolded states. It is proposed in this manuscript that the
cooperative phenomena associated to protein folding has its origin in the
synchronization oscillating phases of the instantaneously induced dipoles that
gives rise to the van der Waals London dispersion interaction. When the
oscillation of the induced dipoles is synchronized an enhanced interaction
regime is triggered and the protein enters a into a folding process. The
propagation of this regime throughout the molecular structure complete the
folding process. If a desynchronizing strong enough perturbation is introduced
at any component cooperative network the system enters in a weakened -- or even
a repulsive -- van der Waals dispersion interaction regime that propagate
itself throughout the structure, that triggers a phase transition that finally
unfolds the protein.


  The special epistemic characteristics of the COVID-19, such as the long
incubation period and the infection through asymptomatic cases, put severe
challenge to the containment of its outbreak. By the end of March 2020, China
has successfully controlled the within-spreading of COVID-19 at a high cost of
locking down most of its major cities, including the epicenter, Wuhan. Since
the low accuracy of outbreak data before the mid of Feb. 2020 forms a major
technical concern on those studies based on statistic inference from the early
outbreak. We apply the supervised learning techniques to identify and train
NP-Net-SIR model which turns out robust under poor data quality condition. By
the trained model parameters, we analyze the connection between population flow
and the cross-regional infection connection strength, based on which a set of
counterfactual analysis is carried out to study the necessity of lock-down and
substitutability between lock-down and the other containment measures. Our
findings support the existence of non-lock-down-typed measures that can reach
the same containment consequence as the lock-down, and provide useful guideline
for the design of a more flexible containment strategy.


  Interactions between biomolecules, electrons and protons are essential to
many fundamental processes sustaining life. It is therefore of interest to
build mathematical models of these bioelectrical processes not only to enhance
understanding but also to enable computer models to complement in vitro and in
vivo experiments.Such models can never be entirely accurate; it is nevertheless
important that the models are compatible with physical principles. Network
Thermodynamics, as implemented with bond graphs, provide one approach to
creating physically compatible mathematical models of bioelectrical systems.
This is illustrated using simple models of ion channels, redox reactions,
proton pumps and electrogenic membrane transporters thus demonstrating that the
approach can be used to build mathematical and computer models of a wide range
of bioelectrical systems.


  The molecular simulations solve the equation of motion of molecular systems,
making 3D shapes of molecules four-dimensional by adding the time coordinate.
These methods have a great potential in drug discovery because they can
realistically model the structures of protein molecules targeted by drugs as
well as the process of binding of a potential drug to its molecular target.
However, routine application of biomolecular simulations is hampered by the
very high computational costs of this method. Several methods have been
developed to address this problem. One of them, metadynamics, disfavors states
of the simulated system that have been already visited and thus forces the
system to explore new and new states. Here we present the package metadynminer
and metadynminer3d to analyze and visualize results from metadynamics, in
particular those produced by a popular metadynamics package Plumed.


  Typically, mathematical simulation studies on COVID-19 pandemic forecasting
are based on deterministic differential equations which assume that both the
number ($n$) of individuals in various epidemiological classes and the time
($t$) on which they depend are quantities that vary continuous. This picture
contrasts with the discrete representation of $n$ and $t$ underlying the real
epidemiological data reported in terms daily numbers of infection cases, for
which a description based on finite difference equations would be more
adequate. Adopting a logistic growth framework, in this paper we present a
quantitative analysis of the errors introduced by the continuous time
description. This analysis reveals that, although the height of the
epidemiological curve maximum is essentially unaffected, the position
$T_{1/2}^{c}$ obtained within the continuous time representation is
systematically shifted backwards in time with respect to the position
$T_{1/2}^{d}$ predicted within the discrete time representation. Rather
counterintuitively, the magnitude of this temporal shift $\tau \equiv
T_{1/2}^{c} - T_{1/2}^{d} < 0$ is basically insensitive to changes in infection
rate $\kappa$. For a broad range of $\kappa$ values deduced from COVID-19 data
at extreme situations (exponential growth in time and complete lockdown), we
found a rather robust estimate $\tau \simeq -2.65\,\mbox{day}^{-1}$. Being
obtained without any particular assumption, the present mathematical results
apply to logistic growth in general without any limitation to a specific real
system.


  In June 2020, Arizona, U.S., emerged as one of the world's worst coronavirus
disease 2019(COVID-19) spots after the stay-at-home order was lifted in the
middle of May. However, with the decisions to reimpose restrictions, the number
of COVID-19 cases has been declining, and Arizona is considered to be a good
model in slowing the epidemic. In this paper, we aimed to examine the COVID-19
situation in Arizona and assess the impact of human mobility change. We
constructed the mobility integrated metapopulation
susceptible-infectious-removed model and fitted to publicly available datasets
on COVID-19 cases and mobility changes in Arizona. Our simulations showed that
by reducing human mobility, the peak time was delayed, and the final size of
the epidemic was decreased in all three regions. Our analysis suggests that
rapid and effective decision making is crucial to control human mobility and,
therefore, COVID-19 epidemics. Until a vaccine is available, reimplementations
of mobility restrictions in response to the increase of new COVID-19 cases
might need to be considered in Arizona and beyond.


  Gene expression (GE) is an inherently random or stochastic or noisy process.
The randomness in different steps of GE, e.g., transcription, translation,
degradation, etc., leading to cell-to-cell variations in mRNA and protein
levels. This variation appears in organisms ranging from microbes to metazoans.
Stochastic gene expression has important consequences for cellular function.
The random fluctuations in protein levels produce variability in cellular
behavior. It is beneficial in some contexts and harmful to others. These
situations include stress response, metabolism, development, cell cycle,
circadian rhythms, and aging. Different model studies e.g., constitutive,
two-state, etc., reveal that the fluctuations in mRNA and protein levels arise
from different steps of gene expression among which the steps in transcription
have the maximum effect. The pulsatile mRNA production through RNAP-II based
reinitiation of transcription is an important part of gene transcription.
Though, the effect of that process on mRNA and protein levels is very little
known. The addition of any biochemical step in the constitutive or two-state
process generally decreases the mean and increases the Fano factor. In this
study, we have shown that the RNAP-II based reinitiation process in gene
transcription can have different effects on both mean and Fano factor at mRNA
levels in different model systems. It decreases the mean and Fano factor both
at the mRNA levels in the constitutive network whereas in other networks it can
simultaneously increase or decrease both quantities or it can have mixed-effect
at mRNA levels. We propose that a constitutive network with reinitiation
behaves like a product independent negative feedback circuit whereas other
networks behave as either product independent positive or negative or mixed
feedback circuit.


  In the framework of homogeneous susceptible-infected-recovered (SIR) models,
we use a control theory approach to identify optimal pandemic mitigation
strategies. We derive rather general conditions for reaching herd immunity
while minimizing the costs incurred by the introduction of societal control
measures (such as closing schools, social distancing, lockdowns, etc.), under
the constraint that the infected fraction of the population does never exceed a
certain maximum corresponding to public health system capacity. Optimality is
derived and verified by variational and numerical methods for a number of model
cost functions. The effects of immune response decay after recovery are taken
into account and discussed in terms of the feasibility of strategies based on
herd immunity.


  Coronavirus Disease of 2019 (COVID-19) created dire consequences globally and
triggered an intense scientific effort from different domains. The resulting
publications created a huge text collection in which finding the studies
related to a biomolecule of interest is challenging for general purpose search
engines because the publications are rich in domain specific terminology. Here,
we present Vapur: an online COVID-19 search engine specifically designed to
find related protein - chemical pairs. Vapur is empowered with a
relation-oriented inverted index that is able to retrieve and group studies for
a query biomolecule with respect to its related entities. The inverted index of
Vapur is automatically created with a BioNLP pipeline and integrated with an
online user interface. The online interface is designed for the smooth
traversal of the current literature by domain researchers and is publicly
available at https://tabilab.cmpe.boun.edu.tr/vapur/ .


  Mechanistic models fit to streaming surveillance data are critical to
understanding the transmission dynamics of an outbreak as it unfolds in
real-time. However, transmission model parameter estimation can be imprecise,
and sometimes even impossible, because surveillance data are noisy and not
informative about all aspects of the mechanistic model. To partially overcome
this obstacle, Bayesian models have been proposed to integrate multiple
surveillance data streams. We devised a modeling framework for integrating
SARS-CoV-2 diagnostics test and mortality time series data, as well as
seroprevalence data from cross-sectional studies, and tested the importance of
individual data streams for both inference and forecasting. Importantly, our
model for incidence data accounts for changes in the total number of tests
performed. We model the transmission rate, infection-to-fatality ratio, and a
parameter controlling a functional relationship between the true case incidence
and the fraction of positive tests as time-varying quantities and estimate
changes of these parameters nonparametrically. We compare our base model
against modified versions which do not use diagnostics test counts or
seroprevalence data to demonstrate the utility of including these often unused
data streams. We apply our Bayesian data integration method to COVID-19
surveillance data collected in Orange County, California between March 2020 and
February 2021 and find that 32--72\% of the Orange County residents experienced
SARS-CoV-2 infection by mid-January, 2021. Despite this high number of
infections, our results suggest that the abrupt end of the winter surge in
January 2021 was due to both behavioral changes and a high level of accumulated
natural immunity.


  The ultimate secret of all lives on earth is hidden in their genomes -- a
totality of DNA sequences. We currently know the whole genome sequence of many
organisms, while our understanding of the genome architecture on a systematic
level remains rudimentary. Applied category theory opens a promising way to
integrate the humongous amount of heterogeneous informations in genomics, to
advance our knowledge regarding genome organization, and to provide us with a
deep and holistic view of our own genomes. In this work we explain why applied
category theory carries such a hope, and we move on to show how it could
actually do so, albeit in baby steps. The manuscript intends to be readable to
both mathematicians and biologists, therefore no prior knowledge is required
from either side.


  Dengue fever impacts populations across the tropics. Dengue is caused by a
mosquito transmitted flavivirus and its burden is projected to increase under
future climate and development scenarios. The transmission process of dengue
virus is strongly moderated by hydro-climatic conditions that impact the
vector's life cycle and behavior. Here, we study the impact of rainfall
seasonality and moisture availability on the monthly distribution of reported
dengue cases in Sri Lanka. Through cluster analysis, we find an association
between seasonal peaks of rainfall and dengue incidence with a two-month lag.
We show that a hydrologically driven epidemiological model (HYSIR), which takes
into account hydrologic memory in addition to the nonlinear dynamics of the
transmission process, captures the two-month lag between rainfall and dengue
cases seasonal peaks. Our analysis reveals a non-monotonic dependence of dengue
cases on moisture, whereby an increase of cases with increasing moisture is
followed by a reduction for very high levels of water availability. Improvement
in prediction of the seasonal peaks in dengue incidence results from a
seasonally varying dependence of transmission rate on water availability.


  Relevance: Laboratory diagnosis of the novel coronavirus (SARS-CoV-2)
infection combined with quarantine for contacts of infected individuals affects
the spread of SARS-CoV-2 infection and levels of related mortality. Practices
for testing for SARS-CoV-2 infection vary geographically in Russia. For
example, in the city of St. Petersburg, where mortality rate for COVID-19 is
the highest in the Russian Federation on Oct. 25, 2020, every death for
COVID-19 corresponds to 15.7 detected cases of COVID-19 in the population,
while the corresponding number for the whole of Russia is 58.1, suggesting
limited detection of mild/moderate cases of COVID-19 in St. Petersburg.
Methods: More active testing for SARS-CoV-2 results in lower case-fatality
ratio (i.e. the proportion of detected COVID-19 cases among all cases of
SARS-CoV-2 infection in the population). We used data on COVID-19 cases and
deaths to examine the correlation between case-fatality ratios and rates of
mortality for COVID-19 in different regions of the Russian Federation. Results:
The correlation between case-fatality ratios and rates of mortality for
COVID-19 in different regions of the Russian Federation on Oct. 25, 2020 is
0.64 (0.50,0.75). For several regions of the Russian Federation, detectability
of SARS-CoV-2 infection is relatively low, while rates of mortality for
COVID-19 are relatively high. Conclusions: Detectability of the SARS-CoV-2
infection is one of the factors that affects the levels of mortality from
COVID-19. To increase detectability, one ought to test all individuals with
respiratory symptoms seeking medical care for SARS-CoV-2 infection, and to
undertake additional measures to increase the volume of testing for SARS-CoV-2.
Such measures, in combination with quarantine for infected cases and their
close contacts help to mitigate the spread of the SARS-CoV-2 infection and
diminish the related mortality.


  The pandemic caused by the novel Coronavirus SARS-CoV2 has been responsible
for life threatening health complications, and extreme pressure on healthcare
systems. While preventive and definite curative medical interventions are yet
to arrive, Non-Pharmaceutical Interventions (NPIs) like physical isolation,
quarantine and drastic social measures imposed by governing agencies are
effective in arresting the spread of infections in a population. In densely
populated countries like India, lockdown interventions are partially effective
due to social and administrative complexities. Using detailed demographic data,
we present an agent based model to imitate the behavior of the population and
its mobility features, even under intervention. We demonstrate the
effectiveness of contact tracing policies and how our model efficiently relates
to empirical findings on testing efficiency. We also present various lockdown
intervention strategies for mitigation - using the bare number of infections,
the effective reproduction rate, as well as using reinforcement learning. Our
analysis can help assess the socio-economic consequences of such interventions,
and provide useful ideas and insights to policy makers for better decision
making.


  Following the spread of the COVID-19 pandemic and pending the establishment
of vaccination campaigns, several non pharmaceutical interventions such as
partial and full lockdown, quarantine and measures of physical distancing have
been imposed in order to reduce the spread of the disease and to lift the
pressure on healthcare system. Mathematical models are important tools for
estimating the impact of these interventions, for monitoring the current
evolution of the epidemic at a national level and for estimating the potential
long-term consequences of relaxation of measures. In this paper, we model the
evolution of the COVID-19 epidemic in Belgium with a deterministic
age-structured extended compartmental model. Our model takes special
consideration for nursing homes which are modelled as separate entities from
the general population in order to capture the specific delay and dynamics
within these entities. The model integrates social contact data and is fitted
on hospitalisations data (admission and discharge), on the daily number of
COVID-19 deaths (with a distinction between general population and nursing home
related deaths) and results from serological studies, with a sensitivity
analysis based on a Bayesian approach. We present the situation as in November
2020 with the estimation of some characteristics of the COVID-19 deduced from
the model. We also present several mid-term and long-term projections based on
scenarios of reinforcement or relaxation of social contacts for different
general sectors, with a lot of uncertainties remaining.


  Recent developments in Omics-technologies revolutionized the investigation of
biology by producing molecular data in multiple dimensions and scale. This
breakthrough in biology raises the crucial issue of their interpretation based
on modelling. In this undertaking, network provides a suitable framework for
modelling the interactions between molecules. Basically a Biological network is
composed of nodes referring to the components such as genes or proteins, and
the edges/arcs formalizing interactions between them. The evolution of the
interactions is then modelled by the definition of a dynamical system. Among
the different categories of network, the Boolean network offers a reliable
qualitative framework for the modelling. Automatically synthesizing a Boolean
network from experimental data therefore remains a necessary but challenging
issue. In this study, we present taboon, an original work-flow for synthesizing
Boolean Networks from biological data. The methodology uses the data in the
form of Boolean profiles for inferring all the potential local formula
inference. They combine to form the model space from which the most truthful
model with regards to biological knowledge and experiments must be found. In
the taboon work-flow the selection of the fittest model is achieved by a
Tabu-search algorithm. taboon is an automated method for Boolean Network
inference from experimental data that can also assist to evaluate and optimize
the dynamic behaviour of the biological networks providing a reliable platform
for further modelling and predictions.


  Many systems across the sciences evolve through a combination of
multiplicative growth and diffusive transport. In the presence of disorder,
these systems tend to form localized structures which alternate between long
periods of relative stasis and short bursts of activity. This behaviour, known
as intermittency in physics and punctuated equilibrium in evolutionary theory,
is difficult to forecast; in particular there is no general principle to locate
the regions where the system will settle, how long it will stay there, or where
it will jump next. Here I introduce a predictive theory of linear intermittency
that closes these gaps. I show that any positive linear system can be mapped
onto a generalization of the "maximal entropy random walk", a Markov process on
graphs with non-local transition rates. This construction reveals the
localization islands as local minima of an effective potential, and
intermittent jumps as barrier crossings in that potential. My results unify the
concepts of intermittency in linear systems and Markovian metastability, and
provide a generally applicable method to reduce, and predict, the dynamics of
disordered linear systems. Applications span physics, evolutionary dynamics and
epidemiology.


  We demonstrate that a simple model based on reaction-diffusion-advection
(RDA) equation forced by realistic surface velocities and nutrients is skilled
in reproducing the distributions of the surface phytoplankton chlorophyll in
the tropical Pacific. We use the low-complexity RDA model to investigate the
scale-relationships in the impact of different drivers (turbulent diffusion,
mean and eddy advection, primary productivity) on the phytoplankton chlorophyll
concentrations. We find that in the 1/4{\deg} (~25km) model, advection has a
substantial impact on the rate of primary productivity, whilst the turbulent
diffusion term has a fairly negligible impact. Turbulent diffusion has an
impact on the phytoplankton variability, with the impact being scale-propagated
and amplified by the larger scale surface currents. We investigate the impact
of a surface nutrient decline and some changes to mesoscale eddy kinetic energy
(climate change projections) on the surface phytoplankton concentrations. The
RDA model suggests that unless mesoscale eddies radically change, phytoplankton
chlorophyll scales sub-linearly with the nutrients, and it is relatively stable
with respect to the nutrient concentrations. Furthermore we explore how a white
multiplicative Gaussian noise introduced into the RDA model on its resolution
scale propagates across spatial scales through the non-linear model dynamics
under different sets of phytoplankton drivers. The unifying message of this
work is that the low complexity (e.g. RDA) models can be successfully used to
realistically model some specific aspects of marine ecosystem dynamics and by
using those models one can explore many questions that would be beyond
computational affordability of the higher-complexity ecosystem models.


  Recent seroprevalence studies have tried to estimate the real number of
asymptomatic cases affected by COVID-19. It is of paramount importance to
understand the impact of these infections in order to prevent a second wave.
This study aims to model the interactions in the population by means of a
complex network and to shed some light on the effectiveness of localised
control measures in Italy in relation to the school opening in mid-September.
The formulation of an epidemiological predictive model is given: the advantage
of using this model lies in that it discriminates between asymptomatic and
symptomatic cases of COVID-19 as the interactions with these two categories of
infected individuals are captured separately, allowing for a study on the
impact of asymptomatic cases. This model is then extended to a structured
nonhomogeneous version by means of the Watts-Strogatz complex network, which is
adopted widely to model societal interactions as it holds the small world
property. Finally, a case study on the situation in Italy is given: first the
homogeneous model is used to compare the official data with the data of the
recent seroprevalence study from Istat; second, in view of the return to school
in mid-September, a study at regional level is conducted. The results of this
study highlight the importance of coordinating the deployment of appropriate
control measures that take into account the role of asymptomatic infections,
especially in younger individuals, and inter-regional connectivity in Italy.


  Optimizing the impact on the economy of control strategies aiming at
containing the spread of COVID-19 is a critical challenge. We use daily new
case counts of COVID-19 patients reported by local health administrations from
different Metropolitan Statistical Areas (MSAs) within the US to parametrize a
model that well describes the propagation of the disease in each area. We then
introduce a time-varying control input that represents the level of social
distancing imposed on the population of a given area and solve an optimal
control problem with the goal of minimizing the impact of social distancing on
the economy in the presence of relevant constraints, such as a desired level of
suppression for the epidemics at a terminal time. We find that with the
exception of the initial time and of the final time, the optimal control input
is well approximated by a constant, specific to each area, which contrasts with
the implemented system of reopening `in phases'. For all the areas considered,
this optimal level corresponds to stricter social distancing than the level
estimated from data. Proper selection of the time period for application of the
control action optimally is important: depending on the particular MSA this
period should be either short or long or intermediate. We also consider the
case that the transmissibility increases in time (due e.g. to increasingly
colder weather), for which we find that the optimal control solution yields
progressively stricter measures of social distancing. {We finally compute the
optimal control solution for a model modified to incorporate the effects of
vaccinations on the population and we see that depending on a number of
factors, social distancing measures could be optimally reduced during the
period over which vaccines are administered to the population.


  The generation of high-quality assemblies, even for large eukaryotic genomes,
has become a routine task for many biologists thanks to recent advances in
sequencing technologies. However, the annotation of these assemblies - a
crucial step towards unlocking the biology of the organism of interest - has
remained a complex challenge that often requires advanced bioinformatics
expertise. Here we present MOSGA, a genome annotation framework for eukaryotic
genomes with a user-friendly web-interface that generates and integrates
annotations from various tools. The aggregated results can be analyzed with a
fully integrated genome browser and are provided in a format ready for
submission to NCBI. MOSGA is built on a portable, customizable, and easily
extendible Snakemake backend, and thus, can be tailored to a wide range of
users and projects. We provide MOSGA as a publicly free available web service
at https://mosga.mathematik.uni-marburg.de and as a docker container at
registry.gitlab.com/mosga/mosga:latest. Source code can be found at
https://gitlab.com/mosga/mosga


  The goal of this work is to consider widespread use of face masks as a
non-pharmaceutical control strategy for the Covid-19 pandemic. A SEIR model
that divides the population into individuals that wear masks and those that do
not is considered. After calculating the basic reproductive number by a next
generation approach, a criterion for determining when an epidemic can be
prevented by the use of masks only and the critical percentage of mask users
for disease prevention in the population are derived. The results are then
applied to real world data from the United States, Brazil and Italy.


  Elucidating physical mechanisms with statistical confidence from molecular
dynamics simulations can be challenging owing to the many degrees of freedom
that contribute to collective motions. To address this issue, we recently
introduced a dynamical Galerkin approximation (DGA) [Thiede et al. J. Phys.
Chem. 150, 244111 (2019)], in which chemical kinetic statistics that satisfy
equations of dynamical operators are represented by a basis expansion. Here, we
reformulate this approach, clarifying (and reducing) the dependence on the
choice of lag time. We present a new projection of the reactive current onto
collective variables and provide improved estimators for rates and committors.
We also present simple procedures for constructing suitable smoothly varying
basis functions from arbitrary molecular features. To evaluate estimators and
basis sets numerically, we generate and carefully validate a dataset of short
trajectories for the unfolding and folding of the trp-cage miniprotein, a
well-studied system. Our analysis demonstrates a comprehensive strategy for
characterizing reaction pathways quantitatively.


  A networked dynamical system is composed of subsystems interconnected through
prescribed interactions. In many engineering applications, however, one
subsystem can also affect others through "unintended" interactions that can
significantly hamper the intended network's behavior. Although unintended
interactions can be modeled as disturbance inputs to the subsystems, these
disturbances depend on the network's states. As a consequence, a disturbance
attenuation property of each isolated subsystem is, alone, insufficient to
ensure that the network behavior is robust to unintended interactions. In this
paper, we provide sufficient conditions on subsystem dynamics and interaction
maps, such that the network's behavior is robust to unintended interactions.
These conditions require that each subsystem attenuates constant external
disturbances, is monotone or "near-monotone", the unintended interaction map is
monotone, and the prescribed interaction map does not contain feedback loops.
We employ this result to guide the design of resource-limited genetic circuits.
More generally, our result provide conditions under which robustness of
constituent subsystems is sufficient to guarantee robustness of the network to
unintended interactions.


  Network theory-based approaches provide valuable insights into the variations
in global structural connectivity between differing dynamical states of
proteins. Our objective is to review network-based analyses to elucidate such
variations, especially in the context of subtle conformational changes. We
present technical details of the construction and analyses of protein structure
networks, encompassing both the non-covalent connectivity and dynamics. We
examine the selection of optimal criteria for connectivity based on the
physical concept of percolation. We highlight the advantages of using
side-chain based network metrics in contrast to backbone measurements. As an
illustrative example, we apply the described network approach to investigate
the global conformational change between the closed and partially open states
of the SARS-CoV-2 spike protein. This conformational change in the spike
protein is crucial for coronavirus entry and fusion into human cells. Our
analysis reveals global structural reorientations between the two states of the
spike protein despite small changes between the two states at the backbone
level. We also observe some differences at strategic locations in the
structures, correlating with their functions, asserting the advantages of the
side-chain network analysis. Finally we present a view of allostery as a subtle
synergistic-global change between the ligand and the receptor, the
incorporation of which would enhance the drug design strategies.


  We propose a mathematical model to analyze the effects of anti-infection
behavior on the equilibrium states of an infectious disease. The anti-infection
behavior is incorporated into a classical epidemiological SIR model, by
considering the behavior adoption rate across the population as an additional
variable. We consider also the effects on the adoption rate produced by the
disease evolution, using a dynamic payoff function and an additional
differential equation. The equilibrium states of the proposed model have
remarkable characteristics: possible coexistence of two locally stable endemic
equilibria, the coexistence of locally stable endemic and disease-free
equilibria, and even the possibility of a stable continuum of endemic
equilibrium points. We show how some of the results obtained may be used to
support strategic planning leading to effective control of the disease in the
long-term.


  This work examines the use of reinforcement learning (RL) to optimize cyclic
lockdowns, which is one of the methods available for control of the COVID-19
pandemic. The problem is structured as an optimal control system for tracking a
reference value, corresponding to the maximum usage level of a critical
resource, such as ICU beds. However, instead of using conventional optimal
control methods, RL is used to find optimal control policies. A framework was
developed to calculate optimal cyclic lockdown timings using an RL-based on-off
controller. The RL-based controller is implemented as an RL agent that
interacts with an epidemic simulator, implemented as an extended SEIR epidemic
model. The RL agent learns a policy function that produces an optimal sequence
of open/lockdown decisions such that goals specified in the RL reward function
are optimized. Two concurrent goals were used: the first one is a public health
goal that minimizes overshoots of ICU bed usage above an ICU bed threshold, and
the second one is a socio-economic goal that minimizes the time spent under
lockdowns. It is assumed that cyclic lockdowns are considered as a temporary
alternative to extended lockdowns when a region faces imminent danger of
overpassing resource capacity limits and when imposing an extended lockdown
would cause severe social and economic consequences due to lack of necessary
economic resources to support its affected population during an extended
lockdown.


  In the last months, due to the emergency of Covid-19, questions related to
the fact of belonging or not to a particular class of individuals (`infected or
not infected'), after being tagged as `positive' or `negative' by a test, have
never been so popular. Similarly, there has been strong interest in estimating
the proportion of a population expected to hold a given characteristics
(`having or having had the virus'). Taking the cue from the many related
discussions on the media, in addition to those to which we took part, we
analyze these questions from a probabilistic perspective (`Bayesian'),
considering several effects that play a role in evaluating the probabilities of
interest. The resulting paper, written with didactic intent, is rather general
and not strictly related to pandemics: the basic ideas of Bayesian inference
are introduced and the uncertainties on the performances of the tests are
treated using the metrological concepts of `systematics', and are propagated
into the quantities of interest following the rules of probability theory; the
separation of `statistical' and `systematic' contributions to the uncertainty
on the inferred proportion of infectees allows to optimize the sample size; the
role of `priors', often overlooked, is stressed, however recommending the use
of `flat priors', since the resulting posterior distribution can be `reshaped'
by an `informative prior' in a later step; details on the calculations are
given, also deriving useful approximated formulae, the tough work being however
done with the help of direct Monte Carlo simulations and Markov Chain Monte
Carlo, implemented in R and JAGS (relevant code provided in appendix).


  Currently, many studies of Alzheimer's disease (AD) are investigating the
neurobiological factors behind the acquisition of beta-amyloid (A), pathologic
tau (T), and neurodegeneration ([N]) biomarkers from neuroimages. However, a
system-level mechanism of how these neuropathological burdens promote
neurodegeneration and why AD exhibits characteristic progression is largely
elusive. In this study, we combined the power of systems biology and network
neuroscience to understand the dynamic interaction and diffusion process of
AT[N] biomarkers from an unprecedented amount of longitudinal Amyloid PET scan,
MRI imaging, and DTI data. Specifically, we developed a network-guided
biochemical model to jointly (1) model the interaction of AT[N] biomarkers at
each brain region and (2) characterize their propagation pattern across the
fiber pathways in the structural brain network, where the brain resilience is
also considered as a moderator of cognitive decline. Our biochemical model
offers a greater mathematical insight to understand the physiopathological
mechanism of AD progression by studying the system dynamics and stability.
Thus, an in-depth system-level analysis allows us to gain a new understanding
of how AT[N] biomarkers spread throughout the brain, capture the early sign of
cognitive decline, and predict the AD progression from the preclinical stage.


  The Sturgis Motorcycle Rally that took place from August 7-16 was one of the
largest public gatherings since the start of the COVID-19 outbreak. Over
460,000 visitors from across the United States travelled to Sturgis, South
Dakota to attend the ten day event. Using anonymous cell phone tracking data we
identify the home counties of visitors to the rally and examine the impact of
the rally on the spread of COVID-19. Our baseline estimate suggests a one
standard deviation increase in Sturgis attendance increased COVID-19 case
growth by 1.1pp in the weeks after the rally.


  COVID-19 is not a universal killer. We study the spread of COVID-19 at the
county level for the United States up until the 15$^{th}$ of August, 2020. We
show that the prevalence of the disease and the death rate are correlated with
the local socio-economic conditions often going beyond local population density
distributions, especially in rural areas. We correlate the COVID-19 prevalence
and death rate with data from the US Census Bureau and point out how the
spreading patterns of the disease show asymmetries in urban and rural areas
separately and are preferentially affecting the counties where a large fraction
of the population is non-white. Our findings can be used for more targeted
policy building and deployment of resources for future occurrence of a pandemic
due to SARS-CoV-2. Our methodology, based on interpretable machine learning and
game theory, can be extended to study the spread of other diseases.


  In this study, we address three important challenges related to disease
transmissions such as the COVID-19 pandemic, namely, (a) providing an early
warning to likely exposed individuals, (b) identifying individuals who are
asymptomatic, and (c) prescription of optimal testing when testing capacity is
limited. First, we present a dynamic-graph based SEIR epidemiological model in
order to describe the dynamics of the disease propagation. Our model considers
a dynamic network that accounts for the interactions between individuals over
time, such as the ones obtained by manual or automated contact tracing, and
uses a diffusion-reaction mechanism to describe the state dynamics. This
dynamic graph model helps identify likely exposed/infected individuals to whom
we can provide early warnings, even before they display any symptoms and/or are
asymptomatic. Moreover, when the testing capacity is limited compared to the
population size, reliable estimation of individual's health state and disease
transmissibility using epidemiological models is extremely challenging. Thus,
estimation of state uncertainty is paramount for both eminent risk assessment,
as well as for closing the tracing-testing loop by optimal testing
prescription. Therefore, we propose the use of arbitrary Polynomial Chaos
Expansion, a popular technique used for uncertainty quantification, to
represent the states, and quantify the uncertainties in the dynamic model. This
design enables us to assign uncertainty of the state of each individual, and
consequently optimize the testing as to reduce the overall uncertainty given a
constrained testing budget. These tools can also be used to optimize vaccine
distribution to curb the disease spread when limited vaccines are available. We
present a few simulation results that illustrate the performance of the
proposed framework, and estimate the impact of incomplete contact tracing data.


  Species living in the extreme cold environment fight against the harsh
conditions using antifreeze proteins (AFPs), that manipulates the freezing
mechanism of water in more than one way. This amazing nature of AFP turns out
to be extremely useful in several industrial and medical applications. The lack
of similarity in their structure and sequence makes their prediction an arduous
task and identifying them experimentally in the wet-lab is time-consuming and
expensive. In this research, we propose a computational framework for the
prediction of AFPs which is essentially based on a sample-specific
classification method using the sparse reconstruction. A linear model and an
over-complete dictionary matrix of known AFPs are used to predict a sparse
class-label vector that provides a sample-association score. Delta-rule is
applied for the reconstruction of two pseudo-samples using lower and upper
parts of the sample-association vector and based on the minimum recovery score,
class labels are assigned. We compare our approach with contemporary methods on
a standard dataset and the proposed method is found to outperform in terms of
Balanced accuracy and Youden's index. The MATLAB implementation of the proposed
method is available at the author's GitHub page
(\{https://github.com/Shujaat123/AFP-SRC}{https://github.com/Shujaat123/AFP-SRC}).


  Covid-19 is raging a devastating trail with the highest mortality-to-infected
ratio ever for a pandemic. Lack of vaccine and therapeutic has rendered social
exclusion through lockdown as the singular mode of containment. Harnessing the
predictive powers of Machine Learning within a 6 dimensional infection kinetic
model, depicting interactive evolution of 6 infection stages - healthy
susceptible ($H$), predisposed comorbid susceptible ($P$), infected ($I$),
recovered ($R$), herd immunized ($V$) and mortality ($D$) - the model, PHIRVD,
provides the first accurate mortality prediction of 18 countries at varying
stages of strategic lockdown, up to 30 days beyond last data training. PHIRVD
establishes mortality-to-infection ratio as the correct pandemic descriptor,
substituting reproduction number, and highlights the importance of early and
prolonged but strategic lockdown to contain secondary relapse.
  Significance Statement:
  1. Accurate prediction of the day-by-day mortality profiles of 18 countries,
30 days beyond the last data of data training.
  2. Precise quantification of the impact of early-vs-later lockdown
impositions.
  3. Accurate prediction of secondary relapse timelines/
  4. Establishment of mortality-to-infected ratio as the correct pandemic
descriptor substituting the popular choice of reproduction number, a proven
failure in predicting future infection kinetics and secondary surge.
  The outcomes have potential to redefine healthy policy landscape,
particularly in light of secondary relapse and possible future SARS-COV/Ebola
group incursion.


  Intent of this research is to explore how mathematical models, specifically
Susceptible-Infected-Removed (SIR) model, can be utilized to forecast peak
outbreak timeline of COVID-19 epidemic amongst a population of interest
starting from the date of first reported case. Till the time of this research,
there was no effective and universally accepted vaccine to control transmission
and spread of this infection. COVID-19 primarily spreads in population through
respiratory droplets from an infected person cough and sneeze which infects
people who are in proximity. COVID-19 is spreading contagiously across the
world. If health policy makers and medical experts could get early and timely
insights into when peak infection rate would occur after first reported case,
they could plan and optimize medical personnel, ventilators supply, and other
medical resources without over-taxing the infrastructure. The predictions may
also help policymakers devise strategies to control the epidemic, potentially
saving many lives. Thus, it can aid in critical decision-making process by
providing actionable insights into COVID-19 outbreak by leveraging available
data.


  By incorporating delayed epidemic information and self-restricted travel
behavior into the SIS model, we have investigated the coupled effects of timely
and accurate epidemic information and people's sensitivity to the epidemic
information on contagion. In the population with only local random movement,
whether the epidemic information is delayed or not has no effect on the spread
of the epidemic. People's high sensitivity to the epidemic information leads to
their risk aversion behavior and the spread of the epidemic is suppressed. In
the population with only global person-to-person movement, timely and accurate
epidemic information helps an individual cut off the connections with the
infected in time and the epidemic is brought under control in no time. A delay
in the epidemic information leads to an individual's misjudgment of who has
been infected and who has not, which in turn leads to rapid progress and a
higher peak of the epidemic. In the population with coexistence of local and
global movement, timely and accurate epidemic information and people's high
sensitivity to the epidemic information play an important role in curbing the
epidemic. A theoretical analysis indicates that people's misjudgment caused by
the delayed epidemic information leads to a higher encounter probability
between the susceptible and the infected and people's self-restricted travel
behavior helps reduce such an encounter probability. A functional relation
between the ratio of infected individuals and the susceptible-infected
encounter probability has been found.


  Observable social traits determine how we interact in society and remain
pervasive even in our globalized world. While a popular hypothesis states that
they may help promote cooperation, the alternative explanation that they
facilitate coordination has gained ground in recent years. Here we explore this
framework and present a model that investigates the role of ethnic markers in
coordination games. We consider fixed markers characterizing agents that use
reinforcement learning to update their strategies in the game. For a wide range
of parameters, we observe the emergence of a collective equilibrium in which
markers play an assorting role. However, if individuals are too conformists or
greedy, markers fail to shape social interactions. These results extend and
complement previous work focused on agent imitation and show that reinforcement
learning is a good candidate to explain many instances of ethnic markers.


  Without a cure, vaccine, or proven long-term immunity against SARS-CoV-2,
test-trace-and-isolate (TTI) strategies present a promising tool to contain its
spread. For any TTI strategy, however, mitigation is challenged by pre- and
asymptomatic transmission, TTI-avoiders, and undetected spreaders, who strongly
contribute to hidden infection chains. Here, we studied a semi-analytical model
and identified two tipping points between controlled and uncontrolled spread:
(1) the behavior-driven reproduction number of the hidden chains becomes too
large to be compensated by the TTI capabilities, and (2) the number of new
infections exceeds the tracing capacity. Both trigger a self-accelerating
spread. We investigated how these tipping points depend on challenges like
limited cooperation, missing contacts, and imperfect isolation. Our model
results suggest that TTI alone is insufficient to contain an otherwise
unhindered spread of SARS-CoV-2, implying that complementary measures like
social distancing and improved hygiene remain necessary.


  Large scale virological testing of SARS-Cov2 is implemented since May 2020 in
France. We assume that the positivity of asymptomatic people not being contact
cases (ANBC) is representative of the positivity of the whole French
population, which allows estimating the real incidence. We estimate, using
Sant\'e Public France reports, that the incidence at the beginning of August
was about 0.8% and rose to about 2.4% of the total population at the beginning
of September. This corresponds to about 1.6 million people simultaneously
infected and 230.000 new infections each day. These evaluations allow to deduce
that intensive care units (ICU) admission rate and infection fatality rate
(IFR) dropped by one order of magnitude since March, and are currently 0.036%
and 0.027% respectively. Basic simulations of the outbreak evolution based on
the hypothesis of negligible reinfection probability are performed for France,
Ile de France, Puy de Dome, Bouches du Rhone and Grand Est. These simulations
are using a reproduction rate (R) constant over time ranging from 1.3 to 1.45
(1.15 in Grand Est). They use an estimation of the number of infection which
occurred during the first wave. These simulations also use the estimated
incidence, which by reducing the susceptible population weeks after weeks
induces a saturation of the otherwise exponential growth of the outbreak. An
incidence peak of 3.5 % is expected at week 39 for France. The calculated total
number of ICU admission and deaths during the second wave are 9000 and 7000
respectively for R=1.3. The cumulative incidence over the two waves is computed
close to 60% for R=1.3 and 70% for R=1.4. This suggests that if individual
immunity exists, herd immunity is likely to be achieved in France by the end of
October 2020. We conclude that Covid-19 is much more spread than previously
thought, but its severity became limited since the end of the first wave.


  Machine learning can precisely identify different cancer tumors at any stage
by classifying cancerous and healthy samples based on their genomic profile. We
have developed novel methods of MLAC (Machine Learning Against Cancer)
achieving perfect results with perfect precision, sensitivity, and specificity.
We have used the whole genome sequencing data acquired by next-generation RNA
sequencing techniques in The Cancer Genome Atlas and Genotype-Tissue Expression
projects for cancerous and healthy tissues respectively. Moreover, we have
shown that unsupervised machine learning clustering has great potential to be
used for cancer diagnosis. Indeed, a creative way to work with data and general
algorithms has resulted in perfect classification i.e. all precision,
sensitivity, and specificity are equal to 1 for most of the different tumor
types even with a modest amount of data, and the same method works well on a
series of cancers and results in great clustering of cancerous and healthy
samples too. Our system can be used in practice because once the classifier is
trained, it can be used to classify any new sample of new potential patients.
One advantage of our work is that the aforementioned perfect precision and
recall are obtained on samples of all stages including very early stages of
cancer; therefore, it is a promising tool for diagnosis of cancers in early
stages. Another advantage of our novel model is that it works with normalized
values of RNA sequencing data, hence people's private sensitive medical data
will remain hidden, protected, and safe. This type of analysis will be
widespread and economical in the future and people can even learn to receive
their RNA sequencing data and do their own preliminary cancer studies
themselves which have the potential to help the healthcare systems. It is a
great step forward toward good health that is the main base of sustainable
societies.


  Motivation: With the aim to amplify and make sense of interactions of
virus-human proteins in the case of SARS-CoV-2, we performed a structural
analysis of the network of protein interactions obtained from the integration
of three sources: 1) proteins of virus SARS-CoV-2, 2)physical interactions
between SARS-CoV-2 and human proteins, 3) known interactions of these human
proteins between them and the dossier of affections in which these proteins are
implicated.
  Results: As a product of this research, we present two networks, one from the
interactions virus-host, and the other restricted to host-host, the last one is
not usually considered for network analysis. We identified the most important
proteins in both networks, those that have the maximal value of calculated
invariants, these proteins are considered as the most: affected, connected or
those that best monitor the flow of information in the network, among them we
find UBC, a human protein related with ubiquination, linked with different
stages of coronavirus disease, and ORF7A a virus protein that induces apoptosis
in infected cells, associated with virion tethering. Using the constructed
networks, we establish the more significant diseases corresponding with human
proteins and their connections with other proteins. It is relevant that the
identified diseases coincide with comorbidities, particularly the subnetwork of
diabetes involves a great quantity of virus and human proteins (56%) and
interactions (60%), this could explain the effect of this condition as an
important cause of disease complications.


  Multisite phosphorylation plays an important role in regulating switchlike
protein activity and has been used widely in mathematical models. With the
development of new experimental techniques and more molecular data, molecular
phosphorylation processes emerge in many systems with increasing complexity and
sizes. These developments call for simple yet valid stochastic models to
describe various multisite phosphorylation processes, especially in large and
complex biochemical networks. To reduce model complexity, this work aims to
simplify the multisite phosphorylation mechanism by a stochastic Hill function
model. Further, this work optimizes regions of parameter space to match
simulation results from the stochastic Hill function with the distributive
multisite phosphorylation process. While traditional parameter optimization
methods have been focusing on finding the best parameter vector, in most
circumstances modelers would like to find a set of parameter vectors that
generate similar system dynamics and results. This paper proposes a general
$\alpha$-$\beta$-$\gamma$ rule to return an acceptable parameter region of the
stochastic Hill function based on a quasi-Newton stochastic optimization
(QNSTOP) algorithm. Different objective functions are investigated
characterizing different features of the simulation-based empirical data, among
which the approximate maximum log-likelihood method is recommended for general
applications. Numerical results demonstrate that with an appropriate parameter
vector value, the stochastic Hill function model depicts the multisite
phosphorylation process well except the initial (transient) period.


  Bone repair using BMP-2 is a promising therapeutic approach in clinical
practices, however, high dosages required to be effective pose issues of cost
and safety. The present study explores the potential of low dose BMP-2
treatment via tissue engineering approach, which amalgamates 3-D
macro/microporous-nanofibrous bacterial cellulose (mNBC) scaffolds and low dose
BMP-2 primed murine mesenchymal stem cells (C3H10T1/2 cells). Initial studies
on cell-scaffold interaction using unprimed C3H10T1/2 cells confirmed that
scaffolds provided a propitious environment for cell adhesion, growth, and
infiltration, owing to its ECM-mimicking nano-micro-macro architecture.
Osteogenic studies were conducted by preconditioning the cells with 50 ng/mL
BMP-2 for 15 minutes, followed by culturing on mNBC scaffolds for up to three
weeks. The results showed an early onset and significantly enhanced bone matrix
secretion and maturation in the scaffolds seeded with BMP-2 primed cells
compared to the unprimed ones. Moreover, mNBC scaffolds alone were able to
facilitate the mineralization of cells to some extent. These findings suggest
that, with the aid of 'osteoinduction' from low dose BMP-2 priming of stem
cells and 'osteoconduction' from nano-macro/micro topography of mNBC scaffolds,
a cost-effective bone tissue engineering strategy can be designed for quick and
excellent in vivo osseointegration.


  The last leg of the year 2019 gave rise to a virus named COVID-19 (Corona
Virus Disease 2019). Since the beginning of this infection in India, the
government implemented several policies and restrictions to curtail its spread
among the population. As the time passed, these restrictions were relaxed and
people were advised to follow precautionary measures by themselves. These
timely decisions taken by the Indian government helped in decelerating the
spread of COVID-19 to a large extent. Despite these decisions, the pandemic
continues to spread and hence, there is an urgent need to plan and control the
spread of this disease. This is possible by finding the future predictions
about the spread. Scientists across the globe are working towards estimating
the future growth of COVID-19. This paper proposes a Modified SEIRD
(Susceptible-Exposed-Infected-Recovered-Deceased) model for projecting COVID-19
infections in India and its five states having the highest number of total
cases. In this model, exposed compartment contains individuals which may be
asymptomatic but infectious. Deep Learning based Long Short-Term Memory (LSTM)
model has also been used in this paper to perform short-term projections. The
projections obtained from the proposed Modified SEIRD model have also been
compared with the projections made by LSTM for next 30 days. The
epidemiological data up to 15th August 2020 has been used for carrying out
predictions in this paper. These predictions will help in arranging adequate
medical infrastructure and providing proper preventive measures to handle the
current pandemic. The effect of different lockdowns imposed by the Indian
government has also been used in modelling and analysis in the proposed
Modified SEIRD model. The results presented in this paper will act as a beacon
for future policy-making to control the COVID-19 spread in India.


  The key to a robust life system is to ensure that each cell population is
maintained in an appropriate state. In this work, a mathematical model was used
to investigate the control of the switching between the migrating and
non-migrating states of the Bacillus subtilis cell population. In this case,
the motile cells and matrix producers were the predominant cell types in the
migrating cell population and non-migrating state, respectively, and could be
suitably controlled according to the environmental conditions and cell density
information. A minimal smooth model consisting of four ordinary differential
equations was used as the mathematical model to control the B. subtilis cell
types. Furthermore, the necessary and sufficient conditions for the hysteresis,
which pertains to the change in the pheromone concentration, were clarified. In
general, the hysteretic control of the cell state enables stable switching
between the migrating and growth states of the B. subtilis cell population,
thereby facilitating the biofilm life cycle. The results of corresponding
culture experiments were examined, and the obtained corollaries were used to
develop a model to input environmental conditions, especially, the external pH.
On this basis, the environmental conditions were incorporated in a simulation
model for the cell type control. In combination with a mathematical model of
the cell population dynamics, a prediction model for colony growth involving
multiple cell states, including concentric circular colonies of B. subtilis,
could be established.


  In the context of epidemiology, policies for disease control are often
devised through a mixture of intuition and brute-force, whereby the set of
logically conceivable policies is narrowed down to a small family described by
a few parameters, following which linearization or grid search is used to
identify the optimal policy within the set. This scheme runs the risk of
leaving out more complex (and perhaps counter-intuitive) policies for disease
control that could tackle the disease more efficiently. In this article, we use
techniques from convex optimization theory and machine learning to conduct
optimizations over disease policies described by hundreds of parameters. In
contrast to past approaches for policy optimization based on control theory,
our framework can deal with arbitrary uncertainties on the initial conditions
and model parameters controlling the spread of the disease, and stochastic
models. In addition, our methods allow for optimization over policies which
remain constant over weekly periods, specified by either continuous or discrete
(e.g.: lockdown on/off) government measures. We illustrate our approach by
minimizing the total time required to eradicate COVID-19 within the
Susceptible-Exposed-Infected-Recovered (SEIR) model proposed by Kissler
\emph{et al.} (March, 2020).


  A biochemical network can be simulated by a set of ordinary differential
equations (ODE) under well stirred reactor conditions, for large numbers of
molecules, and frequent reactions. This is no longer a robust representation
when some molecular species are in small numbers and reactions changing them
are infrequent. In this case, discrete stochastic events trigger changes of the
smooth deterministic dynamics of the biochemical network.
Piecewise-deterministic Markov processes (PDMP) are well adapted for describing
such situations. Although PDMP models are now well established in biology,
these models remain computationally challenging. Previously we have introduced
the push-forward method to compute how the probability measure is spread by the
deterministic ODE flow of PDMPs, through the use of analytic expressions of the
corresponding semigroup. In this paper we provide a more general simulation
algorithm that works also for non-integrable systems. The method can be used
for biochemical simulations with applications in fundamental biology,
biotechnology and biocomputing.This work is an extended version of the work
presented at the conference CMSB2019.


  Protein conformational fluctuations are highly complex and exhibit long-term
correlations. Here, molecular dynamics simulations of small proteins
demonstrate that these conformational fluctuations directly affect the
protein's instantaneous diffusivity $D_I$. We find that the radius of gyration
$R_g$ of the proteins exhibits $1/f$ fluctuations, that are synchronous with
the fluctuations of $D_I$. Our analysis demonstrates the validity of the local
Stokes-Einstein type relation $D_I\propto1/(R_g + R_0)$, where $R_0\sim0.3$ nm
is assumed to be a hydration layer around the protein. From the analysis of
different protein types with both strong and weak conformational fluctuations
the validity of the Stokes-Einstein type relation appears to be a general
property.


  Various coarse-grained models have been proposed to study the spreading
dynamics in the network. A microscopic theory is needed to connect the
spreading dynamics with the individual behaviors. In this letter, we unify the
description of different spreading dynamics on complex networks by decomposing
the microscopic dynamics into two basic processes, the aging process and the
contact process. A microscopic dynamical equation is derived to describe the
dynamics of individual nodes on the network. The hierarchy of a duration
coarse-grained (DCG) approach is obtained to study duration-dependent
processes, where the transition rates depend on the duration of an individual
node on a state. Applied to the epidemic spreading, such formalism is feasible
to reproduce different epidemic models, e.g., the
susceptible-infected-recovered and the susceptible-infected-susceptible models,
and to associate with the corresponding macroscopic spreading parameters with
the microscopic transition rate. The DCG approach enables us to obtain the
steady state of the general SIS model with arbitrary duration-dependent
recovery and infection rates. The current hierarchical formalism can also be
used to describe the spreading of information and public opinions, or to model
a reliability theory in networks.


  Exploring the protein-folding problem has been a long-standing challenge in
molecular biology. Protein folding is highly dependent on folding of secondary
structures as the way to pave a native folding pathway. Here, we demonstrate
that a feature of a large hydrophobic surface area covering most side-chains on
one side or the other side of adjacent $\beta$-strands of a $\beta$-sheet is
prevail in almost all experimentally determined $\beta$-sheets, indicating that
folding of $\beta$-sheets is most likely triggered by multistage hydrophobic
interactions among neighbored side-chains of unfolded polypeptides, enable
$\beta$-sheets fold reproducibly following explicit physical folding codes in
aqueous environments. $\beta$-turns often contain five types of residues
characterized with relatively small exposed hydrophobic proportions of their
side-chains, that is explained as these residues can block hydrophobic effect
among neighbored side-chains in sequence. Temperature dependence of the folding
of $\beta$-sheet is thus attributed to temperature dependence of the strength
of the hydrophobicity. The hydrophobic-effect-based mechanism responsible for
$\beta$-sheets folding is verified by bioinformatics analyses of thousands of
results available from experiments. The folding codes in amino acid sequence
that dictate formation of a $\beta$-hairpin can be deciphered through
evaluating hydrophobic interaction among side-chains of an unfolded polypeptide
from a $\beta$-strand-like thermodynamic metastable state.


  Genome sequence analysis has enabled significant advancements in medical and
scientific areas such as personalized medicine, outbreak tracing, and the
understanding of evolution. Unfortunately, it is currently bottlenecked by the
computational power and memory bandwidth limitations of existing systems, as
many of the steps in genome sequence analysis must process a large amount of
data. A major contributor to this bottleneck is approximate string matching
(ASM).
  We propose GenASM, the first ASM acceleration framework for genome sequence
analysis. We modify the underlying ASM algorithm (Bitap) to significantly
increase its parallelism and reduce its memory footprint, and we design the
first hardware accelerator for Bitap. Our hardware accelerator consists of
specialized compute units and on-chip SRAMs that are designed to match the rate
of computation with memory capacity and bandwidth.
  We demonstrate that GenASM is a flexible, high-performance, and low-power
framework, which provides significant performance and power benefits for three
different use cases in genome sequence analysis: 1) GenASM accelerates read
alignment for both long reads and short reads. For long reads, GenASM
outperforms state-of-the-art software and hardware accelerators by 116x and
3.9x, respectively, while consuming 37x and 2.7x less power. For short reads,
GenASM outperforms state-of-the-art software and hardware accelerators by 111x
and 1.9x. 2) GenASM accelerates pre-alignment filtering for short reads, with
3.7x the performance of a state-of-the-art pre-alignment filter, while
consuming 1.7x less power and significantly improving the filtering accuracy.
3) GenASM accelerates edit distance calculation, with 22-12501x and 9.3-400x
speedups over the state-of-the-art software library and FPGA-based accelerator,
respectively, while consuming 548-582x and 67x less power.


  Situations involving cooperative behaviour are widespread among animals and
humans alike. Game theory and evolutionary dynamics have provided the
theoretical and computational grounds to understand what are the mechanisms
that allow for such cooperation. Studies in this area usually take into
consideration different behavioural strategies and investigate how they can be
fixed in the population under evolving rules. However, how those strategies
emerged from basic evolutionary mechanisms continues to be not fully
understood. To address this issue, here we study the emergence of cooperative
strategies through a model of heuristics selection based on evolutionary
algorithms. In the proposed model, agents interact with other players according
to a heuristic specified by their genetic code and reproduce -- at a longer
time scale -- proportionally to their fitness. We show that the system can
evolve to cooperative regimes for low mutation rates through heuristics
selection while increasing the mutation decreases the level of cooperation. Our
analysis of possible strategies shows that reciprocity and punishment are the
main ingredients for cooperation to emerge, being conditional cooperation the
more frequent strategy. Additionally, we show that if in addition to
behavioural rules, genetic relatedness is included, then kinship plays a
relevant role. Our results illustrate that our evolutionary heuristics model is
a generic and powerful tool to study the evolution of cooperative behaviour.


  Natural isolates of Bacillus subtilis are known for their ability to produce
a large panel of bioactive compounds. Unfortunately, their recalcitrance to
conventional molecular techniques limits their transcript studies. In this
work, diYculties to isolate RNA attributed to the cell wall were overcome,
Wnally authorising powerful RT-PCR's.