  There is growing evidence in favour of the temporal-coding hypothesis that
temporal correlation of neuronal discharges may serve to bind distributed
neuronal activity into unique representations and, in particular, that $\theta$
(3.5-7.5 Hz) and $\delta$ ($0.5<$3.5 Hz) oscillations facilitate information
coding. The $\theta$ and $\delta$ rhythms are shown to be involved in various
sleep stages, and during an{\ae}sthesia, and they undergo changes with the
depth of an{\ae}sthesia. We introduce a thalamocortical model of interacting
neuronal ensembles to describe phase relationships between $\theta$ and
$\delta$ oscillations, especially during deep and light an{\ae}sthesia.
Asymmetric and long range interactions among the thalamocortical neuronal
oscillators are taken into account. The model results are compared with the
experimental observations of Musizza et al. {\it J. Physiol. (London)} 2007
580:315-326. The $\delta$ and $\theta$ activities are found to be separately
generated and are governed by the thalamus and cortex respectively. Changes in
the degree of intra--ensemble and inter--ensemble synchrony imply that the
neuronal ensembles inhibit information coding during deep an{\ae}sthesia and
facilitate it during light an{\ae}sthesia.


  We study a simplified, solvable model of a fully-connected metabolic network
with constrained quenched disorder to mimic the conservation laws imposed by
stoichiometry on chemical reactions. Within a spin-glass type of approach, we
show that in presence of a conserved metabolic pool the flux state
corresponding to maximal growth is stationary independently of the pool size.
In addition, and at odds with the case of unconstrained networks, the volume of
optimal flux configurations remains finite, indicating that the frustration
imposed by stoichiometric constraints, while reducing growth capabilities,
confers robustness and flexibility to the system. These results have a clear
biological interpretation and provide a basic, fully analytical explanation to
features recently observed in real metabolic networks.


  The activation-polymerization-epimerization-depolymerization (APED) model of
Plasson et al. has recently been proposed as a mechanism for the evolution of
homochirality on prebiotic Earth. The dynamics of the APED model in
two-dimensional spatially-extended systems is investigated for various
realistic reaction parameters. It is found that the APED system allows for the
formation of isolated homochiral proto-domains surrounded by a racemate. A
diffusive slowdown of the APED network such as induced through tidal motion or
evaporating pools and lagoons leads to the stabilization of homochiral bounded
structures as expected in the first self-assembled protocells.


  The endemic pine, Pinus canariensis, forms one of the main forest ecosystems
in the Canary Islands. In this archipelago, pine forest is a mosaic of natural
stands (remnants of past forest overexploitation) and artificial stands planted
from the 1940's. The genetic makeup of the artificially regenerated forest is
of some concern. The use of reproductive material with uncontrolled origin or
from a reduced number of parental trees may produce stands ill adapted to local
conditions or unable to adapt in response to environmental change. The genetic
diversity within a transect of reforested stands connecting two natural forest
fragments has been studied with nuclear and chloroplast microsatellites. Little
genetic differentiation and similar levels of genetic diversity to the
surrounding natural stands were found for nuclear markers. However, chloroplast
microsatellites presented lower haplotype diversity in reforested stands, and
this may be a consequence of the lower effective population size of the
chloroplast genome, meaning chloroplast markers have a higher sensitivity to
bottlenecks. Understory natural regeneration within the reforestation was also
analysed to study gene flow from natural forest into artificial stands.
Estimates of immigration rate into artificially regenerated forest were high
(0.68-0.75), producing a significant increase of genetic diversity (both in
chloroplast and nuclear microsatellites), which indicates the capacity for
genetic recovery for P. canariensis reforestations surrounded by larger natural
stands.


  How might systematic left-right (L/R) asymmetry of the body plan originate in
multicellular animals (and plants)? Somehow, the microscopic handedness of
biological molecules must be brought up to macroscopic scales. Basic symmetry
principles suggest that the usual "biological" mechanisms -- diffusion and gene
regulation -- are insufficient to implement the "right-hand rule" defining a
third body axis from the other two. Instead, on the cellular level, "physical"
mechanisms (forces and collective dynamic states) are needed involving the long
stiff fibers of the cytoskeleton. I discuss some possible scenarios; only in
the case of vertebrate internal organs is the answer currently known (and even
that is in dispute).


  We consider a switching rate of a meta-stable reaction scheme, which includes
reactions with arbitrary steps, e.g. $kA\to(k+r)A$. Employing WKB
approximation, controlled by a large system size, we evaluate both the exponent
and the pre-exponential factor for the rate. The results are illustrated on a
number of examples.


  It is an increasingly important problem to study conditions on the structure
of a network that guarantee a given behavior for its underlying dynamical
system. In this paper we report that a Boolean network may fall within the
chaotic regime, even under the simultaneous assumption of several conditions
which in randomized studies have been separately shown to correlate with
ordered behavior. These properties include using at most two inputs for every
variable, using biased and canalyzing regulatory functions, and restricting the
number of negative feedback loops.
  We also prove for n-dimensional Boolean networks that if in addition the
number of outputs for each variable is bounded and there exist periodic orbits
of length c^n for c sufficiently close to 2, any network with these properties
must have a large proportion of variables that simply copy previous values of
other variables. Such systems share a structural similarity to a relatively
small Turing machine acting on one or several tapes.


  We study the large population limit of the Moran process, assuming
weak-selection, and for different scalings. Depending on the particular choice
of scalings, we obtain a continuous model that may highlight the genetic-drift
(neutral evolution) or natural selection; for one precise scaling, both effects
are present. For the scalings that take the genetic-drift into account, the
continuous model is given by a singular diffusion equation, together with two
conservation laws that are already present at the discrete level. For scalings
that take into account only natural selection, we obtain a hyperbolic singular
equation that embeds the Replicator Dynamics and satisfies only one
conservation law. The derivation is made in two steps: a formal one, where the
candidate limit model is obtained, and a rigorous one, where convergence of the
probability density is proved. Additional results on the fixation probabilities
are also presented.


  We have modified the sexual Penna model by introducing the fluctuating
environment and fluctuations representing physiological functions of
individuals. Additionally, we have introduced the mother care corresponding to
the protection against the deleterious influence of the environment, the
learning capacity of individuals corresponding to their immunity and adaptation
to the environment fluctuations and the other risk factors appearing at
puberty. Each of the above mentioned elements influences mainly the survival of
newborns and young individuals while genetic defects accumulated in the genomes
increase the noise of individuals and were responsible for higher mortality of
older individuals. All these modifications enable precise fitting the age
structure of the simulated populations to the age distributions of human
populations.


  Networks have in recent years emerged as an invaluable tool for describing
and quantifying complex systems in many branches of science. Recent studies
suggest that networks often exhibit hierarchical organization, where vertices
divide into groups that further subdivide into groups of groups, and so forth
over multiple scales. In many cases these groups are found to correspond to
known functional units, such as ecological niches in food webs, modules in
biochemical networks (protein interaction networks, metabolic networks, or
genetic regulatory networks), or communities in social networks. Here we
present a general technique for inferring hierarchical structure from network
data and demonstrate that the existence of hierarchy can simultaneously explain
and quantitatively reproduce many commonly observed topological properties of
networks, such as right-skewed degree distributions, high clustering
coefficients, and short path lengths. We further show that knowledge of
hierarchical structure can be used to predict missing connections in partially
known networks with high accuracy, and for more general network structures than
competing techniques. Taken together, our results suggest that hierarchy is a
central organizing principle of complex networks, capable of offering insight
into many network phenomena.


  The interior of cells is crowded thus making it important to assess the
effects of macromolecules on the folding of proteins. Using the Self-Organized
Polymer (SOP) model, which is a coarse-grained representation of polypeptide
chains, we probe the mechanical stability of Ubiquitin (Ub) monomers and
trimers ((Ub)$_3$) in the presence of monodisperse spherical crowding agents.
Crowding increases the volume fraction ($\Phi_c$)-dependent average force
($<f_u(\Phi_c)>$), relative to the value at $\Phi_c = 0$, needed to unfold Ub
and the polyprotein. For a given $\Phi_c$, the values of $<f_u(\Phi_c)>$
increase as the diameter ($\sigma_c$) of the crowding particles decreases. The
average unfolding force $<f_u(\Phi_c)>$ depends on the ratio $\frac{D}{R_g}$,
where $D \approx \sigma_c (\frac{\pi}{6 \Phi_c})^{{1/3}}$ with $R_g$ being the
radius of gyration of Ub (or (Ub)$_3$) in the unfolded state. Examination of
the unfolding pathways shows that, relative to $\Phi_c = 0$, crowding promotes
reassociation of ruptured secondary structural elements. Both the nature of the
unfolding pathways and $<f_u(\Phi_c)>$ for (Ub)$_3$ are altered in the presence
of crowding particles with the effect being most dramatic for the subunit that
unfolds last. We predict, based on SOP simulations and theoretical arguments,
that $<f_u(\Phi_c) > \sim \Phi_c^{\frac{1}{3\nu}}$, where $\nu$ is the Flory
exponent that describes the unfolded (random coil) state of the protein.


  Most of the theoretical models describing the translocation of a polymer
chain through a nanopore use the hypothesis that the polymer is always relaxed
during the complete process. In other words, models generally assume that the
characteristic relaxation time of the chain is small enough compared to the
translocation time that non-equilibrium molecular conformations can be ignored.
In this paper, we use Molecular Dynamics simulations to directly test this
hypothesis by looking at the escape time of unbiased polymer chains starting
with different initial conditions. We find that the translocation process is
not quite in equilibrium for the systems studied, even though the translocation
time tau is about 10 times larger than the relaxation time tau_r. Our most
striking result is the observation that the last half of the chain escapes in
less than ~12% of the total escape time, which implies that there is a large
acceleration of the chain at the end of its escape from the channel.


  Discrete ancestral problems arising in population genetics are investigated.
In the neutral case, the duality concept has proved of particular interest in
the understanding of backward in time ancestral process from the forward in
time branching population dynamics. We show that duality formulae still are of
great use when considering discrete non-neutral Wright-Fisher models. This
concerns a large class of non-neutral models with completely monotone (CM) bias
probabilities. We show that most classical bias probabilities used in the
genetics literature fall within this CM class or are amenable to it through
some `reciprocal mechanism' which we define. Next, using elementary algebra on
CM functions, some suggested novel evolutionary mechanisms of potential
interest are introduced and discussed.


  Small peptides are model molecules for the amino acid residues that are the
constituents of proteins. In any bottom-up approach to understand the
properties of these macromolecules essential in the functioning of every living
being, to correctly describe the conformational behaviour of small peptides
constitutes an unavoidable first step. In this work, we present an study of
several potential energy surfaces (PESs) of the model dipeptide HCO-L-Ala-NH2.
The PESs are calculated using the B3LYP density-functional theory (DFT) method,
with Dunning's basis sets cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, aug-cc-pVTZ, and
cc-pVQZ. These calculations, whose cost amounts to approximately 10 years of
computer time, allow us to study the basis set convergence of the B3LYP method
for this model peptide. Also, we compare the B3LYP PESs to a previous
computation at the MP2/6-311++G(2df,2pd) level, in order to assess their
accuracy with respect to a higher level reference. All data sets have been
analyzed according to a general framework which can be extended to other
complex problems and which captures the nearness concept in the space of model
chemistries (MCs).


  We developed a method for estimating the positional distribution of
transcription fac-tor (TF) binding sites using ChIP-chip data, and applied it
to recently published experiments on binding sites of nine TFs; OCT4, SOX2,
NANOG, HNF1A, HNF4A, HNF6, FOXA2, USF1 and CREB1. The data were obtained from a
genome-wide cov-erage of promoter regions from 8kb upstream of the
Transcription Start Site (TSS) to 2kb downstream. The number of target genes of
each TF ranges from few hundred to several thousand. We found that for each of
the nine TFs the estimated binding site distribution is closely approximated by
a mixture of two components: a narrow peak, localized within 300 base pairs
upstream of the TSS, and a distribution of almost uni-form density within the
tested region. Using Gene Ontology and Enrichment analysis, we were able to
associate (for each of the TFs studied) the target genes of both types of
binding with known biological processes. Most GO terms were enriched either
among the proximal targets or among those with a uniform distribution of
binding sites. For example, the three stemness-related TFs have several hundred
target genes that belong to "development" and "morphogenesis" whose binding
sites belong to the uniform dis-tribution.


  We introduce a novel method to screen the promoters of a set of genes with
shared biological function, against a precompiled library of motifs, and find
those motifs which are statistically over-represented in the gene set. The gene
sets were obtained from the functional Gene Ontology (GO) classification; for
each set and motif we optimized the sequence similarity score threshold,
independently for every location window (measured with respect to the TSS),
taking into account the location dependent nucleotide heterogeneity along the
promoters of the target genes. We performed a high throughput analysis,
searching the promoters (from 200bp downstream to 1000bp upstream the TSS), of
more than 8000 human and 23,000 mouse genes, for 134 functional Gene Ontology
classes and for 412 known DNA motifs. When combined with binding site and
location conservation between human and mouse, the method identifies with high
probability functional binding sites that regulate groups of biologically
related genes. We found many location-sensitive functional binding events and
showed that they clustered close to the TSS. Our method and findings were put
to several experimental tests. By allowing a "flexible" threshold and combining
our functional class and location specific search method with conservation
between human and mouse, we are able to identify reliably functional TF binding
sites. This is an essential step towards constructing regulatory networks and
elucidating the design principles that govern transcriptional regulation of
expression. The promoter region proximal to the TSS appears to be of central
importance for regulation of transcription in human and mouse, just as it is in
bacteria and yeast.


  The search for life elsewhere in the universe is a pivotal question in modern
science. However, to address whether life is common in the universe we must
first understand the likelihood of abiogenesis by studying the origin of life
on Earth. A key missing piece is the origin of biomolecular homochirality:
permeating almost every life-form on Earth is the presence of exclusively
levorotary amino acids and dextrorotary sugars. In this work we discuss recent
results suggesting that life's homochirality resulted from sequential chiral
symmetry breaking triggered by environmental events in a mechanism referred to
as punctuated chirality. Applying these arguments to other potentially
life-bearing platforms has significant implications for the search for
extraterrestrial life: we predict that a statistically representative sampling
of extraterrestrial stereochemistry will be racemic on average.


  The problem of inhibiting viral DNA ejection from bacteriophages by
multivalent counterions, specifically Mg$^{+2}$ counterions, is studied.
Experimentally, it is known that MgSO$_4$ salt has a strong and non-monotonic
effect on the amount of DNA ejected. There exists an optimal concentration at
which the minimum amount of DNA is ejected from the virus. At lower or higher
concentrations, more DNA is ejected from the capsid. We propose that this
phenomenon is the result of DNA overcharging by Mg$^{+2}$ multivalent
counterions. As Mg$^{+2}$ concentration increases from zero, the net charge of
DNA changes from negative to positive. The optimal inhibition corresponds to
the Mg$^{+2}$ concentration where DNA is neutral. At lower/higher
concentrations, DNA genome is charged. It prefers to be in solution to lower
its electrostatic self-energy, which consequently leads to an increase in DNA
ejection. By fitting our theory to available experimental data, the strength of
DNA$-$DNA short range attraction energies, mediated by Mg$^{+2}$, is found to
be $-$0.004 $k_BT$ per nucleotide base. This and other fitted parameters agree
well with known values from other experiments and computer simulations. The
parameters are also in aggreement qualitatively with values for tri- and
tetra-valent counterions.


  We examine genetic statistics used in the study of structured populations. In
a 1999 paper, Wakeley observed that the coalescent process associated with the
finite island model can be decomposed into a scattering phase and a collecting
phase. In this paper, we introduce a class of population structure models,
which we refer to as G/KC models, that obey such a decomposition. In a large
population, large sample limit we derive the distribution of the statistic Fst
for all G/KC models under the assumptions of strong or weak mutation. We show
that in the large population, large sample limit the island and two dimensional
stepping stone models are members of the G/KC class of models, thereby deriving
the distributions of Fst for these two well known models as a special case of a
general formula. We show that our analysis of Fst can be extended to an entire
class of genetic statistics, and we use our approach to examine homozygosity
measures. Our analysis uses coalescent based methods.


  The existence of beyond mean field quasi-cycle oscillations in a simple
spatial model of predator prey interactions is derived from a path integral
formalism. The results agree substantially with those obtained from analysis of
similar models using system size expansions of the master equation. In all of
these analyses, the discrete nature of predator prey populations and finite
size effects lead to persistent oscillations in time, but spatial patterns fail
to form. The path integral formalism goes beyond mean field theory and provides
a focus on individual realizations of the stochastic time evolution of
population not captured in the standard master equation approach.


  We present a self-contained theory for the mechanical response of DNA in
single molecule experiments. Our model is based on a 1D continuum description
of the DNA molecule and accounts both for its elasticity and for DNA-DNA
electrostatic interactions. We consider the classical loading geometry used in
experiments where one end of the molecule is attached to a substrate and the
other one is pulled by a tensile force and twisted by a given number of turns.
We focus on configurations relevant to the limit of a large number of turns,
which are made up of two phases, one with linear DNA and the other one with
superhelical DNA. The model takes into account thermal fluctuations in the
linear phase and electrostatic interactions in the superhelical phase. The
values of the torsional stress, of the supercoiling radius and angle, and key
features of the experimental extension-rotation curves, namely the slope of the
linear region and thermal buckling threshold, are predicted. They are found in
good agreement with experimental data.


  In evolutionary games the fitness of individuals is not constant but depends
on the relative abundance of the various strategies in the population. Here we
study general games among n strategies in populations of large but finite size.
We explore stochastic evolutionary dynamics under weak selection, but for any
mutation rate. We analyze the frequency dependent Moran process in well-mixed
populations, but almost identical results are found for the Wright-Fisher and
Pairwise Comparison processes. Surprisingly simple conditions specify whether a
strategy is more abundant on average than 1/n, or than another strategy, in the
mutation-selection equilibrium. We find one condition that holds for low
mutation rate and another condition that holds for high mutation rate. A linear
combination of these two conditions holds for any mutation rate. Our results
allow a complete characterization of n*n games in the limit of weak selection.


  The fruit fly D. melanogaster is increasingly used as a model organism for
studying acute hypoxia tolerance and for studying aging, but the interactions
between these two factors are not well known. Here we show that hypoxia
tolerance degrades with age in post-hypoxic recovery of whole-body movement,
heart rate and ATP content. We previously used 1H NMR metabolomics and a
constraint-based model of ATP-generating metabolism to discover the end
products of hypoxic metabolism in flies and generate hypotheses for the
biological mechanisms. We expand the reactions in the model using tissue- and
age-specific microarray data from the literature, and then examine metabolomic
profiles of thoraxes after 4 hours at 0.5% O2 and after 5 minutes of recovery
in 40- versus 3-day-old flies. Model simulations were constrained to fluxes
calculated from these data. Simulations suggest that the decreased ATP
production during reoxygenation seen in aging flies can be attributed to
reduced recovery of mitochondrial respiration pathways and concomitant
over-dependence on the acetate production pathway as an energy source.


  Many complex disease syndromes such as asthma consist of a large number of
highly related, rather than independent, clinical phenotypes, raising a new
technical challenge in identifying genetic variations associated simultaneously
with correlated traits. In this study, we propose a new statistical framework
called graph-guided fused lasso (GFlasso) to address this issue in a principled
way. Our approach explicitly represents the dependency structure among the
quantitative traits as a network, and leverages this trait network to encode
structured regularizations in a multivariate regression model over the
genotypes and traits, so that the genetic markers that jointly influence
subgroups of highly correlated traits can be detected with high sensitivity and
specificity. While most of the traditional methods examined each phenotype
independently and combined the results afterwards, our approach analyzes all of
the traits jointly in a single statistical method, and borrow information
across correlated phenotypes to discover the genetic markers that perturbe a
subset of correlated triats jointly rather than a single trait. Using simulated
datasets based on the HapMap consortium data and an asthma dataset, we compare
the performance of our method with the single-marker analysis, and other sparse
regression methods such as the ridge regression and the lasso that do not use
any structural information in the traits. Our results show that there is a
significant advantage in detecting the true causal SNPs when we incorporate the
correlation pattern in traits using our proposed methods.


  Results from molecular systematics and comparative developmental genetics
changed the picture of metazoan and especially bilaterian radiation. According
to this new animal phylogeny (introduced by Adoutte et al. 1999/2000), Grobbens
(1908) widely favoured protostome-deuterostome division of the Bilateria can be
upheld, but only with major rearrangements within these superphyla. On the
cladogenetic level, the Protostomia are split into two unexpected subgroups,
the Lophotrochozoa and Ecdysozoa. The deuterostomes are split into the
subgroups Chordata and Ambulacraria, which is not novel since Grobben (1908)
introduced the Deuterostomia in this way (together with the Chaetognatha as a
third line). However, many details of the new deuterostome phylogeny do not fit
traditional, morphology-based reconstructions. As a consequence, three
relatively unexpected proposals for early deuterostome evolution are favoured
today: An ambulacraria-scenario, a xenoturbellid-scenario, and a
chordate-scenario. The first two proposals are often discussed in the
literature, while the chordate-scenario is almost completely neglected.
Therefore, the paper presented focuses on the chordate scenario, i.e. the
hypothesis of an acrania-like ur-deuterostomian. It is argued that the
acrania-hypothesis is clearly preferable when biomechanic options of a
polysegmented, hydroskeletal body plan are taken into account. The so called
hydroskeleton hypothesis, rooted in the work of W. F. Gutmann, is the most
detailed anagenetic scenario which depicts an acrania-like ur-deuterostome.


  Complex systems are often modeled as Boolean networks in attempts to capture
their logical structure and reveal its dynamical consequences. Approximating
the dynamics of continuous variables by discrete values and Boolean logic gates
may, however, introduce dynamical possibilities that are not accessible to the
original system. We show that large random networks of variables coupled
through continuous transfer functions often fail to exhibit the complex
dynamics of corresponding Boolean models in the disordered (chaotic) regime,
even when each individual function appears to be a good candidate for Boolean
idealization. A suitably modified Boolean theory explains the behavior of
systems in which information does not propagate faithfully down certain chains
of nodes. Model networks incorporating calculated or directly measured transfer
functions reported in the literature on transcriptional regulation of genes are
described by the modified theory.


  We study bifurcations in a spatially extended nonlinear system representing
population dynamics with the help of analytic calculations based on the
time-independent Schr\"{o}dinger equation for a quantum particle subjected to a
uniform gravitational field. Despite the linear character of the
Schr\"{o}dinger equation, the result we obtain helps in the understanding of
the onset of abrupt transitions leading to extinction of biological
populations. The result is expressed in terms of Airy functions and sheds light
on the behavior of bacteria in a Petri dish as well as of large animals such as
rodents moving over a landscape.


  Considerable attention has been paid, in recent years, to the use of networks
in modeling complex real-world systems. Among the many dynamical processes
involving networks, propagation processes -- in which final state can be
obtained by studying the underlying network percolation properties -- have
raised formidable interest. In this paper, we present a bond percolation model
of multitype networks with an arbitrary joint degree distribution that allows
heterogeneity in the edge occupation probability. As previously demonstrated,
the multitype approach allows many non-trivial mixing patterns such as
assortativity and clustering between nodes. We derive a number of useful
statistical properties of multitype networks as well as a general phase
transition criterion. We also demonstrate that a number of previous models
based on probability generating functions are special cases of the proposed
formalism. We further show that the multitype approach, by naturally allowing
heterogeneity in the bond occupation probability, overcomes some of the
correlation issues encountered by previous models. We illustrate this point in
the context of contact network epidemiology.


  Many swimming bacteria are propelled by flagellar motors that stochastically
switch between the clockwise and counterclockwise rotation direction. While the
switching dynamics are one of the most important characteristics of flagellar
motors, the mechanisms that control switching are poorly understood. We present
a statistical-mechanical model of the flagellar rotary motor, which consists of
a number of stator proteins that drive the rotation of a ring of rotor
proteins, which in turn drives the rotation of a flagellar filament. At the
heart of our model is the assumption that the rotor protein complex can exist
in two conformational states corresponding to the two respective rotation
directions, and that switching between these states depends on interactions
with the stator proteins. This naturally couples the switching dynamics to the
rotation dynamics, making the switch sensitive to torque and speed. Another key
element of our model is that after a switching event, it takes time for the
load to build up, due to polymorphic transitions of the filament. Our model
predicts that this slow relaxation dynamics of the filament, in combination
with the load dependence of the switching frequency, leads to a characteristic
switching time, in agreement with recent observations.


  An agent-based model of population dynamics is presented. The model has as
its expected behaviour the population dynamics of the equation-based Webworld
model, within which large communities of species can be grown on evolutionary
time scales. Such communities can be used in the agent-based model without
disrupting the food web structure, and hence a unified model of evolutionary
time and individual-based dynamics can be realised. Individuals encounter
potential prey, and optimal foraging strategies are arrived at through natural
selection.


  The growth of Aeromonas hydrophila and aerobic natural flora (APC) on
gilthead seabream surfaces was evaluated during the refrigerated storage (21
days). The related growth curves were compared with those obtained by a
conventional third order predictive model obtaining a low agreement between
observed and predicted data (Root Mean Squared Error = 1.77 for Aeromonas
hydrophila and 0.64 for APC). The Lotka-Volterra interspecific competition
model was used in order to calculate the degree of interaction between the two
bacterial populations (\beta_{Ah/APC} and \beta{APC/Ah}, respectively, the
interspecific competition coefficients of APC on Aeromonas hydrophila and
vice-versa). Afterwards, the Lotka-Volterra equations were applied as tertiary
predictive model, taking into account, simultaneously, the environmental
fluctuations and the bacterial interspecific competition. This approach allowed
to obtain a best fitting to the observed mean growth curves with a Root Mean
Squared Error of 0.09 for Aeromonas hydrophila and 0.28 for APC. Finally,
authors carry out some considerations about the necessary use of competitive
models in the context of the new trends in predictive microbiology.


  The evolutionary dynamics of a system of cancerous cells in a model of
chronic myeloid leukemia (CML) is investigated by a statistical approach.
Cancer progression is explored by applying a Monte Carlo method to simulate the
stochastic behavior of cell reproduction and death in a population of blood
cells which can experience genetic mutations. In CML front line therapy is
represented by the tyrosine kinase inhibitor imatinib which strongly affects
the reproduction of leukemic cells only. In this work, we analyze the effects
of a targeted therapy on the evolutionary dynamics of normal, first-mutant and
cancerous cell populations. Several scenarios of the evolutionary dynamics of
imatinib-treated leukemic cells are described as a consequence of the efficacy
of the different modeled therapies. We show how the patient response to the
therapy changes when an high value of the mutation rate from healthy to
cancerous cells is present. Our results are in agreement with clinical
observations. Unfortunately, development of resistance to imatinib is observed
in a proportion of patients, whose blood cells are characterized by an
increasing number of genetic alterations. We find that the occurrence of
resistance to the therapy can be related to a progressive increase of
deleterious mutations.


  Many biological and artificial transport channels function without direct
input of metabolic energy during a transport event and without structural
rearrangements involving transitions from a 'closed' to an 'open' state.
Nevertheless, such channels are able to maintain efficient and selective
transport. It has been proposed that attractive interactions between the
transported molecules and the channel can increase the transport efficiency and
that the selectivity of such channels can be based on the strength of the
interaction of the specifically transported molecules with the channel. Herein,
we study the transport through narrow channels in a framework of a general
kinetic theory, which naturally incorporates multi-particle occupancy of the
channel and non-single-file transport. We study how the transport efficiency
and the probability of translocation through the channel are affected by
inter-particle interactions in the confined space inside the channel, and
establish conditions for selective transport. We compare the predictions of the
model with the available experimental data - and find good semi-quantitative
agreement. Finally, we discuss applications of the theory to the design of
artificial nano-molecular sieves.


  We introduce a quantitative measure of the capacity of a small biological
network to evolve. We apply our measure to a stochastic description of the
experimental setup of Guet et al. (Science 296:1466, 2002), treating chemical
inducers as functional inputs to biochemical networks and the expression of a
reporter gene as the functional output. We take an information-theoretic
approach, allowing the system to set parameters that optimize signal processing
ability, thus enumerating each network's highest-fidelity functions. We find
that all networks studied are highly evolvable by our measure, meaning that
change in function has little dependence on change in parameters. Moreover, we
find that each network's functions are connected by paths in the parameter
space along which information is not significantly lowered, meaning a network
may continuously change its functionality without losing it along the way. This
property further underscores the evolvability of the networks.


  Damaged or mismatched bases in DNA can be repaired by Base Excision Repair
(BER) enzymes that replace the defective base. Although the detailed molecular
structures of many BER enzymes are known, how they colocalize to lesions
remains unclear. One hypothesis involves charge transport (CT) along DNA
[Yavin, {\it et al.}, PNAS, {\bf 102}, 3546, (2005)]. In this CT mechanism,
electrons are released by recently adsorbed BER enzymes and travel along the
DNA. The electrons can scatter (by heterogeneities along the DNA) back to the
enzyme, destabilizing and knocking it off the DNA, or, they can be absorbed by
nearby lesions and guanine radicals. We develop a stochastic model to describe
the electron dynamics, and compute probabilities of electron capture by guanine
radicals and repair enzymes. We also calculate first passage times of electron
return, and ensemble-average these results over guanine radical distributions.
Our statistical results provide the rules that enable us to perform
implicit-electron Monte-Carlo simulations of repair enzyme binding and
redistribution near lesions. When lesions are electron absorbing, we show that
the CT mechanism suppresses wasteful buildup of enzymes along intact portions
of the DNA, maximizing enzyme concentration near lesions.


  Enveloped viruses enter host cells either through endocytosis, or by direct
fusion of the viral membrane envelope and the membrane of the host cell.
However, some viruses, such as HIV-1, HSV-1, and Epstein-Barr can enter a cell
through either mechanism, with the choice of pathway often a function of the
ambient physical chemical conditions, such as temperature and pH. We develop a
stochastic model that describes the entry process at the level of binding of
viral glycoprotein spikes to cell membrane receptors and coreceptors. In our
model, receptors attach the cell membrane to the viral membrane, while
subsequent binding of coreceptors enables fusion. The model quantifies the
competition between fusion and endocytotic entry pathways. Relative
probabilities for each pathway are computed numerically, as well as
analytically in the high viral spike density limit. We delineate parameter
regimes in which fusion or endocytosis is dominant. These parameters are
related to measurable and potentially controllable quantities such as membrane
bending rigidity and receptor, coreceptor, and viral spike densities.
Experimental implications of our mechanistic hypotheses are proposed and
discussed.


  The early steps of light response occur in the outer segment of rod and cone
photoreceptor. They involve the hydrolysis of cGMP, a soluble cyclic
nucleotide, that gates ionic channels located in the outer segment membrane. We
shall study here the rate by which cGMP is hydrolyzed by activated
phosphodiesterase (PDE). This process has been characterized experimentally by
two different rate constants $\beta_d$ and $\beta_{sub}$: $\beta_d$ accounts
for the effect of all spontaneously active PDE in the outer segment, and
$\beta_{sub}$ characterizes cGMP hydrolysis induced by a single light-activated
PDE. So far, no attempt has been made to derive the experimental values of
$\beta_d$ and $\beta_{sub}$ from a theoretical model, which is the goal of this
work. Using a model of diffusion in the confined rod geometry, we derive
analytical expressions for $\beta_d$ and $\beta_{sub}$ by calculating the flux
of cGMP molecules to an activated PDE site. We obtain the dependency of these
rate constants as a function of the outer segment geometry, the PDE activation
and deactivation rates and the aqueous cGMP diffusion constant. Our formulas
show good agreement with experimental measurements. Finally, we use our
derivation to model the time course of the cGMP concentration in a
transversally well stirred outer segment.


  The Automated Protein Structure Analysis (APSA) method, which describes the
protein backbone as a smooth line in 3-dimensional space and characterizes it
by curvature kappa and torsion tau as a function of arc length s, was applied
on 77 proteins to determine all secondary structural units via specific
kappa(s) and tau(s) patterns. A total of 533 alpha-helices and 644 beta-strands
were recognized by APSA, whereas DSSP gives 536 and 651 units, respectively.
Kinks and distortions were quantified and the boundaries (entry and exit) of
secondary structures were classified. Similarity between proteins can be easily
quantified using APSA, as was demonstrated for the roll architecture of
proteins ubiquitin and spinach ferridoxin. A twenty-by-twenty comparison of
all-alpha domains showed that the curvature-torsion patterns generated by APSA
provide an accurate and meaningful similarity measurement for secondary,
super-secondary, and tertiary protein structure. APSA is shown to accurately
reflect the conformation of the backbone effectively reducing 3-dimensional
structure information to 2-dimensional representations that are easy to
interpret and understand.


  Myosin-II's rod-like tail drives filament assembly with a head arrangement
that should generate equal and opposite contractile forces on actin--if one
assumes that the filament is a symmetric bipole. Self-assembled myosin
filaments are shown here to be asymmetric in physiological buffer based on
cross-correlated images from both atomic force microscopy (AFM) and total
internal reflection fluorescence (TIRF). Quantitative cross-correlation of
these orthogonal methods produces structural information unavailable to either
method alone in showing that fluorescence intensity along the filament length
is proportional to height. This implies that myosin heads form a shell around
the filament axis, consistent with F-actin binding. A motor density of ~50 -
100 heads/micron is further estimated but with an average of 32% more motors on
one half of any given filament compared to the other, regardless of length. A
purely entropic pyramidal lattice model is developed that qualitatively
captures this lack of length dependence and the distribution of filament
asymmetries. Such strongly asymmetric bipoles are likely to produce an
imbalanced contractile force in cells and in actin-myosin gels, and thereby
contribute to motility as well as cytoskeletal tension.


  We study the evolution of cooperation within the spatial prisoner's dilemma
game on a square lattice where a fraction of players $\mu$ can spread their
strategy more easily than the rest due to a predetermined larger teaching
capability. In addition, players characterized with the larger teaching
capability are allowed to temporarily link with distant opponents of the same
kind with probability $p$, thus introducing shortcut connections among the
distinguished. We show that these additional temporary connections are able to
sustain cooperation throughout the whole range of the temptation to defect.
Remarkably, we observe that as the temptation to defect increases the optimal
$\mu$ decreases, and moreover, only minute values of $p$ warrant the best
promotion of cooperation. Our study thus indicates that influential individuals
must be few and sparsely connected in order for cooperation to thrive in a
defection prone environment.


  The stability properties of models of spontaneous mirror symmetry breaking in
chemistry are characterized algebraically. The models considered here all
derive either from the Frank model or from autocatalysis with limited
enantioselectivity. Emphasis is given to identifying the critical parameter
controlling the chiral symmetry breaking transition from racemic to chiral
steady-state solutions. This parameter is identified in each case, and the
constraints on the chemical rate constants determined from dynamic stability
are derived.


  In fast growing bacteria, ribosomal RNA (rRNA) is required to be transcribed
at very high rates to sustain the high cellular demand on ribosome synthesis.
This results in dense traffic of RNA polymerases (RNAP). We developed a
stochastic model, integrating results of single-molecule and quantitative in
vivo studies of E. coli, to evaluate the quantitative effect of pausing,
termination, and antitermination on rRNA transcription. Our calculations reveal
that in dense RNAP traffic, spontaneous pausing of RNAP can lead to severe
"traffic jams", as manifested in the broad distribution of inter-RNAP distances
and can be a major factor limiting transcription and hence growth. Our results
suggest the suppression of these pauses by the ribosomal antitermination
complex to be essential at fast growth. Moreover, unsuppressed pausing by even
a few non-antiterminated RNAPs can already reduce transcription drastically
under dense traffic. However, the termination factor Rho can remove the
non-antiterminated RNAPs and restore fast transcription. The results thus
suggest an intriguing role by Rho to enhance rather than attenuate rRNA
transcription.


  The classification of life should be based upon the fundamental mechanism in
the evolution of life. We found that the global relationships among species
should be circular phylogeny, which is quite different from the common sense
based upon phylogenetic trees. The genealogical circles can be observed clearly
according to the analysis of protein length distributions of contemporary
species. Thus, we suggest that domains can be defined by distinguished
phylogenetic circles, which are global and stable characteristics of living
systems. The mechanism in genome size evolution has been clarified; hence main
component questions on C-value enigma can be explained. According to the
correlations and quasi-periodicity of protein length distributions, we can also
classify life into three domains.


  In this Supporting information we provide background information on our model
of the bacterial flagellar motor. We also derive the analytical solution of our
coarse-grained model of the switching dynamics and explain the hybrid
stochastic algorithm used for the simulations.


  We generalize the Poland-Scheraga (PS) model to the case of a circular DNA,
taking into account the twisting of the two strains around each other. Guided
by recent single-molecule experiments on DNA strands, we assume that the
torsional stress induced by denaturation enforces formation of supercoils whose
writhe absorbs the linking number expelled by the loops. Our model predicts
that, when the entropy parameter of a loop satisfies $c \le 2$, denaturation
transition does not take place. On the other hand for $c>2$ a first-order
denaturation transition is consistent with our model and may take place in the
actual system, as in the case with no supercoils. These results are in contrast
with other treatments of circular DNA melting where denaturation is assumed to
be accompanied by an increase in twist rather than writhe on the bound
segments.


  The protein backbone is described as a smooth curved and twisted line in
three-dimensional (3D) space and characterized by its curvature $\kappa(s)$ and
torsion $\tau(s)$ both expressed as a function of arc length s. It is shown
that the function $\tau(s)$ is sufficient to analyze the contributions of all
amino acids to the conformation of the protein backbone. The characteristic
peak and trough patterns of the $\tau(s)$ diagrams can be translated into a
16-letter code, which provides a rapid identification of helices, strands, and
turns, specifies entry and exit points of secondary structural units, and
determines their regularity in terms of distortions, kinks or breaks. Via
computer encoding, 3D protein structure is projected into a 1D string of
conformational letters. The 3D-1D-projection procedure represents an extension
of the Automated Protein Structure Analysis (APSA) method. APSA has been
applied to describe 155 super secondary structures from 94 proteins and to
compare results with Efimov's classification system of super secondary
structure. The applicability of the latter is demonstrated.


  Enzyme-mediated reactions may proceed through multiple intermediate
conformational states before creating a final product molecule, and one often
wishes to identify such intermediate structures from observations of the
product creation. In this paper, we address this problem by solving the
chemical master equations for various enzymatic reactions. We devise a
perturbation theory analogous to that used in quantum mechanics that allows us
to determine the first (<n>) and the second (variance) cumulants of the
distribution of created product molecules as a function of the substrate
concentration and the kinetic rates of the intermediate processes. The mean
product flux V=d<n>/dt (or "dose-response" curve) and the Fano factor
F=variance/<n> are both realistically measurable quantities, and while the mean
flux can often appear the same for different reaction types, the Fano factor
can be quite different. This suggests both qualitative and quantitative ways to
discriminate between different reaction schemes, and we explore this
possibility in the context of four sample multistep enzymatic reactions. We
argue that measuring both the mean flux and the Fano factor can not only
discriminate between reaction types, but can also provide some detailed
information about the internal, unobserved kinetic rates, and this can be done
without measuring single-molecule transition events.


  The equilibrium three dimensional shape of relatively short loops of DNA is
studied using an elastic model that takes into account anisotropy in bending
rigidities. Using a reasonable estimate for the anisotropy, it is found that
cyclized DNA with lengths that are not integer multiples of the pitch take on
nontrivial shapes that involve bending out of planes and formation of kinks.
The effect of sequence inhomogeneity on the shape of DNA is addressed, and
shown to enhance the geometrical features. These findings could shed some light
on the role of DNA conformation in protein--DNA interactions.


  At the one-cell stage, the C. elegans embryo becomes polarized along the
anterior-posterior axis. The PAR proteins form complementary anterior and
posterior domains in a dynamic process driven by cytoskeletal rearrangement.
Initially, the PAR proteins are uniformly distributed throughout the embryo.
Following a cue from fertilization, cortical actomyosin contracts towards the
anterior pole. PAR-3/PAR-6/PKC-3 (the anterior PAR proteins) become restricted
to the anterior cortex. PAR-1 and PAR-2 (the posterior PAR proteins) become
enriched in the posterior cortical region. We present a mathematical model of
this polarity establishment process, in which we take a novel approach to
combine reaction-diffusion dynamics of the PAR proteins coupled to a simple
model of actomyosin contraction. We show that known interactions between the
PAR proteins are sufficient to explain many aspects of the observed cortical
PAR dynamics in both wild-type and mutant embryos. However, cytoplasmic PAR
protein polarity, which is vital for generating daughter cells with distinct
molecular components, cannot be properly explained within such a framework. We
therefore consider additional mechanisms that can reproduce the proper
cytoplasmic polarity. In particular we predict that cytoskeletal asymmetry in
the cytoplasm, in addition to the cortical actomyosin asymmetry, is a critical
determinant of PAR protein localization.


  Analyzing the properties of duplicate genes during evolution is useful to
understand the development of new cell functions. The yeast S. cerevisiae is a
useful testing ground for this problem, because its duplicated genes with
different evolutionary birth and destiny are well distinguishable. In
particular, there is a clear detection for the occurrence of a Whole Genome
Duplication (WGD) event in S. cerevisiae, and the genes derived from this event
(WGD paralogs) are known. We studied WGD and non-WGD duplicates by two parallel
analysis based on structural protein domains and on Gene Ontology annotation
scheme respectively. The results show that while a large number of
``duplicable'' structural domains is shared in local and global duplications,
WGD and non-WGD paralogs tend to have different functions. The reason for this
is the existence of WGD and non-WGD specific domains with largely different
functions. In agreement with the recent findings of Wapinski and collaborators
(Nature 449, 2007), WGD paralogs often perform ``core'' cell functions, such as
translation and DNA replication, while local duplications associate with
``peripheral'' functions such as response to stress. Our results also support
the fact that domain architectures are a reliable tool to detect homology, as
the domains of duplicates are largely invariant with date and nature of the
duplication, while their sequences and also their functions might migrate.


  We analyze a system-level model for lytic repression of lambda-phage in E.
coli using reliability theory, showing that the repressor circuit comprises 4
redundant components whose failure mode is prophage induction. Our model
reflects the specific biochemical mechanisms involved in regulation, including
long-range cooperative binding, and its detailed predictions for prophage
induction in E. coli under ultra-violet radiation are in good agreement with
experimental data.


  Background: Transposed elements (TEs) are known to affect transcriptomes,
because either new exons are generated from intronic transposed elements (this
is called exonization), or the element inserts into the exon, leading to a new
transcript. Several examples in the literature show that isoforms generated by
an exonization are specific to a certain tissue (for example the heart muscle)
or inflict a disease. Thus, exonizations can have negative effects for the
transcriptome of an organism. Results: As we aimed at detecting other tissue-
or tumor-specific isoforms in human and mouse genomes which were generated
through exonization of a transposed element, we designed the automated analysis
pipeline SERpredict (SER = Specific Exonized Retroelement) making use of
Bayesian Statistics. With this pipeline, we found several genes in which a
transposed element formed a tissue- or tumor-specific isoform. Conclusion: Our
results show that SERpredict produces relevant results, demonstrating the
importance of transposed elements in shaping both the human and the mouse
transcriptomes. The effect of transposed elements on the human transcriptome is
several times higher than the effect on the mouse transcriptome, due to the
contribution of the primate-specific Alu elements


  Background: Gene duplication and exonization of intronic transposed elements
are two mechanisms that enhance genomic diversity. We examined whether there is
less selection against exonization of transposed elements in duplicated genes
than in single-copy genes. Results: Genome-wide analysis of exonization of
transposed elements revealed a higher rate of exonization within duplicated
genes relative to single-copy genes. The gene for TIF-IA, an RNA polymerase I
transcription initiation factor, underwent a humanoid-specific triplication,
all three copies of the gene are active transcriptionally, although only one
copy retains the ability to generate the TIF-IA protein. Prior to TIF-IA
triplication, an Alu element was inserted into the first intron. In one of the
non-protein coding copies, this Alu is exonized. We identified a single point
mutation leading to exonization in one of the gene duplicates. When this
mutation was introduced into the TIF-IA coding copy, exonization was activated
and the level of the protein-coding mRNA was reduced substantially. A very low
level of exonization was detected in normal human cells. However, this
exonization was abundant in most leukemia cell lines evaluated, although the
genomic sequence is unchanged in these cancerous cells compared to normal
cells. Conclusion: The definition of the Alu element within the TIF-IA gene as
an exon is restricted to certain types of cancers; the element is not exonized
in normal human cells. These results further our understanding of the delicate
interplay between gene duplication and alternative splicing and of the
molecular evolutionary mechanisms leading to genetic innovations. This implies
the existence of purifying selection against exonization in single copy genes,
with duplicate genes free from such constrains.


  Transposable elements may acquire unrelated gene fragments into their
sequences in a process called transduplication. Transduplication of
protein-coding genes is common in plants, but is unknown of in animals. Here,
we report that the Turmoil-1 transposable element in C. elegans has
incorporated two protein-coding sequences into its inverted terminal repeat
(ITR) sequences. The ITRs of Turmoil-1 contain a conserved RNA recognition
motif (RRM) that originated from the rsp- 2 gene and a fragment from the
protein-coding region of the cpg-3 gene. We further report that an open reading
frame specific to C. elegans may have been created as a result of a Turmoil-1
insertion. Mutations at the 5' splice site of this open reading frame may have
reactivated the transduplicated RRM motif


  Background: Transposed elements (TEs) have a substantial impact on mammalian
evolution and are involved in numerous genetic diseases. We compared the impact
of TEs on the human transcriptome and the mouse transcriptome. Results: We
compiled a dataset of all TEs in the human and mouse genomes, identifying
3,932,058 and 3,122,416 TEs, respectively. We than extracted TEs located within
human and mouse genes and, surprisingly, we found that 60% of TEs in both human
and mouse are located in intronic sequences, even though introns comprise only
24% of the human genome. All TE families in both human and mouse can exonize.
TE families that are shared between human and mouse exhibit the same percentage
of TE exonization in the two species, but the exonization level of Alu, a
primatespecific retroelement, is significantly greater than that of other TEs
within the human genome, leading to a higher level of TE exonization in human
than in mouse (1,824 exons compared with 506 exons, respectively). We detected
a primate-specific mechanism for intron gain, in which Alu insertion into an
exon creates a new intron located in the 3' untranslated region (termed
'intronization'). Finally, the insertion of TEs into the first and last exons
of a gene is more frequent in human than in mouse, leading to longer exons in
human. Conclusion: Our findings reveal many effects of TEs on these two
transcriptomes. These effects are substantially greater in human than in mouse,
which is due to the presence of Alu elements in human.


  Transposed elements (TEs) are mobile genetic sequences. During the evolution
of eukaryotes TEs were inserted into active protein-coding genes, affecting
gene structure, expression and splicing patterns, and protein sequences.
Genomic insertions of TEs also led to creation and expression of new functional
non-coding RNAs such as micro- RNAs. We have constructed the TranspoGene
database, which covers TEs located inside proteincoding genes of seven species:
human, mouse, chicken, zebrafish, fruit fly, nematode and sea squirt. TEs were
classified according to location within the gene: proximal promoter TEs,
exonized TEs (insertion within an intron that led to exon creation), exonic TEs
(insertion into an existing exon) or intronic TEs. TranspoGene contains
information regarding specific type and family of the TEs, genomic and mRNA
location, sequence, supporting transcript accession and alignment to the TE
consensus sequence. The database also contains host gene specific data: gene
name, genomic location, Swiss-Prot and RefSeq accessions, diseases associated
with the gene and splicing pattern. In addition, we created microTranspoGene: a
database of human, mouse, zebrafish and nematode TEderived microRNAs. The
TranspoGene and micro- TranspoGene databases can be used by researchers
interested in the effect of TE insertion on the eukaryotic transcriptome.


  Alternative cassette exons are known to originate from two processes
exonization of intronic sequences and exon shuffling. Herein, we suggest an
additional mechanism by which constitutively spliced exons become alternative
cassette exons during evolution. We compiled a dataset of orthologous exons
from human and mouse that are constitutively spliced in one species but
alternatively spliced in the other. Examination of these exons suggests that
the common ancestors were constitutively spliced. We show that relaxation of
the 59 splice site during evolution is one of the molecular mechanisms by which
exons shift from constitutive to alternative splicing. This shift is associated
with the fixation of exonic splicing regulatory sequences (ESRs) that are
essential for exon definition and control the inclusion level only after the
transition to alternative splicing. The effect of each ESR on splicing and the
combinatorial effects between two ESRs are conserved from fish to human. Our
results uncover an evolutionary pathway that increases transcriptome diversity
by shifting exons from constitutive to alternative splicing


  Examination of the human transcriptome reveals higher levels of RNA editing
than in any other organism tested to date. This is indicative of extensive
double-stranded RNA (dsRNA) formation within the human transcriptome. Most of
the editing sites are located in the primate-specific retrotransposed element
called Alu. A large fraction of Alus are found in intronic sequences, implying
extensive Alu-Alu dsRNA formation in mRNA precursors. Yet, the effect of these
intronic Alus on splicing of the flanking exons is largely unknown. Here, we
show that more Alus flank alternatively spliced exons than constitutively
spliced ones; this is especially notable for those exons that have changed
their mode of splicing from constitutive to alternative during human evolution.
This implies that Alu insertions may change the mode of splicing of the
flanking exons. Indeed, we demonstrate experimentally that two Alu elements
that were inserted into an intron in opposite orientation undergo base-pairing,
as evident by RNA editing, and affect the splicing patterns of a downstream
exon, shifting it from constitutive to alternative. Our results indicate the
importance of intronic Alus in influencing the splicing of flanking exons,
further emphasizing the role of Alus in shaping of the human transcriptome


  In this review, we summarize recent developments in stochastic evolutionary
game dynamics of finite populations.


  While slowly turning the ends of a single molecule of DNA at constant applied
force, a discontinuity was recently observed at the supercoiling transition,
when a small plectoneme is suddenly formed. This can be understood as an abrupt
transition into a state in which stretched and plectonemic DNA coexist. We
argue that there should be discontinuities in both the extension and the torque
at the transition, and provide experimental evidence for both. To predict the
sizes of these discontinuities and how they change with the overall length of
DNA, we organize a theory for the coexisting plectonemic state in terms of four
length-independent parameters. We also test plectoneme theories, including our
own elastic rod simulation, finding discrepancies with experiment that can be
understood in terms of the four coexisting state parameters.


  Darwin's hypothesis that all extant life forms are descendants of a last
common ancestor cell and diversification of life forms results from gradual
mutation plus natural selection represents a mainstream view that has
influenced biology and even society for over a century. However, this Darwinian
view on life is contradicted by many observations and lacks a plausible
physico-chemical explanation. Strong evidence suggests that the common ancestor
cell hypothesis is the most fundamental flaw of Darwinism. By contrast, a
totally different perspective on origin and evolution of life claims that
cellular life forms were descendants of already diversified acellular life
forms. Independently originated life forms evolve largely in some parallel ways
even though they also interact with each other. Some evolutionary "gaps"
naturally exist among evolutionary lines. Similarity may not be the only result
of phylogenetic inheritance but may be a result of a convergent mechanism of
origin and evolution. Evolution is not a random process but follows some basic
physico-chemical principles as a result of the interplay of both energy and
entropy on matter.


  Although the existence of robust inverted biomass pyramids seem paradoxical,
they have been observed in planktonic communities, and more recently, in
pristine coral reefs. Understanding the underlying mechanisms which produce
inverted biomass pyramids provides new ecological insights, and for coral
reefs, may help mitigate or restore damaged reefs. We present three classes of
predator-prey models which elucidate mechanisms that generate robust inverted
biomass pyramids. The first class of models exploits well-mixing of predators
and prey, the second class has a refuge (with explicit size) for the prey to
hide, and the third class incorporates the immigration of prey. Our models
indicate that inverted biomass pyramids occur when the prey growth rate, prey
carrying capacity, biomass conversion efficiency, the predator life span, or
the immigration rate of prey fish is sufficiently large. In the second class,
we discuss three hypotheses on how refuge size can impact the amount of prey
available to predators. By explicitly incorporating a refuge size, these can
more realistically model predator-prey interactions than refgue models with
implicit refuge size.


  Demographic oscillators are individual-based systems exhibiting temporal
cycles sustained by the stochastic dynamics of the microscopic interacting
particles. We here use the example of coupled predator-prey oscillators to show
that synchronization to a common frequency can occur between two such systems,
even if they oscillate at different frequencies in the absence of coupling. The
power spectra of the separate and the coupled systems are computed within a van
Kampen expansion in the inverse system size, and it is found that they exhibit
two peaks at separate frequencies at low coupling, but that only one peak is
present at large enough coupling strength. We further make predictions on the
time behaviour of the phases of the two oscillators, and their phase
difference, and confirm the frequency entrainment at sufficiently large
coupling. Theoretical results are verified convincingly in numerical
simulations.


  We introduce a simple "patchy particle" model to study the thermodynamics and
dynamics of self-assembly of homomeric protein complexes. Our calculations
allow us to rationalize recent results for dihedral complexes. Namely, why
evolution of such complexes naturally takes the system into a region of
interaction space where (i) the evolutionarily newer interactions are weaker,
(ii) subcomplexes involving the stronger interactions are observed to be
thermodynamically stable on destabilization of the protein-protein interactions
and (iii) the self-assembly dynamics are hierarchical with these same
subcomplexes acting as kinetic intermediates.


  The past decade has seen great advances in our understanding of the role of
noise in gene regulation and the physical limits to signaling in biological
networks. Here we introduce the spectral method for computation of the joint
probability distribution over all species in a biological network. The spectral
method exploits the natural eigenfunctions of the master equation of
birth-death processes to solve for the joint distribution of modules within the
network, which then inform each other and facilitate calculation of the entire
joint distribution. We illustrate the method on a ubiquitous case in nature:
linear regulatory cascades. The efficiency of the method makes possible
numerical optimization of the input and regulatory parameters, revealing design
properties of, e.g., the most informative cascades. We find, for threshold
regulation, that a cascade of strong regulations converts a unimodal input to a
bimodal output, that multimodal inputs are no more informative than bimodal
inputs, and that a chain of up-regulations outperforms a chain of
down-regulations. We anticipate that this numerical approach may be useful for
modeling noise in a variety of small network topologies in biology.


  Evolution of cooperation in the prisoner's dilemma game is studied where
initially all players are linked via a regular graph, having four neighbors
each. Simultaneously with the strategy evolution, players are allowed to make
new connections and thus permanently extend their neighborhoods, provided they
have been successful in passing their strategy to the opponents. We show that
this simple coevolutionary rule shifts the survival barrier of cooperators
towards high temptations to defect and results in highly heterogeneous
interaction networks with an exponential fit best characterizing their degree
distributions. In particular, there exist an optimal maximal degree for the
promotion of cooperation, warranting the best exchange of information between
influential players.


  We perform molecular dynamics simulations for a simple coarse-grained model
of crambin placed inside of a softly repulsive sphere of radius R. The
confinement makes folding at the optimal temperature slower and affects the
folding scenarios, but both effects are not dramatic. The influence of crowding
on folding are studied by placing several identical proteins within the sphere,
denaturing them, and then by monitoring refolding. If the interactions between
the proteins are dominated by the excluded volume effects, the net folding
times are essentially like for a single protein. An introduction of
inter-proteinic attractive contacts hinders folding when the strength of the
attraction exceeds about a half of the value of the strength of the single
protein contacts. The bigger the strength of the attraction, the more likely is
the occurrence of aggregation and misfolding.


  Non-coding RNAs are crucial regulators of gene expression in prokaryotes and
eukaryotes, but it remains poorly understood how they affect the dynamics of
transcriptional networks. We analyzed the temporal characteristics of the
cyanobacterial iron stress response by mathematical modeling and quantitative
experimental analyses, and focused on the role of a recently discovered small
non-coding RNA, IsrR. We found that IsrR is responsible for a pronounced delay
in the accumulation of isiA mRNA encoding the late-phase stress protein, IsiA,
and that it ensures a rapid decline in isiA levels once external stress
triggers are removed. These kinetic properties allow the system to selectively
respond to sustained (as opposed to transient) stimuli, and thus establish a
temporal threshold, which prevents energetically costly IsiA accumulation under
short-term stress conditions. Biological information is frequently encoded in
the quantitative aspects of intracellular signals (e.g., amplitude and
duration). Our simulations reveal that competitive inhibition and regulated
degradation allow intracellular regulatory networks to efficiently discriminate
between transient and sustained inputs.


  In 1994, I conjectured that insulin-degrading enzyme (IDE) acts as an
inhibitor of malignant transformation by degrading insulin and thus preventing
this major growth-stimulatory hormone from binding and thereby inactivating the
retinoblastoma tumor suppressor protein (RB). Ten years later, I discovered
that a carboxyterminal RB amino acid sequence resembles the catalytic center of
IDE. This structural homology raised the possibility that insulin degradation
is a basic mechanism for tumor suppression shared by RB and IDE. Subsequently,
a first immunohistochemical study on the differential expression of human IDE
in normal tissues, primary tumors and their corresponding lymph node metastases
further corroborated the initial conjecture on IDE being an antineoplastic
molecule. In this report, it is shown that IDE harbors ankyrin repeat-like
amino acid sequences through which it might bind and, as a result, antagonize
the pro-inflammatory factor NF-kappaB as well as cyclin-dependent kinases
(CDKs). As equally revealed here, IDE also contains 2 RXL cyclin-binding motifs
which could contribute to its presumed inhibition of CDKs. These new findings
suggest that IDE is potentially able to suppress both inflammation and
oncogenesis by several mechanisms that ultimately ensure RB function.


  [abridged] Background: The distribution of chemical species in an open system
at metastable equilibrium can be expressed as a function of environmental
variables which can include temperature, oxidation-reduction potential and
others. Calculations of metastable equilibrium for various model systems were
used to characterize chemical transformations among proteins and groups of
proteins found in different compartments of yeast cells.
  Results: With increasing oxygen fugacity, the relative metastability fields
of model proteins for major subcellular compartments go as mitochondrion,
endoplasmic reticulum, cytoplasm, nucleus. In a metastable equilibrium setting
at relatively high oxygen fugacity, proteins making up actin are predominant,
but those constituting the microtubule occur with a low chemical activity. A
reaction sequence involving the microtubule and spindle pole proteins was
predicted by combining the known intercompartmental interactions with a
hypothetical program of oxygen fugacity changes in the local environment. In
further calculations, the most-abundant proteins within compartments generally
occur in relative abundances that only weakly correspond to a metastable
equilibrium distribution. However, physiological populations of proteins that
form complexes often show an overall positive or negative correlation with the
relative abundances of proteins in metastable assemblages.
  Conclusions: This study explored the outlines of a thermodynamic description
of chemical transformations among interacting proteins in yeast cells. The
results suggest that these methods can be used to measure the degree of
departure of a natural biochemical process or population from a local minimum
in Gibbs energy.


  Bifurcations leading to complex dynamical behaviour of non-linear systems are
often encountered when the characteristics of feedback circuits in the system
are varied. In systems with many unknown or varying parameters, it is an
interesting, but difficult problem to find parameter values for which specific
bifurcations occur. In this paper, we develop a loop breaking approach to
evaluate the influence of parameter values on feedback circuit characteristics.
This approach allows a theoretical classification of feedback circuit
characteristics related to possible bifurcations in the system. Based on the
theoretical results, a numerical algorithm for bifurcation search in a possibly
high-dimensional parameter space is developed. The application of the proposed
algorithm is illustrated by searching for a Hopf bifurcation in a model of the
mitogen activated protein kinase (MAPK) cascade, which is a classical example
for biochemical signal transduction.


  The network concept is increasingly used for the description of complex
systems. Here we summarize key aspects of the evolvability and robustness of
the hierarchical network-set of macromolecules, cells, organisms, and
ecosystems. Listing the costs and benefits of cooperation as a necessary
behaviour to build this network hierarchy, we outline the major hypothesis of
the paper: the emergence of hierarchical complexity needs cooperation leading
to the ageing (i.e. gradual deterioration) of the constituent networks. A
stable environment develops cooperation leading to over-optimization, and
forming an "always-old" network, which accumulates damage, and dies in an
apoptosis-like process. A rapidly changing environment develops competition
forming a "forever-young" network, which may suffer an occasional
over-perturbation exhausting system-resources, and causing death in a
necrosis-like process. Giving a number of examples we demonstrate how
cooperation evokes the gradual accumulation of damage typical to ageing.
Finally, we show how various forms of cooperation and consequent ageing emerge
as key elements in all major steps of evolution from the formation of
protocells to the establishment of the globalized, modern human society.


  This paper addresses the theoretical conditions necessary for some subject of
study to survive forever. A probabilistic analysis leads to some prerequisite
conditions for preserving, say, electronic data indefinitely into the future.
The general analysis would also apply to a species, a civilization, or any
subject of study, as long as there is a definition of "survival" available. A
distinction emerges between two approaches to longevity: being many or being
smart. Natural selection relies on the first method, whereas a civilization,
individual, or other singular subject must rely on the latter. A computational
model of survival incorporates the idea of Kolmogorov-type complexity for both
strategies to illustrate the role of data analysis and information processing
that may be required. The survival-through-intelligence strategy has problems
when the subject can self-modify, which is illustrated with a link to Turing's
Halting Problem. The paper concludes with comments on the Fermi Paradox.


  Understanding the time evolution of fragmented animal populations and their
habitats, connected by migration, is a problem of both theoretical and
practical interest. This paper presents a method for calculating the time
evolution of the habitats' population size distribution from a general
stochastic dynamic within each habitat, using a deterministic approximation
which becomes exact for an infinite number of habitats. Fragmented populations
are usually thought to be characterized by a separation of time scale between,
on the one hand, colonization and extinction of habitats and, on the other
hand, the local population dynamics within each habitat. The analysis in this
paper suggests an alternative view: the effective population dynamic stems from
a law of large numbers, where stochastic fluctuations in population size of
single habitats are buffered through the dispersal pool so that the global
population dynamic remains approximately smooth. For illustration, the
deterministic approximation is compared to simulations of a stochastic model
with density dependent local recruitment and mortality. The article is
concluded with a discussion of the general implications of the results, and
possible extensions of the method.


  The extraction of membrane tubes by molecular motors is known to play an
important role for the transport properties of eukaryotic cells. By studying a
generic class of models for the tube extraction, we discover a rich phase
diagram. In particular we show that the density of motors along the tube can
exhibit shocks, inverse shocks and plateaux, depending on parameters which
could in principle be probed experimentally. In addition the phase diagram
exhibits interesting reentrant behavior.


  Gene expression in a cell entails random reaction events occurring over
disparate time scales. Thus, molecular noise that often results in phenotypic
and population-dynamic consequences sets a fundamental limit to biochemical
signaling. While there have been numerous studies correlating the architecture
of cellular reaction networks with noise tolerance, only a limited effort has
been made to understand the dynamic role of protein-protein interactions. Here
we have developed a fully stochastic model for the positive feedback control of
a single gene, as well as a pair of genes (toggle switch), integrating
quantitative results from previous in vivo and in vitro studies. We find that
the overall noise-level is reduced and the frequency content of the noise is
dramatically shifted to the physiologically irrelevant high-frequency regime in
the presence of protein dimerization. This is independent of the choice of
monomer or dimer as transcription factor and persists throughout the multiple
model topologies considered. For the toggle switch, we additionally find that
the presence of a protein dimer, either homodimer or heterodimer, may
significantly reduce its random switching rate. Hence, the dimer promotes the
robust function of bistable switches by preventing the uninduced (induced)
state from randomly being induced (uninduced). The specific binding between
regulatory proteins provides a buffer that may prevent the propagation of
fluctuations in genetic activity. The capacity of the buffer is a non-monotonic
function of association-dissociation rates. Since the protein oligomerization
per se does not require extra protein components to be expressed, it provides a
basis for the rapid control of intrinsic or extrinsic noise.


  In evolutionary game dynamics, reproductive success increases with the
performance in an evolutionary game. If strategy $A$ performs better than
strategy $B$, strategy $A$ will spread in the population. Under stochastic
dynamics, a single mutant will sooner or later take over the entire population
or go extinct. We analyze the mean exit times (or average fixation times)
associated with this process. We show analytically that these times depend on
the payoff matrix of the game in an amazingly simple way under weak selection,
ie strong stochasticity: The payoff difference $\Delta \pi$ is a linear
function of the number of $A$ individuals $i$, $\Delta \pi = u i + v$. The
unconditional mean exit time depends only on the constant term $v$. Given that
a single $A$ mutant takes over the population, the corresponding conditional
mean exit time depends only on the density dependent term $u$. We demonstrate
this finding for two commonly applied microscopic evolutionary processes.


  We study the evolution of cooperation in structured populations within
popular models of social dilemmas, whereby simple coevolutionary rules are
introduced that may enhance players abilities to enforce their strategy on the
opponent. Coevolution thus here refers to an evolutionary process affecting the
teaching activity of players that accompanies the evolution of their
strategies. Particularly, we increase the teaching activity of a player after
it has successfully reproduced, yet we do so depending on the disseminated
strategy. We separately consider coevolution affecting either only the
cooperators or only the defectors, and show that both options promote
cooperation irrespective of the applied game. Opposite to intuitive reasoning,
however, we reveal that the coevolutionary promotion of players spreading
defection is, in the long run, more beneficial for cooperation than the
likewise promotion of cooperators. We explain the contradictory impact of the
two considered coevolutionary rules by examining the differences between
resulting heterogeneities that segregate participating players, and
furthermore, demonstrate that the influential individuals completely determine
the final outcome of the games. Our findings are immune to changes defining the
type of considered social dilemmas and highlight that the heterogeneity of
players, resulting in a positive feedback mechanism, is a fundamental property
promoting cooperation in groups of selfish individuals.


  Simulating the evolution of the Human Immunodeficiency Virus (HIV) epidemic
requires a detailed description of the population network, especially for small
populations in which individuals can be represented in detail and accuracy. In
this paper, we introduce the concept of a Complex Agent Network(CAN) to model
the HIV epidemics by combining agent-based modelling and complex networks, in
which agents represent individuals that have sexual interactions. The
applicability of CANs is demonstrated by constructing and executing a detailed
HIV epidemic model for men who have sex with men (MSM) in Amsterdam, including
a distinction between steady and casual relationships. We focus on MSM contacts
because they play an important role in HIV epidemics and have been tracked in
Amsterdam for a long time. Our experiments show good correspondence between the
historical data of the Amsterdam cohort and the simulation results.


  We present a dynamical model of DNA mechanical unzipping under the action of
a force. The model includes the motion of the fork in the sequence-dependent
landscape, the trap(s) acting on the bead(s), and the polymeric components of
the molecular construction (unzipped single strands of DNA, and linkers).
Different setups are considered to test the model, and the outcome of the
simulations is compared to simpler dynamical models existing in the literature
where polymers are assumed to be at equilibrium.


  We have developed a new simulation method to estimate the distance between
the native state and the first transition state, and the distance between the
intermediate state and the second transition state of a protein which
mechanically unfolds via intermediates. Assuming that the end-to-end extension
$\Delta R$ is a good reaction coordinate to describe the free energy landscape
of proteins subjected to an external force, we define the midpoint extension
$\Delta R^*$ between two transition states from either constant-force or
constant loading rate pulling simulations. In the former case, $\Delta R^*$ is
defined as a middle point between two plateaus in the time-dependent curve of
$\Delta R$, while, in the latter one, it is a middle point between two peaks in
the force-extension curve. Having determined $\Delta R^*$, one can compute
times needed to cross two transition state barriers starting from the native
state. With the help of the Bell and microscopic kinetic theory, force
dependencies of these unfolding times can be used to locate the intermediate
state and to extract unfolding barriers. We have applied our method to the
titin domain I27 and the fourth domain of {\em Dictyostelium discoideum}
filamin (DDFLN4), and obtained reasonable agreement with experiments, using the
C$_{\alpha}$-Go model.


  Using exhaustive Monte Carlo simulations we study the kinetics and mechanism
of fibril formation using lattice models as a function of temperature and the
number of chains. While these models are, at best, caricatures of peptides, we
show that a number of generic features thought to govern fibril assembly are
present in the toy model. The monomer, which contains eight beads made from
three letters (hydrophobic, polar, and charged), adopts a compact conformation
in the native state. The kinetics of fibril assembly occurs in three distinct
stages. In each stage there is a cascade of events that transforms the monomers
and oligomers to ordered structures. In the first "burst" stage highly mobile
oligomers of varying sizes form. The conversion to the aggregation-prone
conformation occurs within the oligomers during the second stage. As time
progresses, a dominant cluster emerges that contains a majority of the chains.
In the final stage, the aggregation-prone conformation particles serve as a
template onto which smaller oligomers or monomers can dock and undergo
conversion to fibril structures. The overall time for growth in the latter
stages is well described by the Lifshitz-Slyazov growth kinetics for
crystallization from super-saturated solutions.


  We present properties of Lotka-Volterra equations describing ecological
competition among a large number of competing species. First we extend to the
case of a non-homogeneous niche space stability conditions for solutions
representing species coexistence. Second, we discuss mechanisms leading to
species clustering and obtain an analytical solution for a lumped state in a
specific instance of the system. We also discuss how realistic ecological
interactions may result in different types of competition coefficients.


  If we follow an asexually reproducing population through time, then the
amount of time that has passed since the most recent common ancestor (MRCA) of
all current individuals lived will change as time progresses. The resulting
"MRCA age" process has been studied previously when the population has a
constant large size and evolves via the diffusion limit of standard
Wright--Fisher dynamics. For any population model, the sample paths of the MRCA
age process are made up of periods of linear upward drift with slope +1
punctuated by downward jumps. We build other Markov processes that have such
paths from Poisson point processes on $\mathbb{R}_{++}\times\mathbb{R}_{++}$
with intensity measures of the form $\lambda\otimes\mu$ where $\lambda$ is
Lebesgue measure, and $\mu$ (the "family lifetime measure") is an arbitrary,
absolutely continuous measure satisfying $\mu((0,\infty))=\infty$ and
$\mu((x,\infty))<\infty$ for all $x>0$. Special cases of this construction
describe the time evolution of the MRCA age in $(1+\beta)$-stable continuous
state branching processes conditioned on nonextinction--a particular case of
which, $\beta=1$, is Feller's continuous state branching process conditioned on
nonextinction. As well as the continuous time process, we also consider the
discrete time Markov chain that records the value of the continuous process
just before and after its successive jumps. We find transition probabilities
for both the continuous and discrete time processes, determine when these
processes are transient and recurrent and compute stationary distributions when
they exist.


  We develop an approach aimed at optimizing the parameters of a network of
biochemical logic gates for reduction of the "analog" noise buildup.
Experiments for three coupled enzymatic AND gates are reported, illustrating
our procedure. Specifically, starch - one of the controlled network inputs - is
converted to maltose by beta-amylase. With the use of phosphate (another
controlled input), maltose phosphorylase then produces glucose. Finally,
nicotinamide adenine dinucleotide (NAD+) - the third controlled input - is
reduced under the action of glucose dehydrogenase to yield the optically
detected signal. Network functioning is analyzed by varying selective inputs
and fitting standardized few-parameters "response-surface" functions assumed
for each gate. This allows a certain probe of the individual gate quality, but
primarily yields information on the relative contribution of the gates to noise
amplification. The derived information is then used to modify our experimental
system to put it in a regime of a less noisy operation.


  The extent of exciton coherence in protein-pigment complexes has significant
implications for the initial light harvesting step in photosynthetic organisms.
In this work we model the main antenna protein of photosystem II, namely light
harvesting complex II (LHC-II), with a single-exciton Hamiltonian with sites
coupled via dipole-dipole interaction, with linear coupling to a dissipative
phonon bath. With appropriate parameters, Monte Carlo path integral (MCPI)
results of the exciton coherence length from 1 K to 500 K show that at
thermodynamic equilibrium, an exciton in LHC-II is localized mostly on 2 single
chlorophyll pigment sites, with persistent short-range coherence over the A2-B2
pair, A3-B3 pair and B1-B5-B6 triplet. Quasi-adiabatic path integral (QUAPI)
calculations of the subsystems mentioned above show a smooth, incoherent
relaxation towards thermodynamic equilibrium. The results obtained imply that
with the exception of small coherent subsystems at cryogenic temperatures,
excitons in LHC-II are more localized than in the analogous light harvesting
complex II (LH-II) of the purple bacterium Rs. molischianum, which may be
expected from the lower symmetry of the former.


  Many types of tumors exhibit chromosomal losses or gains, as well as local
amplifications and deletions. Within any given tumor type, sample specific
amplifications and deletionsare also observed. Typically, a region that is
aberrant in more tumors,or whose copy number change is stronger, would be
considered as a more promising candidate to be biologically relevant to cancer.
We sought for an intuitive method to define such aberrations and prioritize
them. We define V, the volume associated with an aberration, as the product of
three factors: a. fraction of patients with the aberration, b. the aberrations
length and c. its amplitude. Our algorithm compares the values of V derived
from real data to a null distribution obtained by permutations, and yields the
statistical significance, p value, of the measured value of V. We detected
genetic locations that were significantly aberrant and combined them with
chromosomal arm status to create a succint fingerprint of the tumor genome.
This genomic fingerprint is used to visualize the tumors, highlighting events
that are co ocurring or mutually exclusive. We allpy the method on three
different public array CGH datasets of Medulloblastoma and Neuroblastoma, and
demonstrate its ability to detect chromosomal regions that were known to be
altered in the tested cancer types, as well as to suggest new genomic locations
to be tested. We identified a potential new subtype of Medulloblastoma, which
is analogous to Neuroblastoma type 1.


  We propose a simple model of cytochrome c oxidase, including four redox
centers and four protonable sites, to study the time evolution of
electrostatically coupled electron and proton transfers initiated by the
injection of a single electron into the enzyme. We derive a system of master
equations for electron and proton state probabilities and show that an
efficient pumping of protons across the membrane can be obtained for a
reasonable set of parameters. All four experimentally observed kinetic phases
appear naturally from our model. We also calculate the dependence of the
pumping efficiency on the transmembrane voltage at different temperatures and
discuss a possible mechanism of the redox-driven proton translocation.


  The standard pair approximation equations (PA) for the
Susceptible-Infective-Recovered-Susceptible (SIRS) model of infection spread on
a network of homogeneous degree $k$ predict a thin phase of sustained
oscillations for parameter values that correspond to diseases that confer long
lasting immunity. Here we present a study of the dependence of this oscillatory
phase on the parameter $k$ and of its relevance to understand the behaviour of
simulations on networks. For $k=4$, we compare the phase diagram of the PA
model with the results of simulations on regular random graphs (RRG) of the
same degree. We show that for parameter values in the oscillatory phase, and
even for large system sizes, the simulations either die out or exhibit damped
oscillations, depending on the initial conditions. This failure of the standard
PA model to capture the qualitative behaviour of the simulations on large RRGs
is currently being investigated.


  We propose, in this paper, both a presentation and a discussion of Paul
Alsberg's thesis on the supposed specificity principle of human evolution. The
author maintains a difference of nature between Man and Animal relying on an
opposition between "body-adaptation" - that of the Animal - and
"extrabodily-adaptation" - that of Man in which the means of adaptation are
switched outside of the organisms by tool-using. This difference is not a mere
difference of state, but of evolutionary dynamics. Here, Man is not simply
"Homo faber", as in Bergson's view, but produced and made possible by
technique; a technique which then appears as an hominisation factor. Thus, his
"principle of body-liberation" by tool-using is to be retrospectively
understood as a part of the logics of the modification of selection pressure
logics, which reminds us the seminal contemporary niche construction theory (F.
John Odling-Smee). It seems therefore possible to make Paul Alsberg, from his
1922 work, one of the most important precursors of the cultural niche
construction theory.


  The present paper discusses the use of modified Lotka-Volterra equations in
order to stochastically simulate the behaviour of Listeria monocytogenes and
Lactic Acid Bacteria (LAB) during the fermentation period (168 h) of a typical
Sicilian salami. For this purpose, the differential equation system is set
considering T, pH and aw as stochastic variables. Each of them is governed by
dynamics that involve a deterministic linear decrease as a function of the time
t and an "additive noise" term which instantaneously mimics the fluctuations of
T, pH and aw. The choice of a suitable parameter accounting for the interaction
of LAB on L. monocytogenes as well as the introduction of appropriate noise
levels allows to match the observed data, both for the mean growth curves and
for the probability distribution of L. monocytogenes concentration at 168 h.


  Physiological changes which result in changes in bacterial gene expression
are often accompanied by changes in the growth rate for fast adapting enteric
bacteria. Since the availability of RNA polymerase (RNAP) in cells is dependent
on the growth rate, transcriptional control involves not only the regulation of
promoters, but also depends on the available (or free) RNAP concentration which
is difficult to quantify directly. Here we develop a simple physical model
describing the partitioning of cellular RNAP into different classes: RNAPs
transcribing mRNA and ribosomal RNA (rRNA), RNAPs non-specifically bound to
DNA, free RNAP, and immature RNAP. Available experimental data for E. coli
allow us to determine the two unknown parameters of the model and hence deduce
the free RNAP concentration at different growth rates. The results allow us to
predict the growth-rate dependence of the activities of constitutive
(unregulated) promoters, and to disentangle the growth-rate dependent
regulation of promoters (e.g., the promoters of rRNA operons) from changes in
transcription due to changes in the free RNAP concentration at different growth
rates. Our model can quantitatively account for the observed changes in gene
expression patterns in mutant E. coli strains with altered levels of RNAP
expression without invoking additional parameters. Applying our model to the
case of the stringent response following amino acid starvation, we can evaluate
the plausibility of various scenarios of passive transcriptional control
proposed to account for the observed changes in the expression of rRNA and
biosynthetic operons.