  Background Despite its status as a model organism, the development of
Caenorhabditis elegans is not necessarily archetypical for nematodes. The
phylum Nematoda is divided into the Chromadorea (indcludes C. elegans) and the
Enoplea. Compared to C. elegans, enoplean nematodes have very different
patterns of cell division and determination. Embryogenesis of the enoplean
Romanomermis culicivorax has been studied in great detail, but the genetic
circuitry underpinning development in this species is unknown. Results We
created a draft genome of R. culicivorax and compared its developmental gene
content with that of two nematodes, C. elegans and Trichinella spiralis
(another enoplean), and a representative arthropod Tribolium castaneum. This
genome evidence shows that R. culicivorax retains components of the conserved
metazoan developmental toolkit lost in C. elegans. T. spiralis has
independently lost even more of the toolkit than has C. elegans. However, the
C. elegans toolkit is not simply depauperate, as many genes essential for
embryogenesis in C. elegans are unique to this lineage, or have only extremely
divergent homologues in R. culicivorax and T. spiralis. These data imply
fundamental differences in the genetic programmes for early cell specification,
inductive interactions, vulva formation and sex determination. Conclusions Thus
nematodes, despite their apparent phylum-wide morphological conservatism, have
evolved major differences in the molecular logic of their development. R.
culicivorax serves as a tractable, contrasting model to C. elegans for
understanding how divergent genomic and thus regulatory backgrounds can
generate a conserved phenotype. The availability of the draft genome will
promote use of R. culicivorax as a research model.


  Age structure is incorporated in many types of epidemic model. Often it is
convenient to assume that such models converge to early asymptotic behaviour
quickly, before the susceptible population has been appreciably depleted. We
make use of dynamical systems theory to show that for some reasonable parameter
values, this convergence can be slow. Such a possibility should therefore be
considered when parameterising age-structured epidemic models.


  High-throughput transcriptomics experiments have reached the stage where the
count of the number of reads alignable to a given position can be treated as an
almost-continuous signal. This allows to ask questions of
biophysical/biotechnical nature, but which may still have biological
implications. Here we show that when sequencing RNA fragments from one end, as
it is the case on most platforms, an oscillation in the read count is observed
at the other end. We further show that these oscillations can be well described
by Kolmogorov's 1941 broken stick model. We investigate how the model can be
used to improve predictions of gene ends (3' transcript ends) but conclude that
with present data the improvement is only marginal. The results highlight
subtle effects in high-throughput transcriptomics experiments which do not have
a biological origin, but which may still be used to obtain biological
information.


  An evolutionary tree is a cascade of bifurcations starting from a single
common root, generating a growing set of daughter species as time goes by.
Species here is a general denomination for biological species, spoken languages
or any other entity evolving through heredity. From the N currently alive
species within a clade, distances are measured through pairwise comparisons
made by geneticists, linguists, etc. The larger is such a distance for a pair
of species, the older is their last common ancestor. The aim is to reconstruct
the past unknown bifurcations, i.e. the whole clade, from the knowledge of the
N(N-1)/2 quoted distances taken for granted. A mechanical method is presented,
and its applicability discussed.


  Dryland ecosystems commonly exhibit periodic bands of vegetation, thought to
form due to competition between individual plants for heterogeneously
distributed water. In this paper, we develop a Fourier method for locally
identifying the pattern wavenumber and orientation, and apply it to aerial
images from a region of vegetation patterning near Fort Stockton, Texas. We
find that the local pattern wavelength and orientation are typically coherent,
but exhibit both rapid and gradual variation driven by changes in hillslope
gradient and orientation, the potential for water accumulation, or soil type.
Endogenous pattern dynamics, when simulated for spatially homogeneous
topographic and vegetation conditions, predict pattern properties that are much
less variable than the orientation and wavelength observed in natural systems.
Our local pattern analysis, combined with ancillary datasets describing soil
and topographic variation, highlights a largely unexplored correlation between
soil depth, pattern coherence, vegetation cover and pattern wavelength. It
also, surprisingly, suggests that downslope accumulation of water may play a
role in changing vegetation pattern properties.


  Hemagglutinin (HA) and neuraminidase (NA) are highly variable envelope
glycoproteins. Here hydropathic analysis, previously applied to quantify common
flu (H1N1) evolution (1934-), is applied to the evolution of less common but
more virulent (avian derived) H3N2 (1968-), beginning with N2. Whereas N1
exhibited opposing migration and vaccination pressures, the dominant N2 trend
is due to vaccination, with only secondary migration interactions. Separation
and evaluation of these effects is made possible by the use of two distinct
hydropathic scales representing first-order and second-order thermodynamic
interactions. The evolutions of H1 and H3 are more complex, with larger
competing migration and vaccination effects. The linkages of H3 and N2
evolutionary trends are examined on two modular length scales, medium
(glycosidic) and large (corresponding to sialic acid interactions). The
hierarchical hydropathic results complement and greatly extend advanced
phylogenetic results obtained from similarity studies. They exhibit simple
quantitative trends that can be transferred to engineer oncolytic properties of
other viral proteins to treat recalcitrant cancers.


  Genetic interactions pervade every aspect of biology, from evolutionary
theory where they determine the accessibility of evolutionary paths, to
medicine where they contribute to complex genetic diseases. Until very
recently, studies on epistatic interactions have been based on a handful of
mutations, providing at best anecdotal evidence about the frequency and the
typical strength of genetic interactions. In this study we analyze the publicly
available Data Repository of Yeast Genetic INteractions (DRYGIN), which
contains the growth rates of over five million double gene knockout mutants. We
discuss a geometric definition of epistasis which reveals a simple and
surprisingly weak scaling law for the characteristic strength of genetic
interactions as a function of the effects of the mutations being combined. We
then utilize this scaling to quantify the roughness of naturally occurring
fitness landscapes. Finally, we show how the observed roughness differs from
what is predicted by Fisher's geometric model of epistasis and discuss its
consequences on the evolutionary dynamics. Although epistatic interactions
between specific genes remain largely unpredictable, the statistical properties
of an ensemble of interactions can display conspicuous regularities and be
described by simple mathematical laws. By exploiting the amount of data
produced by modern high-throughput techniques it is now possible to thoroughly
test the predictions of theoretical models of genetic interactions and to build
informed computational models of evolution on realistic fitness landscapes.


  We study the epidemic spreading process following contact dynamics with
heavy-tailed waiting time distributions. We show both analytically and
numerically that the temporal heterogeneity of contact dynamics can
significantly suppress the disease's transmissibility, hence the size of
epidemic outbreak, obstructing the spreading process. Furthermore, when the
temporal heterogeneity is strong enough, one obtains the vanishing
transmissibility for any finite recovery time and regardless of underlying
structure of contacts, the condition of which was derived.


  Biological fitness is not an observable quantity and must be inferred from
population dynamics. Bayesian inference applied to the Moran process and
variants yields a robust inference method that can infer fitness in populations
evolving via a Moran dynamic and generalizations. Information about fitness is
derived solely from birth-events in birth-death and death-birth processes in
which selection acts proportionally to fitness, which allows the method to be
applied to populations on a network where the network itself may be changing in
time. Populations may also be allowed to change size while still allowing
estimates for fitness to be inferred.


  We have studied the mobility of the folding catalyst, protein
disulphide-isomerase (PDI), by molecular dynamics and by a rapid approach based
on flexibility. We analysed our simulations using measures of backbone
movement, relative positions and orientations of domains, and distances between
functional sites. Despite their different assumptions, the two methods are
surprisingly consistent. Both methods agree that motion of domains is dominated
by hinge and rotation motion of the a and a' domains relative to the central
b-b' domain core. We identify the a' domain as showing the greatest
intra-domain mobility. The flexibility method, which requires 10^4-fold less
computer power, predicts additional large-scale features of inter-domain motion
that have been observed experimentally. We conclude that the methods offer
complementary insight into the motion of this large protein and provide
detailed structural models that characterise its functionally-significant
conformational changes.


  The budding yeast Saccharomyces cerevisiae is important for human food
production and as a model organism for biological research. The genetic
diversity contained in the global population of yeast strains represents a
valuable resource for a number of fields, including genetics, bioengineering,
and studies of evolution and population structure. Here, we apply a
multiplexed, reduced genome sequencing strategy (known as RAD-seq) to genotype
a large collection of S. cerevisiae strains, isolated from a wide range of
geographical locations and environmental niches. The method permits the
sequencing of the same 1% of all genomes, producing a multiple sequence
alignment of 116,880 bases across 262 strains. We find diversity among these
strains is principally organized by geography, with European, North American,
Asian and African/S. E. Asian populations defining the major axes of genetic
variation. At a finer scale, small groups of strains from cacao, olives and
sake are defined by unique variants not present in other strains. One
population, containing strains from a variety of fermentations, exhibits high
levels of heterozygosity and mixtures of alleles from European and Asian
populations, indicating an admixed origin for this group. In the context of
this global diversity, we demonstrate that a collection of seven strains
commonly used in the laboratory encompasses only one quarter of the genetic
diversity present in the full collection of strains, underscoring the
relatively limited genetic diversity captured by the current set of lab
strains. We propose a model of geographic differentiation followed by
human-associated admixture, primarily between European and Asian populations
and more recently between European and North American populations. The large
collection of genotyped yeast strains characterized here will provide a useful
resource for the broad community of yeast researchers.


  Gene transcription mediated by RNA polymerase II (pol-II) is a key step in
gene expression. The dynamics of pol-II moving along the transcribed region
influence the rate and timing of gene expression. In this work we present a
probabilistic model of transcription dynamics which is fitted to pol-II
occupancy time course data measured using ChIP-Seq. The model can be used to
estimate transcription speed and to infer the temporal pol-II activity profile
at the gene promoter. Model parameters are estimated using either maximum
likelihood estimation or via Bayesian inference using Markov chain Monte Carlo
sampling. The Bayesian approach provides confidence intervals for parameter
estimates and allows the use of priors that capture domain knowledge, e.g. the
expected range of transcription speeds, based on previous experiments. The
model describes the movement of pol-II down the gene body and can be used to
identify the time of induction for transcriptionally engaged genes. By
clustering the inferred promoter activity time profiles, we are able to
determine which genes respond quickly to stimuli and group genes that share
activity profiles and may therefore be co-regulated. We apply our methodology
to biological data obtained using ChIP-seq to measure pol-II occupancy
genome-wide when MCF-7 human breast cancer cells are treated with estradiol
(E2). The transcription speeds we obtain agree with those obtained previously
for smaller numbers of genes with the advantage that our approach can be
applied genome-wide. We validate the biological significance of the pol-II
promoter activity clusters by investigating cluster-specific transcription
factor binding patterns and determining canonical pathway enrichment. We find
that rapidly induced genes are enriched for both estrogen receptor alpha
(ER$\alpha$) and FOXA1 binding in their proximal promoter regions.


  Network analysis provides a powerful framework for the interpretation of
genome-wide data. While static network approaches have proved fruitful, there
is increasing interest in the insights gained from the analysis of cellular
networks under different conditions. In this work, we study the effect of
stress on cellular networks in fission yeast. Stress elicits a sophisticated
and large scale cellular response, involving a shift of resources from cell
growth and metabolism towards protection and maintenance. Previous work has
suggested that these changes can be appreciated at the network level.
  In this paper, we study two types of cellular networks: gene co-regulation
networks and weighted protein interaction networks. We show that in response to
oxidative stress, the co-regulation networks re-organize towards a more
modularised structure: while sets of genes become more tightly co-regulated,
co-regulation between these modules is decreased. This shift translates into
longer average shortest path length, increased transitivity, and decreased
modular overlap in these networks. We also find a similar change in structure
in the weighted protein interaction network in response to both oxidative
stress and nitrogen starvation, confirming and extending previous findings.
  These changes in network structure could represent an increase in network
robustness and/or the emergence of more specialised functional modules.
Additionally, we find stress induces tighter co-regulation of non-coding RNAs,
decreased functional importance of splicing factors, as well as changes in the
centrality of genes involved in chromatin organization, cytoskeleton
organization, cell division, and protein turnover.


  The human Y chromosome exhibits surprisingly low levels of genetic diversity.
This could result from neutral processes if the effective population size of
males is reduced relative to females due to a higher variance in the number of
offspring from males than from females. Alternatively, selection acting on new
mutations, and affecting linked neutral sites, could reduce variability on the
Y chromosome. Here, using genome-wide analyses of X, Y, autosomal and
mitochondrial DNA, in combination with extensive population genetic
simulations, we show that low observed Y chromosome variability is not
consistent with a purely neutral model. Instead, we show that models of
purifying selection are consistent with observed Y diversity. Further, the
number of sites estimated to be under purifying selection greatly exceeds the
number of Y-linked coding sites, suggesting the importance of the highly
repetitive ampliconic regions. While we show that purifying selection removing
deleterious mutations can explain the low diversity on the Y chromosome, we
cannot exclude the possibility that positive selection acting on beneficial
mutations could have also reduced diversity in linked neutral regions, and may
have contributed to lowering human Y chromosome diversity. Because the
functional significance of the ampliconic regions is poorly understood, our
findings should motivate future research in this area.


  In many species, genomic data have revealed pervasive adaptive evolution
indicated by the fixation of beneficial alleles. However, when selection
pressures are highly variable along a species range or through time adaptive
alleles may persist at intermediate frequencies for long periods. So called
balanced polymorphisms have long been understood to be an important component
of standing genetic variation yet direct evidence of the strength of balancing
selection and the stability and prevalence of balanced polymorphisms has
remained elusive. We hypothesized that environmental fluctuations between
seasons in a North American orchard would impose temporally variable selection
on Drosophila melanogaster and consequently maintain allelic variation at
polymorphisms adaptively evolving in response to climatic variation. We
identified hundreds of polymorphisms whose frequency oscillates among seasons
and argue that these loci are subject to strong, temporally variable selection.
We show that these polymorphisms respond to acute and persistent changes in
climate and are associated in predictable ways with seasonally variable
phenotypes. In addition, we show that adaptively oscillating polymorphisms are
likely millions of years old, with some likely predating the divergence between
D. melanogaster and D. simulans. Taken together, our results demonstrate that
rapid temporal fluctuations in climate over generational time promotes adaptive
genetic diversity at loci affecting polygenic phenotypes.


  According to the lipid raft hypothesis, biological lipid membranes are
laterally heterogeneous and filled with nanoscale ordered "raft" domains, which
are believed to play an important role for the organization of proteins in
membranes. However, the mechanisms stabilizing such small rafts are not clear,
and even their existence is sometimes questioned. Here we report the
observation of raft-like structures in a coarse-grained molecular model for
multicomponent lipid bilayers. On small scales, our membranes demix into a
liquid ordered (lo) and a liquid disordered (ld) phase. On large scales, phase
separation is suppressed and gives way to a microemulsion-type state that
contains nanometer size lo domains in a ld environment. Furthermore, we
introduce a mechanism that generates rafts of finite size by a coupling between
monolayer curvature and local composition. We show that mismatch between the
spontaneous curvatures of monolayers in the lo and ld phase induces elastic
interactions, which reduce the line tension between the lo and ld phase and can
stabilize raft domains with a characteristic size of the order of a few
nanometers. Our findings suggest that rafts in multicomponent bilayers might be
closely related to the modulated ripple phase in one-component bilayers.


  In this work we present a flexible tool for tumor progression, which
simulates the evolutionary dynamics of cancer. Tumor progression implements a
multi-type branching process where the key parameters are the fitness
landscape, the mutation rate, and the average time of cell division. The
fitness of a cancer cell depends on the mutations it has accumulated. The input
to our tool could be any fitness landscape, mutation rate, and cell division
time, and the tool produces the growth dynamics and all relevant statistics.


  We study the problem of estimating the origin of an epidemic outbreak --
given a contact network and a snapshot of epidemic spread at a certain time,
determine the infection source. Finding the source is important in different
contexts of computer or social networks. We assume that the epidemic spread
follows the most commonly used susceptible-infected-recovered model. We
introduce an inference algorithm based on dynamic message-passing equations,
and we show that it leads to significant improvement of performance compared to
existing approaches. Importantly, this algorithm remains efficient in the case
where one knows the state of only a fraction of nodes.


  The study of genomic variation has provided key insights into the functional
role of mutations. Predominantly, studies have focused on single nucleotide
variants (SNV), which are relatively easy to detect and can be described with
rich mathematical models. However, it has been observed that genomes are highly
plastic, and that whole regions can be moved, removed or duplicated in bulk.
These structural variants (SV) have been shown to have significant impact on
the phenotype, but their study has been held back by the combinatorial
complexity of the underlying models. We describe here a general model of
structural variation that encompasses both balanced rearrangements and
arbitrary copy-numbers variants (CNV). In this model, we show that the space of
possible evolutionary histories that explain the structural differences between
any two genomes can be sampled ergodically.


  Disease control is of paramount importance in public health with infectious
disease extinction as the ultimate goal. Although diseases may go extinct due
to random loss of effective contacts where the infection is transmitted to new
susceptible individuals, the time to extinction in the absence of control may
be prohibitively long. Thus intervention controls, such as vaccination of
susceptible individuals and/or treatment of infectives, are typically based on
a deterministic schedule, such as periodically vaccinating susceptible children
based on school calendars. In reality, however, such policies are administered
as a random process, while still possessing a mean period. Here, we consider
the effect of randomly distributed intervention as disease control on large
finite populations. We show explicitly how intervention control, based on mean
period and treatment fraction, modulates the average extinction times as a
function of population size and rate of infection spread. In particular, our
results show an exponential improvement in extinction times even though the
controls are implemented using a random Poisson distribution. Finally, we
discover those parameter regimes where random treatment yields an exponential
improvement in extinction times over the application of strictly periodic
intervention. The implication of our results is discussed in light of the
availability of limited resources for control.


  Species do not merely evolve, they also coevolve with other organisms.
Coevolution is a major force driving interacting species to continuously evolve
ex- ploring their fitness landscapes. Coevolution involves the coupling of
species fit- ness landscapes, linking species genetic changes with their
inter-specific ecological interactions. Here we first introduce the Red Queen
hypothesis of evolution com- menting on some theoretical aspects and empirical
evidences. As an introduction to the fitness landscape concept, we review key
issues on evolution on simple and rugged fitness landscapes. Then we present
key modeling examples of coevolution on different fitness landscapes at
different scales, from RNA viruses to complex ecosystems and macroevolution.


  We study adaptive dynamics in games where players abandon the population at a
given rate, and are replaced by naive players characterized by a prior
distribution over the admitted strategies. We demonstrate how such process
leads macroscopically to a variant of the replicator equation, with an
additional term accounting for player turnover. We study how Nash equilibria
and the dynamics of the system are modified by this additional term, for
prototypical examples such as the rock-scissor-paper game and different classes
of two-action games played between two distinct populations. We conclude by
showing how player turnover can account for non-trivial departures from Nash
equilibria observed in data from lowest unique bid auctions.


  Newton set the stage for our view of how science should be done. We remain in
what I will call the `Newtonian Paradigm' in all of physics, including Newton,
Einstein, and Schrodinger. As I will show shortly, Newton invented and
bequeathed to us `efficient cause entailing laws' for the entire becoming of
the universe. With Laplace this became the foundation of contemporary
reductionism in which all that can happen in the world is due to efficient
cause entailing laws. More this framework stands as our dominant way to do
science. The Newtonian Paradigm has done enormous work in science, and helped
lead to the Industrial Revolution, and even our entry into Modernity.
  In this paper I propose to challenge the adequacy of the Newtonian Paradigm
on two ground: 1) For the evolution of the biosphere beyond the watershed of
life, we can formulate no efficient cause entailing laws that allow us to
deduce the evolution of the biosphere. A fortiori, the same holds for the
evolution of the economy, legal systems, social systems, and culture. Because I
have discussed this before with my colleagues Longo and Montevil (1,2) and
elsewhere, (3,4), my discussion of this first point will be rather brief. 2)
What I shall choose to call, after Aristotle's four causes, noted below, Formal
Cause Laws derived from specific `ensemble theories' tell us about the world.
But Formal Cause Laws are not reducible to efficient cause entailing laws of
the Newtonian Paradigm and, critically, have already, unnoticed, crept into
biology concerning the origin of life, and economics concerning economic
growth. Formal cause laws appear to be a new way to do science, independent of
efficient cause entailing laws. Thus Formal Cause laws can be independent of
any specific material substrate. This may bear on the sufficiency of
Materialism in our account of the world.


  We present the phylogeny analysis software SICLE (Sister Clade Extractor), an
easy-to-use, high- throughput tool to describe the nearest neighbors to a node
of interest in a phylogenetic tree as well as the support value for the
relationship. The application is a command line utility that can be embedded
into a phylogenetic analysis pipeline or can be used as a subroutine within
another C++ program. As a test case, we applied this new tool to the published
phylome of Salinibacter ruber, a species of halophilic Bacteriodetes,
identifying 13 unique sister relationships to S. ruber across the 4589 gene
phylogenies. S. ruber grouped with bacteria, most often other Bacteriodetes, in
the majority of phylogenies, but 91 phylogenies showed a branch-supported
sister association between S. ruber and Archaea, an evolutionarily intriguing
relationship indicative of horizontal gene transfer. This test case
demonstrates how SICLE makes it possible to summarize the phylogenetic
information produced by automated phylogenetic pipelines to rapidly identify
and quantify the possible evolutionary relationships that merit further
investigation. SICLE is available for free for noncommercial use at
http://eebweb.arizona.edu/sicle/.


  Moment approximation methods are gaining increasing attention for their use
in the approximation of the stochastic kinetics of chemical reaction systems.
In this paper we derive a general moment expansion method for any type of
propensities and which allows expansion up to any number of moments. For some
chemical reaction systems, more than two moments are necessary to describe the
dynamic properties of the system, which the linear noise approximation (LNA) is
unable to provide. Moreover, also for systems for which the mean does not have
a strong dependence on higher order moments, moment approximation methods give
information about higher order moments of the underlying probability
distribution. We demonstrate the method using a dimerisation reaction,
Michaelis-Menten kinetics and a model of an oscillating p53 system. We show
that for the dimerisation reaction and Michaelis-Menten enzyme kinetics system
higher order moments have limited influence on the estimation of the mean,
while for the p53 system, the solution for the mean can require several moments
to converge to the average obtained from many stochastic simulations. We also
find that agreement between lower order moments does not guarantee that higher
moments will agree. Compared to stochastic simulations our approach is
numerically highly efficient at capturing the behaviour of stochastic systems
in terms of the average and higher moments, and we provide expressions for the
computational cost for different system sizes and orders of approximation. {We
show how the moment expansion method can be employed to efficiently {quantify
parameter sensitivity}.} Finally we investigate the effects of using too few
moments on parameter estimation, and provide guidance on how to estimate if the
distribution can be accurately approximated using only a few moments.


  Background: It is well known that the development of cancer is caused by the
accumulation of somatic mutations within the genome. For oncogenes
specifically, current research suggests that there is a small set of "driver"
mutations that are primarily responsible for tumorigenesis. Further, due to
some recent pharmacological successes in treating these driver mutations and
their resulting tumors, a variety of methods have been developed to identify
potential driver mutations using methods such as machine learning and
mutational clustering. We propose a novel methodology that increases our power
to identify mutational clusters by taking into account protein tertiary
structure via a graph theoretical approach.
  Results: We have designed and implemented GraphPAC (Graph Protein Amino Acid
Clustering) to identify mutational clustering while considering protein spatial
structure. Using GraphPAC, we are able to detect novel clusters in proteins
that are known to exhibit mutation clustering as well as identify clusters in
proteins without evidence of prior clustering based on current methods.
Specifically, by utilizing the spatial information available in the Protein
Data Bank (PDB) along with the mutational data in the Catalogue of Somatic
Mutations in Cancer (COSMIC), GraphPAC identifies new mutational clusters in
well known oncogenes such as EGFR and KRAS. Further, by utilizing graph theory
to account for the tertiary structure, GraphPAC identifies clusters in DPP4,
NRP1 and other proteins not identified by existing methods. The R package is
available at: http://bioconductor.org/packages/release/bioc/html/GraphPAC.html
  Conclusion: GraphPAC provides an alternative to iPAC and an extension to
current methodology when identifying potential activating driver mutations by
utilizing a graph theoretic approach when considering protein tertiary
structure.


  We provide quantitative estimates on the size and the structure of the
epistatic space defined in the main article, "Dynamical Correlations in the
escape strategy of Influenza A virus", EPL 101 68003.


  We apply the "weighted ensemble" (WE) simulation strategy, previously
employed in the context of molecular dynamics simulations, to a series of
systems-biology models that range in complexity from one-dimensional to a
system with 354 species and 3680 reactions. WE is relatively easy to implement,
does not require extensive hand-tuning of parameters, does not depend on the
details of the simulation algorithm, and can facilitate the simulation of
extremely rare events.
  For the coupled stochastic reaction systems we study, WE is able to produce
accurate and efficient approximations of the joint probability distribution for
all chemical species for all time t. WE is also able to efficiently extract
mean first passage times for the systems, via the construction of a
steady-state condition with feedback. In all cases studied here, WE results
agree with independent calculations, but significantly enhance the precision
with which rare or slow processes can be characterized. Speedups over
"brute-force" in sampling rare events via the Gillespie direct Stochastic
Simulation Algorithm range from ~10^12 to ~10^20 for rare states in a
distribution, and ~10^2 to ~10^4 for finding mean first passage times.


  We use a stochastic metapopulation model to study the combined effects of
seasonality and spatial heterogeneity on disease persistence. We find a
pronounced effect of enhanced persistence associated with strong heterogeneity,
intermediate coupling strength and moderate seasonal forcing. Analytic
calculations show that this effect is not related with the phase lag between
epidemic bursts in different patches, but rather with the linear stability
properties of the attractor that describes the steady state of the system in
the large population limit.


  The terrestrial fossil record shows a significant variation in the extinction
and origination rates of species during the past half billion years. Numerous
studies have claimed an association between this variation and the motion of
the Sun around the Galaxy, invoking the modulation of cosmic rays, gamma rays
and comet impact frequency as a cause of this biodiversity variation. However,
some of these studies exhibit methodological problems, or were based on coarse
assumptions (such as a strict periodicity of the solar orbit). Here we
investigate this link in more detail, using a model of the Galaxy to
reconstruct the solar orbit and thus a predictive model of the temporal
variation of the extinction rate due to astronomical mechanisms. We compare
these predictions as well as those of various reference models with
paleontological data. Our approach involves Bayesian model comparison, which
takes into account the uncertainties in the paleontological data as well as the
distribution of solar orbits consistent with the uncertainties in the
astronomical data. We find that various versions of the orbital model are not
favored beyond simpler reference models. In particular, the distribution of
mass extinction events can be explained just as well by a uniform random
distribution as by any other model tested. Although our negative results on the
orbital model are robust to changes in the Galaxy model, the Sun's coordinates
and the errors in the data, we also find that it would be very difficult to
positively identify the orbital model even if it were the true one. (In
contrast, we do find evidence against simpler periodic models.) Thus while we
cannot rule out there being some connection between solar motion and
biodiversity variations on the Earth, we conclude that it is difficult to give
convincing positive conclusions of such a connection using current data.


  It is widely held that a substantial genetic component underlies Bipolar
Disorder (BD) and other neuropsychiatric disease traits. Recent efforts have
been aimed at understanding the genetic basis of disease susceptibility, with
genome-wide association studies (GWAS) unveiling some promising associations.
Nevertheless, the genetic etiology of BD remains elusive with a substantial
proportion of the heritability - which has been estimated to be 80% based on
twin and family studies - unaccounted for by the specific genetic variants
identified by large-scale GWAS. Furthermore, functional understanding of
associated loci generally lags discovery. Studies we report here provide
considerable support to the claim that substantially more remains to be gained
from GWAS on the genetic mechanisms underlying BD susceptibility, and that a
large proportion of the variation in disease risk may be uncovered through
integrative functional genomic approaches. We combine recent analytic advances
in heritability estimation and polygenic modeling and leverage recent
technological advances in the generation of -omics data to evaluate the nature
and scale of the contribution of functional classes of genetic variation to a
relatively intractable disorder. We identified cis eQTLs in cerebellum and
parietal cortex that capture more than half of the total heritability
attributable to SNPs interrogated through GWAS and showed that eQTL-based
heritability estimation is highly tissue-dependent. Our findings show that a
much greater resolution may be attained than has been reported thus far on the
number of common loci that capture a substantial proportion of the heritability
to disease risk and that the functional nature of contributory loci may be
clarified en masse.


  In single molecule laser optical tweezer (LOT) pulling experiments a protein
or RNA is juxtaposed between DNA handles that are attached to beads in optical
traps. The LOT generates folding trajectories under force in terms of
time-dependent changes in the distance between the beads. How to construct the
full intrinsic folding landscape (without the handles and the beads) from the
measured time series is a major unsolved problem. By using rigorous theoretical
methods---which account for fluctuations of the DNA handles, rotation of the
optical beads, variations in applied tension due to finite trap stiffness, as
well as environmental noise and the limited bandwidth of the apparatus---we
provide a tractable method to derive intrinsic free energy profiles. We
validate the method by showing that the exactly calculable intrinsic free
energy profile for a Generalized Rouse Model, which mimics the two-state
behavior in nucleic acid hairpins, can be accurately extracted from simulated
time series in a LOT setup regardless of the stiffness of the handles. We next
apply the approach to trajectories from coarse grained LOT molecular
simulations of a coiled-coil protein based on the GCN4 leucine zipper, and
obtain a free energy landscape that is in quantitative agreement with
simulations performed without the beads and handles. Finally, we extract the
intrinsic free energy landscape from experimental LOT measurements for the
leucine zipper, which is independent of the trap parameters.


  Thermal oscillations of base pairs in the Peyrard-Bishop-Holstein model are
simulated by stochastic fluctuations of base overlap integrals. Numerical
investigation of the model is carried out for a hole transfer in $G_1A_2G_3G_4$
sequence which was previously studied experimentally by F. Lewis et al. A hole
migration between $G_1$ and $G_3G_4$ is determined by the matrix elements of
the charge transition, but presence and amplitude of their stochastic
fluctuations proved to play a key role in reproduction of the experimental
kinetics. Good agreement with the experimental data was obtained for a wide
range of the model parameters' combinations.


  We analyze and reinterpret experimental evidence from the literature to argue
for an ability of tumor cells to self-regulate their aneuploidy rate. We
conjecture that this ability is mediated by a diversification factor that
exploits molecular mechanisms common to embryo stem cells and, to a lesser
extent, adult stem cells, that is eventually reactivated in tumor cells.
Moreover, we propose a direct use of the quasispecies model to cancer cells
based on their significant genomic instability (i.e. aneuploidy rate), by
defining master sequences lengths as the sum of all copy numbers of physically
distinct whole and fragmented chromosomes. We compute an approximate error
threshold such that any aneuploidy rate larger than the threshold would lead to
a loss of fitness of a tumor population, and we confirm that highly aneuploid
cancer populations already function with aneuploidy rates close to the
estimated threshold.


  We investigated the thermodynamic relaxation of giant unilamellar vesicles
(GUVs) which contained small vesicles within their interior. Quenching these
vesicles from their fluid phase (T>Tm) through the phase transition in the gel
state (T<Tm) drives the inner vesicles to be expelled from the larger mother
vesicle via the accompanying decrease in vesicle area by ~25% which forces a
pore to open in the mother vesicle. We demonstrate that the proceeding time
evolution of the resulting efflux follows the relaxation of the membrane area
and describe the entire relaxation process using an Onsager-like nonequilibrium
thermodynamics ansatz. As a consequence of the volume efflux internal vesicles
are expelled from the mother vesicle. Although complete sealing of the pore may
occur during the expulsion, the global relaxation dynamics is conserved.
Finally, comparison of these results to morphological relaxation phenomena
found in earlier studies reveals a universal relaxation behaviour in GUVs. When
quenched from the fluid to gel phase the typical time scale of relaxation shows
little variation and ranges between 4-5 seconds.


  Membrane budding has been extensively studied as an equilibrium process
attributed to the formation of coexisting domains or changes in the vesicle
area to volume ratio (reduced volume). In contrast, non-equilibrium budding
remains experimentally widely unexplored especially when time scales fall well
below the characteristic diffusion time of lipids{\tau} . We show that
localized mechanical perturbations, initiated by driving giant unilamellar
vesicles (GUVs) through their lipid phase transition, leads to the immediate
formation of rapidly growing, multiply localized, non-equilibrium buds, when
the transition takes place at short timescales (<{\tau}). We show that these
buds arise from small fluid-like perturbations and grow as spherical caps in
the third dimension, since in plane spreading is obstructed by the continuous
rigid gel-like matrix. Accounting for both three and two dimensional viscosity,
we demonstrate that dissipation decreases the size scale of the system and
therefore favours the formation of multiple buds as long as the perturbation
takes place above a certain critical rate. This rate depends on membrane and
media viscosity and is qualitatively and quantitatively correctly predicted by
our theoretical description.


  It has been known that epidemic outbreaks in the SIR model on networks are
described by phase transitions. Despite the similarity with percolation
transitions, whether an epidemic outbreak occurs or not cannot be predicted
with probability one in the thermodynamic limit. We elucidate its mechanism by
deriving a simple Langevin equation that captures an essential aspect of the
phenomenon. We also calculate the probability of epidemic outbreaks near the
transition point.


  We consider mechanical stability of dimeric and monomeric proteins with the
cystine knot motif. A structure based dynamical model is used to demonstrate
that all dimeric and some monomeric proteins of this kind should have
considerable resistance to stretching that is significantly larger than that of
titin. The mechanisms of the large mechanostability are elucidated. In most
cases, it originates from the induced formation of one or two cystine
slipknots. Since there are four termini in a dimer, there are several ways of
selecting two of them to pull by. We show that in the cystine knot systems,
there is strong anisotropy in mechanostability and force patterns related to
the selection. We show that the thermodynamic stability of the dimers is
enhanced compared to the constituting monomers whereas machanostability is
either lower or higher.


  Population pharmacokinetic (PK) modeling methods can be statistically
classified as either parametric or nonparametric (NP). Each classification can
be divided into maximum likelihood (ML) or Bayesian (B) approaches. In this
paper we discuss the nonparametric case using both maximum likelihood and
Bayesian approaches. We present two nonparametric methods for estimating the
unknown joint population distribution of model parameter values in a
pharmacokinetic/pharmacodynamic (PK/PD) dataset. The first method is the NP
Adaptive Grid (NPAG). The second is the NP Bayesian (NPB) algorithm with a
stick-breaking process to construct a Dirichlet prior. Our objective is to
compare the performance of these two methods using a simulated PK/PD dataset.
Our results showed excellent performance of NPAG and NPB in a realistically
simulated PK study. This simulation allowed us to have benchmarks in the form
of the true population parameters to compare with the estimates produced by the
two methods, while incorporating challenges like unbalanced sample times and
sample numbers as well as the ability to include the covariate of patient
weight. We conclude that both NPML and NPB can be used in realistic PK/PD
population analysis problems. The advantages of one versus the other are
discussed in the paper. NPAG and NPB are implemented in R and freely available
for download within the Pmetrics package from www.lapk.org.


  Mathematical models are increasingly being used to understand complex
biochemical systems, to analyze experimental data and make predictions about
unobserved quantities. However, we rarely know how robust our conclusions are
with respect to the choice and uncertainties of the model. Using algebraic
techniques we study systematically the effects of intermediate, or transient,
species in biochemical systems and provide a simple, yet rigorous mathematical
classification of all models obtained from a core model by including
intermediates. Main examples include enzymatic and post-translational
modification systems, where intermediates often are considered insignificant
and neglected in a model, or they are not included because we are unaware of
their existence. All possible models obtained from the core model are
classified into a finite number of classes. Each class is defined by a
mathematically simple canonical model that characterizes crucial dynamical
properties, such as mono- and multistationarity and stability of steady states,
of all models in the class. We show that if the core model does not have
conservation laws, then the introduction of intermediates does not change the
steady-state concentrations of the species in the core model, after suitable
matching of parameters. Importantly, our results provide guidelines to the
modeler in choosing between models and in distinguishing their properties.
Further, our work provides a formal way of comparing models that share a common
skeleton.


  Transcription factor (TF) molecules translocate by facilitated diffusion (a
combination of 3D diffusion around and 1D random walk on the DNA). Despite the
attention this mechanism received in the last 40 years, only a few studies
investigated the influence of the cellular environment on the facilitated
diffusion mechanism and, in particular, the influence of `other' DNA binding
proteins competing with the TF molecules for DNA space. Molecular crowding on
the DNA is likely to influence the association rate of TFs to their target site
and the steady state occupancy of those sites, but it is still not clear how it
influences the search in a genome-wide context, when the model includes
biologically relevant parameters (such as: TF abundance, TF affinity for DNA
and TF dynamics on the DNA). We performed stochastic simulations of TFs
performing the facilitated diffusion mechanism, and considered various
abundances of cognate and non-cognate TFs. We show that, for both obstacles
that move on the DNA and obstacles that are fixed on the DNA, changes in search
time are not statistically significant in case of biologically relevant
crowding levels on the DNA. In the case of non-cognate proteins that slide on
the DNA, molecular crowding on the DNA always leads to statistically
significant lower levels of occupancy, which may confer a general mechanism to
control gene activity levels globally. When the `other' molecules are immobile
on the DNA, we found a completely different behaviour, namely: the occupancy of
the target site is always increased by higher molecular crowding on the DNA.
Finally, we show that crowding on the DNA may increase transcriptional noise
through increased variability of the occupancy time of the target sites.


  Transcription factors (TFs) are proteins that bind to specific sites on the
DNA and regulate gene activity. Identifying where TF molecules bind and how
much time they spend on their target sites is key for understanding
transcriptional regulation. It is usually assumed that the free energy of
binding of a TF to the DNA (the affinity of the site) is highly correlated to
the amount of time the TF remains bound (the occupancy of the site). However,
knowing the binding energy is not sufficient to infer actual binding site
occupancy. This mismatch between the occupancy predicted by the affinity and
the observed occupancy may be caused by various factors, such as TF abundance,
competition between TFs or the arrangement of the sites on the DNA. We
investigated the relationship between the affinity of a TF for a set of binding
sites and their occupancy. In particular, we considered the case of lac
repressor (lacI) in E.coli and performed stochastic simulations of the TF
dynamics on the DNA for various combinations of lacI abundance in competition
with TFs that contribute to macromolecular crowding. Our results showed that
for medium and high affinity sites, TF competition does not play a significant
role in genomic occupancy, except in cases when the abundance of lacI is
significantly increased or when a low-information content PWM was used.
Nevertheless, for medium and low affinity sites, an increase in TF abundance
(for both lacI or other molecules) leads to an increase in occupancy at several
sites. Keywords: facilitated diffusion, Position Weight Matrix, thermodynamic
equilibrium, motif information content, molecular crowding


  A model with six mutually-exclusive compartments related to dengue is
studied. Three vector control tools are considered: insecticides (larvicide and
adulticide) and mechanical control. The basic reproduction number associated to
the model is presented. The problem is studied using an optimal control
approach. The human data is based on the dengue outbreak that occurred in Cape
Verde. Control measures are simulated in different scenarios and their
consequences analyzed.


  The Red cell-Hb-CSF functions as a sensor adapting response to Hb
heterotropic equilibriums. At the lungs O2 and Mg2+, each one increasing
affinity for the other stabilize the relax (R) form [(O2)4Hb(Mg)2].(H2O)R. At
tissue level, the inclusion of H+ and 2,3-DPG excludes O2 and Mg2+ to stabilize
the tense (T) form 2,3-DPG-deoxyHb-(H2O)T. Both senses are integrated into a
cycle T into R and R into T, without involving a direct reversal. The
dissipative potential of water cluster (H2O)n interacts with the hydrophilic
asymmetries of Hb, to restrict randomness of the kinetic sense implicated in a
single peak for activation energy (Ea). The hydration shells could sequence an
enhanced Ea into several peaks, to sequentially activate transitions states.
Hence, changes in dipole state, sliding, pKa, n-H-bonds, etc., could became
concatenated for vectoriality. (H2O)n by the loss of H-bonds couple with to the
hydration turnover of proteins and ions to result in incomplete water cluster
(H2O)n*, with a lower-n. (H2O)n* became a carrier of heat/entropy into the
cerebrospinal fluid (CSF) which has to be replaced 3.7 times per day. OxyHb
formation involves sliding-down of alpha vs beta chains, to shift alpha1 and
alpha2 Pro 44 into allowing the entrance of a fully hydrated
[Mg.(H2O)6](H2O)12-14(2+) (or Zn2+) into the hydrophilic beta2-alpha1 and
beta1-alpha2 interfaces. OxyHb pKa of 6.4 leads to H+-dissociation increasing
negative charge of R-groups. This at beta2-alpha1 sequence two tetradentate
chelates, first an Mg2+, bonding with beta2 His 92 and a second Mg2+ with
alpha1 His 87, to cooperatively release hindrance. The interconversion of
oxy-to-deoxyHb, pKa=8, leads to the amphoteric imidazole to became positively
charged and proximal histidines return into hindrance position, releasing the
incomplete hydrated Mg.(H2O)inc(2+) and O2 into CSF.


  Speciation is of fundamental importance to understanding the huge diversity
of life on Earth. In contrast to current phenomenological models, we develop a
biophysically motivated approach to study speciation involving the co-evolution
of protein binding DNA for two geographically isolated populations. Our results
predict that, despite neutral diffusion of hybrids in trait space, smaller
populations have a higher rate of speciation, due to sequence entropy poising
populations more closely to incompatible regions of phenotype space. A key
lesson of this work is that non-trivial contributions of sequence entropy give
rise to a strong population size dependence on speciation rates.


  Collective decision-making in biological systems requires all individuals in
the group to go through a behavioural change of state. During this transition,
the efficiency of information transport is a key factor to prevent cohesion
loss and preserve robustness. The precise mechanism by which natural groups
achieve such efficiency, though, is currently not fully understood. Here, we
present an experimental study of starling flocks performing collective turns in
the field. We find that the information to change direction propagates across
the flock linearly in time with negligible attenuation, hence keeping group
decoherence to a minimum. This result contrasts with current theories of
collective motion, which predict a slower and dissipative transport of
directional information. We propose a novel theory whose cornerstone is the
existence of a conserved spin current generated by the gauge symmetry of the
system. The theory turns out to be mathematically identical to that of
superfluid transport in liquid helium and it explains the dissipationless
propagating mode observed in turning flocks. Superfluidity also provides a
quantitative expression for the speed of propagation of the information,
according to which transport must be swifter the stronger the group's
orientational order. This prediction is verified by the data. We argue that the
link between strong order and efficient decision-making required by
superfluidity may be the adaptive drive for the high degree of behavioural
polarization observed in many living groups. The mathematical equivalence
between superfluid liquids and turning flocks is a compelling demonstration of
the far-reaching consequences of symmetry and conservation laws across
different natural systems.


  Analyzing coexistence and survival scenarios of Lotka-Volterra (LV) networks
in which the total biomass is conserved is of vital importance for the
characterization of long-term dynamics of ecological communities. Here, we
introduce a classification scheme for coexistence scenarios in these
conservative LV models and quantify the extinction process by employing the
Pfaffian of the network's interaction matrix. We illustrate our findings on
global stability properties for general systems of four and five species and
find a generalized scaling law for the extinction time.


  The famous series of Fibonacci numbers is defined by a recursive equation
saying that each number is the sum of its two predecessors, with the initial
condition that the first two numbers are equal to unity. Here, we show that the
numbers of fatty acids (straight-chain aliphatic monocarboxylic acids) with n
carbon atoms is exactly given by the Fibonacci numbers. Thus, by investing one
more carbon atom into extending a fatty acid, an organism can increase the
variability of the fatty acids approximately by the factor of the Golden
section, 1.618. As the Fibonacci series grows asymptotically exponentially, our
results are in line with combinatorial complexity found generally in biology.
We also outline potential extensions of the calculations to modified (e.g.,
hydroxylated) fatty acids. The presented enumeration method may be of interest
for lipidomics, combinatorial chemistry, synthetic biology and the theory of
evolution (including prebiotic evolution).


  The aim of this work was to show few examples and few perspective of modeling
in epidemiology. We began with differential equations which were a first tool
to describe and predict that phenomena. Wroclaw as a cite was very important,
because statistics from smallpox epidemic were used by Bernoulli to estimate
parameters of first mathematical model of epidemic. Next step were SIR models
and those also appeared first as differential equations. They were very popular
in begin of XX century. When computer simulation changed the world of
mathematical modeling agent-based models gave more possibilities in
epidemiology. That models have a big privilege on differential equation,
because of information of social network people habits and reaction on
infections, which can me involved in agent-based models as well as governmental
intervention. We showed in this work how that human relations are important in
transmitting diseases and there is example, where it is possible to conduct
experiments of significant policy relevance (vaccinating), such as
investigating the initial growth of an epidemic on a real-world network.
Presented H1N1 model could be observed in real time (prediction was made in
September 2009 for winter season 2009/2010) what make it more exiting and also
practice.


  We show that the resolution of social dilemmas on random graphs and
scale-free networks is facilitated by imitating not the strategy of better
performing players but rather their emotions. We assume sympathy and envy as
the two emotions that determine the strategy of each player by any given
interaction, and we define them as probabilities to cooperate with players
having a lower and higher payoff, respectively. Starting with a population
where all possible combinations of the two emotions are available, the
evolutionary process leads to a spontaneous fixation to a single emotional
profile that is eventually adopted by all players. However, this emotional
profile depends not only on the payoffs but also on the heterogeneity of the
interaction network. Homogeneous networks, such as lattices and regular random
graphs, lead to fixations that are characterized by high sympathy and high
envy, while heterogeneous networks lead to low or modest sympathy but also low
envy. Our results thus suggest that public emotions and the propensity to
cooperate at large depend, and are in fact determined by the properties of the
interaction network.


  To explain the differentiation of stem cells in terms of dynamical systems
theory, models of interacting cells with intracellular protein expression
dynamics are analyzed and simulated. Simulations were carried out for all
possible protein expression networks consisting of two genes under cell--cell
interactions mediated by the diffusion of a protein. Networks that show cell
differentiation are extracted and two forms of symmetric differentiation based
on Turing's mechanism and asymmetric differentiation are identified. In the
latter network, the intracellular protein levels show oscillatory dynamics at a
single-cell level, while cell-to-cell synchronicity of the oscillation is lost
with an increase in the number of cells. Differentiation to a fixed-point type
behavior follows with a further increase in the number of cells. The cell type
with oscillatory dynamics corresponds to a stem cell that can both proliferate
and differentiate, while the latter fixed-point type only proliferates. This
differentiation is analyzed as a saddle-node bifurcation on an invariant
circle, while the number ratio of each cell type is shown to be robust against
perturbations due to self-consistent determination of the effective bifurcation
parameter as a result of the cell--cell interaction. Complex cell
differentiation is designed by combing these simple two-gene networks. The
generality of the present differentiation mechanism, as well as its biological
relevance, is discussed.


  We investigate the properties of large random conservative chemical reaction
networks composed of elementary reactions endowed with either mass-action or
saturating kinetics, assigning kinetic parameters in a
thermodynamically-consistent manner. We find that such complex networks exhibit
qualitatively similar behavior when fed with external nutrient flux. The
nutrient is preferentially transformed into one specific chemical that is an
intrinsic property of the network. We propose a self-consistent proto-cell toy
model in which the preferentially synthesized chemical is a precursor for the
cell membrane, and show that such proto-cells can exhibit sustainable
homeostatic growth when fed with any nutrient diffusing through the membrane,
provided that nutrient is metabolized at a sufficient rate.


  We propose a method that uses genetic data to test for the occurrence of a
recent range expansion and to infer the location of the origin of the
expansion. We introduce a statistic for pairs of populations $\psi$ (the
directionality index) that detects asymmetries in the two-dimensional allele
frequency spectrum caused by the series of founder events that happen during an
expansion. Such asymmetry arises because low frequency alleles tend to be lost
during founder events, thus creating clines in the frequencies of surviving
low-frequency alleles. Using simulations, we further show that $\psi$ is more
powerful for detecting range expansions than both $F_{ST}$ and clines in
heterozygosity. We illustrate the utility of $\psi$ by applying it to a data
set from modern humans and show how we can include more complicated scenarios
such as multiple expansion origins or barriers to migration in the model.


  Adenosine triphosphate (ATP)-binding cassette (ABC) transporters form a
family of molecular motor proteins that couple ATP hydrolysis to substrate
translocation across cell membranes. Each nucleotide binding domain of
ABC-transporters contains a highly conserved H-loop Histidine residue, whose
precise mechanistic role to motor functions has remained elusive. By using
combined quantum mechanical and molecular mechanical calculations, we showed
that the conserved H-loop residue H662 in E. coli HlyB, a bacterial ABC
transporter, can act first as a general acid and then as a general base to
facilitate proton transfers in ATP hydrolysis. Without the assistance of H662,
direct proton transfer from the lytic water to ATP results in a greatly
elevated barrier height. Our findings suggest that the essential function of
the H-loop residue H662 is to provide a "chemical linchpin" that shuttles
protons between reactants through a relay mechanism, thereby catalyzing ATP
hydrolysis in HlyB.


  Understanding the population structure and patterns of gene flow within
species is of fundamental importance to the study of evolution. In the fields
of population and evolutionary genetics, measures of genetic differentiation
are commonly used to gather this information. One potential caveat is that
these measures assume gene flow to be symmetric. However, asymmetric gene flow
is common in nature, especially in systems driven by physical processes such as
wind or water currents. Since information about levels of asymmetric gene flow
among populations is essential for the correct interpretation of the
distribution of contemporary genetic diversity within species, this should not
be overlooked. To obtain information on asymmetric migration patterns from
genetic data, complex models based on maximum likelihood or Bayesian approaches
generally need to be employed, often at great computational cost. Here, a new
simpler and more efficient approach for understanding gene flow patterns is
presented. This approach allows the estimation of directional components of
genetic divergence between pairs of populations at low computational effort,
using any of the classical or modern measures of genetic differentiation. These
directional measures of genetic differentiation can further be used to
calculate directional relative migration and to detect asymmetries in gene flow
patterns. This can be done in a user-friendly web application called
divMigrate-online introduced in this paper. Using simulated data sets with
known gene flow regimes, we demonstrate that the method is capable of resolving
complex migration patterns under a range of study designs.


  Most ectotherms show an inverse relationship between developmental
temperature and body size, a phenomenon known as the temperature size rule
(TSR). Several competing hypotheses have been proposed to explain its
occurrence. According to one set of views, the TSR results from inevitable
biophysical effects of temperature on the rates of growth and differentiation,
whereas other views suggest the TSR is an adaptation that can be achieved by a
diversity of mechanisms in different taxa. Our data reveal that the fruit fly,
Drosophila melanogaster, obeys the TSR using a novel mechanism: reduction of
critical size at higher temperatures. In holometabolous insects, attainment of
critical size initiates the hormonal cascade that terminates growth, and hence,
Drosophila larvae appear to instigate the signal to stop growth at a smaller
size at higher temperatures. This is in contrast to findings from another
holometabolous insect, Manduca sexta, in which the TSR results from the effect
of temperature on the rate and duration of growth. This contrast suggests that
there is no single mechanism that accounts for the TSR. Instead, the TSR
appears to be an adaptation that is achieved at a proximate level through
different mechanisms in different taxa.


  Building circuits and studying their behavior in cells is a major goal of
systems and synthetic biology. Synthetic biology enables the precise control of
cellular states for systems studies, the discovery of novel parts, control
strategies, and interactions for the design of robust synthetic systems. To the
best of our knowledge,there are no literature reports for the synthetic circuit
construction for protozoan parasites. This paper describes the construction of
genetic circuit for the targeted enzyme inositol phosphorylceramide synthase
belonging to the protozoan parasite Leishmania. To explore the dynamic nature
of the circuit designed, simulation was done followed by circuit validation by
qualitative and quantitative approaches. The genetic circuit designed for
inositol phosphorylceramide synthase shows responsiveness, oscillatory and
bistable behavior, together with intrinsic robustness.


  In this thesis we present a mathematical model describing the population
dynamics of molecules in an artificial chemistry where large molecules can be
produced by successive ligation of pairs of smaller molecules. The chemistry
contains a large number of spontaneous reactions of which a small subset could
be catalyzed by molecules produced in the chemistry with varying catalytic
strengths. We show ACSs, if present in the catalytic network, can focus the
resources of the system into a sparse set of molecules. ACSs can produce a
bistability in the population dynamics and, in particular, steady states
wherein the ACS molecules dominate the population, i.e., have higher
concentrations compared to the rest of molecules in the chemistry (background).
In this thesis we attempt to address two main questions: First, under what
circumstances do molecules belonging to the ACSs dominate over the background,
and second, starting from an initial condition that does not contain good
catalysts, can a sparse set of large molecules (containing several tens or a
few hundred monomers) that are good catalysts arise and be maintained in the
system at concentrations significantly above the background? We show that if an
ACS catalyzed by large molecules contains within it (or partially overlaps
with) a smaller ACS catalyzed by smaller molecules (referred to as a `nested
ACS' structure), the catalytic strength required for the large ACS to dominate
comes down significantly. We show that when the network contains a cascade of
nested ACSs with the catalytic strengths of molecules increasing gradually with
their size, a sparse subset of molecules including some very large molecules
can come to dominate the system.


  Confirmed cases during the early stage of the 2009 H1N1 pdm in various
countries showed an age shift between importations and local transmission
cases, with adults mainly responsible for seeding unaffected regions and
children most frequently driving community outbreaks. We introduce a multi-host
stochastic metapopulation model with two age classes to analytically address
the role of a heterogeneously mixing population and its associated
non-homogeneous travel behaviors on the risk of a major epidemic. We inform the
model with statistics on demography, mixing and travel behavior for Europe and
Mexico, and calibrate it to the 2009 H1N1 pdm early outbreak. We varied model
parameters to explore the invasion conditions under different scenarios. We
derive the expression for the global invasion potential of the epidemic that
depends on disease transmissibility, transportation network and mobility
features, demographic profile and mixing pattern. Highly assortative mixing
favor the spatial containment of the epidemic, this effect being contrasted by
an increase in the social activity of adults vs. children. Heterogeneity of the
mobility network topology and traffic flows strongly favor the disease
invasion, as also a larger fraction of children traveling. Variations in the
demography and mixing habits across countries lead to heterogeneous outbreak
situations. Results are compatible with the H1N1 spatial spread observed. The
work illustrates the importance of age-dependent mixing profiles and mobility
features in the study of the conditions for the spatial invasion of an emerging
influenza pandemic. Its results allow the immediate assessment of the risk of a
major epidemic for a specific scenario upon availability of data, and the
evaluation of the effectiveness of public health interventions targeting
specific age groups, their interactions and mobility behaviors.


  In this paper we introduce a novel parallel pipeline for fast and accurate
mapping of RNA sequences on servers equipped with multicore processors. Our
software, named HPG-Aligner, leverages the speed of the Burrows-Wheeler
Transform to map a large number of RNA fragments (reads) rapidly, as well as
the accuracy of the Smith-Waterman algorithm, that is employed to deal with
conflictive reads. The aligner is complemented with a careful strategy to
detect splice junctions based on the division of RNA reads into short segments
(or seeds), which are then mapped onto a number of candidate alignment
locations, providing useful information for the successful alignment of the
complete reads.
  Experimental results on platforms with AMD and Intel multicore processors
report the remarkable parallel performance of HPG-Aligner, on short and long
RNA reads, which excels in both execution time and sensitivity to an
state-of-the-art aligner such as TopHat 2 built on top of Bowtie and Bowtie 2.


  We propose a simple approach, based on the minimization of the total
(entropic plus unfolding) energy of a two-state system, describing the
stretch-induced unfolding of macromolecules (proteins, silks, nanopolymers,
DNA/RNA). The model is fully analytical and enlightens the role of the
different energetic components regulating the unfolding evolution. As an
explicit example of application we compare the analytical results with the
titin Atomic Force Microscopy experiments showing the ability of the model to
quantitatively reproduce the mechanical behavior of macromolecules unfolding.


  The study of functional genomics--particularly in non-model organisms has
been dramatically improved over the last few years by use of transcriptomes and
RNAseq. While these studies are potentially extremely powerful, a
computationally intensive procedure--the de novo construction of a reference
transcriptome must be completed as a prerequisite to further analyses. The
accurate reference is critically important as all downstream steps, including
estimating transcript abundance are critically dependent on the construction of
an accurate reference. Though a substantial amount of research has been done on
assembly, only recently have the pre-assembly procedures been studied in
detail. Specifically, several stand-alone error correction modules have been
reported on, and while they have shown to be effective in reducing errors at
the level of sequencing reads, how error correction impacts assembly accuracy
is largely unknown. Here, we show via use of a simulated dataset, that applying
error correction to sequencing reads has significant positive effects on
assembly accuracy, by reducing assembly error by nearly 50%, and therefore
should be applied to all datasets.


  Chemical reaction networks (CRNs) formally model chemistry in a well-mixed
solution. CRNs are widely used to describe information processing occurring in
natural cellular regulatory networks, and with upcoming advances in synthetic
biology, CRNs are a promising programming language for the design of artificial
molecular control circuitry. Due to a formal equivalence between CRNs and a
model of distributed computing known as population protocols, results transfer
readily between the two models.
  We show that if a CRN respects finite density (at most O(n) additional
molecules can be produced from n initial molecules), then starting from any
dense initial configuration (all molecular species initially present have
initial count Omega(n), where n is the initial molecular count and volume),
then every producible species is produced in constant time with high
probability.
  This implies that no CRN obeying the stated constraints can function as a
timer, able to produce a molecule, but doing so only after a time that is an
unbounded function of the input size. This has consequences regarding an open
question of Angluin, Aspnes, and Eisenstat concerning the ability of population
protocols to perform fast, reliable leader election and to simulate arbitrary
algorithms from a uniform initial state.


  Recent advances in non-equilibrium statistical mechanics and single molecule
technologies make it possible to extract free energy differences from
irreversible work measurements in pulling experiments. To date, free energy
recovery has been focused on native or equilibrium molecular states, whereas
free energy measurements of kinetic states (i.e. finite lifetime states that
are generated dynamically and are metastable) have remained unexplored. Kinetic
states can play an important role in various domains of physics, such as
nanotechnology or condensed matter physics. In biophysics, there are many
examples where they determine the fate of molecular reactions: protein and
peptide-nucleic acid binding, specific cation binding, antigen-antibody
interactions, transient states in enzymatic reactions or the formation of
transient intermediates and non-native structures in molecular folders. Here we
demonstrate that it is possible to obtain free energies of kinetic states by
applying extended fluctuation relations. This is shown by using optical
tweezers to mechanically unfold and refold DNA structures exhibiting
intermediate and misfolded kinetic states.


  The Histone-like Nucleoid Structuring protein (H-NS) is a nucleoid-associated
protein, which is involved in both gene regulation and DNA compaction. H-NS can
bind to DNA in two different ways: in trans, by binding to two separate DNA
duplexes, or in cis, by binding to different sites on the same duplex. Based on
scanning force microscopy imaging and optical trap-driven unzipping assays, it
has recently been suggested that DNA compaction may result from the
antagonistic effects of H-NS binding to DNA in trans and cis configurations. In
order to get more insight into the compaction mechanism, we constructed a
coarse-grained model description of the compaction of bacterial DNA by H-NS.
These simulations highlight the fact that DNA compaction indeed results from
the subtle equilibrium between several competing factors, which include the
deformation dynamics of the plasmid and the several binding modes of protein
dimers to DNA, i.e. dangling configurations, cis- and trans-binding. In
particular, the degree of compaction is extremely sensitive to the difference
in binding energies of the cis and trans configurations. Our simulations also
point out that the conformations of the DNA-protein complexes are significantly
different in bulk and in planar conditions, suggesting that conformations
observed on mica surfaces may differ significantly from those that prevail in
living cells.


  We study systems with two symmetric absorbing states, such as the voter model
and variations of it, which have been broadly used as minimal neutral models in
genetics, population ecology, sociology, etc. We analyze the effects of a key
ingredient ineluctably present in most real applications: random-field-like
quenched disorder. In accord with simulations and previous findings,
coexistence between the two competing states/opinions turns out to be strongly
favored by disorder in the standard voter model; actually, a disorder-induced
phase transition is generated for any finite system size in the presence of an
arbitrary small spontaneous-inversion rate (preventing absorbing states from
being stable). For non-linear versions of the voter model a general theory (by
AlHammal et al.) explains that the spontaneous breaking of the up/down symmetry
and an absorbing state phase transition can occur either together or
separately, giving raise to two different scenarios. Here, we show that he
presence of quenched disorder in non-linear voter models does not allow the
separation of the up-down (Ising-like) symmetry breaking from the
active-to-absorbing phase transition in low-dimensional systems: both phenomena
can occur only simultaneously, as a consequence of the well-known Imry-Ma
argument generalized to these non-equilibrium problems. When the two phenomena
occur at unison, resulting into a genuinely non-equilibrium ("Generalized
Voter") transition, the Imry-Ma argument is violated and the symmetry can be
spontaneously broken even in low dimensions.


  Chemical mapping methods probe RNA structure by revealing and leveraging
correlations of a nucleotide's structural accessibility or flexibility with its
reactivity to various chemical probes. Pioneering work by Lucks and colleagues
has expanded this method to probe hundreds of molecules at once on an Illumina
sequencing platform, obviating the use of slab gels or capillary
electrophoresis on one molecule at a time. Here, we describe optimizations to
this method from our lab, resulting in the MAP-seq protocol (Multiplexed
Accessibility Probing read out through sequencing), version 1.0. The protocol
permits the quantitative probing of thousands of RNAs at once, by several
chemical modification reagents, on the time scale of a day using a table-top
Illumina machine. This method and a software package MAPseeker
(http://simtk.org/home/map_seeker) address several potential sources of bias,
by eliminating PCR steps, improving ligation efficiencies of ssDNA adapters,
and avoiding problematic heuristics in prior algorithms. We hope that the
step-by-step description of MAP-seq 1.0 will help other RNA mapping
laboratories to transition from electrophoretic to next-generation sequencing
methods and to further reduce the turnaround time and any remaining biases of
the protocol.


  We discuss two cases that can be connected to the dynamics of interacting
populations: (I.) density waves for the case or negligible random fluctuations
of the populations densities, and (II.) probability distributions connected to
the model equations for of spatially averaged populations densities for the
case of significant random fluctuations of the independent quantity that can be
associated with the population density. For the case (I.) we consider model
equations containing polynomial nonlinearities. Such nonlinearities can arise
as a consequence of interaction among the populations (for the case of large
population densities) or as a result of a Taylor series expansion (for the case
of small density of interacting populations). In the both cases we can apply
the modified method of simplest equation to obtain exact traveling-wave
solutions connected to migration of population members. Such solutions are
obtained for systems consisting of 1 or 3 populations respectively. For the
case (II.) we discuss model equations of the Fokker-Planck kind for the
evolution of the statistical distributions of population densities. We derive
several stationary distributions for the population density and calculate the
expected exit time connected to the extinction of the studied population.


  We show that a normalized alignment free similarity measure, called D2z, can
be used to detect potential regulatory relations for gene sets when little is
known about the regulatory elements involved. One scenario where such gene sets
arise is genome wide association studies (GWAS). In this work we consider a
gene set from a GWAS on childhood general cognitive ability. We build a
co-regulation network for the GWAS genes based on the D2z scores, which shows
potential co-regulatory relationships between the genes as well as predict
additional genes that are likely to be part of the network. We found that the
set of the predicted genes is enriched in genes associated with mental
retardation and GO terms such as synapse and neuron development. In particular,
we found strong evidence of regulatory connection between the GWAS genes and
CHL1, a gene known to be involved in mental retardation.


  Graph rewrite systems are powerful tools to model and study complex problems
in various fields of research. Their successful application to chemical
reaction modelling on a molecular level was shown but no appropriate and simple
system is available at the moment.
  The presented Graph Grammar Library (GGL) implements a generic Double Push
Out approach for general graph rewrite systems. The framework focuses on a high
level of modularity as well as high performance, using state-of-the-art
algorithms and data structures, and comes with extensive documentation. The
large GGL chemistry module enables extensive and detailed studies of chemical
systems. It well meets the requirements and abilities envisioned by Yadav et
al. (2004) for such chemical rewrite systems. Here, molecules are represented
as undirected labeled graphs while chemical reactions are described by
according graph grammar rules. Beside the graph transformation, the GGL offers
advanced cheminformatics algorithms for instance to estimate energies
ofmolecules or aromaticity perception. These features are illustrated using a
set of reactions from polyketide chemistry a huge class of natural compounds of
medical relevance.
  The graph grammar based simulation of chemical reactions offered by the GGL
is a powerful tool for extensive cheminformatics studies on a molecular level.
The GGL already provides rewrite rules for all enzymes listed in the KEGG
LIGAND database is freely available at
http://www.tbi.univie.ac.at/software/GGL/.


  The proper biological functioning of proteins often relies on the occurrence
of coordinated fluctuations around their native structure, or of wider and
sometimes highly elaborated motions. Coarse-grained elastic-network
descriptions are known to capture essential aspects of conformational dynamics
in proteins, but have so far remained mostly phenomenological, and unable to
account for the chemical specificities of amino acids. Here, we propose a
method to derive residue- and distance-specific effective harmonic potentials
from the statistical analysis of an extensive dataset of NMR conformational
ensembles. These potentials constitute dynamical counterparts to the mean-force
statistical potentials commonly used for static analyses of protein structures.
In the context of the elastic network model, they yield a strongly improved
description of the cooperative aspects of residue motions, and give the
opportunity to systematically explore the influence of sequence details on
protein dynamics.


  Motivated by empirical observations of algebraic duplicated sequence length
distributions in a broad range of natural genomes, we analytically formulate
and solve a class of simple discrete duplication/substitution models that
generate steady-states sharing this property. Continuum equations are derived
for arbitrary time-independent duplication length source distribution, a limit
that we show can be mapped directly onto certain fragmentation models that have
been intensively studied by physicists in recent years. Quantitative agreement
with simulation is demonstrated. These models account for the algebraic form
and exponent of naturally occuring duplication length distributions without the
need for fine-tuning of parameters.


  Most of the common used models of epidemic spreading allow contaminating many
neighbors of a particular node in the network. They are usually analyzed by
differential equations on probability vectors. We propose a model of epidemic
spreading, which restricts to at most one contamination per time step and
analyze it by discrete approach, working on vectors of possible states of the
system. Theoretical predictions of epidemic treshold, stationary state and time
needed to reach it are given and appear to be perfectly consistent with
computer simulations. We also point ou that the model appears to be well suited
to mimic epidemic spreading within student communes.


  An allometric height-mass exponent $\gamma$ gives an approximative power-law
relation $< M> \propto H^\gamma$ between the average mass $< M>$ and the height
$H$, for a sample of individuals. The individuals in the present study are
humans but could be any biological organism. The sampling can be for a specific
age of the individuals or for an age-interval. The body-mass index (BMI) is
often used for practical purposes when characterizing humans and it is based on
the allometric exponent $\gamma=2$. It is here shown that the actual value of
$\gamma$ is to large extent determined by the degree of correlation between
mass and height within the sample studied: no correlation between mass and
height means $\gamma=0$, whereas if there was a precise relation between mass
and height such that all individuals had the same shape and density then
$\gamma=3$. The connection is demonstrated by showing that the value of
$\gamma$ can be obtained directly from three numbers characterizing the spreads
of the relevant random Gaussian statistical distributions: the spread of the
height and mass distributions together with the spread of the mass distribution
for the average height. Possible implications for allometric relations in
general are discussed.


  Rapidly growing public gene expression databases contain a wealth of data for
building an unprecedentedly detailed picture of human biology and disease. This
data comes from many diverse measurement platforms that make integrating it all
difficult. Although RNA-sequencing (RNA-seq) is attracting the most attention,
at present the rate of new microarray studies submitted to public databases far
exceeds the rate of new RNA-seq studies. There is clearly a need for methods
that make it easier to combine data from different technologies. In this paper,
we propose a new method for processing RNA-seq data that yields gene expression
estimates that are much more similar to corresponding estimates from microarray
data, hence greatly improving cross-platform comparability. The method we call
PREBS is based on estimating the expression from RNA-seq reads overlapping the
microarray probe regions, and processing these estimates with microarray
summarisation algorithm RMA. Using paired microarray and RNA-seq samples from
TCGA LAML data set we show that PREBS expression estimates derived from RNA-seq
are more similar to microarray-based expression estimates than those from other
RNA-seq processing methods. In an experiment to retrieve paired microarray
samples from a database using an RNA-seq query sample, gene signatures defined
based on PREBS expression estimates were found to be much more accurate than
those from other methods. PREBS also allows new ways of using RNA-seq data,
such as expression estimation for microarray probe sets. An implementation of
the proposed method is available in the Bioconductor package `prebs'.


  Biological competition is widely believed to result in the evolution of
selfish preferences. The related concept of the `homo economicus' is at the
core of mainstream economics. However, there is also experimental and empirical
evidence for other-regarding preferences. Here we present a theory that
explains both, self-regarding and other-regarding preferences. Assuming
conditions promoting non-cooperative behaviour, we demonstrate that
intergenerational migration determines whether evolutionary competition results
in a `homo economicus' (showing self-regarding preferences) or a `homo
socialis' (having other-regarding preferences). Our model assumes spatially
interacting agents playing prisoner's dilemmas, who inherit a trait determining
`friendliness', but mutations tend to undermine it. Reproduction is ruled by
fitness-based selection without a cultural modification of reproduction rates.
Our model calls for a complementary economic theory for `networked minds' (the
`homo socialis') and lays the foundations for an evolutionarily grounded theory
of other-regarding agents, explaining individually different utility functions
as well as conditional cooperation.


  Adaptive social networks, in which nodes and network structure co-evolve, are
often described using a mean-field system of equations for the density of node
and link types. These equations constitute an open system due to dependence on
higher order topological structures. We propose a new approach to moment
closure based on the analytical description of the system in an asymptotic
regime. We apply the proposed approach to two examples of adaptive networks:
recruitment to a cause model and adaptive epidemic model. We show a good
agreement between the improved mean-field prediction and simulations of the
full network system.


  We present a two-layer hidden Markov model to detect structure of haplotypes
for unrelated individuals. This allows modeling two scales of linkage
disequilibrium (one within a group of haplotypes and one between groups),
thereby taking advantage of rich haplotype information to infer local ancestry
for admixed individuals. Our method outperforms competing state-of-art methods,
particularly for regions of small ancestral track lengths. Applying our method
to Mexican samples in HapMap3, we found five coding regions, ranging from $0.3
-1.3$ megabase (Mb) in lengths, that exhibit excessive European ancestry
(average dosage > 1.6). A particular interesting region of 1.1Mb (with average
dosage 1.95) locates on Chromosome 2p23 that harbors two genes, PXDN and MYT1L,
both of which are associated with autism and schizophrenia. In light of the low
prevalence of autism in Hispanics, this region warrants special attention. We
confirmed our findings using Mexican samples from the 1000 genomes project. A
software package implementing methods described in the paper is freely
available at \url{http://bcm.edu/cnrc/mcmcmc}


  In order to formulate the Fundamental Theorem of Natural Selection, Fisher
defined the average excess and average effect of a gene substitution. Finding
these notions to be somewhat opaque, some authors have recommended
reformulating Fisher's ideas in terms of covariance and regression, which are
classical concepts of statistics. We argue that Fisher intended his two
averages to express a distinction between correlation and causation. On this
view the average effect is a specific weighted average of the actual phenotypic
changes that result from physically changing the allelic states of homologous
genes. We show that the statistical and causal conceptions of the average
effect, perceived as inconsistent by Falconer, can be reconciled if certain
relationships between the genotype frequencies and non-additive residuals are
conserved. There are certain theory-internal considerations favoring Fisher's
original formulation in terms of causality; for example, the frequency-weighted
mean of the average effects equaling zero at each locus becomes a derivable
consequence rather than an arbitrary constraint. More broadly, Fisher's
distinction between correlation and causation is of critical importance to
gene-trait mapping studies and the foundations of evolutionary biology.


  Individuals within any species exhibit differences in size, developmental
state, or spatial location. These differences coupled with environmental
fluctuations in demographic rates can have subtle effects on population
persistence and species coexistence. To understand these effects, we provide a
general theory for coexistence of structured, interacting species living in a
stochastic environment. The theory is applicable to nonlinear, multi species
matrix models with stochastically varying parameters. The theory relies on
long-term growth rates of species corresponding to the dominant Lyapunov
exponents of random matrix products. Our coexistence criterion requires that a
convex combination of these long-term growth rates is positive with probability
one whenever one or more species are at low density. When this condition holds,
the community is stochastically persistent: the fraction of time that a species
density goes below $\delta>0$ approaches zero as $\delta$ approaches zero.
Applications to predator-prey interactions in an autocorrelated environment, a
stochastic LPA model, and spatial lottery models are provided. These
applications demonstrate that positive autocorrelations in temporal
fluctuations can disrupt predator-prey coexistence, fluctuations in
log-fecundity can facilitate persistence in structured populations, and
long-lived, relatively sedentary competing populations are likely to coexist in
spatially and temporally heterogenous environments.


  Proofreading/editing in protein synthesis is essential for accurate
translation of information from the genetic code. In this article we present a
theoretical investigation of efficiency of a kinetic proofreading mechanism
that employs hydrolysis of the wrong substrate as the discriminatory step in
enzyme catalytic reactions. We consider aminoacylation of tRNA^{Ile} which is a
crucial step in protein synthesis and for which experimental results are now
available. We present an augmented kinetic scheme and then employ methods of
stochastic simulation algorithm to obtain time dependent concentrations of
different substances involved in the reaction and their rates of formation. We
obtain the rates of product formation and ATP hydrolysis for both correct and
wrong substrates (isoleucine and valine in our case), in single molecular
enzyme as well as ensemble enzyme kinetics. The present theoretical scheme
correctly reproduces (i) the amplitude of the discrimination factor in the
overall rates between isoleucine and valine which is obtained as (1.8 \times
10^2).(4.33 \times 10^2) = 7.8 \times 10^4, (ii) the rates of ATP hydrolysis
for both Ile and Val at different substrate concentrations in the
aminoacylation of tRNA^{Ile}. The present study shows a non-michaelis type
dependence of rate of reaction on tRNA^{Ile} concentration in case of valine.
The editing in steady state is found to be independent of amino acid
concentration. Interestingly, the computed ATP hydrolysis rate for valine at
high substrate concentration is same as the rate of formation of Ile-tRNA^{Ile}
whereas at intermediate substrate concentration the ATP hydrolysis rate is
relatively low.


  In March and early April 2013, a new avian-origin influenza A (H7N9) virus
(A-OIV) emerged in the eastern China. During the first week of April, the 18
infection cases have been confirmed and 6 people have died since March. This
virus has caused global concern as a potential pandemic threat. Here we use
evolutionary analysis to reconstruct the origins and early development of the
A-OIV viruses. We found that A-OIV was derived from a reassortment of three
avian flu virus strains, and substantial mutations have been detected. Our
results highlight the need for systematic surveillance of influenza in birds,
and provide evidence that the mixing of new genetic elements in avian can
result in the emergence of viruses with pandemic potential in humans.


  The process controlling the diferentiation of stem, or progenitor, cells into
one specific functional direction is called lineage specification. An important
characteristic of this process is the multi-lineage priming, which requires the
simultaneous expression of lineage-specific genes. Prior to commitment to a
certain lineage, it has been observed that these genes exhibit intermediate
values of their expression levels. Multi-lineage differentiation has been
reported for various progenitor cells, and it has been explained through the
bifurcation of a metastable state. During the differentiation process the
dynamics of the core regulatory network follows a bifurcation, where the
metastable state, corresponding to the progenitor cell, is destabilized and the
system is forced to choose between the possible developmental alternatives.
While this approach gives a reasonable interpretation of the cell fate decision
process, it fails to explain the multi-lineage priming characteristic. Here, we
describe a new multi-dimensional switch-like model that captures both the
process of cell fate decision and the phenomenon of multi-lineage priming. We
show that in the symmetrical interaction case, the system exhibits a new type
of degenerate bifurcation, characterized by a critical hyperplane, containing
an infinite number of critical steady states. This critical hyperplane may be
interpreted as the support for the multi-lineage priming states of the
progenitor. Also, the cell fate decision (the multi-stability and switching
behavior) can be explained by a symmetry breaking in the parameter space of
this critical hyperplane. These analytical results are confirmed by Monte-Carlo
simulations of the corresponding chemical master equations.