  Cooperative behavior, where one individual incurs a cost to help another, is
a wide spread phenomenon. Here we study direct reciprocity in the context of
the alternating Prisoner's Dilemma. We consider all strategies that can be
implemented by one and two-state automata. We calculate the payoff matrix of
all pairwise encounters in the presence of noise. We explore deterministic
selection dynamics with and without mutation. Using different error rates and
payoff values, we observe convergence to a small number of distinct equilibria.
Two of them are uncooperative strict Nash equilibria representing always-defect
(ALLD) and Grim. The third equilibrium is mixed and represents a cooperative
alliance of several strategies, dominated by a strategy which we call Forgiver.
Forgiver cooperates whenever the opponent has cooperated; it defects once when
the opponent has defected, but subsequently Forgiver attempts to re-establish
cooperation even if the opponent has defected again. Forgiver is not an
evolutionarily stable strategy, but the alliance, which it rules, is
asymptotically stable. For a wide range of parameter values the most commonly
observed outcome is convergence to the mixed equilibrium, dominated by
Forgiver. Our results show that although forgiving might incur a short-term
loss it can lead to a long-term gain. Forgiveness facilitates stable
cooperation in the presence of exploitation and noise.


  Despite the importance of a thermodynamically stable structure with a
conserved fold for protein function, almost all evolutionary models neglect
site-site correlations that arise from physical interactions between
neighboring amino acid sites. This is mainly due to the difficulty in
formulating a computationally tractable model since rate matrices can no longer
be used. Here we introduce a general framework, based on factor graphs, for
constructing probabilistic models of protein evolution with site
interdependence. Conveniently, efficient approximate inference algorithms, like
Belief Propagation, can be used to calculate likelihoods for these models. We
fit an amino acid substitution model of this type that accounts for both
solvent accessibility and site-site correlations. Comparisons of the new model
with rate matrix models and a model accounting only for solvent accessibility
demonstrate that it better fits the sequence data. We also examine evolution
within a family of homohexameric enzymes and find that site-site correlations
between most contacting subunits contribute to a higher likelihood. In
addition, we show that the new substitution model has a similar mathematical
form to the one introduced in (Rodrigue et al. 2005), although with different
parameter interpretations and values. We also perform a statistical analysis of
the effects of amino acids at neighboring sites on substitution probabilities
and find a significant perturbation of most probabilities, further supporting
the significant role of site-site interactions in protein evolution and
motivating the development of new evolutionary models like the one described
here. Finally, we discuss possible extensions and applications of the new
substitution models.


  We present an experimentally guided, bi-phasic, multi-species ionic gel model
to compare and make quantitative predictions on the viscoelastic properties of
healthy mucus versus Cystic-Fibrosis infected mucus. The mixture theory
consists of the mucus (polymer phase) and water (solvent phase) as well as
several different ions: H$^+$, Na$^+$, Cl$^-$ and Ca$^{2+}$. The model is
calibrated with the experimental data of mucus diffusivities. Numerical studies
reveal that the Donnan forces are the dominating mechanism driving the mucus
swelling/deswelling transition. However, these transitions could be continuous
or discontinuous, depending on the complex interplay between ionic bath
concentrations, pH, polymer mass in the solvent and the average charge per
monomer.


  Containing the spreading of crime is a major challenge for society. Yet,
since thousands of years, no effective strategy has been found to overcome
crime. To the contrary, empirical evidence shows that crime is recurrent, a
fact that is not captured well by rational choice theories of crime. According
to these, strong enough punishment should prevent crime from happening. To gain
a better understanding of the relationship between crime and punishment, we
consider that the latter requires prior discovery of illicit behavior and study
a spatial version of the inspection game. Simulations reveal the spontaneous
emergence of cyclic dominance between ''criminals'', ''inspectors'', and
''ordinary people'' as a consequence of spatial interactions. Such cycles
dominate the evolutionary process, in particular when the temptation to commit
crime or the cost of inspection are low or moderate. Yet, there are also
critical parameter values beyond which cycles cease to exist and the population
is dominated either by a stable mixture of criminals and inspectors or one of
these two strategies alone. Both continuous and discontinuous phase transitions
to different final states are possible, indicating that successful strategies
to contain crime can be very much counter-intuitive and complex. Our results
demonstrate that spatial interactions are crucial for the evolutionary outcome
of the inspection game, and they also reveal why criminal behavior is likely to
be recurrent rather than evolving towards an equilibrium with monotonous
parameter dependencies.


  Computational methods are needed to differentiate the small fraction of
missense mutations that contribute to disease by disrupting protein function
from neutral variants. We describe several complementary methods using
large-scale homology modeling of human protein complexes to detect non-neutral
mutations. Importantly, unlike sequence conservation-based methods, this
structure-based approach provides experimentally testable biochemical
mechanisms for mutations in disease. Specifically, we infer metal ion, small
molecule, protein-protein, and nucleic acid binding sites by homology and find
that disease-associated missense mutations are more prevalent in each class of
binding site than are neutral mutations. Importantly, our approach identifies
considerably more binding sites than those annotated in the RefSeq database.
Furthermore, an analysis of metal ion and protein-protein binding sites
predicted by machine learning shows a similar preponderance of
disease-associated mutations in these sites. We also derive a statistical score
for predicting how mutations affect metal ion binding and find many dbSNP
mutations that likely disrupt ion binding but were not previously considered
deleterious. We also cluster mutations in the protein structure to discover
putative functional regions. Finally, we develop a machine learning predictor
for detecting disease-associated missense mutations and show that it
outperforms two other prediction methods on an independent test set.


  Artificial nanowalkers are inspired by biomolecular counterparts from living
cells, but remain far from comparable to the latter in design principles. The
walkers reported to date mostly rely on chemical mechanisms to gain a
direction; they all produce chemical wastes. Here we report a light-powered DNA
bipedal walker based on a design principle derived from cellular walkers. The
walker has two identical feet and the track has equal binding sites; yet the
walker gains a direction by pure physical mechanisms that autonomously amplify
an intra-site asymmetry into a ratchet effect. The nanowalker is free of any
chemical waste. It has a distinct thermodynamic feature that it possesses the
same equilibrium before and after operation, but generates a truly
non-equilibrium distribution during operation. The demonstrated design
principle exploits mechanical effects and is adaptable for use in other
nanomachines.


  Metapopulation (multipatch) models are widely used to study the patterns of
spatial spread of epidemics. In this paper we study the impact of inter-patch
connection weights on the predictions of these models. We contrast arbitrary,
uniform link weights with link weights predicted using a gravity model based on
patch populations and distance. In a synthetic system with one large driver
city and many small follower cities, we show that under uniform link weights,
epidemics in the follower regions are perfectly synchronized. In contrast,
gravity-based links allow a more realistic, less synchronized distribution of
epidemic peaks in the follower regions. We then fit a three-patch
metapopulation model to regional dengue fever data from Peru -- a country
experiencing yearly, spatially defined epidemics. We use data for 2002-2008
(studying the seasonal disease patterns in the country and the yearly
reinfection patterns from jungle to the coast) and 2000-2001 (one large
epidemic of a new disease strain across the country). We present numerical
results.


  Models of biochemical networks are frequently high-dimensional and complex.
Reduction methods that preserve important dynamical properties are therefore
essential in their study. Interactions between the nodes in such networks are
frequently modeled using a Hill function, $x^n/(J^n+x^n)$. Reduced ODEs and
Boolean networks have been studied extensively when the exponent $n$ is large.
However, the case of small constant $J$ appears in practice, but is not well
understood. In this paper we provide a mathematical analysis of this limit, and
show that a reduction to a set of piecewise linear ODEs and Boolean networks
can be mathematically justified. The piecewise linear systems have closed form
solutions that closely track those of the fully nonlinear model. On the other
hand, the simpler, Boolean network can be used to study the qualitative
behavior of the original system. We justify the reduction using geometric
singular perturbation theory and compact convergence, and illustrate the
results in networks modeling a genetic switch and a genetic oscillator.


  In a combinatorial communication system, some signals consist of the
combinations of other signals. Such systems are more efficient than equivalent,
non-combinatorial systems, yet despite this they are rare in nature. Why?
Previous explanations have focused on the adaptive limits of combinatorial
communication, or on its purported cognitive difficulties, but neither of these
explains the full distribution of combinatorial communication in the natural
world. Here we present a nonlinear dynamical model of the emergence of
combinatorial communication that, unlike previous models, considers how
initially non-communicative behaviour evolves to take on a communicative
function. We derive three basic principles about the emergence of combinatorial
communication. We hence show that the interdependence of signals and responses
places significant constraints on the historical pathways by which
combinatorial signals might emerge, to the extent that anything other than the
most simple form of combinatorial communication is extremely unlikely. We also
argue that these constraints can be bypassed if individuals have the
socio-cognitive capacity to engage in ostensive communication. Humans, but
probably no other species, have this ability. This may explain why language,
which is massively combinatorial, is such an extreme exception to nature's
general trend for non-combinatorial communication.


  Mutualistic networks are formed when the interactions between two classes of
species are mutually beneficial. They are important examples of cooperation
shaped by evolution. Mutualism between animals and plants plays a key role in
the organization of ecological communities. Such networks in ecology have
generically evolved a nested architecture independent of species composition
and latitude - specialists interact with proper subsets of the nodes with whom
generalists interact. Despite sustained efforts to explain observed network
structure on the basis of community-level stability or persistence, such
correlative studies have reached minimal consensus. Here we demonstrate that
nested interaction networks could emerge as a consequence of an optimization
principle aimed at maximizing the species abundance in mutualistic communities.
Using analytical and numerical approaches, we show that because of the
mutualistic interactions, an increase in abundance of a given species results
in a corresponding increase in the total number of individuals in the
community, as also the nestedness of the interaction matrix. Indeed, the
species abundances and the nestedness of the interaction matrix are correlated
by an amount that depends on the strength of the mutualistic interactions.
Nestedness and the observed spontaneous emergence of generalist and specialist
species occur for several dynamical implementations of the variational
principle under stationary conditions. Optimized networks, while remaining
stable, tend to be less resilient than their counterparts with randomly
assigned interactions. In particular, we analytically show that the abundance
of the rarest species is directly linked to the resilience of the community.
Our work provides a unifying framework for studying the emergent structural and
dynamical properties of ecological mutualistic networks.


  In this study, we investigate interactions of extended conformations of
homodimeric peptides made of small (glycine or alanine) and large hydrophobic
(valine or leucine) sidechains using all-atom molecular dynamics simulations to
decipher driving forces for \beta-sheet formation. We make use of a periodic
boundary condition setup in which individual peptides are infinitely long and
stretched. Dimers adopt \beta-sheet conformations at short interpeptide
distances (\xi ~ 0.5 nm) and at intermediate distances (~ 0.8 nm), valine and
leucine homodimers assume cross-\beta-like conformations with side chains
interpenetrating each other. These two states are identified as minima in the
Potential of Mean Force (PMF). While the number of interpeptide hydrogen bonds
increases with decreasing interpeptide distance, the total hydrogen bond number
in the system does not change significantly, suggesting that formation of
\beta-sheet structures from extended conformations is not driven by hydrogen
bonds. This is supported by an increase in electrostatic energy at short
interpeptide distances. A remarkable correlation between the volume of the
system and the total electrostatic energy is observed, further reinforcing the
idea that excluding water in proteins comes with an enthalpic penalty. We also
discuss microscopic mechanisms accounting for \beta-sheet formation based on
computed enthalpy and entropy and we show that they are different for peptides
with small and large side chains.


  An ongoing challenge in the mathematical modeling of biological aggregations
is to strengthen the connection between models and biological data by
quantifying the rules that individuals follow. We model aggregation of the pea
aphid, Acyrthosiphon pisum. Specifically, we conduct experiments to track the
motion of aphids walking in a featureless circular arena in order to deduce
individual-level rules. We observe that each aphid transitions stochastically
between a moving and a stationary state. Moving aphids follow a correlated
random walk. The probabilities of motion state transitions, as well as the
random walk parameters, depend strongly on distance to an aphid's nearest
neighbor. For large nearest neighbor distances, when an aphid is essentially
isolated, its motion is ballistic with aphids moving faster, turning less, and
being less likely to stop. In contrast, for short nearest neighbor distances,
aphids move more slowly, turn more, and are more likely to become stationary;
this behavior constitutes an aggregation mechanism. From the experimental data,
we estimate the state transition probabilities and correlated random walk
parameters as a function of nearest neighbor distance. With the
individual-level model established, we assess whether it reproduces the
macroscopic patterns of movement at the group level. To do so, we consider
three distributions, namely distance to nearest neighbor, angle to nearest
neighbor, and percentage of population moving at any given time. For each of
these three distributions, we compare our experimental data to the output of
numerical simulations of our nearest neighbor model, and of a control model in
which aphids do not interact socially. Our stochastic, social nearest neighbor
model reproduces salient features of the experimental data that are not
captured by the control.


  The role of positive selection in human evolution remains controversial. On
the one hand, scans for positive selection have identified hundreds of
candidate loci and the genome-wide patterns of polymorphism show signatures
consistent with frequent positive selection. On the other hand, recent studies
have argued that many of the candidate loci are false positives and that most
apparent genome-wide signatures of adaptation are in fact due to reduction of
neutral diversity by linked recurrent deleterious mutations, known as
background selection. Here we analyze human polymorphism data from the 1,000
Genomes project (Abecasis et al. 2012) and detect signatures of pervasive
positive selection once we correct for the effects of background selection. We
show that levels of neutral polymorphism are lower near amino acid
substitutions, with the strongest reduction observed specifically near
functionally consequential amino acid substitutions. Furthermore, amino acid
substitutions are associated with signatures of recent adaptation that should
not be generated by background selection, such as the presence of unusually
long and frequent haplotypes and specific distortions in the site frequency
spectrum. We use forward simulations to show that the observed signatures
require a high rate of strongly adaptive substitutions in the vicinity of the
amino acid changes. We further demonstrate that the observed signatures of
positive selection correlate more strongly with the presence of regulatory
sequences, as predicted by ENCODE (Gerstein et al. 2012), than the positions of
amino acid substitutions. Our results establish that adaptation was frequent in
human evolution and provide support for the hypothesis of King and Wilson (King
and Wilson 1975) that adaptive divergence is primarily driven by regulatory
changes.


  Recent research has identified interactions between networks as crucial for
the outcome of evolutionary games taking place on them. While the consensus is
that interdependence does promote cooperation by means of organizational
complexity and enhanced reciprocity that is out of reach on isolated networks,
we here address the question just how much interdependence there should be.
Intuitively, one might assume the more the better. However, we show that in
fact only an intermediate density of sufficiently strong interactions between
networks warrants an optimal resolution of social dilemmas. This is due to an
intricate interplay between the heterogeneity that causes an asymmetric
strategy flow because of the additional links between the networks, and the
independent formation of cooperative patterns on each individual network.
Presented results are robust to variations of the strategy updating rule, the
topology of interdependent networks, and the governing social dilemma, thus
suggesting a high degree of universality.


  A fitness landscape is a genetic space -- with two genotypes adjacent if they
differ in a single locus -- and a fitness function. Evolutionary dynamics
produce a flow on this landscape from lower fitness to higher; reaching
equilibrium only if a local fitness peak is found. I use computational
complexity to question the common assumption that evolution on static fitness
landscapes can quickly reach a local fitness peak. I do this by showing that
the popular NK model of rugged fitness landscapes is PLS-complete for K >= 2;
the reduction from Weighted 2SAT is a bijection on adaptive walks, so there are
NK fitness landscapes where every adaptive path from some vertices is of
exponential length. Alternatively -- under the standard complexity theoretic
assumption that there are problems in PLS not solvable in polynomial time --
this means that there are no evolutionary dynamics (known, or to be discovered,
and not necessarily following adaptive paths) that can converge to a local
fitness peak on all NK landscapes with K = 2. Applying results from the
analysis of simplex algorithms, I show that there exist single-peaked
landscapes with no reciprocal sign epistasis where the expected length of an
adaptive path following strong selection weak mutation dynamics is
$e^{O(n^{1/3})}$ even though an adaptive path to the optimum of length less
than n is available from every vertex. The technical results are written to be
accessible to mathematical biologists without a computer science background,
and the biological literature is summarized for the convenience of
non-biologists with the aim to open a constructive dialogue between the two
disciplines.


  Protein aggregation is an important field of investigation because it is
closely related to the problem of neurodegenerative diseases, to the
development of biomaterials, and to the growth of cellular structures such as
cyto-skeleton. Self-aggregation of protein amyloids, for example, is a
complicated process involving many species and levels of structures. This
complexity, however, can be dealt with using statistical mechanical tools, such
as free energies, partition functions, and transfer matrices. In this article,
we review general strategies for studying protein aggregation using statistical
mechanical approaches and show that canonical and grand canonical ensembles can
be used in such approaches. The grand canonical approach is particularly
convenient since competing pathways of assembly and dis-assembly can be
considered simultaneously. Another advantage of using statistical mechanics is
that numerically exact solutions can be obtained for all of the thermodynamic
properties of fibrils, such as the amount of fibrils formed, as a function of
initial protein concentration. Furthermore, statistical mechanics models can be
used to fit experimental data when they are available for comparison.


  The evolution of RNA viruses such as HIV, Hepatitis C and Influenza virus
occurs so rapidly that the viruses' genomes contain information on past
ecological dynamics. Hence, we develop a phylodynamic method that enables the
joint estimation of epidemiological parameters and phylogenetic history. Based
on a compartmental susceptible-infected-removed (SIR) model, this method
provides separate information on incidence and prevalence of infections.
Detailed information on the interaction of host population dynamics and
evolutionary history can inform decisions on how to contain or entirely avoid
disease outbreaks.
  We apply our Birth-Death SIR method (BDSIR) to two viral data sets. First,
five human immunodeficiency virus type 1 clusters sampled in the United Kingdom
between 1999 and 2003 are analyzed. The estimated basic reproduction ratios
range from 1.9 to 3.2 among the clusters. All clusters show a decline in the
growth rate of the local epidemic in the middle or end of the 90's.
  The analysis of a hepatitis C virus (HCV) genotype 2c data set shows that the
local epidemic in the C\'ordoban city Cruz del Eje originated around 1906
(median), coinciding with an immigration wave from Europe to central Argentina
that dates from 1880--1920. The estimated time of epidemic peak is around 1970.


  We propose a kinetic model for the self-aggregation by amyloid proteins. By
extending several well-known models for protein aggregation, the time evolution
of aggregate concentrations containing $r$ proteins, denoted $c_r(t)$, can be
written in terms of generalized Smoluchowski kinetics. With this approach we
take into account all possible aggregation and fragmentation reactions
involving clusters of any size. Correspondingly, an aggregate of size x+y could
be formed by or break-up into two smaller constituent aggregates of sizes x and
y. The rates of each aggregation or fragmentation reaction, called kernels, are
specified in terms of the aggregate size, and we solve $c_r(t)$ for large
cluster sizes using numerical techniques. We show that by using Smoluchowski
kinetics many pathways to fibrillation are possible and quantities, such as the
aggregate length distribution at an arbitrary time, can be calculated. We show
that the predicted results of the model are in agreement with the experimental
observations.


  We present a new experimental-computational technology of inferring network
models that predict the response of cells to perturbations and that may be
useful in the design of combinatorial therapy against cancer. The experiments
are systematic series of perturbations of cancer cell lines by targeted drugs,
singly or in combination. The response to perturbation is measured in terms of
levels of proteins and phospho-proteins and of cellular phenotype such as
viability. Computational network models are derived de novo, i.e., without
prior knowledge of signaling pathways, and are based on simple non-linear
differential equations. The prohibitively large solution space of all possible
network models is explored efficiently using a probabilistic algorithm, belief
propagation, which is three orders of magnitude more efficient than Monte Carlo
methods. Explicit executable models are derived for a set of perturbation
experiments in Skmel-133 melanoma cell lines, which are resistant to the
therapeutically important inhibition of Raf kinase. The resulting network
models reproduce and extend known pathway biology. They can be applied to
discover new molecular interactions and to predict the effect of novel drug
perturbations, one of which is verified experimentally. The technology is
suitable for application to larger systems in diverse areas of molecular
biology.


  The var genes of the human malaria parasite Plasmodium falciparum present a
challenge to population geneticists due to their extreme diversity, which is
generated by high rates of recombination. These genes encode a primary antigen
protein called PfEMP1, which is expressed on the surface of infected red blood
cells and elicits protective immune responses. Var gene sequences are
characterized by pronounced mosaicism, precluding the use of traditional
phylogenetic tools that require bifurcating tree-like evolutionary
relationships. We present a new method that identifies highly variable regions
(HVRs), and then maps each HVR to a complex network in which each sequence is a
node and two nodes are linked if they share an exact match of significant
length. Here, networks of var genes that recombine freely are expected to have
a uniformly random structure, but constraints on recombination will produce
network communities that we identify using a stochastic block model. We
validate this method on synthetic data, showing that it correctly recovers
populations of constrained recombination, before applying it to the Duffy
Binding Like-{\alpha} (DBL{\alpha}) domain of var genes. We find nine HVRs
whose network communities map in distinctive ways to known DBL{\alpha}
classifications and clinical phenotypes. We show that the recombinational
constraints of some HVRs are correlated, while others are independent. These
findings suggest that this micromodular structuring facilitates independent
evolutionary trajectories of neighboring mosaic regions, allowing the parasite
to retain protein function while generating enormous sequence diversity. Our
approach therefore offers a rigorous method for analyzing evolutionary
constraints in var genes, and is also flexible enough to be easily applied more
generally to any highly recombinant sequences.


  More than any other species, humans form social ties to individuals who are
neither kin nor mates, and these ties tend to be with similar people. Here, we
show that this similarity extends to genotypes. Across the whole genome,
friends' genotypes at the SNP level tend to be positively correlated
(homophilic); however, certain genotypes are negatively correlated
(heterophilic). A focused gene set analysis suggests that some of the overall
correlation can be explained by specific systems; for example, an olfactory
gene set is homophilic and an immune system gene set is heterophilic. Finally,
homophilic genotypes exhibit significantly higher measures of positive
selection, suggesting that, on average, they may yield a synergistic fitness
advantage that has been helping to drive recent human evolution.


  Chemical crosslinking of collagen is a general strategy to reproduce
macroscale tissue properties in physiological environment. However,
simultaneous control of protein conformation, material properties and
biofunctionality is highly challenging with current synthetic strategies.
Consequently, the potentially-diverse clinical applications of collagen-based
biomaterials cannot be fully realised. In order to establish defined
biomacromolecular systems for mineralised tissue applications, type I collagen
was functionalised with 1,3-Phenylenediacetic acid (Ph) and investigated at the
molecular, macroscopic and functional levels. Preserved triple helix
conformation was observed in obtained covalent networks via ATR-FTIR and WAXS,
while network crosslinking degree could be adjusted based on specific reaction
conditions. Decreased swelling ratio and increased thermo-mechanical properties
were observed compared to state of-the-art carbodiimide (EDC)-crosslinked
collagen controls, likely related to the intermolecular covalent incorporation
of the aromatic segment. Ph crosslinked hydrogels displayed nearly intact
material integrity and only a slight mass decrease following 1-week incubation
in either PBS or simulated body fluid (SBF), in contrast to EDC crosslinked
collagen. Furthermore, FTIR, SEM and EDS revealed deposition of a
calcium-phosphate phase on SBF-retrieved samples, whereby an increased calcium
phosphate ratio was observed in hydrogels with higher Ph content. 72-hour
material extracts were well tolerated by L929 mouse fibroblasts, whereby cell
confluence and metabolic activity (MTS assay) were comparable to those of cells
cultured in cell culture medium (positive control). In light of their
controlled structure-function properties, these biocompatible collagen
hydrogels represent attractive material systems for potential mineralised
tissue applications.


  We present a hybrid model of a biological filter, a genetic circuit which
removes fast fluctuations in the cell's internal representation of the extra
cellular environment. The model takes the classic feed-forward loop (FFL) motif
and represents it as a network of continuous protein concentrations and binary,
unobserved gene promoter states. We address the problem of statistical
inference and parameter learning for this class of models from partial,
discrete time observations. We show that the hybrid representation leads to an
efficient algorithm for approximate statistical inference in this circuit, and
show its effectiveness on a simulated data set.


  Comparative genetic studies of non-model organisms are transforming rapidly
due to major advances in sequencing technology. A limiting factor in these
studies has been the identification and screening of orthologous loci across an
evolutionarily distant set of taxa. Here, we evaluate the efficacy of genomic
markers targeting ultraconserved DNA elements (UCEs) for analyses at shallow
evolutionary timescales. Using sequence capture and massively parallel
sequencing to generate UCE data for five co-distributed Neotropical rainforest
bird species, we recovered 776-1,516 UCE loci across the five species. Across
species, 53-77 percent of the loci were polymorphic, containing between 2.0 and
3.2 variable sites per polymorphic locus, on average. We performed species tree
construction, coalescent modeling, and species delimitation, and we found that
the five co-distributed species exhibited discordant phylogeographic histories.
We also found that species trees and divergence times estimated from UCEs were
similar to those obtained from mtDNA. The species that inhabit the understory
had older divergence times across barriers, contained a higher number of
cryptic species, and exhibited larger effective population sizes relative to
species inhabiting the canopy. Because orthologous UCEs can be obtained from a
wide array of taxa, are polymorphic at shallow evolutionary time scales, and
can be generated rapidly at low cost, they are effective genetic markers for
studies investigating evolutionary patterns and processes at shallow time
scales.


  We develop a model to describe the development of dominance relations between
social animals as they use past experiences to inform future interactions.
Using the game-theoretic framework of a Hawk-Dove game with asymmetric
resource-holding potentials (RHPs), we derive a simple model that describes the
social interactions of animals that compete for resources. We then consider a
game-playing strategy in which animals acquire information about their RHP
asymmetry from the results of their previous contests and subsequently use
their asymmetry assessment to inform their behavior in future contests. We
examine how directly incorporating the fact that animals have incomplete
information in their interactions can lead to outcomes that differ from what
would be expected if one considers the situation as a static game in which the
animals have perfect information about the asymmetry size. We thereby obtain
results that are consistent with observations of dominance-hierarchy formation
in social animals. We also discuss how increased interactions between animals
can speed up the asymmetry assessment process and how this can explain why
aggression can sometimes decrease as the number of interactions between animals
increases over time.


  Seasonal influenza is a significant public health concern in the United
States and globally. While influenza vaccines are the single most effective
intervention to reduce influenza morbidity and mortality, there is considerable
debate surrounding the merits and consequences of repeated seasonal
vaccination. Here, we describe a two-season influenza epidemic contact network
model and use it to demonstrate that increasing the level of continuity in
vaccination across seasons reduces the burden on public health. We show that
revaccination reduces the influenza attack rate not only because it reduces the
overall number of susceptible individuals, but also because it better protects
highly-connected individuals, who would otherwise make a disproportionately
large contribution to influenza transmission. Our work thus contributes a
population-level perspective to debates about the merits of repeated influenza
vaccination and advocates for public health policy to incorporate individual
vaccine histories.


  It is likely that the strength of selection acting upon a mutation varies
through time due to changes in the environment. However, most population
genetic theory assumes that the strength of selection remains constant. Here we
investigate the consequences of fluctuating selection pressures on the
quantification of adaptive evolution using McDonald-Kreitman (MK) style
approaches. In agreement with previous work, we show that fluctuating selection
can generate evidence of adaptive evolution even when the expected strength of
selection on a mutation is zero. However, we also find that the mutations,
which contribute to both polymorphism and divergence tend, on average, to be
positively selected during their lifetime, under fluctuating selection models.
This is because mutations that fluctuate, by chance, to positive selected
values, tend to reach higher frequencies in the population than those that
fluctuate towards negative values. Hence the evidence of positive adaptive
evolution detected under a fluctuating selection model by MK type approaches is
genuine since fixed mutations tend to be advantageous on average during their
lifetime. Never-the-less we show that methods tend to underestimate the rate of
adaptive evolution when selection fluctuates.


  The membrane components of cellular organelles have been shown to segregate
into domains as the result of biochemical maturation. We propose that the
dynamical competition between maturation and lateral segregation of membrane
components regulates domain formation. We study a two- component fluid membrane
in which enzymatic reaction irreversibly converts one component into another,
and phase separation triggers the formation of transient membrane domains. The
maximum domains size is shown to depend on the maturation rate as a power-law
similar to the one observed for domain growth with time in the absence of
maturation, despite this time dependence not being verified in the case of
irreversible maturation. This control of domain size by enzymatic activity
could play a critical role in intra-organelle dynamics.


  A reaction--diffusion replicator equation is studied. A novel method to apply
the principle of global regulation is used to write down the model with
explicit spatial structure. Properties of stationary solutions together with
their stability are analyzed analytically, and relationships between stability
of the rest points of the non-distributed replicator equation and distributed
system are shown. A numerical example is given to show that the spatial
variable in this particular model promotes the system's permanence.


  A theoretical framework for the description of Susceptible-Infected-Removed
(SIR) spreading processes with synergistic transmission of infection on a
lattice is developed. The model incorporates explicitly the effects of a
time-dependent environment on the transmission of infection between hosts.
Exact solution of the model shows that time-dependence of the neighbourhood of
recipient hosts is a key factor for synergistic spreading processes. It is
demonstrated that the higher the connectivity of a lattice, the more prominent
is the effect of synergy on spread.


  The transmission dynamics of some infectious diseases is related to the
contact structure between individuals in a network. We used five algorithms to
generate contact networks with different topological structure but with the
same scale-free degree distribution. We simulated the spread of acute and
chronic infectious diseases on these networks, using SI (Susceptible -
Infected) and SIS (Susceptible - Infected - Susceptible) epidemic models. In
the simulations, our objective was to observe the effects of the topological
structure of the networks on the dynamics and prevalence of the simulated
diseases. We found that the dynamics of spread of an infectious disease on
different networks with the same degree distribution may be considerably
different.


  Protein function depends on both protein structure and amino acid (aa)
sequence. Here we show that modular features of both structure and function can
be quantified from the aa sequence alone for the amyloid 770 aa precursor
protein A4. Both the new second order hydropathicity scale, based on
evolutionary optimization (self-organized criticality), and the standard first
order scale, based on complete protein (water-air) unfolding, yield
thermodynamically significant features of amyloid function related to spinodal
decomposition into modules. Spinodal crossings are marked by breaks in slope of
algebraic functions of aa sequences. The new scale is associated with a lower
effective temperature than the standard scale.


  The disappearance of the guanosine triphosphate (GTP)-tubulin cap is widely
believed to be the forerunner event for the growth-shrinkage transition
(`catastrophe') in microtubule filaments in eukaryotic cells. We study a
discrete version of a stochastic model of the GTP cap dynamics, originally
proposed by Flyvbjerg, Holy and Leibler (Flyvbjerg, Holy and Leibler, Phys.
Rev. Lett. 73, 2372, 1994). Our model includes both spontaneous and vectorial
hydrolysis, as well as dissociation of a non-hydrolyzed dimer from the filament
after incorporation. In the first part of the paper, we apply this model to a
single protofilament of a microtubule. A catastrophe transition is defined for
each protofilament, similar to the earlier one-dimensional models, the
frequency of occurrence of which is then calculated under various conditions,
but without explicit assumption of steady state conditions. Using a
perturbative approach, we show that the leading asymptotic behavior of the
protofilament catastrophe in the limit of large growth velocities is remarkably
similar across different models. In the second part of the paper, we extend our
analysis to the entire filament by making a conjecture that a minimum number of
such transitions are required to occur for the onset of microtubule
catastrophe. The frequency of microtubule catastrophe is then determined using
numerical simulations, and compared with analytical/semi-analytical estimates
made under steady state/quasi-steady state assumptions respectively for the
protofilament dynamics. A few relevant experimental results are analyzed in
detail, and compared with predictions from the model. Our results indicate that
loss of GTP cap in 2-3 protofilaments is necessary to trigger catastrophe in a
microtubule.


  Iron is essential for all known life due to its redox properties, however
these same properties can also lead to its toxicity in overload through the
production of reactive oxygen species. Robust systemic and cellular control are
required to maintain safe levels of iron and the liver seems to be where this
regulation is mainly located. Iron misregulation is implicated in many diseases
and as our understanding of iron metabolism improves the list of iron-related
disorders grows. Recent developments have resulted in greater knowledge of the
fate of iron in the body and have led to a detailed map of its metabolism,
however a quantitative understanding at the systems level of how its components
interact to produce tight regulation remains elusive.
  A mechanistic computational model of human liver iron metabolism, which
includes the core regulatory components, was constructed based on known
mechanisms of regulation and on their kinetic properties, obtained from several
publications. The model was then quantitatively validated by comparing its
results with previously published physiological data, and it is able to
reproduce multiple experimental findings. A time course simulation following an
oral dose of iron was compared to a clinical time course study and the
simulation was found to recreate the dynamics and time scale of the systems
response to iron challenge. A disease simulation of haemochromatosis was
created by altering a single reaction parameter that mimics a human
haemochromatosis gene (HFE) mutation. The simulation provides a quantitative
understanding of the liver iron overload that arises in this disease.
  This model supports and supplements understanding of the role of the liver as
an iron sensor and provides a framework for further modelling, including
simulations to identify valuable drug targets and design of experiments to
improve further our knowledge of this system.


  Molecular dynamics simulations have a prominent role in biophysics and drug
discovery due to the atomistic information they provide on the structure,
energetics and dynamics of biomolecules. Specialized software packages are
required to analyze simulated trajectories, either interactively or via
scripts, to derive quantities of interest and provide insight for further
experiments. This paper presents the Density Profile Tool, a package that
enhances the Visual Molecular Dynamics environment with the ability to
interactively compute and visualize 1-D projections of various density
functions of molecular models. We describe how the plugin is used to perform
computations both via a graphical interface and programmatically. Results are
presented for realistic examples, all-atom bilayer models, showing how mass and
electron densities readily provide measurements such as membrane thickness,
location of structural elements, and how they compare to X-ray diffraction
experiments.


  Motivation: The mapping of RNA-seq reads to their transcripts of origin is a
fundamental task in transcript expression estimation and differential
expression scoring. Where ambiguities in mapping exist due to transcripts
sharing sequence, e.g. alternative isoforms or alleles, the problem becomes an
instance of non-trivial probabilistic inference. Bayesian inference in such a
problem is intractable and approximate methods must be used such as Markov
chain Monte Carlo (MCMC) and Variational Bayes. Standard implementations of
these methods can be prohibitively slow for large datasets and complex gene
models.
  Results: We propose an approximate inference scheme based on Variational
Bayes applied to an existing model of transcript expression inference from
RNA-seq data. We apply recent advances in Variational Bayes algorithmics to
improve the convergence of the algorithm beyond the standard variational
expectation-maximisation approach. We apply our algorithm to simulated and
biological datasets, demonstrating that the increase in speed requires only a
small trade-off in accuracy of expression level estimation.
  Availability: The methods were implemented in R and C++, and are available as
part of the BitSeq project at https://code.google.com/p/bitseq/. The methods
will be made available through the BitSeq Bioconductor package at the next
stable release.


  Tumors are defined by their intense proliferation, but sometimes cancer cells
turn senescent and stop replicating. In the stochastic cancer model in which
all cells are tumorigenic, senescence is seen as the result of random
mutations, suggesting that it could represent a barrier to tumor growth. In the
hierarchical cancer model a subset of the cells, the cancer stem cells, divide
indefinitely while other cells eventually turn senescent. Here we formulate
cancer growth in mathematical terms and obtain predictions for the evolution of
senescence. We perform experiments in human melanoma cells which are compatible
with the hierarchical model and show that senescence is a reversible process
controlled by survivin. We conclude that enhancing senescence is unlikely to
provide a useful therapeutic strategy to fight cancer, unless the cancer stem
cells are specifically targeted


  The identification of cancer stem cells in vivo and in vitro relies on
specific surface markers that should allow to sort cancer cells in
phenotypically distinct subpopulations. Experiments report that sorted cancer
cell populations after some time tend to express again all the original
markers, leading to the hypothesis of phenotypic switching, according to which
cancer cells can transform stochastically into cancer stem cells. Here we
explore an alternative explanation based on the hypothesis that markers are not
perfect and are thus unable to identify all cancer stem cells. Our analysis is
based on a mathematical model for cancer cell proliferation that takes into
account phenotypic switching, imperfect markers and error in the sorting
process. Our conclusion is that the observation of reversible expression of
surface markers after sorting does not provide sufficient evidence in support
of phenotypic switching.


  Background: AQP1 belongs to aquaporins family, water-specific,
membrane-channel proteins expressed in diverse tissues. Recent papers showed
that during angiogenesis, AQP1 is expressed preferentially by microvessels,
favoring angiogenesis via the increase of permeability In particular, in AQP1
null mice, endothelial cell migration is impaired without altering their
proliferation or adhesion. Therefore, AQP1 has been proposed as a novel
promoter of tumor angiogenesis. Methods/Findings: Using targeted silencing of
AQP1 gene expression, an impairment in the organization of F-actin and a
reduced migration capacity was demonstrated in human endothelial and melanoma
cell lines. Interestingly, we showed, for the first time, that AQP1
co-immunoprecipitated with Lin-7. Lin7-GFP experiments confirmed
co-immunoprecipitation. In addition, the knock down of AQP1 decreased the level
of expression of Lin-7 and b-catenin and the inhibition of proteasome
contrasted partially such a decrease. Conclusions/Significance: All together,
our findings show that AQP1 plays a role inside the cells through
Lin-7/b-catenin interaction. Such a role of AQP1 is the same in human melanoma
and endothelial cells, suggesting that AQP1 plays a global physiological role.
A model is presented.


  Background: A fundamental problem in cancer research is identifying the cell
type that is capable of sustaining neoplastic growth and its origin from normal
tissue cells. Recent investigations of a variety of tumor types have shown that
phenotypically identifiable and isolable subfractions of cells possess the
tumor-forming ability. In the present paper, using two lineage-related human
melanoma cell lines, primary melanoma line IGR39 and its metastatic derivative
line IGR37, two main observations are reported. The first one is the first
phenotypic evidence to support the origin of melanoma cancer stem cells (CSCs)
from mutated tissue-specific stem cells; and the second one is the
identification of a more aggressive subpopulation of CSCs in melanoma that are
CXCR6+. Conclusions/Significance: The association of a more aggressive tumor
phenotype with asymmetric self-renewal phenotype reveals a previously
unrecognized aspect of tumor cell physiology. Namely, the retention of some
tissue-specific stem cell attributes, like the ability to asymmetrically
self-renew, impacts the natural history of human tumor development. Knowledge
of this new aspect of tumor development and progression may provide new targets
for cancer prevention and treatment.


  Over the recent years, the field of whole metagenome shotgun sequencing has
witnessed significant growth due to the high-throughput sequencing technologies
that allow sequencing genomic samples cheaper, faster, and with better coverage
than before. This technical advancement has initiated the trend of sequencing
multiple samples in different conditions or environments to explore the
similarities and dissimilarities of the microbial communities. Examples include
the human microbiome project and various studies of the human intestinal tract.
With the availability of ever larger databases of such measurements, finding
samples similar to a given query sample is becoming a central operation. In
this paper, we develop a content-based exploration and retrieval method for
whole metagenome sequencing samples. We apply a distributed string mining
framework to efficiently extract all informative sequence $k$-mers from a pool
of metagenomic samples and use them to measure the dissimilarity between two
samples. We evaluate the performance of the proposed approach on two human gut
metagenome data sets as well as human microbiome project metagenomic samples.
We observe significant enrichment for diseased gut samples in results of
queries with another diseased sample and very high accuracy in discriminating
between different body sites even though the method is unsupervised. A software
implementation of the DSM framework is available at
https://github.com/HIITMetagenomics/dsm-framework


  We study the patterns generated in finite-time sweeps across
symmetry-breaking bifurcations in individual-based models. Similar to the
well-known Kibble-Zurek scenario of defect formation, large-scale patterns are
generated when model parameters are varied slowly, whereas fast sweeps produce
a large number of small domains. The symmetry breaking is triggered by
intrinsic noise, originating from the discrete dynamics at the micro-level.
Based on a linear-noise approximation, we calculate the characteristic length
scale of these patterns. We demonstrate the applicability of this approach in a
simple model of opinion dynamics, a model in evolutionary game theory with a
time-dependent fitness structure, and a model of cell differentiation. Our
theoretical estimates are confirmed in simulations. In further numerical work,
we observe a similar phenomenon when the symmetry-breaking bifurcation is
triggered by population growth.


  The attempt to determine the population growth rate from field data reveals
several ambiguities in its definition(s), which seem to throw into question the
very concept itself. However, an alternative point of view is proposed that not
only preserves the identity of the concept, but also helps discriminate between
competing models for capturing the data.


  An autocatalytic pattern matching polymer system is studied as an abstract
model for chemical ecosystem evolution. Highly ordered populations with
particular sequence patterns appear spontaneously out of a vast number of
possible states. The interplay between the selected microscopic sequence
patterns and the macroscopic cooperative structures is examined. Stability,
fluctuations, and evolutionary selection mechanisms are investigated for the
involved self-organizing processes.


  Using the growth of population in Australia in the past 10,000 years it is
illustrated here how an illusion created by hyperbolic distributions may lead
easily to incorrect conclusions. Contrary to the published claim, there was no
change in the mechanism of growth of the ancient human population in Australia
around 5000 years before present (BP). Data for the number of rock-shelter
sites, interpreted as reflecting the growth of human population, are analysed.
They are shown to be monotonically increasing with time without any sign of
intensification in the past 10,000 years. The growth of human population in
Australia is in excellent agreement with the similar pattern describing the
growth of the world population.


  The amount of non-unique sequence (non-singletons) in a genome directly
affects the difficulty of read alignment to a reference assembly for high
throughput-sequencing data. Although a greater length increases the chance for
reads being uniquely mapped to the reference genome, a quantitative analysis of
the influence of read lengths on mappability has been lacking. To address this
question, we evaluate the k-mer distribution of the human reference genome. The
k-mer frequency is determined for k ranging from 20 to 1000 basepairs. We use
the proportion of non-singleton k-mers to evaluate the mappability of reads for
a corresponding read length. We observe that the proportion of non-singletons
decreases slowly with increasing k, and can be fitted by piecewise power-law
functions with different exponents at different k ranges. A faster decay at
smaller values for k indicates more limited gains for read lengths > 200
basepairs. The frequency distributions of k-mers exhibit long tails in a
power-law-like trend, and rank frequency plots exhibit a concave Zipf's curve.
The location of the most frequent 1000-mers comprises 172 kilobase-ranged
regions, including four large stretches on chromosomes 1 and X, containing
genes with biomedical implications. Even the read length 1000 would be
insufficient to reliably sequence these specific regions.


  In this comment, we analyze the mathematical model by Holst et al that
estimates the increased mortality of Inachis io larvae caused by Bt-maize
pollen in European farmaland. In particular, we focus our attention on the part
of Holst et al's work that models the Bt-maize pollen loss on plant leaves by a
differential equation. We provide the exact solution of the differential
equation, showing that the behaviour of the solution, when the time goes to
infinity, is different from what claimed in Holst et al's work.


  Our current understanding of the taxonomic and phylogenetic diversity of
cellular organisms, especially the bacteria and archaea, is mostly based upon
studies of sequences of the small- subunit rRNAs (ssu-rRNAs). To address the
limitation of ssu-rRNA as a phylogenetic marker, such as copy number variation
among organisms and complications introduced by horizontal gene transfer,
convergent evolution, or evolution rate variations, we have identified protein-
coding gene families as alternative Phylogenetic and Phylogenetic Ecology
markers (PhyEco). Current nucleotide sequence similarity based Operational
Taxonomic Unit (OTU) classification methods are not readily applicable to amino
acid sequences of PhyEco markers. We report here the development of TreeOTU, a
phylogenetic tree structure based OTU classification method that takes into
account of differences in rates of evolution between taxa and between genes.
OTU sets built by TreeOTU are more faithful to phylogenetic tree structures
than sequence clustering (non phylogenetic) methods for ssu-rRNAs. OTUs built
from phylogenetic trees of protein coding PhyEco markers are comparable to our
current taxonomic classification at different levels. With the included OTU
comparing tools, the TreeOTU is robust in phylogenetic referencing with
different phylogenetic markers and trees.


  The illness-death model of a chronic disease consists of the states 'Normal',
'Disease' and 'Death'. In general, the transition rates between the states
depend on three time scales: calendar time, age and duration of the chronic
disease. Previous works have shown that the age-specific prevalence of the
chronic disease can be described by differential equations if the duration is
negligible. This article derives a partial differential equation (PDE) in the
presence of duration dependency. As an important application, the PDE allows
the calculation of the age-specific incidence from cross-sectional surveys.


  We have investigated the role that the mutation rate and the structure of
genetic variation at a locus play in determining whether a gene is involved in
disease. We predict that the mutation rate and its genetic diversity should be
higher in genes associated with disease, unless all genes that could cause
disease have already been identified. Consistent with our predictions we find
that genes associated with Mendelian and complex disease are substantially
longer than non-disease genes. However, we find that both Mendelian and complex
disease genes are found in regions of the genome with relatively low mutation
rates, as inferred from intron divergence between humans and chimpanzees.
Complex disease gene are predicted to have higher rates of non-synonymous
mutation than non-disease genes, but the opposite pattern is found in Mendelian
disease genes. Finally, we find that disease genes are in regions of
significantly elevated genetic diversity, even when variation in the rate of
mutation is controlled for. The effect is small nevertheless. Our results
suggest that variation in the genic mutation rate and the genetic architecture
of the locus play a minor role in determining whether a gene is associated with
disease.


  I discuss the so-called SuperSpreader epidemic, for which SARS is the
canonical examples (and, perhaps, MERS will be another). I use simulation by an
agent-based model as well as the mathematics of multi-type branching-processes
to illustrate how the SS epidemic differs from the more familiar uniform
epidemic (e.g., caused by influenza). The conclusions may surprise the reader:
(a) The SS epidemic must be described by at least two numbers, such as the mean
reproductive number (of "secondary" cases caused by a "primary case"), R0, and
the variance of same, call it V0; (b) Even if R0 > 1, if V0 >> R0 the
probability that an infection-chain caused by one primary case goes extinct
without intervention may be close to one (e.g., 0.97); (c) The SS epidemic may
have a long "kindling period" in which sporadic cases appear (transmitted from
some unknown host) and generate a cluster of cases, but the chains peter out,
perhaps generating a false sense of security that a pandemic will not occur;
(d) Interventions such as isolation (or contact-tracing and secondary case
isolation) may prove efficacious even without driving R0 below one; (e) The
efficacy of such interventions diminishes, but slowly, with increasing V0 at
fixed R0. From these considerations, I argue that the SS epidemic has dynamics
sufficiently distinct from the uniform case that efficacious public-health
interventions can be designed even in the absence of a vaccine or other form of
treatment.


  In this article, we introduce a novel backward method to model stochastic
gene expression and protein level dynamics. The protein amount is regarded as a
diffusion process and is described by a backward stochastic differential
equation (BSDE). Unlike many other SDE techniques proposed in the literature,
the BSDE method is backward in time; that is, instead of initial conditions it
requires the specification of endpoint ("final") conditions, in addition to the
model parametrization. To validate our approach we employ Gillespie's
stochastic simulation algorithm (SSA) to generate (forward) benchmark data,
according to predefined gene network models. Numerical simulations show that
the BSDE method is able to correctly infer the protein level distributions that
preceded a known final condition, obtained originally from the forward SSA.
This makes the BSDE method a powerful systems biology tool for time reversed
simulations, allowing, for example, the assessment of the biological conditions
(e.g. protein concentrations) that preceded an experimentally measured event of
interest (e.g. mitosis, apoptosis, etc.).


  It has long been accepted that the multiple-ion single-file transport model
is appropriate for many kinds of ion channels. However, most of the purely
theoretical works in this field did not capture all of the important features
of the realistic systems. Nowadays, large-scale atomic-level simulations are
more feasible. Discrepancy between theories, simulations and experiments are
getting obvious, enabling people to carefully examine the missing parts of the
theoretical models and methods. In this work, it is attempted to find out the
essential features that such kind of models should possess, in order that the
physical properties of an ion channel be adequately reflected.


  Accumulating evidence suggests that many tumors have a hierarchical
organization, with the bulk of the tumor composed of relatively differentiated
short-lived progenitor cells that are maintained by a small population of
undifferentiated long-lived cancer stem cells. It is unclear, however, whether
cancer stem cells originate from normal stem cells or from dedifferentiated
progenitor cells. To address this, we mathematically modeled the effect of
dedifferentiation on carcinogenesis. We considered a hybrid
stochastic-deterministic model of mutation accumulation in both stem cells and
progenitors, including dedifferentiation of progenitor cells to a stem
cell-like state. We performed exact computer simulations of the emergence of
tumor subpopulations with two mutations, and we derived semi-analytical
estimates for the waiting time distribution to fixation. Our results suggest
that dedifferentiation may play an important role in carcinogenesis, depending
on how stem cell homeostasis is maintained. If the stem cell population size is
held strictly constant (due to all divisions being asymmetric), we found that
dedifferentiation acts like a positive selective force in the stem cell
population and thus speeds carcinogenesis. If the stem cell population size is
allowed to vary stochastically with density-dependent reproduction rates
(allowing both symmetric and asymmetric divisions), we found that
dedifferentiation beyond a critical threshold leads to exponential growth of
the stem cell population. Thus, dedifferentiation may play a crucial role, the
common modeling assumption of constant stem cell population size may not be
adequate, and further progress in understanding carcinogenesis demands a more
detailed mechanistic understanding of stem cell homeostasis.


  Psammophilous plants are special plants that flourish in sand moving
environments. There are two main mechanisms by which the wind affects these
plants: (i) sand drift exposes roots and covers branches--the exposed roots
turn into new plants and the covered branches turn into new roots; both
mechanisms result in an enhanced growth rate of the psammophilous plant cover
of the dunes; (ii) strong winds, often associated with sand movement, tear
branches and seed them in nearby locations, resulting in new plants and an
enhanced growth rate of the psammophilous plant cover of the dunes. Despite
their important role in dune dynamics, to our knowledge, psammophilous plants
have never been incorporated into mathematical models of sand dunes. Here, we
attempt to model the effects of these plants on sand dune dynamics. We
construct a set of three ordinary differential equations for the fractions of
surface cover of regular vegetation, biogenic soil crust and psammophilous
plants. The latter reach their optimal growth under (i) specific sand drift or
(ii) specific wind power. We show that psammophilous plants enrich the sand
dune dynamics. Depending on the climatological conditions, it is possible to
obtain one, two, or three steady dune states. The activity of the dunes can be
associated with the surface cover--bare dunes are active, and dunes with
significant cover of vegetation, biogenic soil crust, or psammophilous plants
are fixed. Our model shows that under suitable precipitation rates and wind
power, the dynamics of the different cover types is in accordance with the
common view that dunes are initially stabilized by psammophilous plants that
reduce sand activity, thus enhancing the growth of regular vegetation that
eventually dominates the cover of the dunes and determines their activity.


  We consider a general mathematical model of a within-host viral infection
with $n$ virus strains and explicit age-since-infection structure for infected
cells. In the model, multiple virus strains compete for a population of target
cells. Cells infected with virus strain $i\in\left\{1,...,n\right\}$ die at
per-capita rate $\delta_i(a)$ and produce virions at per-capita rate $p_i(a)$,
where $\delta_i(a)$ and $p_i(a)$ are functions of the age-since-infection of
the cell. Viral strain $i$ has a basic reproduction number, $\mathcal{R}_i$,
and a corresponding positive single strain equilibrium, $E_i$, when
$\mathcal{R}_i>1$. If $\mathcal{R}_i<1$, then the total concentration of virus
strain $i$ will converge to 0 asymptotically. The main result is that when
$\max_i \mathcal{R}_i>1$ and all of the reproduction numbers are distinct, i.e.
$\mathcal{R}_i\neq \mathcal{R}_j \ \forall i\neq j$, the viral strain with the
maximal basic reproduction number competitively excludes the other strains. As
an application of the model, HIV evolution is considered and simulations are
provided.


  The emerging threat of a human pandemic caused by the H5N1 avian influenza
virus strain magnifies the need for controlling the incidence of H5N1 infection
in domestic bird populations. Culling is one of the most widely used control
measures and has proved effective for isolated outbreaks. However, the
socio-economic impacts of mass culling, in the face of a disease which has
become endemic in many regions of the world, can affect the implementation and
success of culling as a control measure. We use mathematical modeling to
understand the dynamics of avian influenza under different culling approaches.
We incorporate culling into an SI model by considering the per capita culling
rates to be general functions of the number of infected birds. Complex dynamics
of the system, such as backward bifurcation and forward hysteresis, along with
bi-stability, are detected and analyzed for two distinct culling scenarios. In
these cases, employing other control measures temporarily can drastically
change the dynamics of the solutions to a more favorable outcome for disease
control.


  Mass-vaccination campaigns are an important strategy in the global fight
against poliomyelitis and measles. The large-scale logistics required for these
mass immunisation campaigns magnifies the need for research into the
effectiveness and optimal deployment of pulse vaccination. In order to better
understand this control strategy, we propose a mathematical model accounting
for the disease dynamics in connected regions, incorporating seasonality,
environmental reservoirs and independent periodic pulse vaccination schedules
in each region. The effective reproduction number, $R_e$, is defined and proved
to be a global threshold for persistence of the disease. Analytical and
numerical calculations show the importance of synchronising the pulse
vaccinations in connected regions and the timing of the pulses with respect to
the pathogen circulation seasonality. Our results indicate that it may be
crucial for mass-vaccination programs, such as national immunisation days, to
be synchronised across different regions. In addition, simulations show that a
migration imbalance can increase $R_e$ and alter how pulse vaccination should
be optimally distributed among the patches, similar to results found with
constant-rate vaccination. Furthermore, contrary to the case of constant-rate
vaccination, the fraction of environmental transmission affects the value of
$R_e$ when pulse vaccination is present.


  Translating metabolic networks into dynamic models is difficult if kinetic
constants are unknown. Structural Kinetic Modelling (SKM) replaces reaction
elasticities by independent random numbers. Here I propose a variant that
accounts for reversible reactions and thermodynamics: in Structural
Thermokinetic Modelling (STM), correlated elasticities are computed from enzyme
saturation values and thermodynamic forces, which are physically independent.
STM relies on a dependency schema in which basic variables can be sampled,
fitted to data, or optimised, while all other variables are computed from them.
Probability distributions in the dependency schema define a model ensemble,
which leads to probabilistic predictions even if data are scarce. STM
highlights the importance of variabilities, dependencies and covariances of
biological variables. By choosing or sampling the basic variables, we can
convert metabolic networks into kinetic models with consistent reversible rate
laws. Metabolic control coefficients obtained from these models can tell us
about metabolic dynamics, including responses and optimal adaptations to
perturbations as well as enzyme synergies, metabolite correlations, and
metabolic fluctuations arising from chemical noise. By comparing model variants
with different network structures, fluxes, thermodynamic forces, regulation, or
types of rate laws, we can quantify the effects of these model features. To
showcase STM, I study metabolic control, metabolic fluctuations, and enzyme
synergies, and how they are shaped by thermodynamic forces. Thermodynamics can
be used to obtain more precise predictions of flux control, enzyme synergies,
correlated flux and metabolite variations, and of the emergence and propagation
of metabolic noise.


  Genetic drift at the frontiers of two-dimensional range expansions of
microorganisms can frustrate local cooperation between different genetic
variants, demixing the population into distinct sectors. In a biological
context, mutualistic or antagonistic interactions will typically be asymmetric
between variants. By taking into account both the asymmetry and the interaction
strength, we show that the much weaker demixing in three dimensions allows for
a mutualistic phase over a much wider range of asymmetric cooperative benefits,
with mutualism prevailing for any positive, symmetric benefit. We also
demonstrate that expansions with undulating fronts roughen dramatically at the
boundaries of the mutualistic phase, with severe consequences for the
population genetics along the transition lines.


  The typical process for classifying and submitting a newly sequenced virus to
the NCBI database involves two steps. First, a BLAST search is performed to
determine likely family candidates. That is followed by checking the candidate
families with the Pairwise Sequence Alignment tool for similar species. The
submitter's judgement is then used to determine the most likely species
classification. The aim of this paper is to show that this process can be
automated into a fast, accurate, one-step process using the proposed
alignment-free method and properly implemented machine learning techniques.
  We present a new family of alignment-free vectorizations of the genome, the
generalized vector, that maintains the speed of existing alignment-free methods
while outperforming all available methods. This new alignment-free
vectorization uses the frequency of genomic words (k-mers), as is done in the
composition vector, and incorporates descriptive statistics of those k-mers'
positional information, as inspired by the natural vector.
  We analyze 5 different characterizations of genome similarity using
$k$-nearest neighbor classification, and evaluate these on two collections of
viruses totaling over 10,000 viruses. We show that our proposed method performs
better than, or as well as, other methods at every level of the phylogenetic
hierarchy.
  The data and R code is available upon request.


  The breeder's equation is a cornerstone of quantitative genetics and is
widely used in evolutionary modeling. The equation which reads R=h^{2}S relates
response to selection R (the mean phenotype of the progeny) to the selection
differential S (mean phenotype of selected parents) through a simple
proportionality relation. The validity of this relation however relies strongly
on the normal (Gaussian) distribution of parent's genotype which is an
unobservable quantity and cannot be ascertained. In contrast, we show here that
if the fitness (or selection) function is Gaussian, an alternative, exact
linear equation in the form of R'=j^{2}S' can be derived, regardless of the
parental genotype distribution. Here R' and S' stand for the mean phenotypic
lag behind the mean of the fitness function in the offspring and selected
populations. To demonstrate this relation, we derive the exact functional
relation between the mean phenotype in the selected and the offspring
population and deduce all cases that lead to a linear relation between these
quantities. These computations, which are confirmed by individual based
numerical simulations, generalize naturally to the multivariate Lande's
equation \Delta\mathbf{\bar{z}}=GP^{-1}\mathbf{S} .


  We demonstrate that the coupling of excitonic and vibrational motion in
biological complexes can provide mechanisms to explain the long-lived
oscillations that have been obtained in non linear spectroscopic signals of
different photosynthetic pigment protein complexes and we discuss the
contributions of excitonic versus purely vibrational components to these
oscillatory features. Considering a dimer model coupled to a structured
spectral density we exemplify the fundamental aspects of the electron-phonon
dynamics, and by analyzing separately the different contributions to the non
linear signal, we show that for realistic parameter regimes purely electronic
coherence is of the same order as purely vibrational coherence in the
electronic ground state. Moreover, we demonstrate how the latter relies upon
the excitonic interaction to manifest. These results link recently proposed
microscopic, non-equilibrium mechanisms to support long lived coherence at
ambient temperatures with actual experimental observations of oscillatory
behaviour using 2D photon echo techniques to corroborate the fundamental
importance of the interplay of electronic and vibrational degrees of freedom in
the dynamics of light harvesting aggregates.


  The phenotypic consequences of individual mutations are modulated by the wild
type genetic background in which they occur.Although such background dependence
is widely observed, we do not know whether general patterns across species and
traits exist, nor about the mechanisms underlying it. We also lack knowledge on
how mutations interact with genetic background to influence gene expression,
and how this in turn mediates mutant phenotypes. Furthermore, how genetic
background influences patterns of epistasis remains unclear. To investigate the
genetic basis and genomic consequences of genetic background dependence of the
scallopedE3 allele on the Drosophila melanogaster wing, we generated multiple
novel genome level datasets from a mapping by introgression experiment and a
tagged RNA gene expression dataset. In addition we used whole genome
re-sequencing of the parental lines two commonly used laboratory strains to
predict polymorphic transcription factor binding sites for SD. We integrated
these data with previously published genomic datasets from expression
microarrays and a modifier mutation screen. By searching for genes showing a
congruent signal across multiple datasets, we were able to identify a robust
set of candidate loci contributing to the background dependent effects of
mutations in sd. We also show that the majority of background-dependent
modifiers previously reported are caused by higher-order epistasis, not
quantitative non-complementation. These findings provide a useful foundation
for more detailed investigations of genetic background dependence in this
system, and this approach is likely to prove useful in exploring the genetic
basis of other traits as well.


  The DevRS two component system of Mycobacterium tuberculosis is responsible
for its dormancy in host and becomes operative under hypoxic condition. It is
experimentally known that phosphorylated DevR controls the expression of
several downstream genes in a complex manner. In the present work we propose a
theoretical model to show role of binding sites in DevR mediated gene
expression. Individual and collective role of binding sites in regulating DevR
mediated gene expression has been shown via modeling. Objective of the present
work is two fold. First, to describe qualitatively the temporal dynamics of
wild type genes and their known mutants. Based on these results we propose that
DevR controlled gene expression follows a specific pattern which is efficient
in describing other DevR mediated gene expression. Second, to analyze behavior
of the system from information theoretical point of view. Using the tools of
information theory we have calculated molecular efficiency of the system and
have shown that it is close to the maximum limit of isothermal efficiency.


  Genotype networks are a method used in systems biology to study the
"innovability" of a set of genotypes having the same phenotype. In the past
they have been applied to determine the genetic heterogeneity, and stability to
mutations, of systems such as metabolic networks and RNA folds. Recently, they
have been the base for re-conciliating the two neutralist and selectionist
schools on evolution.
  Here, we adapted the concept of genotype networks to the study of population
genetics data, applying them to the 1000 Genomes dataset. We used networks
composed of short haplotypes of Single Nucleotide Variants (SNV), and defined
phenotypes as the presence or absence of a haplotype in a human population. We
used coalescent simulations to determine if the number of samples in the 1000
Genomes dataset is large enough to represent the genetic variation of real
populations. The result is a scan of how properties related to the genetic
heterogeneity and stability to mutations are distributed along the human
genome. We found that genes involved in acquired immunity, such as some HLA and
MHC genes, tend to have the most heterogeneous and connected networks; and we
have also found that there is a small, but significant difference between
networks of coding regions and those of non-coding regions, suggesting that
coding regions are both richer in genotype diversity, and more stable to
mutations. Together, the work presented here may constitute a starting point
for applying genotype networks to study genome variation, as larger datasets of
next-generation data will become availa


  Non-Poissonian bursty processes are ubiquitous in natural and social
phenomena, yet little is known about their effects on the large-scale spreading
dynamics. In order to characterize these effects we devise an analytically
solvable model of Susceptible-Infected (SI) spreading dynamics in infinite
systems for arbitrary inter-event time distributions and for the whole time
range. Our model is stationary from the beginning, and the role of lower bound
of inter-event times is explicitly considered. The exact solution shows that
for early and intermediate times the burstiness accelerates the spreading as
compared to a Poisson-like process with the same mean and same lower bound of
inter-event times. Such behavior is opposite for late time dynamics in finite
systems, where the power-law distribution of inter-event times results in a
slower and algebraic convergence to fully infected state in contrast to the
exponential decay of the Poisson-like process. We also provide an intuitive
argument for the exponent characterizing algebraic convergence.


  Here we propose a new approach to modeling gene expression based on the
theory of random dynamical systems (RDS) that provides a general coupling
prescription between the nodes of any given regulatory network given the
dynamics of each node is modeled by a RDS. The main virtues of this approach
are the following: (i) it provides a natural way to obtain arbitrarily large
networks by coupling together simple basic pieces, thus revealing the
modularity of regulatory networks; (ii) the assumptions about the stochastic
processes used in the modeling are fairly general, in the sense that the only
requirement is stationarity; (iii) there is a well developed mathematical
theory, which is a blend of smooth dynamical systems theory, ergodic theory and
stochastic analysis that allows one to extract relevant dynamical and
statistical information without solving the system; (iv) one may obtain the
classical rate equations form the corresponding stochastic version by averaging
the dynamic random variables (small noise limit). It is important to emphasize
that unlike the deterministic case, where coupling two equations is a trivial
matter, coupling two RDS is non-trivial, specially in our case, where the
coupling is performed between a state variable of one gene and the switching
stochastic process of another gene and, hence, it is not a priori true that the
resulting coupled system will satisfy the definition of a random dynamical
system. We shall provide the necessary arguments that ensure that our coupling
prescription does indeed furnish a coupled regulatory network of random
dynamical systems. Finally, the fact that classical rate equations are the
small noise limit of our stochastic model ensures that any validation or
prediction made on the basis of the classical theory is also a validation or
prediction of our model.


  Close analysis of the published interpretation of the number of rock-shelter
sites in Australia provides further evidence that there was no intensification
in the growth of human population between 1000 and 10,000 years BP. An
alternative way of determining the time-dependent distribution of the size of
human population between 1000 and 10,000 years BP is discussed.


  We study the evolution of allele frequencies in a large population where
random mating is violated in a particular way that is related to recent works
on speciation. Specifically, we consider non-random encounters in haploid
organisms described by biallelic genes at two loci and assume that individuals
whose alleles differ at both loci are incompatible. We show that evolution
under these conditions leads to the disappearance of one of the alleles and
substantially reduces the diversity of the population. The allele that
disappears, and the other allele frequencies at equilibrium, depend only on
their initial values, and so does the time to equilibration. However, certain
combinations of allele frequencies remain constant during the process,
revealing the emergence of strong correlation between the two loci promoted by
the epistatic mechanism of incompatibility. We determine the geometrical
structure of the haplotype frequency space and solve the dynamical equations,
obtaining a simple rule to determine equilibrium solution from the initial
conditions. We show that our results are equivalent to selection against double
heterozigotes for a population of diploid individuals and discuss the relevance
of our findings to speciation.


  Vascular endothelial growth factor (VEGF) signaling is involved in the
process of blood vessel development and maintenance. Signaling is initiated by
binding of the bivalent VEGF ligand to the membrane-bound receptors (VEGFR),
which in turn stimulates receptor dimerization. Herein, we discuss experimental
evidence that VEGF receptors localize in caveloae and other regions of the
plasma membrane, and for other receptors, it has been shown that receptor
clustering has an impact on dimerization and thus also on signaling. Overall,
receptor clustering is part of a complex ecosystem of interactions and how
receptor clustering impacts dimerization is not well understood. To address
these questions, we have formulated the simplest possible model. We have
postulated the existence of a single high affinity region in the cell membrane,
which acts as a transient trap for receptors. We have defined an ODE model by
introducing high- and low-density receptor variables and introduce the
corresponding reactions from a realistic model of VEGF signal initiation.
Finally, we use the model to investigate the relation between the degree of
VEGFR concentration, ligand availability, and signaling. In conclusion, our
simulation results provide a deeper understanding of the role of receptor
clustering in cell signaling.


  Hybrid modeling provides an effective solution to cope with multiple time
scales dynamics in systems biology. Among the applications of this method, one
of the most important is the cell cycle regulation. The machinery of the cell
cycle, leading to cell division and proliferation, combines slow growth,
spatio-temporal re-organisation of the cell, and rapid changes of regulatory
proteins concentrations induced by post-translational modifications. The
advancement through the cell cycle comprises a well defined sequence of stages,
separated by checkpoint transitions. The combination of continuous and discrete
changes justifies hybrid modelling approaches to cell cycle dynamics. We
present a piecewise-smooth version of a mammalian cell cycle model, obtained by
hybridization from a smooth biochemical model. The approximate hybridization
scheme, leading to simplified reaction rates and binary event location
functions, is based on learning from a training set of trajectories of the
smooth model. We discuss several learning strategies for the parameters of the
hybrid model.


  This paper presents a method for modeling biological systems which combines
formal techniques on intervals, numerical simulations and satisfaction of
Signal Temporal Logic (STL) formulas. The main modeling challenge addressed by
this approach is the large uncertainty in the values of the parameters due to
the experimental difficulties of getting accurate biological data. This method
considers intervals for each parameter and a formal description of the expected
behavior of the model. In a first step, it produces reduced intervals of
possible parameter values. Then by performing a systematic search in these
intervals, it defines sets of parameter values used in the next step. This
procedure aims at finding a sub-space where the model robustly behaves as
expected. We apply this method to the modeling of the cellular iron homeostasis
network in erythroid progenitors. The produced model describes explicitly the
regulation mechanism which acts at the translational level.


  Myopia in human and animals is caused by the axial elongation of the eye and
is closely linked to the thinning of the sclera which supports the eye tissue.
This thinning has been correlated with the overproduction of matrix
metalloproteinase (MMP-2), an enzyme that degrades the collagen structure of
the sclera. In this short paper, we propose a descriptive model of a regulatory
network with hysteresis, which seems necessary for creating oscillatory
behavior in the hybrid model between MMP-2, MT1-MMP and TIMP-2. Numerical
results provide insight on the type of equilibria present in the system.


  Predicting the ancestral sequences of a group of homologous sequences related
by a phylogenetic tree has been the subject of many studies, and numerous
methods have been proposed to this purpose. Theoretical results are available
that show that when the mutation rate become too large, reconstructing the
ancestral state at the tree root is no longer feasible. Here, we also study the
reconstruction of the ancestral changes that occurred along the tree edges. We
show that, depending on the tree and branch length distribution, reconstructing
these changes (i.e. reconstructing the ancestral state of all internal nodes in
the tree) may be easier or harder than reconstructing the ancestral root state.
However, results from information theory indicate that for the standard Yule
tree, the task of reconstructing internal node states remains feasible, even
for very high substitution rates. Moreover, computer simulations demonstrate
that for more complex trees and scenarios, this result still holds. For a large
variety of counting, parsimony-based and likelihood-based methods, the
predictive accuracy of a randomly selected internal node in the tree is indeed
much higher than the accuracy of the same method when applied to the tree root.
Moreover, parsimony- and likelihood-based methods appear to be remarkably
robust to sampling bias and model mis-specification.


  Messenger RNA encodes a sequence of amino acids by using codons. For most
amino acids there are multiple synonymous codons that can encode the amino
acid. The translation speed can vary from one codon to another, thus there is
room for changing the ribosome speed while keeping the amino acid sequence and
hence the resulting protein. Recently, it has been noticed that the choice of
the synonymous codon, via the resulting distribution of slow- and
fast-translated codons, affects not only on the average speed of one ribosome
translating the messenger RNA (mRNA) but also might have an effect on nearby
ribosomes by affecting the appearance of "traffic jams" where multiple
ribosomes collide and form queues. To test this "context effect" further, we
here investigate the effect of the sequence of synonymous codons on the
ribosome traffic by using a ribosome traffic model with codon-dependent rates,
estimated from experiments. We compare the ribosome traffic on wild type
sequences and sequences where the synonymous codons were swapped randomly. By
simulating translation of 87 genes, we demonstrate that the wild type
sequences, especially those with a high bias in codon usage, tend to have the
ability to reduce ribosome collisions, hence optimizing the cellular investment
in the translation apparatus. The magnitude of such reduction of the
translation time might have a significant impact on the cellular growth rate
and thereby have importance for the survival of the species.


  In bacteria, regulatory proteins search for a specific DNA binding target via
"facilitated diffusion": a series of rounds of 3D diffusion in the cytoplasm,
and 1D linear diffusion along the DNA contour. Using large scale Brownian
dynamics simulations we find that each of these steps is affected differently
by crowding proteins, which can either be bound to the DNA acting as a road
block to the 1D diffusion, or freely diffusing in the cytoplasm. Macromolecular
crowding can strongly affect mechanistic features such as the balance between
3D and 1D diffusion, but leads to surprising robustness of the total search
time.


  Quantifying convergence and sufficient sampling of macromolecular molecular
dynamics simulations is more often than not a source of controversy (and of
various ad hoc solutions) in the field. Clearly, the only reasonable,
consistent and satisfying way to infer convergence (or otherwise) of a
molecular dynamics trajectory must be based on probability theory. Ideally, the
question we would wish to answer is the following : "What is the probability
that a molecular configuration important for the analysis in hand has not yet
been observed ?". Here we propose a method for answering a variant of this
question by using the Good-Turing formalism for frequency estimation of
unobserved species in a sample. Although several approaches may be followed in
order to deal with the problem of discretizing the configurational space, for
this work we use the classical RMSD matrix as a means to answering the
following question: "What is the probability that a molecular configuration
with an RMSD (from all other already observed configurations) higher than a
given threshold has not actually been observed ?". We apply the proposed method
to several different trajectories and show that the procedure appears to be
both computationally stable and internally consistent. A free, open-source
program implementing these ideas is immediately available for download via
public repositories.


  This paper presents an accurate short-read mapper for next-generation
sequencing data which is widely used in the 1000 Genomes Project, and human
clinical and other species genome studies.


  Studies on the genetics of adaptation typically neglect the possibility that
a deleterious mutation might fix. Nonetheless, here we show that, in many
regimes, the first substitution is most often deleterious, even when fitness is
expected to increase in the long term. In particular, we prove that this
phenomenon occurs under weak mutation for any house-of-cards model with an
equilibrium distribution. We find that the same qualitative results hold under
Fisher's geometric model. We also provide a simple intuition for the surprising
prevalence of unconditionally deleterious substitutions during early
adaptation. Importantly, the phenomenon we describe occurs on fitness
landscapes without any local maxima and is therefore distinct from
"valley-crossing". Our results imply that the common practice of ignoring
deleterious substitutions leads to qualitatively incorrect predictions in many
regimes. Our results also have implications for the substitution process at
equilibrium and for the response to a sudden decrease in population size.


  Primary tumors often emerge within genetically altered fields of premalignant
cells that appear histologically normal but have a high chance of progression
to malignancy. Clinical observations have suggested that these premalignant
fields pose high risks for emergence of secondary recurrent tumors if left
behind after surgical removal of the primary tumor. In this work, we develop a
spatio-temporal stochastic model of epithelial carcinogenesis, combining
cellular reproduction and death dynamics with a general framework for
multi-stage genetic progression to cancer. Using this model, we investigate how
macroscopic features (e.g. size and geometry of premalignant fields) depend on
microscopic cellular properties of the tissue (e.g.\ tissue renewal rate,
mutation rate, selection advantages conferred by genetic events leading to
cancer, etc). We develop methods to characterize how clinically relevant
quantities such as waiting time until emergence of second field tumors and
recurrence risk after tumor resection. We also study the clonal relatedness of
recurrent tumors to primary tumors, and analyze how these phenomena depend upon
specific characteristics of the tissue and cancer type. This study contributes
to a growing literature seeking to obtain a quantitative understanding of the
spatial dynamics in cancer initiation.


  In this paper, I tentatively put forward a new hypothesis that the emergence
of a single chiral form of biomolecules in living organisms is specifically
determined by the electron spin state during their enzyme-catalyzed synthesis
processes. Specifically speaking, the electrons released from the coenzyme
NAD(P)H of amino acid synthase are heterogeneous in spin states; however, when
they pass through the chiral alpha-helix structure of the enzymes to the site
of amino acid synthesis at the other end of the helix, their spin states are
filtered and polarized, producing only spin up electrons; once the
spin-polarized electrons participate in the reductive reaction between
alpha-oxo acid and ammonia, only L-amino acids are formed according to the
Pauli exclusion principle.


  The functional diversity in factors identified in association with
non-alcoholic fatty liver disease (NAFLD) necessitates the utilization of
holistic approaches to investigate the network properties of NAFLD
pathogenesis. We describe the generation of the first mouse-specific
multi-tissue metabolic model, the SteatoNet, to investigate NAFLD-related
steatosis, by utilising an object-oriented modelling approach and incorporating
numerous hepatic metabolic pathways, their interaction with peripheral tissues
and hierarchical feedback regulation at the transcriptional and
post-translational level. The unique validated model is based on steady-state
analysis of ordinary differential equations that allows investigation of
systemic behaviour in the absence of estimated kinetic parameters. Model
analysis indicated the importance of the metabolic flux distribution parameter
in controlling metabolite concentration and identified critical focal points in
the network that may play a pivotal role in initiating NAFLD. Albeit requiring
experimental validation, the candidate triggers identified by SteatoNet provide
insight into the network properties of NAFLD and emphasize the promising scope
and potential of SteatoNet as a primary investigative tool for hypotheses
generation and exploring related metabolic diseases.