  In the companion paper of this set (Capitan and Cuesta, 2010) we have
developed a full analytical treatment of the model of species assembly
introduced in Capitan et al. (2009). This model is based on the construction of
an assembly graph containing all viable configurations of the community, and
the definition of a Markov chain whose transitions are the transformations of
communities by new species invasions. In the present paper we provide an
exhaustive numerical analysis of the model, describing the average time to the
recurrent state, the statistics of avalanches, and the dependence of the
results on the amount of available resource. Our results are based on the fact
that the Markov chain provides an asymptotic probability distribution for the
recurrent states, which can be used to obtain averages of observables as well
as the time variation of these magnitudes during succession, in an exact
manner. Since the absorption times into the recurrent set are found to be
comparable to the size of the system, the end state is quickly reached (in
units of the invasion time). Thus, the final ecosystem can be regarded as a
fluctuating complex system where species are continually replaced by newcomers
without ever leaving the set of recurrent patterns. The assembly graph is
dominated by pathways in which most invasions are accepted, triggering small
extinction avalanches. Through the assembly process, communities become less
resilient (e.g., have a higher return time to equilibrium) but become more
robust in terms of resistance against new invasions.


  The interplay of seasonality, the system's nonlinearities and intrinsic
stochasticity is studied for a seasonally forced
susceptible-exposed-infective-recovered stochastic model. The model is explored
in the parameter region that corresponds to childhood infectious diseases such
as measles. The power spectrum of the stochastic fluctuations around the
attractors of the deterministic system that describes the model in the
thermodynamic limit is computed analytically and validated by stochastic
simulations for large system sizes. Size effects are studied through additional
simulations. Other effects such as switching between coexisting attractors
induced by stochasticity often mentioned in the literature as playing an
important role in the dynamics of childhood infectious diseases are also
investigated. The main conclusion is that stochastic amplification, rather than
these effects, is the key ingredient to understand the observed incidence
patterns.


  Experimental and theoretical analyses of the lactate dehydrogenase and
glutathione reductase based enzymatic AND logic gates in which the enzymes and
their substrates serve as logic inputs are performed. These two systems are
examples of the novel, previously unexplored, class of biochemical logic gates
that illustrate potential biomedical applications of biochemical logic. They
are characterized by input concentrations at logic 0 and 1 states corresponding
to normal and abnormal physiological conditions. Our analysis shows that the
logic gates under investigation have similar noise characteristics. Both
significantly amplify random noise present in inputs, however we establish that
for realistic widths of the input noise distributions, it is still possible to
differentiate between the logic 0 and 1 states of the output. This indicates
that reliable detection of abnormal biomedical conditions is indeed possible
with such enzyme-logic systems.


  A model for the transmission of dengue disease is presented. It consists of
eight mutually-exclusive compartments representing the human and vector
dynamics. It also includes a control parameter (adulticide spray) in order to
combat the mosquito. The model presents three possible equilibria: two
disease-free equilibria (DFE) --- where humans, with or without mosquitoes,
live without the disease --- and another endemic equilibrium (EE). In the
literature it has been proved that a DFE is locally asymptotically stable,
whenever a certain epidemiological threshold, known as the basic reproduction
number, is less than one. We show that if a minimum level of insecticide is
applied, then it is possible to maintain the basic reproduction number below
unity. A case study, using data of the outbreak that occured in 2009 in Cape
Verde, is presented.


  The folding pathway and rate coefficients of the folding of a knotted protein
are calculated for a potential energy function with minimal energetic
frustration. A kinetic transition network is constructed using the discrete
path sampling approach, and the resulting potential energy surface is
visualized by constructing disconnectivity graphs. Owing to topological
constraints, the low-lying portion of the landscape consists of three distinct
regions, corresponding to the native knotted state and to configurations where
either the N- or C-terminus is not yet folded into the knot. The fastest
folding pathways from denatured states exhibit early formation of the
N-terminus portion of the knot and a rate-determining step where the C-terminus
is incorporated. The low-lying minima with the N-terminus knotted and the
C-terminus free therefore constitute an off-pathway intermediate for this
model. The insertion of both the N- and C-termini into the knot occur late in
the folding process, creating large energy barriers that are the rate limiting
steps in the folding process. When compared to other protein folding proteins
of a similar length, this system folds over six orders of magnitude more
slowly.


  We present a simple model of protein folding dynamics that captures key
qualitative elements recently seen in all-atom simulations. The goals of this
theory are to serve as a simple formalism for gaining deeper insight into the
physical properties seen in detailed simulations as well as to serve as a model
to easily compare why these simulations suggest a different kinetic mechanism
than previous simple models. Specifically, we find that non-native contacts
play a key role in determining the mechanism, which can shift dramatically as
the energetic strength of non-native interactions is changed. For protein-like
non-native interactions, our model finds that the native state is a kinetic
hub, connecting the strength of relevant interactions directly to the nature of
folding kinetics.


  We study large fluctuations in evolutionary games belonging to the
coordination and anti-coordination classes. The dynamics of these games,
modeling cooperation dilemmas, is characterized by a coexistence fixed point
separating two absorbing states. We are particularly interested in the problem
of fixation that refers to the possibility that a few mutants take over the
entire population. Here, the fixation phenomenon is induced by large
fluctuations and is investigated by a semi-classical WKB
(Wentzel-Kramers-Brillouin) theory generalized to treat stochastic systems
possessing multiple absorbing states. Importantly, this method allows us to
analyze the combined influence of selection and random fluctuations on the
evolutionary dynamics \textit{beyond} the weak selection limit often considered
in previous works. We accurately compute, including pre-exponential factors,
the probability distribution function in the long-lived coexistence state and
the mean fixation time necessary for a few mutants to take over the entire
population in anti-coordination games, and also the fixation probability in the
coordination class. Our analytical results compare excellently with extensive
numerical simulations. Furthermore, we demonstrate that our treatment is
superior to the Fokker-Planck approximation when the selection intensity is
finite.


  Gene duplications are one of major primary driving forces for evolutionary
novelty. We took population genetics models of genes duplicate to study how
evolutionary forces acting during the fixation of mutant allele at duplicate
loci. We study the fixation time of mutant allele at duplicate loci under
double null recessive model (DNR) and haploinsufficient model (HI). And we also
investigate how selection coefficients with other evolutionary force influence
the fixation frequency of mutant allele at duplicate loci. Our results suggest
that the selection plays a role in the evolutionary fate of duplicate genes,
and tight linkage would help the mutant allele preserved at duplicate loci. Our
theoretical simulation agree with the genomics data analysis result well, that
selection, rather than drift, plays a important role in the establishment of
duplicate loci, and recombination have a great opportunity to be acted upon
selection.


  The approach for calculation of the mode intensities of DNA conformational
vibrations in the Raman spectra is developed. It is based on the valence-optic
theory and the model for description of conformational vibrations of DNA with
counterions. The calculations for Na- and Cs-DNA low-frequency Raman spectra
show that the vibrations of DNA backbone chains near 15 cm-1 have the greatest
intensity. In the spectrum of Na-DNA at frequency range upper than 40 cm-1 the
modes of H-bond stretching in base pairs have the greatest intensities, while
the modes of ion-phosphate vibrations have the lowest intensity. In Cs-DNA
spectra at this frequency range the mode of ion-phosphate vibrations is
prominent. Its intensity is much higher than the intensities of Na-DNA modes of
this spectra range. Other modes of Cs-DNA have much lower intensities than in
the case of Na-DNA. The comparison of our calculations with the experimental
data shows that developed approach gives the understanding of the sensitivity
of DNA low-frequency Raman bands to the neutralization of the double helix by
light and heavy counterions.


  We study the evolution of cooperation in spatial public goods games where,
besides the classical strategies of cooperation (C) and defection (D), we
consider punishing cooperators (PC) or punishing defectors (PD) as an
additional strategy. Using a minimalist modeling approach, our goal is to
separately clarify and identify the consequences of the two punishing
strategies. Since punishment is costly, punishing strategies loose the
evolutionary competition in case of well-mixed interactions. When spatial
interactions are taken into account, however, the outcome can be strikingly
different, and cooperation may spread. The underlying mechanism depends on the
character of the punishment strategy. In case of cooperating punishers,
increasing the fine results in a rising cooperation level. In contrast, in the
presence of the PD strategy, the phase diagram exhibits a reentrant transition
as the fine is increased. Accordingly, the level of cooperation shows a
non-monotonous dependence on the fine. Remarkably, punishing strategies can
spread in both cases, but based on largely different mechanisms, which depend
on the cooperativeness (or not) of punishers.


  The hypothesis of a companion object (Nemesis) orbiting the Sun was motivated
by the claim of a terrestrial extinction periodicity, thought to be mediated by
comet showers. The orbit of a distant companion to the Sun is expected to be
perturbed by the Galactic tidal field and encounters with passing stars, which
will induce variation in the period. We examine the evidence for the previously
proposed periodicity, using two modern, greatly improved paleontological
datasets of fossil biodiversity. We find that there is a narrow peak at 27 My
in the cross-spectrum of extinction intensity time series between these
independent datasets. This periodicity extends over a time period nearly twice
that for which it was originally noted. An excess of extinction events are
associated with this periodicity at 99% confidence. In this sense we confirm
the originally noted feature in the time series for extinction. However, we
find that it displays extremely regular timing for about 0.5 Gy. The regularity
of the timing compared with earlier calculations of orbital perturbation would
seem to exclude the Nemesis hypothesis as a causal factor.


  Components of biological systems interact with each other in order to carry
out vital cell functions. Such information can be used to improve estimation
and inference, and to obtain better insights into the underlying cellular
mechanisms. Discovering regulatory interactions among genes is therefore an
important problem in systems biology. Whole-genome expression data over time
provides an opportunity to determine how the expression levels of genes are
affected by changes in transcription levels of other genes, and can therefore
be used to discover regulatory interactions among genes.
  In this paper, we propose a novel penalization method, called truncating
lasso, for estimation of causal relationships from time-course gene expression
data. The proposed penalty can correctly determine the order of the underlying
time series, and improves the performance of the lasso-type estimators.
Moreover, the resulting estimate provides information on the time lag between
activation of transcription factors and their effects on regulated genes. We
provide an efficient algorithm for estimation of model parameters, and show
that the proposed method can consistently discover causal relationships in the
large $p$, small $n$ setting. The performance of the proposed model is
evaluated favorably in simulated, as well as real, data examples. The proposed
truncating lasso method is implemented in the R-package grangerTlasso and is
available at http://www.stat.lsa.umich.edu/~shojaie.


  This study investigates the effect of competition between individuals on
population dynamics when they compete for different resources during different
seasons or during different growth stages. Individuals are assumed to compete
for a single resource during each of these periods according to one of the
following competition types: scramble, contest, or an intermediate between the
two. The effect of two successive competition periods is determined to be
expressed by simple relations on products of two "transition matrices" for
various sets of competition types for the two periods. In particular, for the
scramble and contest competition combination, results vary widely depending on
the order of the two competition types. Furthermore, the stability properties
of derived population models as well as the effect of more than two successive
competition periods are discussed.


  Spatial evolutionary games model individuals who are distributed in a spatial
domain and update their strategies upon playing a normal form game with their
neighbors. We derive integro-differential equations as deterministic
approximations of the microscopic updating stochastic processes. This
generalizes the known mean-field ordinary differential equations and provide a
powerful tool to investigate the spatial effects in populations evolution. The
deterministic equations allow to identify many interesting features of the
evolution of strategy profiles in a population, such as standing and traveling
waves, and pattern formation, especially in replicator-type evolutions.


  Studies into the statistical properties of biological networks have led to
important biological insights, such as the presence of hubs and hierarchical
modularity. There is also a growing interest in studying the statistical
properties of networks in the context of cancer genomics. However, relatively
little is known as to what network features differ between the cancer and
normal cell physiologies, or between different cancer cell phenotypes. Based on
the observation that frequent genomic alterations underlie a more aggressive
cancer phenotype, we asked if such an effect could be detectable as an increase
in the randomness of local gene expression patterns. Using a breast cancer gene
expression data set and a model network of protein interactions we derive
constrained weighted networks defined by a stochastic information flux matrix
reflecting expression correlations between interacting proteins. Based on this
stochastic matrix we propose and compute an entropy measure that quantifies the
degree of randomness in the local pattern of information flux around single
genes. By comparing the local entropies in the non-metastatic versus metastatic
breast cancer networks, we here show that breast cancers that metastasize are
characterised by a small yet significant increase in the degree of randomness
of local expression patterns. We validate this result in three additional
breast cancer expression data sets and demonstrate that local entropy better
characterises the metastatic phenotype than other non-entropy based measures.
We show that increases in entropy can be used to identify genes and signalling
pathways implicated in breast cancer metastasis. Further exploration of such
integrated cancer expression and protein interaction networks will therefore be
a fruitful endeavour.


  Super-spreading events for infectious diseases occur when some infected
individuals infect more than the average number of secondary cases. Several
super-spreading individuals have been identified for the 2003 outbreak of
severe acute respiratory syndrome (SARS). We develop a model for
super-spreading events of infectious diseases, which is based on the outbreak
of SARS. Using this model we describe two methods for estimating the parameters
of the model, which we demonstrate with the small-scale SARS outbreak at the
Amoy Gardens, Hong Kong, and the large-scale outbreak in the entire Hong Kong
Special Administrative Region. One method is based on parameters calculated for
the classical susceptible - infected - removed (SIR) disease model. The second
is based on parameter estimates found in the literature. Using the parameters
calculated for the SIR model, our model predicts an outcome similar to that for
the SIR model. On the other hand, using parameter estimates from SARS
literature our model predicts a much more serious epidemic.


  Much has been studied on the spreading speed and traveling wave solutions for
cooperative reaction-diffusion systems. In this paper, we shall establish the
spreading speed for a large class of non-cooperative reaction-diffusion systems
and characterize the spreading speed as the slowest speed of a family of
non-constant traveling wave solutions. As an application, our results are
applied to a partially cooperative system describing interactions between
ungulates and grass.


  We present a step towards the metabolome-wide computational inference of
cellular metabolic reaction networks from metabolic profiling data, such as
mass spectrometry. The reconstruction is based on identification of irreducible
statistical interactions among the metabolite activities using the ARACNE
reverse-engineering algorithm and on constraining possible metabolic
transformations to satisfy the conservation of mass. The resulting algorithms
are validated on synthetic data from an abridged computational model of
Escherichia coli metabolism. Precision rates upwards of 50% are routinely
observed for identification of full metabolic reactions, and recalls upwards of
20% are also seen.


  We propose a theoretical formalism, molecular finite automata (MFA), to
describe individual proteins as rule-based computing machines. The MFA
formalism provides a framework for modeling individual protein behaviors and
systems-level dynamics via construction of programmable and executable
machines. Models specified within this formalism explicitly represent the
context-sensitive dynamics of individual proteins driven by external inputs and
represent protein-protein interactions as synchronized machine
reconfigurations. Both deterministic and stochastic simulations can be applied
to quantitatively compute the dynamics of MFA models. We apply the MFA
formalism to model and simulate a simple example of a signal transduction
system that involves a MAP kinase cascade and a scaffold protein.


  Males are often the "sicker" sex with male biased parasitism found in a
taxonomically diverse range of species. There is considerable interest in the
processes that could underlie the evolution of sex-biased parasitism. Mating
system differences along with differences in lifespan may play a key role. We
examine whether these factors are likely to lead to male-biased parasitism
through natural selection taking into account the critical role that ecological
feedbacks play in the evolution of defence. We use a host-parasite model with
two-sexes and the techniques of adaptive dynamics to investigate how mating
system and sexual differences in competitive ability and longevity can select
for a bias in the rates of parasitism. Male-biased parasitism is selected for
when males have a shorter average lifespan or when males are subject to greater
competition for resources. Male-biased parasitism evolves as a consequence of
sexual differences in life history that produce a greater proportion of
susceptible females than males and therefore reduce the cost of avoiding
parasitism in males. Different mating systems such as monogamy, polygamy or
polyandry did not produce a bias in parasitism through these ecological
feedbacks but may accentuate an existing bias.


  We examine birth--death processes with state dependent transition
probabilities and at least one absorbing boundary. In evolution, this describes
selection acting on two different types in a finite population where
reproductive events occur successively. If the two types have equal fitness the
system performs a random walk. If one type has a fitness advantage it is
favored by selection, which introduces a bias (asymmetry) in the transition
probabilities. How long does it take until advantageous mutants have invaded
and taken over? Surprisingly, we find that the average time of such a process
can increase, even if the mutant type always has a fitness advantage. We
discuss this finding for the Moran process and develop a simplified model which
allows a more intuitive understanding. We show that this effect can occur for
weak but non--vanishing bias (selection) in the state dependent transition
rates and infer the scaling with system size. We also address the Wright-Fisher
model commonly used in population genetics, which shows that this stochastic
slowdown is not restricted to birth-death processes.


  The large reservoir of antibiotic resistant bacteria in raw and treated water
supplies is a matter of public health concern. Currently, the National
Antimicrobial Resistance Monitoring Systems, a collaborative effort of the
Centers for Disease Control, the US Department of Agriculture, and the US Food
and Drug Administration, does not monitor antimicrobial resistance in surface
waters. Given the serious nature of antibiotic resistance in clinical settings,
and the likelihood that antibiotic resistant bacteria can be transmitted to
humans from large environmental reservoirs via drinking water, explanations for
the distribution of antibiotic resistant bacteria and tools for studying this
distribution must be found. Here we focus on mathematical modeling of
cultivable bacteria in a river, which will be used to study the distribution of
antibiotic resistant bacteria in the environment. We consider both antibiotic
resistant and non-antibiotic resistant bacteria in the model, and, taking into
account the strong correlation between land use and antibiotic resistant
bacteria in rivers, we include a function for the influx of bacteria into the
river from the shore. We simulate the model for two different time scales and
show that if there is too many bacteria from the land entering the river, the
river entirely fills with antibiotic resistant bacteria, while less frequent
influxes allows time for the bacteria to lose the antibiotic resistant gene.
This mathematically verifies that reduction in antibiotic use near the banks of
rivers, will reduce the counts of antibiotic resistant bacteria in rivers.


  Studying the topology of so-called real networks, that is networks obtained
from sociological or biological data for instance, has become a major field of
interest in the last decade. One way to deal with it is to consider that
networks are built from small functional units called motifs, which can be
found by looking for small subgraphs whose numbers of occurrences in the whole
network are surprisingly high. In this article, we propose to define motifs
through a local overrepresentation in the network and develop a statistic to
detect them without relying on simulations. We then illustrate the performance
of our procedure on simulated and real data, recovering already known
biologically relevant motifs. Moreover, we explain how our method gives some
information about the respective roles of the vertices in a motif.


  There are a great many proteins that localize to and collectively generate
curvature in biological fluid membranes. We study changes in the topology of
fluid membranes due to the presence of highly anisotropic, curvature-inducing
proteins. Generically, we find a surprisingly rich phase diagram with phases of
both positive and negative Gaussian curvature. As a concrete example modeled on
experiments, we find that a lamellar phase in a negative Gaussian curvature
regime exhibits a propensity to form screw dislocations of definite burgers
scalar but of both chirality. The induced curvature depends strongly on the
membrane rigidity, suggesting membrane composition can be a factor regulating
membrane sculpting to to curvature-inducing proteins.


  Both community ecology and conservation biology seek further understanding of
factors governing the advance of an invasive species. We model biological
invasion as an individual-based, stochastic process on a two-dimensional
landscape. An ecologically superior invader and a resident species compete for
space preemptively. Our general model includes the basic contact process and a
variant of the Eden model as special cases. We employ the concept of a
"roughened" front to quantify effects of discreteness and stochasticity on
invasion; we emphasize the probability distribution of the front-runner's
relative position. That is, we analyze the location of the most advanced
invader as the extreme deviation about the front's mean position. We find that
a class of models with different assumptions about neighborhood interactions
exhibit universal characteristics. That is, key features of the invasion
dynamics span a class of models, independently of locally detailed demographic
rules. Our results integrate theories of invasive spatial growth and generate
novel hypotheses linking habitat or landscape size (length of the invading
front) to invasion velocity, and to the relative position of the most advanced
invader.


  Our understanding of the physics of biological molecules, such as proteins
and DNA, is limited because the approximations we usually apply to model inert
materials are not in general applicable to soft, chemically inhomogeneous
systems. The configurational complexity of biomolecules means the entropic
contribution to the free energy is a significant factor in their behaviour,
requiring detailed dynamical calculations to fully evaluate. Computer
simulations capable of taking all interatomic interactions into account are
therefore vital. However, even with the best current supercomputing facilities,
we are unable to capture enough of the most interesting aspects of their
behaviour to properly understand how they work. This limits our ability to
design new molecules, to treat diseases, for example. Progress in biomolecular
simulation depends crucially on increasing the computing power available.
Faster classical computers are in the pipeline, but these provide only
incremental improvements. Quantum computing offers the possibility of
performing huge numbers of calculations in parallel, when it becomes available.
We discuss the current open questions in biomolecular simulation, how these
might be addressed using quantum computation and speculate on the future
importance of quantum-assisted biomolecular modelling.


  The intrinsic stochasticity of gene expression can lead to large variations
in protein levels across a population of cells. To explain this variability,
different sources of mRNA fluctuations ('Poisson' and 'Telegraph' processes)
have been proposed in stochastic models of gene expression. Both Poisson and
Telegraph scenario models explain experimental observations of noise in protein
levels in terms of 'bursts' of protein expression. Correspondingly, there is
considerable interest in establishing relations between burst and steady-state
protein distributions for general stochastic models of gene expression. In this
work, we address this issue by considering a mapping between stochastic models
of gene expression and problems of interest in queueing theory. By applying a
general theorem from queueing theory, Little's Law, we derive exact relations
which connect burst and steady-state distribution means for models with
arbitrary waiting-time distributions for arrival and degradation of mRNAs and
proteins. The derived relations have implications for approaches to quantify
the degree of transcriptional bursting and hence to discriminate between
different sources of intrinsic noise in gene expression. To illustrate this, we
consider a model for regulation of protein expression bursts by small RNAs. For
a broad range of parameters, we derive analytical expressions (validated by
stochastic simulations) for the mean protein levels as the levels of regulatory
small RNAs are varied. The results obtained show that the degree of
transcriptional bursting can, in principle, be determined from changes in mean
steady-state protein levels for general stochastic models of gene expression.


  A common survival strategy of microorganisms subjected to stress involves the
generation of phenotypic heterogeneity in the isogenic microbial population
enabling a subset of the population to survive under stress. In a recent study,
a mycobacterial population of M. smegmatis was shown to develop phenotypic
heterogeneity under nutrient depletion. The observed heterogeneity is in the
form of a bimodal distribution of the expression levels of the Green
Fluorescent Protein (GFP) as reporter with the gfp fused to the promoter of the
rel gene. The stringent response pathway is initiated in the subpopulation with
high rel activity.In the present study, we characterize quantitatively the
single cell promoter activity of the three key genes, namely, mprA, sigE and
rel, in the stringent response pathway with gfp as the reporter. The origin of
bimodality in the GFP distribution lies in two stable expression states, i.e.,
bistability. We develop a theoretical model to study the dynamics of the
stringent response pathway. The model incorporates a recently proposed
mechanism of bistability based on positive feedback and cell growth retardation
due to protein synthesis. Based on flow cytometry data, we establish that the
distribution of GFP levels in the mycobacterial population at any point of time
is a linear superposition of two invariant distributions, one Gaussian and the
other lognormal, with only the coefficients in the linear combination depending
on time. This allows us to use a binning algorithm and determine the time
variation of the mean protein level, the fraction of cells in a subpopulation
and also the coefficient of variation, a measure of gene expression noise.The
results of the theoretical model along with a comprehensive analysis of the
flow cytometry data provide definitive evidence for the coexistence of two
subpopulations with overlapping protein distributions.


  Recent work suggests a link between chiral asymmetry in the amino acid
iso-valine extracted from the Murchison meteorite and the extent of hydrous
alteration. We present the results of neutron scattering experiments on an
exchanged, 1-dimensionally ordered n-propyl ammonium vermiculite clay. The
vermiculite gel has a (001) d-spacing of order 5nm at the temperature and
concentration of the experiments and the d-spacing responds sensitively to
changes in concentration, temperature and electronic environment. The data show
that isothermal addition of D-histidine or L-histidine solutions produces
shifts in the d-spacing that are different for each enantiomer. This chiral
specificity is of interest for the question of whether clays could have played
an important role in the origin of biohomochirality.


  Transport through nano-channels plays an important role in many biological
processes and industrial applications. Gaining insights into the functioning of
biological transport processes and the design of man-made nano-devices requires
an understanding of the basic physics of such transport. A simple exclusion
process has proven to be very useful in ex- plaining the properties of several
artificial and biological nano-channels. It is particularly useful for modeling
the influence of inter-particle interactions on transport characteristics. In
this paper, we explore several models of the exclusion process using a mean
field approach and computer simulations. We examine the effects of crowding
inside the channel and its immediate vicinity on the mean flux and the
transport times of single molecules. Finally, we discuss the robustness of the
theory's predictions with respect to the crucial characteristics of the
hindered diffusion in nano-channels that need to be included in the model.


  Site-specific recombination on supercoiled circular DNA molecules can yield a
variety of knots and catenanes. Twist knots are some of the most common
conformations of these products and they can act as substrates for further
rounds of site-specific recombination. They are also one of the simplest
families of knots and catenanes. Yet, our systematic understanding of their
implication in DNA and important cellular processes like site-specific
recombination is very limited. Here we present a topological model of
site-specific recombination characterising all possible products of this
reaction on twist knot substrates, extending previous work of Buck and Flapan.
We illustrate how to use our model to examine previously uncharacterised
experimental data. We also show how our model can help determine the sequence
of products in multiple rounds of processive recombination and distinguish
between products of processive and distributive recombination. This model
studies generic site- specific recombination on arbitrary twist knot
substrates, a subject for which there is limited global understanding. We also
provide a systematic method of applying our model to a variety of different
recombination systems.


  Systems with the coexistence of different stable attractors are widely
exploited in systems biology in order to suitably model the differentiating
processes arising in living cells. In order to describe genetic regulatory
networks several deterministic models based on systems of nonlinear ordinary
differential equations have been proposed. Few studies have been developed to
characterize how either an external input or the coupling can drive systems
with different coexisting states. For the sake of simplicity, in this
manuscript we focus on systems belonging to the class of radial isochron clocks
that exhibits hard excitation, in order to investigate their complex dynamics,
local and global bifurcations arising in presence of constant external inputs.
In particular the occurrence of saddle node on limit cycle bifurcations is
detected.


  Integrins are allosteric cell adhesion receptors that control many important
processes, including cell migration, proliferation, and apoptosis. Ligand
binding activates integrins by stabilizing an integrin conformation with
separated cytoplasmic tails, thus enabling the binding of proteins that mediate
cytoplasmic signaling. Experiments demonstrate a high sensitivity of integrin
signaling to ligand density and this has been accounted mainly to avidity
effects. Based on experimental data we have developed a quantitative Monte
Carlo model for integrin signal initiation. We show that within the
physiological ligand density range avidity effects cannot explain the
sensitivity of cellular signaling to small changes in ligand density. Src
kinases are among the first proteins to be activated, possibly by trans
auto-phosphorylation. We calculate the extent of integrin and ligand clustering
as well as the speed and extent of Src kinase activation by trans
auto-phosphorylation or direct binding at different experimentally monitored
ligand densities. We find that the experimentally observed ligand density
dependency can be reproduced if Src kinases are activated by trans
auto-phosphorylation or some other mechanism limits integrin-dependent Src
kinase activation. We propose that Src kinase and thus cell activation by trans
auto-phosphorylation may provide a mechanism to enable ligand-density dependent
responses at physiological ligand densities. The capacity to detect small
differences in ligand density at a ligand density that is large enough to
permit cell adhesion is likely to be important for haptotaxis.


  How do living cells achieve sufficient abundances of functional protein
complexes while minimizing promiscuous non-functional interactions? Here we
study this problem using a first-principle model of the cell whose phenotypic
traits are directly determined from its genome through biophysical properties
of protein structures and binding interactions in crowded cellular environment.
The model cell includes three independent prototypical pathways, whose
topologies of Protein-Protein Interaction (PPI) sub-networks are different, but
whose contributions to the cell fitness are equal. Model cells evolve through
genotypic mutations and phenotypic protein copy number variations. We found a
strong relationship between evolved physical-chemical properties of protein
interactions and their abundances due to a "frustration" effect: strengthening
of functional interactions brings about hydrophobic interfaces, which make
proteins prone to promiscuous binding. The balancing act is achieved by
lowering concentrations of hub proteins while raising solubilities and
abundances of functional monomers. Based on these principles we generated and
analyzed a possible realization of the proteome-wide PPI network in yeast. In
this simulation we found that high-throughput affinity capture - mass
spectroscopy experiments can detect functional interactions with high fidelity
only for high abundance proteins while missing most interactions for low
abundance proteins.


  Acute myeloid leukemia (AML) involves a block in terminal differentiation of
the myeloid lineage and uncontrolled proliferation of a progenitor state. Using
phorbol myristate acetate (PMA), it is possible to overcome this block in THP-1
cells (an M5-AML containing the MLL-MLLT3 fusion), resulting in differentiation
to an adherent monocytic phenotype. As part of FANTOM4, we used microarrays to
identify 23 microRNAs that are regulated by PMA. We identify four PMA-induced
micro- RNAs (mir-155, mir-222, mir-424 and mir-503) that when overexpressed
cause cell-cycle arrest and partial differentiation and when used in
combination induce additional changes not seen by any individual microRNA. We
further characterize these prodifferentiative microRNAs and show that mir-155
and mir-222 induce G2 arrest and apoptosis, respectively. We find mir-424 and
mir-503 are derived from a polycistronic precursor mir-424-503 that is under
repression by the MLL-MLLT3 leukemogenic fusion. Both of these microRNAs
directly target cell-cycle regulators and induce G1 cell-cycle arrest when
overexpressed in THP-1. We also find that the pro-differentiative mir-424 and
mir-503 downregulate the anti-differentiative mir-9 by targeting a site in its
primary transcript. Our study highlights the combinatorial effects of multiple
microRNAs within cellular systems.


  The melting behavior of long, heterogeneous DNA chains is examined within the
framework of the nonlinear lattice dynamics based Peyrard-Bishop-Dauxois (PBD)
model. Data for the pBR322 plasmid and the complete T7 phage have been used to
obtain model fits and determine parameter dependence on salt content. Melting
curves predicted for the complete fd phage and the Y1 and Y2 fragments of the
$\phi$X174 phage without any adjustable parameters are in good agreement with
experiment. The calculated probabilities for single base-pair opening are
consistent with values obtained from imino proton exchange experiments.


  A dynamic mathematical model has been developed and validated to describe the
synthesis of pectate lyases (Pels), the major virulence factors in Dickeya
dadantii. This work focuses on the simultaneous modeling of the metabolic
degradation of pectin by Pel enzymes and the genetic regulation of pel genes by
2-keto-3-deoxygluconate (KDG), a catabolite product of pectin which inactivates
KdgR, one of the main repressors of pel genes. This modeling scheme takes into
account the fact that the system is composed of two time-varying compartments:
the extracellular medium, where Pel enzymes cleave pectin into oligomers, and
the bacterial cytoplasm where, after internalization, oligomers are converted
to KDG. Using the quasi-stationary state approximations, the model consists of
some nonlinear differential equations for which most of the parameters could be
estimated from the literature or from independent experiments. The few
remaining unknown parameters were obtained by fitting the model equations
against a set of Pel activity data. Model predictions were verified by
measuring the time courses of bacterial growth, Pel production, pel mRNA
accumulation and pectin consumption under various growth conditions. This work
reveals that pectin is almost totally consumed before the burst of Pel
production. This paradoxical behaviour can be interpreted as an evolutionary
strategy to control the diffusion process so that as soon as a small amount of
pectin is detected by the bacteria in its surroundings it anticipates more
pectin to come. The model also predicts the possibility of bistable steady
states in the presence of constant pectin compounds.


  We propose a degree-based coarse graining approach that not just accelerates
the evaluation of dynamics on complex networks, but also satisfies the
consistency conditions for both equilibrium statistical distributions and
nonequilibrium dynamical flows. For the Ising model and
susceptible-infected-susceptible epidemic model, we introduce these required
conditions explicitly and further prove that they are satisfied by our
coarse-grained network construction within the annealed network approximation.
Finally, we numerically show that the phase transitions and fluctuations on the
coarse-grained network are all in good agreements with those on the original
one.


  Biological networks of interacting agents exhibit similar topological
properties for a wide range of scales, from cellular to ecological levels,
suggesting the existence of a common evolutionary origin. A general
evolutionary mechanism based on global stability has been proposed recently [J
I Perotti, O V Billoni, F A Tamarit, D R Chialvo, S A Cannas, Phys. Rev. Lett.
103, 108701 (2009)]. This mechanism is incorporated into a model of a growing
network of interacting agents in which each new agent's membership in the
network is determined by the agent's effect on the network's global stability.
We show that, out of this stability constraint, several topological properties
observed in biological networks emerge in a self organized manner. The
influence of the stability selection mechanism on the dynamics associated to
the resulting network is analyzed as well.


  Conventional wisdom holds that influenza A and B are such genetically
dissimilar viruses that infection with one cannot confer cross-immunity to the
other. However, our examination of the records of the past 25 influenza seasons
in the U.S. reveals that almost every time there is an early and severe
influenza A outbreak, the annual influenza B epidemic is almost entirely
suppressed (and is never suppressed otherwise). Temporary broad-spectrum (aka
"heterologous") immunity in the aftermath of influenza infection is the most
direct explanation for this phenomenon. We find a remarkably weak degree of
temporary cross-immunity is needed to explain these patterns, and that indeed
influenza B provides an ideal setting for the observation of heterologous
immune effects outside of the carefully controlled environment of a laboratory.


  Understanding how monomeric proteins fold under in vitro conditions is
crucial to describing their functions in the cellular context. Significant
advances both in theory and experiments have resulted in a conceptual framework
for describing the folding mechanisms of globular proteins. The experimental
data and theoretical methods have revealed the multifaceted character of
proteins. Proteins exhibit universal features that can be determined using only
the number of amino acid residues (N) and polymer concepts. The sizes of
proteins in the denatured and folded states, cooperativity of the folding
transition, dispersions in the melting temperatures at the residue level, and
time scales of folding are to a large extent determined by N. The consequences
of finite N especially on how individual residues order upon folding depends on
the topology of the folded states. Such intricate details can be predicted
using the Molecular Transfer Model that combines simulations with measured
transfer free energies of protein building blocks from water to the desired
concentration of the denaturant. By watching one molecule fold at a time, using
single molecule methods, the validity of the theoretically anticipated
heterogeneity in the folding routes, and the N-dependent time scales for the
three stages in the approach to the native state have been established. Despite
the successes of theory, of which only a few examples are documented here, we
conclude that much remains to be done to solve the "protein folding problem" in
the broadest sense.


  We study the role of demographic fluctuations in typical endemics as
exemplified by the stochastic SIRS model. The birth-death master equation of
the model is simulated using exact numerics and analysed within the linear
noise approximation. The endemic fixed point is unstable to internal
demographic noise, and leads to sustained oscillations. This is ensured when
the eigenvalues ($\lambda$) of the linearised drift matrix are complex, which
in turn, is possible only if detailed balance is violated. In the oscillatory
state, the phases decorrelate asymptotically, distinguishing such oscillations
from those produced by external periodic forcing. These so-called quasicycles
are of sufficient strength to be detected reliably only when the ratio
$|Im(\lambda)/Re(\lambda)|$ is of order unity. The coherence or regularity of
these oscillations show a maximum as a function of population size, an effect
known variously as stochastic coherence or coherence resonance. We find that
stochastic coherence can be simply understood as resulting from a non-monotonic
variation of $|Im(\lambda)/Re(\lambda)|$ with population size. Thus, within the
linear noise approximation, stochastic coherence can be predicted from a purely
deterministic analysis. The non-normality of the linearised drift matrix,
associated with the violation of detailed balance, leads to enhanced
fluctuations in the population amplitudes.


  We review the existing mathematical models which describe physicochemical
mechanisms capable of producing a symmetry-breaking transition to a state in
which one chirality dominates the other. A new model is proposed, with the aim
of elucidating the fundamental processes at work in the crystal grinding
systems of Viedma [Phys Rev Lett 94, 065504, (2005)] and Noorduin [J Am Chem
Soc 130, 1158, (2008)]. We simplify the model as far as possible to uncover the
fundamental competitive process which causes the symmetry-breaking, and analyse
other simplifications which might be expected to show symmetry-breaking.


  We present a differential equation-based mathematical model of nectar
foraging by the honey bee Apis mellifera. The model focuses on two behavioural
classes; nectar foragers and nectar receivers. Results generated from the model
are used to demonstrate how different classes within a collective can
collaborate to combine information and produce finely tuned decisions through
simple interactions. In particular we show the importance of the `search time'
- the time a returning forager takes to find an available nectar receiver - in
restricting the forager population to a level consistent with colony-wide
needs.


  We discuss two different ways of chromosomes' and genomes' evolution.
Purifying selection dominates in large panmictic populations, where Mendelian
law of independent gene assortment is valid. If the populations are small,
recombination processes are not effective enough to ensure an independent
assortment of linked genes and larger clusters of genes could be inherited as
the genetic units. There are whole clusters of genes which tend to complement
in such conditions instead of single pairs of alleles like in the case of
purifying selection. Computer simulations have shown that switching in-between
complementation and purification strategies has a character of a phase
transition. This is also responsible for specific distribution of recombination
events observed along eukaryotic chromosomes - higher recombination rate is
observed in subtelomeric regions than in central parts of chromosomes - for
sympatric speciation and probably for non-monotonous relation between
reproduction potential and genetic distance between parents.


  Bacteria are the unseen majority on our planet, with millions of species and
comprising most of the living protoplasm. While current methods enable in-depth
study of a small number of communities, a simple tool for breadth studies of
bacterial population composition in a large number of samples is lacking. We
propose a novel approach for reconstruction of the composition of an unknown
mixture of bacteria using a single Sanger-sequencing reaction of the mixture.
This method is based on compressive sensing theory, which deals with
reconstruction of a sparse signal using a small number of measurements.
Utilizing the fact that in many cases each bacterial community is comprised of
a small subset of the known bacterial species, we show the feasibility of this
approach for determining the composition of a bacterial mixture. Using
simulations, we show that sequencing a few hundred base-pairs of the 16S rRNA
gene sequence may provide enough information for reconstruction of mixtures
containing tens of species, out of tens of thousands, even in the presence of
realistic measurement noise. Finally, we show initial promising results when
applying our method for the reconstruction of a toy experimental mixture with
five species. Our approach may have a potential for a practical and efficient
way for identifying bacterial species compositions in biological samples.


  In the past few years, overwhelming evidence has accrued that a high level of
expression of the protein neuroglobin protects neurons in vitro, in animal
models, and in humans, against cell death associated with hypoxic and amyloid
insult. However, until now, the exact mechanism of neuroglobin's protective
action has not been determined. Using cell biology and biochemical approaches
we demonstrate that neuroglobin inhibits the intrinsic pathway of apoptosis in
vitro and intervenes in activation of pro-caspase 9 by interaction with
cytochrome c. Using systems level information of the apoptotic signalling
reactions we have developed a quantitative model of neuroglobin inhibition of
apoptosis, which simulates neuroglobin blocking of apoptosome formation at a
single cell level. Furthermore, this model allows us to explore the effect of
neuroglobin in conditions not easily accessible to experimental study. We found
that the protection of neurons by neuroglobin is very concentration sensitive.
The impact of neuroglobin may arise from both its binding to cytochrome c and
its subsequent redox reaction, although the binding alone is sufficient to
block pro-caspase 9 activation. These data provides an explanation the action
of neuroglobin in the protection of nerve cells from unwanted apoptosis.


  BH3 mimetics have been proposed as new anticancer therapeutics. They target
anti-apoptotic Bcl-2 proteins, up-regulation of which has been implicated in
the resistance of many cancer cells, particularly leukemia and lymphoma cells,
to apoptosis. Using probabilistic computational modeling of the mitochondrial
pathway of apoptosis, verified by single-cell experimental observations, we
develop a model of Bcl-2 inhibition of apoptosis. Our results clarify how Bcl-2
imparts its anti-apoptotic role by increasing the time-to-death and
cell-to-cell variability. We also show that although the commitment to death is
highly impacted by differences in protein levels at the time of stimulation,
inherent stochastic fluctuations in apoptotic signaling are sufficient to
induce cell-to-cell variability and to allow single cells to escape death. This
study suggests that intrinsic cell-to-cell stochastic variability in apoptotic
signaling is sufficient to cause fractional killing of cancer cells after
exposure to BH3 mimetics. This is an unanticipated facet of cancer
chemoresistance.


  In this work we study a weak Prisoner\^as Dilemma game in which both
strategies and update rules are subjected to evolutionary pressure.
Interactions among agents are specified by complex topologies, and we consider
both homogeneous and heterogeneous situations. We consider deterministic and
stochastic update rules for the strategies, which in turn may consider single
links or full context when selecting agents to copy from. Our results indicate
that the co-evolutionary process preserves heterogeneous networks as a suitable
framework for the emergence of cooperation. Furthermore, on those networks, the
update rule leading to a larger fraction, which we call replicator dynamics, is
selected during co-evolution. On homogeneous networks we observe that even if
replicator dynamics turns out again to be the selected update rule, the
cooperation level is larger than on a fixed update rule framework. We conclude
that for a variety of topologies, the fact that the dynamics coevolves with the
strategies leads in general to more cooperation in the weak Prisoner's Dilemma
game.


  The adhesion zone of immune cells, the 'immunological synapse', exhibits
characteristic domains of receptor-ligand complexes. The domain formation is
likely caused by a length difference of the receptor-ligand complexes, and has
been investigated in experiments in which T cells adhere to supported membranes
with anchored ligands. For supported membranes with two types of anchored
ligands, MHCp and ICAM1, that bind to the receptors TCR and LFA1 in the cell
membrane, the coexistence of domains of TCR-MHCp and LFA1-ICAM1 complexes in
the cell adhesion zone has been observed for a wide range of ligand
concentrations and affinities. For supported membranes with long and short
ligands that bind to the same cell receptor CD2, in contrast, domain
coexistence has been observed for a rather narrow ratio of ligand
concentrations. In this article, we determine detailed phase diagrams for cells
adhering to supported membranes with a statistical-physical model of cell
adhesion. We find a characteristic difference between the adhesion scenarios in
which two types of ligands in a supported membrane bind (i) to the same cell
receptor or (ii) to two different cell receptors, which helps to explain the
experimental observations. Our phase diagrams fully include thermal shape
fluctuations of the cell membranes on nanometer scales, which lead to a
critical point for the domain formation and to a cooperative binding of the
receptors and ligands.


  The genetic code maps the sixty-four nucleotide triplets (codons) to twenty
amino-acids. While the biochemical details of this code were unraveled long
ago, its origin is still obscure. We review information-theoretic approaches to
the problem of the code's origin and discuss the results of a recent work that
treats the code in terms of an evolving, error-prone information channel. Our
model - which utilizes the rate-distortion theory of noisy communication
channels - suggests that the genetic code originated as a result of the
interplay of the three conflicting evolutionary forces: the needs for diverse
amino-acids, for error-tolerance and for minimal cost of resources. The
description of the code as an information channel allows us to mathematically
identify the fitness of the code and locate its emergence at a second-order
phase transition when the mapping of codons to amino-acids becomes nonrandom.
The noise in the channel brings about an error-graph, in which edges connect
codons that are likely to be confused. The emergence of the code is governed by
the topology of the error-graph, which determines the lowest modes of the
graph-Laplacian and is related to the map coloring problem.


  In this paper, we investigate a conjecture by von Haeseler concerning the
Maximum Parsimony method for phylogenetic estimation, which was published by
the Newton Institute in Cambridge on a list of open phylogenetic problems in
2007. This conjecture deals with the question whether Maximum Parsimony trees
are hereditary. The conjecture suggests that a Maximum Parsimony tree for a
particular (DNA) alignment necessarily has subtrees of all possible sizes which
are most parsimonious for the corresponding subalignments. We answer the
conjecture affirmatively for binary alignments on five taxa but also show how
to construct examples for which Maximum Parsimony trees are not hereditary.
Apart from showing that a most parsimonious tree cannot generally be reduced to
a most parsimonious tree on fewer taxa, we also show that compatible most
parsimonious quartets do not have to provide a most parsimonious supertree.
Last, we show that our results can be generalized to Maximum Likelihood for
certain nucleotide substitution models.


  We describe the structure of the graphs with the smallest average distance
and the largest average clustering given their order and size. There is usually
a unique graph with the largest average clustering, which at the same time has
the smallest possible average distance. In contrast, there are many graphs with
the same minimum average distance, ignoring their average clustering. The form
of these graphs is shown with analytical arguments. Finally, we measure the
sensitivity to rewiring of this architecture with respect to the clustering
coefficient, and we devise a method to make these networks more robust with
respect to vertex removal.


  The three dimensional structure of a protein is an outcome of the
interactions of its constituent amino acids in 3D space. Considering the amino
acids as nodes and the interactions among them as edges we have constructed and
analyzed protein contact networks at different length scales, long and
short-range. While long and short-range interactions are determined by the
positions of amino acids in primary chain, the contact networks are constructed
based on the 3D spatial distances of amino acids. We have further divided these
networks into sub-networks of hydrophobic, hydrophilic and charged residues.
Our analysis reveals that a significantly higher percentage of assortative
sub-clusters of long-range hydrophobic networks helps a protein in
communicating the necessary information for protein folding in one hand; on the
other hand the higher values of clustering coefficients of hydrophobic
sub-clusters play a major role in slowing down the process so that necessary
local and global stability can be achieved through intra connectivities of the
amino acid residues. Further, higher degrees of hydrophobic long-range
interactions suggest their greater role in protein folding and stability. The
small-range all amino acids networks have signature of hierarchy. The present
analysis with other evidences suggest that in a protein's 3D conformational
space, the growth of connectivity is not evolved either through preferential
attachment or through random connections; rather, it follows a specific
structural necessity based guiding principle - where some of the interactions
are primary while the others, generated as a consequence of these primary
interactions are secondary.


  Strategy changes are an essential part of evolutionary games. Here we
introduce a simple rule that, depending on the value of a single parameter $w$,
influences the selection of players that are considered as potential sources of
the new strategy. For positive $w$ players with high payoffs will be considered
more likely, while for negative $w$ the opposite holds. Setting $w$ equal to
zero returns the frequently adopted random selection of the opponent. We find
that increasing the probability of adopting the strategy from the fittest
player within reach, i.e. setting $w$ positive, promotes the evolution of
cooperation. The robustness of this observation is tested against different
levels of uncertainty in the strategy adoption process and for different
interaction network. Since the evolution to widespread defection is tightly
associated with cooperators having a lower fitness than defectors, the fact
that positive values of $w$ facilitate cooperation is quite surprising. We show
that the results can be explained by means of a negative feedback effect that
increases the vulnerability of defectors although initially increasing their
survivability. Moreover, we demonstrate that the introduction of $w$
effectively alters the interaction network and thus also the impact of
uncertainty by strategy adoptions on the evolution of cooperation.


  We discuss, in terms of rate-distortion theory, the fitness of molecular
codes as the problem of designing an optimal information channel. The fitness
is governed by an interplay between the cost and quality of the channel, which
induces smoothness in the code. By incorporating this code fitness into
population dynamics models, we suggest that the emergence and evolution of
molecular codes may be explained by simple channel design considerations.


  We consider systems where a cycle born via the Hopf bifurcation blows up to
infinity as a parameter ranges over a finite interval. Two examples
demonstrating this effect are presented: planar Lotka-Volterra type systems
with a competition-cooperation term and quasi-linear higher order equations.


  To perform recognition, molecules must locate and specifically bind their
targets within a noisy biochemical environment with many look-alikes. Molecular
recognition processes, especially the induced-fit mechanism, are known to
involve conformational changes. This arises a basic question: does molecular
recognition gain any advantage by such conformational changes? By introducing a
simple statistical-mechanics approach, we study the effect of conformation and
flexibility on the quality of recognition processes. Our model relates
specificity to the conformation of the participant molecules and thus suggests
a possible answer: Optimal specificity is achieved when the ligand is slightly
off target, that is a conformational mismatch between the ligand and its main
target improves the selectivity of the process. This indicates that
deformations upon binding serve as a conformational proofreading mechanism,
which may be selected for via evolution.


  Rubisco, probably the most abundant protein in the biosphere, performs an
essential part in the process of carbon fixation through photosynthesis thus
facilitating life on earth. Despite the significant effect that Rubisco has on
the fitness of plants and other photosynthetic organisms, this enzyme is known
to have a remarkably low catalytic rate and a tendency to confuse its
substrate, carbon dioxide, with oxygen. This apparent inefficiency is puzzling
and raises questions regarding the roles of evolution versus biochemical
constraints in shaping Rubisco. Here we examine these questions by analyzing
the measured kinetic parameters of Rubisco from various organisms in various
environments. The analysis presented here suggests that the evolution of
Rubisco is confined to an effectively one-dimensional landscape, which is
manifested in simple power law correlations between its kinetic parameters.
Within this one dimensional landscape, which may represent biochemical and
structural constraints, Rubisco appears to be tuned to the intracellular
environment in which it resides such that the net photosynthesis rate is nearly
optimal. Our analysis indicates that the specificity of Rubisco is not the main
determinant of its efficiency but rather the tradeoff between the carboxylation
velocity and CO2 affinity. As a result, the presence of oxygen has only
moderate effect on the optimal performance of Rubisco, which is determined
mostly by the local CO2 concentration. Rubisco appears as an experimentally
testable example for the evolution of proteins subject both to strong selection
pressure and to biochemical constraints which strongly confine the evolutionary
plasticity to a low dimensional landscape.


  Molecular recognition, which is essential in processing information in
biological systems, takes place in a crowded noisy biochemical environment and
requires the recognition of a specific target within a background of various
similar competing molecules. We consider molecular recognition as a
transmission of information via a noisy channel and use this analogy to gain
insights on the optimal, or fittest, molecular recognizer. We focus on the
optimal structural properties of the molecules such as flexibility and
conformation. We show that conformational changes upon binding, which often
occur during molecular recognition, may optimize the detection performance of
the recognizer. We thus suggest a generic design principle termed
'conformational proofreading' in which deformation enhances detection. We
evaluate the optimal flexibility of the molecular recognizer, which is
analogous to the stochasticity in a decision unit. In some scenarios, a
flexible recognizer, i.e., a stochastic decision unit, performs better than a
rigid, deterministic one. As a biological example, we discuss conformational
changes during homologous recombination, the process of genetic exchange
between two DNA strands.


  The function of the organism hinges on the performance of its
information-processing networks, which convey information via molecular
recognition. Many paths within these networks utilize molecular codebooks, such
as the genetic code, to translate information written in one class of molecules
into another molecular "language" . The present paper examines the emergence
and evolution of molecular codes in terms of rate-distortion theory and reviews
recent results of this approach. We discuss how the biological problem of
maximizing the fitness of an organism by optimizing its molecular coding
machinery is equivalent to the communication engineering problem of designing
an optimal information channel. The fitness of a molecular code takes into
account the interplay between the quality of the channel and the cost of
resources which the organism needs to invest in its construction and
maintenance. We analyze the dynamics of a population of organisms that compete
according to the fitness of their codes. The model suggests a generic mechanism
for the emergence of molecular codes as a phase transition in an information
channel. This mechanism is put into biological context and demonstrated in a
simple example.


  Transcription factor proteins bind specific DNA sequences to control the
expression of genes. They contain DNA binding domains which belong to several
super-families, each with a specific mechanism of DNA binding. The total number
of transcription factors encoded in a genome increases with the number of genes
in the genome. Here, we examined the number of transcription factors from each
super-family in diverse organisms.
  We find that the number of transcription factors from most super-families
appears to be bounded. For example, the number of winged helix factors does not
generally exceed 300, even in very large genomes. The magnitude of the maximal
number of transcription factors from each super-family seems to correlate with
the number of DNA bases effectively recognized by the binding mechanism of that
super-family. Coding theory predicts that such upper bounds on the number of
transcription factors should exist, in order to minimize cross-binding errors
between transcription factors. This theory further predicts that factors with
similar binding sequences should tend to have similar biological effect, so
that errors based on mis-recognition are minimal. We present evidence that
transcription factors with similar binding sequences tend to regulate genes
with similar biological functions, supporting this prediction.
  The present study suggests limits on the transcription factor repertoire of
cells, and suggests coding constraints that might apply more generally to the
mapping between binding sites and biological function.


  The control of gene expression involves complex mechanisms that show large
variation in design. For example, genes can be turned on either by the binding
of an activator (positive control) or the unbinding of a repressor (negative
control). What determines the choice of mode of control for each gene? This
study proposes rules for gene regulation based on the assumption that free
regulatory sites are exposed to nonspecific binding errors, whereas sites bound
to their cognate regulators are protected from errors. Hence, the selected
mechanisms keep the sites bound to their designated regulators for most of the
time, thus minimizing fitness-reducing errors. This offers an explanation of
the empirically demonstrated Savageau demand rule: Genes that are needed often
in the natural environment tend to be regulated by activators, and rarely
needed genes tend to be regulated by repressors; in both cases, sites are bound
for most of the time, and errors are minimized. The fitness advantage of error
minimization appears to be readily selectable. The present approach can also
generate rules for multi-regulator systems. The error-minimization framework
raises several experimentally testable hypotheses. It may also apply to other
biological regulation systems, such as those involving protein-protein
interactions.


  Numerous biological functions-such as enzymatic catalysis, the immune
response system, and the DNA-protein regulatory network-rely on the ability of
molecules to specifically recognize target molecules within a large pool of
similar competitors in a noisy biochemical environment. Using the basic
framework of signal detection theory, we treat the molecular recognition
process as a signal detection problem and examine its overall performance.
Thus, we evaluate the optimal properties of a molecular recognizer in the
presence of competition and noise. Our analysis reveals that the optimal design
undergoes a "phase transition" as the structural properties of the molecules
and interaction energies between them vary. In one phase, the recognizer should
be complementary in structure to its target (like a lock and a key), while in
the other, conformational changes upon binding, which often accompany molecular
recognition, enhance recognition quality. Using this framework, the abundance
of conformational changes may be explained as a result of increasing the
fitness of the recognizer. Furthermore, this analysis may be used in future
design of artificial signal processing devices based on biomolecules.


  Single cell experiments of simple regulatory networks can markedly differ
from cell population experiments. Such differences arise from stochastic events
in individual cells that are averaged out in cell populations. For instance,
while individual cells may show sustained oscillations in the concentrations of
some proteins, such oscillations may appear damped in the population average.
In this paper we investigate the role of RNA stochastic fluctuations as a
leading force to produce a sustained excitatory behavior at the single cell
level. Opposed to some previous models, we build a fully stochastic model of a
negative feedback loop that explicitly takes into account the RNA stochastic
dynamics. We find that messenger RNA random fluctuations can be amplified
during translation and produce sustained pulses of protein expression.
Motivated by the recent appreciation of the importance of non--coding
regulatory RNAs in post--transcription regulation, we also consider the
possibility that a regulatory RNA transcript could bind to the messenger RNA
and repress translation. Our findings show that the regulatory transcript helps
reduce gene expression variability both at the single cell level and at the
cell population level.


  A large amount of population models use the concept of a carrying capacity.
Simulated populations are bounded by invoking finite resources through a
survival probability, commonly referred to as the Verhulst factor. The fact,
however, that resources are not easily accounted for in actual biological
systems makes the carrying capacity parameter ill-defined. Henceforth, we deem
it essential to consider cases for which the parameter is unnecessary. This
work demonstrates the possibility of Verhulst-free steady states using the
Penna aging model, with one semelparous birth per adult. Stable populations are
obtained by setting a mutation threshold that is higher than the reproduction
age.


  We investigate a population dynamics model that exhibits a Neimark Sacker
bifurcation with a period that is naturally close to 4. Beyond the bifurcation,
the period becomes soon locked at 4 due to a strong resonance, and a second
attractor of period 2 emerges, which coexists with the first attractor over a
considerable parameter range. A linear stability analysis and a numerical
investigation of the second attractor reveal that the bifurcations producing
the second attractor occur naturally in this type of system.


  Our main tenet argues that the primary role of positive thymic selection and
the resulting T cell population is the maintenance of a homeostatic equilibrium
with self MHC-self peptide complexes. The homeostatic T cell repertoire can
recognize infections non-specifically and this is an indirect (negative)
recognition: the whole homeostatic T cell population together "holds a mirror"
to the whole self, and any MHC-peptide complex that is "not reflected in the
mirror" can be perceived by surrounding homeostatic T cells as a signal of the
presence of a foreign entity. On the other hand, infection-specific T cell
clones arise in a different pathway in the periphery, do not enter the thymus,
and form a functionally different population. Here we summarize the basic
assumptions and consequences of a logic-based new model, which differs from
conventional models in many respects.


  We developed a low-cost, high-throughput microbiome profiling method that
uses combinatorial sequence tags attached to PCR primers that amplify the rRNA
V6 region. Amplified PCR products are sequenced using an Illumina paired-end
protocol to generate millions of overlapping reads. Combinatorial sequence
tagging can be used to examine hundreds of samples with far fewer primers than
is required when sequence tags are incorporated at only a single end. The
number of reads generated permitted saturating or near-saturating analysis of
samples of the vaginal microbiome. The large number of reads al- lowed an
in-depth analysis of errors, and we found that PCR-induced errors composed the
vast majority of non-organism derived species variants, an ob- servation that
has significant implications for sequence clustering of similar high-throughput
data. We show that the short reads are sufficient to assign organisms to the
genus or species level in most cases. We suggest that this method will be
useful for the deep sequencing of any short nucleotide region that is
taxonomically informative; these include the V3, V5 regions of the bac- terial
16S rRNA genes and the eukaryotic V9 region that is gaining popularity for
sampling protist diversity.


  A model for the transmission of dengue disease is presented. It consists of
eight mutually-exclusive compartments representing the human and vector
dynamics. It also includes a control parameter (insecticide) in order to fight
the mosquitoes. The main goal of this work is to investigate the best way to
apply the control in order to effectively reduce the number of infected humans
and mosquitoes. A case study, using data of the outbreak that occurred in 2009
in Cape Verde, is presented.


  Using the Dominant Reaction Pathways method, we perform an ab-initio
quantum-mechanical simulation of a conformational transition of a peptide
chain. The method we propose makes it possible to investigate the
out-of-equilibrium dynamics of these systems, without resorting to an empirical
representation of the molecular force field. It also allows to study rare
transitions involving rearrangements in the electronic structure. By comparing
the results of the ab-initio simulation with those obtained employing a
standard force field, we discuss its capability to describe the non-equilibrium
dynamics of conformational transitions.


  Organisms modulate their fitness in heterogeneous environments by dispersing.
Prior work shows that there is selection against "unconditional" dispersal in
spatially heterogeneous environments. "Unconditional" means individuals
disperse at a rate independent of their location. We prove that if within-patch
fitness varies spatially and between two values temporally, then there is
selection for unconditional dispersal: any evolutionarily stable strategy (ESS)
or evolutionarily stable coalition (ESC) includes a dispersive phenotype.
Moreover, at this ESS or ESC, there is at least one sink patch (i.e. geometric
mean of fitness less than one) and no sources patches (i.e. geometric mean of
fitness greater than one). These results coupled with simulations suggest that
spatial-temporal heterogeneity due to abiotic forcing result in either an ESS
with a dispersive phenotype or an ESC with sedentary and dispersive phenotypes.
In contrast, spatial-temporal heterogeneity due to biotic interactions can
select for higher dispersal rates that ultimately spatially synchronize
population dynamics.


  With the recent developments in proteomic technologies, a complete human
proteome project (HPP) appears feasible for the first time. However, there is
still debate as to how it should be designed and what it should encompass. In
"proteomics speak", the debate revolves around the central question as to
whether a gene-centric or a protein-centric proteomics approach is the most
appropriate way forward. In this paper, we try to shed light on what these
definitions mean, how large-scale proteomics such as a HPP can insert into the
larger omics chorus, and what we can reasonably expect from a HPP in the way it
has been proposed so far.


  We consider a branching process with Poissonian immigration where individuals
have inheritable types. At rate theta, new individuals singly enter the total
population and start a new population which evolves like a supercritical,
homogeneous, binary Crump-Mode-Jagers process: individuals have i.i.d.
lifetimes durations (non necessarily exponential) during which they give birth
independently at constant rate b. First, using spine decomposition, we relax
previously known assumptions required for a.s. convergence of total population
size. Then, we consider three models of structured populations: either all
immigrants have a different type, or types are drawn in a discrete spectrum or
in a continuous spectrum. In each model, the vector (P_1,P_2,...) of relative
abundances of surviving families converges a.s. In the first model, the limit
is the GEM distribution with parameter theta/b.


  We investigate how the pattern of contacts between species in mutualistic
ecosystems is affected by the phylogenetic proximity between the species of
each guild. We develop several theoretical tools to measure that effect and we
use them to examine some real mutualistic sytems. We aim at establishing the
role of such proximity in the emergence of a nested pattern of contacts. We
conclude that although phylogenetic proximity is compatible with nestedness it
can not be claimed to determine it. We find that nestedness can instead be
attributed to a general rule by which species tend to behave as generalists
holding contacts with counterparts that already have a large number of
contacts. A nested ecosystem generated by this rule, shows high phylogenetic
diversity. This is to say, the counterparts of species having similar degrees
are not phylogenetic neighbours.


  Pedigrees are directed acyclic graphs that represent ancestral relationships
between individuals in a population. Based on a schematic recombination
process, we describe two simple Markov models for sequences evolving on
pedigrees - Model R (recombinations without mutations) and Model RM
(recombinations with mutations). For these models, we ask an identifiability
question: is it possible to construct a pedigree from the joint probability
distribution of extant sequences? We present partial identifiability results
for general pedigrees: we show that when the crossover probabilities are
sufficiently small, certain spanning subgraph sequences can be counted from the
joint distribution of extant sequences. We demonstrate how pedigrees that
earlier seemed difficult to distinguish are distinguished by counting their
spanning subgraph sequences.


  Smectite clays are hydrated layer silicates that, like micas, occur naturally
in abundance. Importantly, they have readily modifiable interlayer spaces that
provide excellent sites for nanochemistry. Vermiculite is one such smectite
clay and in the presence of small chain-length alkyl-NH3Cl ions, forms
sensitive, 1-D ordered model clay systems with expandable nano-pore inter-layer
regions. These inter-layers readily adsorb organic molecules. N-propyl NH3Cl
vermiculite clay gels were used to determine the adsorption of alanine, lysine
and histidine by chiral HPLC. The results show that during reaction with fresh
vermiculite interlayers, significant chiral enrichment of either L- and
D-enantiomers occurs depending on the amino acid. Chiral enrichment of the
supernatant solutions is up to about 1% per pass. In contrast, addition to clay
interlayers already reacted with amino acid solutions resulted in little or no
change in D/L ratio during the time of the experiment. Adsorption of small
amounts of amphiphilic organic molecules in clay inter-layers is known to
produce Layer-by-Layer or Langmuir-Blodgett films. Moreover atomistic
simulations show that self-organization of organic species in clay interlayers
is important. These non-centrosymmetric, chirally active nanofilms may cause
clays to act subsequently as chiral amplifiers, concentrating organic material
from dilute solution and having different adsorption energetics for D- and
L-enantiomers. The additional role of clays in RNA oligimerization already
postulated by Ferris and others, together with the need for the organization of
amphiphilic molecules and lipids noted by Szostak and others, suggests that
such chiral separation by clays in lagoonal environments at normal biological
temperatures might also have played a significant role in the origin of
biochirality.


  We obtain a Fokker-Planck equation describing experimental data on the
collective motion of locusts. The noise is of internal origin and due to the
discrete character and finite number of constituents of the swarm. The
stationary probability distribution shows a rich phenomenology including
non-monotonic behavior of several order/disorder transition indicators in noise
intensity. This complex behavior arises naturally as a result of the randomness
in the system. Its counterintuitive character challenges standard
interpretations of noise induced transitions and calls for an extension of this
theory in order to capture the behavior of certain classes of biologically
motivated models. Our results suggest that the collective switches of the
group's direction of motion might be due to a random ergodic effect and, as
such, they are inherent to group formation.


  The protein folding is regarded as a quantum transition between torsion
states on polypeptide chain. The deduction of the folding rate formula in our
previous studies is reviewed. The rate formula is generalized to the case of
frequency variation in folding. Then the following problems about the
application of the rate theory are discussed: 1) The unified theory on the
two-state and multi-state protein folding is given based on the concept of
quantum transition. 2) The relationship of folding and unfolding rates vs
denaturant concentration is studied. 3) The temperature dependence of folding
rate is deduced and the non-Arrhenius behaviors of temperature dependence are
interpreted in a natural way. 4) The inertial moment dependence of folding rate
is calculated based on the model of dynamical contact order and consistent
results are obtained by comparison with one-hundred-protein experimental
dataset. 5) The exergonic and endergonic foldings are distinguished through the
comparison between theoretical and experimental rates for each protein. The
ultrafast folding problem is viewed from the point of quantum folding theory
and a new folding speed limit is deduced from quantum uncertainty relation. And
finally, 6) since only the torsion-accessible states are manageable in the
present formulation of quantum transition how the set of torsion-accessible
states can be expanded by using statistical energy landscape approach is
discussed. All above discussions support the view that the protein folding is
essentially a quantum transition between conformational states.


  Ribosome is a molecular machine that polymerizes a protein where the sequence
of the amino acid residues, the monomers of the protein, is dictated by the
sequence of codons (triplets of nucleotides) on a messenger RNA (mRNA) that
serves as the template. The ribosome is a molecular motor that utilizes the
template mRNA strand also as the track. Thus, in each step the ribosome moves
forward by one codon and, simultaneously, elongates the protein by one amino
acid. We present a theoretical model that captures most of the main steps in
the mechano-chemical cycle of a ribosome. The stochastic movement of the
ribosome consists of an alternating sequence of pause and translocation; the
sum of the durations of a pause and the following translocation is the time of
dwell of the ribosome at the corresponding codon. We derive the analytical
expression for the distribution of the dwell times of a ribosome in our model.
Whereever experimental data are available, our theoretical predictions are
consistent with those results. We suggest appropriate experiments to test the
new predictions of our model, particularly, the effects of the quality control
mechanism of the ribosome and that of their crowding on the mRNA track.


  Markov models are extensively used in the analysis of molecular evolution. A
recent line of research suggests that pairs of proteins with functional and
physical interactions co-evolve with each other. Here, by analyzing hundreds of
orthologous sets of three fungi and their co-evolutionary relations, we
demonstrate that co-evolutionary assumption may violate the Markov assumption.
Our results encourage developing alternative probabilistic models for the cases
of extreme co-evolution.


  Ecosystems often undergo abrupt regime shifts in response to gradual external
changes. These shifts are theoretically understood as a regime switch between
alternative stable states of the ecosystem dynamical response to smooth changes
in external conditions. Usual models introduce nonlinearities in the
macroscopic dynamics of the ecosystem that lead to different stable attractors
among which the shift takes place. Here we propose an alternative explanation
of catastrophic regime shifts based on a recent model that pictures ecological
communities as systems in continuous fluctuation, according to certain
transition probabilities, between different micro-states in the phase space of
viable communities. We introduce a spontaneous extinction rate that accounts
for gradual changes in external conditions, and upon variations on this control
parameter the system undergoes a regime shift with similar features to those
previously reported. Under our microscopic viewpoint we recover the main
results obtained in previous theoretical and empirical work (anomalous
variance, hysteresis cycles, trophic cascades). The model predicts a gradual
loss of species in trophic levels from bottom to top near the transition. But
more importantly, the spectral analysis of the transition probability matrix
allows us to rigorously establish that we are observing the fingerprints, in a
finite size system, of a true phase transition driven by background
extinctions.


  We define a general mathematical framework for studying post-translational
modification processes under the assumption of mass action kinetics.


  Posttranslational modification of proteins is key in transmission of signals
in cells. Many signaling pathways contain several layers of modification cycles
that mediate and change the signal through the pathway. Here, we study a simple
signaling cascade consisting of n layers of modification cycles, such that the
modified protein of one layer acts as modifier in the next layer. Assuming
mass-action kinetics and taking the formation of intermediate complexes into
account, we show that the steady states are solutions to a polynomial in one
variable, and in fact that there is exactly one steady state for any given
total amounts of substrates and enzymes.
  We demonstrate that many steady state concentrations are related through
rational functions, which can be found recursively. For example,
stimulus-response curves arise as inverse functions to explicit rational
functions. We show that the stimulus-response curves of the modified substrates
are shifted to the left as we move down the cascade. Further, our approach
allows us to study enzyme competition, sequestration and how the steady state
changes in response to changes in the total amount of substrates.
  Our approach is essentially algebraic and follows recent trends in the study
of posttranslational modification systems.


  Homologous recombination plays a key role in generating genetic diversity,
while maintaining protein functionality. The mechanisms by which RecA enables a
single-stranded segment of DNA to recognize a homologous tract within a whole
genome are poorly understood. The scale by which homology recognition takes
place is of a few tens of base pairs, after which the quest for homology is
over. To study the mechanism of homology recognition, RecA-promoted homologous
recombination between short DNA oligomers with different degrees of heterology
was studied in vitro, using fluorescence resonant energy transfer. RecA can
detect single mismatches at the initial stages of recombination, and the
efficiency of recombination is strongly dependent on the location and
distribution of mismatches. Mismatches near the 5' end of the incoming strand
have a minute effect, whereas mismatches near the 3' end hinder strand exchange
dramatically. There is a characteristic DNA length above which the sensitivity
to heterology decreases sharply. Experiments with competitor sequences with
varying degrees of homology yield information about the process of homology
search and synapse lifetime. The exquisite sensitivity to mismatches and the
directionality in the exchange process support a mechanism for homology
recognition that can be modeled as a kinetic proofreading cascade.


  We present a comparison between Fe K-edge x-ray absorption spectra of
carbonmonoxy-myoglobin and its simulation based on density-functional theory
determination of the structure and vibrations and spectral simulation with
multiple-scattering theory. An excellent comparison is obtained for the main
part of the molecular structure without any structural fitting parameters. The
geometry of the CO ligand is reliably determined using a synergic approach to
data analysis. The methodology underlying this approach is expected to be
especially useful in similar situations in which high-resolution data for
structure and vibrations are available.


  In this paper we show that living cells can multiplex biochemical signals,
i.e. transmit multiple signals through the same signaling pathway
simultaneously, and yet respond to them very specifically. We demonstrate how
two binary input signals can be encoded in the concentration of a common
signaling protein, which is then decoded such that each of the two output
signals provides reliable information about one corresponding input. Under
biologically relevant conditions the network can reach the maximum amount of
information that can be transmitted, which is 2 bits.


  We have studied nucleation dynamics of the Ising model in scale-free networks
with degree distribution $P(k)\sim k^{-\gamma}$ by using forward flux sampling
method, focusing on how the network topology would influence the nucleation
rate and pathway. For homogeneous nucleation, the new phase clusters grow from
those nodes with smaller degree, while the cluster sizes follow a power-law
distribution. Interestingly, we find that the nucleation rate $R_{Hom}$ decays
exponentially with the network size $N$, and accordingly the critical nucleus
size increases linearly with $N$, implying that homogeneous nucleation is not
relevant in the thermodynamic limit. These observations are robust to the
change of $\gamma$ and also present in random networks. In addition, we have
also studied the dynamics of heterogeneous nucleation, wherein $w$ impurities
are initially added, either to randomly selected nodes or to targeted ones with
largest degrees. We find that targeted impurities can enhance the nucleation
rate $R_{Het}$ much more sharply than random ones. Moreover, $\ln
(R_{Het}/R_{Hom})$ scales as $w^{\gamma-2/\gamma-1}$ and $w$ for targeted and
random impurities, respectively. A simple mean field analysis is also present
to qualitatively illustrate above simulation results.


  We propose a biophysical model of Escherichia coli that predicts growth rate
and an effective cellular composition from an effective, coarse-grained
representation of its genome. We assume that E. coli is in a state of balanced
exponential steadystate growth, growing in a temporally and spatially constant
environment, rich in resources. We apply this model to a series of past
measurements, where the growth rate and rRNA-to-protein ratio have been
measured for seven E. coli strains with an rRNA operon copy number ranging from
one to seven (the wild-type copy number). These experiments show that growth
rate markedly decreases for strains with fewer than six copies. Using the
model, we were able to reproduce these measurements. We show that the model
that best fits these data suggests that the volume fraction of macromolecules
inside E. coli is not fixed when the rRNA operon copy number is varied.
Moreover, the model predicts that increasing the copy number beyond seven
results in a cytoplasm densely packed with ribosomes and proteins. Assuming
that under such overcrowded conditions prolonged diffusion times tend to weaken
binding affinities, the model predicts that growth rate will not increase
substantially beyond the wild-type growth rate, as indicated by other
experiments. Our model therefore suggests that changing the rRNA operon copy
number of wild-type E. coli cells growing in a constant rich environment does
not substantially increase their growth rate. Other observations regarding
strains with an altered rRNA operon copy number, such as nucleoid compaction
and the rRNA operon feedback response, appear to be qualitatively consistent
with this model. In addition, we discuss possible design principles suggested
by the model and propose further experiments to test its validity.


  Migration between different habitats is ubiquitous among biological
populations. In this Letter, we study a simple quasispecies model for evolution
in two different habitats, with different fitness landscapes, coupled through
one-way migration. Our model applies to asexual, rapidly evolving organisms
such as microbes. Our key finding is a dynamical phase transition at a critical
value of the migration rate. The time to reach steady state diverges at this
critical migration rate. Above the transition, the population is dominated by
immigrants from the primary habitat. Below the transition, the genetic
composition of the population is highly non-trivial, with multiple coexisting
quasispecies which are not native to either habitat. Using results from
localization theory, we show that the critical migration rate may be very small
--- demonstrating that evolutionary outcomes can be very sensitive to even a
small amount of migration.


  The assembly of RecA on single-stranded DNA is measured and interpreted as a
stochastic finite-state machine that is able to discriminate fine differences
between sequences, a basic computational operation. RecA filaments efficiently
scan DNA sequence through a cascade of random nucleation and disassembly events
that is mechanistically similar to the dynamic instability of microtubules.
This iterative cascade is a multistage kinetic proofreading process that
amplifies minute differences, even a single base change. Our measurements
suggest that this stochastic Turing-like machine can compute certain integral
transforms.


  One of the major functions of RecA protein in the cell is to bind
single-stranded DNA exposed upon damage, thereby triggering the SOS repair
response.We present fluorescence anisotropy measurements at the binding onset,
showing enhanced DNA length discrimination induced by adenosine triphosphate
consumption. Our model explains the observed DNA length sensing as an outcome
of out-of equilibrium binding fluctuations, reminiscent of microtubule dynamic
instability. The cascade architecture of the binding fluctuations is a
generalization of the kinetic proofreading mechanism. Enhancement of precision
by an irreversible multistage pathway is a possible design principle in the
noisy biological environment.


  Individual Based Models (IBMs) and Agent Based Models (ABMs) have become
widely used tools to understand complex biological systems. However, general
methods of parameter inference for IBMs are not available. In this paper we
show that it is possible to address this problem with a traditional
likelihood-based approach, using an example of an IBM developed to describe the
spread of Chytridiomycosis in a population of frogs as a case study. We show
that if the IBM satisfies certain criteria we can find the likelihood (or
posterior) analytically, and use standard computational techniques, such as
Markov Chain Monte Carlo (MCMC), for parameter inference.