  Separate artificial neural network (ANN) models were developed from data in
two geographical regions and years apart for a marsh-nesting bird, the
red-winged blackbird Agelaius phoeniceus. Each model was independently tested
on the spatially and temporally distinct data from the other region to
determine how generalizable it was. The first model was developed to predict
occurrence of nests in two wetlands on Lake Erie, Ohio in 1995 and 1996. The
second model was developed to predict breeding success in two marshes in
Connecticut, USA in 1969 and 1970. Independent variables were vegetation
durability, stem density, stem/nest height, distance to open water, distance to
edge, and water depth. With input variable relevances, sensitivity analyses and
neural interpretation diagrams we were able to understand how the different
models predicted nest occurrence and breeding success and compare their
differences and similarities. Both models also predicted increasing nest
occurrence/breeding success with increasing water depth under the nest and
increasing distance to edge. However, relationships for prediction differed in
the models. Generalizability of the models was poor except when the marshes had
similar values of important variables in the model. ANN models performed better
than generalized linear models (GLM) on marshes with similar structures.
Generalizability of the models did not differ in nest occurrence and breeding
success data. Extensive testing also showed that the GLMs were not necessarily
more generalizable than ANNs, suggesting that ANN models make good definitions
of a study system but are too specific to generalize well to other ecologically
complex systems unless input variable distributions are very similar.


  We review the use of artificial neural networks, particularly the feedforward
multilayer perceptron with back-propagation for training (MLP), in ecological
modelling. Overtraining on data or giving vague references to how it was
avoided is the major problem. Various methods can be used to determine when to
stop training in artificial neural networks: 1) early stopping based on
cross-validation, 2) stopping after a analyst defined error is reached or after
the error levels off, 3) use of a test data set. We do not recommend the third
method as the test data set is then not independent of model development. Many
studies used the testing data to optimize the model and training. Although this
method may give the best model for that set of data it does not give
generalizability or improve understanding of the study system. The importance
of an independent data set cannot be overemphasized as we found dramatic
differences in model accuracy assessed with prediction accuracy on the training
data set, as estimated with bootstrapping, and from use of an independent data
set. The comparison of the artificial neural network with a general linear
model (GLM) as a standard procedure is recommended because a GLM may perform as
well or better than the MLP. MLP models should not be treated as black box
models but instead techniques such as sensitivity analyses, input variable
relevances, neural interpretation diagrams, randomization tests, and partial
derivatives should be used to make the model more transparent, and further our
ecological understanding which is an important goal of the modelling process.
Based on our experience we discuss how to build a MLP model and how to optimize
the parameters and architecture.


  This paper focuses on a mathematical model interpreting the prime number life
cycle of periodical cicadas.Changed the viewpoint to predators rather than the
prey,this model fits reality very well by utilizing some principles and
assumption.With the definition of the predator income,natural selection from
predators seems to be the main reason for such a long life cycle.Consequent
solution of this model is exactly the fact of real nature.


  For a fixed number of food levels on a lattice, the Chowdhury model for
ecosystem is simulated and found to give a number of different species which
approaches exponentially a plateau if the considered area is increased.


  It is shown that all simple ecological, i.e. population dynamical equations
(unlimited exponential population growth (or decrease) dynamics, logistic or
Verhulst equation, usual and generalized Lotka-Volterra equations) hold a
symmetry, called $C_{E}PT$ symmetry. Namely, all simple ecological dynamical
equations are invariant (symmetric) in respect to successive application of the
time reversal transformation - $T$, space coordinates reversal or parity
transformation - $P$, and predator-prey reversal transformation - $C_{E}$ that
changes preys in the predators or pure (healthy) in the impure (fatal)
environment, and vice versa. It is deeply conceptually analogous to remarkable
$CPT$ symmetry of the fundamental physical dynamical equations. Further, it is
shown that by more accurate, "microscopic" analysis, given $C_{E}PT$ symmetry
becomes explicitly broken.


  We study the primary DNA structure of four of the most completely sequenced
human chromosomes (including chromosome 19 which is the most dense in coding),
using Non-extensive Statistics. We show that the exponents governing the decay
of the coding size distributions vary between $5.2 \le r \le 5.7$ for the short
scales and $1.45 \le q \le 1.50$ for the large scales. On the contrary, the
exponents governing the decay of the non-coding size distributions in these
four chromosomes, take the values $2.4 \le r \le 3.2$ for the short scales and
$1.50 \le q \le 1.72$ for the large scales. This quantitative difference, in
particular in the tail exponent $q$, indicates that the non-coding (coding)
size distributions have long (short) range correlations. This non-trivial
difference in the DNA statistics is attributed to the non-conservative
(conservative) evolution dynamics acting on the non-coding (coding) DNA
sequences.


  We construct a minimal content-based realization of the duplication and
divergence model of genomic networks introduced by Wagner [A. Wagner, Proc.
Natl. Acad. Sci. {\bf 91}, 4387 (1994)] and investigate the scaling properties
of the directed degree distribution and clustering coefficient. We find that
the content based network exhibits crossover between two scaling regimes, with
log-periodic oscillations for large degrees. These features are not present in
the original gene duplication model, but inherent in the content based model of
Balcan and Erzan. The scaling exponents $\gamma_1$ and $\gamma_2=\gamma_1-1/2$
of the Balcan-Erzan model turn out to be robust under duplication and point
mutations, but get modified in the presence of splitting and merging of
strings. The clustering coefficient as a function of the degree, $C(d)$, is
found, for the Balcan-Erzan model, to behave in a way qualitatively similar to
the out-degree distribution, however with a very small exponent $\alpha_1=
1-\gamma_1$ and an envelope for the oscillatory part, which is essentially
flat, thus $\alpha_2= 0$. Under duplication and mutations including splitting
and merging of strings, $C(d)$ is found to decay exponentially.


  Apparent biodiversity on earth exists only if we compare different species
separated from their environments. Meanwhile coexisting species have to be
identical in terms of energetic interactions. Consider the biosphere as a
network of chemical reactions. This leads to the conclusion that forces which
drive life evolution may be found inside the process of the transformations of
molecules. From the thermodynamic point of view the system of reacting
particles reaches a steady state when the chemical potentials of reagents
become equal. Interactions between chemical compounds taking part in concurrent
reactions and the equilibration of their chemical potentials is the essence of
biological evolution and its only driving force.


  Many biomedical experiments are carried out by pooling individual biological
samples. However, pooling samples can potentially hide biological variance and
give false confidence concerning the data significance. In the context of
microarray experiments for detecting differentially expressed genes, recent
publications have addressed the problem of the efficiency of sample-pooling,
and some approximate formulas were provided for the power and sample size
calculations. It is desirable to have exact formulas for these calculations and
have the approximate results checked against the exact ones. We show that the
difference between the approximate and exact results can be large. In this
study, we have characterized quantitatively the effect of pooling samples on
the efficiency of microarray experiments for the detection of differential gene
expression between two classes. We present exact formulas for calculating the
power of microarray experimental designs involving sample pooling and technical
replications. The formulas can be used to determine the total numbers of arrays
and biological subjects required in an experiment to achieve the desired power
at a given significance level. The conditions under which pooled design becomes
preferable to non-pooled design can then be derived given the unit cost
associated with a microarray and that with a biological subject. This paper
thus serves to provide guidance on sample pooling and cost effectiveness. The
formulation in this paper is outlined in the context of performing microarray
comparative studies, but its applicability is not limited to microarray
experiments. It is also applicable to a wide range of biomedical comparative
studies where sample pooling may be involved.


  We consider some mathematical issues raised by the modelling of gene
networks. The expression of genes is governed by a complex set of regulations,
which is often described symbolically by interaction graphs. Once such a graph
has been established, there remains the difficult task to decide which
dynamical properties of the gene network can be inferred from it, in the
absence of precise quantitative data about their regulation. In this paper we
discuss a rule proposed by R.Thomas according to which the possibility for the
network to have several stationary states implies the existence of a positive
circuit in the corresponding interaction graph. We prove that, when properly
formulated in rigorous terms, this rule becomes a theorem valid for several
different types of formal models of gene networks. This result is already known
for models of differential or boolean type. We show here that a stronger
version of it holds in the differential setup when the decay of protein
concentrations is taken into account. This allows us to verify also the
validity of Thomas' rule in the context of piecewise-linear models and the
corresponding discrete models. We discuss open problems as well.


  A vital constituent of a virus is its protein shell, called the viral capsid,
that encapsulates and hence provides protection for the viral genome. Viral
capsids are usually spherical, and for a significant number of viruses exhibit
overall icosahedral symmetry. The corresponding surface lattices, that encode
the locations of the capsid proteins and intersubunit bonds, can be modelled by
Viral Tiling Theory.
  It has been shown in vitro that under a variation of the experimental
boundary conditions, such as the pH value and salt concentration, tubular
particles may appear instead of, or in addition to, spherical ones. In order to
develop models that describe the simultaneous assembly of both spherical and
tubular variants, and hence study the possibility of triggering tubular
malformations as a means of interference with the replication mechanism, Viral
Tiling Theory has to be extended to include tubular lattices with end caps.
This is done here for the case of Papovaviridae, which play a distinguished
role from the viral structural point of view as they correspond to all pentamer
lattices, i.e. lattices formed from clusters of five protein subunits
throughout. These results pave the way for a generalisation of recently
developed assembly models.


  A new N-dimensional vector space of DNA sequences over the Galois field of
the 64 codons (GF(64)) was recently presented. Now, in order to include
deletions and insertions (indel mutations), we have defined a new Galois field
over the set of elements X1X2X3 (C125), where Xi belong to {O, A, C, G, C}. We
have called this set, the extended triplet set and the elements X1X2X3, the
extended triplets. The order of the bases is derived from the Z64-algebra of
the genetic code -recently published-. Starting from the natural bijection phi:
GF(5^3)-> C125 between the polynomial representation of elements from GF(5^3)
and the elements X1X2X3, a novel Galois field over the set of elements X1X2X3
is defined. Taking the polynomial coefficients a0, a1, a2 belong to GF(5) and
the bijective function f: GF(5) ->{O, A, C, G, C}, where f(0) = O, f(1) = A,
f(2) = C, f(3) = G, f(4) = U, bijection phi is induced such that phi(a0 + a1x +
a2x^2) = (f(a1), f(a2), f(a0)) = (X1X2X3). Next, by means of the bijection phi
we define sum "+" and product "*" operations in the set of codons C125, in such
a way that the resultant field (C125, +, *) turns isomorphic to the Galois
Field GF(5^3). This field allows the definition of a novel N-dimensional vector
space (S) over the field GF (5^3) on the set of all 125^N sequences of extended
triplets in which all possible DNA sequence alignments of length N are
included. Here the "classical gap" produced by alignment algorithms corresponds
to the neutral element "O". It is verified that the homologous (generalized)
recombination between two homologous DNA duplexes involving a reciprocal
exchange of DNA sequences -e.g. between two chromosomes that carry the same
genetic loci- algebraically corresponds to the action of two automorphism pairs
(or two translation pairs) over two paired DNA duplexes.


  Transcriptional interactions in a cell are modulated by a variety of
mechanisms that prevent their representation as pure pairwise interactions
between a transcription factor and its target(s). These include, among others,
transcription factor activation by phosphorylation and acetylation, formation
of active complexes with one or more co-factors, and mRNA/protein degradation
and stabilization processes.
  This paper presents a first step towards the systematic, genome-wide
computational inference of genes that modulate the interactions of specific
transcription factors at the post-transcriptional level. The method uses a
statistical test based on changes in the mutual information between a
transcription factor and each of its candidate targets, conditional on the
expression of a third gene. The approach was first validated on a synthetic
network model, and then tested in the context of a mammalian cellular system.
By analyzing 254 microarray expression profiles of normal and tumor related
human B lymphocytes, we investigated the post transcriptional modulators of the
MYC proto-oncogene, an important transcription factor involved in
tumorigenesis. Our method discovered a set of 100 putative modulator genes,
responsible for modulating 205 regulatory relationships between MYC and its
targets. The set is significantly enriched in molecules with function
consistent with their activities as modulators of cellular interactions,
recapitulates established MYC regulation pathways, and provides a notable
repertoire of novel regulators of MYC function. The approach has broad
applicability and can be used to discover modulators of any other transcription
factor, provided that adequate expression profile data are available.


  Studies of species distributions on elevational gradients are essential to
understand principles of community organisation as well as to conserve species
in montane regions. This study examined the patterns of species richness,
abundance, composition, range sizes, and distribution of rainforest birds at 14
sites along an elevational gradient (500-1400 m) in the Kalakad-Mundanthurai
Tiger Reserve (KMTR) of the Western Ghats, India. In contrast to theoretical
expectation, resident bird species richness did not change significantly with
elevation although the species composition changed substantially (<10%
similarity) between the lowest and highest elevation sites. Constancy in
species richness was possibly due to relative constancy in productivity and
lack of elevational trends in vegetation structure. Elevational range size of
birds, expected to increase with elevation according to Rapoport's rule, was
found to show a contrasting inverse U-shaped pattern because species with
narrow elevational distributions, including endemics, occurred at both ends of
the gradient (below 800 m and above 1,200 m). Bird species composition also did
not vary randomly along the gradient as assessed using a hierarchy of null
models of community assembly, from completely unconstrained models to ones with
species richness and range-size distribution restrictions. Instead, bird
community composition was significantly correlated with elevation and tree
species composition of sites, indicating the influence of deterministic factors
on bird community structure. Conservation of low- and high-elevation areas and
maintenance of tree species composition against habitat alteration are
important for bird conservation in the southern Western Ghats rainforests.


  Distance based algorithms are a common technique in the construction of
phylogenetic trees from taxonomic sequence data. The first step in the
implementation of these algorithms is the calculation of a pairwise distance
matrix to give a measure of the evolutionary change between any pair of the
extant taxa. A standard technique is to use the log det formula to construct
pairwise distances from aligned sequence data. We review a distance measure
valid for the most general models, and show how the log det formula can be used
as an estimator thereof. We then show that the foundation upon which the log
det formula is constructed can be generalized to produce a previously unknown
estimator which improves the consistency of the distance matrices constructed
from the log det formula. This distance estimator provides a consistent
technique for constructing quartets from phylogenetic sequence data under the
assumption of the most general Markov model of sequence evolution.


  Experimental investigations of the biosynthesis of a number of proteins have
pointed out that part of the native structure can be acquired already during
translation. We carried out a comprehensive statistical analysis of some
average structural properties of proteins that have been put forward as
possible signatures of this progressive buildup process. Contrary to a
widespread belief, it is found that there is no major propensity of the amino
acids to form contacts with residues that are closer to the N terminus.
Moreover, it is found that the C terminus is significantly more compact and
locally-organized than the N one. Also this bias, though, is unlikely to be
related to vectorial effects, since it correlates with subtle differences in
the primary sequence. These findings indicate that even if proteins aquire
their structure vectorially no signature of this seems to be detectable in
their average structural properties.


  The aim of this article is to present a developed method that decomposes the
autofluorescence spectrum into the spectra of naturally occurring biochemical
components of biotissue. It requires knowledge of detailed spectrum behaviour
of different endogenous fluorophores. We have studied the main bio-markers in
human tissue and proposed a simple modelling algorithm for their spectra
shapes. The empirical method was tested theoretically by quantum-mechanical
calculations of the spectra in the unharmonic Morse potential approach.


  Techniques such as detrended fluctuation analysis (DFA) and its extensions
have been widely used to determine the nature of scaling in nucleotide
sequences. In this brief communication we show that tandem repeats which are
ubiquitous in nucleotide sequences can prevent reliable estimation of possible
long-range correlations. Therefore, it is important to investigate the presence
of tandem repeats prior to scaling exponent estimation.


  It is known since the early days of molecular biology that proteins locate
their specific targets on DNA up to two orders of magnitude faster than the
Smoluchowski 3D diffusion rate. It was the idea due to Delbruck that they are
non-specifically adsorbed on DNA, and sliding along DNA provides for the faster
1D search. Surprisingly, the role of DNA conformation was never considered in
this context. In this article, we explicitly address the relative role of 3D
diffusion and 1D sliding along coiled or globular DNA and the possibility of
correlated re-adsorbtion of desorbed proteins. We have identified a wealth of
new different scaling regimes. We also found the maximal possible acceleration
of the reaction due to sliding, we found that the maximum on the
rate-versus-ionic strength curve is asymmetric, and that sliding can lead not
only to acceleration, but in some regimes to dramatic deceleration of the
reaction.


  We present a method to investigate the kinetics of protein folding on a long
time-scale and the dynamics underlying the formation of secondary and tertiary
structures during the entire reaction. The approach is based on the formal
analogy between thermal and quantum diffusion: by writing the solution of the
Fokker-Planck equation for the time-evolution of a protein in a viscous
heat-bath in terms of a path integral, we derive a Hamilton-Jacobi variational
principle from which we are able to compute the most probable pathway of
folding. The method is applied to the folding of the Villin Headpiece
Subdomain, in the framework of a Go-model. We have found that, in this model,
the transition occurs through an initial collapsing phase driven by the
starting coil configuration and a later rearrangement phase, in which secondary
structures are formed and all computed paths display strong similarities. This
method is completely general, does not require the prior knowledge of any
reaction coordinate and represents an efficient tool to perfom ab-initio
simulations of the entire folding process with available computers.


  Gene expression data matrices often contain missing expression values. In
this paper, we describe a new algorithm, named improved fixed rank
approximation algorithm (IFRAA), for missing values estimations of the large
gene expression data matrices. We compare the present algorithm with the two
existing and widely used methods for reconstructing missing entries for DNA
microarray gene expression data: the Bayesian principal component analysis
(BPCA) and the local least squares imputation method (LLS). The three
algorithms were applied to four microarray data sets and two synthetic low-rank
data matrices. Certain percentages of the elements of these data sets were
randomly deleted, and the three algorithms were used to recover them. In
conclusion IFRAA appears to be the most reliable and accurate approach for
recovering missing DNA microarray gene expression data, or any other noisy data
matrices that are effectively low rank.


  Experimental evidence suggests that the folding and aggregation of the
amyloid $\beta$-protein (A$\beta$) into oligomers is a key pathogenetic event
in Alzheimer's disease (AD). Inhibiting the pathologic folding and
oligomerization of A$\beta$ could be effective in the prevention and treatment
of AD. Here, using all-atom molecular dynamics simulations in explicit solvent,
we probe the initial stages of folding of a decapeptide segment of A$\beta$,
A$\beta_{21-30}$, shown experimentally to nucleate the folding process. In
addition, we examine the folding of a homologous decapeptide containing an
amino acid substitution linked to hereditary cerebral hemorrhage with
amyloidosis--Dutch type, [Gln22]A$\beta_{21-30}$. We find that: (i) when the
decapeptide is in water, hydrophobic interactions and transient salt bridges
between Lys28 and either Glu22 or Asp23 are important in the formation of a
loop in the Val24--Lys28 region of the wild type decapeptide; (ii) in the
presence of salt ions, salt bridges play a more prominent role in the
stabilization of the loop; (iii) in water with a reduced density, the
decapeptide forms a helix, indicating the sensitivity of folding to different
aqueous environments; (iv) the ``Dutch'' peptide in water, in contrast to the
wild type peptide, fails to form a long-lived Val24--Lys28 loop, suggesting
that loop stability is a critical factor in determining whether A$\beta$ folds
into pathologic structures. Our results are relevant to understand the
mechanism of A$\beta$ peptide folding in different environments, such as intra-
and extracellular milieus or cell membranes, and how amino acid substitutions
linked to familial forms of amyloidosis cause disease.


  We study the mechanical properties of semiflexible polymers when the contour
length of the polymer is comparable to its persistence length. We compute the
exact average end-to-end distance and shape of the polymer for different
boundary conditions, and show that boundary effects can lead to significant
deviations from the well-known long-polymer results. We also consider the case
of stretching a uniformly charged biopolymer by an electric field, for which we
compute the average extension and the average shape, which is shown to be
trumpetlike. Our results also apply to long biopolymers when thermal
fluctuations have been smoothed out by a large applied field or force.


  Motivated by the formation of ring-like filament structures in the cortex of
plant and animal cells, we study the dynamics of a two-dimensional layer of
cytoskeletal filaments and motor proteins near a surface by a general continuum
theory. As a result of active processes, dynamic patterns of filament
orientation and density emerge via instabilities. We show that
self-organization phenomena can lead to the formation of stationary and
oscillating rings. We present state diagrams which reveal a rich scenario of
asymptotic behaviors and discuss the role of boundary conditions.


  Transfer-RNA genes in archaea often have introns intervening between exon
sequences. The structural motif at the boundary between exon and intron is the
bulge-helix-bulge. Computational investigations of these boundary structures in
H. marismortui lead us to propose that tRNA-isoleucine and tRNA-tryptophan
genes are co-located. Precise insilico identification of the splice-sites on
the bulges at the exon-intron boundaries conduce us to infer that a single
intron-containing composite tRNA-gene can give rise to more than one gene
produc.


  A cell's behavior is a consequence of the complex interactions between its
numerous constituents, such as DNA, RNA, proteins and small molecules. Cells
use signaling pathways and regulatory mechanisms to coordinate multiple
processes, allowing them to respond to and adapt to an ever-changing
environment. The large number of components, the degree of interconnectivity
and the complex control of cellular networks are becoming evident in the
integrated genomic and proteomic analyses that are emerging. It is increasingly
recognized that the understanding of properties that arise from whole-cell
function require integrated, theoretical descriptions of the relationships
between different cellular components. Recent theoretical advances allow us to
describe cellular network structure with graph concepts, and have revealed
organizational features shared with numerous non-biological networks. How do we
quantitatively describe a network of hundreds or thousands of interacting
components? Does the observed topology of cellular networks give us clues about
their evolution? How does cellular networks' organization influence their
function and dynamical responses? This article will review the recent advances
in addressing these questions.


  The structure of a stochastic nonlinear gene regulatory network is uncovered
by studying its response to input signal generators. Four applications are
studied in detail: a nonlinear connection of two linear systems, the design of
a logic pulse, a molecular amplifier and the interference of three signal
generators in E2F1 regulatory element. The gene interactions are presented
using molecular diagrams that have a precise mathematical structure and retain
the biological meaning of the processes.


  Interpretation of empirical results based on a taxa's lifetime distribution
shows apparently conflicting results. Species' lifetime is reported to be
exponentially distributed, whereas higher order taxa, such as families or
genera, follow a broader distribution, compatible with power law decay. We show
that both these evidences are consistent with a simple evolutionary model that
does not require specific assumptions on species interaction. The model
provides a zero-order description of the dynamics of ecological communities and
its species lifetime distribution can be computed exactly. Different behaviors
are found: an initial $t^{-3/2}$ power law, emerging from a random walk type of
dynamics, which crosses over to a steeper $t^{-2}$ branching process-like
regime and finally is cutoff by an exponential decay which becomes weaker and
weaker as the total population increases. Sampling effects can also be taken
into account and shown to be relevant: if species in the fossil record were
sampled according to the Fisher log-series distribution, lifetime should be
distributed according to a $t^{-1}$ power law. Such variability of behaviors in
a simple model, combined with the scarcity of data available, cast serious
doubts on the possibility to validate theories of evolution on the basis of
species lifetime data.


This paper was withdrawn by the authors.


  A set of rules is defined to systematically number the groups and the atoms
of organic molecules and, particularly, of polypeptides in a modular manner.
Supported by this numeration, a set of internal coordinates is defined. These
coordinates (termed Systematic, Approximately Separable and Modular Internal
Coordinates, SASMIC) are straightforwardly written in Z-matrix form and may be
directly implemented in typical Quantum Chemistry packages. A number of Perl
scripts that automatically generate the Z-matrix files for polypeptides are
provided as supplementary material. The main difference with other
Z-matrix-like coordinates normally used in the literature is that normal
dihedral angles (``principal dihedrals'' in this work) are only used to fix the
orientation of whole groups and a somewhat non-standard type of dihedrals,
termed ``phase dihedrals'', are used to describe the covalent structure inside
the groups. This physical approach allows to approximately separate soft and
hard movements of the molecule using only topological information and to
directly implement constraints. As an application, we use the coordinates
defined and ab initio quantum mechanical calculations to assess the commonly
assumed approximation of the free energy, obtained from ``integrating out'' the
side chain degree of freedom chi, by the Potential Energy Surface (PES) in the
protected dipeptide HCO-L-Ala-NH2. We also present a sub-box of the Hessian
matrix in two different sets of coordinates to illustrate the approximate
separation of soft and hard movements when the coordinates defined in this work
are used.


  It has long been hypothesized that futile cycles in cellular metabolism are
involved in the regulation of biochemical pathways. Following the work of
Newsholme and Crabtree, we develop a quantitative theory for this idea based on
open-system thermodynamics and metabolic control analysis. It is shown that the
{\it stoichiometric sensitivity} of an intermediary metabolite concentration
with respect to changes in steady-state flux is governed by the effective
equilibrium constant of the intermediate formation, and the equilibrium can be
regulated by a futile cycle. The direction of the shift in the effective
equilibrium constant depends on the direction of operation of the futile cycle.
High stoichiometric sensitivity corresponds to ultrasensitivity of an
intermediate concentration to net flow through a pathway; low stoichiometric
sensitivity corresponds to super-robustness of concentration with respect to
changes in flux. Both cases potentially play important roles in metabolic
regulation. Futile cycles actively shift the effective equilibrium by expending
energy; the magnitude of changes in effective equilibria and sensitivities is a
function of the amount of energy used by a futile cycle. This proposed
mechanism for control by futile cycles works remarkably similarly to kinetic
proofreading in biosynthesis. The sensitivity of the system is also intimately
related to the rate of concentration fluctuations of intermediate metabolites.
The possibly different roles of the two major mechanisms for cellular
biochemical regulation, namely reversible chemical modifications via futile
cycles and shifting equilibrium by macromolecular binding, are discussed.


  We develop a class of models with which we simulate the assembly of particles
into T1 capsid-like objects using Newtonian dynamics. By simulating assembly
for many different values of system parameters, we vary the forces that drive
assembly. For some ranges of parameters, assembly is facile, while for others,
assembly is dynamically frustrated by kinetic traps corresponding to malformed
or incompletely formed capsids. Our simulations sample many independent
trajectories at various capsomer concentrations, allowing for statistically
meaningful conclusions. Depending on subunit (i.e., capsomer) geometries,
successful assembly proceeds by several mechanisms involving binding of
intermediates of various sizes. We discuss the relationship between these
mechanisms and experimental evaluations of capsid assembly processes.


  It is increasingly evident that there is more to biological evolution than
natural selection; moreover, the concept of evolution is not limited to
biology. We propose an integrative framework for characterizing how entities
evolve, in which evolution is viewed as a process of context-driven
actualization of potential (CAP). Processes of change differ according to the
degree of nondeterminism, and the degree to which they are sensitive to,
internalize, and depend upon a particular context. The approach enables us to
embed phenomena across disciplines into a broad conceptual framework. We give
examples of insights into physics, biology, culture and cognition that derive
from this unifying framework.


  The conditions of the chromosomes inside the nucleus in the Rabl
configuration have been modelled as self-avoiding polymer chains under
restraining conditions. To ensure that the chromosomes remain stretched out and
lined up, we fixed their end points to two opposing walls. The numbers of
segments $N$, the distances $d_1$ and $d_2$ between the fixpoints, and the
wall-to-wall distance $z$ (as measured in segment lengths) determine an
approximate value for the Kuhn segment length $k_l$. We have simulated the
movement of the chromosomes using molecular dynamics to obtain the expected
distance distribution between the genetic loci in the absence of further
attractive or repulsive forces. A comparison to biological experiments on
\textit{Drosophila Melanogaster} yields information on the parameters for our
model. With the correct parameters it is possible to draw conclusions on the
strength and range of the attraction that leads to pairing.


  For a model of molecular evolution to be useful for phylogenetic inference,
the topology of evolutionary trees must be identifiable. That is, from a joint
distribution the model predicts, it must be possible to recover the tree
parameter. We establish tree identifiability for a number of phylogenetic
models, including a covarion model and a variety of mixture models with a
limited number of classes. The proof is based on the introduction of a more
general model, allowing more states at internal nodes of the tree than at
leaves, and the study of the algebraic variety formed by the joint
distributions to which it gives rise. Tree identifiability is first established
for this general model through the use of certain phylogenetic invariants.


  Biological genomes are divided into coding and non-coding regions. Introns
are non-coding parts within genes, while the remaining non-coding parts are
intergenic sequences. To study the evolutionary significance of recombination
inside introns we have used two models based on the Monte Carlo method. In our
computer simulations we have implemented the internal structure of genes by
declaring the probability of recombination between exons. One situation when
inside-intron recombination is advantageous is recovering functional genes by
combining proper exons dispersed in the genetic pool of the population after a
long period without selection for the function of the gene. Populations have to
pass through the bottleneck, then. These events are rather rare and we have
expected that there should be other phenomena giving profits from the
inside-intron recombination. In fact we have found that inside-intron
recombination is advantageous only in the case when after recombination,
besides the recombinant forms, parental haplotypes are available and selection
is set already on gametes.


  Integrating biological information from different sources to understand
cellular processes is an important problem in systems biology. We use data from
mRNA expression arrays and chemical kinetics to formulate a metabolic model
relevant to K562 erythroleukemia cells. MAP kinase pathway activation alters
the expression of metabolic enzymes in K562 cells. Our array data show changes
in expression of lactate dehydrogenase (LDH) isoforms after treatment with
phorbol 12-myristate 13-acetate (PMA), which activates MAP kinase signaling. We
model the change in lactate production which occurs when the MAP kinase pathway
is activated, using a non-equilibrium, chemical-kinetic model of homolactic
fermentation. In particular, we examine the role of LDH isoforms, which
catalyze the conversion of pyruvate to lactate. Changes in the isoform ratio
are not the primary determinant of the production of lactate. Rather, the total
concentration of LDH controls the lactate concentration.


  We show that the hydrophobicity of sequences is the leading term in
Miyazawa-Jernigan interactions. Being the source of additive (solvation) terms
in pair-contact interactions, they were used to reduce the energy parameters
while resulting in a clear vector manipulation of energy. The reduced
(additive) potential performs considerably successful in predicting the
statistical properties of arbitrary structures. The evaluated designabilities
of the structures by both models are highly correlated. Suggesting
geometrically non-degenerate vectors (structures) as protein-like structures,
the additive model is a powerful tool for protein design. Moreover, a crossing
point in the log-linear diagram of designability-ranking shows that about 1/e
of the structures have designabilities above the average, independent on the
used model.


  On the basis of epidemiological studies, infection was suggested to play a
role in the etiology of human cancer. While for some cancers such a role was
indeed demonstrated, there is no direct biological support for the role of
viral pathogens in the pathogenesis of childhood leukemia. Using a novel
bioinformatic tool, that alternates between clustering and standard statistical
methods of analysis, we performed a "double blind" search of published gene
expression data of subjects with different childhood ALL subtypes, looking for
unanticipated partitions of patients, induced by unexpected groups of genes
with correlated expression. We discovered a group of about thirty genes,
related to the interferon response pathway, whose expression levels divide the
ALL samples into two subgroups; high in 50, low in 285 patients. Leukemic
subclasses prevalent in early childhood (the age most susceptible to infection)
are over-represented in the high expression subgroup. Similar partitions,
induced by the same genes, were found also in breast and ovarian cancer but not
in lung cancer, prostate cancer and lymphoma. About 40% of breast cancer
samples expressed the "interferon- related" signature. It is of interested that
several studies demonstrated MMTV-like sequences in about 40% of breast cancer
samples. Our discovery of an unanticipated strong signature of an interferon
induced pathway provides molecular support for a role for either inflammation
or viral infection in the pathogenesis of childhood leukemia as well as breast
and ovarian cancer.


  r-selection refers to evolutionary competition in the rate of a population's
exponential increase. This is contrasted with K-selection, in which populations
in steady-state compete in efficiency of resource conversion. Evolution in
nature is thought to combine these two in various proportions. But in modeling
the evolution of life histories, theorists have used r-selection exclusively;
up until now, there has not been a practical algorithm for computing the target
function of K-selection. The Malthusian parameter, as computed from the
Euler-Lotka equation, is a quantitative rendering of the r in r-selection,
computed from the fundamental life history variables mortality and fertility.
Herein, a quantitative formulation of K is derived in similar terms. The basis
for our model is the logistic equation which, we argue, applies more generally
than is commonly appreciated. Support is offered for the utility of this
paradigm, and one example computation is exhibited, in which K-selection
appears to support pleiotropic explanations for senescence only one fourth as
well as r-selection.


  Deciphering regulatory events that drive malignant transformation represents
a major challenge for systems biology. Here we analyzed genome-wide
transcription profiling of an in-vitro transformation process. We focused on a
cluster of genes whose expression levels increased as a function of p53 and
p16INK4A tumor suppressors inactivation. This cluster predominantly consists of
cell cycle genes and constitutes a signature of a diversity of cancers. By
linking expression profiles of the genes in the cluster with the dynamic
behavior of p53 and p16INK4A, we identified a promoter architecture that
integrates signals from the two tumor suppressive channels and that maps their
activity onto distinct levels of expression of the cell cycle genes, which in
turn, correspond to different cellular proliferation rates. Taking components
of the mitotic spindle as an example, we experimentally verified our
predictions that p53-mediated transcriptional repression of several of these
novel targets is dependent on the activities of p21, NFY and E2F. Our study
demonstrates how a well-controlled transformation process allows linking
between gene expression, promoter architecture and activity of upstream
signaling molecules.


  A generalized model of the spread of the Hantavirus in mice populations is
presented on the basis of recent observational findings concerning the movement
characteristics of the mice that carry the infection. The factual information
behind the generalization is based on mark-recapture observations reported in
Giuggioli et al. [Bull. Math. Biol. 67, 1135 (2005)] that have necessitated the
introduction of home ranges in the simple model of Hantavirus spread presented
by Abramson and Kenkre [Phys. Rev. E 66, 11912 (2002)]. The essential feature
of the model presented here is the existence of adult mice that remain largely
confined to locations near their home ranges, and itinerant juvenile mice that
are not so confined, and, during their search for their own homes, move and
infect both other juveniles and adults that they meet during their movement.
The model is presented at three levels of description: mean field, kinetic and
configuration. Results of calculations are shown explicitly from the mean field
equations and the simulation rules, and are found to agree in some respects and
to differ in others. The origin of the differences is shown to lie in spatial
correlations. It is indicated how mark-recapture observations in the field may
be employed to verify the applicability of the theory.


  Simple random walk considerations are used to interpret rodent population
data collected in Hantavirus-related investigations in Panama regarding the
short-tailed cane mouse, \emph{Zygodontomys brevicauda}. The diffusion constant
of mice is evaluated to be of the order of (and larger than) 200 meters squared
per day. The investigation also shows that the rodent mean square displacement
saturates in time, indicating the existence of a spatial scale which could, in
principle, be the home range of the rodents. This home range is concluded to be
of the order of 70 meters. Theoretical analysis is provided for interpreting
animal movement data in terms of an interplay of the home ranges, the diffusion
constant, and the size of the grid used to monitor the movement. The study
gives impetus to a substantial modification of existing theory of the spread of
the Hantavirus epidemic which has been based on simple diffusive motion of the
rodents, and additionally emphasizes the importance for developing more
accurate techniques for the measurement of rodent movement.


  A theory is provided for the estimation of home ranges of animals from the
standard mark-recapture technique in which data are collected by capturing,
tagging and recapturing the animals. The theoretical tool used is the
Fokker-Planck equation, its characteristic quantities being the diffusion
constant which describes the motion of the animals, and the attractive
potential which addresses their tendency to live near their burrows. The
measurement technique is shown to correspond to the calculation of a certain
kind of mean square displacement of the animals relevant to the specific
probing window in space corresponding to the trapping region. The output of the
theory is a sigmoid curve of the observable mean square displacement as a
function of the ratio of distances characteristic of the home range and the
trapping region, along with an explicit prescription to extract the home range
form observations. Applications of the theory to rodent movement in Panama and
New Mexico are pointed out. An analysis is given of the sensitivity of our
theory to the choice of the confining potential via the use of various
representative cases. A comparison is provided between home range size inferred
from our method and from other procedures employed in the literature.
Consequences of home range overlap are also discussed.


  We determine the persistence length, $l_p$, for a bacterial group I ribozyme
as a function of concentration of monovalent and divalent cations by fitting
the distance distribution functions $P(r)$ obtained from small angle X-ray
scattering intensity data to the asymptotic form of the calculated $P_{WLC}(r)$
for a worm-like chain (WLC). The $l_p$ values change dramatically over a narrow
range of \Mg concentration from $\sim$21 \AA in the unfolded state (\textbf{U})
to $\sim$10 \AA in the compact ($\mathrm{I_C}$) and native states. Variations
in $l_p$ with increasing \Na concentration are more gradual. In accord with the
predictions of polyelectrolyte theory we find $l_p \propto 1/ \kappa^2$ where
$\kappa$ is the inverse Debye-screening length.


  We analyze data from a long term field project in New Mexico, consisting of
repeated sessions of mark-recaptures of Peromyscus maniculatus (Rodentia:
Muridae), the host and reservoir of Sin Nombre Virus (Bunyaviridae:
Hantavirus). The displacements of the recaptured animals provide a means to
study their movement from a statistical point of view. We extract two
parameters from the data with the help of a simple model: the diffusion
constant of the rodents, and the size of their home range. The short time
behavior shows the motion to be approximately diffusive and the diffusion
constant to be 470+/-50m^2/day. The long time behavior provides an estimation
of the diameter of the rodent home ranges, with an average value of 100+/-25m.
As in previous investigations directed at Zygodontomys brevicauda observations
in Panama, we use a box model for home range estimation. We also use a harmonic
model in the present investigation to study the sensitivity of the conclusions
to the model used and find that both models lead to similar estimates.


  We present a novel framework for inferring regulatory and sequence-level
information from gene co-expression networks. The key idea of our methodology
is the systematic integration of network inference and network topological
analysis approaches for uncovering biological insights. We determine the gene
co-expression network of Bacillus subtilis using Affymetrix GeneChip time
series data and show how the inferred network topology can be linked to
sequence-level information hard-wired in the organism's genome. We propose a
systematic way for determining the correlation threshold at which two genes are
assessed to be co-expressed by using the clustering coefficient and we expand
the scope of the gene co-expression network by proposing the slope ratio metric
as a means for incorporating directionality on the edges. We show through
specific examples for B. subtilis that by incorporating expression level
information in addition to the temporal expression patterns, we can uncover
sequence-level biological insights. In particular, we are able to identify a
number of cases where (i) the co-expressed genes are part of a single
transcriptional unit or operon and (ii) the inferred directionality arises due
to the presence of intra-operon transcription termination sites.


  In this report, correlation of the pixels comprising a microarray spot is
investigated. Subsequently, correlation statistics namely: Pearson correlation
and Spearman rank correlation are used to segment the foreground and background
intensity of microarray spots. The performance of correlation-based
segmentation is compared to clustering-based (PAM, k-means) and seeded-region
growing techniques (SPOT). It is shown that correlation-based segmentation is
useful in flagging poorly hybridized spots, thus minimizes false-positives. The
present study also raises the intriguing question of whether a change in
correlation can be an indicator of differential gene expression.


  We present some results of simulations of population growth and evolution,
using the standard asexual Penna model, with individuals characterized by a
string of bits representing a genome containing some possible mutations. After
about 20000 simulation steps, when only a few genetic families are still
present from among rich variety of families at the beginning of the simulation
game, strong peaks in mutation distribution functions are observed. This known
effect is due to evolution rules with hereditary mechanism. The birth and death
balance in the simulation game also leads to elimination of families specified
by different genomes. Number of families $G(t)$ versus time $t$ follow the
power law, $G \propto t^n$. Our results show the power coefficient exponent $n$
is changing as the time goes. Starting from about --1, smoothly achieves about
--2 after hundreds of steps, and finally has semi-smooth transition to 0, when
only one family exists in the environment. This is in contrast with constant
$n$ about --1 as found, for example, in \cite{bib:evolution}. We suspect that
this discrepancy may be due to two different time scales in simulations -
initial stages follow the $n\approx-1$ law, yet for large number of simulation
steps we get $n\approx-2$, providing random initial population was sufficiently
big to allow for still reliable statistical analysis. The $n\approx-1$
evolution stage seems to be associated with the Verhulst mechanism of
population elimination due to the limited environmental capacity - when the
standard evolution rules were modified, we observed a plateau ($n=0$) in the
power law in short time scale, again followed by $n\approx -2$ law for longer
times. The modified model uses birth rate controlled by the current population
instead of the standard Verhulst death factor.


  We study secondary structures of random RNA molecules by means of a
renormalized field theory based on an expansion in the sequence disorder. We
show that there is a continuous phase transition from a molten phase at higher
temperatures to a low-temperature glass phase. The primary freezing occurs
above the critical temperature, with local islands of stable folds forming
within the molten phase. The size of these islands defines the correlation
length of the transition. Our results include critical exponents at the
transition and in the glass phase.


  In this paper the heat transport in microtubules (MT) is investigated. When
the dimension of the structure is of the order of the de Broglie wave length
the transport phenomena must be analyzed within quantum mechanics. In this
paper we developed the Dirac type thermal equation for MT .The solution of the
equation-the temperature fields for electrons can be wave type or diffusion
type depending on the dynamics of the scattering. Key words: Microtubules
ultrashort laser pulses, Dirac thermal equation, temperature fields.


  Mouse prion protein PrP106-126 is a peptide corresponding to the residues
107-127 of human prion protein. It has been shown that PrP106-126 can reproduce
the main neuropathological features of prionrelated transmissible spongiform
encephalopathies and can form amyloid-like fibrils in vitro. The conformational
characteristics of PrP106-126 fibril have been investigated by electron
microscopy, CD spectroscopy, NMR and molecular dynamics simulations. Recent
researches have found out that PrP106-126 in water assumes a stable structure
consisting of two parallel beta-sheets that are tightly packed against each
other. In this work we perform molecular dynamics simulation to reveal the
elongation mechanism of PrP106-126 fibril. Influenced by the edge strands of
the fibril which already adopt beta-sheets conformation, single PrP106-126
peptide forms beta-structure and becomes a new element of the fibril. Under
acidic condition, single PrP106-126 peptide adopts a much larger variety of
conformations than it does under neural condition, which makes a peptide easier
to be influenced by the edge strands of the fibril. However, acidic condition
dose not largely affect the stability of PrP106-126 peptide fibril. Thus, the
speed of fibril elongation can be dramatically increased by lowering the pH
value of the solution. The pH value was adjusted by either changing the
protonation state of the residues or adding hydronium ions (acidic solution) or
hydroxyl ions (alkaline solution). The differences between these two approaches
are analyzed here.


  Spontaneous switching events in most characterized genetic switches are rare,
resulting in extremely stable epigenetic properties. We show how simple
arguments lead to theories of the rate of such events much like the absolute
rate theory of chemical reactions corrected by a transmission factor. Both the
probability of the rare cellular states that allow epigenetic escape, and the
transmission factor, depend on the rates of DNA binding and unbinding events
and on the rates of protein synthesis and degradation. Different mechanisms of
escape from the stable attractors occur in the nonadiabatic, weakly adiabatic
and strictly adiabatic regimes, characterized by the relative values of those
input rates.


  We consider the directed evolution of a population after an intervention that
has significantly altered the underlying fitness landscape. We model the space
of genotypes as a distributive lattice; the fitness landscape is a real-valued
function on that lattice. The risk of escape from intervention, i.e., the
probability that the population develops an escape mutant before extinction, is
encoded in the risk polynomial. Tools from algebraic combinatorics are applied
to compute the risk polynomial in terms of the fitness landscape. In an
application to the development of drug resistance in HIV, we study the risk of
viral escape from treatment with the protease inhibitors ritonavir and
indinavir.


  Precisely how ecological factors influence animal social structure is far
from clear. We explore this question using an agent-based model inspired by the
fission-fusion society of Spider monkeys (Ateles spp). Our model introduces a
realistic, complex foraging environment composed of many resource patches with
size varying according to an inverse power-law frequency distribution with
exponent $\beta$. Foragers do not interact among them and start from random
initial locations. They have either a complete or partial knowledge of the
environment and maximize the ratio between the size of the next visited patch
and the distance traveled to it, ignoring previously visited patches. At
intermediate values of $\beta$, when large patches are neither too scarce nor
too abundant, foragers form groups (coincide at the same patch) with a similar
size frequency distribution as the spider monkey's subgroups. Fission-fusion
events create a network of associations that contains weak bonds among foragers
that meet only rarely and strong bonds among those that repeat associations
more frequently than would be expected by chance. The latter form sub-networks
with the highest number of bonds and a high clustering coefficient at
intermediate values of $\beta$. The weak bonds enable the whole social network
to percolate. Some of our results are similar to those found in long-term field
studies of spider monkeys and other fission-fusion species. We conclude that
hypotheses about the ecological causes of fission-fusion and the origin of
complex social structures should consider the heterogeneity and complexity of
the environment in which social animals live.


  Cellular life requires the presence of a set of biochemical mechanisms in
order to maintain a predictable process of growth and division. Several
attempts have been made towards the building of minimal protocells from a
top-down approach, i.e. by using available biomolecules This type of synthetic
approach has so far been unsuccesful, and the lack of appropriate models of the
synthetic protocell cycle might be needed to guide future experiments. In this
paper we present a simple biochemically and physically feasible model of cell
replication involving a discrete semi-permeable vesicle with an internal
minimal metabolism involving two reactive centers. It is shown that such a
system can effectively undergo a whole cell replication cycle. The model can be
used as a basic framework to model whole protocell dynamics including more
complex sets of reactions. The possible implementation of our design in future
synthetic protocells is outlined.


  At an early stage in pre-biotic evolution, groups of replicating molecules
must coordinate their reproduction to form aggregated units of selection.
Mechanisms that enable this to occur are currently not well understood. In this
paper we introduce a deterministic model of primitive replicating aggregates,
proto-organisms, that host populations of replicating information carrying
molecules. Some of the molecules promote the reproduction of the proto-organism
at the cost of their individual replication rate. A situation resembling that
of group selection arises. We derive and analytically solve a partial
differential equation that describes the system. We find that the relative
prevalence of fast and slow replicators is determined by the relative strength
of selection at the aggregate level to the selection strength at the molecular
level. The analysis is concluded by a preliminary treatment of finite
population size effects.


  The relationship between the design and functionality of molecular networks
is now a key issue in biology. Comparison of regulatory networks performing
similar tasks can give insights into how network architecture is constrained by
the functions it directs. We here discuss methods of network comparison based
on network architecture and signaling logic. Introducing local and global
signaling scores for the difference between two networks we quantify
similarities between evolutionary closely and distantly related bacteriophages.
Despite the large evolutionary separation between phage $\lambda$ and 186 their
networks are found to be similar when difference is measured in terms of global
signaling. We finally discuss how network alignment can be used to to pinpoint
protein similarities viewed from the network perspective.


  The Penna model is a model of evolutionary ageing through mutation
accumulation where traditionally time and the age of an organism are treated as
discrete variables and an organism's genome by a binary bit string. We
reformulate the asexual Penna model and show that, a universal scale invariance
emerges as we increase the number of discrete genome bits to the limit of a
continuum. The continuum model, introduced by Almeida and Thomas in
[Int.J.Mod.Phys.C, 11, 1209 (2000)] can be recovered from the discrete model in
the limit of infinite bits coupled with a vanishing mutation rate per bit.
Finally, we show that scale invariant properties may lead to the ubiquitous
Gompertz Law for mortality rates for early ages, which is generally regarded as
being empirical.


  The Min system in Escherichia coli directs division to the centre of the cell
through pole-to-pole oscillations of the MinCDE proteins. We present a one
dimensional stochastic model of these oscillations which incorporates membrane
polymerisation of MinD into linear chains. This model reproduces much of the
observed phenomenology of the Min system, including pole-to-pole oscillations
of the Min proteins. We then apply this model to investigate the Min system
during cell division. Oscillations continue initially unaffected by the closing
septum, before cutting off rapidly. The fractions of Min proteins in the
daughter cells vary widely, from 50%-50% up to 85%-15% of the total from the
parent cell, suggesting that there may be another mechanism for regulating
these levels in vivo.


  Recombination is introduced into Eigen's theory of quasispecies evolution.
Comparing numerical simulations of the rate equations in the non-recombining
and recombining cases show that recombination has a strong effect on the error
threshold and, for a wide range of mutation rates, gives rise to two stable
fixed points in the dynamics. This bi-stability results in the existence of two
error thresholds. We prove that, under some assumptions on the fitness
landscape but for general crossover probability, a fixed point localized about
the sequence with superior fitness is globally stable for low mutation rates.


  SELEX (Systematic Evolution of Ligands by Exponential Enrichment) is an
experimental procedure that allows extracting, from an initially random pool of
DNA, those oligomers with high affinity for a given DNA-binding protein. We
address what is a suitable experimental and computational procedure to infer
parameters of transcription factor-DNA interaction from SELEX experiments. To
answer this, we use a biophysical model of transcription factor-DNA
interactions to quantitatively model SELEX. We show that a standard procedure
is unsuitable for obtaining accurate interaction parameters. However, we
theoretically show that a modified experiment in which chemical potential is
fixed through different rounds of the experiment allows robust generation of an
appropriate data set. Based on our quantitative model, we propose a novel
bioinformatic method of data analysis for such modified experiment and apply it
to extract the interaction parameters for a mammalian transcription factor
CTF/NFI. From a practical point of view, our method results in a significantly
improved false positive/false negative trade-off, as compared to both the
standard information theory based method and a widely used empirically
formulated procedure.


  We use a combination of analytic models and computer simulations to gain
insight into the dynamics of evolution. Our results suggest that certain
interesting phenomena should eventually emerge from the fossil record. For
example, there should be a ``tortoise and hare effect'': Those genera with the
smallest species death rate are likely to survive much longer than genera with
large species birth and death rates. A complete characterization of the
behavior of a branch of the phylogenetic tree corresponding to a genus and
accurate mathematical representations of the various stages are obtained. We
apply our results to address certain controversial issues that have arisen in
paleontology such as the importance of punctuated equilibrium and whether
unique Cambrian phyla have survived to the present.


  Some microbial organisms are known to randomly slip into and out of
hibernation, irrespective of environmental conditions [1]. In a (genetically)
uniform population a typically very small subpopulation becomes metabolically
inactive whereas the majority subpopulation remains active and grows. Bacteria
such as E. coli, Staphylococcus aureus (MRSA-superbug), Mycobacterium
tuberculosis, and Pseudomonas aeruginosa [1-3] show persistence. It can render
bacteria less vulnerable in adverse environments [1, 4, 5] and their effective
eradication through medication more difficult [2, 3, 6]. Here we show that
medication treatment regimes may have to be modified when persistence is taken
into account and characterize optimal approaches assuming that the total
medication dose is constrained. The determining factors are cumulative
toxicity, eradication power of the medication and bacterial response
timescales. Persistent organisms have to be fought using tailored eradication
strategies which display two fundamental characteristics. Ideally, the
treatment time should be significantly longer than in the case of persistence
with the medication uniformly spread out over time; however, if treatment time
has to be limited, then the application of medication has to be concentrated
towards the beginning and end of the treatment. These findings deviate from
current clinical practice, and may therefore help to optimize and simplify
treatments. Our use of multi-objective optimization [7] to map out the optimal
strategies can be generalized to other related problems.


  Mechanically induced protein unfolding in the force-clamp apparatus is shown,
in a coarse-grained model of ubiquitin, to have lognormal statistics above a
treshold force and exponential below it. Correspondingly, the mean unfolding
time is slowly varying and exponentially decreases as a function of the force.
The time dependencies of the end-to-end distances are also distinct. The time
sequence of unfolding events weakly depends on force and much of it resembles
that for stretching at constant speed. A more complicated time dependence
arises for integrin.


  Evidence is presented suggesting, for the first time, that the protein
foldability metric sigma=(T_theta - T_f)/T_theta, where T_theta and T_f are,
respectively, the collapse and folding transition temperatures, could be used
also to measure the foldability of RNA sequences. The importance of sigma is
discussed in the context of the in silico design of rapidly foldable RNA
sequences.


  We study a spatially inhomogeneous model of cancer growth based on
Michaelis--Menten kinetics, subjected to additive Gaussian noise and
multiplicative dichotomous noise. In presence of the latter, we can observe a
transition between two stationary states of the system. The transient behaviour
generates a spatial pattern of two phases, where cancer cells or immune cells
predominate.


  Phage lambda is one of the most studied biological models in modern molecular
biology. Over the past 50 years quantitative experimental knowledge on this
biological model has been accumulated at all levels: physics, chemistry,
genomics, proteomics, functions, and more. All its components have been known
to a great detail. The theoretical task has been to integrate its components to
make the organism working quantitatively in a harmonic manner. This would test
our biological understanding and would lay a solid fundamental for further
explorations and applications, an obvious goal of systems biology. One of the
outstanding challenges in doing so has been the so-called stability puzzle of
lambda switch: the biologically observed robustness and its difficult
mathematical reconstruction based on known experimental values. In this chapter
we review the recent theoretical and experimental efforts on tackling this
problem. An emphasis is put on the minimum quantitative modeling where a
successful numerical agreement between experiments and modeling has been
achieved. A novel method tentatively named stochastic dynamical structure
analysis emerged from such study is also discussed within a broad modeling
perspective.


  The classic algorithms of Needleman--Wunsch and Smith--Waterman find a
maximum a posteriori probability alignment for a pair hidden Markov model
(PHMM). In order to process large genomes that have undergone complex genome
rearrangements, almost all existing whole genome alignment methods apply fast
heuristics to divide genomes into small pieces which are suitable for
Needleman--Wunsch alignment. In these alignment methods, it is standard
practice to fix the parameters and to produce a single alignment for subsequent
analysis by biologists.
  Our main result is the construction of a whole genome parametric alignment of
Drosophila melanogaster and Drosophila pseudoobscura. Parametric alignment
resolves the issue of robustness to changes in parameters by finding all
optimal alignments for all possible parameters in a PHMM. Our alignment draws
on existing heuristics for dividing whole genomes into small pieces for
alignment, and it relies on advances we have made in computing convex polytopes
that allow us to parametrically align non-coding regions using biologically
realistic models. We demonstrate the utility of our parametric alignment for
biological inference by showing that cis-regulatory elements are more conserved
between Drosophila melanogaster and Drosophila pseudoobscura than previously
thought. We also show how whole genome parametric alignment can be used to
quantitatively assess the dependence of branch length estimates on alignment
parameters.
  The alignment polytopes, software, and supplementary material can be
downloaded at http://bio.math.berkeley.edu/parametric/.


  This article presents a new way to understand the descriptive ability of tree
shape statistics. Where before tree shape statistics were chosen by their
ability to distinguish between macroevolutionary models, the ``resolution''
presented in this paper quantifies the ability of a statistic to differentiate
between similar and different trees. We term this a ``geometric'' approach to
differentiate it from the model-based approach previously explored. A distinct
advantage of this perspective is that it allows evaluation of multiple tree
shape statistics describing different aspects of tree shape. After developing
the methodology, it is applied here to make specific recommendations for a
suite of three statistics which will hopefully prove useful in applications.
The article ends with an application of the tree shape statistics to clarify
the impact of omission of taxa on tree shape.


  There are several common ways to encode a tree as a matrix, such as the
adjacency matrix, the Laplacian matrix (that is, the infinitesimal generator of
the natural random walk), and the matrix of pairwise distances between leaves.
Such representations involve a specific labeling of the vertices or at least
the leaves, and so it is natural to attempt to identify trees by some feature
of the associated matrices that is invariant under relabeling. An obvious
candidate is the spectrum of eigenvalues (or, equivalently, the characteristic
polynomial). We show for any of these choices of matrix that the fraction of
binary trees with a unique spectrum goes to zero as the number of leaves goes
to infinity. We investigate the rate of convergence of the above fraction to
zero using numerical methods. For the adjacency and Laplacian matrices, we show
that that the {\em a priori} more informative immanantal polynomials have no
greater power to distinguish between trees.


  The transport of cargo particles which are pulled by several molecular motors
in a cooperative manner is studied theoretically. The transport properties
depend primarily on the maximal number, $N$, of motor molecules that may pull
simultaneously on the cargo particle. Since each motor must unbind from the
filament after a finite number of steps but can also rebind to it again, the
actual number of pulling motors is not constant but varies with time between
zero and $N$. An increase in the maximal number $N$ leads to a strong increase
of the average walking distance (or run length) of the cargo particle. If the
cargo is pulled by up to $N$ kinesin motors, e.g., the walking distance is
estimated to be $5^{N-1}/N$ micrometers which implies that seven or eight
kinesin molecules are sufficient to attain an average walking distance in the
centimeter range. If the cargo particle is pulled against an external load
force, this force is shared between the motors which provides a nontrivial
motor-motor coupling and a generic mechanism for nonlinear force-velocity
relationships. With increasing load force, the probability distribution of the
instantenous velocity is shifted towards smaller values, becomes broader, and
develops several peaks. Our theory is consistent with available experimental
data and makes quantitative predictions that are accessible to systematic in
vitro experiments.


  For two genotypes that have the same mean number of offspring but differ in
the variance in offspring number, natural selection will favor the genotype
with lower variance. The concept of fitness becomes cloudy under these
conditions because the outcome of evolution is not deterministic. However, the
effect of variance in offspring number on the fixation probability of mutant
strategies has been calculated under several scenarios with the general
conclusion that variance in offspring number reduces fitness but only in
proportion to the inverse of the population size (Gillespie 1974, Proulx 2000).
This relationship becomes more complicated under a metapopulation scenario
where the "effective" population size depends on migration rate, population
structure, and life cycle. We show that under hard selection and weak migration
fitness in a metapopulation composed of equal sized demes is determined by deme
size. Conversely, for high migration rates and hard selection the effective
fitness depends on the total size of the metapopulation. Interestingly, under
soft selection there is no effect of migration or neighboring population
structure on effective fitness, and fitness depends only on deme size. We use
individual based simulations in developed in Shpak (2005) to validate our
analytical approximations and investigate deviations of our assumption of equal
deme size.


  We model a system of networking agents that seek to optimize their centrality
in the network while keeping their cost, the number of connections they are
participating in, low. Unlike other game-theory based models for network
evolution, the success of the agents is related only to their position in the
network. The agents use strategies based on local information to improve their
chance of success. Both the evolution of strategies and network structure are
investigated. We find a dramatic time evolution with cascades of strategy
change accompanied by a change in network structure. On average the network
self-organizes to a state close to the transition between a fragmented state
and a state with a giant component. Furthermore, with increasing system size
both the average degree and the level of fragmentation decreases. We also
observe that the network keeps on actively evolving, although it does not have
to, thus suggesting a Red Queen-like situation where agents have to keep on
networking and responding to the moves of the others in order to stay
successful.


  The collective dynamics of $N$ weakly coupled processive molecular motors are
considered theoretically. We show, using a discrete lattice model, that the
velocity-force curves strongly depend on the effective dynamic interactions
between motors and differ significantly from a simple mean field prediction.
They become essentially independent of $N$ if it is large enough. For strongly
biased motors such as kinesin this occurs if $N\gtrsim 5$. The study of a
two-state model shows that the existence of internal states can induce
effective interactions.


  Recombination is an important event in the evolution of HIV. It affects the
global spread of the pandemic as well as evolutionary escape from host immune
response and from drug therapy within single patients. Comprehensive
computational methods are needed for detecting recombinant sequences in large
databases, and for inferring the parental sequences.
  We present a hidden Markov model to annotate a query sequence as a
recombinant of a given set of aligned sequences. Parametric inference is used
to determine all optimal annotations for all parameters of the model. We show
that the inferred annotations recover most features of established hand-curated
annotations. Thus, parametric analysis of the hidden Markov model is feasible
for HIV full-length genomes, and it improves the detection and annotation of
recombinant forms.
  All computational results, reference alignments, and C++ source code are
available at http://bio.math.berkeley.edu/recombination/.


  In order to elucidate the role of the native state topology and the stability
of subdomains in protein folding, we investigate free energy landscape of human
lysozyme, which is composed of two subdomains, by Monte Carlo simulations. A
realistic lattice model with Go-like interaction is used. We take the relative
interaction strength (stability, in other word) of two subdomains as a variable
parameter and study the folding process. A variety of folding process is
observed and we obtained a phase diagram of folding in terms of temperature and
the relative stability. Experimentally-observed diversity in folding process of
c-type lysozimes is thus understood as a consequence of the difference in the
relative stability of subdomains.


  A complete parametric analysis of dynamic regimes of a conceptual model of
anti-tumor virus therapy is presented. The role and limitations of mass-action
kinetics are discussed. A functional response, which is a function of the ratio
of uninfected to infected tumor cells, is proposed to describe the spread of
the virus infection in the tumor. One of the main mathematical features of
ratio-dependent models is that the origin is a complicated equilibrium point
whose characteristics crucially determine the main properties of the model. It
is shown that, in a certain area of parameter values, the trajectories of the
model form a family of homoclinics to the origin (so-called elliptic sector).
Biologically, this means that both infected and uninfected tumor cells can be
eliminated with time, and complete recovery is possible as a result of the
virus therapy within the framework of deterministic models


  Reserved territories seem to be best reference sites of wildnature, where the
long-term observations are carried out. Simulation model of spatially
distributed processes of contamination of the state reservation is developed,
and the dynamics of some pollutants is studied. An issue of the generalized
evaluation of an ecological system status is discussed.


  By considering the energy cost of messages carried by proteins as
proportional to their information content we found experimental proof that
proteins from all living organisms tend to have their estimated semantic
content of information per unit mass, statistically, close to a constant. Thus,
in the message carried by proteins -to achieve minimum energy waste- the rate
of information content per unit mass tends to be optimized in living organisms.
The experimental evidence of this new information law resembles a marathon
where highly optimized proteins correspond to advanced runners followed by a
main bunch and the stragglers -lowly optimized proteins. Our results suggest
the existence of a continuous optimization process that living organisms had to
face, in which a compromise between biological functionality, economic
feasibility and the survival requirements is established.


  During the eukaryotic cell cycle, chromatin undergoes several conformational
changes, which are believed to play key roles in gene expression regulation
during interphase, and in genome replication and division during mitosis. In
this paper, we propose a scenario for chromatin structural reorganization
during mitosis, which bridges all the different scales involved in chromatin
architecture, from nucleosomes to chromatin loops. We build a model for
chromatin, based on available data, taking into account both physical and
topological constraints DNA has to deal with. Our results suggest that the
mitotic chromosome condensation/decondensation process is induced by a
structural change at the level of the nucleosome itself.


  We cloned, expressed and purified the Escherichia coli yhbO gene product,
which is homolog to the Bacillus subtilis general stress protein 18 (the yfkM
gene product), the Pyrococcus furiosus intracellular protease PfpI, and the
human Parkinson disease protein DJ-1. The gene coding for YhbO was generated by
amplifying the yhbO gene from E. coli by polymerase chain reaction. It was
inserted in the expression plasmid pET-21a, under the transcriptional control
of the bacteriophage T7 promoter and lac operator. A BL21(DE3) E. coli strain
transformed with the YhbO-expression vector pET-21a-yhbO, accumulates large
amounts of a soluble protein of 20 kDa in SDS-PAGE that matches the expected
YhbO molecular weight. YhbO was purified to homogeneity by HPLC DEAE ion
exchange chromatography and hydroxylapatite chromatography and its identity was
confirmed by N-terminal sequencing and mass spectrometry analysis. The native
protein exists in monomeric, trimeric and hexameric forms.


  The thermodynamics of the small SH3 protein domain is studied by means of a
simplified model where each bead-like amino acid interacts with the others
through a contact potential controlled by a 20x20 random matrix. Good folding
sequences, characterized by a low native energy, display three main
thermodynamical phases, namely a coil-like phase, an unfolded globule and a
folded phase (plus other two phases, namely frozen and random coil, populated
only at extremes temperatures). Interestingly, the unfolded globule has some
regions already structured. Poorly designed sequences, on the other hand,
display a wide transition from the random coil to a frozen state. The
comparison with the analytic theory of heteropolymers is discussed.


  Effects of predators of juvenile mice on the spread of the Hantavirus are
analyzed in the context of a recently proposed model. Two critical values of
the predation probability are identified. When the smaller of them is exceeded,
the hantavirus infection vanishes without extinguishing the mice population.
When the larger is exceeded, the entire mice population vanishes. These results
suggest the possibility of control of the spread of the epidemic by introducing
predators in areas of mice colonies in a suitable way so that such control does
not kill all the mice but lowers the epidemic spread.


  Background: Approximately 40% of adults in Botswana are HIV-infected. The
Botswana antiretroviral program began in 2002 and currently treats 34,000
patients with a goal of treating 85,000 patients (~30% of HIV-infected adults)
by 2009. We predict the evolution of drug-resistant strains of HIV that may
emerge as a consequence of this treatment program. We discuss the implications
of our results for the World Health Organization's (WHO's) proposed
surveillance system for detecting drug-resistant strains of HIV in Africa.
Methods: We use a mathematical model of the emergence of drug resistance. We
incorporate demographic and treatment data to make specific predictions as to
when the WHO surveillance threshold is likely to be exceeded. Results: Our
results show - even if rates of acquired resistance are high, but the
drug-resistant strains that evolve are only half as transmissible as wild-type
strains - that transmission of drug-resistant strains will remain low (< 5% by
2009) and are unlikely to exceed the WHO's surveillance threshold. However,our
results show that transmission of drug-resistant strains in Botswana could
increase to ~15% by 2009 if resistant strains are as transmissible as wild-type
strains. Conclusion: The WHO's surveillance system is designed to detect
transmitted resistance that exceeds a threshold level of 5%. Whether this
system will detect drug-resistant strains in Botswana by 2009 will depend upon
the transmissibility of the strains that emerge. Our results imply that it
could be many years before the WHO detects transmitted resistance in other
sub-Saharan African countries with less ambitious treatment programs than
Botswana.


  We model the spontaneous assembly of a capsid (a virus's closed outer shell)
from many copies of identical units, using entirely irreversible steps and only
information local to the growing edge. Our model is formulated in terms of (i)
an elastic Hamiltonian with stretching and bending stiffness and a spontaneous
curvature, and (ii) a set of rate constants for addition of new units or bonds.
An ensemble of highly irregular capsids is generated, unlike the well-known
icosahedrally symmetric viruses, but (we argue) plausible as a way to model the
irregular capsids of retroviruses such as HIV. We found that (i) the
probability of successful capsid completion decays exponentially with capsid
size; (ii) capsid size depends strongly on spontaneous curvature and weakly on
the ratio of the bending and stretching elastic stiffnesses of the shell; (iii)
the degree of localization of Gaussian curvature (a measure of facetedness)
depends heavily on the ratio of elastic stiffnesses.


  The free energy of looping DNA by proteins and protein complexes determines
to what extent distal DNA sites can affect each other. We inferred its in vivo
value through a combined computational-experimental approach for different
lengths of the loop and found that, in addition to the intrinsic periodicity of
the DNA double helix, the free energy has an oscillatory component with about
half the helical period. Moreover, the oscillations have such an amplitude that
the effects of regulatory molecules become strongly dependent on their precise
DNA positioning and yet easily tunable by their cooperative interactions. These
unexpected results can confer to the physical properties of DNA a more
prominent role at shaping the properties of gene regulation than previously
thought.


  Many real world networks are characterized by adaptive changes in their
topology depending on the dynamic state of their nodes. Here we study epidemic
dynamics in an adaptive network, where susceptibles are able to avoid contact
with infected by rewiring their network connections. We demonstrate that
adaptive rewiring has profound consequences for the emerging network structure,
giving rise to assortative degree correlation and a separation into two loosely
connected sub-compartments. This leads to dynamics such as oscillations,
hysteresis and 1st order transitions. We describe the system in terms of a
simple model using a pair-approximation and present a full local bifurcation
analysis. Our results indicate that the interplay between dynamics and topology
can have important consequences for the spreading of infectious diseases and
related applications.


  We suggest a fast method to find possibly overlapping network communities of
a desired size and link density. Our method is a natural generalization of the
finite-$T$ superparamegnetic Potts clustering introduced by Blatt, Wiseman, and
Domany (Phys. Rev. Lett. v.76, 3251 (1996) and the recently suggested by
Reichard and Bornholdt (Phys. Rev. Lett. v.93, 21870 (2004)) annealing of Potts
model with global antiferromagnetic term. Similarly to both preceding works,
the proposed generalization is based on ordering of ferromagnetic Potts model;
the novelty of the proposed approach lies in the adjustable dependence of the
antiferromagnetic term on the population of each Potts state, which
interpolates between the two previously considered cases. This adjustability
allows to empirically tune the algorithm to detect the maximum number of
communities of the given size and link density. We illustrate the method by
detecting protein complexes in high-throughput protein binding networks.


  We study the general properties of stochastic two-species models for
predator-prey competition and coexistence with Lotka-Volterra type interactions
defined on a $d$-dimensional lattice. Introducing spatial degrees of freedom
and allowing for stochastic fluctuations generically invalidates the classical,
deterministic mean-field picture. Already within mean-field theory, however,
spatial constraints, modeling locally limited resources, lead to the emergence
of a continuous active-to-absorbing state phase transition. Field-theoretic
arguments, supported by Monte Carlo simulation results, indicate that this
transition, which represents an extinction threshold for the predator
population, is governed by the directed percolation universality class. In the
active state, where predators and prey coexist, the classical center
singularities with associated population cycles are replaced by either nodes or
foci. In the vicinity of the stable nodes, the system is characterized by
essentially stationary localized clusters of predators in a sea of prey. Near
the stable foci, however, the stochastic lattice Lotka-Volterra system displays
complex, correlated spatio-temporal patterns of competing activity fronts.
Correspondingly, the population densities in our numerical simulations turn out
to oscillate irregularly in time, with amplitudes that tend to zero in the
thermodynamic limit. Yet in finite systems these oscillatory fluctuations are
quite persistent, and their features are determined by the intrinsic
interaction rates rather than the initial conditions. We emphasize the
robustness of this scenario with respect to various model perturbations.


  We present a simple method for determining the exact noise power spectra in
linear chemical reaction networks. We apply the method to networks which are
representative of biochemical processes such as gene expression and signal
detection. Our results clarify how noise is transmitted by signal detection
motifs, and indicate how to coarse-grain networks by the elimination of fast
reactions.


  The Z_64-algebra of the genetic code and DNA sequences of length N was
recently stated. In order to beat the limits of this structure such as the
impossibility of non-coding region analysis in genomes and the impossibility of
the insertions and deletions analysis (indel mutations), we have develop a
cycle group structure over the of extended base triplets of DNA X_1X_2X_3, X_i
belong to {O, A, C, G, U}, where the letter O denote the base omission
(deletion) in the codon. The obtained group is isomorphic to the abelian
5-group Z_125 of integer module 125. Next, it is defined the abelian finite
group S over a set of DNA alignment sequences of length N. The group S could be
represented as the direct sum of homocyclic groups: 2-group and 5-group. In
particular, DNA subsequences without indel mutation could be considered
building block of genes represented by homocyclic 2-groups (described in the
previous Z_64-algebra). While those DNA subsequences affected by indel
mutations are described by means of homocyclic 5-groups. This representation
suggests identify genome block structures by way of a regular grammar capable
of recognize it. In addition, this novel structure allows us a general analysis
of the mutational pathways follow by genes and isofunctional genome regions by
means of the automorphism group on S.


  A general mechanism is proposed by which small intrinsic fluctuations in a
system far from equilibrium can result in nearly deterministic dynamical
behaviors which are markedly distinct from those realized in the meanfield
limit. The mechanism is demonstrated for the kinetic Monte-Carlo version of the
Schnakenberg reaction where we identified a scaling limit in which the global
deterministic bifurcation picture is fundamentally altered by fluctuations.
Numerical simulations of the model are found to be in quantitative agreement
with theoretical predictions.


  We have performed individual-based lattice simulations of SIR and SEIR
dynamics to investigate both the short and long-term dynamics of childhood
epidemics. In our model, infection takes place through a combination of local
and long-range contacts, in practice generating a dynamic small-world network.
Sustained oscillations emerge with a period much larger than the duration of
infection. We found that the network topology has a strong impact on the
amplitude of oscillations and in the level of persistence. Diseases do not
spread very effectively through local contacts. This can be seen by measuring
an {\em effective} transmission rate $\beta_{\mbox {\scriptsize eff}}$ as well
as the basic reproductive rate $R_0$. These quantities are lower in the
small-world network than in an homogeneously mixed population, whereas the
average age at infection is higher.


  How cooperation emerges in human societies is still a puzzle. Evolutionary
game theory has been the standard framework to address this issue. In most
models, every individual plays with all others, and then reproduce and die
according to what they earn. This amounts to assuming that selection takes
place at a slow pace with respect to the interaction time scale. We show that,
quite generally, if selection speeds up, the evolution outcome changes
dramatically. Thus, in games such as Harmony, where cooperation is the only
equilibrium and the only rational outcome, rapid selection leads to dominance
of defectors. Similar non trivial phenomena arise in other binary games and
even in more complicated settings such as the Ultimatum game. We conclude that
the rate of selection is a key element to understand and model the emergence of
cooperation, and one that has so far been overlooked.


  Samoilov, Plyasunov, and Arkin provide an example of a chemical reaction
whose full stochastic (Master Equation) model exhibits bistable behavior, but
for which the deterministic (mean field) version has a unique steady state at
least for special parameter values. In this short note, we provide a proof of
uniqueness valid for all possible parameter values.


  In a seminal paper Caspar and Klug established a theory that provides a
family of polyhedra as blueprints for the structural organisation of viral
capsids. In particular, they encode the locations of the proteins in the shells
that encapsulate, and hence provide protection for, the viral genome. Despite
of its huge success and numerous applications in virology experimental results
have provided evidence for the fact that the theory is too restrictive to
describe all known viruses. Especially, the family of Papovaviridae, which
contains cancer-causing viruses, falls out of the scope of this theory.
  In a recent paper we have shown that certain members of the family of
Papovaviridae can be described via tilings. In this paper, we develop a
comprehensive mathematical framework for the derivation of all surface
structures of viral particles in this family. We show that this formalism fixes
the structure and relative sizes of all particles collectively so that there
exists only one scaling factor that relates the sizes of all particles with
their biological counterparts.
  The series of polyhedra derived here complements the Caspar-Klug family of
polyhedra. It is the first mathematical result that provides a common
organisational principle for different types of viral particles in the family
of Papovaviridae and paves the way for an understanding of Papovaviridae
polymorphism. Moreover, it provides crucial input for the construction of
assembly models.


  We propose a formalism for describing two-strand systems of a DNA type by
means of soliton von Neumann equations, and illustrate how it works on a simple
example exactly solvably by a Darboux transformation. The main idea behind the
construction is the link between solutions of von Neumann equations and
entangled states of systems consisting of two subsystems evolving in time in
opposite directions. Such a time evolution has analogies in realistic DNA where
the polymerazes move on leading and lagging strands in opposite directions.


  Cooperative interactions among the binding of multiple signaling molecules is
a common mechanism for enhancing the sensitivity of biological signaling
systems. It is widely assumed that this increase in sensitivity of the mean
response implies the ability to detect smaller signals. We show that, quite
generally, there is a component of the noise in such systems that can be traced
to the random arrival of the signaling molecules at their receptor sites, and
this diffusive noise is not reduced by cooperativity. Cooperativity makes it
easier for real systems to reach this physical limit, but cannot reduce the
limit itself.


  Linear discrimination, from the point of view of numerical linear algebra,
can be treated as solving an ill-posed system of linear equations. In order to
generate a solution that is robust in the presence of noise, these problems
require regularization. Here, we examine the ill-posedness involved in the
linear discrimination of cancer gene expression data with respect to outcome
and tumor subclasses. We show that a filter factor representation, based upon
Singular Value Decomposition, yields insight into the numerical ill-posedness
of the hyperplane-based separation when applied to gene expression data. We
also show that this representation yields useful diagnostic tools for guiding
the selection of classifier parameters, thus leading to improved performance.