  Organisms show a remarkable range of sizes, yet the dimensions of a single
cell rarely exceed $100$ $\mu$m. While the physical and biological origins of
this constraint remain poorly understood, exceptions to this rule give valuable
insights. A well-known counterexample is the aquatic plant $Chara$, whose cells
can exceed $10$ cm in length and $1$ mm in diameter. Two spiraling bands of
molecular motors at the cell periphery drive the cellular fluid up and down at
speeds up to $100$ $\mu$m/s, motion that has been hypothesized to mitigate the
slowness of metabolite transport on these scales and to aid in homeostasis.
This is the most organized instance of a broad class of continuous motions
known as "cytoplasmic streaming", found in a wide range of eukaryotic organisms
- algae, plants, amoebae, nematodes, and flies - often in unusually large
cells. In this overview of the physics of this phenomenon, we examine the
interplay between streaming, transport and cell size, and discuss the possible
role of self-organization phenomena in establishing the observed patterns of
streaming.


  The emergence and maintenance of cooperative behavior is a fascinating topic
in evolutionary biology and social science. The public goods game (PGG) is a
paradigm for exploring cooperative behavior. In PGG, the total resulting payoff
is divided equally among all participants. This feature still leads to the
dominance of defection without substantially magnifying the public good by a
multiplying factor. Much effort has been made to explain the evolution of
cooperative strategies, including a recent model in which only a portion of the
total benefit is shared by all the players through introducing a new strategy
named persistent cooperation. A persistent cooperator is a contributor who is
willing to pay a second cost to retrieve the remaining portion of the payoff
contributed by themselves. In a previous study, this model was analyzed in the
framework of well-mixed populations. This paper focuses on discussing the
persistent cooperation in lattice-structured populations. The evolutionary
dynamics of the structured populations consisting of three types of competing
players (pure cooperators, defectors and persistent cooperators) are revealed
by theoretical analysis and numerical simulations. In particular, the
approximate expressions of fixation probabilities for strategies are derived on
one-dimensional lattices. The phase diagrams of stationary states, the
evolution of frequencies and spatial patterns for strategies are illustrated on
both one-dimensional and square lattices by simulations. Our results are
consistent with the general observation that, at least in most situations, a
structured population facilitates the evolution of cooperation. Specifically,
here we find that the existence of persistent cooperators greatly suppresses
the spreading of defectors under more relaxed conditions in structured
populations compared to that obtained in well-mixed population.


  Existing elastic network models are typically parametrized at a given cutoff
distance and often fail to properly predict the thermal fluctuation of many
macromolecules that involve multiple characteristic length scales. We introduce
a multiscale flexibility-rigidity index (mFRI) method to resolve this problem.
The proposed mFRI utilizes two or three correlation kernels parametrized at
different length scales to capture protein interactions at corresponding
scales. It is about 20% more accurate than the Gaussian network model (GNM) in
the B-factor prediction of a set of 364 proteins. Additionally, the present
method is able to delivery accurate predictions for multiscale macromolecules
that fail GNM. Finally, or a protein of $N$ residues, mFRI is of linear scaling
(O(N)) in computational complexity, in contrast to the order of O(N^3) for GNM.


  (Abridged) We present palaeoeconomy reconstructions for pre-modern
agriculture, with the Cucuteni-Trypillia Cultural unity (5,400-2,700 BC, modern
Ukraine, Moldova and Romania) as example. The starting point of our analysis is
the palaeodiet structure suggested by archaeological data, stable isotope
analyses of human remains, and palynology. We allow for the archeologically
attested contributions of domesticated and wild animal products to the diet,
develop plausible estimates of the yield of ancient cereal varieties cultivated
with ancient techniques, and quantify the yield dependence on the time after
initial planting and on rainfall (as a climate proxy). Our conclusions involve
analysis of the labour costs of the agricultural cycle of both an individual
and a farmer's family. Finally, we put our results into the context of the
exploitation territory and catchment analysis. The simplest economic complex
based on cereals, domestic and wild animal products, with fallow cropping,
appears to be capable of supporting an isolated, relatively small farming
community of 50-300 people (2-10 ha). Our results strongly suggest that dairy
products played a significant role in the dietary and labour balance. The
smaller settlements are typical of the earliest Trypillia A but remain
predominant at the later stages. A larger settlement of several hundred people
could function in isolation only with technological innovations, such as manure
fertiliser and ard tillage. Very large settlements of a few hundred hectares
could function only if supported by satellite farming villages and stable
exchange networks. We also discuss, quantify and assess some strategies to
mitigate the risks of arable agriculture associated with strong temporal
fluctuations in the cereal yield, such as manure fertilisation, increased
fraction of cereals in the diet combined with producing grain surplus for
emergency storage.


  Enterotoxigenic Escherichia coli (ETEC) are a major cause of diarrhea
worldwide, and infection of children in underdeveloped countries often leads to
high mortality rates. Isolated ETEC express a plethora of colonization factors
(fimbriae/pili), of which CFA/I and CFA/II that are assembled via the alternate
chaperone pathway (ACP), are amongst the most common. Fimbriae are filamentous
structures, whose shafts are primarily composed of helically arranged single
pilin-protein subunits, with a unique biomechanical capability allowing them to
unwind and rewind. A sustained ETEC infection, under adverse conditions of
dynamic shear forces, is primarily attributed to this biomechanical feature of
ETEC fimbriae. Recent understandings about the role of fimbriae as virulence
factors are pointing to an evolutionary adaptation of their structural and
biomechanical features. In this work, we investigated the biophysical
properties of CS2 fimbriae from the CFA/II group. Homology modelling its major
structural subunit CotA reveals structural clues and these are related to the
niche in which they are expressed. Using optical tweezers force spectroscopy we
found that CS2 fimbriae unwind at a constant force of 10 pN and have a corner
velocity of 1300 nm/s, i.e., the velocity at which the force required for
unwinding rises exponentially with increased speed. The biophysical properties
of CS2 fimbriae assessed in this work classify them into a low-force unwinding
group of fimbriae together with the CFA/I and CS20 fimbriae expressed by ETEC
strains. The three fimbriae are expressed by ETEC, colonize in similar gut
environments, and exhibit similar biophysical features, but differ in their
biogenesis. Our observation suggests that the environment has a strong impact
on the biophysical characteristics of fimbriae expressed by ETEC.


  Multiple-objective optimization is common in biological systems. In the
mammalian olfactory system, each sensory neuron stochastically expresses only
one out of up to thousands of olfactory receptor (OR) gene alleles; at organism
level the types of expressed ORs need to be maximized. Existing models focus
only on monoallele activation, and cannot explain recent observations in
mutants, especially the reduced global diversity of expressed ORs in G9a/GLP
knockouts. In this work we integrated existing information on OR expression,
and constructed a comprehensive model that has all its components based on
physical interactions. Analyzing the model reveals an evolutionarily optimized
three-layer regulation mechanism, which includes zonal segregation, epigenetic
barrier crossing coupled to a negative feedback loop that mechanistically
differs from previous theoretical proposals, and a previously unidentified
enhancer competition step. This model not only recapitulates monoallelic OR
expression, but also elucidates how the olfactory system maximizes and
maintains the diversity of OR expression, and has multiple predictions
validated by existing experimental results. Through making analogy to a
physical system with thermally activated barrier crossing and comparative
reverse engineering analyses, the study reveals that the olfactory receptor
selection system is optimally designed, and particularly underscores
cooperativity and synergy as a general design principle for multi-objective
optimization in biology.


  Recent experimental advances in biology allow researchers to obtain gene
expression profiles at single-cell resolution over hundreds, or even thousands
of cells at once. These single-cell measurements provide snapshots of the
states of the cells that make up a tissue, instead of the population-level
averages provided by conventional high-throughput experiments. This new data
therefore provides an exciting opportunity for computational modelling. In this
paper we introduce the idea of viewing single-cell gene expression profiles as
states of an asynchronous Boolean network, and frame model inference as the
problem of reconstructing a Boolean network from its state space. We then give
a scalable algorithm to solve this synthesis problem. We apply our technique to
both simulated and real data. We first apply our technique to data simulated
from a well established model of common myeloid progenitor differentiation. We
show that our technique is able to recover the original Boolean network rules.
We then apply our technique to a large dataset taken during embryonic
development containing thousands of cell measurements. Our technique
synthesises matching Boolean networks, and analysis of these models yields new
predictions about blood development which our experimental collaborators were
able to verify.


  Crystallography may be the gold standard of protein structure determination,
but obtaining the necessary high-quality crystals is also in some ways akin to
prospecting for the precious metal. The tools and models developed in soft
matter physics to understand colloidal assembly offer some insights into the
problem of crystallizing proteins. This topical review describes the various
analogies that have been made between proteins and colloids in that context. We
highlight the explanatory power of patchy particle models, but also the
challenges of providing guidance for crystallizing specific proteins. We
conclude with a presentation of possible future research directions. This
article is intended for soft matter scientists interested in protein
crystallization as a self-assembly problem, and as an introduction to the
pertinent physics literature for protein scientists more generally.


  Complexes of physically interacting proteins constitute fundamental
functional units responsible for driving biological processes within cells. A
faithful reconstruction of the entire set of complexes is therefore essential
to understand the functional organization of cells. In this review, we discuss
the key contributions of computational methods developed till date
(approximately between 2003 and 2015) for identifying complexes from the
network of interacting proteins (PPI network). We evaluate in depth the
performance of these methods on PPI datasets from yeast, and highlight
challenges faced by these methods, in particular detection of sparse and small
or sub- complexes and discerning of overlapping complexes. We describe methods
for integrating diverse information including expression profiles and 3D
structures of proteins with PPI networks to understand the dynamics of complex
formation, for instance, of time-based assembly of complex subunits and
formation of fuzzy complexes from intrinsically disordered proteins. Finally,
we discuss methods for identifying dysfunctional complexes in human diseases,
an application that is proving invaluable to understand disease mechanisms and
to discover novel therapeutic targets. We hope this review aptly commemorates a
decade of research on computational prediction of complexes and constitutes a
valuable reference for further advancements in this exciting area.


  Cells sense external concentrations and, via biochemical signaling, respond
by regulating the expression of target proteins. Both in signaling networks and
gene regulation there are two main mechanisms by which the concentration can be
encoded internally: amplitude modulation (AM), where the absolute concentration
of an internal signaling molecule encodes the stimulus, and frequency
modulation (FM), where the period between successive bursts represents the
stimulus. Although both mechanisms have been observed in biological systems,
the question of when it is beneficial for cells to use either AM or FM is
largely unanswered. Here, we first consider a simple model for a single
receptor (or ion channel), which can either signal continuously whenever a
ligand is bound, or produce a burst in signaling molecule upon receptor
binding. We find that bursty signaling is more accurate than continuous
signaling only for sufficiently fast dynamics. This suggests that modulation
based on bursts may be more common in signaling networks than in gene
regulation. We then extend our model to multiple receptors, where continuous
and bursty signaling are equivalent to AM and FM respectively, finding that AM
is always more accurate. This implies that the reason some cells use FM is
related to factors other than accuracy, such as the ability to coordinate
expression of multiple genes or to implement threshold crossing mechanisms.


  Peptide Nucleic Acids (PNA) are non-natural oligonucleotides mimics.
{\gamma}-PNA backbone are formed by standard nucleic acids nucleobases
connected by a neutral N-(2-aminoethyl)glycine backbone linked by a peptide
bond. In this study, we use Nuclear Magnetic Resonance (NMR) and Circular
Dichroism (CD) to explore the properties of the supramolecular duplexes formed
by these species. We show that standard Watson-Crick base pair as well as
non-standard ones are formed in solution. The duplexes thus formed present
marked melting transition temperatures substantially higher than their nucleic
acid homologs. Moreover, the presence of a chiral group on the
{\gamma}-peptidic backbone increases further more this transition temperature,
leading to very stable duplexes.


  Two endemic plant species in the Mediterranean island of Cyprus, Crocus
cyprius and Ophrys kotschyi, were used as a case study. We have coupled climate
change scenarios, and land use change models with species distribution models.
Future land use scenarios were modelled by initially calculating the rate of
current land use changes between two time snapshots (2000 and 2006) on the
island, and based on these transition probabilities markov-chain cellular
automata were used to generate future land use changes for 2050. Climate change
scenarios A1B, A2, B1 and B2A were derived from the IPCC reports. Species
climatic preferences were derived from their current distributions using
classification trees while habitats preferences were derived from the Red Data
Book of the Flora of Cyprus. A bioclimatic model for Crocus cyprius was built
using mean temperature of wettest quarter, max temperature of warmest month and
precipitation seasonality, while for Ophrys kotchyi the bioclimatic model was
built using precipitation of wettest month, mean temperature of warmest
quarter, isothermality, precipitation of coldest quarter, and annual
precipitation. Sequentially, simulation scenarios were performed regarding
future species distributions by accounting climate alone and both climate and
land use changes. The distribution of the two species resulting from the
bioclimatic models was then filtered by future land use changes, providing the
species projected potential distribution. The species projected potential
distribution varies depending on the type and scenario used, but many of both
species current sites/locations are projected to be outside their future
potential distribution. Our results demonstrate the importance of including
both land use and climatic changes in predictive species modeling.


  Many proteins carry out their biological functions by forming the
characteristic tertiary structures. Therefore, the search of the stable states
of proteins by molecular simulations is important to understand their functions
and stabilities. However, getting the stable state by conformational search is
difficult, because the energy landscape of the system is characterized by many
local minima separated by high energy barriers. In order to overcome this
difficulty, various sampling and optimization methods for conformations of
proteins have been proposed. In this study, we propose a new conformational
search method for proteins by using genetic crossover and Metropolis criterion.
We applied this method to an $\alpha$-helical protein. The conformations
obtained from the simulations are in good agreement with the experimental
results.


  Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely
used in studies of enzymatic catalysis. Until recently, it has been cost
prohibitive to determine the asymptotic limit of key energetic and structural
properties with respect to increasingly large QM regions. Leveraging recent
advances in electronic structure efficiency and accuracy, we investigate
catalytic properties in catechol O-methyltransferase, a representative example
of a methyltransferase critical to human health. Using QM regions ranging in
size from reactants-only (64 atoms) to nearly one-third of the entire protein
(940 atoms), we show that properties such as the activation energy approach
within chemical accuracy of the large-QM asymptotic limits rather slowly,
requiring approximately 500-600 atoms if the QM residues are chosen simply by
distance from the substrate. This slow approach to asymptotic limit is due to
charge transfer from protein residues to the reacting substrates. Our large
QM/MM calculations enable identification of charge separation for fragments in
the transition state as a key component of enzymatic methyl transfer rate
enhancement. We introduce charge shift analysis that reveals the minimum number
of protein residues (ca. 11-16 residues or 200-300 atoms for COMT) needed for
quantitative agreement with large-QM simulations. The identified residues are
not those that would be typically selected using criteria such as chemical
intuition or proximity. These results provide a recipe for a more careful
determination of QM region sizes in future QM/MM studies of enzymes.


  Symbiotic interactions between microbes and their multicellular hosts have
manifold impacts on molecular, cellular and organismal biology. To identify
candidate bacterial genes involved in maintaining endosymbiotic associations
with insect hosts, we analyzed genome-wide patterns of gene expression in the
alpha-proteobacteria Wolbachia pipientis across the life cycle of Drosophila
melanogaster using public data from the modENCODE project that was generated in
a Wolbachia-infected version of the ISO1 reference strain. We find that the
majority of Wolbachia genes are expressed at detectable levels in D.
melanogaster across the entire life cycle, but that only 7.8% of 1195 Wolbachia
genes exhibit robust stage- or sex-specific expression differences when studied
in the "holo-organism" context. Wolbachia genes that are differentially
expressed during development are typically up-regulated after D. melanogaster
embryogenesis, and include many bacterial membrane, secretion system and
ankyrin-repeat containing proteins. Sex-biased genes are often organised as
small operons of uncharacterised genes and are mainly up-regulated in adult
males D. melanogaster in an age-dependent manner suggesting a potential role in
cytoplasmic incompatibility. Our results indicate that large changes in
Wolbachia gene expression across the Drosophila life-cycle are relatively rare
when assayed across all host tissues, but that candidate genes to understand
host-microbe interaction in facultative endosymbionts can be successfully
identified using holo-organism expression profiling. Our work also shows that
mining public gene expression data in D. melanogaster provides a rich set of
resources to probe the functional basis of the Wolbachia-Drosophila symbiosis
and annotate the transcriptional outputs of the Wolbachia genome.


  The recent explosion of genomic data has underscored the need for
interpretable and comprehensive analyses that can capture complex phylogenetic
relationships within and across species. Recombination, reassortment and
horizontal gene transfer constitute examples of pervasive biological phenomena
that cannot be captured by tree-like representations. Starting from hundreds of
genomes, we are interested in the reconstruction of potential evolutionary
histories leading to the observed data. Ancestral recombination graphs
represent potential histories that explicitly accommodate recombination and
mutation events across orthologous genomes. However, they are computationally
costly to reconstruct, usually being infeasible for more than few tens of
genomes. Recently, Topological Data Analysis (TDA) methods have been proposed
as robust and scalable methods that can capture the genetic scale and frequency
of recombination. We build upon previous TDA developments for detecting and
quantifying recombination, and present a novel framework that can be applied to
hundreds of genomes and can be interpreted in terms of minimal histories of
mutation and recombination events, quantifying the scales and identifying the
genomic locations of recombinations. We implement this framework in a software
package, called TARGet, and apply it to several examples, including small
migration between different populations, human recombination, and horizontal
evolution in finches inhabiting the Gal\'apagos Islands.


  Inference with population genetic data usually treats the population pedigree
as a nuisance parameter, the unobserved product of a past history of random
mating. However, the history of genetic relationships in a given population is
a fixed, unobserved object, and so an alternative approach is to treat this
network of relationships as a complex object we wish to learn about, by
observing how genomes have been noisily passed down through it. This paper
explores this point of view, showing how to translate questions about
population genetic data into calculations with a Poisson process of mutations
on all ancestral genomes. This method is applied to give a robust
interpretation to the $f_4$ statistic used to identify admixture, and to design
a new statistic that measures covariances in mean times to most recent common
ancestor between two pairs of sequences. The method more generally interprets
population genetic statistics in terms of sums of specific functions over
ancestral genomes, thereby providing concrete, broadly interpretable
interpretations for these statistics. This provides a method for describing
demographic history without simplified demographic models. More generally, it
brings into focus the population pedigree, which is averaged over in
model-based demographic inference.


  We combined the genetic crossover, which is one of the operations of genetic
algorithm, and replica-exchange method in parallel molecular dynamics
simulations. The genetic crossover and replica-exchange method can search the
global conformational space by exchanging the corresponding parts between a
pair of conformations of a protein. In this study, we applied this method to an
$\alpha$-helical protein, Trp-cage mini protein, which has 20 amino-acid
residues. The conformations obtained from the simulations are in good agreement
with the experimental results.


  Divergence date estimates are central to understand evolutionary processes
and depend, in the case of molecular phylogenies, on tests of molecular clocks.
Here we propose two non-parametric tests of strict and relaxed molecular clocks
built upon a framework that uses the empirical cumulative distribution (ECD) of
branch lengths obtained from an ensemble of Bayesian trees and well known
non-parametric (one-sample and two-sample) Kolmogorov-Smirnov (KS)
goodness-of-fit test. In the strict clock case, the method consists in using
the one-sample Kolmogorov-Smirnov (KS) test to directly test if the phylogeny
is clock-like, in other words, if it follows a Poisson law. The ECD is computed
from the discretized branch lengths and the parameter $\lambda$ of the expected
Poisson distribution is calculated as the average branch length over the
ensemble of trees. To compensate for the auto-correlation in the ensemble of
trees and pseudo-replication we take advantage of thinning and effective sample
size, two features provided by Bayesian inference MCMC samplers. Finally, it is
observed that tree topologies with very long or very short branches lead to
Poisson mixtures and in this case we propose the use of the two-sample KS test
with samples from two continuous branch length distributions, one obtained from
an ensemble of clock-constrained trees and the other from an ensemble of
unconstrained trees. Moreover, in this second form the test can also be applied
to test for relaxed clock models. The use of a statistically equivalent
ensemble of phylogenies to obtain the branch lengths ECD, instead of one
consensus tree, yields considerable reduction of the effects of small sample
size and provides again of power.


  We propose and study {a generalized Beverton-Holt competition model} subject
to Allee effects to obtain insights on how the interplay of Allee effects and
contest competition affects the persistence and the extinction of two competing
species. By using {the theory of monotone dynamics} and the properties of
critical curves for non-invertible maps, our analysis shows that our model has
simple equilibrium dynamics for most conditions. The coexistence of two
competing species occurs only if the system has four interior equilibria.} We
provide an {approximation to} the basins of the boundary attractors (i.e., the
extinction of one or both species) where our results suggests that {contest
species are more prone to extinction than scramble ones are at low densities}.
In addition, {in comparison to the} dynamics of two species scramble
competition models subject to Allee effects, our study suggests that (i) Both
contest and scramble competition models can have only three boundary attractors
without the coexistence equilibria, or four attractors among which only one is
the {persistent attractor}, {whereas} scramble competition {models} may have
the extinction of both species as its only attractor under certain conditions,
i.e., \emph{the essential extinction} of two species due to strong Allee
effects; (ii) {Scramble competition models like Ricker type models can} have
much more complicated dynamical structure of interior attractors than contest
ones like Beverton-Holt type models have; and (iii) Scramble competition models
{like Ricker type competition models} may be more likely to promote the
coexistence of two species at low and high densities under certain conditions.


  Recent experiments on the bacterial flagellar motor have shown that the
structure of this nanomachine, which drives locomotion in a wide range of
bacterial species, is more dynamic than previously believed. Specifically, the
number of active torque-generating units (stators) was shown to vary across
applied loads. This finding invalidates the experimental evidence reporting
that limiting (zero-torque) speed is independent of the number of active
stators. Here, we propose that, contrary to previous assumptions, the maximum
speed of the motor is not universal, but rather increases as additional
torque-generators are recruited. This result arises from our assumption that
stators disengage from the motor for a significant portion of their
mechanochemical cycles at low loads. We show that this assumption is consistent
with current experimental evidence and consolidate our predictions with
arguments that a processive motor must have a high duty ratio at high loads.


  Evolution in finite populations is often modelled using the classical Moran
process. Over the last ten years this methodology has been extended to
structured populations using evolutionary graph theory. An important question
in any such population, is whether a rare mutant has a higher or lower chance
of fixating (the fixation probability) than the Moran probability, i.e. that
from the original Moran model, which represents an unstructured population. As
evolutionary graph theory has developed, different ways of considering the
interactions between individuals through a graph and an associated matrix of
weights have been considered, as have a number of important dynamics. In this
paper we revisit the original paper on evolutionary graph theory in light of
these extensions to consider these developments in an integrated way. In
particular we find general criteria for when an evolutionary graph with general
weights satisfies the Moran probability for the set of six common evolutionary
dynamics.


  Although the role of lateral gene transfer is well recognized in the
evolution of bacteria, it is generally assumed that it has had less influence
among eukaryotes. To explore this hypothesis we compare the dynamics of genome
evolution in two groups of organisms: Cyanobacteria and Fungi. Ancestral
genomes are inferred in both clades using two types of methods. First, Count, a
gene tree unaware method that models gene duplications, gains and losses to
explain the observed numbers of genes present in a genome. Second, ALE, a more
recent gene tree-aware method that reconciles gene trees with a species tree
using a model of gene duplication, loss, and transfer. We compare their merits
and their ability to quantify the role of transfers, and assess the impact of
taxonomic sampling on their inferences. We present what we believe is
compelling evidence that gene transfer plays a significant role in the
evolution of Fungi.


  The Luria-Delbr{\"u}ck experiment is a cornerstone of evolutionary theory,
demonstrating the randomness of mutations before selection. The distribution of
the number of mutants in this experiment has been the subject of intense
investigation during the last 70 years. Despite this considerable effort, most
of the results have been obtained under the assumption of constant growth rate,
which is far from the experimental condition. We derive here the properties of
this distribution for arbitrary growth function, for both the deterministic and
stochastic growth of the mutants. The derivation we propose uses the number of
wild type bacteria as the independent variable instead of time. The derivation
is surprisingly simple and versatile, allowing many generalizations to be taken
easily into account.


  In this work, we report the heat rectifying capability of a-helices. Using
molecular dynamics simulations we show an increased thermal diffusivity in the
C-Terminal to N-Terminal direction of propagation. The origin of this effect
seems to be a function of the particular orientation of the hydrogen bonds
stabilizing these a-helices. Our results may be relevant for the design of
thermal rectification devices for materials science and lend support to the
role of normal length hydrogen bonds in the asymmetrical energy flow in
proteins.


  The Set Covering Machine (SCM) is a greedy learning algorithm that produces
sparse classifiers. We extend the SCM for datasets that contain a huge number
of features. The whole genetic material of living organisms is an example of
such a case, where the number of feature exceeds 10^7. Three human pathogens
were used to evaluate the performance of the SCM at predicting antimicrobial
resistance. Our results show that the SCM compares favorably in terms of
sparsity and accuracy against L1 and L2 regularized Support Vector Machines and
CART decision trees. Moreover, the SCM was the only algorithm that could
consider the full feature space. For all other algorithms, the latter had to be
filtered as a preprocessing step.


  Detecting the interactions of genetic compounds like genes, SNPs, proteins,
metabolites, etc. can potentially unravel the mechanisms behind complex traits
and common genetic disorders. Several methods have been taken into
consideration for the analysis of different types of genetic data, regression
being one of the most widely adopted. Without any doubt, a common data type is
represented by gene expression profiles, from which gene regulatory networks
have been inferred with different approaches. In this work we review nine
penalised regression methods applied to microarray data to infer the topology
of the network of interactions. We evaluate each method with respect to the
complexity of biological data. We analyse the limitations of each of them in
order to suggest a number of precautions that should be considered to make
their predictions more significant and reliable.


  Comment on "Fluctuations in the DNA double helix: A critical review" by
Frank-Kamenetskii and Prakash


  Biological systems are driven by intricate interactions among the complex
array of molecules that comprise the cell. Many methods have been developed to
reconstruct network models of those interactions. These methods often draw on
large numbers of samples with measured gene expression profiles to infer
connections between genes (or gene products). The result is an aggregate
network model representing a single estimate for the likelihood of each
interaction, or "edge," in the network. While informative, aggregate models
fail to capture the heterogeneity that is represented in any population. Here
we propose a method to reverse engineer sample-specific networks from aggregate
network models. We demonstrate the accuracy and applicability of our approach
in several data sets, including simulated data, microarray expression data from
synchronized yeast cells, and RNA-seq data collected from human lymphoblastoid
cell lines. We show that these sample-specific networks can be used to study
changes in network topology across time and to characterize shifts in gene
regulation that may not be apparent in expression data. We believe the ability
to generate sample-specific networks will greatly facilitate the application of
network methods to the increasingly large, complex, and heterogeneous
multi-omic data sets that are currently being generated, and ultimately support
the emerging field of precision network medicine.


  A major challenge in next-generation genome sequencing (NGS) is to assemble
massive overlapping short reads that are randomly sampled from DNA fragments.
To complete assembling, one needs to finish a fundamental task in many leading
assembly algorithms: counting the number of occurrences of k-mers (length-k
substrings in sequences). The counting results are critical for many components
in assembly (e.g. variants detection and read error correction). For large
genomes, the k-mer counting task can easily consume a huge amount of memory,
making it impossible for large-scale parallel assembly on commodity servers.
  In this paper, we develop MSPKmerCounter, a disk-based approach, to
efficiently perform k-mer counting for large genomes using a small amount of
memory. Our approach is based on a novel technique called Minimum Substring
Partitioning (MSP). MSP breaks short reads into multiple disjoint partitions
such that each partition can be loaded into memory and processed individually.
By leveraging the overlaps among the k-mers derived from the same short read,
MSP can achieve astonishing compression ratio so that the I/O cost can be
significantly reduced. For the task of k-mer counting, MSPKmerCounter offers a
very fast and memory-efficient solution. Experiment results on large real-life
short reads data sets demonstrate that MSPKmerCounter can achieve better
overall performance than state-of-the-art k-mer counting approaches.
  MSPKmerCounter is available at http://www.cs.ucsb.edu/~yangli/MSPKmerCounter


  In recent years experiments have demonstrated that living cells can measure
low chemical concentrations with high precision, and much progress has been
made in understanding what sets the fundamental limit to the precision of
chemical sensing. Chemical concentration measurements start with the binding of
ligand molecules to receptor proteins, which is an inherently noisy process,
especially at low concentrations. The signaling networks that transmit the
information on the ligand concentration from the receptors into the cell have
to filter this noise extrinsic to the cell as much as possible. These networks,
however, are also stochastic in nature, which means that they will also add
noise to the transmitted signal. In this review, we will first discuss how the
diffusive transport and binding of ligand to the receptor sets the receptor
correlation time, and then how downstream signaling pathways integrate the
noise in the receptor state; we will discuss how the number of receptors, the
receptor correlation time, and the effective integration time together set a
fundamental limit on the precision of sensing. We then discuss how cells can
remove the receptor noise while simultaneously suppressing the intrinsic noise
in the signaling network. We describe why this mechanism of time integration
requires three classes of resources---receptors and their integration time,
readout molecules, energy---and how each resource class sets a fundamental
sensing limit. We also briefly discuss the scheme of maximum-likelihood
estimation, the role of receptor cooperativity, and how cellular copy protocols
differ from canonical copy protocols typically considered in the computational
literature, explaining why cellular sensing systems can never reach the
Landauer limit on the optimal trade-off between accuracy and energetic cost.


  Mating preferences of many biological species are not constant but
season-dependent. Within the framework of evolutionary game theory this can be
modeled with two finite opposite-sex populations playing against each other
following the rules that are periodically changing. By combining Floquet theory
and the concept of quasi-stationary distributions, we reveal existence of
metastable time-periodic states in the evolution of finite game-driven
populations. The evolutionary Floquet states correspond to time-periodic
probability flows in the strategy space which cannot be resolved within the
mean-field framework. The lifetime of metastable Floquet states increases with
the size $N$ of populations so that they become attractors in the limit $N
\rightarrow \infty$.


  Maximum parsimony is one of the most frequently-discussed tree reconstruction
methods in phylogenetic estimation. However, in recent years it has become more
and more apparent that phylogenetic trees are often not sufficient to describe
evolution accurately. For instance, processes like hybridization or lateral
gene transfer that are commonplace in many groups of organisms and result in
mosaic patterns of relationships cannot be represented by a single phylogenetic
tree. This is why phylogenetic networks, which can display such events, are
becoming of more and more interest in phylogenetic research. It is therefore
necessary to extend concepts like maximum parsimony from phylogenetic trees to
networks. Several suggestions for possible extensions can be found in recent
literature, for instance the softwired and the hardwired parsimony concepts. In
this paper, we analyze the so-called big parsimony problem under these two
concepts, i.e. we investigate maximum parsimonious networks and analyze their
properties. In particular, we show that finding a softwired maximum parsimony
network is possible in polynomial time. We also show that the set of maximum
parsimony networks for the hardwired definition always contains at least one
phylogenetic tree. Lastly, we investigate some parallels of parsimony to
different likelihood concepts on phylogenetic networks.


  Metagenomics characterizes the taxonomic diversity of microbial communities
by sequencing DNA directly from an environmental sample. One of the main
challenges in metagenomics data analysis is the binning step, where each
sequenced read is assigned to a taxonomic clade. Due to the large volume of
metagenomics datasets, binning methods need fast and accurate algorithms that
can operate with reasonable computing requirements. While standard
alignment-based methods provide state-of-the-art performance, compositional
approaches that assign a taxonomic class to a DNA read based on the k-mers it
contains have the potential to provide faster solutions. In this work, we
investigate the potential of modern, large-scale machine learning
implementations for taxonomic affectation of next-generation sequencing reads
based on their k-mers profile. We show that machine learning-based
compositional approaches benefit from increasing the number of fragments
sampled from reference genome to tune their parameters, up to a coverage of
about 10, and from increasing the k-mer size to about 12. Tuning these models
involves training a machine learning model on about 10 8 samples in 10 7
dimensions, which is out of reach of standard soft-wares but can be done
efficiently with modern implementations for large-scale machine learning. The
resulting models are competitive in terms of accuracy with well-established
alignment tools for problems involving a small to moderate number of candidate
species, and for reasonable amounts of sequencing errors. We show, however,
that compositional approaches are still limited in their ability to deal with
problems involving a greater number of species, and more sensitive to
sequencing errors. We finally confirm that compositional approach achieve
faster prediction times, with a gain of 3 to 15 times with respect to the
BWA-MEM short read mapper, depending on the number of candidate species and the
level of sequencing noise.


  Neutral models for the dynamics of a system of competing species are used,
nowadays, to describe a wide variety of empirical communities. These models are
used in many situations, ranging from population genetics and ecological
biodiversity to macroevolution and cancer tumors. One of the main issues
discussed within this framework is the relationships between the abundance of a
species and its age. Here we provide a comprehensive analysis of the
age-abundance relationships for fixed-size and growing communities. Explicit
formulas for the average and the most likely age of a species with abundance
$n$ are given, together with the full probability distribution function. We
further discuss the universality of these results and their applicability to
the tropical forest community.


  The 2014 Ebola outbreak in west Africa raised many questions about the
control of infectious disease in an increasingly connected global society.
Limited availability of contact information made contact tracing difficult or
impractical in combating the outbreak. We consider the development of
multi-scale public health strategies and simulate policies for community-level
response aimed at early screening of communities rather than individuals, as
well as travel restrictions to prevent community cross-contamination. Our
analysis shows the policies to be effective even at a relatively low level of
compliance. In our simulations, 40% of individuals conforming to these policies
is enough to stop the outbreak. Simulations with a 50% compliance rate are
consistent with the case counts in Liberia during the period of rapid decline
after mid September, 2014. We also find the travel restriction to be effective
at reducing the risks associated with compliance substantially below the 40%
level, shortening the outbreak and enabling efforts to be focused on affected
areas. Our results suggest that the multi-scale approach can be used to further
evolve public health strategy for defeating emerging epidemics.


  Chromosomal rearrangements, which shuffle DNA throughout the genome, are an
important source of divergence across taxa. Using a paired-end read approach
with Illumina sequence data for archaic humans, I identify changes in genome
structure that occurred recently in human evolution. Hundreds of rearrangements
indicate genomic trafficking between the sex chromosomes and autosomes, raising
the possibility of sex-specific changes. Additionally, genes adjacent to genome
structure changes in Neanderthals are associated with testis-specific
expression, consistent with evolutionary theory that new genes commonly form
with expression in the testes. I identify one case of new-gene creation through
transposition from the Y chromosome to chromosome 10 that combines the 5' end
of the testis-specific gene Fank1 with previously untranscribed sequence. This
new transcript experienced copy number expansion in archaic genomes, indicating
rapid genomic change. Among rearrangements identified in Neanderthals, 13% are
transposition of selfish genetic elements, while 32% appear to be ectopic
exchange between repeats. In Denisovan, the pattern is similar but numbers are
significantly higher with 18% of rearrangements reflecting transposition and
40% ectopic exchange between distantly related repeats. There is an excess of
divergent rearrangements relative to polymorphism in Denisovan, which might
result from non-uniform rates of mutation, possibly reflecting a burst of TE
activity in the lineage that led to Denisovan. Finally, loci containing genome
structure changes show diminished rates of introgression from Neanderthals into
modern humans, consistent with the hypothesis that rearrangements serve as
barriers to gene flow during hybridization. Together, these results suggest
that this previously unidentified source of genomic variation has important
biological consequences in human evolution.


  The cellular basis for the adaptive immune response during antigen
recognition relies on a specialized protein interface known as the
immunological synapse (IS). Understanding the biophysical basis for protein
patterning by deciphering the quantitative rules for their formation and motion
is an important aspect of characterizing immune cell recognition and thence the
rules for immune system activation. We propose a minimal mathematical model for
the physical basis of membrane protein patterning in the IS, which encompass
membrane mechanics, protein binding kinetics and motion, and fluid flow in the
synaptic cleft. Our theory leads to simple predictions for the spatial and
temporal scales of protein cluster formation, growth and arrest as a function
of membrane stiffness, rigidity and kinetics of the adhesive proteins, and the
fluid in the synaptic cleft. Numerical simulations complement these scaling
laws by quantifying the nucleation, growth and stabilization of proteins
domains on the size of the cell. Direct comparison with experiment shows that
passive elastohydrodynamics and kinetics of protein binding in the synaptic
cleft can describe the short-time formation and organization of protein
clusters, without evoking any active processes in the cytoskeleton. Despite the
apparent complexity of the process, our analysis highlights the role of just
two dimensionless parameters that characterize the spatial and temporal
evolution of the protein pattern: a ratio of membrane elasticity to protein
elasticity, and the ratio of a hydrodynamic time scale for fluid flow relative
to the protein binding rate, and we present a simple phase diagram that
encompasses the variety of patterns that can arise.


  SweepFinder is a popular program that implements a powerful likelihood-based
method for detecting recent positive selection, or selective sweeps. Here, we
present SweepFinder2, an extension of SweepFinder with increased sensitivity
and robustness to the confounding effects of mutation rate variation and
background selection, as well as increased flexibility that enables the user to
examine genomic regions in greater detail and to specify a fixed distance
between test sites. Moreover, SweepFinder2 enables the use of invariant sites
for sweep detection, increasing both its power and precision relative to
SweepFinder.


  Gene regulatory networks (GRNs) are increasingly used for explaining
biological processes with complex transcriptional regulation. A GRN links the
expression levels of a set of genes via regulatory controls that gene products
exert on one another. Boolean networks are a common modeling choice since they
balance between detail and ease of analysis. However, even for Boolean networks
the problem of fitting a given network model to an expression dataset is
NP-Complete. Previous methods have addressed this issue heuristically or by
focusing on acyclic networks and specific classes of regulation functions. In
this paper we introduce a novel algorithm for this problem that makes use of
sampling in order to handle large datasets. Our algorithm can handle time
series data for any network type and steady state data for acyclic networks.
Using in-silico time series data we demonstrate good performance on large
datasets with a significant level of noise.


  Using molecular dynamics simulations, the structure of model mini-protein was
thoroughly characterized in the imidazolium-based amino acid ionic liquids and
their aqueous solutions. We report that the mini-protein is more stable when
AAIL is added as a cosolvent. Complete substitution of water by organic cations
and anions further results in hindered conformational flexibility of the
mini-protein. This observation suggests that AAILs are able to defend proteins
from thermally induced denaturation. We show by means of radial distributions
that the mini-protein is efficiently solvated by both solvents due to agood
mutual miscibility. However, amino acid based anions prevail in the first
coordination sphere of the mini-protein.


  Living systems need to be highly responsive, and also to keep fluctuations
low. These goals are incompatible in equilibrium systems due to the Fluctuation
Dissipation Theorem (FDT). Here, we show that biological sensory systems,
driven far from equilibrium by free energy consumption, can reduce their
intrinsic fluctuations while maintaining high responsiveness. By developing a
continuum theory of the E. coli chemotaxis pathway, we demonstrate that
adaptation can be understood as a non-equilibrium phase transition controlled
by free energy dissipation, and it is characterized by a breaking of the FDT.
We show that the maximum response at short time is enhanced by free energy
dissipation. At the same time, the low frequency fluctuations and the
adaptation error decrease with the free energy dissipation algebraically and
exponentially, respectively.


  Understanding cell-fate decisions during tumorigenesis and metastasis is a
major challenge in modern cancer biology. One canonical cell-fate decision that
cancer cells undergo is Epithelial-to-Mesenchymal Transition (EMT) and its
reverse Mesenchymal-to-Epithelial Transition (MET). While transitioning between
these two phenotypes - epithelial and mesenchymal - cells can also attain a
hybrid epithelial/mesenchymal (i.e. partial or intermediate EMT) phenotype.
Cells in this phenotype have mixed epithelial (e.g. adhesion) and mesenchymal
(e.g. migration) properties, thereby allowing them to move collectively as
clusters of Circulating Tumor Cells (CTCs). If these clusters enter the
circulation, they can be more apoptosis-resistant and more capable of
initiating metastatic lesions than cancer cells moving individually with wholly
mesenchymal phenotypes, having undergo a complete EMT. Here, we review the
operating principles of the core regulatory network for EMT/MET that acts as a
three-way switch giving rise to three distinct phenotypes - epithelial,
mesenchymal and hybrid epithelial/mesenchymal. We further characterize this
hybrid E/M phenotype in terms of its capabilities in terms of collective cell
migration, tumor-initiation, cell-cell communication, and drug resistance. We
elucidate how the highly interconnected coupling between these modules
coordinates cell-fate decisions among a population of cancer cells in the
dynamic tumor, hence facilitating tumor-stoma interactions, formation of CTC
clusters, and consequently cancer metastasis. Finally, we discuss the multiple
advantages that the hybrid epithelial/mesenchymal phenotype have as compared to
a complete EMT phenotype and argue that these collectively migrating cells are
the primary 'bad actors' of metastasis.


  Self-assembly in the laboratory can now yield `information-rich'
nanostructures in which each component is of a distinct type and has a defined
spatial position. Ensuring the thermodynamic stability of such structures
requires inter-component interaction energies to increase logarithmically with
structure size, in order to counter the entropy gained upon mixing component
types in solution. However, self-assembly in the presence of strong
interactions results in general in kinetic trapping, so suggesting a limit to
the size of an (equilibrium) structure that can be self-assembled from
distinguishable components. Here we study numerically a two-dimensional
hierarchical assembly scheme already considered in experiment. We show that
this scheme is immune to the kinetic traps associated with strong `native'
interactions (interactions designed to stabilize the intended structure), and
so, in principle, offers a way to make large information-rich structures. In
this scheme the size of an assembled structure scales exponentially with the
stage of assembly, and assembly can continue as long as random motion is able
to bring structures into contact. The resulting superstructure could provide a
template for building in the third dimension. The chief drawback of this scheme
is that it is particularly susceptible to kinetic traps that result from
`non-native' interactions (interactions not required to stabilize the intended
structure); the scale on which such a scheme can be realized therefore depends
upon how effectively this latter kind of interaction can be suppressed.


  The ability to monitor nutrient and other environmental conditions with high
sensitivity is crucial for cell growth and survival. Sensory adaptation allows
a cell to recover its sensitivity after a transient response to a shift in the
strength of extracellular stimulus. The working principles of adaptation have
been established previously based on rate equations which do not consider
fluctuations in a thermal environment. Recently, G. Lan et al. (Nature Phys.,
8:422-8, 2012) performed a detailed analysis of a stochastic model for the E.
coli sensory network. They showed that accurate adaptation is possible only
when the system operates in a nonequilibrium steady-state (NESS). They further
proposed an energy-speed-accuracy (ESA) trade-off relation. We present here
analytic results on the NESS of the model through a mapping to a
one-dimensional birth-death process. An exact expression for the entropy
production rate is also derived. Based on these results, we are able to discuss
the ESA relation in a more general setting. Our study suggests that the
adaptation error can be reduced exponentially as the methylation range
increases. Finally, we show that a nonequilibrium phase transition exists in
the infinite methylation range limit, despite the fact that the model contains
only two discrete variables.


  The mechanisms of cold- and pressure-denaturation of proteins are matter of
debate and are commonly understood as due to water-mediated interactions. Here
we study several cases of proteins, with or without a unique native state, with
or without hydrophilic residues, by means of a coarse-grain protein model in
explicit solvent. We show, using Monte Carlo simulations, that taking into
account how water at the protein interface changes its hydrogen bond properties
and its density fluctuations is enough to predict protein stability regions
with elliptic shapes in the temperature-pressure plane, consistent with
previous theories. Our results clearly identify the different mechanisms with
which water participates to denaturation and open the perspective to develop
advanced computational design tools for protein engineering.


  Mutational neighbourhoods in genotype-phenotype (GP) maps are widely believed
to be more likely to share characteristics than expected from random chance.
Such genetic correlations should, as John Maynard Smith famously pointed out,
strongly influence evolutionary dynamics. We explore and quantify these
intuitions by comparing three GP maps - RNA SS, HP for tertiary, Polyominoes
for protein quaternary structure - to a simple random null model that maintains
the number of genotypes mapping to each phenotype, but assigns genotypes
randomly. The mutational neighbourhood of a genotype in these GP maps is much
more likely to contain (mutationally neutral) genotypes mapping to the same
phenotype than in the random null model. These neutral correlations can
increase the robustness to mutations by orders of magnitude over that of the
null model, raising robustness above the critical threshold for the formation
of large neutral networks that enhance the capacity for neutral exploration. We
also study {\em non-neutral correlations}: Compared to the null model, i) If a
particular (non-neutral) phenotype is found once in the 1-mutation
neighbourhood of a genotype, then the chance of finding that phenotype multiple
times in this neighbourhood is larger than expected; ii) If two genotypes are
connected by a single neutral mutation, then their respective non-neutral
1-mutation neighbourhoods are more likely to be similar; iii) If a genotype
maps to a folding or self-assembling phenotype, then its non-neutral neighbours
are less likely to be a potentially deleterious non-folding or non-assembling
phenotype. Non-neutral correlations of type i) and ii) reduce the rate at which
new phenotypes can be found by neutral exploration, and so may diminish
evolvability, while non-neutral correlations of type iii) may instead
facilitate evolutionary exploration and so increase evolvability.


  BarraCUDA is a C program which uses the BWA algorithm in parallel with nVidia
CUDA to align short next generation DNA sequences against a reference genome.
The genetically improved (GI) code is up to three times faster on short paired
end reads from The 1000 Genomes Project and 60percent more accurate on a short
BioPlanet.com GCAT alignment benchmark. GPGPU Barracuda running on a single K80
Tesla GPU can align short paired end nextgen sequences up to ten times faster
than bwa on a 12 core CPU.


  Most phylogenetic models assume that the evolutionary process is stationary
and reversible. As a result, the root of the tree cannot be inferred as part of
the analysis because the likelihood of the data does not depend on the position
of the root. Yet defining the root of a phylogenetic tree is a key component of
phylogenetic inference because it provides a point of reference for polarising
ancestor/descendant relationships and therefore interpreting the tree. In this
paper we investigate the effect of relaxing the reversibility assumption and
allowing the position of the root to be another unknown in the model. We
propose two hierarchical models that are centred on a reversible model but
perturbed to allow non-reversibility. The models differ in the degree of
structure imposed on the perturbations. The analysis is performed in the
Bayesian framework using Markov chain Monte Carlo methods. We illustrate the
performance of the two non-reversible models in analyses of simulated data sets
using two types of topological priors. We then apply the models to a real
biological data set, the radiation of polyploid yeasts, for which there is a
robust biological opinion about the root position. Finally we apply the models
to a second biological data set for which the rooted tree is controversial: the
ribosomal tree of life. We compare the two non-reversible models and conclude
that both are useful in inferring the position of the root from real biological
data sets.


  It is shown how to take advantage of the frequent occurrence of double-strand
breaks in the genome of prokaryotic cells, in order to reduce their high
efficient reproductive capability. The analysis examines the physical status of
the free ends of each break and considers how this status can interfere with an
external physical apparatus, with the aim of undermining the repair processes.
We indicate the biological consequences of this interaction and we give an
approximate evaluation of the topological and dynamical effects that arise on
the genomic material involved. The overall result suggests a significant
reduction of the dynamics of the repair.


  Major discrepancies have been noted for some time between fossil ages and
molecular divergence dates for a variety of taxa. Recently, systematic trends
within avian clades have been uncovered. The trends show that the disparity is
much larger for mitochondrial DNA than for nuclear DNA, also that it is larger
for crown fossil dates than stem fossil dates. It has been argued that this
pattern is largely inconsistent with incompleteness of the fossil record as the
principal driver of the disparity. A case is presented that given the expected
mutations from a fluctuating background of astrophysical radiation from such
sources as supernovae, the rate of molecular clocks is variable and should
increase back to a few Ma, before returning to the long-term average rate. This
is a possible explanation for the disparity. One test of this hypothesis is to
look for an acceleration of molecular clocks 2 to 2.5 Ma due to one or more
moderately nearby supernovae known to have happened at that time. Another is to
look for reduced disparity in benthic organisms of the deep ocean. In addition,
due to the importance of highly penetrating muon irradiation, the disparity
should be magnified for megafauna.


  Since the arrival of genetic typing methods in the late 1960's, researchers
have puzzled at the clinical consequence of observed strain mixtures within
clinical isolates of Plasmodium falciparum. We present a new statistical model
that infers the number of strains present and the amount of admixture with the
local population (panmixia) using whole-genome sequence data. The model
provides a rigorous statistical approach to inferring these quantities as well
as the proportions of the strains within each sample. Applied to 168 samples of
whole-genome sequence data from northern Ghana, the model provides
significantly improvement fit over models implementing simpler approaches to
mixture for a large majority (129/168) of samples. We discuss the possible uses
of this model as a window into within-host selection for clinical and
epidemiological studies and outline possible means for experimental validation.


  Motivation: The consistent amount of different types of omics data requires
novel methods of analysis and data integration. In this work we describe
Regression2Net, a computational approach to analyse gene expression and
methylation profiles via regression analysis and network-based techniques.
  Results: We identified 284 and 447 unique candidate genes potentially
associated to the Glioblastoma pathology from two networks inferred from mixed
genetic datasets. In-depth biological analysis of these networks reveals genes
that are related to energy metabolism, cell cycle control (AATF), immune system
response and several types of cancer. Importantly, we observed significant
over- representation of cancer related pathways including glioma especially in
the methylation network. This confirms the strong link between methylation and
glioblastomas. Potential glioma suppressor genes ACCN3 and ACCN4 linked to
NBPF1 neuroblastoma breakpoint family have been identified in our expression
network. Numerous ABC transporter genes (ABCA1, ABCB1) present in the
expression network suggest drug resistance of glioblastoma tumors.


  DNA nanostructures are made using synthetic DNA strands, the sequences of
which are designed such that they will self-assemble into the desired form by
hybridization of complementary domains. Various structures and devices have
been presented, including DNA tweezers, nanorobots and a range of linear motors
such as bipedal walkers. Inspiration for the latter is drawn from naturally
occurring molecular motors like kinesin.
  This paper describes a concept for an autonomous rotary nanomotor made from
DNA, which utilizes the well-known and widely-studied phenomenon of
toehold-mediated DNA strand displacement. The motor is to be driven by a series
of strand displacement reactions, the order of which is controlled by steric
constraints arising from the secondary structure of the DNA strands comprising
the motor mechanism.
  The capabilities of DNA motors would be extended significantly if autonomous
rotary motion could be achieved. The device has a range of potential
applications, including molecular computation and single-molecule manipulation.


  Recent development of high-resolution mass spectrometry (MS) instruments
enables chemical cross-linking (XL) to become a high-throughput method for
obtaining structural information about proteins. Restraints derived from XL-MS
experiments have been used successfully for structure refinement and
protein-protein docking. However, one formidable question is under which
circumstances XL-MS data might be sufficient to determine a protein's tertiary
structure de novo? Answering this question will not only include understanding
the impact of XL-MS data on sampling and scoring within a de novo protein
structure prediction algorithm, it must also determine an optimal cross-linker
type and length for protein structure determination. While a longer
cross-linker will yield more restraints, the value of each restraint for
protein structure prediction decreases as the restraint is consistent with a
larger conformational space.
  In this study, the number of cross-links and their discriminative power was
systematically analyzed in silico on a set of 2,055 non-redundant protein folds
considering Lys-Lys, Lys-Asp, Lys-Glu, Cys-Cys, and Arg-Arg reactive
cross-linkers between 1 {\AA} and 60 {\AA}. Depending on the protein size a
heuristic was developed that determines the optimal cross-linker length. Next,
simulated restraints of variable length were used to de novo predict the
tertiary structure of fifteen proteins using the BCL::Fold algorithm. The
results demonstrate that a distinct cross-linker length exists for which
information content for de novo protein structure prediction is maximized. The
sampling accuracy improves on average by 1.0 {\AA} and up to 2.2 {\AA} in the
most prominent example. XL-MS restraints enable consistently an improved
selection of native-like models with an average enrichment of 2.1.


  As costs of genome sequencing have dropped precipitously, development of
efficient bioinformatic methods to analyze genome structure and evolution have
become ever more urgent. For example, most published phylogenomic studies
involve either massive concatenation of sequences, or informal comparisons of
phylogenies inferred on a small subset of orthologous genes, neither of which
provides a comprehensive overview of evolution or systematic identification of
genes with unusual and interesting evolution (e.g. horizontal gene transfers,
gene duplication and subsequent neofunctionalization). We are interested in
identifying such "outlying" gene trees from the set of gene trees and
estimating the distribution of the tree over the "tree space".
  This paper describes an improvement to the kdetrees algorithm, an adaptation
of classical kernel density estimation to the metric space of phylogenetic
trees (Billera-Holmes-Vogtman treespace), whereby the kernel normalizing
constants, are estimated through the use of the novel holonomic gradient
methods. As the original kdetrees paper, we have applied kdetrees to a set of
Apicomplexa genes and it identified several unreliable sequence alignments
which had escaped previous detection, as well as a gene independently reported
as a possible case of horizontal gene transfer.


  Transcriptional and post-transcriptional regulation shape
tissue-type-specific proteomes, but their relative contributions remain
contested. Estimates of the factors determining protein levels in human tissues
do not distinguish between (i) the factors determining the variability between
the abundances of different proteins, i.e., mean-level-variability and, (ii)
the factors determining the physiological variability of the same protein
across different tissue types, i.e., across-tissues variability. We sought to
estimate the contribution of transcript levels to these two orthogonal sources
of variability, and found that scaled mRNA levels can account for most of the
mean-level-variability but not necessarily for across-tissues variability. The
reliable quantification of the latter estimate is limited by substantial
measurement noise. However, protein-to-mRNA ratios exhibit substantial
across-tissues variability that is functionally concerted and reproducible
across different datasets, suggesting extensive post-transcriptional
regulation. These results caution against estimating protein fold-changes from
mRNA fold-changes between different cell-types, and highlight the contribution
of post-transcriptional regulation to shaping tissue-type-specific proteomes.


  Cells use surface receptors to estimate the concentration of external
ligands. Limits on the accuracy of such estimations have been well studied for
pairs of ligand and receptor species. However, the environment typically
contains many ligands, which can bind to the same receptors with different
affinities, resulting in cross-talk. In traditional rate models, such
cross-talk prevents accurate inference of individual ligand concentrations. In
contrast, here we show that knowing the precise timing sequence of stochastic
binding and unbinding events allows one receptor to provide information about
multiple ligands simultaneously and with a high accuracy. We argue that such
high-accuracy estimation of multiple concentrations can be realized by the
familiar kinetic proofreading mechanism.


  Species evolution is essentially a random process of interaction between
biological populations and their environments. As a result, some physical
parameters in evolution models are subject to statistical fluctuations. In this
paper, two important parameters in the Eigen model, the fitness and mutation
rate, are treated as Gaussian distributed random variables simultaneously to
examine the property of the error threshold. Numerical simulation results show
that the error threshold in the fully random model appears as a crossover
region instead of a phase transition point, and as the fluctuation strength
increases the crossover region becomes smoother and smoother. Furthermore, it
is shown that the randomization of the mutation rate plays a dominant role in
changing the error threshold in the fully random model, which is consistent
with the existing experimental data. The implication of the threshold change
due to the randomization for antiviral strategies is discussed.


  The stochastic Eigen model proposed by Feng et al. (Journal of theoretical
biology, 246 (2007) 28) showed that error threshold is no longer a phase
transition point but a crossover region whose width depends on the strength of
the random fluctuation in an environment. The underlying cause of this
phenomenon has not yet been well examined. In this article, we adopt a single
peak Gaussian distributed fitness landscape instead of a constant one to
investigate and analyze the change of the error threshold and the statistical
property of the quasi-species population. We find a roughly linear relation
between the width of the error threshold and the fitness fluctuation strength.
For a given quasi-species, the fluctuation of the relative concentration has a
minimum with a normal distribution of the relative concentration at the maximum
of the averaged relative concentration, it has however a largest value with a
bimodal distribution of the relative concentration near the error threshold.
The above results deepen our understanding of the quasispecies and error
threshold and are heuristic for exploring practicable antiviral strategies.


  Due to rapid advancement in high-throughput techniques, such as microarrays
and next generation sequencing technologies, biological data are increasing
exponentially. The current challenge in computational biology and
bioinformatics research is how to analyze these huge raw biological data to
extract biologically meaningful knowledge. This review paper presents the
applications of formal concept analysis for the analysis and knowledge
discovery from biological data, including gene expression discretization, gene
co-expression mining, gene expression clustering, finding genes in gene
regulatory networks, enzyme/protein classifications, binding site
classifications, and so on. It also presents a list of FCA-based software tools
applied in biological domain and covers the challenges faced so far.


  Adequate read filtering is critical when processing high-throughput data in
marker-gene-based studies. Sequencing errors can cause the mis-clustering of
otherwise similar reads, artificially increasing the number of retrieved
Operational Taxonomic Units (OTUs) and therefore leading to the overestimation
of microbial diversity. Sequencing errors will also result in OTUs that are not
accurate reconstructions of the original biological sequences. Herein we
present a novel and sensitive sequence filtering algorithm that minimizes both
problems by calculating the exact error-probability distribution of a sequence
from its quality scores. In order to validate our method, we quality-filtered
thirty-seven publicly available datasets obtained by sequencing mock and
environmental microbial communities with the Roche 454, Illumina MiSeq and
IonTorrent PGM platforms, and compared our results to those obtained with
previous approaches such as the ones included in mothur, QIIME and UPARSE. Our
algorithm retained substantially more reads than its predecessors, while
resulting in fewer and more accurate OTUs. This improved sensitiveness produced
more faithful representations, both quantitatively and qualitatively, of the
true microbial diversity present in the studied samples. Furthermore, the
method introduced in this work is computationally inexpensive and can be
readily applied in conjunction with any existent analysis pipeline.


  Cancer has been characterized as a constellation of hundreds of diseases
differing in underlying mutations and depending on cellular environments.
Carcinogenesis as a stochastic physical process has been studied for over sixty
years, but there is no accepted standard model. We show that the hazard rates
of all cancers are characterized by a simple dynamic stochastic process on a
half-line, with a universal linear restoring force balancing a universal simple
Brownian motion starting from a universal initial distribution. Only a critical
radius defining the transition from normal to tumorigenic genomes distinguishes
between different cancer types when time is measured in cell--cycle units.
Reparametrizing to chronological time units introduces two additional
parameters: the onset of cellular senescence with age and the time interval
over which this cessation in replication takes place. This universality implies
that there may exist a finite separation between normal cells and tumorigenic
cells in all tissue types that may be a viable target for both early detection
and preventive therapy.


  Predicting the spread of vector-borne diseases in response to incursions
requires knowledge of both host and vector demographics in advance of an
outbreak. Whereas host population data is typically available, for novel
disease introductions there is a high chance of the pathogen utilising a vector
for which data is unavailable. This presents a barrier to estimating the
parameters of dynamical models representing host-vector-pathogen interaction,
and hence limits their ability to provide quantitative risk forecasts. The
Theileria orientalis (Ikeda) outbreak in New Zealand cattle demonstrates this
problem: even though the vector has received extensive laboratory study, a high
degree of uncertainty persists over its national demographic distribution.
Addressing this, we develop a Bayesian data assimilation approach whereby
indirect observations of vector activity inform a seasonal spatio-temporal risk
surface within a stochastic epidemic model. We provide quantitative predictions
for the future spread of the epidemic, quantifying uncertainty in the model
parameters, case infection times, and the disease status of undetected
infections. Importantly, we demonstrate how our model learns sequentially as
the epidemic unfolds, and provides evidence for changing epidemic dynamics
through time. Our approach therefore provides a significant advance in rapid
decision support for novel vector-borne disease outbreaks.


  The Persistent-Phylogeny Model is an extension of the widely studied
Perfect-Phylogeny Model, encompassing a broader range of evolutionary
phenomena. Biological and algorithmic questions concerning persistent phylogeny
have been intensely investigated in recent years. In this paper, we explore two
alternative approaches to the persistent-phylogeny problem that grow out of our
previous work on perfect phylogeny, and on galled trees. We develop an integer
programming solution to the Persistent-Phylogeny Problem; empirically explore
its efficiency; and empirically explore the utility of using fast algorithms
that recognize galled trees, to recognize persistent phylogeny. The empirical
results identify parameter ranges where persistent phylogeny are galled trees
with high frequency, and show that the integer programming approach can
efficiently identify persistent phylogeny of much larger size than has been
previously reported.


  Recent experiments showing scaling of the intrachromosomal contact
probability, $P(s)\sim s^{-1}$ with the genomic distance $s$, are interpreted
to mean a self-similar fractal-like chromosome organization. However, scaling
of $P(s)$ varies across organisms, requiring an explanation. We illustrate
dynamical arrest in a highly confined space as a discriminating marker for
genome organization, by modeling chromosome inside a nucleus as a homopolymer
confined to a sphere of varying sizes. Brownian dynamics simulations show that
the chain dynamics slows down as the polymer volume fraction ($\phi$) inside
the confinement approaches a critical value $\phi_c$. The universal value of
$\phi_c^{\infty}\approx 0.44$ for a sufficiently long polymer ($N\gg 1$) allows
us to discuss genome dynamics using $\phi$ as a single parameter. Our study
shows that the onset of glassy dynamics is the reason for the segregated
chromosome organization in human ($N\approx 3\times 10^9$,
$\phi\gtrsim\phi_c^{\infty}$), whereas chromosomes of budding yeast ($N\approx
10^8$, $\phi<\phi_c^{\infty}$) are equilibrated with no clear signature of such
organization.


  The Verhulst model is probably the best known macroscopic rate equation in
population ecology. It depends on two parameters, the intrinsic growth rate and
the carrying capacity. These parameters can be estimated for different
populations and are related to the reproductive fitness and the competition for
limited resources, respectively. We investigate analytically and numerically
the simplest possible microscopic scenarios that give rise to the logistic
equation in the deterministic mean-field limit. We provide a definition of the
two parameters of the Verhulst equation in terms of microscopic parameters. In
addition, we derive the conditions for extinction or persistence of the
population by employing either the "momentum-space" spectral theory or the
"real-space" Wentzel-Kramers-Brillouin (WKB) approximation to determine the
probability distribution function and the mean time to extinction of the
population. Our analytical results agree well with numerical simulations.


  Induced effects by direct exposure to ionizing radiation (IR) are a central
issue in many fields like radiation protection, clinic diagnosis and
oncological therapies. Direct irradiation at certain doses induce cell death,
but similar effects can also occur in cells no directly exposed to IR, a
mechanism known as bystander effect. Non-IR (radiofrequency waves) can induce
the death of cells loaded with MNPs in a focused oncological therapy known as
magnetic hyperthermia. Indirect mechanisms are also able to induce the death of
unloaded MNPs cells. Using in vitro cell models, we found that colocalization
of the MNPs at the lysosomes and the non-increase of the temperature induces
bystander effect under non-IR. Our results provide a landscape in which
bystander effects are a more general mechanism, up to now only observed and
clinically used in the field of radiotherapy.


  Photosynthesis is a mechanism developed by terrestrial life to utilize the
energy from photons of solar origin for biological use. Subsurface regions are
isolated from the photosphere, and consequently are incapable of utilizing this
energy. This opens up the opportunity for life to evolve alternate mechanisms
for harvesting available energy. Bacterium Candidatus Desulforudis audaxviator,
found 2.8 km deep in a South African mine, harvests energy from radiolysis,
induced by particles emitted from radioactive U, Th and K present in
surrounding rock. Another radiation source in the subsurface environments is
secondary particles generated by Galactic Cosmic Rays (GCRs). Using Monte Carlo
simulations, it is shown that it is a steady source of energy comparable to
that produced by radioactive substances, and the possibility of a slow
metabolizing life flourishing on it cannot be ruled out. Two mechanisms are
proposed through which GCR-induced secondary particles can be utilized for
biological use in subsurface environments: (1) GCRs injecting energy in the
environment through particle-induced radiolysis, and (2) organic synthesis from
GCR secondaries interacting with the medium. Laboratory experiments to test
these hypotheses are also proposed. Implications of these mechanisms on finding
life in the Solar System and elsewhere in the Universe are discussed.


  Joint peak detection is a central problem when comparing samples in genomic
data analysis, but current algorithms for this task are unsupervised and
limited to at most 2 sample types. We propose PeakSegJoint, a new constrained
maximum likelihood segmentation model for any number of sample types. To select
the number of peaks in the segmentation, we propose a supervised penalty
learning model. To infer the parameters of these two models, we propose to use
a discrete optimization heuristic for the segmentation, and convex optimization
for the penalty learning. In comparisons with state-of-the-art peak detection
algorithms, PeakSegJoint achieves similar accuracy, faster speeds, and a more
interpretable model with overlapping peaks that occur in exactly the same
positions across all samples.


  Termites, like many social insects, build nests of complex architecture.
These constructions have been proposed to optimize different structural
features. Here we describe the nest network of the termite Nasutitermes
ephratae, which is among the largest nest-network reported for termites and
show that it optimizes diverse parameters defining the network architecture.
The network structure avoids multiple crossing of galleries and minimizes the
overlap of foraging territories. Thus, these termites are able to minimize the
number of galleries they build, while maximizing the foraging area available at
the nest mounds. We present a simple computer algorithm that reproduces the
basics characteristics of this termite nest network, showing that simple rules
can produce complex architectural designs efficiently.


  This work by Feng and papers which published its conclusions do not cite nor
do they deal with objections by the author published in 2013-2014. There are
many fundamental problems. We will summarize here the principal problems, as
published by Melott and Bambach which render irrelevant most of the work
presented by Feng.


  In the face of uncertain biological response to climate change and the many
critiques concerning model complexity it is increasingly important to develop
predictive mechanistic frameworks that capture the dominant features of
ecological communities and their dependencies on environmental factors. This is
particularly important for critical global processes such as biomass changes,
carbon export, and biogenic climate feedback. Past efforts have successfully
understood a broad spectrum of plant and community traits across a range of
biological diversity and body size, including tree size distributions and
maximum tree height, from mechanical, hydrodynamic, and resource constraints.
Recently it was shown that global scaling relationships for net primary
productivity are correlated with local meteorology and the overall biomass
density within a forest. Along with previous efforts, this highlights the
connection between widely observed allometric relationships and predictive
ecology. An emerging goal of ecological theory is to gain maximum predictive
power with the least number of parameters. Here we show that the explicit
dependence of such critical quantities can be systematically predicted knowing
just the size of the largest tree. This is supported by data showing that
forests converge to our predictions as they mature. Since maximum tree size can
be calculated from local meteorology this provides a general framework for
predicting the generic structure of forests from local environmental parameters
thereby addressing a range of critical Earth-system questions.


  The simplest theories often have much merit and many limitations, and in this
vein, the value of Neutral Theory (NT) has been the subject of much debate over
the past 15 years. NT was proposed at the turn of the century by Stephen
Hubbell to explain pervasive patterns observed in the organization of
ecosystems. Its originally tepid reception among ecologists contrasted starkly
with the excitement it caused among physicists and mathematicians. Indeed, NT
spawned several theoretical studies that attempted to explain empirical data
and predicted trends of quantities that had not yet been studied. While there
are a few reviews of NT oriented towards ecologists, our goal here is to review
the quantitative results of NT and its extensions for physicists who are
interested in learning what NT is, what its successes are and what important
problems remain unresolved. Furthermore, we hope that this review could also be
of interest to theoretical ecologists because many potentially interesting
results are buried in the vast NT literature. We propose to make these more
accessible by extracting them and presenting them in a logical fashion. We
conclude the review by discussing how one might introduce realistic non-neutral
elements into the current models.


  The twin processes of viral evolutionary escape and reversion in response to
host immune pressure, in particular the cytotoxic T-lymphocyte (CTL) response,
shape Human Immunodeficiency Virus-1 sequence evolution in infected host
populations. The tempo of CTL escape and reversion is known to differ between
CTL escape variants in a given host population. Here, we ask: are rates of
escape and reversion comparable across infected host populations? For three
cohorts taken from three continents, we estimate escape and reversion rates at
23 escape sites in optimally defined Gag epitopes. We find consistent escape
rate estimates across the examined cohorts. Reversion rates are also consistent
between a Canadian and South African infected host population. Certain Gag
escape variants that incur a large replicative fitness cost are known to revert
rapidly upon transmission. However, the relationship between escape/reversion
rates and viral replicative capacity across a large number of epitopes has not
been interrogated. We investigate this relationship by examining $in$ $vitro$
replicative capacities of viral sequences with minimal variation: point escape
mutants induced in a lab strain. Remarkably, despite the complexities of
epistatic effects exemplified by pathways to escape in famous epitopes, and the
diversity of both hosts and viruses, CTL escape mutants which escape rapidly
tend to be those with the highest replicative capacity when applied as a single
point mutation. Similarly, mutants inducing the greatest costs to viral
replicative capacity tend to revert more quickly. These data suggest that
escape rates in Gag are consistent across host populations, and that in general
these rates are dominated by site specific effects upon viral replicative
capacity.


  We develop a latent variable model and an efficient spectral algorithm
motivated by the recent emergence of very large data sets of chromatin marks
from multiple human cell types. A natural model for chromatin data in one cell
type is a Hidden Markov Model (HMM); we model the relationship between multiple
cell types by connecting their hidden states by a fixed tree of known
structure. The main challenge with learning parameters of such models is that
iterative methods such as EM are very slow, while naive spectral methods result
in time and space complexity exponential in the number of cell types. We
exploit properties of the tree structure of the hidden states to provide
spectral algorithms that are more computationally efficient for current
biological datasets. We provide sample complexity bounds for our algorithm and
evaluate it experimentally on biological data from nine human cell types.
Finally, we show that beyond our specific model, some of our algorithmic ideas
can be applied to other graphical models.


  Recent work draws attention to community-community encounters ("coalescence")
as likely an important factor shaping natural ecosystems. This work builds on
MacArthur's classic model of competitive coexistence to investigate such
community-level competition in a minimal theoretical setting. It is shown that
the ability of a species to survive a coalescence event is best predicted by a
community-level "fitness" of its native community rather than the intrinsic
performance of the species itself. The model presented here allows formalizing
a macroscopic perspective whereby a community harboring organisms at varying
abundances becomes equivalent to a single organism expressing genes at
different levels. While most natural communities do not satisfy the strict
criteria of multicellularity developed by multi-level selection theory, the
effective cohesion described here is a generic consequence of division of
labor, requires no cooperative interactions, and can be expected to be
widespread in microbial ecosystems.


  Vibrio cholerae is a globally distributed water-borne pathogen that causes
severe diarrheal disease and mortality, with current outbreaks as part of the
seventh pandemic. Further understanding of the role of environmental factors in
potential pathogen distribution and corresponding V. cholerae disease
transmission over time and space is urgently needed to target surveillance of
cholera and other climate and water-sensitive diseases. We used an ecological
niche model (ENM) to identify environmental variables associated with V.
cholerae presence in marine environments, to project a global model of V.
cholerae distribution in ocean waters under current and future climate
scenarios. We generated an ENM using published reports of V. cholerae in
seawater and freely available remotely sensed imagery. Models indicated that
factors associated with V. cholerae presence included chlorophyll-a, pH, and
sea surface temperature (SST), with chlorophyll-a demonstrating the greatest
explanatory power from variables selected for model calibration. We identified
specific geographic areas for potential V. cholerae distribution. Coastal
Bangladesh, where cholera is endemic, was found to be environmentally similar
to coastal areas in Latin America. In a conservative climate change scenario,
we observed a predicted increase in areas with environmental conditions
suitable for V. cholerae. Findings highlight the potential for vulnerability
maps to inform cholera surveillance, early warning systems, and disease
prevention and control.


  Identifying the entire set of complexes is essential not only to understand
complex formations, but also to map the high level organisation of the cell.
Computational prediction of protein complexes faces several challenges
including the lack of sufficient protein interactions, presence of noise in
protein interaction datasets and difficulty in predicting small and sparse
complexes. These challenges are covered in most reviews of complex prediction
methods. However, an important challenge that needs to be addressed is the
prediction of membrane complexes. These are often ignored because existing
protein interaction detection techniques do not detect interactions between
membrane proteins. But, recently there have been several new experimental
techniques including MY2H that are capable of detecting membrane protein
interactions. In the light of this new data, we discuss here new challenges and
the kind of open problems that need to be solved to effectively detect membrane
complexes.


  Transcription factor binding to the surface of DNA regulatory regions is one
of the primary causes of regulating gene expression levels. A probabilistic
approach to model protein-DNA interactions at the sequence level is through
Position Weight Matrices (PWMs) that estimate the joint probability of a DNA
binding site sequence by assuming positional independence within the DNA
sequence.
  Here we construct conditional PWMs that depend on the motif signatures in the
flanking DNA sequence, by conditioning known binding site loci on the presence
or absence of additional binding sites in the flanking sequence of each site's
locus. Pooling known sites with similar flanking sequence patterns allows for
the estimation of the conditional distribution function over the binding site
sequences.
  We apply our model to the Dorsal transcription factor binding sites active in
patterning the Dorsal-Ventral axis of Drosophila development. We find that
those binding sites that cooperate with nearby Twist sites on average contain
about 0.5 bits of information about the presence of Twist transcription factor
binding sites in the flanking sequence. We also find that Dorsal binding site
detectors conditioned on flanking sequence information make better predictions
about what is a Dorsal site relative to background DNA than detection without
information about flanking sequence features.


  We study the mechanical unfolding of a simple model protein. The Langevin
dynamics results are analyzed using Markov-model methods which allow to
describe completely the configurational space of the system. Using transition
path theory we also provide a quantitative description of the unfolding
pathways followed by the system. Our study shows a complex dynamical scenario.
In particular, we see that the usual one-dimensional picture: free-energy vs
end-to-end distance representation, gives a misleading description of the
process. Unfolding can occur following different pathways and configurations
which seem to play a central role in one-dimensional pictures are not the
intermediate states of the unfolding dynamics.


  The problem of unicellular-multicellular transition is one of the main issues
that is discussing in evolutionary biology. In [1] the fitness of a colony of
cells is considered in terms of its two basic components, viability and
fecundity. Intrinsic trade-off function of each cell defines a type of cell. We
elaborate models providing in [1]. Assuming that all intrinsic trade-off
functions are linear, we construct a model with different cell types and show
that the differentiation of these types tends to full specialization. In
addition, we attempt to consider the fact that environmental factors influence
on the fitness of the colony. Thus, we introduce an energy restriction to the
model and show that in optimum we get situations in which there exists a set of
states, each of them allowing colony to achieve the same maximum level of
fitness. In some states arbitrary chosen cell may be specialized, in some -
unspecialized, but fecundity and viability of each cell belong to limited
ranges (which are unique for each cell). It is worth pointing out that the
models from [1] are not robust. We try to overcome this disadvantage.


  Black cumin (Nigella sativa) is a spice having medicinal properties with
pungent and bitter odour. It is used since thousands of years to treat various
ailments, including cancer mainly in South Asia and Middle Eastern regions.
Substantial evidence in multiple research studies emphasizes about the
therapeutic importance of bioactive principles of N. sativa in cancer
bioassays; however, the exact mechanism of their anti-tumour action is still to
be fully comprehended. The current study makes an attempt in this direction by
exploiting the advancements in the Insilico reverse screening technology. In
this study, three different Insilico Reverse Screening approaches have been
employed for identifying the putative molecular targets of the bioactive
principles in Black cumin (thymoquinone, alpha-hederin, dithymoquinone and
thymohydroquinone) relevant to its anti-tumour functionality. The identified
set of putative targets is further compared with the existing set of
experimentally validated targets, so as to estimate the performance of insilico
platforms. Subsequently, molecular docking simulations studies were performed
to elucidate the molecular interactions between the bioactive compounds & their
respective identified targets. The molecular interactions of one such target
identified i.e. VEGF2 along with thymoquinone depicted one H-bond formed at the
catalytic site. The molecular targets identified in this study need further
confirmatory tests on cancer bioassays, in order to justify the research
findings from Insilico platforms. This study has brought to light the
effectiveness of usage of Insilico Reverse Screening protocols to characterise
the un-identified target-ome of poly pharmacological bioactive agents in
spices.