  Protein crystal production is a major bottleneck for the structural
characterisation of proteins. To advance beyond large-scale screening, rational
strategies for protein crystallization are crucial. Understanding how chemical
anisotropy (or patchiness) of the protein surface due to the variety of amino
acid side chains in contact with solvent, contributes to protein protein
contact formation in the crystal lattice is a major obstacle to predicting and
optimising crystallization. The relative scarcity of sophisticated theoretical
models that include sufficient detail to link collective behaviour, captured in
protein phase diagrams, and molecular level details, determined from
high-resolution structural information is a further barrier. Here we present
two crystals structures for the P23TR36S mutant of gamma D-crystallin, each
with opposite solubility behaviour, one melts when heated, the other when
cooled. When combined with the protein phase diagram and a tailored patchy
particle model we show that a single temperature dependent interaction is
sufficient to stabilise the inverted solubility crystal. This contact, at the
P23T substitution site, relates to a genetic cataract and reveals at a
molecular level, the origin of the lowered and retrograde solubility of the
protein. Our results show that the approach employed here may present an
alternative strategy for the rationalization of protein crystallization.


  Bacterial DNA gyrase introduces negative supercoils into chromosomal DNA and
relaxes positive supercoils introduced by replication and transiently by
transcription. Removal of these positive supercoils is essential for
replication fork progression and for the overall unlinking of the two duplex
DNA strands, as well as for ongoing transcription. To address how gyrase copes
with these topological challenges, we used high-speed single-molecule
fluorescence imaging in live Escherichia coli cells. We demonstrate that at
least 300 gyrase molecules are stably bound to the chromosome at any time, with
~12 enzymes enriched near each replication fork. Trapping of reaction
intermediates with ciprofloxacin revealed complexes undergoing catalysis. Dwell
times of ~2 s were observed for the dispersed gyrase molecules, which we
propose maintain steady-state levels of negative supercoiling of the
chromosome. In contrast, the dwell time of replisome-proximal molecules was ~8
s, consistent with these catalyzing processive positive supercoil relaxation in
front of the progressing replisome.


  We propose an active walker model for the motion of individual ants
communicating via chemical signals. It is assumed that communication takes the
form of a time-dependent pheromone field that feedbacks into the motion ants
through tropotaxis: individuals can sense the gradient of the pheromone
concentration field and adjust their orientation accordingly. The individual
model takes the form of a Langevin system of equations in polar coordinates
driven by two-dimensional Gaussian fluctuations and with orientation changes in
response to two pheromone fields: one emanating from the nest, and other
actively produced by ants in their nest-bound journey after finding a food
source. We explicitly track the evolution of both fields in three dimensions.
The proposed tropotaxis model relating the pheromone field to the orientation
changes is similar to Weber's law, but depends explicitly only on the gradient
of the pheromone concentration. We identify ranges of values for the model
parameters that yield the emergence of two key foraging patterns: successful
recruitment to newly found sources, and colony-wide trail networks.


  This paper describes a method to efficiently retrieve protein database
sequences similar to a query sequence, while allowing for significant numbers
of mutations. We call this method SEQR for SEQuence Retrieval. This approach
increases the speed of sequence similarity searches by an order of magnitude
compared to conventional algorithms at the expense of sensitivity. Furthermore,
retrieval time increases less than linearly with the number of sequences, a
desirable property during an era when next generation sequencing technologies
have yielded greater than exponential increases in sequence records. The lower
sensitivity of the algorithm for distantly related sequences compared to
benchmarks is not intrinsic to the method itself, but rather due to the
procedure used to construct the indexing terms, and may be improved. The
indexing terms themselves can be added to standard information retrieval
engines, enabling complex queries that include sequence similarity and other
descriptors such as taxonomy and text descriptions.


  We consider self-assembly of proteins into a virus capsid by the methods of
molecular dynamics. The capsid corresponds either to SPMV or CCMV and is
studied with and without the RNA molecule inside. The proteins are flexible and
described by the structure-based coarse-grained model augmented by
electrostatic interactions. Previous studies of the capsid self-assembly
involved solid objects of a supramolecular scale, e.g. corresponding to
capsomeres, with engineered couplings and stochastic movements. In our
approach, a single capsid is dissociated by an application of a high
temperature for a variable period and then the system is cooled down to allow
for self-assembly. The restoration of the capsid proceeds to various extent,
depending on the nature of the dissociated state, but is rarely complete
because some proteins depart too far unless the process takes place in a
confined space.


  As of October 12, 2018, 211 cases of Ebola virus disease (EVD) were reported
in North Kivu Province, Democratic Republic of Congo. Since the beginning of
October the outbreak has largely shifted into regions in which active armed
conflict is occurring, and in which EVD cases and their contacts are difficult
for health workers to reach. We modeled EVD transmission using a branching
process with gradually quenching transmission estimated from past EVD
outbreaks, with outbreak trajectories conditioned on agreement with the course
of the current outbreak, and with multiple levels of vaccination coverage. We
used an autoregression for short-term projections, a regression model for final
sizes, and a simple Gott's law rule as an ensemble of forecasts. Short-term
model projections were validated against actual case counts. During validation
of short-term projections, models consistently scored higher on shorter-term
forecasts. Based on case counts as of October 13, the stochastic model
projected a median case count of 226 by October 27 (95% prediction interval:
205-268) and 245 by November 10 (95% PI: 208-315), while the auto-regression
model projected median case counts of 240 (95% PI: 215-307) and 259 (95% PI:
216-395) for those dates, respectively. Projected median final counts range
from 274 to 421. Except for Gott's law, the projected probability of an
outbreak surpassing 2013-2016 is exceedingly small. The stochastic model
estimates that vaccine coverage in this outbreak is lower than reported in its
trial. Based on our projections we believe that the epidemic had not yet peaked
at the time of these estimates, though an outbreak like 2013-2016 is not
likely. We estimate that transmission rates are higher than under target levels
of vaccine coverage, and this model estimate may offer a surrogate indicator
for the outbreak response challenges.


  The elastic network (EN) is a prime model that describes the long-time
dynamics of biomolecules. However, the use of harmonic potentials renders this
model insufficient for studying large conformational changes of proteins (e.g.
stretching of proteins, folding and thermal unfolding). Here, we extend the
capabilities of the EN model by using a harmonic approximation described by
Lennard-Jones (LJ) interactions for far contacts and native contacts obtained
from the standard overlap criterion as in the case of Go-like models. While our
model is validated against the EN model by reproducing the equilibrium
properties for a number of proteins, we also show that the model is suitable
for the study of large conformation changes by providing various examples. In
particular, this is illustrated on the basis of pulling simulations that
predict with high accuracy the experimental data on the rupture force of the
studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling
simulations highlight the advantages of our model with respect to Go-like
approaches, where the latter fail to reproduce previous results obtained by
all-atom simulations that predict an additional characteristic peak for this
protein. In addition, folding simulations of small peptides yield different
folding times for alpha-helix and beta-hairpin, in agreement with experiment,
in this way providing further opportunities for the application of our model in
studying large conformational changes of proteins. In contrast to the EN model,
our model is suitable for both normal mode analysis and molecular dynamics
simulation. We anticipate that the proposed model will find applications in a
broad range of problems in biology, including, among others, protein folding
and thermal unfolding.


  Amino acid sequence portrays most intrinsic form of a protein and expresses
primary structure of protein. The order of amino acids in a sequence enables a
protein to acquire a particular stable conformation that is responsible for the
functions of the protein. This relationship between a sequence and its function
motivates the need to analyse the sequences for predicting protein functions.
Early generation computational methods using BLAST, FASTA, etc. perform
function transfer based on sequence similarity with existing databases and are
computationally slow. Although machine learning based approaches are fast, they
fail to perform well for long protein sequences (i.e., protein sequences with
more than 300 amino acid residues). In this paper, we introduce a novel method
for construction of two separate feature sets for protein sequences based on
analysis of 1) single fixed-sized segments and 2) multi-sized segments, using
bi-directional long short-term memory network. Further, model based on proposed
feature set is combined with the state of the art Multi-lable Linear
Discriminant Analysis (MLDA) features based model to improve the accuracy.
Extensive evaluations using separate datasets for biological processes and
molecular functions demonstrate promising results for both single-sized and
multi-sized segments based feature sets. While former showed an improvement of
+3.37% and +5.48%, the latter produces an improvement of +5.38% and +8.00%
respectively for two datasets over the state of the art MLDA based classifier.
After combining two models, there is a significant improvement of +7.41% and
+9.21% respectively for two datasets compared to MLDA based classifier.
Specifically, the proposed approach performed well for the long protein
sequences and superior overall performance.


  Intracellular transport of organelles is fundamental to cell function and
health. The mounting evidence suggests that this transport is in fact
anomalous. However, the reasons for the anomaly is still under debate. We
examined experimental trajectories of organelles inside a living cell and
propose a mathematical model that describes the previously reported transition
from sub-diffusive to super-diffusive motion. In order to explain
super-diffusive behaviour at long times, we introduce non-Markovian detachment
kinetics of the cargo: the rate of detachment is inversely proportional to the
time since the last attachment. Recently, we observed the non-Markovian
detachment rate experimentally in eukaryotic cells. Here we further discuss
different scenarios of how this effective non-Markovian detachment rate could
arise. The non-Markovian model is successful in simultaneously describing the
time-averaged variance (the time-averaged mean squared displacement corrected
for directed motion), the mean first passage time of trajectories and the
multiple peaks observed in the distributions of cargo velocities. We argue that
non-Markovian kinetics could be biologically beneficial compared to the
Markovian kinetics commonly used for modelling, by increasing the average
distance the cargoes travel when a microtubule is blocked by other filaments.
In turn, sub-diffusion allows cargoes to reach neighbouring filaments with
higher probability, which promotes active motion along the microtubules.


  Beta-diversity has been repeatedly shown to decline with increasing
elevation, but the causes of this pattern remain unclear, partly because they
are confounded by coincident variation in alpha- and gamma-diversity. We used
8,795 forest vegetation-plot records from the Czech National Phytosociological
Database to compare the observed patterns of beta diversity to null-model
expectations (beta-deviation) controlling for the effects of alpha- and
gamma-diversity. We tested whether \b{eta}-diversity patterns along a 1,200 m
elevation gradient exclusively depend on the effect of varying species pool
size, or also on the variation of the magnitude of community assembly
mechanisms determining the distribution of species across communities (e.g.,
environmental filtering, dispersal limitation). The null model we used is a
novel extension of an existing null-model designed for presence/absence data
and was specifically designed to disrupt the effect of community assembly
mechanisms, while retaining some key features of observed communities such as
average species richness and species abundance distribution. Analyses were
replicated in ten subregions with comparable elevation ranges. Beta-diversity
declined along the elevation gradient due to a decrease in gamma-diversity,
which was steeper than the decrease in alpha-diversity. This pattern persisted
after controlling for alpha- and gamma-diversity variation, and the results
were robust when different resampling schemes and diversity metrics were used.
We conclude that in temperate forests the pattern of decreasing beta-diversity
with elevation does not exclusively depend on variation in species pool size,
as has been hypothesized, but also on variation in community assembly
mechanisms. The results were consistent across resampling schemes and diversity
measures, thus supporting the use of vegetation plot databases for
understanding...


  The 3D organisation of the genome in interphase cells is not a randomly
folded polymer. Rather, experiments show that chromosomes arrange into a
network of 3D compartments that correlate with biological processes, such as
transcription, chromatin modifications, and protein binding. However, these
compartments do not exist during cell division when the DNA is condensed, and
it is unclear how and when they emerge. In this paper, we focus on the early
stages after cell-division as the chromosomes start to decondense. We use a
simple polymer model to understand the types of 3D structures that emerge from
local unfolding of a compact initial state. From simulations, we recover 3D
compartments, such as TADs and A/B compartments, that are consistently detected
in Chromosome Capture Experiments across cell types and organisms. This
suggests that the large-scale 3D organisation is a result of an inflation
process.


  Binding of neutral Pheophorbide-a methyl ester (MePheo-a) to various
synthetic polynucleotides, double-stranded poly(A)-poly(U), poly(G)-poly(C) and
four-stranded poly(G), as well as to calf thymus DNA, was studied using the
methods of absorption and polarized fluorescent spectroscopy. Measurements were
performed in aqueous bufferes solutions (pH 6.9) of low ionic strength (2 mM
Na+) in a wide range of molar phosphate-to-dye ratios (P/D). Absorption and
fluorescence characteristics of complexes formed between the dye and
biopolymers were determined. Binding of MePheo-a to four-stranded poly(G) is
shown to be accompanied by the most significant spectral transformations:
hypochromism of dye absorption,large batochromic shift of Soret absorption band
(~26 nm) and fluorescence band (~9 nm) maxima, 48-fold enhancement of the dye
emission intensity. In contrast, its binding to the double-stranded
polynucleotides and native DNA induces only small shifts of absorption and
fluorescence bands, as well as no more than 4-fold rise of fluorescence
intensity. Substantial spectral changes and high value of fluorescence
polarization degree (0.26) observed upon binding of MePheo-a to quadruplex
poly(G) allow us to suggest the intercalation of the dye chromophore between
guanine tetrads. At the same time, small spectral changed and insignificant
increase of MePheo-a fluorescence polarization degree (0.12) upon binding of
the dye to double-stranded biopolymers point to another binding type.
Incorporation of MePheo-a to a helix groove is supposed to occur, presumably in
the dimeric form. Substantial enhancement of MePheo-a emission upon binding to
four-stranded poly(G) allows us to propose this compouns as a new fluorescent
probe G-quadruplex structures.


  Forecasting elections -- a challenging, high-stakes problem -- is the subject
of much uncertainty, subjectivity, and media scrutiny. To shed light on this
process, we develop a method for forecasting elections from the perspective of
dynamical systems. Our model borrows ideas from epidemiology, and we use
polling data from United States elections to determine its parameters.
Surprisingly, our general model performs as well as popular forecasters for the
2012 and 2016 U.S. races for president, senators, and governors. Although
contagion and voting dynamics differ, our work suggests a valuable approach to
elucidate how elections are related across states. It also illustrates the
effect of accounting for uncertainty in different ways, provides an example of
data-driven forecasting using dynamical systems, and suggests avenues for
future research on political elections. We conclude with our forecasts for the
senatorial and gubernatorial races on 6~November 2018, which we posted on 5
November 2018.


  Ecological networks originating as a result of three different ecological
processes are examined and cross-compared to assess if the underlying
ecological processes in these systems produce considerable difference in the
structure of the networks. Absence of any significant difference in the
structure of the networks may indicate towards the possibility of a universal
structural pattern in these ecological networks. The underlying graphs of the
networks derived by the ecological processes, namely host-parasite interaction,
plant pollination and seed dispersion are all bipartite graphs and thus several
algebraic structural measures fail to distinguish between the structure of
these networks. In this work we use weighted spectral distribution (WSD) of
normalized graph Laplacian, which have been effectively used earlier to
discriminate graphs with different topologies, to investigate the possibility
of existence of structural dissimilarity in these networks. Graph spectrum is
often considered a signature of the graph and WSD of the graph Laplacian is
shown to be related to the distribution of some small subgraphs in a graph and
hence represent the global structure of a network.We use random projections of
WSD to $\mathbb{R}^{2}$ and $\mathbb{R}^{3}$ and establish that the structure
of plant pollinator networks is significantly different as compared to
host-parasite and seed dispersal networks. The structures of host parasite
networks and seed dispersal networks are found to be identical. Furthermore, we
use some algebraic structural measures in order to quantify the differences as
well as similarities observed in the structure of the three kinds of networks.
We thus infer that our work suggests an absence of a universal structural
pattern in these three different kinds of networks.


  Biochemical reactions involving three or more reactants, called
higher-molecular reactions, play an important role in theoretical systems and
synthetic biology. In particular, such reactions underpin a variety of
important bio-dynamical phenomena, such as multi-stability/multi-modality,
oscillations, bifurcations, and noise-induced effects. However, only reactions
with at most two reactants, called bi-molecular reactions, are experimentally
feasible. To bridge the gap, in this paper we put forward an algorithm for
systematically approximating arbitrary higher-molecular reactions with
bi-molecular ones, while preserving the underlying stochastic dynamics.
Properties of the algorithm and convergence are established via singular
perturbation theory. The algorithm is applied to a variety of higher-molecular
biochemical networks, and is shown to play an important role in
nucleic-acid-based synthetic biology.


  A population is considered stationary if the growth rate is zero and the age
structure is constant. It thus follows that a population is considered
non-stationary if either its growth rate is non-zero and/or its age structure
is non-constant. We propose three properties that are related to the stationary
population identity (SPI) of population biology by connecting it with
stationary populations and non-stationary populations which are approaching
stationarity. One of these important properties is that SPI can be applied to
partition a population into stationary and non-stationary components. These
properties provide deeper insights into cohort formation in real-world
populations and the length of the duration for which stationary and
non-stationary conditions hold. The new concepts are based on the time gap
between the occurrence of stationary and non-stationary populations within the
SPI framework that we refer to as Oscillatory SPI and the Amplitude of SPI.
This article will appear in Bulletin of Mathematical Biology (Springer)


  Study of network dynamics is very active area in biological and social
sciences. However, the relationship between the network structure and the
attractors of the dynamics has not been fully understood yet. In this study, we
numerically investigated the role of degenerate self-loops on the attractors
and its basin size using the budding yeast cell-cycle network model. In the
network, all self-loops negatively surpress the node (self-inhibition loops)
and the attractors are only fixed points, i.e. point attractors. It is found
that there is a simple division rule of the state space by removing the
self-loops when the attractors consist only of point attractors. The point
attractor with largest basin size is robust against the change of the
self-inhibition loop. Furthermore, some limit cycles of period 2 appear as new
attractor when a self-activation loop is added to the original network. It is
also shown that even in that case, the point attractor with largest basin size
is robust.


  The Minimal Ancestral Deviation (MAD) method is a recently introduced
procedure for estimating the root of a phylogenetic tree, based only on the
shape and branch lengths of the tree. The method is loosely derived from the
midpoint rooting method, but, unlike its predecessor, makes use of all pairs of
OTUs when positioning the root. In this note we establish properties of this
method and then describe a fast and memory efficient algorithm. As a proof of
principle, we use our algorithm to determine the MAD roots for simulated
phylogenies with up to 100,000 OTUs. The calculations take a few minutes on a
standard laptop.


  Reconstruction of population histories is a central problem in population
genetics. Existing coalescent-based methods, like the seminal work of Li and
Durbin (Nature, 2011), attempt to solve this problem using sequence data but
have no rigorous guarantees. Determining the amount of data needed to correctly
reconstruct population histories is a major challenge. Using a variety of tools
from information theory, the theory of extremal polynomials, and approximation
theory, we prove new sharp information-theoretic lower bounds on the problem of
reconstructing population structure -- the history of multiple subpopulations
that merge, split and change sizes over time. Our lower bounds are exponential
in the number of subpopulations, even when reconstructing recent histories. We
demonstrate the sharpness of our lower bounds by providing algorithms for
distinguishing and learning population histories with matching dependence on
the number of subpopulations. Along the way and of independent interest, we
essentially determine the optimal number of samples needed to learn an
exponential mixture distribution information-theoretically, proving the upper
bound by analyzing natural (and efficient) algorithms for this problem.


  Estimating the human longevity and computing of life expectancy are central
to the population dynamics. These aspects were studied seriously by scientists
since fifteenth century, including renowned astronomer Edmund Halley. From
basic principles of population dynamics, we propose a method to compute life
expectancy from incomplete data.


  The phase transition of $(PolyA/PolyT)_{100}$ duplex into the denaturated
state is studied in the Peyrard-Bishop-Dauxois model by the method of direct
molecular-dynamical modeling. The temperature dependencies of the total energy
and heat capacity of the duplex are calculated. The approach applied can be
used to calculate the statistical properties of the duplexes of any length and
nucleotide composition.


  Contact-assisted protein folding has made very good progress, but two
challenges remain. One is accurate contact prediction for proteins lack of many
sequence homologs and the other is that time-consuming folding simulation is
often needed to predict good 3D models from predicted contacts. We show that
protein distance matrix can be predicted well by deep learning and then
directly used to construct 3D models without folding simulation at all. Using
distance geometry to construct 3D models from our predicted distance matrices,
we successfully folded 21 of the 37 CASP12 hard targets with a median family
size of 58 effective sequence homologs within 4 hours on a Linux computer of 20
CPUs. In contrast, contacts predicted by direct coupling analysis (DCA) cannot
fold any of them in the absence of folding simulation and the best CASP12 group
folded 11 of them by integrating predicted contacts into complex,
fragment-based folding simulation. The rigorous experimental validation on 15
CASP13 targets show that among the 3 hardest targets of new fold our
distance-based folding servers successfully folded 2 large ones with <150
sequence homologs while the other servers failed on all three, and that our ab
initio folding server also predicted the best, high-quality 3D model for a
large homology modeling target. Further experimental validation in CAMEO shows
that our ab initio folding server predicted correct fold for a membrane protein
of new fold with 200 residues and 229 sequence homologs while all the other
servers failed. These results imply that deep learning offers an efficient and
accurate solution for ab initio folding on a personal computer.


  CRISPR-Cas mediated immunity in bacteria allows bacterial populations to
protect themselves against pathogens. However, it also exposes them to the
dangers of auto-immunity by developing protection that targets its own genome.
Using a simple model of the coupled dynamics of phage and bacterial
populations, we explore how acquisition rates affect the survival rate of the
bacterial colony. We find that the optimal strategy depends on the initial
population sizes of both viruses and bacteria. Additionally, certain
combinations of acquisition and dynamical rates and initial population sizes
guarantee protection, due to a dynamical balance between the evolving
population sizes, without relying on acquisition of viral spacers. Outside this
regime, the high cost of auto-immunity limits the acquisition rate. We discuss
these optimal survival strategies in terms of recent experiments.


  The lncRNA HOTAIR has been implicated in several human cancers. Here, we
evaluated the molecular alterations and upstream regulatory mechanisms of
HOTAIR in glioma, the most common primary brain tumors, and its clinical
relevance. HOTAIR gene expression, methylation, copy-number and prognostic
value were investigated in human gliomas integrating data from online datasets
and our cohorts. High levels of HOTAIR were associated with higher grades of
glioma, particularly IDH wild-type cases. Mechanistically, HOTAIR was
overexpressed in a gene dosage-independent manner, while DNA methylation levels
of particular CpGs in HOTAIR locus were associated with HOTAIR expression
levels in GBM clinical specimens and cell lines. Concordantly, the
demethylating agent 5-Aza-2'-deoxycytidine affected HOTAIR transcriptional
levels in a cell line-dependent manner. Importantly, HOTAIR was frequently
co-expressed with HOXA9 in high-grade gliomas from TCGA, Oncomine, and our
Portuguese and French datasets. Integrated in silico analyses, chromatin
immunoprecipitation, and qPCR data showed that HOXA9 binds directly to the
promoter of HOTAIR. Clinically, GBM patients with high HOTAIR expression had a
significantly reduced overall survival, independently of other prognostic
variables. In summary, this work reveals HOXA9 as a novel direct regulator of
HOTAIR, and establishes HOTAIR as an independent prognostic marker, providing
new therapeutic opportunities to treat this highly aggressive cancer.


  Parental allele-specific expression of imprinted genes is mediated by
imprinting control regions (ICRs) that are constitutively marked by DNA
methylation imprints on the maternal or paternal allele. Mono-allelic DNA
methylation is strictly required for the process of imprinting and has to be
faithfully maintained during the entire lifespan. While the regulation of DNA
methylation itself is well understood, the mechanisms whereby the opposite
allele remains unmethylated are unclear. Here, we show that in the mouse, at
maternally methylated ICRs, the paternal allele, which is constitutively
associated with H3K4me2/3, is marked by default by H3K27me3 when these ICRs are
transcriptionally inactive, leading to the formation of a bivalent chromatin
signature. Our data suggest that at ICRs, chromatin bivalency has a protective
role by ensuring that DNA on the paternal allele remains unmethylated and
protected against spurious and unscheduled gene expression. Moreover , they
provide the proof of concept that, beside pluripotent cells, chromatin
bivalency is the default state of transcriptionally inactive CpG island
promoters , regardless of the developmental stage, thereby contributing to
protect cell identity.


  CpG islands (CGI) marked by bivalent chromatin in stem cells are believed to
be more prone to aberrant DNA methylation in tumor cells. The robustness and
genome-wide extent of this instructive program in different cancer types remain
to be determined. To address this issue we developed a user-friendly approach
to integrate the stem cell chromatin signature in customized DNA methylation
analyses. We used publicly available ChIP-sequencing datasets of several human
embryonic stem cell (hESC) lines to determine the extent of bivalent chromatin
genome-wide. We then created annotated lists of high-confidence bivalent,
H3K4me3-only and H3K27me3-only chromatin regions. The main features of bivalent
regions included localization in CGI/promoters, depletion in retroelements and
enrichment in specific histone modifications, including the poorly
characterized H3K23me2 mark. Moreover, bivalent promoters could be classified
in three clusters based on PRC2 and PolII complexes occupancy. Genes with
bivalent promoters of the PRC2-defined cluster displayed the lowest expression
upon differentiation. As proof-of-concept, we assessed the DNA methylation
pattern of eight types of tumors and confirmed that aberrant cancer-associated
DNA hypermethylation preferentially targets CGI characterized by bivalent
chromatin in hESCs. We also found that such aberrant DNA hypermethylation
affected particularly bivalent CGI/promoters associated with genes that tend to
remain repressed upon differentiation. Strikingly, bivalent CGI were the most
affected by aberrant DNA hypermethylation in both CpG Island Methylator
Phenotype-positive (CIMP+) and CIMP-negative tumors, suggesting that, besides
transcriptional silencing in the pre-tumorigenic cells, the bivalent chromatin
signature in hESCs is a key determinant of the instructive program for aberrant
DNA methylation.


  Current methods to optimize vaccine dose are purely empirically based,
whereas in the drug development field, dosing determinations use far more
advanced quantitative methodology to accelerate decision-making. Applying these
established methods in the field of vaccine development may reduce the
currently large clinical trial sample sizes, long time frames, high costs, and
ultimately have a better potential to save lives. We propose the field of
immunostimulation/immunodynamic (IS/ID) modelling, which aims to translate
mathematical frameworks used for drug dosing towards optimizing vaccine dose
decision-making. Analogous to PK/PD modelling, IS/ID modelling approaches apply
mathematical models to describe the underlying mechanisms by which the immune
response is stimulated by vaccination (IS) and the resulting measured immune
response dynamics (ID). To move IS/ID modelling forward, existing datasets and
further data on vaccine allometry and dose-dependent dynamics need to be
generated and collate, requiring a collaborative environment with input from
academia, industry, regulators, governmental and non-governmental agencies to
share modelling expertise, and connect modellers to vaccine data.


  It has recently been suggested that the fundamental haploid-diploid cycle of
eukaryotic sex exploits a rudimentary form of the Baldwin effect. Thereafter
the other associated phenomena can be explained as evolution tuning the amount
and frequency of learning experienced by an organism. Using the well-known NK
model of fitness landscapes it is here shown that the emergence of dominance
can also be explained under this view of eukaryotic evolution.


  We use topological data analysis as a tool to analyze the fit of mathematical
models to experimental data. This study is built on data obtained from motion
tracking groups of aphids in [Nilsen et al., PLOS One, 2013] and two random
walk models that were proposed to describe the data. One model incorporates
social interactions between the insects, and the second model is a control
model that excludes these interactions. We compare data from each model to data
from experiment by performing statistical tests based on three different sets
of measures. First, we use time series of order parameters commonly used in
collective motion studies. These order parameters measure the overall
polarization and angular momentum of the group, and do not rely on a priori
knowledge of the models that produced the data. Second, we use order parameter
time series that do rely on a priori knowledge, namely average distance to
nearest neighbor and percentage of aphids moving. Third, we use computational
persistent homology to calculate topological signatures of the data. Analysis
of the a priori order parameters indicates that the interactive model better
describes the experimental data than the control model does. The topological
approach performs as well as these a priori order parameters and better than
the other order parameters, suggesting the utility of the topological approach
in the absence of specific knowledge of mechanisms underlying the data.


  We propose the chemlambda artificial chemistry, whose behavior strongly
suggests that real molecules which embed Interaction Nets patterns and real
chemical reactions which resemble Interaction Nets graph rewrites could be a
realistic path towards molecular computers, in the sense explained in the
article.


  Alzheimer's disease is the most common cause of dementia. It is the
fifth-leading cause of death among elderly people. With high genetic
heritability (79%), finding disease causal genes is a crucial step in find
treatment for AD. Following the International Genomics of Alzheimer's Project
(IGAP), many disease-associated genes have been identified; however, we don't
have enough knowledge about how those disease-associated genes affect gene
expression and disease-related pathways. We integrated GWAS summary data from
IGAP and five different expression level data by using TWAS method and
identified 15 disease causal genes under strict multiple testing (alpha<0.05),
4 genes are newly identified; identified additional 29 potential disease causal
genes under false discovery rate(alpha < 0.05), 21 of them are newly
identified. Many genes we identified are also associated with some autoimmune
disorder.


  Most protocols for the high-throughput directed evolution of enzymes rely on
random encapsulation to link phenotype and genotype. In order to optimize these
approaches, or compare one to another, one needs a measure of their performance
at extracting the best variants. We introduce here a new metric named the
Selection Quality Index (SQI), which can be computed from a simple mock
experiment with a known initial fraction of active variants. As opposed to
previous approaches, our index integrates the random co-encapsulation of
entities in compartments and comes with a straightforward experimental
interpretation. We further show how this new metric can be used to extract
general trends of protocol efficiency, or reveal hidden mechanisms such as a
counterintuitive form of beneficial poisoning in the Compartmentalized
Self-Replication protocol.


  Eukaryotic cells transmit information by signaling through complex networks
of interacting proteins. Here we develop a theoretical and computational
framework that relates the biophysics of protein-protein interactions (PPIs)
within a signaling network to its information processing properties. To do so,
we generalize statistical physics-inspired models for protein binding to
account for interactions that depend on post-translational state (e.g.
phosphorylation). By combining these models with information-theoretic methods,
we find that PPIs are a key determinant of information transmission within a
signaling network, with weak interactions giving rise to "noise" that
diminishes information transmission. While noise can be mitigated by increasing
interaction strength, the accompanying increase in transmission comes at the
expense of a slower dynamical response. This suggests that the biophysics of
signaling protein interactions give rise to a fundamental "speed-information"
trade-off. Surprisingly, we find that cross-talk between pathways in complex
signaling networks do not significantly alter information capacity--an
observation that may partially explain the promiscuity and ubiquity of weak
PPIs in heavily interconnected networks. We conclude by showing how our
framework can be used to design synthetic biochemical networks that maximize
information transmission, a procedure we dub "InfoMax" design.


  Currently, third-generation sequencing techniques, which allow to obtain much
longer DNA reads compared to the next-generation sequencing technologies, are
becoming more and more popular. There are many possibilities to combine data
from next-generation and third-generation sequencing.
  Herein, we present a new application called dnaasm-link for linking contigs,
a result of \textit{de novo} assembly of second-generation sequencing data,
with long DNA reads. Our tool includes an integrated module to fill gaps with a
suitable fragment of appropriate long DNA read, which improves the consistency
of the resulting DNA sequences. This feature is very important, in particular
for complex DNA regions, as presented in the paper. Finally, our implementation
outperforms other state-of-the-art tools in terms of speed and memory
requirements, which may enable the usage of the presented application for
organisms with a large genome, which is not possible in~existing applications.
  The presented application has many advantages as (i) significant memory
optimization and reduction of computation time (ii) filling the gaps through
the appropriate fragment of a specified long DNA read (iii) reducing number of
spanned and unspanned gaps in the existing genome drafts.
  The application is freely available to all users under GNU Library or Lesser
General Public License version 3.0 (LGPLv3). The demo application, docker image
and source code are available at http://dnaasm.sourceforge.net.


  Networks are fundamental building blocks for representing data, and
computations. Remarkable progress in learning in structurally defined (shallow
or deep) networks has recently been achieved. Here we introduce evolutionary
exploratory search and learning method of topologically flexible networks under
the constraint of producing elementary computational steady-state input-output
operations.
  Our results include; (1) the identification of networks, over four orders of
magnitude, implementing computation of steady-state input-output functions,
such as a band-pass filter, a threshold function, and an inverse band-pass
function. Next, (2) the learned networks are technically controllable as only a
small number of driver nodes are required to move the system to a new state.
Furthermore, we find that the fraction of required driver nodes is constant
during evolutionary learning, suggesting a stable system design. (3), our
framework allows multiplexing of different computations using the same network.
For example, using a binary representation of the inputs, the network can
readily compute three different input-output functions. Finally, (4) the
proposed evolutionary learning demonstrates transfer learning. If the system
learns one function A, then learning B requires on average less number of steps
as compared to learning B from tabula rasa.
  We conclude that the constrained evolutionary learning produces large robust
controllable circuits, capable of multiplexing and transfer learning. Our study
suggests that network-based computations of steady-state functions,
representing either cellular modules of cell-to-cell communication networks or
internal molecular circuits communicating within a cell, could be a powerful
model for biologically inspired computing. This complements conceptualizations
such as attractor based models, or reservoir computing.


  It is well known that the kinetics of an intracellular biochemical network is
stochastic. This is due to intrinsic noise arising from the random timing of
biochemical reactions in the network as well as due to extrinsic noise stemming
from the interaction of unknown molecular components with the network and from
the cell's changing environment. While there are many methods to study the
effect of intrinsic noise on the system dynamics, few exist to study the
influence of both types of noise. Here we show how one can extend the
conventional linear-noise approximation to allow for the rapid evaluation of
the molecule numbers statistics of a biochemical network influenced by
intrinsic noise and by slow lognormally distributed extrinsic noise. The theory
is applied to simple models of gene regulatory networks and its validity
confirmed by comparison with exact stochastic simulations. In particular we
show how extrinsic noise modifies the dependence of the variance of the
molecule number fluctuations on the rate constants, the mutual information
between input and output signalling molecules and the robustness of
feed-forward loop motifs.


  Attempting to recognize a tree inside a phylogenetic network is a fundamental
undertaking in evolutionary analysis. In the last few years, therefore,
tree-based phylogenetic networks, which are defined by a spanning tree called a
subdivision tree, have attracted attention of theoretical biologists. However,
the application of such networks is still not easy, due to many problems whose
time complexities are not clearly understood. In this paper, we provide a
general framework for solving those various old or new problems from a coherent
perspective, rather than analyzing the complexity of each individual problem or
developing an algorithm one by one. More precisely, we establish a structure
theorem that gives a way to canonically decompose any rooted binary
phylogenetic network N into maximal zig-zag trails that are uniquely
determined, and use it to characterize the set of subdivision trees of N in the
form of a direct product, in a way reminiscent of the structure theorem for
finitely generated Abelian groups. From the main results, we derive a series of
linear time and linear time delay algorithms for the following problems: given
a rooted binary phylogenetic network N, 1) determine whether or not N has a
subdivision tree and find one if there exists any; 2) measure the deviation of
N from being tree-based; 3) compute the number of subdivision trees of N; 4)
list all subdivision trees of N; and 5) find a subdivision tree to maximize or
minimize a prescribed objective function. All algorithms proposed here are
optimal in terms of time complexity. Our results do not only imply and unify
various known results, but also answer many open questions and moreover enable
novel applications, such as the estimation of a maximum likelihood tree
underlying a tree-based network. The results and algorithms in this paper still
hold true for a special class of rooted non-binary phylogenetic networks.


  We formulate an age-structured nonlinear partial differential equation model
that features short-term immigration effects in a population. Individuals can
immigrate into the population as any of the three stages in the model:
susceptible, infected or recovered}. Global stability of the immigration-free
and infection-free equilibria is discussed. A generalized numerical framework
is established and specific short-term immigration scenarios are explored.


  We use the McKendrick equation with variable ageing rate and randomly
distributed maturation time to derive a state dependent distributed delay
differential equation. We show that the resulting delay differential equation
preserves non-negativity of initial conditions and we characterise local
stability of equilibria. By specifying the distribution of maturation age, we
recover state dependent discrete, uniform and gamma distributed delay
differential equations. We show how to reduce the uniform case to a system of
state dependent discrete delay equations and the gamma distributed case to a
system of ordinary differential equations. To illustrate the benefits of these
reductions, we convert previously published transit compartment models into
equivalent distributed delay differential equations.


  Among the different biomarkers of aging based on omics and clinical data, DNA
methylation clocks stand apart providing unmatched accuracy in assessing the
biological age of both humans and animal models of aging. Here, we discuss
robustness of DNA methylation clocks and bounds on their out-of-sample
performance and review computational strategies for development of the clocks.


  In here presented in silico study we suggest a way how to implement the
evolutionary principles into anti-cancer therapy design. We hypothesize that
instead of its ongoing supervised adaptation, the therapy may be constructed as
a self-sustaining evolutionary process in a dynamic fitness landscape
established implicitly by evolving cancer cells, microenvironment and the
therapy itself. For these purposes, we replace a unified therapy with the
`therapy species', which is a population of heterogeneous elementary therapies,
and propose a way how to turn the toxicity of the elementary therapy into its
fitness in a way conforming to evolutionary causation. As a result, not only
the therapies govern the evolution of different cell phenotypes, but the cells'
resistances govern the evolution of the therapies as well. We illustrate the
approach by the minimalistic ad hoc evolutionary model. Its results indicate
that the resistant cells could bias the evolution towards more toxic elementary
therapies by inhibiting the less toxic ones. As the evolutionary causation of
cancer drug resistance has been intensively studied for a few decades, we refer
to cancer as a special case to illustrate purely theoretical analysis.


  Quantifying and estimating wildlife population sizes is a foundation of
wildlife management. However, many carnivore species are cryptic, leading to
innate difficulties in estimating their populations. We evaluated the potential
for using more rigorous statistical models to estimate the populations of black
bears (Ursus americanus) in Wisconisin, and their applicability to other
furbearers such as bobcats (Lynx rufus). To estimate black bear populations, we
developed an AAH state-space model in a Bayesian framework based on Norton
(2015) that can account for variation in harvest and population demographics
over time. Our state-space model created an accurate estimate of the black bear
population in Wisconsin based on age-at-harvest data and improves on previous
models by using little demographic data, no auxiliary data, and not being
sensitive to initial population size. The increased accuracy of the AAH
state-space models should allow for better management by being able to set
accurate quotas to ensure a sustainable harvest for the black bear population
in each zone. We also evaluated the demography and annual harvest data of
bobcats in Wisconsin to determine trends in harvest, method, and hunter
participation in Wisconsin. Each trait of harvested bobcats that we tested
(mass, male:female sex ratio, and age) changed over time, and because these
were interrelated, and we inferred that harvest selection for larger size
biased harvests in favor of a) larger, b) older, and c) male bobcats by hound
hunters. We also found an increase in the proportion of bobcats that were
harvested by hound hunting compared to trapping, and that hound hunters were
more likely to create taxidermy mounts than trappers. We also found that
decreasing available bobcat tags and increasing success have occurred over
time, and correlate with substantially increasing both hunter populations and
hunter interest.


  We present the epithelial-to-mesenchymal transition (EMT) from two
perspectives: experimental/technological and theoretical. We review the state
of the current understanding of the regulatory networks that underlie EMT in
three physiological contexts: embryonic development, wound healing, and
metastasis. We describe the existing experimental systems and manipulations
used to better understand the molecular participants and factors that influence
EMT and metastasis. We review the mathematical models of the regulatory
networks involved in EMT, with a particular emphasis on the network motifs
(such as coupled feedback loops) that can generate intermediate hybrid states
between the epithelial and mesenchymal states. Ultimately, the understanding
gained about these networks should be translated into methods to control
phenotypic outcomes, especially in the context of cancer therapeutic
strategies. We present emerging theories of how to drive the dynamics of a
network toward a desired dynamical attractor (e.g. an epithelial cell state)
and emerging synthetic biology technologies to monitor and control the state of
cells.


  We introduce a simple physical picture to explain the process of molecular
sorting, whereby specific proteins are concentrated and distilled into
submicrometric lipid vesicles in eukaryotic cells. To this purpose, we
formulate a model based on the coupling of spontaneous molecular aggregation
with vesicle nucleation. Its implications are studied by means of a
phenomenological theory describing the diffusion of molecules towards multiple
sorting centers that grow due to molecule absorption and are extracted when
they reach a sufficiently large size. The predictions of the theory are
compared with numerical simulations of a lattice-gas realization of the model
and with experimental observations. The efficiency of the distillation process
is found to be optimal for intermediate aggregation rates, where the density of
sorted molecules is minimal and the process obeys simple scaling laws.
Quantitative measures of endocytic sorting performed in primary endothelial
cells are compatible with the hypothesis that these optimal conditions are
realized in living cells.


  Cyclic predator-prey systems have been shown to give rise to rich, and novel,
space-time patterns, as for example coarsening domains with non-trivial
in-domain dynamics. In this work we study numerically the responses of a cyclic
six-species model, characterized by the formation of spirals inside coarsening
domains, to two different types of perturbations: changing the values of the
predation and reproduction rates as well as changing the interaction scheme.
For both protocols we monitor the time evolution of the system after the onset
of the perturbation through the measurement of dynamical correlation functions
and time-dependent densities of empty sites. In this way we gain insights into
the complex responses to different perturbations in a system where spirals,
which are due to the formation of cyclic alliances, dominate the dynamics
inside the coarsening domains.


  We present a novel approach for the prediction of anticancer compound
sensitivity by means of multi-modal attention-based neural networks (PaccMann).
In our approach, we integrate three key pillars of drug sensitivity, namely,
the molecular structure of compounds, transcriptomic profiles of cancer cells
as well as prior knowledge about interactions among proteins within cells. Our
models ingest a drug-cell pair consisting of SMILES encoding of a compound and
the gene expression profile of a cancer cell and predicts an IC50 sensitivity
value. Gene expression profiles are encoded using an attention-based encoding
mechanism that assigns high weights to the most informative genes. We present
and study three encoders for SMILES string of compounds: 1) bidirectional
recurrent 2) convolutional 3) attention-based encoders. We compare our devised
models against a baseline model that ingests engineered fingerprints to
represent the molecular structure. We demonstrate that using our
attention-based encoders, we can surpass the baseline model. The use of
attention-based encoders enhance interpretability and enable us to identify
genes, bonds and atoms that were used by the network to make a prediction.


  Whether or not cooperation is favored in evolutionary games on graphs depends
on the population structure and spatial properties of the interaction network.
Population structures can be expressed as configurations. Such configurations
extend scenarios with a single cooperator among defectors to any number of
cooperators and any arrangement of cooperators and defectors. Thus, as a single
cooperator can be interpreted as a lone mutant, the discussion about fixation
properties based on configurations also applies to multiple mutants. For
interaction networks modeled as regular graphs and for weak selection, the
emergence of cooperation can be assessed by structure coefficients, which are
specific for a configuration and a graph. We analyze these structure
coefficients and particularly show that under certain conditions the
coefficients strongly correlate to the average shortest path length between
cooperators on the evolutionary graph. Thus,for multiple cooperators fixation
properties on regular evolutionary graphs can be linked to cooperator path
length.


  In 2001 Allen and Steel showed that, if subtree and chain reduction rules
have been applied to two unrooted phylogenetic trees, the reduced trees will
have at most 28k taxa where k is the TBR (Tree Bisection and Reconnection)
distance between the two trees. Here we reanalyse Allen and Steel's
kernelization algorithm and prove that the reduced instances will in fact have
at most 15k-9 taxa. Moreover we show, by describing a family of instances which
have exactly 15k-9 taxa after reduction, that this new bound is tight. These
instances also have no common clusters, showing that a third
commonly-encountered reduction rule, the cluster reduction, cannot further
reduce the size of the kernel in the worst case. To achieve these results we
introduce and use "unrooted generators" which are analogues of rooted
structures that have appeared earlier in the phylogenetic networks literature.
Using similar argumentation we show that, for the minimum hybridization problem
on two rooted trees, 9k-2 is a tight bound (when subtree and chain reduction
rules have been applied) and 9k-4 is a tight bound (when, additionally, the
cluster reduction has been applied) on the number of taxa, where k is the
hybridization number of the two trees.


  Drug resistance is still a major challenge in cancer therapy. Drug
combination is expected to overcome drug resistance. However, the number of
possible drug combinations is enormous, and thus it is infeasible to
experimentally screen all effective drug combinations considering the limited
resources. Therefore, computational models to predict and prioritize effective
drug combinations is important for combinatory therapy discovery in cancer. In
this study, we proposed a novel deep learning model, AuDNNsynergy, to
prediction drug combinations by integrating multi-omics data and chemical
structure data. In specific, three autoencoders were trained using the gene
expression, copy number and genetic mutation data of all tumor samples from The
Cancer Genome Atlas. Then the physicochemical properties of drugs combined with
the output of the three autoencoders, characterizing the individual cancer
cell-lines, were used as the input of a deep neural network that predicts the
synergy value of given pair-wise drug combinations against the specific cancer
cell-lines. The comparison results showed the proposed AuDNNsynergy model
outperforms four state-of-art approaches, namely DeepSynergy, Gradient Boosting
Machines, Random Forests, and Elastic Nets. Moreover, we conducted the
interpretation analysis of the deep learning model to investigate potential
vital genetic predictors and the underlying mechanism of synergistic drug
combinations on specific cancer cell-lines.


  In this manuscript we investigate the long-term behavior of a single-species
fishery, which is harvested by several fleets. The time evolution of this
population is modeled by a discrete time stochastic age-structured model. We
assume that incertitude only affects the recruitment. First, for the
deterministic version of this model, we characterize the equilibrium yields in
terms of the fishing mortality. Then, for the stochastic version, we introduce
the concepts of maximum expected, log expected and harmonic expected
sustainable yield, and we analyze how the incertitude affects the behavior of
these yields and their stationary distribution. All the numerical simulations
are performed with data obtained from Patagonian Tooth-fish fishery, which is
harvested by four different type of fleets: Chilean Industrial fleet, Chilean
Artisanal fleet, Argentinean longline fleet, and Argentinean Artisanal fleet.


  Taking advantages of high-throughput genotyping technology of single
nucleotide polymorphism (SNP), large genome-wide association studies (GWASs)
have been considered as the promise to unravel the complex relationships
between genotypes and phenotypes, in particularly common diseases. However,
current multi-locus-based methods are insufficient, in terms of computational
cost and discrimination power, to detect statistically significant interactions
and they are lacking in the ability of finding diverse genetic effects on
multifarious diseases. Especially, multiple statistic tests for high-order
epistasis ($ \geq $ 2 SNPs) will raise huge analytical challenges because the
computational cost increases exponentially as the growth of the cardinality of
SNPs in an epistatic module. In this paper, we develop a simple, fast and
powerful method, named JS-MA, using the Jensen-Shannon divergence and a
high-dimensional $ k $-mean clustering algorithm for mapping the genome-wide
multi-locus epistatic interactions on multiple diseases. Compared with some
state-of-the-art association mapping tools, our method is demonstrated to be
more powerful and efficient from the experimental results on the systematical
simulations. We also applied JS-MA to the GWAS datasets from WTCCC for two
common diseases, i.e. Rheumatoid Arthritis and Type 1 Diabetes. JS-MA not only
confirms some recently reported biologically meaningful associations but also
identifies some novel findings. Therefore, we believe that our method is
suitable and efficient for the full-scale analysis of multi-disease-related
interactions in the large GWASs.


  Microscopic models of reaction-diffusion processes on the cell membrane can
link local spatiotemporal effects to macroscopic self-organized patterns often
observed on the membrane. Simulation schemes based on the microscopic lattice
method (MLM) can model these processes at the microscopic scale by tracking
individual molecules, represented as hard-spheres, on fine lattice voxels.
Although MLM is simple to implement and is generally less computationally
demanding than off-lattice approaches, its accuracy and consistency in modeling
surface reactions have not been fully verifed. Using the Spatiocyte scheme, we
study the accuracy of MLM in diffusion-influenced surface reactions. We derive
the lattice-based bimolecular association rates for two-dimensional
surface-surface reaction and one-dimensional volume-surface adsorption
according to the Smoluchowski-Collins-Kimball model and random walk theory. We
match the time-dependent rates on lattice with off-lattice counterparts to
obtain the correct expressions for MLM parameters in terms of physical
constants. The expressions indicate that the voxel size needs to be at least
0.6% larger than the molecule to accurately simulate surface reactions on
triangular lattice. On square lattice, the minimum voxel size should be even
larger, at 5%. We also demonstrate the ability of MLM-based schemes such as
Spatiocyte to simulate a reaction-diffusion model that involves all dimensions:
three-dimensional diffusion in the cytoplasm, two-dimensional diffusion on the
cell membrane and one-dimensional cytoplasm-membrane adsorption. With the
model, we examine the contribution of the 2D reaction pathway to the overall
reaction rate at different reactant diffusivity, reactivity and concentrations.


  The sampling of interaction partners depends on often implicit modelling
assumptions, yet has marked effects on the dynamics in evolutionary games. One
particularly important aspect is whether or not competitors also interact.
Population structures naturally affect sampling such that in a microscopic
interpretation of the replicator dynamics in well-mixed populations competing
individuals do not interact but do interact in structured populations. In
social dilemmas interactions with competitors invariably inhibit cooperation
while limited local interactions in structured populations support cooperation
by reducing exploitation through cluster formation. These antagonistic effects
of population structures on cooperation affect interpretations and the
conclusions depend on the details of the comparison. For example, in the
Snowdrift game, spatial structure may inhibit cooperation when compared to the
replicator dynamics. However, modifying the replicator dynamics to include
interactions between competitors lowers the equilibrium frequency of
cooperators, which changes the conclusions and space is invariably beneficial,
just as in the Prisoner's Dilemma. These conclusions are confirmed by
comparisons with random matching models, which mimic population structures but
randomly reshuffle individuals to inhibit spatial correlations. Finally, the
differences in the dynamics with and without interactions among competing
individuals underlies the differences between death-birth and birth-death
updating in the spatial Moran process: death-birth supports cooperation because
competitors tend not to interact whereas they tend to do for birth-death
updating and hence cooperators provide direct support to competitors to their
own detriment.


  It is important to study on genes affecting to growing environment of main
crops. Especially the recognition problem of promoter region, which is the
problem to predict whether DNA sequences contain promoter regions or not, is
prior to find abiotic stress-inducible genes. Studies on predicting promoter
sequences in DNA sequences have been studied by traditional pattern matching
methods and machine learning methods in biology and computer science.


  Gene expression and its regulation is a nonequilibrium stochastic process.
Different molecules are involved in several biochemical steps in this process
with low copies. It is observed that the stochasticity in biochemical processes
is mainly due to the low copy number of the molecules present in the system.
Several studies also show that the nonequilibrium biochemical processes require
energy cost. But cellular system has developed itself through natural evolution
by minimizing energy cost for optimum output. Here we study the role of
stochasticity qualitatively in a network of two genes using stochastic
simulation method and approximately measure the energy consumption for the gene
expression process. We find that the noise in gene expression process reduces
the energy cost of protein synthesis. Therefore, we argued that the
stochasticity in gene expression may be a choice of cellular system for protein
synthesis with minimum energy cost.


  Although Darwinian models are rampant in the social sciences, social
scientists do not face the problem that motivated Darwin's theory of natural
selection: the problem of explaining how lineages evolve despite that any
traits they acquire are regularly discarded at the end of the lifetime of the
individuals that acquired them. While the rationale for framing culture as an
evolutionary process is correct, it does not follow that culture is a Darwinian
or selectionist process, or that population genetics provides viable starting
points for modeling cultural change. This paper lays out step-by-step arguments
as to why a selectionist approach to cultural evolution is inappropriate,
focusing on the lack of randomness, and lack of a self-assembly code. It
summarizes an alternative evolutionary approach to culture: self-other
reorganization via context-driven actualization of potential.


  One of the main properties of biological systems is modularity, which
manifests itself at all levels of their organization, starting with the level
of molecular genetics, ending with the level of whole organisms and their
communities. In a simplified form, these basic principles were transferred from
the genetics of populations to the field of evolutionary computations, in order
to solve applied optimization problems. Over almost half a century of
development in this field of computer science, considerable practical
experience has been gained and interesting theoretical results have been
obtained. In this survey, the phenomena and patterns associated with modularity
in genetics and evolutionary computations are compared. An analysis of
similarities and differences in the results obtained in these areas is carried
out from the modularity view point. The possibilities for knowledge transfer
between the areas are discussed.


  The evolutionary and ecological processes behind the origin of species are
among the most fundamental problems in biology. In fact, many theoretical
hypothesis on different type of speciation have been proposed. In particular,
models of sympatric speciation leading to the formation of new species without
geographical isolation, are based on the niche hypothesis: the diversification
of the population is induced by the competition for a limited set of the
available resources. On the other hand, neutral models of evolution have shown
that stochastic forces are sufficient to generate coexistence of different
species. In this work, we bring this dichotomy to the context of species
formation, and we study how neutral and niche forces contribute to sympatric
speciation in a model ecosystem. In particular, we study the evolution of a
population of individuals with asexual reproduction whose inherited characters
or phenotypes are specified by both niche-based and neutral traits. We analyse
the stationary state of the dynamics, and study the distribution of individuals
in the whole space of possible phenotypes. We show, both by numerical
simulations and analytics, that there is a non-trivial coupling between neutral
and niche forces induced by stochastic effects in the evolution of the
population that allows the formation of clusters (i.e., species) in the
phenotypic space. Our framework can be generalised also to sexual reproduction
or other type of population dynamics.


  The ribosome flow model with input and output (RFMIO) is a deterministic
dynamical system that has been used to study the flow of ribosomes during mRNA
translation. The RFMIO and its variants encapsulate important properties that
are relevant to modeling ribosome flow such as the possible evolution of
"traffic jams" and non-homogeneous elongation rates along the mRNA molecule,
and can also be used for studying additional intracellular processes such as
transcription, transport, and more. Here we consider networks of interconnected
RFMIOs as a fundamental tool for modeling, analyzing and re-engineering the
complex mechanisms of protein production. In these networks, the output of each
RFMIO may be divided, using connection weights, between several inputs of other
RFMIOs. We show that under quite general feedback connections the network has
two important properties: (1) it admits a unique steady-state and every
trajectory converges to this steady-state, and (2) the problem of how to
determine the connection weights so that the network steady-state output is
maximized is a convex optimization problem. These mathematical properties make
these networks highly suitable as models of various phenomena: property (1)
means that the behavior is predictable and ordered, and property (2) means that
determining the optimal weights is numerically tractable even for large-scale
networks. For the specific case of a feed-forward network of RFMIOs we prove an
additional useful property, namely, that there exists a spectral representation
for the network steady-state, and thus it can be determined without any
numerical simulations of the dynamics. We describe the implications of these
result to several fundamental biological phenomena and biotechnological
objectives.


  Heart failure is a major cause for premature death. Given heterogeneity of
the heart failure syndrome, identifying genetic determinants of cardiac
function and structure may provide greater insights into heart failure. Despite
progress in understanding the genetic basis of heart failure through genome
wide association studies, heritability of heart failure is not well understood.
Gaining further insights into mechanisms that contribute to heart failure
requires systematic approaches that go beyond single trait analysis. We
integrated Bayesian multi-trait approach and Bayesian networks for the analysis
of 10 correlated traits of cardiac structure and function measured for 3387
individuals with whole exome sequence data. While using single-trait based
approaches did not find any significant genetic variant, applying the
integrative Bayesian multi-trait approach, we identified 3 novel variants
located in genes, RGS3, CHD3, and MRPL38 with significant impact on the cardiac
traits such as left ventricular volume index, parasternal long axis
interventricular septum thickness, and mean left ventricular wall thickness.
Among these, the rare variant NC_000009.11:g.116346115C>A (rs144636307) in RGS3
showed pleiotropic effect on left ventricular mass index, left ventricular
volume index and Maximum left atrial anterior-posterior diameter while RGS3 can
inhibit TGF-beta signaling associated with left ventricle dilation and systolic
dysfunction.


  Development of fast methods to conduct in silico experiments using
computational models of cellular signaling is a promising approach toward
advances in personalized medicine. However, software-based cellular network
simulation has run-times plagued by wasted CPU cycles and unnecessary
processes. Hardware-based simulation affords substantial speedup, but prior
attempts at hardware-based biological simulation have been limited in scope and
have suffered from inaccuracies due to poor random number generation. In this
work, we propose several hardware-based simulation schemes utilizing novel
random update index generation techniques for step-based and round-based
stochastic simulations of cellular networks. Our results show improved runtimes
while maintaining simulation accuracy compared to software implementations.


  We present a novel mechanism for resolving the mechanical rigidity of
nanoscopic circular polymers that flow in a complex environment. The emergence
of a regime of negative differential mobility induced by topological
interactions between the rings and the substrate is the key mechanism for
selective sieving of circular polymers with distinct flexibilities. A simple
model accurately describes the sieving process observed in molecular dynamics
simulations and yields experimentally verifiable analytical predictions, which
can be used as a reference guide for improving filtration procedures of
circular filaments. The topological sieving mechanism we propose ought to be
relevant also in probing the microscopic details of complex substrates.


  Sequential data is being generated at an unprecedented pace in various forms,
including text and genomic data. This creates the need for efficient
compression mechanisms to enable better storage, transmission and processing of
such data. To solve this problem, many of the existing compressors attempt to
learn models for the data and perform prediction-based compression. Since
neural networks are known as universal function approximators with the
capability to learn arbitrarily complex mappings, and in practice show
excellent performance in prediction tasks, we explore and devise methods to
compress sequential data using neural network predictors. We combine recurrent
neural network predictors with an arithmetic coder and losslessly compress a
variety of synthetic, text and genomic datasets. The proposed compressor
outperforms Gzip on the real datasets and achieves near-optimal compression for
the synthetic datasets. The results also help understand why and where neural
networks are good alternatives for traditional finite context models


  With the advent of the big data era, generative models of complex networks
are becoming elusive from direct computational simulation. We present an exact,
linear-algebraic reduction scheme of generative models of networks. By
exploiting the bilinear structure of the matrix representation of the
generative model, we separate its null eigenspace, and reduce the exact
description of the generative model to a smaller vector space. After reduction,
we group generative models in universality classes according to their rank and
metric signature, and work out, in a computationally affordable way, their
relevant properties (e.g., spectrum). The reduction also provides the
environment for a simplified computation of their properties. The proposed
scheme works for any generative model admitting a matrix representation, and
will be very useful in the study of dynamical processes on networks, as well as
in the understanding of generative models to come, according to the provided
classification.


  Under mass-action kinetics, biochemical reaction networks give rise to
polynomial autonomous dynamical systems whose parameters are often difficult to
estimate. We deal in this paper with the problem of identifying the kinetic
parameters of a class of biochemical networks which are abundant, such as
multisite phosphorylation systems and phosphorylation cascades (for example,
MAPK cascades). For any system of this class we explicitly exhibit a single
species for each connected component of the associated digraph such that the
successive total derivatives of its concentration allow us to identify all the
parameters occurring in the component. The number of derivatives needed is
bounded essentially by the length of the corresponding connected component of
the digraph. Moreover, in the particular case of the cascades, we show that the
parameters can be identified from a bounded number of successive derivatives of
the last product of the last layer. This theoretical result induces also a
heuristic interpolation-based identifiability procedure to recover the values
of the rate constants from exact measurements.


  Recent experimental results have shown that active enzymes can diffuse faster
when they are in the presence of their substrates. Fluorescence correlation
spectroscopy (FCS), which relies on analyzing the fluctuations in fluorescence
intensity signal to measure the diffusion coefficient of particles, has
typically been employed in most of the prior studies. However, flaws in the FCS
method, due to its high sensitivity to the environment, have recently been
evaluated, calling the prior diffusion results into question. It behooves us to
adopt complimentary and direct methods to measure the mobility of enzymes in
solution. Herein, we use a novel technique of direct single-molecule imaging to
observe the diffusion of single enzymes. This technique is less sensitive to
intensity fluctuations and gives the diffusion coefficient directly based on
the trajectory of the enzymes. Our measurements recapitulate that enzyme
diffusion is enhanced in the presence of its substrate and find that the
relative increase in diffusion of a single enzyme is even higher than those
previously reported using FCS. We also use this complementary method to test if
the total enzyme concentration affects the relative increase in diffusion and
if enzyme oligomerization state changes during catalytic turnover. We find that
the diffusion increase is independent of the total background concentration of
enzyme and the catalysis of substrate does not change the oligomerization state
of enzymes.


  Cooperation is ubiquitous across all levels of biological systems ranging
from microbial communities to human societies. It, however, seemingly
contradicts the evolutionary theory, since cooperators are exploited by
free-riders and thus are disfavored by natural selection. Many studies based on
evolutionary game theory have tried to solve the puzzle and figure out the
reason why cooperation exists and how it emerges. Network reciprocity is one of
the mechanisms to promote cooperation, where nodes refer to individuals and
links refer to social relationships. The spatial arrangement of mutant
individuals, which refers to the clustering of mutants, plays a key role in
network reciprocity. Besides, many other mechanisms supporting cooperation
suggest that the clustering of mutants plays an important role in the expansion
of mutants. However, the clustering of mutants and the game dynamics are
typically coupled. It is still unclear how the clustering of mutants alone
alters the evolutionary dynamics. To this end, we employ a minimal model with
frequency independent fitness on a circle. It disentangles the clustering of
mutants from game dynamics. The distance between two mutants on the circle is
adopted as a natural indicator for the clustering of mutants or assortment. We
find that the assortment is an amplifier of the selection for the connected
mutants compared with the separated ones. Nevertheless, as mutants are
separated, the more dispersed mutants are, the greater the chance of invasion
is. It gives rise to the non-monotonic effect of clustering, which is
counterintuitive. On the other hand, we find that less assortative mutants
speed up fixation. Our model shows that the clustering of mutants plays a
non-trivial role in fixation, which has emerged even if the game interaction is
absent.


  Use of a combination of statistical thermodynamics and the Gershgorin theorem
enable us to guess, in the thermodynamic limit, a plausible value for the upper
bound free-energy difference between native-like structures of monomeric
globular proteins. Support to our result in light of both the observed
free-energy change between the native and denatured states and the
microstability free-energy values obtained from the observed micro-unfolding
tendency of nine globular proteins, will be here discussed.


  We built a two-state model of an asexually reproducing organism in a periodic
environment endowed with the capability to anticipate an upcoming environmental
change and undergo pre-emptive switching. By virtue of these anticipatory
transitions, the organism oscillates between its two states that is a time
$\theta$ out of sync with the environmental oscillation. We show that an
anticipation-capable organism increases its long-term fitness over an organism
that oscillates in-sync with the environment, provided $\theta$ does not exceed
a threshold. We also show that the long-term fitness is maximized for an
optimal anticipation time that decreases approximately as $1/n$, $n$ being the
number of cell divisions in time $T$. Furthermore, we demonstrate that optimal
"anticipators" outperforms "bet-hedgers" in the range of parameters considered.
For a sub-optimal ensemble of anticipators, anticipation performs better to
bet-hedging only when the variance in anticipation is small compared to the
mean and the rate of pre-emptive transition is high. Taken together, our work
suggests that anticipation increases overall fitness of an organism in a
periodic environment and it is a viable alternative to bet-hedging provided the
error in anticipation is small.


  Delays are an important phenomenon arising in a wide variety of real world
systems. They occur in biological models because of diffusion effects or as
simplifying modeling elements. We propose here to consider delayed stochastic
reaction networks. The difficulty here lies in the fact that the state-space of
a delayed reaction network is infinite-dimensional, which makes their analysis
more involved. We demonstrate here that a particular class of stochastic
time-varying delays, namely those that follow a phase-type distribution, can be
exactly implemented in terms of a chemical reaction network. Hence, any
delay-free network can be augmented to incorporate those delays through the
addition of delay-species and delay-reactions. Hence, for this class of
stochastic delays, which can be used to approximate any delay distribution
arbitrarily accurately, the state-space remains finite-dimensional and,
therefore, standard tools developed for standard reaction network still apply.
In particular, we demonstrate that for unimolecular mass-action reaction
networks that the delayed stochastic reaction network is ergodic if and only if
the non-delayed network is ergodic as well. Bimolecular reactions are more
difficult to consider but an analogous result is also obtained. These results
tell us that delays that are phase-type distributed, regardless of their
distribution, are not harmful to the ergodicity property of reaction networks.
We also prove that the presence of those delays adds convolution terms in the
moment equation but does not change the value of the stationary means compared
to the delay-free case. Finally, the control of a certain class of delayed
stochastic reaction network using a delayed antithetic integral controller is
considered. It is proven that this controller achieves its goal provided that
the delay-free network satisfy the conditions of ergodicity and
output-controllability.


  Inspired by recent experiments on chromosomal dynamics, we introduce an
exactly solvable model for the interaction between a flexible polymer and a set
of motor-like enzymes. The enzymes can bind and unbind to specific sites of the
polymer and when bound produce a dipolar force on two neighboring monomers. We
study the resulting non-equilibrium dynamics of the polymer and find that the
motion of the monomers has several properties that were observed experimentally
for chromosomal loci: a subdiffusive mean squared displacement and the
appearance of regions of correlated motion. We also determine the velocity
autocorrelation of the monomers and find that the underlying stochastic process
is not fractional Brownian motion. Finally, we show that the active forces
swell the polymer by an amount that becomes constant for large polymers.


  Prediction of metabolism in cytochrome P450s remains to be a crucial yet
challenging topic in discovering and designing drugs, agrochemicals and
nutritional supplements. The problem is challenging because the rate of P450
metabolism depends upon both the intrinsic chemical reactivity of the site and
the protein-ligand geometry that is energetically accessible in the active site
of a given P450 isozyme. We have addressed this problem using a two-level
screening system. The first level implements an empirical QSAR-based scoring
function employing the local chemical motifs to characterize the intrinsic
reactivity. The second level uses molecular docking and molecular mechanics to
account for the geometrical effects, including induced-fit effects in the
protein which can be very important in P450 interactions with ligands. This
approach has achieved high accuracy for both the P450 3A4 and 2D6 isoforms. In
identifying at least one metabolic site in the top two ranked positions, the
prediction rate can reach as high as 92.7% for the test set of isoform 3A4. For
the 2D6 isoform, 100% accuracy is achieved on this basic evaluation metric,
and, because this active site is considerably smaller and more selective than
3A4, very high precision is attained for full prediction of all metabolic
sites. The method also requires considerably less CPU time than our previous
efforts, which involved a large number of expensive simulations for each ligand
to be evaluated. After screening using the empirical score function, only a few
best candidates are left for each ligand, making the number of necessary
estimations in the second level very small, which significantly reduces the
computation time.


  The biological processes that execute complex multiple functions, such as
cell cycle, must ensure the order of sequential events and keep the dynamic
robustness against various fluctuations. Here, we examine the dynamic mechanism
and the fundamental structure to achieve these properties in the cell-cycle
process of budding yeast Saccharomyces cerevisiae. We show that the budding
yeast cell-cycle process behaves like an excitable system containing three
well-coupled saddle-node bifurcations to execute DNA replication and mitosis
events. The yeast cell-cycle regulatory network can be separated into G1/S
phase module, early M module and late M phase module, where the positive
feedbacks in each module and the interactions among the modules play important
role. If the cell-cycle process operates near the critical points of the
saddle-node bifurcations, there is a critical slowing down or ghost effect.
This can provide the cell-cycle process with a sufficient duration for each
event and an attractive manifold for the state checking of the completion of
DNA replication and mitosis; moreover, the fluctuation in the early
module/event is forbidden to transmit to the latter module/event. Our results
suggest both a fundamental structure of cell-cycle regulatory network and a
hint for the evolution of eukaryotic cell-cycle processes, from the dynamic
checking mechanism to the molecule checkpoint pathway.


  In spite of decades of research, much remains to be discovered about folding:
the detailed structure of the initial (unfolded) state, vestigial folding
instructions remaining only in the unfolded state, the interaction of the
molecule with the solvent, instantaneous power at each point within the
molecule during folding, the fact that the process is stable in spite of myriad
possible disturbances, potential stabilization of trajectory by chaos, and, of
course, the exact physical mechanism (code or instructions) by which the
folding process is specified in the amino acid sequence. Simulations based upon
microscopic physics have had some spectacular successes and continue to
improve, particularly as super-computer capabilities increase. The simulations,
exciting as they are, are still too slow and expensive to deal with the
enormous number of molecules of interest. In this paper, we introduce an
approximate model based upon physics, empirics, and information science which
is proposed for use in machine learning applications in which very large
numbers of sub-simulations must be made. In particular, we focus upon machine
learning applications in the learning phase and argue that our model is
sufficiently close to the physics that, in spite of its approximate nature, can
facilitate stepping through machine learning solutions to explore the mechanics
of folding mentioned above. We particularly emphasize the exploration of energy
flow (power) within the molecule during folding, the possibility of energy
scale invariance (above a threshold), vestigial information in the unfolded
state as attractive targets for such machine language analysis, and statistical
analysis of an ensemble of folding micro-steps.


  The actin cytoskeleton is composed of multiple networks which are specialized
for several processes such as cell motility or cell division. Each of these
networks are composed of organized actin microfilaments which are decorated
with specific sets of actin binding proteins (ABPs). The molecular mechanisms
guiding ABPs to specific actin networks are still poorly understood, but
cooperativity, the mechanism by which the binding of an ABP is positively
influenced by proximal bound ABPs, plays a crucial role in generating locally
dense stretches of ABPs. Cooperative binding is characterized by its amplitude,
but also by the range at which its effects are propagated along an actin
filament through long-range allosteric interactions. The range of these
allosteric effects is still debated, but is likely to be significant at the
lengthscale of actin filaments in cells. Here, we investigated how
cooperativity influences the clustering of ABPs, using a stochastic
computational model of binding of ABPs to actin filaments. The model reproduces
the formation of ABP clusters observed experimentally at the single filament
scale, and provides a theoretical estimation of the range of cooperativity for
proteins such as ADF/cofilin. We found that both the amplitude and the spatial
range of cooperativity dramatically impact the properties of clustering.
However, the parameters of cooperativity modulate differently the rate of
assembly, size and dynamics of the ABP clusters, suggesting that cooperativity
is an efficient mechanism to regulate precisely the recruitment of ABPs in
cells. This work provides a more general framework for future understanding of
how actin networks acquire distinct and specific protein compositions from a
common cytoplasm.


  We study fixation probabilities for the Moran stochastic process for the
evolution of a population with three or more types of individuals and
frequency-dependent fitnesses. Contrarily to the case of populations with two
types of individuals, in which fixation probabilities may be calculated by an
exact formula, here we must solve a large system of linear equations. We first
show that this system always has a unique solution. Other results are upper and
lower bounds for the fixation probabilities obtained by coupling the Moran
process with three strategies with birth-death processes with only two
strategies. We also apply our bounds to the problem of evolution of cooperation
in a population with three types of individuals already studied in a
deterministic setting by N\'u\~nez Rodr\'iguez and Neves (J. Math. Biol. (2016)
73:1665--1690). We argue that cooperators will be fixated in the population
with probability arbitrarily close to 1 for a large region of initial
conditions and large enough population sizes.


  The ideas that we forge creatively as individuals and groups build on one
another in a manner that is cumulative and adaptive, forming open-ended
lineages across space and time. Thus, human culture is believed to evolve. The
pervasiveness of cross-domain creativity--as when a song inspires a
painting--would appear indicative of discontinuities in cultural lineages.
However, if what evolves through culture is our worldviews--the webs of
thoughts, ideas, and attitudes that constitutes our way of seeing being in the
world--then the problem of discontinuities is solved. The state of a worldview
can be affected by information assimilated in one domain, and this
change-of-state can be expressed in another domain. In this view, the gesture,
narrative, or artifact that constitutes a specific creative act is not what is
evolving; it is merely the external manifestation of the state of an evolving
worldview. Like any evolutionary process, cultural evolution requires a balance
between novelty, via the generation of variation, and continuity, via the
preservation of variants that are adaptive. In cultural evolution, novelty is
generated through creativity, and continuity is provided by social learning
processes, e.g., imitation. Both the generative and imitative aspects of
cultural evolution are affected by social media. We discuss the trajectory from
social ideation to social innovation, focusing on the role of
self-organization, renewal, and perspective-taking at the individual and social
group level.


  Understanding the three-dimensional (3D) structure of the genome is essential
for elucidating vital biological processes and their links to human disease. To
determine how the genome folds within the nucleus, chromosome conformation
capture methods such as HiC have recently been employed. However, computational
methods that exploit the resulting high-throughput, high-resolution data are
still suffering from important limitations. In this work, we explore the idea
of manifold learning for the 3D chromatin structure inference and present a
novel method, REcurrent Autoencoders for CHromatin 3D structure prediction
(REACH-3D). Our framework employs autoencoders with recurrent neural units to
reconstruct the chromatin structure. In comparison to existing methods,
REACH-3D makes no transfer function assumption and permits dynamic analysis.
Evaluating REACH-3D on synthetic data indicated high agreement with the ground
truth. When tested on real experimental HiC data, REACH-3D recovered most
faithfully the expected biological properties and obtained the highest
correlation coefficient with microscopy measurements. Last, REACH-3D was
applied to dynamic HiC data, where it successfully modeled chromatin
conformation during the cell cycle.


  De novo design seeks to generate molecules with required property profiles by
virtual design-make-test cycles. With the emergence of deep learning and neural
generative models in many application areas, models for molecular design based
on neural networks appeared recently and show promising results. However, the
new models have not been profiled on consistent tasks, and comparative studies
to well-established algorithms have only seldom been performed.
  To standardize the assessment of both classical and neural models for de novo
molecular design, we propose an evaluation framework, GuacaMol, based on a
suite of standardized benchmarks. The benchmark tasks encompass measuring the
fidelity of the models to reproduce the property distribution of the training
sets, the ability to generate novel molecules, the exploration and exploitation
of chemical space, and a variety of single and multi-objective optimization
tasks. The benchmarking open-source Python code, and a leaderboard can be found
on https://benevolent.ai/guacamol


  Complex systems may contain heterogeneous types of variables that interact in
a multi-level and multi-scale manner. In this context, high-level layers may
considered as groups of variables interacting in lower-level layers. This is
particularly true in biology, where, for example, genes are grouped in pathways
and two types of interactions are present: pathway-pathway interactions and
gene-gene interactions. However, from data it is only possible to measure the
expression of genes while it is impossible to directly measure the activity of
pathways. Nevertheless, the knowledge on the inter-dependence between the
groups and the variables allows for a multi-layer network inference, on both
observed variables and groups, even if no direct information on the latter is
present in the data (hence groups are considered as latent). In this paper, we
propose an extension of the latent graphical lasso method that leverages on the
knowledge of the inter-links between the hidden (groups) and observed layers.
The method exploits the knowledge of group structure that influence the
behaviour of observed variables to retrieve a two layers network. Its efficacy
was tested on synthetic data to check its ability in retrieving the network
structure compared to the ground truth. We present a case study on
Neuroblastoma, which shows how our multi-level inference is relevant in real
contexts to infer biologically meaningful connections.


  Research on kinase-targeting drugs has made great strides over the last 30
years and is attracting greater attention for the treatment of yet more
kinase-related diseases. Currently, 42 kinase drugs have been approved by the
FDA, most of which (39) are Type I/II inhibitors. Notwithstanding these
advances, it is desirable to target additional kinases for drug development as
more than 200 diseases, particularly cancers, are directly associated with
aberrant kinase regulation and signaling. Here, we review the extant Type I/II
drugs systematically to obtain insights into the binding pocket
characteristics, the associated features of Type I/II drugs, and the mechanism
of action to facilitate future kinase drug design and discovery. We conclude by
summarizing the main successes and limitations of targeting kinase for the
development of drugs.