  The ejection of DNA from a bacterial virus (``phage'') into its host cell is
a biologically important example of the translocation of a macromolecular chain
along its length through a membrane. The simplest mechanism for this motion is
diffusion, but in the case of phage ejection a significant driving force
derives from the high degree of stress to which the DNA is subjected in the
viral capsid. The translocation is further sped up by the ratcheting and
entropic forces associated with proteins that bind to the viral DNA in the host
cell cytoplasm. We formulate a generalized diffusion equation that includes
these various pushing and pulling effects and make estimates of the
corresponding speed-ups in the overall translocation process. Stress in the
capsid is the dominant factor throughout early ejection, with the pull due to
binding particles taking over at later stages. Confinement effects are also
investigated, in the case where the phage injects its DNA into a volume
comparable to the capsid size. Our results suggest a series of in vitro
experiments involving the ejection of DNA into vesicles filled with varying
amounts of binding proteins from phage whose state of stress is controlled by
ambient salt conditions or by tuning genome length.


  Phi-values are experimental measures of how the kinetics of protein folding
is changed by single-site mutations. Phi-values measure energetic quantities,
but are often interpreted in terms of the structures of the transition state
ensemble. Here we describe a simple analytical model of the folding kinetics in
terms of the formation of protein substructures. The model shows that
Phi-values have both structural and energetic components. In addition, it
provides a natural and general interpretation of "nonclassical" Phi-values
(i.e., less than zero, or greater than one). The model reproduces the
Phi-values for 20 single-residue mutations in the alpha-helix of the protein
CI2, including several nonclassical Phi-values, in good agreement with
experiments.


  In this review, we discuss the applications of the theory of birth-and-death
processes to problems in biology, primarily, those of evolutionary genomics.
The mathematical principles of the theory of these processes are briefly
described. Birth-and-death processes, with some straightforward additions such
as innovation, are a simple, natural formal framework for modeling a vast
variety of biological processes such as population dynamics, speciation, genome
evolution, including growth of paralogous gene families and horizontal gene
transfer, and somatic evolution of cancers. We further describe how empirical
data, e.g., distributions of paralogous gene family size, can be used to choose
the model that best reflects the actual course of evolution among different
versions of birth-death-and-innovation models. It is concluded that
birth-and-death processes, thanks to their mathematical transparency,
flexibility and relevance to fundamental biological process, are going to be an
indispensable mathematical tool for the burgeoning field of systems biology.


  Most population models assume that individuals within a given population are
identical, that is, the fundamental role of variation is ignored. Inhomogeneous
models of populations and communities allow for birth and death rates to vary
among individuals; recently, theorems of existence and asymptotic of solutions
of such models were investigated. Here we develop another approach to modeling
heterogeneous populations by reducing the model to the Cauchy problem for a
special system of ODEs. As a result, the total population size and current
distribution of the vector-parameter can be found in explicit analytical form
or computed effectively. The developed approach is extended to the models of
inhomogeneous communities.


  The rank ordered distribution of the codon usage frequencies for 123
bacteriae is best fitted by a three parameters function that is the sum of a
constant, an exponential and a linear term in the rank n. The parameters depend
(two parabolically) from the total GC content. The rank ordered distribution of
the amino acids is fitted by a straight line. The Shannon entropy computed over
all the codons is well fitted by a parabola in the GC content, while the
partial entropies computed over subsets of the codons show peculiar different
behavior, exhibiting therefore a first conspiracy effect. Moreover the sum of
the codon usage frequencies over particular sets, e.g. with C and A
(respectively G and U) as i-th nucleotide, shows a clear linear dependence from
the GC content, exhibiting another conspiracy effect.


  We propose a novel variant of the algorithm by Simson et al. [R. Simson, E.D.
Sheets, K. Jacobson, Biophys. J. 69, 989 (1995)]. Their algorithm was developed
to detect transient confinement zones in experimental single particle tracking
trajectories of diffusing membrane proteins or lipids. We show that our
algorithm is able to detect confinement in a wider class of confining potential
shapes than Simson et al.'s one. Furthermore it enables to detect not only
temporary confinement but also jumps between confinement zones. Jumps are
predicted by membrane skeleton fence and picket models. In the case of
experimental trajectories of $\mu$-opioid receptors, which belong to the family
of G-protein-coupled receptors involved in a signal transduction pathway, this
algorithm confirms that confinement cannot be explained solely by rigid fences.


  The evolution of the full repertoire of proteins encoded in a given genome is
mostly driven by gene duplications, deletions, and sequence modifications of
existing proteins. Indirect information about relative rates and other
intrinsic parameters of these three basic processes is contained in the
proteome-wide distribution of sequence identities of pairs of paralogous
proteins. We introduce a simple mathematical framework based on a stochastic
birth-and-death model that allows one to extract some of this information and
apply it to the set of all pairs of paralogous proteins in seven model
organisms. It was found that the histogram of sequence identities p generated
by an all-to-all alignment of all protein sequences encoded in a genome is well
fitted with a power-law form ~p^(-gamma) with the value of the exponent gamma
around 4 for the majority of organisms used in this study. This implies that
the intra-protein variability of substitution rates is best described by the
Gamma-distribution with the exponent alpha ~ 0.33. We separately measure the
short-term (``raw'') duplication and deletion rates r*_dup, r*_del which
include gene copies that will be removed soon after the duplication event and
their dramatically reduced long-term counterparts r_dup, r_del. Systematic
trends of each of the four duplication/deletion rates with the total number of
genes in the genome were analyzed. All but the deletion rate of recent
duplicates r*_del were shown to systematically increase with N_genes.
Abnormally flat shapes of sequence identity histograms observed for yeast and
human are consistent with lineages leading to these organisms undergoing one or
more whole-genome duplications.


  Rapoport's rule claims that latitudinal ranges of plant and animal species
are generally smaller at low than at high latitudes. However, doubts as to the
generality of the rule have been expressed, because studies providing evidence
against the rule are more numerous than those in support of it. In groups for
which support has been provided, the trend of increasing latitudinal ranges
with latitude is restricted to or at least most distinct at high latitudes,
suggesting that the effect may be a local phenomenon, for example the result of
glaciations. Here we test the rule using two models, a simple one-dimensional
one with a fixed number of animals expanding in a northern or southerly
direction only, and the evolutionary/ecological Chowdhury model using birth,
ageing, death, mutation, speciation, prey-predator relations and food levels.
Simulations with both models gave results contradicting Rapoport's rule. In the
first, latitudinal ranges were roughly independent of latitude, in the second,
latitudinal ranges were greatest at low latitudes, as also shown empirically
for some well studied groups of animals.


  An inhomogeneous helicoidal nearest-neighbor model with continuous degrees of
freedom is shown to predict the same DNA melting properties as traditional
long-range Ising models, for free DNA molecules in solution, as well as
superhelically stressed DNA with a fixed linking number constraint. Without
loss of accuracy, the continuous degrees of freedom can be discretized using a
minimal number of discretization points, yielding an effective transfer matrix
model of modest dimension (d=36). The resulting algorithms to compute DNA
melting profiles are both simple and efficient.


  Apparent competition between species is believed to be one of the principle
driving forces that structure ecological communities, although the precise
mecha nisms have yet to be characterized. Here we develop a model system that
isolates phage-mediated interactions by neutralizing resource competition using
two genetically identical Bordetella bronchiseptica strains that differ only in
that one is the carrier of a phage and the other is susceptible to the phage.
We observe and quantify the competitive advantage of the bacterial strain
bearing the prophage in both invading and in resisting invasion by bacteria
susceptible to the phage, and use our measurements to develop a mathematical
model of phage-mediated competition. The model predicts, and experimental
evidence confirms, that the competitive advantage conferred by the phage
depends only on the relative phage pathology and is independent of other phage
and host parameters. This work combines experimental and mathematical
approaches to the study of phage-driven competition, and provides an
experimentally tested framework for evaluation of the effects of
pathogens/parasites on interspecific competition.


  The kinetic folding of RNA sequences into secondary structures is modeled as
a complex adaptive system, the components of which are possible RNA structural
rearrangements (SRs) and their associated bases and base pairs. RNA bases and
base pairs engage in local stacking interactions that determine the
probabilities (or fitnesses) of possible SRs. Meanwhile, selection operates at
the level of SRs; an autonomous stochastic process periodically (i.e., from one
time step to another) selects a subset of possible SRs for realization based on
the fitnesses of the SRs. Using examples based on selected natural and
synthetic RNAs, the model is shown to qualitatively reproduce characteristic
(nonlinear) RNA folding dynamics such as the attainment by RNAs of alternative
stable states. Possible applications of the model to the analysis of properties
of fitness landscapes, and of the RNA sequence to structure mapping are
discussed.


  Models of coevolution have recently been studied, in which a species is
defined by a genome in the form of a bitstring, and the interactions between
species i and j are given by a fixed matrix with independent, randomly
distributed elements M_{ij}. A consequence of the stochastic independence is
that species whose genotypes differ even by a single bit may have completely
different phenotypes. This is clearly unrealistic, as closely related species
should be similar in their interactions with the rest of the ecosystem. Here we
therefore study a model, in which the M_{ij} are correlated to a controllable
degree. We calculate, both analytically and numerically, the correlation
function for matrix elements M_{ij} and M_{kl} versus the Hamming distance
between the bitstrings representing the species. We compare Monte Carlo
simulations of coevolution models with uncorrelated and correlated
interactions. In particular, we consider the lifetimes of individual species.
The species-lifetime distribution is close to a power law with an exponent near
-2 in both uncorrelated and correlated cases. The durations of quasi-steady
states and power spectral densities for the diversity indices display
noticeable differences. However, some qualitative features, like 1/f behaviour
in power spectral densities for the diversity indices, are not affected by the
correlations in the interaction matrix.


  We present a base-pairing model of oligonuleotide duplex formation and show
in detail its equivalence to the Nearest-Neighbour dimer methods from fits to
free energy of duplex formation data for short DNA-DNA and DNA-RNA hybrids
containing only Watson Crick pairs. In this approach the connection between
rank-deficient polymer and rank-determinant oligonucleotide parameter, sets for
DNA duplexes is transparent. The method is generalised to include RNA/DNA
hybrids where the rank-deficient model with 11 dimer parameters in fact
provides marginally improved predictions relative to the standard method with
16 independent dimer parameters ($\Delta G$ mean errors of 4.5 and 5.4 %
respectively).


  The so-called Mixed Feedback Loop (MFL) is a small two-gene network where
protein A regulates the transcription of protein B and the two proteins form a
heterodimer. It has been found to be statistically over-represented in
statistical analyses of gene and protein interaction databases and to lie at
the core of several computer-generated genetic networks. Here, we propose and
mathematically study a model of the MFL and show that, by itself, it can serve
both as a bistable switch and as a clock (an oscillator) depending on kinetic
parameters. The MFL phase diagram as well as a detailed description of the
nonlinear oscillation regime are presented and some biological examples are
discussed. The results emphasize the role of protein interactions in the
function of genetic modules and the usefulness of modelling RNA dynamics
explicitly.


  Successful cell division requires a tight regulation of chromosome motion via
the activity of molecular motors. Many of the key players at the origin of the
forces generating the movement have been identified, but their spatial and
temporal organization remains elusive. The protein complex Kinetochore on the
chromosome associates with microtubules emanating from one of the spindle poles
and drives the chromosome toward the pole. Chromokinesin motors on the
chromosome arms also interact with microtubules, ejecting the chromosome away
from the pole. In animal cells, a monooriented chromosome (associated to a
single pole) periodically switches between phases of poleward and away from the
pole movement[, a behavior tentatively explained so far by the existence of a
complex switching mechanism within the kinetochore itself. Here we show that
the interplay between the morphology of the mitotic spindle and the collective
kinetics of chromokinesins can account for the highly non-linear periodic
chromosome motion. Our analysis provides a natural explanation for the origin
of chromosome directional instability and for the mechanism by which
chromosomes feel their position in space.


  The question what determines the structure of natural food webs has been
listed among the nine most important unanswered questions in ecology. It arises
naturally from many problems related to ecosystem stability and resilience. The
traditional view is that population-dynamical stability is crucial for
understanding the observed structures. But phylogeny (evolutionary history) has
also been suggested as the dominant mechanism. Here we show that observed
topological features of predatory food webs can be reproduced to unprecedented
accuracy by a mechanism taking into account only phylogeny, size constraints,
and the heredity of the trophically relevant traits of prey and predators. The
analysis reveals a tendency to avoid resource competition rather than apparent
competition. In food webs with many parasites this pattern is reversed.


  Background: Caenorhabditis elegans is a major model system in biology, yet
very little is known about its biology outside the laboratory. Especially, its
unusual mode of reproduction with self-fertile hermaphrodites and facultative
males raises the question of its frequency of outcrossing in natural
populations. Results: We describe the first analysis of C. elegans individuals
sampled directly from natural populations. C. elegans is found predominantly in
the dauer stage and with a very low frequency of males compared with
hermaphrodites. While C. elegans was previously shown to display a low
worldwide genetic diversity, we find by comparison a surprisingly high local
genetic diversity of C. elegans populations; this local diversity is
contributed in great part by immigration of new alleles rather than by
mutation. Our results on heterozygote frequency, male frequency and linkage
disequilibrium furthermore show that selfing is the predominant mode of
reproduction in C. elegans natural populations, yet that infrequent outcrossing
events occur, at a rate of approximately 1%. Conclusions: Our results give a
first insight in the biology of C. elegans in the natural populations. They
demonstrate that local populations of C. elegans are genetically diverse and
that a low frequency of outcrossing allows for the recombination of these
locally diverse genotypes.


  The translocation of structured RNA or DNA molecules through narrow pores
necessitates the opening of all base pairs. Here, we study the interplay
between the dynamics of translocation and base-pairing theoretically, using
kinetic Monte Carlo simulations and analytical methods. We find that the
transient formation of basepairs that do not occur in the ground state can
significantly speed up translocation.


  Signal transduction pathways are largely conserved throughout the animal
kingdom. The repertoire of pathways is limited and each pathway is used in
different intercellular signaling events during the development of a given
animal. For example, Wnt signaling is recruited, sometimes redundantly with
other molecular pathways, in four cell specification events during
Caenorhabditis elegans vulva development, including the activation of vulval
differentiation. Strikingly, a recent study finds that Wnts act to repress
vulval differentiation in the nematode Pristionchus pacificus, demonstrating
evolutionary flexibility in the use of intercellular signaling pathways.


  Various biological system models have been proposed in systems biology, which
are based on the complex biological reactions kinetic of various components.
These models are not practical because we lack of kinetic information. In this
paper, it is found that the enzymatic reaction and multi-order reaction rate is
often controlled by the transport of the reactants in biological systems. A
Computational Fluid Dynamic (CFD) approach, which is based on transport of the
components and kinetics of biological reactions, is introduced for biological
system modeling. We apply this approach to a biological wastewater treatment
system for the study of metabolism of organic carbon substrates and the
population of microbial. The results show that CFD model coupled with reaction
kinetics is more accurate and more feasible than kinetic models for biological
system modeling.


  This is an introductory review of deterministic mutation-selection models for
asexual populations (i.e., quasispecies theory) and related topics. First, the
basic concepts of fitness, mutations, and sequence space are introduced.
Different types of mutation-selection dynamics are defined and their relation
to problems of statistical physics are outlined. Then the stationary population
distribution in simple, single peak fitness landscapes is discussed at length,
with particular emphasis on the error threshold phenomenon. Extensions of the
theory covering e.g. epistatic interactions, diploid organisms,
semiconservative replication and time-dependent fitness peaks are briefly
described. A further section is devoted to randomly rugged fitness landscapes,
which may display fitness correlations of various degree as well as extended
neutral networks. The final two sections address evolutionary dynamics in both
rugged and smooth fitness landscapes, and provide a brief overview of pertinent
experiments.


  Proteins, chain molecules of amino acids, behave in ways which are similar to
each other yet quite distinct from standard compact polymers. We demonstrate
that the Flory theorem, derived for polymer melts, holds for compact protein
native state structures and is not incompatible with the existence of
structured building blocks such as $\alpha$-helices and $\beta$-strands. We
present a discussion on how the notion of the thickness of a polymer chain,
besides being useful in describing a chain molecule in the continuum limit,
plays a vital role in interpolating between conventional polymer physics and
the phase of matter associated with protein structures.


  Motivated by experiments in which single-stranded DNA with a short hairpin
loop at one end undergoes unforced diffusion through a narrow pore, we study
the first passage times for a particle, executing one-dimensional brownian
motion in an asymmetric sawtooth potential, to exit one of the boundaries. We
consider the first passage times for the case of classical diffusion,
characterized by a mean-square displacement of the form $< (\Delta x)^2> \sim
t$, and for the case of anomalous diffusion or subdiffusion, characterized by a
mean-square displacement of the form $< (\Delta x)^2> \sim t^{\gamma}$ with
$0<\gamma<1$. In the context of classical diffusion, we obtain an expression
for the mean first passage time and show that this quantity changes when the
direction of the sawtooth is reversed or, equivalently, when the reflecting and
absorbing boundaries are exchanged. We discuss at which numbers of `teeth' $N$
(or number of DNA nucleotides) and at which heights of the sawtooth potential
this difference becomes significant. For large $N$, it is well known that the
mean first passage time scales as $N^2$. In the context of subdiffusion, the
mean first passage time does not exist. Therefore we obtain instead the
distribution of first passage times in the limit of long times. We show that
the prefactor in the power relation for this distribution is simply the
expression for the mean first passage time in classical diffusion. We also
describe a hypothetical experiment to calculate the average of the first
passage times for a fraction of passage events that each end within some time
$t^*$. We show that this average first passage time scales as $N^{2/\gamma}$ in
subdiffusion.


  In this work we consider time reversal of the increasing geometrical
progression of the population of a simple biological species without any
enemies (predators) in the appropriate environment with unlimited resources
(food, territory, etc.). It is shown that such time reversal corresponds to
appearance of the cannibalism, i.e. self-predaciousness or self-damping
phenomena which can be described by a type of difference Verhulst equation.


  We present several new results pertaining to haplotyping. These results
concern the combinatorial problem of reconstructing haplotypes from incomplete
and/or imperfectly sequenced haplotype fragments. We consider the complexity of
the problems Minimum Error Correction (MEC) and Longest Haplotype
Reconstruction (LHR) for different restrictions on the input data.
Specifically, we look at the gapless case, where every row of the input
corresponds to a gapless haplotype-fragment, and the 1-gap case, where at most
one gap per fragment is allowed. We prove that MEC is APX-hard in the 1-gap
case and still NP-hard in the gapless case. In addition, we question earlier
claims that MEC is NP-hard even when the input matrix is restricted to being
completely binary. Concerning LHR, we show that this problem is NP-hard and
APX-hard in the 1-gap case (and thus also in the general case), but is
polynomial time solvable in the gapless case.


  This paper develops a quasispecies model where cells can adopt a two-cell
survival strategy. Within this strategy, pairs of cells join together, at which
point one of the cells sacrifices its own replicative ability for the sake of
the other cell. We develop a simplified model for the evolutionary dynamics of
this process, allowing us to solve for the steady-state using standard
approaches from quasispecies theory. We find that our model exhibits two
distinct regimes of behavior: At low concentrations of limiting resource, the
two-cell strategy outcompetes the single-cell survival strategy, while at high
concentrations of limiting resource, the single-cell survival strategy
dominates. Associated with the two solution regimes of our model is a
localization to delocalization transition over the portion of the genome coding
for the multicell strategy, analogous to the error catastrophe in standard
quasispecies models. The existence of such a transition indicates that
multicellularity can emerge because natural selection does not act on specific
cells, but rather on replicative strategies. Within this framework, individual
cells become the means by which replicative strategies are propagated. Such a
framework is therefore consistent with the concept that natural selection does
not act on individuals, but rather on populations.


  A formalism is developed which allows to determine the locations of all local
symmetry axes of three-dimensional particles with overall icosahedral symmetry.
It relies on the fact that the root system of the non-crystallographic Coxeter
group H_3 encodes the locations of the planes of reflection that generate the
discrete rotational symmetries of the particles. Via an appropriate extension
of the root system, new planes of reflection are introduced which determine
local axes of rotational symmetry. An easy-to-implement formalism is derived
that allows to compute the surface structure of any three-dimensional
icosahedral particle with local symmetries. It can be used also for particles
with overall octahedral and tetrahedral symmetry in conjunction with the root
systems of the corresponding reflection groups.
  Applications to viruses are discussed explicitly. It is shown that the
concept of quasi-equivalence in Caspar-Klug Theory corresponds to the special
case of local six-fold symmetry axes contained in the theory developed here,
and the corresponding geometries can hence be obtained with this formalism
based on the root system of H_3.
  Moreover, as a by-product, the theory answers the long-standing open question
why only certain types of capsomeres, i.e. clusters of protein subunits, are
observed in the surface structures of viruses. Since the types of the
capsomeres are determined by the orders of the local symmetry axes on which
they are located, the possible types of capsomeres are restricted by the
spectrum of local symmetry axes allowed by the theory. Based on this we
determine the spectrum of all capsomere types that may occur in viral capsids
and give explicit examples for the lower-order cases.


  We introduce an analytical model for population dynamics with intra-specific
competition, mutation and assortative mating as basic ingredients. The set of
equations that describes the time evolution of population size in a mean-field
approximation may be decoupled. We find a phase transition leading to sympatric
speciation as a parameter that quantifies competition strength is varied. This
transition, previously found in a computational model, occurs to be of first
order.


  Parapatric speciation is studied using an individual--based model with sexual
reproduction. We combine the theory of mutation accumulation for biological
ageing with an environmental selection pressure that varies according to the
individuals geographical positions and phenotypic traits. Fluctuations and
genetic diversity of large populations are crucial ingredients to model the
features of evolutionary branching and are intrinsic properties of the model.
Its implementation on a spatial lattice gives interesting insights into the
population dynamics of speciation on a geographical landscape and the
disruptive selection that leads to the divergence of phenotypes. Our results
suggest that assortative mating is not an obligatory ingredient to obtain
speciation in large populations at low gene flow.


  Based upon the correlation matrix of the human promoter sequences, a complex
network is constructed to capture the principal relationships between these
promoters. It is a complex network has the properties of the right-skewed
degree distribution and the clustering simultaneously, i.e., a hierarchical
structure. An eigenvector centrality (EC) based method is used to reconstruct
this hierarchical structure.


  We precisely position a noncanonical intron in the odd second copy of
tRNAAsp(GTC) gene in the newly sequenced crenarchaea S.acidocaldarius. The
uniform assortment of some features from normal aspartate tDNA and some from
those corresponding to non-standard amino acids conduce us to conjecture it to
be a novel tRNA gene, probably coding for a modified aspartate residue. Further
we reposition intron in tRNAHis(GUG) gene in P.aerophilum.The BHB motif at the
exon-intron boundaries are re-analyzed and found to support our conjectures.


  The immune system normally protects the human host against death by
infection. However, when an immune response is mistakenly directed at self
antigens, autoimmune disease can occur. We describe a model of protein
evolution to simulate the dynamics of the adaptive immune response to antigens.
Computer simulations of the dynamics of antibody evolution show that different
evolutionary mechanisms, namely gene segment swapping and point mutation, lead
to different evolved antibody binding affinities. Although a combination of
gene segment swapping and point mutation can yield a greater affinity to a
specific antigen than point mutation alone, the antibodies so evolved are
highly cross-reactive and would cause autoimmune disease, and this is not the
chosen dynamics of the immune system. We suggest that in the immune system a
balance has evolved between binding affinity and specificity in the mechanism
for searching the amino acid sequence space of antibodies.


  We explore the interplay between the protein-protein interactions network and
the expression of the interacting proteins. It is shown that interacting
proteins are expressed in significantly more similar cellular concentrations.
This is largely due to interacting pairs which are part of protein complexes.
We solve a generic model of complex formation and show explicitly that
complexes form most efficiently when their members have roughly the same
concentrations. Therefore, the observed similarity in interacting protein
concentrations could be attributed to optimization for efficiency of complex
formation.


  We investigate the Volterra ecosystem driven by correlated noises. The
competition of the predators induces an increasing in population density of the
predators. The competition of the preys, however, leads the predators to decay.
The predators may have better stability under strong correlated noises. The
predators undergo a sensitivity to a random environment, whereas the preys
exhibit a surprising endurance to this stochasticity.


  We present a statistical analysis of biological evolution processes.
Specifically, we study the stochastic replication-mutation-death model where
the population of a species may grow or shrink by birth or death, respectively,
and additionally, mutations lead to the creation of new species. We rank the
various species by the chronological order by which they originate. The average
population N_k of the kth species decays algebraically with rank, N_k ~ M^{mu}
k^{-mu}, where M is the average total population. The characteristic exponent
mu=(alpha-gamma)/(alpha+beta-gamma)$ depends on alpha, beta, and gamma, the
replication, mutation, and death rates. Furthermore, the average population P_k
of all descendants of the kth species has a universal algebraic behavior, P_k ~
M/k.


  We compare and contrast the long-time dynamical properties of two
individual-based models of biological coevolution. Selection occurs via
multispecies, stochastic population dynamics with reproduction probabilities
that depend nonlinearly on the population densities of all species resident in
the community. New species are introduced through mutation. Both models are
amenable to exact linear stability analysis, and we compare the analytic
results with large-scale kinetic Monte Carlo simulations, obtaining the
population size as a function of an average interspecies interaction strength.
Over time, the models self-optimize through mutation and selection to
approximately maximize a community fitness function, subject only to
constraints internal to the particular model. If the interspecies interactions
are randomly distributed on an interval including positive values, the system
evolves toward self-sustaining, mutualistic communities. In contrast, for the
predator-prey case the matrix of interactions is antisymmetric, and a nonzero
population size must be sustained by an external resource. Time series of the
diversity and population size for both models show approximate 1/f noise and
power-law distributions for the lifetimes of communities and species. For the
mutualistic model, these two lifetime distributions have the same exponent,
while their exponents are different for the predator-prey model. The difference
is probably due to greater resilience toward mass extinctions in the food-web
like communities produced by the predator-prey model.


  We consider a nonlinear partial differential equation of conservation type to
describe the dynamics of vesicular stomatitis virus observed in aliquots of
fixed particle number taken from an evolving clone at periodic intervals of
time \cite{novella 95}. The changes in time behavior of fitness function
noticed in experimental data are related to a crossover exhibited by the
solutions to this equation in the transient regime for pulse-like initial
conditions. As a consequence, the average replication rate of the population is
predicted to reach a plateau as a power t^{-(1/2)}.


  Understanding cell function requires an accurate description of how noise is
transmitted through biochemical networks. We present an analytical result for
the power spectrum of the output signal of a biochemical network that takes
into account the correlations between the noise in the input signal (the
extrinsic noise) and the noise in the reactions that constitute the network
(the intrinsic noise). These correlations arise from the fact that the
reactions by which biochemical signals are detected necessarily affect the
signaling molecules and the detection components of the network simultaneously.
We show that anti-correlation between the extrinsic and intrinsic noise
enhances the robustness of zero-order ultrasensitive networks to biochemical
noise. We discuss the consequences of the correlation between extrinsic and
intrinsic noise for a modular description of noise transmission through large
biochemical networks in the context of the mitogen-activated protein kinase
cascade.


  The tethered-particle method is a single-molecule technique that has been
used to explore the dynamics of a variety of macromolecules of biological
interest. We give a theoretical analysis of the particle motions in such
experiments. Our analysis reveals that the proximity of the tethered bead to a
nearby surface (the microscope slide) gives rise to a volume-exclusion effect,
resulting in an entropic force on the molecule. This force stretches the
molecule, changing its statistical properties. In particular, the proximity of
bead and surface brings about intriguing scaling relations between key
observables (statistical moments of the bead) and parameters such as the bead
size and contour length of the molecule. We present both approximate analytic
solutions and numerical results for these effects in both flexible and
semiflexible tethers. Finally, our results give a precise,
experimentally-testable prediction for the probability distribution of the
distance between the polymer attachment point and the center of the mobile
bead.


  We investigate the translocation of a stiff polymer consisting of M monomers
through a nanopore in a membrane, in the presence of binding particles
(chaperones) that bind onto the polymer, and partially prevent backsliding of
the polymer through the pore. The process is characterized by the rates: k for
the polymer to make a diffusive jump through the pore, q for unbinding of a
chaperone, and the rate q kappa for binding (with a binding strength kappa);
except for the case of no binding kappa=0 the presence of the chaperones give
rise to an effective force that drives the translocation process. Based on a
(2+1) variate master equation, we study in detail the coupled dynamics of
diffusive translocation and (partial) rectification by the binding proteins. In
particular, we calculate the mean translocation time as a function of the
various physical parameters.


  The distribution of inequivalent geometries occurring during self-assembly of
the major capsid protein in thermodynamic equilibrium is determined based on a
master equation approach. These results are implemented to characterize the
assembly of SV40 virus and to obtain information on the putative pathways
controlling the progressive build-up of the SV40 capsid. The experimental
testability of the predictions is assessed and an analysis of the geometries of
the assembly intermediates on the dominant pathways is used to identify targets
for antiviral drug design.


  A vital constituent of a virus is its protein shell, called the viral capsid,
that encapsulates and hence provides protection for the viral genome. Assembly
models are developed for viral capsids built from protein building blocks that
can assume different local bonding structures in the capsid. This situation
occurs, for example, for viruses in the family of Papovaviridae, which are
linked to cancer and are hence of particular interest for the health sector.
More specifically, the viral capsids of the (pseudo-) T=7 particles in this
family consist of pentamers that exhibit two different types of bonding
structures. While this scenario cannot be described mathematically in terms of
Caspar-Klug Theory (Caspar and Klug 1962), it can be modelled via tiling theory
(Twarock 2004). The latter is used to encode the local bonding environment of
the building blocks in a combinatorial structure, called the assembly tree,
which is a basic ingredient in the derivation of assembly models for
Papovaviridae along the lines of the equilibrium approach of Zlotnick (Zlotnick
1994). A phase space formalism is introduced to characterize the changes in the
assembly pathways and intermediates triggered by the variations in the
association energies characterizing the bonds between the building blocks in
the capsid. Furthermore, the assembly pathways and concentrations of the
statistically dominant assembly intermediates are determined. The example of
Simian Virus 40 is discussed in detail.


  We present an analytically tractable variant of the voter model that provides
a quantitatively accurate description of beta-diversity (two-point correlation
function) in two tropical forests. The model exhibits novel scaling behavior
that leads to links between ecological measures such as relative species
abundance and the species area relationship.


  We discuss recent progress towards the establishment of important
structure-property-function relationships in eumelanins - key functional
bio-macromolecular systems responsible for photo-protection and immune response
in humans, and implicated in the development of melanoma skin cancer. We focus
on the link between eumelanin's secondary structure and optical properties such
as broad band UV-visible absorption and strong non-radiative relaxation; both
key features of the photo-protective function. We emphasise the insights gained
through a holistic approach combining optical spectroscopy with first
principles quantum chemical calculations, and advance the hypothesis that the
robust functionality characteristic of eumelanin is related to extreme chemical
and structural disorder at the secondary level. This inherent disorder is a low
cost natural resource, and it is interesting to speculate as to whether it may
play a role in other functional bio-macromolecular systems.


  In this work, we discovered a fundamental connection between selection for
protein stability and emergence of preferred structures of proteins. Using
standard exact 3-dimensional lattice model we evolve sequences starting from
random ones and determining exact native structure after each mutation.
Acceptance of mutations is biased to select for stable proteins. We found that
certain structures, wonderfold, are independently discovered numerous times as
native states of stable proteins in many unrelated runs of selection. Diversity
of sequences that fold into wonderfold structures gives rise to superfamilies,
i.e. sets of dissimilar sequences that fold into the same or very similar
structures. Wonderfolds appear to be the most designable structures out of
complete set of compact lattice proteins. Furthermore, proteins having
wondefolds as their native structure tend to be most thermostable among all
evolved proteins. This effect is purely due to the favorable geometric
properties of wonderfolds and, thus, dominates any dependence on sequences. The
present work establishes a model of prebiotic structure selection, which
identifies dominant structural patterns emerging upon optimization of proteins
for survival in hot environment. Convergently discovered prebiotic initial
superfamilies with wonderfold structures could have served as a seed for
subsequent biological evolution involving gene duplications and divergence.


  A method is presented that, when used in conjunction with single molecule
experimental techniques, allows for the extraction of rates and mechanical
properties of a biomolecule undergoing transitions between mechanically
distinct states. This analysis enables the exploration of systems where the
lifetimes of survival are of order of the intrinsic time constant of the
experimental apparatus; permitting the study of kinetic events whose transition
rates are an order of magnitude (or two) larger than those that can be studied
with traditional averaging methods. Using current experimental techniques, this
allows for the study of biomolecules whose lifetime of survival in a given
state are as low as milliseconds down to microseconds. The relevant observable
is the auto-correlation function of the experimental probe that is attached to
the biomolecule of interest. General solutions are expressed in terms of a
series. Closed form solutions are found for two physically opposing limits:
where transitions between the mechanically distinct states of the biomolecule
occur on either much faster or much slower time scales than those governing the
motion of the experimental probe. Motivated by the derivation of these opposing
bounds, two series solutions for the general case are then presented. We
present an error analysis for truncating each series at arbitrary order and
obtain a range of parameters for which this method could be applicable to the
study of the two state biomolecular problem. Finally, both series (up to third
order) are expressed for the two state problem when the system obeys Markov
statistics. These solutions should be amiable to the analysis of experimental
data, expanding temporal analysis of data from single molecule experiments.


  We consider the folding of a self-avoiding homopolymer on a lattice, with
saturating hydrogen bond interactions. Our goal is to numerically evaluate the
statistical distribution of the topological genus of pseudoknotted
configurations. The genus has been recently proposed for classifying
pseudoknots (and their topological complexity) in the context of RNA folding.
We compare our results on the distribution of the genus of pseudoknots, with
the theoretical predictions of an existing combinatorial model for an
infinitely flexible and stretchable homopolymer. We thus obtain that steric and
geometric constraints considerably limit the topological complexity of
pseudoknotted configurations, as it occurs for instance in real RNA molecules.
We also analyze the scaling properties at large homopolymer length, and the
genus distributions above and below the critical temperature between the
swollen phase and the compact-globule phase, both in two and three dimensions.


  Including spatial structure and stochastic noise invalidates the classical
Lotka-Volterra picture of stable regular population cycles emerging in models
for predator-prey interactions. Growth-limiting terms for the prey induce a
continuous extinction threshold for the predator population whose critical
properties are in the directed percolation universality class. Here, we discuss
the robustness of this scenario by considering an ecologically inspired
stochastic lattice predator-prey model variant where the predation process
includes next-nearest-neighbor interactions. We find that the corresponding
stochastic model reproduces the above scenario in dimensions 1< d \leq 4, in
contrast with mean-field theory which predicts a first-order phase transition.
However, the mean-field features are recovered upon allowing for
nearest-neighbor particle exchange processes, provided these are sufficiently
fast.


  We give an exact solution to the Kolmogorov equation describing genetic drift
for an arbitrary number of alleles at a given locus. This is achieved by
finding a change of variable which makes the equation separable, and therefore
reduces the problem with an arbitrary number of alleles to the solution of a
set of equations that are essentially no more complicated than that found in
the two-allele case. The same change of variable also renders the Kolmogorov
equation with the effect of mutations added separable, as long as the mutation
matrix has equal entries in each row. Thus this case can also be solved exactly
for an arbitrary number of alleles. The general solution, which is in the form
of a probability distribution, is in agreement with the previously known
results--which were for the cases of two and three alleles only. Results are
also given for a wide range of other quantities of interest, such as the
probabilities of extinction of various numbers of alleles, mean times to these
extinctions, and the means and variances of the allele frequencies. To aid
dissemination, these results are presented in two stages: first of all they are
given without derivations and too much mathematical detail, and then
subsequently derivations and a more technical discussion are provided.


  It has recently been demonstrated that many biological networks exhibit a
scale-free topology where the probability of observing a node with a certain
number of edges (k) follows a power law: i.e. p(k) ~ k^-g. This observation has
been reproduced by evolutionary models. Here we consider the network of
protein-protein interactions and demonstrate that two published independent
measurements of these interactions produce graphs that are only weakly
correlated with one another despite their strikingly similar topology. We then
propose a physical model based on the fundamental principle that (de)solvation
is a major physical factor in protein-protein interactions. This model
reproduces not only the scale-free nature of such graphs but also a number of
higher-order correlations in these networks. A key support of the model is
provided by the discovery of a significant correlation between number of
interactions made by a protein and the fraction of hydrophobic residues on its
surface. The model presented in this paper represents the first physical model
for experimentally determined protein-protein interactions that comprehensively
reproduces the topological features of interaction networks. These results have
profound implications for understanding not only protein-protein interactions
but also other types of scale-free networks.


  We investigate the evolutionary dynamics of a finite population of RNA
sequences adapting to a neutral fitness landscape. Despite the lack of
differential fitness between viable sequences, we observe typical properties of
adaptive evolution, such as increase of mean fitness over time and punctuated
equilibrium transitions. We discuss the implications of these results for
understanding evolution at high mutation rates, and extend the relevance of the
quasispecies concept to finite populations and time scales. Our results imply
that the quasispecies concept and neutral drift are not complementary concepts,
and that the relative importance of each is determined by a combination of
population size and mutation rate.


  We calculate the equation of state of DNA under tension for the case that the
DNA features loops. Such loops occur transiently during DNA condensation in the
presence of multivalent ions or sliding cationic protein linkers. The
force-extension relation of such looped DNA modelled as a wormlike chain is
calculated via path integration in the semiclassical limit. This allows us to
determine rigorously the high stretching asymptotics. Notably the functional
form of the force-extension curve resembles that of straight DNA, yet with a
strongly renormalized apparent persistence length. That means that the
experimentally extracted single molecule elasticity does not necessarily
reflect the bare DNA stiffness only, but can also contain additional
contributions that depend on the overall chain conformation and length.


  In this paper we present a network model to study the impact of spatial
distribution of constituents, coupling between them and diffusive processes in
the context of biological situations. The model is in terms of network of
mobile elements whose internal dynamics is given by differential equations
exhibiting switching and/or oscillatory behaviour. To make the model more
consistent with the underlying biological phenomena we incorporate properties
like growth and decay into the network.
  We characterise this network by calculating the usual network measures like
network efficiency, entropy growth, vertex degree distribution, geodesic
lengths, centrality as well as fractal dimensions and generalised entropy. It
is seen that the model can demonstrate the features of both scale free networks
as well as small worlds network in different parameter domains. The formation
of spatio-temporal patterns is another feature of such networks which makes
them appealing for understanding broad qualitative aspects of problems like
cell differentiation and synchronization. The response of the network to
various attack strategies(isolated and multiple) is also studied.


  The coexistence of different viral strains (quasispecies) within the same
host are nowadays observed for a growing number of viruses, most notably HIV,
Marburg and Ebola, but the conditions for the formation and survival of new
strains have not yet been understood. We present a model of HIV quasispecies
competition, that describes the conditions of viral quasispecies coexistence
under different immune system conditions. Our model incorporates both T and B
cells responses, and we show that the role of B cells is important and additive
to that of T cells. Simulations of coinfection (simultaneous infection) and
superinfection (delayed secondary infection) scenarios in the early stages
(days) and in the late stages of the infection (years) are in agreement with
emerging molecular biology findings. The immune response induces a competition
among similar phenotypes, leading to differentiation (quasi-speciation), escape
dynamics and complex oscillations of viral strain abundance. We found that the
quasispecies dynamics after superinfection or coinfection has time scales of
several months and becomes even slower when the immune system response is weak.
Our model represents a general framework to study the speed and distribution of
HIV quasispecies during disease progression, vaccination and therapy.


  A generalized computational method for folding proteins with a fully
transferable potential and geometrically realistic all-atom model is presented
and tested on seven different helix bundle proteins. The protocol, which
includes graph-theoretical analysis of the ensemble of resulting folded
conformations, was systematically applied and consistently produced structure
predictions of approximately 3 Angstroms without any knowledge of the native
state. To measure and understand the significance of the results, extensive
control simulations were conducted. Graph theoretic analysis provides a means
for systematically identifying the native fold and provides physical insight,
conceptually linking the results to modern theoretical views of protein
folding. In addition to presenting a method for prediction of structure and
folding mechanism, our model suggests that a accurate all-atom amino acid
representation coupled with a physically reasonable atomic interaction
potential (that does not require optimization to the test set) and hydrogen
bonding are essential features for a realistic protein model.


  A stochastic evolutionary dynamics of two strategies given by 2 x 2 matrix
games is studied in finite populations. We focus on stochastic properties of
fixation: how a strategy represented by a single individual wins over the
entire population. The process is discussed in the framework of a random walk
with arbitrary hopping rates. The time of fixation is found to be identical for
both strategies in any particular game. The asymptotic behavior of the fixation
time and fixation probabilities in the large population size limit is also
discussed. We show that fixation is fast when there is at least one pure
evolutionary stable strategy (ESS) in the infinite population size limit, while
fixation is slow when the ESS is the coexistence of the two strategies.


  Most population models assume that individuals within a given population are
identical, that is, the fundamental role of variation is ignored. Here we
develop a general approach to modeling heterogeneous populations with discrete
evolutionary time step. The theory is applied to population dynamics of natural
rotifer populations. We show that under particular conditions the behavior of
the inhomogeneous model possesses complex transition regimes, which depends
both on the mean and the variance of the initial parameter distribution and the
final state of the population depends on the least possible value from the
domain of the parameter. The question of persistence of the population is
discussed.


  We describe a stochastic birth-and-death model of evolution of horizontally
transferred genes in microbial populations. The model is a generalization of
the stochastic model described by Berg and Kurland and includes five
parameters: the rate of mutational inactivation, selection coefficient,
immigration rate (i.e., rate of arrival of a novel sequence from outside of the
recipient population), within-population horizontal transmission rate, and
population size. The model of Berg and Kurland included four parameters,
namely, mutational inactivation, selection coefficient, population size, and
transmission rate. However, the effect of transmission was disregarded in the
interpretation of the results, and the overall conclusion was that horizontally
acquired sequences can be fixed in a population only when they confer a
substantial selective advantage onto the recipient and therefore are subject to
strong positive selection. By contrast, analysis of the present model in
different domains of parameter values shows that, as long as the rate of
within-population horizontal transmission is comparable to the mutational
inactivation rate and there is even a low rate of immigration, horizontally
acquired sequences can be fixed in the population or at least persist for a
long time even when they are neutral or slightly deleterious. The available
biological data strongly suggest that intense within-population and even
between-populations gene flows are realistic for at least some prokaryotic
species and environments. Therefore our modeling results are compatible with
the notion of a pivotal role of horizontal gene transfer in the evolution of
prokaryotes.


  A simple weakly frequency dependent model for the dynamics of a population
with a finite number of types is proposed, based upon an advantage of being
rare. In the infinite population limit, this model gives rise to a non-smooth
dynamical system that reaches its globally stable equilibrium in finite time.
This dynamical system is sufficiently simple to permit an explicit solution,
built piecewise from solutions of the logistic equation in continuous time. It
displays an interesting tree-like structure of coalescing components.


  The simultaneous growth of multiple microbial species is a problem of
fundamental ecological interest. In media containing more than one
growth-limiting substrate, multiple species can coexist. The question then
arises: Can single-species data predict the existence and stability of
mixed-culture steady states in mixed-substrate environments? This question has
been extensively studied with the help of resource-based models. These studies
have shown that the single-species data required to predict mixed-culture
behavior consists of the growth isoclines and consumption vectors, which in
turn are determined from single-substrate data by making specific assumptions
about the kinetics of mixed-substrate growth. Here, we show that these
assumptions are not valid for microbial growth on mixtures of substitutable
substrates. However, the theory can be developed by determining the growth
isoclines and consumption vectors directly from the mixed-substrate data, thus
obviating the need for specific assumptions about the kinetics of
mixed-substrate growth. We show furthermore that in addition to the growth
isoclines and consumption vectors, the single-species, mixed-substrate data
yields a new family of curves, which we call the consumption curves.
Consideration of the growth isoclines and the consumption curves yields deeper
insights into the behavior of the mixed cultures. It yields a priori bounds on
the substrate concentrations achieved during coexistence, permits the extension
of the theory to systems in which the growth isoclines are non-monotonic, and
clarifies earlier results obtained by considering only the growth isoclines.


  When steady state chemostat cultures are abruptly exposed to substrate-excess
conditions, they exhibit long lags before adjusting to the new environment. The
identity of the rate-limiting step for this slow response can be inferred from
the initial yields and specific growth rates measured by exposing steady state
cultures at various dilution rates to substrate-excess conditions. We measured
these parameters for glucose-limited cultures of E. coli ML308 growing at
various dilution rates between 0.03 and 0.6 1/hr. In all the cases, the initial
yields were 20-30% less than the steady state yields. The decline of the yield
implies that overflow metabolism is triggered in response to excess glucose. It
is therefore unlikely that the initial response of the cells is limited by
substrate uptake. The initial specific growth rates of cultures growing at low
dilution rates (D = 0.03, 0.05, 0.075, 0.1, 0.3 1/hr) were significantly higher
than the steady state specific growth rates. However, the increment in the
specific growth rate decreased with the dilution rate, and at D=0.6 1/hr, there
was no improvement in the specific growth rate. The initial specific growth
rates varied hyperbolically with the dilution, decreasing sharply at dilution
rates below 0.1 1/hr and saturating at D=0.6 1/hr. This is consistent with a
picture in which the initial response is limited by the activity of glutamate
dehydrogenase.


  Fermentation is generally modelled by kinetic equations giving the time
evolutions for biomass, substrate, and product concentrations. Although these
equations can be solved analytically in simple cases if substrate/product
inhibition and biomass death are included, they are typically solved
numerically. We propose an analytical treatment of the kinetic equations
--including cell death and an arbitrary number of inhibitions-- in which
constant yield needs not be assumed. Equations are solved in phase space, i.e.
the biomass concentration is written explicitly as a function of the substrate
concentration.


  We predict a novel temperature-driven phase transition of DNA below the
melting transition. The additional, intermediate phase exists for repetitive
sequences, when the two strands have different lengths. In this phase, the
excess bases of the longer strand are completely absorbed as bulge loops inside
the helical region. When the temperature is lowered, the excess bases desorb
into overhanging ends, resulting in a contour length change. This continuous
transition is in many aspects analogous to Bose Einstein condensation. Weak
sequence disorder renders the transition discontinuous.


  The TGF-beta pathway plays a central role in tissue homeostasis and
morphogenesis. It transduces a variety of extracellular signals into
intracellular transcriptional responses that control a plethora of cellular
processes, including cell growth, apoptosis, and differentiation. We use
computational modeling to show that coupling of signaling with receptor
trafficking results in a highly versatile signal-processing unit, able to sense
by itself absolute levels of ligand, temporal changes in ligand concentration,
and ratios of multiple ligands. This coupling controls whether the response of
the receptor module is transient or permanent and whether or not different
signaling channels behave independently of each other. Our computational
approach unifies seemingly disparate experimental observations and suggests
specific changes in receptor trafficking patterns that can lead to phenotypes
that favor tumor progression.


  The NF-kB signalling system is involved in a variety of cellular processes
including immune response, inflammation, and apoptosis. Recent experiments have
found oscillations in the nuclear-cytoplasmic translocation of the NF-kB
transcription factor. How the cell uses the oscillations to differentiate input
conditions and send specific signals to downstream genes is an open problem. We
shed light on this issue by examining the small core network driving the
oscillations, which, we show, is designed to produce periodic spikes in nuclear
NF-kB concentration. The oscillations can be used to regulate downstream genes
in a variety of ways. In particular, we show that genes to whose operator sites
NF-kB binds and dissociates fast can respond very sensitively to changes in the
input signal, with effective Hill coefficients in excess of 20.


  The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the
`Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic
neurotransmission, and whose other members include receptors for glycine,
gamma-aminobutyric acid, and serotonin. Cryo-electron microscopy has yielded a
three dimensional structure of the nAChR in its closed state. However, the
exact nature and location of the channel gate remains uncertain. Although the
transmembrane pore is constricted close to its center, it is not completely
occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A.
Model calculations suggest that such a constriction may form a hydrophobic
gate, preventing movement of ions through a channel. We present a detailed and
quantitative simulation study of the hydrophobic gating model of the nicotinic
receptor, in order to fully evaluate this hypothesis. We demonstrate that the
hydrophobic constriction of the nAChR pore indeed forms a closed gate.
Potential of mean force (PMF) calculations reveal that the constriction
presents a barrier of height ca. 10 kT to the permeation of sodium ions,
placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A
radius hydrophobic pore can form a functional barrier to the permeation of a 1
A radius Na+ ion. Using a united atom force field for the protein instead of an
all atom one retains the qualitative features but results in differing
conductances, showing that the PMF is sensitive to the detailed molecular
interactions.


  Being the HIV-1 Protease (HIV-1-PR) an essential enzyme in the viral life
cycle, its inhibition can control AIDS. The folding of single domain proteins,
like each of the monomers forming the HIV-1-PR homodimer, is controlled by
local elementary structures (LES, folding units stabilized by strongly
interacting, highly conserved, as a rule hydrophobic, amino acids). These LES
have evolved over myriad of generations to recognize and strongly attract each
other, so as to make the protein fold fast and be stable in its native
conformation. Consequently, peptides displaying a sequence identical to those
segments of the monomers associated with LES are expected to act as competitive
inhibitors and thus destabilize the native structure of the enzyme. These
inhibitors are unlikely to lead to escape mutants as they bind to the protease
monomers through highly conserved amino acids which play an essential role in
the folding process. The properties of one of the most promising inhibitors of
the folding of the HIV-1-PR monomers found among these peptides is demonstrated
with the help of spectrophotometric assays and CD spectroscopy.


  It is shown that a small subset of modes which are likely to be involved in
protein functional motions of large amplitude can be determined by retaining
the most robust normal modes obtained using different protein models. This
result should prove helpful in the context of several applications proposed
recently, like for solving difficult molecular replacement problems or for
fitting atomic structures into low-resolution electron density maps. Moreover,
it may also pave the way for the development of methods allowing to predict
such motions accurately.


  The study of Planktonic Foraminifera abundances permits to obtain climatic
curves on the basis of percentage ratio between tropical and temperate/polar
forms. Climatic changes were controlled by several phenomena as: (i)
Milankovitch's cycles, produced by variations of astronomical parameters such
as precession, obliquity and eccentricity; (ii) continental geodynamic
evolution and orogenic belt; (iii) variations of atmospheric and oceanic
currents; (iv) volcanic eruptions; (v) meteor impacts. But while astronomical
parameters have a quasi-regular periodicity, the other phenomena can be
considered as "noise signal" in natural systems. The interplay between cyclical
astronomical variations, the "noise signal" and the intrinsic nonlinearity of
the ecologic system produces strong glacial or interglacial period according to
the stochastic resonance phenomenon.


  We present Monte-Carlo (MC) simulations of the stretching of a single 30 nm
chromatin fiber. The model approximates the DNA by a flexible polymer chain
with Debye-H\"uckel electrostatics and uses a two-angle zig-zag model for the
geometry of the linker DNA connecting the nucleosomes. The latter are
represented by flat disks interacting via an attractive Gay-Berne potential.
Our results show that the stiffness of the chromatin fiber strongly depends on
the linker DNA length. Furthermore, changing the twisting angle between
nucleosomes from 90 deg to 130 deg increases the stiffness significantly. An
increase in the opening angle from 22 deg to 34 deg leads to softer fibers for
small linker lengths. We observe that fibers containing a linker histone at
each nucleosome are stiffer compared to those without the linker histone. The
simulated persistence lengths and elastic moduli agree with experimental data.
Finally, we show that the chromatin fiber does not behave as an isotropic
elastic rod, but its rigidity depends on the direction of deformation:
chromatin is much more resistant to stretching than to bending.


  The information regarding the structure of a single protein is encoded in the
network of interacting amino acids. Considering each protein as a weighted and
unweighted network of amino acids we have analyzed a total of forty nine
protein structures that covers the three branches of life on earth. Our results
show that the probability degree distribution of network connectivity follows
Poisson's distribution; whereas the probability strength distribution does not
follow any known distribution. However, the average strength of amino acid node
depends on its degree (k). For some of the proteins, the strength of a node
increases linearly with k. On the other hand, for a set of other proteins,
although the strength increases linaerly with k for smaller values of k, we
have not obtained any clear functional relationship of strength with degree at
higher values of k. The results also show that the weight of the amino acid
nodes belonging to the highly connected nodes tend to have a higher value. The
result that the average clustering coefficient of weighted network is less than
that of unweighted network implies that the topological clustering is generated
by edges with low weights. The ratio of average clustering coefficients of
protein network to that of the corresponding classical random network varies
linearly with the number (N) of amino acids of a protein; whereas the ratio of
characteristic path lengths varies logarithmically with N. The power law
behaviour of clustering coefficients of weighted and unweighted network as a
function of degree k indicates that the network has a signature of hierarchical
network. It has also been observed that the network is of assortative type.


  In euryarchaeal methanogen M.kandleri and in Nanoarchaea N. equitans some of
the missing tRNA genes are embedded in others. We argue from bioinformatic
evidence that position specific intron splicing is the key behind co-location
of these tRNA genes.


  We study the stochastic dynamics of sequences evolving by single site
mutations, segmental duplications, deletions, and random insertions. These
processes are relevant for the evolution of genomic DNA. They define a
universality class of non-equilibrium 1D expansion-randomization systems with
generic stationary long-range correlations in a regime of growing sequence
length. We obtain explicitly the two-point correlation function of the sequence
composition and the distribution function of the composition bias in sequences
of finite length. The characteristic exponent $\chi$ of these quantities is
determined by the ratio of two effective rates, which are explicitly calculated
for several specific sequence evolution dynamics of the universality class.
Depending on the value of $\chi$, we find two different scaling regimes, which
are distinguished by the detectability of the initial composition bias. All
analytic results are accurately verified by numerical simulations. We also
discuss the non-stationary build-up and decay of correlations, as well as more
complex evolutionary scenarios, where the rates of the processes vary in time.
Our findings provide a possible example for the emergence of universality in
molecular biology.


  Molecular domestication of a transposable element is defined as its
functional recruitment by the host genome. To date, two independent events of
molecular domestication of the P transposable element have been described: in
the Drosophila obscura species group and in the Drosophila montium species
subgroup. These P neogenes consist to stationary, non repeated sequences,
potentially encoding 66 kDa repressor-like proteins (RLs). Here we investigate
the function of the montium P neogenes. We provide evidence for the presence of
RLs proteins in two montium species (D. tsacasi and D. bocqueti) specifically
expressed in adult and larval brain and gonads. We tested the hypothesis that
the montium P neogenes function is related to the repression of the
transposition of distant related mobile P elements which coexist in the genome.
Our results strongly suggest that the montium P neogenes are not recruited to
down regulate the P element transposition. Given that all the proteins encoded
by mobile or stationary P homologous sequences show a strong conservation of
the DNA Binding Domain, we tested the capacity of the RLs proteins to bind DNA
in vivo. Immunstaining of polytene chromosomes in D. melanogaster transgenic
lines strongly suggest that montium P neogenes encode proteins that bind DNA in
vivo. RLs proteins show multiple binding to the chromosomes. We suggest that
the property recruited in the case of the montium P neoproteins is their DNA
binding property. The possible functions of these neogenes are discussed.


  Recent work has shown that expression level is the main predictor of a
gene’s evolutionary rate, and that more highly expressed genes evolve
slower. A possible explanation for this observation is selection for proteins
which fold properly despite mistranslation, in short selection for
translational robustness. Translational robustness leads to the somewhat
paradoxical prediction that highly expressed genes are extremely tolerant to
missense substitutions but nevertheless evolve very slowly. Here, we study a
simple theoretical model of translational robustness that allows us to gain
analytic insight into how this paradoxical behavior arises.


  We introduce a simple computational model that, with a microscopic dynamics
driven by natural selection and mutation alone, allows the description of true
speciation events. A statistical analysis of the so generated evolutionary tree
captures realistic features showing power laws for frequency distributions in
time and size. Albeit these successful predictions, the difficulty in obtaining
punctuated dynamics with mass extinctions suggests the necessity of decoupling
micro and macro-evolutionary mechanisms in agreement with some ideas of Gould's
and Eldredge's theory of punctuated equilibrium.


  The current understanding of supernova and gamma-ray burst events suggests
important effects on the biosphere if one of more of them happened to strike
the earth in the past. In this paper we evaluate the possibility that life
extinctions which probably occurred due to excess of radiation occurring in the
geologic past might have left a genetic signature on surviving species. We
emphasize the signatures of these extinctions, proposing a quantitative model
to evaluate the surviving probability of the species, based on kinetic aspects
of the frequency of mutations and the DNA repair rate.


  Background: The relative contribution of different aging-associated processes
to the age phenotype may differ among individuals, creating variability in
aging manifestations among age-peers. Capturing this variability can
significantly advance understanding the aging and mortality. An index of
age-associated health disorders (deficits), called a "frailty index" (FI),
appears to be a promising characteristic of such processes. In this study we
address the connections of the FI with age focusing on disabled individuals who
might be at excessive risk of frailty. Methods: The National Long Term Care
Survey (NLTCS) assessed health and functioning of the U.S. elderly in 1982,
1984, 1989, 1994, and 1999. Detailed information for our sample was assessed
from about 26,700 interviews. The individual FI is defined as a proportion of
deficits for a given person. We perform cross-sectional empirical analysis of
the FI age-patterns. Results: FI in the NLTCS exhibits accelerated (quadratic)
increase with age. Deficits might accumulate faster among the elderly who, at
younger ages, had a low mean FI ("healthy" group) than a high FI ("disabled"
group). Age-patterns for "healthy" and "disabled" groups converge at advanced
ages. The rate of deficit accumulation is sex-sensitive. Convergence of the
(sex-specific) FI for "healthy" and "disabled" groups in later ages determines
biological age limits, associated with given levels of health-maintenance in
the society, which correspond to 109.4 years for females and 92.5 years for
males. Conclusions: The FI can be employed as a measure of biological age and
population heterogeneity for modeling aging processes and mortality in elderly
individuals.


  We present an approximate analytic study of our previously introduced model
of evolution including the effects of genetic exchange. This model is motivated
by the process of bacterial transformation. We solve for the velocity, the rate
of increase of fitness, as a function of the fixed population size, $N$. We
find the velocity increases with $\ln N$, eventually saturated at an $N$ which
depends on the strength of the recombination process. The analytical treatment
is seen to agree well with direct numerical simulations of our model equations.


  We introduce a Markov model for the evolution of a gene family along a
phylogeny. The model includes parameters for the rates of horizontal gene
transfer, gene duplication, and gene loss, in addition to branch lengths in the
phylogeny. The likelihood for the changes in the size of a gene family across
different organisms can be calculated in O(N+hM^2) time and O(N+M^2) space,
where N is the number of organisms, $h$ is the height of the phylogeny, and M
is the sum of family sizes. We apply the model to the evolution of gene content
in Preoteobacteria using the gene families in the COG (Clusters of Orthologous
Groups) database.


  To better understand mortality change with age capturing the variability in
individuals' rates of aging, we performed comprehensive analysis of statistical
properties of a cumulative index of age-associated disorders (deficits), called
a "frailty index" (FI). This index is calculated as the proportion of the
health deficits in an individual. It is found, first, that frequency,
time-to-death, mortality-rate, and relative-risk-of-death exhibit remarkably
similar FI- and age- patterns. Second, the FI, on the one hand, and mortality
rate and relative risk, on the other hand, also exhibit similar age patterns
with accelerated increase up to oldest-old ages and with subsequent
deceleration and even decline. Third, distribution of the FI with time-to-death
is sharper than that of age with time-to-death. These and related findings
support the conclusion that the FI can describe aging processes and population
heterogeneity. We also discuss the ability of the FI to capture physiological
processes underlying aging both on individual and population levels.


  The frequencies of A, C, G and T in mitochondrial DNA vary among species due
to unequal rates of mutation between the bases. The frequencies of bases at
four-fold degenerate sites respond directly to mutation pressure. At 1st and
2nd positions, selection reduces the degree of frequency variation. Using a
simple evolutionary model, we show that 1st position sites are less constrained
by selection than 2nd position sites, and therefore that the frequencies of
bases at 1st position are more responsive to mutation pressure than those at
2nd position. We define a similarity measure between amino acids that is a
function of 8 measured physical properties. We define a proximity measure for
each amino acid, which is the average similarity between an amino acid and all
others that are accessible via single point mutations in the genetic code. We
also define a responsiveness for each amino acid, which measures how rapidly an
amino acid frequency changes as a result of mutation pressure acting on the
base frequencies. There is a strong correlation between responsiveness and
proximity, and both these quantities are also correlated with the mutability of
amino acids estimated from the mtREV substitution rate matrix. We also consider
the variation of base frequencies between strands and between genes on a
strand. These trends are consistent with the patterns expected from analysis of
the variation among genomes


  A useful approach to the mathematical analysis of large-scale biological
networks is based upon their decompositions into monotone dynamical systems.
This paper deals with two computational problems associated to finding
decompositions which are optimal in an appropriate sense. In graph-theoretic
language, the problems can be recast in terms of maximal sign-consistent
subgraphs. The theoretical results include polynomial-time approximation
algorithms as well as constant-ratio inapproximability results. One of the
algorithms, which has a worst-case guarantee of 87.9 percent from optimality,
is based on the semidefinite programming relaxation approach of
Goemans-Williamson. The algorithm was implemented and tested on a Drosophila
segmentation network and an Epidermal Growth Factor Receptor pathway model, and
it was found to perform close to optimally.


  This paper develops a simplified model for sexual replication within the
quasispecies formalism. We assume that the genomes of the replicating organisms
are two-chromosomed and diploid, and that the fitness is determined by the
number of chromosomes that are identical to a given master sequence. We also
assume that there is a cost to sexual replication, given by a characteristic
time $ \tau_{seek} $ during which haploid cells seek out a mate with which to
recombine. If the mating strategy is such that only viable haploids can mate,
then when $ \tau_{seek} = 0 $, it is possible to show that sexual replication
will always outcompete asexual replication. However, as $ \tau_{seek} $
increases, sexual replication only becomes advantageous at progressively higher
mutation rates. Once the time cost for sex reaches a critical threshold, the
selective advantage for sexual replication disappears entirely. The results of
this paper suggest that sexual replication is not advantageous in small
populations per se, but rather in populations with low replication rates. In
this regime, the cost for sex is sufficiently low that the selective advantage
obtained through recombination leads to the dominance of the strategy. In fact,
at a given replication rate and for a fixed environment volume, sexual
replication is selected for in high populations because of the reduced time
spent finding a reproductive partner.


  Mixed-substrate microbial growth is among the most intensely studied systems
in molecular microbiology. Several mathematical models have been developed to
account for the genetic regulation of such systems, especially those resulting
in diauxic growth. In this work, we compare the dynamics of three such models
(Narang, Biotech. Bioeng., 59, 116, 1998; Thattai & Shraiman, Biophys. J, 85,
744, 2003; Brandt et al, Water Research, 38, 1004, 2004). We show that these
models are dynamically similar - the initial motion of the inducible enzymes in
all the models is described by Lotka-Volterra equations for competing species.
The dynamic similarity occurs because in all the models, the inducible enzymes
possess properties characteristic of competing species: Their synthesis is
autocatalytic, and they inhibit each other. Despite this dynamic similarity,
the models vary with respect to the range of dynamics captured. The Brandt et
al model captures only the diauxic growth pattern, whereas the remaining two
models capture both diauxic and non-diauxic growth patterns. The models also
differ with respect to the mechanisms that generate the mutual inhibition
between the enzymes. In the Narang model, the mutual inhibition occurs because
the enzymes for each substrate enhance the dilution of the enzymes for the
other substrate. In the Thattai & Shraiman model, the mutual inhibition is
entirely due to competition for the phosphoryl groups.


  We analyze gene expression time-series data of yeast S. cerevisiae measured
along two full cell-cycles. We quantify these data by using q-exponentials,
gene expression ranking and a temporal mean-variance analysis. We construct
gene interaction networks based on correlation coefficients and study the
formation of the corresponding giant components and minimum spanning trees. By
coloring genes according to their cell function we find functional clusters in
the correlation networks and functional branches in the associated trees. Our
results suggest that a percolation point of functional clusters can be
identified on these gene expression correlation networks.


  For many infectious diseases, a small-world network on an underlying regular
lattice is a suitable simplified model for the contact structure of the host
population. It is well known that the contact network, described in this
setting by a single parameter, the small-world parameter $p$, plays an
important role both in the short term and in the long term dynamics of epidemic
spread. We have studied the effect of the network structure on models of immune
for life diseases and found that in addition to the reduction of the effective
transmission rate, through the screening of infectives, spatial correlations
may strongly enhance the stochastic fluctuations. As a consequence, time series
of unforced Susceptible-Exposed-Infected-Recovered (SEIR) models provide
patterns of recurrent epidemics with realistic amplitudes, suggesting that
these models together with complex networks of contacts are the key ingredients
to describe the prevaccination dynamical patterns of diseases such as measles
and pertussis. We have also studied the role of the host contact strucuture in
pathogen antigenic variation, through its effect on the final outcome of an
invasion by a viral strain of a population where a very similar virus is
endemic. Similar viral strains are modelled by the same infection and
reinfection parameters, and by a given degree of cross immunity that represents
the antigenic distance between the competing strains. We have found, somewhat
surprisingly, that clustering on the network decreases the potential to sustain
pathogen diversity.


  Motivated by recent research of aging in E. coli, we explore the effects of
aging on bacterial fitness. The disposable soma theory of aging was developed
to explain how differences in lifespans and aging rates could be linked to life
history trade-offs. Although generally applied for multicellular organisms, it
is also useful for exploring life history strategies of single celled organisms
such as bacteria. Starting from the Euler-Lotka equation, we propose a
mathematical model to explore how a finite lifespan effects fitness of
bacteria. We find that that there is surprisingly little loss of fitness when
the bacterium has limited opportunities to reproduce. Instead, the fitness
gained each time the bacteria reproduces decreases rapidly, and investing
resources to survive to reproduce the first few times is likely more
advantageous than investing additional resources to try to maintain cell
integrity longer.


  The transcription factor p53 is an important regulator of cell fate.
Mutations in p53 gene are associated with many cancers. In response to signals
such as DNA damage, p53 controls the transcription of a series of genes that
cause cell cycle arrest during which DNA damage is repaired, or triggers
programmed cell death that eliminates possibly cancerous cells wherein DNA
damage might have remained unrepaired. Previous experiments showed oscillations
in p53 level in response to DNA damage, but the mechanism of oscillation
remained unclear. Here we examine a model where the concentrations of p53
isoforms are regulated by Mdm22, Arf, Siah, and beta-catenin. The extent of DNA
damage is signalled through the switch-like activity of a DNA damage sensor,
the DNA-dependent protein kinase Atm. This switch is responsible for initiating
and terminating oscillations in p53 concentration. The strength of the DNA
damage signal modulates the number of oscillatory waves of p53 and Mdm22 but
not the frequency or amplitude of oscillations{a result that recapitulates
experimental findings. A critical fnding was that the phosphorylated form of
Nbs11, a member of the DNA break repair complex Mre11-Rad50-Nbs11 (MRN), must
augment the activity of Atm kinase. While there is in vitro support for this
assumption, this activity appears essential for p53 dynamics. The model
provides several predictions concerning, among others, degradation of the
phosphorylated form of p53, the rate of DNA repair as a function of DNA damage,
the sensitivity of p53 oscillation to transcription rates of SIAH, beta-CATENIN
and ARF, and the hysteretic behavior of active Atm kinase levels with respect
to the DNA damage signal


  The simplest model of a smart spatial redistribution of individuals is
proposed. A single-species population is considered, to be composed of two
discrete subpopulations inhabiting two stations; migration is a transfer
between them. The migration is not random and yields the maximization of net
reproduction, with respect to the transaction costs. The organisms are supposed
to be globally informed. Discrete time model is studied, since it shows all the
features of a smart migrations, while the continuous time case brings no new
knowledge but the technical problems. Some properties of the model are studied
and discussed.


  We consider a Susceptible-Infective-Recovered (SIR) model, where the
mechanism for the renewal of susceptibles is demographic, on a ring with next
nearest neighbour interactions, and a family of correlated pair approximations
(CPA), parametrized by a measure of the relative contributions of loops and
open triplets of the sites involved in the infection process. We have found
that the phase diagram of the CPA, at fixed coordination number, changes
qualitatively as the relative weight of the loops increases, from the phase
diagram of the uncorrelated pair approximation to phase diagrams typical of
one-dimensional systems. In addition, we have performed computer simulations of
the same model and shown that while the CPA with a constant correlation
parameter cannot describe the global behaviour of the model, a reasonable
description of the endemic equilibria as well as of the phase diagram may be
obtained by allowing the parameter to depend on the demographic rate.


  Protein crystallization in vivo provides some fascinating examples of
biological self-assembly. Here, we provide a selective survey to show the
diversity of functions for which protein crystals are used, and the physical
properties of the crystals thatare exploited. Where known, we emphasize how the
nature of the protein-protein interactions leads to control of the
crystallization behaviour.


  Biologists measure information in different ways. Neurobiologists and
researchers in bioinformatics often measure information using
information-theoretic measures such as Shannon's entropy or mutual information.
Behavioral biologists and evolutionary ecologists more commonly use
decision-theoretic measures, such the value of information, which assess the
worth of information to a decision maker. Here we show that these two kinds of
measures are intimately related in the context of biological evolution. We
present a simple model of evolution in an uncertain environment, and calculate
the increase in Darwinian fitness that is made possible by information about
the environmental state. This fitness increase -- the fitness value of
information -- is a composite of both Shannon's mutual information and the
decision-theoretic value of information. Furthermore, we show that in certain
cases the fitness value of responding to a cue is exactly equal to the mutual
information between the cue and the environment. In general the Shannon entropy
of the environment, which seemingly fails to take anything about organismal
fitness into account, nonetheless imposes an upper bound on the fitness value
of information.


  We study a version of the Tangled Nature model of evolutionary ecology
redefined in a phenotype space where mutants have properties correlated to
their parents. The model has individual-based dynamics whilst incorporating
species scale competitive constraints and a system scale resource constraint.
Multiple species arise that coexist in a species interaction network with
evolving global properties. Both the mean interaction strength and the network
connectance increase relative to the null system as mutualism becomes more
extensive. From a study of the dependence of average degree on the resource
level we extract the diversity-connectance relationship which conforms to the
hyperbolic form seen in field data. This is adjudged to arise as a consequence
of the evolutionary pressure to achieve positive interactions. The network
degree distributions conform more strongly to exponential than to the null
binomial distributions in all cases. This effect is believed to be caused by
correlations in the reproductive process. We also study how resource
availability influences the phenotypical lifetime distribution which is
approximately of power law form. We observe that the mean lifetime is inversely
related to the resource level.


  We study a class of growth algorithms for directed graphs that are candidate
models for the evolution of genetic regulatory networks. The algorithms involve
partial duplication of nodes and their links, together with innovation of new
links, allowing for the possibility that input and output links from a newly
created node may have different probabilities of survival. We find some
counterintuitive trends as parameters are varied, including the broadening of
indegree distribution when the probability for retaining input links is
decreased. We also find that both the scaling of transcription factors with
genome size and the measured degree distributions for genes in yeast can be
reproduced by the growth algorithm if and only if a special seed is used to
initiate the process.


  The compartmental model is a basic tool for studying signal propagation in
neurons, and, if the model parameters are adequately defined, it can also be of
help in the study of electrical or fluid transport. Here we show that the input
resistance, in different networks which simulate the passive properties of
neurons, is the result of an interplay between the relevant conductances,
morphology and size. These results suggest that neurons must grow in such a way
that facilitates the current flow. We propose that power consumption is an
important factor by which neurons attain their final morphological appearance.


  A new formalism for calculation of the partition function of single stranded
nucleic acids is presented. Secondary structures and the topology of structure
elements are the level of resolution that is used. The folding model deals with
matches, mismatches, symmetric and asymmetric interior loops, stacked pairs in
loop and dangling end regions, multi-branched loops, bulges and single base
stacking that might exist at duplex ends or at the ends of helices.
Calculations on short and long sequences show, that for short oligonucleotides,
a duplex formation often displays a two-state transition. However, for longer
oligonucleotides, the thermodynamic properties of the single self-folding
transition affects the transition nature of the duplex formation, resulting in
a population of intermediate hairpin species in the solution. The role of
intermediate hairpin species is analyzed in the case when a short
oligonucleotides (molecular beacons) have to reliably identify and hybridize to
accessible nucleotides within their targeted mRNA sequences. It is shown that
the enhanced specificity of the molecular beacons is a result of their
constrained conformational flexibility and the all-or-none mechanism of their
hybridization to the target sequence.


  The cell cycle is an orderly sequence of events which ultimately lead to the
division of a single cell into two daughter cells. In the case of DNA damage by
radiation or chemicals, the damage checkpoints in the $G_{1}$ and $G_{2}$
phases of the cell cycle are activated. This results in an arrest of the cell
cycle so that the DNA damage can be repaired. Once this is done, the cell
continues with its usual cycle of activity. We study a mathematical model of
the DNA damage checkpoint in the $G_{2}$ phase which arrests the transition
from the $G_{2}$ to the $M$ (mitotic) phase of the cell cycle. The tumor
suppressor protein p53 plays a key role in activating the pathways leading to
cell cycle arrest in mammalian systems. If the DNA damage is severe, the p53
proteins activate other pathways which bring about apoptosis, i.e., programmed
cell death. Loss of the p53 gene results in the proliferation of cells
containing damaged DNA, i.e., in the growth of tumors which may ultimately
become cancerous. There is some recent experimental evidence which suggests
that the mutation of a single copy of the p53 gene (in the normal cell each
gene has two identical copies) is sufficient to trigger the formation of
tumors. We study the effect of reducing the gene copy number of the p53 and two
other genes on cell cycle arrest and obtain results consistent with
experimental observations.


  The unicellular cyanobacterium Synechococcus PCC 7942 has been used as a
model organism for studies of prokaryotic circadian rhythms,
carbon-concentrating mechanisms, response to a variety of nutrient and
environmental stresses, and cell division. This paper presents the results of
the first proteomic exploratory study of Synechococcus PCC 7942. The proteome
was analyzed using two-dimensional gel electrophoresis followed by MALDI-TOF
mass spectroscopy, and database searching. Of 140 analyzed protein spots, 110
were successfully identified as 62 different proteins, many of which occurred
as multiple spots on the gel. The identified proteins were organized into 18
different functional categories reflecting the major metabolic and cellular
processes occurring in the cyanobacterial cells in the exponential growth
phase. Among the identified proteins, 14 previously unknown or considered to be
hypothetical are here shown to be true gene products in Synechococcus sp. PCC
7942, and may be helpful for annotation of the newly sequenced genome.


  Although ARTMAP and ART-based models were introduced in early 70's they were
not used in characterizing and classifying ecological observations. ART-based
models have been extensively used for classification models based on satellite
imagery. This report, to our knowledge, is the first application of ART-based
methods and specifically ARTMAP for predicting habitat selection and spatial
distribution of species. We compare the performance of ARTMAP to assess the
breeding success of three bird species (Lanius senator, Hippolais pallida, and
Calandrella brachydactyla) based on multi-spectral satellite imagery and
environmental variables. ARTMAP is superior both in terms of performance
(percent correctly classified - pcc = 1.00) and generalizability (pcc >0.96) to
those of feedforward multilayer backpropogation (>0.87, >0.65), linear and
quadratic discriminant analysis (>0.48, >0.46) and k-nearest neighbor (>0.82,
>0.66) methods. Compared to other methods, ARTMAP is able to incorporate new
observations with far less computational effort and can easily add data to
already trained models.