  We review the key role played by mathematical modelling in elucidating two
centre-finding patterning systems in E. coli: midcell division positioning by
the MinCDE system and DNA partitioning by the ParABS system. We focus
particularly on how, despite much experimental effort, these systems were
simply too complex to unravel by experiments alone, and instead required key
injections of quantitative, mathematical thinking. We conclude the review by
analysing the frequency of modelling approaches in microbiology over time. We
find that while such methods are increasing in popularity, they are still
probably heavily under-utilised for optimal progress on complex biological
questions.


  Stochastic birth-death processes are described as continuous-time Markov
processes in models of population dynamics. A system of infinite, coupled
ordinary differential equations (the so-called master equation) describes the
time-dependence of the probability of each system state. Using a generating
function, the master equation can be transformed into a partial differential
equation. In this contribution we analyze the integrability of two types of
stochastic birth-death processes (with polynomial birth and death rates) using
standard differential Galois theory. We discuss the integrability of the PDE
via a Laplace transform acting over the temporal variable. We show that the PDE
is not integrable except for the (trivial) case in which rates are linear
functions of the number of individuals.


  The understanding of mechanisms that control epigenetic changes is an
important research area in modern functional biology. Epigenetic modifications
such as DNA methylation are in general very stable over many cell divisions.
DNA methylation can however be subject to specific and fast changes over a
short time scale even in non-dividing (i.e. not-replicating) cells. Such
dynamic DNA methylation changes are caused by a combination of active
demethylation and de novo methylation processes which have not been
investigated in integrated models. Here we present a hybrid (hidden) Markov
model to describe the cycle of methylation and demethylation over (short) time
scales. Our hybrid model decribes several molecular events either happening at
deterministic points (i.e. describing mechanisms that occur only during cell
division) and other events occurring at random time points. We test our model
on mouse embryonic stem cells using time-resolved data. We predict methylation
changes and estimate the efficiencies of the different modification steps
related to DNA methylation and demethylation.


  Many cellular responses to surrounding cues require temporally concerted
transcriptional regulation of multiple genes. In prokaryotic cells, a
single-input-module motif with one transcription factor regulating multiple
target genes can generate coordinated gene expression. In eukaryotic cells,
transcriptional activity of a gene is affected by not only transcription
factors but also the epigenetic modifications and three-dimensional chromosome
structure of the gene. To examine how local gene environment and transcription
factor regulation are coupled, we performed a combined analysis of time-course
RNA-seq data of TGF-\b{eta} treated MCF10A cells and related epigenomic and
Hi-C data. Using Dynamic Regulatory Events Miner (DREM), we clustered
differentially expressed genes based on gene expression profiles and associated
transcription factors. Genes in each class have similar temporal gene
expression patterns and share common transcription factors. Next, we defined a
set of linear and radial distribution functions, as used in statistical
physics, to measure the distributions of genes within a class both spatially
and linearly along the genomic sequence. Remarkably, genes within the same
class despite sometimes being separated by tens of million bases (Mb) along
genomic sequence show a significantly higher tendency to be spatially close
despite sometimes being separated by tens of Mb along the genomic sequence than
those belonging to different classes do. Analyses extended to the process of
mouse nervous system development arrived at similar conclusions. Future studies
will be able to test whether this spatial organization of chromosomes
contributes to concerted gene expression.


  A tragedy of the commons (TOC) occurs when individuals acting in their own
self-interest deplete commonly-held resources, leading to a worse outcome than
had they cooperated. Over time, the depletion of resources can change
incentives for subsequent actions. Here, we investigate long-term feedback
between game and environment across a continuum of incentives in an
individual-based framework. We identify payoff-dependent transition rules that
lead to oscillatory TOC-s in stochastic simulations and the mean field limit.
Further extending the stochastic model, we find that spatially explicit
interactions can lead to emergent, localized dynamics, including the
propagation of cooperative wave fronts and cluster formation of both social
context and resources. These dynamics suggest new mechanisms underlying how
TOCs arise and how they might be averted.


  In the context of cancer, internal "checkerboard" structures are normally
found in the matrices of gene expression data, which correspond to genes that
are significantly up- or down-regulated in patients with specific types of
tumors. In this paper, we propose a novel method, called dual
graph-regularization principal component analysis (DGPCA). The main innovation
of this method is that it simultaneously considers the internal geometric
structures of the condition manifold and the gene manifold. Specifically, we
obtain principal components (PCs) to represent the data and approximate the
cluster membership indicators through Laplacian embedding. This new method is
endowed with internal geometric structures, such as the condition manifold and
gene manifold, which are both suitable for bi-clustering. A closed-form
solution is provided for DGPCA. We apply this new method to simultaneously
cluster genes and conditions (e.g., different samples) with the aim of finding
internal "checkerboard" structures on gene expression data, if they exist.
Then, we use this new method to identify regulatory genes under the particular
conditions and to compare the results with those of other state-of-the-art
PCA-based methods. Promising results on gene expression data have been verified
by extensive experiments


  We inquire to what extent protein peptide plane and side chain dynamics can
be reconstructed from knowledge of C-alpha dynamics. Due to lack of
experimental data we analyze all atom molecular dynamics trajectories from
Anton supercomputer, and for clarity we limit our attention to the peptide
plane O atoms and side chain C-beta atoms. We try and reconstruct their
dynamics using four different approaches. Three of these are the publicly
available reconstruction programs Pulchra, Remo Scwrl4. The fourth, Statistical
Method, builds entirely on statistical analysis of Protein Data Bank (PDB)
structures. All four methods place the O and C-beta atoms accurately along the
Anton trajectories. However, the Statistical Method performs best. The results
suggest that under physiological conditions, the all atom dynamics is slaved to
that of C-alpha atoms. The results can help improve all atom force fields, and
advance reconstruction and refinement methods for reduced protein structures.
The results provide impetus for development of effective coarse grained force
fields in terms of reduced coordinates.


  Using a population dynamics inspired by an ensemble of growing cells, a set
of fluctuation theorems linking observables measured at the lineage and
population levels are derived. One of these relations implies specific
inequalities comparing the population doubling time with the mean generation
time at the lineage or population levels. While these inequalities have been
derived before for age controlled models with negligible mother-daughter
correlations, we show that they also hold for a broad class of size-controlled
models. We discuss the implications of this result for the interpretation of a
recent experiment in which the growth of bacteria strains has been probed at
the single cell level.


  For centuries isolation has been the main control strategy of unforeseen
epidemic outbreaks. When implemented in full and without delay, isolation is
very effective. However, flawless implementation is seldom feasible in
practice. We present an epidemic model called SIQ with an isolation protocol,
focusing on the consequences of delays and incomplete identification of
infected hosts. The continuum limit of this model is a system of Delay
Differential Equations, the analysis of which reveals clearly the dependence of
epidemic evolution on model parameters including disease reproductive number,
isolation probability, speed of identification of infected hosts and recovery
rates. Our model offers estimates on minimum response capabilities needed to
curb outbreaks,and predictions of endemic states when containment fails.
Critical response capability is expressed explicitly in terms of parameters
that are easy to obtain, to assist in the evaluation of funding priorities
involving preparedness and epidemics management.


  Summary: Human alpha satellite and satellite 2/3 contribute to several
percent of the human genome. However, identifying these sequences with
traditional algorithms is computationally intensive. Here we develop dna-brnn,
a recurrent neural network to learn the sequences of the two classes of
centromeric repeats. It achieves high similarity to RepeatMasker and is times
faster. Dna-brnn explores a novel application of deep learning and may
accelerate the study of the evolution of the two repeat classes.
  Availability and implementation: https://github.com/lh3/dna-nn
  Contact: hli@jimmy.harvard.edu


  Living in groups brings benefits to many animals, such as a protection
against predators and an improved capacity for sensing and making decisions
while searching for resources in uncertain environments. A body of studies has
shown how collective behaviors within animal groups on the move can be useful
for pooling information about the current state of the environment. The effects
of interactions on collective motion have been mostly studied in models of
agents with no memory. Thus, whether coordinated behaviors can emerge from
individuals with memory and different foraging experiences is still poorly
understood. By means of an agent based model, we quantify how individual memory
and information fluxes can contribute to improving the foraging success of a
group in complex environments. In this context, we define collective learning
as a coordinated change of behavior within a group resulting from individual
experiences and information transfer. We show that an initially scattered
population of foragers visiting dispersed resources can gradually achieve
cohesion and become selectively localized in space around the most salient
resource sites. Coordination is lost when memory or information transfer among
individuals is suppressed. The present modelling framework provides predictions
for empirical studies of collective learning and could also find applications
in swarm robotics and motivate new search algorithms based on reinforcement.


  Cooperation is the cornerstone of human evolutionary success. Like no other
species, we champion the sacrifice of personal benefits for the common good,
and we work together to achieve what we are unable to achieve alone. Knowledge
and information from past generations is thereby often instrumental in ensuring
we keep cooperating rather than deteriorating to less productive ways of
coexistence. Here we present a mathematical model based on evolutionary game
theory that shows how using the past as the benchmark for evolutionary success,
rather than just current performance, significantly improves cooperation in the
future. Interestingly, the details of just how the past is taken into account
play only second-order importance, whether it be a weighted average of past
payoffs or just a single payoff value from the past. Cooperation is promoted
because information from the past disables fast invasions of defectors, thus
enhancing the long-term benefits of cooperative behavior.


  Logical models have been successfully used to describe regulatory and
signaling networks without requiring quantitative data. However, existing data
is insufficient to adequately define a unique model, rendering the
parametrization of a given model a difficult task.
  Here, we focus on the characterization of the set of Boolean functions
compatible with a given regulatory structure, i.e. the set of all monotone
nondegenerate Boolean functions. We then propose an original set of rules to
locally explore the direct neighboring functions of any function in this set,
without explicitly generating the whole set. Also, we provide relationships
between the regulatory functions and their corresponding dynamics.
  Finally, we illustrate the usefulness of this approach by revisiting
Probabilistic Boolean Networks with the model of T helper cell differentiation
from Mendoza & Xenarios.


  Accurate protein structural ensembles can be determined with metainference, a
Bayesian inference method that integrates experimental information with prior
knowledge of the system and deals with all sources of uncertainty and errors as
well as with system heterogeneity. Furthermore, metainference can be
implemented using the metadynamics approach, which enables the computational
study of complex biological systems requiring extensive conformational
sampling. In this chapter, we provide a step-by-step guide to perform and
analyse metadynamic metainference simulations using the ISDB module of the
open-source PLUMED library, as well as a series of practical tips to avoid
common mistakes. Specifically, we will guide the reader in the process of
learning how to model the structural ensemble of a small disordered peptide by
combining state-of-the-art molecular mechanics force fields with nuclear
magnetic resonance data, including chemical shifts, scalar couplings and
residual dipolar couplings.


  The spread of an epidemic process is considered in the context of a spatial
SIR stochastic model that includes a parameter $0\le p\le 1$ that assigns
weights $p$ and $1- p$ to global and local infective contacts respectively. The
model was previously studied by other authors in different contexts. In this
work we characterized the behavior of the system around the threshold for
epidemic spreading. We first used a deterministic approximation of the
stochastic model and checked the existence of a threshold value of $p$ for
exponential epidemic spread. An analytical expression, which defines a function
of the quotient $\alpha$ between the transmission and recovery rates, is
obtained to approximate this threshold. We then performed different analyses
based on intensive stochastic simulations and found that this expression is
also a good estimate for a similar threshold value of $p$ obtained in the
stochastic model. The dynamics of the average number of infected individuals
and the average size of outbreaks show a behavior across the threshold that is
well described by the deterministic approximation. The distributions of the
outbreak sizes at the threshold present common features for all the cases
considered corresponding to different values of $\alpha>1$. These features are
otherwise already known to hold for the standard stochastic SIR model at its
threshold, $\alpha=1$: (i) the probability of having an outbreak of size $n$
goes asymptotically as $n^{-3/2}$ for an infinite system, (ii) the maximal size
of an outbreak scales as $N^{2/3}$ for a finite system of size $N$.


  We investigate the dynamics of the voter model in which the population itself
changes endogenously via the birth-death process. There are two species of
voters, labeled A and B, and the population of each species can grow or shrink
by the birth-death process at equal rates $b$. Individuals of opposite species
also undergo voter model dynamics in which an AB pair can equiprobably become
AA or BB with rate $v$---neutral evolution. In the limit $b/v\to\infty$, the
distribution of consensus times varies as $t^{-3}$ and the probability that the
population size equals $n$ at the moment of consensus varies as $n^{-3}$. As
the birth/death rate $b$ is increased, fixation occurs more more quickly; that
is, population fluctuations promote consensus.


  The actin cortex is a thin layer of actin filaments and myosin motors beneath
the outer membrane of animal cells. It determines the cells' mechanical
properties and forms important morphological structures. Physical descriptions
of the cortex as a contractile active gel suggest that these structures can
result from dynamic instabilities. However, in these analyses the cortex is
described as a two-dimensional layer. Here, we show that the dynamics of the
cortex is qualitatively different when gel fluxes in the direction
perpendicular to the membrane are taken into account. In particular, an
isotropic cortex is then stable for arbitrarily large active stresses. If
lateral contractility exceeds vertical contractility, the system can either
from protrusions with an apparently chaotic dynamics or a periodic static
pattern of protrusions.


  Summary statistics of genome-wide association studies (GWAS) teach causal
relationship between millions of genetic markers and tens and thousands of
phenotypes. However, underlying biological mechanisms are yet to be elucidated.
We can achieve necessary interpretation of GWAS in a causal mediation
framework, looking to establish a sparse set of mediators between genetic and
downstream variables, but there are several challenges. Unlike existing methods
rely on strong and unrealistic assumptions, we tackle practical challenges
within a principled summary-based causal inference framework. We analyzed the
proposed methods in extensive simulations generated from real-world genetic
data. We demonstrated only our approach can accurately redeem causal genes,
even without knowing actual individual-level data, despite the presence of
competing non-causal trails.


  Phenotypic heterogeneity in cancer cells is widely observed and is often
linked to drug resistance. In several cases, such heterogeneity in drug
sensitivity of tumors is driven by stochastic and reversible acquisition of a
drug tolerant phenotype by individual cells even in an isogenic population.
Accumulating evidence further suggests that cell-fate transitions such as the
epithelial to mesenchymal transition (EMT) are associated with drug resistance.
In this study, we analyze stochastic models of phenotypic switching to provide
a framework for analyzing cell-fate transitions such as EMT as a source of
phenotypic variability in drug sensitivity. Motivated by our cell-culture based
experimental observations connecting phenotypic switching in EMT and drug
resistance, we analyze a coarse-grained model of phenotypic switching between
two states in the presence of cytotoxic stress from chemotherapy. We derive
analytical results for time-dependent probability distributions that provide
insights into the rates of phenotypic switching and characterize initial
phenotypic heterogeneity of cancer cells. The results obtained can also shed
light on fundamental questions relating to adaptation and selection scenarios
in tumor response to cytotoxic therapy.


  Cells possess non-membrane-bound bodies, many of which are now understood as
phase-separated condensates. One class of such condensates is composed of two
polymer species, where each consists of repeated binding sites that interact in
a one-to-one fashion with the binding sites of the other polymer. Previous
biologically-motivated modeling of such a two-component system surprisingly
revealed that phase separation is suppressed for certain combinations of
numbers of binding sites. This phenomenon, dubbed the "magic-number effect",
occurs if the two polymers can form fully-bonded small oligomers by virtue of
the number of binding sites in one polymer being an integer multiple of the
number of binding sites of the other. Here we use lattice-model simulations and
analytical calculations to show that this magic-number effect can be greatly
enhanced if one of the polymer species has a rigid shape that allows for
multiple distinct bonding conformations. Moreover, if one species is rigid, the
effect is robust over a much greater range of relative concentrations of the
two species. Our findings advance our understanding of the fundamental physics
of two-component polymer-based phase-separation and suggest implications for
biological and synthetic systems.


  Fifty years ago, Robert MacArthur showed that stable equilibria optimize
quadratic functions of the population sizes in several important ecological
models. Here, we generalize this finding to a broader class of systems within
the framework of contemporary niche theory, and precisely state the conditions
under which an optimization principle (not necessarily quadratic) can be
obtained. We show that conducting the optimization in the space of
environmental states instead of population sizes leads to a universal and
transparent physical interpretation of the objective function. Specifically,
the equilibrium state minimizes the perturbation of the environment induced by
the presence of the competing species, subject to the constraint that no
species has a positive net growth rate. We use this "minimum environmental
perturbation principle" to make new predictions for eco-evolution and community
assembly, and describe a simple experimental setting where its conditions of
validity have been empirically tested.


  Considering that life on earth evolved about 3.7 billion years ago,
vertebrates are young, appearing in the fossil record during the Cambrian
explosion about 542 to 515 million years ago. Results from sequence analyses of
genomes from bacteria, yeast, plants, invertebrates and vertebrates indicate
that receptors for adrenal steroids (aldosterone, cortisol), and sex steroids
(estrogen, progesterone, testosterone) also are young, with receptors for
estrogens and 3-ketosteroids first appearing in basal chordates
(cephalochordates: amphioxus), which are close ancestors of vertebrates. An
ancestral progesterone receptor and an ancestral corticoid receptor, the common
ancestor of the glucocorticoid and mineralocorticoid receptors, evolved in
jawless vertebrates (cyclostomes: lampreys, hagfish). This was followed by
evolution of an androgen receptor and distinct glucocorticoid and
mineralocorticoid receptors in cartilaginous fishes (gnathostomes: sharks).
Adrenal and sex steroid receptors are not found in echinoderms: and
hemichordates, which are ancestors in the lineage of cephalochordates and
vertebrates. The presence of steroid receptors in vertebrates, in which these
steroid receptors act as master switches to regulate differentiation,
development, reproduction, immune responses, electrolyte homeostasis and stress
responses, argues for an important role for steroid receptors in the
evolutionary success of vertebrates, considering that the human genome contains
about 22,000 genes, which is not much larger than genomes of invertebrates,
such as Caenorhabditis elegans (~18,000 genes) and Drosophila (~14,000 genes).


  In this paper, we present the global analysis of a HCV model under therapy.
We prove that the solutions with positive initial values are global, positive,
bounded and not display periodic orbits. In addition, we show that the model is
globally asymptotically stable, by using appropriate Lyapunov functions.


  Dynamical mean field theory (DMFT) is a tool that allows to analyze the
stochastic dynamics of $N$ interacting degrees of freedom in terms of a
self-consistent $1$-body problem. In this work, focusing on models of
ecosystems, we present the derivation of DMFT through the dynamical cavity
method, and we develop a method for solving it numerically. Our numerical
procedure can be applied to a large variety of systems for which DMFT holds. We
implement and test it for the generalized random Lotka-Volterra model, and show
that complex dynamical regimes characterized by chaos and aging can be captured
and studied by this framework.


  Phenotypic variation is a hallmark of cellular physiology. Metabolic
heterogeneity, in particular, underpins single-cell phenomena such as microbial
drug tolerance and growth variability. Much research has focussed on
transcriptomic and proteomic heterogeneity, yet it remains unclear if such
variation permeates to the metabolic state of a cell. Here we propose a
stochastic model to show that complex forms of metabolic heterogeneity emerge
from fluctuations in enzyme expression and catalysis. The analysis predicts
clonal populations to split into two or more metabolically distinct
subpopulations. We reveal mechanisms not seen in deterministic models, in which
enzymes with unimodal expression distributions lead to metabolites with a
bimodal or multimodal distribution across the population. Based on published
data, the results suggest that metabolite heterogeneity may be more pervasive
than previously thought. Our work casts light on links between gene expression
and metabolism, and provides a theory to probe the sources of metabolite
heterogeneity.


  We introduce GeNet, a method for shotgun metagenomic classification from raw
DNA sequences that exploits the known hierarchical structure between labels for
training. We provide a comparison with state-of-the-art methods Kraken and
Centrifuge on datasets obtained from several sequencing technologies, in which
dataset shift occurs. We show that GeNet obtains competitive precision and good
recall, with orders of magnitude less memory requirements. Moreover, we show
that a linear model trained on top of representations learned by GeNet achieves
recall comparable to state-of-the-art methods on the aforementioned datasets,
and achieves over 90% accuracy in a challenging pathogen detection problem.
This provides evidence of the usefulness of the representations learned by
GeNet for downstream biological tasks.


  The Iberian Mediterranean Basin is home to one of the largest groups of
prehistoric rock art sites in Europe. Despite the cultural relevance of
prehistoric Spanish Levantine rock art, pigment composition remains partially
unknown, and the nature of the binders used for painting has yet to be
disclosed. In this work, we present the first omic analysis applied to one of
the flagship Levantine rock art sites: the Valltorta ravine (Castell{\'o}n,
Spain). We used high-throughput sequencing to provide the first description of
the bacterial communities colonizing the rock art patina, which proved to be
dominated by Firmicutes species and might have a protective effect on the
paintings. Proteomic analysis was also performed on rock art microsamples in
order to determine the organic binders present in Levantine prehistoric rock
art pigments. This information could shed light on the controversial dating of
this UNESCO Cultural Heritage, and contribute to defining the chrono-cultural
framework of the societies responsible for these paintings.


  The human neurotensin receptor 1 (NTSR1) is a G protein-coupled receptor that
can be expressed in HEK293T cells by stable transfection. Its ligand is a
13-amino-acid peptide that binds with nanomolar affinity from the extracellular
side to NTSR1. Ligand binding induces conformational changes that trigger the
intracellular signaling processes. Recent single-molecule studies revealed a
dynamic monomer - dimer equilibrium of the receptor in vitro. Here we report on
the oligomerization state of the human NTSR1 in the plasma membrane of HEK293T
cells in vivo. We fused different fluorescent marker proteins mRuby3 or
mNeonGreen to the C-terminus of NTSR1 and mutated a tetracysteine motif into
intracellular loop 3 (ICL3) for subsequent FlAsH labeling. Oligomerization of
NTSR1 was monitored before and after stimulation of the receptor with its
ligand by FLIM and homoFRET microscopy (i.e. Forster resonance energy transfer
between identical fluorophores detected by fluorescence anisotropy), by
colocalization microscopy and by time-lapse imaging using structured
illumination microscopy (SIM).


  We study a mathematical model inspired by a common scenario in northern
Patagonia consisting of humans and their livestock living together with native
wildlife. The main production system in this area is sheep farming, which due
to wrong historical management has led to desertification of the habitat, with
impact on both native wildlife and livestock. Particularly the largest native
herbivore, the guanaco, has reflected in their numbers this alteration induced
in the environment. We analyze a mathematical model that captures the main
characteristics of the interaction between sheep and guanaco: the hierarchical
competition and the advantage granted by humans to the herds, and also
incorporates a dynamic for the habitat. Using the metapopulation formalism, the
trophic web of two herbivores is extended over a patchy landscape that
considers two characteristic times for the dynamics of the resource. Our study
stresses the dependence of the metapopulations dynamics on the recovery time of
the resource. These results are backed up by a deterministic mean field model
which shares some similarities with the stochastic and spatially extended one.
We find different regimes depending on the parameters considered: coexistence
of both species, survival of a single species and extinction of the other, and
extinction of both. Remarkably, in some regions of parameters space we detect
the presence of periodical spatio-temporal patterns, with persistent
oscillations of constant amplitude. Based on these results, we perform a
characterization of the observed scenarios in order to gain insight about the
system.


  Mechanisms to use nanoparticles to separate sunlight into photovoltaic useful
range and thermally useful range to increase the efficiency of solar cells and
to dissipate heat radiatively are discussed based upon lessons we learnt from
photosynthesis. We show that the dual-band maxima in the absorption spectrum of
bacterial light harvestors not only are due to the bacteriochlorophylls
involved but also come from the geometry of the light harvestor. Being able to
manipulate these two bands arbitrarily enables us to fabricate the
nanoparticles required. Such mechanisms are also useful for the design of
remote power charging and light sensors.


  Boolean networks are a popular modeling framework in computational biology to
capture the dynamics of molecular networks, such as gene regulatory networks.
It has been observed that many published models of such networks are defined by
regulatory rules driving the dynamics that have certain so-called canalizing
properties. In this paper, we investigate the dynamics of a random Boolean
network with such properties using analytical methods and simulations.
  From our simulations, we observe that Boolean networks with higher canalizing
depth have generally fewer attractors, the attractors are smaller, and the
basins are larger, with implications for the stability and robustness of the
models. These properties are relevant to many biological applications.
Moreover, our results show that, from the standpoint of the attractor
structure, high canalizing depth, compared to relatively small positive
canalizing depth, has a very modest impact on dynamics.
  Motivated by these observations, we conduct mathematical study of the
attractor structure of a random Boolean network of canalizing depth one (i.e.,
the smallest positive depth). For every positive integer $\ell$, we give an
explicit formula for the limit of the expected number of attractors of length
$\ell$ in an $n$-state random Boolean network as $n$ goes to infinity.


  We train a neural network to predict chemical toxicity based on gene
expression data. The input to the network is a full expression profile
collected either in vitro from cultured cells or in vivo from live animals. The
output is a set of fine grained predictions for the presence of a variety of
pathological effects in treated animals. When trained on the Open TG-GATEs
database it produces good results, outperforming classical models trained on
the same data. This is a promising approach for efficiently screening chemicals
for toxic effects, and for more accurately evaluating drug candidates based on
preclinical data.


  Motivation: Ab initio protein docking represents a major challenge for
optimizing a noisy and costly "black box"-like function in a high-dimensional
space. Despite progress in this field, there is no docking method available for
rigorous uncertainty quantification (UQ) of its solution quality (e.g.
interface RMSD or iRMSD).
  Results: We introduce a novel algorithm, Bayesian Active Learning (BAL), for
optimization and UQ of such black-box functions and flexible protein docking.
BAL directly models the posterior distribution of the global optimum (or native
structures for protein docking) with active sampling and posterior estimation
iteratively feeding each other. Furthermore, we use complex normal modes to
represent a homogeneous Euclidean conformation space suitable for
high-dimension optimization and construct funnel-like energy models for
encounter complexes. Over a protein docking benchmark set and a CAPRI set
including homology docking, we establish that BAL significantly improve against
both starting points by rigid docking and refinements by particle swarm
optimization, providing for one third targets a top-3 near-native prediction.
BAL also generates tight confidence intervals with half range around 25% of
iRMSD and confidence level at 85%. Its estimated probability of a prediction
being native or not achieves binary classification AUROC at 0.93 and AUPRC over
0.60 (compared to 0.14 by chance); and also found to help ranking predictions.
To the best of our knowledge, this study represents the first uncertainty
quantification solution for protein docking, with theoretical rigor and
comprehensive assessment.
  Source codes are available at https://github.com/Shen-Lab/BAL.


  Cancers are mainly caused by somatic genomic alterations (SGAs) that perturb
cellular signaling systems and eventually activate oncogenic processes.
Therefore, understanding the functional impact of SGAs is a fundamental task in
cancer biology and precision oncology. Here, we present a deep neural network
model with encoder-decoder architecture, referred to as genomic impact
transformer (GIT), to infer the functional impact of SGAs on cellular signaling
systems through modeling the statistical relationships between SGA events and
differentially expressed genes (DEGs) in tumors. The model utilizes a
multi-head self-attention mechanism to identify SGAs that likely cause DEGs, or
in other words, differentiating potential driver SGAs from passenger ones in a
tumor. GIT model learns a vector (gene embedding) as an abstract representation
of functional impact for each SGA-affected gene. Given SGAs of a tumor, the
model can instantiate the states of the hidden layer, providing an abstract
representation (tumor embedding) reflecting characteristics of perturbed
molecular/cellular processes in the tumor, which in turn can be used to predict
multiple phenotypes. We apply the GIT model to 4,468 tumors profiled by The
Cancer Genome Atlas (TCGA) project. The attention mechanism enables the model
to better capture the statistical relationship between SGAs and DEGs than
conventional methods, and distinguishes cancer drivers from passengers. The
learned gene embeddings capture the functional similarity of SGAs perturbing
common pathways. The tumor embeddings are shown to be useful for tumor status
representation, and phenotype prediction including patient survival time and
drug response of cancer cell lines.


  Monte Carlo (MC) permutation test is considered the gold standard for
statistical hypothesis testing, especially when standard parametric assumptions
are not clear or likely to fail. However, in modern data science settings where
a large number of hypothesis tests need to be performed simultaneously, it is
rarely used due to its prohibitive computational cost. In genome-wide
association studies, for example, the number of hypothesis tests $m$ is around
$10^6$ while the number of MC samples $n$ for each test could be greater than
$10^8$, totaling more than $nm$=$10^{14}$ samples. In this paper, we propose
Adaptive MC multiple Testing (AMT) to estimate MC p-values and control false
discovery rate in multiple testing. The algorithm outputs the same result as
the standard full MC approach with high probability while requiring only
$\tilde{O}(\sqrt{n}m)$ samples. This sample complexity is shown to be optimal.
On a Parkinson GWAS dataset, the algorithm reduces the running time from 2
months for full MC to an hour. The AMT algorithm is derived based on the theory
of multi-armed bandits.


  Tandem-duplication-random-loss (TDRL) is an important genome rearrangement
operation studied in evolutionary biology. This paper investigates some of the
formal properties of TDRL operations on the symmetric group (the space of
permutations over an $ n $-set). In particular, the cardinality of `balls' of
radius one in the TDRL metric, as well as the cardinality of the maximum
intersection of two such balls, are determined. The corresponding problems for
the so-called mirror (or palindromic) TDRL rearrangement operations are also
solved. The results represent an initial step in the study of error correction
and reconstruction problems in this context and are of potential interest in
DNA-based data storage applications.


  Rapid progress in deep learning has spurred its application to bioinformatics
problems including protein structure prediction and design. In classic machine
learning problems like computer vision, progress has been driven by
standardized data sets that facilitate fair assessment of new methods and lower
the barrier to entry for non-domain experts. While data sets of protein
sequence and structure exist, they lack certain components critical for machine
learning, including high-quality multiple sequence alignments and insulated
training / validation splits that account for deep but only weakly detectable
homology across protein space. We have created the ProteinNet series of data
sets to provide a standardized mechanism for training and assessing data-driven
models of protein sequence-structure relationships. ProteinNet integrates
sequence, structure, and evolutionary information in programmatically
accessible file formats tailored for machine learning frameworks. Multiple
sequence alignments of all structurally characterized proteins were created
using substantial high-performance computing resources. Standardized data
splits were also generated to emulate the difficulty of past CASP (Critical
Assessment of protein Structure Prediction) experiments by resetting protein
sequence and structure space to the historical states that preceded six prior
CASPs. Utilizing sensitive evolution-based distance metrics to segregate
distantly related proteins, we have additionally created validation sets
distinct from the official CASP sets that faithfully mimic their difficulty.
ProteinNet thus represents a comprehensive and accessible resource for training
and assessing machine-learned models of protein structure.


  Genetic regulatory networks control ontogeny. For fifty years Boolean
networks have served as models of such systems, ranging from ensembles of
random Boolean networks as models for generic properties of gene regulation to
working dynamical models of a growing number of sub-networks of real cells. At
the same time, their statistical mechanics has been thoroughly studied. Here we
recapitulate their original motivation in the context of current theoretical
and empirical research. We discuss ensembles of random Boolean networks whose
dynamical attractors model cell types. A sub-ensemble is the critical ensemble.
There is now strong evidence that genetic regulatory networks are dynamically
critical, and that evolution is exploring the critical sub-ensemble. The
generic properties of this sub-ensemble predict essential features of cell
differentiation. In particular, the number of attractors in such networks
scales as the DNA content raised to the 0.63 power. Data on the number of cell
types as a function of the DNA content per cell shows a scaling relationship of
0.88. Thus, the theory correctly predicts a power law relationship between the
number of cell types and the DNA contents per cell, and a comparable slope. We
discuss these new scaling values and show prospects for new research lines for
Boolean networks as a base model for systems biology.


  We study several lattice random walk models with stochastic resetting to
previously visited sites which exhibit a phase transition between an anomalous
diffusive regime and a localization regime where diffusion is suppressed. The
localized phase settles above a critical resetting rate, or rate of memory use,
and the probability density asymptotically adopts in this regime a
non-equilibrium steady state similar to that of the well known problem of
diffusion with resetting to the origin. The transition occurs because of the
presence of a single impurity site where the resetting rate is lower than on
other sites, and around which the walker spontaneously localizes. Near
criticality, the localization length diverges with a critical exponent that
falls in the same class as the self-consistent theory of Anderson localization
of waves in random media. The critical dimensions are also the same in both
problems. Our study provides analytically tractable examples of localization
transitions in path-dependent, reinforced stochastic processes, which can be
also useful for understanding spatial learning by living organisms.


  Proteomics is the large-scale analysis of the proteins. The common method for
identifying proteins and characterising their amino acid sequences is to digest
the proteins into peptides, analyse the peptides using mass spectrometry and
assign the resulting tandem mass spectra (MS/MS) to peptides using database
search tools. However, database search algorithms are highly dependent on a
reference protein database and they cannot identify peptides and proteins not
included in the database. Therefore, de novo sequencing algorithms are
developed to overcome the problem by directly reconstructing the peptide
sequence of an MS/MS spectrum without using any protein database. Current de
novo sequencing algorithms often fail to construct the completely matched
sequences, and produce partial matches. In this study, we propose a genetic
algorithm based method, GA-Novo, to solve the complex optimisation task of de
novo peptide sequencing, aiming at constructing full length sequences. Given an
MS/MS spectrum, GA-Novo optimises the amino acid sequences to best fit the
input spectrum. On the testing dataset, GA-Novo outperforms PEAKS, the most
commonly used software for this task, by constructing 8% higher number of fully
matched peptide sequences, and 4% higher recall at partially matched sequences.


  Artificial RNA molecules with novel functionality have many applications in
synthetic biology, pharmacy and white biotechnology. The de novo design of such
devices using computational methods and prediction tools is a
resource-efficient alternative to experimental screening and selection
pipelines. In this review, we describe methods common to many such
computational approaches, thoroughly dissect these methods and highlight open
questions for the individual steps. Initially, it is essential to investigate
the biological target system, the regulatory mechanism that will be exploited,
as well as the desired components in order to define design objectives.
Subsequent computational design is needed to combine the selected components
and to obtain novel functionality. This process can usually be split into
constrained sequence sampling, the formulation of an optimization problem and
an in silico analysis to narrow down the number of candidates with respect to
secondary goals. Finally, experimental analysis is important to check whether
the defined design objectives are indeed met in the target environment and
detailed characterization experiments should be performed to improve the
mechanistic models and detect missing design requirements.


  Studying evolutionary correlations in alignments of homologous sequences by
means of an inverse Potts model has proven useful to obtain residue-residue
contact energies and to predict contacts in proteins. The quality of the
results depend much on several choices of the detailed model and on the
algorithms used. We built, in a very controlled way, synthetic alignments with
statistical properties similar to those of real proteins, and used them to
assess the performance of different inversion algorithms and of their variants.
Realistic synthetic alignments display typical features of low--temperature
phases of disordered systems, a feature that affects the inversion algorithms.
We showed that a Boltzmann--learning algorithm is computationally feasible and
performs well in predicting the energy of native contacts. However, all
algorithms suffer of false positives quite equally, making the quality of the
prediction of native contacts with the different algorithm much
system--dependent.


  A question in evolutionary biology is why the number of males is
approximately equal to that of females in many species, and Fisher's theory of
equal investment answers that it is the evolutionarily stable state. The
Fisherian mechanism can be given a concrete form by a genetic model based on
the following assumptions: (1) Males and females mate at random. (2) An allele
acts on the father to determine the expected progeny sex ratio. (3) The
offspring inherits the allele from either side of the parents with equal
probability. The model is known to achieve the 1:1 sex ratio due to the
invasion of mutant alleles with different progeny sex ratios. In this study,
however, we argue that mutation plays a more subtle role in that fluctuations
caused by mutation renormalize the sex ratio and thereby keep it away from 1:1
in general. This finding shows how the sex ratio is affected by mutation in a
systematic way, whereby the effective mutation rate can be estimated from an
observed sex ratio.


  G protein-coupled receptors (GPCRs) are a large superfamily of membrane
proteins that are activated by extracellular small molecules or photons.
Neurotensin receptor 1 (NTSR1) is a GPCR that is activated by neurotensin, i.e.
a 13 amino acid peptide. Binding of neurotensin induces conformational changes
in the receptor that trigger the intracellular signaling processes. While
recent single-molecule studies have reported a dynamic monomer - dimer
equilibrium of NTSR1 in vitro, a biophysical characterization of the
oligomerization status of NTSR1 in living mammalian cells is complicated. Here
we report on the oligomerization state of the human NTSR1 tagged with mRuby3 by
dissolving the plasma membranes of living HEK293T cells into 10 nm-sized
soluble lipid nanoparticles by addition of styrene-maleic acid copolymers
(SMALPs). Single SMALPs were analyzed one after another in solution by
multi-parameter single molecule spectroscopy including brightness, fluorescence
lifetime and anisotropy for homoFRET. Brightness analysis was improved using
single SMALP detection in a confocal ABELtrap for extended observation times in
solution. A bimodal brightness distribution indicated a significant fraction of
dimeric NTSR1 in SMALPs or in the plasma membrane, respectively, before
addition of neurotensin.


  Vaccination is crucial for the control of epidemics. Yet it is a social
dilemma since non-vaccinators can benefit from the herd immunity created by the
vaccinators. Thus the optimum vaccination level is not reached via voluntary
vaccination at times. Intensive studies incorporate social networks to study
vaccination behavior, and it is shown that vaccination can be promoted on some
networks. The underlying network, however, is often assumed to be static,
neglecting the dynamical nature of social networks. We investigate the
vaccination behavior on dynamical social networks using both simulations and
mean-field approximations. We find that the more robust the
vaccinator-infected-non-vaccinator links are or the more fragile the
vaccinator-healthy-non-vaccinator links are, the higher the final vaccination
level is. This result is true for arbitrary rationality. Furthermore, we show
that, under strong selection, the vaccination level can be higher than that in
the well-mixed population. In addition, we show that vaccination on evolving
social network is equivalent to the vaccination in well mixed population with a
rescaled basic reproductive ratio. Our results highlight the dynamical nature
of social network on the vaccination behavior, and can be insightful for the
epidemic control.


  Hsp70 molecular chaperones are abundant ATP-dependent nanomachines that
actively reshape non-native, misfolded proteins and assist a wide variety of
essential cellular processes. Here we combine complementary
computational/theoretical approaches to elucidate the structural and
thermodynamic details of the chaperone-induced expansion of a substrate
protein, with a particular emphasis on the critical role played by ATP
hydrolysis. We first determine the conformational free-energy cost of the
substrate expansion due to the binding of multiple chaperones using
coarse-grained molecular simulations. We then exploit this result to implement
a non-equilibrium rate model which estimates the degree of expansion as a
function of the free energy provided by ATP hydrolysis. Our results are in
quantitative agreement with recent single-molecule FRET experiments and
highlight the stark non-equilibrium nature of the process, showing that Hsp70s
are optimized to convert effectively chemical energy into mechanical work close
to physiological conditions.


  Multiple-merger coalescents, e.g. $\Lambda$-$n$-coalescents, have been
proposed as models of the genealogy of $n$ sampled individuals for a range of
populations whose genealogical structures are not captured well by Kingman's
$n$-coalescent. $\Lambda$-$n$-coalescents can be seen as the limit process of
the discrete genealogies of Cannings models with fixed population size, when
time is rescaled and population size $N\to\infty$. As established for Kingman's
$n$-coalescent, moderate population size fluctuations in the discrete
population model should be reflected by a time-change of the limit coalescent.
For $\Lambda$-$n$-coalescents, this has been explicitly shown for only a
limited subclass of $\Lambda$-$n$-coalescents and exponentially growing
populations. This article gives a general construction of time-changed
$\Lambda$-$n$-coalescents as limits of specific Cannings models with rather
arbitrary time changes.


  Membrane Proteins (MPs) account for around 15-39% of the human proteome and
assume a critical role in a vast set of cellular and physiological mechanisms,
including molecular transport, nutrient uptake, toxin and waste product
clearance, respiration, and signaling. While roughly 60% of all FDA-approved
drugs target MPs, there is a shortage of structural and biochemical data on
them mainly hindered by their localization in the lipid bilayer. We present
here MEmbrane protein dimer Novel Structure Analyser database (MENSAdb), a real
time web-application exposing a broad array of fundamental features about MPs
surface and their interfacial regions. In particular, we present conservation,
four distinctive Accessible Solvent Area (ASA) descriptors, average and
environment-specific B-factors, intermolecular contacts at 2.5 and 4.0
angstroms distance cutoffs, salt-bridges, hydrogen-bonds, hydrophobic, pi-pi
interactions, t-stacking and cation-pi interactions. Additionally, users can
closely inspect differences in values between three distinctive residues
classes: i) non-surface, ii) surface and non-interfacial and iii) interfacial.
The database is freely available at www.moreiralab.com/resources/mensadb.


  Near a bifurcation point, the response time of a system is expected to
diverge due to the phenomenon of critical slowing down. We investigate critical
slowing down in well-mixed stochastic models of biochemical feedback by
exploiting a mapping to the mean-field Ising universality class. This mapping
allows us to quantify critical slowing down in experiments where we measure the
response of T cells to drugs. Specifically, the addition of a drug is
equivalent to a sudden quench in parameter space, and we find that quenches
that take the cell closer to its critical point result in slower responses. We
further demonstrate that our class of biochemical feedback models exhibits the
Kibble-Zurek collapse for continuously driven systems, which predicts the
scaling of hysteresis in cellular responses to more gradual perturbations. We
discuss the implications of our results in terms of the tradeoff between a
precise and a fast response.


  Phylogenetic diversity indices provide a formal way to apportion
'evolutionary heritage' across species. Two natural diversity indices are Fair
Proportion (FP) and Equal Splits (ES). FP is also called 'evolutionary
distinctiveness' and, for rooted trees, is identical to the Shapley Value (SV),
which arises from cooperative game theory. In this paper, we investigate the
extent to which FP and ES can differ, characterise tree shapes on which the
indices are identical, and study the equivalence of FP and SV and its
implications in more detail. We also define and investigate analogues of these
indices on unrooted trees (where SV was originally defined), including an index
that is closely related to the Pauplin representation of phylogenetic
diversity.


  Agricultural production is affected by climate extremes, which are increasing
because of global warming. This motivates the need of a proper evaluation of
the agricultural production systems resilience to enhance food security, market
stability, and the general ability of society to cope with the effects of
climate change. Resilience is generally assessed through holistic approaches
involving a large number of indicators for the environmental, social and
economic factors that influence food availability, access and utilization.
Here, we investigate the problem of measuring resilience in a simplified
framework, focusing on the crop production component of the agricultural
system. For an idealized production system composed of a single crop, using the
original definition of resilience, we identify the best combination of the mean
and the variance of annual crop production data to estimate crop resilience
i.e. the crop resilience indicator. Through numerical experiments conducted
with a conceptual crop model, we show the general properties of this indicator
applied to production systems for different levels of adaptation to climate
variability and in case of increasing frequencies of extreme events. Finally,
we discuss the applicability of the proposed approach to real agricultural
production systems and the expected effects of crop diversity on the resilience
of crop production systems, following directly from the mathematical definition
of the crop resilience indicator.


  Individual-based model (IBM) has been used to simulate and to design control
strategies for dynamic systems that are subject to stochasticity and
heterogeneity, such as infectious diseases. In the IBM, an individual is
represented by a set of specific characteristics that may change dynamically
over time. This feature allows a more realistic analysis of the spread of an
epidemic. This paper presents a literature survey of IBM applied to biomedical
and epidemiology research. The main goal is to present existing techniques,
advantages and future perspectives in the development of the model. We
evaluated 89 articles, which mostly analyze interventions aimed at endemic
infections. In addition to the review, an overview of IBM is presented as an
alternative to complement or replace compartmental models, such as the SIR
(Susceptible-Infected-Recovered) model. Numerical simulations also illustrate
the capabilities of IBM, as well as some limitations regarding the effects of
discretization. We show that similar side-effects of discretization scheme for
compartmental models may also occur in IBM, which requires careful attention.


  In this brief, we study epidemic spreading dynamics taking place in complex
networks. We specifically investigate the effect of synergy, where multiple
interactions between nodes result in a combined effect larger than the simple
sum of their separate effects. Although synergistic effects play key roles in
various biological and social phenomena, their analyses have been often
performed by means of approximation techniques and for limited types of
networks. In order to address this limitation, this paper proposes a rigorous
approach to quantitatively understand the effect of synergy in the
Susceptible-Infected-Susceptible model taking place in an arbitrary complex
network. We derive an upper bound on the growth rate of the synergistic
Susceptible-Infected-Susceptible model in terms of the eigenvalues of a matrix
whose size grows quadratically with the number of the nodes in the network. We
confirm the effectiveness of our result by numerical simulations on empirically
observed human and animal social networks.


  Eco-evolutionary frameworks can explain certain features of communities in
which ecological and evolutionary processes occur over comparable timescales.
Here, we investigate whether an evolutionary dynamics may interact with the
spatial structure of a prey-predator community in which both species show
limited mobility and predator perceptual ranges are subject to natural
selection. In these conditions, our results unveil an eco-evolutionary feedback
between species spatial mixing and predators perceptual range: different levels
of mixing select for different perceptual ranges, which in turn reshape the
spatial distribution of prey and its interaction with predators. This emergent
pattern of interspecific interactions feeds back to the efficiency of the
various perceptual ranges, thus selecting for new ones. Finally, since
prey-predator mixing is the key factor that regulates the intensity of
predation, we explore the community-level implications of such feedback and
show that it controls both coexistence times and species extinction
probabilities.


  Recent experimental findings on a protein showed the diffusion of vibrational
energy does not occur along the backbone interaction,as it might be expected,
but prevalently on non-bonded contacts. These results are explained presenting
a theoretical picture, supported by computational calculations, that accounts
for these different behaviors in vibrational energy exchange process showing
the collective motions on the backbone present a $superlocalized$ nature as
their decay with the distance $r$ is $exp(-r^{d})$ with $d \sim 1.8$, whereas
those associated to non-bonded contacts result simply localized with $d \sim
1$.


  We find for the first time the asymptotic representation of the solution to
the space dependent variable order fractional diffusion and Fokker-Planck
equations. We identify a new advection term that causes ultra-slow spatial
aggregation of subdiffusive particles due to dominance over the standard
advection and diffusion terms, in the long-time limit. This uncovers the
anomalous mechanism by which non-uniform distributions can occur. We perform
experiments on intracellular lysosomal distributions and Monte Carlo
simulations and find excellent agreement between the asymptotic solution,
particle histograms and experiments.


  Single-cell gene expression measurements offer opportunities in deriving
mechanistic understanding of complex diseases, including cancer. However, due
to the complex regulatory machinery of the cell, gene regulatory network (GRN)
model inference based on such data still manifests significant uncertainty. The
goal of this paper is to develop optimal classification of single-cell
trajectories accounting for potential model uncertainty. Partially-observed
Boolean dynamical systems (POBDS) are used for modeling gene regulatory
networks observed through noisy gene-expression data. We derive the exact
optimal Bayesian classifier (OBC) for binary classification of single-cell
trajectories. The application of the OBC becomes impractical for large GRNs,
due to computational and memory requirements. To address this, we introduce a
particle-based single-cell classification method that is highly scalable for
large GRNs with much lower complexity than the optimal solution. The
performance of the proposed particle-based method is demonstrated through
numerical experiments using a POBDS model of the well-known T-cell large
granular lymphocyte (T-LGL) leukemia network with noisy time-series
gene-expression data.


  The sensitivity analysis of biological system models can significantly
contribute to identifying and explaining influences of internal or external
changes on model and its elements. We propose here a comprehensive framework to
study sensitivity of intra-cellular networks and to identify key intervention
pathways, by performing both static and dynamic sensitivity analysis. While the
static sensitivity analysis focuses on the impact of network topology and
update functions, the dynamic analysis accounts for context-dependent transient
state distributions. To study sensitivity, we use discrete models, where each
element is represented as a discrete variable and assigned an update rule,
which is a function of element's known direct and indirect regulators. Our
sensitivity analysis framework allows for assessing the effect of context on
individual element sensitivity, as well as on element criticality in reaching
preferred outcomes. The framework also enables discovery of most influential
pathways in the model that are essential for satisfying important system
properties, and thus, could be used for interventions. We discuss the role of
nine different network attributes in identifying key elements and intervention
pathways, and evaluate their performance using model checking method. Finally,
we apply our methods on the model of naive T cell differentiation, and further
demonstrate the importance of context-based sensitivity analysis in identifying
most influential elements and pathways.


  Biological network alignment (NA) aims to find a node mapping between
species' molecular networks that uncovers similar network regions, thus
allowing for transfer of functional knowledge between the aligned nodes.
However, current NA methods do not end up aligning functionally related nodes.
A likely reason is that they assume it is topologically similar nodes that are
functionally related. However, we show that this assumption does not hold well.
So, a paradigm shift is needed with how the NA problem is approached. We
redefine NA as a data-driven framework, TARA (daTA-dRiven network Alignment),
which attempts to learn the relationship between topological relatedness and
functional relatedness without assuming that topological relatedness
corresponds to topological similarity, like traditional NA methods do. TARA
trains a classifier to predict whether two nodes from different networks are
functionally related based on their network topological patterns. We find that
TARA is able to make accurate predictions. TARA then takes each pair of nodes
that are predicted as related to be part of an alignment. Like traditional NA
methods, TARA uses this alignment for the across-species transfer of functional
knowledge. Clearly, TARA as currently implemented uses topological but not
protein sequence information for this task. We find that TARA outperforms
existing state-of-the-art NA methods that also use topological information,
WAVE and SANA, and even outperforms or complements a state-of-the-art NA method
that uses both topological and sequence information, PrimAlign. Hence, adding
sequence information to TARA, which is our future work, is likely to further
improve its performance.


  We introduce a notion of finite sampling consistency for phylogenetic trees
and show that the set of finitely sampling consistent and exchangeable
distributions on n leaf phylogenetic trees is a polytope. We use this polytope
to show that the set of all exchangeable and infinite sampling consistent
distributions on 4 leaf phylogenetic trees is exactly Aldous' beta-splitting
model and give a description of some of the vertices for the polytope of
distributions on 5 leaves. We also introduce a new semialgebraic set of
exchangeable and sampling consistent models we call the multinomial model and
use it to characterize the set of exchangeable and sampling consistent
distributions.


  Unsupervised clustering has broad applications in data stratification,
pattern investigation and new discovery beyond existing knowledge. In
particular, clustering of bioactive molecules facilitates chemical space
mapping, structure-activity studies, and drug discovery. These tasks,
conventionally conducted by similarity-based methods, are complicated by data
complexity and diversity. We ex-plored the superior learning capability of deep
autoencoders for unsupervised clustering of 1.39 mil-lion bioactive molecules
into band-clusters in a 3-dimensional latent chemical space. These
band-clusters, displayed by a space-navigation simulation software, band
molecules of selected bioactivity classes into individual band-clusters
possessing unique sets of common sub-structural features beyond structural
similarity. These sub-structural features form the frameworks of the
literature-reported pharmacophores and privileged fragments. Within each
band-cluster, molecules are further banded into selected sub-regions with
respect to their bioactivity target, sub-structural features and molecular
scaffolds. Our method is potentially applicable for big data clustering tasks
of different fields.


  Classical stochastic demography predicts that environmental stochasticity
reduces population growth rates and, thereby, can increase extinction risk. In
contrast, the SAS-CFF model demonstrates that environmental stochasticity can
promote genetic diversity. Extending the SAS-CFF to account for demography, I
examine the simultaneous effects of environmental stochasticity on genetic
diversity and population persistence. Consistent with Gillespie's analysis, if
the log-fitness function is concave and allelic responses to the environment
are not perfectly correlated, then per-capita growth rates of rare alleles are
positive and genetic diversity is maintained in the sense of stochastic
persistence i.e. allelic frequencies tend to stay away from zero almost-surely
and in probability. Alternatively, if the log-fitness function is convex, then
per-capita growth rates of rare alleles are negative and an allele
asymptotically fixates with probability one. If the population's low-density,
per-capita growth rate is positive, then the population persists in the sense
of stochastic persistence, else it goes asymptotically extinct with probability
one. In contrast to per-capita growth rates of rare alleles, the population's
per-capita growth rate is a decreasing function of the concavity of the
log-fitness function. Moreover, when the log-fitness function is concave,
allelic diversity increases the population's per-capita growth rate while
decreasing the per-capita growth rate of rare alleles, and environmental
stochasticity increases the per-capita growth rate of rare alleles but
decreases the population's per-capita growth rate. Collectively, these results
(i) highlight how mechanisms promoting population persistence may be at odds
with mechanisms promoting genetic diversity, and (ii) provide conditions under
which population persistence relies on existing standing genetic variation.


  Shape had been intuitively recognized to play a dominant role in determining
the global motion patterns of bio-molecular assemblies. However, it is not
clear exactly how shape determines the motion patterns. What about the local
interactions that hold a structure together to a certain shape? The
contributions of global shape and local interactions usually mix together and
are difficult to tease part. In this work, we use symmetry to elucidate the
distinct roles of global shape and local interactions in protein dynamics.
Symmetric complexes provide an ideal platform for this task since in them the
effects of local interactions and global shape are separable, allowing their
distinct roles to be identified. Our key findings based on symmetric assemblies
are: (i) the motion patterns of each subunit are determined primarily by
intra-subunit interactions (IRSi), and secondarily by inter-subunit
interactions (IESi); (ii) the motion patterns of the whole assembly are fully
dictated by the global symmetry/shape and have nothing to do with local iESi or
IRSi. This is followed by a discussion on how the findings may be generalized
to complexes in any shape, with or without symmetry.


  Despite the recognised effectiveness of networks of Marine Protected Areas
(MPAs) as a biodiversity conservation instrument, nowadays MPA network design
frequently disregards the importance of connectivity patterns. In the case of
sedentary marine populations, connectivity stems not only from the stochastic
nature of the physical environment that affects early-life stages dispersal,
but also from the spawning stock attributes that affect the reproductive output
(e.g., passive eggs and larvae) and its survivorship. Early-life stages are
virtually impossible to track in the ocean. Therefore, numerical ocean current
simulations coupled to egg and larval Lagrangian transport models remain the
most common approach for the assessment of marine larval connectivity. Inferred
larval connectivity may be different depending on the type of connectivity
considered; consequently, the prioritisation of sites for marine populations'
conservation might also differ. Here, we introduce a framework for evaluating
and designing MPA networks based on the identification of connectivity hotspots
using graph theoretic analysis. We use as a case of study a network of
open-access areas and MPAs, off Mallorca Island (Spain), and test its
effectiveness for the protection of the painted comber Serranus scriba. Outputs
from network analysis are used to: (1) identify critical areas for improving
overall larval connectivity; (2) assess the impact of species' biological
parameters in network connectivity; and (3) explore alternative MPA
configurations to improve average network connectivity. Results demonstrate the
potential of graph theory to identify non-trivial egg/larval dispersal patterns
and emerging collective properties of the MPA network which are relevant for
increasing protection efficiency.


  Algorithms for community detection are usually stochastic, leading to
different partitions for different choices of random seeds. Consensus
clustering has proven to be an effective technique to derive more stable and
accurate partitions than the ones obtained by the direct application of the
algorithm. However, the procedure requires the calculation of the consensus
matrix, which can be quite dense if (some of) the clusters of the input
partitions are large. Consequently, the complexity can get dangerously close to
quadratic, which makes the technique inapplicable on large graphs. Here we
present a fast variant of consensus clustering, which calculates the consensus
matrix only on the links of the original graph and on a comparable number of
additional node pairs, suitably chosen. This brings the complexity down to
linear, while the performance remains comparable as the full technique.
Therefore, our fast consensus clustering procedure can be applied on networks
with millions of nodes and links.


  Understanding the dynamics of infectious disease spread in a heterogeneous
population is an important factor in designing control strategies. Here, we
develop a novel tensor-driven multi-compartment version of the classic
Susceptible-Infected-Recovered (SIR) model and apply it to Internet data to
reveal information about the complex spatial structure of disease spread. The
model is used to analyze state-level Google search data from the US pertaining
to two viruses, Respiratory Syncytial Virus (RSV), and West Nile Virus (WNV).
We fit the data with correlations of $R^2=0.70$, and $0.52$ for RSV and WNV,
respectively. Although no prior assumptions on spatial structure are made,
human movement patterns in the US explain 27-30\% of the estimated inter-state
transmission rates. The transmission rates within states are correlated with
known demographic indicators, such as population density and average age.
Finally, we show that the patterns of disease load for subsequent seasons can
be predicted using the model parameters estimated for previous seasons and as
few as $7$ weeks of data from the current season. Our results are applicable to
other countries and similar viruses, allowing the identification of disease
spread parameters and prediction of disease load for seasonal viruses earlier
in season.


  The ultimatum game has been a prominent paradigm in studying the evolution of
fairness. It predicts that responders should accept any nonzero offer and
proposers should offer the smallest possible amount according to orthodox game
theory. However, the prediction strongly contradicts with experimental findings
where responders usually reject low offers below $20\%$ and proposers usually
make higher offers than expected. To explain the evolution of such fair
behaviors, we here introduce empathy in group-structured populations by
allowing a proportion $\alpha$ of the population to play empathetic strategies.
Interestingly, we find that for high mutation probabilities, the mean offer
decreases with $\alpha$ and the mean demand increases, implying empathy
inhibits the evolution of fairness. For low mutation probabilities, the mean
offer and demand approach to the fair ones with increasing $\alpha$, implying
empathy promotes the evolution of fairness. Furthermore, under both weak and
strong intensities of natural selection, we analytically calculate the mean
offer and demand for different levels of $\alpha$. Counterintuitively, we
demonstrate that although a higher mutation probability leads to a higher level
of fairness under weak selection, an intermediate mutation probability
corresponds to the lowest level of fairness under strong selection. Our study
provides systematic insights into the evolutionary origin of fairness in
group-structured populations with empathetic strategies.


  One way of investigating how genes affect human traits would be with a
genome-wide association study (GWAS). Genetic markers, known as
single-nucleotide polymorphism (SNP), are used in GWAS. This raises privacy and
security concerns as these genetic markers can be used to identify individuals
uniquely. This problem is further exacerbated by a large number of SNPs needed,
which produce reliable results at a higher risk of compromising the privacy of
participants.
  We describe a method using homomorphic encryption (HE) to perform GWAS in a
secure and private setting. This work is based on a proposed algorithm. Our
solution mainly involves homomorphically encrypted matrix operations and
suitable approximations that adapts the semi-parallel GWAS algorithm for HE. We
leverage the complex space of the CKKS encryption scheme to increase the number
of SNPs that can be packed within a ciphertext. We have also developed a cache
module that manages ciphertexts, reducing the memory footprint.
  We have implemented our solution over two HE open source libraries, HEAAN and
SEAL. Our best implementation took $24.70$ minutes for a dataset with $245$
samples, over $4$ covariates and $10643$ SNPs.
  We demonstrate that it is possible to achieve GWAS with homomorphic
encryption with suitable approximations.


  Long reads produced by third-generation sequencing technologies are used to
construct an assembly (i.e., the subject's genome), which is further used in
downstream genome analysis. Unfortunately, long reads have high sequencing
error rates and a large proportion of bps in these long reads are incorrectly
identified. These errors propagate to the assembly and affect the accuracy of
genome analysis. Assembly polishing algorithms minimize such error propagation
by polishing or fixing errors in the assembly by using information from
alignments between reads and the assembly (i.e., read-to-assembly alignment
information). However, assembly polishing algorithms can only polish an
assembly using reads either from a certain sequencing technology or from a
small assembly. Such technology-dependency and assembly-size dependency require
researchers to 1) run multiple polishing algorithms and 2) use small chunks of
a large genome to use all available read sets and polish large genomes. We
introduce Apollo, a universal assembly polishing algorithm that scales well to
polish an assembly of any size (i.e., both large and small genomes) using reads
from all sequencing technologies (i.e., second- and third-generation). Our goal
is to provide a single algorithm that uses read sets from all available
sequencing technologies to improve the accuracy of assembly polishing and that
can polish large genomes. Apollo 1) models an assembly as a profile hidden
Markov model (pHMM), 2) uses read-to-assembly alignment to train the pHMM with
the Forward-Backward algorithm, and 3) decodes the trained model with the
Viterbi algorithm to produce a polished assembly. Our experiments with real
read sets demonstrate that Apollo is the only algorithm that 1) uses reads from
any sequencing technology within a single run and 2) scales well to polish
large assemblies without splitting the assembly into multiple parts.


  Calcium (Ca2+) plays a central role in mediating both contractile function
and hypertrophic signalling in ventricular cardiomyocytes. L-type Ca2+ channels
trigger release of Ca2+ from ryanodine receptors (RyRs) for cellular
contraction, while signalling downstream of Gq coupled receptors stimulates
Ca2+ release via inositol 1,4,5-trisphosphate receptors (IP3Rs), engaging
hypertrophic signalling pathways. Modulation of the amplitude, duration, and
duty cycle of the cytosolic Ca2+ contraction signal, and spatial localisation,
have all been proposed to encode this hypertrophic signal. Given current
knowledge of IP3Rs, we develop a model describing the effect of functional
interaction (cross-talk) between RyR and IP3R channels on the Ca2+ transient,
and examine the sensitivity of the Ca2+ transient shape to properties of IP3R
activation. A key result of our study is that IP3R activation increases Ca2+
transient duration for a broad range of IP3R properties, but the effect of IP3R
activation on Ca2+ transient amplitude is dependent on IP3 concentration.
Furthermore we demonstrate that IP3-mediated Ca2+ release in the cytosol
increases the duty cycle of the Ca2+ transient, the fraction of the cycle for
which [Ca2+] is elevated, across a broad range of parameter values and IP3
concentrations. When coupled to a model of downstream transcription factor
(NFAT) activation, we demonstrate that there is a high correspondence between
the Ca transient duty cycle and the proportion of activated NFAT in the
nucleus. These findings suggest increased cytosolic Ca2+ duty cycle as a
plausible mechanism for IP3-dependent hypertrophic signalling via
Ca2+-sensitive transcription factors such as NFAT in ventricular
cardiomyocytes.


  A biomimetic model of cell-cell communication was developed to probe the
passive molecular transport across ion channels inserted in synthetic lipid
bilayers formed between contacting droplets arranged in a linear array.
Diffusion of a fluorescent probe across the array was measured for different
pore concentrations. The diffusion characteristic time scale is found to vary
non-linearly with the pore concentration. Our measurements are successfully
modeled by a continuous time random walk description, whose waiting time is the
first exit time from a droplet through a cluster of pores. The size of the
cluster of pores is found to increase with their concentration. Our results
provide a direct link between the mesoscopic permeation properties and the
microscopic characteristics of the pores such as their number, size and spatial
arrangement.


  Many developmental processes in biology utilize Notch-Delta signaling to
construct an ordered pattern of cellular differentiation. This signaling
modality is based on nearest-neighbor contact, as opposed to the more familiar
mechanism driven by the release of diffusible ligands. Here, exploiting this
"juxtacrine" property, we present an exact treatment of the pattern formation
problem via a system of nine coupled ordinary differential equations. The
possible patterns that are realized for realistic parameters can be analyzed by
considering a co-dimension 2 pitchfork bifurcation of this system. This
analysis explains the observed prevalence of hexagonal patterns with high Delta
at their center, as opposed to those with central high Notch levels. We show
that outside this range of parameters, in particular for low cis-coupling, a
novel kind of pattern is produced, where high Delta cells have high Notch as
well. It also suggests that the biological system is only weakly first order,
so that an additional mechanism is required to generate the observed
defect-free patterns. We construct a simple strategy for producing such
defect-free patterns.


  Selective inference is considered for testing trees and edges in phylogenetic
tree selection from molecular sequences. This improves the previously proposed
approximately unbiased test by adjusting the selection bias when testing many
trees and edges at the same time. The newly proposed selective inference
$p$-value is useful for testing selected edges to claim that they are
significantly supported if $p>1-\alpha$, whereas the non-selective $p$-value is
still useful for testing candidate trees to claim that they are rejected if
$p<\alpha$. The selective $p$-value controls the type-I error conditioned on
the selection event, whereas the non-selective $p$-value controls it
unconditionally. The selective and non-selective approximately unbiased
$p$-values are computed from two geometric quantities called signed distance
and mean curvature of the region representing tree or edge of interest in the
space of probability distributions. These two geometric quantities are
estimated by fitting a model of scaling-law to the non-parametric multiscale
bootstrap probabilities. Our general method is applicable to a wider class of
problems; phylogenetic tree selection is an example of model selection, and it
is interpreted as the variable selection of multiple regression, where each
edge corresponds to each predictor. Our method is illustrated in a previously
controversial phylogenetic analysis of human, rabbit and mouse.


  As cancer advances, cells often spread from the primary tumor to other parts
of the body and form metastases. This is the main cause of cancer related
mortality. Here we investigate a conceptually simple model of metastasis
formation where metastatic lesions are initiated at a rate which depends on the
size of the primary tumor. The evolution of each metastasis is described as an
independent branching process. We assume that the primary tumor is resected at
a given size and study the earliest time at which any metastasis reaches a
minimal detectable size. The parameters of our model are estimated
independently for breast, colorectal, headneck, lung and prostate cancers. We
use these estimates to compare predictions from our model with values reported
in clinical literature. For some cancer types, we find a remarkably wide range
of resection sizes such that metastases are very likely to be present, but none
of them are detectable. Our model predicts that only very early resections can
prevent recurrence, and that small delays in the time of surgery can
significantly increase the recurrence probability.


  Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs which play a
significant role in several biological processes. RNA-seq based transcriptome
sequencing has been extensively used for identification of lncRNAs. However,
accurate identification of lncRNAs in RNA-seq datasets is crucial for exploring
their characteristic functions in the genome as most coding potential
computation (CPC) tools fail to accurately identify them in transcriptomic
data. Well-known CPC tools such as CPC2, lncScore, CPAT are primarily designed
for prediction of lncRNAs based on the GENCODE, NONCODE and CANTATAdb
databases. The prediction accuracy of these tools often drops when tested on
transcriptomic datasets. This leads to higher false positive results and
inaccuracy in the function annotation process. In this study, we present a
novel tool, PLIT, for the identification of lncRNAs in plants RNA-seq datasets.
PLIT implements a feature selection method based on L1 regularization and
iterative Random Forests (iRF) classification for selection of optimal
features. Based on sequence and codon-bias features, it classifies the RNA-seq
derived FASTA sequences into coding or long non-coding transcripts. Using L1
regularization, 31 optimal features were obtained based on lncRNA and
protein-coding transcripts from 8 plant species. The performance of the tool
was evaluated on 7 plant RNA-seq datasets using 10-fold cross-validation. The
analysis exhibited superior accuracy when evaluated against currently available
state-of-the-art CPC tools.


  Statistical methods for genomewide association studies (GWAS) continue to
improve. However, the increasing volume and variety of genetic and genomic data
make computational speed and ease of data manipulation mandatory in future
software. In our view, a collaborative effort of statistical geneticists is
required to develop open source software targeted to genetic epidemiology. Our
attempt to meet this need is called the OPENMENDELproject
(https://openmendel.github.io). It aims to (1) enable interactive and
reproducible analyses with informative intermediate results, (2) scale to big
data analytics, (3) embrace parallel and distributed computing, (4) adapt to
rapid hardware evolution, (5) allow cloud computing, (6) allow integration of
varied genetic data types, and (7) foster easy communication between
clinicians, geneticists, statisticians, and computer scientists. This article
reviews and makes recommendations to the genetic epidemiology community in the
context of the OPENMENDEL project.


  High throughput sequencing of RNA (RNA-Seq) can provide us with millions of
short fragments of RNA transcripts from a sample. How to better recover the
original RNA transcripts from those fragments (RNA-Seq assembly) is still a
difficult task. For example, RNA-Seq assembly tools typically require
hyper-parameter tuning to achieve good performance for particular datasets.
This kind of tuning is usually unintuitive and time-consuming. Consequently,
users often resort to default parameters, which do not guarantee consistent
good performance for various datasets.
  Here we propose BOAssembler (https://github.com/olivomao/boassembler), a
framework that enables end-to-end automatic tuning of RNA-Seq assemblers, based
on Bayesian Optimization principles. Experiments show this data-driven approach
is effective to improve the overall assembly performance. The approach would be
helpful for downstream (e.g. gene, protein, cell) analysis, and more broadly,
for future bioinformatics benchmark studies.


  Statistical inference of evolutionary parameters from molecular sequence data
relies on coalescent models to account for the shared genealogical ancestry of
the samples. However, inferential algorithms do not scale to available data
sets. A strategy to improve computational efficiency is to rely on simpler
coalescent and mutation models, resulting in smaller hidden state spaces. An
estimate of the cardinality of the state-space of genealogical trees at
different resolutions is essential to decide the best modeling strategy for a
given dataset. To our knowledge, there is neither an exact nor approximate
method to determine these cardinalities. We propose a sequential importance
sampling algorithm to estimate the cardinality of the space of genealogical
trees under different coalescent resolutions. Our sampling scheme proceeds
sequentially across the set of combinatorial constraints imposed by the data.
We analyse the cardinality of different genealogical tree spaces on simulations
to study the settings that favor coarser resolutions. We estimate the
cardinality of genealogical tree spaces from mtDNA data from the 1000 genomes
and a sample from a Melanesian population to illustrate the settings in which
it is advantageous to employ coarser resolutions.


  Non-uniform rates of morphological evolution and evolutionary increases in
organismal complexity, captured in metaphors like "adaptive zones", "punctuated
equilibrium" and "blunderbuss patterns", require more elaborate explanations
than a simple gradual accumulation of mutations. Here we argue that non-uniform
evolutionary increases in phenotypic complexity can be caused by a
threshold-like response to growing ecological pressures resulting from
evolutionary diversification at a given level of complexity. Acquisition of a
new phenotypic feature allows an evolving species to escape this pressure but
can typically be expected to carry significant physiological costs. Therefore,
the ecological pressure should exceed a certain level to make such an
acquisition evolutionarily successful. We present a detailed quantitative
description of this process using a microevolutionary competition model as an
example. The model exhibits sequential increases in phenotypic complexity
driven by diversification at existing levels of complexity and the resulting
increase in competitive pressure, which can push an evolving species over the
barrier of physiological costs of new phenotypic features.


  Overstretching of B-DNA is currently understood as force-induced melting.
Depending on the geometry of the stretching experiment, the force threshold for
the overstretching transition is around 65 or 110 pN. Although the mechanisms
behind force-induced melting have been correctly described by Rouzina and
Bloomfield \cite{RouzinaBloomfield2001a}, neither force threshold has been
exactly calculated by theory. In this work, a detailed analysis of the
force-extension curve is presented, based on a description of single-stranded
DNA in terms of the discrete Kratky-Porod model, consistent with (i) the
contour length expected from the crystallographically determined monomer
distance, and (ii) a high value of the elastic stretch modulus arising from
covalent bonding. The value estimated for the ss-DNA persistence length,
$\lambda = 1.0 $ nm, is at the low end of currently known estimates and
reflects the intrinsic stiffness of the partially, or fully stretched state,
where electrostatic repulsion effects are expected to be minimal. A detailed
analysis of single- and double-stranded DNA free energies provides estimates of
the overstretching force thresholds. In the unconstrained geometry, the
predicted threshold is 64 pN. In the constrained geometry, after allowing for
the entropic penalty of the plectonemic topology of the molten state, the
predicted threshold is 111 pN.


  Homochirality is a generic and unique property of all biochemical life and
the fractional circular polarization of light it induces therefore constitutes
a potentially unambiguous biosignature.} However, while high-quality circular
polarimetric spectra can be easily and quickly obtained in the laboratory,
accurate measurements in the field are much more challenging due to large
changes in illumination and target movement. In this study we have measured
various targets in the field, up to distances of a few kilometers, using the
dedicated circular spectropolarimeter TreePol. We show how photosynthetic life
can readily be distinguished from abiotic matter. We underline the potential of
circular polarization signals as a remotely accessible means to characterize
and monitor terrestrial vegetation, e.g. for agriculture and forestry.
Additionally, we discuss the potential of circular polarization for the remote
detection of extraterrestrial life.


  The highly charged RNA molecules, with each phosphate carrying a single
negative charge, cannot fold into well-defined architectures with tertiary
interactions, in the absence of ions. For ribozymes, divalent cations are known
to be more efficient than monovalent ions in driving them to a compact state
although Mg$^{2+}$ ions are needed for catalytic activities. Therefore, how
ions interact with RNA is relevant in understanding RNA folding. It is often
thought that most of the ions are territorially and non-specifically bound to
the RNA, as predicted by the counterion condensation (CIC) theory. Here, we
show using simulations of ${\it Azoarcus}$ ribozyme, based on an accurate
coarse-grained Three Site Interaction (TIS) model, with explicit divalent and
monovalent cations, that ion condensation is highly specific and depends on the
nucleotide position. The regions with high coordination between the phosphate
groups and the divalent cations are discernible even at very low Mg$^{2+}$
concentrations when the ribozyme does not form tertiary interactions.
Surprisingly, these regions also contain the secondary structural elements that
nucleate subsequently in the self-assembly of RNA, implying that ion
condensation is determined by the architecture of the folded state. These
results are in sharp contrast to interactions of ions (monovalent and divalent)
with rigid charged rods in which ion condensation is uniform and position
independent. The differences are explained in terms of the dramatic
non-monotonic shape fluctuations in the ribozyme as it folds with increasing
Mg$^{2+}$ or Ca$^{2+}$ concentration.


  One challenge in gene network inference is distinguishing between direct and
indirect regulation. Some algorithms, including ARACNE and Phixer, approach
this problem by using pruning methods to eliminate redundant edges in an
attempt to explain the observed data with the simplest possible network
structure. However, we hypothesize that there may be a cost in accuracy to
simplifying the predicted networks in this way, especially due to the
prevalence of redundant connections, such as feed forward loops, in gene
networks. In this paper, we evaluate the pruning methods of ARACNE and Phixer,
and score their accuracy using receiver operating characteristic curves and
precision-recall curves. Our results suggest that while pruning can be useful
in some situations, it may have a negative effect on overall accuracy that has
not been previously studied. Researchers should be aware of both the advantages
and disadvantages of pruning when inferring networks, in order to choose the
best inference strategy for their experimental context.


  We present molecular phylogenetic investigation of Thuidiaceae, especially on
Thudium and Pelekium. Three chloroplast sequences (trnL-F, rps4, and atpB-rbcL)
and one nuclear sequence (ITS) were analyzed. Data partitions were analyzed
separately and in combination by employing MP (maximum parsimony) and Bayesian
methods. The influence of data conflict in combined analyses was further
explored by two methods: the incongruence length difference (ILD) test and the
partition addition bootstrap alteration approach (PABA). Based on the results,
ITS 1& 2 had crucial effect in phylogenetic reconstruction in this study, and
more chloroplast sequences should be combinated into the analyses since their
stability for reconstructing within genus of pleurocarpous mosses. We supported
that Helodiaceae including Actinothuidium, Bryochenea, and Helodium still
attributed to Thuidiaceae, and the monophyletic Thuidiaceae s. lat. should also
include several genera (or species) from Leskeaceae such as Haplocladium and
Leskea. In the Thuidiaceae, Thuidium and Pelekium were resolved as two
monophyletic groups separately. The results from molecular phylogeny were
supported by the crucial morphological characters in Thuidiaceae s. lat.,
Thuidium and Pelekium.


  Native shortcut networks (SCN0) are sub-graphs of native protein residue
networks (PRN0). In this paper, we propose the Network Dynamics (ND) model,
which reconstructs a PRN0 by adding back its edges according to some recipe,
while its nodes are fixed at native (PDB) locations. A PRN0 reconstruction will
eventually reconstruct its SCN0, but only after exploring several non-native
shortcut network (SCN) configurations. It is these other SCN configurations
that are of interest to us, as they produce the statistics to evaluate the
different recipes. The recipes vary from each other slightly to investigate the
effect of different edge orderings on protein folding. An edge ordering is
deemed more successful if it produces a stronger correlation with experimental
folding rates. Over proteins of different chain fold types, this basic
requirement is best satisfied when a recipe favours earlier restoration of
edges with smaller sequence separation, and earlier restoration of edges
incident on nodes with larger remaining degree within the 1-hop neighborhood of
a previously selected node. Further, this recipe generated more route-like
trajectories over SCN space and a wider range of simulated folding rates, both
of which signal cooperative two-state folding behavior. It also produced better
correspondence between ND calculated phi-values and experimental phi-values
obtained from a set of ten transition state ensembles (TSE). Also introduced is
the local centrality measure which uses centrality of initial fold
substructures to yield calculated phi-values from PRN0s with even better
overall correspondence than ND calculated phi-values. Preliminary investigation
with mixed ND recipes by chain fold type suggests that concocting further
context-sensitive recipes by chain fold type and combining them could improve
ND. This search task may be better undertaken by machine learning.