  Ni\'eb\'e is a food leguminous plant cultivated in tropical Africa for its
seeds rich in proteins. The main problem setted by its production is the
conservation of harvests. In the fields as in the stocks, the seeds are
destroyed by pests (bruchids). These bruchids are always associated with
several entomophagous species of hymenoptera. Four entomophagous species were
listed : an egg parasitoid (U lariophaga Stephan), and three solitary larval
and pupal ectoparasitoids (D. Basalis Rondoni, Pteromalidae; E. vuilleti
Crawford and E. orientalis Crawford, Eupelmidae). The survey of the populations
shows that at the beginning of storage, E orientalis is the most abundant
specie (72 %) whereas E. vuilleti and D. Basalis respectively represent 12 %
and 16 % of the hymenoptera. During storage, the E orientalis population
decreases gradually and it disappears completely in less than two months after
the beginning of storage. E. Vuilleti population becomes gradually more
important than D. basalis population which regress until less than 10 % of the
emerging parasitoids. E vuilleti adopts ovicide and larvicide behaviour against
D. Basalis. This behaviour explains its population regression inside granaries.
If the aggressive behaviour of this Eupelmidae is a constant, that could also
explain the disappearance of E orientalis. However if this species is
maintained in stocks, it would be an effective control agent of bruchids
according to their parasitic capacities. This study shows that ovicide and
larvicide behaviour of E vuilleti is not expressed against E orientalis. When
the females have exclusively the hosts already parasitized by E orientalis,
they do not lay eggs. The disappearance of E orientalis could not thus be
explained by the presence of E. vuilleti.


  Bacterial chemotaxis in Escherichia coli is a canonical system for the study
of signal transduction. A remarkable feature of this system is the coexistence
of precise adaptation in population with large fluctuating cellular behavior in
single cells (Korobkova et al. 2004, Nature, 428, 574). Using a stochastic
model, we found that the large behavioral variability experimentally observed
in non-stimulated cells is a direct consequence of the architecture of this
adaptive system. Reversible covalent modification cycles, in which methylation
and demethylation reactions antagonistically regulate the activity of
receptor-kinase complexes, operate outside the region of first-order kinetics.
As a result, the receptor-kinase that governs cellular behavior exhibits a
sigmoidal activation curve. This curve simultaneously amplifies the inherent
stochastic fluctuations in the system and lengthens the relaxation time in
response to stimulus. Because stochastic fluctuations cause large behavioral
variability and the relaxation time governs the average duration of runs in
response to small stimuli, cells with the greatest fluctuating behavior also
display the largest chemotactic response. Finally, Large-scale simulations of
digital bacteria suggest that the chemotaxis network is tuned to simultaneously
optimize the random spread of cells in absence of nutrients and the cellular
response to gradients of attractant.


  Biological networks have evolved to be highly functional within uncertain
environments while remaining extremely adaptable. One of the main contributors
to the robustness and evolvability of biological networks is believed to be
their modularity of function, with modules defined as sets of genes that are
strongly interconnected but whose function is separable from those of other
modules. Here, we investigate the in silico evolution of modularity and
robustness in complex artificial metabolic networks that encode an increasing
amount of information about their environment while acquiring ubiquitous
features of biological, social, and engineering networks, such as scale-free
edge distribution, small-world property, and fault-tolerance. These networks
evolve in environments that differ in their predictability, and allow us to
study modularity from topological, information-theoretic, and gene-epistatic
points of view using new tools that do not depend on any preconceived notion of
modularity. We find that for our evolved complex networks as well as for the
yeast protein-protein interaction network, synthetic lethal pairs consist
mostly of redundant genes that lie close to each other and therefore within
modules, while knockdown suppressor pairs are farther apart and often straddle
modules, suggesting that knockdown rescue is mediated by alternative pathways
or modules. The combination of network modularity tools together with genetic
interaction data constitutes a powerful approach to study and dissect the role
of modularity in the evolution and function of biological networks.


  Structural fluctuations in the thermal equilibrium of the kinesin motor
domain are studied using a lattice protein model with Go interactions. By means
of the multi-self-overlap ensemble (MSOE) Monte Carlo method and the principal
component analysis (PCA), the free-energy landscape is obtained. It is shown
that kinesins have two subdomains that exhibit partial folding/unfolding at
functionally important regions: one is located around the nucleotide binding
site and the other includes the main microtubule binding site. These subdomains
are consistent with structural variability that was reported recently based on
experimentally-obtained structures. On the other hand, such large structural
fluctuations have not been captured by B-factor or normal mode analyses. Thus,
they are beyond the elastic regime, and it is essential to take into account
chain connectivity for studying the function of kinesins.


  It is common to use reciprocal best hits, also known as a boomerang
criterion, for determining orthology between sequences. The best hits may be
found by blast, or by other more recently developed algorithms. Previous work
seems to have assumed that reciprocal best hits is a sufficient but not
necessary condition for orthology. In this article, I explain why reciprocal
best hits cannot logically be a sufficient condition for orthology. If
reciprocal best hits is neither sufficient nor necessary for orthology, it
would seem worthwhile to examine further the logical foundations of some
unsupervised algorithms that are used to identify orthologs.


  The persistent pods of the tree, Acacia caven, that do not fall from the tree
provide opportunities for the appearance of a diverse group of insects the
following season. Such pods collected during the spring of 1999 in Chile were
indehiscent with highly sclerified pod walls. In contrast, persistent pods
collected in Uruguay after a wet winter and spring (2002) were partially
dehiscent, inducing the deterioration of the woody pods, and consequently
exposing the seeds. These persistent pods are a natural refuge for insect
species, namely two bruchid beetles (Pseudopachymeria spinipes, Stator
furcatus), one scolytidae (Dendroctonus sp), lepidopterous larvae, ant colonies
(Camponotus sp),one species of oophagous parasitoid (Uscana espinae group
senex), the gregarious larval-pupae parasitoid Monoksa dorsiplana
(Pteromalidae) and two species of Horismenus spp. (Eulophidae). The patriline
of M. dorsiplana is frequently formed by 1 son +7 daughters.


  A coarse-grained variational model is used to investigate the polymer
dynamics of barrier crossing for a diverse set of two-state folding proteins.
The model gives reliable folding rate predictions provided excluded volume
terms that induce minor structural cooperativity are included in the
interaction potential. In general, the cooperative folding routes have sharper
interfaces between folded and unfolded regions of the folding nucleus and
higher free energy barriers. The calculated free energy barriers are strongly
correlated with native topology as characterized by contact order. Increasing
the rigidity of the folding nucleus changes the local structure of the
transition state ensemble non-uniformly across the set of protein studied.
Neverthless, the calculated prefactors k0 are found to be relatively uniform
across the protein set, with variation in 1/k0 less than a factor of five. This
direct calculation justifies the common assumption that the prefactor is
roughly the same for all small two-state folding proteins. Using the barrier
heights obtained from the model and the best fit monomer relaxation time 30ns,
we find that 1/k0 (1-5)us (with average 1/k0 4us). This model can be extended
to study subtle aspects of folding such as the variation of the folding rate
with stability or solvent viscosity, and the onset of downhill folding.


  Pathways are integral to systems biology. Their classical representation has
proven useful but is inconsistent in the meaning assigned to each arrow (or
edge) and inadvertently implies the isolation of one pathway from another.
Conversely, modern high-throughput experiments give rise to standardized
networks facilitating topological calculations. Combining these perspectives,
we can embed classical pathways within large-scale networks and thus
demonstrate the crosstalk between them. As more diverse types of
high-throughput data become available, we can effectively merge both
perspectives, embedding pathways simultaneously in multiple networks. However,
the original problem still remains - the current edge representation is
inadequate to accurately convey all the information in pathways. Therefore, we
suggest that a standardized, well-defined, edge ontology is necessary and
propose a prototype here, as a starting point for reaching this goal.


  In this paper, we consider the statistical analysis of a protein interaction
network. We propose a Bayesian model that uses a hierarchy of probabilistic
assumptions about the way proteins interact with one another in order to: (i)
identify the number of non-observable functional modules; (ii) estimate the
degree of membership of proteins to modules; and (iii) estimate typical
interaction patterns among the functional modules themselves. Our model
describes large amount of (relational) data using a relatively small set of
parameters that we can reliably estimate with an efficient inference algorithm.
We apply our methodology to data on protein-to-protein interactions in
saccharomyces cerevisiae to reveal proteins' diverse functional roles. The case
study provides the basis for an overview of which scientific questions can be
addressed using our methods, and for a discussion of technical issues.


  We consider an asexual population evolving on rugged fitness landscapes which
are defined on the multi-dimensional genotypic space and have many local
optima. We track the most populated genotype as it changes when the population
jumps from a fitness peak to a better one during the process of adaptation.
This is done using the dynamics of the shell model which is a simplified
version of the quasispecies model for infinite populations and standard
Wright-Fisher dynamics for large finite populations. We show that the
population fraction of a genotype obtained within the quasispecies model and
the shell model match for fit genotypes and at short times, but the dynamics of
the two models are identical for questions related to the most populated
genotype. We calculate exactly several properties of the jumps in infinite
populations some of which were obtained numerically in previous works. We also
present our preliminary simulation results for finite populations. In
particular, we measure the jump distribution in time and find that it decays as
$t^{-2}$ as in the quasispecies problem.


  A dynamical system model for tumor -- immune system interaction together with
a method to mimic radiation therapy are proposed. A large population of virtual
patients is simulated following an ideal radiation treatment. A characteristic
parameter, the Immune System -- Tumor Efficiency Rate (ISTER), is introduced.
ISTER dependence of treatment success and other features is studied.
Statistical results allow us to give a patient classification scheme.
Radiotherapy treatment biological effective dose (BED) is thus optimized based
on the patient physical condition, following the ALARA (As Low As Reasonably
Achievable) criterion.


  We present a model for diffusion in a molecularly crowded environment. The
model consists of random barriers in percolation network. Random walks in the
presence of slowly moving barriers show normal diffusion for long times, but
anomalous diffusion at intermediate times. The effective exponents for square
distance versus time usually are below one at these intermediate times, but can
be also larger than one for high barrier concentrations. Thus we observe sub-
as well as super-diffusion in a crowded environment.


  We study the problem of the emergence of cooperation in the spatial
Prisoner's Dilemma. The pioneering work by Nowak and May showed that large
initial populations of cooperators can survive and sustain cooperation in a
square lattice with imitate-the-best evolutionary dynamics. We revisit this
problem in a cost-benefit formulation suitable for a number of biological
applications. We show that if a fixed-amount reward is established for
cooperators to share, a single cooperator can invade a population of defectors
and form structures that are resilient to re-invasion even if the reward
mechanism is turned off. We discuss analytically the case of the invasion by a
single cooperator and present agent-based simulations for small initial
fractions of cooperators. Large cooperation levels, in the sustainability
range, are found. In the conclusions we discuss possible applications of this
model as well as its connections with other mechanisms proposed to promote the
emergence of cooperation.


  We discuss various models of ion transport through cell membrane channels.
Recent experimental data shows that sizes of ion channels are compared to those
of ions and that only few ions may be simultaneously in any single channel.
Theoretical description of ion transport in such channels should therefore take
into account interactions between ions and between ions and channel proteins.
This is not satisfied by macroscopic continuum models based on
Poisson-Nernst-Planck equations. More realistic descriptions of ion transport
are offered by microscopic Brownian and molecular dynamics. One should also
take into account a dynamical character of the channel structure. This is not
yet addressed in the literature


  In this paper we present a novel framework for sequence to shape maps. These
combinatorial maps realize exponentially many shapes, and have preimages which
contain extended connected subgraphs of diameter n (neutral networks). We prove
that all basic properties of RNA folding maps also hold for combinatorial maps.
Our construction is as follows: suppose we are given a graph $H$ over the $\{1
>...,n\}$ and an alphabet of nucleotides together with a symmetric relation
$\mathcal{R}$, implied by base pairing rules. Then the shape of a sequence of
length n is the maximal H subgraph in which all pairs of nucleotides incident
to H-edges satisfy $\mathcal{R}$. Our main result is to prove the existence of
at least $\sqrt{2}^{n-1}$ shapes with extended neutral networks, i.e. shapes
that have a preimage with diameter $n$ and a connected component of size at
least $(\frac{1+\sqrt{5}}{2})^n+(\frac{1-\sqrt{5}}{2})^n$. Furthermore, we show
that there exists a certain subset of shapes which carries a natural graph
structure. In this graph any two shapes are connected by a path of shapes with
respective neutral networks of distance one. We finally discuss our results and
provide a comparison with RNA folding maps.


  We introduce a topology-based nonlinear network model of protein dynamics
with the aim of investigating the interplay of spatial disorder and
nonlinearity. We show that spontaneous localization of energy occurs
generically and is a site-dependent process. Localized modes of nonlinear
origin form spontaneously in the stiffest parts of the structure and display
site-dependent activation energies. Our results provide a straightforward way
for understanding the recently discovered link between protein local stiffness
and enzymatic activity. They strongly suggest that nonlinear phenomena may play
an important role in enzyme function, allowing for energy storage during the
catalytic process.


  We have investigated the migration behaviors of uranium (U) and thorium (Th)
in Earth and other terrestrial planets. Theoretical models of U and Th
migration have been proposed. These models suggest that the unique features of
Earth are closely connected with its unique U and Th migration models and
distribution patterns. In the Earth, U and Th can combine with oxidative
volatile components and water, migrate up to the asthenosphere position to form
an enrichment zone (EZ) of U and Th first, and then migrate up further to the
crusts through magmatism and metamorphism. We emphasize that the formation of
an EZ of U, Th and other heat-producing elements is a prerequisite for the
formation of a plate tectonic system. The heat-producing elements, currently
mainly U and Th, in the EZ are also the energy sources that drive the formation
and evolution of the crust of Earth and create special granitic continental
crusts. In other terrestrial planets, including Mercury, Venus, and Mars, an EZ
can not be formed because of a lack of oxidative volatile components and water.
For this reason, a plate tectonic system can not been developed in these
planets. We also emphasize the influence of U and Th in EZ on the development
and evolution of life on Earth. We propose that since the Earth and planets
were born in a united solar system, there should be some common mechanisms to
create the similarities and differences between them. We have tried to develop
an integrated view to explain some problems in the tectonics of Earth and
evolution, bio-evolution, and planetary dynamics through U and Th geochemistry.
We believe that a comprehensive exploration on energy sources and their
evolution is a good way to build bridges between different disciplines of
science in order to better understand the Earth and planets.


  We incorporate hydrodynamic interactions in a structure-based model of
ubiquitin and demonstrate that the hydrodynamic coupling may reduce the peak
force when stretching the protein at constant speed, especially at larger
speeds. Hydrodynamic interactions are also shown to facilitate unfolding at
constant force and inhibit stretching by fluid flows.


  Lieberman and Melott built their recent arXiv preprint 0704.2896 on my
published paper and (a preprint of) a subsequent comment by Liebermans
associate Cornette. But had this group waited for the Cornette comment to
actually appear in print together with the expected Reply, they would have
learned that his comment exposes Cornettes confusion that likely was due to
journal misprint of my figure. Thus 0704.2896 is baseless. Despite receiving
the extended Reply with Errata, these authors still fail to recognize that
detrending of paleontological records-which they erroneously promote as a
must-is an arbitrary rather than a universal operation.


  The influence of three well-known disaccharides, namely trehalose, maltose
and sucrose, on some structural and dynamical properties of lysozyme has been
investigated by means of molecular dynamics computer simulations in the 37-60
wt % concentration range. The effects of sugars on the protein conformation are
found relatively weak, in agreement with the preferential hydration of
lysozyme. Conversely, sugars seem to increase significantly the relaxation
times of the protein. These effects are shown to be correlated to the
fractional solvent accessibilities of lysozyme residues and further support the
slaving of protein dynamics. Moreover, a significant increase in the relaxation
times of lysozyme, sugars and water molecules is observed within the studied
concentration range and may result from the percolation of the hydrogen-bond
network of sugar molecules. This percolation appears to be of primary
importance to explain the influence of sugars on the dynamical properties of
lysozyme and water.


  Initial reaction rate data for lactic dehydrogenase / pyruvate, lactic
dehydrogenase / lactate and malic dehydrogenase / malate enzyme reactions were
analyzed to obtain activation free energy changes of -329, -195 and -221
cal/mole, respectively, for rate increases associated with time-specific
irradiation of the crystalline substrates prior to dissolution and
incorporation in the reaction solutions. These energies, presumably, correspond
to conformational or vibrational changes in the reactants or the activated
complex. For the lactic dehydrogenase / pyruvate reaction, it is estimated that
on the order of 10% of the irradiation energy (546 nm, 400 footcandles for 5
seconds) would be required to produce the observed reaction rate increase if a
presumed photoproduct is consumed stoichiometrically with the pyruvate
substrate. These findings are consistent with the proposition that the observed
reaction rate enhancement involves photoproducts derived from oscillatory
atmospheric gas reactions at the crystalline enzyme substrate surfaces rather
than photo-excitations of the substrate molecules, per se.


  Persistent activity is postulated to drive neural network plasticity and
learning. To investigate its underlying cellular mechanisms, we developed a
biophysically tractable model that explains the emergence, sustenance, and
eventual termination of short-term persistent activity. Using the model, we
reproduced the features of reverberating activity that were observed in small
(50-100 cells) networks of cultured hippocampal neurons, such as the appearance
of polysynaptic current clusters, the typical inter-cluster intervals, the
typical duration of reverberation, and the response to changes in
extra-cellular ionic composition. The model relies on action
potential-triggered residual presynaptic calcium, which we suggest plays an
important role in sustaining reverberations. We show that reverberatory
activity is maintained by enhanced asynchronous transmitter release from
pre-synaptic terminals, which in itself depends on the dynamics of residual
presynaptic calcium. Hence, asynchronous release, rather than being a "synaptic
noise", can play an important role in network dynamics. Additionally, we found
that a fast timescale synaptic depression is responsible for oscillatory
network activation during reverberations, whereas the onset of a slow timescale
depression leads to the termination of reverberation. The simplicity of our
model enabled a number of predictions that were confirmed by additional
analyses of experimental manipulations.


  In 1968, Sorin Comorosan first reported a phenomenon wherein irradiation of
the substrate of an enzyme reaction, in the crystalline state, for a specific
number of seconds could lead to an enhanced aqueous solution reaction rate for
the enzyme(up to 30%). Dependence on crystal irradiation time was found to be
oscillatory with a fixed period. The basis for this unusual phenomenon has
remained a mystery. Previously unreported experimental results are presented
which demonstrate, for the LDH / pyruvate reaction, that the identity of the
crystalline material irradiated is, largely, inconsequential. It is proposed
here that the irradiation procedure drives oscillatory reactions involving
atmospheric gases adsorbed on the crystals and that these photoproducts, or
related dark-reaction species, when dissolved, function as enzyme cofactors.


  We demonstrate a new algorithm for finding protein conformations that
minimize a non-bonded energy function. The new algorithm, called the difference
map, seeks to find an atomic configuration that is simultaneously in two
constraint spaces. The first constraint space is the space of atomic
configurations that have a valid peptide geometry, while the second is the
space of configurations that have a non-bonded energy below a given target.
These two constraint spaces are used to define a deterministic dynamical
system, whose fixed points produce atomic configurations in the intersection of
the two constraint spaces. The rate at which the difference map produces low
energy protein conformations is compared with that of a contemporary search
algorithm, parallel tempering. The results indicate the difference map finds
low energy protein conformations at a significantly higher rate then parallel
tempering.


  The Epstein-Barr virus (EBV) infects more than 90% of the human population,
and is the cause of several both serious and mild diseases. It is a
tumorivirus, and has been widely studied as a model system for gene
(de)regulation in human. A central feature of the EBV life cycle is its ability
to persist in human B cells in states denoted latency I, II and III. In latency
III the host cell is driven to cell proliferation and hence expansion of the
viral population, but does not enter the lytic pathway, and no new virions are
produced, while the latency I state is almost completely dormant. In this paper
we study a physico-chemical model of the switch between latency I and latency
III in EBV. We show that the unusually large number of binding sites of two
competing transcription factors, one viral and one from the host, serves to
make the switch sharper (higher Hill coefficient), either by cooperative
binding between molecules of the same species when they bind, or by competition
between the two species if there is sufficient steric hindrance.


  Vibrational energy transfer of the amide I mode of N-methylacetamide (NMA) is
studied theoretically using the vibrational configuration interaction method. A
quartic force field of NMA is constructed at the B3LYP/6-31G+(d) level of
theory and its accuarcy is checked by comparing the resulting anharmonic
frequencies with available theoretical and experimental values. Quantum
dynamics calculations for the amide I mode excitation clarify the dominant
energy transfer pathways, which sensitively depend on the anharmonic couplings
among vibrational modes. A ratio of the anharmonic coupling to the frequency
mismatch is employed to predict and interpret the dominant energy flow
pathways.


  Recent studies have reported that T cells can integrate signals between
interrupted encounters with Antigen Presenting Cells (APCs) in such a way that
the process of signal integration exhibits a form of memory. Here, we carry out
a computational study using a simple mathematical model of T cell activation to
investigate the ramifications of interrupted T cell-APC contacts on signal
integration. We consider several mechanisms of how signal integration at these
time scales may be achieved and conclude that feedback control of immediate
early gene products (IEGs) appears to be a highly plausible mechanism that
allows for effective signal integration and cytokine production from multiple
exposures to APCs. Analysis of these computer simulations provides an
experimental roadmap involving several testable predictions.


  Fifteen years ago, a structural analysis of the hormone insulin and the
retinoblastoma tumor suppressor protein (RB) revealed that they may physically
interact with one another. Subsequently, an RB peptide corresponding to the
proposed RB binding site for insulin was found to recognize full-length insulin
in vitro. As part of efforts aimed at developing this RB peptide into an
anti-cancer drug, this molecule was chemically coupled to a cellular
internalization signal and termed "MCR peptide". Meanwhile, several such MCR
peptide variants have been demonstrated to restrain the proliferation of
different human cancer cells in vitro and in vivo. Moreover, one of the MCR
peptides coined MCR-10 was shown to be capable of interfering with the complex
formation between insulin and RB in HepG2 human hepatoma cells, as monitored by
immunofluorescence. This latter result indicating an in vivo association
between insulin and RB was confirmed by a follow-up study combining the methods
of co-immunoprecipitation and immunoblotting. Here, we provide evidence for the
existence of the insulin-RB complex in A549 human non-small cell lung cancer
cells. Specifically, we demonstrate this heterodimer by means of a magnetic
beads-based immunoprecipitation approach and equally show that this dimer can
be disrupted by MCR-4 or MCR-10 each of which is known to possess
antiproliferative properties, yet to a much lesser extent by a control peptide.
Thus, this investigation has yielded another important proof for the occurrence
of the insulin-RB dimer and, furthermore, its validity as a target for
antineoplastic MCR peptides.


  As I compress on the canvas of a few pages here major results of my research
on the retinoblastoma tumor suppressor protein (RB) spreading over the past 15
years, an exciting picture emerges on this unique host molecule which surpasses
in its complexity even that of the most capable viral proteins known to date.
Accordingly, RB has the potential to bind not only growth-promoting proteins
such as insulin, but also to attach itself to calcium and oxygen, as well as to
be secreted into the extracellular environment. Moreover, RB may exert
proteolytic, antimicrobial and anti-aging activities. These condensed
structure-based insights on RB are the substance of a scientific revolution I
have initiated a long time ago, yet likely to gain even further speed in the
years to come, thus expanding both our understanding of life at the molecular
level and the possibilities for pharmacological modulation of fundamental
biological phenomena, particularly in oncology and gerontology.


  We investigate the ability of algorithms developed for reverse engineering of
transcriptional regulatory networks to reconstruct metabolic networks from
high-throughput metabolite profiling data. For this, we generate synthetic
metabolic profiles for benchmarking purposes based on a well-established model
for red blood cell metabolism. A variety of data sets is generated, accounting
for different properties of real metabolic networks, such as experimental
noise, metabolite correlations, and temporal dynamics. These data sets are made
available online. We apply ARACNE, a mainstream transcriptional networks
reverse engineering algorithm, to these data sets and observe performance
comparable to that obtained in the transcriptional domain, for which the
algorithm was originally designed.


  The model of smart migration for two-species community is developed, where
the individuals implement reflexive strategy of spatial redistribution.
Simulations have been used to figure out the situations where reflexy gives an
advantage over a non-reflexive spatial behaviour, and vice versa.


  Codon usage bias measure is defined through the mutual entropy calculation of
real codon frequency distribution against the quasi-equilibrium one. This
latter is defined in three manners: (1) the frequency of synonymous codons is
supposed to be equal (i.e., the arithmetic mean of their frequencies); (2) it
coincides to the frequency distribution of triplets; and, finally, (3) the
quasi-equilibrium frequency distribution is defined as the expected frequency
of codons derived from the dinucleotide frequency distribution. The measure of
bias in codon usage is calculated for 125 bacterial genomes.


  RNA polymerase (RNAP) is an enzyme that synthesizes a messenger RNA (mRNA)
strand which is complementary to a single-stranded DNA template. From the
perspective of physicists, an RNAP is a molecular motor that utilizes chemical
energy input to move along the track formed by a ssDNA. In some circumstances,
which are described in this paper, a large number of RNAPs move simultaneously
along the same track. We refer to such collective movements of the RNAPs as
RNAP traffic because of the similarities between the collective dynamics of the
RNAPs on ssDNA track and that of vehicles in highway traffic. In this paper we
develop a theoretical model for RNAP traffic by incorporating the steric
interactions between RNAPs as well as the mechano-chemical cycle of individual
RNAPs during the elongation of the mRNA. By a combination of analytical and
numerical techniques, we calculate the rates of mRNA synthesis and the average
density profile of the RNAPs on the ssDNA track. We also suggest novel
experiments for testing our theoretical predictions.


  We study the microrheology of nanoparticle shells [Dinsmore et al. Science
298, 1006 (2002)] and viral capsids [Ivanovska et al. PNAS 101, 7600 (2004)] by
computing the mechanical response function and thermal fluctuation spectrum of
a viscoelastic spherical shell that is permeable to the surrounding solvent. We
determine analytically the damped dynamics of the shear, bend, and compression
modes of the shell coupled to the solvent both inside and outside the sphere in
the zero Reynolds number limit. We identify fundamental length and time scales
in the system, and compute the thermal correlation function of displacements of
antipodal points on the sphere and the mechanical response to pinching forces
applied at these points. We describe how such a frequency-dependent antipodal
correlation and/or response function, which should be measurable in new
AFM-based microrheology experiments, can probe the viscoelasticity of these
synthetic and biological shells constructed of nanoparticles.


  We study the length distribution of a particular class of DNA sequences known
as 5'UTR exons. These exons belong to the messanger RNA of protein coding
genes, but they are not coding (they are located upstream of the coding portion
of the mRNA) and are thus less constrained from an evolutionary point of view.
We show that both in mouse and in human these exons show a very clean power law
decay in their length distribution and suggest a simple evolutionary model
which may explain this finding. We conjecture that this power law behaviour
could indeed be a general feature of higher eukaryotes.


  Ongoing sub-cellular oscillation of Min proteins is required to block
minicelling in E. coli. Experimentally, Min oscillations are seen in newly
divided cells and no minicells are produced. In model Min systems many daughter
cells do not oscillate following septation because of unequal partitioning of
Min proteins between the daughter cells. Using the 3D model of Huang et al., we
investigate the septation process in detail to determine the cause of the
asymmetric partitioning of Min proteins between daughter cells. We find that
this partitioning problem arises at certain phases of the MinD and MinE
oscillations with respect to septal closure and it persists independently of
parameter variation. At most 85% of the daughter cells exhibit Min oscillation
following septation. Enhanced MinD binding at the static polar and dynamic
septal regions, consistent with cardiolipin domains, does not substantially
increase this fraction of oscillating daughters. We believe that this problem
will be shared among all existing Min models and discuss possible biological
mechanisms that may minimize partitioning errors of Min proteins following
septation.


  We perform theoretical studies of stretching of 20 proteins with knots within
a coarse grained model. The knot's ends are found to jump to well defined
sequential locations that are associated with sharp turns whereas in
homopolymers they diffuse around and eventually slide off. The waiting times of
the jumps are increasingly stochastic as the temperature is raised. Larger
knots do not return to their native locations when a protein is released after
stretching.


  The utilization of multiple phosphorylation sites in regulating a biological
response is ubiquitous in cell signaling. If each site contributes an
additional, equivalent binding site, then one consequence of an increase in the
number of phosphorylations may be to increase the probability that, upon
disassociation, a ligand immediately rebinds to its receptor. How such effects
may influence cell signaling systems has been less studied. Here, a
self-consistent integral equation formalism for ligand rebinding, in
conjunction with Monte Carlo simulations, is employed to further investigate
the effects of multiple, equivalent binding sites on shaping biological
responses. Multiple regimes that characterize qualitatively different physics
due to the differential prevalence of rebinding effects are predicted.
Calculations suggest that when ligand rebinding contributes significantly to
the dose response, a purely allovalent model can influence the binding curves
nonlinearly. The model also predicts that ligand rebinding in itself appears
insufficient to generative a highly cooperative biological response.


  A pedigree is a directed graph that describes how individuals are related
through ancestry in a sexually-reproducing population. In this paper we explore
the question of whether one can reconstruct a pedigree by just observing
sequence data for present day individuals. This is motivated by the increasing
availability of genomic sequences, but in this paper we take a more theoretical
approach and consider what models of sequence evolution might allow pedigree
reconstruction (given sufficiently long sequences). Our results complement
recent work that showed that pedigree reconstruction may be fundamentally
impossible if one uses just the degrees of relatedness between different extant
individuals. We find that for certain stochastic processes, pedigrees can be
recovered up to isomorphism from sufficiently long sequences.


  The important task of determining the connectivity of gene networks, and at a
more detailed level even the kind of interaction existing between genes, can
nowadays be tackled by microarraylike technologies. Yet, there is still a large
amount of unknowns with respect to the amount of data provided by a single
microarray experiment, and therefore reliable gene network retrieval procedures
must integrate all of the available biological knowledge, even if coming from
different sources and of different nature. In this paper we present a reverse
engineering algorithm able to reveal the underlying gene network by using
time-series dataset on gene expressions considering the system response to
different perturbations. The approach is able to determine the sparsity of the
gene network, and to take into account possible {\it a priori} biological
knowledge on it. The validity of the reverse engineering approach is
highlighted through the deduction of the topology of several {\it simulated}
gene networks, where we also discuss how the performance of the algorithm
improves enlarging the amount of data or if any a priori knowledge is
considered. We also apply the algorithm to experimental data on a nine gene
network in {\it Escherichia coli


  We work in the context of the infinitely many alleles model. The allelic
partition associated with a coalescent process started from n individuals is
obtained by placing mutations along the skeleton of the coalescent tree; for
each individual, we trace back to the most recent mutation affecting it and
group together individuals whose most recent mutations are the same. The number
of blocks of each of the different possible sizes in this partition is the
allele frequency spectrum. The celebrated Ewens sampling formula gives precise
probabilities for the allele frequency spectrum associated with Kingman's
coalescent. This (and the degenerate star-shaped coalescent) are the only
Lambda coalescents for which explicit probabilities are known, although they
are known to satisfy a recursion due to Moehle. Recently, Berestycki,
Berestycki and Schweinsberg have proved asymptotic results for the allele
frequency spectra of the Beta(2-alpha,alpha) coalescents with alpha in (1,2).
In this paper, we prove full asymptotics for the case of the
Bolthausen-Sznitman coalescent.


  A simple lattice model, recently introduced as a generalization of the
Wako--Sait\^o model of protein folding, is used to investigate the properties
of widely studied molecules under external forces. The equilibrium properties
of the model proteins, together with their energy landscape, are studied on the
basis of the exact solution of the model. Afterwards, the kinetic response of
the molecules to a force is considered, discussing both force clamp and dynamic
loading protocols and showing that theoretical expectations are verified. The
kinetic parameters characterizing the protein unfolding are evaluated by using
computer simulations and agree nicely with experimental results, when these are
available. Finally, the extended Jarzynski equality is exploited to investigate
the possibility of reconstructing the free energy landscape of proteins with
pulling experiments.


  We found theoretically that competition between ~Kq^4 and ~Qq^2 terms in the
Fourier transformed conformational energy of a single lipid chain, in
combination with inter-chain entropic repulsion in the hydrophobic part of the
lipid (bi)layer, may cause a crossover on the bilayer pressure-area isotherm
P(A)~(A-A_0)^{-n}. The crossover manifests itself in the transition from n=5/3
to n=3. Our microscopic model represents a single lipid molecule as a worm-like
chain with finite irreducible cross-section area A_0, flexural rigidity K and
stretching modulus Q in a parabolic potential with self-consistent curvature
B(A) formed by entropic interactions between hydrocarbon chains in the lipid
layer. The crossover area per lipid A* obeys relation Q^2/(KB(A*))~1 . We
predict a peculiar possibility to deduce effective elastic moduli K and Q of
the individual hydrocarbon chain from the analysis of the isotherm possessing
such crossover. Also calculated is crossover-related behavior of the area
compressibility modulus K_a, equilibrium area per lipid A_t, and chain order
parameter S.


  This review deals with computer simulation of biological ageing, particularly
with the Penna model of 1995.


  In this paper we present the asymptotic enumeration of RNA structures with
pseudoknots. We develop a general framework for the computation of exponential
growth rate and the sub exponential factors for $k$-noncrossing RNA structures.
Our results are based on the generating function for the number of
$k$-noncrossing RNA pseudoknot structures, ${\sf S}_k(n)$, derived in
\cite{Reidys:07pseu}, where $k-1$ denotes the maximal size of sets of mutually
intersecting bonds. We prove a functional equation for the generating function
$\sum_{n\ge 0}{\sf S}_k(n)z^n$ and obtain for $k=2$ and $k=3$ the analytic
continuation and singular expansions, respectively. It is implicit in our
results that for arbitrary $k$ singular expansions exist and via transfer
theorems of analytic combinatorics we obtain asymptotic expression for the
coefficients. We explicitly derive the asymptotic expressions for 2- and
3-noncrossing RNA structures. Our main result is the derivation of the formula
${\sf S}_3(n) \sim \frac{10.4724\cdot 4!}{n(n-1)...(n-4)}
(\frac{5+\sqrt{21}}{2})^n$.


  It previously has been discovered that visible light irradiation of
crystalline substrates can lead to enhancement of subsequent enzymatic reaction
rates as sharply peaked oscillatory functions of irradiation time. The
particular activating irradiation times can vary with source of a given enzyme
and thus, presumably, its molecular structure. The experiments reported here
demonstrate that the potential for this anomalous enzyme reaction rate
enhancement can be transferred from one bacterial species to another coincident
with transfer of the genetic determinant for the relevant enzyme. In
particular, the effect of crystal-irradiated chloramphenicol on growth of
bacterial strains in which a transferable R-factor DNA plasmid coding for
chloramphenicol resistance was or was not present (S. panama R+, E. coli R+,
and E. coli R-) was determined. Chloramphenicol samples irradiated 10, 35 and
60 sec produced increased growth rates (diminished inhibition) for the
resistant S. panama and E. coli strains, while having no such effect on growth
rate of the sensitive E. coli strain. Consistent with past findings,
chloramphenicol samples irradiated 5, 30 and 55 sec produced decreased growth
rates (increased inhibition) for all three strains.


  With the advent of high-throughput profiling methods, interest in reverse
engineering the structure and dynamics of biochemical networks is high.
Recently an algorithm for reverse engineering of biochemical networks was
developed by Laubenbacher and Stigler. It is a top-down approach using time
discrete dynamical systems. One of its key steps includes the choice of a term
order. The aim of this paper is to identify minimal requirements on data sets
to be used with this algorithm and to characterize optimal data sets. We found
minimal requirements on a data set based on how many terms the functions to be
reverse engineered display. Furthermore, we identified optimal data sets, which
we characterized using a geometric property called "general position".
Moreover, we developed a constructive method to generate optimal data sets,
provided a codimensional condition is fulfilled. In addition, we present a
generalization of their algorithm that does not depend on the choice of a term
order. For this method we derived a formula for the probability of finding the
correct model, provided the data set used is optimal. We analyzed the
asymptotic behavior of the probability formula for a growing number of
variables n (i.e. interacting chemicals). Unfortunately, this formula converges
to zero as fast as r^(q^n), where q is a natural number and 0<r<1. Therefore,
even if an optimal data set is used and the restrictions in using term orders
are overcome, the reverse engineering problem remains unfeasible, unless
prodigious amounts of data are available. Such large data sets are
experimentally impossible to generate with today's technologies.


  Computational analyses of, e.g., genomic, proteomic, or metabolomic data,
commonly result in one or more sets of candidate genes, proteins, or enzymes.
These sets are often the outcome of clustering algorithms. Subsequently, it has
to be tested if, e.g., the candidate gene-products are members of known
metabolic processes. With OrfMapper we provide a powerful but easy-to-use,
web-based database application, that supports such analyses. All services
provided by OrfMapper are freely available at http://www.orfmapper.com


  The spatial structure of populations is a key element in the understanding of
the large scale spreading of epidemics. Motivated by the recent empirical
evidence on the heterogeneous properties of transportation and commuting
patterns among urban areas, we present a thorough analysis of the behavior of
infectious diseases in metapopulation models characterized by heterogeneous
connectivity and mobility patterns. We derive the basic reaction-diffusion
equation describing the metapopulation system at the mechanistic level and
derive an early stage dynamics approximation for the subpopulation invasion
dynamics. The analytical description uses degree block variables that allows us
to take into account arbitrary degree distribution of the metapopulation
network. We show that along with the usual single population epidemic threshold
the metapopulation network exhibits a global threshold for the subpopulation
invasion. We find an explicit analytic expression for the invasion threshold
that determines the minimum number of individuals traveling among
subpopulations in order to have the infection of a macroscopic number of
subpopulations. The invasion threshold is a function of factors such as the
basic reproductive number, the infectious period and the mobility process and
it is found to decrease for increasing network heterogeneity. We provide
extensive mechanistic numerical Monte Carlo simulations that recover the
analytical finding in a wide range of metapopulation network connectivity
patterns. The results can be useful in the understanding of recent data driven
computational approaches to disease spreading in large transportation networks
and the effect of containment measures such as travel restrictions.


  In this paper, we use a variety of mathematical techniques to explore
existence, local stability, and global stability of equilibria in abstract
models of mitochondrial metabolism. The class of models constructed is defined
by the biological description of the system, with minimal mathematical
assumptions. The key features are an electron transport chain coupled to a
process of charge translocation across a membrane. In the absence of charge
translocation these models have previously been shown to behave in a very
simple manner with a single, globally stable equilibrium. We show that with
charge translocation the conclusion about a unique equilibrium remains true,
but local and global stability do not necessarily follow. In sufficiently low
dimensions - i.e. for short electron transport chains - it is possible to make
claims about local and global stability of the equilibrium. On the other hand,
for longer chains, these general claims are no longer valid. Some particular
conditions which ensure stability of the equilibrium for chains of arbitrary
length are presented.


  Cells store information in DNA and in stable programs of gene expression,
which thereby implement forms of long-term cellular memory. Cells must also
possess short-term forms of information storage, implemented
post-translationally, to transduce and interpret external signals. CaMKII, for
instance, is thought to implement a one-bit (bistable) short-term memory
required for learning at post-synaptic densities. Here we show by mathematical
analysis that multisite protein phosphorylation, which is ubiquitous in all
eukaryotic signalling pathways, exhibits multistability for which the maximal
number of steady states increases with the number of sites. If there are n
sites, the maximal information storage capacity is at least log_2 (n+2)/2 bits
when n is even and log_2 (n+1)/2 bits when n is odd. Furthermore, when
substrate is in excess, enzyme saturation together with an alternating low/high
pattern in the site-specific relative catalytic efficiencies, enriches for
multistability. That is, within physiologically plausible ranges for
parameters, multistability becomes more likely than monostability. We discuss
the experimental challenges in pursuing these predictions and in determining
the biological role of short-term information storage.


  Inherent structure theory is used to discover strong connections between
simple characteristics of protein structure and the energy landscape of a Go
model. The potential energies and vibrational free energies of inherent
structures are highly correlated, and both reflect simple measures of networks
of native contacts. These connections have important consequences for models of
protein dynamics and thermodynamics.


  We study the passage times of a translocating polymer of length $N$ in three
dimensions, while it is pulled through a narrow pore with a constant force $F$
applied to one end of the polymer. At small to moderate forces, satisfying the
condition $FN^{\nu}/k_BT\lesssim1$, where $\nu\approx0.588$ is the Flory
exponent for the polymer, we find that $\tau_N$, the mean time the polymer
takes to leave the pore, scales as $N^{2+\nu}$ independent of $F$, in agreement
with our earlier result for F=0. At strong forces, i.e., for
$FN^{\nu}/k_BT\gg1$, the behaviour of the passage time crosses over to
$\tau_N\sim N^2/F$. We show here that these behaviours stem from the polymer
dynamics at the immediate vicinity of the pore -- in particular, the memory
effects in the polymer chain tension imbalance across the pore.


  OBJECTIVES: Changes in the way dietary fat is metabolized can be considered
causative in obesity. The role of sedentary behavior in this defect has not
been determined. We hypothesized that physical inactivity partitions dietary
fats toward storage and that a resistance exercise training program mitigates
storage.


  We solve the Unanimity Rule on networks with exponential, uniform and
scalefree degree distributions. In particular we arrive at equations relating
the asymptotic number of nodes in one of two states to the initial fraction of
nodes in this state. The solutions for exponential and uniform networks are
exact, the approximation for the scalefree case is in perfect agreement with
simulation results. We use these solutions to provide a theoretical
understanding for experimental data on biodiversity loss in foodwebs, which is
available for the three network types discussed. The model allows in principle
to estimate the critical value of species that have to be removed from the
system to induce its complete collapse.


  Silver staining is used to detect proteins after electrophoretic separation
on polyacrylamide gels. It combines excellent sensitivity (in the low nanogram
range) with the use of very simple and cheap equipment and chemicals. It is
compatible with downstream processing, such as mass spectrometry analysis after
protein digestion. The sequential phases of silver staining are protein
fixation, then sensitization, then silver impregnation and finally image
development. Several variants of silver staining are described here, which can
be completed in a time range from 2 h to 1 d after the end of the
electrophoretic separation. Once completed, the stain is stable for several
weeks.


  A probabilistic structure on sequential dynamical systems is introduced here,
the new model will be called Probabilistic Sequential Network, PSN. The
morphisms of Probabilistic Sequential Networks are defined using two algebraic
conditions. It is proved here that two homomorphic Probabilistic Sequential
Networks have the same equilibrium or steady state probabilities if the
morphism is either an epimorphism or a monomorphism. Additionally, the proof of
the set of PSN with its morphisms form the category PSN, having the category of
sequential dynamical systems SDS, as a full subcategory is given. Several
examples of morphisms, subsystems and simulations are given.


  The diversity of virus populations within single infected hosts presents a
major difficulty for the natural immune response as well as for vaccine design
and antiviral drug therapy. Recently developed pyrophosphate based sequencing
technologies (pyrosequencing) can be used for quantifying this diversity by
ultra-deep sequencing of virus samples. We present computational methods for
the analysis of such sequence data and apply these techniques to pyrosequencing
data obtained from HIV populations within patients harboring drug resistant
virus strains. Our main result is the estimation of the population structure of
the sample from the pyrosequencing reads. This inference is based on a
statistical approach to error correction, followed by a combinatorial algorithm
for constructing a minimal set of haplotypes that explain the data. Using this
set of explaining haplotypes, we apply a statistical model to infer the
frequencies of the haplotypes in the population via an EM algorithm. We
demonstrate that pyrosequencing reads allow for effective population
reconstruction by extensive simulations and by comparison to 165 sequences
obtained directly from clonal sequencing of four independent, diverse HIV
populations. Thus, pyrosequencing can be used for cost-effective estimation of
the structure of virus populations, promising new insights into viral
evolutionary dynamics and disease control strategies.


  The problem of intermediates in the fossil record has been frequently
discussed ever since Darwin. The extent of `gaps' (missing transitional stages)
has been used to argue against gradual evolution from a common ancestor.
Traditionally, gaps have often been explained by the improbability of
fossilization and the discontinuous selection of found fossils. Here we take an
analytical approach and demonstrate why, under certain sampling conditions, we
may not expect intermediates to be found. Using a simple null model, we show
mathematically that the question of whether a taxon sampled from some time in
the past is likely to be morphologically intermediate to other samples (dated
earlier and later) depends on the shape and dimensions of the underlying
phylogenetic tree that connects the taxa, and the times from which the fossils
are sampled.


  3' untranslated regions (3' UTRs) contain binding sites for many regulatory
elements, and in particular for microRNAs (miRNAs). The importance of
miRNA-mediated post-transcriptional regulation has become increasingly clear in
the last few years.
  We propose two complementary approaches to the statistical analysis of
oligonucleotide frequencies in mammalian 3' UTRs aimed at the identification of
candidate binding sites for regulatory elements. The first method is based on
the identification of sets of genes characterized by evolutionarily conserved
overrepresentation of an oligonucleotide. The second method is based on the
identification of oligonucleotides showing statistically significant strand
asymmetry in their distribution in 3' UTRs.
  Both methods are able to identify many previously known binding sites located
in 3'UTRs, and in particular seed regions of known miRNAs. Many new candidates
are proposed for experimental verification.


  Xenopus tropicalis is an anuran amphibian species used as model in vertebrate
comparative genomics. It provides the same advantages as Xenopus laevis but is
diploid and has a smaller genome of 1.7 Gbp. Therefore X. tropicalis is more
amenable to systematic transcriptome surveys. We initiated a large-scale
partial cDNA sequencing project to provide a functional genomics resource on
genes expressed in the nervous system during early embryogenesis and
metamorphosis in X. tropicalis. A gene index was defined and analysed after the
collection of over 48,785 high quality sequences. Partial cDNA sequences were
obtained from an embryonic head and retina library (30,272 sequences) and from
a metamorphic brain and spinal cord library (27,602 sequences). These ESTs are
estimated to represent 9,693 transcripts derived from an estimated 6,000 genes.
An estimated 46% of these cDNA sequences contain their start codon. Further
annotation included Gene Ontology functional classification, InterPro domain
analysis, alternative splicing and non-coding RNA identification. Gene
expression profiles were derived from EST counts and used to define transcripts
specific to metamorphic stages of development. Moreover, these ESTs allowed
identification of a set of 225 polymorphic microsatellites that can be used as
genetic markers. These cDNA sequences permit in silico cloning of numerous
genes and will facilitate studies aimed at deciphering the roles of cognate
genes expressed in the nervous system during neural development and
metamorphosis. The genomic resources developed to study X. tropicalis biology
will accelerate exploration of amphibian physiology and genetics.


  Experimental variables of optical tweezers instrumentation that affect RNA
folding/unfolding kinetics were investigated. A model RNA hairpin, P5ab, was
attached to two micron-sized beads through hybrid RNA/DNA handles; one bead was
trapped by dual-beam lasers and the other was held by a micropipette. Several
experimental variables were changed while measuring the unfolding/refolding
kinetics, including handle lengths, trap stiffness, and modes of force applied
to the molecule. In constant-force mode where the tension applied to the RNA
was maintained through feedback control, the measured rate coefficients varied
within 40% when the handle lengths were changed by 10 fold (1.1 to 10.2 Kbp);
they increased by two- to three-fold when the trap stiffness was lowered to one
third (from 0.1 to 0.035 pN/nm). In the passive mode, without feedback control
and where the force applied to the RNA varied in response to the end-to-end
distance change of the tether, the RNA hopped between a high-force folded-state
and a low-force unfolded-state. In this mode, the rates increased up to
two-fold with longer handles or softer traps. Overall, the measured rates
remained with the same order-of-magnitude over the wide range of conditions
studied. In the companion paper (1), we analyze how the measured kinetics
parameters differ from the intrinsic molecular rates of the RNA, and thus how
to obtain the molecular rates.


  By exerting mechanical force it is possible to unfold/refold RNA molecules
one at a time. In a small range of forces, an RNA molecule can hop between the
folded and the unfolded state with force-dependent kinetic rates. Here, we
introduce a mesoscopic model to analyze the hopping kinetics of RNA hairpins in
an optical tweezers setup. The model includes different elements of the
experimental setup (beads, handles and RNA sequence) and limitations of the
instrument (time lag of the force-feedback mechanism and finite bandwidth of
data acquisition). We investigated the influence of the instrument on the
measured hopping rates. Results from the model are in good agreement with the
experiments reported in the companion article (1). The comparison between
theory and experiments allowed us to infer the values of the intrinsic
molecular rates of the RNA hairpin alone and to search for the optimal
experimental conditions to do the measurements. We conclude that long handles
and soft laser traps represent the best conditions to extract rate estimates
that are closest to the intrinsic molecular rates. The methodology and
rationale presented here can be applied to other experimental setups and other
molecules.


  Part 1 of the study intends to show that the universal trend of amino acid
gain and loss discovered by Jordan et al. (2005) can be accounted for by the
spontaneity of DNA typical damages. These damages lead to replacements of
guanine and cytosine by thymine. Part 2 proposes a hypothesis of the evolution
of the genetic code, the leading mechanism of which is the nucleotide
spontaneous damage. The hypothesis accounts for the universal trend of amino
acid gain and loss, stability of the genetic code towards point mutations, the
presence of code dialects, and the symmetry of the genetic code table.


  Degeneracy of the genetic code is a biological way to minimize effects of the
undesirable mutation changes. Degeneration has a natural description on the
5-adic space of 64 codons $\mathcal{C}_5 (64) = \{n_0 + n_1 5 + n_2 5^2
  : n_i = 1, 2, 3, 4 \} ,$ where $n_i$ are digits related to nucleotides as
follows: C = 1, A = 2, T = U = 3, G = 4. The smallest 5-adic distance between
codons joins them into 16 quadruplets, which under 2-adic distance decay into
32 doublets. p-Adically close codons are assigned to one of 20 amino acids,
which are building blocks of proteins, or code termination of protein
synthesis. We shown that genetic code multiplets are made of the p-adic nearest
codons.


  In birds, the link between parental care behaviour and prolactin release
during incubation persists after hatching in altricial birds, but has never
been precisely studied during the whole rearing period in precocial species,
such as ducks. The present study aims to understand how changes in parental
care after hatching are related to circulating prolactin levels in mallard hens
rearing ducklings. Blood was sampled in hens over at least 13 post-hatching
weeks and the behaviour of the hens and the ducklings was recorded daily until
fledging. Contacts between hens and the ducklings, leadership of the ducklings
and gathering of them steadily decreased over post-hatching time. Conversely,
resting, preening and agonistic behaviour of hens towards ducklings increased.
Plasma prolactin concentrations remained at high levels after hatching and then
fell after week 6 when body mass and structural size of the young were close to
those of the hen. Parental care behaviour declined linearly with brood age,
showed a disruption of the hen-brood bond at week 6 post-hatching and was
related to prolactin concentration according to a sigmoid function. Our results
suggest that a definite threshold in circulating prolactin is necessary to
promote and/or to maintain post-hatching parental care in ducks.


This paper has been withdrawn by the author


  Amyotrophic lateral sclerosis (ALS) is the most common adult motor neuron
disease, causing motor neuron degeneration, muscle atrophy, paralysis, and
death. Despite this degenerative process, a stable hypermetabolic state has
been observed in a large subset of patients. Mice expressing a mutant form of
Cu/Zn-superoxide dismutase (mSOD1 mice) constitute an animal model of ALS that,
like patients, exhibits unexpectedly increased energy expenditure.
Counterbalancing for this increase with a high-fat diet extends lifespan and
prevents motor neuron loss. Here, we investigated whether lipid metabolism is
defective in this animal model. Hepatic lipid metabolism was roughly normal,
whereas gastrointestinal absorption of lipids as well as peripheral clearance
of triglyceride-rich lipoproteins were markedly increased, leading to decreased
postprandial lipidemia. This defect was corrected by the high-fat regimen that
typically induces neuroprotection in these animals. Together, our findings show
that energy metabolism in mSOD1 mice shifts toward an increase in the
peripheral use of lipids. This metabolic shift probably accounts for the
protective effect of dietary lipids in this model.


  The success of an infectious disease to invade a population is strongly
controlled by the population's specific connectivity structure. Here a network
model is presented as an aid in understanding the role of social behavior and
heterogeneous connectivity in determining the spatio-temporal patterns of
disease dynamics. We explore the controversial origins of long-term recurrent
oscillations believed to be characteristic to diseases that have a period of
temporary immunity after infection. In particular, we focus on sexually
transmitted diseases such as syphilis, where this controversy is currently
under review. Although temporary immunity plays a key role, it is found that in
realistic small-world networks, the social and sexual behavior of individuals
also has great influence in generating long-term cycles. The model generates
circular waves of infection with unusual spatial dynamics that depend on focal
areas that act as pacemakers in the population. Eradication of the disease can
be efficiently achieved by eliminating the pacemakers with a targeted
vaccination scheme. A simple difference equation model is derived, that
captures the infection dynamics of the network model and gives insights into
their origins and their eradication through vaccination.


  Sequence alignment forms the basis of many methods for functional annotation
by phylogenetic comparison, but becomes unreliable in the `twilight' regions of
high sequence divergence and short gene length. Here we perform a cross-species
comparison of two herpesviruses, VZV and KSHV, with a hybrid method called
graph alignment. The method is based jointly on the similarity of protein
interaction networks and on sequence similarity. In our alignment, we find open
reading frames for which interaction similarity concurs with a low level of
sequence similarity, thus confirming the evolutionary relationship. In
addition, we find high levels of interaction similarity between open reading
frames without any detectable sequence similarity. The functional predictions
derived from this alignment are consistent with genomic position and gene
expression data.


  The free-energy landscape of the alpha-helix of protein G is studied by means
of metadynamics coupled with a solute tempering algorithm. Metadynamics allows
to overcome large energy barriers, whereas solute tempering improves the
sampling with an affordable computational effort. From the sampled free-energy
surface we are able to reproduce a number of experimental observations, such as
the fact that the lowest minimum corresponds to a globular conformation
displaying some degree of beta-structure, that the helical state is metastable
and involves only 65% of the chain. The calculations also show that the system
populates consistently a pi-helix state and that the hydrophobic staple motif
is present only in the free-energy minimum associated with the helices, and
contributes to their stabilization. The use of metadynamics coupled with solute
tempering results then particularly suitable to provide the thermodynamics of a
short peptide, and its computational efficiency is promising to deal with
larger proteins.


  The lac operon of Escherichia coli exhibits bistability. Early studies showed
that bistability occurs during growth on TMG/succinate and lactose + glucose,
but not during growth on lactose. More recent studies with lacGFP-transfected
cells show bistability with TMG/succinate, but not with lactose and lactose +
glucose. In the literature, these results are attributed to variations of the
positive feedback generated by induction. Specifically, during growth on
TMG/succinate, induction generates positive feedback because the permease
stimulates the accumulation of TMG, which, in turn, promotes the synthesis of
more permease. This positive feedback is attenuated during growth on lactose
because hydrolysis of lactose by galactosidase suppresses the stimulatory
effect of the permease. But the stabilizing effect of dilution also changes
dramatically as a function of the medium composition. For instance, during
growth on TMG/succinate, the dilution rate of the permease is proportional to
its activity, $e$, because the specific growth rate is independent of $e$.
However, during growth on lactose, the permease dilution rate is proportional
to $e^2$ because the specific growth rate is proportional to the specific
lactose uptake rate, which in turn, proportional to $e$. Here, we show that:
(a) This dependence on $e^2$ creates such a strong stabilizing effect that
bistability is virtually impossible during growth on lactose, even in the face
of positive feedback. (b) This stabilizing effect is weakened during growth on
lactose + glucose because the specific growth rate on glucose is independent of
$e$, so that the dilution rate once again contains a term that is proportional
to $e$. We discuss the experimental data in the light of these results.


  Control of the living cell functions with remarkable reliability despite the
stochastic nature of the underlying molecular networks -- a property presumably
optimized by biological evolution. We here ask to what extent the property of a
stochastic dynamical network to produce reliable dynamics is an evolvable
trait. Using an evolutionary algorithm based on a deterministic selection
criterion for the reliability of dynamical attractors, we evolve dynamical
networks of noisy discrete threshold nodes. We find that, starting from any
random network, reliability of the attractor landscape can often be achieved
with only few small changes to the network structure. Further, the evolvability
of networks towards reliable dynamics while retaining their function is
investigated and a high success rate is found.


  As in many other areas of science, systems biology makes extensive use of
statistical association and significance estimates in contingency tables, a
type of categorical data analysis known in this field as enrichment (also
over-representation or enhancement) analysis. In spite of efforts to create
probabilistic annotations, especially in the Gene Ontology context, or to deal
with uncertainty in high throughput-based datasets, current enrichment methods
largely ignore this probabilistic information since they are mainly based on
variants of the Fisher Exact Test. We developed an open-source R package to
deal with probabilistic categorical data analysis, ProbCD, that does not
require a static contingency table. The contingency table for the enrichment
problem is built using the expectation of a Bernoulli Scheme stochastic process
given the categorization probabilities. An on-line interface was created to
allow usage by non-programmers and is available at:
http://xerad.systemsbiology.net/ProbCD/ . We present an analysis framework and
software tools to address the issue of uncertainty in categorical data
analysis. In particular, concerning the enrichment analysis, ProbCD can
accommodate: (i) the stochastic nature of the high-throughput experimental
techniques and (ii) probabilistic gene annotation.


  Apoptosis, or genetically programmed cell death, is a crucial cellular
process that maintains the balance between life and death in cells. The precise
molecular mechanism of apoptosis signaling and how these two pathways are
differentially activated under distinct apoptotic stimuli is poorly understood.
We developed a Monte Carlo-based stochastic simulation model that can
characterize distinct signaling behaviors in the two major pathways of
apoptotic signaling using a novel probability distribution-based approach.
Specifically, we show that for a weak death signal, such as low levels of death
ligand Fas (CD95) binding or under stress conditions, the type 2 mitochondrial
pathway dominates apoptotic signaling. Our results also show signaling in the
type 2 pathway is stochastic, where the population average over many cells does
not capture the cell-to-cell fluctuations in the time course (~1 - 10 hours) of
downstream caspase-3 activation. On the contrary, the probability distribution
of caspase-3 activation for the mitochondrial pathway shows a distinct bimodal
behavior that can be used to characterize the stochastic signaling in type 2
apoptosis. Interestingly, such stochastic fluctuations in apoptosis signaling
happen even in the presence of large numbers of signaling molecules. In a
fluctuating environment, such stochasticity in the timecourse of caspase-3
activation may be an adaptive mechanism for allowing a competing survival
signal to win over a weak death trigger before the critical cell fate decision
is made and thus minimizes the risk of pathologies.


  Chemical reactions in cell are subject to intense stochastic fluctuations. An
important question is how the fundamental physiological behavior of cell is
kept stable against those noisy perturbations. In this paper a stochastic model
of cell cycle of budding yeast is constructed to analyze the effects of noise
on the cell cycle oscillation. The model predicts intense noise in levels of
mRNAs and proteins, and the simulated protein levels explain the observed
statistical tendency of noise in populations of synchronous and asynchronous
cells. In spite of intense noise in levels of proteins and mRNAs, cell cycle is
stable enough to bring the largely perturbed cells back to the physiological
cyclic oscillation. The model shows that consecutively appearing fixed points
are the origin of this stability of cell cycle.


  The possible coexistence of one host, one aggressive parasite and one
non-lethal parasite is simulated using the Penna model of biological ageing. If
the aggressive parasites survive the difficult initial times where they have to
adjust genetically to the proper host age, all three species may survive,
though the host number may be diminished by increasing parasite aggressivity.


  We present a computer simulation, which is inspired by Penna model, to help
understanding the effect of genetic coding tables on population dynamics. To
represent populations we used real and artificial gene sequences in this model.
We coded these sequences using different amino acid tables in Nature, the
standard table as well as the tables which are used by mithocondria and some
eukaryotes. Contrary to common belief we find that the standard code table
which is used in most organisms in Nature, does not give the most resilient
coding against point mutations.


  We study a stochastic predator-prey model on a square lattice, where each of
the six species has two superior and two inferior partners. The invasion
probabilities between species depend on the predator-prey pair and are
supplemented by Gaussian noise. Conditions are identified that warrant the
largest impact of noise on the evolutionary process, and the results of Monte
Carlo simulations are qualitatively reproduced by a four-point cluster
dynamical mean-field approximation. The observed noise-guided evolution is
deeply routed in short-range spatial correlations, which is supported by
simulations on other host lattice topologies. Our findings are conceptually
related to the coherence resonance phenomenon in dynamical systems via the
mechanism of threshold duality. We also show that the introduced concept of
noise-guided evolution via the exploitation of threshold duality is not limited
to predator-prey cyclical interactions, but may apply to models of evolutionary
game theory as well, thus indicating its applicability in several different
fields of research.


  Magnetic tweezers are used to study the mechanical response under torsion of
single nucleosome arrays reconstituted on tandem repeats of 5S positioning
sequences. Regular arrays are extremely resilient and can reversibly
accommodate a large amount of supercoiling without much change in length. This
behavior is quantitatively described by a molecular model of the chromatin 3-D
architecture. In this model, we assume the existence of a dynamic equilibrium
between three conformations of the nucleosome, which are determined by the
crossing status of the entry/exit DNAs (positive, null or negative). Torsional
strain, in displacing that equilibrium, extensively reorganizes the fiber
architecture. The model explains a number of long-standing topological
questions regarding DNA in chromatin, and may provide the ground to better
understand the dynamic binding of most chromatin-associated proteins.


  We propose a biologically motivated quantity, twinness, to evaluate local
similarity between nodes in a network. The twinness of a pair of nodes is the
number of connected, labeled subgraphs of size n in which the two nodes possess
identical neighbours. The graph animal algorithm is used to estimate twinness
for each pair of nodes (for subgraph sizes n=4 to n=12) in four different
protein interaction networks (PINs). These include an Escherichia coli PIN and
three Saccharomyces cerevisiae PINs -- each obtained using state-of-the-art
high throughput methods. In almost all cases, the average twinness of node
pairs is vastly higher than expected from a null model obtained by switching
links. For all n, we observe a difference in the ratio of type A twins (which
are unlinked pairs) to type B twins (which are linked pairs) distinguishing the
prokaryote E. coli from the eukaryote S. cerevisiae. Interaction similarity is
expected due to gene duplication, and whole genome duplication paralogues in S.
cerevisiae have been reported to co-cluster into the same complexes. Indeed, we
find that these paralogous proteins are over-represented as twins compared to
pairs chosen at random. These results indicate that twinness can detect
ancestral relationships from currently available PIN data.


  Using a time-dependent perturbation theory, vibrational energy relaxation
(VER) of isotopically labeled amide I modes in cytochrome c solvated with water
is investigated. Contributions to the VER are decomposed into two contributions
from the protein and water. The VER pathways are visualized using radial and
angular excitation functions for resonant normal modes. Key differences of VER
among different amide I modes are demonstrated, leading to a detailed picture
of the spatial anisotropy of the VER. The results support the experimental
observation that amide I modes in proteins relax with sub picosecond
timescales, while the relaxation mechanism turns out to be sensitive to the
environment of the amide I mode.


  Phylogenetic networks are a generalization of phylogenetic trees that allow
for the representation of non-treelike evolutionary events, like recombination,
hybridization, or lateral gene transfer. In a recent series of papers devoted
to the study of reconstructibility of phylogenetic networks, Moret, Nakhleh,
Warnow and collaborators introduced the so-called {tripartition metric for
phylogenetic networks. In this paper we show that, in fact, this tripartition
metric does not satisfy the separation axiom of distances (zero distance means
isomorphism, or, in a more relaxed version, zero distance means
indistinguishability in some specific sense) in any of the subclasses of
phylogenetic networks where it is claimed to do so. We also present a subclass
of phylogenetic networks whose members can be singled out by means of their
sets of tripartitions (or even clusters), and hence where the latter can be
used to define a meaningful metric.


  The distribution of genetic polymorphisms in a population contains
information about the mutation rate and the strength of natural selection at a
locus. Here, we show that the Poisson Random Field (PRF) method of
population-genetic inference suffers from systematic biases that tend to
underestimate selection pressures and mutation rates, and that erroneously
infer positive selection. These problems arise from the infinite-sites
approximation inherent in the PRF method. We introduce three new inference
techniques that correct these problems. We present a finite-site modification
of the PRF method, as well as two new methods for inferring selection pressures
and mutation rates based on diffusion models. Our methods can be used to infer
not only a "weighted average" of selection pressures acting on a gene sequence,
but also the distribution of selection pressures across sites. We evaluate the
accuracy of our methods, as well that of the original PRF approach, by
comparison with Wright-Fisher simulations.


  Surface modified amorphous nanoporous silica molecules with hydrophobic as
well as hydrophilic character can be effectively used as therapeutic drug for
combating chicken malaria in poultry industry. The amorphous nanosilica was
developed by top-down approach using volcanic soil derived silica as source
material. Amorphous silica has long been used as feed additive for poultry
industry and considered to be safe for human consumption by WHO and USDA. The
basic mechanism of action of these nanosilica molecules is mediated by the
physical absorption of VLDL, serum triglycerides and other serum cholesterol
components in the lipophilic nanopores of nanosilica. This reduces the supply
of the host derived cholesterol, thus limiting the growth of the malarial
parasite in vivo.


  Biochemistry and mechanics are closely coupled in cell adhesion. At sites of
cell-matrix adhesion, mechanical force triggers signaling through the
Rho-pathway, which leads to structural reinforcement and increased
contractility in the actin cytoskeleton. The resulting force acts back to the
sites of adhesion, resulting in a positive feedback loop for mature adhesion.
Here we model this biochemical-mechanical feedback loop for the special case
when the actin cytoskeleton is organized in stress fibers, which are
contractile bundles of actin filaments. Activation of myosin II molecular
motors through the Rho-pathway is described by a system of reaction-diffusion
equations, which are coupled into a viscoelastic model for a contractile actin
bundle. We find strong spatial gradients in the activation of contractility and
in the corresponding deformation pattern of the stress fiber, in good agreement
with experimental findings.


  One of the most direct human mechanisms of promoting cooperation is rewarding
it. We study the effect of sharing a reward among cooperators in the most
stringent form of social dilemma, namely the Prisoner's Dilemma. Specifically,
for a group of players that collect payoffs by playing a pairwise Prisoner's
Dilemma game with their partners, we consider an external entity that
distributes a fixed reward equally among all cooperators. Thus, individuals
confront a new dilemma: on the one hand, they may be inclined to choose the
shared reward despite the possibility of being exploited by defectors; on the
other hand, if too many players do that, cooperators will obtain a poor reward
and defectors will outperform them. By appropriately tuning the amount to be
shared a vast variety of scenarios arises, including traditional ones in the
study of cooperation as well as more complex situations where unexpected
behavior can occur. We provide a complete classification of the equilibria of
the $n$-player game as well as of its evolutionary dynamics.


  Malaria and other parasites, including virus often induce an increase in host
lipids which the invaders use to their own advantage. We obtained encouraging
results in our investigations on bird malaria with a new approach namely the
use of nanosilica to mop up excess host lipids. While this project is
continuing we have investigated another, simpler system namely silkworms which
suffer from a deadly baculovirus, BmNPV. This virus decimates the infected
population within 24 hours or so and no known antibiotic antidote or
genetically resistant strain of silkworm3 exists. We report here a partial
success, which is worth following up. Our rationale, we believe, has a broad
and interdisciplinary appeal, for, this nanosilica treatment might be used
together with other arsenals on all sorts of virus which take advantage of
enhanced host lipids. It has not escaped our notice that Ebola and HIV also
belong to this category. Nanoparticles are being preferentially harnessed,
because they offer a greater surface area, circulate more easily and in
lepidopteran system4 they are removed within 24 hours from the body. Lawry
surmised, on cogent theoretical grounds that particles significantly smaller
than micron order would be less harmful in the hemocoele. Furthermore, Hui-peng
et al. pointed out that lipase treatment, the only viable option for
controlling BmNPV interferes in hormonal balance and cannot be applied to pre
molting stage.


  Various types of surface functionalized nanosilica (50-60 nm size with 3-10
nm inner pore size range) have been used to kill insect pests by sucking up
cuticular lipids and breaking the water barrier. We have also utilized
nanosilica for mopping up host lipids induced by the malarial parasite, P.
gallinaceum in poultry birds; VLDL cholesterol and serum triglycerides are
brought back to the normal level with a concomitant check in parasite growth.
While this work continues, we have explored another more convenient system,
silkworm (Bombyx mori) that is frequently decimated by a baculovirus, NPV for
which no antidote is known so far. Here, too, viral infection enhances host
lipids. Eight different types of nanosilica were injected in the virus infected
silkworm (batches of 10 worms) after ensuring 100% survival up to cocoon
formation in control larvae (injected with the same volume of ethanol, the
medium of nanosilica). Of these 8, AL60102 and AL60106, have the most marked
effect on infected silkworm, both as prophylactic and pharmaceutical agents.
Normal larvae injected with these nanosilica survive up to cocoon formation.


  Wild animals, pets, zoo animals and mammals of veterinary importance heavily
suffer from trypanosomiasis. Drugs with serious side effects are currently
mainstay of therapies used by veterinarians. Trypanosomiasis is caused by
Trypanosoma sp. leading to sleeping sickness in humans. Surface modified
(hydrophobic and lipophilic) amorphous nanoporous silica molecules could be
effectively used as therapeutic drug for combating trypanosomiasis. The
amorphous nanosilica was developed by top-down approach using volcanic soil
derived silica (Advasan; 50- 60 nm size with 3-10 nm inner pore size range) and
diatomaceous earth (FS; 60-80 nm size with 3-5 nm inner pore size range) as
source materials. According to WHO and USDA standards amorphous silica has long
been used as feed additives for several veterinary industries and considered to
be safe for human consumption. The basic mechanism of action of these
nanosilica molecules is mediated by the physical absorption of HDL components
in the lipophilic nanopores of nanosilica. This reduces the supply of the host
derived cholesterol, thus limiting the growth of the Trypanosoma sp. in vivo.


  We investigate the cycles in the transcription network of S. cerevisiae.
Unlike a similar network of E. coli, it contains many cycles. We characterize
properties of these cycles and their place in the regulatory mechanism of the
cell. Almost all cycles in the transcription network of S. cerevisiae are
contained in a single strongly connected component, which we call LSCC (L for
``largest''), except for a single cycle of two transcription factors. Among
different physiological conditions, cell cycle has the most significant
relationship with LSCC, as the set of 64 transcription interactions that are
active in all phases of the cell cycle has overlap of 27 with the interactions
of LSCC (of which there are 49). Conversely, if we remove the interactions that
are active in all phases of the cell cycle (fewer than 1% of the total), the
LSCC would have only three nodes and 5 edges, 4 of which are active only in the
stress response subnetwork. LSCC has a special place in the topology of the
network and it can be used to define a natural hierarchy in the network; in
every physiological subnetwork LSCC plays a pivotal role. Apart from those
well-defined conditions, the transcription network of S. cerevisiae is devoid
of cycles. It was observed that two conditions that were studied and that have
no cycles of their own are exogenous: diauxic shift and DNA repair, while cell
cycle, sporulation are endogenous.


  Jansson and Sung showed that, given a dense set of input triplets T
(representing hypotheses about the local evolutionary relationships of triplets
of species), it is possible to determine in polynomial time whether there
exists a level-1 network consistent with T, and if so to construct such a
network. They also showed that, unlike in the case of trees (i.e. level-0
networks), the problem becomes NP-hard when the input is non-dense. Here we
further extend this work by showing that, when the set of input triplets is
dense, the problem is even polynomial-time solvable for the construction of
level-2 networks. This shows that, assuming density, it is tractable to
construct plausible evolutionary histories from input triplets even when such
histories are heavily non-tree like. This further strengthens the case for the
use of triplet-based methods in the construction of phylogenetic networks. We
also show that, in the non-dense case, the level-2 problem remains NP-hard.


  Metapopulation models describing cities with different populations coupled by
the travel of individuals are of great importance in the understanding of
disease spread on a large scale. An important example is the Rvachev-Longini
model [{\it Math. Biosci.} {\bf 75}, 3-22 (1985)] which is widely used in
computational epidemiology. Few analytical results are however available and in
particular little is known about paths followed by epidemics and disease
arrival times. We study the arrival time of a disease in a city as a function
of the starting seed of the epidemics. We propose an analytical Ansatz, test it
in the case of a spreading on the world wide air transportation network, and
show that it predicts accurately the arrival order of a disease in world-wide
cities.


  Long range charge transfer experiments in DNA oligomers and the subsequently
measured -- and very diverse -- transport response of DNA wires in solid state
experiments exemplifies the need for a thorough theoretical understanding of
charge migration in DNA-based natural and artificial materials. Here we present
a review of tight-binding models for DNA conduction which have the intrinsic
merit of containing more structural information than plain rate-equation models
while still retaining sufficient detail of the electronic properties. This
allows for simulations of transport properties to be more manageable with
respect to density functional theory methods or correlated first principle
algorithms.


  Estimates of molecular divergence times when rates of evolution vary require
the assumption of a model of rate change. Brownian motion is one such model,
and since rates cannot become negative, a log Brownian model seems appropriate.
Divergence time estimates can then be made using weighted least squares
penalties. As sequences become long, this approach effectively becomes
equivalent to penalized likelihood or Bayesian approaches. Different forms of
the least squares penalty are considered to take into account correlation due
to shared ancestors. It is shown that a scale parameter is also needed since
the sum of squares changes with the scale of time. Errors or uncertainty on
fossil calibrations, may be folded in with errors due to the stochastic nature
of Brownian motion and ancestral polymorphism, giving a total sum of squares to
be minimized. Applying these methods to placental mammal data the estimated age
of the root decreases from 125 to about 94 mybp. However, multiple fossil
calibration points and relative molecular divergence times inflate the sum of
squares more than expected. If fossil data are also bootstrapped, then the
confidence interval for the root of placental mammals varies widely from ~70 to
130 mybp. Such a wide interval suggests that more and better fossil calibration
data is needed and/or better models of rate evolution are needed and/or better
molecular data are needed. Until these issues are thoroughly investigated, it
is premature to declare either the old molecular dates frequently obtained
(e.g. > 110 mybp) or the lack of identified placental fossils in the
Cretaceous, more indicative of when crown-group placental mammals evolved.


  This paper studies the mutation-selection balance in three simplified
replication models. The first model considers a population of organisms
replicating via the production of asexual spores. The second model considers a
sexually replicating population that produces identical gametes. The third
model considers a sexually replicating population that produces distinct sperm
and egg gametes. All models assume diploid organisms whose genomes consist of
two chromosomes, each of which is taken to be functional if equal to some
master sequence, and defective otherwise. In the asexual population, the
asexual diploid spores develop directly into adult organisms. In the sexual
populations, the haploid gametes enter a haploid pool, where they may fuse with
other haploids. The resulting immature diploid organisms then proceed to
develop into mature organisms. Based on an analysis of all three models, we
find that, as organism size increases, a sexually replicating population can
only outcompete an asexually replicating population if the adult organisms
produce distinct sperm and egg gametes. A sexual replication strategy that is
based on the production of large numbers of sperm cells to fertilize a small
number of eggs is found to be necessary in order to maintain a sufficiently low
cost for sex for the strategy to be selected for over a purely asexual
strategy. We discuss the usefulness of this model in understanding the
evolution and maintenance of sexual replication as the preferred replication
strategy in complex, multicellular organisms.


  In a recent article, Desai and Fisher (2007) proposed that the speed of
adaptation in an asexual population is determined by the dynamics of the
stochastic edge of the population, that is, by the emergence and subsequent
establishment of rare mutants that exceed the fitness of all sequences
currently present in the population. Desai and Fisher perform an elaborate
stochastic calculation of the mean time $\tau$ until a new class of mutants has
been established, and interpret $1/\tau$ as the speed of adaptation. As they
note, however, their calculations are valid only for moderate speeds. This
limitation arises from their method to determine $\tau$: Desai and Fisher
back-extrapolate the value of $\tau$ from the best-fit class' exponential
growth at infinite time. This approach is not valid when the population adapts
rapidly, because in this case the best-fit class grows non-exponentially during
the relevant time interval. Here, we substantially extend Desai and Fisher's
analysis of the stochastic edge. We show that we can apply Desai and Fisher's
method to high speeds by either exponentially back-extrapolating from finite
time or using a non-exponential back-extrapolation. Our results are compatible
with predictions made using a different analytical approach (Rouzine et al.
2003, 2007), and agree well with numerical simulations.


  Feedback loops in a dynamic network play an important role in determining the
dynamics of that network. Through a computational study, in this paper we show
that networks with fewer independent negative feedback loops tend to exhibit
more regular behavior than those with more negative loops. To be precise, we
study the relationship between the number of independent feedback loops and the
number and length of the limit cycles in the phase space of dynamic Boolean
networks. We show that, as the number of independent negative feedback loops
increases, the number (length) of limit cycles tends to decrease (increase).
These conclusions are consistent with the fact, for certain natural biological
networks, that they on the one hand exhibit generally regular behavior and on
the other hand show less negative feedback loops than randomized networks with
the same numbers of nodes and connectivity.


  We use traveling-wave theory to derive expressions for the rate of
accumulation of deleterious mutations under Muller's ratchet and the speed of
adaptation under positive selection in asexual populations. Traveling-wave
theory is a semi-deterministic description of an evolving population, where the
bulk of the population is modeled using deterministic equations, but the class
of the highest-fitness genotypes, whose evolution over time determines loss or
gain of fitness in the population, is given proper stochastic treatment. We
derive improved methods to model the highest-fitness class (the stochastic
edge) for both Muller's ratchet and adaptive evolution, and calculate analytic
correction terms that compensate for inaccuracies which arise when treating
discrete fitness classes as a continuum. We show that traveling wave theory
makes excellent predictions for the rate of mutation accumulation in the case
of Muller's ratchet, and makes good predictions for the speed of adaptation in
a very broad parameter range. We predict the adaptation rate to grow
logarithmically in the population size until the population size is extremely
large.