  The aim of this paper is to study the qualitative behaviour predicted by a
mathematical model for the initial stage of T cell activation. The state
variables in the model are the concentrations of phosphorylation states of the
T cell receptor complex and the phosphatase SHP-1 in the cell. It is shown that
these quantities cannot approach zero and that the model possesses more than
one positive steady state for certain values of the parameters. It can also
exhibit damped oscillations. It is proved that the chemical concentration which
represents the degree of activation of the cell, that of the maximally
phosphorylated form of the T cell receptor complex, is in general a
non-monotone function of the activating signal. In particular there are cases
where there is a value of the dissociation constant of the ligand from the
receptor which produces an optimal activation of the T cell. In this way the
results of certain simulations in the literature have been confirmed rigorously
and some important features which had not previously been seen have been
discovered.


  In an adapted population of mutators in which most mutations are deleterious,
a nonmutator that lowers the mutation rate is under indirect selection and can
sweep to fixation. Using a multitype branching process, we calculate the
fixation probability of a rare nonmutator in a large population of asexual
mutators. We show that when beneficial mutations are absent, the fixation
probability is a nonmonotonic function of the mutation rate of the mutator: it
first increases sublinearly and then decreases exponentially. We also find that
beneficial mutations can enhance the fixation probability of a nonmutator. Our
analysis is relevant to an understanding of recent experiments in which a
reduction in the mutation rates has been observed.


  Technical signs of progress during the last decades has led to a situation in
which the accumulation of genome sequence data is increasingly fast and cheap.
The huge amount of molecular data available nowadays can help addressing new
and essential questions in Evolution. However, reconstructing evolution of DNA
sequences requires models, algorithms, statistical and computational methods of
ever increasing complexity. Since most dramatic genomic changes are caused by
genome rearrangements (gene duplications, gain/loss events), it becomes crucial
to understand their mechanisms and reconstruct ancestors of the given genomes.
This problem is known to be NP-complete even in the "simplest" case of three
genomes. Heuristic algorithms are usually executed to provide approximations of
the exact solution.
  We state that, even if the ancestral reconstruction problem is NP-hard in
theory, its exact resolution is feasible in various situations, encompassing
organelles and some bacteria. Such accurate reconstruction, which identifies
too some highly homoplasic mutations whose ancestral status is undecidable,
will be initiated in this work-in-progress, to reconstruct ancestral genomes of
two Mycobacterium pathogenetic bacterias. By mixing automatic reconstruction of
obvious situations with human interventions on signaled problematic cases, we
will indicate that it should be possible to achieve a concrete, complete, and
really accurate reconstruction of lineages of the Mycobacterium tuberculosis
complex. Thus, it is possible to investigate how these genomes have evolved
from their last common ancestors.


  We first revisit the multi-allelic mutation-fitness balance problem,
especially when mutations obey a house of cards condition, where the
discrete-time deterministic evolutionary dynamics of the allelic frequencies
derives from a Shahshahani potential. We then consider multi-allelic
Wright-Fisher stochastic models whose deviation to neutrality is from the
Shahsha-hani mutation/selection potential. We next focus on the weak selection,
weak mutation cases and, making use of a Gamma calculus, we compute the
normalizing partition functions of the invariant probability densities
appearing in their Wright-Fisher diffusive approximations. Using these results,
Generalized Ewens sampling formulae (ESF) from the equilibrium distributions
are derived. We start treating the ESF in the mixed mutation/selection
potential case and then we restrict ourselves to the ESF in the simpler
house-of-cards mutations only situation. We also address some issues concerning
sampling problems from infinitely-many alleles weak limits.


  A fundamental question in systems biology is what combinations of mean and
variance of the species present in a stochastic biochemical reaction network
are attainable by perturbing the system with an external signal. To address
this question, we show that the moments evolution in any generic network can be
either approximated or, under suitable assumptions, computed exactly as the
solution of a switched affine system. Motivated by this application, we propose
a new method to approximate the reachable set of switched affine systems. A
remarkable feature of our approach is that it allows one to easily compute
projections of the reachable set for pairs of moments of interest, without
requiring the computation of the full reachable set, which can be prohibitive
for large networks. As a second contribution, we also show how to select the
external signal in order to maximize the probability of reaching a target set.
To illustrate the method we study a renown model of controlled gene expression
and we derive estimates of the reachable set, for the protein mean and
variance, that are more accurate than those available in the literature and
consistent with experimental data.


  Divergence between populations for a given trait can be driven by natural or
sexual selection, interacting with migration behaviour. Mating preference for
different phenotypes can lead to the emergence and persistence of
differentiated populations. Dominance between alleles encoding for divergent
phenotypes can interfere in such processes. Using a diploid model of trait
determining both mating success and migration rate, we explored differentiation
between two connected populations, assuming either co-dominance or strict
dominance between alleles. The model assumes that individuals prefer mating
with partners displaying the same phenotype and therefore tend to move to the
other population when their own phenotype is rare. We show that the emergence
of differentiated populations in this diploid model is limited as compared to
results obtained with the same model assuming haploidy. When assuming
co-dominance, differentiation arises only when migration is limited as compared
to preference. Such differentiation is less dependent on migration when
assuming strict dominance between haplotypes. Dominant alleles frequently
invade populations because their phenotype is more frequently expressed,
resulting in higher mating success and rapid decrease in migration. However,
depending on the initial distribution of alleles, this advantage associated
with dominance (i.e. Haldane's sieve) may lead to fixation of the dominant
allele throughout both populations. Depending on the initial distribution of
heterozygotes, persistence of polymorphisms within populations can also occur
because heterozygotes displaying the predominant phenotype benefit from mating
preferences. Altogether, our results highlight that heterozygotes' behaviour
has a strong impact on population differentiation and stress out the need of
diploid models of differentiation and speciation driven by natural and sexual
selection.


  Information is increasingly being viewed as a resource used by organisms to
increase their fitness. Indeed, it has been formally shown that there is a
sensible way to assign a reproductive value to information and it is
non-negative. However, all of this work assumed that information collection is
cost-free. Here, we account for such a cost and provide conditions for when the
reproductive value of information will be negative. In these instances,
counter-intuitively, it is in the interest of the organism to remain ignorant.
We link our results to empirical studies where Bayesian behaviour appears to
break down in complex environments and provide an alternative explanation of
lowered arousal thresholds in the evolution of sleep.


  Ebola virus is one of the most virulent pathogens for humans. We present a
mathematical description of different Susceptible-Exposed-Infectious-Recovered
(SEIR) models. By using mathematical modeling and analysis, the latest major
outbreak of Ebola virus in West Africa is described. Our aim is to study and
discuss the properties of SEIR models with respect to Ebola virus, the
information they provide, and when the models make sense. We added to the basic
SEIR model demographic effects in order to analyze the equilibria with vital
dynamics. Numerical simulations confirm the theoretical analysis. The control
of the propagation of the virus through vaccination is investigated and the
case study of Liberia is discussed in detail.


  The host immune response can often efficiently suppress a virus infection,
which may lead to selection for immune-resistant viral variants within the
host. For example, during HIV infection, an array of CTL immune response
populations recognize specific epitopes (viral proteins) presented on the
surface of infected cells to effectively mediate their killing. However HIV can
rapidly evolve resistance to CTL attack at different epitopes, inducing a
dynamic network of interacting viral and immune response variants. We consider
models for the network of virus and immune response populations, consisting of
Lotka-Volterra-like systems of ordinary differential equations. Stability of
feasible equilibria and corresponding uniform persistence of distinct variants
are characterized via a Lyapunov function. We specialize the model to a "binary
sequence" setting, where for $n$ epitopes there can be $2^n$ distinct viral
variants mapped on a hypercube graph. The dynamics in several cases are
analyzed and sharp polychotomies are derived characterizing persistent
variants. In particular, we prove that if the viral fitness costs for gaining
resistance to each epitope are equal, then the system of $2^n$ virus strains
converges to a "perfectly nested network" with less than or equal to $n+1$
persistent virus strains. Overall, our results suggest that immunodominance,
i.e. relative strength of immune response to an epitope, is the most important
factor determining the persistent network structure.


  Sewall Wright's adaptive landscape metaphor penetrates a significant part of
evolutionary thinking. Supplemented with Fisher's fundamental theorem of
natural selection and Kimura's maximum principle, it provides a unifying and
intuitive representation of the evolutionary process under the influence of
natural selection as the hill climbing on the surface of mean population
fitness. On the other hand, it is also well known that for many more or less
realistic mathematical models this picture is a sever misrepresentation of what
actually occurs. Therefore, we are faced with two questions. First, it is
important to identify the cases in which adaptive landscape metaphor actually
holds exactly in the models, that is, to identify the conditions under which
system's dynamics coincides with the process of searching for a (local) fitness
maximum. Second, even if the mean fitness is not maximized in the process of
evolution, it is still important to understand the structure of the mean
fitness manifold and see the implications of this structure on the system's
dynamics. Using as a basic model the classical replicator equation, in this
note we attempt to answer these two questions and illustrate our results with
simple well studied systems.


  We modify the Tangled Nature Model of Christensen et. al. so that the agents
affect the carrying capacity. This leads to a model of species-environment
co-evolution where the system tends to have a larger carrying capacity with
life than without. We discuss the model as an example of an entropic hierarchy
and some implications for Gaia theory.


  We have earlier reported the MOLSDOCK technique to perform rigid
receptor/flexible ligand docking. The method uses the MOLS method, developed in
our laboratory. In this paper we report iMOLSDOCK, the 'flexible receptor'
extension we have carried out to the algorithm MOLSDOCK. iMOLSDOCK uses
mutually orthogonal Latin squares (MOLS) to sample the conformation and the
docking pose of the ligand and also the flexible residues of the receptor
protein. The method then uses a variant of the mean field technique to analyze
the sample to arrive at the optimum. We have benchmarked and validated
iMOLSDOCK with a dataset of 44 peptide-protein complexes with peptides. We have
also compared iMOLSDOCK with other flexible receptor docking tools GOLD v5.2.1
and AutoDock Vina. The results obtained show that the method works better than
these two algorithms, though it consumes more computer time.


  RNA sequencing allows one to study allelic imbalance of gene expression,
which may be due to genetic factors or genomic imprinting. It is desirable to
model both genetic and parent-of-origin effects simultaneously to avoid
confounding and to improve the power to detect either effect. In a study of
experimental cross, separation of genetic and parent-of-origin effects can be
achieved by studying reciprocal cross of two inbred strains. In contrast, this
task is much more challenging for an outbred population such as human
population. To address this challenge, we propose a new framework to combine
experimental strategies and novel statistical methods. Specifically, we propose
to collect genotype data from family trios as well as RNA-seq data from the
children of family trios. We have developed a new statistical method to
estimate both genetic and parent-of-origin effects from such data sets. We
demonstrated this approach by studying 30 trios of HapMap samples. Our results
support some of previous finding of imprinted genes and also recover new
candidate imprinted genes.


  Local Ca Releases (LCRs) are crucial events involved in cardiac pacemaker
cell function. However, specific algorithms for automatic LCR detection and
analysis have not been developed in live, spontaneously beating pacemaker
cells. Here we measured LCRs using a high-speed 2D-camera in spontaneously
contracting sinoatrial (SA) node cells isolated from rabbit and guinea pig and
developed a new algorithm capable of detecting and analyzing the LCRs spatially
in two-dimensions, and in time. Our algorithm tracks points along the midline
of the contracting cell. It uses these points as a coordinate system for affine
transform, producing a transformed image series where the cell does not
contract. Action potential-induced Ca transients and LCRs were thereafter
isolated from recording noise by applying a series of spatial filters. The LCR
birth and death events were detected by a differential (frame-to-frame)
sensitivity algorithm. An LCR was detected when its signal changes sufficiently
quickly within a sufficiently large area. The LCR is considered to have died
when its amplitude decays substantially, or when it merges into the rising
whole cell Ca transient. Our algorithm provides major LCR parameters such as
period, signal mass, duration, and path area. As LCRs propagate within cells,
the algorithm identifies splitting and merging behaviors, indicating the
importance of Ca-induced-Ca-release for the fate of LCRs and for generating a
powerful ensemble Ca signal. Thus, our new computer algorithms eliminate motion
artifacts and detect 2D local spatiotemporal Ca release events from recording
noise and global signals. While the algorithms detect LCRs in sinoatrial nodal
cells, they have the potential to be used in other applications in biophysics
and cell physiology, for example, to detect Ca wavelets (abortive waves),
sparks and embers in muscle cells and Ca puffs and syntillas in neurons.


  We introduce a powerful iterative algorithm to compute protein folding
pathways, with realistic all-atom force fields. Using the path integral
formalism, we explicitly derive a modified Langevin equation which samples
directly the ensemble of reactive pathways, exponentially reducing the cost of
simulating thermally activated transitions. The algorithm also yields a
rigorous stochastic estimate of the reaction coordinate. After illustrating
this approach on a simple toy model, we successfully validate it against the
results of ultra-long plain MD protein folding simulations for a fast folding
protein (Fip35), which were performed on the Anton supercomputer. Using our
algorithm, computing a folding trajectory for this protein requires only 1000
core hours, a computational load which could be even carried out on a desktop
workstation.


  Phenotypical variability in the absence of genetic variation often reflects
complex energetic landscapes associated with underlying gene regulatory
networks (GRNs). In this view, different phenotypes are associated with
alternative states of complex nonlinear systems: stable attractors in
deterministic models or modes of stationary distributions in stochastic
descriptions. We provide theoretical and practical characterizations of these
landscapes, specifically focusing on stochastic slow promoter kinetics, a time
scale relevant when transcription factor binding and unbinding are affected by
epigenetic processes like DNA methylation and chromatin remodeling. In this
case, largely unexplored except for numerical simulations, adiabatic
approximations of promoter kinetics are not appropriate. In contrast to the
existing literature, we provide rigorous analytic characterizations of multiple
modes. A general formal approach gives insight into the influence of parameters
and the prediction of how changes in GRN wiring, for example through mutations
or artificial interventions, impact the possible number, location, and
likelihood of alternative states. We adapt tools from the mathematical field of
singular perturbation theory to represent stationary distributions of Chemical
Master Equations for GRNs as mixtures of Poisson distributions and obtain
explicit formulas for the locations and probabilities of metastable states as a
function of the parameters describing the system. As illustrations, the theory
is used to tease out the role of cooperative binding in stochastic models in
comparison to deterministic models, and applications are given to various model
systems, such as toggle switches in isolation or in communicating populations,
a synthetic oscillator and a trans-differentiation network.


  The raw data from study of variation of D-loop mtDNA of white sturgeon
Acipenser transmontanus (Mol. Biol. Evol. 1993. 10: 326-341) was re-analyzed.
Re-calculated nucleotide substitution rate ({\mu}) was 0.782-0.939 x 10-7
substitutions/site/year/lineage, which was 1.4 times less than the estimate
given in above-mentioned publication. The use of new {\mu} has led to an
increase in estimates of long-term female effective population size (Nef) and
coalescence times for mtDNA haplotypes, previously calculated for samples of
Amur sturgeon A. schrenckii and kaluga A. dauricus, but it did not affect the
conclusions on the critical state of their natural populations.


  We report here new electrical laws, derived from nonlinear electro-diffusion
theory, about the effect of the local geometrical structure, such as curvature,
on the electrical properties of a cell. We adopt the Poisson-Nernst-Planck
(PNP) equations for charge concentration and electric potential as a model of
electro-diffusion. In the case at hand, the entire boundary is impermeable to
ions and the electric field satisfies the compatibility condition of Poisson's
equation. We construct an asymptotic approximation for certain singular limits
to the steady-state solution in a ball with an attached cusp-shaped funnel on
its surface. As the number of charge increases, they concentrate at the end of
cusp-shaped funnel. These results can be used in the design of nano-pipettes
and help to understand the local voltage changes inside dendrites and axons
with heterogenous local geometry.


  How important is resource depletion (e.g. depletion of target cells) in
controlling infection within a host? And how can we distinguish between
resource depletion and other mechanisms that may contribute to decline of
pathogen load or lead to pathogen clearance? In this paper we examine data from
a previously published experiment. In this experiment, mice were infected with
influenza virus carrying a green fluorescent protein reporter gene, and the
proportion of lung epithelial cells that were influenza infected was measured
as a function of time. Three inoculum dose groups - 10^4 PFU, 10^6 PFU and 10^7
PFU - were used. The proportion of cells infected was estimated to be about 21
(95% confidence interval: 14-32) fold higher in the highest dose group than in
the lowest dose group with the middle dose group in between. We show that this
pattern is highly inconsistent with a model where target cell depletion is the
principal means of controlling infection, and we argue that such a pattern
constitutes a reasonable criterion for rejecting many resource depletion
models. A model with an innate interferon response that renders susceptible
cells resistant fits the data reasonably well. This model suggests that target
cell depletion is only a minor factor in controlling natural influenza
infection.


  Recent outbreaks of Ebola, H1N1 and other infectious diseases have shown that
the assumptions underlying the established theory of epidemics management are
too idealistic. For an improvement of procedures and organizations involved in
fighting epidemics, extended models of epidemics management are required. The
necessary extensions consist in a representation of the management loop and the
potential frictions influencing the loop. The effects of the non-deterministic
frictions can be taken into account by including the measures of robustness and
risk in the assessment of management options. Thus, besides of the increased
structural complexity resulting from the model extensions, the computational
complexity of the task of epidemics management - interpreted as an optimization
problem - is increased as well. This is a serious obstacle for analyzing the
model and may require an additional pre-processing enabling a simplification of
the analysis process. The paper closes with an outlook discussing some
forthcoming problems.


  We present DiSH-Sim, a simulator for large discrete models of biological
signal transduction pathways, capable of simulating networks with multi-valued
elements in both deterministic and stochastic manner. We focus on order of
update and thus incorporate information about timing, taking into account that
biological processes are not synchronized and certain biochemical changes occur
slower than others. Another feature of our simulator is the use of grouped
rules to model multi-valued elements and delays. The DiSH-Sim is publicly
available and is being used to validate discrete cancer microenvironment and
infectious disease models. It is also incorporated within a large architecture
that includes natural language processing tools that read biological literature
to assemble logical models. This paper demonstrates the functionalities and
ease of use of DiSH-Sim, making it a very useful tool for discrete modeling.


  The spread of an infectious disease can, in some cases, promote the
propagation of other pathogens favouring violent outbreaks, which cause a
discontinuous transition to an endemic state. The topology of the contact
network plays a crucial role in these cooperative dynamics. We consider a
susceptible--infected--removed (SIR) type model with two mutually cooperative
pathogens: an individual already infected with one disease has an increased
probability of getting infected by the other. We present an heterogeneous
mean-field theoretical approach to the co--infection dynamics on generic
uncorrelated power-law degree-distributed networks and validate its results by
means of numerical simulations. We show that, when the second moment of the
degree distribution is finite, the epidemic transition is continuous for low
cooperativity, while it is discontinuous when cooperativity is sufficiently
high. For scale-free networks, i.e. topologies with diverging second moment,
the transition is instead always continuous. In this way we clarify the effect
of heterogeneity and system size on the nature of the transition and we
validate the physical interpretation about the origin of the discontinuity.


  The increasing availability of sequencing data of cancer samples is fueling
the development of algorithmic strategies to investigate tumor heterogeneity
and infer reliable models of cancer evolution. We here build up on previous
works on cancer progression inference from genomic alteration data, to deliver
two distinct Cytoscape-based applications, which allow to produce, visualize
and manipulate cancer evolution models, also by interacting with public genomic
and proteomics databases. In particular, we here introduce cyTRON, a
stand-alone Cytoscape app, and cyTRON/JS, a web application which employs the
functionalities of Cytoscape/JS.
  cyTRON was developed in Java; the code is available at
https://github.com/BIMIB-DISCo/cyTRON and on the Cytoscape App Store
http://apps.cytoscape.org/apps/cytron. cyTRON/JS was developed in JavaScript
and R; the source code of the tool is available at
https://github.com/BIMIB-DISCo/cyTRON-js and the tool is accessible from
https://bimib.disco.unimib.it/cytronjs/welcome.


  Life science is entering a new era of petabyte-level sequencing data.
Converting such big data to biological insights represents a huge challenge for
computational analysis. To this end, we developed DeepMetabolism, a
biology-guided deep learning system to predict cell phenotypes from
transcriptomics data. By integrating unsupervised pre-training with supervised
training, DeepMetabolism is able to predict phenotypes with high accuracy
(PCC>0.92), high speed (<30 min for >100 GB data using a single GPU), and high
robustness (tolerate up to 75% noise). We envision DeepMetabolism to bridge the
gap between genotype and phenotype and to serve as a springboard for
applications in synthetic biology and precision medicine.


  The spreading and regulation of epigenetic marks on chromosomes is crucial to
establish and maintain cellular identity. Nonetheless, the dynamical mechanism
leading to the establishment and maintenance of a given, cell-line specific,
epigenetic pattern is still poorly understood. In this work we propose, and
investigate in silico, a possible experimental strategy to illuminate the
interplay between 3D chromatin structure and epigenetic dynamics. We consider a
set-up where a reconstituted chromatin fibre is stretched at its two ends
(e.g., by laser tweezers), while epigenetic enzymes (writers) and
chromatin-binding proteins (readers) are flooded into the system. We show that,
by tuning the stretching force and the binding affinity of the readers for
chromatin, the fibre undergoes a sharp transition between a stretched,
epigenetically disordered, state and a crumpled, epigenetically coherent, one.
We further investigate the case in which a knot is tied along the chromatin
fibre, and find that the knotted segment enhances local epigenetic order,
giving rise to "epigenetic solitons" which travel and diffuse along chromatin.
Our results point to an intriguing coupling between 3D chromatin topology and
epigenetic dynamics, which may be investigated via single molecule experiments.


  The interaction between proteins and DNA is a key driving force in a
significant number of biological processes such as transcriptional regulation,
repair, recombination, splicing, and DNA modification. The identification of
DNA-binding sites and the specificity of target proteins in binding to these
regions are two important steps in understanding the mechanisms of these
biological activities. A number of high-throughput technologies have recently
emerged that try to quantify the affinity between proteins and DNA motifs.
Despite their success, these technologies have their own limitations and fall
short in precise characterization of motifs, and as a result, require further
downstream analysis to extract useful and interpretable information from a
haystack of noisy and inaccurate data. Here we propose MotifMark, a new
algorithm based on graph theory and machine learning, that can find binding
sites on candidate probes and rank their specificity in regard to the
underlying transcription factor. We developed a pipeline to analyze
experimental data derived from compact universal protein binding microarrays
and benchmarked it against two of the most accurate motif search methods. Our
results indicate that MotifMark can be a viable alternative technique for
prediction of motif from protein binding microarrays and possibly other related
high-throughput techniques.


  Archaeological settlement systems are usually analysed in terms of the
relation between the rank of a settlement and its size (Zipf's law). We argue
that this approach is unreliable, and can be misleading, in application to
archaeological data where the recovery rate of settlements can be low and their
size estimates are often approximate at best. An alternative framework for the
settlement data interpretation, better connected with the theoretical concepts
of the stochastic evolution of settlement systems, is applied to the evolving
system of settlements of the Late Neolithic-Bronze Age Trypillia cultural
complex (5,400-2,800 BC) in modern Ukraine. The stochastic evolution model
provides a consistent and accurate explanation of the frequency of occurrence
of Trypillian settlement areas in the range from 0.05 to 500 ha. Thus
validated, the model leads to reliable estimates of the typical size of a newly
formed settlement as well as the growth rates of the total number of
settlements and their areas. The parameters of the settlement system thus
revealed are consistent with palaeoeconomy reconstructions for the Trypillia
area.


  Bats are unique mammals. This note discusses some questions regarding bat
evolution including why they are nocturnal and why they can echolocate. It is
hypothesized that echolocation was necessary for bats to survive the period of
limited visibility that followed the Cretaceous-Paleogene (K-Pg) extinction
event.


  The recent development of single-cell transcriptomics has enabled gene
expression to be measured in individual cells instead of being
population-averaged. Despite this considerable precision improvement, inferring
regulatory networks remains challenging because stochasticity now proves to
play a fundamental role in gene expression. In particular, mRNA synthesis is
now acknowledged to occur in a highly bursty manner. We propose to view the
inference problem as a fitting procedure for a mechanistic gene network model
that is inherently stochastic and takes not only protein, but also mRNA levels
into account. We first explain how to build and simulate this network model
based upon the coupling of genes that are described as piecewise-deterministic
Markov processes. Our model is modular and can be used to implement various
biochemical hypotheses including causal interactions between genes. However, a
naive fitting procedure would be intractable. By performing a relevant
approximation of the stationary distribution, we derive a tractable procedure
that corresponds to a statistical hidden Markov model with interpretable
parameters. This approximation turns out to be extremely close to the
theoretical distribution in the case of a simple toggle-switch, and we show
that it can indeed fit real single-cell data. As a first step toward inference,
our approach was applied to a number of simple two-gene networks simulated in
silico from the mechanistic model and satisfactorily recovered the original
networks. Our results demonstrate that functional interactions between genes
can be inferred from the distribution of a mechanistic, dynamical stochastic
model that is able to describe gene expression in individual cells. This
approach seems promising in relation to the current explosion of single-cell
expression data.


  Important effects of noise on a one-dimensional gene expression model
involving a single gene have recently been discussed. However, few works have
been devoted to the transition in two-dimensional models which include the
interaction of genes. Therefore, we investigate here, a quantitative
two-dimensional model (MeKS network) of gene expression dynamics describing the
competence development in the B. subtilis under the influence of L\'evy as well
as Brownian motions, where noises can do the B. subtilis a favor in nutrient
depletion. To analyze the transitions between the vegetative and the competence
regions therein, two deterministic quantities, the mean first exit time (MFET)
and the first escape probability (FEP) from a microscopic perspective, as well
as their averaged versions from a macroscopic perspective, are applied. The
relative contribution factor (RCF), the ratio of non-Gaussian and Gaussian
noise strengths, is adopted to implement optimal control in these transitions.
Schematic representations indicate that there exists an optimum choice that
makes the transition occurring at the highest probability. Additionally, we use
a geometric concept, the stochastic basin of attraction, to exhibit a pictorial
comprehension about the influence of the L\'{e}vy motion on the basin stability
of the competence state.


  Molecular dynamics (MD) simulations allow the exploration of the phase space
of biopolymers through the integration of equations of motion of their
constituent atoms. The analysis of MD trajectories often relies on the choice
of collective variables (CVs) along which the dynamics of the system is
projected. We developed a graphical user interface (GUI) for facilitating the
interactive choice of the appropriate CVs. The GUI allows: defining
interactively new CVs; partitioning the configurations into microstates
characterized by similar values of the CVs; calculating the free energies of
the microstates for both unbiased and biased (metadynamics) simulations;
clustering the microstates in kinetic basins; visualizing the free energy
landscape as a function of a subset of the CVs used for the analysis. A simple
mouse click allows one to quickly inspect structures corresponding to specific
points in the landscape.


  According to the competitive exclusion principle, in a finite ecosystem,
extinction occurs naturally when two or more species compete for the same
resources. An important question that arises is: when coexistence is not
possible, which mechanisms confer an advantage to a given species against the
other(s)? In general, it is expected that the species with the higher
reproductive/death ratio will win the competition, but other mechanisms, such
as asymmetry in interspecific competition or unequal diffusion rates, have been
found to change this scenario dramatically. In this work, we examine
competitive advantage in the context of quasi-neutral population models,
including stochastic models with spatial structure as well as macroscopic
(mean-field) descriptions. We employ a two-species contact process in which the
"biological clock" of one species is a factor of $\alpha$ slower than that of
the other species. Our results provide new insights into how stochasticity and
competition interact to determine extinction in finite spatial systems. We find
that a species with a slower biological clock has an advantage if resources are
limited, winning the competition against a species with a faster clock, in
relatively small systems. Periodic or stochastic environmental variations also
favor the slower species, even in much larger systems.


  Mathematical models of real life phenomena are highly nonlinear involving
multiple parameters and often exhibiting complex dynamics. Experimental data
sets are typically small and noisy, rendering estimation of parameters from
such data unreliable and difficult. This paper presents a study of the Bayesian
posterior distribution for unknown parameters of two chaotic discrete dynamical
systems conditioned on observations of the system. The study shows how the
qualitative properties of the posterior are affected by the intrinsic noise
present in the data, the representation of this noise in the parameter
estimation process, and the length of the data-set. The results indicate that
increasing length of dataset does not significantly increase the precision of
the estimate, and this is true for both periodic and chaotic data. On the other
hand, increasing precision of the measurements leads to significant increase in
precision of the estimated parameter in case of periodic data, but not in the
case of chaotic data. These results are highly useful in designing laboratory
and field-based studies in biology in general, and ecology and conservation in
particular.


  Background: Tree reconciliation problems have long been studied in
phylogenetics. A particular variant of the reconciliation problem for a gene
tree T and a species tree S assumes that for each interior vertex x of T it is
known whether x represents a speciation or a duplication. This problem appears
in the context of analyzing orthology data.
  Results: We show that S is a species tree for T if and only if S displays all
rooted triples of T that have three distinct species as their leaves and are
rooted in a speciation vertex. A valid reconciliation map can then be found in
polynomial time. Simulated data shows that the event-labeled gene trees convey
a large amount of information on underlying species trees, even for a large
percentage of losses.
  Conclusions: The knowledge of event labels in a gene tree strongly constrains
the possible species tree and, for a given species tree, also the possible
reconciliation maps. Nevertheless, many degrees of freedom remain in the space
of feasible solutions. In order to disambiguate the alternative solutions
additional external constraints as well as optimization criteria could be
employed.


  Evolutionary dynamics experienced by mixed microbial populations of
cooperators and cheaters has been examined in experiments in the literature
using a protocol of periodic dilution to investigate the properties of
resilience and adaptability to environmental changes. Data depicted on an
appropriate phase diagram indicate, among other features, a stable equilibrium
point at which cooperators and cheaters coexist [A. Sanchez, J. Gore, PLOs
Biology, 11 (4), e1001547 (2013)]. We present here a phenomenological analysis
of these data focusing on an eco-evolutionary-game perspective. To that end, we
work on an extension of the model proposed by Tao and Cressman Y. Tao, R.
Cressman, Bull. Math. Biol. 69, 1377 - 1399 (2007). It's original version takes
into account changes of the total population density while the individuals
experience pairwise Prisoner's Dilemma game. The extension devised here
contains a dilution factor to be conform with the experimental procedure, in
addition of a term accounting for Allee effects. Differently from other
descriptions proposed in similar contexts, however, the model here does not
account for assortative encounters, group or kin selection. Nonetheless, it
describes surprisingly well both qualitatively and quantitatively the features
of the observed phase diagram. We discuss these results in terms of the
behavior of an effective payoff matrix defined accordingly.


  Models of biological processes are often subject to different sources of
noise. Developing an understanding of the combined effects of different types
of uncertainty is an open challenge. In this paper, we study a variant of the
susceptible-infective-recovered model of epidemic spread, which combines both
agent-to-agent heterogeneity and intrinsic noise. We focus on epidemic cycles,
driven by the stochasticity of infection and recovery events, and study in
detail how heterogeneity in susceptibilities and propensities to pass on the
disease affects these quasi-cycles. While the system can only be described by a
large hierarchical set of equations in the transient regime, we derive a
reduced closed set of equations for population-level quantities in the
stationary regime. We analytically obtain the spectra of quasi-cycles in the
linear-noise approximation. We find that the characteristic frequency of these
cycles is typically determined by population averages of susceptibilities and
infectivities, but that their amplitude depends on higher-order moments of the
heterogeneity. We also investigate the synchronisation properties and phase lag
between different groups of susceptible and infected individuals.


  Hfq is a bacterial protein that is involved in several aspects of nucleic
acids metabolism. It has been described as one of the nucleoid associated
proteins shaping the bacterial chromosome, although it is better known to
influence translation and turnover of cellular RNAs. Here, we explore the role
of Escherichia coli Hfq C-terminal domain in the compaction of double stranded
DNA. Various experimental methodologies, including fluorescence microscopy
imaging of single DNA molecules confined inside nanofluidic channels, atomic
force microscopy, isothermal titration microcalorimetry, and electrophoretic
mobility assays have been used to follow the assembly of the C-terminal and
N-terminal regions of Hfq on DNA. Results highlight the role of Hfq C-terminal
arms in DNA binding, change in mechanical properties of the double helix and
compaction of DNA into a condensed form. The propensity for bridging and
compaction of DNA by the C-terminal domain might be related to aggregation of
bound protein and may have implications for protein binding related gene
regulation.


  The potential number of drug like small molecules is estimated to be between
10^23 and 10^60 while current databases of known compounds are orders of
magnitude smaller with approximately 10^8 compounds. This discrepancy has led
to an interest in generating virtual libraries using hand crafted chemical
rules and fragment based methods to cover a larger area of chemical space and
generate chemical libraries for use in in silico drug discovery endeavors. Here
it is explored to what extent a recurrent neural network with long short term
memory cells can figure out sensible chemical rules and generate synthesizable
molecules by being trained on existing compounds encoded as SMILES. The
networks can to a high extent generate novel, but chemically sensible
molecules. The properties of the molecules are tuned by training on two
different datasets consisting of fragment like molecules and drug like
molecules. The produced molecules and the training databases have very similar
distributions of molar weight, predicted logP, number of hydrogen bond
acceptors and donors, number of rotatable bonds and topological polar surface
area when compared to their respective training sets. The compounds are for the
most cases synthesizable as assessed with SA score and Wiley ChemPlanner.


  Experiments measuring DNA extension in nanochannels are at odds with even the
most basic predictions of current scaling arguments for the conformations of
confined semiflexible polymers such as DNA. We show that a theory based on a
weakly self-avoiding, one-dimensional "telegraph" process collapses
experimental data and simulation results onto a single master curve throughout
the experimentally relevant region of parameter space and explains the
mechanisms at play.


  Whether or not biodiversity dynamics tend toward stable equilibria remains an
unsolved question in ecology and evolution with important implications for our
understanding of diversity and its conservation. Phylo/population genetic
models and macroecological theory represent two primary lenses through which we
view biodiversity. While phylo/population genetics provide an averaged view of
changes in demography and diversity over timescales of generations to
geological epochs, macroecology provides an ahistorical description of
commonness and rarity across contemporary co-occurring species. Our goal is to
combine these two approaches to gain novel insights into the non-equilibrium
nature of biodiversity. We help guide near future research with a call for
bioinformatic advances and an outline of quantitative predictions made possible
by our approach.


  In this review, we discuss recent advances on the plasticity of cancer stem
cells and highlight their relevance to understand the metastatic process and to
guide therapeutic interventions. Recent results suggest that the strict
hierarchical structure of cancer cell populations advocated by the cancer stem
cell model must be reconsidered since the depletion of cancer stem cells leads
the other tumor cells to switch back into the cancer stem cell phenotype. This
plasticity has important implications for metastasis since migrating cells do
not need to be cancer stem cells in order to seed a metastasis. We also discuss
the important role of the immune system and the microenvironment in modulating
phenotypic switching and suggest possible avenues to exploit our understanding
of this process to develop an effective strategy for precision medicine.


  In many dryland environments, vegetation self-organizes into bands that can
be clearly identified in remotely-sensed imagery. The status of individual
bands can be tracked over time, allowing for a detailed remote analysis of how
human populations affect the vital balance of dryland ecosystems. In this
study, we characterize vegetation change in areas of the Horn of Africa where
imagery taken in the early 1950s is available. We find that substantial change
is associated with steep increases in human activity, which we infer primarily
through the extent of road and dirt track development. A seemingly paradoxical
signature of human impact appears as an increase in the widths of the
vegetation bands, which effectively increases the extent of vegetation cover in
many areas. We show that this widening occurs due to altered rates of
vegetation colonization and mortality at the edges of the bands, and conjecture
that such changes are driven by human-induced shifts in plant species
composition. Our findings suggest signatures of human impact that may aid in
identifying and monitoring vulnerable drylands in the Horn of Africa.


  Two basic features of assemblages of unicellular plankton: (1) their high
biodiversity and (2) the power-law structure of their abundance, can be
explained by an allometric scaling of cell growth and mortality with respect to
cell size. To show this, we describe a numerical study of a size-structured,
multispecies, population-dynamic model; the model has a single resource,
supporting an arbitrary number of phytoplankton and zooplankton species. If the
number of plankton species is large enough, the death rate of prey and cell
growth rate of predators have approximate allometric scalings with cell size.
Together, these scalings give rise to an equilibrium distribution of abundance
near the power law, on which many species can coexist. Scalings of this kind
cannot be achieved if the number of species is small. This suggests that the
conjunction of species-richness and power-law structures in plankton
communities is more than a coincidence. Although the exact allometric scalings
used here should not be expected in practice, exclusion of species should be
relatively slow if they lie close to the power law. Thus the forces needed to
achieve coexistence could be effective, even if they are relatively weak.


  The cryptocurrency market surpassed the barrier of \$100 billion market
capitalization in June 2017, after months of steady growth. Despite its
increasing relevance in the financial world, however, a comprehensive analysis
of the whole system is still lacking, as most studies have focused exclusively
on the behaviour of one (Bitcoin) or few cryptocurrencies. Here, we consider
the history of the entire market and analyse the behaviour of 1,469
cryptocurrencies introduced between April 2013 and June 2017. We reveal that,
while new cryptocurrencies appear and disappear continuously and their market
capitalization is increasing (super-)exponentially, several statistical
properties of the market have been stable for years. These include the number
of active cryptocurrencies, the market share distribution and the turnover of
cryptocurrencies. Adopting an ecological perspective, we show that the
so-called neutral model of evolution is able to reproduce a number of key
empirical observations, despite its simplicity and the assumption of no
selective advantage of one cryptocurrency over another. Our results shed light
on the properties of the cryptocurrency market and establish a first formal
link between ecological modelling and the study of this growing system. We
anticipate they will spark further research in this direction.


  This paper studies invasion waves in the diffusive
Competitor-Competitor-Mutualist model generalizing the 2-species Lotka-Volterra
model studied by Weinberger et al. The mutualist may benefit the invading or
the resident species producing two different types of invasions. Sufficient
conditions for linear determinacy are derived in both cases, and when they
hold, explicit formulas for linear spreading speeds of the invasions are
obtained by linearizing the model. While in the first case the linear speed is
increased by the mutualist, it is unaffected in the second case. Mathematical
methods are based on converting the model into a cooperative reaction-diffusion
system.


  Crowded environments modify the diffusion of macromolecules, generally
slowing their movement and inducing transient anomalous subdiffusion. The
presence of obstacles also modifies the kinetics and equilibrium behavior of
tracers. While previous theoretical studies of particle diffusion have
typically assumed either impenetrable obstacles or binding interactions that
immobilize the particle, in many cellular contexts bound particles remain
mobile. Examples include membrane proteins or lipids with some entry and
diffusion within lipid domains and proteins that can enter into membraneless
organelles or compartments such as the nucleolus. Using a lattice model, we
studied the diffusive movement of tracer particles which bind to soft
obstacles, allowing tracers and obstacles to occupy the same lattice site. For
sticky obstacles, bound tracer particles are immobile, while for slippery
obstacles, bound tracers can hop without penalty to adjacent obstacles. In both
models, binding significantly alters tracer motion. The type and degree of
motion while bound is a key determinant of the tracer mobility: slippery
obstacles can allow nearly unhindered diffusion, even at high obstacle filling
fraction. To mimic compartmentalization in a cell, we examined how obstacle
size and a range of bound diffusion coefficients affect tracer dynamics. The
behavior of the model is similar in two and three spatial dimensions. Our work
has implications for protein movement and interactions within cells.


  A large number of protein sequences are becoming available through the
application of novel high-throughput sequencing technologies. Experimental
functional characterization of these proteins is time-consuming and expensive,
and is often only done rigorously for few selected model organisms.
Computational function prediction approaches have been suggested to fill this
gap. The functions of proteins are classified using the Gene Ontology (GO),
which contains over 40,000 classes. Additionally, proteins have multiple
functions, making function prediction a large-scale, multi-class, multi-label
problem.
  We have developed a novel method to predict protein function from sequence.
We use deep learning to learn features from protein sequences as well as a
cross-species protein-protein interaction network. Our approach specifically
outputs information in the structure of the GO and utilizes the dependencies
between GO classes as background information to construct a deep learning
model. We evaluate our method using the standards established by the
Computational Assessment of Function Annotation (CAFA) and demonstrate a
significant improvement over baseline methods such as BLAST, with significant
improvement for predicting cellular locations.


  Genetic alterations initiate tumors and enable the evolution of drug
resistance. The pro-cancer view of mutations is however incomplete, and several
studies show that mutational load can reduce tumor fitness. Given its negative
effect, genetic load should make tumors more sensitive to anticancer drugs.
Here, we test this hypothesis across all major types of cancer from the Cancer
Cell Line Encyclopedia, which provides genetic and expression data of 496 cell
lines together with their response to 24 common anticancer drugs. We found that
the efficacy of 9 out of 24 drugs showed significant association with genetic
load in a pan-cancer analysis. The associations for some tissue-drug
combinations were remarkably strong, with genetic load explaining up to 83% of
the variance in the drug response. Overall, the role of genetic load depended
on both the drug and the tissue type with 10 tissues being particularly
vulnerable to genetic load. We also identified changes in gene expression
associated with increased genetic load, which included cell-cycle checkpoints,
DNA damage and apoptosis. Our results show that genetic load is an important
component of tumor fitness and can predict drug sensitivity. Beyond being a
biomarker, genetic load might be a new, unexplored vulnerability of cancer.


  A mathematico-physically valid formulation is required to infer properties of
disordered protein conformations from single-molecule F\"orster resonance
energy transfer (smFRET). Conformational dimensions inferred by conventional
approaches that presume a homogeneous conformational ensemble can be
unphysical. When all possible---heterogeneous as well as
homogeneous---conformational distributions are taken into account without
prejudgement, a single value of average transfer efficiency $\langle E\rangle$
between dyes at two chain ends is generally consistent with highly diverse,
multiple values of the average radius of gyration $\langle R_{\rm g}\rangle$.
Here we utilize unbiased conformational statistics from a coarse-grained
explicit-chain model to establish a general logical framework to quantify this
fundamental ambiguity in smFRET inference. As an application, we address the
long-standing controversy regarding the denaturant dependence of $\langle
R_{\rm g}\rangle$ of unfolded proteins, focusing on Protein L as an example.
Conventional smFRET inference concluded that $\langle R_{\rm g}\rangle$ of
unfolded Protein L is highly sensitive to [GuHCl], but data from small-angle
X-ray scattering (SAXS) suggested a near-constant $\langle R_{\rm g}\rangle$
irrespective of [GuHCl]. Strikingly, the present analysis indicates that
although the reported $\langle E\rangle$ values for Protein L at [GuHCl] = 1 M
and 7 M are very different at 0.75 and 0.45, respectively, the Bayesian
$R^2_{\rm g}$ distributions consistent with these two $\langle E\rangle$ values
overlap by as much as $75\%$. Our findings suggest, in general, that the
smFRET-SAXS discrepancy regarding unfolded protein dimensions likely arise from
highly heterogeneous conformational ensembles at low or zero denaturant, and
that additional experimental probes are needed to ascertain the nature of this
heterogeneity.


  To cluster sequences given only their read-set representations, one may try
to reconstruct each one from the corresponding read set, and then employ
conventional (dis)similarity measures such as the edit distance on the
assembled sequences. This approach is however problematic and we propose
instead to estimate the similarities directly from the read sets. Our approach
is based on an adaptation of the Monge-Elkan similarity known from the field of
databases. It avoids the NP-hard problem of sequence assembly. For low coverage
data it results in a better approximation of the true sequence similarities and
consequently in better clustering, in comparison to the
first-assemble-then-cluster approach.


  At 49 C erythrocytes undergo morphological changes due to an internal force,
but the origin of the force that drives changes is not clear. Here we point out
that our recent experiments on thermally induced force-release in hemoglobin
can provide an explanation for the morphological changes of erythrocytes.


  Neurotransmitter receptor molecules, concentrated in synaptic membrane
domains along with scaffolds and other kinds of proteins, are crucial for
signal transmission across chemical synapses. In common with other membrane
protein domains, synaptic domains are characterized by low protein copy numbers
and protein crowding, with rapid stochastic turnover of individual molecules.
We study here in detail a stochastic lattice model of the receptor-scaffold
reaction-diffusion dynamics at synaptic domains that was found previously to
capture, at the mean-field level, the self-assembly, stability, and
characteristic size of synaptic domains observed in experiments. We show that
our stochastic lattice model yields quantitative agreement with mean-field
models of nonlinear diffusion in crowded membranes. Through a combination of
analytic and numerical solutions of the master equation governing the reaction
dynamics at synaptic domains, together with kinetic Monte Carlo simulations, we
find substantial discrepancies between mean-field and stochastic models for the
reaction dynamics at synaptic domains. Based on the reaction and diffusion
properties of synaptic receptors and scaffolds suggested by previous
experiments and mean-field calculations, we show that the stochastic
reaction-diffusion dynamics of synaptic receptors and scaffolds provide a
simple physical mechanism for collective fluctuations in synaptic domains, the
molecular turnover observed at synaptic domains, key features of the observed
single-molecule trajectories, and spatial heterogeneity in the effective rates
at which receptors and scaffolds are recycled at the cell membrane. Our work
sheds light on the physical mechanisms and principles linking the collective
properties of membrane protein domains to the stochastic dynamics that rule
their molecular~components.


  Much of contemporary systems biology owes its success to the abstraction of a
network, the idea that diverse kinds of molecular, cellular, and organismal
species and interactions can be modeled as relational nodes and edges in a
graph of dependencies. Since the advent of high-throughput data acquisition
technologies in fields such as genomics, metabolomics, and neuroscience, the
automated inference and reconstruction of such interaction networks directly
from large sets of activation data, commonly known as reverse-engineering, has
become a routine procedure. Whereas early attempts at network
reverse-engineering focused predominantly on producing maps of system
architectures with minimal predictive modeling, reconstructions now play
instrumental roles in answering questions about the statistics and dynamics of
the underlying systems they represent. Many of these predictions have clinical
relevance, suggesting novel paradigms for drug discovery and disease treatment.
While other reviews focus predominantly on the details and effectiveness of
individual network inference algorithms, here we examine the emerging field as
a whole. We first summarize several key application areas in which inferred
networks have made successful predictions. We then outline the two major
classes of reverse-engineering methodologies, emphasizing that the type of
prediction that one aims to make dictates the algorithms one should employ. We
conclude by discussing whether recent breakthroughs justify the computational
costs of large-scale reverse-engineering sufficiently to admit it as a mainstay
in the quantitative analysis of living systems.


  Understanding the genetic basis of phenotypic plasticity is crucial for
predicting and managing climate change effects on wild plants and crops. Here,
we combined crop modeling and quantitative genetics to study the genetic
control of oil yield plasticity for multiple abiotic stresses in sunflower.
  First we developed stress indicators to characterize 14 environments for
three abiotic stresses (cold, drought and nitrogen) using the SUNFLO crop model
and phenotypic variations of three commercial varieties. The computed plant
stress indicators better explain yield variation than descriptors at the
climatic or crop levels. In those environments, we observed oil yield of 317
sunflower hybrids and regressed it with three selected stress indicators. The
slopes of cold stress norm reaction were used as plasticity phenotypes in the
following genome-wide association study.
  Among the 65,534 tested SNP, we identified nine QTL controlling oil yield
plasticity to cold stress. Associated SNP are localized in genes previously
shown to be involved in cold stress responses: oligopeptide transporters, LTP,
cystatin, alternative oxidase, or root development. This novel approach opens
new perspectives to identify genomic regions involved in
genotype-by-environment interaction of a complex traits to multiple stresses in
realistic natural or agronomical conditions.


  This work is devoted to investigating the evolution of concentration in a
genetic regulation system, when the synthesis reaction rate is under additive
and multiplicative asymmetric stable L\'evy fluctuations. By focusing on the
impact of skewness (i.e., non-symmetry) in the probability distributions of
noise, we find that via examining the mean first exit time (MFET) and the first
escape probability (FEP), the asymmetric fluctuations, interacting with
nonlinearity in the system, lead to peculiar likelihood for transcription. This
includes, in the additive noise case, realizing higher likelihood of
transcription for larger positive skewness (i.e., asymmetry) index $\beta$,
causing a stochastic bifurcation at the non-Gaussianity index value $\alpha=1$
(i.e., it is a separating point or line for the likelihood for transcription),
and achieving a turning point at the threshold value $\beta \approx -0.5$
(i.e., beyond which the likelihood for transcription suddenly reversed for
$\alpha$ values). The stochastic bifurcation and turning point phenomena do not
occur in the symmetric noise case ($\beta =0$). While in the multiplicative
noise case, non-Gaussianity index value $\alpha=1$ is a separating point or
line for both the mean first exit time (MFET) and the first escape probability
(FEP). We also investigate the noise enhanced stability phenomenon.
Additionally, we are able to specify the regions in the whole parameter space
for the asymmetric noise, in which we attain desired likelihood for
transcription. We have conducted a series of numerical experiments in
`regulating' the likelihood of gene transcription by tuning asymmetric stable
L\'evy noise indexes. This work offers insights for possible ways of achieving
gene regulation in experimental research.


  Natural phenomena frequently involve a very large number of interacting
molecules moving in confined regions of space. Cellular transport by motor
proteins is an example of such collective behavior. We derive a deterministic
compartmental model for the unidirectional flow of particles along a
one-dimensional lattice of sites with nearest-neighbor interactions between the
particles. The flow between consecutive sites is governed by a soft simple
exclusion principle and by attracting or repelling forces between neighboring
particles. Using tools from contraction theory, we prove that the model admits
a unique steady-state and that every trajectory converges to this steady-state.
Analysis and simulations of the effect of the attracting and repelling forces
on this steady-state highlight the crucial role that these forces may play in
increasing the steady-state flow, and reveal that this increase stems from the
alleviation of traffic jams along the lattice. Our theoretical analysis
clarifies microscopic aspects of complex multi-particle dynamic processes.


  The edit distance under the DCJ model can be computed in linear time for
genomes with equal content or with Indels. But it becomes NP-Hard in the
presence of duplications, a problem largely unsolved especially when Indels are
considered. In this paper, we compare two mainstream methods to deal with
duplications and associate them with Indels: one by deletion, namely
DCJ-Indel-Exemplar distance; versus the other by gene matching, namely
DCJ-Indel-Matching distance. We design branch-and-bound algorithms with set of
optimization methods to compute exact distances for both. Furthermore, median
problems are discussed in alignment with both of these distance methods, which
are to find a median genome that minimizes distances between itself and three
given genomes. Lin-Kernighan (LK) heuristic is leveraged and powered up by
sub-graph decomposition and search space reduction technologies to handle
median computation. A wide range of experiments are conducted on synthetic data
sets and real data sets to show pros and cons of these two distance metrics per
se, as well as putting them in the median computation scenario.


  Telomeres, the nucleoprotein complexes at the termini of linear chromosomes,
are essential for the processes of end replication, end-protection, and
chromatin segregation. The Mre11 complex is involved in multiple cellular roles
in DNA repair and structure in the regulation and function of telomere size
homeostasis. In this study, we characterize yeast telomere chromatin structure,
phenotypic heritability, and chromatin segregation in both wild-type [MRE11]
and A470 motif alleles. MRE11 strains confer a telomere size of 300 base pairs
of G+T irregular simple sequence repeats. This DNA and a portion of
subtelomeric DNA is embedded in a telosome: an MNase-resistant non-nucleosomal
particle. Chromatin immunoprecipitation shows a three to four-fold lower
occupancy of Mre11A470T proteins than wild-type proteins in telosomes.
Telosomes containing the Mre11A470T protein confer a greater resistance to
MNase digestion than wild-type telosomes. The integration of a wild-type MRE11
allele into an ectopic locus in the genome of a mre11A470T mutant and the
introduction of a mre11A470T allele at an ectopic site in a wild-type strain
lead to unexpectedly differing results. In each case, the replicated sister
chromatids inherit telosomes containing only the protein encoded by the genomic
mre11 locus, even in the presence of protein encoded by the opposing ectopic
allele. We hypothesize that the telosome segregates by a conservative
mechanism. These data support a mechanism for the linkage between sister
chromatid replication and maintenance of either identical mutant or identical
wild-type telosomes after replication of sister chromatids. These data suggest
the presence of an active mechanism for chromatin segregation in yeast.


  Stochastic reaction networks that exhibit bistability are common in many
fields such as systems biology and materials science. Sampling of the
stationary distribution is crucial for understanding and characterizing the
long term dynamics of bistable stochastic dynamical systems. However, this is
normally hindered by the insufficient sampling of the rare transitions between
the two metastable regions. In this paper, we apply the parallel replica
(ParRep) method for continuous time Markov chain to accelerate the stationary
distribution sampling of bistable stochastic reaction networks. The proposed
method uses parallel computing to accelerate the sampling of rare transitions
and it is very easy to implement. We combine ParRep with the path space
information bounds for parametric sensitivity analysis. We demonstrate the
efficiency and accuracy of the method by studying the Schl\"{o}gl model and the
genetic switches network.


  Proteins perform their functions usually by interacting with other proteins.
Predicting which proteins interact is a fundamental problem. Experimental
methods are slow, expensive, and have a high rate of error. Many computational
methods have been proposed among which sequence-based ones are very promising.
However, so far no such method is able to predict effectively the entire human
interactome: they require too much time or memory. We present SPRINT (Scoring
PRotein INTeractions), a new sequence-based algorithm and tool for predicting
protein-protein interactions. We comprehensively compare SPRINT with
state-of-the-art programs on seven most reliable human PPI datasets and show
that it is more accurate while running orders of magnitude faster and using
very little memory. SPRINT is the only program that can predict the entire
human interactome. Our goal is to transform the very challenging problem of
predicting the entire human interactome into a routine task. The source code of
SPRINT is freely available from github.com/lucian-ilie/SPRINT/ and the datasets
and predicted PPIs from www.csd.uwo.ca/faculty/ilie/SPRINT/.


  Distances between sequences based on their $k$-mer frequency counts can be
used to reconstruct phylogenies without first computing a sequence alignment.
Past work has shown that effective use of k-mer methods depends on 1)
model-based corrections to distances based on $k$-mers and 2) breaking long
sequences into blocks to obtain repeated trials from the sequence-generating
process. Good performance of such methods is based on having many high-quality
blocks with many homologous sites, which can be problematic to guarantee a
priori.
  Nature provides natural blocks of sequences into homologous regions---namely,
the genes. However, directly using past work in this setting is problematic
because of possible discordance between different gene trees and the underlying
species tree. Using the multispecies coalescent model as a basis, we derive
model-based moment formulas that involve the divergence times and the
coalescent parameters. From this setting, we prove identifiability results for
the tree and branch length parameters under the Jukes-Cantor model of sequence
mutations.


  Extensive cooperation among unrelated individuals is unique to humans, who
often sacrifice personal benefits for the common good and work together to
achieve what they are unable to execute alone. The evolutionary success of our
species is indeed due, to a large degree, to our unparalleled other-regarding
abilities. Yet, a comprehensive understanding of human cooperation remains a
formidable challenge. Recent research in social science indicates that it is
important to focus on the collective behavior that emerges as the result of the
interactions among individuals, groups, and even societies. Non-equilibrium
statistical physics, in particular Monte Carlo methods and the theory of
collective behavior of interacting particles near phase transition points, has
proven to be very valuable for understanding counterintuitive evolutionary
outcomes. By studying models of human cooperation as classical spin models, a
physicist can draw on familiar settings from statistical physics. However,
unlike pairwise interactions among particles that typically govern solid-state
physics systems, interactions among humans often involve group interactions,
and they also involve a larger number of possible states even for the most
simplified description of reality. The complexity of solutions therefore often
surpasses that observed in physical systems. Here we review experimental and
theoretical research that advances our understanding of human cooperation,
focusing on spatial pattern formation, on the spatiotemporal dynamics of
observed solutions, and on self-organization that may either promote or hinder
socially favorable states.


  Realistic description of patient-specific mechanical properties of clotting
dynamics presents a major challenge. Available patient-specific data falls
short of robustly characterizing myriads of complex dynamic interactions that
happen during clotting. We propose a simplified switching model for a key part
of the coagulation cascade that describes dynamics of just four variables. The
model correctly predicts prolonged activity of thrombin, an important enzyme in
the clotting process, in certain plasma factor compositions. The activity
sustains beyond the time which is conventionally considered to be the end of
clotting. This observation along with the simplified model is hypothesized as a
necessary step towards effectively studying patient-specific properties of
clotting dynamics in realistic geometries.


  We present a proof of principle for the phenomenon of the tragedy of the
commons that is at the center of many theories on the evolution of cooperation.
We establish the tragedy in the context of a general chemostat model with two
species, the cooperator and the cheater. Both species have the same growth rate
function and yield constant, but the cooperator allocates a portion of the
nutrient uptake towards the production of a public good -the "Commons" in the
Tragedy- which is needed to digest the externally supplied nutrient. The
cheater on the other hand does not produce this enzyme, and allocates all
nutrient uptake towards its own growth. We prove that when the cheater is
present initially, both the cooperator and the cheater will eventually go
extinct, hereby confirming the occurrence of the tragedy. We also show that
without the cheater, the cooperator can survive indefinitely, provided that at
least a low level of public good or processed nutrient is available initially.
Our results provide a predictive framework for the analysis of
cooperator-cheater dynamics in a powerful model system of experimental
evolution.


  We study here global and local entanglements of open protein chains by
implementing the concept of knotoids. Knotoids have been introduced in 2012 by
Vladimir Turaev as a generalization of knots in 3-dimensional space. More
precisely, knotoids are diagrams representing projections of open curves in 3D
space, in contrast to knot diagrams which represent projections of closed
curves in 3D space. The intrinsic difference with classical knot theory is that
the generalization provided by knotoids admits non-trivial topological
entanglement of the open curves provided that their geometry is frozen as it is
the case for crystallized proteins. Consequently, our approach doesn't require
the closure of chains into loops which implies that the geometry of analysed
chains does not need to be changed by closure in order to characterize their
topology. Our study revealed that the knotoid approach detects protein regions
that were classified earlier as knotted and also new, topologically interesting
regions that we classify as pre-knotted.


  Essential genes constitute the core of genes which cannot be mutated too much
nor lost along the evolutionary history of a species. Natural selection is
expected to be stricter on essential genes and on conserved (highly shared)
genes, than on genes that are either nonessential or peculiar to a single or a
few species. In order to further assess this expectation, we study here how
essentiality of a gene is connected with its degree of conservation among
several unrelated bacterial species, each one characterised by its own codon
usage bias. Confirming previous results on E. coli, we show the existence of a
universal exponential relation between gene essentiality and conservation in
bacteria. Moreover we show that, within each bacterial genome, there are at
least two groups of functionally distinct genes, characterised by different
levels of conservation and codon bias: i) a core of essential genes, mainly
related to cellular information processing; ii) a set of less conserved
nonessential genes with prevalent functions related to metabolism. In
particular, the genes in the first group are more retained among species, are
subject to a stronger purifying conservative selection and display a more
limited repertoire of synonymous codons. The core of essential genes is close
to the minimal bacterial genome, which is in the focus of recent studies in
synthetic biology, though we confirm that orthologs of genes that are essential
in one species are not necessarily essential in other species. We also list a
set of highly shared genes which, reasonably, could constitute a reservoir of
targets for new anti-microbial drugs.


  One of the first successes of neutral ecology was to predict
realistically-broad distributions of rare and abundant species. However, it has
remained an outstanding theoretical challenge to describe how this distribution
of abundances changes with spatial scale, and this gap has hampered attempts to
use observed species abundances as a way to quantify what non-neutral processes
are needed to fully explain observed patterns. To address this, we introduce a
new formulation of spatial neutral biodiversity theory and derive analytical
predictions for the way abundance distributions change with scale. For tropical
forest data where neutrality has been extensively tested before now, we apply
this approach and identify an incompatibility between neutral fits at regional
and local scales. We use this approach derive a sharp quantification of what
remains to be explained by non-neutral processes at the local scale, setting a
quantitative target for more general models for the maintenance of
biodiversity.


  Sex hormone-binding globulin (SHBG) is a binding protein that regulates
availability of steroids hormones in the plasma. Although best known as steroid
carrier, studies have associated SHBG in modulating behavioral aspects related
to sexual receptivity. Among steroids, estradiol (17\b{eta}-estradiol,
oestradiol or E2) is well recognized as the most active endogenous female
hormone, exerting important roles in reproductive and nonreproductive
functions. Thus, in this study we aimed to employ molecular dynamics (MD) and
docking techniques for quantifying the interaction energy between a complex
aqueous solution, composed by different salts, SHBG and E2. Due to glucose
concentration resembles those observed in diabetic levels, special emphasis was
devoted to uncover the main consequences of this carbohydrate on the SHBG and
E2 molecules. We also examined possible energetic changes due to solution on
the binding energy of SHBG-E2 complex. In this framework, our calculations
uncovered a remarkable interaction energy between glucose and SHBG surface.
Surprisingly, we also observed solute components movement toward SHBG yielding
clusters surrounding the protein. This finding, corroborated by the higher
energy and shorter distance found between glucose and SHBG, suggests a scenario
in favor of a detainment state. In addition, in spite of protein superficial
area increment it does not exerted modification on binding site area nor over
binding energy SHBG-E2 complex. Finally, our calculations also highlighted an
interaction between E2 and glucose when the hormone was immersed in the
solution. In summary, our findings contribute for a better comprehension of
both SHBG and E2 interplay with aqueous solution components.


  The origins of life stands among the great open scientific questions of our
time. While a number of proposals exist for possible starting points in the
pathway from non-living to living matter, these have so far not achieved states
of complexity that are anywhere near that of even the simplest living systems.
A key challenge is identifying the properties of living matter that might
distinguish living and non-living physical systems such that we might build new
life in the lab. This review is geared towards covering major viewpoints on the
origin of life for those new to the origin of life field, with a forward look
towards considering what it might take for a physical theory that universally
explains the phenomenon of life to arise from the seemingly disconnected array
of ideas proposed thus far. The hope is that a theory akin to our other
theories in fundamental physics might one day emerge to explain the phenomenon
of life, and in turn finally permit solving its origins.


  We will give a simple, unified, possible explanation of several debated
genetic issues on today's humans, Neandertals and Denisovans. In particular it
is shown by means of a simple mathematical model why there is little genetic
variation in todays's human population or in Western Neandertal population, why
all mtDNA and y-chromosomes in today's humans seem to have African origin with
no trace of Neandertal nor Denosovan mtDNA or y-chromosomes, why a big part of
the European gene pool is young (from Neolitic time), and why today's East
Asians have mode Neandertal genes than today's Europeans.


  High throughput immune repertoire sequencing is promising to lead to new
statistical diagnostic tools for medicine and biology. Successful
implementations of these methods require a correct characterization, analysis
and interpretation of these datasets. We present IGoR -- a new comprehensive
tool that takes B or T-cell receptors sequence reads and quantitatively
characterizes the statistics of receptor generation from both cDNA and gDNA. It
probabilistically annotates sequences and its modular structure can investigate
models of increasing biological complexity for different organisms. For B-cells
IGoR returns the hypermutation statistics, which we use to reveal
co-localization of hypermutations along the sequence. We demonstrate that IGoR
outperforms existing tools in accuracy and estimate the sample sizes needed for
reliable repertoire characterization.


  Controlling the self-assembly of supramolecular structures is vital for
living cells, and a central challenge for engineering at the nano- and
microscales. Nevertheless, even particles without optimized shapes can robustly
form well-defined morphologies. This is the case in numerous medical conditions
where normally soluble proteins aggregate into fibers. Beyond the diversity of
molecular mechanisms involved, we propose that fibers generically arise from
the aggregation of irregular particles with short-range interactions. Using a
minimal model of ill-fitting, sticky particles, we demonstrate robust fiber
formation for a variety of particle shapes and aggregation conditions.
Geometrical frustration plays a crucial role in this process, and accounts for
the range of parameters in which fibers form as well as for their metastable
character.


  Current understanding of how contractility emerges in disordered actomyosin
networks of non-muscle cells is still largely based on the intuition derived
from earlier works on muscle contractility. This view, however, largely
overlooks the free energy gain following passive cross-linker binding, which,
even in the absence of active fluctuations, provides a thermodynamic drive
towards highly overlapping filamentous states. In this work, we shed light on
this phenomenon, showing that passive cross-linkers, when considered in the
context of two anti-parallel filaments, generate noticeable contractile forces.
However, as binding free energy of cross-linkers is increased, a sharp onset of
kinetic arrest follows, greatly diminishing effectiveness of this contractility
mechanism, allowing the network to contract only with weakly resisting tensions
at its boundary. We have carried out stochastic simulations elucidating this
mechanism, followed by a mean-field treatment that predicts how contractile
forces asymptotically scale at small and large binding energies, respectively.
Furthermore, when considering an active contractile filament pair, based on
non-muscle myosin II, we found that the non-processive nature of these motors
leads to highly inefficient force generation, due to recoil slippage of the
overlap during periods when the motor is dissociated. However, we discovered
that passive cross-linkers can serve as a structural ratchet during these
unbound motor time spans, resulting in vast force amplification. Our results
shed light on the non-equilibrium effects of transiently binding proteins in
biological active matter, as observed in the non-muscle actin cytoskeleton,
showing that highly efficient contractile force dipoles result from synergy of
passive cross-linker and active motor dynamics, via a ratcheting mechanism on a
funneled energy landscape.


  Biochemical networks play a crucial role in biological systems, implementing
a broad range of vital functions. They normally operate at low copy numbers and
in spatial settings, but this is often ignored and well-stirred conditions are
assumed. Yet, it is increasingly becoming clear that even microscopic spatial
inhomogeneities oftentimes can induce significant differences on the
macroscopic level. Since experimental observation of single-molecule behavior
is extremely challenging, theoretical modeling of biochemical reactions on the
single-particle level is an important tool for understanding spatial effects in
biochemical systems. While purely analytical models quickly become intractable
here, spatial-stochastic simulations can capture a wide range of biochemical
processes with the necessary levels of detail. Here we discuss different
techniques for spatial-stochastic simulation of reaction-diffusion systems, and
explain important precautions required to make them biochemically accurate and
efficient. We illustrate non-negligible accuracy issues arising even in the
most simple approaches to biochemical simulation, and present methods to deal
with them. We first explain how Brownian Dynamics, a widely used particle-based
diffusion simulation technique with fixed propagation time, can be adapted to
incorporate chemical reactions, and portray a range of schemes that elaborate
on this idea. We then introduce event-driven spatial-stochastic simulation
methods, in which system updates are performed asynchronously with
situation-dependent, varying time steps; here we particularly focus on eGFRD, a
computationally efficient particle-based algorithm that makes use of analytical
functions to accurately sample interparticle reactions and diffusive motion
with large jumps in time and space. We end by briefly presenting recent
developments in the field of spatial-stochastic biochemical simulation.


  The analysis of benthic assemblages is a valuable tool to describe the
ecological status of transitional water ecosystems, but species are extremely
sensitive and respond to both microhabitat and seasonal differences. The
identification of changes in the composition of the macrobenthic community in
specific microhabitats can then be used as an "early warning" for environmental
changes which may affect the economic and ecological importance of lagoons,
through their provision of Ecosystem Services. From a conservational point of
view, the appropriate definition of the spatial aggregation level of
microhabitats or local communities is of crucial importance. The main objective
of this work is to assess the role of the spatial scale in the analysis of
lagoon biodiversity. First, we analyze the variation in the sample coverage for
alternative aggregations of the monitoring stations in three lagoons of the Po
River Delta. Then, we analyze the variation of a class of entropy indices by
mixed effects models, properly accounting for the fixed effects of biotic and
abiotic factors and random effects ruled by nested sources of variability
corresponding to alternative definitions of local communities. Finally, we
address biodiversity partitioning by a generalized diversity measure, namely
the Tsallis entropy, and for alternative definitions of the local communities.
The main results obtained by the proposed statistical protocol are presented,
discussed and framed in the ecological context.


  To understand how neurons and nervous systems first evolved, we need an
account of the origins of neural elongations: Why did neural elongations (axons
and dendrites) first originate, such that they could become the central
component of both neurons and nervous systems? Two contrasting conceptual
accounts provide different answers to this question. Braitenberg's vehicles
provide the iconic illustration of the dominant input-output (IO) view. Here
the basic role of neural elongations is to connect sensors to effectors, both
situated at different positions within the body. For this function, neural
elongations are thought of as comparatively long and specific connections,
which require an articulated body involving substantial developmental processes
to build. Internal coordination (IC) models stress a different function for
early nervous systems. Here the coordination of activity across extended parts
of a multicellular body is held central, in particular for the contractions of
(muscle) tissue. An IC perspective allows the hypothesis that the earliest
proto-neural elongations could have been functional even when they were
initially simple short and random connections, as long as they enhanced the
patterning of contractile activity across a multicellular surface. The present
computational study provides a proof of concept that such short and random
neural elongations can play this role. While an excitable epithelium can
generate basic forms of patterning for small body-configurations, adding
elongations allows such patterning to scale up to larger bodies. This result
supports a new, more gradual evolutionary route towards the origins of the very
first full neurons and nervous systems.


  Phylogenetic models have polynomial parametrization maps. For symmetric
group-based models, Matsen studied the polynomial inequalities that
characterize the joint probabilities in the image of these parametrizations. We
employ this description for maximum likelihood estimation via numerical
algebraic geometry. In particular, we explore an example where the maximum
likelihood estimate does not exist, which would be difficult to discover
without using algebraic methods. We also study the embedding problem for
symmetric group-based models, i.e. we identify which mutation matrices are
matrix exponentials of rate matrices that are invariant under a group action.


  We study a deterministic mechanistic model for the flow of ribosomes along
the mRNA molecule, called the ribosome flow model with extended objects
(RFMEO). This model encapsulates many realistic features of translation
including non-homogeneous transition rates along the mRNA, the fact that every
ribosome covers several codons, and the fact that ribosomes cannot overtake one
another.
  The RFMEO is a mean-field approximation of an important model from
statistical mechanics called the totally asymmetric simple exclusion process
with extended objects (TASEPEO). We demonstrate that the RFMEO describes
biophysical aspects of translation better than previous mean-field
approximations, and that its predictions correlate well with those of TASEPEO.
However, unlike TASEPEO, the RFMEO is amenable to rigorous analysis using tools
from systems and control theory. We show that the ribosome density profile
along the mRNA in the RFMEO converges to a unique steady-state density that
depends on the length of the mRNA, the transition rates along it, and the
number of codons covered by every ribosome, but not on the initial density of
ribosomes along the mRNA. In particular, the protein production rate also
converges to a unique steady-state. Furthermore, if the transition rates along
the mRNA are periodic with a common period T then the ribosome density along
the mRNA and the protein production rate converge to a unique periodic pattern
with period T, that is, the model entrains to periodic excitations in the
transition rates.
  We believe that the RFMEO could be useful for modeling, understanding, and
re-engineering translation as well as other important biological processes.


  The ribosome flow model on a ring (RFMR) is a deterministic model for
translation of a circularized mRNA. We derive a new spectral representation for
the optimal steady-state production rate and the corresponding optimal
steady-state ribosomal density in the RFMR. This representation has several
important advantages. First, it provides a simple and numerically stable
algorithm for determining the optimal values even in very long rings. Second,
it enables efficient computation of the sensitivity of the optimal production
rate to small changes in the transition rates along the mRNA. Third, it implies
that the optimal steady-state production rate is a strictly concave function of
the transition rates. Thus maximizing the optimal steady-state production rate
with respect to the rates, under an affine constraint on the rates becomes a
convex optimization problem that admits a unique solution, which can be
determined numerically using highly efficient algorithms. This optimization
problem is important, for example, when re-engineering heterologous genes in a
host organism. We describe the implications of our results to this and other
aspects of translation.


  Synthetic biology is a growing interdisciplinary field, with far-reaching
applications, which aims to design biochemical systems that behave in a desired
manner. With the advancement of strand-displacement DNA computing, a large
class of abstract biochemical networks may be physically realized using DNA
molecules. Methods for systematic design of the abstract systems with
prescribed behaviors have been predominantly developed at the (less-detailed)
deterministic level. However, stochastic effects, neglected at the
deterministic level, are increasingly found to play an important role in
biochemistry. In such circumstances, methods for controlling the intrinsic
noise in the system are necessary for a successful network design at the
(more-detailed) stochastic level. To bridge the gap, the noise-control
algorithm for designing biochemical networks is developed in this paper. The
algorithm structurally modifies any given reaction network under mass-action
kinetics, in such a way that (i) controllable state-dependent noise is
introduced into the stochastic dynamics, while (ii) the deterministic dynamics
are preserved. The capabilities of the algorithm are demonstrated on a
production-decay reaction system, and on an exotic system displaying
bistability. For the production-decay system, it is shown that the algorithm
may be used to redesign the network to achieve noise-induced multistability.
For the exotic system, the algorithm is used to redesign the network to control
the stochastic switching, and achieve noise-induced oscillations.


  In this paper we propose a method to study a general vector-hosts
mathematical model in order to explain how the changes in biodiversity could
influence the dynamics of vector-borne diseases. We find that under the
assumption of frequency-dependent transmission, i.e. the assumption that the
number of contacts are diluted by the total population of hosts, the presence
of a competent host is a necessary condition for the existence of an endemic
state. In addition, we obtain that in the case of an endemic disease with a
unique competent and resilient host, an increase in its density amplifies the
disease.


  We consider inference about the history of a sample of DNA sequences,
conditional upon the haplotype counts and the number of segregating sites
observed at the present time. After deriving some theoretical results in the
coalescent setting, we implement rejection sampling and importance sampling
schemes to perform the inference. The importance sampling scheme addresses an
extension of the Ewens Sampling Formula for a configuration of haplotypes and
the number of segregating sites in the sample. The implementations include both
constant and variable population size models. The methods are illustrated by
two human Y chromosome data sets.


  Chiral heteropolymers such as larger globular proteins can simultaneously
support multiple length scales. The interplay between different scales brings
about conformational diversity, and governs the structure of the energy
landscape. Multiple scales produces also complex dynamics, which in the case of
proteins sustains live matter. However, thus far no clear understanding exist,
how to distinguish the various scales that determine the structure and dynamics
of a complex protein. Here we propose a systematic method to identify the
scales in chiral heteropolymers such as a protein. For this we introduce a
novel order parameter, that not only reveals the scales but also probes the
phase structure. In particular, we argue that a chiral heteropolymer can
simultaneously display traits of several different phases, contingent on the
length scale at which it is scrutinized. Our approach builds on a variant of
Kadanoff's block-spin transformation that we employ to coarse grain piecewise
linear chains such as the C$\alpha$ backbone of a protein. We derive
analytically and then verify numerically a number of properties that the order
parameter can display. We demonstrate how, in the case of crystallographic
protein structures in Protein Data Bank, the order parameter reveals the
presence of different length scales, and we propose that a relation must exist
between the scales, phases, and the complexity of folding pathways.


  Molecular and genetic definitions of human-specific changes to genomic
regulatory networks (GRNs) contributing to development of unique to human
phenotypes remain a highly significant challenge. Genome-wide proximity
placement analysis of diverse families of human-specific genomic regulatory
loci (HSGRL) identified topologically-associating domains (TADs) that are
significantly enriched for HSGRL and designated rapidly-evolving in humans TADs
(Genome Biol Evol. 2016 8; 2774-88). Here, the analysis of HSGRL, hESC-enriched
enhancers, super-enhancers (SEs), and specific sub-TAD structures termed
super-enhancer domains (SEDs) has been performed. Markedly distinct features of
the principal regulatory structures of interphase chromatin evolved in the hESC
genome compared to mouse: the SED quantity is 3-fold higher and the median SED
size is significantly larger. Concomitantly, the overall TAD quantity is
increased by 42% while the median TAD size is significantly decreased (p =
9.11E-37) in the hESC genome. Present analyses illustrate a putative global
role for HSGRL in shaping the human-specific features of the interphase
chromatin organization and functions, which are facilitated by accelerated
creation of new enhancers associated with targeted placement of HSGRL at
defined genomic coordinates. A trend toward the convergence of TAD and SED
architectures of interphase chromatin in the hESC genome may reflect changes of
3D-folding patterns of linear chromatin fibers designed to enhance both
regulatory complexity and functional precision of GRNs by creating
predominantly a single gene per regulatory domain structures.