  This paper is to show that most discrete models used for population dynamics
in ecology are inherently pathological that their predications cannot be
independently verified by experiments because they violate a fundamental
principle of physics. The result is used to tackle an on-going controversy
regarding ecological chaos. Another implication of the result is that all
continuous dynamical systems must be modeled by differential equations. As a
result it suggests that researches based on discrete modeling must be closely
scrutinized and the teaching of calculus and differential equations must be
emphasized for students of biology.


  The general Markov plus invariable sites (GM+I) model of biological sequence
evolution is a two-class model in which an unknown proportion of sites are not
allowed to change, while the remainder undergo substitutions according to a
Markov process on a tree. For statistical use it is important to know if the
model is identifiable; can both the tree topology and the numerical parameters
be determined from a joint distribution describing sequences only at the leaves
of the tree? We establish that for generic parameters both the tree and all
numerical parameter values can be recovered, up to clearly understood issues of
`label swapping.' The method of analysis is algebraic, using phylogenetic
invariants to study the variety defined by the model. Simple rational formulas,
expressed in terms of determinantal ratios, are found for recovering numerical
parameters describing the invariable sites.


  In this paper we present a discussion of a phenomenological model of the MAPK
cascade which was originally proposed by Angeli et al. (PNAS 101, 1822 (2004)).
The model and its solution are extended in several respects: a) an analytical
solution is given for the cascade equilibria, exploiting a parameter-based
symmetry of the rate equations; b) we discuss the cooperativity (Hill
coefficients) of the cascade and show that a feedforward loop within the
cascade increases its cooperativity. The relevance of this result for the
notion of modularity is discussed; c) the feedback model for cascade
bistability by Angeli et al. is reconsidered. We argue that care must be taken
in modeling the interactions and a biologically realistic phenomenological
model cannot be too reductionist. The inclusion of a time-dependent degradation
rate is needed to account for a switching of the cascade.


  A large number of theoretical studies predict that the dynamics of spatially
structured populations (metapopulations) can be altered by constant
perturbations to local population size. However, these studies presume large
metapopulations inhabiting noise-free, zero-extinction environments, and their
predictions have never been empirically verified.
  Here we report an empirical study on the effects of localized perturbations
on global dynamics and stability, using fruitfly metapopulations in the
laboratory. We find that constant addition of individuals to a particular
subpopulation in every generation stabilizes that subpopulation locally, but
does not have any detectable effect on the dynamics and stability of the
metapopulation. Simulations of our experimental system using a simple but
widely applicable model of population dynamics were able to recover the
empirical findings, indicating the generality of our results. We then simulated
the possible consequences of perturbing more subpopulations, increasing the
strength of perturbations, and varying the rate of migration, but found that
none of these conditions were expected to alter the outcomes of our
experiments. Finally, we show that our main results are robust to the presence
of local extinctions in the metapopulation.
  Our study shows that localized perturbations are unlikely to affect the
dynamics of real metapopulations, a finding that has cautionary implications
for ecologists and conservation biologists faced with the problem of
stabilizing unstable metapopulations in nature.


  Finding the dominant direction of flow of information in densely
interconnected regulatory or signaling networks is required in many
applications in computational biology and neuroscience. This is achieved by
first identifying and removing links which close up feedback loops in the
original network and hierarchically arranging nodes in the remaining network.
In mathematical language this corresponds to a problem of making a graph
acyclic by removing as few links as possible and thus altering the original
graph in the least possible way. Practically in all applications the exact
solution of this problem requires an enumeration of all combinations of removed
links, which is computationally intractable. We introduce and compare two
algorithms: the deterministic, 'greedy' algorithm that preferentially cuts the
links that participate in the largest number of feedback cycles, and the
probabilistic one based on a simulated annealing of a hierarchical layout of
the network which minimizes the number of ``backward'' links going from lower
to higher hierarchical levels. We find that the annealing algorithm outperforms
the deterministic one in terms of speed, memory requirement, and the actual
number of removed links. Implications for system biology and directions for
further research are discussed.


  We study the dynamics of gene activities in relatively small size biological
networks (up to a few tens of nodes), e.g. the activities of cell-cycle
proteins during the mitotic cell-cycle progression. Using the framework of
deterministic discrete dynamical models, we characterize the dynamical
modifications in response to structural perturbations in the network
connectivities. In particular, we focus on how perturbations affect the set of
fixed points and sizes of the basins of attraction. Our approach uses two
analytical measures: the basin entropy $H$ and the perturbation size $\Delta$,
a quantity that reflects the distance between the set of fixed points of the
perturbed network to that of the unperturbed network. Applying our approach to
the yeast-cell cycle network introduced by Li \textit{et al.} provides a low
dimensional and informative fingerprint of network behavior under large classes
of perturbations. We identify interactions that are crucial for proper network
function, and also pinpoints functionally redundant network connections.
Selected perturbations exemplify the breadth of dynamical responses in this
cell-cycle model.


  The mechanism of propulsion of host bacteria under the action of actin gel
networks is examined by means of a continuum model of the dynamics of F-actin
concentration. The model includes the elasticity of the network, its attachment
to the host and the polymerization at the interface with it. A formula for the
cruise velocity is derived wherefrom the contributions of elasticity and
polymerization are made explicit.
  The velocity is proportional to the size of the active portion of the gel,
and in turn it elongates its shape. It is shown how a finite latency time is
needed for the onset of motion: the related threshold condition is explicitly
given. It is numerically found that at steady state the velocity has a
pulsating behavior, with sudden decelerations and subsequent recoveries: this
unveils the dominant role of the elastic forces in propulsion.


  We investigate a simple model of microtubule dynamics in which a microtubule
evolves by: (i) attachment of guanosine triphosphate (GTP) to its end at rate
lambda, (ii) GTP converting irreversibly to guanosine diphosphate (GDP) at rate
1, and (iii) detachment of GDP from the end of a microtubule at rate mu. As a
function of these elemental rates, the microtubule can grow steadily or its
length can fluctuate wildly. A master equation approach is developed to
characterize these intriguing features. For mu=0, we find exact expressions for
tubule and GTP cap length distributions, as well as a power-law length
distributions of GTP and GDP islands. For mu=oo, we find the average time
between catastrophes, where the microtubule shrinks to zero length, and extend
this approach to also determine the size distribution of avalanches (sequence
of consecutive GDP detachment events). We obtain the phase diagram for general
rates and verify our predictions by numerical simulations.


  The results of Brownian dynamics simulations of a single DNA molecule in
shear flow are presented taking into account the effect of internal viscosity.
The dissipative mechanism of internal viscosity is proved necessary in the
research of DNA dynamics. A stochastic model is derived on the basis of the
balance equation for forces acting on the chain. The Euler method is applied to
the solution of the model. The extensions of DNA molecules for different
Weissenberg numbers are analyzed. Comparison with the experimental results
available in the literature is carried out to estimate the contribution of the
effect of internal viscosity.


  Phenotype of biological systems needs to be robust against mutation in order
to sustain themselves between generations. On the other hand, phenotype of an
individual also needs to be robust against fluctuations of both internal and
external origins that are encountered during growth and development. Is there a
relationship between these two types of robustness, one during a single
generation and the other during evolution? Could stochasticity in gene
expression have any relevance to the evolution of these robustness? Robustness
can be defined by the sharpness of the distribution of phenotype; the variance
of phenotype distribution due to genetic variation gives a measure of `genetic
robustness' while that of isogenic individuals gives a measure of
`developmental robustness'. Through simulations of a simple stochastic gene
expression network that undergoes mutation and selection, we show that in order
for the network to acquire both types of robustness, the phenotypic variance
induced by mutations must be smaller than that observed in an isogenic
population. As the latter originates from noise in gene expression, this
signifies that the genetic robustness evolves only when the noise strength in
gene expression is larger than some threshold. In such a case, the two
variances decrease throughout the evolutionary time course, indicating increase
in robustness. The results reveal how noise that cells encounter during growth
and development shapes networks' robustness to stochasticity in gene
expression, which in turn shapes networks' robustness to mutation. The
condition for evolution of robustness as well as relationship between genetic
and developmental robustness is derived through the variance of phenotypic
fluctuations, which are measurable experimentally.


  The analysis of secreted proteins represents a challenge for current
proteomics techniques. Proteins are usually secreted at low concentrations in
the culture media, which makes their recovery difficult. In addition, culture
media are rich in salts and other compounds interfering with most proteomics
techniques, which makes selective precipitation of proteins almost mandatory
for a correct subsequent proteomics analysis. Last but not least, the
non-secreted proteins liberated in the culture medium upon lysis of a few dead
cells heavily contaminate the so-called secreted proteins preparations. Several
techniques have been used in the past for concentration of proteins secreted in
culture media. These techniques present several drawbacks, such as
coprecipitation of salts or poor yields at low protein concentrations. Improved
techniques based on carrier-assisted trichloroacetic acid precipitation are
described and discussed in this paper. These techniques have been used to
analyse the secretome of myeloid cells (macrophages, dendritic cells) and
enabled to analyze proteins secreted at concentrations close to 1 ng/ml,
thereby allowing to detect some of the cytokines (TNF, IL-12) secreted by the
myeloid cells upon activation by bacterial products.


  Viral infection requires the binding of receptors on the target cell membrane
to glycoproteins, or ``spikes,'' on the viral membrane. The initial entry is
usually classified as fusogenic or endocytotic. However, binding of viral
spikes to cell surface receptors not only initiates the viral adhesion and the
wrapping process necessary for internalization, but can simultaneously initiate
direct fusion with the cell membrane. Both fusion and internalization have been
observed to be viable pathways for many viruses. We develop a stochastic model
for viral entry that incorporates a competition between receptor mediated
fusion and endocytosis. The relative probabilities of fusion and endocytosis of
a virus particle initially nonspecifically adsorbed on the host cell membrane
are computed as functions of receptor concentration, binding strength, and
number of spikes. We find different parameter regimes where the entry pathway
probabilities can be analytically expressed. Experimental tests of our
mechanistic hypotheses are proposed and discussed.


  In a recent series of ground-breaking experiments, Nakajima et al. [Science
308, 414-415 (2005)] showed that the three cyanobacterial clock proteins KaiA,
KaiB, and KaiC are sufficient in vitro to generate circadian phosphorylation of
KaiC. Here, we present a mathematical model of the Kai system. At its heart is
the assumption that KaiC can exist in two conformational states, one favoring
phosphorylation and the other dephosphorylation. Each individual KaiC hexamer
then has a propensity to be phosphorylated in a cyclic manner. To generate
macroscopic oscillations, however, the phosphorylation cycles of the different
hexamers must be synchronized. We propose a novel synchronisation mechanism
based on differential affinity: KaiA stimulates KaiC phosphorylation, but the
limited supply of KaiA dimers binds preferentially to those KaiC hexamers that
are falling behind in the oscillation. KaiB sequesters KaiA and stabilizes the
dephosphorylating KaiC state. We show that our model can reproduce a wide range
of published data, including the observed insensitivity of the oscillation
period to variations in temperature, and that it makes nontrivial predictions
about the effects of varying the concentrations of the Kai proteins.


  Despite improved rational drug design and a remarkable progress in genomic,
proteomic and high-throughput screening methods, the number of novel,
single-target drugs fell much behind expectations during the past decade.
Multi-target drugs multiply the number of pharmacologically relevant target
molecules by introducing a set of indirect, network-dependent effects. Parallel
with this the low-affinity binding of multi-target drugs eases the constraints
of druggability, and significantly increases the size of the druggable
proteome. These effects tremendously expand the number of potential drug
targets, and will introduce novel classes of multi-target drugs with smaller
side effects and toxicity. Here we review the recent progress in this field,
compare possible network attack strategies, and propose several methods to find
target-sets for multi-target drugs.


  The paper represents three supplements to the source paper, q-bio/0610044
[q-bio.OT], with three new series of harmonic structures of the genetic code,
determined by Gauss arithmetical algorithm; by Table of Minimal Adding, as in
(Rakocevic, 2011a: Table 4; 2011b: Table 4); all structures in relation to
Binary-code tree (Rakocevic, 1998). The determination itself is realized
through atom and nucleon number balancing and nuancing of molekular polarity.
In the first supplement the word is about some additional harmonic structures
in relation to a previous our paper (Rakocevic, 2004); in the second one about
the relation that structures with the polarity of protein amino acids. In the
third supplement we give new ideas about the genetic code by an inclusion of
the notions cipher of the genetic code and the key of that cipher.


  In this paper it is shown that within a Combined Genetic Code Table, realized
through a combination of Watson-Crick Table and Codon Path Cube it exists,
without an exception, a strict distinction between two classes of enzymes
aminoacyl-tRNA synthetases, corresponding two classes of amino acids and
belonging codons. By this, the distinction itself is followed by a strict
balance of atom number within two subclasses of class I as well as two
subclasses of class II of amino acids.


  Myosin motor is the machine, which performs mechanical work in the course of
adenosine triphosphate molecule hydrolysis and myosin head conformations
accompanying this process. For displacement of individual fragments of large
molecule relative to each other to arise and work to be performed, force must
be born inside protein. What kind of interaction generates this force?
  Models based on Huxley 1957 theory ascertain relations between chemical
reactions rate constants and energies of crossbridge conformations.
Nevertheless, understand in the framework of thermodynamics how myosin motor
works in principle is impossible: it is smoothly heated device cyclically
producing mechanical work (second law). Furthermore, in every working cycle
myosin head captures and splits a single molecule. Hence, ordinary dynamic laws
rather than stochastic laws govern this process.
  The simple mechanism of chemomechanical transduction is proposed. The two
products of adenosine triphosphate hydrolysis, adenosine diphosphate and
inorganic phosphate, have charges of the same sense and Coulomb interaction of
these charges produces the force pushing backdoor and rotating converter
domain. The velocity of filaments sliding becomes the principal parameter in
the model and new mechanism of indirect interaction between the cross-bridges
radically different from one suggested by Huxley and Simmons in 1971 appears.
The working stroke duration is inversely proportional to the velocity now.
Therefore Hill equation appears and the parameter values obtained are in
reasonable agreement with experiment.


  We construct an individual-based kinematic model of rolling migratory locust
swarms. The model incorporates social interactions, gravity, wind, and the
effect of the impenetrable boundary formed by the ground. We study the model
using numerical simulations and tools from statistical mechanics, namely the
notion of H-stability. For a free-space swarm (no wind and gravity), as the
number of locusts increases, it approaches a crystalline lattice of fixed
density if it is H-stable, and in contrast becomes ever more dense if it is
catastrophic. Numerical simulations suggest that whether or not a swarm rolls
depends on the statistical mechanical properties of the corresponding
free-space swarm. For a swarm that is H-stable in free space, gravity causes
the group to land and form a crystalline lattice. Wind, in turn, smears the
swarm out along the ground until all individuals are stationary. In contrast,
for a swarm that is catastrophic in free space, gravity causes the group to
land and form a bubble-like shape. In the presence of wind, the swarm migrates
with a rolling motion similar to natural locust swarms. The rolling structure
is similar to that observed by biologists, and includes a takeoff zone, a
landing zone, and a stationary zone where grounded locusts can rest and feed.


  Recent studies on the evolutionary dynamics of the Prisoner's Dilemma game in
scale-free networks have demonstrated that the heterogeneity of the network
interconnections enhances the evolutionary success of cooperation. In this
paper we address the issue of how the characterization of the asymptotic states
of the evolutionary dynamics depends on the initial concentration of
cooperators. We find that the measure and the connectedness properties of the
set of nodes where cooperation reaches fixation is largely independent of
initial conditions, in contrast with the behavior of both the set of nodes
where defection is fixed, and the fluctuating nodes. We also check for the
robustness of these results when varying the degree heterogeneity along a
one-parametric family of networks interpolating between the class of
Erdos-Renyi graphs and the Barabasi-Albert networks.


  We study the genetic behaviour of a population formed by haploid individuals
which reproduce asexually. The genetic information for each individual is
stored along a bit-string (or chromosome) with L bits, where 0-bits represent
the wild-type allele and 1-bits correspond to harmful mutations. Each newborn
inherits this chromosome from its parent with some few random mutations: on
average a fixed number m of bits are flipped. Selection is implemented
according to the number N of 1-bits counted along the individual's chromosome:
the smaller N the higher the probability an individual has to survive a new
time step. Such a population evolves, with births and deaths, and its genetic
distribution becomes stabilised after many enough generations have passed.
  The question we pose concerns the procedure of increasing L. The aim is to
get the same distribution of relative genetic loads N/L among the equilibrated
population, in spite of a larger L. Should we keep the same mutation rate m/L
for different values of L? The answer is yes, which intuitively seems to be
plausible. However, this conclusion is not trivial, according to our
simulational results: the question involves also the population size.


  Sexual reproduction in Nature requires two sexes, which raises the question
why the reproductive scheme did not evolve to have three or more sexes. Here we
construct a constrained optimization model based on the communication theory to
analyze trade-offs among reproductive schemes with arbitrary number of sexes.
More sexes on one hand lead to higher reproductive diversity, but on the other
hand incur greater cost in time and energy for reproductive success. Our model
shows that the two-sexes reproduction scheme maximizes the recombination
entropy-to-cost ratio, and hence is the optimal solution to the problem.


  We study a genetic regulatory network model developed to demonstrate that
genetic robustness can evolve through stabilizing selection for optimal
phenotypes. We report preliminary results on whether such selection could
result in a reorganization of the state space of the system. For the chosen
parameters, the evolution moves the system slightly toward the more ordered
part of the phase diagram. We also find that strong memory effects cause the
Derrida annealed approximation to give erroneous predictions about the model's
phase diagram.


  The population dynamics of mixed-culture and sequential-cultures in
fruit-must fermentation is a very interesting problem from both the theoretical
and technological stand point. By mixed-culture we refer to those fermentations
in which more that one strain/species are present from the beginning of the
process; and by sequential to those in which different microorganisms are added
along the process. These kinds of fermentations play a key role in industry of
fermented beverages. The mathematical models of such processes should represent
the dynamics of multiple species growing in a mixture of substrates as the
fruit must composition includes an important proportion of hexoses and
pentoses. Since the flavor is a basic aspect of beverages, and in general of
all fermented foods, it is fundamental to study the population dynamics and its
impact in the organoleptic properties, through the influence in the volatile
compound profile. The goal of this paper is to present a panorama of the
investigations of these fermentative ecosystems from a biological,
mathematical, and technological stand point.


  The network paradigm is increasingly used to describe the topology and
dynamics of complex systems. Here we review the results of the topological
analysis of protein structures as molecular networks describing their
small-world character, and the role of hubs and central network elements in
governing enzyme activity, allosteric regulation, protein motor function,
signal transduction and protein stability. We summarize available data how
central network elements are enriched in active centers and ligand binding
sites directing the dynamics of the entire protein. We assess the feasibility
of conformational and energy networks to simplify the vast complexity of rugged
energy landscapes and to predict protein folding and dynamics. Finally, we
suggest that modular analysis, novel centrality measures, hierarchical
representation of networks and the analysis of network dynamics will soon lead
to an expansion of this field.


  We modify the classical Shimada-Yamakawa theory of DNA looping by
generalizing the form of boundary constraints. This generalization is important
in the context of DNA cyclization experiments since it mimics the reduced local
rigidity of the "nicked" DNA loop. Our results indicate that the non-trivial
boundary constraints may be responsible for the existing dramatic discrepancy
between various DNA cyclization experiments. The developed effective
Hamiltonian method may be further extended to a broad class of DNA looping
problems.


  We present a simple model which describes elastic response of single-stranded
DNA (ssDNA) to stretching, including the regime of very high force (up to 1000
pN). ssDNA is modelled as a discreet persistent chain, whose ground state is a
zigzag rather than a straight line configuration. This mimics the underlying
molecular architecture and helps to explain the experimentally observed
staturation of the stretching curve at very high force.


  We present a simplified model of the dynamics of translocation of RNA through
a nanopore which only allows the passage of unbound nucleotides. In particular,
we consider the disorder averaged translocation dynamics of random,
two-component, single-stranded nucleotides, by reducing the dynamics to the
motion of a random walker on a one-dimensional free energy landscape of
translocation. These translocation landscapes are calculated from the folds of
the RNA sequences and the voltage bias applied across the nanopore. We compute
these landscapes for 1500 randomly drawn two-letter sequences of length 4000.
Simulations of the dynamics on these landscapes display anomalous
characteristics, similar to random forcing energy landscapes, where the
translocation process proceeds slower than linearly in time for sufficiently
small voltage biases across the nanopore, but moves linearly in time at large
voltage biases. We argue that our simplified model provides an upper bound to
the more realistic translocation dynamics, and thus we expect that all RNA
translocation models will exhibit anomalous regimes.


  We analyze three simple genetic circuits which involve transcriptional
regulation and feedback: the autorepressor, the switch and the repressilator,
that consist of one, two and three genes, respectively. Such systems are
commonly simulated using rate equations, that account for the concentrations of
the mRNAs and proteins produced by these genes. Rate equations are suitable
when the concentrations of the relevant molecules in a cell are large and
fluctuations are negligible. However, when some of the proteins in the circuit
appear in low copy numbers, fluctuations become important and the rate
equations fail. In this case stochastic methods, such as direct numerical
integration of the master equation or Monte Carlo simulations are required.
Here we present deterministic and stochastic simulations of the autorepressor,
the switch and the repressilator. We show that fluctuations give rise to
quantitative and qualitative changes in the dynamics of these systems. In
particular, we demonstrate a fluctuations-induced bistability in a variant of
the genetic switch and and noisy oscillations obtained in the repressilator
circuit.


  Computer simulations prove that we should spend more money on red wine than
on medication. We optimise the problem of how to distribute a fixed amount of
money between medication and food spendings, using the Penna ageing model.


  We investigate the dynamics of an idealized model of microtubule growth that
evolves by: (i) attachment of guanosine triphosphate (GTP) at rate lambda, (ii)
conversion of GTP to guanosine diphosphate (GDP) at rate 1, and (iii)
detachment of GDP at rate mu. As a function of these rates, a microtubule can
grow steadily or its length can fluctuate wildly. For mu=0, we find the exact
tubule and GTP cap length distributions, and power-law length distributions of
GTP and GDP islands. For mu=infinity, we argue that the time between
catastrophes, where the microtubule shrinks to zero length, scales as
exp(lambda). We also find the phase boundary between a growing and shrinking
microtubule.


  The role of the selection pressure and mutation amplitude on the behavior of
a single-species population evolving on a two-dimensional lattice, in a
periodically changing environment, is studied both analytically and
numerically. The mean-field level of description allows to highlight the
delicate interplay between the different time-scale processes in the resulting
complex dynamics of the system. We clarify the influence of the amplitude and
period of the environmental changes on the critical value of the selection
pressure corresponding to a phase-transition "extinct-alive" of the population.
However, the intrinsic stochasticity and the dynamically-built in correlations
among the individuals, as well as the role of the mutation-induced variety in
population's evolution are not appropriately accounted for. A more refined
level of description, which is an individual-based one, has to be considered.
The inherent fluctuations do not destroy the phase transition "extinct-alive",
and the mutation amplitude is strongly influencing the value of the critical
selection pressure. The phase diagram in the plane of the population's
parameters -- selection and mutation is discussed as a function of the
environmental variation characteristics. The differences between a smooth
variation of the environment and an abrupt, catastrophic change are also
addressesd.


  We introduce new methods for phylogenetic tree quartet construction by using
machine learning to optimize the power of phylogenetic invariants. Phylogenetic
invariants are polynomials in the joint probabilities which vanish under a
model of evolution on a phylogenetic tree. We give algorithms for selecting a
good set of invariants and for learning a metric on this set of invariants
which optimally distinguishes the different models. Our learning algorithms
involve linear and semidefinite programming on data simulated over a wide range
of parameters. We provide extensive tests of the learned metrics on simulated
data from phylogenetic trees with four leaves under the Jukes-Cantor and Kimura
3-parameter models of DNA evolution. Our method greatly improves on other uses
of invariants and is competitive with or better than neighbor-joining. In
particular, we obtain metrics trained on trees with short internal branches
which perform much better than neighbor joining on this region of parameter
space.


  Immunity is hypothesized to share limited resources with other physiological
functions and may mediate life history trade-offs, for example between
reproduction and survival. However, vertebrate immune defense is a complex
system that consists of three components. To date, no study has assessed all of
these components for the same animal model and within a given situation.
Previous studies have determined that the acquired immunity of common eiders
(Somateria mollissima) is suppressed during incubation. The present paper aims
to assess the innate immune response in fasting eiders in relation to their
initial body condition. Innate immunity was assessed by measuring plasma nitric
oxide (NO) levels, prior to and after injection of lipopolysaccharides (LPS), a
method which is easily applicable to many wild animals. Body condition index
and corticosterone levels were subsequently determined as indicators of body
condition and stress level prior to LPS injection. The innate immune response
in eiders did not vary significantly throughout the incubation period. The
innate immune response of eiders did not vary significantly in relation to
their initial body condition but decreased significantly when corticosterone
levels increased. However, NO levels after LPS injection were significantly and
positively related to initial body condition, while there was a significant
negative relationship with plasma corticosterone levels. Our study suggests
that female eiders preserve an effective innate immune response during
incubation and this response might be partially determined by the initial body
condition.


  Methods to effectively detect multi-locus genetic association are becoming
increasingly relevant in the genetic dissection of complex trait in humans.
Current approaches typically consider a limited number of hypotheses, most of
which are related to the effect of a single locus or of a relatively small
number of neighboring loci on a chromosomal region. We have developed a novel
method that is specifically designed to detect genetic association involving
multiple disease-susceptibility loci, possibly on different chromosomes. Our
approach relies on the efficient discovery of patterns comprising spatially
unrestricted polymorphic markers and on the use of appropriate test statistics
to evaluate pattern-trait association. Power calculations using multi-locus
disease models demonstrate significant gain of power by using this method in
detecting multi-locus genetic association when compared to a standard single
marker analysis method. When analyzing a Schizophrenia dataset, we confirmed a
previously identified gene-gene interaction. In addition, a less conspicuous
association involving different markers on the same two genes was also
identified, implicating genetic heterogeneity.


  The FKPP equation with a variable growth rate and advection by an
incompressible velocity field is considered as a model for plankton dispersed
by ocean currents. If the average growth rate is negative then the model has a
survival-extinction transition; the location of this transition in the
parameter space is constrained using variational arguments and delimited by
simulations. The statistical steady state reached when the system is in the
survival region of parameter space is characterized by integral constraints and
upper and lower bounds on the biomass and productivity that follow from
variational arguments and direct inequalities. In the limit of
zero-decorrelation time the velocity field is shown to act as Fickian diffusion
with an eddy diffusivity much larger than the molecular diffusivity and this
allows a one-dimensional model to predict the biomass, productivity and
extinction transitions. All results are illustrated with a simple growth and
stirring model.


  Spatially resolved genetic data is increasingly used to reconstruct the
migrational history of species. To assist such inference, we study, by means of
simulations and analytical methods, the dynamics of neutral gene frequencies in
a population undergoing a continual range expansion in one dimension. During
such a colonization period, lineages can fix at the wave front by means of a
``surfing'' mechanism [Edmonds C.A., Lillie A.S. & Cavalli-Sforza L.L. (2004)
Proc Natl Acad Sci USA 101: 975-979]. We quantify this phenomenon in terms of
(i) the spatial distribution of lineages that reach fixation and, closely
related, (ii) the continual loss of genetic diversity (heterozygosity) at the
wave front, characterizing the approach to fixation. Our simulations show that
an effective population size can be assigned to the wave that controls the
(observable) gradient in heterozygosity left behind the colonization process.
This effective population size is markedly higher in pushed waves than in
pulled waves, and increases only sub-linearly with deme size. To explain these
and other findings, we develop a versatile analytical approach, based on the
physics of reaction-diffusion systems, that yields simple predictions for any
deterministic population dynamics.


  We present a spatial, individual-based predator-prey model in which dispersal
is dependent on the local community. We determine species suitability to the
biotic conditions of their local environment through a time and space varying
fitness measure. Dispersal of individuals to nearby communities occurs whenever
their fitness falls below a predefined tolerance threshold. The spatiotemporal
dynamics of the model is described in terms of this threshold. We compare this
dynamics with the one obtained through density-independent dispersal and find
marked differences. In the community-driven scenario, the spatial correlations
in the population density do not vary in a linear fashion as we increase the
tolerance threshold. Instead we find the system to cross different dynamical
regimes as the threshold is raised. Spatial patterns evolve from disordered, to
scale-free complex patterns, to finally becoming well-organized domains. This
model therefore predicts that natural populations, the dispersal strategies of
which are likely to be influenced by their local environment, might be subject
to complex spatiotemporal dynamics.


  In this Supporting Information, we provide background information on our
model of the in vitro Kai system and the calculations that we have performed.
We will closely follow the outline of the main text.


  In this paper, we study a two-lane totally asymmetric simple exclusion
process (TASEP) coupled with random attachment and detachment of particles
(Langmuir kinetics) in both lanes under open boundary conditions. Our model can
describe the directed motion of molecular motors, attachment and detachment of
motors, and free inter-lane transition of motors between filaments. In this
paper, we focus on some finite-size effects of the system because normally the
sizes of most real systems are finite and small (e.g., size $\leq 10,000$). A
special finite-size effect of the two-lane system has been observed, which is
that the density wall moves left first and then move towards the right with the
increase of the lane-changing rate. We called it the jumping effect. We find
that increasing attachment and detachment rates will weaken the jumping effect.
We also confirmed that when the size of the two-lane system is large enough,
the jumping effect disappears, and the two-lane system has a similar density
profile to a single-lane TASEP coupled with Langmuir kinetics. Increasing
lane-changing rates has little effect on density and current after the density
reaches maximum. Also, lane-changing rate has no effect on density profiles of
a two-lane TASEP coupled with Langmuir kinetics at a large
attachment/detachment rate and/or a large system size. Mean-field approximation
is presented and it agrees with our Monte Carlo simulations.


  A recently developed treatment of stochastic processes leads to the
construction of a potential landscape for the dynamical evolution of complex
systems. Since the existence of a potential function in generic settings has
been frequently questioned in literature,herewe study several related
theoretical issues that lie at core of the construction. We showthat the novel
treatment,via a transformation,is closely related to the symplectic structure
that is central in many branches of theoretical physics. Using this insight, we
demonstrate an invariant under the transformation. We further explicitly
demonstrate, in one-dimensional case, the contradistinction among the new
treatment to those of Ito and Stratonovich, as well as others.Our results
strongly suggest that the method from statistical physics can be useful in
studying stochastic, complex systems in general.


  The death of a biological population is an extreme event which we investigate
here for a host-parasitoid system. Our simulations using the Penna ageing model
show how biological evolution can ``teach'' the parasitoids to avoid extinction
by waiting for the right age of the host. We also show the dependence of
extinction time on the population size.


  The development of new techniques to quantitatively measure gene expression
in cells has shed light on a number of systems that display oscillations in
protein concentration. Here we review the different mechanisms which can
produce oscillations in gene expression or protein concentration, using a
framework of simple mathematical models. We focus on three eukaryotic genetic
regulatory networks which show "ultradian" oscillations, with time period of
the order of hours, and involve, respectively, proteins important for
development (Hes1), apoptosis (p53) and immune response (NFkB). We argue that
underlying all three is a common design consisting of a negative feedback loop
with time delay which is responsible for the oscillatory behaviour.


  The shape of semiflexible polymer rings is studied over their whole range of
flexibility. Investigating the joint distribution of asphericity and nature of
asphericity as well as their respective averages we find two distinct shape
regimes depending on the flexibility of the polymer. For small perimeter to
persistence length the fluctuating rings exhibit only planar, elliptical
configurations. At higher flexibilities three dimensional, crumpled structures
arise. Analytic calculations for tight polymer rings confirm an elliptical
shape in the stiff regime.


  In this paper a quantum mechanical description of the assembly/disassembly
process for microtubules is proposed. We introduce creation and annihilation
operators that raise or lower the microtubule length by a tubulin layer.
Following that, the Hamiltonian and corresponding equations of motion are
derived that describe the dynamics of microtubules. These Heisenberg-type
equations are then transformed to semi-classical equations using the method of
coherent structures. The latter equations are very similar to the
phenomenological equations that describe dynamic instability of microtubules in
a tubulin solution.


  Some biologists accept Wright's adaptive landscape idea, believing it is one
of most profound concepts in evolutionary dynamics. Some wouldn't, believing
that "the idea that there is such a quantity remains one of the most widely
held popular misconceptions about evolution." The two groups usually have very
limited communication with each other. Sometimes such isolation can be good,
because it protects budding ideas in a harsh environment. Thanks to the
theoretical and experimental progress during past few years, time may have
arrived to consider both sides seriously. Present letter is a small attempt to
bridge this gap.


  We attempted to test biogeographic hypotheses proposed for the evolution of
tropical faunas using mitochondrial DNA sequences of selected African mammalian
taxa (Insectivora, Rodentia and Primates). 1/ we constructed a molecular
phylogeny of taxa in order to ascertain their monophyly and calibrate a
molecular clock; 2/ we analysed and compared the phylogeographic patterns of
five forest-dwelling taxa; 3/ we evaluated the evolutionary processes
potentially involved in the speciation of cercopithecine Primates. Phylogenetic
results confirm that gene history is not necessarily the same as organism
history. Phylogeographic analyses reveal distinct patterns for each model
species, suggesting differences in initial distributions or different responses
to the same events. They indicate a role of Plio-Pleistocene vicariance events
in the intra-specific diversification of small mammals. Thus, genetic
divergence would be much older than the last glacial cycles. In cercopithecine
Primates, speciation would have been predominantly allopatric and driven by
Miocene and Pliocene vicariance events. Altogether, results give support to the
refuge hypothesis, without excluding the riverine barrier nor the
paleogrographic ones. They emphasize the role of paleo-ecological changes in
generating diversity and of the main riverine barriers in shaping the present
distribution of that diversity.


  The similarity between protein sequences is a directly and easly computed
quantity from which to deduce information about their evolutionary distance and
to detect homologous proteins. The SIMAP database -- Similarity Matrix of
Proteins -- provides a pre-computed similarity matrix covering the similarity
space formed by about all publicly available amino acid sequences from public
databases and completely sequenced genomes. From SIMAP we construct the protein
homology network, where the proteins are the nodes and the links represent
homology relationships. With more than 5 million nodes and about 70 10^9 edges
it is the greatest protein homology network ever been builded. We describe the
basic features and we perform a global statistical analysis of the network.
Starting from the Smith-Waterman similarity score, we define for each edge a
weight w to measure the similarity distance between two nodes. Keeping only
edges with a weigth greater than a minimal w_min, and by varying w_min we build
a family of networks with different degree of similarity. We investigate the
distribution of connected components (clusters) of the networks at different
w_min and in particular we find a behaviour similar to a phase transition
guided by the formation of a giant component. Moreover we study selected
sequence features and protein domains of protein pairs that connect different
clusters in the networks at different level of similarity. We observed
specific, non-random distributions of the protein features and domains for
proteins connecting clusters at certain weight intervals.


  We analyse a simple discrete-time stochastic process for the theoretical
modeling of the evolution of protein lengths. At every step of the process a
new protein is produced as a modification of one of the proteins already
existing and its length is assumed to be random variable which depends only on
the length of the originating protein. Thus a Random Recursive Trees (RRT) is
produced over the natural integers. If (quasi) scale invariance is assumed, the
length distribution in a single history tends to a lognormal form with a
specific signature of the deviations from exact gaussianity. Comparison with
the very large SIMAP protein database shows good agreement.


  I point out the similarity between the microtubule experiment reported by
Priel et al [Biophys. J. 90, 4639 (2006)] and the ZnO nanowire experiment of
Wang et al [Nanolett. 6, 2768 (2006)]. It is quite possible that MTs are
similar to a piezoelectric field effect transistor for which the role of the
control gate electrode is played by the piezo-induced electric field across the
width of the MT walls and their elastic bending features


  It is the goal of systems biology to understand the behavior of the whole in
terms of the knowledge of the parts. This is hard to achieve in many cases due
to the difficulty of characterizing the many constituents and their complex web
of interactions involved in a biological system. The lac promoter of E. coli
offers a possibility of confronting system-leve properties of transcriptional
regulation with the known biochemistry of the molecular constituents and their
mutual interactions. Such confrontations can reveal previously unknown
constituents and interactions, as well as offering new insight into how the
components work together as a whole. Here we study the combinatorial control of
the lac promoter by the regulators LacR and CRP. A previous in vivo study
[Setty et al., PNAS 100: 7702-7 (2003)] found gross disagreement between the
observed promoter activites and the expected behavior based on the known
molecular mechanisms. We repeated the study by identifying and removing several
extraneous factors which significantly modulated the expression of the lac
promoter. Through quantitative, systematic characterization of promoter
activity for a number of key mutants and guided by the thermodynamic model of
transcriptional gene regulation, we are able to account for the combinatorial
control of the lac promoter quantitatively, in terms of a cooperative
interaction between CRP and LacR-mediated DNA looping. Specifically, our
analysis indicates that the sensitivity of the inducer response results from
LacR-mediated DNA looping, which is significantly enhanced by CRP.


  We propose a game-theoretic dynamics of a population of replicating
individuals. It consists of two parts: the standard replicator one and a
migration between two different habitats. We consider symmetric two-player
games with two evolutionarily stable strategies: the efficient one in which the
population is in a state with a maximal payoff and the risk-dominant one where
players are averse to risk. We show that for a large range of parameters of our
dynamics, even if the initial conditions in both habitats are in the basin of
attraction of the risk-dominant equilibrium (with respect to the standard
replication dynamics without migration), in the long run most individuals play
the efficient strategy.


  In this paper we consider the evolution of the HIV infection under a highly
effective treatment based on a combination of RTI and/or PI drugs. This is
usually modelled with a set of ODEs where the drug efficacy is a constant. In
this paper we consider the treatment efficacy as a dynamic variable which
evolves during the treatment. To model the dynamics of the drug efficacy we use
an evolutive type assumption. The resulting set of ODEs are able to reproduce
the typical pattern of the illness evolution under treatment: a period of
remission, which can last several years, followed by a progressive
recrudescence and a viral rebound.


  Elucidating the architecture and dynamics of large scale genetic regulatory
networks of cells is an important goal in systems biology. We study the system
level dynamical properties of the genetic network of Escherichia coli that
regulates its metabolism, and show how its design leads to biologically useful
cellular properties. Our study uses the database (Covert et al., Nature 2004)
containing 583 genes and 96 external metabolites which describes not only the
network connections but also the boolean rule at each gene node that controls
the switching on or off of the gene as a function of its inputs. We have
studied how the attractors of the boolean dynamical system constructed from
this database depend on the initial condition of the genes and on various
environmental conditions corresponding to buffered minimal media. We find that
the system exhibits homeostasis in that its attractors, that turn out to be
fixed points or low period cycles, are highly insensitive to initial conditions
or perturbations of gene configurations for any given fixed environment. At the
same time the attractors show a wide variation when external media are varied
implying that the system mounts a highly flexible response to changed
environmental conditions. The regulatory dynamics acts to enhance the cellular
growth rate under changed media. Our study shows that the reconstructed genetic
network regulating metabolism in {\it E. coli} is hierarchical, modular, and
largely acyclic, with environmental variables controlling the root of the
hierarchy. This architecture makes the cell highly robust to perturbations of
gene configurations as well as highly responsive to environmental changes. The
twin properties of homeostasis and response flexibility are achieved by this
dynamical system even though it is not close to the edge of chaos.


  The capacity of proteins to interact specifically with one another underlies
our conceptual understanding of how living systems function. Systems-level
study of specificity in protein-protein interactions is complicated by the fact
that the cellular environment is crowded and heterogeneous; interaction pairs
may exist at low relative concentrations and thus be presented with many more
opportunities for promiscuous interactions compared to specific interaction
possibilities. Here we address these questions using a simple computational
model that includes specifically designed interacting model proteins immersed
in a mixture containing hundreds of different unrelated ones; all of them
undergo simulated diffusion and interaction. We find that specific complexes
are quite robust to interference from promiscuous interaction partners, only in
the range of temperatures Tdesign>T>Trand. At T>Tdesign specific complexes
become unstable, while at T<Trand formation of specific complexes is suppressed
by promiscuous interactions. Specific interactions can form only if
Tdesign>Trand. This condition requires an energy gap between binding energy in
a specific complex and set of binding energies between randomly associating
proteins, providing a general physical constraint on evolutionary selection or
design of specific interacting protein interfaces. This work has implications
for our understanding of how the protein repertoire functions and evolves
within the context of cellular systems.


  Recently there has been growing interest in the use of Maximum Relative
Entropy (MaxREnt) as a tool for statistical inference in ecology. In contrast,
here we propose MaxREnt as a tool for applying statistical mechanics to
ecology. We use MaxREnt to explain and predict species abundance patterns in
ecological communities in terms of the most probable behaviour under given
environmental constraints, in the same way that statistical mechanics explains
and predicts the behaviour of thermodynamic systems. We show that MaxREnt
unifies a number of different ecological patterns: (i) at relatively local
scales a unimodal biodiversity-productivity relationship is predicted in good
agreement with published data on grassland communities, (ii) the predicted
relative frequency of rare vs. abundant species is very similar to the
empirical lognormal distribution, (iii) both neutral and non-neutral species
abundance patterns are explained, (iv) on larger scales a monotonic
biodiversity-productivity relationship is predicted in agreement with the
species-energy law, (v) energetic equivalence and power-law self-thinning
behaviour are predicted in resource-rich communities. We identify mathematical
similarities between these ecological patterns and the behaviour of
thermodynamic systems, and conclude that the explanation of ecological patterns
is not unique to ecology but rather reflects the generic statistical behaviour
of complex systems with many degrees of freedom under very general types of
environmental constraints.


  Many socio-economic and biological processes can be modeled as systems of
interacting individuals. The behaviour of such systems can be often described
within game-theoretic models. In these lecture notes, we introduce fundamental
concepts of evolutionary game theory and review basic properties of
deterministic replicator dynamics and stochastic dynamics of finite
populations. We discuss stability of equilibria in deterministic dynamics with
migration, time-delay, and in stochastic dynamics of well-mixed populations and
spatial games with local interactions. We analyze the dependence of the
long-run behaviour of a population on various parameters such as the time
delay, the noise level, and the size of the population.


  Transcription regulation is largely governed by the profile and the dynamics
of transcription factors' binding to DNA. Stochastic effects are intrinsic to
this dynamics and the binding to functional sites must be controled with a
certain specificity for living organisms to be able to elicit specific cellular
responses. Specificity stems here from the interplay between binding affinity
and cellular abundancy of transcription factor proteins and the binding of such
proteins to DNA is thus controlled by their chemical potential.
  We combine large-scale protein abundance data in the budding yeast with
binding affinities for all transcription factors with known DNA binding site
sequences to assess the behavior of their chemical potentials. A sizable
fraction of transcription factors is apparently bound non-specifically to DNA
and the observed abundances are marginally sufficient to ensure high
occupations of the functional sites. We argue that a biological cause of this
feature is related to its noise-filtering consequences: abundances below
physiological levels do not yield significant binding of functional targets and
mis-expressions of regulated genes are thus tamed.


  Immunity is believed to share limited resources with other physiological
functions and this may partly account for the fitness costs of reproduction.
Previous studies have shown that the acquired immunity of female common eiders
(Somateria mollissima) is suppressed during the incubation fast. To save
energy, triiodothyronine (T3) is adaptively decreased during fasting in most
bird species, despite T3 levels are maintained throughout incubation in female
eiders. However, the relationship between thyroid hormones and the immune
system is not fully understood. The current study aimed to determine the
endocrine mechanisms that underlie immunosuppression in incubating female
eiders. ...


  Many cases of non-standard genetic codes are known in mitochondrial genomes.
We carry out analysis of phylogeny and codon usage of organisms for which the
complete mitochondrial genome is available, and we determine the most likely
mechanism for codon reassignment in each case. Reassignment events can be
classified according to the gain-loss framework. The gain represents the
appearance of a new tRNA for the reassigned codon or the change of an existing
tRNA such that it gains the ability to pair with the codon. The loss represents
the deletion of a tRNA or the change in a tRNA so that it no longer translates
the codon. One possible mechanism is Codon Disappearance, where the codon
disappears from the genome prior to the gain and loss events. In the
alternative mechanisms the codon does not disappear. In the Unassigned Codon
mechanism, the loss occurs first, whereas in the Ambiguous Intermediate
mechanism, the gain occurs first. Codon usage analysis gives clear evidence of
cases where the codon disappeared at the point of the reassignment and also
cases where it did not disappear. Codon disappearance is the probable
explanation for stop to sense reassignments and a small number of reassignments
of sense codons. However, the majority of sense to sense reassignments cannot
be explained by codon disappearance. In the latter cases, by analysis of the
presence or absence of tRNAs in the genome and of the changes in tRNA
sequences, it is sometimes possible to distinguish between the Unassigned Codon
and Ambiguous Intermediate mechanisms. We emphasize that not all reassignments
follow the same scenario and that it is necessary to consider the details of
each case carefully.


  Evolutionary relationships between species are usually represented in
phylogenies, i.e. evolutionary trees, which are a type of networks. The
terminal nodes of these trees represent species, which are made of individuals
and populations among which gene flow occurs. This flow can also be represented
as a network. In this paper we briefly show some properties of these complex
networks of evolutionary and ecological relationships. First, we characterize
large scale evolutionary relationships in the Tree of Life by a degree
distribution. Second, we represent genetic relationships between individuals of
a Mediterranean marine plant, Posidonia oceanica, in terms of a Minimum
Spanning Tree. Finally, relationships among plant shoots inside populations are
represented as networks of genetic similarity.


  Replication of DNA and synthesis of proteins are studied from the view-point
of quantum database search. Identification of a base-pairing with a quantum
query gives a natural (and first ever) explanation of why living organisms have
4 nucleotide bases and 20 amino acids. It is amazing that these numbers arise
as solutions to an optimisation problem. Components of the DNA structure which
implement Grover's algorithm are identified, and a physical scenario is
presented for the execution of the quantum algorithm. It is proposed that
enzymes play a crucial role in maintaining quantum coherence of the process.
Experimental tests that can verify this scenario are pointed out.


  The biological hierarchy and the differences between living and non-living
matter are considered from the standpoint of quantum mechanics.


  Does quantum dynamics play a role in DNA replication? What type of tests
would reveal that? Some statistical checks that distinguish classical and
quantum dynamics in DNA replication are proposed.


  Information is often encoded as an aperiodic chain of building blocks. Modern
digital computers use bits as the building blocks, but in general the choice of
building blocks depends on the nature of the information to be encoded. What
are the optimal building blocks to encode structural information? This can be
analysed by substituting the operations of addition and multiplication of
conventional arithmetic with translation and rotation. It is argued that at the
molecular level, the best component for encoding discretised structural
information is carbon. Living organisms discovered this billions of years ago,
and used carbon as the back-bone for constructing proteins that function
according to their structure. Structural analysis of polypeptide chains shows
that an efficient and versatile structural language of 20 building blocks is
needed to implement all the tasks carried out by proteins. Properties of amino
acids indicate that the present triplet genetic code was preceded by a more
primitive one, coding for 10 amino acids using two nucleotide bases.


  Living organisms are not just random collections of organic molecules. There
is continuous information processing going on in the apparent bouncing around
of molecules of life. Optimisation criteria in this information processing can
be searched for using the laws of physics. Quantum dynamics can explain why
living organisms have 4 nucleotide bases and 20 amino acids, as optimal
solutions of the molecular assembly process. Experiments should be able to tell
whether evolution indeed took advantage of quantum dynamics or not.


  Possible explanation of the 64/20 redundancy of the triplet genetic code
based on the assumption of quantum nature of genetic information is proposed.


  It is a fascinating subject to explore how well we can understand the
processes of life on the basis of fundamental laws of physics. It is emphasised
that viewing biological processes as manipulation of information extracts their
essential features. This information processing can be analysed using
well-known methods of computer science. The lowest level of biological
information processing, involving DNA and proteins, is the easiest one to link
to physical properties. Physical underpinnings of the genetic information that
could have led to the universal language of 4 nucleotide bases and 20 amino
acids are pointed out. Generalisations of Boolean logic, especially features of
quantum dynamics, play a crucial role.


  Competition for available resources is natural amongst coexisting species,
and the fittest contenders dominate over the rest in evolution. The dynamics of
this selection is studied using a simple linear model. It has similarities to
features of quantum computation, in particular conservation laws leading to
destructive interference. Compared to an altruistic scenario, competition
introduces instability and eliminates the weaker species in a finite time.


  We describe the evolution of macromolecules as an information transmission
process and apply tools from Shannon information theory to it. This allows us
to isolate three independent, competing selective pressures that we term
compression, transmission, and neutrality selection. The first two affect
genome length: the pressure to conserve resources by compressing the code, and
the pressure to acquire additional information that improves the channel,
increasing the rate of information transmission into each offspring. Noisy
transmission channels (replication with mutations) gives rise to a third
pressure that acts on the actual encoding of information; it maximizes the
fraction of mutations that are neutral with respect to the phenotype. This
neutrality selection has important implications for the evolution of
evolvability. We demonstrate each selective pressure in experiments with
digital organisms.


  We examine the classical contents of quantum games. It is shown that a
quantum strategy can be interpreted as a classical strategies with effective
density-dependent game matrices composed of transposed matrix elements. In
particular, successful quantum strategies in dilemma games are interpreted in
terms of a symmetrized game matrix that corresponds to an altruistic game plan.


  A quantum spin chain is identified by the labels of a vector state of a
Kashiwara crystal basis. The intensity of the one-spin flip is assumed to
depend from the variation of the labels. The rank ordered plot of the
numerically computed, averaged in time, transition probabilities is nicely
fitted by a Yule distribution, which is the observed distribution of the ranked
short oligonucleotides frequency in DNA.


  We show that the quantum decoherence of Forster resonant energy transfer
between two optically active molecules can be described by a spin-boson model.
This allows us to give quantitative criteria, in terms of experimentally
measurable system parameters, that are necessary for coherent Bloch
oscillations of excitons between the chromophores. Experimental tests of our
results should be possible with Flourescent Resonant Energy Transfer (FRET)
spectroscopy. Although we focus on the case of protein-pigment complexes our
results are also relevant to quantum dots and organic molecules in a dielectric
medium.


  We formulate and study a quantum field theory of a microtubule, a basic
element of living cells. Following the quantum theory of consciousness by
Hameroff and Penrose, we let the system to reduce to one of the classical
states without measurement if certain conditions are
satisfied(self-reductions), and calculate the self-reduction time $\tau_N$ (the
mean interval between two successive self-reductions) of a cluster consisting
of more than $N$ neighboring tubulins (basic units composing a microtubule).
$\tau_N$ is interpreted there as an instance of the stream of consciousness. We
analyze the dependence of $\tau_N$ upon $N$ and the initial conditions, etc.
For relatively large electron hopping amplitude, $\tau_N$ obeys a power law
$\tau_N \sim N^b$, which can be explained by the percolation theory. For
sufficiently small values of the electron hopping amplitude, $\tau_N$ obeys an
exponential law, $\tau_N \sim \exp(c' N)$. By using this law, we estimate the
condition for $\tau_N $ to take realistic values $\tau_N$
\raisebox{-0.5ex}{$\stackrel{>}{\sim}$} $10^{-1}$ sec as $N$ \raisebox{-0.5ex}
{$\stackrel{>}{\sim}$} 1000.


Amelogenins are a group of protein isoforms produced by alternative splicing or proteolysis from the AMELX gene, on the X chromosome, and also the AMELY gene in males, on the Y chromosome. They are involved in amelogenesis, the development of enamel. Amelogenins are type of extracellular matrix protein, which, together with ameloblastins, enamelins and tuftelins, direct the mineralization of enamel to form a highly organized matrix of rods, interrod crystal and proteins.
Although the precise role of amelogenin(s) in regulating the mineralization process is unknown, it is known that amelogenins are abundant during amelogenesis. Developing human enamel contains about 70% protein, 90% of which are amelogenins.

Function
Amelogenins are believed to be involved in the organizing of enamel rods during tooth development. The latest research indicates that these proteins regulate the initiation and growth of hydroxyapatite crystals during the mineralization of enamel. In addition, amelogenins appear to aid in the development of cementum by directing cementoblasts to the tooth's root surface.

Variants
The amelogenin gene has been most widely studied in humans, where it is a single copy gene, located on the X and Y chromosomes at Xp22.1–Xp22.3 and Yp 11.2 [5]. The amelogenin gene's location on sex chromosomes has implications for variability both between the X chromosome form (AMELX) and the Y chromosome form (AMELY), and between alleles of AMELY among different populations. This is because AMELY exists in the non-recombining region of chromosome Y, effectively isolating it from normal selection pressures. Other sources of amelogenin variation arise from the various isoforms of AMELX obtained from alternative splicing of mRNA transcripts. Specific roles for isoforms have yet to be established. Among other organisms, amelogenin is well conserved among eutherians, and has homologs in monotremes, reptiles and amphibians.

Application in sex determination
Differences between the X chromosome and Y chromosome versions of the amelogenin gene (AMELX and AMELY respectively) enable it to be used in sex determination of unknown human samples. AMELX’s intron 1 contains a 6-base-pair deletion relative to intron 1 of AMELY. This can be detected at low cost using polymerase chain reaction (PCR) of intron 1, followed by gel electrophoresis. Two bands of DNA, at 555 bps and 371 bps, are resolved if both the AMELX and AMELY versions of the gene are present (i.e. the sample is from a male) or one band of DNA, at 555 bps, if the AMELX version only is present (i.e. the sample is from a female).However, because of AMELY variation among individuals and populations, this method of sex determination is not 100% accurate. Mutation in regions of AMELY intron 1 commonly used as primer annealing sites may disable PCR amplification. A 6bp insertion to AMELY intron 1 results in an amplicon identical in length to that of AMELX. In some males AMELY may be deleted entirely. In any of these cases only one band is visualized during gel electrophoresis of PCR products, causing misidentification of the sample as female. The misidentification rate may vary among populations, but in general appears to be low. In one study in Spain, the amelogenin sex determination test using AMELX (977bps) and AMELY (790bps) bands was performed for 1224 individuals of known gender with a 99.84% (1222/1224) accuracy rate. Another study in India, however, found 5 of its 270 men studied (1.85%) possessed an AMELY deletion, terming them "deleted-amelogenin males" (DAMs). In response the authors suggested that while the amelogenin sex test may be accurate in general, other Y chromosome markers such as SRY, STR, or 50f2 can be used for less ambiguous gender identification.

Clinical significance
Mutations in AMELX can cause amelogenesis imperfecta, a disorder of tooth enamel development.

== References ==


Amely was an American rock band from Orlando, Florida, United States, formed in 2008. The band comprised four members; Petie Pizarro (Vocals/Guitar), Brandon Walden (Guitar), Patrick Ridgen (Bass) and Nate Parsell (Drums). The sound of the band was a mix of rock with power pop elements. Having been a band for a short period of time, Amely managed to be signed to an independent record label, for this genre of music, Fearless Records.

History
Amely did not start out as a fully formed band. It started in 2008 with a studio project for the singer and rhythm guitarist Petie Pizarro (24 years old), who was moving forward as a soloist after the breakup of his previous band. He began writing songs that he would afterwards choose from, for the recording of an  EP. Later on, he dubbed "On My Own," at Vision Sound Studios Vision Sound Studios in Jacksonville, Florida with producer Daryl Phennager (Red Jumpsuit Apparatus). Even though Pizarro played most of the instruments on the recording, including guitar, bass, piano and synths, he brought Nate Parsell in to record some complicated drum parts, since they had previously been together in a band. Upon the completion of this EP, in June 2008, Pizarro and Parsell realized they had something worth further exploring musically and agreed to form a band. Bassist, Patrick Ridgen, was soon found and lead guitarist, Brandon Walden, joined them a little later and so Amely was formed. Parsell started as a drummer in the basement of his old house in rural Ohio. He studied percussion privately throughout his high school career, while also playing in marching band, concert band, jazz ensemble, pep band, and pit orchestra for four years, as well as a handful of fruitless "rock bands." He then studied in a music conservatory and earned a degree in recording arts and music business from Full Sail University. Some of his early influences were John Bonham and Charlie Watts along with some more modern influences such as Dave Grohl, Travis Barker, Carter Beauford, and Zac Farro from Paramore.

Touring, unsigned
A blend of rock combined with a dose of power pop influences, Amely  started  making their mark on the music scene, despite only having been formed after Petie Pizarro's project/EP, for  a year. Since playing their first show in March 2009, Amely had completed three national tours, opened for All Time Low, played at the Unsilent Night Festival, with A Day to Remember, Red Jumpsuit Apparatus and Breathe Carolina, and performed on the Kevin Says Stage at the Florida dates of Vans Warped Tour 2009. The band had such a success, that they played in American Samoa in the South Pacific, in front of 5000 people, for Blue Sky's 10 Year Bash. Regarding the geared used in their live shows, the drummer, Parsell, used a four-piece Tama Starclassic Birch kit: 10×12 rack, 16×16 floor, and 18×22 kick—one of the best-sounding kick drums out there, according to him. Also, he's been using Zildjian cymbals, a 24″ K Light ride, a 21″ Brilliant Sweet ride (as a left-side crash), a 20″ A Custom Projection crash, and 14″ A Custom Mastersound hats. In addition to their touring, the band reached nearly a million plays on their MySpace page and was named one of PureVolume's Top 20 Unsigned bands of 2009.

Recording
After touring, and their success, they caught a major recording label's eye. Although they had only been formed as a band for less than a year, in January 2010, Amely was signed by Fearless Records. After being signed, the band headed immediately into the studio with Zack Odem, Kenneth Mount (All Time Low, Mayday Parade), and Mike Green (Paramore, Breathe Carolina) to record their debut  EP Hello World, which was released on May 25, 2010. The EP consisted of eight catchy tracks that introduced the band's melodic, pop rock-influenced sound. At this point, lead singer Pizzaro, compared their music to pop-rock bands such as The All American Rejects and Jimmy Eat World and wanted his band to have mass appeal without any age range, gender, or mindset restrictions. "We want to be a mature pop/rock band with a sound that kids can get into, but adults will also enjoy," Pizarro notes. Apart from their EP, the first holiday-themed compilation album series created by Fearless Records called 'Tis the Season to Be Fearless featuring eight Christmas-themed songs was released on November 22, 2010, including Amely's  cover of the song "It's Christmas Time Again" originally by Backstreet Boys. The band has also covered the famous song "The Climb" by Miley Cyrus. The song was originally in the  soundtrack of the Walt Disney Pictures's movie Hannah Montana and the album was released by Walt Disney Records, but the purpose of the band recording it was to be in the Fearless Records compilation album Punk Goes Pop 3. However, their cover was not on the album and, instead, one their original songs ("Come Back to Me") was included on the bonus CD of Punk Goes Pop 3. Apart from recording, the band went on "Let's Be Animals Tour"  with The Downtown Fiction, He Is We and Cady Groves in spring 2011.

Break-up
On November 1, 2011 Amely announced via webchat that they would be breaking up. In addition, before the official break up, the band announced the release of another EP titled Raleigh Sessions, which was self recorded, produced, and mixed. The entire process from start to finish was created by the band alone in Raleigh, North Carolina. The EP was released on November 21, 2011. However, it was released only in a digital form and can be found on iTunes and Spotify. Despite the sad news for the fans, Pizarro, Parsell and Ridgen also announced that they will be forming a new band called Sovereign Soldiers, and Walden would be moving forward as a solo artist. The band's announcement stated:

This decision did not come lightly and was the effect of a great number of things that seemed beyond our control. It has been a more than amazing journey with our best friends and this destination we've now reached is not an end, it is a new beginning. Although we have been extremely discouraged by the discourse that fell upon us, we have not let it knock us down and we hope that you can stay strong with all of us. Since the start of this year and even before that we have been writing for what we intended to be our first full-length album. As the chance to professionally produce our music appeared to grow further in distance we took it upon ourselves to record the songs that we had poured our hearts into. After combining the diverse and plenty skills of all 4 of us, these home recordings are all that's left of Amely. We believe this is some of our greatest work and we want to be able to share it with you. Though it may come as a parting gift, we hope this collection of songs finds you well. It may be the remains of Amely, but we hope it reminds you, this is not the end.

Band members
Petie Pizarro - Lead Vocalist/Rhythm Guitarist
Brandon Walden – Lead Guitarist
Patrick Ridgen - Bassist
Nate Parsell - Drummer
Jason Seife Guitarist

Discography
The Raleigh Sessions (2011)
Live Under Lights (Live) EP (2011)
'Tis The Season To Be Fearless (2010)
Hello World (2010)
Amely EP (2009)

References
External links
"Amely | Listen and Stream Free Music, Albums, New Releases, Photos, Videos". Myspace.
"Amely". Facebook.com.
"Amely". YouTube.
"Page Title". Nomadnarrative.com.
"Petie Pizarro Music". YouTube.
"TheRecoveryFL". YouTube.
"Fearless Records | Join The Fearless Family". Fearless Records.
"Fearless Records". YouTube.
"SovereignSoldiers". YouTube.


Aniridia is the absence of the iris, a muscular structure that opens and closes the pupil to allow light into the eye. It is also responsible for eye color. Without it, the central eye appears all black. It can be congenital, in which both eyes are usually involved, or caused by a penetrant injury.  Isolated aniridia is a congenital disorder that is not limited to a defect in iris development,  but is a  panocular condition with macular and optic nerve hypoplasia, cataract,  and corneal changes. Vision may be severely compromised and the disorder is frequently associated with some ocular complications: nystagmus, amblyopia, buphthalmos, and cataract. Aniridia in some individuals occurs as part of a syndrome, such as WAGR syndrome (kidney nephroblastoma (Wilms tumour), genitourinary anomalies and intellectual disability), or Gillespie syndrome (cerebellar ataxia).

PAX6
The AN2 region of the short  arm of chromosome 11 (11p13) includes the PAX6 gene (named for its PAired boX status), whose gene product helps regulate a cascade of other genetic processes involved in the development of the eye (as well as other non-ocular structures).  This PAX6 gene is around 95% similar to the pax gene found in zebrafish, a creature whose ancestors diverged from human evolutionary development around 400 million years ago. Thus the PAX6 gene is  highly conserved across evolutionary  lineages.
Defects in the PAX6 gene cause aniridia-like ocular defects in mice (as well as Drosophila). Aniridia is a heterozygous disorder, meaning that only one of the two chromosome 11 copies is affected. When both copies are altered (homozygous condition), the result is a uniformly fatal condition with near complete failure of entire eye formation. In 2001, two cases of  homozygous aniridia patients were reported; the fetuses died prior to birth and had severe brain damage. In mice, homozygous small eye defect (mouse Pax-6) leads to loss of the eyes and nose and the murine fetuses sustain severe brain damage.

Types
Aniridia may be broadly divided  into hereditary and sporadic forms. Hereditary aniridia is usually transmitted in an autosomal dominant manner (each offspring has a 50% chance of being affected), although rare autosomal recessive forms (such as Gillespie syndrome) have also been reported. Sporadic aniridia mutations may affect the WT1 region adjacent to the AN2 aniridia region, causing a kidney cancer called nephroblastoma (Wilms tumor). These patients often also have genitourinary abnormalities and intellectual disability (WAGR syndrome).
Several different mutations may affect the PAX6 gene. Some mutations appear to inhibit gene function more than others, with subsequent variability in the severity of the disease. Thus, some aniridic individuals are only missing a relatively small amount of iris, do not have foveal hypoplasia, and retain relatively normal vision. Presumably, the genetic defect in these individuals causes less "heterozygous insufficiency," meaning they retain enough gene function to yield a milder phenotype.

Online Mendelian Inheritance in Man (OMIM): 106210 AN
Online Mendelian Inheritance in Man (OMIM): 106220 Aniridia and absent patella
Online Mendelian Inheritance in Man (OMIM): 106230 Aniridia, microcornea, and  spontaneously reabsorbed cataract
Online Mendelian Inheritance in Man (OMIM): 206700 Aniridia, cerebellar ataxia, and mental deficiency (Gillespie syndrome)

Mutational analysis
Molecular (DNA) testing for PAX6 gene mutations (by sequencing of the entire coding region and deletion/duplication analysis) is available for isolated aniridia and the Gillespie syndrome. For the WAGR syndrome, high-resolution cytogenetic analysis and fluorescence in situ hybridization (FISH) can be utilized to identify deletions within chromosome band 11p13, where both the PAX6 and WT1 genes are located.

Symptoms
Aniridia can cause many symptoms, such as:

Poor vision (not always present)
More sensitivity to light
Fast, uncontrolled, shaking "to and from" eye movements (nystagmus)
Eyes don't line up (strabismus)

Treatment
In May 2018, the U.S. Food and Drug Administration approved the CustomFlex Artificial Iris, the first synthetic iris for use in adults and children with congenital aniridia or iris defects related to other conditions, such as albinism, traumatic injury, or surgical removal due to ocular melanoma. The artificial iris is a surgically implanted device made of thin, foldable, medical-grade silicone and is custom-sized and colored for each individual patient. The prosthetic iris is held in place by the anatomical structures of the eye or, if needed, by sutures.

See also
WAGR syndrome
Scleral lenses

References
External links

GeneReviews/NCBI/NIH/UW entry on Aniridia
NCBI/Molecular diagnosis of aniridia
OMIM entries on Aniridia
GeneReviews/NIH/NCBI/UW entry on Wilms Tumor Overview


The antagonistic pleiotropy hypothesis was first proposed by George C. Williams in 1957 as an evolutionary explanation for senescence. Pleiotropy is the phenomenon where one gene controls more than one phenotypic trait in an organism. A gene is considered to possess antagonistic pleiotropy if it controls more than one trait, where at least one of these traits is beneficial to the organism's fitness early on in life and at least one is detrimental to the organism's fitness later on due to a decline in the force of natural selection. The theme of G. C. William's idea about antagonistic pleiotropy was that if a gene caused both increased reproduction in early life and aging in later life, then senescence would be adaptive in evolution. For example, one study suggests that since follicular depletion in human females causes both more regular cycles in early life and loss of fertility later in life through menopause, it can be selected for by having its early benefits outweigh its late costs.

As a constraint on perfection
Antagonistic pleiotropy is one of the several reasons evolutionary biologists give for organisms never being able to reach perfection through natural selection. Antagonistically pleiotropic genes are the explanation for fitness trade-offs. This means that genes that are pleiotropic control for some beneficial traits and some detrimental traits; thus, if they happen to persist through natural selection, this will prevent organisms from reaching perfection because if they possess the benefits of the gene, they must also possess the imperfections or faults. An example of this would be female rodents that live in a nest with other females and may end up feeding young that are not theirs due to their intense parental drive. This strong parental drive will be selected for, but the organisms will still make the mistake of feeding young that are not theirs and misallocating their resources.

Benefits and costs
Antagonistic pleiotropy has several negative consequences.  It results in delayed adaptation, an altered path of evolution, and reduced adaptation of other traits. In addition, the overall benefit of alleles is cut down significantly (by about half) by pleiotropy. Still, antagonistic pleiotropy has some evolutionary benefits.  In fact, the conservation of genes is directly related to the pleiotropic character of an organism. This implies that genes that control for multiple traits, even if the traits have different implications for the organism's fitness, have more staying power in an evolutionary context.

Role in sexual selection
It is generally accepted that the evolution of secondary sexual characteristics persists until the relative costs of survival outweigh the benefits of reproductive success. At the level of genes, this means a trade-off between variation and expression of selected traits. Strong, persistent sexual selection should result in decreased genetic variation for these traits. However, higher levels of variation have been reported in sexually-selected traits compared to non-sexually selected traits. This phenomenon is especially clear in lek species, where males' courtship behavior confers no immediate advantage to the female. Female choice presumably depends on correlating male displays (secondary sexual characteristics) with overall genetic quality. If such directional sexual selection depletes variation in males, why would female choice continue to exist? Rowe and Houle answer this question (the lek paradox) using the notion of genetic capture, which couples the sexually-selected traits with the overall condition of the organism. They posit that the genes for secondary sexual characteristics must be pleiotropically linked to condition, a measure of the organism's fitness. In other words, the genetic variation in secondary sexual characteristics is maintained due to variation in the organism's condition.

Role in disease
The survival of many serious genetic disorders in our long evolutionary history has led researchers to reassess the role of antagonistic pleiotropy in disease. If genetic disorders are defined by the existence of deleterious alleles, then natural selection acting over evolutionary time would result in a lower frequency of mutations than are currently observed. In a recent article, Carter and Nguyen identify several genetic disorders, arguing that far from being a rare phenomenon, antagonistic pleiotropy might be a fundamental mechanism for the survival of these non-optimal alleles.
In one of these studies, 99 individuals with Laron syndrome (a rare form of dwarfism) were monitored alongside their non-dwarf kin for a period of ten years. Patients with Laron syndrome possess one of three genotypes for the growth hormone receptor gene (GHR). Most patients have an A->G splice site mutation in position 180 in exon 6. Some others possess a nonsense mutation (R43X), while the rest are heterozygous for the two mutations. Laron syndrome patients experienced a lower incidence of cancer mortality and diabetes compared to their non-dwarf kin. This suggests a role for antagonistic pleiotropy, whereby a deleterious mutation is preserved in a population because it still confers some survival benefit.
Another instance of antagonistic pleiotropy is manifested in Huntington's disease, a rare neurodegenerative disorder characterized by a high number of CAG repeats within the Huntingtin gene. The onset of Huntington's is usually observed post-reproductive age and generally involves involuntary muscle spasms, cognitive difficulties and psychiatric problems. Incidentally, the high number of CAG repeats is associated with increased activity of p53, a tumor suppressing protein that participates in apoptosis. It has been hypothesized that this explains the lower rates of cancer among Huntington's patients. Huntington's disease is also correlated with high fecundity.Additionally, it was found that individuals with a higher pro-inflammatory ratio TNFα/IL-10 had a significantly higher incidence of death due to cardiovascular disease in old age. Yet, it was hypothesized that this genotype was prevalent because higher ratios of TNFα/IL-10 allow individuals to more effectively combat infection during reproductive years.Sickle cell anemia, Beta-thalassemia, and cystic fibrosis are some other examples of the role antagonistic pleiotropy may play in genetic disorders.

Ubiquity
Although there are so many negative effects related to genes that are antagonistically pleiotropic, it is still present among most forms of life. Indeed, pleiotropy is one of the most common traits possessed by genes overall. In addition to that, pleiotropy is under strong stabilizing selection. In one experiment with mice and the morphology of the mandible, 1/5 of the loci had effects of pleiotropy for the entire mandible. One other example was in the Russian biologist Dmitry K. Belyaev's study on the domestication of the fox. In Dmitry K. Belyaev's farm-fox experiment, wild foxes were bred for docile behavior alone. After 40 generations, other physiological changes had surfaced including shortened tails, floppy ears, a white star in the forehead, rolled tails, shorter legs. Since the only thing being selected for was behavior, this leads scientists to believe that these secondary characteristics were controlled by the same gene or genes as docile behavior.

Adaptivity and senescence
An antagonistically pleiotropic gene can be selected for if it has beneficial effects in early life while having its negative effects in later life because genes tend to have larger impacts on fitness in an organism's prime than in their old age. An example of this is testosterone levels in male humans.  Higher levels of this hormone lead to increased fitness in early life, while causing decreased fitness in later life due to a higher risk for prostate cancer. This is an example of antagonistic pleiotropy being an explanation for senescence. Senescence is the act of ageing in individuals; it's the failure over time of the individual's life processes by natural causes. Williams's theory has been the motivation for many of the experimental studies on the reasons for aging in the last 25 years. However, there is more than one theory out there for aging. The competing model to explain senescence is Medawar's "mutation accumulation" hypothesis, saying that "over evolutionary time, late-acting mutations will accumulate at a much faster rate than early-acting mutation. These late-acting mutations will thus lead to declining viability and/or fertility as an organism ages." Medawar's theory is based around the older concept of selection shadow that had been discussed throughout the early 1900s and led to Medawar's theory after discussions with J. B. S. Haldane in the 1940s.

Potential Examples
DNA repair
DNA Damage Theory of Aging
A prominent explanation for aging at the molecular level is the DNA damage theory of aging.  It has been proposed that genetic elements that regulate DNA repair in somatic cells may constitute an important example of age-dependent pleiotropic "genes".  As pointed out by Vijg, genome repair and maintenance is beneficial early in life by swiftly eliminating DNA damage or damaged cells. However, studies of DNA repair in the brain and in muscle indicate that the transition from mitotic cell division to the post-mitotic condition that occurs early in life is accompanied by a reduction in DNA repair.  The reduced expression of DNA repair is presumably part of an evolutionary adaptation for diverting the resources of the cell that were previously used for DNA repair, as well as for replication and cell division, to more essential neuronal and muscular functions.The harmful effect of this genetically controlled reduction in expression is to allow increased accumulation of DNA damage.  Reduced DNA repair causes increased impairment of transcription and progressive loss of cell and tissue function.  However, these harmful effects of DNA damage are cumulative and most severe in chronologically older individuals whose numbers diminish with time (by causes of death that can be independent of senescence).  As a consequence, the beneficial effects of the genetic elements that control the reduction of DNA repair early in life would predominate.  Thus regulatory genetic elements that reduce expression of DNA repair genes in post-mitotic cells appear to be important examples of the postulated pleiotropic "genes" that are beneficial in youth but deleterious at an older age.

Telomere Theory
Another example related to aging is the Telomere theory. Telomere theory proposes that telomeres shorten with repeated cell division which attribute to cell senescence and tissue damage. The end replication problem explains the mechanism behind the inability of DNA polymerase to commence the RNA primer to perform its function in completing the lagging strand due to the shortening of DNA. Telomere shortening is common in somatic cells. However, germ line and stem cells prevent the end replication problem with the help of telomerase. Telomerase elongates the 3’ end that is then formed into a t-loop to prevent the cell from entering the G0 phase and cell senescence.Inflammation and damage to tissue are the underlying problems due to increased senescent cells. In several studies shortened telomeres have been associated with age related sarcopenia, atherosclerotic cardiovascular disease, and cancer. However, there is still the question whether telomere length causes these diseases or if the diseases cause shortened telomeres. Hence, the shortening of telomeres complies with antagonistic pleiotropy theory. The trade-off exists as the cell benefits from telomerase which prevents permanent growth arrest but telomere shortening is associated with functional loss.

Free Radical Theory
Another example related to aging is the Free Radical theory. Free Radical Theory suggests that the free radicals, which are being produced by aerobic respiration, are causing oxidative stress to be put on the body. This oxidative stress will result in aging and lead to death. Oxygen centered radicals are very reactive and can cause the accumulation of damage on lipids, nucleic acids as well as proteins within the body. This accumulation of damage on the biological molecules changes the framework and leads to a reduction in the molecules' activity levels. Lipid peroxides accumulate in the membrane phospholipids, which in turn diminishes the mitochondrial membrane's effectiveness as a barrier. The process of DNA transcription and translation also acquires oxidative damage. The result is alterations in the base pairings of the DNA sequence. Research has found that DNA mutations from free radical damage are highly uncommon but would still lead to the build up of damaged proteins as well as decreased biological activity.However, a large body of research has shown that antioxidants (compounds that defuse free radicals) do not slow aging or age-related diseases, putting the popular free radical theory into doubt.

See also
Evolution of ageing

== References ==