  Reputation plays a major role in human societies, and it has been proposed as
an explanation for the evolution of cooperation. While the majority of previous
studies equates reputation with a transparent and complete history of players'
past decisions, in real life, reputations are often ambiguous and opaque. Using
web-based experiments, we explore the extent to which opaque reputation works
in isolating defectors, with and without partner selection opportunities. Our
results show that low reputation works as a signal of untrustworthiness,
whereas medium or high reputation are not taken into account by participants
for orienting their choices. We also find that reputation without partner
selection does not promote cooperative behavior; that is, defectors do not turn
into cooperators only for the sake of getting a positive reputation. Finally,
in a third study, we find that, when reputation is pivotal to selection, then a
substantial proportion of would-be-defectors turn into cooperators. Taken
together, these results provide insights on the characteristics of reputation
and on the way in which humans make use of it when selecting partners but also
when knowing that they will be selected.


  After birth setting up an effective immune system is a major challenge for
all living organisms. In this paper we show that this process can be explored
by using the age-specific infant death rate as a kind of sensor. This is made
possible because, as shown by the authors in Berrut et al. (2016), between
birth and a critical age t_c, for all mammals the death rate decreases with age
as an hyperbolic function. For humans t_c is equal to 10 years. At some ages
the hyperbolic fall displays spikes which, it is assumed, correspond to
specific events in the organism's response to exogenous factors. One of these
spikes occurs 10 days after birth and there is another at the age of about 300
days. It is shown that the first spike is related to viral infections whereas
the second is related to bacterial diseases. By going back to former time
periods during which infant mortality was much higher than currently, it is
possible to get a magnified view of these peaks which in turn may give us
useful information about how an organism adapts to new conditions. Apart from
pathogens, the same methodology can be used to study the response to changes in
other external conditions, e.g. supply of food, temperature or oxygen level.


  MicroRNAs (miRNAs) are non-coding RNAs with approximately 22 nucleotides (nt)
that are derived from precursor molecules. These precursor molecules or
pre-miRNAs often fold into stem-loop hairpin structures. However, a large
number of sequences with pre-miRNA-like hairpins can be found in genomes. It is
a challenge to distinguish the real pre-miRNAs from other hairpin sequences
with similar stem-loops (referred to as pseudo pre-miRNAs). Several
computational methods have been developed to tackle this challenge. In this
paper we propose a new method, called MirID, for identifying and classifying
microRNA precursors. We collect 74 features from the sequences and secondary
structures of pre-miRNAs; some of these features are taken from our previous
studies on non-coding RNA prediction while others were suggested in the
literature. We develop a combinatorial feature mining algorithm to identify
suitable feature sets. These feature sets are then used to train support vector
machines to obtain classification models, based on which classifier ensemble is
constructed. Finally we use an AdaBoost algorithm to further enhance the
accuracy of the classifier ensemble. Experimental results on a variety of
species demonstrate the good performance of the proposed method, and its
superiority over existing tools.


  The advent of modern genome sequencing techniques allows for a more stringent
test of the neutrality hypothesis of Darwinian evolution, where all individuals
have the same fitness. Using the individual based model of Wright and Fisher,
we compute the amplitude of neutral aggregation in the genome space, i.e., the
probability of finding two individuals at genetic (hamming) distance k as a
function of genome size L, population size N and mutation probability per base
\nu. In well mixed populations, we show that for N\nu\textless{}1/L, neutral
aggregation is the dominant force and most individuals are found at short
genetic distances from each other. For N\nu\textgreater{}1 on the contrary,
individuals are randomly dispersed in genome space. The results are extended to
geographically dispersed population, where the controlling parameter is shown
to be a combination of mutation and migration probability. The theory we
develop can be used to test the neutrality hypothesis in various ecological and
evolutionary systems.


  Nucleation processes are at the heart of a large number of phenomena, from
cloud formation to protein crystallization. A recently emerging area where
nucleation is highly relevant is the initiation of filamentous protein
self-assembly, a process that has broad implications from medicine to
nanotechnology. As such, spontaneous nucleation of protein fibrils has received
much attention in recent years with many theoretical and experimental studies
focussing on the underlying physical principles. In this paper we make a step
forward in this direction and explore the early time behaviour of filamentous
protein growth in the context of nucleation theory. We first provide an
overview of the thermodynamics and kinetics of spontaneous nucleation of
protein filaments in the presence of one relevant degree of freedom, namely the
cluster size. In this case, we review how key kinetic observables, such as the
reaction order of spontaneous nucleation, are directly related to the physical
size of the critical nucleus. We then focus on the increasingly prominent case
of filament nucleation that includes a conformational conversion of the
nucleating building-block as an additional slow step in the nucleation process.
Using computer simulations, we study the concentration dependence of the
nucleation rate. We find that, under these circumstances, the reaction order of
spontaneous nucleation with respect to the free monomer does no longer relate
to the overall physical size of the nucleating aggregate but rather to the
portion of the aggregate that actively participates in the conformational
conversion. Our results thus provide a novel interpretation of the common
kinetic descriptors of protein filament formation, including the reaction order
of the nucleation step or the scaling exponent of lag times, and put into
perspective current theoretical descriptions of protein aggregation.


  Genetic oscillators play important roles in cell life regulation. The
regulatory efficiency usually depends strongly on the emergence of stable
collective dynamic modes, which requires designing the interactions between
genetic networks. We investigate the dynamics of two identical synthetic
genetic repressilators coupled by an additional plasmid which implements quorum
sensing (QS) in each network thereby supporting global coupling. In a basic
genetic ring oscillator network in which three genes inhibit each other in
unidirectional manner, QS stimulates the transcriptional activity of chosen
genes providing for competition between inhibitory and stimulatory activities
localized in those genes. The "promoter strength", the Hill cooperativity
coefficient of transcription repression, and the coupling strength, i.e.,
parameters controlling the basic rates of genetic reactions, were chosen for
extensive bifurcation analysis. The results are presented as a map of dynamic
regimes. We found that the remarkable multistability of the anti-phase limit
cycle and stable homogeneous and inhomogeneous steady states exists over broad
ranges of control parameters. We studied the anti-phase limit cycle stability
and the evolution of irregular oscillatory regimes in the parameter areas where
the anti-phase cycle loses stability. In these regions we observed developing
complex oscillations, collective chaos, and multistability between regular
limit cycles and complex oscillations over uncommonly large intervals of
coupling strength. QS-coupling stimulates the appearance of intrachaotic
periodic windows with spatially symmetric and asymmetric partial limit cycles
which, in turn, change the type of chaos from a simple anti-phase character
into chaos composed of pieces of the trajectories having alternating polarity.
Abstract truncated...


  Epidemic spreading has been studied for a long time and most of them are
focused on the growing aspect of a single epidemic outbreak. Recently, we
extended the study to the case of recurrent epidemics (Sci. Rep. {\bf 5}, 16010
(2015)) but limited only to a single network. We here report from the real data
of coupled regions or cities that the recurrent epidemics in two coupled
networks are closely related to each other and can show either synchronized
outbreak phase where outbreaks occur simultaneously in both networks or mixed
outbreak phase where outbreaks occur in one network but do not in another one.
To reveal the underlying mechanism, we present a two-layered network model of
coupled recurrent epidemics to reproduce the synchronized and mixed outbreak
phases. We show that the synchronized outbreak phase is preferred to be
triggered in two coupled networks with the same average degree while the mixed
outbreak phase is preferred for the case with different average degrees.
Further, we show that the coupling between the two layers is preferred to
suppress the mixed outbreak phase but enhance the synchronized outbreak phase.
A theoretical analysis based on microscopic Markov-chain approach is presented
to explain the numerical results. This finding opens a new window for studying
the recurrent epidemics in multi-layered networks.


  We present a new computational approach, Action-CSA, to sample multiple
reaction pathways with fixed initial and final states through global
optimization of the Onsager-Machlup action using the conformational space
annealing method. This approach successfully samples not only the most dominant
pathway but also many other possible paths without initial guesses on reaction
pathways. Pathway space is efficiently sampled by crossover operations of a set
of paths and preserving the diversity of sampled pathways. The sampling ability
of the approach is assessed by finding pathways for the conformational changes
of alanine dipeptide and hexane. The benchmarks demonstrate that the rank order
and the transition time distribution of multiple pathways identified by the new
approach are in good agreement with those of long molecular dynamics
simulations. We also show that the lowest action folding pathway of the
mini-protein FSD-1 identified by the new approach is consistent with previous
molecular dynamics simulations and experiments.


  Regulatory networks consist of interacting molecules with a high degree of
mutual chemical specificity. How can these molecules evolve when their function
depends on maintenance of interactions with cognate partners and simultaneous
avoidance of deleterious "crosstalk" with non-cognate molecules? Although
physical models of molecular interactions provide a framework in which
co-evolution of network components can be analyzed, most theoretical studies
have focused on the evolution of individual alleles, neglecting the network. In
contrast, we study the elementary step in the evolution of gene regulatory
networks: duplication of a transcription factor followed by selection for TFs
to specialize their inputs as well as the regulation of their downstream genes.
We show how to coarse grain the complete, biophysically realistic
genotype-phenotype map for this process into macroscopic functional outcomes
and quantify the probability of attaining each. We determine which evolutionary
and biophysical parameters bias evolutionary trajectories towards fast
emergence of new functions and show that this can be greatly facilitated by the
availability of "promiscuity-promoting" mutations that affect TF specificity.


  In phylogenetic analysis, for non-molecular data, particularly morphology,
parsimony optimization is the most commonly employed approach. In the past and
present application of the parsimony principle, extra step numbers have been
added across different characters without proper justification. This practice,
however, has caused the impacts of characters to be inflated or deflated
without a valid reason. To resolve this methodological deficiency, here I
present a corrected parsimony criterion for reconstructing phylogenies. In
essence, character rather than step is the most fundamental unit. Accordingly,
the most parsimonious tree should maximize the sum or average of the
phylogenetic signals, quantified by retention index, contributed by each
character. Assigning proper weights to characters is a separate task that
requires information other than the intra-character step number changing range.


  For single-cell or metagenomic sequencing projects, it is necessary to
sequence with a very high mean coverage in order to make sure that all parts of
the sample DNA get covered by the reads produced. This leads to huge datasets
with lots of redundant data. A filtering of this data prior to assembly is
advisable. Titus Brown et al. (2012) presented the algorithm Diginorm for this
purpose, which filters reads based on the abundance of their $k$-mers. We
present Bignorm, a faster and quality-conscious read filtering algorithm. An
important new feature is the use of phred quality scores together with a
detailed analysis of the $k$-mer counts to decide which reads to keep. With
recommended parameters, in terms of median we remove 97.15% of the reads while
keeping the mean phred score of the filtered dataset high. Using the SDAdes
assembler, we produce assemblies of high quality from these filtered datasets
in a fraction of the time needed for an assembly from the datasets filtered
with Diginorm. We conclude that read filtering is a practical method for
reducing read data and for speeding up the assembly process. Our Bignorm
algorithm allows assemblies of competitive quality in comparison to Diginorm,
while being much faster. Bignorm is available for download at
https://git.informatik.uni-kiel.de/axw/Bignorm.git


  Understanding overfishing phenomenon and regulating fishing quotas is a major
global challenge for the 21st Century both in terms of providing food for
humankind and to preserve the oceans ecosystems. However, fishing is a complex
economic activity, affected not just by overfishing but also by such factors as
pollution, technology, financial factors and more. For this reason, it is often
difficult to state with complete certainty that overfishing is the cause of the
decline of a fishery. In this study, we developed a simple dynamic model based
on the earlier, well-known Lotka-Volterra model or Prey-Predator model. To
describe exploitation patterns, we assume that the fish stock and the fishing
industry are coupled stock variables in the model and they dynamically affect
each other, with the fishing yield proportional to both the fishing capital and
the fish stock. The model is based on the concept that the fishing industry
acts as the predator of the resource and that its growth and subsequent decline
is directly related to the abundance of the fish stock. If the model can be fit
historical data relative to specific fisheries, then it is a strong indication
that the fishing industry is strongly affected by the magnitude of the fish
stock and that, in particular, the decline of the yield and the decline of the
stock are linked to each other. The model does not pretend to be a general
description of the fishing industry in all its varied forms; however, the data
reported here show that the model can indeed qualitatively describe several
historical case of the collapse of fisheries. The model can also be used as a
qualitative guide to understand the behavior of several other fisheries. These
result indicate that one of the main factors causing the present crisis of the
world's fisheries is the overexploitation of the fish stocks.


  The present study revealed that populations of climbing perch (Anabas
testudineus) inhabiting river, backwater, shallow water channel, ponds with and
without vegetation cover, marsh, sewage canal and aquaculture tank varied
significantly in the number of food items consumed. Chironomous larvae, organic
debris and filamentous algae were the common ingredients of the menu of this
species across focal ecosystems, whereas sewage canal population was found
surviving solely on insect larvae and organic debris.


  In the traditional setup of public goods game all players are involved in
every available groups and the mutual benefit is shared among competing
cooperator and defector strategies. But in real life situations the group
formation of players could be more sophisticated because not all players are
attractive enough for others to participate in a joint venture. What if only
those players can initiate a group formation and establish a game who are
successful enough to the neighbors? To elaborate this idea we employ a modified
protocol and demonstrate that a carefully chosen threshold to establish joint
venture could efficiently improve the cooperation level even if the synergy
factor would suggest a full defector state otherwise. The microscopic mechanism
which is responsible for this effect is based on the asymmetric consequences of
competing strategies: while the success of a cooperator provides a long-time
well-being for the neighborhood, the temporary advantage of defection cannot be
maintained if the protocol is based on the success of leaders.


  Influenza is a contagious respiratory illness that causes significant human
morbidity and mortality, affecting 5-15% of the population in a typical
epidemic season. Human influenza epidemics are caused by types A and B, with
roughly 25% of human cases due to influenza B. Influenza B is a single-stranded
RNA virus with a high mutation rate, and both prior immune history and
vaccination put significant pressure on the virus to evolve. Due to the high
rate of viral evolution, the influenza B vaccine component of the annual
influenza vaccine is updated, roughly every other year in recent years. To
predict when an update to the vaccine is needed, an estimate of expected
vaccine effectiveness against a range of viral strains is required. We here
introduce a method to measure antigenic distance between the influenza B
vaccine and circulating viral strains. The measure correlates well with
effectiveness of the influenza B component of the annual vaccine in humans
between 1979 and 2014. We discuss how this measure of antigenic distance may be
used in the context of annual influenza vaccine design and prediction of
vaccine effectiveness.


  The males of the specie of frogs Engystomops pustulosus produce simple and
com- plex calls to lure females, as a way of Intersexual selection. Complex
calls lead males to a greater reproductive success than simple calls do.
However, the complex calls are also more attractive to the main predator of
these amphibians, the bat Trachops cirrhosus. Therefore, as M. Ryan suggests,
the complexity of the calls let the frogs keep a trade off between reproductive
success and predation. In this paper, we first propose to model the proportion
of simple to complex calls as a symmetric game of two strategies. We also
propose a model with three strategies (simple callers, complex callers and
quiet males), where we assess the effect of a male that keeps quiet and
intercepts females, which would play a role of Intrasexual selection. We
analyze the stable points of the replicator equations of the models that we
propose. Under the assumption that the decision of the males takes into account
this trade off between reproductive success and predation, our model reproduces
the observed behavior reported in the literature with minimal assumption on the
parameters. From the three strategies model, we verify that the quiet strategy
could only coexists with the simple and complex strategies as long as the rate
at which quiet males intercept females is high. We conclude that the
reproductive strategy of the male frog Engystomops pustulosus is rational.


  This is a method for discrete event simulation specified by survival
analysis. It presents a sequence of steps. First, hazard rates from survival
analysis specify the rates of a set of counting processes. Second, those
counting processes define a transition kernel. Third, there are four different
ways to sample that transition kernel, including a first-principles derivation
of exact stochastic simulation algorithms (SSA) in continuous time. This
simulation allows time-dependent intensities which include both continuous and
atomic components. Separating the steps involved makes a clear correspondence
between mathematical formulation and algorithmic implementation.


  In order to better understand the process of gene translation, the ribosome
flow model (RFM) with pool was introduced recently. This model describes the
movement of several ribosomes along an mRNA template and simultaneously
captures the dynamics of the finite pool of ribosomes. Studying this system
with respect to the number and stability of its equilibria was so far based on
monotone systems theory (Margaliot and Tuller, 2012). We extend the results
obtained therein by using a geometric approach, showing that the equilibria of
the system constitute a normally hyperbolic invariant submanifold.
Subsequently, we analyze the Jacobi linearization of the system evaluated at
the equilibria in order to show that the equilibria are asymptotically stable
relative to certain affine subspaces. As this approach does not require any
monotonicity features of the system, it may also be applied for more complex
systems of the same kind such as bi-directional ribosome flows or time-varying
template numbers.


  Reliable in-silico design of synthetic gene networks necessitates novel
approaches to model the process of protein synthesis under the influence of
limited resources. We present such a novel protein synthesis model which
originates from the Ribosome Flow Model and among other things describes the
movement of RNA-polymerase and Ribosomes on mRNA and DNA templates
respectively. By analyzing the convergence properties of this model based upon
geometric considerations we present additional insights into the dynamic
mechanisms of the process of protein synthesis. Further, we exemplarily show
how this model can be used to evaluate the performance of synthetic gene
circuits under different loading scenarios.


  McNamara, Trimmer, and Houston (2012) claim to provide an explanation of
certain systematic deviations from rational behavior using a mechanism that
could arise through natural selection. We provide an arguably much simpler
mechanism in terms of computational limitations, that performs better in the
environment described by McNamara, Trimmer, and Houston (2012). To argue
convincingly that animals' use of state-dependent valuation is adaptive and is
likely to be selected for by natural selection, one must argue that, in some
sense, it is a better approach than the simple strategies that we propose.


  The RNA-sequencing (RNA-seq) is becoming increasingly popular for quantifying
gene expression levels. Since the RNA-seq measurements are relative in nature,
between-sample normalization of counts is an essential step in differential
expression (DE) analysis. The normalization of existing DE detection algorithms
is ad hoc and performed once for all prior to DE detection, which may be
suboptimal since ideally normalization should be based on non-DE genes only and
thus coupled with DE detection. We propose a unified statistical model for
joint normalization and DE detection of log-transformed RNA-seq data.
Sample-specific normalization factors are modeled as unknown parameters in the
gene-wise linear models and jointly estimated with the regression coefficients.
By imposing sparsity-inducing L1 penalty (or mixed L1/L2-norm for multiple
treatment conditions) on the regression coefficients, we formulate the problem
as a penalized least-squares regression problem and apply the augmented
lagrangian method to solve it. Simulation studies show that the proposed model
and algorithms outperform existing methods in terms of detection power and
false-positive rate when more than half of the genes are differentially
expressed and/or when the up- and down-regulated genes among DE genes are
unbalanced in amount.


  A simple theory of ion permeation through a channel is presented, in which
diffusion occurs according to Fick's law and drift according to Ohm's law, in
the electric field determined by all the charges present. This theory accounts
for permeation in the channels studied to date in a wide range of solutions.
Interestingly, the theory works because the shape of the electric field is a
sensitive function of experimental conditions, e.g., ion concentration. Rate
constants for flux are sensitive functions of ionic concentration because the
fixed charge of the channel protein is shielded by the ions in an near it. Such
shielding effects are not included in traditional theories of ionic channels,
or other proteins, for that matter.


  Pan-genome analysis is a standard procedure to decipher genome heterogeneity
and diversification of bacterial species. Specie evolution is traced by
defining and comparing the core (conserved), accessory (dispensable) and unique
(strain-specific) gene pool with other strains of interest. Here, we present
pan-genome analysis of the genus Serratia, comprising of a dataset of 100
genomes. The isolates have clinical to environmental origin and consist of ten
different species from the genus, along with two subspecies of the
representative strain Serratia marcescens. Out of 19430 non-redundant coding
DNA sequences (CDS) from the dataset, 972 (5%) belonged to the core genome.
Majority of these genes were linked to metabolic function, followed by cellular
processes/signalling, information storage/processing while rest of them were
poorly characterized. 10,135 CDSs (52.16%) were associated with dispensible
genome while 8,321 CDSs (42.82%) were singletons or strain specific. The
Pan-genome orthologs indicated a positive correlation to the number of genomes
whereas negative correlation was obtained for core genome. Genomes were aligned
to obtain information about synteny, insertion/inversion, deletion and
duplications. This study provides insights into variation of Serratia species
and paves way for pan-genome analysis of other bacterial species at genus
level.


  Knots in proteins have been proposed to resist proteasomal degradation. Ample
evidence associates proteasomal degradation with neurodegeneration. One
interesting possibility is that indeed knotted conformers stall this machinery
leading to toxicity. However, although the proteasome is known to unfold
mechanically its substrates, at present there are no experimental methods to
emulate this particular traction geometry. Here, we consider several dynamical
models of the proteasome in which the complex is represented by an effective
potential with an added pulling force. This force is meant to induce
translocation of a protein or a polypeptide into the catalytic chamber. The
force is either constant or applied periodically. The translocated proteins are
modelled in a coarse-grained fashion. We do comparative analysis of several
knotted globular proteins and the transiently knotted polyglutamine tracts of
length 60 alone and fused in exon 1 of the huntingtin protein. Huntingtin is
associated with Huntington disease, a well-known genetically-determined
neurodegenerative disease. We show that the presence of a knot hinders and
sometimes even jams translocation. We demonstrate that the probability to do so
depends on the protein, the model of the proteasome, the magnitude of the
pulling force, and the choice of the pulled terminus. In any case, the net
effect would be a hindrance in the proteasomal degradation process in the cell.
This would then yield toxicity \textit{via} two different mechanisms: one
through toxic monomers compromising degradation and another by the formation of
toxic oligomers. Our work paves the way to the mechanistic investigation of the
mechanical unfolding of knotted structures by the proteasome and its relation
to toxicity and disease.


  The influenza pandemic of 1918-1919 killed at least 50 million people. The
reasons why this pandemic was so deadly remain largely unknown. However, It has
been shown that the 1918 viral hemagglutinin allows to reproduce the hallmarks
of the illness observed during the original pandemic. Thanks to the wealth of
hemagglutinin sequences accumulated over the last decades, amino-acid
substitutions that are found in the 1918-1919 sequences but rare otherwise can
be identified with high confidence. Such an analysis reveals that Gly188, which
is located within a key motif of the receptor binding site, is so rarely found
in hemagglutinin sequences that, taken alone, it is likely to be deleterious.
Monitoring this singular mutation in viral sequences may help prevent another
dramatic pandemic.


  We discuss the population dynamics with selection and random diffusion,
keeping the total population constant, in a fitness landscape associated with
Constraint Satisfaction, a paradigm for difficult optimization problems. We
obtain a phase diagram in terms of the size of the population and the diffusion
rate, with a glass phase inside which the dynamics keeps searching for better
configurations, and outside which deleterious `mutations' spoil the
performance. The phase diagram is analogous to that of dense active matter in
terms of temperature and drive.


  We develop a stochastic kinetic model of a pre-formed attachment of a
mictrotuble (MT) with a cell cortex, in which the MT is tethered to the cell by
a group of active motor proteins. Such an attachment is a particularly unique
case of ligand-receptor bonds: The MT ligand changes its length (and thus
binding sites) with time by polymerization-depolymerization kinetics, while
multiple motor receptors tend to walk actively along the MT length. These
processes, combined with force-mediated unbinding of the motors, result in an
elaborate behavior of the MT connection to the cell cortex. We present results
for the strength and lifetime of the system through the well-established
force-clamp and force-ramp protocols when external tension is applied to the
MT. The simulation results reveal that the MT-cell attachment behaves as a
catch-bond or slip-bond depending on system parameters. We provide analytical
approximations of the lifetime and discuss implications of our results on
in-vitro experiments.


  In this paper we explore the very important case of finding a health measure
in the lines of the survival curve but independent of the standard deviation
parameter. This is done by estimating the health state curve and calculating
the total area under the curve. An interesting comparison with the survival
curves is done. The health status measure (HSM or HS) resulting is of the form
of the life expectancy (LE) as it is expressed in terms of years of age. The HS
is independent of the standard deviation. When a perfect rectangularization of
the survival curve appear, the LE is equal to the age where the health state
curve approaches zero. The provided HS form gives rise to an interesting
classification of various countries. The results are close to those obtained by
the LE estimator with an interesting reorganization of country ranking. We also
provide illustrations of our estimated Health Status along with comparative
presentations of the Healthy Life Expectancy and the HALE measures in several
countries. Theoretical issues are provided and important stochastic simulations
are done along with reproduction of the death probability density forms.


  Understanding the evolutionary relationship among species is of fundamental
importance to the biological sciences. The location of the root in any
phylogenetic tree is critical as it gives an order to evolutionary events. None
of the popular models of nucleotide evolution used in likelihood or Bayesian
methods are able to infer the location of the root without exogenous
information. It is known that the most general Markov models of nucleotide
substitution can also not identify the location of the root or be fitted to
multiple sequence alignments with less than three sequences. We prove that the
location of the root and the full model can be identified and statistically
consistently estimated for a non-stationary, strand-symmetric substitution
model given a multiple sequence alignment with two or more sequences. We also
generalise earlier work to provide a practical means of overcoming the
computationally intractable problem of labelling hidden states in a
phylogenetic model.


  Population dynamics is constrained by the environment, which needs to obey
certain conditions to support population growth. We consider a standard model
for the evolution of a single species population density, that includes
reproduction, competition for resources and spatial spreading, while subject to
an external harmful effect. The habitat is spatially heterogeneous, there
existing a refuge where the population can be protected. Temporal variability
is introduced by the intermittent character of the refuge. This scenario can
apply to a wide range of situations, from a lab setting where bacteria can be
protected by a blinking mask from ultraviolet radiation, to large scale
ecosystems, like a marine reserve where there can be seasonal fishing
prohibitions. Using analytical and numerical tools, we investigate the
asymptotic behavior of the total population as a function of the size and
characteristic time scales of the refuge. We obtain expressions for the minimal
size required for population survival, in the slow and fast time scale limits.


  Humans, like all organisms, are subject to fundamental biophysical laws. Van
Valen predicted that, because of zero-sum dynamics, all populations of all
species in a given environment flux the same amount of energy on average.
Damuth's 'energetic equivalence rule' supported Van Valen's conjecture by
showing a trade off between few big animals per area with high individual
metabolic rates compared to abundant small species with low energy
requirements. We use established metabolic scaling theory to compare variation
in densities and individual energy use in human societies to other land
mammals. We show that hunter-gatherers occurred at lower densities than a
mammal of our size. Most modern humans, in contrast, concentrate in large
cities at densities that are up to four orders of magnitude greater than
hunter-gatherers yet cities consume up to two orders of magnitude greater
energy per capita. Today, cities across the globe flux greater energy than net
primary productivity on a per area basis. This is possible through enormous
fluxes of energy and materials across urban boundaries to sustain hyper-dense,
modern humans. The metabolic rift with nature created by hyper-dense cities
supported by fossil fuel energy poses formidable challenges for establishing a
sustainable relationship on a rapidly urbanizing, yet finite planet.


  Summary: Protein quality assessment is a long-standing problem in
bioinformatics. For more than a decade we have developed state-of-art
predictors by carefully selecting and optimising inputs to a machine learning
method. The correlation has increased from 0.60 in ProQ to 0.81 in ProQ2 and
0.85 in ProQ3 mainly by adding a large set of carefully tuned descriptions of a
protein. Here, we show that a substantial improvement can be obtained using
exactly the same inputs as in ProQ2 or ProQ3 but replacing the support vector
machine by a deep neural network. This improves the Pearson correlation to 0.90
(0.85 using ProQ2 input features).
  Availability: ProQ3D is freely available both as a webserver and a
stand-alone program at http://proq3.bioinfo.se/


  The spatiotemporal oscillation patterns of the proteins MinD and MinE are
used by the bacterium E. coli to sense its own geometry. Strikingly, both
computer simulations and experiments have recently shown that for the same
geometry of the reaction volume, different oscillation patterns can be stable,
with stochastic switching between them. Here we use particle-based Brownian
dynamics simulations to predict the relative frequency of different oscillation
patterns over a large range of three-dimensional compartment geometries, in
excellent agreement with experimental results. Fourier analyses as well as
pattern recognition algorithms are used to automatically identify the different
oscillation patterns and the switching rates between them. We also identify
novel oscillation patterns in three-dimensional compartments with
membrane-covered walls and identify a linear relation between the bound
Min-protein densities and the volume-to-surface ratio. In general, our work
shows how geometry sensing is limited by multistability and stochastic
fluctuations.


  Gene regulatory networks with dynamics characterized by multiple stable
states underlie cell fate-decisions. Quantitative models that can link
molecular-level knowledge of gene regulation to a global understanding of
network dynamics have the potential to guide cell-reprogramming strategies.
Networks are often modeled by the stochastic Chemical Master Equation, but
methods for systematic identification of key properties of the global dynamics
are currently lacking. We present a method for analyzing global dynamics of
gene networks using the Markov State Model (MSM) framework, which utilizes a
separation-of-timescales based clustering method to obtain a coarse-grained
state transition graph that approximates global gene network dynamics. The
method identifies the number, phenotypes, and lifetimes of long-lived network
states. Application of transition path theory to the constructed MSM decomposes
global dynamics into a set of dominant transition paths and associated relative
probabilities for stochastic state-switching. In this proof-of-concept study,
we found that the MSM provides a general framework for analyzing and
visualizing stochastic multistability and state-transitions in gene networks.
Our results suggest that the MSM framework, adopted from the field of atomistic
Molecular Dynamics, can be a useful tool for quantitative Systems Biology at
the network scale.


  Reliability of kinetic parameters are crucial in understanding enzyme
kinetics within cellular system. The present study suggests a few cautions that
need introspection for estimation of parameters like K(M), V(max) and K(I)
using Lineweaver-Burk plots. The quality of IC(50) too needs a thorough
reinvestigation because of its direct link with K(I) and K(M) values.
Inhibition kinetics under both steady-state and non-steady-state conditions are
studied and errors in estimated parameters are compared against actual values
to settle the question of their adequacy.


  Although a number of studies have shown that natural and laboratory
populations initially well-adapted to their environment can evolve rapidly when
conditions suddenly change, the dynamics of rapid adaptation are not well
understood. Here a population genetic model of polygenic selection is analyzed
to describe the short-term response of a quantitative trait after a sudden
shift of the phenotypic optimum. We provide explicit analytical expressions for
the time scales over which the trait mean approaches the new optimum. We find
that when the effect sizes are small relative to a scaled mutation rate, the
genomic signatures of polygenic selection are small to moderate allele
frequency changes that occur in the short-term phase in a synergistic fashion.
In contrast, selective sweeps, i.e., dramatic changes in the allele frequency
may occur provided the size of the effect is sufficiently large. Applications
of our theoretical results to the relationship between QTL and selective sweep
mapping and to tests of fast polygenic adaptation are discussed.


  We study the binding of plant hormone IAA on its receptor TIR1 introducing a
novel computational method that we call tomographic docking and that accounts
for interactions occurring along the depth of the binding pocket. Our results
suggest that selectivity is related to constraints that potential ligands
encounter on their way from the surface of the protein to their final position
at the pocket bottom. Tomographic docking helps develop specific hypotheses
about ligand binding, distinguishing binders from non-binders, and suggests
that binding is a three-step mechanism, consisting of engagement with a niche
in the back wall of the pocket, interaction with a molecular filter which
allows or precludes further descent of ligands, and binding on the pocket base.
Only molecules that are able to descend the pocket and bind at its base allow
the co-receptor IAA7 to bind on the complex, thus behaving as active auxins.
Analyzing the interactions at different depths, our new method helps in
identifying critical residues that constitute preferred future study targets
and in the quest for safe and effective herbicides. Also, it has the potential
to extend the utility of docking from ligand searches to the study of processes
contributing to selectivity.


  Over the years, advances in experimental and computational methods have
helped us to understand the role of thermodynamic, kinetic and active
(chaperone-aided) effects in coordinating the folding steps required to
achieving a knotted native state. Here, we review such developments by paying
particular attention to the complementarity of experimental and computational
studies. Key open issues that could be tackled with either or both approaches
are finally pointed out.


  We introduce phylodyn, an R package for phylodynamic analysis based on gene
genealogies. The package main functionality is Bayesian nonparametric
estimation of effective population size fluctuations over time. Our
implementation includes several Markov chain Monte Carlo-based methods and an
integrated nested Laplace approximation-based approach for phylodynamic
inference that have been developed in recent years. Genealogical data describe
the timed ancestral relationships of individuals sampled from a population of
interest. Here, individuals are assumed to be sampled at the same point in time
(isochronous sampling) or at different points in time (heterochronous
sampling); in addition, sampling events can be modeled with preferential
sampling, which means that the intensity of sampling events is allowed to
depend on the effective population size trajectory. We assume the coalescent
and the sequentially Markov coalescent processes as generative models of
genealogies. We include several coalescent simulation functions that are useful
for testing our phylodynamics methods via simulation studies. We compare the
performance and outputs of various methods implemented in phylodyn and outline
their strengths and weaknesses. R package phylodyn is available at
https://github.com/mdkarcher/phylodyn.


  Motivated by issues arising in population dynamics, we consider the problem
of iterating a given analytic function a number of times. We use the celebrated
technique known as Carleman linearization that turns (for a certain class of
functions) this problem into simply taking the power of a real number. We
expand this method, showing in particular that it can be used for population
models with immigration, and we also apply it to the famous logistic map. We
also are able to give a number of results for the invariant density of this
map, some being related to the Carleman linearization.


  We study the dynamics of the out-of-equilibrium nonlinear q-voter model with
two types of susceptible voters and zealots, introduced in [EPL 113, 48001
(2016)]. In this model, each individual supports one of two parties and is
either a susceptible voter of type $q_1$ or $q_2$, or is an inflexible zealot.
At each time step, a $q_i$-susceptible voter ($i = 1,2$) consults a group of
$q_i$ neighbors and adopts their opinion if all group members agree, while
zealots are inflexible and never change their opinion. This model violates
detailed balance whenever $q_1 \neq q_2$ and is characterized by two distinct
regimes of low and high density of zealotry. Here, by combining analytical and
numerical methods, we investigate the non-equilibrium stationary state of the
system in terms of its probability distribution, non-vanishing currents and
unequal-time two-point correlation functions. We also study the switching times
properties of the model by exploiting an approximate mapping onto the model of
[Phys. Rev. E 92, 012803 (2015)] that satisfies the detailed balance, and also
outline some properties of the model near criticality.


  The two primary causal dimensions of age-related disease are rate and
function. Change in rate of disease development shifts the age of onset. Change
in physiological function provides necessary steps in disease progression. A
causal factor may alter the rate of physiological change, but that causal
factor itself may have no direct physiological role. Alternatively, a causal
factor may provide a necessary physiological function, but that causal factor
itself may not alter the rate of disease onset. The rate-function duality
provides the basis for solving puzzles of age-related disease. Causal factors
of cancer illustrate the duality between rate processes of discovery, such as
somatic mutation, and necessary physiological functions, such as invasive
penetration across tissue barriers. Examples from cancer suggest general
principles of age-related disease.


  The motor neuron disease amyotrophic lateral sclerosis (ALS) typically begins
with localized muscle weakness. Progressive, widespread paralysis often follows
over a few years. Does the disease begin with local changes in a small piece of
neural tissue and then spread? Or does neural decay happen independently across
diverse spatial locations? The distinction matters, because local initiation
may arise by local changes in a tissue microenvironment, by somatic mutation,
or by various epigenetic or regulatory fluctuations in a few cells. A local
trigger must be coupled with a mechanism for spread. By contrast, independent
decay across spatial locations cannot begin by a local change, but must depend
on some global predisposition or spatially distributed change that leads to
approximately synchronous decay. This article outlines the conceptual frame by
which one contrasts local triggers and spread versus parallel spatially
distributed decay. Various neurodegenerative diseases differ in their
mechanistic details, but all can usefully be understood as falling along a
continuum of interacting local and global processes. Cancer provides an example
of disease progression by local triggers and spatial spread, setting a
conceptual basis for clarifying puzzles in neurodegeneration. Heart disease
also has crucial interactions between global processes, such as circulating
lipid levels, and local processes in the development of atherosclerotic
plaques. The distinction between local and global processes helps to understand
these various age-related diseases.


  Recent methodological developments have facilitated the integration of
high-throughput data into genome-scale models to obtain context-specific
metabolic reconstructions. A unique solution to this data integration problem
often may not be guaranteed, leading to a multitude of context-specific
predictions equally concordant with the integrated data. Yet, little attention
has been paid to the alternative optima resulting from the integration of
context-specific data. Here we present computational approaches to analyze
alternative optima for different context-specific data integration instances.
By using these approaches on metabolic reconstructions for the leaf of
Arabidopsis thaliana and the human liver, we show that the analysis of
alternative optima is key to adequately evaluating the specificity of the
predictions in particular cellular contexts. While we provide several ways to
reduce the ambiguity in the context-specific predictions, our findings indicate
that the existence of alternative optimal solutions warrant caution in detailed
context-specific analyses of metabolism.


  We develop a theoretical framework on the mechanism of combinatorial binding
of transcription factors (TFs) with their specific binding sites on DNA. We
consider three possible mechanisms viz. monomer, hetero-oligomer and
coordinated recruitment pathways. In the monomer pathway, combinatorial TFs
search for their targets in an independent manner and the protein-protein
interactions among them will be insignificant. The protein-protein interactions
are very strong so that the hetero-oligomer complex of TFs as a whole searches
for the cognate sites in case of hetero-oligomer pathway. The TF which arrived
first will recruit the adjacent TFs in a sequential manner in the recruitment
pathway. The free energy released from the protein-protein interactions among
TFs will be in turn utilized to stabilize the TFs-DNA complex. Such coordinated
binding of TFs in fact emerges as the cooperative effect. Monomer and
hetero-oligomer pathways are efficient only when few TFs are involved in the
combinatorial regulation. Detailed random walk simulations suggest that when
the number of TFs in a combination increases then the searching efficiency of
TFs in these pathways decreases with the increasing number of TFs in a power
law manner. The power law exponent associated with the monomer pathway seems to
be strongly dependent on the number of TFs, distance between the initial
position of TFs from their specific binding sites and the hop size associated
with the dynamics of TFs on DNA.


  The Sherman-Woodbury-Morrison (SWM) formula gives an explicit formula for the
inverse perturbation of a matrix in terms of the inverse of the original matrix
and the perturbation. This formula is useful for numerical applications. We
have produced similar results, giving an expression for the inverse of a matrix
when the $i$th row and column are removed. However, our expression involves
taking a limit, which inhibits use in similar applications as the SWM formula.
However, using our expression to find an analytical result on the spectral
radius of a special product of two matrices leads to an application. In
particular, we find a way to compute the fundamental reproductive ratio of a
relapsing disease being spread by a vector among two species of host that
undergo a different number of relapses.


  Markov state models (MSMs) and Master equation models are popular approaches
to approximate molecular kinetics, equilibria, metastable states, and reaction
coordinates in terms of a state space discretization usually obtained by
clustering. Recently, a powerful generalization of MSMs has been introduced,
the variational approach (VA) of molecular kinetics and its special case the
time-lagged independent component analysis (TICA), which allow us to
approximate slow collective variables and molecular kinetics by linear
combinations of smooth basis functions or order parameters. While it is known
how to estimate MSMs from trajectories whose starting points are not sampled
from an equilibrium ensemble, this has not yet been the case for TICA and the
VA. Previous estimates from short trajectories, have been strongly biased and
thus not variationally optimal. Here, we employ Koopman operator theory and
ideas from dynamic mode decomposition (DMD) to extend the VA and TICA to
non-equilibrium data. The main insight is that the VA and TICA provide a
coefficient matrix that we call Koopman model, as it approximates the
underlying dynamical (Koopman) operator in conjunction with the basis set used.
This Koopman model can be used to compute a stationary vector to reweight the
data to equilibrium. From such a Koopman-reweighted sample, equilibrium
expectation values and variationally optimal reversible Koopman models can be
constructed even with short simulations. The Koopman model can be used to
propagate densities, and its eigenvalue decomposition provide estimates of
relaxation timescales and slow collective variables for dimension reduction.
Koopman models are generalizations of Markov state models, TICA and the linear
VA and allow molecular kinetics to be described without a cluster
discretization.


  Cyanobacteria are an integral part of the Earth's biogeochemical cycles and a
promising resource for the synthesis of renewable bioproducts from atmospheric
CO2 . Growth and metabolism of cyanobacteria are inherently tied to the diurnal
rhythm of light availability. As yet, however, insight into the stoichiometric
and energetic constraints of cyanobacterial diurnal growth is limited. Here, we
develop a computational platform to evaluate the optimality of diurnal
phototrophic growth using a high-quality genome-scale metabolic reconstruction
of the cyanobacterium Synechococcus elongatus PCC 7942. We formulate
phototrophic growth as a self-consistent autocatalytic process and evaluate the
resulting time-dependent resource allocation problem using constraint-based
analysis. Based on a narrow and well defined set of parameters, our approach
results in an ab initio prediction of growth properties over a full diurnal
cycle. In particular, our approach allows us to study the optimality of
metabolite partitioning during diurnal growth. The cyclic pattern of glycogen
accumulation, an emergent property of the model, has timing characteristics
that are shown to be a trade-off between conflicting cellular objectives. The
approach presented here provides insight into the time-dependent resource
allocation problem of phototrophic diurnal growth and may serve as a general
framework to evaluate the optimality of metabolic strategies that evolved in
photosynthetic organisms under diurnal conditions.


  Studying metabolic networks is vital for many areas such as novel drugs and
bio-fuels. For biologists, a key challenge is that many reactions are
impractical or expensive to be found through experiments. Our task is to
recover the missing reactions. By exploiting the problem structure, we model
reaction recovery as a hyperlink prediction problem, where each reaction is
regarded as a hyperlink connecting its participating vertices (metabolites).
Different from the traditional link prediction problem where two nodes form a
link, a hyperlink can involve an arbitrary number of nodes. Since the
cardinality of a hyperlink is variable, existing classifiers based on a fixed
number of input features become infeasible. Traditional methods, such as common
neighbors and Katz index, are not applicable either, since they are restricted
to pairwise similarities. In this paper, we propose a novel hyperlink
prediction algorithm, called Matrix Boosting (MATBoost). MATBoost conducts
inference jointly in the incidence space and adjacency space by performing an
iterative completion-matching optimization. We carry out extensive experiments
to show that MATBoost achieves state-of-the-art performance. For a metabolic
network with 1805 metabolites and 2583 reactions, our algorithm can
successfully recover nearly 200 reactions out of 400 missing reactions.


  We consider a model for an infectious disease in the onset of an outbreak. We
introduce contact tracing incorporating a tracing delay. The effect of
randomness in the delay and the effect of the length of this delay in
comparison to the infectious period of the disease respectively to a latency
period on the effect of tracing, given e.g. by the change of the reproduction
number, is analyzed. We focus particularly on the effect of randomness in the
tracing delay.


  Evolutionary game theory is a mathematical approach to studying how social
behaviors evolve. In many recent works, evolutionary competition between
strategies is modeled as a stochastic process in a finite population. In this
context, two limits are both mathematically convenient and biologically
relevant: weak selection and large population size. These limits can be
combined in different ways, leading to potentially different results. We
consider two orderings: the wN limit, in which weak selection is applied before
the large population limit, and the Nw limit, in which the order is reversed.
Formal mathematical definitions of the Nw and wN limits are provided. Applying
these definitions to the Moran process of evolutionary game theory, we obtain
asymptotic expressions for fixation probability and conditions for success in
these limits. We find that the asymptotic expressions for fixation probability,
and the conditions for a strategy to be favored over a neutral mutation, are
different in the Nw and wN limits. However, the ordering of limits does not
affect the conditions for one strategy to be favored over another.


  We construct two ordinary-differential-equation models of a predator feeding
adaptively on two prey types, and we evaluate the models' ability to fit data
on freshwater plankton. We model the predator's switch from one prey to the
other in two different ways: (1) smooth switching using a hyperbolic tangent
function; and (2) by incorporating a parameter that changes abruptly across the
switching boundary as a system variable that is coupled to the population
dynamics. We conduct linear stability analyses, use approximate Bayesian
computation (ABC) combined with a population Monte Carlo (PMC) method to fit
model parameters, and compare model predictions quantitatively to data for
ciliate predators and their two algal prey groups collected from Lake Constance
on the German--Swiss--Austrian border. We show that the two models fit the data
well when the smooth transition is steep, supporting the simplifying assumption
of a discontinuous prey switching behavior for this scenario. We thus conclude
that prey switching is a possible mechanistic explanation for the observed
ciliate--algae dynamics in Lake Constance in spring, but that these data cannot
distinguish between the details of prey switching that are encoded in these
different models.


  Real food web data available in the literature presents us with the relations
between various species, sizes of these species, metabolic types of the species
and other useful information, which allows us to define parameters for the
mathematical dynamical models of these food webs, and perform theoretical
studies of these models. Unfortunately, the researches often face the problem
of the extinction of the species in such situations, which could be an
important limiting factor. In this paper, we propose a simple algorithm of
parameterisation that leads to the existence of positive steady state and
improves persistence of the species in multi-species models of real food webs.


  Intrinsic noise, the stochastic cell-to-cell fluctuations in mRNAs and
proteins, has been observed and proved to play important roles in cellular
systems. Due to the recent development in single-cell-level measurement
technology, the studies on intrinsic noise are becoming increasingly popular
among scholars. The chemical master equation (CME) has been used to model the
evolutions of complex chemical and biological systems since 1940, and are often
served as the standard tool for modeling intrinsic noise in gene regulation
system. A CME-based model can capture the discrete, stochastic, and dynamical
nature of gene regulation system, and may offer casual and physical explanation
of the observed data at single-cell level. Nonetheless, the complexity of CME
also pose serious challenge for researchers in proposing practical modeling and
inference frameworks. In this article, we will review the existing works on the
modelings and inference of intrinsic noise in gene regulation system within the
framework of CME model. We will explore the principles in constructing a CME
model for studying gene regulation system and discuss the popular
approximations of CME. Then we will study the simulation simulation methods as
well as the analytical and numerical approaches that can be used to obtain
solution to a CME model. Finally we will summary the exiting statistical
methods that can be used to infer the unknown parameters or structures in CME
model using single-cell-level gene expression data.


  To address the large gap between time scales that can be easily reached by
molecular simulations and those required to understand protein dynamics, we
propose a rapid self-consistent approximation of the side chain free energy at
every integration step. In analogy with the adiabatic Born-Oppenheimer
approximation for electronic structure, the protein backbone dynamics are
simulated as preceding according to the dictates of the free energy of an
instantaneously-equilibrated side chain potential. The side chain free energy
is computed on the fly, allowing the protein backbone dynamics to traverse a
greatly smoothed energetic landscape. This results in extremely rapid
equilibration and sampling of the Boltzmann distribution. Because our method
employs a reduced model involving single-bead side chains, we also provide a
novel, maximum-likelihood method to parameterize the side chain model using
input data from high resolution protein crystal structures. We demonstrate
state-of-the-art accuracy for predicting $\chi_1$ rotamer states while
consuming only milliseconds of CPU time. We also show that the resulting free
energies of side chains is sufficiently accurate for de novo folding of some
proteins.


  Mixed-species growth models are needed as a synthesis of ecological knowledge
and for guiding forest management. Individual-tree models have been commonly
used, but the difficulties of reliably scaling from the individual to the stand
level are often underestimated. Emergent properties and statistical issues
limit their effectiveness. A more holistic modelling of aggregates at the whole
stand level is a potentially attractive alternative. This work explores
methodology for developing biologically consistent dynamic mixture models where
the state is described by aggregate stand-level variables for species or
age/size cohorts. The methods are demonstrated and tested with a two-cohort
model for spruce-aspen mixtures named SAM. The models combine single-species
submodels and submodels for resource partitioning among the cohorts. The
partitioning allows for differences in competitive strength among species and
size classes, and for complementarity effects. Height growth reduction in
suppressed cohorts is also modelled. SAM fits well the available data, and
exhibits behaviors consistent with current ecological knowledge. The general
framework can be applied to any number of cohorts, and should be useful as a
basis for modelling other mixed-species or uneven-aged stands.


  Genetic variation in human populations is influenced by geographic ancestry
due to spatial locality in historical mating and migration patterns. Spatial
population structure in genetic datasets has been traditionally analyzed using
either model-free algorithms, such as principal components analysis (PCA) and
multidimensional scaling, or using explicit spatial probabilistic models of
allele frequency evolution. We develop a general probabilistic model and an
associated inference algorithm that unify the model-based and data-driven
approaches to visualizing and inferring population structure. Our algorithm,
Geographic Ancestry Positioning (GAP), relates local genetic distances between
samples to their spatial distances, and can be used for visually discerning
population structure as well as accurately inferring the spatial origin of
individuals on a two-dimensional continuum. On both simulated and several real
datasets from diverse human populations, GAP exhibits substantially lower error
in reconstructing spatial ancestry coordinates compared to PCA.
  Our spatial inference algorithm can also be effectively applied to the
problem of population stratification in genome-wide association studies (GWAS),
where hidden population structure can create fictitious associations when
population ancestry is correlated with both the genotype and the trait. We
develop an association test that uses the ancestry coordinates inferred by GAP
to accurately account for ancestry-induced correlations in GWAS. Based on
simulations and analysis of a dataset of 10 metabolic traits measured in a
Northern Finland cohort, which is known to exhibit significant population
structure, we find that our method has superior power to current approaches.


  Determination of protein concentration in often an absolute pre-requisite in
preparing samples for biochemical and proteomic analyses. However, current
protein assay methods are not compatible with both reducers and detergents,
which are however present simultaneously in most denaturing extraction buffers
used in proteomics and electrophoresis, and in particular in SDS
electrophoresis. We found that inclusion of cyclodextrins in a Coomassie
blue-based assay made it compatible with detergents, as cyclodextrins complex
detergents in a 1:1 molecular ratio. As this type of assay is intrinsically
resistant to reducers, we have thus developed a single step assay that is both
detergent and reducer compatible. Depending on the type and concentration of
detergents present in the sample buffer, either beta-cyclodextrin or
alpha-cyclodextrin can be used, the former being able to complex a wider range
of detergents and the latter being able to complex higher amounts of detergents
due to its greater solubility in water. Cyclodextrins are used at final
concentrations of 2-10 mg/mL in the assay mix. This typically allows to measure
samples containing as little as 0.1 mg/mL protein, in the presence of up to 2%
detergent and reducers such as 5 % mercaptoethanol or 50 mM DTT in a single
step with a simple spectrophotometric assay. This article is protected by
copyright. All rights reserved.


  One of the most important questions in game theory concerns how mutual
cooperation can be achieved and maintained in a social dilemma. In Axelrod's
tournaments of the iterated prisoner's dilemma, Tit-for-Tat (TFT) demonstrated
the role of reciprocity in the emergence of cooperation. However, the stability
of TFT does not hold in the presence of implementation error, and a TFT
population is prone to neutral drift to unconditional cooperation, which
eventually invites defectors. We argue that a combination of TFT and anti-TFT
(ATFT) overcomes these difficulties in a noisy environment, provided that ATFT
is defined as choosing the opposite to the opponent's last move. According to
this TFT-ATFT strategy, a player normally uses TFT; turns to ATFT upon
recognizing his or her own error; returns to TFT either when mutual cooperation
is recovered or when the opponent unilaterally defects twice in a row. The
proposed strategy provides simple and deterministic behavioral rules for
correcting implementation error in a way that cannot be exploited by the
opponent, and suppresses the neutral drift to unconditional cooperation.


  Given a gene tree and a species tree, ancestral configurations represent the
combinatorially distinct sets of gene lineages that can reach a given node of
the species tree. They have been introduced as a data structure for use in the
recursive computation of the conditional probability under the multispecies
coalescent model of a gene tree topology given a species tree, the cost of this
computation being affected by the number of ancestral configurations of the
gene tree in the species tree. For matching gene trees and species trees, we
obtain enumerative results on ancestral configurations. We study ancestral
configurations in balanced and unbalanced families of trees determined by a
given seed tree, showing that for seed trees with more than one taxon, the
number of ancestral configurations increases for both families exponentially in
the number of taxa $n$. For fixed $n$, the maximal number of ancestral
configurations tabulated at the species tree root node and the largest number
of labeled histories possible for a labeled topology occur for trees with
precisely the same unlabeled shape. For ancestral configurations at the root,
the maximum increases with $k_0^n$, where $k_0 \approx 1.5028$ is a quadratic
recurrence constant. Under a uniform distribution over the set of labeled trees
of given size, the mean number of root ancestral configurations grows with
$\sqrt{3/2}(4/3)^n$ and the variance with approximately $1.4048(1.8215)^n$. The
results provide a contribution to the combinatorial study of gene trees and
species trees.


  We consider the goal of predicting how complex networks respond to chronic
(press) perturbations when characterizations of their network topology and
interaction strengths are associated with uncertainty. Our primary result is
the derivation of exact formulas for the expected number and probability of
qualitatively incorrect predictions about a system's responses under
uncertainties drawn form arbitrary distributions of error. These formulas
obviate the current use of simulations, algorithms, and qualitative modeling
techniques. Additional indices provide new tools for identifying which links in
a network are most qualitatively and quantitatively sensitive to error, and for
determining the volume of errors within which predictions will remain
qualitatively determinate (i.e. sign insensitive). Together with recent
advances in the empirical characterization of uncertainty in ecological
networks, these tools bridge a way towards probabilistic predictions of network
dynamics.


  Phylogenetics, the inference of evolutionary trees from molecular sequence
data such as DNA, is an enterprise that yields valuable evolutionary
understanding of many biological systems. Bayesian phylogenetic algorithms,
which approximate a posterior distribution on trees, have become a popular if
computationally expensive means of doing phylogenetics. Modern data collection
technologies are quickly adding new sequences to already substantial databases.
With all current techniques for Bayesian phylogenetics, computation must start
anew each time a sequence becomes available, making it costly to maintain an
up-to-date estimate of a phylogenetic posterior. These considerations highlight
the need for an \emph{online} Bayesian phylogenetic method which can update an
existing posterior with new sequences.
  Here we provide theoretical results on the consistency and stability of
methods for online Bayesian phylogenetic inference based on Sequential Monte
Carlo (SMC) and Markov chain Monte Carlo (MCMC). We first show a consistency
result, demonstrating that the method samples from the correct distribution in
the limit of a large number of particles. Next we derive the first reported set
of bounds on how phylogenetic likelihood surfaces change when new sequences are
added. These bounds enable us to characterize the theoretical performance of
sampling algorithms by bounding the effective sample size (ESS) with a given
number of particles from below. We show that the ESS is guaranteed to grow
linearly as the number of particles in an SMC sampler grows. Surprisingly, this
result holds even though the dimensions of the phylogenetic model grow with
each new added sequence.


  The unwelcome evolution of malignancy during cancer progression emerges
through a selection process in a complex heterogeneous population structure. In
the present work, we investigate evolutionary dynamics in a phenotypically
heterogeneous population of stem cells (SCs) and their associated progenitors.
The fate of a malignant mutation is determined not only by overall stem cell
and differentiated cell growth rates but also differentiation and
dedifferentiation rates. We investigate the effect of such a complex population
structure on the evolution of malignant mutations. We derive exact analytic
results for the fixation probability of a mutant arising in each of the
subpopulations. The analytic results are in almost perfect agreement with the
numerical simulations. Moreover, a condition for evolutionary advantage of a
mutant cell versus the wild type population is given in the present study. We
also show that microenvironment-induced plasticity in invading mutants leads to
more aggressive mutants with higher fixation probability. Our model predicts
that decreasing polarity between stem and differentiated cells turnover would
raise the survivability of non-plastic mutants; while it would suppress the
development of malignancy for plastic mutants. We discuss our model in the
context of colorectal/intestinal cancer (at the epithelium). This novel
mathematical framework can be applied more generally to a variety of problems
concerning selection in heterogeneous populations, in other contexts such as
population genetics, and ecology.


  Corynebacterium glutamicum is a Gram-positive, anaerobic, rod-shaped soil
bacterium able to grow on a diversity of carbon sources like sugars and organic
acids. It is a biotechnological relevant organism because of its highly
efficient ability to biosynthesize amino acids, such as L-glutamic acid and
L-lysine. Here, we reconstructed the most complete C. glutamicum regulatory
network to date and comprehensively analyzed its global organizational
properties, systems-level features and functional architecture. Our analyses
show the tremendous power of Abasy Atlas to study the functional organization
of regulatory networks. We created two models of the C. glutamicum regulatory
network: all-evidences (containing both weak and strong supported interactions,
genomic coverage = 73%) and strongly-supported (only accounting for strongly
supported evidences, genomic coverage = 71%). Using state-of-the-art
methodologies, we prove that power-law behaviors truly govern the connectivity
and clustering coefficient distributions. We found a non-previously reported
circuit motif that we named complex feed-forward motif. We highlighted the
importance of feedback loops for the functional architecture, beyond whether
they are statistically over-represented or not in the network. We show that the
previously reported top-down approach is inadequate to infer the hierarchy
governing a regulatory network because feedback bridges different hierarchical
layers, and the top-down approach disregards the presence of intermodular genes
shaping the integration layer. Our findings all together further support a
diamond-shaped, three-layered hierarchy exhibiting some feedback between
processing and coordination layers, which is shaped by four classes of
systems-level elements: global regulators, locally autonomous modules, basal
machinery and intermodular genes.


  Investigations of topological uniqueness of gene interaction networks in
cancer cells are essential for understanding this disease. Based on the random
matrix theory, we study the distribution of the nearest neighbor level spacings
$P(s)$ of interaction matrices for gene networks in human cancer cells. The
interaction matrices are computed using the Cancer Network Galaxy (TCNG)
database, which is a repository of gene interactions inferred by a Bayesian
network model. 256 NCBI GEO entries regarding gene expressions in human cancer
cells have been selected for the Bayesian network calculations in TCNG. We
observe the Wigner distribution of $P(s)$ when the gene networks are dense
networks that have more than $\sim 38,000$ edges. In the opposite case, when
the networks have smaller numbers of edges, the distribution $P(s)$ becomes the
Poisson distribution. We investigate relevance of $P(s)$ both to the size of
the networks and to edge frequencies that manifest reliance of the inferred
gene interactions.


  Standard models assign disease progression to discrete categories or stages
based on well-characterized clinical markers. However, such a system is
potentially at odds with our understanding of the underlying biology, which in
highly complex systems may support a (near-)continuous evolution of disease
from inception to terminal state. To learn such a continuous disease score one
could infer a latent variable from dynamic "omics" data such as RNA-seq that
correlates with an outcome of interest such as survival time. However, such
analyses may be confounded by additional data such as clinical covariates
measured in electronic health records (EHRs). As a solution to this we
introduce covariate latent variable models, a novel type of latent variable
model that learns a low-dimensional data representation in the presence of two
(asymmetric) views of the same data source. We apply our model to TCGA
colorectal cancer RNA-seq data and demonstrate how incorporating
microsatellite-instability (MSI) status as an external covariate allows us to
identify genes that stratify patients on an immune-response trajectory.
Finally, we propose an extension termed Covariate Gaussian Process Latent
Variable Models for learning nonparametric, nonlinear representations. An R
package implementing variational inference for covariate latent variable models
is available at http://github.com/kieranrcampbell/clvm.


  Dynamical systems driven by nonlinear delay SDEs with small noise can exhibit
important rare events on long timescales. When there is no delay, classical
large deviations theory quantifies rare events such as escapes from metastable
fixed points. Near such fixed points, one can approximate nonlinear delay SDEs
by linear delay SDEs. Here, we develop a fully explicit large deviations
framework for (necessarily Gaussian) processes $X_t$ driven by linear delay
SDEs with small diffusion coefficients. Our approach enables fast numerical
computation of the action functional controlling rare events for $X_t$ and of
the most likely paths transiting from $X_0 = p$ to $X_T=q$. Via linear noise
local approximations, we can then compute most likely routes of escape from
metastable states for nonlinear delay SDEs. We apply our methodology to the
detailed dynamics of a genetic regulatory circuit, namely the co-repressive
toggle switch, which may be described by a nonlinear chemical Langevin SDE with
delay.


  In this study, spatial stochastic and logistic model (SSLM) describing
dynamics of a population of a certain species was analysed. The behaviour of
the extinction threshold as a function of model parameters was studied. More
specifically, we studied how the critical values for the model parameters that
separate the cases of extinction and persistence depend on the spatial scales
of the competition and dispersal kernels. We compared the simulations and
analytical results to examine if and how the mathematical approximations break
down at the vicinity of the extinction threshold, and found a functional form
of the naive approximation for which higher-order term of analytical
approximation converges.


  Frameshift translation is an important phenomenon that contributes to the
appearance of novel Coding DNA Sequences (CDS) and functions in gene evolution,
by allowing alternative amino acid translations of genes coding regions.
Frameshift translations can be identified by aligning two CDS, from a same gene
or from homologous genes, while accounting for their codon structure. Two main
classes of algorithms have been proposed to solve the problem of aligning CDS,
either by amino acid sequence alignment back-translation, or by simultaneously
accounting for the nucleotide and amino acid levels. The former does not allow
to account for frameshift translations and up to now, the latter exclusively
accounts for frameshift translation initiation, not accounting for the length
of the translation disruption caused by a frameshift.
  Here, we introduce a new scoring scheme with an algorithm for the pairwise
alignment of CDS accounting for frameshift translation initiation and length,
while simultaneously accounting for nucleotide and amino acid sequences. We
compare the method to other CDS alignment methods based on an application to
the comparison of pairs of CDS from homologous \emph{human}, \emph{mouse} and
\emph{cow} genes of ten mammalian gene families from the Ensembl-Compara
database. The results show that our method is particularly robust to parameter
changes as compared to existing methods. It also appears to be a good
compromise, performing well both in the presence and absence of frameshift
translations between the CDS. An implementation of the method is available at
https://github.com/UdeS-CoBIUS/FsePSA.


  Is it possible to reverse the low total fertility rate (TFR) in the developed
world? Using a hypothetical model of population we have analysed the decline of
the TFR which have took place in the background of ongoing global economic
changes, and a liberalization process after the end of the Cold War. These
phenomena have affected more than 110 millions of inhabitants of Central Europe
and the Baltics and approximately 80 millions of inhabitants in Germany. The
model has features of complex and evolving system of interacting individuals,
and it enables to investigate a broad spectrum of input factors on individual
decisions to limit the offspring. In the case of the TFR <1.5, our initial
analysis show a need of radical changes of the global economy that will
stimulate series of self-regulations of demographic processes and evolution
toward the safe TFR>2.1. The changes should stimulate more uniform spatial
distribution of wealth, capital and usage. They will increase a number of
self-sufficient and cooperative territories, to decrease the income inequality,
to decrease labour and social mobilities. Societies should investigate the
impacts of economic regulations and actions on the TFR trends in advance and
take into account a biological nature of women more responsible.


  The value of research containing novel combinations of molecules can be seen
in many innovative and award-winning research programs. Despite calls to use
innovative approaches to address common diseases, an increasing majority of
research funding goes toward "safe" incremental research. Counteracting this
trend by nurturing novel and potentially transformative scientific research is
challenging, it must be supported in competition with established research
programs. Therefore, we propose a tool that helps to resolve the tension
between safe but fundable research vs. high-risk but potentially
transformational research. It does this by identifying hidden overlapping
interest around novel molecular research topics. Specifically, it identifies
paths of molecular interactions that connect research topics and hypotheses
that would not typically be associated, as the basis for scientific
collaboration. Because these collaborations are related to the scientists'
present trajectory, they are low risk and can be initiated rapidly. Unlike most
incremental steps, these collaborations have the potential for leaps in
understanding, as they reposition research for novel disease applications. We
demonstrate the use of this tool to identify scientists who could contribute to
understanding the cellular role of genes with novel associations with
Alzheimer's disease, which have not been thoroughly characterized, in part due
to the funding emphasis on established research.


  We propose a new deterministic growth model which captures certain features
of both the Gompertz and Korf laws. We investigate its main properties, with
special attention to the correction factor, the relative growth rate, the
inflection point, the maximum specific growth rate, the lag time and the
threshold crossing problem. Some data analytic examples and their performance
are also considered. Furthermore, we study a stochastic counterpart of the
proposed model, that is a linear time-inhomogeneous birth-death process whose
mean behaves as the deterministic one. We obtain the transition probabilities,
the moments and the population ultimate extinction probability for this
process. We finally treat the special case of a simple birth process, which
better mimics the proposed growth model.


  ATP-driven proton pumps, which are critical to the operation of a cell,
maintain cytosolic and organellar pH levels within a narrow functional range.
These pumps employ two very different mechanisms: an elaborate rotary mechanism
used by V-ATPase H+ pumps, and a simpler alternating access mechanism used by
P-ATPase H+ pumps. Why are two different mechanisms used to perform the same
function? Systematic analysis, without parameter fitting, of kinetic models of
the rotary, alternating access and other possible mechanisms suggest that, when
the ratio of protons transported per ATP hydrolyzed exceeds one, the
one-at-a-time proton transport by the rotary mechanism is faster than other
possible mechanisms across a wide range of driving conditions. When the ratio
is one, there is no intrinsic difference in the free energy landscape between
mechanisms, and therefore all mechanisms can exhibit the same kinetic
performance. To our knowledge all known rotary pumps have an H+:ATP ratio
greater than one, and all known alternating access ATP-driven proton pumps have
a ratio of one. Our analysis suggests a possible explanation for this apparent
relationship between coupling ratio and mechanism. When the conditions under
which the pump must operate permit a coupling ratio greater than one, the
rotary mechanism may have been selected for its kinetic advantage. On the other
hand, when conditions require a coupling ratio of one or less, the alternating
access mechanism may have been selected for other possible advantages resulting
from its structural and functional simplicity.


  The advent of large scale, high-throughput genomic screening has introduced a
wide range of tests for diagnostic purposes. Prominent among them are tests
using miRNA expression levels. Genomics and proteomics now provide expression
levels of hundreds of miRNAs at a time. However, for actual diagnostic tools to
become reality requires the simultaneous development of methods to interpret
the large amounts of miRNA expression data that can be generated from a single
patient sample. Because these data are in numeric form, quantitative methods
must be developed. Statistics such as p-values and log fold change give some
insight, but the diagnostic effectiveness of each miRNA test must first be
evaluated. Here, the author has developed a traditional, sensitivity- and
specificity-based algorithm, as well as a modern machine learning algorithm,
and evaluated their diagnostic potential for lung cancer against a publicly
available database. The findings suggest that the machine learning algorithm
achieves higher accuracy (97% for cancerous and 73% for normal samples), in
addition to providing confidence intervals that could provide valuable
diagnostic support. The machine learning algorithm also has significant
potential for expansion to more complex diagnoses of lung cancer sub-types, to
other cancers as well diseases beyond cancer. Both algorithms are available on
the Github repo: https://github.com/neerja-g/machine-learning-miRNA.


  Motivated by single-molecule experiments on synaptic membrane protein
domains, we use a stochastic lattice model to study protein reaction and
diffusion processes in crowded membranes. We find that the stochastic
reaction-diffusion dynamics of synaptic proteins provide a simple physical
mechanism for collective fluctuations in synaptic domains, the molecular
turnover observed at synaptic domains, key features of the single-molecule
trajectories observed for synaptic proteins, and spatially inhomogeneous
protein lifetimes at the cell membrane. Our results suggest that central
aspects of the single-molecule and collective dynamics observed for membrane
protein domains can be understood in terms of stochastic reaction-diffusion
processes at the cell membrane.


  The living cell is an open nonequilibrium biochemical system, where ATP
hydrolysis serves as the energy source for a wide range of intracellular
processes including the assurance for decision-making. In the fission yeast
cell cycle, the transition from G2 phase to M phase is triggered by the
activation of Cdc13/Cdc2 and Cdc25, and the deactivation of Wee1. Each of these
three events involves a phosphorylation-dephosphorylation (PdP) cycle, and
together they form a regulatory circuit with feedback loops. Almost all
quantitative models for cellular networks in the past have invalid
thermodynamics due to the assumption of irreversible enzyme kinetics. We
constructed a thermodynamically realistic kinetic model of the G2/M circuit,
and show that the phosphorylation energy ($\Delta G$), which is determined by
the cellular ATP/ADP ratio, critically controls the dynamics and the bistable
nature of Cdc2 activation. Using fission yeast nucleoplasmic extract (YNPE), we
are able to experimentally verify our model prediction that increased , being
synergistic to the accumulation of Cdc13, drives the activation of Cdc2.
Furthermore, Cdc2 activation exhibits bistability and hysteresis in response to
changes in phosphorylation energy. These findings suggest that adequate
maintenance of phosphorylation energy ensures the bistability and robustness of
the activation of Cdc2 in the G2/M transition. Free energy might play a
widespread role in biological decision-making processes, connecting
thermodynamics with information processing in biology.


  We have performed all atom molecular dynamics simulations on the projection
domain of the intrinsically disordered htau40 protein. After generating a
suitable ensemble of starting conformations at high temperatures, at room
temperature in an adaptive box algorithm we have generated histograms for the
radius of gyration, secondary structure time series, generated model small
angle x-ray scattering intensities, and model chemical shift plots for
comparison to nuclear magnetic resonance data for solvated and filamentous tau.
Significantly, we find that the chemical shift spectrum is more consistent with
filamentous tau than full length solution based tau. We have also carried out
principle component analysis and find three basics groups: compact globules,
tadpoles, and extended hinging structures. To validate the adaptive box and our
force field choice, we have run limited simulations in a large conventional box
with varying force fields and find that our essential results are unchanged. We
also performed two simulations with the TIP4P-D water model, the effects of
which depended on whether the initial configuration was compact or extended.


  One of the properties that make ecological systems so unique is the range of
complex behavioural patterns that can be exhibited by even the simplest
communities with only a few species. Much of this complexity is commonly
attributed to stochastic factors which have very high-degrees of freedom.
Orthodox study of the evolution of these simple networks has generally been
limited in its ability to explain complexity, since it restricts evolutionary
adaptation to an inertia-free process with few degrees of freedom in which only
gradual, moderately complex behaviours are possible. We propose a model
inspired by particle mediated field phenomena in classical physics in
combination with fundamental concepts in adaptation, that suggests that small
but high-dimensional chaotic dynamics near to the adaptive trait optimum could
help explain complex properties shared by most ecological datasets, such as
aperiodicity and pink, fractal noise spectra. By examining a simple
predator-prey model and appealing to real ecological data, we show that this
type of complexity could be easily confused for or confounded by stochasticity,
especially when spurred on or amplified by stochastic factors that share
variational and spectral properties with the underlying dynamics.


  The architecture of eukaryotic coding genes allows the production of several
different protein isoforms by genes. Current gene phylogeny reconstruction
methods make use of a single protein product per gene, ignoring information on
alternative protein isoforms. These methods often lead to inaccurate gene tree
reconstructions that require to be corrected before being used in phylogenetic
tree reconciliation analyses or gene products phylogeny reconstructions. Here,
we propose a new approach for the reconstruction of accurate gene trees and
protein trees accounting for the production of alternative protein isoforms by
the genes of a gene family. We extend the concept of reconciliation to protein
trees, and we define a new reconciliation problem called MinDRGT that consists
in finding a gene tree that minimizes a double reconciliation cost with a given
protein tree and a given species tree. We define a second problem called
MinDRPGT that consists in finding a protein tree and a gene tree minimizing a
double reconciliation cost, given a species tree and a set of protein subtrees.
We provide algorithmic exact and heuristic solutions for some versions of the
problems, and we present the results of an application to the correction of
gene trees from the Ensembl database. An implementation of the heuristic method
is available at https://github.com/UdeS-CoBIUS/Protein2GeneTree.


  A Supertree synthesizes the topologies of a set of phylogenetic trees
carrying overlapping taxa set. In process, conflicts in the tree topologies are
aimed to be resolved with the consensus clades. Such a problem is proved to be
NP-hard. Various heuristics on local search, maximum parsimony, graph cut, etc.
lead to different supertree approaches, of which the most popular methods are
based on analyzing fixed size subtree topologies (such as triplets or
quartets). Time and space complexities of these methods, however, depend on the
subtree size considered. Our earlier proposed supertree method COSPEDTree, uses
evolutionary relationship among individual couplets (taxa pair), to produce
slightly conservative (not fully resolved) supertrees. Here we propose its
improved version COSPEDTree-II, which produces better resolved supertree with
lower number of missing branches, and incurs much lower running time. Results
on biological datasets show that COSPEDTree-II belongs to the category of high
performance and computationally efficient supertree methods.


  Microtubule asters - radial arrays of microtubules organized by centrosomes -
play a fundamental role in the spatial coordination of animal cells. The
standard model of aster growth assumes a fixed number of microtubules
originating from the centrosomes. However, aster morphology in this model does
not scale with cell size, and we recently found evidence for non-centrosomal
microtubule nucleation. Here, we combine autocatalytic nucleation and
polymerization dynamics to develop a biophysical model of aster growth. Our
model predicts that asters expand as traveling waves and recapitulates all
major aspects of aster growth. As the nucleation rate increases, the model
predicts an explosive transition from stationary to growing asters with a
discontinuous jump of the growth velocity to a nonzero value. Experiments in
frog egg extract confirm the main theoretical predictions. Our results suggest
that asters observed in large frog and amphibian eggs are a meshwork of short,
unstable microtubules maintained by autocatalytic nucleation and provide a
paradigm for the assembly of robust and evolvable polymer networks.


  The fast blood stream of animals is associated with large shear stresses.
Consequently, blood cells have evolved a special morphology and a specific
internal architecture allowing them to maintain their integrity over several
weeks. For instance, non-mammalian red blood cells, mammalian erythroblasts and
platelets have a peripheral ring of microtubules, called the marginal band,
that flattens the overall cell morphology by pushing on the cell cortex. In
this article, we model how the shape of these cells stems from the balance
between marginal band elasticity and cortical tension. We predict that the
diameter of the cell scales with the total microtubule polymer, and verify the
predicted law across a wide range of species. Our analysis also shows that the
combination of the marginal band rigidity and cortical tension increases the
ability of the cell to withstand forces without deformation. Finally, we model
the marginal band coiling that occurs during the disc-to-sphere transition
observed for instance at the onset of blood platelet activation. We show that
when cortical tension increases faster than crosslinkers can unbind, the
marginal band will coil, whereas if the tension increases slower, the marginal
band may shorten as microtubules slide relative to each other.


  We consider two approaches to modelling of cell metabolism of
6-mercaptopurine, which is one of the important chemotherapy drugs used for
treating of acute lymphocytic leukemia: kinetic ordinary differential equations
and random Boolean networks, and analyse their interplay with respect to taking
into account ATP concentration as a key parameter of switching between
different pathways. It is shown that Boolean networks, which allow for avoiding
complexity of general kinetic modelling, preserve an opportunity to the
reproduction of the principal switching mechanism. To keep a detailed
quantitative measure of the control parameter, a combined Boolean-ODE method is
proposed.


  A complex network approach to protein folding is proposed. The graph object
is the network of shortcut edges present in a native-state protein (SCN0).
Although SCN0s are found via an intuitive message passing algorithm (S.
Milgram, Psychology Today, May 1967 pp. 61-67), they are meaningful enough that
the logarithm form of their contact order (SCN0_lnCO) correlates significantly
with protein kinetic rates, regardless of protein size. Further, the clustering
coefficient of a SCN0 (CSCN0) can be used to combine protein segments
iteratively within the Restricted Binary Collision model to form the whole
native structure. This simple yet surprisingly effective strategy identified
reasonable folding pathways for 12 small single-domain two-state folders, and
three non-canonical proteins: ACBP (non-two-state), Top7 (non-cooperative) and
DHFR (non-single-domain, > 100 residues). For two-state folders, CSCN0 is
relatable to folding rates, transition-state placement and stability. The
influence of CSCN0 on folding extends to non-native structures. Moreover, SCN
analysis of non-native structures could suggest three fold success factors for
the fast folding Villin headpiece peptide. These results support the view of
protein folding as a bottom-up hierarchical process guided from above by
native-state topology, and could facilitate future constructive demonstrations
of this long held hypothesis for larger proteins.