  The conjunction of insights from structural biology, solution biochemistry,
genetics and single molecule biophysics has provided a renewed impetus for the
construction of quantitative models of biological processes. One area that has
been a beneficiary of these experimental techniques is the study of viruses. In
this paper we describe how the insights obtained from such experiments can be
utilized to construct physical models of processes in the viral life cycle. We
focus on dsDNA bacteriophages and show that the bending elasticity of DNA and
its electrostatics in solution can be combined to determine the forces
experienced during packaging and ejection of the viral genome. Furthermore, we
quantitatively analyze the effect of fluid viscosity and capsid expansion on
the forces experienced during packaging. Finally, we present a model for DNA
ejection from bacteriophages based on the hypothesis that the energy stored in
the tightly packed genome within the capsid leads to its forceful ejection. The
predictions of our model can be tested through experiments in vitro where DNA
ejection is inhibited by the application of external osmotic pressure.


  We used a metacommunity of 49 discrete communities of aquatic invertebrates
to analyze the dynamical relationship between community and metacommunity
species distributions as a test of the neutral theory of biodiversity and
biogeography. At the community scale, observed variation in species richness
and relative abundance was greater than predicted by neutral models, and
revealed important differences among species in competitive ability and
tolerance for predation. At the metacommunity scale, species with metacommunity
proportions of less than 0.01% (38% of the observed metacommunity) were
consistently more abundant than predicted by models. Our results are at
variance with the neutral theory, and suggest that the use of an identical
survival probability for all species in neutral models misrepresents
substantial aspects of community assembly. Nevertheless, building and testing
neutral models can provide valuable insights into the processes that determine
species distributions.


  A newly developed statistical pair potential based on Distance-scaled Finite
Ideal-gas REference (DFIRE) state is applied to unbound protein-protein docking
structure selections. The performance of the DFIRE energy function is compared
to those of the well-established ZDOCK energy scores and RosettaDock energy
function using the comprehensive decoy sets generated by ZDOCK and RosettaDock.
Despite significant difference in the functional forms and complexities of the
three energy scores, the differences in overall performance for docking
structure selections are small between DFIRE and ZDOCK2.3 and between DFIRE and
RosettaDock. This result is remarkable considering that a single-term DFIRE
energy function was originally designed for monomer proteins while
multiple-term energy functions of ZDOCK and RosettaDock were specifically
optimized for docking. This provides hope that the accuracy of the existing
energy functions for docking can be improved.


  The Penna ageing model is based on mutation accumulation theory. We show that
it also allows for self-organization of antagonistic pleiotropy which helps at
young age at the expense of old age. This can be interpreted as emergence of
altruism.


  Amyloid fibers are aggregates of proteins. They are built out of a peptide
called $\beta$--amyloid (A$\beta$) containing between 41 and 43 residues,
produced by the action of an enzyme which cleaves a much larger protein known
as the Amyloid Precursor Protein (APP). X-ray diffraction experiments have
shown that these fibrils are rich in $\beta$--structures, whereas the shape of
the peptide displays an $\alpha$--helix structure within the APP in its
biologically active conformation. A realistic model of fibril formation is
developed based on the seventeen residues A$\beta$12--28 amyloid peptide, which
has been shown to form fibrils structurally similar to those of the whole
A$\beta$ peptide. With the help of physical arguments and in keeping with
experimental findings, the A$\beta$12--28 monomer is assumed to be in four
possible states (i.e., native helix conformation, $\beta$--hairpin, globular
low--energy state and unfolded state). Making use of these monomeric states,
oligomers (dimers, tertramers and octamers) were constructed. With the help of
short, detailed Molecular Dynamics (MD) calculations of the three monomers and
of a variety of oligomers, energies for these structures were obtained. Making
use of these results within the framework of a simple yet realistic model to
describe the entropic terms associated with the variety of amyloid
conformations, a phase diagram can be calculated of the whole many--body
system, leading to a thermodynamical picture in overall agreement with the
experimental findings. In particular, the existence of micellar metastable
states seem to be a key issue to determine the thermodynamical properties of
the system.


  The possibility of deriving the contact potentials between amino acids from
their frequencies of occurence in proteins is discussed in evolutionary terms.
This approach allows the use of traditional thermodynamics to describe such
frequencies and, consequently, to develop a strategy to include in the
calculations correlations due to the spatial proximity of the amino acids and
to their overall tendency of being conserved in proteins. Making use of a
lattice model to describe protein chains and defining a "true" potential, we
test these strategies by selecting a database of folding model sequences,
deriving the contact potentials from such sequences and comparing them with the
"true" potential. Taking into account correlations allows for a markedly better
prediction of the interaction potentials.


  While all the information required for the folding of a protein is contained
in its amino acid sequence, one has not yet learned how to extract this
information to predict the three--dimensional, biologically active, native
conformation of a protein whose sequence is known. Using insight obtained from
simple model simulations of the folding of proteins, in particular of the fact
that this phenomenon is essentially controlled by conserved (native) contacts
among (few) strongly interacting ("hot"), as a rule hydrophobic, amino acids,
which also stabilize local elementary structures (LES, hidden, incipient
secondary structures like $\alpha$--helices and $\beta$--sheets) formed early
in the folding process and leading to the postcritical folding nucleus (i.e.,
the minimum set of native contacts which bring the system pass beyond the
highest free--energy barrier found in the whole folding process) it is possible
to work out a succesful strategy for reading the native structure of designed
proteins from the knowledge of only their amino acid sequence and of the
contact energies among the amino acids. Because LES have undergone millions of
years of evolution to selectively dock to their complementary structures, small
peptides made out of the same amino acids as the LES are expected to
selectively attach to the newly expressed (unfolded) protein and inhibit its
folding, or to the native (fluctuating) native conformation and denaturate it.
These peptides, or their mimetic molecules, can thus be used as effective
non--conventional drugs to those already existing (and directed at neutralizing
the active site of enzymes), displaying the advantage of not suffering from the
uprise of resistance.


  Methods for alignment of protein sequences typically measure similarity by
using substitution matrix with scores for all possible exchanges of one amino
acid with another. Although widely used, the matrices derived from homologous
sequence segments, such as Dayhoff's PAM matrices and Henikoff's BLOSUM
matrices, are not specific for protein conformation identification. Using a
different approach, we got many amino acid segment blocks. For each of them,
the protein secondary structure is identical. Based on these blocks, we have
derived new amino acid substitution matrices. The application of these matrices
led to marked improvements in conformation segment search and homologues
detection in twilight zone.


  The formation of bubbles in nucleic acids (NAs) are fundamental in many
biological processes such as DNA replication, recombination, telomeres
formation, nucleotide excision repair, as well as RNA transcription and
splicing. These precesses are carried out by assembled complexes with enzymes
that separate selected regions of NAs. Within the frame of a nonlinear dynamics
approach we model the structure of the DNA duplex by a nonlinear network of
coupled oscillators. We show that in fact from certain local structural
distortions there originate oscillating localized patterns, that is radial and
torsional breathers, which are associated with localized H-bond deformations,
being reminiscent of the replication bubble. We further study the temperature
dependence of these oscillating bubbles. To this aim the underlying nonlinear
oscillator network of the DNA duplex is brought in contact with a heat bath
using the Nos$\rm{\acute{e}}$-Hoover-method. Special attention is paid to the
stability of the oscillating bubbles under the imposed thermal perturbations.
It is demonstrated that the radial and torsional breathers, sustain the impact
of thermal perturbations even at temperatures as high as room temperature.
Generally, for nonzero temperature the H-bond breathers move coherently along
the double chain whereas at T=0 standing radial and torsional breathers result.


  We study evolutionary processes induced by spatio-temporal dynamics in
prebiotic evolution. Using numerical simulations we demonstrate that
hypercycles emerge from complex interaction structures in multispecies systems.
In this work we also find that `hypercycle hybrid' protects the hypercycle from
its environment during the growth process. There is little selective advantage
for one hypercycle to maintain coexistence with others. This brings the
possibility of the outcompetition between hypercycles resulting in the negative
effect on information diversity. To enrich the information in hypercycles,
symbiosis with parasites is suggested. It is shown that symbiosis with
parasites can play an important role in the prebiotic immunology.


  We use Monte Carlo simulations to study attractive and excluded volume
interactions between nucleosome core particles in 30 nm-chromatin fibers. The
nucleosomes are treated as disk-like objects having an excluded volume and
short range attraction modelled by a variant of the Gay-Berne potential. The
nucleosomes are connected via bendable and twistable linker DNA in the crossed
linker fashion. We investigate the influence of the nucleosomal excluded volume
on the stiffness of the fiber. For parameter values that correspond to chicken
erythrocyte chromatin we find that the persistence length is governed to a
large extent by that excluded volume whereas the soft linker backbone
elasticity plays only a minor role. We further find that internucleosomal
attraction can induce the formation of hairpin configurations. Tension-induced
opening of such configurations into straight fibers manifests itself in a
quasi-plateau in the force-extension curve that resembles results from recent
micromanipulation experiments. Such hairpins may play a role in the formation
of higher order structures in chromosomes like chromonema fibers.


  Histone octamers show a heat-induced mobility along DNA. Recent theoretical
studies have established two mechanisms that are qualitatively and
quantitatively compatible with in vitro experiments on nucleosome sliding:
Octamer repositiong through one-basepair twist defects and through ten-basepair
bulge defects. A recent experiment demonstrated that the repositioning is
strongly suppressed in the presence of minor-groove binding DNA ligands. In the
present study we give a quantitative theory for nucleosome repositioning in the
presence of such ligands. We show that the experimentally observed octamer
mobilities are consistent with the picture of bound ligands blocking the
passage of twist defects through the nucleosome. This strongly supports the
model of twist defects inducing a corkscrew motion of the nucleosome as the
underlying mechanism of nucleosome sliding. We provide a theoretical estimate
of the nucleosomal mobility without adjustable parameters, as a function of
ligand concentration, binding affinity, binding site orientiation, temperature
and DNA anisotropy. Having this mobility at hand we speculate about the
interaction between a nucleosome and a transcribing RNA polymerase and suggest
a novel mechanism that might account for polymerase induced nucleosome
repositioning.


  A square lattice is introduced into the Penna model for biological aging in
order to study the evolution of diploid sexual populations under certain
conditions when one single locus in the individual's genome is considered as
identifier of species. The simulation results show, after several generations,
the flourishing and coexistence of two separate species in the same
environment, i.e., one original species splits up into two on the same
territory (sympatric speciation). As well, the mortalities obtained are in a
good agreement with the Gompertz law of exponential increase of mortality with
age.


  The ability of oligonucleotide microarrays to measure gene expression has
been hindered by an imperfect understanding of the relationship between input
RNA concentrations and output signals. We argue that this relationship can be
understood based on the underlying statistical mechanics of these devices. We
present a model that includes the relevant interactions between the molecules.
Our model for the first time accounts for partially zippered probe-target
hybrids in a physically realistic manner, and also includes target-target
binding in solution. Large segments of the target molecules are not bound to
the probes, often in an asymmetric pattern, emphasizing the importance of
modeling zippering properly. The resultant fit between the model and training
data using optimized parameters is excellent, and it also does well at
predicting test data.


  Biological networks have been recently found to exhibit many topological
properties of the so-called complex networks. It has been reported that they
are, in general, both highly skewed and directed. In this paper, we report on
the dynamics of a Michaelis-Menten like model when the topological features of
the underlying network resemble those of real biological networks.
Specifically, instead of using a random graph topology, we deal with a complex
heterogeneous network characterized by a power-law degree distribution coupled
to a continuous dynamics for each network's component. The dynamics of the
model is very rich and stationary, periodic and chaotic states are observed
upon variation of the model's parameters. We characterize these states
numerically and report on several quantities such as the system's phase diagram
and size distributions of clusters of stationary, periodic and chaotic nodes.
The results are discussed in view of recent debate about the ubiquity of
complex networks in nature and on the basis of several biological processes
that can be well described by the dynamics studied.


  A four-state mutation-selection model for the evolution of populations of
DNA-sequences is investigated with particular interest in the phenomenon of
error thresholds. The mutation model considered is the Kimura 3ST mutation
scheme, fitness functions, which determine the selection process, come from the
permutation-invariant class. Error thresholds can be found for various fitness
functions, the phase diagrams are more interesting than for equivalent
two-state models. Results for (small) finite sequence lengths are compared with
those for infinite sequence length, obtained via a maximum principle that is
equivalent to the principle of minimal free energy in physics.


  High-throughput protein interaction detection methods are strongly affected
by false positive and false negative results. Focused experiments are needed to
complement the large-scale methods by validating previously detected
interactions but it is often difficult to decide which proteins to probe as
interaction partners. Developing reliable computational methods assisting this
decision process is a pressing need in bioinformatics. We show that we can use
the conserved properties of the protein network to identify and validate
interaction candidates. We apply a number of machine learning algorithms to the
protein connectivity information and achieve a surprisingly good overall
performance in predicting interacting proteins. Using a 'leave-one-out'
approach we find average success rates between 20-50% for predicting the
correct interaction partner of a protein. We demonstrate that the success of
these methods is based on the presence of conserved interaction motifs within
the network. A reference implementation and a table with candidate interacting
partners for each yeast protein are available at http://www.protsuggest.org


  It is becoming increasingly appreciated that the signal transduction systems
used by eukaryotic cells to achieve a variety of essential responses represent
highly complex networks rather than simple linear pathways. While significant
effort is being made to experimentally measure the rate constants for
individual steps in these signaling networks, many of the parameters required
to describe the behavior of these systems remain unknown, or at best,
estimates. With these goals and caveats in mind, we use methods of statistical
mechanics to extract useful predictions for complex cellular signaling
networks. To establish the usefulness of our approach, we have applied our
methods towards modeling the nerve growth factor (NGF)-induced differentiation
of neuronal cells. Using our approach, we are able to extract predictions that
are highly specific and accurate, thereby enabling us to predict the influence
of specific signaling modules in determining the integrated cellular response
to the two growth factors. We show that extracting biologically relevant
predictions from complex signaling models appears to be possible even in the
absence of measurements of all the individual rate constants. Our methods also
raise some interesting insights into the design and possible evolution of
cellular systems, highlighting an inherent property of these systems wherein
particular ''soft'' combinations of parameters can be varied over wide ranges
without impacting the final output and demonstrating that a few ''stiff''
parameter combinations center around the paramount regulatory steps of the
network. We refer to this property -- which is distinct from robustness -- as
''sloppiness.''


  We determine conditions under which a random biochemical system is likely to
contain a subsystem that is both autocatalytic and able to survive on some
ambient `food' source. Such systems have previously been investigated for their
relevance to origin-of-life models. In this paper we extend earlier work, by
finding precisely the order of catalysation required for the emergence of such
self-sustaining autocatalytic networks. This answers questions raised in
earlier papers, yet also allows for a more general class of models. We also
show that a recently-described polynomial time algorithm for determining
whether a catalytic reaction system contains an autocatalytic, self-sustaining
subsystem is unlikely to adapt to allow inhibitory catalysation - in this case
we show that the associated decision problem is NP-complete.


  A new approach to the modular, complex systems analysis of nonlinear dynamics
in cell cycling network transformations involved in carcinogenesis is proposed.
Carcinogenesis is a complex process that involves dynamically inter-connected
biomolecules in the intercellular, membrane, cytosolic, nuclear and nucleolar
compartments that form numerous inter-related pathways. One such family of
pathways contains the cell cyclins. Cyclins are proteins that link several
critical pro-apoptotic and other cell cycling/division components, including
the tumor suppressor gene TP53 and its product, the Thomsen-Friedenreich
antigen (T antigen), Rb, mdm2, c-Myc, p21, p27, Bax, Bad and Bcl-2, which all
play major roles in neoplastic transformation of many tissues. This novel
theoretical analysis based on recently published studies of cyclin signaling,
with special emphasis placed on the roles of cyclins D1 and E, suggests novel
clinical trials and rational therapies of cancer through reestablishment of
cell cycling inhibition in metastatic cancer cells.


  In this paper we review some recent results that shed light on a fundamental
question in molecular systematics: how much phylogenetic `signal' can we expect
from characters that have evolved under some Markov process? There are many
sides to this question and we begin by describing some explicit bounds on the
probability of correctly reconstructing an ancestral state from the states
observed at the tips. We show how this bound sets upper limits on the
probability of tree reconstruction from aligned sequences, and we provide some
new extensions that allow site-to-site rate variation or a covarion mechanism.
We then explore the relationship between the number of sites required for
accurate tree reconstruction and other model parameters - such as the number of
species, and substitution probabilities, and we describe a phase transition
that occurs when substitution probabilities exceed a critical value. In the
remainder of this paper we turn to models of character evolution where the
state space is assumed to be either infinite or very large. These models have
some relevance to certain types of genomic data (such as gene order) and here
we again investigate how many characters are required for accurate tree
reconstruction.


  We define and completely solve a content-based directed network whose nodes
consist of random words and an adjacency rule involving perfect or approximate
matches, for an alphabet with an arbitrary number of letters. The analytic
expression for the out-degree distribution shows a crossover from a leading
power law behavior to a log-periodic regime bounded by a different power law
decay. The leading exponents in the two regions have a weak dependence on the
mean word length, and an even weaker dependence on the alphabet size. The
in-degree distribution, on the other hand, is much narrower and does not show
scaling behavior. The results might be of interest for understanding the
emergence of genomic interaction networks, which rely, to a large extent, on
mechanisms based on sequence matching, and exhibit similar global features to
those found here.


  The advent of new experimental genomic technologies and the massive increase
of DNA sequence information is helping researchers better understand how our
genes work. Recently, experiments on mRNA abundance (gene expression) have
revealed that gene expression shows a stationary organization described by a
power-law distribution (scale-free organization) (i.e., gene expression $k$
decays as $k^{-\gamma}$), which is highly conserved in all the major five
kingdoms of life, from Bacteria to Human. An underlying gene expression
dynamics "rich-travel-more" was suggested to recover that evolutional
conservation of transcriptional organization. Here we propose a constructive
approach to gene expression dynamics with larger scope. Our gene expression
construction restores the stationary state, predicts the power-law exponent for
different organisms with natural explanation for small correction at high and
low expression levels, describes the intermediate state dynamics (time finite)
and elucidates the gene expression stability. This approach requires only one
assumption: Markov property.


  We study the effects of the sequence on the propagation of nonlinear
excitations in simple models of DNA in which we incorporate actual DNA
sequences obtained from human genome data. We show that kink propagation
requires forces over a certain threshold, a phenomenon already found for
aperiodic sequences [F. Dom\'\i nguez-Adame {\em et al.}, Phys. Rev. E {\bf
52}, 2183 (1995)]. For forces below threshold, the final stop positions are
highly dependent on the specific sequence. The results of our model are
consistent with the stick-slip dynamics of the unzipping process observed in
experiments. We also show that the effective potential, a collective coordinate
formalism introduced by Salerno and Kivshar [Phys. Lett. A {\bf 193}, 263
(1994)] is a useful tool to identify key regions in DNA that control the
dynamical behavior of large segments. Additionally, our results lead to further
insights in the phenomenology observed in aperiodic systems.


  A preliminary analysis of the temporal evolution of a population of
\emph{Peromyscus maniculatus} infected with Hantavirus Sin Nombre is made.
Ecological and epidemiological parameters are derived from the data, and they
are used as inputs for the analytical model presented in [Abramson and Kenkre,
Phys. Rev. E \textbf{66}, 011912 (2002)]. A prediction of the critical carrying
capacity and its associated critical mouse density is made, and the time series
is analyzed under the light of these. It is found that the sporadic
disappearances and reappearances of the infected phase correspond to the
bifurcation predicted by the model.


  This paper extends the semiconservative quasispecies equations to account for
arbitrary post-replication lesion repair efficiency. Such an extension could be
an important tool for understanding processes such as cancer development and
stem cell growth. Starting from the quasispecies dynamics over the space of
genomes, we derive an equivalent dynamics over the space of ordered sequence
pairs. From this set of equations, we are able to derive the infinite sequence
length form of the dynamics for a class of ``master-genome''-based fitness
landscapes. We use these equations to solve for a ``generalized''
single-fitness-peak landscape, where the master genome can sustain a maximum
number of lesions and remain viable. The central pattern that emerges from our
studies is that imperfect lesion repair often leads to increased mutational
robustness over semiconservative replication with completely efficient lesion
repair. The reason for this is that imperfect lesion repair breaks some of the
correlation between the parent and daughter strands, thereby preventing
replication errors from destroying the information in the original genome. The
result is a delayed error catastrophe over that expected from the original
semiconservative quasispecies model. In particular, we show that when only of
the strands is necessary for conferring viability, then, when lesion repair is
turned off, a semiconservatively replicating system becomes an effectively
conservatively replicating one.


  Using a simple model for the trail formation of ants, the relation between
i)the schedule of feeding which represents the unsteady natural environment,
ii)emerging patterns of trails connecting a nest with food resources, and
iii)the foraging efficiency is studied. Simulations and a simple analysis show
that the emergent trail pattern flexibly varies depending on the feeding
schedule by which ants can make an efficient foraging according to the
underlying unsteady environment.


  Reaction-diffusion equations based on a polymerization model are solved to
simulate the spreading of hypothetic left and right handed life forms on the
Earth's surface. The equations exhibit front-like behavior as is familiar from
the theory of the spreading of epidemics. It is shown that the relevant time
scale for achieving global homochirality is not, however, the time scale of
front propagation, but the much longer global diffusion time. The process can
be sped up by turbulence and large scale flows. It is speculated that, if the
deep layers of the early ocean were sufficiently quiescent, there may have been
the possibility of competing early life forms with opposite handedness.


  The internal cell wall structure of the bacterium Lactobacillus helveticus
has been observed in situ in aqueous solution using an atomic force microscope
(AFM). The AFM tip was used not only for imaging but presumably to remove
mechanically large patches of the outer cell wall after appropriate chemical
treatment, which typically leaves the bacteria alive. The surface exposed after
such a surgery revealed ca. 26 nm thick twisted strands within the cell wall.
The structure and location of the observed strands are consistent with the
glycan backbone of peptidoglycan fibers that give strength to the cell wall.
The found structural organization of these fibers has not been observed
previously.


  By observing reconstituted chromatin by fluorescence microscopy (FM) and
atomic force microscopy (AFM), we found that the density of nucleosomes
exhibits a bimodal profile, i.e., there is a large transition between the dense
and dispersed states in reconstituted chromatin. Based on an analysis of the
spatial distribution of nucleosome cores, we deduced an effective thermodynamic
potential as a function of the nucleosome-nucleosome distance. This enabled us
to interpret the folding transition of chromatin in terms of a first-order
phase transition. This mechanism for the condensation of chromatin is discussed
in terms of its biological significance.


  Concomitant with the evolution of biological diversity must have been the
evolution of mechanisms that facilitate evolution, due to the essentially
infinite complexity of protein sequence space. We describe how evolvability can
be an object of Darwinian selection, emphasizing the collective nature of the
process. We quantify our theory with computer simulations of protein evolution.
These simulations demonstrate that rapid or dramatic environmental change leads
to selection for greater evolvability. The selective pressure for large scale
genetic moves, such as DNA exchange, becomes increasingly strong as the
environmental conditions become more uncertain. Our results demonstrate that
evolvability is a selectable trait and allow for the explanation of a large
body of experimental results.


  Receptor-ligand binding is a critical first step in signal transduction and
the duration of the interaction can impact signal generation. In mammalian
cells, clustering of receptors may be facilitated by heterogeneous zones of
lipids, known as lipid rafts. In vitro experiments show that disruption of
rafts significantly alters the dissociation of fibroblast growth factor-2
(FGF-2) from heparan sulfate proteoglycans, co-receptors for FGF-2. In this
paper, we develop a continuum stochastic formalism in order to (i) study how
rebinding affects the dissociation of ligands from a planar substrate, and (ii)
address the question of how receptor clustering influences ligand rebinding. We
find that clusters reduce the effective dissociation rate dramatically when the
clusters are dense and the overall surface density of receptors is low. The
effect is much less pronounced in the case of high receptor density and shows
non-monotonic behavior with time. These predictions are verified via lattice
Monte Carlo simulations. Comparison with experimental results suggests that the
theory does not capture the complete biological system. We speculate that
additional co-operative mechanisms might be present in order to increase ligand
retention, and present one possible ``internal diffusion'' model.


  We study the basepairing dynamics of DNA with repetitive sequences where
local strand slippage can create, annihilate, and move bulge loops. Using an
explicit theoretical model, we find a rich dynamical behavior as a function of
an applied shear force f: reptation-like dynamics at f=f_c with a rupture time
scaling as N^3 with its length N, drift-diffusion dynamics for f_c < f < f^*,
and a dynamical transition to an unraveling mode of strand separation at f =
f^*. We predict a viscoelastic behavior for periodic DNA with time and force
scales that can be programmed into its sequence.


  Here we study the effects of many-body interactions on rate and mechanism in
protein folding, using the results of molecular dynamics simulations on
numerous coarse-grained C-alpha-model single-domain proteins. After adding
three-body interactions explicitly as a perturbation to a Go-like Hamiltonian
with native pair-wise interactions only, we have found 1) a significantly
increased correlation with experimental phi-values and folding rates, 2) a
stronger correlation of folding rate with contact order, matching the
experimental range in rates when the fraction of three-body energy in the
native state is ~ 20%, and 3) a considerably larger amount of 3-body energy
present in Chymotripsin inhibitor than other proteins studied.


  We study the spread of Hantavirus over a host population of deer mice using a
population dynamics model. We show that taking into account the internal
fluctuations in the mouse population due to its discrete character strongly
alters the behaviour of the system. In addition to the familiar transition
present in the deterministic model, the inclusion of internal fluctuations
leads to the emergence of an additional deterministically hidden transition. We
determine parameter values that lead to maximal propagation of the disease, and
discuss some implications for disease prevention policies.


  We analyze a mathematical model of the epizootic of Hantavirus in mice
population, including the effect of species which compete with the host. We
show that the existence of the second species has an important consequence for
the prevalence of the infectious agent in the host. When the two mice species
survive in the ecosystem, the competitive pressure of the second species may
lead to reduction or complete elimination of the prevalence of infection. The
transition between the disappearance of the infection and its presence occurs
at a critical value of the competitor's population, resembling a second order
phase transition in a statistical system. The results provide a rigorous
framework for the study of the impact of biodiversity in the propagation of
infectious diseases, and further lends itself to future experimental
verification.


  By observing trends in the folding kinetics of experimental 2-state proteins
at their transition midpoints, and by observing trends in the barrier heights
of numerous simulations of coarse grained, C-alpha model, Go proteins, we show
that folding rates correlate with the degree of heterogeneity in the formation
of native contacts. Statistically significant correlations are observed between
folding rates and measures of heterogeneity inherent in the native topology, as
well as between rates and the variance in the distribution of either
experimentally measured or simulated phi-values.


  Proteins are minimally frustrated polymers. However, for realistic protein
models non-native interactions must be taken into account. In this paper we
analyze the effect of non-native interactions on the folding rate and on the
folding free energy barrier. We present an analytic theory to account for the
modification on the free energy landscape upon introduction of non-native
contacts, added as a perturbation to the strong native interactions driving
folding. Our theory predicts a rate-enhancement regime at fixed temperature,
under the introduction of weak, non-native interactions. We have thoroughly
tested this theoretical prediction with simulations of a coarse-grained protein
model, by employing an off-lattice $C_\alpha$ model of the src-SH3 domain. The
strong agreement between results from simulations and theory confirm the non
trivial result that a relatively small amount of non-native interaction energy
can actually assist the folding to the native structure.


  The prediction of the biologically active native conformation of a protein is
one of the fundamental challenges of structural biology. This problem remains
yet unsolved mainly due to three factors: the partial knowledge of the
effective free energy function that governs the folding process, the enormous
size of the conformational space of a protein and, finally, the relatively
small differences of energy between conformations, in particular, between the
native one and the ones that make up the unfolded state.
  Herein, we recall the importance of taking into account, in a detailed
manner, the many interactions involved in the protein folding problem (such as
steric volume exclusion, Ramachandran forces, hydrogen bonds, weakly polar
interactions, coulombic energy or hydrophobic attraction) and we propose a
strategy to effectively construct a free energy function that, including the
effects of the solvent, could be numerically tractable. It must be pointed out
that, since the internal free energy function that is mainly described does not
include the constraints of the native conformation, it could only help to reach
the 'molten globule' state. We also discuss about the limits and the lacks from
which suffer the simple models that we, physicists, love so much.


  Can unicellular organisms survive a drastic temperature change, and adapt to
it after many generations? In simulations of the Penna model of biological
ageing, both extinction and adaptation were found for asexual and sexual
reproduction as well as for parasex. These model investigations are the basis
for the design of evolution experiments with heterotrophic flagellates.


  It has recently been discovered that many biological systems, when
represented as graphs, exhibit a scale-free topology. One such system is the
set of structural relationships among protein domains. The scale-free nature of
this and other systems has previously been explained using network growth
models that, while motivated by biological processes, do not explicitly
consider the underlying physics or biology. In the present work we explore a
sequence-based model for the evolution protein structures and demonstrate that
this model is able to recapitulate the scale-free nature observed in graphs of
real protein structures. We find that this model also reproduces other
statistical feature of the protein domain graph. This represents, to our
knowledge, the first such microscopic, physics-based evolutionary model for a
scale-free network of biological importance and as such has strong implications
for our understanding of the evolution of protein structures and of other
biological networks.


  Stacking interactions in single stranded nucleic acids give rise to
configurations of an annealed rod-coil multiblock copolymer. Theoretical
analysis identifies the resulting signatures for long homopolynucleotides: A
non monotonous dependence of size on temperature, corresponding effects on
cyclization and a plateau in the extension force law. Explicit numerical
results for poly(dA) and poly(rU) are presented.


  An effective potential function is critical for protein structure prediction
and folding simulation. For simplified models of proteins where coordinates of
only $C_\alpha$ atoms need to be specified, an accurate potential function is
important. Such a simplified model is essential for efficient search of
conformational space. In this work, we present a formulation of potential
function for simplified representations of protein structures. It is based on
the combination of descriptors derived from residue-residue contact and
sequence-dependent local geometry. The optimal weight coefficients for contact
and local geometry is obtained through optimization by maximizing margins among
native and decoy structures. The latter are generated by chain growth and by
gapless threading. The performance of the potential function in blind test of
discriminating native protein structures from decoys is evaluated using several
benchmark decoy sets. This potential function have comparable or better
performance than several residue-based potential functions that require in
addition coordinates of side chain centers or coordinates of all side chain
atoms.


  A four-node network consisting of a negative loop controlling a positive one
is studied. It models some of the features of the p53 gene network. Using
piecewise linear dynamics with thresholds, the allowed dynamical classes are
fully characterized and coded. The biologically relevant situations are
identified and conclusions drawn concerning the effectiveness of the p53
network as a tumour inhibitor mechanism.


  The observed correlations between pairs of homologous protein sequences are
typically explained in terms of a Markovian dynamic of amino acid substitution.
This model assumes that every location on the protein sequence has the same
background distribution of amino acids, an assumption that is incompatible with
the observed heterogeneity of protein amino acid profiles and with the success
of profile multiple sequence alignment. We propose an alternative model of
amino acid replacement during protein evolution based upon the assumption that
the variation of the amino acid background distribution from one residue to the
next is sufficient to explain the observed sequence correlations of homologs.
The resulting dynamical model of independent replacements drawn from
heterogeneous backgrounds is simple and consistent, and provides a unified
homology match score for sequence-sequence, sequence-profile and
profile-profile alignment.


  Recently we have extended our the "unified" model of evolutionary ecology to
incorporate the {\it spatial inhomogeneities} of the eco-system and the {\it
migration} of individual organisms from one patch to another within the same
eco-system. In this paper an extension of our recent model is investigated so
as to describe the {\it migration} and {\it speciation} in a more realistic
way.


  Increasingly used high throughput experimental techniques, like DNA or
protein microarrays give as a result groups of interesting, e.g. differentially
regulated genes which require further biological interpretation. With the
systematic functional annotation provided by the Gene Ontology the information
required to automate the interpretation task is now accessible. However, the
determination of statistical significant e.g. molecular functions within these
groups is still an open question. In answering this question, multiple testing
issues must be taken into account to avoid misleading results. Here we present
a statistical framework that tests whether functions, processes or locations
described in the Gene Ontology are significantly enriched within a group of
interesting genes when compared to a reference group. First we define an exact
analytical expression for the expected number of false positives that allows us
to calculate adjusted p-values to control the false discovery rate. Next, we
demonstrate and discuss the capabilities of our approach using publicly
available microarray data on cell-cycle regulated genes. Further, we analyze
the robustness of our framework with respect to the exact gene group
composition and compare the performance with earlier approaches. The software
package GOSSIP implements our method and is made freely available at
http://gossip.gene-groups.net/


  Changing base composition during the evolution of biological sequences can
mislead some of the phylogenetic inference techniques in current use. However,
detecting whether such a process has occurred may be difficult, since
convergent evolution may lead to similar base frequencies emerging from
different lineages.
  To study this situation, algebraic models of biological sequence evolution
are introduced in which the base composition is fixed throughout evolution.
Basic properties of the associated algebraic varieties are investigated,
including the construction of some phylogenetic invariants.


  The Poland-Scheraga (PS) model for the helix-coil transition of DNA considers
the statistical mechanics of the thermally induced binding of two complementary
strands of DNA. In this paper, we show how to modify the PS model when a torque
is applied to the extremities of DNA: We propose a simple model for the energy
of twisted DNA and compute the entropy of a loop, subject to angular
constraints (supercoiling). The denaturation curves are shifted towards lower
or higher temperatures depending on the sign of the torque, and the UV
absorption peaks are softened. The properties of supercoiled DNA can be deduced
through the use of a numerical Legendre transform. In the homogeneous case, we
find that for weak supercoiling, the phenomenological quadratic law relating
the torsional energy to the number of unpaired bases is recovered.


  Motivation:Microarray experiments result in large scale data sets that
require extensive mining and refining to extract useful information. We
demonstrate the usefulness of (nonmetric) multidimensional scaling (MDS) method
in analyzing a large number of genes. Applying MDS to the microarray data is
certainly not new, but the existing works are all on small numbers
  (< 100) of points to be analyzed. We have been developing an efficient novel
algorithm for nonmetric multidimensional scaling (nMDS) analysis for very large
data sets as a maximally unsupervised data mining device. We wish to
demonstrate its usefulness in the context of bioinformatics (unraveling
relational patterns among genes from time series data in this paper).
  Results: The Pearson correlation coefficient with its sign flipped is used to
measure the dissimilarity of the gene activities in transcriptional response of
cell-cycle-synchronized human fibroblasts to serum [Iyer {\it et al}., Science
{\bf 283}, 83 (1999)]. These dissimilarity data have been analyzed with our
nMDS algorithm to produce an almost circular relational pattern of the genes.
The obtained pattern expresses a temporal order in the data in this example;
the temporal expression pattern of the genes rotates along this circular
arrangement and is related to the cell cycle. For the data we analyze in this
paper we observe the following. If an appropriate preparation procedure is
applied to the original data set, linear methods such as the principal
component analysis (PCA) could achieve reasonable results, but without data
preprocessing linear methods such as PCA cannot achieve a useful picture.
Furthermore, even with an appropriate data preprocessing, the outcomes of
linear procedures are not as clearcut as those by nMDS without preprocessing.


  Circular permutation connects the N and C termini of a protein and
concurrently cleaves elsewhere in the chain, providing an important mechanism
for generating novel protein fold and functions. However, their in genomes is
unknown because current detection methods can miss many occurances, mistaking
random repeats as circular permutation. Here we develop a method for detecting
circularly permuted proteins from structural comparison. Sequence order
independent alignment of protein structures can be regarded as a special case
of the maximum-weight independent set problem, which is known to be
computationally hard. We develop an efficient approximation algorithm by
repeatedly solving relaxations of an appropriate intermediate integer
programming formulation, we show that the approximation ratio is much better
then the theoretical worst case ratio of $r = 1/4$. Circularly permuted
proteins reported in literature can be identified rapidly with our method,
while they escape the detection by publicly available servers for structural
alignment.


  Motivation. Protein design aims to identify sequences compatible with a given
protein fold but incompatible to any alternative folds. To select the correct
sequences and to guide the search process, a design scoring function is
critically important. Such a scoring function should be able to characterize
the global fitness landscape of many proteins simultaneously.
  Results. To find optimal design scoring functions, we introduce two geometric
views and propose a formulation using mixture of nonlinear Gaussian kernel
functions. We aim to solve a simplified protein sequence design problem. Our
goal is to distinguish each native sequence for a major portion of
representative protein structures from a large number of alternative decoy
sequences, each a fragment from proteins of different fold. Our scoring
function discriminate perfectly a set of 440 native proteins from 14 million
sequence decoys. We show that no linear scoring function can succeed in this
task. In a blind test of unrelated proteins, our scoring function misclassfies
only 13 native proteins out of 194. This compares favorably with about 3-4
times more misclassifications when optimal linear functions reported in
literature are used. We also discuss how to develop protein folding scoring
function.


  We present a self-consistent field approximation to the problem of the
genetic switch composed of two mutually repressing/activating genes. The
protein and DNA state dynamics are treated stochastically and on equal footing.
In this approach the mean influence of the proteomic cloud created by one gene
on the action of another is self-consistently computed. Within this
approximation a broad range of stochastic genetic switches may be solved
exactly in terms of finding the probability distribution and its moments. A
much larger class of problems, such as genetic networks and cascades also
remain exactly solvable with this approximation. We discuss in depth certain
specific types of basic switches, which are used by biological systems and
compare their behavior to the expectation for a deterministic switch.


  The statistics of outbreaks in a model for the propagation of meningococcal
diseases is analyzed, taking into account the possibility that the population
is fragmented into weakly connected patches. It is shown that, depending on the
size of of the sample studied, the ratio between the variance and the mean of
infected cases can vary from one (Poisson statistics) to the inverse of the
infection rate.


  The effect of spatial correlations on the spread of infectious diseases was
investigated using a stochastic SIR (Susceptible-Infective-Recovered) model on
complex networks. It was found that in addition to the reduction of the
effective transmission rate, through the screening of infectives, spatial
correlations may have another major effect through the enhancement of
stochastic fluctuations. As a consequence large populations will have to become
even larger to 'average out' significant differences from the solution of
deterministic models. Furthermore, time series of the (unforced) model provide
patterns of recurrent epidemics with slightly irregular periods and realistic
amplitudes, suggesting that stochastic models together with complex networks of
contacts may be sufficient to describe the long term dynamics of some diseases.
The spatial effects were analysed quantitatively by modelling measles and
pertussis, using a SEIR (Susceptible-Exposed-Infective-Recovered) model. Both
the period and the spatial coherence of the epidemic peaks of pertussis are
well described by the unforced model for realistic values of the parameters.


  We study a set of six-species ecological models where each species has two
predators and two preys. On a square lattice the time evolution is governed by
iterated invasions between the neighboring predator-prey pairs chosen at random
and by a site exchange with a probability Xs between the neutral pairs. These
models involve the possibility of spontaneous formation of different defensive
alliances whose members protect each other from the external invaders. The
Monte Carlo simulations show a surprisingly rich variety of the stable spatial
distributions of species and subsequent phase transitions when tuning the
control parameter Xs. These very simple models are able to demonstrate that the
competition between these associations influences their composition. Sometimes
the dominant association is developed via a domain growth. In other cases
larger and larger invasion processes preceed the prevalence of one of the
stable asociations. Under some conditions the survival of all the species can
be maintained by the cyclic dominance occuring between these associations.


  Understanding the design of the universe of protein structures may provide
insights into protein evolution. We study the architecture of the protein
domain universe, which has been found to poses peculiar scale-free properties
(Dokholyan et al., Proc. Natl. Acad. Sci. USA 99: 14132-14136 (2002)). We
examine the origin of these scale-free properties of the graph of protein
domain structures (PDUG) and determine that that the PDUG is not modular, i.e.
it does not consist of modules with uniform properties. Instead, we find the
PDUG to be self-similar at all scales. We further characterize the PDUG
architecture by studying the properties of the hub nodes that are responsible
for the scale-free connectivity of the PDUG. We introduce a measure of the
betweenness centrality of protein domains in the PDUG and find a power-law
distribution of the betweenness centrality values. The scale-free distribution
of hubs in the protein universe suggests that a set of specific statistical
mechanics models, such as the self-organized criticality model, can potentially
identify the principal driving forces of molecular evolution. We also find a
gatekeeper protein domain, removal of which partitions the largest cluster into
two large sub-clusters. We suggest that the loss of such gatekeeper protein
domains in the course of evolution is responsible for the creation of new fold
families.


  In this work we employ various methods of analysis (unfolding simulations and
comparative analysis of structures and sequences of proteomes of thermophilic
organisms) to show that organisms can follow two major strategies of
thermophilic adaptation: (i) General, non-specific, structure-based, when
proteomes of certain thermophilic organisms show significant structural bias
toward proteins of higher compactness. In this case thermostability is achieved
by greater overall number of stabilizing contacts, none of which may be
especially strong, and (ii) Specific, sequence-based, whereby sequence
variations aimed at strengthening specific types of interactions (e.g.
electrostatics) are applied without significantly changing structures of
proteins. The choice of a certain strategy is a direct consequence of
evolutionary history and environmental conditions of particular (hyper)
thermophilic species: ancient hyperthermophilic organisms that directly evolved
in hot environment, pursued mostly structure-based strategy, while later
evolved organisms whose thermophilic adaptation was a consequence of their
recolonization of hot environment, pursued specific, sequence-based strategy of
thermophilic adaptation.


  Protein structures can be studied as complex networks of interacting amino
acids. We study proteins of different structural classes from the network
perspective. Our results indicate that proteins, regardless of their structural
class, show small-world network property. Various network parameters offer
insight into the structural organisation of proteins and provide indications of
modularity in protein networks.


  Naturally occurring networks exhibit quantitative features revealing
underlying growth mechanisms. Numerous network mechanisms have recently been
proposed to reproduce specific properties such as degree distributions or
clustering coefficients. We present a method for inferring the mechanism most
accurately capturing a given network topology, exploiting discriminative tools
from machine learning. The Drosophila melanogaster protein network is
confidently and robustly (to noise and training data subsampling) classified as
a duplication-mutation-complementation network over preferential attachment,
small-world, and other duplication-mutation mechanisms. Systematic
classification, rather than statistical study of specific properties, provides
a discriminative approach to understand the design of complex networks.


  The rates-across-sites assumption in phylogenetic inference posits that the
rate matrix governing the Markovian evolution of a character on an edge of the
putative phylogenetic tree is the product of a character-specific scale factor
and a rate matrix that is particular to that edge. Thus, evolution follows
basically the same process for all characters, except that it occurs faster for
some characters than others. To allow estimation of tree topologies and edge
lengths for such models, it is commonly assumed that the scale factors are not
arbitrary unknown constants, but rather unobserved, independent, identically
distributed draws from a member of some parametric family of distributions. A
popular choice is the gamma family. We consider an example of a clock-like tree
with three taxa, one unknown edge length, and a parametric family of scale
factor distributions that contain the gamma family. This model has the property
that, for a generic choice of unknown edge length and scale factor
distribution, there is another edge length and scale factor distribution which
generates data with exactly the same distribution, so that even with infinitely
many data it will be typically impossible to make correct inferences about the
unknown edge length.


  Being HIV-1-PR an essential enzyme in the viral life cycle, its inhibition
can control AIDS. Because the folding of single domain proteins, like HIV-1-PR
is controlled by local elementary structures (LES, folding units stabilized by
strongly interacting, highly conserved amino acids) which have evolved over
myriads of generations to recognize and strongly attract each other so as to
make the protein fold fast, we suggest a novel type of HIV-1-PR inhibitors
which interfere with the folding of the protein: short peptides displaying the
same amino acid sequence of that of LES. Theoretical and experimental evidence
for the specificity and efficiency of such inhibitors are presented.


  Herein it is shown that in order to study the statistical properties of DNA
sequences in bacterial chromosomes it suffices to consider only one half of the
chromosome because they are similar to its corresponding complementary sequence
in the other half. This is due to the inverse bilateral symmetry of bacterial
chromosomes. Contrary to the classical result that DNA coding regions of
bacterial genomes are purely uncorrelated random sequences, here it is shown,
via a renormalization group approach, that DNA random fluctuations of single
bases are modulated by log-periodic variations. Distance series of triplets
display long-range correlations in each half of the intact chromosome and in
intronless protein-coding sequences, or both long-range correlations and
log-periodic modulations along the whole chromosome. Hence scaling analyses of
distance series of DNA sequences have to consider the functional units of
bacterial chromosomes.


  Recently, it has been confirmed that long duplex DNA molecules with sizes
larger than several tens of kilo-base pairs (kbp), exhibit a discrete
conformational transition from an elongated coil state to a compact globule
state upon the addition of various kinds of chemical species that usually
induce DNA condensation. In this study, we performed a single-molecule
observation on a large DNA, Lambda ZAP II DNA (ca. 41 kbp), in a solution
containing RNA polymerase and substrates along with spermine, a tetravalent
cation, at different concentrations, by use of fluorescence staining of both
DNA and RNA. We found that transcription, or RNA production, is completely
inhibited in the compact state, but is actively performed in the unfolded coil
state. Such an all-or-none effect on transcriptional activity induced by the
discrete conformational transition of single DNA molecules is discussed in
relation to the mechanism of the regulation of large-scale genetic activity.


  Until now, design of the annual influenza vaccine has relied on phylogenetic
or whole-sequence comparisons of the viral coat proteins hemagglutinin and
neuraminidase, with vaccine effectiveness assumed to correlate monotonically to
the vaccine-influenza sequence difference. We use a theory from statistical
mechanics to quantify the non-monotonic immune response that results from
antigenic drift in the epitopes of the hemagglutinin and neuraminidase
proteins. The results explain the ineffectiveness of the 2003--2004 influenza
vaccine in the United States and provide an accurate measure by which to
optimize the effectiveness of future annual influenza vaccines.


  The rules that specify how the information contained in DNA codes amino
acids, is called "the genetic code". Using a simplified version of the Penna
nodel, we are using computer simulations to investigate the importance of the
genetic code and the number of amino acids in Nature on population dynamics. We
find that the genetic code is not a random pairing of codons to amino acids and
the number of amino acids in Nature is an optimum under mutations.


  We report first principles density functional calculations for
5,6-dihydroxyindole-2-carboxylic acid (DHICA) and several reduced forms. DHICA
and 5,6-dihydroxyindole (DHI) are believed to be the basic building blocks of
the eumelanins. Our results show that carboxylation has a significant effect on
the physical properties of the molecules. In particular, the relative
stabilities and the HOMO-LUMO gaps (calculated with the $\Delta$SCF method) of
the various redox forms are strongly affected. We predict that, in contrast to
DHI, the density of unpaired electrons, and hence the ESR signal, in DHICA is
negligibly small.


  Collective behavior of active elements inspired by mass of biological
organisms is addressed. Especially, two topics are focused on among amazing
behaviors performed by colony of ants. First, task allocation phenomena are
treated from the viewpoint of proportion regulation of population between
different states. Using a dynamical model consisting of elements and external
``stock materials'', adaptability against various disturbances is numerically
studied. In addition, a dynamical model for a colony ants interacting via two
kind of pheromones is studied, in which simulated ants, as a mass, are shown to
make an efficient foraging flexibly varying the foraging tactics according to
feeding schedules. Finally, experiments are performed with robots moving in
virtual pheromone fields simulated by CG and CCD camera feedback system. Trail
formation processes are demonstrated by this multi-robot system.


  We study the evolution of recombination using a microscopic model developed
within the frame of the theory of quantitative traits. Two components of
fitness are considered: a static one that describes adaptation to environmental
factors not related to the population itself, and a dynamic one that accounts
for interactions between organisms e.g. competition. We focus on the dynamics
of colonization of an empty niche. As competition is a function of the
population, selection pressure rapidly changes in time. The simulations show
that both in the case of flat and steep static fitness landscapes,
recombination provides a high velocity of movement in the phenotypic space thus
allowing recombinants to colonize the highest fitness regions earlier than non
recombinants that are often driven to extinction. The stabilizing effects of
competition and assortativity are also discussed. Finally, the analysis of
phase diagrams shows that competition is the key factor for the evolution of
recombination, while assortativity plays a significant role only in small
populations.


  The n-person Prisoner's Dilemma is a widely used model for populations where
individuals interact in groups. The evolutionary stability of populations has
been analysed in the literature for the case where mutations in the population
may be considered as isolated events. For this case, and assuming simple
trigger strategies and many iterations per game, we analyse the rate of
convergence to the evolutionarily stable populations. We find that for some
values of the payoff parameters of the Prisoner's Dilemma this rate is so low
that the assumption, that mutations in the population are infrequent on that
timescale, is unreasonable. Furthermore, the problem is compounded as the group
size is increased. In order to address this issue, we derive a deterministic
approximation of the evolutionary dynamics with explicit, stochastic mutation
processes, valid when the population size is large. We then analyse how the
evolutionary dynamics depends on the following factors: mutation rate, group
size, the value of the payoff parameters, and the structure of the initial
population. In order to carry out the simulations for groups of more than just
a few individuals, we derive an efficient way of calculating the fitness
values. We find that when the mutation rate per individual and generation is
very low, the dynamics is characterised by populations which are evolutionarily
stable. As the mutation rate is increased, other fixed points with a higher
degree of cooperation become stable. For some values of the payoff parameters,
the system is characterised by (apparently) stable limit cycles dominated by
cooperative behaviour. The parameter regions corresponding to high degree of
cooperation grow in size with the mutation rate, and in number with the group
size.


  Vibrational energy relaxation (VER) of a selected mode in cytochrome c
(hemeprotein) in vacuum is studied using two theoretical approaches: One is the
equilibrium simulation approach with quantum correction factors, and the other
is the reduced model approach which describes the protein as an ensemble of
normal modes coupled with nonlinear coupling elements. Both methods result in
estimates of VER time (sub ps) for a CD stretching mode in the protein at room
temperature, that are in accord with the experimental data of Romesberg's
group. The applicability of the two methods is examined through a discussion of
the validity of Fermi's golden rule on which the two methods are based.


  We study the medium effects on the selection of sequences in protein folding
by taking account of the surface potential in HP-model. Our analysis on the
proportion of H and P monomers in the sequences gives a direct interpretation
that the lowly designable structures possess small average gap. The numerical
calculation by means of our model exhibits that the surface potential enhances
the average gap of highly designable structures. It also shows that a most
stable structure may be no longer the most stable one if the medium parameters
changed.


  Dynamics of the polymorphic A<->B transitions in DNA is compared for two
polypurine sequences, poly(dA).poly(dT) and poly(dG).poly(dC), long known to
exhibit contrasting properties in experiments. In free molecular dynamics
simulations reversible transitions are induced by changing the size of a water
drop around DNA neutralized by Na ions. In poly(dG).poly(dC) the B<->A
transitions are easy, smooth and perfectly reversible. In contrast, a B->A
transition in poly(dA).poly(dT) dodecamer fragment could not be obtained even
though its A-form is stable under low hydration. Normal B->A transitions are
observed, however, in long poly(dA).poly(dT) stretches flanked by GC pairs. An
intermediate range of hydration numbers is identified where opposite
transitions are observed in the two dodecamer fragments, namely, A->B in
poly(dA).poly(dT) and B->A in poly(dG).poly(dC). With hydration numbers close
to the stability limit of the B-form, the two sequences exhibit qualitatively
different counterion distributions, with a characteristic accumulation of Na
ions next to the opening of the minor groove in poly(dA).poly(dT). This
difference can explain an increased persistence of poly(dA).poly(dT) DNA
towards A-form in crystalline and amorphous fibers as compared to solution
conditions. The good overall agreement with experimental data corroborates the
general role of the electrostatic condensation mechanism in the A/B
polymorphism in DNA.


  ... By additionally considering a two-mutation model for carcinogenesis and
indefinite proliferation by the activation of telomerase, we demonstrate that
the risk of dying doe to cancer can outweigh the positive effect of longer
telomeres on the longevity.


  Recent studies of in vitro evolution of DNA via protein binding indicate that
the evolution behavior is qualitatively different in different parameter
regimes. I here present a general theory that is valid for a wide range of
parameters, and which reproduces and extends previous results. Specifically,
the mean-field theory of a general translation-invariant model can be reduced
to the basic diffusion equation with a dynamic boundary condition. The simple
analytical form yields both quantitatively accurate predictions and valuable
insight into the principles involved. In particular, I introduce a cutoff
criterion for finite populations that illustrates both of these qualities.


  The 16-22 amino acid fragment of the beta-amyloid peptide associated with the
Alzheimer's disease, Abeta, is capable of forming amyloid fibrils. Here we
study the aggregation mechanism of Abeta(16-22) peptides by unbiased
thermodynamic simulations at the atomic level for systems of one, three and six
Abeta(16-22) peptides. We find that the isolated Abeta(16-22) peptide is mainly
a random coil in the sense that both the alpha-helix and beta-strand contents
are low, whereas the three- and six-chain systems form aggregated structures
with a high beta-sheet content. Furthermore, in agreement with experiments on
Abeta(16-22) fibrils, we find that large parallel beta-sheets are unlikely to
form. For the six-chain system, the aggregated structures can have many
different shapes, but certain particularly stable shapes can be identified.


  In E. coli, accurate cell division depends upon the oscillation of Min
proteins from pole to pole. We provide a model for the polar localization of
MinD based only on diffusion, a delay for nucleotide exchange, and different
rates of attachment to the bare membrane and the occupied membrane. We derive
analytically the probability density, and correspondingly the length scale, for
MinD attachment zones. Our simple analytical model illustrates the processes
giving rise to the observed localization of cellular MinD zones.


  The charges in live cells interact with or produce electric fields, which
results in enormous dielectric responses, flexoelectricity, and related
phenomena. Here we report on a contraction of schizosacchraoymces pombe
(fission yeast) cells induced by magnetic fields, as observed using a phase
sensitive projection image technique. Unlike electric fields, magnetic fields
only act on moving charges. The observed behavior is quite remarkable, and may
result from a contractile Lorentz force acting on diamagnetic screening
currents. This would indicate extremely high intracellular charge mobilities.
Besides, we observed a large electro - optical response from fission yeast
cells.


  We identified latent periodicity in catalytic domains of approximately 85% of
serine/threonine and tyrosine protein kinases. Similar results were obtained
for other 22 protein domains. We also designed the method of noise
decomposition, which is aimed to distinguish between different periodicity
types of the same period length. The method is to be used in conjunction with
the cyclic profile alignment, and this combination is able to reveal
structure-related or function-related patterns of latent periodicity. Possible
origins of the periodic structure of protein kinase active sites are discussed.
Summarizing, we presume that latent periodicity is the common property of many
catalytic protein domains.


  We found latent periodicity of 150 protein families now. We suppose that
latent periodicity can determine a spectrum of resonance oscillations in
proteins.


  Small Ubiquitin-related modifier (SUMO) proteins are widely expressed in
eukaryotic cells, which are reversibly coupled to their substrates by motif
recognition, called sumoylation. Two interesting questions are 1) how many
potential SUMO substrates may be included in mammalian proteomes, such as human
and mouse, 2) and given a SUMO substrate, can we recognize its sumoylation
sites? To answer these two questions, previous prediction systems of SUMO
substrates mainly adopted the pattern recognition methods, which could get high
sensitivity with relatively too many potential false positives. So we use
phylogenetic conservation between mouse and human to reduce the number of
potential false positives.


  We present a simple theory that uses thermodynamic parameters to predict the
probability that a protein retains the wildtype structure after one or more
random amino acid substitutions. Our theory predicts that for large numbers of
substitutions the probability that a protein retains its structure will decline
exponentially with the number of substitutions, with the severity of this
decline determined by properties of the structure. Our theory also predicts
that a protein can gain extra robustness to the first few substitutions by
increasing its thermodynamic stability. We validate our theory with simulations
on lattice protein models and by showing that it quantitatively predicts
previously published experimental measurements on subtilisin and our own
measurements on variants of TEM1 beta-lactamase. Our work unifies observations
about the clustering of functional proteins in sequence space, and provides a
basis for interpreting the response of proteins to substitutions in protein
engineering applications.


  Proteins have regular tertiary structures but irregular amino acid sequences.
This made it very difficult to decode the structural information in the protein
sequences. Here we demonstrate that many small alpha protein domains have
hidden sequence symmetries characteristic of their pseudo-symmetric tertiary
structures. We also present a modified method of recurrent plot to reveal this
kind of the hidden sequence symmetry. The results may enable us understand
parts of the relations between protein sequences and their tertiary structures,
i.e, how the primary sequence of a protein determines its tertiary structure.


  Motivation: Predicting the metastatic potential of primary malignant tissues
has direct bearing on choice of therapy. Several microarray studies yielded
gene sets whose expression profiles successfully predicted survival (Ramaswamy
et al 2003; Sorlie et al 2001; van't Veer et al 2003). Nevertheless, the
overlap between these gene sets is almost zero. Such small overlaps were
observed also in other complex diseases (Lossos et al 2003; Miklos and Maleszka
2004), and the variables that could account for the differences had evoked a
wide interest. One of the main open questions in this context is whether the
disparity can be attributed only to trivial reasons such as different
technologies, different patients and different types of analysis. Results: To
answer this question we concentrated on one single breast cancer dataset, and
analyzed it by one single method, the one which was used by van't Veer et al to
produce a set of outcome predictive genes. We showed that in fact the resulting
set of genes is not unique; it is strongly influenced by the subset of patients
used for gene selection. Many equally predictive lists could have been produced
from the same analysis. Three main properties of the data explain this
sensitivity: (a) many genes are correlated with survival; (b) the differences
between these correlations are small; (c) the correlations fluctuate strongly
when measured over different subsets of patients. A possible biological
explanation for these properties is discussed.


  Summary: The F2CS server provides access to the software, F2CS2.00, that
implements an automated prediction method of SCOP and CATH classifications of
proteins, based on their FSSP Z-scores (Getz et al., 2002), Availability: Free,
at http://www.weizmann.ac.il/physics/complex/compphys/f2cs/. Contact:
eytan.domany@weizmann.ac.il Supplementary information: The site contains links
to additional figures and tables.


  We explore some aspects of the relationship between biological evolution
processes and the mathematical theory of records. For Eigen's quasispecies
model with an uncorrelated fitness landscape, we show that the evolutionary
trajectories traced out by a population initially localized at a randomly
chosen point in sequence space can be described in close analogy to record
dynamics, with two complications. First, the increasing number of genotypes
that become available with increasing distance from the starting point implies
that fitness records are more frequent than for the standard case of
independent, identically distributed random variables. Second, fitness records
can be bypassed, which strongly reduces the number of genotypes that take part
in an evolutionary trajectory. For exponential and Gaussian fitness
distributions, this number scales with sequence length $N$ as $\sqrt{N}$, and
it is of order unity for distributions with a power law tail. This is in strong
contrast to the number of records, which is of order $N$ for any fitness
distribution.


  Immune cells coordinate their efforts for the correct and efficient
functioning of the immune system (IS). Each cell type plays a distinct role and
communicates with other cell types through mediators such as cytokines,
chemokines and hormones, among others, that are crucial for the functioning of
the IS and its fine tuning. Nevertheless, a quantitative analysis of the
topological properties of an immunological network involving this complex
interchange of mediators among immune cells is still lacking. Here we present a
method for quantifying the relevance of different mediators in the immune
network, which exploits a definition of centrality based on the concept of
efficient communication. The analysis, applied to the human immune system,
indicates that its mediators significantly differ in their network relevance.
We found that cytokines involved in innate immunity and inflammation and some
hormones rank highest in the network, revealing that the most prominent
mediators of the IS are molecules involved in these ancestral types of defence
mechanisms highly integrated with the adaptive immune response, and at the
interplay among the nervous, the endocrine and the immune systems.


  The replica exchange molecular dynamics (REMD) approach is applied to four
oligomeric peptide systems. At physiologically relevant temperature values REMD
samples conformation space and aggregation transitions more efficiently than
constant temperature molecular dynamics (CTMD). During the aggregation process
the energetic and structural properties are essentially the same in REMD and
CTMD. A condensation stage toward disordered aggregates precedes the
$\beta$-sheet formation. Two order parameters, borrowed from anisotropic fluid
analysis, are used to monitor the aggregation process. The order parameters do
not depend on the peptide sequence and length and therefore allow to compare
the amyloidogenic propensity of different peptides.


  Biological systems rely on robust internal information processing: Survival
depends on highly reproducible dynamics of regulatory processes. Biological
information processing elements, however, are intrinsically noisy (genetic
switches, neurons, etc.). Such noise poses severe stability problems to system
behavior as it tends to desynchronize system dynamics (e.g. via fluctuating
response or transmission time of the elements). Synchronicity in parallel
information processing is not readily sustained in the absence of a central
clock. Here we analyze the influence of topology on synchronicity in networks
of autonomous noisy elements. In numerical and analytical studies we find a
clear distinction between non-reliable and reliable dynamical attractors,
depending on the topology of the circuit. In the reliable cases, synchronicity
is sustained, while in the unreliable scenario, fluctuating responses of single
elements can gradually desynchronize the system, leading to non-reproducible
behavior. We find that the fraction of reliable dynamical attractors strongly
correlates with the underlying circuitry. Our model suggests that the observed
motif structure of biological signaling networks is shaped by the biological
requirement for reproducibility of attractors.


  We define the ``shift-match number'' for a binary string and we compute the
probability of occurrence of a given string as a subsequence in longer strings
in terms of its shift-match number. We thus prove that the string matching
probabilities depend not only on the length of shorter strings, but also on the
equivalence class of the shorter string determined by its shift-match number.


  We construct two models of discrete-time replicator dynamics with time delay.
In the social-type model, players imitate opponents taking into account average
payoffs of games played some units of time ago. In the biological-type model,
new players are born from parents who played in the past. We consider
two-player games with two strategies and a unique mixed evolutionarily stable
strategy. We show that in the first type of dynamics, it is asymptotically
stable for small time delays and becomes unstable for big ones when the
population oscillates around its stationary state. In the second type of
dynamics, however, evolutionarily stable strategy is asymptotically stable for
any size of a time delay.


  Activated processes such as protein unfolding are highly sensitive to
heterogeneity in the environment. We study a highly simplified model of a
protein in a random heterogeneous environment, a model of the in vivo
environment. It is found that if the heterogeneity is sufficiently large the
total rate of the process is essentially a random variable; this may be the
cause of the species-to-species variability in the rate of prion protein
conversion found by Deleault et al. [Nature, 425 (2003) 717].


  We discuss stochastic dynamics of populations of individuals playing games.
Our models possess two evolutionarily stable strategies: an efficient one,
where a population is in a state with the maximal payoff (fitness) and a
risk-dominant one, where players are averse to risks. We assume that
individuals play with randomly chosen opponents (they do not play against
average strategies as in the standard replicator dynamics). We show that the
long-run behavior of a population depends on its size and the mutation level.


  Intermittent density fluctuations of nucleotide molecules (adenine, guanine,
cytosine and thymine) along DNA sequences are studied in the framework of a
hierarchical structure (HS) model originally proposed for the study of fully
developed turbulence [She and Leque, Phys. Rev. Lett. 72}, 336 (1994)]. Large
scale (10^3 < \ell < 10^5 bp) base density fluctuation is shown to satisfy the
HS similarity. The derived values of a HS parameter $\beta$ from a large number
of genome data (including Bacteria, Archaea, human chromosomes and viruses)
characterize different biological properties such as strand symmetry,
phylogenetic relations and horizontal gene transfer. It is suggested that the
HS analysis offers a useful quantitative description for heterogeneity,
sequence complexity and large scale structures of genomes.


  We provide a classification of symmetric three-player games with two
strategies and investigate evolutionary and asymptotic stability (in the
replicator dynamics) of their Nash equilibria. We discuss similarities and
differences between two-player and multi-player games. In particular, we
construct examples which exhibit a novel behavior not found in two-player
games.


  We present a statistical approach to protein structure by introducing a
representation of protein folds based on simple observables defined as
frequencies of oriented cycles in contact graphs. Motivated by the idea that
these cycles may form a code of the entire protein structure we investigate its
Shannon entropy. The latter shows a characteristic transitory behavior when
traced over different values of the geometric threshold $t$ which defines
protein residues in contact. To account for this observation, we propose a
non-linear mechanical model-- formulated as a Hamiltonian system in which $t$
is regarded as a continuous parameter-- that allows to identify the behavior of
the Shannon entropy as a first-order phase transition between a disordered and
an ordered phase of the proposed mechanical system. The transition itself
reflects the formation of protein structure when represented in terms of
contact graph cycles, and it is identified as an example of a fragmentation
transition known from several other statistical systems without a proper
thermodynamic limit. Some interesting implications follow from our model
including chirality, broken $t$-symmetry and one-dimensionality. Although
defined from purely structural considerations, we further show that for native
protein structures these observables follow some of the quantitative rules that
are typically valid for amino acid types in protein sequences. Moreover, this
relation between polypeptide structures and their amino acid sequences suggests
a specific protein design alphabet.


  Water molecules and molecular chaperones efficiently help the protein folding
process. Here we describe their action in the context of the energy and
topological networks of proteins. In energy terms water and chaperones were
suggested to decrease the activation energy between various local energy minima
smoothing the energy landscape, rescuing misfolded proteins from conformational
traps and stabilizing their native structure. In kinetic terms water and
chaperones may make the punctuated equilibrium of conformational changes less
punctuated and help protein relaxation. Finally, water and chaperones may help
the convergence of multiple energy landscapes during protein-macromolecule
interactions. We also discuss the possibility of the introduction of protein
games to narrow the multitude of the energy landscapes when a protein binds to
another macromolecule. Both water and chaperones provide a diffuse set of
rapidly fluctuating weak links (low affinity and low probability interactions),
which allow the generalization of all these statements to a multitude of
networks.


  In the context of complex systems and, particularly, of protein folding, a
physically meaningful distance is defined which allows to make useful
statistical statements about the way in which energy differences are modified
when two different instances of the same potential-energy function are used.
When the two instances arise from the fact that different algorithms or
different approximations are used, the distance herein defined may be used to
evaluate the relative accuracy of the two methods. When the difference is due
to a change in the free parameters of which the potential depends on, the
distance can be used to quantify, in each region of parameter space, the
robustness of the modeling to such a change and this, in turn, may be used to
assess the significance of a parameters' fit. Both cases are illustrated with a
practical example: the study of the Poisson-based solvation energy in the
Trp-Cage protein (PDB code 1L2Y).


  As mathematical model for the evolutionary equations of species the
masterequation is choiced. Two formulations will be demonstrated to include the
changes of parameters into the masterequation - that is, on the one hand, the
formation of a second masterequation for the development of parameters, and, on
the other hand, the use of the Wigner-distribution to describe the development
of parameters. Moreover, the Wigner-distribution is used to describe morphic
fields and involved in the theorie of selforganization.


  The genes of eukaryotes are characterized by protein coding fragments, the
exons, interrupted by introns, i.e. stretches of DNA which do not carry any
useful information for the protein synthesis. We have analyzed the melting
behavior of randomly selected human cDNA sequences obtained from the genomic
DNA by removing all introns. A clear correspondence is observed between exons
and melting domains. This finding may provide new insights in the physical
mechanisms underlying the evolution of genes.


  We describe an new algorithm for visualizing an alignment of biological
sequences according to a probabilistic model of evolution. The resulting data
array is readily interpreted by the human eye and amenable to digital image
techniques. We present examples using mRNA sequences from mouse and rat: three
cytochromes, Rattus norvegicus Cyp2a1, Cyp2a2, (Medline: 90212624) and Mus
musculus Cyp2a12 (Medline: 93249380); and two zinc finger proteins, Mus
musculus zfp111 and zfp235 (Medline: 22683274). The underlying evolutionary
model is derived from one proposed by Thorne, Kishino, and Felsenstein and
improved by Hein and others. The demonstration implementation aligns two
sequences using time and memory quadratic in the mean sequence length. The
algorithm is extensible, after Hein, to multiple sequences. We mention a basic
method to reduce time and memory demands.


  High-gain nonlinear observers occur in the nonlinear automatic control theory
and are in standard usage in chemical engineering processes. We apply such a
type of analysis in the context of a very simple one-gene regulation circuit.
In general, an observer combines an analytical differential-equation-based
model with partial measurement of the system in order to estimate the
non-measured state variables. We use one of the simplest observers, that of
Gauthier et al., which is a copy of the original system plus a correction term
which is easy to calculate. For the illustration of this procedure, we employ a
biological model, recently adapted from Goodwin's old book by De Jong, in which
one plays with the dynamics of the concentrations of the messenger RNA coding
for a given protein, the protein itself, and a single metabolite. Using the
observer instead of the metabolite, it is possible to rebuild the non-measured
concentrations of the mRNA and the protein


  Protein-fragment seqlets typically feature about 10 amino acid residue
positions that are fixed to within conservative substitutions but usually
separated by a number of prescribed gaps with arbitrary residue content. By
quantifying a general amino acid residue sequence in terms of the associated
codon number sequence, we have found a precise modular Fibonacci sequence in a
continuous gap-free 10-residue seqlet with either 3 or 4 conservative amino
acid substitutions. This modular Fibonacci sequence is genuinely biophysical,
for it occurs nine times in the SWISS-Prot/TrEMBL database of natural proteins.


  The mammalian cortex is divided into architectonic and functionally distinct
areas. There is growing experimental evidence that their emergence and
development is controlled by both epigenetic and genetic factors. The latter
were recently implicated as dominating the early cortical area specification.
In this paper, we present a theoretical model that explicitly considers the
genetic factors and that is able to explain several sets of experiments on
cortical area regulation involving transcription factors Emx2 and Pax6, and
fibroblast growth factor FGF8. The model consists of the dynamics of thalamo-
cortical connections modulated by signaling molecules that are regulated
genetically, and by axonal competition for neocortical space. The model can
make predictions and provides a basic mathematical framework for the early
development of the thalamo-cortical connections and area patterning that can be
further refined as more experimental facts become known.


  We study by mean-field analysis and stochastic simulations chemical models
for genetic toggle switches formed from pairs of genes that mutually repress
each other. In order to determine the stability of the genetic switches, we
make a connection with reactive flux theory and transition state theory. The
switch stability is characterised by a well defined lifetime $\tau$. We find
that $\tau$ grows exponentially with the mean number $\Nmean$ of transcription
factor molecules involved in the switching. In the regime accessible to direct
numerical simulations, the growth law is well characterised by
$\tau\sim\Nmean{}^{\alpha}\exp(b\Nmean)$, where $\alpha$ and $b$ are
parameters. The switch stability is decreased by phenomena that increase the
noise in gene expression, such as the production of multiple copies of a
protein from a single mRNA transcript (shot noise), and fluctuations in the
number of proteins produced per transcript. However, robustness against
biochemical noise can be drastically enhanced by arranging the transcription
factor binding domains on the DNA such that competing transcription factors
mutually exclude each other on the DNA. We also elucidate the origin of the
enhanced stability of the exclusive switch with respect to that of the general
switch: while the kinetic prefactor is roughly the same for both switches, the
`barrier' for flipping the switch is significantly higher for the exclusive
switch than for the general switch.


  We investigate 738 complete genomes of viruses to detect the presence of
short inverted repeats. The number of inverted repeats found is compared with
the prediction obtained for a Bernoullian and for a Markovian control model. We
find as a statistical regularity that the number of observed inverted repeats
is often greater than the one expected in terms of a Bernoullian or Markovian
model in several of the viruses and in almost all those with a genome longer
than 30,000 bp.


  By driven to extinction species less or poorly adapted, the Darwinian
evolutionary theory is intrinsically an optimization theory. We investigate two
optimization algorithms with such evolutionary characteristics: the Bak-Sneppen
and the Extremal Optimization. By comparing their mean fitness in the steady
state regime, we conclude that the Bak-Sneppen dynamics is more efficient than
the Extremal Optimization if the parameter $\tau$ is in the interval
$[0,0.86]$. The determination of the spatial correlation and the probability
distribution of the avalanches show that the Extremal Optimization dynamics
does not lead the system into a critical self-organized state. Trough a
discrete form of the Bak-Sneppen model we argument that biodiversity is an
essential requisite to preserve the self-organized criticality.


  With the sponsorship of ``Fundacio La Caixa'' we met in Barcelona, November
21st and 22nd, to analyze the reasons why, after the completion of the human
genome sequence, the identification all protein coding genes and their variants
remains a distant goal. Here we report on our discussions and summarize some of
the major challenges that need to be overcome in order to complete the human
gene catalog.


  Recent progress has clarified many features of the global architecture of
biological metabolic networks, which have highly organized and optimized
tolerances and tradeoffs (HOT) for functional requirements of flexibility,
efficiency, robustness, and evolvability, with constraints on conservation of
energy, redox, and many small moieties. One consequence of this architecture is
a highly structured modularity that is self-dissimilar and scale-rich, with
extremes in low and high variability, including power laws, in both metabolite
and reaction degree distributions. This paper illustrates these features using
the well-understood stoichiometry of metabolic networks in bacteria, and a
simple model of an abstract metabolism.


  It is proposed that the co-expression of statistically significant motifs
among the sequences of a proteome is a phylogenetic trait. From the
co-expression matrix of such motifs in a group of prokaryotic proteomes a
suitable definition of a phylogenetic distance is introduced and the
corresponding distance matrix between proteomes is constructed. From the
distance matrix a phylogenetic tree is inferred, following a standard
procedure. It compares well with a reference tree deduced from a distance
matrix obtained from the alignment of ribosomal RNA sequences. Our results are
consistent with the hypothesis that biological evolution manifests itself with
a modulation of basic correlations between shared peptides of short length,
present in protein sequences. Moreover, the simple procedure we propose
reconfirms that it is possible, sampling entire proteomes, to average the
effects of lateral gene transfer and infer reasonable phylogenies.


  Finding analytically the statistics of the longest common subsequence (LCS)
of a pair of random sequences drawn from c alphabets is a challenging problem
in computational evolutionary biology. We present exact asymptotic results for
the distribution of the LCS in a simpler, yet nontrivial, variant of the
original model called the Bernoulli matching (BM) model which reduces to the
original model in the large c limit. We show that in the BM model, for all c,
the distribution of the asymptotic length of the LCS, suitably scaled, is
identical to the Tracy-Widom distribution of the largest eigenvalue of a random
matrix whose entries are drawn from a Gaussian unitary ensemble. In particular,
in the large c limit, this provides an exact expression for the asymptotic
length distribution in the original LCS problem.


  Selection pressures on proteins are usually measured by comparing homologous
nucleotide sequences (Zuckerkandl and Pauling 1965). Recently we introduced a
novel method, termed `volatility', to estimate selection pressures on protein
sequences from their synonymous codon usage (Plotkin and Dushoff 2003, Plotkin
et al 2004a). Here we provide a theoretical foundation for this approach. We
derive the expected frequencies of synonymous codons as a function of the
strength of selection, the mutation rate, and the effective population size. We
analyze the conditions under which we can expect to draw inferences from biased
codon usage, and we estimate the time scales required to establish and maintain
such a signal. Our results indicate that, over a broad range of parameters,
synonymous codon usage can reliably distinguish between negative selection,
positive selection, and neutrality. While the power of volatility to detect
negative selection depends on the population size, there is no such dependence
for the detection of positive selection. Furthermore, we show that phenomena
such as transient hyper-mutators in microbes can improve the power of
volatility to detect negative selection, even when the typical observed neutral
site heterozygosity is low.


  We study the role of recombination, as practiced by genetically-competent
bacteria, in speeding up Darwinian evolution. This is done by adding a new
process to a previously-studied Markov model of evolution on a smooth fitness
landscape; this new process allows alleles to be exchanged with those in the
surrounding medium. Our results, both numerical and analytic, indicate that for
a wide range of intermediate population sizes, recombination dramatically
speeds up the evolutionary advance.