  It is shown that the density of modes of the vibrational spectrum of globular
proteins is universal, i.e., regardless of the protein in question it closely
follows one universal curve. The present study, including 135 proteins analyzed
with a full atomic empirical potential (CHARMM22) and using the full complement
of all atoms Cartesian degrees of freedom, goes far beyond previous claims of
universality, confirming that universality holds even in the high-frequency
range (300- 4000 1/cm), where peaks and turns in the density of states are
faithfully reproduced from one protein to the next. We also characterize
fluctuations of the spectral density from the average, paving the way to a
meaningful discussion of rare, unusual spectra and the structural reasons for
the deviations in such "outlier" proteins. Since the method used for the
derivation of the vibrational modes (potential energy formulation, set of
degrees of freedom employed, etc.) has a dramatic effect on the spectral
density, another significant implication of our findings is that the
universality can provide an exquisite tool for assessing and improving the
quality of various models used for NMA computations. Finally, we show that the
input configuration too affects the density of modes, thus emphasizing the
importance of simplified potential energy formulations that are minimized at
the outset.


  An ecosystem is a nonlinear dynamical system, its orbits giving rise to the
observed complexity in the system. The diverse components of the ecosystem
interact in discrete time to give rise to emergent features that determine the
trajectory of system's time evolution. The paper studies the evolutionary
dynamics of a toy two species ecosystem modelled as a discrete time Kolmogorov
system. It is assumed that only the two species comprise the ecosystem and
compete with each other for obtaining growth resources, mediated through inter
as well as intraspecific coupling constants to obtain resources for growth.
Numerical simulations reveal the transition from regular to irregular dynamics
and emergence of chaos during the process of evolution of these populations. We
find that the presence or absence of chaotic dynamics is being determined by
the interspecific interaction coefficients. For values of the interspecific
interaction constants widely apart, the system emerges to regular dynamics,
implying coexistence of the competing populations on a long term evolutionary
scenario.


  Synonymous codons, i.e., DNA nucleotide triplets coding for the same amino
acid, are used differently across the variety of living organisms. The
biological meaning of this phenomenon, known as codon usage bias, is still
controversial. In order to shed light on this point, we propose a new codon
bias index, $CompAI$, that is based on the competition between cognate and
near-cognate tRNAs during translation, without being tuned to the usage bias of
highly expressed genes. We perform a genome-wide evaluation of codon bias for
$E.coli$, comparing $CompAI$ with other widely used indices: $tAI$, $CAI$, and
$Nc$. We show that $CompAI$ and $tAI$ capture similar information by being
positively correlated with gene conservation, measured by ERI, and
essentiality, whereas, $CAI$ and $Nc$ appear to be less sensitive to
evolutionary-functional parameters. Notably, the rate of variation of $tAI$ and
$CompAI$ with ERI allows to obtain sets of genes that consistently belong to
specific clusters of orthologous genes (COGs). We also investigate the
correlation of codon bias at the genomic level with the network features of
protein-protein interactions in $E.coli$. We find that the most densely
connected communities of the network share a similar level of codon bias (as
measured by $CompAI$ and $tAI$). Conversely, a small difference in codon bias
between two genes is, statistically, a prerequisite for the corresponding
proteins to interact. Importantly, among all codon bias indices, $CompAI$ turns
out to have the most coherent distribution over the communities of the
interactome, pointing to the significance of competition among cognate and
near-cognate tRNAs for explaining codon usage adaptation.


  The birth of new genes is an important motor of evolutionary innovation.
Whereas many new genes arise by gene duplication, others originate at genomic
regions that do not contain any gene or gene copy. Some of these newly
expressed genes may acquire coding or non-coding functions and be preserved by
natural selection. However, it is yet unclear which is the prevalence and
underlying mechanisms of de novo gene emergence. In order to obtain a
comprehensive view of this process we have performed in-depth sequencing of the
transcriptomes of four mammalian species, human, chimpanzee, macaque and mouse,
and subsequently compared the assembled transcripts and the corresponding
syntenic genomic regions. This has resulted in the identification of over five
thousand new transcriptional multiexonic events in human and/or chimpanzee that
are not observed in the rest of species. By comparative genomics we show that
the expression of these transcripts is associated with the gain of regulatory
motifs upstream of the transcription start site (TSS) and of U1 snRNP sites
downstream of the TSS. We also find that the coding potential of the new genes
is higher than expected by chance, consistent with the presence of
protein-coding genes in the dataset. Using available human tissue proteomics
and ribosome profiling data we identify several de novo genes with translation
evidence. These genes show significant purifying selection signatures,
indicating that they are probably functional. Taken together, the data supports
a model in which frequently-occurring new transcriptional events in the genome
provide the raw material for the evolution of new proteins.


  Non-selective effects, like genetic drift, are an important factor in modern
conceptions of evolution, and have been extensively studied for constant
population sizes. Here, we consider non-selective evolution in the case of
growing populations that are of small size and have varying trait compositions
(e.g. after a population bottleneck). We find that, in these conditions,
populations never fixate to a trait, but tend to a random limit composition,
and that the distribution of compositions 'freezes' to a steady state This
final state is crucially influenced by the initial conditions. We obtain these
findings from a combined theoretical and experimental approach, using multiple
mixed subpopulations of two Pseudomonas putida strains in non-selective growth
conditions as model system. The experimental results for the population
dynamics match the theoretical predictions based on the P\'olya urn model for
all analyzed parameter regimes. In summary, we show that exponential growth
stops genetic drift. This result contrasts with previous theoretical analyses
of non-selective evolution (e.g. genetic drift), which investigated how traits
spread and eventually take over populations (fixate). Moreover, our work
highlights how deeply growth influences non-selective evolution, and how it
plays a key role in maintaining genetic variability. Consequently, it is of
particular importance in life-cycles models of periodically shrinking and
expanding populations.


  In the Monte desert, increasing population density, changing land rights and
infrastructure may encourage livestock activity, with unknown consequences on
ecosystems. Factors that influence livestock settlement distribution may affect
ecosystem degradation. We hypothesize that surface and groundwater availability
influence livestock settlements distribution. We evaluated this hypothesis with
a Monte Carlo based model to simulate Settlement Dynamics in Drylands (SeDD),
which calculates probabilities on a gridded region based on six environmental
factors: groundwater depth, vegetation type, proximity to rivers, paved roads,
old river beds, and existing settlements. A parameter sweep, including millions
of simulations, was run with combinations of parameters related to these
factors. The sets of parameter values that minimized the residuals between
simulations and observations indicated the relative importance of each factors
on settlement distribution. Distances to rivers and old river beds were
critical to explain the current distribution of settlements, while vegetation,
paved roads, and water table depths were not important. Spatial distribution of
simulated vegetation, which included degradation around livestock settlements,
generally agreed with remotely sensed vegetation classes. The model could be a
useful tool to evaluate the effects of land use changes (water provision, river
flows), on settlement distribution and vegetation degradation in arid
environments.


  Reaction-diffusion models are used to describe systems in fields as diverse
as physics, chemistry, ecology and biology. The fundamental quantities in such
models are individual entities such as atoms and molecules, bacteria, cells or
animals, which move and/or react in a stochastic manner. If the number of
entities is large, accounting for each individual is inefficient, and often
partial differential equation (PDE) models are used in which the stochastic
behaviour of individuals is replaced by a description of the averaged, or mean
behaviour of the system. In some situations the number of individuals is large
in certain regions and small in others. In such cases, a stochastic model may
be inefficient in one region, and a PDE model inaccurate in another. To
overcome this problem, we develop a scheme which couples a stochastic
reaction-diffusion system in one part of the domain with its mean field
analogue, i.e. a discretised PDE model, in the other part of the domain. The
interface in between the two domains occupies exactly one lattice site and is
chosen such that the mean field description is still accurate there. This way
errors due to the flux between the domains are small. Our scheme can account
for multiple dynamic interfaces separating multiple stochastic and
deterministic domains, and the coupling between the domains conserves the total
number of particles. The method preserves stochastic features such as
extinction not observable in the mean field description, and is significantly
faster to simulate on a computer than the pure stochastic model.


  The three-compartmental demographic model by Korotaeyv-Malkov-Khaltourina,
connecting population size, economic surplus, and educational level, is
considered from the point of view of dynamical systems theory. It is shown that
there exist two integrals of motion, which enable the system to be reduced to
one non-linear ordinary differential equation. The study of its structure
provides analytical criteria for the dominance ranges of the dynamics of
Malthus and Kremer. Additionally, the particular ranges of parameters enable
the derived general ordinary differential equations to be reduced to the models
of Gompertz and Thoularis-Wallace.


  Structural and dynamical fingerprints of evolutionary optimization in
biological networks are still unclear. We here analyze the dynamics of genetic
regulatory networks responsible for the regulation of cell cycle and cell
differentiation in three organisms or cell types each, and show that they obey
a version of Hebb's rule which we term as coherence. More precisely, we find
that simultaneously expressed genes with a common target are less likely to
conflict at the attractors of the regulatory dynamics. We then investigate the
dependence of coherence on structural parameters, such as the mean number of
inputs per node and the activatory/repressory interaction ratio, as well as on
dynamically determined quantities, such as the basin size and the number of
expressed genes.


  A small portion of a tissue defines a microstate in gene expression space.
Mutations, epigenetic events or external factors cause microstate displacements
which are modeled by combining small independent gene expression variations and
large Levy jumps, resulting from the collective variations of a set of genes.
The risk of cancer in a tissue is estimated as the microstate probability to
transit from the normal to the tumor region in gene expression space. The
formula coming from the contribution of large Levy jumps seems to provide a
qualitatively correct description of the lifetime risk of cancer, and reveals
an interesting connection between the risk and the way the tissue is protected
against infections.


  The fitness landscape defines the relationship between genotypes and fitness
in a given environment, and underlies fundamental quantities such as the
distribution of selection coefficient, or the magnitude and type of epistasis.
A better understanding of variation of landscape structure across species and
environments is thus necessary to understand and predict how populations will
adapt. An increasing number of experiments investigates the properties of
fitness landscapes by identifying mutations, constructing genotypes with
combinations of these mutations, and measuring the fitness of these genotypes.
Yet these empirical landscapes represent a very small sample of the vast space
of all possible genotypes, and this sample is often biased by the protocol used
to identify mutations. Here we develop a rigorous statistical framework based
on Approximate Bayesian Computation to address these concerns, and use this
flexible framework to fit a broad class of phenotypic fitness models (including
Fisher's model) to 26 empirical landscapes representing 9 diverse biological
systems. In spite of uncertainty due to the small size of most published
empirical landscapes, the inferred landscapes have similar structure in similar
biological systems. Surprisingly, goodness of fit tests reveal that this class
of phenotypic models, which has been successful so far in interpreting
experimental data, is a plausible model in only 3 out of 9 biological systems.
More precisely, although Fisher's model was able to explain several statistical
properties of the landscapes - including mean and standard deviation of
selection and epistasis coefficients -, it was often unable to explain the full
structure of fitness landscapes.


  Undergraduate biochemistry laboratory courses often do not provide students
with an authentic research experience, particularly when the express purpose of
the laboratory is purely instructional. However, an instructional laboratory
course that is inquiry- and research-based could simultaneously impart
scientific knowledge and foster a student's research expertise and confidence.
We have developed a year-long undergraduate biochemistry laboratory curriculum
wherein students determine, via experiment and computation, the function of a
protein of known three-dimensional structure. The first half of the course is
inquiry-based and modular in design; students learn general biochemical
techniques while gaining preparation for research experiments in the second
semester. Having learned standard biochemical methods in the first semester,
students independently pursue their own (original) research projects in the
second semester. This new curriculum has yielded an improvement in student
performance and confidence as assessed by various metrics. To disseminate
teaching resources to students and instructors alike, a freely accessible
Biochemistry Laboratory Education resource is available at
http://biochemlab.org.


  What are the structural determinants of protein sequence evolution? A number
of site-specific structural characteristics have been proposed, most of which
are broadly related to either the density of contacts or the solvent
accessibility of individual residues. Most importantly, there has been
disagreement in the literature over the relative importance of solvent
accessibility and local packing density for explaining site-specific sequence
variability in proteins. We show here that this discussion has been confounded
by the definition of local packing density. The most commonly used measures of
local packing, such as the contact number and the weighted contact number,
represent by definition the combined effects of local packing density and
longer-range effects. As an alternative, we here propose a truly local measure
of packing density around a single residue, based on the Voronoi cell volume.
We show that the Voronoi cell volume, when calculated relative to the geometric
center of amino-acid side chains, behaves nearly identically to the relative
solvent accessibility, and both can explain, on average, approximately 34\% of
the site-specific variation in evolutionary rate in a data set of 209 enzymes.
An additional 10\% of variation can be explained by non-local effects that are
captured in the weighted contact number. Consequently, evolutionary variation
at a site is determined by the combined action of the immediate amino-acid
neighbors of that site and of effects mediated by more distant amino acids. We
conclude that instead of contrasting solvent accessibility and local packing
density, future research should emphasize the relative importance of immediate
contacts and longer-range effects on evolutionary variation.


  The 2014-2015 Ebola Virus Disease (EVD) epidemic in West Africa was the
largest ever recorded, representing a fundamental shift in Ebola epidemiology
with unprecedented spatiotemporal complexity. We developed spatial transmission
models using a gravity-model framework to explain spatiotemporal dynamics of
EVD in West Africa at both the national and district-level scales, and to
compare effectiveness of local interventions (e.g. local quarantine) and
long-range interventions (e.g. border-closures). Incorporating spatial
interactions, the gravity model successfully captures the multiple waves of
epidemic growth observed in Guinea. Model simulations indicate that
local-transmission reductions were most effective in Liberia, while long-range
transmission was dominant in Sierra Leone. The model indicates the presence of
spatial herd protection, wherein intervention in one region has a protective
effect on surrounding regions. The district-level intervention analysis
indicates the presence of intervention-amplifying regions, which provide
above-expected levels of reduction in cases and deaths beyond their borders.
The gravity-modeling approach accurately captured the spatial spread patterns
of EVD at both country and district levels, and helps to identify the most
effective locales for intervention. This model structure and intervention
analysis provides information that can be used by public health policymakers to
assist planning and response efforts for future epidemics.


  Each and every biological function in living organism happens as a result of
protein-protein interactions.The diseases are no exception to this. Identifying
one or more proteins for a particular disease and then designing a suitable
chemical compound (known as drug) to destroy these proteins has been an
interesting topic of research in bio-informatics. In previous methods, drugs
were designed using only seven chemical components and were represented as a
fixed-length tree. But in reality, a drug contains many chemical groups
collectively known as pharmacophore. Moreover, the chemical length of the drug
cannot be determined before designing the drug.In the present work, a Particle
Swarm Optimization (PSO) based methodology has been proposed to find out a
suitable drug for a particular disease so that the drug-protein interaction
becomes stable. In the proposed algorithm, the drug is represented as a
variable length tree and essential functional groups are arranged in different
positions of that drug. Finally, the structure of the drug is obtained and its
docking energy is minimized simultaneously. Also, the orientation of chemical
groups in the drug is tested so that it can bind to a particular active site of
a target protein and the drug fits well inside the active site of target
protein. Here, several inter-molecular forces have been considered for accuracy
of the docking energy. Results show that PSO performs better than the earlier
methods.


  We propose a new multi-network-based strategy to integrate different layers
of genomic information and use them in a coordinate way to identify driving
cancer genes. The multi-networks that we consider combine transcription factor
co-targeting, microRNA co-targeting, protein-protein interaction and gene
co-expression networks. The rationale behind this choice is that gene
co-expression and protein-protein interactions require a tight coregulation of
the partners and that such a fine tuned regulation can be obtained only
combining both the transcriptional and post-transcriptional layers of
regulation. To extract the relevant biological information from the
multi-network we studied its partition into communities. To this end we applied
a consensus clustering algorithm based on state of art community detection
methods. Even if our procedure is valid in principle for any pathology in this
work we concentrate on gastric, lung, pancreas and colorectal cancer and
identified from the enrichment analysis of the multi-network communities a set
of candidate driver cancer genes. Some of them were already known oncogenes
while a few are new. The combination of the different layers of information
allowed us to extract from the multi-network indications on the regulatory
pattern and functional role of both the already known and the new candidate
driver genes.


  The Shapley Value and the Fair Proportion Index of phylogenetic trees have
been frequently discussed as prioritization tools in conservation biology. Both
indices rank species according to their contribution to total phylogenetic
diversity, allowing for a simple conservation criterion. While both indices
have their specific advantages and drawbacks, it has recently been shown that
both values are closely related. However, as different authors use different
definitions of the Shapley Value, the specific degree of relatedness depends on
the specific version of the Shapley Value - it ranges from a high correlation
index to equality of the indices. In this note, we first give an overview of
the different indices. Then we turn our attention to the mere ranking order
provided by either of the indices. We compare the rankings obtained from
different versions of the Shapley Value for a phylogenetic tree of European
amphibians and illustrate their differences. We then undertake further analyses
on simulated data and show that even though the chance of two rankings being
exactly identical (when obtained from different versions of the Shapley Value)
decreases with an increasing number of taxa, the distance between the two
rankings converges to zero, i.e., the rankings are becoming more and more
alike. Moreover, we introduce our freely available software package
FairShapley, which was implemented in Perl and with which all calculations have
been performed.


  A distance-based method to reconstruct a phylogenetic tree with $n$ leaves
takes a distance matrix, $n \times n$ symmetric matrix with $0$s in the
diagonal, as its input and reconstructs a tree with $n$ leaves using tools in
combinatorics. A safety radius is a radius from a tree metric (a distance
matrix realizing a true tree) within which the input distance matrices must all
lie in order to satisfy a precise combinatorial condition under which the
distance-based method is guaranteed to return a correct tree. A stochastic
safety radius is a safety radius under which the distance-based method is
guaranteed to return a correct tree within a certain probability. In this paper
we investigated stochastic safety radii for the neighbor-joining (NJ) method
and balanced minimal evolution (BME) method for $n = 5$.


  Cell reproduction involves replication of diverse molecule species, in
contrast to simple replication system with fewer components. Here, we address
why such diversity is sustained despite the efficiency of simple replication
systems, using a cell model with catalytic reaction dynamics that grew by
uptake of environmental resources. Limited resources led to increased diversity
of components within the system, and the number of coexisting species increased
with a negative power of the resource abundances. The diversity was explained
from the optimum growth speed of the cell, determined by a tradeoff between the
utility of diverse resources and the concentration onto fewer components to
increase the reaction rate.


  Dispersal between different habitats influences the dynamics and stability of
populations considerably. Furthermore, these effects depend on the local
interactions of a population with other species. Here, we perform a general and
comprehensive study of the simplest possible system that includes dispersal and
local interactions, namely a 2-patch 2-species system. We evaluate the impact
of dispersal on stability and on the occurrence of bifurcations, including
pattern forming bifurcations that lead to spatial heterogeneity, in 19
different classes of models with the help of the generalized modelling
approach. We find that dispersal often destabilizes equilibria, but it can
stabilize them if it increases population losses. If dispersal is nonrandom,
i.e. if emigration or immigration rates depend on population densities, the
correlation of stability with migration rates is positive in part of the
models. We also find that many systems show all four types of bifurcations and
that antisynchronous oscillations occur mostly with nonrandom dispersal.


  This paper takes a two-pronged approach to investigate the phenomenon of
cross-domain influence on creativity. We present a study in which creative
individuals were asked to list influences on their creative work. More than
half the listed influences were unrelated to their creative domain, thus
demonstrating empirically that cross-domain influence is widespread. We then
present a preliminary model of exaptation, a form of cross-domain influence on
creativity in which a different context suggests a new use for an existing
item, using an example from the study.


  Data obtained from ISSR amplification may readily be extracted but only
allows us to know, for each gene, if a specific allele is present or not. From
this partial information we provide a probabilistic method to reconstruct the
pedigree corresponding to some families of diploid cultivars. This method
consists in determining for each individual what is the most likely couple of
parent pair amongst all older individuals, according to some probability
measure. The construction of this measure bears on the fact that the
probability to observe the specific alleles in the child, given the status of
the parents does not depend on the generation and is the same for each gene.
This assumption is then justified from a convergence result of gene frequencies
which is proved here. Our reconstruction method is applied to a family of 85
living accessions representing the common broom {\it Cytisus scoparius}.


  Comparative and evolutive ecologists are interested in the distribution of
quantitative traits among related species. The classical framework for these
distributions consists of a random process running along the branches of a
phylogenetic tree relating the species. We consider shifts in the process
parameters, which reveal fast adaptation to changes of ecological niches. We
show that models with shifts are not identifiable in general. Constraining the
models to be parsimonious in the number of shifts partially alleviates the
problem but several evolutionary scenarios can still provide the same joint
distribution for the extant species. We provide a recursive algorithm to
enumerate all the equivalent scenarios and to count the effectively different
scenarios. We introduce an incomplete-data framework and develop a maximum
likelihood estimation procedure based on the EM algorithm. Finally, we propose
a model selection procedure, based on the cardinal of effective scenarios, to
estimate the number of shifts and prove an oracle inequality.


  The structure and dynamics of negatively charged nucleic acids strongly
correlate with the concentration and charge of the oppositely charged
counter-ions. It is well known that the structural collapse of DNA is favored
in the presence of additional salt, a source of excess oppositely charged ions.
Under such conditions single stranded DNA adopts a collapsed coil like
conformation, typically characterized by stacking base pairs. Using atomistic
molecular dynamics simulation, we demonstrate that in the presence of
additional divalent salt (MgCl2) single stranded DNA (Dickerson Drew dodecamer)
initially collapses and then expands with increasing salt concentration. This
is due to the overcharging induced DNA chain swelling, a dominant factor at a
higher divalent salt concentration. In a nutshell, our simulations show how in
the presence of divalent salt, non-sequential base stacking and overcharging
competes and affect single stranded DNA dynamics unlike a monovalent salt.


  We consider a model for an epidemic in a population that occupies
geographically distinct locations. The disease is spread within subpopulations
by contacts between infective and susceptible individuals, and is spread
between subpopulations by the migration of infected individuals. We show how
susceptible individuals can act collectively to limit the spread of disease
during the initial phase of an epidemic, by specifying the distribution that
minimises the growth rate of the epidemic when the infectives are migrating so
as to maximise the growth rate. We also give an explicit strategy that
minimises the basic reproduction number, which is also shown be optimal in
terms of the probability of extinction and total size of the epidemic.


  Discovery of microRNAs (miRNAs) relies on predictive models for
characteristic features from miRNA precursors (pre-miRNAs). The short length of
miRNA genes and the lack of pronounced sequence features complicate this task.
To accommodate the peculiarities of plant and animal miRNAs systems, tools for
both systems have evolved differently. However, these tools are biased towards
the species for which they were primarily developed and, consequently, their
predictive performance on data sets from other species of the same kingdom
might be lower. While these biases are intrinsic to the species, the
characterization of their occurrence can lead to computational approaches able
to diminish their negative effect on the accuracy of pre-miRNAs predictive
models. Here, we investigate in this study how 45 predictive models induced for
data sets from 45 species, distributed in eight subphyla, perform when applied
to a species different from the species used in its induction. Our
computational experiments show that the separability of pre-miRNAs and pseudo
pre-miRNAs instances is species-dependent and no feature set performs well for
all species, even within the same subphylum. Mitigating this species
dependency, we show that an ensemble of classifiers reduced the classification
errors for all 45 species. As the ensemble members were obtained using
meaningful, and yet computationally viable feature sets, the ensembles also
have a lower computational cost than individual classifiers that rely on energy
stability parameters, which are of prohibitive computational cost in large
scale applications. In this study, the combination of multiple pre-miRNAs
feature sets and multiple learning biases enhanced the predictive accuracy of
pre-miRNAs classifiers of 45 species. This is certainly a promising approach to
be incorporated in miRNA discovery tools towards more accurate and less
species-dependent tools.


  With different genomes available, unsupervised learning algorithms are
essential in learning genome-wide biological insights. Especially, the
functional characterization of different genomes is essential for us to
understand lives. In this book chapter, we review the state-of-the-art
unsupervised learning algorithms for genome informatics from DNA to MicroRNA.
  DNA (DeoxyriboNucleic Acid) is the basic component of genomes. A significant
fraction of DNA regions (transcription factor binding sites) are bound by
proteins (transcription factors) to regulate gene expression at different
development stages in different tissues. To fully understand genetics, it is
necessary of us to apply unsupervised learning algorithms to learn and infer
those DNA regions. Here we review several unsupervised learning methods for
deciphering the genome-wide patterns of those DNA regions.
  MicroRNA (miRNA), a class of small endogenous non-coding RNA (RiboNucleic
acid) species, regulate gene expression post-transcriptionally by forming
imperfect base-pair with the target sites primarily at the 3$'$ untranslated
regions of the messenger RNAs. Since the 1993 discovery of the first miRNA
\emph{let-7} in worms, a vast amount of studies have been dedicated to
functionally characterizing the functional impacts of miRNA in a network
context to understand complex diseases such as cancer. Here we review several
representative unsupervised learning frameworks on inferring miRNA regulatory
network by exploiting the static sequence-based information pertinent to the
prior knowledge of miRNA targeting and the dynamic information of miRNA
activities implicated by the recently available large data compendia, which
interrogate genome-wide expression profiles of miRNAs and/or mRNAs across
various cell conditions.


  Within-host genetic sequencing from samples collected over time provides a
dynamic view of how viruses evade host immunity. Immune-driven mutations might
stimulate neutralization breadth by selecting antibodies adapted to cycles of
immune escape that generate within-subject epitope diversity. Comprehensive
identification of immune-escape mutations is experimentally and computationally
challenging. With current technology, many more viral sequences can readily be
obtained than can be tested for binding and neutralization, making
down-selection necessary. Typically, this is done manually, by picking variants
that represent different time-points and branches on a phylogenetic tree. Such
strategies are likely to miss many relevant mutations and combinations of
mutations, and to be redundant for other mutations. Longitudinal Antigenic
Sequences and Sites from Intrahost Evolution (LASSIE) uses transmitted-founder
loss to identify virus "hot-spots" under putative immune selection and chooses
sequences that represent recurrent mutations in selected sites. LASSIE favors
earliest sequences in which mutations arise. With well-characterized
longitudinal Env sequences, we confirmed selected sites were concentrated in
antibody contacts and selected sequences represented diverse antigenic
phenotypes. Practical applications include rapidly identifying immune targets
under selective pressure within a subject, selecting minimal sets of reagents
for immunological assays that characterize evolving antibody responses, and for
immunogens in polyvalent "cocktail" vaccines.


  Since genetic algorithm was proposed by John Holland (Holland J. H., 1975) in
the early 1970s, the study of evolutionary algorithm has emerged as a popular
research field (Civicioglu & Besdok, 2013). Researchers from various scientific
and engineering disciplines have been digging into this field, exploring the
unique power of evolutionary algorithms (Hadka & Reed, 2013). Many applications
have been successfully proposed in the past twenty years. For example,
mechanical design (Lampinen & Zelinka, 1999), electromagnetic optimization
(Rahmat-Samii & Michielssen, 1999), environmental protection (Bertini, Felice,
Moretti, & Pizzuti, 2010), finance (Larkin & Ryan, 2010), musical orchestration
(Esling, Carpentier, & Agon, 2010), pipe routing (Furuholmen, Glette, Hovin, &
Torresen, 2010), and nuclear reactor core design (Sacco, Henderson,
Rios-Coelho, Ali, & Pereira, 2009). In particular, its function optimization
capability was highlighted (Goldberg & Richardson, 1987) because of its high
adaptability to different function landscapes, to which we cannot apply
traditional optimization techniques (Wong, Leung, & Wong, 2009). Here we review
the applications of evolutionary algorithms in bioinformatics.


  The dynamics of short-lived mRNA results in bursts of protein production in
gene regulatory networks. We investigate the propagation of bursting noise
between different levels of mathematical modelling, and demonstrate that
conventional approaches based on diffusion approximations can fail to capture
bursting noise. An alternative coarse-grained model, the so-called piecewise
deterministic Markov process (PDMP), is seen to outperform the diffusion
approximation in biologically relevant parameter regimes. We provide a
systematic embedding of the PDMP model into the landscape of existing
approaches, and we present analytical methods to calculate its stationary
distribution and switching frequencies.


  Protein-DNA interactions are crucial for all biological processes. One of the
most important fundamental aspects of these interactions is the process of
protein searching and recognizing specific binding sites on DNA. A large number
of experimental and theoretical investigations have been devoted to uncovering
the molecular description of these phenomena, but many aspects of the
mechanisms of protein search for the targets on DNA remain not well understood.
One of the most intriguing problems is the role of multiple targets in protein
search dynamics. Using a recently developed theoretical framework we analyze
this question in detail. Our method is based on a discrete-state stochastic
approach that takes into account most relevant physical-chemical processes and
leads to fully analytical description of all dynamic properties. Specifically,
systems with two and three targets have been explicitly investigated. It is
found that multiple targets in most cases accelerate the search in comparison
with a single target situation. However, the acceleration is not always
proportional to the number of targets. Surprisingly, there are even situations
when it takes longer to find one of the multiple targets in comparison with the
single target. It depends on the spatial position of the targets, distances
between them, average scanning lengths of protein molecules on DNA, and the
total DNA lengths. Physical-chemical explanations of observed results are
presented. Our predictions are compared with experimental observations as well
as with results from a continuum theory for the protein search. Extensive Monte
Carlo computer simulations fully support our theoretical calculations.


  Characterization of pluripotent states, in which cells can both self-renew
and differentiate, and the irreversible loss of pluripotency are important
research areas in developmental biology. In particular, an understanding of
these processes is essential to the reprogramming of cells for biomedical
applications, i.e., the experimental recovery of pluripotency in differentiated
cells. Based on recent advances in dynamical-systems theory for gene
expression, we propose a gene-regulatory-network model consisting of several
pluripotent and differentiation genes. Our results show that cellular-state
transition to differentiated cell types occurs as the number of cells
increases, beginning with the pluripotent state and oscillatory expression of
pluripotent genes. Cell-cell signaling mediates the differentiation process
with robustness to noise, while epigenetic modifications affecting gene
expression dynamics fix the cellular state. These modifications ensure the
cellular state to be protected against external perturbation, but they also
work as an epigenetic barrier to recovery of pluripotency. We show that
overexpression of several genes leads to the reprogramming of cells, consistent
with the methods for establishing induced pluripotent stem cells. Our model,
which involves the inter-relationship between gene expression dynamics and
epigenetic modifications, improves our basic understanding of cell
differentiation and reprogramming.


  Dengue is one of the most important infectious diseases in the world, in
terms of death and economic cost. Hence, the modeling of dengue is of great
importance to help us understand the dynamics disease, and interfering with its
spreading mathematical by the proposition of control methods. In this work,
control strategies in an attempt to eliminate the Aedes aegypti mosquito, as
well as proposals for the vaccination campaign are evaluated. In our
mathematical model, the mechanical control is accomplished through the
environmental support capacity affected by a discrete function that represents
removal of breeding. Chemical control is carried out from the use of
insecticide and larvicidal. The efficacious of vaccination is studied through
the removal of a fraction of individuals, proportional to the rate of
vaccination, from the susceptible compartment and its transfer to the recovered
compartment. Our major find is that the dengue fever epidemic is only
eradicated with the use of an immunizing vaccine because control measures
directed against the vector are not enough to halt disease spreading. Even
where the infected mosquitoes are eliminated from the system, susceptible
mosquitoes are still present, and infected humans cause dengue fever to
reappear in the human population.


  Viral capsids are structurally constrained by interactions amongst the amino
acids of the constituting proteins. Therefore, epistasis is expected to evolve
amongst sites engaged in physical interactions, and to influence their
substitution rates. In order to study the distribution of structural epistasis,
we modeled \emph{in silico} the capsid of 18 species of the \phix \, family,
including the wild type. \phix \, is amongst the simplest organisms, making it
suitable for experimental evolution and \emph{in silico} modeling. We found
nearly 40 variable amino acid sites in the main capsid protein across the 18
species. To study how epistasis evolved in this group, we reconstructed the
ancestral sequences using a Bayesian phylogenetic framework. The ancestral
states include 8 variable amino acids, for a total of 256 possible haplotypes.
The $dN/dS$ ratio is low, suggesting strong purifying selection, consistent
with the idea that the structure is constrained by some form of stabilizing
selection. For each haplotype in the ancestral node and for the extant species
we estimated \emph{in silico} the distribution of free energies and of
epistasis. We found that free energy has not significantly increased but
epistasis has. We decomposed epistasis up to fifth order and found that
high-order epistasis can sometimes compensate pairwise interactions, often
making the free energy seem additive. By synthesizing some of the ancestral
haplotypes of the capsid gene, we measured their fitness experimentally, and
found that the predicted deviations in the coat protein free energy do not
significantly affect fitness, which is consistent with the stabilizing
selection hypothesis.


  Developing an accurate scoring function is essential for successfully
predicting protein structures. In this study, we developed a statistical
potential function, called OPUS-Beta, for energetically evaluating beta-sheet
contact pattern (the entire residue-residue beta-contacts of a protein)
independent of the atomic coordinate information. The OPUS-Beta potential
contains five terms, i.e., a self-packing term, a pairwise inter-strand packing
term, a pairwise intra-strand packing term, a lattice term and a
hydrogen-bonding term. The results show that, in recognizing the native
beta-contact pattern from decoys, OPUS-Beta potential outperforms the existing
methods in literature, especially in combination with a method using
2D-recursive neural networks (about 5% and 23% improvements in top-1 and top-5
selections). We expect OPUS-Beta potential to be useful in beta-sheet modeling
for proteins.


  The Wright-Fisher Fokker-Planck equation describes the stochastic dynamics of
self-reproducing, competing variants at fixed population size. We use Fisher's
angular transformation, which defines a natural length for this stochastic
process, to remove the co-ordinate dependence of it's diffusive dynamics,
resulting in simple Brownian motion in an unstable potential, driving variants
to extinction or fixation. This insight allows calculation of very accurate
asymptotic formula for the Green's function under neutrality and selection,
using a novel heuristic Gaussian approximation.


  Human immunodeficiency virus (HIV) evolves with extraordinary rapidity.
However, its evolution is constrained by interactions between mutations in its
fitness landscape. Here we show that an Ising model describing these
interactions, inferred from sequence data obtained prior to the use of
antiretroviral drugs, can be used to identify clinically significant sites of
resistance mutations. Successful predictions of the resistance sites indicate
progress in the development of successful models of real viral evolution at the
single residue level, and suggest that our approach may be applied to help
design new therapies that are less prone to failure even where resistance data
is not yet available.


  We study the classical Kimura diffusion operator defined on the n-simplex,
$$L^{Kim}=\sum_{1\leq i,j\leq n+1}x_ix_j\partial_{x_i}\partial_{x_j}$$ We give
novel constructions for the basis of eigenpolynomials, and the solution to the
inhomogeneous Dirichlet problem, which are well adapted to numerical
applications. Our solution of the Dirichlet problem is quite explicit and
provides a precise description of the singularities that arise along the
boundary.


  Introduction The focus of analyzing data from microarray experiments and
extracting biological insight from such data has experienced a shift from
identification of individual genes in association with a phenotype to that of
biological pathways or gene sets. Meanwhile, feature selection algorithm
becomes imperative to cope with the high dimensional nature of many modeling
tasks in bioinformatics. Many feature selection algorithms use information
contained within a gene set as a biological priori, and select relevant
features by incorporating such information. Thus, an integration of gene set
analysis with feature selection is highly desired. Significance analysis of
microarray to gene-set reduction analysis (SAM-GSR) algorithm is a novel
direction of gene set analysis, aiming at further reduction of gene set into a
core subset. Here, we explore the feature selection trait possessed by SAM-GSR
and then modify SAM-GSR specifically to better fulfill this role. Results and
Conclusions Training on a multiple sclerosis (MS) microarray data using both
SAM-GSR and our modification of SAM-GSR, excellent discriminative performance
on an independent test set was achieved. To conclude, absorbing biological
information from a gene set may be helpful for classification and feature
selection. Discussion Given the fact the complete pathway information is far
from completeness, a statistical method capable of constructing biologically
meaningful gene networks is in demand. The basic requirement is that interplay
among genes must be taken into account.


  Conventional analysis of in vitro assays of motor proteins rests on the
assumption that all proteins with the same chemical composition function
identically; however molecule-to-molecule variation is often seen even in
well-controlled experiments. In an effort to obtain a statistically meaningful
set of time traces that simultaneously avoid any experimental artifacts, we
performed quantum-dot labeled kinesin experiments on both surface and levitated
microtubules. Similar to glassy systems, we found that mean velocities of
individual kinesin motors vary widely from one motor to another, the variation
of which is greater than that expected from the stochastic variation of
stepping times. In the presence of heterogeneity, an ensemble-averaged quantity
such as diffusion constant or randomness parameter is ill-defined. We propose
to analyze heterogeneous data from single molecule measurements by decomposing
them into homogeneous subensembles.


  In recent years moment-closure approximations (MA) of the chemical master
equation have become a popular method for the study of stochastic effects in
chemical reaction systems. Several different MA methods have been proposed and
applied in the literature, but it remains unclear how they perform with respect
to each other. In this paper we study the normal, Poisson, log-normal and
central-moment-neglect MAs by applying them to understand the stochastic
properties of chemical systems whose deterministic rate equations show the
properties of bistability, ultrasensitivity and oscillatory behaviour. Our
results suggest that the normal MA is favourable over the other studied MAs. In
particular we found that (i) the size of the region of parameter space where a
closure gives physically meaningful results, e.g. positive mean and variance,
is considerably larger for the normal closure than for the other three
closures; (ii) the accuracy of the predictions of the four closures (relative
to simulations using the stochastic simulation algorithm) is comparable in
those regions of parameter space where all closures give physically meaningful
results; (iii) the Poisson and log-normal MAs are not uniquely defined for
systems involving conservation laws in molecule numbers. We also describe the
new software package MOCA which enables the automated numerical analysis of
various MA methods in a graphical user interface and which was used to perform
the comparative analysis presented in this paper. MOCA allows the user to
develop novel closure methods and can treat polynomial, non-polynomial, as well
as time-dependent propensity functions, thus being applicable to virtually any
chemical reaction system.


  Methods for patterning biomolecules on a substrate at the single molecule
level have been studied as a route to sensors with single-molecular sensitivity
or as a way to probe biological phenomena at the single-molecule level.
However, the arrangement and orientation of single biomolecules on substrates
has been less investigated. Here, we examined the arrangement and orientation
of two rod-like coiled-coil proteins, cortexillin and tropomyosin, around
patterned gold nanostructures. The high aspect ratio of the coiled coils made
it possible to study their orientations and to pursue a strategy of protein
orientation via two-point attachment. The proteins were anchored to the
surfaces using thiol groups, and the number of cysteine residues in tropomyosin
was varied to test how this variation affected the structure and arrangement of
the surface-attached proteins. Molecular dynamics studies were used to
interpret the observed positional distributions. Based on initial studies of
protein attachment to gold post structures, two 31-nm-long tropomyosin
molecules were aligned between the two sidewalls of a trench with a width of 68
nm. Because the approach presented in this study uses one of twenty natural
amino acids, this method provides a convenient way to pattern biomolecules on
substrates using standard chemistry.


  Statistics of Poincar\'e recurrences is studied for the base-pair breathing
dynamics of an all-atom DNA molecule in realistic aqueous environment with
thousands of degrees of freedom. It is found that at least over five decades in
time the decay of recurrences is described by an algebraic law with the
Poincar\'e exponent close to $\beta=1.2$. This value is directly related to the
correlation decay exponent $\nu = \beta -1$, which is close to $\nu\approx
0.15$ observed in the time resolved Stokes shift experiments. By applying the
virial theorem we analyse the chaotic dynamics in polynomial potentials and
demonstrate analytically that exponent $\beta=1.2$ is obtained assuming the
dominance of dipole-dipole interactions in the relevant DNA dynamics. Molecular
dynamics simulations also reveal the presence of strong low frequency noise
with the exponent $\eta=1.6$. We trace parallels with the chaotic dynamics of
symplectic maps with a few degrees of freedom characterized by the Poincar\'e
exponent $\beta \sim 1.5$.


  Reconstructing the tree of life from molecular sequences is a fundamental
problem in computational biology. Modern data sets often contain a large number
of genes, which can complicate the reconstruction problem due to the fact that
different genes may undergo different evolutionary histories. This is the case
in particular in the presence of horizontal genetic transfer (HGT), where a
gene is inherited from a distant species rather than an immediate ancestor.
Such an event produces a gene tree which is distinct from, but related to, the
species phylogeny.
  In previous work, a natural stochastic models of HGT was introduced and
studied. It was shown, both in simulation and theoretical studies, that a
species phylogeny can be reconstructed from gene trees despite surprisingly
high rates of HGT under this model. Rigorous lower and upper bounds on this
achievable rate were also obtained, but a large gap remained. Here we close
this gap, up to a constant. Specifically we show that a species phylogeny can
be reconstructed correctly from gene trees even when, on each gene, each edge
of the species tree has a constant probability of being the location of an HGT
event. Our new reconstruction algorithm, which relies only on unrooted gene
tree topologies, builds the tree recursively from the leaves and runs in
polynomial time.
  We also provide a matching bound in the negative direction (up to a constant)
and extend our results to some cases where gene trees are not perfectly known.


  Reconstructing evolutionary trees from molecular sequence data is a
fundamental problem in computational biology. Stochastic models of sequence
evolution are closely related to spin systems that have been extensively
studied in statistical physics and that connection has led to important
insights on the theoretical properties of phylogenetic reconstruction
algorithms as well as the development of new inference methods. Here, we study
maximum likelihood, a classical statistical technique which is perhaps the most
widely used in phylogenetic practice because of its superior empirical
accuracy.
  At the theoretical level, except for its consistency, that is, the guarantee
of eventual correct reconstruction as the size of the input data grows, much
remains to be understood about the statistical properties of maximum likelihood
in this context. In particular, the best bounds on the sample complexity or
sequence-length requirement of maximum likelihood, that is, the amount of data
required for correct reconstruction, are exponential in the number, $n$, of
tips---far from known lower bounds based on information-theoretic arguments.
Here we close the gap by proving a new upper bound on the sequence-length
requirement of maximum likelihood that matches up to constants the known lower
bound for some standard models of evolution.
  More specifically, for the $r$-state symmetric model of sequence evolution on
a binary phylogeny with bounded edge lengths, we show that the sequence-length
requirement behaves logarithmically in $n$ when the expected amount of mutation
per edge is below what is known as the Kesten-Stigum threshold. In general, the
sequence-length requirement is polynomial in $n$. Our results imply moreover
that the maximum likelihood estimator can be computed efficiently on randomly
generated data provided sequences are as above.


  Geophysical and geochemical anomalies may have a mutagenic effect on plants
growing in active fault zones being the factors of evolutionary transformation
of plant populations. To test this assumption we evaluated the mutation
activity of a Lonicera caerulea natural population in one of the active fault
zones in the Altai Mountains. We derived principal cytogenetic indices (i.e.,
mitotic, prophase, metaphase, anaphase, and telophase indices as well as
proportion and range of abnormal mitoses) for meristematic cells of Lonicera
caerulea seedlings. We found that the local geological and geophysical
environment (i.e., mineralogical composition of rocks and anomalies of the
magnetic field) increases the mitotic activity and the number of abnormal
mitoses in the meristematic cells. The results may help to clarify the role of
environmental conditions of tectonically active regions in microevolutionary
processes.


  Discoveries of endogenous human stem cell-associated retroviruses (SCARs)
revealed consistent activation of specific endogenous retroviral elements in
human preimplantation embryos and documented the essential role of the
sustained retroviral activities in the maintenance of pluripotency, functional
identity and integrity of naive-state embryonic stem cells, and anti-viral
resistance of the early-stage human embryos. Activation of specific SCARs,
namely LTR7.HERVH and LTR5Hs.HERVK, has been demonstrated in patients diagnosed
with multiple types of cancer, autoimmune diseases, neurodegenerative disorders
and it is likely associated with the emergence of clinically lethal therapy
resistant death-from-cancer phenotypes in a sub-set of cancer patients
diagnosed with different types of malignant tumors.


  Signal transduction in biological cells is effected by signaling pathways
that typically include multiple feedback loops. Here we analyze information
transfer through a prototypical signaling module with biochemical feedback. The
module switches stochastically between an inactive and active state; the input
to the module governs the activation rate while the output (i.e., the product
concentration) perturbs the inactivation rate. Using a novel perturbative
approach, we compute the rate with which information about the input is gained
from observation of the output. We obtain an explicit analytical result valid
to first order in feedback strength and to second order in the strength of
input. The total information gained during an extended time interval is found
to depend on the feedback strength only through the total number of
activation/inactivation events.


  The problem of DNA$-$DNA interaction mediated by divalent counterions is
studied using a generalized Grand-canonical Monte-Carlo simulation for a system
of two salts. The effect of the divalent counterion size on the condensation
behavior of the DNA bundle is investigated. Experimentally, it is known that
multivalent counterions have strong effect on the DNA condensation phenomenon.
While tri- and tetra-valent counterions are shown to easily condense free DNA
molecules in solution into toroidal bundles, the situation with divalent
counterions are not as clear cut. Some divalent counterions like Mg$^{+2}$ are
not able to condense free DNA molecules in solution, while some like Mn$^{+2}$
can condense them into disorder bundles. In restricted environment such as in
two dimensional system or inside viral capsid, Mg$^{+2}$ can have strong effect
and able to condense them, but the condensation varies qualitatively with
different system, different coions. It has been suggested that divalent
counterions can induce attraction between DNA molecules but the strength of the
attraction is not strong enough to condense free DNA in solution. However, if
the configuration entropy of DNA is restricted, these attractions are enough to
cause appreciable effects. The variations among different divalent salts might
be due to the hydration effect of the divalent counterions. In this paper, we
try to understand this variation using a very simple parameter, the size of the
divalent counterions. We investigate how divalent counterions with different
sizes can leads to varying qualitative behavior of DNA condensation in
restricted environments. Additionally a Grand canonical Monte-Carlo method for
simulation of systems with two different salts is presented in detail.


  When a system has more than one stable state, how can the stability of these
states be compared? This deceptively simple question has important consequences
for ecosystems, because systems with alternative stable states can undergo
dramatic regime shifts. The probability, frequency, duration, and dynamics of
these shifts will all depend on the relative stability of the stable states.
Unfortunately, the concept of stability in ecology has suffered from
substantial confusion and this is particularly problematic for systems where
stochastic perturbations can cause shifts between coexisting alternative stable
states. A useful way to visualize stable states in stochastic systems is with a
ball-in-cup diagram, in which the state of the system is represented as the
position of a ball rolling on a surface, and the random perturbations can push
the ball from one basin of attraction to another. The surface is determined by
a potential function, which provides a natural stability metric. However,
systems amenable to this representation, called gradient systems, are quite
rare. As a result, the potential function is not widely used and other
approaches based on linear stability analysis have become standard. Linear
stability analysis is designed for local analysis of deterministic systems and,
as we show, can produce a highly misleading picture of how the system will
behave under continual, stochastic perturbations. In this paper, we show how
the potential function can be generalized so that it can be applied broadly,
employing a concept from stochastic analysis called the quasi-potential. Using
three classic ecological models, we demonstrate that the quasi-potential
provides a useful way to quantify stability in stochastic systems.


  The movement of a particle described by Brownian motion is quantified by a
single parameter, $D$, the diffusion constant. The estimation of $D$ from a
discrete sequence of noisy observations is a fundamental problem in biological
single particle tracking experiments since it can report on the environment
and/or the state of the particle itself via hydrodynamic radius. Here we
present a method to estimate $D$ that takes into account several effects that
occur in practice, that are important for correct estimation of $D$, and that
have hitherto not been combined together for estimation of $D$. These effects
are motion blur from finite integration time of the camera, intermittent
trajectories, and time-dependent localization uncertainty. Our estimation
procedure, a maximum likelihood estimation, follows directly from the
likelihood expression for a discretely observed Brownian trajectory that
explicitly includes these effects. The manuscript begins with the formulation
of the likelihood expression and then presents three methods to find the exact
solution. Each method has its own advantages in either computational
robustness, theoretical insight, or the estimation of hidden variables. We then
compare our estimator to previously published estimators using a squared log
loss function to demonstrate the benefit of including these effects.


  Large-scale simultaneous mRNA translation and the resulting competition for
the available ribosomes has important implications to the cell's functioning
and evolution. Developing a better understanding of the intricate correlations
between these simultaneous processes, rather than focusing on the translation
of a single isolated transcript, should help in gaining a better understanding
of mRNA translation regulation and the way elongation rates affect organismal
fitness. A model of simultaneous translation is specifically important when
dealing with highly expressed genes, as these consume more resources. In
addition, such a model can lead to more accurate predictions that are needed in
the interconnection of translational modules in synthetic biology. We develop
and analyze a general model for large-scale simultaneous mRNA translation and
competition for ribosomes. This is based on combining several ribosome flow
models (RFMs) interconnected via a pool of free ribosomes. We prove that the
compound system always converges to a steady-state and that it always entrains
or phase locks to periodically time-varying transition rates in any of the mRNA
molecules. We use this model to explore the interactions between the various
mRNA molecules and ribosomes at steady-state. We show that increasing the
length of an mRNA molecule decreases the production rate of all the mRNAs.
Increasing any of the codon translation rates in a specific mRNA molecule
yields a local effect: an increase in the translation rate of this mRNA, and
also a global effect: the translation rates in the other mRNA molecules all
increase or all decrease. These results suggest that the effect of codon
decoding rates of endogenous and heterologous mRNAs on protein production is
more complicated than previously thought.


  We investigate the dynamics of two models of biological networks with purely
suppressive interactions between the units; species interacting via niche
competition and neurons via inhibitory synaptic coupling. In both of these
cases, power-law scaling of the density of states with probability arises
without any fine-tuning of the model parameters. These results argue against
the increasingly popular notion that non-equilibrium living systems operate at
special critical points, driven by there by evolution so as to enable adaptive
processing of input data.


  The capacity to predict and control bioprocesses is perhaps one of the most
important objectives of biotechnology. Computational simulation is an
established methodology for the design and optimization of bioprocesses, where
the finite elements method (FEM) is at the state-of-art engineering
multi-physics simulation system, with tools such as Finite Element Analysis
(FEA) and Computational Fluid Dynamics (CFD). Although FEA and CFD are
currently applied to bioreactor design, most simulations are restricted to the
multi-physics capabilities of the existing sofware packages. This manuscript is
a contribution for the consolidation of FEM in computational biotechnology, by
presenting a comprehensive review of finite element procedures of the most
common enzymatic mechanisms found in biotechnological processes, such as,
enzyme activation, Michaelis Menten, competitive inhibition, non-competitive
inhibition, anti-competitive inhibition, competition by substrate, sequential
random mechanism, ping-pong bi-bi and Theorel-Chance. Most importantly, the
manuscript opens the possibility for the use of FEM in conjunction with
{\guillemotleft}in-silico{\guillemotright} models of metabolic networks, as
well as, chemical networks in order to simulate complex bioprocesses in
biotechnology, putting emphasis into flux balance analysis, pheno-metabolomics
space exploration in time and space, overcoming the limitations of assuming
chemostat conditions in systems biology computations.


  Phylogenetic trees are used to model evolution: leaves are labelled to
represent contemporary species ("taxa") and interior vertices represent extinct
ancestors. Informally, convex characters are measurements on the contemporary
species in which the subset of species (both contemporary and extinct) that
share a given state, form a connected subtree. Given an unrooted, binary
phylogenetic tree T on a set of n >= 2 taxa, a closed (but fairly opaque)
expression for the number of convex characters on T has been known since 1992,
and this is independent of the exact topology of T. In this note we prove that
this number is actually equal to the (2n-1)th Fibonacci number. Next, we define
g_k(T) to be the number of convex characters on T in which each state appears
on at least k taxa. We show that, somewhat curiously, g_2(T) is also
independent of the topology of T, and is equal to to the (n-1)th Fibonacci
number. As we demonstrate, this topological neutrality subsequently breaks down
for k >= 3. However, we show that for each fixed k >= 1, g_k(T) can be computed
in O(n) time and the set of characters thus counted can be efficiently listed
and sampled. We use these insights to give a simple but effective exact
algorithm for the NP-hard maximum parsimony distance problem that runs in time
$\Theta( \phi^{n} \cdot n^2 )$, where $\phi \approx 1.618...$ is the golden
ratio, and an exact algorithm which computes the tree bisection and
reconnection distance (equivalently, a maximum agreement forest) in time
$\Theta( \phi^{2n}\cdot \text{poly}(n))$, where $\phi^2 \approx 2.619$.


  miRNAs serve as crucial post-transcriptional regulators in various essential
cell fate decision. However, the contribution of the mRNA-miRNA mutual
regulation to bistability is not fully understood. Here, we built a set of
mathematical models of mRNA-miRNA interactions and systematically analyzed the
sensitivity of response curves under various conditions. First, we found that
mRNA-miRNA reciprocal regulation could manifest ultrasensitivity to subserve
the generation of bistability when equipped with a positive feedback loop.
Second, the region of bistability is expanded by a stronger competing mRNA
(ceRNA). Interesting, bistability can be emerged without feedback loop if
multiple miRNA binding sites exist on a target mRNA. Thus, we demonstrated the
importance of simple mRNA-miRNA reciprocal regulation in cell fate decision.


  The following manuscript reviews various theories of bipedalism and provides
a holistic answer to human evolution. There are two questions regarding
bipedalism: i) why were the earliest hominins partially bipedal? and ii) why
did hominins become increasingly bipedal over time and replace their less
bipedal ancestors? To answer these questions, the prominent theories in the
field, such as the savanna-based theory, the postural feeding hypotheses, and
the provisioning model, are collectively examined. Because biological evolution
is an example of trial and error and not a simple causation, there may be
multiple answers to the evolution of bipedalism. The postural feeding
hypothesis (reaching for food/balancing) provides an explanation for the
partial bipedalism of the earliest hominins. The savannah-based theory
describes how the largely bipedal hominins that started to settle on the ground
became increasingly bipedal. The provisioning model (food-gathering/monogamy)
explains questions arising after the postural feeding hypothesis and before the
savannah theory in an evolutionary timeline. Indeed, there are no straight
lines between the theories, and multiple forces could have pushed the evolution
of bipedalism at different points. Finally, this manuscript states that the
arboreal hominins that possessed ambiguous traits of bipedalism were eliminated
through choice and selection. Using the biological analogy of the okapi and
giraffe, I explain how one of the branches (Homo) became increasingly bipedal
while the other (Pan) adapted to locomotion for forest life by narrowing the
anatomical/biological focus in evolution.


  During cell migration, cells become polarized, change their shape, and move
in response to various internal and external cues. Cell polarization is defined
through the spatio-temporal organization of molecules such as PI3K or small
GTPases, and is determined by intracellular signaling networks. It results in
directional forces through actin polymerization and myosin contractions. Many
existing mathematical models of cell polarization are formulated in terms of
reaction-diffusion systems of interacting molecules, and are often defined in
one or two spatial dimensions. In this paper, we introduce a 3D
reaction-diffusion model of interacting molecules in a single cell, and find
that cell geometry has an important role affecting the capability of a cell to
polarize, or change polarization when an external signal changes direction. Our
results suggest a geometrical argument why more roundish cells can repolarize
more effectively than cells which are elongated along the direction of the
original stimulus, and thus enable roundish cells to turn faster, as has been
observed in experiments. On the other hand, elongated cells preferentially
polarize along their main axis even when a gradient stimulus appears from
another direction. Furthermore, our 3D model can accurately capture the effect
of binding and unbinding of important regulators of cell polarization to and
from the cell membrane. This spatial separation of membrane and cytosol, not
possible to capture in 1D or 2D models, leads to marked differences of our
model from comparable lower-dimensional models.


  A protein residue network (PRN) is a network of interacting amino acids
within a protein. We describe characteristics of a sparser, highly central and
more volatile sub-network of a PRN called the short-cut network (SCN), as a
protein folds under molecular dynamics (MD) simulation with the goal of
understanding how proteins form navigable small-world networks within
themselves. The edges of an SCN are found via a local greedy search on a PRN.
SCNs grow in size and transitivity strength as a protein folds, and SCNs from
successful MD trajectories are better formed in these terms. Findings from an
investigation on how to model the formation of SCNs using dynamic graph theory,
and suggestions to move forward are presented. A SCN is enriched with
short-range contacts and its formation correlates positively with secondary
structure formation. Thus our approach to modeling PRN formation, in essence
protein folding from a graph theoretic view point, is more in tune with the
notion of increasing order to a random graph than the other way around, and
this increase in order coincides with improved navigability of PRNs.


  Molecular dynamics (MD) simulations of ions (K$^+$, Na$^+$, Ca$^{2+}$ and
Cl$^-$) in aqueous solutions are investigated. Water is described using the
SPC/E model. A stochastic coarse-grained description for ion behaviour is
presented and parameterized using MD simulations. It is given as a system of
coupled stochastic and ordinary differential equations, describing the ion
position, velocity and acceleration. The stochastic coarse-grained model
provides an intermediate description between all-atom MD simulations and
Brownian dynamics (BD) models. It is used to develop a multiscale method which
uses all-atom MD simulations in parts of the computational domain and (less
detailed) BD simulations in the remainder of the domain.


  We analyse the effect of intrinsic fluctuations on the properties of bistable
stochastic systems with time scale separation operating under1 quasi-steady
state conditions. We first formulate a stochastic generalisation of the
quasi-steady state approximation based on the semi-classical approximation of
the partial differential equation for the generating function associated with
the Chemical Master Equation. Such approximation proceeds by optimising an
action functional whose associated set of Euler-Lagrange (Hamilton) equations
provide the most likely fluctuation path. We show that, under appropriate
conditions granting time scale separation, the Hamiltonian can be re-scaled so
that the set of Hamilton equations splits up into slow and fast variables,
whereby the quasi-steady state approximation can be applied. We analyse two
particular examples of systems whose mean-field limit has been shown to exhibit
bi-stability: an enzyme-catalysed system of two mutually-inhibitory proteins
and a gene regulatory circuit with self-activation. Our theory establishes that
the number of molecules of the conserved species are order parameters whose
variation regulates bistable behaviour in the associated systems beyond the
predictions of the mean-field theory. This prediction is fully confirmed by
direct numerical simulations using the stochastic simulation algorithm. This
result allows us to propose strategies whereby, by varying the number of
molecules of the three conserved chemical species, cell properties associated
to bistable behaviour (phenotype, cell-cycle status, etc.) can be controlled.


  Cubical complexes are metric spaces constructed by gluing together unit cubes
in an analogous way to the construction of simplicial complexes. We construct
Brownian motion on such spaces, define random walks, and prove that the
transition kernels of the random walks converge to that for Brownian motion.
The proof involves pulling back onto the complex the distribution of Brownian
sample paths on the standard cube, and combining this with a distribution on
walks between cubes in the complex. The main application lies in analysing sets
of evolutionary trees: several tree spaces are cubical complexes and we briefly
describe our results and some applications in this context. Our results extend
readily to a class of polyhedral complex in which every cell of maximal
dimension is isometric to a given fixed polyhedron.


  The eukaryotic Cys2His2 zinc finger proteins bind to DNA ubiquitously at
highly conserved domains, responsible for gene regulation and the spatial
organization of DNA. To study and understand the zinc finger DNA-protein
interaction, we use the extended ladder in the DNA model proposed by Zhu,
Rasmussen, Balatsky \& Bishop (2007) \cite{Zhu-2007}. Considering one single
spinless electron in each nucleotide $\pi$-orbital along a double DNA chain
(dDNA), we find a typical pattern for the {\color{black} bottom of the occupied
molecular orbital (BOMO)}, highest occupied molecular orbital (HOMO) and lowest
unoccupied orbital (LUMO) along the binding sites. We specifically looked at
two members of zinc finger protein family: specificity protein 1 (SP1) and
early grown response 1 transcription factors (EGR1). When the valence band is
filled, we find electrons in the purines along the nucleotide sequence,
compatible with the electric charges of the binding amino acids in SP1 and EGR1
zinc finger.


  At the core of high throughput DNA sequencing platforms lies a bio-physical
surface process that results in a random geometry of clusters of homogenous
short DNA fragments typically hundreds of base pairs long - bridge
amplification. The statistical properties of this random process and length of
the fragments are critical as they affect the information that can be
subsequently extracted, i.e., density of successfully inferred DNA fragment
reads. The ensemble of overlapping DNA fragment reads are then used to
computationally reconstruct the much longer target genome sequence, e.g,
ranging from hundreds of thousands to billions of base pairs. The success of
the reconstruction in turn depends on having a sufficiently large ensemble of
DNA fragments that are sufficiently long. In this paper using stochastic
geometry we model and optimize the end-to-end process linking and partially
controlling the statistics of the physical processes to the success of the
computational step. This provides, for the first time, a framework capturing
salient features of such sequencing platforms that can be used to study cost,
performance or sensitivity of the sequencing process.


  We present a theoretical framework to analyze the dynamics of gene expression
with stochastic bursts. Beginning with an individual-based model which fully
accounts for the messenger RNA (mRNA) and protein populations, we propose a
novel expansion of the master equation for the joint process. The resulting
coarse-grained model reduces the dimensionality of the system, describing only
the protein population while fully accounting for the effects of discrete and
fluctuating mRNA population. Closed form expressions for the stationary
distribution of the protein population and mean first-passage times of the
coarse-grained model are derived and large-scale Monte Carlo simulations show
that the analysis accurately describes the individual-based process accounting
for mRNA population, in contrast to the failure of commonly proposed
diffusion-type models.


  Although proteins have been recognized as small-world networks and their
small-world network properties of clustering and short paths have been
exploited computationally to produce biologically relevant information, they
have not been truly explored as such, i.e. as navigable small-world networks in
the original spirit of Milgram's work. This research seeks to fill this gap by
exploring local search on a network representation of proteins and to probe the
source of navigability in proteins. Previously, we confirmed that proteins are
navigable small-world networks and observed that local search paths exhibit
different characteristics from global search paths. In this paper, we
investigate the biological relevance of the differences in path characteristics
on a type III receptor tyrosine kinase (KIT). A chief difference that works in
favour of local search paths as intra-protein communication pathways is their
weaker proclivity, compared to global search paths, for traversing long-range
edges. Long-range edges tend to be less stable and their inclusion tends to
decrease the communication propensity of a path. The source of protein
navigability is traced to clustering provided by short-range edges. The
majority of a protein's short-range edges reside within structures deemed
important for long-range energy transport and modulation of allosteric
communication in proteins. Therefore, the disruption of intra-protein
communication as a result of the destruction of these structures via random
rewiring is expected. A local search perspective leads us to this expected
conclusion while a global search perspective does not. These findings initiate
the compilation of a list of path properties that are characteristic of
intra-protein pathways and could suggest fresh avenues for evolving and
regulating navigable (small-world) networks.


  Boolean networks are an important class of computational models for molecular
interaction networks. Boolean canalization, a type of hierarchical clustering
of the inputs of a Boolean function, has been extensively studied in the
context of network modeling where each layer of canalization adds a degree of
stability in the dynamics of the network. Recently, dynamic network control
approaches have been used for the design of new therapeutic interventions and
for other applications such as stem cell reprogramming. This work studies the
role of canalization in the control of Boolean molecular networks. It provides
a method for identifying the potential edges to control in the wiring diagram
of a network for avoiding undesirable state transitions. The method is based on
identifying appropriate input-output combinations on undesirable transitions
that can be modified using the edges in the wiring diagram of the network.
Moreover, a method for estimating the number of changed transitions in the
state space of the system as a result of an edge deletion in the wiring diagram
is presented. The control methods of this paper were applied to a mutated
cell-cycle model and to a p53-mdm2 model to identify potential control targets.


  We employ a recently developed coarse-grained model for peptides and proteins
where the effect of pH is automatically included. We explore the effect of pH
in the aggregation process of the amyloidogenic peptide KTVIIE and two related
sequences, using three different pH environments. Simulations using large
systems (24 peptides chains per box) allow us to correctly account for the
formation of realistic peptide aggregates. We evaluate the thermodynamic and
kinetic implications of changes in sequence and pH upon peptide aggregation,
and we discuss how a minimalistic coarse-grained model can account for these
details.


  The Kets, an ethnic group in the Yenisei River basin, Russia, are considered
the last nomadic hunter-gatherers of Siberia, and Ket language has no
transparent affiliation with any language family. We investigated connections
between the Kets and Siberian and North American populations, with emphasis on
the Mal'ta and Paleo-Eskimo ancient genomes using original data from 46
unrelated samples of Kets and 42 samples of their neighboring ethnic groups
(Uralic-speaking Nganasans, Enets, and Selkups). We genotyped over 130,000
autosomal SNPs, determined mitochondrial and Y-chromosomal haplogroups, and
performed high-coverage genome sequencing of two Ket individuals. We
established that the Kets belong to the cluster of Siberian populations related
to Paleo-Eskimos. Unlike other members of this cluster (Nganasans, Ulchi,
Yukaghirs, and Evens), Kets and closely related Selkups have a high degree of
Mal'ta ancestry. Implications of these findings for the linguistic hypothesis
uniting Ket and Na-Dene languages into a language macrofamily are discussed.


  In population and evolutionary biology, hypotheses about micro-evolutionary
and macro-evolutionary processes are commonly tested by comparing the shape
indices of empirical evolutionary trees with those predicted by neutral models.
A key ingredient in this approach is the ability to compute and quantify
distributions of various tree shape indices under random models of interest. As
a step to meet this challenge, in this paper we investigate the joint
distribution of cherries and pitchforks (that is, subtrees with two and three
leaves) under two widely used null models: the Yule-Harding-Kingman (YHK) model
and the proportional to distinguishable arrangements (PDA) model. Based on two
novel recursive formulae, we propose a dynamic approach to numerically compute
the exact joint distribution (and hence the marginal distributions) for trees
of any size. We also obtained insights into the statistical properties of trees
generated under these two models, including a constant correlation between the
cherry and the pitchfork distributions under the YHK model, the log-concavity
and unimodality of cherry distributions under both models. In particular, we
show the existence of a unique change point for cherry distribution between the
two models, that is, there exists a critical value $\tau_n$ for each $n\geq 4$
such that the probability that a random tree with $n$ leaves generated under
the YHK model contains $k$ cherries is lower than that under the PDA model if
$1<k< \tau_n$, and higher if $\tau_n<k\le n/2$.


  Enzymes have been recently found to exhibit enhanced diffusion due to their
catalytic activities. A recent experiment [C. Riedel et al., Nature 517, 227
(2015)] has found evidence that suggests this phenomenon might be controlled by
the degree of exothermicity of the catalytic reaction involved. Four mechanisms
that can lead to this effect, namely, self-thermophoresis, boost in kinetic
energy, stochastic swimming, and collective heating, are critically discussed,
and it is shown that only the last two could be strong enough to account for
the observations. The resulting quantitative description is used to examine the
biological significance of the effect.


  Inspired by the fact that people have diverse propensities to punish
wrongdoers, we study a spatial public goods game with defectors and different
types of punishing cooperators. During the game, cooperators punish defectors
with class-specific probabilities and subsequently share the associated costs
of sanctioning. We show that in the presence of different punishing cooperators
the highest level of public cooperation is always attainable through a
selection mechanism. Interestingly, the selection not necessarily favors the
evolution of punishers who would be able to prevail on their own against the
defectors, nor does it always hinder the evolution of punishers who would be
unable to prevail on their own. Instead, the evolutionary success of punishing
strategies depends sensitively on their invasion velocities, which in turn
reveals fascinating examples of both competition and cooperation among them.
Furthermore, we show that under favorable conditions, when punishment is not
strictly necessary for the maintenance of public cooperation, the less
aggressive, mild form of sanctioning is the sole victor of selection process.
Our work reveals that natural strategy selection can not only promote, but
sometimes also hinder competition among prosocial strategies.


  Advances in de novo synthesis of DNA and computational gene design methods
make possible the customization of genes by direct manipulation of features
such as codon bias and mRNA secondary structure. Codon context is another
feature significantly affecting mRNA translational efficiency, but existing
methods and tools for evaluating and designing novel optimized protein coding
sequences utilize untested heuristics and do not provide quantifiable
guarantees on design quality. In this study we examine statistical properties
of codon context measures in an effort to better understand the phenomenon. We
analyze the computational complexity of codon context optimization and design
exact and efficient heuristic gene recoding algorithms under reasonable
constraint models. We also present a web-based tool for evaluating codon
context bias in the appropriate context.


  The properties of a front between two different phases in the presence of a
smoothly inhomogeneous external field that takes its critical value at the
crossing point is analyzed. Two generic scenarios are studied. In the first,
the system admits a bistable solution and the external field governs the rate
in which one phase invades the other. The second mechanism corresponds to a
second order transition that, in the case of reactive systems, takes the form
of a transcritical bifurcation at the crossing point. We solve for the front
shape and its response to external white noise, showing that static properties
and also some of the dynamics features cannot distinguish between the two
scenarios. The only reliable indicator turns out to be the fluctuation
statistics. These take a Gaussian form in the bifurcation case and a
double-peak shape in a bistable system. The results of a recent analysis of the
morphogenesis process in Drosophila embryos are reanalyzed and we show, in
contrast with the interpretation suggested by Krotov et. al., that the
plausible underlying dynamics is bistable and not bifurcational.


  Quantitative Fitness Analysis (QFA) is a high-throughput experimental and
computational methodology for measuring the growth of microbial populations.
QFA screens can be used to compare the health of cell populations with and
without a mutation in a query gene in order to infer genetic interaction
strengths genome-wide, examining thousands of separate genotypes. We introduce
Bayesian, hierarchical models of population growth rates and genetic
interactions that better reflect QFA experimental design than current
approaches. Our new approach models population dynamics and genetic interaction
simultaneously, thereby avoiding passing information between models via a
univariate fitness summary. Matching experimental structure more closely,
Bayesian hierarchical approaches use data more efficiently and find new
evidence for genes which interact with yeast telomeres within a published
dataset.


  This study investigates the influence of lattice structure in evolutionary
games. The snowdrift games is considered in networks with high clustering
coefficients, that use four different strategy-updating. Analytical conjectures
using pair approximation were compared with the numerical results. Results
indicate that general statements asserting that the lattice structure enhances
cooperation are misleading.


  Exaggerate radical-induced DNA damage under magnetic fields is of great
concerns to medical biosafety and to bio-molecular device based upon DNA
electronic conductivity. In this report, the effect of applying an external
magnetic field (MF) on DNA-mediated charge transport (CT) was investigated by
studying guanine oxidation by a kinetics trap (8CPG) via photoirradiation of
anthraquinone (AQ) in the presence of an external MF. Positive enhancement in
CT efficiencies was observed in both the proximal and distal 8CPG after
applying a static MF of 300 mT. MF assisted CT has shown sensitivities to
magnetic field strength, duplex structures, and the integrity of base pair
stacking. MF effects on spin evolution of charge injection upon AQ irradiation
and alignment of base pairs to CT-active conformation during radical
propagation were proposed to be the two major factors that MF attributed to
facilitate DNA-mediated CT. Herein, our results suggested that the electronic
conductivity of duplex DNA can be enhanced by applying an external MF. MF
effects on DNA-mediated CT may offer a new avenue for designing DNA-based
electronic device, and unraveled MF effects on redox and radical relevant
biological processes.


  Background. Dynamical models of gene regulatory networks (GRNs) are highly
effective in describing complex biological phenomena and processes, such as
cell differentiation and cancer development. Yet, the topological and
functional characterization of real GRNs is often still partial and an
exhaustive picture of their functioning is missing.
  Motivation. We here introduce CABeRNET, a Cytoscape app for the generation,
simulation and analysis of Boolean models of GRNs, specifically focused on
their augmentation when a only partial topological and functional
characterization of the network is available. By generating large ensembles of
networks in which user-defined entities and relations are added to the original
core, CABeRNET allows to formulate hypotheses on the missing portions of real
networks, as well to investigate their generic properties, in the spirit of
complexity science.
  Results. CABeRNET offers a series of innovative simulation and modeling
functions and tools, including (but not being limited to) the dynamical
characterization of the gene activation patterns ruling cell types and
differentiation fates, and sophisticated robustness assessments, as in the case
of gene knockouts. The integration within the widely used Cytoscape framework
for the visualization and analysis of biological networks, makes CABeRNET a new
essential instrument for both the bioinformatician and the computational
biologist, as well as a computational support for the experimentalist. An
example application concerning the analysis of an augmented T-helper cell GRN
is provided.


  We propose a working hypothesis supported by numerical simulations that brain
networks evolve based on the principle of the maximization of their internal
information flow capacity. We find that synchronous behavior and capacity of
information flow of the evolved networks reproduce well the same behaviors
observed in the brain dynamical networks of Caenorhabditis elegans and humans,
networks of Hindmarsh-Rose neurons with graphs given by these brain networks.
We make a strong case to verify our hypothesis by showing that the neural
networks with the closest graph distance to the brain networks of
Caenorhabditis elegans and humans are the Hindmarsh-Rose neural networks
evolved with coupling strengths that maximize information flow capacity.
Surprisingly, we find that global neural synchronization levels decrease during
brain evolution, reflecting on an underlying global no Hebbian-like evolution
process, which is driven by no Hebbian-like learning behaviors for some of the
clusters during evolution, and Hebbian-like learning rules for clusters where
neurons increase their synchronization.


  Gene Regulatory Networks (GRNs) consisting of combinations of transcription
factors (TFs) and their cis promoters are assumed to be sufficient to direct
the development of organisms. Mutations in GRNs are assumed to be the primary
drivers for the evolution of multicellular life. Here it is proven that neither
of these assumptions is correct. They are inconsistent with fundamental
principles of combinatorics of bounded encoded networks. It is shown there are
inherent complexity and control capacity limits for any gene regulatory network
that is based solely on protein coding genes such as transcription factors.
This result has significant practical consequences for understanding
development, evolution, the Cambrian Explosion, as well as multi-cellular
diseases such as cancer. If the arguments are sound, then genes cannot explain
the development of complex multicellular organisms and genes cannot explain the
evolution of complex multicellular life.


  When relaxation towards an equilibrium or steady state is exponential at
large times, one usually considers that the associated relaxation time $\tau$,
i.e., the inverse of that decay rate, is the longest characteristic time in the
system. However that need not be true, and in particular other times such as
the lifetime of an infinitesimal perturbation can be much longer. In the
present work we demonstrate that this paradoxical property can arise even in
quite simple systems such as a chain of reactions obeying mass action kinetics.
By mathematical analysis of simple reaction networks, we pin-point the reason
why the standard relaxation time does not provide relevant information on the
potentially long transient times of typical infinitesimal perturbations.
Overall, we consider four characteristic times and study their behavior in both
simple chains and in more complex reaction networks taken from the publicly
available database "Biomodels." In all these systems involving mass action
rates, Michaelis-Menten reversible kinetics, or phenomenological laws for
reaction rates, we find that the characteristic times corresponding to
lifetimes of tracers and of concentration perturbations can be much longer than
$\tau$.


  In differential expression (DE) analysis of RNA-seq count data, it is known
that genes with a larger read number are more likely to be differentially
expressed. This bias has a profound effect on the subsequent Gene Ontology (GO)
analysis by perturbing the ranks of gene-sets. Another known bias is that the
commonly used parametric DE analysis methods (e.g., edgeR, DESeq and baySeq)
tend to yield more DE genes as the sequencing depth is increased. We
nevertheless show that these biases are in fact confined to data of the
technical replicate type. We also show the GO or gene-set enrichment analysis
methods applied to technical replicate data result in considerable number of
false positives. In conclusion, the current DE and enrichment analysis methods
can be confidently used for biological replicate count data, while caution
should be exercised when analysing technical replicate data.


  We recognise that Stegosaurus exhibited exterior chirality and could,
therefore, have assumed either of two distinct, mirror-image forms. Our
preliminary investigations suggest that both existed. Stegosaurus's exterior
chirality raises new questions such as the validity of well-known exhibits
whilst offering new insights into long-standing questions such as the function
of the plates. We inform our discussions throughout with examples of modern-day
animals that exhibit exterior chirality.