  Protein functions in cells may be activated or modified by the attachment of
several kinds of chemical groups. While protein phosphorylation, i.e. the
attachment of a phosphoryl (PO$_3^-$) group, is the most studied form of
protein modification, and is known to regulate the functions of many proteins,
protein behavior can also be modified by nitrosylation, acetylation,
methylation, etc. A protein can have multiple modification sites, and display
some form of transition only when enough sites are modified. In a previous
paper we have modeled the generic equilibrium properties of multisite protein
modification (R.Chignola, C. Dalla Pellegrina, A. Del Fabbro, E.Milotti,
Physica A {\bf 371}, 463 (2006)) and we have shown that it can account both for
sharp, robust thresholds and for information transfer between processes with
widely separated timescales. Here we use the same concepts to expand that
analysis starting from a dynamical description of multisite modification: we
give analytical results for the basic dynamics and numerical results in an
example where the modification chain is cascaded with a Michaelis-Menten step.
We modify the dynamics and analyze an example with realistic
phosphorylation/dephosphorylation steps, and give numerical evidence of the
independence of the allosteric effect from the details of the
attachment-detachment processes. We conclude that multisite protein
modification is dynamically equivalent to the classic allosteric effect.


  In eukaryotic cells, many motor proteins can move simultaneously on a single
microtubule track. This leads to interesting collective phenomena like jamming.
Recently we reported ({\it Phys. Rev. Lett. {\bf 95}, 118101 (2005)}) a
lattice-gas model which describes traffic of unconventional (single-headed)
kinesins KIF1A. Here we generalize this model, introducing a novel interaction
parameter $c$, to account for an interesting mechano-chemical process which has
not been considered in any earlier model. We have been able to extract all the
parameters of the model, except $c$, from experimentally measured quantities.
In contrast to earlier models of intra-cellular molecular motor traffic, our
model assigns distinct ``chemical'' (or, conformational) states to each kinesin
to account for the hydrolysis of ATP, the chemical fuel of the motor. Our model
makes experimentally testable theoretical predictions. We determine the phase
diagram of the model in planes spanned by experimentally controllable
parameters, namely, the concentrations of kinesins and ATP. Furthermore, the
phase-separated regime is studied in some detail using analytical methods and
simulations to determine e.g. the position of shocks. Comparison of our
theoretical predictions with experimental results is expected to elucidate the
nature of the mechano-chemical process captured by the parameter $c$.


  To study the fluctuations and dynamics in chemical reaction processes,
stochastic differential equations based on the rate equation involving chemical
concentrations are often adopted. When the number of molecules is very small,
however, the discreteness in the number of molecules cannot be neglected since
the number of molecules must be an integer. This discreteness can be important
in biochemical reactions, where the total number of molecules is not
significantly larger than the number of chemical species. To elucidate the
effects of such discreteness, we study autocatalytic reaction systems
comprising several chemical species through stochastic particle simulations.
The generation of novel states is observed; it is caused by the extinction of
some molecular species due to the discreteness in their number. We demonstrate
that the reaction dynamics are switched by a single molecule, which leads to
the reconstruction of the acting network structure. We also show the strong
dependence of the chemical concentrations on the system size, which is caused
by transitions to discreteness-induced novel states.


  We study the evolutionary Prisoner's Dilemma on two social networks obtained
from actual relational data. We find very different cooperation levels on each
of them that can not be easily understood in terms of global statistical
properties of both networks. We claim that the result can be understood at the
mesoscopic scale, by studying the community structure of the networks. We
explain the dependence of the cooperation level on the temptation parameter in
terms of the internal structure of the communities and their interconnections.
We then test our results on community-structured, specifically designed
artificial networks, finding perfect agreement with the observations in the
real networks. Our results support the conclusion that studies of evolutionary
games on model networks and their interpretation in terms of global properties
may not be sufficient to study specific, real social systems. In addition, the
community perspective may be helpful to interpret the origin and behavior of
existing networks as well as to design structures that show resilient
cooperative behavior.


  Multisite protein modification is a ubiquitous mechanism utilized by cells to
control protein functions. We have recently proposed a dynamical description of
multisite protein modification which embodies all the essential features of the
process (E. Milotti, A. Del Fabbro, C. Dalla Pellegrina, and R. Chignola,
Physica A, in press), and we have used this model to analyze the stability and
the time-scales of this mechanism. The same model can be used to understand how
the system responds to stimuli: here we show that it displays
frequency-dependent dynamical hysteresis. This behavior closely parallels --
with the due differences -- what is observed in magnetic systems. By selecting
model parameters that span the known biological ranges, we find that the
frequency-dependent features cover the band of the observed oscillations of
molecular intracellular signals, and this suggests that this mechanism may have
an important role in cellular information processing.


  Proteins are polymers of amino acids. These macromolecules are synthesized by
intracellular machines called ribosomes. Although the experimental
investigation of protein synthesis has been a traditional area of research in
molecular cell biology, important quantitative models of protein synthesis have
been reported in research journals devoted to statistical physics and related
interdisciplinary topics. From the perspective of a physicist, protein
synthesis is the classical transport of interacting ribosomes on a messenger
RNA (mRNA) template that dictates the sequence of the amino acids on the
protein. We discuss appropriate simplification of the models and methods. In
particular, we develop and analyze a simple toy model using some elementary
techniques of non-equilibrium statistical mechanics and predict the average
rate of protein synthesis and the spatial organization of the ribosomes in the
steady state.


  We propose the n-clique network as a powerful tool for understanding global
structures of combined highly-interconnected subgraphs, and provide theoretical
predictions for statistical properties of the n-clique networks embedded in a
complex network using the degree distribution and the clustering spectrum.
Furthermore, using our theoretical predictions, we find that the statistical
properties are invariant between 3-clique networks and original networks for
several observable real-world networks with the scale-free connectivity and the
hierarchical modularity. The result implies that structural properties are
identical between the 3-clique networks and the original networks.


  We introduce a toy model of the "rat race" in which individuals try to better
themselves relative to the rest of the population. An individual is
characterized by a real-valued fitness and each advances at a constant rate by
an amount that depends on its standing in the population. The leader advances
to remain ahead of its nearest neighbor, while all others advance by an amount
that is set by the distance to the leader. A rich dynamics occurs as a function
of the mean jump size of the trailing particles. For small jumps, the leader
maintains its position, while for large jumps, there are long periods of stasis
that are punctuated by episodes of explosive advancement and many lead changes.
Intermediate to these two regimes, agents reach a common fitness and evolution
grinds to a halt.


  Allosteric regulation is often viewed as thermodynamic in nature. However
protein internal motions during an enzymatic reaction cycle can be slow hopping
processes over numerous potential barriers. We propose that regulating
molecules may function by modifying the nonequilibrium protein dynamics. The
theory predicts that an enzyme under the new mechanism has different
temperature dependence, waiting time distribution of the turnover cycle, and
dynamic fluctuation patterns with and without effector. Experimental tests of
the theory are proposed.


  Due to the time-scale limitations of all-atom simulation of proteins, there
has been substantial interest in coarse-grained approaches. Some methods, like
"Resolution Exchange," [E. Lyman et al., Phys. Rev. Lett. 96, 028105 (2006)]
can accelerate canonical all-atom sampling, but require properly distributed
coarse ensembles. We therefore demonstrate that full sampling can indeed be
achieved in a sufficiently simplified protein model, as verified by a recently
developed convergence analysis. The model accounts for protein backbone
geometry in that rigid peptide planes rotate according to atomistically defined
dihedral angles, but there are only two degrees of freedom (phi and psi
dihedrals) per residue. Our convergence analysis indicates that small proteins
(up to 89 residues in our tests) can be simulated for more than 50 "structural
decorrelation times" in less than a week on a single processor. We show that
the fluctuation behavior is reasonable, as well as discussing applications,
limitations, and extensions of the model.


  Boolean networks may be viewed as idealizations of biological genetic
networks, where each node is represented by an on-off switch which is a
function of the binary output from some other nodes. We evolve connectivity in
a single Boolean network, and demonstrate how the sole requirement of
sequential matching of attractors may open for an evolution that exhibits
punctuated equilibrium.


  In a recent paper, Newman surveys the literature on power law spectra in
evolution, self-organised criticality and presents a model of his own to arrive
at a conclusion that self-organised criticality is not necessary for evolution.
Not only did he miss a key model (Ecolab) that has a clear self-organised
critical mechanism, but also Newman's model exhibits the same mechanism that
gives rise to power law behaviour as does Ecolab. Newman's model is, in fact, a
``mean field'' approximation of a self-organised critical system. In this
paper, I have also implemented Newman's model using the Ecolab software,
removing the restriction that the number of species remains constant. It turns
out that the requirement of constant species number is non-trivial, leading to
a global coupling between species that is similar in effect to the species
interactions seen in Ecolab. In fact, the model must self-organise to a state
where the long time average of speciations balances that of the extinctions,
otherwise the system either collapses or explodes. In view of this, Newman's
model does not provide the hoped-for counter example to the presence of
self-organised criticality in evolution, but does provide a simple, almost
analytic model that can used to understand more intricate models such as
Ecolab.


  Internal coordinate molecular dynamics (ICMD) is an efficient method for
studying biopolymers, but it is readily applicable only to molecules with tree
topologies, that is with no internal flexible rings. Common examples violating
this condition are prolines and loops closed by S-S bridges in proteins. The
most important such case, however, is nucleic acids because the flexibility of
the furanose rings always plays an essential role in conformational transitions
both in DNA and RNA. There are a few long-known theoretical approaches to this
problem, but, in practice, rings with fixed bond lengths are closed by adding
appropriate harmonic distance restraints, which is not always acceptable
especially in dynamics. This paper tries to overcome this handicap of ICMD and
proposes a rational strategy which results in practical numerical algorithms.
It gives a unified analytical treatment which shows that this problem is very
close to the difficulties encountered by the method of constraints in Cartesian
coordinate dynamics, and certain ideas of the latter appear helpful in the
context of ICMD. The method is affordable for large molecules and generally
applicable to all kinds of rings. A specific implementation for five-membered
rings is described and tested for a proline-rich polypeptide and a decamer DNA
duplex. In both cases conditions are found which make possible time steps
around 10 fsec in ICMD calculations.


  This paper considers the effect of electrostatics on the stability of a
charged membrane. We show that at low ionic strength and high surface charge
density, repulsion between charges on the membrane renders it unstable to the
formation of holes. A free straight edge is unstable to modulations with
wavelength longer than the Debye screening length. Hence at low ionic strength,
membranes will disintegrate into vesicles. We use these results to interpret
stable holes in red blood cell ghosts (Steck et. al., Science 168:255(1970)).


  Bacterial chemotaxis is controlled by the signalling of a cluster of
receptors. A cooperative model is presented, in which coupling between
neighbouring receptor dimers enhances the sensitivity with which stimuli can be
detected, without diminishing the range of chemoeffector concentration over
which chemotaxis can operate. Individual receptor dimers have two stable
conformational states: one active, one inactive. Noise gives rise to a
distribution between these states, with the probability influenced by ligand
binding, and also by the conformational states of adjacent receptor dimers. The
two-state model is solved, based on an equivalence with the Ising model in a
randomly distributed magnetic field. The model has only two effective
parameters, and unifies a number of experimental findings. According to the
value of the parameter comparing coupling and noise, the signal can be
arbitrarily sensitive to changes in the fraction of receptor dimers to which
ligand is bound. The counteracting effect of a change of methylation level is
mapped to an induced field in the Ising model. By returning the activity to the
pre-stimulus level, this adapts the receptor cluster to a new ambient
concentration of chemoeffector and ensures that a sensitive response can be
maintained over a wide range of concentrations.


  We present a so-called adaptive Ising model (AIM) to provide a unifying
explanation for sensitivity and perfect adaptation in bacterial chemotactic
signalling, based on coupling among receptor dimers. In an AIM, an external
field, representing ligand binding, is randomly applied to a fraction of spins,
representing the states of the receptor dimers, and there is a delayed negative
feedback from the spin value on the local field. This model is solved in an
adiabatic approach. If the feedback is slow and weak enough, as indeed in
chemotactic signalling, the system evolves through quasi-equilibrium states and
the ``magnetization'', representing the signal, always attenuates towards zero
and is always sensitive to a subsequent stimulus.


  Motivated by experiments in which a polynucleotide is driven through a
proteinaceous pore by an electric field, we study the diffusive motion of a
polymer threaded through a narrow channel with which it may have strong
interactions. We show that there is a range of polymer lengths in which the
system is approximately translationally invariant, and we develop a
coarse-grained description of this regime. From this description, general
features of the distribution of times for the polymer to pass through the pore
may be deduced. We also introduce a more microscopic model. This model provides
a physically reasonable scenario in which, as in experiments, the polymer's
speed depends sensitively on its chemical composition, and even on its
orientation in the channel. Finally, we point out that the experimental
distribution of times for the polymer to pass through the pore is much broader
than expected from simple estimates, and speculate on why this might be.


  The idea of adaptive walks on fitness landscapes as a means of studying
evolutionary processes on large time scales is extended to fitness landscapes
that are slowly changing over time. The influence of ruggedness and of the
amount of static fitness contributions are investigated for model landscapes
derived from Kauffman's $NK$ landscapes. Depending on the amount of static
fitness contributions in the landscape, the evolutionary dynamics can be
divided into a percolating and a non-percolating phase. In the percolating
phase, the walker performs a random walk over the regions of the landscape with
high fitness.


  Recently it has been found that stable and accurate molecular dynamics (MD)
of B-DNA duplexes can be obtained in relatively inexpensive computational
conditions with the bulk solvent represented implicitly, but the minor groove
filled with explicit water (J. Am. Chem. Soc. 1998, 120, 10928). The present
paper further explores these simulation conditions in order to understand the
main factors responsible for the observed surprisingly good agreement with
experimental data. It appears that in the case of the EcoRI dodecamer certain
sequence specific hydration patterns in the minor groove earlier known from
experimental data are formed spontaneously in the course of MD simulations. The
effect is reliably reproduced in several independent computational experiments
in different simulation conditions. With all major groove water removed,
closely similar results are obtained, with even better reproducibility. On the
other hand, without explicit hydration, metastable dynamics around a B-DNA like
state can be obtained which, however, only poorly compares with experimental
data. It appears, therefore, that a right-handed DNA helix with explicitly
hydrated minor groove is a minimal model system where the experimental
properties of B-DNA can be reproduced with non-trivial sequence-dependent
effects. Its small size makes possible virtually exhaustive sampling, which is
a major advantage with respect to alternative approaches. An appendix is
included with a correction to the implicit leapfrog integrator used in internal
coordinate MD.


  We study analytically the steady-state regime of a network of n error-prone
self-replicating templates forming an asymmetric hypercycle and its error tail.
We show that the existence of a master template with a higher non-catalyzed
self-replicative productivity, a, than the error tail ensures the stability of
chains in which m<n-1 templates coexist with the master species. The stability
of these chains against the error tail is guaranteed for catalytic coupling
strengths (K) of order of a. We find that the hypercycle becomes more stable
than the chains only for K of order of a2. Furthermore, we show that the
minimal replication accuracy per template needed to maintain the hypercycle,
the so-called error threshold, vanishes like sqrt(n/K) for large K and n<=4.


  We present a simple model based on a reaction-diffusion equation to explain
pattern formation in a multicellular bacterium (Streptomyces). We assume
competition for resources as the basic mechanism that leads to pattern
formation; in particular we are able to reproduce the spatial pattern formed by
bacterial aerial mycelium in case of growth in minimal (low resources) and
maximal (large resources) culture media.


  The Monte Carlo simulation of the dynamics of complex molecules produces
trajectories with a large number of different configurations to sample
configuration space. It is expected that these configurations can be classified
into a small number of conformations representing essential changes in the
shape of the molecule. We present a method to visualize these conformations by
point sets in the plane based on a geometrical distance measure between
individual configurations. It turns out that different conformations appear as
well-separated point sets. The method is further improved by performing a
cluster analysis of the data set. The point-cluster representation is used to
control a three-dimensional molecule viewer application to show individual
configurations and conformational changes. The extraction of essential
coordinates and visualization of molecular shape is discussed.


  We study self-replicating molecules under externally varying conditions.
Changing conditions such as temperature variations and/or alterations in the
environment's resource composition lead to both non-constant replication and
decay rates of the molecules. In general, therefore, molecular evolution takes
place in a dynamic rather than a static fitness landscape. We incorporate
dynamic replication and decay rates into the standard quasispecies theory of
molecular evolution, and show that for periodic time-dependencies, a system of
evolving molecules enters a limit cycle for $t\to\infty$. For fast periodic
changes, we show that molecules adapt to the time-averaged fitness landscape,
whereas for slow changes they track the variations in the landscape arbitrarily
closely. We derive a general approximation method that allows us to calculate
the attractor of time-periodic landscapes, and demonstrate using several
examples that the results of the approximation and the limiting cases of very
slow and very fast changes are in perfect agreement. We also discuss landscapes
with arbitrary time dependencies, and show that very fast changes again lead to
a system that adapts to the time-averaged landscape. Finally, we analyze the
dynamics of a finite population of molecules in a dynamic landscape, and
discuss its relation to the infinite population limit.


  Monotone systems constitute one of the most important classes of dynamical
systems used in mathematical biology modeling. The objective of this paper is
to extend the notion of monotonicity to systems with inputs and outputs, a
necessary first step in trying to understand interconnections, especially
including feedback loops, built up out of monotone components. Basic
definitions and theorems are provided, as well as an application of a theorem
regarding negative feedback loops to the study of a model of one of the cell's
most important subsystems (MAPK cascades) .


  This paper studies the emergence of multi-stability and hysteresis in those
systems that arise, under positive feedback, starting from monotone systems
with well-defined steady-state responses. Such feedback configurations appear
routinely in several fields of application, and especially in biology.
Characterizations of global stability behavior are stated in terms of easily
checkable graphical conditions. An example of a signaling cascade under
positive feedback is presented.


  This note provides a theorem showing, under suitable technical assumptions,
that if a system S adapts to a class of external signals U, in the sense of
egulation against disturbances or tracking signals in U, then S must ecessarily
contain a subsystem which is capable of generating all the signals in U. It is
not assumed that regulation is robust, nor is there a prior requirement for the
system to be partitioned into separate plant and controller components.
Instead, one assumes that a ``signal detection'' property holds. The result was
motivated by questions of adaptation in bacterial chemotaxis, but the general
mathematical principle is of wide applicability.


  Mobility of taxol inside microtubules was investigated using fluorescence
recovery after photobleaching (FRAP) on flow-aligned bundles. Bundles were made
of microtubules with either GMPCPP or GTP at the exchangeable site on the
tubulin dimer. Recovery times were sensitive to bundle thickness and packing,
indicating that taxol molecules are able to move laterally through the bundle.
The density of open binding sites along a microtubule was varied by controlling
the concentration of taxol in solution for GMPCPP samples. With > 63% sites
occupied, recovery times were independent of taxol concentration and,
therefore, inversely proportional to the microscopic dissociation rate,
k_{off}. It was found that 10*k_{off} (GMPCPP) ~ k_{off} (GTP), consistent
with, but not fully accounting for, the difference in equilibrium constants for
taxol on GMPCPP and GTP microtubules. With < 63% sites occupied, recovery times
decreased as ~ [Tax]^{-1/5} for both types of microtubules. We conclude that
the diffusion of taxol along the microtubule interior is hindered by rebinding
events when open sites are within ~7 nm of each other.


  Upon the covalent-bonding hybrid of the nitrogen atoms taken as a measure for
the structural regularity in nucleobases, it can be identified that the
internal relation within the 20 amino acids follows a cooperative
vector-in-space addition principle based on the spherical and rotational
symmetry of a quasi-28-gon (quasi-icosikaioctagon), with two evolutionary axes.


  We consider the regime in which the bands of the torsional acoustic (TA) and
the hydrogen-bond-stretch (HBS) modes of the DNA interpenetrate each other.
Within the framework of a model that accommodates the structure of the double
helix, we find the three-wave interaction between the TA- and the HBS-modes,
and show that microwave radiation could bring about torsional vibrations that
could serve as a pump mode for maintaining the HBS-one. Rayleigh's threshold
condition for the parametric resonance provides an estimate for the power
density of the mw-field necessary for generating the HBS-mode.


  New analyses of the organization of the genetic code system together with
their relation to the two classes of aminoacyl-tRNA synthetases are reported in
this work. A closer inspection revealed how the enzymes and the 20 amino acids
of the genetic code are intertwined on a polyhedron model. Complimentarily and
cooperative symmetry between class I and class II aminoacyl-tRNA synthetases
displayed by a 28-gon model are discussed, and we found that the two previously
suggested evolutionary axes of the genetic code overlap the three two-fold
symmetry axes within the two classes of aminoacyl-tRNA synthetases. Moreover,
it is identified that the amino-acid side-chain carbon-atom numbers (1, 3, 4
and 7) in the overwhelming majority of the amino acids recognized by each of
the two classes of aminoacyl-tRNA synthetases meet a mathematical relationship,
the Lucas series.


  Biophysicists and structural biologists increasingly acknowledge the role
played by the mechanical properties of macromolecules as a critical element in
many biological processes. This change has been brought about, in part, by the
advent of single molecule biophysics techniques that have made it possible to
exert piconewton forces on key macromolecules and observe their deformations at
nanometer length scales, as well as to observe the mechanical action of
macromolecules such as molecular motors. This has opened up immense
possibilities for a new generation of mechanical investigations that will
respond to such measurements in an attempt to develop a coherent theory for the
mechanical behavior of macromolecules under conditions where thermal and
chemical effects are on an equal footing with deterministic forces. This paper
presents an application of the principles of mechanics to the problem of DNA
packaging, one of the key events in the life cycle of bacterial viruses with
special reference to the nature of the internal forces that are built up during
the DNA packaging process.


  The representation of a biochemical network as a graph is the coarsest level
of description in cellular biochemistry. By studying the network structure one
can draw conclusions on the large scale organisation of the biochemical
processes. We describe methods how one can extract hierarchies of subnetworks,
how these can be interpreted and further deconstructed to find autonomous
subnetworks. The large-scale organisation we find is characterised by a tightly
connected core surrounded by increasingly loosely connected substrates.


  We present a simple model for the underlying structure of protein-protein
pairwise interaction graphs that is based on the way in which proteins attach
to each other in experiments such as yeast two-hybrid assays. We show that data
on the interactions of human proteins lend support to this model. The frequency
of the number of connections per protein under this model does not follow a
power law, in contrast to the reported behaviour of data from large scale yeast
two-hybrid screens of yeast protein-protein interactions. Sampling sub-graphs
from the underlying graphs generated with our model, in a way analogous to the
sampling performed in large scale yeast two-hybrid searches, gives degree
distributions that differ subtly from the power law and that fit the observed
data better than the power law itself. Our results show that the observation of
approximate power law behaviour in a sampled sub-graph does not imply that the
underlying graph follows a power law.


  We compare asynchronous vs. synchronous update of discrete dynamical networks
and find that a simple time delay in the nodes may induce a reproducible
deterministic dynamics even in the case of asynchronous update in random order.
In particular we observe that the dynamics under synchronous parallel update
can be reproduced accurately under random asynchronous serial update for a
large class of networks. This mechanism points at a possible general principle
of how computation in gene regulation networks can be kept in a
quasi-deterministic ``clockwork mode'' in spite of the absence of a central
clock. A delay similar to the one occurring in gene regulation causes
synchronization in the model. Stability under asynchronous dynamics disfavors
topologies containing loops, comparing well with the observed strong
suppression of loops in biological regulatory networks.


  In the present work, 16 genetic code doublets and their cognate amino acids
in the genetic code are fitted into a polyhedron model. Based on the structural
regularity in nucleobases, and by using a series of common-sense topological
approaches to rearranging the Hamiltonian-type graph of the codon map, it is
identified that the degeneracy of codons and the internal relation of the 20
amino acids within the genetic code are in agreement with the spherical and
polyhedral symmetry of a quasi-28-gon, i.e., icosikaioctagon. Hence, a
quasi-central, quasi-polyhedral and rotational symmetry within the genetic code
is described. Accordingly, the rotational symmetry of the numerical
distribution of side-chain carbon atoms of the 20 amino acids and the
side-chain skeleton atoms (carbon, nitrogen, oxygen and sulfur) of the 20 amino
acids are presented in the framework of this quasi-28-gon model. Two
evolutionary axes within the 20 standard amino acids are suggested.


  A population evolving in an inhomogeneous environment will adapt differently
to different regions. We study the conditions under which such a population can
maintain adaptations to a particular region when that region is not stationary,
but can move. In particular, we study a quasispecies living near a favorable
patch ("oasis") in the middle of a large "desert." The population has two
genetic states, one of which which conveys a relative advantage while in the
oasis at the cost of a disadvantage in the desert. We consider the population
dynamics when the oasis is moving, or equivalently some form of "wind" is
blowing the population away from the oasis. We find that the ratio of the two
types of individuals exhibits sharp transitions at particular oasis velocities.
We calculate an extinction velocity, and a switching velocity above which the
dominance switches from the oasis-adapted genotype to the desert-adapted one.
This switching velocity is analagous to the quasispecies mutational error
threshold. Above this velocity, the population cannot maintain adaptations to
the properties of the oasis.


  Identifying the driving forces and the mechanism of association of
huntingtin-exon1, a close marker for the progress of Huntington's disease, is
an important prerequisite towards finding potential drug targets, and
ultimately a cure. We introduce here a modelling framework based on a key
analogy of the physico-chemical properties of the exon1 fragment to block
copolymers. We use a systematic mesoscale methodology, based on Dissipative
Particle Dynamics, which is capable of overcoming kinetic barriers, thus
capturing the dynamics of significantly larger systems over longer times than
considered before. Our results reveal that the relative hydrophobicity of the
poly-glutamine block as compared to the rest of the (proline-based) exon1
fragment, ignored to date, constitutes a major factor in the initiation of the
self-assembly process. We find that the assembly is governed by both the
concentration of exon1 and the length of the poly-glutamine stretch, with a low
length threshold for association even at the lowest volume fractions we
considered. Moreover, this self-association occurs irrespective of whether the
glutamine stretch is in random coil or hairpin configuration, leading to
spherical or cylindrical assemblies, respectively. We discuss the implications
of these results for reinterpretation of existing research within this context,
including that the routes towards aggregation of exon1 may be distinct to those
of the widely studied homopolymeric poly-glutamine peptides.


  We study counter expressed gene networks constructed from gene-expression
data obtained from many types of cancers. The networks are synthesized by
connecting vertices belonging to each others' list of K-farthest-neighbors,
with K being an a priori selected non-negative integer. In the range of K
corresponding to minimum homogeneity, the degree distribution of the networks
shows scaling. Clustering in these networks is smaller than that in equivalent
random graphs and remains zero till significantly large K. Their small
diameter, however, implies small-world behavior which is corroborated by their
eigenspectrum. We discuss implications of these findings in several contexts.


  We identify a set of 575 human genes that are expressed in all conditions
tested in a publicly available database of microarray results. Based on this
common occurrence, the set is expected to be rich in "housekeeping" genes,
showing constitutive expression in all tissues. We compare selected aspects of
their genomic structure with a set of background genes. We find that the
introns, untranslated regions and coding sequences of the housekeeping genes
are shorter, indicating a selection for compactness in these genes.


  The defense response in plants challenged with pathogens is characterized by
the activation of a diverse set of genes. Many of the same genes are induced in
the defense responses of a wide range of plant species. How plant defense gene
families evolve may therefore provide an important clue to our understanding of
how disease resistance evolves. Because studies usually focus on a single host
species, little data are available regarding changes in defense gene expression
patterns as species diverge. The expression of defense-induced genes PR10,
chitinase and chalcone synthase was assayed in four pea species (Pisum sativum,
P. humile, P. elatius and P. fulvum) and two Lathyrus species (L. sativus and
L. tingitanus) which exhibited a range of infection phenotypes with Fusarium
solani . In P. sativum, resistance was accompanied by a strong induction of
defense genes at 8 hr. post-inoculation. Weaker induction was seen in
susceptible interactions in wild species. Divergence in the timing of PR10
expression was most striking between P. sativum and its closest realtive, P.
humile. Two members of this multigene family, designated PR10.1 and PR10.2, are
strongly-expressed in response to Fusarium, while the PR10.3 gene is more
weakly expressed, among Pisum species. The rapidity with which PR10 expression
evolves raises the question, is divergence of defense gene expression a part of
the phenotypic diversity underlying plant/pathogen coevolution?


  Hierarchical structure is an essential part of complexity, important notion
relevant for a wide range of applications ranging from biological population
dynamics through robotics to social sciences. In this paper we propose a simple
cellular-automata tool for study of hierarchical population dynamics.


  We investigate the process of fixation of advantageous mutations in an
asexual population. We assume that the effect of each beneficial mutation is
exponentially distributed with mean value $\omega_{med}=1/\beta$. The model
also considers that the effect of each new deleterious mutation reduces the
fitness of the organism independent on the previous number of mutations. We use
the branching process formulation and also extensive simulations to study the
model. The agreement between the analytical predictions and the simulational
data is quite satisfactory. Surprisingly, we observe that the dependence of the
probability of fixation $P_{fix}$ on the parameter $\omega_{med}$ is precisely
described by a power-law relation, $P_{fix} \sim \omega_{med}^{\gamma}$. The
exponent $\gamma$ is an increase function of the rate of deleterious mutations
$U$, whereas the probability $P_{fix}$ is a decreasing function of $U$. The
mean value $\omega_{fix}$ of the beneficial mutations which reach ultimate
fixation depends on $U$ and $\omega_{med}$. The ratio
$\omega_{fix}/\omega_{med}$ increases as we consider higher values of mutation
value $U$ in the region of intermediate to large values of $\omega_{med}$,
whereas for low $\omega_{med}$ we observe the opposite behavior.


  Many biological processes involve one dimensional diffusion over a correlated
inhomogeneous energy landscape with a correlation length $\xi_c$. Typical
examples are specific protein target location on DNA, nucleosome repositioning,
or DNA translocation through a nanopore, in all cases with $\xi_c\approx$ 10
nm. We investigate such transport processes by the mean first passage time
(MFPT) formalism, and find diffusion times which exhibit strong sample to
sample fluctuations. For a a displacement $N$, the average MFPT is diffusive,
while its standard deviation over the ensemble of energy profiles scales as
$N^{3/2}$ with a large prefactor. Fluctuations are thus dominant for
displacements smaller than a characteristic $N_c \gg \xi_c$: typical values are
much less than the mean, and governed by an anomalous diffusion rule. Potential
biological consequences of such random walks, composed of rapid scans in the
vicinity of favorable energy valleys and occasional jumps to further valleys,
is discussed.


  We expand the current understanding of ecological resilience to include the
nested hierarchy of cognitive submodules that particularly, if not uniquely,
delineates human ecosystems. These modules, ranging from the immune system to
the local social network, are embedded in a cultural milieu which, to take the
perspective of the evolutionary anthropologist Robert Boyd, ''is as much a part
of human biology as the enamel on our teeth''. We begin by extending recent
treatments of cognitive process as associated with characteristic information
sources to a certain class of ecosystems through a generalization of
coarse-graining. In the spirit of the Large Deviations Program, we then import
renormalization formalism via the Asymptotic Equipartition Theorem to obtain
punctuated response to parameters of increasing habitat degradation. A Legendre
transform of an appropriate joint information permits analysis away from
critical points, and generates the expected quasi-stability in a highly natural
manner. We interweave the discussion with applications to the public health
impacts of the massive deurbanization and deindustrialization presently
afflicting the United States.


  The interactions between proteins, DNA, and RNA in living cells constitute
molecular networks that govern various cellular functions. To investigate the
global dynamical properties and stabilities of such networks, we studied the
cell-cycle regulatory network of the budding yeast. With the use of a simple
dynamical model, it was demonstrated that the cell-cycle network is extremely
stable and robust for its function. The biological stationary state--the G1
state--is a global attractor of the dynamics. The biological pathway--the
cell-cycle sequence of protein states--is a globally attracting trajectory of
the dynamics. These properties are largely preserved with respect to small
perturbations to the network. These results suggest that cellular regulatory
networks are robustly designed for their functions.


  Many important cellular functions are implemented by protein complexes that
act as sophisticated molecular machines of varying size and temporal stability.
Here we demonstrate quantitatively that protein complexes in the yeast,
Saccharomyces cerevisiae, are comprised of a core in which subunits are highly
coexpressed, display the same deletion phenotype (essential or non-essential)
and share identical functional classification and cellular localization. This
core is surrounded by a functionally mixed group of proteins, which likely
represent short-lived- or spurious attachments. The results allow us to define
the deletion phenotype and cellular task of most known complexes, and to
identify with high confidence the biochemical role of hundreds of proteins with
yet unassigned functionality.


  The molecular mechanism of the solvent motion that is required to instigate
the protein structural relaxation above a critical hydration level or
transition temperature has yet to be determined. In this work we use
quasi-elastic neutron scattering (QENS) and molecular dynamics simulation to
investigate hydration water dynamics near a greatly simplified protein surface.
We consider the hydration water dynamics near the completely deuterated
N-acetyl-leucine-methylamide (NALMA) solute, a hydrophobic amino acid side
chain attached to a polar blocked polypeptide backbone, as a function of
concentration between 0.5M-2.0M, under ambient conditions. In this
Communication, we focus our results of hydration dynamics near a model protein
surface on the issue of how enzymatic activity is restored once a critical
hydration level is reached, and provide a hypothesis for the molecular
mechanism of the solvent motion that is required to trigger protein structural
relaxation when above the hydration transition.


  We study a four species ecological system with cyclic dominance whose
individuals are distributed on a square lattice. Randomly chosen individuals
migrate to one of the neighboring sites if it is empty or invade this site if
occupied by their prey. The cyclic dominance maintains the coexistence of all
the four species if the concentration of vacant sites is lower than a threshold
value. Above the treshold, a symmetry breaking ordering occurs via growing
domains containing only two neutral species inside. These two neutral species
can protect each other from the external invaders (predators) and extend their
common territory. According to our Monte Carlo simulations the observed phase
transition is equivalent to those found in spreading models with two equivalent
absorbing states although the present model has continuous sets of absorbing
states with different portions of the two neutral species. The selection
mechanism yielding symmetric phases is related to the domain growth process
whith wide boundaries where the four species coexist.


  We study a recent model for calcium signal transduction. This model displays
spiking, bursting and chaotic oscillations in accordance with experimental
results. We calculate bifurcation diagrams and study the bursting behaviour in
detail. This behaviour is classified according to the dynamics of separated
slow and fast subsystems. It is shown to be of the Fold-Hopf type, a type which
was previously only described in the context of neuronal systems, but not in
the context of signal transduction in the cell.


  The immune system is a real-time example of an evolving system that navigates
the essentially infinite complexity of protein sequence space. How this system
responds to disease and vaccination is discussed. Of particular focus is the
case when vaccination leads to increased susceptibility to disease, a
phenomenon termed original antigenic sin. A physical theory of protein
evolution to explain limitations in the immune system response to vaccination
and disease is discussed, and original antigenic sin is explained as stemming
from localization of the immune system response in antibody sequence space.
This localization is a result of the roughness in sequence space of the evolved
antibody affinity constant for antigen and is observed for diseases with high
year-to-year mutation rates, such as influenza.


  We analyze the dependence of thermal denaturation transition and folding
rates of globular proteins on the number of amino acid residues, N. Using
lattice Go models we show that DeltaT/T_F ~ N^-1, where T_F is the folding
transition temperature and DeltaT is the folding transition width. This finding
is consistent with finite size effects expected for the systems undergoing a
phase transition from a disordered to an ordered phase. The dependence of the
folding rates k_F on N for lattice models and the dataset of 57 proteins and
peptides shows that k_F = k_F^0 exp(-CN^beta) provides a good fit, if 0 < beta
<= 2/3 and C is a constant. We find that k_F = k_F^0 exp(-1.1N^0.5) with k_F^0
=(0.4x10^-6 s)^-1 can estimate optimal protein folding rates to within an order
of magnitude in most cases. By using this fit for a set of proteins with
beta-sheet topology we find that k_F^0 is approximately equal to k_U^0, the
prefactor for unfolding rates. The maximum ratio of k_U^0/k_F^0 is 10 for this
class of proteins.


  Former work on an application of order-disorder theory is recalled as a
vehicle to add further development and significance to the recent paper on
motifs in protein interactions.


  We propose a new kernel for biological sequences which borrows ideas and
techniques from information theory and data compression. This kernel can be
used in combination with any kernel method, in particular Support Vector
Machines for protein classification. By incorporating prior biological
assumptions on the properties of amino-acid sequences and using a Bayesian
averaging framework, we compute the value of this kernel in linear time and
space, benefiting from previous achievements proposed in the field of universal
coding. Encouraging classification results are reported on a standard protein
homology detection experiment.


  The asymmetry in the shapes of folded and unfolded states are probed using
two parameters, one being a measure of the sphericity and the other that
describes the shape. For the folded states, whose interiors are densely packed,
the radii of gyration (Rg) and these two parameters are calculated using the
coordinates of the experimentally determined structures. Although Rg scales as
expected for maximally compact structures, the distributions of the shape
parameters show that there is considerable asymmetry in the shapes of folded
structures. The degree of asymmetry is greater for proteins that form
oligomers. Analysis of the two- and three-body contacts in the native
structures shows that the presence of near equal number of contacts between
backbone and side-chains and between side-chains gives rise to dense packing.
We suggest that proteins with relatively large values of shape parameters can
tolerate volume mutations without greatly affecting the network of contacts or
their stability. To probe shape characteristics of denatured states we have
developed a model of a WW-like domain. The shape parameters, which are
calculated using Langevin simulations, change dramatically in the course of
coil to globule transition. Comparison of the values of shape parameters
between the globular state and the folded state of WW domain shows that both
energetic (especially dispersion in the hydrophobic interactions) and steric
effects are important in determining packing in proteins.


  Understanding why some cellular components are conserved across species,
while others evolve rapidly is a key question of modern biology. Here we
demonstrate that in S. cerevisiae proteins organized in cohesive patterns of
interactions are conserved to a significantly higher degree than those that do
not participate in such motifs. We find that the conservation of proteins
within distinct topological motifs correlates with the motif's
inter-connectedness and function and also depends on the structure of the
overall interactome topology. These findings indicate that motifs may represent
evolutionary conserved topological units of cellular networks molded in
accordance with the specific biological function in which they participate.


  Instead of conformation states of single residues, refined conformation
states of quintuplets are proposed to reflect conformation correlation. Simple
hidden Markov models combining with sliding window scores are used for
predicting secondary structure of a protein from its amino acid sequence. Since
the length of protein conformation segments varies in a narrow range, we ignore
the duration effect of the length distribution. The window scores for residues
are a window version of the Chou-Fasman propensities estimated under an
approximation of conditional independency. Different window widths are
examined, and the optimal width is found to be 17. A high accuracy about 70% is
achieved.


  We introduce a random bit-string model of post-transcriptional genetic
regulation based on sequence matching. The model spontaneously yields a scale
free network with power law scaling with $ \gamma=-1$ and also exhibits
log-periodic behaviour. The in-degree distribution is much narrower, and
exhibits a pronounced peak followed by a Gaussian distribution. The network is
of the smallest world type, with the average minimum path length independent of
the size of the network, as long as the network consists of one giant cluster.
The percolation threshold depends on the system size.


  Evolving biomolecular networks have to combine the stability against
perturbations with flexibility allowing their constituents to assume new roles
in the cell. Gene duplication followed by functional divergence of associated
proteins is a major force shaping molecular networks in living organisms.
Recent availability of system-wide data for yeast S. Cerevisiae allow us to
access the effects of gene duplication on robustness and plasticity of
molecular networks. We demonstrate that the upstream transcriptional regulation
of duplicated genes diverges fast, losing on average 4% of their common
transcription factors for every 1% divergence of their amino acid sequences. In
contrast, the set of physical interaction partners of their protein products
changes much slower. The relative stability of downstream functions of
duplicated genes, is further corroborated by their ability to substitute for
each other in gene knockout experiments. We believe that the combination of the
upstream plasticity and the downstream robustness is a general feature
determining the evolvability of molecular networks.


  We have performed a statistical analysis of the spatial distribution of
operons in the transcriptional regulation network of Escherichia coli. The
analysis reveals that operons that regulate each other and operons that are
coregulated tend to lie next to each other on the genome. Moreover, these pairs
of operons tend to be transcribed in diverging directions. This spatial
arrangement of operons allows the upstream regulatory regions to interfere with
each other. This affords additional regulatory control, as illustrated by a
mean-field analysis of a feed-forward loop. Our results suggest that regulatory
control can provide a selection pressure that drives operons together in the
course of evolution.


  Animal mitochondrial genomes usually have two transfer RNAs for Leucine: one,
with anticodon UAG, translates the four-codon family CUN, whilst the other,
with anticodon UAA, translates the two-codon family UUR. These two genes must
differ at the third anticodon position, but in some species the genes differ at
many additional sites, indicating that these genes have been independent for a
long time. Duplication and deletion of genes in mitochondrial genomes occurs
frequently during the evolution of the Metazoa. If a tRNA-Leu gene were
duplicated and a substitution occurred in the anticodon, this would effectively
turn one type of tRNA into the other. The original copy of the second tRNA type
might then be lost by a deletion elsewhere in the genome. There are several
groups of species in which the two tRNA-Leu genes occur next to one another (or
very close) on the genome, which suggests that tandem duplication has occurred.
Here we use RNA-specific phylogenetic methods to determine evolutionary trees
for both genes. We present evidence that the process of duplication, anticodon
mutation and deletion of tRNA-Leu genes has occurred at least five times during
the evolution of the Metazoa - once in the common ancestor of all Protostomes,
once in the common ancestor of Echinoderms and Hemichordates, once in the
hermit crab, and twice independently in Molluscs.


  The transition from localized to systemic spreading of bacteria, viruses and
other agents is a fundamental problem that spans medicine, ecology, biology and
agriculture science. We have conducted experiments and simulations in a simple
one-dimensional system to determine the spreading of bacterial populations that
occurs for an inhomogeneous environment under the influence of external
convection. Our system consists of a long channel with growth inhibited by
uniform UV illumination except in a small ``oasis'', which is shielded from the
UV light. To mimic blood flow or other flow past a localized infection, the
oasis is moved with a constant velocity through the UV-illuminated ``desert''.
The experiments are modeled with a convective reaction-diffusion equation. In
both the experiment and model, localized or extinct populations are found to
develop, depending on conditions, from an initially localized population. The
model also yields states where the population grows everywhere. Further, the
model reveals that the transitions between localized, extended, and extinct
states are continuous and non-hysteretic. However, it does not capture the
oscillations of the localized population that are observed in the experiment.


  Background
  Nucleotide sequences contain multiple codes responsible for organism's
functioning and structure. They can be investigated by various signal
processing methods. These techniques are well suited for indication of
frequently encountered sequence motifs (i.e., repeats). However, if there are
two or more codes containing the same motif, the local nucleotide distribution
(i.e., profile), resulting from sequence alignment by the motif position, will
represent overlapping of the code patterns.
  Results
  The novel algorithm for decomposition of pattern overlapping is proposed. It
is capable to work with dispersed repeats as well. The algorithm is based on
cross-correlation procedure applied locally in a cumulative fashion. Its
sensitivity was tested on human genomic sequences.
  Conclusions
  Cumulative Local Cross-Correlation was successfully used to decompose
overlapping of nucleotide patterns in human genomic sequences. Being very
general technique (as general as cross-correlation), it can be easily adopted
in other signal processing applications and naturally extended for
multidimensional cases.
  Software implementation of the algorithm is available on request from the
authors.


  Is protein secondary structure primarily determined by local interactions
between residues closely spaced along the amino acid backbone, or by non-local
tertiary interactions? To answer this question we have measured the entropy
densities of primary structure and secondary structure sequences, and the local
inter-sequence mutual information density. We find that the important
inter-sequence interactions are short ranged, that correlations between
neighboring amino acids are essentially uninformative, and that only 1/4 of the
total information needed to determine the secondary structure is available from
local inter-sequence correlations. Since the remaining information must come
from non-local interactions, this observation supports the view that the
majority of most proteins fold via a cooperative process where secondary and
tertiary structure form concurrently. To provide a more direct comparison to
existing secondary structure prediction methods, we construct a simple hidden
Markov model (HMM) of the sequences. This HMM achieves a prediction accuracy
comparable to other single sequence secondary structure prediction algorithms,
and can extract almost all of the inter-sequence mutual information. This
suggests that these algorithms are almost optimal, and that we should not
expect a dramatic improvement in prediction accuracy. However, local
correlations between secondary and primary structure are probably of
under-appreciated importance in many tertiary structure prediction methods,
such as threading.


  Mortality is an instrument of natural selection. Evolutionary motivated
theories imply its irreversibility and life history dependence. This is
inconsistent with mortality data for protected populations. Accurate analysis
yields mortality law, which is specific for their evolutionary unprecedented
conditions, yet universal for species as evolutionary remote as humans and
flies. The law is exact, instantaneous, reversible, stepwise, and allows for a
rapid (within less than two years for humans) and significant mortality
decrease at any (but very old) age.


  We propose that the evolutionary purpose of junk DNA is to protect the gene.
Mutation agents, such as retro-viruses, hit ``decoy'' DNA most of the time.
Although the argument is far from general, we propose that the percentage of
junk DNA should correlate with the number of retroviruses attacking a given
species. It should also anti-correlate with the ideal mutation rates (higher in
insects than in mammals).


  An analysis of the Kimura 3ST model of DNA sequence evolution is given on the
basis of its continuous Lie symmetries. The rate matrix commutes with a
U(1)xU(1)xU(1) phase subgroup of the group GL(4) of 4x4x4 invertible complex
matrices acting on a linear space spanned by the 4 nucleic acid base letters.
The diagonal `branching operator' representing speciation is defined, and shown
to intertwine the U(1)xU(1)xU(1) action. Using the intertwining property, a
general formula for the probability density on the leaves of a binary tree
under the Kimura model is derived, which is shown to be equivalent to
established phylogenetic spectral transform methods.


  While it is not possible to directly the observe evolution of multigene
families, the best alternative is to compare orthologous family members among
several closely-related species with varying degrees of reproductive isolation.
Using RT-PCR we show that in pea (Pisum sativum) each member of the
pathogenesis-related PR10 family has a distinct pattern of expression in
response to the fungus Fusarium solani, and in treatment with salicylic acid,
chitosan and abcisic acid. Sequencing reveals that PR10.1, PR10.2 and PR10.3
exist in P. humile, P. elatius and P. fulvum, except that no PR10.2 orthologue
was identified in P. elatius. PR10.1, PR10.2 and PR10.3 appear to have diverged
from a single gene in the common Pisum ancestor. For the recently diverged
PR10.1 and PR10.2, the timing of fungal-induced expression differs greatly
among species. For example, PR10.1 was strongly induced in P. sativum by F.
solani within 8 hours postinoculation (h.p.i.), whereas little PR10.1
expression was seen in pea's closest relative, P. humile, and in the more
distantly-related P. elatius. In P. fulvum, expression did not peak until 48
h.p.i. Expression of the more ancient PR10.4 and PR10.5 genes is more tightly
conserved among Pisum species. These data indicate that expression, as well as
sequence, can evolve rapidly. We hypothesize that changes in differential
expression of multigene family members could provide a source of phenotypic
diversity in populations, which may be of particular importance to
plant/pathogen coevolution.


  BACKGROUND: Transcriptional regulation is a key mechanism in the functioning
of the cell, and is mostly effected through transcription factors binding to
specific recognition motifs located upstream of the coding region of the
regulated gene. The computational identification of such motifs is made easier
by the fact that they often appear several times in the upstream region of the
regulated genes, so that the number of occurrences of relevant motifs is often
significantly larger than expected by pure chance. RESULTS: To exploit this
fact, we construct sets of genes characterized by the statistical
overrepresentation of a certain motif in their upstream regions. Then we study
the functional characterization of these sets by analyzing their annotation to
Gene Ontology terms. For the sets showing a statistically significant specific
functional characterization, we conjecture that the upstream motif
characterizing the set is a binding site for a transcription factor involved in
the regulation of the genes in the set. CONCLUSIONS: The method we propose is
able to identify many known binding sites in S. cerevisiae and new candidate
targets of regulation by known transcription factors. Its application to less
well studied organisms is likely to be valuable in the exploration of their
regulatory interaction network.


  We study genetic switches formed from pairs of mutually repressing operons.
The switch stability is characterised by a well defined lifetime which grows
sub-exponentially with the number of copies of the most-expressed transcription
factor, in the regime accessible by our numerical simulations. The stability
can be markedly enhanced by a suitable choice of overlap between the upstream
regulatory domains. Our results suggest that robustness against biochemical
noise can provide a selection pressure that drives operons, that regulate each
other, together in the course of evolution.


  We propose a generic model of eco-systems, with a {\it hierarchical} food web
structure. In our computer simulations we let the eco-system evolve
continuously for so long that that we can monitor extinctions as well as
speciations over geological time scales. {\it Speciation} leads not only to
horizontal diversification of species at any given trophic level but also to
vertical bio-diversity that accounts for the emergence of complex species from
simpler forms of life. We find that five or six trophic levels appear as the
eco-system evolves for sufficiently long time, starting initially from just one
single level. Moreover, the time intervals between the successive collections
of ecological data is so short that we could also study ``micro''-evolution of
the eco-system, i.e., the birth, ageing and death of individual organisms.


  The recently proposed compressed backbone theory suggested that the intrinsic
curvature in DNA can result from a geometric mismatch between the specific
backbone length and optimal base stacking orientations. It predicted that the
curvature in A-tract repeats can be relaxed by introducing single stranded
breaks (nicks). This effect has not been tested earlier and it would not be
accounted for by alternative models of DNA bending. Here the curvature in a
specifically designed series of nicked DNA fragments is tested experimentally
by gel mobility assays and, simultaneously, by free molecular dynamics
simulations. Single stranded breaks produce virtually no effect upon the gel
mobility of the random sequence DNA. In contrast, nicked A-tract fragments
reveal a regular modulation of curvature depending upon the position of the
strand break with respect to the overall bend. Maximal relaxation is observed
when nicks occur inside A-tracts. The results are partially reproduced in
simulations. Analysis of computed curved DNA conformations reveals a group of
sugar atoms that exhibit reduced backbone length within A-tracts, which can
correspond to the compression hypothesis.


  We present a solvable model that predicts the folding kinetics of two-state
proteins from their native structures. The model is based on conditional chain
entropies. It assumes that folding processes are dominated by small-loop
closure events that can be inferred from native structures. For CI2, the src
SH3 domain, TNfn3, and protein L, the model reproduces two-state kinetics, and
it predicts well the average Phi-values for secondary structures. The barrier
to folding is the formation of predominantly local structures such as helices
and hairpins, which are needed to bring nonlocal pairs of amino acids into
contact.


  With a generalized language-of-thought argument for immune cognition, we
model how population-directed structured psychosocial stress can impose an
image of itself on the coevolutionary conflict between a highly adaptive
chronic infection and the immune system. We raise the possibility that, for
successful 'evolution machines' like HIV and malaria, simplistic
individual-oriented magic bullet drug treatments, vaccines, behavior
modification, or other interventions that do not address the fundamental living
and working conditions underlying disease ecology will fail to control current
epidemics and may select for holistic pathogenic life histories which increase
disease virulence.


  This paper summarizes and re-evaluates Prigogine's evolution of thought, from
a more classical view of thermodynamics with negative implications for the
evolution and persistence of life, through to a far more general and open
formulation of thermodynamics, which does not require assumption of a Big Bang.
The paper also proposes an encoding scheme for the statistics of Lorentzian
systems into mixed Fermi-Bose density matrices, to complete the program given
in quant-ph/0309087, which forcibly points towards symmetry in time at the
microscopic level. The paper suggests future directions to reconcile
microscopic time symmetry and macroscopic asymmetry and evolution, and
concludes that the space-time continuum may be larger and more interesting than
we are able yet to pin down.


  The determination of the folding mechanisms of proteins is critical to
understand the topological change that can propagate Alzheimer and
Creutzfeld-Jakobs diseases, among others. The computational community has paid
considerable attention to this problem; however, the associated time scale,
typically on the order of milliseconds or more, represents a formidable
challenge. Ab initio protein folding from long molecular dynamics (MD)
simulations or ensemble dynamics is not feasible with ordinary computing
facilities and new techniques must be introduced. Here we present a detailed
study of the folding of a 16-residue beta-hairpin, described by a generic
energy model and using the activation-relaxation technique. From a total of 90
trajectories at 300 K, three folding pathways emerge. All involve a
simultaneous optimization of the complete hydrophobic and hydrogen bonding
interactions. The first two follow closely those observed by previous
theoretical studies. The third pathway, never observed by previous all-atom
folding, unfolding and equilibrium simulations, can be described as a reptation
move of one strand of the beta-sheet with respect to the other. This reptation
move indicates that non-native interactions can play a dominant role in the
folding of secondary structures. These results point to a more complex folding
picture than expected for a simple beta-hairpin.


  This paper presents a unified mathematical framework for inference in
graphical models, building on the observation that graphical models are
algebraic varieties.
  From this geometric viewpoint, observations generated from a model are
coordinates of a point in the variety, and the sum-product algorithm is an
efficient tool for evaluating specific coordinates. The question addressed here
is how the solutions to various inference problems depend on the model
parameters. The proposed answer is expressed in terms of tropical algebraic
geometry. A key role is played by the Newton polytope of a statistical model.
Our results are applied to the hidden Markov model and to the general Markov
model on a binary tree.


  Analytic estimates for the forces and free energy generated by bilayer
deformation reveal a compelling and intuitive model for MscL channel gating
analogous to the nucleation of a second phase. We argue that the competition
between hydrophobic mismatch and tension results in a surprisingly rich story
which can provide both a quantitative comparison to measurements of opening
tension for MscL when reconstituted in bilayers of different thickness and
qualitative insights into the function of the MscL channel and other
transmembrane proteins.


  It is important to understand how protein folding and evolution influences
each other. Several studies based on entropy calculation correlating
experimental measurement of residue participation in folding nucleus and
sequence conservation have reached different conclusions. Here we report
analysis of conservation of folding nucleus using an evolutionary model
alternative to entropy based approaches. We employ a continuous time Markov
model of codon substitution to distinguish mutation fixed by evolution and
mutation fixed by chance. This model takes into account bias in codon
frequency, bias favoring transition over transversion, as well as explicit
phylogenetic information. We measure selection pressure using the ratio
$\omega$ of synonymous vs. non-synonymous substitution at individual residue
site. The $\omega$-values are estimated using the {\sc Paml} method, a
maximum-likelihood estimator. Our results show that there is little correlation
between the extent of kinetic participation in protein folding nucleus as
measured by experimental $\phi$-value and selection pressure as measured by
$\omega$-value. In addition, two randomization tests failed to show that
folding nucleus residues are significantly more conserved than the whole
protein. These results suggest that at the level of codon substitution, there
is no indication that folding nucleus residues are significantly more conserved
than other residues. We further reconstruct candidate ancestral residues of the
folding nucleus and suggest possible test tube mutation studies of ancient
folding nucleus.


  Biological phenomena differ significantly from physical phenomena. At the
heart of this distinction is the fact that biological entities have
computational abilities and thus they are inherently difficult to predict. This
is the reason why simplified models that provide the minimal requirements for
computation turn out to be very useful to study networks of many components. In
this chapter, we briefly review the dynamical aspects of models of regulatory
networks, discussing their most salient features, and we also show how these
models can give clues about the way in which networks may organize their
capacity to evolve, by providing simple examples of the implementation of
robustness and modularity.


  The emergence of complex patterns of organization close to the Cambrian
boundary is known to have happened over a (geologically) short period of time.
It involved the rapid diversification of body plans and stands as one of the
major transitions in evolution. How it took place is a controversial issue.
Here we explore this problem by considering a simple model of pattern formation
in multicellular organisms. By modeling gene network-based morphogenesis and
its evolution through adaptive walks, we explore the question of how
combinatorial explosions might have been actually involved in the Cambrian
event. Here we show that a small amount of genetic complexity including both
gene regulation and cell-cell signaling allows one to generate an extraordinary
repertoire of stable spatial patterns of gene expression compatible with
observed anteroposterior patterns in early development of metazoans. The
consequences for the understanding of the tempo and mode of the Cambrian event
are outlined.


  A new version of DNA walks, where nucleotides are regarded unequal in their
contribution to a walk is introduced, which allows us to study thoroughly the
"fine structure" of nucleotide sequences. The approach is based on the
assumption that nucleotides have an inner abstract characteristics, the
determinative degree, which reflects phenomenological properties of genetic
code and is adjusted to nucleotides physical properties. We consider each
position in codon independently, which gives three separate walks being
characterized by different angles and lengths, and such an object is called
triander which reflects the "strength" of branch. A general method of
identification of DNA sequence "by triander", which can be treated as a unique
"genogram", "gene passport" is proposed. The two- and three-dimensional
trianders are considered. The difference of sequence fine structure in genes
and the intergenic space is shown. A clear triplet signal in coding locuses is
found which is absent in the intergenic space and is independent from the
sequence length. The topological classification of trianders is presented which
can allow us to provide a detail working out signatures of functionally
different genomic regions.


  Mechanical stretching of six proteins is studied through molecular dynamics
simulations. The model is Go-like, with Lennard-Jones interactions at native
contacts. Low temperature unfolding scenarios are remarkably complex and
sensitive to small structural changes. Thermal fluctuations reduce the peak
forces and the number of metastable states during unfolding. The unfolding
pathways also simplify as temperature rises. In the entropic limit, all
proteins show a monotonic decrease of the extension where bonds rupture with
their separation along the backbone (contact order).


  The use of gene microchips has enabled a rapid accumulation of
gene-expression data. One of the major challenges of analyzing this data is the
diversity, in both size and signal strength, of the various modules in the gene
regulatory networks of organisms. Based on the Iterative Signature Algorithm
[Bergmann, S., Ihmels, J. and Barkai, N. (2002) Phys. Rev. E 67, 031902], we
present an algorithm - the Progressive Iterative Signature Algorithm (PISA) -
that, by sequentially eliminating modules, allows unsupervised identification
of both large and small regulatory modules. We applied PISA to a large set of
yeast gene-expression data, and, using the Gene Ontology annotation database as
a reference, found that our algorithm is much better able to identify
regulatory modules than methods based on high-throughput transcription-factor
binding experiments or on comparative genomics.


  We describe a new global multiple alignment program capable of aligning a
large number of genomic regions. Our progressive alignment approach
incorporates the following ideas: maximum-likelihood inference of ancestral
sequences, automatic guide-tree construction, protein based anchoring of
ab-initio gene predictions, and constraints derived from a global homology map
of the sequences. We have implemented these ideas in the MAVID program, which
is able to accurately align multiple genomic regions up to megabases long.
MAVID is able to effectively align divergent sequences, as well as incomplete
unfinished sequences. We demonstrate the capabilities of the program on the
benchmark CFTR region which consists of 1.8Mb of human sequence and 20
orthologous regions in marsupials, birds, fish, and mammals. Finally, we
describe two large MAVID alignments: an alignment of all the available HIV
genomes and a multiple alignment of the entire human, mouse and rat genomes.


  Expression of the Drosophila segment polarity genes is initiated by a
prepattern of pair-rule gene products and maintained by a network of regulatory
interactions throughout several stages of embryonic development. Analysis of a
model of gene interactions based on differential equations showed that
wild-type expression patterns of these genes can be obtained for a wide range
of kinetic parameters, which suggests that the steady states are determined by
the topology of the network and the type of regulatory interactions between
components, not the detailed form of the rate laws. To investigate this, we
propose and analyze a Boolean model of this network which is based on a binary
ON/OFF representation of transcription and protein levels, and in which the
interactions are formulated as logical functions. In this model the spatial and
temporal patterns of gene expression are determined by the topology of the
network and whether components are present or absent, rather than the absolute
levels of the mRNAs and proteins and the functional details of their
interactions. The model is able to reproduce the wild type gene expression
patterns, as well as the ectopic expression patterns observed in
over-expression experiments and various mutants. Furthermore, we compute
explicitly all steady states of the network and identify the basin of
attraction of each steady state. The model gives important insights into the
functioning of the segment polarity gene network, such as the crucial role of
the wingless and sloppy paired genes, and the network's ability to correct
errors in the prepattern.


  Recent work on mutation-selection models has revealed that, under specific
assumptions on the fitness function and the mutation rates, asymptotic
estimates for the leading eigenvalue of the mutation-reproduction matrix may be
obtained through a low-dimensional maximum principle in the limit N to infinity
(where N is the number of types). In order to extend this variational principle
to a larger class of models, we consider here a family of reversible N by N
matrices and identify conditions under which the high-dimensional Rayleigh-Ritz
variational problem may be reduced to a low-dimensional one that yields the
leading eigenvalue up to an error term of order 1/N. For a large class of
mutation-selection models, this implies estimates for the mean fitness, as well
as a concentration result for the ancestral distribution of types.


  Increased biological complexity is generally associated with the addition of
new genetic information, which must be integrated into the existing regulatory
network that operates within the cell. General arguments on network control, as
well as several recent genomic observations, indicate that regulatory gene
number grows disproportionally fast with increasing genome size. We present two
models for the growth of regulatory networks. Both predict that the number of
transcriptional regulators will scale quadratically with total gene number.
This appears to be in good quantitative agreement with genomic data from 89
fully sequenced prokaryotes. Moreover, the empirical curve predicts that any
new non-regulatory gene will be accompanied by more than one additional
regulator beyond a genome size of about 20,000 genes, within a factor of two of
the observed ceiling. Our analysis places transcriptional regulatory networks
in the class of accelerating networks. We suggest that prokaryotic complexity
may have been limited throughout evolution by regulatory overhead, and
conversely that complex eukaryotes must have bypassed this constraint by novel
strategies.


  Many signalling functions in molecular biology require proteins bind to
substrates such as DNA in response to environmental signals such as the
simultaneous binding to a small molecule. Examples are repressor proteins which
may transmit information via a conformational change in response to the ligand
binding. An alternative entropic mechanism of ``allostery'' suggests that the
inducer ligand changes the intramolecular vibrational entropy not just the
static structure. We present a quantitative, coarse-grained model of entropic
allostery that suggests design rules for internal cohesive potentials in
proteins employing this effect. It also addresses the issue of how the signal
information to bind or unbind is transmitted through the protein. The model may
be applicable to a wide range of repressors and also to signalling in
transmembrane proteins.


  Due to the conventional distinction between ecological (rapid) and
evolutionary (slow)timescales, ecological and population models to date have
typically ignored the effects of evolution. Yet the potential for rapid
evolutionary change has been recently established and may be critical to
understanding how populations adapt to changing environments. In this paper we
examine the relationship between ecological and evolutionary dynamics, focusing
on a well-studied experimental aquatic predator-prey system (Fussmann et al.
2000; Shertzer et al. 2002; Yoshida et al. 2003). Major properties of
predator-prey cycles in this system are determined by ongoing evolutionary
dynamics in the prey population. Under some conditions, however, the
populations tend to apparently stable steady-state densities. These are the
subject of the present paper. We examine a previously developed model for the
system, to determine how evolution shapes properties of the equilibria, in
particular the number and identity of coexisting prey genotypes. We then apply
these results to explore how evolutionary dynamics can shape the responses of
the system to "management": externally imposed alterations in conditions.
Specifically, we compare the behavior of the system including evolutionary
dynamics, with predictions that would be made if the potential for rapid
evolutionary change is negelected. Finally, we posit some simple experiments to
verify our prediction that evolution can have significant qualitative effects
on observed population-level responses to changing conditions.


  Examining elevated rates of systemic lupus erythematosus in African-American
women from perspectives of immune cognition suggests the disease constitutes an
internalized physiological image of external structured psychosocial stress, a
'pathogenic social hierarchy' involving the synergism of racism and gender
discrimination in the context of policy-driven social disintegration which has
particularly affected ethnic minorities in the USA. The disorder represents the
punctuated resetting of normal immune self-image to a self-attacking excited
state, a process formally analogous to models of punctuated equilibrium in
evolutionary theory. Both onset and progression of disease may be stratified by
a relation to cyclic physiological responses which are long in comparison with
heartbeat period: circadian, hormonal, and annual light/termperature cycles.


  We consider a model of large regulatory gene expression networks where the
thresholds activating the sigmoidal interactions between genes and the signs of
these interactions are shuffled randomly. Such an approach allows for a
qualitative understanding of network dynamics in a lack of empirical data
concerning the large genomes of living organisms. Local dynamics of network
nodes exhibits the multistationarity and oscillations and depends crucially
upon the global topology of a "maximal" graph (comprising of all possible
interactions between genes in the network). The long time behavior observed in
the network defined on the homogeneous "maximal" graphs is featured by the
fraction of positive interactions ($0\leq \eta\leq 1$) allowed between genes.
There exists a critical value $\eta_c<1$ such that if $\eta<\eta_c$, the
oscillations persist in the system, otherwise, when $\eta>\eta_c,$ it tends to
a fixed point (which position in the phase space is determined by the initial
conditions and the certain layout of switching parameters). In networks defined
on the inhomogeneous directed graphs depleted in cycles, no oscillations arise
in the system even if the negative interactions in between genes present
therein in abundance ($\eta_c=0$). For such networks, the bidirectional edges
(if occur) influence on the dynamics essentially. In particular, if a number of
edges in the "maximal" graph is bidirectional, oscillations can arise and
persist in the system at any low rate of negative interactions between genes
($\eta_c=1$). Local dynamics observed in the inhomogeneous scalable regulatory
networks is less sensitive to the choice of initial conditions. The scale free
networks demonstrate their high error tolerance.


  In Xenopus early embryos, replication origins neither require specific DNA
sequences nor is there an efficient S/M checkpoint, even though the whole
genome (3 billion bases) is completely duplicated within 10-20 minutes. This
leads to the"random-completion problem" of DNA replication in embryos, where
one needs to find a mechanism that ensures complete, faithful, timely
reproduction of the genome without any sequence dependence of replication
origins. We analyze recent DNA replication data in Xenopus laevis egg extracts
and find discrepancies with models where replication origins are distributed
independently of chromatin structure. Motivated by these discrepancies, we have
investigated the role that chromatin looping may play in DNA replication. We
find that the loop-size distribution predicted from a wormlike-chain model of
chromatin can account for the spatial distribution of replication origins in
this system quantitatively. Together with earlier findings of increasing
frequency of origin firings, our results can explain the random-completion
problem. The agreement between experimental data (molecular combing) and
theoretical predictions suggests that the intrinsic stiffness of chromatin
loops plays a fundamental biological role in DNA replication in early-embryo
Xenopus in regulating the origin spacing.


  How DNA repair enzymes find the relatively rare sites of damage is not known
in great detail. Recent experiments and molecular data suggest that the
individual repair enzymes do not work independently of each other, but rather
interact with each other through currents exchanged along DNA. A damaged site
in DNA hinders this exchange and this makes it possible to quickly free up
resources from error free stretches of DNA. Here the size of the speedup gained
from this current exchange mechanism is calculated and the characteristic
length and time scales are identified. In particular for Escherichia coli we
estimate the speedup to be 50000/N, where N is the number of repair enzymes
participating in the current exchange mechanism. Even though N is not exactly
known a speedup of order 10 is not entirely unreasonable. Furthermore upon over
expression of repair enzymes the detection time only varies as one over the
squareroot of N and not as 1/N. This behavior is of interest in assessing the
impact of stress full and radioactive environments on individual cell mutation
rates.


  A simulation approach to the stochastic growth of bacterial towers is
presented, in which a non-uniform and finite nutrient supply essentially
determines the emerging structure through elementary chemotaxis. The method is
based on cellular automata and we use simple, microscopic, local rules for
bacterial division in nutrient-rich surroundings. Stochastic nutrient
diffusion, while not crucial to the dynamics of the total population, is
influential in determining the porosity of the bacterial tower and the
roughness of its surface. As the bacteria run out of food, we observe an
exponentially rapid saturation to a carrying capacity distribution, similar in
many respects to that found in a recently proposed phenomenological
hierarchical population model, which uses heuristic parameters and macroscopic
rules. Complementary to that phenomenological model, the simulation aims at
giving more microscopic insight into the possible mechanisms for one of the
recently much studied bacterial morphotypes, known as "towering biofilm",
observed experimentally using confocal laser microscopy. A simulation
suggesting a mechanism for biofilm resistance to antibiotics is also shown.


  We investigate the transport of proteins inside the proteasome and propose an
active transport mechanism based on a spatially asymmetric interaction
potential of peptide chains. The transport is driven by fluctuations which are
always present in such systems. We compute the peptide-size dependent transport
rate which is essential for the functioning of the proteasome. In agreement
with recent experiments, varying temperature changes the transport mechanism
qualitatively.


  We study DNA adsorption and renaturation in a water-phenol two-phase system,
with or without shaking. In very dilute solutions, single-stranded DNA is
adsorbed at the interface in a salt-dependent manner. At high salt
concentrations the adsorption is irreversible. The adsorption of the
single-stranded DNA is specific to phenol and relies on stacking and hydrogen
bonding. We establish the interfacial nature of a DNA renaturation at a high
salt concentration. In the absence of shaking, this reaction involves an
efficient surface diffusion of the single-stranded DNA chains. In the presence
of a vigorous shaking, the bimolecular rate of the reaction exceeds the
Smoluchowski limit for a three-dimensional diffusion-controlled reaction. DNA
renaturation in these conditions is known as the Phenol Emulsion Reassociation
Technique or PERT. Our results establish the interfacial nature of PERT. A
comparison of this interfacial reaction with other approaches shows that PERT
is the most efficient technique and reveals similarities between PERT and the
renaturation performed by single-stranded nucleic acid binding proteins. Our
results lead to a better understanding of the partitioning of nucleic acids in
two-phase systems, and should help design improved extraction procedures for
damaged nucleic acids. We present arguments in favor of a role of phenol and
water-phenol interface in prebiotic chemistry. The most efficient renaturation
reactions (in the presence of condensing agents or with PERT) occur in
heterogeneous systems. This reveals the limitations of homogeneous approaches
to the biochemistry of nucleic acids. We propose a heterogeneous approach to
overcome the limitations of the homogeneous viewpoint.


  We reexamine the generalized cognitive homunculus, an organism's internalized
image of its physiological, psychological, and social state, which, when
properly adjusted, can quickly detect subtle deviations from a reference
configuration. We particularly seek to extend the treatment beyond
'language-of-thought' systems modeled as ergodic information sources. Such
extension would generate an exceedingly rich response repertoire, not limited
by fixed patterns of grammar and syntax. Rather, these would themselves be
tunable according to the changing short-term contextual demands faced by the
organism, possibly providing significant long-term adaptive advantage.


  We have computed the widely neglected first-order interaction between neutral
atoms. At interatomic separations typical of condensed media, it is nearly
equal to the 1/r^6 second-order London energy inferred from interactions in
gasses. Our results, without the exchange forces that lead to covalent bonding,
suggest that the quality of non-bonding attraction between neutral atoms of
molecules in condensed media differs from the 1/r^6 form usually ascribed to
it. If we add first-order and all second-order terms, we obtain energies nearly
five times the 1/r^6 London energies which dominate only at the atomic
separations found in gasses. For computation, we propose a practical, accurate
form of energy to replace the qualitatively inaccurate Lennard-Jones and
harmonic forms casually assumed to hold at the interatomic separations found in
condensed media.


  The Relevance Vector Machine (RVM) is a recently developed machine learning
framework capable of building simple models from large sets of candidate
features. Here, we describe a protocol for using the RVM to explore very large
numbers of candidate features, and a family of models which apply the power of
the RVM to classifying and detecting interesting points and regions in
biological sequence data. The models described here have been used successfully
for predicting transcription start sites and other features in genome
sequences.


  Background: In addition to known protein-coding genes, large amount of
apparently non-coding sequence are conserved between the human and mouse
genomes. It seems reasonable to assume that these conserved regions are more
likely to contain functional elements than less-conserved portions of the
genome. Here we used a motif-oriented machine learning method to extract the
strongest signal from a set of non-coding conserved sequences.
  Results: We successfully fitted models to reflect the non-coding sequences,
and showed that the results were quite consistent for repeated training runs.
Using the learned model to scan genomic sequence, we found that it often made
predictions close to the start of annotated genes. We compared this method with
other published promoter-prediction systems, and show that the set of promoters
which are detected by this method seems to be substantially similar to that
detected by existing methods.
  Conclusions: The results presented here indicate that the promoter signal is
the strongest single motif-based signal in the non-coding functional fraction
of the genome. They also lend support to the belief that there exists a
substantial subset of promoter regions which share common features and are
detectable by a variety of computational methods.


  A theoretical model of stabilization of a microtubule assembly due to
microtubule-associated-proteins(MAP) is presented. MAPs are assumed to bind to
the microtubule filaments, thus preventing their disintegration following
hydrolysis and enhancing further polymerization. Using mean-field rate
equations and explicit numerical simulations, we show that the density of MAP
(number of MAP per tubulin in the microtubule) has to exceed a critical value
$\rho_{c}$ to stabilize the structure against de-polymerization. At lower
densities $\rho < \rho_{c}$, the microtubule population consists mostly of
short polymers with exponentially decaying length distribution, whereas at
$\rho > \rho_{c}$ the average length increases linearly with time and the
microtubules ultimately extend to the cell boundary. Using experimentally
measured values of various parameters, the critical ratio of MAP to tubulin
required for unlimited growth is seen to be of the order of 1:100 or even
smaller.


  Understanding of how protein interaction networks (PIN) of living organisms
have evolved or are organized can be the first stepping stone in unveiling how
life works on a fundamental ground. Here we introduce a hybrid network model
composed of the yeast PIN and the protein family interaction network. The
essential ingredient of the model includes the protein family identity and its
robustness under evolution, as well as the three previously proposed ones: gene
duplication, divergence, and mutation. We investigate diverse structural
properties of our model with parameter values relevant to yeast, finding that
the model successfully reproduces the empirical data.


  Hydrophobicity is thought to be one of the primary forces driving the folding
of proteins. On average, hydrophobic residues occur preferentially in the core,
whereas polar residues tends to occur at the surface of a folded protein. By
analyzing the known protein structures, we quantify the degree to which the
hydrophobicity sequence of a protein correlates with its pattern of surface
exposure. We have assessed the statistical significance of this correlation for
several hydrophobicity scales in the literature, and find that the computed
correlations are significant but far from optimal. We show that this less than
optimal correlation arises primarily from the large degree of mutations that
naturally occurring proteins can tolerate. Lesser effects are due in part to
forces other than hydrophobicity and we quantify this by analyzing the surface
exposure distributions of all amino acids. Lastly we show that our database
findings are consistent with those found from an off-lattice hydrophobic-polar
model of protein folding.


  Evolution has provided many organisms with sophisticated sensory systems that
enable them to respond to signals in their environment. The response frequently
involves alteration in the pattern of movement, such as the chemokinesis of the
bacterium Escherichia coli, which swims by rotating its flagella. When rotated
counterclockwise (CCW) the flagella coalesce into a propulsive bundle,
producing a relatively straight ``run'', and when rotated clockwise (CW) they
fly apart, resulting in a ``tumble'' which reorients the cell with little
translocation. A stochastic process generates the runs and tumbles, and in a
chemoeffector gradient runs that carry the cell in a favorable direction are
extended. The overall structure of the signal transduction pathways is
well-characterized in E. coli, but important details are still not understood.
Only recently has a source of gain in the signal transduction network been
identified experimentally, and here we present a mathematical model based on
dynamic assembly of receptor teams that can explain this observation.