  A blood cell lineage consists of several consecutive developmental stages
from the pluri- or multipotent stem cell to a state of terminal
differentiation. Despite their importance for human biology, the regulatory
pathways and gene networks that govern these differentiation processes are not
yet fully understood. This is in part due to challenges associated with
delineating the interactions between transcription factors (TFs) and their
target genes. A possible path forward in this issue is provided by increasingly
available expression data as a basis for linking differentiation stages and
gene activities. Here, we present a novel hierarchical approach to identify
characteristic expression peak patterns that global regulators expose along the
differentiation path of cell lineages. Based on such simple patterns, we
identify cell state-specific marker genes and extract TFs that likely drive
their differentiation. Integration of the mean expression values of
stage-specific key player genes yields a distinct peaking pattern for each
lineage that is used to identify further genes in the dataset behaving
similarly. Incorporating the set of TFs which regulate these genes incurred at
a set of stage-specific regulators controlling the biological process of cell
fate. As proof of concept, we consider two expression datasets covering key
differentiation events in blood cell formation of mice.


  GTPases are molecular switches that regulate a wide range of cellular
processes, such as organelle biogenesis, position, shape, and signal
transduction. These enzymes operate by toggling between an active ("ON")
guanosine triphosphate (GTP)-bound state and an inactive ("OFF") guanosine
diphosphate (GDP)-bound state; such a toggle is regulated by GEFs (guanine
nucleotide exchange factors) and GAPs (GTPase activating proteins). Here we
dissect a network motif between monomeric (m) and trimeric (t) GTPases
assembled exclusively in eukaryotic cells of multicellular organisms. To this
end, we develop a system of ordinary differential equations in which these two
classes of GTPases are interlinked conditional to their ON/OFF states within a
motif through feedforward and feedback loops. We provide formulas for the
steady states of the system and perform local stability analysis to investigate
the role of the different connections between the GTPase switches. A
feedforward from the active mGTPase to the GEF of the tGTPase was sufficient to
provide two locally stable states: one where both active/inactive forms of the
mGTPase can be interpreted as having low concentrations and the other where
both m- and tGTPase have high concentrations. When a feedback loop from the GEF
of the tGTPase to the GAP of the mGTPase was added to the feedforward system,
two other locally stable states emerged, both having the tGTPase inactivated
and being interpreted as having low active tGTPase concentrations. Finally, the
addition of a second feedback loop, from the active tGTPase to the GAP of the
mGTPase, gives rise to a family of steady states parametrized by the inactive
tGTPase concentrations. Our findings reveal that the coupling of these two
different GTPase motifs can dramatically change their steady state behaviors
and shed light on how such coupling may impact information processing in
eukaryotic cells.


  An interesting inference drawn by some Covid-19 epidemiological models is
that there exists a proportion of the population who are not susceptible to
infection -- even at the start of the current pandemic. This paper introduces a
model of the immune response to a virus. This is based upon the same sort of
mean-field dynamics as used in epidemiology. However, in place of the location,
clinical status, and other attributes of people in an epidemiological model, we
consider the state of a virus, B and T-lymphocytes, and the antibodies they
generate. Our aim is to formalise some key hypotheses as to the mechanism of
resistance. We present a series of simple simulations illustrating changes to
the dynamics of the immune response under these hypotheses. These include
attenuated viral cell entry, pre-existing cross-reactive humoral
(antibody-mediated) immunity, and enhanced T-cell dependent immunity. Finally,
we illustrate the potential application of this sort of model by illustrating
variational inversion (using simulated data) of this model to illustrate its
use in testing hypotheses. In principle, this furnishes a fast and efficient
immunological assay--based on sequential serology--that provides a (i)
quantitative measure of latent immunological responses and (ii) a Bayes optimal
classification of the different kinds of immunological response (c.f., glucose
tolerance tests used to test for insulin resistance). This may be especially
useful in assessing SARS-CoV-2 vaccines.


  This work is related to the setup of overflowing exponential fed-batch
cultures (O-EFBC) derived from carbon limited EFBC dedicated to the production
of mycosubtilin, an antifungal lipopeptide belonging to the iturin family.
O-EFBC permits the continuous removal of the product from the bioreactor
achieving a complete extraction of mycosubtilin. This paper also provides a
dynamical Monod-based growth model of this process that is accurate enough to
simulate the evolution of the specific growth rate and to correlate it to the
mycosubtilin specific productivity. Two particular and dependant phenomena
related to the foam overflow are taken into account by the model: the outgoing
flow rate of a broth volume and the loss of biomass. Interestingly, the biomass
concentration in the foam was found to be lower than the biomass concentration
in the bioreactor relating this process to a recycling one. Parameters of this
model are the growth yield on substrate and the maximal specific growth rate
estimated from experiments led at feed rates of 0.062, 0.071 and 0.086 h --1.
The model was extrapolated to five additional experiments carried out at feed
rates of 0.008, 0.022, 0.040, 0.042 and 0.062 h --1 enabling the correlation of
the mean specific growth rates with productivity results. Finally, a feed rate
of 0.086 h --1 corresponding to a mean specific growth rate of 0.070 h --1
allowed a specific productivity of 1.27 mg of mycosubtilin g --1 of dried
biomass h --1 .


  Secondary structure plays an important role in determining the function of
non-coding RNAs. Hence, identifying RNA secondary structures is of great value
to research. Computational prediction is a mainstream approach for predicting
RNA secondary structure. Unfortunately, even though new methods have been
proposed over the past 40 years, the performance of computational prediction
methods has stagnated in the last decade. Recently, with the increasing
availability of RNA structure data, new methods based on machine-learning
technologies, especially deep learning, have alleviated the issue. In this
review, we provide a comprehensive overview of RNA secondary structure
prediction methods based on machine-learning technologies and a tabularized
summary of the most important methods in this field. The current pending issues
in the field of RNA secondary structure prediction and future trends are also
discussed.


  Accurately determining a change in protein binding affinity upon mutations is
important for the discovery and design of novel therapeutics and to assist
mutagenesis studies. Determination of change in binding affinity upon mutations
requires sophisticated, expensive, and time-consuming wet-lab experiments that
can be aided with computational methods. Most of the computational prediction
techniques require protein structures that limit their applicability to protein
complexes with known structures. In this work, we explore the sequence-based
prediction of change in protein binding affinity upon mutation. We have used
protein sequence information instead of protein structures along with machine
learning techniques to accurately predict the change in protein binding
affinity upon mutation. Our proposed sequence-based novel change in protein
binding affinity predictor called PANDA gives better accuracy than existing
methods over the same validation set as well as on an external independent test
dataset. On an external test dataset, our proposed method gives a maximum
Pearson correlation coefficient of 0.52 in comparison to the state-of-the-art
existing protein structure-based method called MutaBind which gives a maximum
Pearson correlation coefficient of 0.59. Our proposed protein sequence-based
method, to predict a change in binding affinity upon mutations, has wide
applicability and comparable performance in comparison to existing protein
structure-based methods. A cloud-based webserver implementation of PANDA and
its python code is available at
https://sites.google.com/view/wajidarshad/software and
https://github.com/wajidarshad/panda.


  We study mechanisms for reopening economic activities that explore the trade
off between containing the spread of COVID-19 and maximizing economic impact.
This is of current importance as many organizations, cities, and states are
formulating reopening strategies. Our mechanisms, referred to as group
scheduling, are based on partitioning the population into groups and scheduling
each group on appropriate days with possible gaps (when all are quarantined).
Each group interacts with no other group and, importantly, any person who is
symptomatic in a group is quarantined. Specifically, our mechanisms are
characterized by three parameters $(g,d,t)$, where $g$ is the number of groups,
$d$ is the number of days a group is continuously scheduled, and $t$ is the gap
between cycles. We show that our mechanisms effectively trade off economic
activity for more effective control of the COVID-19 virus. In particular, we
show that the $(2,5,0)$ mechanism, which partitions the population into two
groups that alternatively work for five days each, flat lines the number of
COVID-19 cases quite effectively, while still maintaining economic activity at
70% of pre-COVID-19 level. We also study mechanisms such as $(2,3,2)$ and
$(3,3,0)$ that achieve a somewhat lower economic output (about 50%) at the cost
of more aggressive control of the virus; these could be applicable in
situations when the disease spread is more rampant in the population. We
demonstrate the efficacy of our mechanisms by theoretical analysis and
extensive experimental simulations on various epidemiological models. Our
mechanisms prove beneficial just by regulating human interactions. Moreover,
our results show that if the disease transmission (reproductive) rate is made
lower by following social distancing, mask wearing, and other public health
guidelines, it can further increase the efficacy of our mechanisms.


  We investigate how the initial number of infected individuals affects the
behavior of the critical susceptible-infected-recovered process. We analyze the
outbreak size distribution, duration of the outbreaks, and the role of
fluctuations.


  Susceptibility governs the dynamics of contagion. The classical SIR model is
one of the simplest compartmental models of contagion spread, assuming a single
shared susceptibility level. However, variation in susceptibility over a
population can fundamentally affect the dynamics of contagion and thus the
ultimate outcome of a pandemic. We develop mathematical machinery which
explicitly considers susceptibility variation, illuminates how the
susceptibility distribution is sculpted by contagion, and thence how such
variation affects the SIR differential questions that govern contagion. Our
methods allow us to derive closed form expressions for herd immunity thresholds
as a function of initial susceptibility distributions and suggests an
intuitively satisfying approach to inoculation when only a fraction of the
population is accessible to such intervention. Of particular interest, if we
assume static susceptibility of individuals in the susceptible pool, ignoring
susceptibility diversity {\em always} results in overestimation of the herd
immunity threshold and that difference can be dramatic. Therefore, we should
develop robust measures of susceptibility variation as part of public health
strategies for handling pandemics.


  The basis of several recent methods for drug repurposing is the key principle
that an efficacious drug will reverse the disease molecular 'signature' with
minimal side-effects. This principle was defined and popularized by the
influential 'connectivity map' study in 2006 regarding reversal relationships
between disease- and drug-induced gene expression profiles, quantified by a
disease-drug 'connectivity score.' Over the past 15 years, several studies have
proposed variations in calculating connectivity scores towards improving
accuracy and robustness in light of massive growth in reference drug profiles.
However, these variations have been formulated inconsistently using various
notations and terminologies even though they are based on a common set of
conceptual and statistical ideas. Therefore, we present a systematic
reconciliation of multiple disease-drug similarity metrics (ES, css, Sum,
Cosine, XSum, XCor, XSpe, XCos, EWCos) and connectivity scores (CS, RGES, NCS,
WCS, Tau, CSS, EMUDRA) by defining them using consistent notation and
terminology. In addition to providing clarity and deeper insights, this
coherent definition of connectivity scores and their relationships provides a
unified scheme that newer methods can adopt, enabling the computational
drug-development community to compare and investigate different approaches
easily. To facilitate the continuous and transparent integration of newer
methods, this article will be available as a live document
(https://jravilab.github.io/connectivity_scores) coupled with a GitHub
repository (https://github.com/jravilab/connectivity_scores) that any
researcher can build on and push changes to.


  Consider a population evolving from year to year through three seasons:
spring, summer and winter. Every spring starts with $N$ dormant individuals
waking up independently of each other according to a given distribution. Once
an individual is awake, it starts reproducing at a constant rate. By the end of
spring, all individuals are awake and continue reproducing independently as
Yule processes during the whole summer. In the winter, $N$ individuals chosen
uniformly at random go to sleep until the next spring, and the other
individuals die. We show that because an individual that wakes up unusually
early can have a large number of surviving descendants, for some choices of
model parameters the genealogy of the population will be described by a
$\Lambda$-coalescent. In particular, the beta coalescent can describe the
genealogy when the rate at which individuals wake up increases exponentially
over time. We also characterize the set of all $\Lambda$-coalescents that can
arise in this framework.


  Contact-tracing is an essential tool in order to mitigate the impact of
pandemic such as the COVID-19. In order to achieve efficient and scalable
contact-tracing in real time, digital devices can play an important role. While
a lot of attention has been paid to analyzing the privacy and ethical risks of
the associated mobile applications, so far much less research has been devoted
to optimizing their performance and assessing their impact on the mitigation of
the epidemic. We develop Bayesian inference methods to estimate the risk that
an individual is infected. This inference is based on the list of his recent
contacts and their own risk levels, as well as personal information such as
results of tests or presence of syndromes. We propose to use probabilistic risk
estimation in order to optimize testing and quarantining strategies for the
control of an epidemic. Our results show that in some range of epidemic
spreading (typically when the manual tracing of all contacts of infected people
becomes practically impossible, but before the fraction of infected people
reaches the scale where a lock-down becomes unavoidable), this inference of
individuals at risk could be an efficient way to mitigate the epidemic. Our
approaches translate into fully distributed algorithms that only require
communication between individuals who have recently been in contact. Such
communication may be encrypted and anonymized and thus compatible with privacy
preserving standards. We conclude that probabilistic risk estimation is capable
to enhance performance of digital contact tracing and should be considered in
the currently developed mobile applications.


  The spread of pandemics such as COVID-19 is strongly linked to human
activities. The objective of this paper is to specify and examine early
indicators of disease spread risk in cities during the initial stages of
outbreak based on patterns of human activities obtained from digital trace
data. In this study, the Venables distance (D_v), and the activity density
(D_a) are used to quantify and evaluate human activities for 193 US counties,
whose cumulative number of confirmed cases was greater than 100 as of March 31,
2020. Venables distance provides a measure of the agglomeration of the level of
human activities based on the average distance of human activities across a
city or a county (less distance could lead to a greater contact risk). Activity
density provides a measure of level of overall activity level in a county or a
city (more activity could lead to a greater risk). Accordingly, Pearson
correlation analysis is used to examine the relationship between the two human
activity indicators and the basic reproduction number in the following weeks.
The results show statistically significant correlations between the indicators
of human activities and the basic reproduction number in all counties, as well
as a significant leader-follower relationship (time lag) between them. The
results also show one to two weeks' lag between the change in activity
indicators and the decrease in the basic reproduction number. This result
implies that the human activity indicators provide effective early indicators
for the spread risk of the pandemic during the early stages of the outbreak.
Hence, the results could be used by the authorities to proactively assess the
risk of disease spread by monitoring the daily Venables distance and activity
density in a proactive manner.


  The emerging field of hybrid DNA - protein nanotechnology brings with it the
potential for many novel materials which combine the addressability of DNA
nanotechnology with versatility of protein interactions. However, the design
and computational study of these hybrid structures is difficult due to the
system sizes involved. To aid in the design and in silico analysis process, we
introduce here a coarse-grained DNA/RNA-protein model that extends the
oxDNA/oxRNA models of DNA/RNA with a coarse-grained model of proteins based on
an anisotropic network model representation. Fully equipped with analysis
scripts and visualization, our model aims to facilitate hybrid nanomaterial
design towards eventual experimental realization, as well as enabling study of
biological complexes. We further demonstrate its usage by simulating
DNA-protein nanocage, DNA wrapped around histones, and a nascent RNA in
polymerase.


  This paper shows that the generalized logistic distribution model is derived
from the well-known compartment model, consisting of susceptible, infected and
recovered compartments, abbreviated as the SIR model, under certain conditions.
In the SIR model, there are uncertainties in predicting the final values for
the number of infected population and the infectious parameter. However, by
utilizing the information obtained from the generalized logistic distribution
model, we can perform the SIR numerical computation more stably and more
accurately. Applications to severe acute respiratory syndrome (SARS) and
Coronavirus disease 2019 (COVID-19) using this combined method are also
introduced.


  It is well known that the classical recombination equation for two parent
individuals is equivalent to the law of mass action of a strongly reversible
chemical reaction network, and can thus be reformulated as a generalised
gradient system. Here, this is generalised to the case of an arbitrary number
of parents. Furthermore, the gradient structure of the backward-time
partitioning process is investigated.


  Lotka's seminal work (A.J. Lotka A., Proc. Natl. Acad. Sci. U.S.A. 6 (1920)
410) "on certain rhythmic relations'' is already one hundred years old, but the
research activity about pattern formations due to cyclical dominance is more
vibrant than ever. It is because non-transitive interactions have paramount
role on maintaining biodiversity and adequate human intervention into
ecological systems requires deeper understanding of related dynamical
processes. In this perspective article we overview different aspects of
biodiversity, with focus on how it can be maintained based on mathematical
modeling of last years. We also briefly discuss the potential links to
evolutionary game models of social systems, and finally, give an overview about
potential prospects for future research.


  Interactions between COVID-19 and other pathogens may change their dynamics.
Specifically, this may hinder the modelling of empirical data when the symptoms
of both infections are hard to distinguish. We introduce a model coupling the
dynamics of COVID-19 and seasonal influenza, simulating cooperation,
competition and asymmetric interactions. We find that the coupling synchronizes
both infections, with a strong influence on the dynamics of influenza, reducing
its time extent to a half.


  Single-molecule F\"orster resonance energy transfer (smFRET) of molecular
motors provides transformative insights into their dynamics and conformational
changes both at high temporal and spatial resolution simultaneously. However, a
key challenge of such FRET investigations is to observe a molecule in action
for long enough without restricting its natural function. The Anti-Brownian
ELectrokinetic Trap (ABEL trap) sets out to combine smFRET with molecular
confinement to enable observation times of up to several seconds while removing
any requirement of tethered surface attachment of the molecule in question. In
addition, the ABEL trap's inherent ability to selectively capture FRET active
molecules accelerates the data acquisition process. Here we exemplify the
capabilities of the ABEL trap in performing extended timescale smFRET
measurements on the molecular motor Rep, which is crucial for removing protein
blocks ahead of the advancing DNA replication machinery and for restarting
stalled DNA replication. We are able to monitor single Rep molecules up to 6 s
with 1 ms time resolution capturing multiple conformational switching events
during the observation time. Here we provide a step-by-step guide for the
rational design, construction and implementation of the ABEL trap for smFRET
detection of Rep in vitro. We include details of how to model the electric
potential at the trap site and use Hidden Markov analysis of the smFRET
trajectories.


  COVID-19 pandemic represents an unprecedented global health crisis in the
last 100 years. Its economic, social and health impact continues to grow and is
likely to end up as one of the worst global disasters since the 1918 pandemic
and the World Wars. Mathematical models have played an important role in the
ongoing crisis; they have been used to inform public policies and have been
instrumental in many of the social distancing measures that were instituted
worldwide.
  In this article we review some of the important mathematical models used to
support the ongoing planning and response efforts. These models differ in their
use, their mathematical form and their scope.


  Some of the key questions of interest during the COVID-19 pandemic (and all
outbreaks) include: where did the disease start, how is it spreading, who is at
risk, and how to control the spread. There are a large number of complex
factors driving the spread of pandemics, and, as a result, multiple modeling
techniques play an increasingly important role in shaping public policy and
decision making. As different countries and regions go through phases of the
pandemic, the questions and data availability also changes. Especially of
interest is aligning model development and data collection to support response
efforts at each stage of the pandemic. The COVID-19 pandemic has been
unprecedented in terms of real-time collection and dissemination of a number of
diverse datasets, ranging from disease outcomes, to mobility, behaviors, and
socio-economic factors. The data sets have been critical from the perspective
of disease modeling and analytics to support policymakers in real-time. In this
overview article, we survey the data landscape around COVID-19, with a focus on
how such datasets have aided modeling and response through different stages so
far in the pandemic. We also discuss some of the current challenges and the
needs that will arise as we plan our way out of the pandemic.


  The COVID-19 pandemic has resulted in more than 30.35 million infections and
9, 50, 625 deaths in 212 countries over the last few months. Different drug
intervention acting at multiple stages of pathogenesis of COVID-19 can
substantially reduce the infection induced mortality. The current within-host
mathematical modeling studies deals with the optimal drug regimen and the
efficacy of combined therapy in treatment of COVID-19. The drugs/interventions
considered include Arbidol, Remdesivir, Inteferon (INF) and
Lopinavir/Ritonavir. It is concluded that these drug interventions when
administered individually or in combination reduce the infected cells and viral
load. Four scenarios involving administration of single drug intervention, two
drug interventions, three drug interventions and all the four have been
discussed. In all these scenarios the optimal drug regimen is proposed based on
two methods. In the first method these medical interventions are modeled as
control interventions and a corresponding objective function and optimal
control problem is formulated. In this setting the optimal drug regimen is
proposed. Later using the the comparative effectiveness method the optimal drug
regimen is proposed based on basic reproduction number and viral load. The
average infected cell count and viral load decreased the most when all the four
interventions were applied together. On the other hand the average susceptible
cell count decreased the best when Arbidol alone was administered. The basic
reproduction number and viral count decreased the best when all the four
interventions were applied together reinstating the fact obtained earlier in
the optimal control setting. These findings may help physicians with decision
making in treatment of life-threatening COVID-19 pneumonia.


  Bacillus subtilis ATCC6633 produces mycosubtilin, a non-ribosomally
synthesized lipopeptide of the iturin family which presents antagonistic
activities toward various phytopathogens. Different homologues with fatty acid
moiety varying from C 15 to C 17 are usually co-produced, with their biological
activities increasing with the number of carbons in the fatty acid chain. In
the present report, we highlight that growth temperature modulates both the
extent of mycosubtilin production and the relative abundance of the different
homologues. A 30-fold increase in mycosubtilin production was observed when the
temperature was decreased from 37 C to 25 C for both strain ATCC6633 and its
derivative BBG100, a constitutive mycosubtilin overproducer. However, no
significant difference in either the expression of the mycosubtilin synthetase
encoding genes or in the intracellular synthetase concentration could be found,
suggesting that the observed phenotype originated from a higher mycosubtilin
synthetase turnover at lower temperature. We also point out that lower growth
temperature leads to an increased proportion of odd-numbered fatty acid
homologues as a consequence of de novo synthesis of C 17 anteiso fatty acid
following cell adaptation to low temperatures.


  A Bacillus subtilis derivative was obtained from strain ATCC 6633 by
replacement of the native promoter of the mycosubtilin operon by a constitutive
promoter originating from the replication gene repU of the Staphylococcus
aureus plasmid pUB110. The recombinant strain, designated BBG100, produced up
to 15-fold more mycosubtilin than the wild type produced. The overproducing
phenotype was related to enhancement of the antagonistic activities against
several yeasts and pathogenic fungi. Hemolytic activities were also clearly
increased in the modified strain. Mass spectrometry analyses of enriched
mycosubtilin extracts showed similar patterns of lipopeptides for BBG100 and
the wild type. Interestingly, these analyses also revealed a new form of
mycosubtilin which was more easily detected in the BBG100 sample. When tested
for its biocontrol potential, wild-type strain ATCC 6633 was almost ineffective
for reducing a Pythium infection of tomato seedlings. However, treatment of
seeds with the BBG100 overproducing strain resulted in a marked increase in the
germination rate of seeds. This protective effect afforded by mycosubtilin
overproduction was also visualized by the significantly greater fresh weight of
emerging seedlings treated with BBG100 compared to controls or seedlings
inoculated with the wild-type strain.


  Genomes may be analyzed from an information viewpoint as very long strings,
containing functional elements of variable length, which have been assembled by
evolution. In this work an innovative information theory based algorithm is
proposed, to extract significant (relatively small) dictionaries of genomic
words. Namely, conceptual analyses are here combined with empirical studies, to
open up a methodology for the extraction of variable length dictionaries from
genomic sequences, based on the information content of some factors. Its
application to human chromosomes highlights an original inter-chromosomal
similarity in terms of factor distributions.


  We analyse the dynamics of fishing vessels with different home ports in an
area where these vessels, in choosing where to fish, are influenced by their
own experience in the past and by their current observation of the locations of
other vessels in the fleet. Empirical data from the boats near Ancona and
Pescara shows stylized statistical properties that are reminiscent of Kirman
and F\"ollmer's ant recruitment model, although with two ant colonies
represented by the two ports. From the point of view of a fisherman, the two
fishing areas are not equally attractive, and he tends to prefer the one closer
to where he is based. This piece of evidence led us to extend the original ants
model to a situation with two asymmetric zones and finite resources. We show
that, in the mean-field regime, our model exhibits the same properties as the
empirical data. We obtain a phase diagram that separates high and low herding
regimes, but also fish population extinction. Our analysis has interesting
policy implications for the ecology of fishing areas. It also suggests that
herding behaviour here, just as in financial markets, will lead to significant
fluctuations in the amount of fish landed, as the boat concentration on one
area at a given point in time will diminish the overall catch, such loss not
being compensated by the reproduction of fish in the other area. In other
terms, individually rational behaviour will not lead to collectively optimal
results.


  Motivated by recent experiments on an antibiotic resistance gene, we
investigate genetic interactions between synonymous mutations in the framework
of exclusion models of translation. We show that the range of possible
interactions is markedly different depending on whether translation efficiency
is assumed to be proportional to ribosome current or ribosome speed. In the
first case every mutational effect has a definite sign that is independent of
genetic background, whereas in the second case the effect-sign can vary
depending on the presence of other mutations. The latter result is demonstrated
using configurations of multiple translational bottlenecks induced by slow
codons.


  Resources are rarely distributed uniformly within a population. Heterogeneity
in the concentration of a drug, the quality of breeding sites, or wealth can
all affect evolutionary dynamics. In this study, we represent a collection of
properties affecting the fitness at a given location using a color. A green
node is rich in resources while a red node is poorer. More colors can represent
a broader spectrum of resource qualities. For a population evolving according
to the birth-death Moran model, the first question we address is which
structures, identified by graph connectivity and graph coloring, are
evolutionarily equivalent. We prove that all properly two-colored, undirected,
regular graphs are evolutionarily equivalent (where "properly colored" means
that no two neighbors have the same color). We then compare the effects of
background heterogeneity on properly two-colored graphs to those with
alternative schemes in which the colors are permuted. Finally, we discuss
dynamic coloring as a model for spatiotemporal resource fluctuations, and we
illustrate that random dynamic colorings often diminish the effects of
background heterogeneity relative to a proper two-coloring.


  We introduce the notion of non-oscillation, propose a constructive method for
its robust verification, and study its application to biological interaction
networks (also known as, chemical reaction networks). We begin by revisiting
Muldowney's result on non-existence of periodic solutions based on the study of
the variational system of the second additive compound of the Jacobian of a
nonlinear system. We show that exponential stability of the latter rules out
limit cycles, quasi-periodic solutions, and broad classes of oscillatory
behavior. We focus then on nonlinear equations arising in biological
interaction networks with general kinetics, and we show that the dynamics of
the aforementioned variational system can be embedded in a linear differential
inclusion. We then propose algorithms for constructing piecewise linear
Lyapunov functions to certify global robust non-oscillatory behavior. Finally,
we apply our techniques to study several regulated enzymatic cycles where
available methods are not able to provide any information about their
qualitative global behavior.


  Recent work has proven the existence of extreme inbreeding in a European
ancestry sample taken from the contemporary UK population \cite{nature_01}.
This result brings our attention again to a math problem related to inbreeding
family trees and diversity. Groups with a finite number of individuals could
give a variety of genetic relationships. { In previous works
\cite{PhysRevE.92.052132, PhysRevE.90.022125, JARNE20191}, we have addressed
the issue of building inbreeding trees for biparental reproduction using
Markovian models. Here, we extend these studies by presenting an algorithm to
generate and represent inbreeding trees with no overlapping generations. We
explicitly assume a two-gender reproductory scheme, and we pay particular
attention to the links between nodes. We show that even for a simple case with
a relatively small number of nodes in the tree, there are a large number of
possible ways to rearrange the links between generations. We present an
open-source python code to generate the tree graph, the adjacency matrix, and
the histogram of the links for each different tree representation. We show how
this mapping reflects the difference between tree realizations, and how
valuable information may be extracted upon inspection of these matrices. The
algorithm includes a feature to average several tree realizations, obtain the
connectivity distribution, and calculate the average and mean value. We used
this feature to compare trees with a different number of generations and nodes.
The code presented here, available in Git-Hub, may be easily modified to be
applied to other areas of interest involving connections between individuals,
extend the study to add more characteristics of the different nodes, etc.


  We analyse the stability of linear dynamical systems defined on sparse,
random graphs with predator-prey, competitive, and mutualistic interactions.
These systems are aimed at modelling the stability of fixed points in large
systems defined on complex networks, such as, ecosystems consisting of a large
number of species that interact through a food-web. We develop an exact theory
for the spectral distribution and the leading eigenvalue of the corresponding
sparse Jacobian matrices. This theory reveals that the nature of local
interactions have a strong influence on system's stability. We show that, in
general, linear dynamical systems defined on random graphs with a prescribed
degree distribution of unbounded support are unstable if they are large enough,
implying a tradeoff between stability and diversity. Remarkably, in contrast to
the generic case, antagonistic systems that only contain interactions of the
predator-prey type can be stable in the infinite size limit. This qualitatively
feature for antagonistic systems is accompanied by a peculiar oscillatory
behaviour of the dynamical response of the system after a perturbation, when
the mean degree of the graph is small enough. Moreover, for antagonistic
systems we also find that there exist a dynamical phase transition and critical
mean degree above which the response becomes non-oscillatory.


  A theory on the conformation transition for SARS-CoV-2 spike protein (S) is
proposed. The conformation equilibrium between open (up) and closed (down)
conformations of receptor binding domain (RBD) of the spike is studied from the
first-principle. The conformational state population is deduced from the free
energy change in conformation transition of S protein. We demonstrated that the
free energy includes two parts, one from the multi-minima of conformational
potential and another from the variation of structural elasticity. Both factors
are dependent of amino acid mutation. The former is related to the change of
affinity of RBD to ACE 2 due to the mutation in the subdomain RBM (receptor
binding motif) of RBD. The latter is caused by the change of elastic energy of
S protein. When the affinity has increased significantly and/or the elastic
energy has been reduced substantially the equilibrium is biased to the open
conformation. Only then can the virus infection process continue. Possible new
SARS-CoV-2 variants from amino acid mutations in 5-9 sites on RBD interface are
predicted. The elastic energy variation needed for conformation transition is
estimated quantitatively. Taking the elastic-structural change into account
more virus variants are possible. Both the D614G mutation, the K986P mutation
and the new variants 501Y in current SARS-CoV-2 pandemic can be interpreted
from the presented theory. The comparison of the infectivity of SARS-CoV-2 with
SARS-CoV-1 is made from the point of conformation equilibrium. Why the virus
entrance takes priority in lower temperature and higher humidity is interpreted
by the present theory. The conformational transition influenced by
electromagnetic field is also discussed briefly.


  Objective: Ghrelin is an orexigenic peptide hormone involved in the
regulation of energy homeostasis, food intake and glucose metabolism. Serum
levels increase anticipating a meal and fall afterwards. Underlying genetic
mechanisms of the ghrelin secretion are unknown. Methods: Total serum ghrelin
was measured in 1501 subjects selected from the population-based
LIFE-ADULT-sample after an overnight fast. A genome-wide association study
(GWAS) was performed. Gene-based expression association analyses
(transcriptome-wide association study (TWAS)) were done using MetaXcan.
Results: In the GWAS, three loci reached genome-wide significance: the
WW-domain containing the oxidoreductase-gene (WWOX; p=1.80E-10) on chromosome
16q23.3-24.1 (SNP: rs76823993); the Contactin-Associated Protein-Like 2 gene
(CNTNAP2; p=9.0E-9) on chromosome 7q35-q36 (SNP: rs192092592) and the Ghrelin
And Obestatin Prepropeptide gene (GHRL; p=2.72E-8) on chromosome 3p25.3 (SNP:
rs143729751). In the TWAS, serum ghrelin was negatively associated with RNA
expression of the GDNF Family Receptor Alpha Like (GFRAL), receptor of the
anorexigenic Growth Differentiation Factor-15 (GDF15), (z-score=-4.288,
p=1.81E-05). Furthermore, ghrelin was positively associated with Ribosomal
Protein L36 (RPL36; z-score=4.848, p=1.25E-06). Conclusions: Our findings
provide evidence of a functional link between two major players of weight
regulation, the ghrelin system and the GDF15/GFRAL-pathway.


  The COVID-19 pandemic has resulted in more than 14.5 million infections and
6,04,917 deaths in 212 countries over the last few months. Different drug
intervention acting at multiple stages of pathogenesis of COVID-19 can
substantially reduce the infection induced,thereby decreasing the mortality.
Also population level control strategies can reduce the spread of the COVID-19
substantially. Motivated by these observations, in this work we propose and
study a multi scale model linking both within-host and between-host dynamics of
COVID-19. Initially the natural history dealing with the disease dynamics is
studied. Later, comparative effectiveness is performed to understand the
efficacy of both the within-host and population level interventions. Findings
of this study suggest that a combined strategy involving treatment with drugs
such as Arbidol, remdesivir, Lopinavir/Ritonavir that inhibits viral
replication and immunotherapies like monoclonal antibodies, along with
environmental hygiene and generalized social distancing proved to be the best
and optimal in reducing the basic reproduction number and environmental spread
of the virus at the population level.


  In phylogenetics it is of interest for rate matrix sets to satisfy closure
under matrix multiplication as this makes finding the set of corresponding
transition matrices possible without having to compute matrix exponentials. It
is also advantageous to have a small number of free parameters as this, in
applications, will result in a reduction of computation time. We explore a
method of building a rate matrix set from a rooted tree structure by assigning
rates to internal tree nodes and states to the leaves, then defining the rate
of change between two states as the rate assigned to the most recent common
ancestor of those two states. We investigate the properties of these matrix
sets from both a linear algebra and a graph theory perspective and show that
any rate matrix set generated this way is closed under matrix multiplication.
The consequences of setting two rates assigned to internal tree nodes to be
equal are then considered. This methodology could be used to develop
parameterised models of amino acid substitution which have a small number of
parameters but convey biological meaning.


  The dense social contact networks and high mobility in congested urban areas
facilitate the rapid transmission of infectious diseases. Typical mechanistic
epidemiological models are either based on uniform mixing with ad-hoc contact
processes or need real-time or archived population mobility data to simulate
the social networks. However, the rapid and global transmission of the novel
coronavirus (SARS-CoV-2) has led to unprecedented lockdowns at global and
regional scales, leaving the archived datasets to limited use. While it is
often hypothesized that population density is a significant driver in disease
propagation, the disparate disease trajectories and infection rates exhibited
by the different cities with comparable densities require a high-resolution
description of the disease and its drivers. In this study, we explore the
impact of the creation of containment zones on travel patterns within the city.
Further, we use a dynamical network-based infectious disease model to
understand the key drivers of disease spread at sub-kilometer scales
demonstrated in the city of Ahmedabad, India, which has been classified as a
SARS-CoV-2 hotspot. We find that in addition to the contact network and
population density, road connectivity patterns and ease of transit are strongly
correlated with the rate of transmission of the disease. Given the limited
access to real-time traffic data during lockdowns, we generate road
connectivity networks using open-source imageries and travel patterns from
open-source surveys and government reports. Within the proposed framework, we
then analyze the relative merits of social distancing, enforced lockdowns, and
enhanced testing and quarantining mitigating the disease spread.


  We develop a novel hybrid epidemiological model and a specific methodology
for its calibration to distinguish and assess the impact of mobility
restrictions (given by Apple's mobility trends data) from other complementary
non-pharmaceutical interventions (NPIs) used to control the spread of COVID-19.
Using the calibrated model, we estimate that mobility restrictions contribute
to 47 % (US States) and 47 % (worldwide) of the overall suppression of the
disease transmission rate using data up to 13/08/2020. The forecast capacity of
our model was evaluated doing four-weeks ahead predictions. Using data up to
30/06/20 for calibration, the mean absolute percentage error (MAPE) of the
prediction of cumulative deceased individuals was 5.0 % for the United States
(51 states) and 6.7 % worldwide (49 countries). This MAPE was reduced to 3.5%
for the US and 3.8% worldwide using data up to 13/08/2020. We find that the
MAPE was higher for the total confirmed cases at 11.5% worldwide and 10.2% for
the US States using data up to 13/08/2020. Our calibrated model achieves an
average R-Squared value for cumulative confirmed and deceased cases of 0.992
using data up to 30/06/20 and 0.98 using data up to 13/08/20.


  We present a modified age-structured SIR model based on known patterns of
social contact and distancing measures within Washington, USA. We find that
population age-distribution has a significant effect on disease spread and
mortality rate, and contribute to the efficacy of age-specific contact and
treatment measures. We consider the effect of relaxing restrictions across less
vulnerable age-brackets, comparing results across selected groups of varying
population parameters. Moreover, we analyze the mitigating effects of
vaccinations and examine the effectiveness of age-targeted distributions.
Lastly, we explore how our model can applied to other states to reflect
social-distancing policy based on different parameters and metrics.


  The interaction between actin filaments and microtubules is crucial for many
eukaryotic cellular processes, such as, among others, cell polarization, cell
motility and cellular wound healing. The importance of this interaction has
long been recognised, yet very little is understood about both the underlying
mechanisms and the consequences for the spatial (re)organization of the
cellular cytoskeleton. At the same time, understanding the causes and the
consequences of the interaction between different biomolecular components are
key questions for \emph{in vitro} research involving reconstituted biomolecular
systems, especially in the light of current interest in creating minimal
synthetic cells. In this light, recent \emph{in vitro} experiments have shown
that the actin-microtubule interaction mediated by the cytolinker TipAct, which
binds to actin lattice and microtubule tip, causes the directed transport of
actin filaments. We develop an analytical theory of dynamically unstable
microtubules, nucleated from the center of a spherical cell, in interaction
with actin filaments. We show that, depending on the balance between the
diffusion of unbound actin filaments and propensity to bind microtubules, actin
is either concentrated in the center of the cell, where the density of
microtubules is highest, or becomes localized to the cell cortex.


  Drug discovery using deep learning has attracted a lot of attention of late
as it has obvious advantages like higher efficiency, less manual guessing and
faster process time. In this paper, we present a novel neural network for
generating small molecules similar to the ones in the training set. Our network
consists of an encoder made up of bi-GRU layers for converting the input
samples to a latent space, predictor for enhancing the capability of encoder
made up of 1D-CNN layers and a decoder comprised of uni-GRU layers for
reconstructing the samples from the latent space representation. Condition
vector in latent space is used for generating molecules with the desired
properties. We present the loss functions used for training our network,
experimental details and property prediction metrics. Our network outperforms
previous methods using Molecular weight, LogP and Quantitative Estimation of
Drug-likeness as the evaluation metrics.


  Tree-based phylogenetic networks, which may be roughly defined as
leaf-labeled networks built by adding arcs only between the original tree
edges, have elegant properties for modeling evolutionary histories. We answer
an open question of Francis, Semple, and Steel about the complexity of
determining how far a phylogenetic network is from being tree-based, including
non-binary phylogenetic networks. We show that finding a phylogenetic tree
covering the maximum number of nodes in a phylogenetic network can be be
computed in polynomial time via an encoding into a minimum-cost maximum flow
problem.


  Motivation: With the development of third-generation sequencing technologies,
people are able to obtain DNA sequences with lengths from 10s to 100s of kb.
These long reads allow protein domain annotation without assembly, thus can
produce important insights into the biological functions of the underlying
data. However, the high error rate in third-generation sequencing data raises a
new challenge to established domain analysis pipelines. The state-of-the-art
methods are not optimized for noisy reads and have shown unsatisfactory
accuracy of domain classification in third-generation sequencing data. New
computational methods are still needed to improve the performance of domain
prediction in long noisy reads. Results: In this work, we introduce ProDOMA, a
deep learning model that conducts domain classification for third-generation
sequencing reads. It uses deep neural networks with 3-frame translation
encoding to learn conserved features from partially correct translations. In
addition, we formulate our problem as an open-set problem and thus our model
can reject unrelated DNA reads such as those from noncoding regions. In the
experiments on simulated reads of protein coding sequences and real reads from
the human genome, our model outperforms HMMER and DeepFam on protein domain
classification. In summary, ProDOMA is a useful end-to-end protein domain
analysis tool for long noisy reads without relying on error correction.
Availability: The source code and the trained model are freely available at
https://github.com/strideradu/ProDOMA. Contact: yannisun@cityu.edu.hk


  Time delay has been incorporated in models to reflect certain physical or
biological meaning. The theory of delay differential equations (DDEs), which
has seen extensive growth in the last seventy years or so, can be used to
examine the effects of time delay in the dynamical behavior of systems being
considered. Numerical tools to study DDEs have played a significant role not
only in illustrating theoretical results but also in discovering interesting
dynamics of the model. DDE-Biftool, which is a Matlab package for numerical
continuation and numerical bifurcation analysis of DDEs, is one of the most
utilized and popular numerical tools for DDEs. In this paper, we present a
guide to using the latest version of DDE-Biftool targeted to researchers who
are new to the study of time delay systems. A short discussion of an example
application, which is a harvested predator-prey model with a single discrete
time delay, will be presented first. We then implement this example model in
DDE-Biftool, pointing out features where beginners need to be cautious. We end
with a comparison of our theoretical and numerical results.


  Rapid climate change has wide-ranging implications for the Arctic region,
including sea ice loss, increased geopolitical attention, and expanding
economic activity, including a dramatic increase in shipping activity. As a
result, the risk of harmful non-native marine species being introduced into
this critical region will increase unless policy and management steps are
implemented in response. Using big data about shipping, ecoregions, and
environmental conditions, we leverage network analysis and data mining
techniques to assess, visualize, and project ballast water-mediated species
introductions into the Arctic and dispersal of non-native species within the
Arctic. We first identify high-risk connections between the Arctic and
non-Arctic ports that could be sources of non-native species over 15 years
(1997-2012) and observe the emergence of shipping hubs in the Arctic where the
cumulative risk of non-native species introduction is increasing. We then
consider how environmental conditions can constrain this Arctic introduction
network for species with different physiological limits, thus providing a
species-level tool for decision-makers. Next, we focus on within-Arctic
ballast-mediated species dispersal where we use higher-order network analysis
to identify critical shipping routes that may facilitate species dispersal
within the Arctic. The risk assessment and projection framework we propose
could inform risk-based assessment and management of ship-borne invasive
species in the Arctic.


  Due to SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) being a
novel virus, there are currently no known effective antiviral drugs capable of
slowing its progress. To accelerate the discovery of potential drug candidates,
bioinformatics based in silico drug discovery can be applied as a very robust
tool. In the present study, more than 60 antiviral drugs already available on
the market, were chosen after literature survey. These can be used in clinical
trials for the treatment of COVID-19. In this study, these candidate drugs were
ranked based on their potential to interact with the Spike protein and RdRp
(RNA-dependent RNA polymerase) of SARS-CoV-2. Additionally, the mechanism of
their action as well as how the virus infection can utilize Hemoglobin to
decrease the oxygen level in blood is explained. Moreover, multiple sequence
alignments of the Spike protein with 75 sequences of different viruses from the
Orthocoronavirinae subfamily were performed. This gives insight into the
evolutionarily conserved domains that can be targeted using drug or antibody
treatment. This multidimensional study opens a new layer of understanding about
the most effective drug-targetable sites on the Spike protein of SARS-CoV-2.


  We present a highly efficient and accurate hybrid stochastic simulation
algorithm (HSSA) for the purpose of simulating a subset of biochemical
reactions of large gene regulatory networks (GRN). The algorithm relies on the
separability of a GRN into two groups of reactions, A and B, such that the
reactions in A can be simulated via a stochastic simulation algorithm (SSA),
while those in group B can yield to a deterministic description via ordinary
differential equations. First, we derive exact expressions needed to sample the
next reaction time and reaction type, and then give two examples of how a GRN
can be partitioned. Although the methods presented here can be applied to a
variety of different stochastic systems within GRN, we focus on simulating
mRNAs in particular. To demonstrate the accuracy and efficiency of this
algorithm, we apply it to a three-gene oscillator, first in one cell, and then
in an array of cells (up to 64 cells) interacting via molecular diffusion, and
compare its performance to the Gillespie algorithm (GA). Depending on the
particular numerical values of the system parameters, and the partitioning
itself, we show that our algorithm is between 11 and 445 times faster than the
GA.


  In the study of infectious diseases on networks, researchers calculate
epidemic thresholds to help forecast whether a disease will eventually infect a
large fraction of a population. Because network structure typically changes in
time, which fundamentally influences the dynamics of spreading processes on
them and in turn affects epidemic thresholds for disease propagation, it is
important to examine epidemic thresholds in temporal networks. Most existing
studies of epidemic thresholds in temporal networks have focused on models in
discrete time, but most real-world networked systems evolve continuously in
time. In our work, we encode the continuous time-dependence of networks into
the evaluation of the epidemic threshold of a
susceptible--infected--susceptible (SIS) process by studying an SIS model on
tie-decay networks. We derive the epidemic-threshold condition of this model,
and we perform numerical experiments to verify it. We also examine how
different factors---the decay coefficients of the tie strengths in a network,
the frequency of interactions between nodes, and the sparsity of the underlying
social network in which interactions occur---lead to decreases or increases of
the critical values of the threshold and hence contribute to facilitating or
impeding the spread of a disease. We thereby demonstrate how the features of
tie-decay networks alter the outcome of disease spread.


  We study how international flights can facilitate the spread of an epidemic
to a worldwide scale. We combine an infrastructure network of flight
connections with a population density dataset to derive the mobility network,
and then we define an epidemic framework to model the spread of the disease.
Our approach combines a compartmental SEIRS model with a graph diffusion model
to capture the clusteredness of the distribution of the population. The
resulting model is characterised by the dynamics of a metapopulation SEIRS,
with amplification or reduction of the infection rate which is determined also
by the mobility of individuals. We use simulations to characterise and study a
variety of realistic scenarios that resemble the recent spread of COVID-19.
Crucially, we define a formal framework that can be used to design epidemic
mitigation strategies: we propose an optimisation approach based on genetic
algorithms that can be used to identify an optimal airport closure strategy,
and that can be employed to aid decision making for the mitigation of the
epidemic, in a timely manner.


  Evolutionary game theory has proven to be an elegant framework providing many
fruitful insights in population dynamics and human behaviour. Here, we focus on
the aspect of behavioural plasticity and its effect on the evolution of
populations. We consider games with only two strategies in both well-mixed
infinite and finite populations settings. We assume that individuals might
exhibit behavioural plasticity referred to as incompetence of players. We study
the effect of such heterogeneity on the outcome of local interactions and,
ultimately, on global competition. For instance, a strategy that was dominated
before can become desirable from the selection perspective when behavioural
plasticity is taken into account. Furthermore, it can ease conditions for a
successful fixation in infinite populations' invasions. We demonstrate our
findings on the examples of Prisoners' Dilemma and Snowdrift game, where we
define conditions under which cooperation can be promoted.


  In this paper we consider the global qualitative properties of a
stochastically perturbed logistic model of population growth. In this model,
the stochastic perturbations are assumed to be of the white noise type and are
proportional to the current population size. Using the direct Lyapunov method,
we established the global properties of this stochastic differential equation.
In particular, we found that solutions of the equation oscillate around an
interval, and explicitly found the end points of this interval. Moreover, we
found that, if the magnitude of the noise exceeds a certain critical level
(which is also explicitly found), then the stochastic stabilisation
("stabilisation by noise") of the zero solution occurs. In this case, (i) the
origin is the lower boundary of the interval, and (ii) the extinction of the
population due to stochasticity occurs almost sure (a.s.) for a finite time.


  As per the records of theWorld Health Organization, the first formally
reported incidence of Zika virus occurred in Brazil in May 2015. The disease
then rapidly spread to other countries in Americas and East Asia, affecting
more than 1,000,000 people. Zika virus is primarily transmitted through bites
of infected mosquitoes of the species Aedes (Aedes aegypti and Aedes
albopictus). The abundance of mosquitoes and, as a result, the prevalence of
Zika virus infections are common in areas which have high precipitation, high
temperature, and high population density.Nonlinear spatio-temporal dependency
of such data and lack of historical public health records make prediction of
the virus spread particularly challenging. In this article, we enhance Zika
forecasting by introducing the concepts of topological data analysis and,
specifically, persistent homology of atmospheric variables, into the virus
spread modeling. The topological summaries allow for capturing higher order
dependencies among atmospheric variables that otherwise might be unassessable
via conventional spatio-temporal modeling approaches based on geographical
proximity assessed via Euclidean distance. We introduce a new concept of
cumulative Betti numbers and then integrate the cumulative Betti numbers as
topological descriptors into three predictive machine learning models: random
forest, generalized boosted regression, and deep neural network. Furthermore,
to better quantify for various sources of uncertainties, we combine the
resulting individual model forecasts into an ensemble of the Zika spread
predictions using Bayesian model averaging. The proposed methodology is
illustrated in application to forecasting of the Zika space-time spread in
Brazil in the year 2018.


  The rate of mRNA translation depends on the initiation, elongation, and
termination rates of ribosomes along the mRNA. These rates depend on many
"local" factors like the abundance of free ribosomes and tRNA molecules in the
vicinity of the mRNA molecule. All these factors are stochastic and their
experimental measurements are also noisy. An important question is how protein
production in the cell is affected by this considerable variability. We develop
a new theoretical framework for addressing this question by modeling the rates
as identically and independently distributed random variables and using tools
from random matrix theory to analyze the steady-state production rate. The
analysis reveals a principle of universality: the average protein production
rate depends only on the of the set of possible values that the random variable
may attain. This explains how total protein production can be stabilized
despite the overwhelming stochasticticity underlying cellular processes.


  This document compares SRI's Pathway Tools (PTools) software with the
Reactome software. Both software systems serve the pathway bioinformatics area,
including representation and analysis of metabolic pathways and signaling
pathways. The comparison covers pathway bioinformatics capabilities, but does
not cover other major facets of Pathway Tools that are completely absent from
the Reactome software: Pathway Tools genome-informatics capabilities,
regulatory informatics capabilities, and table-based analysis tools
(SmartTables).
  Our overall findings are as follows. (1) PTools is significantly ahead of
Reactome in its basic information pages. For example, PTools pathway layout
algorithms have been developed to an advanced state over several decades,
whereas Reactome pathway layouts are illegible, omit important information, and
are created manually and therefore cannot scale to thousands of genomes. (2)
PTools is far ahead of Reactome in omics analysis. PTools includes all of the
omics-analysis methods that Reactome provides, and includes multiple methods
that Reactome lacks. (3) PTools contains a metabolic route search tool
(searching for paths through the metabolic network), which Reactome lacks. (4)
PTools is significantly ahead of Reactome in inference of metabolic pathways
from genome information to create new metabolic databases. (5) PTools has an
extensive complement of metabolic-modeling tools whereas Reactome has none. (6)
PTools is more scalable than Reactome, handling 18,000 genomes versus 90
genomes for Reactome. (7) PTools has a larger user base than Reactome. PTools
powers 17 websites versus two for Reactome. PTools has been licensed by 10,800
users (Reactome licensed user count is unknown).


  Recreational fishing is a highly socio-ecological process. Although
recreational fisheries are self-regulating and resilient, changing
anthropogenic pressure drives these fisheries to overharvest and collapse.
Here, we evaluate the effect of demographic and environmental stochasticity for
a social-ecological two-species fish model. In the presence of noise, we find
that an increase in harvesting rate drives a critical transition from
high-yield/low-price fisheries to low-yield/high-price fisheries. To calculate
stochastic trajectories for demographic noise, we derive the master equation
corresponding to the model and perform Monte-Carlo simulation. Moreover, the
analysis of probabilistic potential and mean first-passage time reveals the
resilience of alternative steady states. We also describe the efficacy of a few
generic indicators in forecasting sudden transitions. Furthermore, we show that
incorporating social norms on the model allows moderate fish density to
maintain despite higher harvesting rates. Overall, our study highlights the
occurrence of critical transitions in a stochastic social-ecological model and
suggests ways to mitigate them.


  Proteins are translated from the N- to the C-terminus, raising the basic
question of how this innate directionality affects their evolution. To explore
this question, we analyze 16,200 structures from the protein data bank (PDB).
We find remarkable enrichment of $\alpha$-helices at the C terminus and
$\beta$-strands at the N terminus. Furthermore, this $\alpha$-$\beta$ asymmetry
correlates with sequence length and contact order, both determinants of folding
rate, hinting at possible links to co-translational folding (CTF). Hence, we
propose the 'slowest-first' scheme, whereby protein sequences evolved
structural asymmetry to accelerate CTF: the slowest of the
cooperatively-folding segments are positioned near the N terminus so they have
more time to fold during translation. A phenomenological model predicts that
CTF can be accelerated by asymmetry, up to double the rate, when folding time
is commensurate with translation time; analysis of the PDB reveals that
structural asymmetry is indeed maximal in this regime. This correspondence is
greater in prokaryotes, which generally require faster protein production.
Altogether, this indicates that accelerating CTF is a substantial evolutionary
force whose interplay with stability and functionality is encoded in sequence
asymmetry.


  In order to design a more potent and effective chemical entity, it is
essential to identify molecular structures with the desired chemical
properties. Recent advances in generative models using neural networks and
machine learning are being widely used by many emerging startups and
researchers in this domain to design virtual libraries of drug-like compounds.
Although these models can help a scientist to produce novel molecular
structures rapidly, the challenge still exists in the intelligent exploration
of the latent spaces of generative models, thereby reducing the randomness in
the generative procedure. In this work we present a manifold traversal with
heuristic search to explore the latent chemical space. Different heuristics and
scores such as the Tanimoto coefficient, synthetic accessibility, binding
activity, and QED drug-likeness can be incorporated to increase the validity
and proximity for desired molecular properties of the generated molecules. For
evaluating the manifold traversal exploration, we produce the latent chemical
space using various generative models such as grammar variational autoencoders
(with and without attention) as they deal with the randomized generation and
validity of compounds. With this novel traversal method, we are able to find
more unseen compounds and more specific regions to mine in the latent space.
Finally, these components are brought together in a simple platform allowing
users to perform search, visualization and selection of novel generated
compounds.


  Mass spectrometry (MS) based omics data analysis require significant time and
resources. To date, few parallel algorithms have been proposed for deducing
peptides from mass spectrometry-based data. However, these parallel algorithms
were designed, and developed when the amount of data that needed to be
processed was smaller in scale. In this paper, we prove that the communication
bound that is reached by the \emph{existing} parallel algorithms is
$\Omega(mn+2r\frac{q}{p})$, where $m$ and $n$ are the dimensions of the
theoretical database matrix, $q$ and $r$ are dimensions of spectra, and $p$ is
the number of processors. We further prove that communication-optimal strategy
with fast-memory $\sqrt{M} = mn + \frac{2qr}{p}$ can achieve
$\Omega({\frac{2mnq}{p}})$ but is not achieved by any existing parallel
proteomics algorithms till date. To validate our claim, we performed a
meta-analysis of published parallel algorithms, and their performance results.
We show that sub-optimal speedups with increasing number of processors is a
direct consequence of not achieving the communication lower-bounds. We further
validate our claim by performing experiments which demonstrate the
communication bounds that are proved in this paper. Consequently, we assert
that next-generation of \emph{provable}, and demonstrated superior parallel
algorithms are urgently needed for MS based large systems-biology studies
especially for meta-proteomics, proteogenomic, microbiome, and proteomics for
non-model organisms. Our hope is that this paper will excite the parallel
computing community to further investigate parallel algorithms for highly
influential MS based omics problems.


  Protein aggregates in the brain play a central role in cognitive decline and
structural damage associated with neurodegenerative diseases. For instance, in
Alzheimer's disease the formation of Amyloid-beta plaques and tau proteins
neurofibrillary tangles follows from the accumulation of different proteins
into large aggregates through specific mechanisms such as nucleation and
elongation. These mechanisms have been studied in vitro where total protein
mass is conserved. However, in vivo, clearance mechanisms may play an important
role in limiting the formation of aggregates. Here, we generalise classical
models of protein aggregation to take into account both production of monomers
and the clearance of protein aggregates. Depending on the clearance model, we
show that there may be a critical clearance value above which aggregation does
not take place. Our result offers further evidence in support of the hypotheses
that clearance mechanisms play a potentially crucial role in neurodegenerative
disease initiation and progression; and as such, are a possible therapeutic
target.


  Identifying similar protein sequences is a core step in many computational
biology pipelines such as detection of homologous protein sequences, generation
of similarity protein graphs for downstream analysis, functional annotation and
gene location. Performance and scalability of protein similarity searches have
proven to be a bottleneck in many bioinformatics pipelines due to increases in
cheap and abundant sequencing data. This work presents a new distributed-memory
software, PASTIS. PASTIS relies on sparse matrix computations for efficient
identification of possibly similar proteins. We use distributed sparse matrices
for scalability and show that the sparse matrix infrastructure is a great fit
for protein similarity searches when coupled with a fully-distributed
dictionary of sequences that allows remote sequence requests to be fulfilled.
Our algorithm incorporates the unique bias in amino acid sequence substitution
in searches without altering the basic sparse matrix model, and in turn,
achieves ideal scaling up to millions of protein sequences.


  We study the evolutionary dynamics of the prisoner's dilemma game in which
cooperators and defectors interact with another actor type called exiters.
Rather than being exploited by defectors, exiters exit the game in favour of a
small payoff. We find that this simple extension of the game allows cooperation
to flourish in well-mixed populations when iterations or reputation are added.
In networked populations, however, the exit option is less conducive to
cooperation. Instead, it enables the coexistence of cooperators, defectors, and
exiters through cyclic dominance. Other outcomes are also possible as the exit
payoff increases or the network structure changes, including network-wide
oscillations in actor abundances that may cause the extinction of exiters and
the domination of defectors, although game parameters should favour exiting.
The complex dynamics that emerges in the wake of a simple option to exit the
game implies that nuances matter even if our analyses are restricted to
incentives for rational behaviour.


  Translation is one of the main steps in the synthesis of proteins. It
consists of ribosomes that translate sequences of nucleotides encoded on mRNA
into polypeptide sequences of amino acids. Ribosomes bound to mRNA move
unidirectionally, while unbound ribosomes diffuse in the cytoplasm. It has been
hypothesized that finite diffusion of ribosomes plays an important role in
ribosome recycling and that mRNA circularization enhances the efficiency of
translation. In order to estimate the effect of cytoplasmic diffusion on the
rate of translation, we consider a Totally Asymmetric Simple Exclusion Process
(TASEP) coupled to a finite diffusive reservoir, which we call the Ribosome
Transport model with Diffusion (RTD). In this model, we derive an analytical
expression for the rate of protein synthesis as a function of the diffusion
constant of ribosomes, which is corroborated with results from continuous-time
Monte Carlo simulations. Using a wide range of biological relevant parameters,
we conclude that diffusion in biological cells is fast enough so that it does
not play a role in controlling the rate of translation initiation.


  Population-level scaling in ecological systems arises from individual growth
and death with competitive constraints. We build on a minimal dynamical model
of metabolic growth where the tension between individual growth and mortality
determines population size distribution. We include resource competition based
on shared capture area separately. By varying relative rates of growth, death,
and competitive attrition, we connect regular and random spatial patterns
across sessile organisms from forests to ants, termites, and fairy circles.
Then, we consider transient temporal dynamics in the context of asymmetric
competition that primarily weakens the smaller of two competitors such as
canopy shading or large colony dominance. When such competition couples slow
timescales of growth with fast competitive death, it generates population shock
waves similar to those observed in forest demographic data. Our minimal
quantitative theory unifies spatiotemporal patterns across sessile organisms
through local competition mediated by the laws of metabolic growth which in
turn result from long-term evolutionary dynamics.


  We study a p-adic reaction-diffusion system and the associated Turing
patterns. We establish an instability criteria and show that the Turing
patterns are not classical patterns consisting of alternating domains. Instead
of this, a Turing pattern consists of several domains (clusters), each of them
supporting a different pattern but with the same parameter values. This type of
patterns are typically produced by reaction-diffusion equations on large
networks.


  In comparison with individual testing, group testing (also known as pooled
testing) is more efficient in reducing the number of tests and potentially
leading to tremendous cost reduction. As indicated in the recent article posted
on the US FDA website, the group testing approach for COVID-19 has received a
lot of interest lately. There are two key elements in a group testing
technique: (i) the pooling matrix that directs samples to be pooled into
groups, and (ii) the decoding algorithm that uses the group test results to
reconstruct the status of each sample. In this paper, we propose a new family
of pooling matrices from packing the pencil of lines (PPoL) in a finite
projective plane. We compare their performance with various pooling matrices
proposed in the literature, including 2D-pooling, P-BEST, and Tapestry, using
the two-stage definite defectives (DD) decoding algorithm. By conducting
extensive simulations for a range of prevalence rates up to 5%, our numerical
results show that there is no pooling matrix with the lowest relative cost in
the whole range of the prevalence rates. To optimize the performance, one
should choose the right pooling matrix, depending on the prevalence rate. The
family of PPoL matrices can dynamically adjust their column weights according
to the prevalence rates and could be a better alternative than using a fixed
pooling matrix.


  Liquid-liquid phase separated (LLPS) states are key to compartmentalise
components in the absence of membranes, however it is unclear whether LLPS
condensates are actively and specifically organized in the sub-cellular space
and by which mechanisms. Here, we address this question by focusing on the
ParABS DNA segregation system, composed of a centromeric-like sequence (parS),
a DNA-binding protein (ParB) and a motor (ParA). We show that parS-ParB
associate to form nanometer-sized, round condensates. ParB molecules diffuse
rapidly within the nucleoid volume, but display confined motions when trapped
inside ParB condensates. Single ParB molecules are able to rapidly diffuse
between different condensates, and nucleation is strongly favoured by parS.
Notably, the ParA motor is required to prevent the fusion of ParB condensates.
These results describe a novel active mechanism that splits, segregates and
localises non-canonical LLPS condensates in the sub-cellular space.


  Mathematical models in epidemiology are an indispensable tool to determine
the dynamics and important characteristics of infectious diseases. Apart from
their scientific merit, these models are often used to inform political
decisions and intervention measures during an ongoing outbreak. However,
reliably inferring the dynamics of ongoing outbreaks by connecting complex
models to real data is still hard and requires either laborious manual
parameter fitting or expensive optimization methods which have to be repeated
from scratch for every application of a given model. In this work, we address
this problem with a novel combination of epidemiological modeling with
specialized neural networks. Our approach entails two computational phases: In
an initial training phase, a mathematical model describing the epidemic is used
as a coach for a neural network, which acquires global knowledge about the full
range of possible disease dynamics. In the subsequent inference phase, the
trained neural network processes the observed data of an actual outbreak and
infers the parameters of the model in order to realistically reproduce the
observed dynamics and reliably predict future progression. With its flexible
framework, our simulation-based approach is applicable to a variety of
epidemiological models. Moreover, since our method is fully Bayesian, it is
designed to incorporate all available prior knowledge about plausible parameter
values and returns complete joint posterior distributions over these
parameters. Application of our method to the early Covid-19 outbreak phase in
Germany demonstrates that we are able to obtain reliable probabilistic
estimates for important disease characteristics, such as generation time,
fraction of undetected infections, likelihood of transmission before symptom
onset, and reporting delays using a very moderate amount of real-world
observations.


  After the first lockdown in response to the COVID-19 outbreak, many countries
faced difficulties in balancing infection control with economics. Due to
limited prior knowledge, economists began researching this issue using
cost-benefit analysis and found that infection control processes significantly
affect economic efficiency. A UK study used economic parameters to numerically
demonstrate an optimal balance in the process, including keeping the infected
population stationary. However, universally applicable knowledge, which is
indispensable for the guiding principles of infection control, has not yet been
clearly developed because of the methodological limitations of simulation
studies. Here, we propose a simple model and theoretically prove the universal
result of economic irreversibility by applying the idea of thermodynamics to
pandemic control. This means that delaying infection control measures is more
expensive than implementing infection control measures early while keeping
infected populations stationary. This implies that once the infected population
increases, society cannot return to its previous state without extra
expenditures. This universal result is analytically obtained by focusing on the
infection-spreading phase of pandemics, and is applicable not just to COVID-19,
regardless of "herd immunity." It also confirms the numerical observation of
stationary infected populations in its optimally efficient process. Our
findings suggest that economic irreversibility is a guiding principle for
balancing infection control with economic effects.


  We propose a method for gene expression based analysis of cancer phenotypes
incorporating network biology knowledge through unsupervised construction of
computational graphs. The structural construction of the computational graphs
is driven by the use of topological clustering algorithms on protein-protein
networks which incorporate inductive biases stemming from network biology
research in protein complex discovery. This structurally constrains the
hypothesis space over the possible computational graph factorisation whose
parameters can then be learned through supervised or unsupervised task
settings. The sparse construction of the computational graph enables the
differential protein complex activity analysis whilst also interpreting the
individual contributions of genes/proteins involved in each individual protein
complex. In our experiments analysing a variety of cancer phenotypes, we show
that the proposed methods outperform SVM, Fully-Connected MLP, and
Randomly-Connected MLPs in all tasks. Our work introduces a scalable method for
incorporating large interaction networks as prior knowledge to drive the
construction of powerful computational models amenable to introspective study.


  The field of Artificial Intelligence (AI) is going through a period of great
expectations, introducing a certain level of anxiety in research, business and
also policy. This anxiety is further energised by an AI race narrative that
makes people believe they might be missing out. Whether real or not, a belief
in this narrative may be detrimental as some stake-holders will feel obliged to
cut corners on safety precautions, or ignore societal consequences just to
"win". Starting from a baseline model that describes a broad class of
technology races where winners draw a significant benefit compared to others
(such as AI advances, patent race, pharmaceutical technologies), we investigate
here how positive (rewards) and negative (punishments) incentives may
beneficially influence the outcomes. We uncover conditions in which punishment
is either capable of reducing the development speed of unsafe participants or
has the capacity to reduce innovation through over-regulation. Alternatively,
we show that, in several scenarios, rewarding those that follow safety measures
may increase the development speed while ensuring safe choices. Moreover, in
{the latter} regimes, rewards do not suffer from the issue of over-regulation
as is the case for punishment. Overall, our findings provide valuable insights
into the nature and kinds of regulatory actions most suitable to improve safety
compliance in the contexts of both smooth and sudden technological shifts.


  Since the first reports of a coronavirus (CoV) disease 2019 (COVID-19) caused
by severe acute respiratory syndrome virus (SARS-CoV-2) in Wuhan, Hubei
province, China, scientists are working around the clock to find sound answers
to the issue of its origin. While the number of scientific articles on
SARS-CoV-2 is increasing, there are still many gaps as to its origin. All
studies failed to find a coronavirus in other animals that is more similar to
human SARS-COV2 than the bat virus, considered to be the primary reservoir. In
this paper we address a new hypothesis, based on a possible recombination
between a DNA and SARS-CoV viruses, to explain the rise of SRAS-CoV-2. By
comparing SARS-CoV-2 and related CoVs with circoviruses (CVs), we found strong
sequence similarity of the genomic region at the 3-end of Bat-CoV ORF1a and the
origin of replication (Ori) of porcine CV type 2 (PCV2), as well as similar RNA
secondary structures of the region encompassing the cleavage site of CoV S gene
with the PCV2 Ori. This constitutes a primary evidence that supports a possible
recombination, which occurrence might explain the origin of SARS-CoV-2.


  Current models for the folding of the human genome see a hierarchy stretching
down from chromosome territories, through A/B compartments and TADs
(topologically-associating domains), to contact domains stabilized by cohesin
and CTCF. However, molecular mechanisms underlying this folding, and the way
folding affects transcriptional activity, remain obscure. Here we review
physical principles driving proteins bound to long polymers into clusters
surrounded by loops, and present a parsimonious yet comprehensive model for the
way the organization determines function. We argue that clusters of active RNA
polymerases and their transcription factors are major architectural features;
then, contact domains, TADs, and compartments just reflect one or more loops
and clusters. We suggest tethering a gene close to a cluster containing
appropriate factors -- a transcription factory -- increases the firing
frequency, and offer solutions to many current puzzles concerning the actions
of enhancers, super-enhancers, boundaries, and eQTLs (expression quantitative
trait loci). As a result, the activity of any gene is directly influenced by
the activity of other transcription units around it in 3D space, and this is
supported by Brownian-dynamics simulations of transcription factors binding to
cognate sites on long polymers.


  The hyperbolic dependence of catalytic rate on substrate concentration is a
classical result in enzyme kinetics, quantified by the celebrated
Michaelis-Menten equation. The ubiquity of this relation in diverse chemical
and biological contexts has recently been rationalized by a graph-theoretic
analysis of deterministic reaction networks. Experiments, however, have
revealed that "molecular noise" - intrinsic stochasticity at the molecular
scale - leads to significant deviations from classical results and to
unexpected effects like "molecular memory", i.e., the breakdown of statistical
independence between turnover events. Here we show, through a new method of
analysis, that memory and non-hyperbolicity have a common source in an initial,
and observably long, transient peculiar to stochastic reaction networks of
multiple enzymes. Networks of single enzymes do not admit such transients. The
transient yields, asymptotically, to a steady-state in which memory vanishes
and hyperbolicity is recovered. We propose new statistical measures, defined in
terms of turnover times, to distinguish between the transient and steady states
and apply these to experimental data from a landmark experiment that first
observed molecular memory in a single enzyme with multiple binding sites. Our
study shows that catalysis at the molecular level with more than one enzyme
always contains a non-classical regime and provides insight on how the
classical limit is attained.


  This work demonstrates how a complete statistical model quantifying the
evolution of pairs of aligned proteins can be constructed from a
time-parameterised substitution matrix and a time-parameterised 3-state
alignment machine. All parameters of such a model can be inferred from any
benchmark data-set of aligned protein sequences. This allows us to examine nine
well-known substitution matrices on six benchmarks curated using various
structural alignment methods; any matrix that does not explicitly model a
"time"-dependent Markov process is converted to a corresponding base-matrix
that does. In addition, a new optimal matrix is inferred for each of the six
benchmarks. Using Minimum Message Length (MML) inference, all 15 matrices are
compared in terms of measuring the Shannon information content of each
benchmark. This has resulted in a new and clear overall best performed
time-dependent Markov matrix, MMLSUM, and its associated 3-state machine, whose
properties we have analysed in this work. For standard use, the MMLSUM series
of (log-odds) \textit{scoring} matrices derived from the above Markov matrix,
are available at https://lcb.infotech.monash.edu.au/mmlsum.


  Real-world examples of periods of periodical organisms range from cicadas
whose life-cycles are larger prime numbers, like 13 or 17, to bamboos whose
periods are large multiples of small primes, like 40 or even 120. The
periodicity is caused by interaction of species, be it a predator-prey
relationship, symbiosis, commensialism, or competition exclusion principle. We
propose a simple mathematical model which explains and models all those
principles, including listed extremel cases. This, rather universal,
qualitative model is based on the concept of a local fitness function, where a
randomly chosen new period is selected if the value of the global fitness
function of the species increases. Arithmetically speaking, the different
observed interactions are related to only four principles: given a couple of
integer periods either (1) their greatest common divisor is one, (2) one of the
periods is prime, (3) both periods are equal, or (4) one period is an integer
multiple of the other.


  The death toll for Covid-19 may be reduced by dividing the population into
two classes, the vulnerable and the fit, with different lockdown regimes.
Instead of one reproduction number there now are four parameters. These make it
possible to quantify the effect of the social distancing measures. There is a
simple stochastic model for epidemics in a two type population. Apart from the
size of the population of the vulnerable and the fit, and the initial number of
infected in the two classes, only the four reproduction parameters are needed
to run the two type Reed-Frost model. The program is simple and fast. On a pc
it takes less than five minutes to do a hundred thousand simulations of the
epidemic for a population of the size of the US. Epidemics are non-linear
processes. Results may be counterintuitive. The average number of vulnerable
persons infected by an infectious fit person is a crucial parameter of the
epidemic in the two type population. Intuitively this parameter should be
small. However simulations show that even if this parameter is small the death
toll may be higher than without shielding. Under certain conditions increasing
the value of the parameter may reduce the death toll. The article addresses
these blind spots in our intuition.


  In response to the worldwide outbreak of the coronavirus disease COVID-19, a
variety of nonpharmaceutical interventions such as face masks and social
distancing have been implemented. A careful assessment of the effects of such
containment strategies is required to avoid exceeding social and economical
costs as well as a dangerous "second wave" of the pandemic. In this work, we
combine a recently developed dynamical density functional theory model and an
extended SIRD model with hysteresis to study effects of various measures and
strategies using realistic parameters. Depending on intervention thresholds, a
variety of phases with different numbers of shutdowns and deaths are found.
Spatiotemporal simulations provide further insights into the dynamics of a
second wave. Our results are of crucial importance for public health policy.


  This study explored how population mobility flows form commuting networks
across US counties and influence the spread of COVID-19. We utilized 3-level
mixed effects negative binomial regression models to estimate the impact of
network COVID-19 exposure on county confirmed cases and deaths over time. We
also conducted weighting-based analyses to estimate the causal effect of
network exposure. Results showed that commuting networks matter for COVID-19
deaths and cases, net of spatial proximity, socioeconomic, and demographic
factors. Different local racial and ethnic concentrations are also associated
with unequal outcomes. These findings suggest that commuting is an important
causal mechanism in the spread of COVID-19 and highlight the significance of
interconnected of communities. The results suggest that local level mitigation
and prevention efforts are more effective when complemented by similar efforts
in the network of connected places. Implications for research on inequality in
health and flexible work arrangements are discussed.


  Pandemic outbreak creates a life threatening situation around the world and
for a while it feels as if it slows down the world. A lot of effort is being
taken to mitigate the spread of the pandemics. The main objective of this paper
is to build a mathematical model of the pandemic spread in India based on the
mass gathering using distribution functions and Bayesian approach. Subsequently
modeled the pandemic spread based on the flights arrived in India using
generalized linear regression. We validated the effect of these events in the
number of confirmed cases in India and formulated the number of confirmed cases
in the positive sense. Subsequently, studied the progression of the infective
using progression series and flattened the curve using an appropriate
convergence criterion. We validated the theoretical aspects with the statistics
released by the Government.


  The mechanical properties of DNA are typically described by elastic theories
with purely local couplings (on-site models). We discuss and analyze
coarse-grained (oxDNA) and all-atom simulations, which indicate that in DNA
distal sites are coupled. Hence, off-site models provide a more realistic
description of the mechanics of the double helix. We show that off-site
interactions are responsible for a length scale dependence of the elasticity,
and we develop an analytical framework to estimate bending and torsional
persistence lengths in models including these interactions. Our simulations
indicate that off-site couplings are particularly strong for certain degrees of
freedom, while they are very weak for others. If stiffness parameters obtained
from DNA data are used, the theory predicts large length scale dependent
effects for torsional fluctuations and a modest effect in bending fluctuations,
which is in agreement with experiments.


  Human societies include diverse social relationships. Friends, family,
business colleagues, and online contacts can all contribute to one's social
life. Individuals may behave differently in different domains, but success in
one domain may engender success in another. Here, we study this problem using
multilayer networks to model multiple domains of social interactions, in which
individuals experience different environments and may express different
behaviors. We provide a mathematical analysis and find that coupling between
layers tends to promote prosocial behavior. Even if prosociality is disfavored
in each layer alone, multilayer coupling can promote its proliferation in all
layers simultaneously. We apply this analysis to six real-world multilayer
networks, ranging from the socio-emotional and professional relationships in a
Zambian community, to the online and offline relationships within an academic
University. We discuss the implications of our results, which suggest that
small modifications to interactions in one domain may catalyze prosociality in
a different domain.


  Identifying one or more biologically-active/native decoys from millions of
non-native decoys is one of the major challenges in computational structural
biology. The extreme lack of balance in positive and negative samples (native
and non-native decoys) in a decoy set makes the problem even more complicated.
Consensus methods show varied success in handling the challenge of decoy
selection despite some issues associated with clustering large decoy sets and
decoy sets that do not show much structural similarity. Recent investigations
into energy landscape-based decoy selection approaches show promises. However,
lack of generalization over varied test cases remains a bottleneck for these
methods. We propose a novel decoy selection method, ML-Select, a machine
learning framework that exploits the energy landscape associated with the
structure space probed through a template-free decoy generation. The proposed
method outperforms both clustering and energy ranking-based methods, all the
while consistently offering better performance on varied test-cases. Moreover,
ML-Select shows promising results even for the decoy sets consisting of mostly
low-quality decoys. ML-Select is a useful method for decoy selection. This work
suggests further research in finding more effective ways to adopt machine
learning frameworks in achieving robust performance for decoy selection in
template-free protein structure prediction.