  Microarray data are often used to determine which genes are differentially
expressed between groups, for example, between treatment and control groups.
There are methods of determining which genes have a high probability of
differential expression, but those methods depend on the estimation of
probability densities. Theoretical results have shown such estimation to be
unreliable when high-probability genes are identified.
  The genes that are probably differentially expressed can be found using
decision theory instead of density estimation. Simulations show that the
proposed decision-theoretic method is much more reliable than a
density-estimation method. The proposed method is used to determine which genes
to consider differentially expressed between patients with different types of
cancer.
  The proposed method determines which genes have a high probability of
differential expression. It can be applied to data sets that have replicate
microarrays in each of two or more groups of patients or experiments.


  Many genomic experiments, notably microarray experiments seeking to detect
differential gene expression, involve calculating a large number of p-values.
This leads to the multiple testing problem: when the number of null hypotheses
is large, the probability of accepting at least one false alternative
hypothesis is often much greater than the significance level of the tests,
which tends to mislead investigators. Software called HighProbability provides
a simple, fast, reliable solution to the multiple testing problem, with
applications to many areas of bioinformatics. For example, in a microarray
study, HighProbability can determine which genes are probably differentially
expressed. Given a set of p-values not adjusted for multiple testing,
HighProbability determines which ones are low enough to imply a high
probability of the truth of their alternative hypotheses. The set of p-values
may be determined by conventional hypothesis testing or by random permutations
using existing R or S-PLUS software. HighProbability is freely available under
license through http://www.davidbickel.com . Coded in S, HighProbability
currently requires an installation of R or S-PLUS, but the algorithm is short
enough for fast implementation in non-S languages as well.


  Cellular metabolism, the integrated interconversion of thousands of metabolic
substrates through enzyme-catalyzed biochemical reactions, is the most
investigated complex intercellular web of molecular interactions. While the
topological organization of individual reactions into metabolic networks is
increasingly well understood, the principles governing their global functional
utilization under different growth conditions pose many open questions. We
implement a flux balance analysis of the E. coli MG1655 metabolism, finding
that the network utilization is highly uneven: while most metabolic reactions
have small fluxes, the metabolism's activity is dominated by several reactions
with very high fluxes. E. coli responds to changes in growth conditions by
reorganizing the rates of selected fluxes predominantly within this high flux
backbone. The identified behavior likely represents a universal feature of
metabolic activity in all cells, with potential implications to metabolic
engineering.


  A probability model is presented for the dynamics of mutation-selection
balance in a haploid infinite-population infinite-sites setting sufficiently
general to cover mutation-driven changes in full age-specific demographic
schedules. The model accommodates epistatic as well as additive selective
costs. Closed form characterizations are obtained for solutions in finite time,
along with proofs of convergence to stationary distributions and a proof of the
uniqueness of solutions in a restricted case. Examples are given of
applications to the biodemography of aging, including instabilities in current
formulations of mutation accumulation.


  We study the effects of the sequence on the propagation of nonlinear
excitations in simple models of DNA, and how those effects are modified by
noise. Starting from previous results on soliton dynamics on lattices defined
by aperiodic potentials, [F. Dom\'\i nguez-Adame {\em et al.}, Phys. Rev. E
{\bf 52}, 2183 (1995)], we analyze the behavior of lattices built from real DNA
sequences obtained from human genome data. We confirm the existence of
threshold forces, already found in Fibonacci sequences, and of stop positions
highly dependent on the specific sequence. Another relevant conclusion is that
the effective potential, a collective coordinate formalism introduced by
Salerno and Kivshar [Phys. Lett. A {\bf 193}, 263 (1994)] is a useful tool to
identify key regions that control the behaviour of a larger sequence. We then
study how the fluctuations can assist the propagation process by helping the
excitations to escape the stop positions. Our conclusions point out to
improvements of the model which look promising to describe mechanical
denaturation of DNA. Finally, we also consider how randomly distributed energy
focus on the chain as a function of the sequence.


  Noise, through its interaction with the nonlinearity of the living systems,
can give rise to counter-intuitive phenomena such as stochastic resonance,
noise-delayed extinction, temporal oscillations, and spatial patterns. In this
paper we briefly review the noise-induced effects in three different
ecosystems: (i) two competing species; (ii) three interacting species, one
predator and two preys, and (iii) N-interacting species. The transient dynamics
of these ecosystems are analyzed through generalized Lotka-Volterra equations
in the presence of multiplicative noise, which models the interaction between
the species and the environment. The interaction parameter between the species
is random in cases (i) and (iii), and a periodical function, which accounts for
the environmental temperature, in case (ii). We find noise-induced phenomena
such as quasi-deterministic oscillations, stochastic resonance, noise-delayed
extinction, and noise-induced pattern formation with nonmonotonic behaviors of
patterns areas and of the density correlation as a function of the
multiplicative noise intensity. The asymptotic behavior of the time average of
the \emph{$i^{th}$} population when the ecosystem is composed of a great number
of interacting species is obtained and the effect of the noise on the
asymptotic probability distributions of the populations is discussed.


  Gene regulatory network of lambda phage is one the best studied model systems
in molecular biology. More 50 years of experimental study has provided a
tremendous amount of data at all levels: physics, chemistry, DNA, protein, and
function. However, its stability and robustness for both wild type and mutants
has been a notorious theoretical/mathematical problem. In this paper we report
our successful calculation on the properties of this gene regulatory network.
We believe it is of its first kind. Our success is of course built upon
numerous previous theoretical attempts, but following 3 features make our
modeling uniqu:
  1) A new modeling method particular suitable for stability and robustness
study;
  2) Paying a close attention to the well-known difference of in vivo and in
vitro;
  3) Allowing more important role for noise and stochastic effect to play.
  The last two points have been discussed by two of us (Ao and Yin,
cond-mat/0307747), which we believe would be enough to make some of previous
theoretical attempts successful, too. We hope the present work would stimulate
a further interest in the emerging field of gene regulatory network.


  Demographic data and recent experiments verify earlier predictions that
mortality has short (few percent of the life span) memory of the previous life
history, may be significantly decreased, reset to its value at a much younger
age, and (until certain age) eliminated. Such mortality dynamics is
demonstrated to be characteristic only of evolutionary unprecedented protected
populations. When conditions improve, their mortality decreases stepwise. At
crossovers the rate of decrease rapidly changes. The crossovers manifest the
edges of the stairs in the universal ladder of rapid mortality adjustment to
changing conditions. Mortality is dominated by the established universal law
which reduces it to few biologically explicit parameters and which is verified
with human and fly mortality data. Specific experiments to test universality of
the law for other animals, and to unravel the mechanism of stepwise life
extension, are suggested.


  Well protected human and laboratory animal populations with abundant
resources are evolutionary unprecedented. Physical approach, which takes
advantage of their extensively quantified mortality, establishes that its
dominant fraction yields the exact law, whose universality from yeast to humans
is unprecedented, and suggests its unusual mechanism. Singularities of the law
demonstrate new kind of stepwise adaptation. The law proves that universal
mortality is an evolutionary byproduct, which at any age is reversible,
independent of previous life history, and may be disposable. Recent experiments
verify these predictions. Life expectancy may be extended, arguably to
immortality, by relatively small and universal biological amendments in the
animals. Indeed, it doubled with improving conditions in humans; increased
2.4-fold with genotype change in Drosophila, and 6-fold (to 430 years in human
terms), with no apparent loss in health and vitality, in nematodes with a small
number of perturbed genes and tissues. The law suggests a physical mechanism of
the universal mortality and its regulation.


  We propose a model based on coupled multiplicative stochastic processes to
understand the dynamics of competing species in an ecosystem. This process can
be conveniently described by a Fokker-Planck equation. We provide an analytical
expression for the marginalized stationary distribution. Our solution is found
in excellent agreement with numerical simulations and compares rather well with
observational data from tropical forests.


  In recent we introduced, developed and established a new concept, model,
methodology and principle for studying human longevity in terms of demographic
basis. We call the new model the "Weon model", which is a general model
modified from the Weibull model with an age-dependent shape parameter to
describe human survival and mortality curves. We demonstrate the application of
the Weon model to the mortality dynamics and the mathematical limit of
longevity (the mortality rate to be mathematically zero, implying a maximum
longevity) in the Section I. The mathematical limit of longevity can be induced
by the mortality dynamics in nature. As a result, we put forward the
complementarity principle, which explains the recent paradoxical trends that
the mathematical limit decreases as the longevity increases, in the Section II.
Our findings suggest that the human longevity can be limited by the
complementarity principle.


  We perform a reverse engineering from the ``extended Spellman data'',
consisting of 6178 mRNA levels measured by microarrays at 73 instances in four
time series during the cell cycle of the yeast Saccharomyces cerevisae. By
assuming a linear model of the genetic regulatory network, and imposing an
extra constraint (the Lasso), we obtain a unique inference of coupling
parameters. These parameters are transfered into an adjacent matrix, which is
analyzed with respect to topological properties and biological relevance. We
find a very broad distribution of outdegrees in the network, compatible with
earlier findings for biological systems and totally incompatible with a random
graph, and also indications of modules in the network. Finally, we show there
is an excess of genes coding for transcription factors among the genes of
highest outdegrees, a fact which indicates that our approach has biological
relevance.


  We present a model for the growth of West Nile virus in mosquito and bird
populations based on observations of the initial epidemic in the U.S. Increase
of bird mortality as a result of infection, which is a feature of the epidemic,
is found to yield an effect which is observable in principle, viz., periodic
variations in the extent of infection. The vast difference between mosquito and
bird lifespans, another peculiarity of the system, is shown to lead to
interesting consequences regarding delay in the onset of the steady-state
infection. An outline of a framework is provided to treat mosquito diffusion
and bird migration.


  Biological networks, such as cellular metabolic pathways or networks of
corticocortical connections in the brain, are intricately organized, yet
remarkably robust toward structural damage. Whereas many studies have
investigated specific aspects of robustness, such as molecular mechanisms of
repair, this article focuses more generally on how local structural features in
networks may give rise to their global stability. In many networks the failure
of single connections may be more likely than the extinction of entire nodes,
yet no analysis of edge importance (edge vulnerability) has been provided so
far for biological networks. We tested several measures for identifying
vulnerable edges and compared their prediction performance in biological and
artificial networks. Among the tested measures, edge frequency in all shortest
paths of a network yielded a particularly high correlation with vulnerability,
and identified inter-cluster connections in biological but not in random and
scale-free benchmark networks. We discuss different local and global network
patterns and the edge vulnerability resulting from them.


  Based on our physical and biological studies we have recently developed a
mathematical framework for the analysis of nonlnear dynamics. We call this
framework the dynamical structure analysis. It has four dynamical elements:
potential landscape, transverse matrix, descendant matrix, and stochastic
drive. In particular, the importance and the existence of the potential
landscape is emphasized.
  The dynamical structure analysis is illustrated in detail by the study of
stability, robustness, and efficiency of the simplest gene regulatory network
of phage lambda.


  Accurate demographic functions help scientists define and understand
longevity. We summarize a new demographic model, the Weon model, and show the
application to the demographic data for Switzerland (1876-2002). Particularly,
the Weon model simply defines the maximum longevity, which is induced in nature
by the mortality dynamics. In this study, we reconsider the definition of the
maximum longevity and the effectiveness for longevity by the combined effect of
the survival and mortality functions. The results suggest that the mortality
function should be zero at the maximum longevity, since the density function is
zero but the survival function is not zero. Furthermore, the effectiveness for
longevity can be maximized at the characteristic life by the complementarity
between the survival and mortality functions, which suggests that there may be
two parts of rectangularization for longevity. The historical trends for
Switzerland (1876-2002) implies that there may be a fundamental limiting force
to restrict the increase of the effectiveness. As a result, it seems that the
density function is essential to define and understand the mortality dynamics,
the maximum longevity, the effectiveness for longevity, the paradigm of
rectangularization and the historical trends of the effectiveness by the
complementarity between the survival and mortality functions.


  The transition from a normal to cancerous cell requires a number of highly
specific mutations that affect cell cycle regulation, apoptosis,
differentiation, and many other cell functions. One hallmark of cancerous
genomes is genomic instability, with mutation rates far greater than those of
normal cells. In microsatellite instability (MIN tumors), these are often
caused by damage to mismatch repair genes, allowing further mutation of the
genome and tumor progression. These mutation rates may lie near the error
catastrophe found in the quasispecies model of adaptive RNA genomes, suggesting
that further increasing mutation rates will destroy cancerous genomes. However,
recent results have demonstrated that DNA genomes exhibit an error threshold at
mutation rates far lower than their conservative counterparts. Furthermore,
while the maximum viable mutation rate in conservative systems increases
indefinitely with increasing master sequence fitness, the semiconservative
threshold plateaus at a relatively low value. This implies a paradox, wherein
inaccessible mutation rates are found in viable tumor cells. In this paper, we
address this paradox, demonstrating an isomorphism between the conservatively
replicating (RNA) quasispecies model and the semiconservative (DNA) model with
post-methylation DNA repair mechanisms impaired. Thus, as DNA repair becomes
inactivated, the maximum viable mutation rate increases smoothly to that of a
conservatively replicating system on a transformed landscape, with an upper
bound that is dependent on replication rates. We postulate that inactivation of
post-methylation repair mechanisms are fundamental to the progression of a
tumor cell and hence these mechanisms act as a method for prevention and
destruction of cancerous genomes.


  We first review how to determine the rate of vibrational energy relaxation
(VER) using perturbation theory. We then apply those theoretical results to the
problem of VER of a CD stretching mode in the protein cytochrome c. We model
cytochrome c in vacuum as a normal mode system with the lowest-order anharmonic
coupling elements. We find that, for the ``lifetime'' width parameter $\gamma=3
\sim 30$ cm$^{-1}$, the VER time is $0.2 \sim 0.3$ ps, which agrees rather well
with the previous classical calculation using the quantum correction factor
method, and is consistent with spectroscopic experiments by Romesberg's group.
We decompose the VER rate into separate contributions from two modes, and find
that the most significant contribution, which depends on the ``lifetime'' width
parameter, comes from those modes most resonant with the CD vibrational mode.


  Here we postulate three laws which form a mathematical framework to capture
the essence of Darwinian evolutionary dynamics. The second law is most
quantitative and is explicitly expressed by a unique form of stochastic
differential equation. A precise definition of Wright's adaptive landscape is
given and a new and consistent interpretation of Fisher's fundamental theorem
of natural selection is provided. Based on a recently discovered theorem the
generality of the proposed laws is illustrated by an explicit demonstration of
their equivalence to a general conventional non-equilibrium dynamics
formulation. The proposed laws provide a coherence framework to discuss several
current evolutionary problems, such as speciation and stability, and gives a
firm base for the application of statistical physics tools in Darwinian
dynamics.


  We study how a fraction of a population should be vaccinated to most
efficiently top epidemics. We argue that only local information (about the
neighborhood of specific vertices) is usable in practice, and hence we consider
only local vaccination strategies. The efficiency of the vaccination strategies
is investigated with both static and dynamical measures. Among other things we
find that the most efficient strategy for many real-world situations is to
iteratively vaccinate the neighbor of the previous vaccinee that has most links
out of the neighborhood.


  The fact that humans cooperate with non-kin in large groups, or with people
they will never meet again, is a long-standing evolutionary puzzle with
profound implications. Cooperation is linked to altruism, the capacity to
perform costly acts that confer benefits on others. Theoretical approaches had
so far disregarded costly acts that do not yield future benefits for the
altruist, either directly or indirectly. Recently, strong reciprocity, i.e.,
the predisposition to cooperate with others and to punish non-cooperators at
personal cost, has been proposed as a schema for understanding altruism in
humans. While behavioral experiments support the existence of strong
reciprocity, its evolutionary origins remain unclear: group and cultural
selection are generally invoked to compensate for the negative effects that
reciprocity is assumed to have on individuals. Here we show, by means of an
agent-based model inspired on the Ultimatum Game, that selection acting on
individuals capable of other-regarding behavior can give rise to strong
reciprocity. The results, consistent with the existence of neural correlates of
fairness, are in good agreement with observations on humans and monkeys.


  Background: Exonic splice enhancers are sequences embedded within exons which
promote and regulate the splicing of the transcript in which they are located.
A class of exonic splice enhancers are the SR proteins, which are thought to
mediate interactions between splicing factors bound to the 5' and 3' splice
sites. Method and results: We present a novel strategy for analysing
protein-coding sequence by first randomizing the codons used at each position
within the coding sequence, then applying a motif-based machine learning
algorithm to compare the true and randomized sequences. This strategy
identified a collection of motifs which can successfully discriminate between
real and randomized coding sequence, including -- but not restricted to --
several previously reported splice enhancer elements. As well as successfully
distinguishing coding exons from randomized sequences, we show that our model
is able to recognize non-coding exons. Conclusions: Our strategy succeeded in
detecting signals in coding exons which seem to be orthogonal to the sequences'
primary function of coding for proteins. We believe that many of the motifs
detected here may represent binding sites for previously unrecognized proteins
which influence RNA splicing. We hope that this development will lead to
improved knowledge of exonic splice enhancers, and new developments in the
field of computational gene prediction.


  We prove the global asymptotic stability of a well-known delayed
negative-feedback model of testosterone dynamics, which has been proposed as a
model of oscillatory behavior. We establish stability (and hence the
impossibility of oscillations) even in the presence of delays of arbitrary
length.


  The possibility that the sliding motion of proteins on DNA is influenced by
the base sequence through a base pair reading interaction, is considered.
Referring to the case of the T7 RNA-polymerase, we show that the protein should
follow a noise-influenced sequence-dependent motion which deviate from the
standard random walk usually assumed. The general validity and the implications
of the results are discussed.


  Techniques recently developed for the analysis of human social networks are
applied to the social network of bottlenose dolphins living in Doubtful Sound,
New Zealand. We identify communities and subcommunities within the dolphin
population and present evidence that sex- and age-related homophily play a role
in the formation of clusters of preferred companionship. We also identify
brokers who act as links between subcommunities and who appear to be crucial to
the social cohesion of the population as a whole. The network is found to be
similar to human social networks in some respects but different in some others
such as the level of assortative mixing by degree within the population. This
difference elucidates some of the means by which the network formed and
evolves.


  In this paper, we investigate fitness landscapes (under point mutation and
recombination) from the standpoint of whether the induced evolutionary dynamics
have a "fast-slow" time scale associated with the differences in relaxation
time between local quasi-equilibria and the global equilibrium. This dynamical
behavior has been formally described in the econometrics literature in terms of
the spectral properties of the appropriate operator matrices by Simon and Ando
(1961), and we use the relations they derive to ask which fitness functions and
mutation/recombination operators satisfy these properties. It turns out that
quite a wide range of landscapes satisfy the condition (at least trivially)
under point mutation given a sufficiently low mutation rate, while the property
appears to be difficult to satisfy under genetic recombination. In spite of the
fact that Simon-Ando decomposability can be realized over fairly wide range of
parameters, it imposes a number of restrictions on which landscape
partitionings are possible. For these reasons, the Simon-Ando formalism doesn't
appear to be applicable to other forms of decomposition and aggregation of
variables that are important in evolutionary systems.


  The genetic code is considered to be universal. In order to test if some
statistical properties of the coding bacterial genome were due to inherent
properties of the genetic code, we compared the autocorrelation function, the
scaling properties and the maximum entropy of the distribution of distances of
amino acids in sequences obtained by translating protein-coding regions from
the genome of Borrelia burgdorferi, under different genetic codes. Overall our
results indicate that these properties are very stable to perturbations made by
altering the genetic code. We also discuss the evolutionary likely implications
of the present results.


  We discuss properties which must be satisfied by a genetic network in order
for it to allow differentiation.
  These conditions are expressed as follows in mathematical terms. Let $F$ be a
differentiable mapping from a finite dimensional real vector space to itself.
The signs of the entries of the Jacobian matrix of $F$ at a given point $a$
define an interaction graph, i.e. a finite oriented finite graph $G(a)$ where
each edge is equipped with a sign. Ren\'e Thomas conjectured twenty years ago
that, if $F$ has at least two non degenerate zeroes, there exists $a$ such that
$G(a)$ contains a positive circuit. Different authors proved this in special
cases, and we give here a general proof of the conjecture. In particular, we
get this way a necessary condition for genetic networks to lead to
multistationarity, and therefore to differentiation.
  We use for our proof the mathematical literature on global univalence, and we
show how to derive from it several variants of Thomas' rule, some of which had
been anticipated by Kaufman and Thomas.


  The conformation space of a 20-residue antiparallel $\beta$-sheet peptide,
sampled by molecular dynamics simulations, is mapped to a network.
Conformations are nodes of the network, and the transitions between them are
links. The conformation space network describes the significant free energy
minima and their dynamic connectivity without projections into arbitrarily
chosen reaction coordinates. As previously found for the Internet and the
World-Wide Web as well as for social and biological networks, the conformation
space network is scale-free and contains highly connected hubs like the native
state which is the most populated free energy basin. Furthermore, the native
basin exhibits a hierarchical organization which is not found for a random
heteropolymer lacking a predominant free-energy minimum. The network topology
is used to identify conformations in the folding transition state ensemble, and
provides a basis for understanding the heterogeneity of the transition state
and denaturated state ensemble as well as the existence of multiple pathway


  A fundamental question in aging research concerns the demographic
trajectories at the highest ages, especially for supercentenarians (persons
aged 110 or more). We wish to demonstrate that the Weon model enables
scientists to describe the demographic trajectories for supercentenarians. We
evaluate the average survival data from the modern eight countries and the
valid and complete data for supercentenarians from the International Database
on Longevity (Robine and Vaupel, (2002) North American Actuarial Journal 6,
54-63). The results suggest that the Weon model predicts the maximum longevity
to exist around ages 120-130, which indicates that there is an intrinsic limit
to human longevity, and that the Weon model allows the best possible
description of the demographic trajectories for supercentenarians.


  Information theoretic analysis of genetic languages indicates that the
naturally occurring 20 amino acids and the triplet genetic code arose by
duplication of 10 amino acids of class-II and a doublet genetic code having
codons NNY and anticodons $\overleftarrow{\rm GNN}$. Evidence for this scenario
is presented based on the properties of aminoacyl-tRNA synthetases, amino acids
and nucleotide bases.


  One important issue commonly encountered in the analysis of microarray data
is to decide which and how many genes should be selected for further studies.
For discriminant microarray data analyses based on statistical models, such as
the logistic regression models, gene selection can be accomplished by a
comparison of the maximum likelihood of the model given the real data,
$\hat{L}(D|M)$, and the expected maximum likelihood of the model given an
ensemble of surrogate data with randomly permuted label, $\hat{L}(D_0|M)$.
Typically, the computational burden for obtaining $\hat{L}(D_0|M)$ is immense,
often exceeding the limits of computing available resources by orders of
magnitude. Here, we propose an approach that circumvents such heavy
computations by mapping the simulation problem to an extreme-value problem. We
present the derivation of an asymptotic distribution of the extreme-value as
well as its mean, median, and variance. Using this distribution, we propose two
gene selection criteria, and we apply them to two microarray datasets and three
classification tasks for illustration.


  Alternative splicing allows an organism to make different proteins in
different cells at different times, all from the same gene. In a cell that uses
alternative splicing, the total length of all the exons is much shorter than in
a cell that encodes the same set of proteins without alternative splicing. This
economical use of exons makes genes more stable during reproduction and
development because a genome with a shorter exon length is more resistant to
harmful mutations. Genomic stability may be the reason why higher vertebrates
splice alternatively. For a broad class of alternatively spliced genes, a
formula is given for the increase in their stability.


  We have generalized our ``unified'' model of evolutionary ecology by taking
into account the possible movements of the organisms from one ``patch'' to
another within the same eco-system. We model the spatial extension of the
eco-system (i.e., the geography) by a square lattice where each site
corresponds to a distinct ``patch''. A self-organizing hierarchical food web
describes the prey-predator relations in the eco-system. The same species at
different patches have identical food habits but differ from each other in
their reproductive characteristic features. By carrying out computer
simulations up to $10^9$ time steps, we found that, depending on the values of
the set of parameters, the distribution of the lifetimes of the species can be
either exponential or a combination of power laws. Some of the other features
of our ``unified'' model turn out to be robust against migration of the
organisms.


  We study the joint effect of thermal bath fluctuations and an external noise
tuning activity of cytotoxic cells on the triggered immune response in a
growing cancerous tissue. The immune response is assumed to be primarily
mediated by effector cells that develop a cytotoxic activity against the
abnormal tissue. The kinetics of such a reaction is represented by an
enzymatic-like Michaelis-Menten two step process. Effective free-energy surface
for the process is further parametrized by the fluctuating energy barrier
between the states of high and low concentration of cancerous cells. By
analyzing the far from equilibrium escape problem across the fluctuating
potential barrier, we determine condtions of the most efficient decay kinetics
of the cancer cell-population in the presence of dichotomously fluctuating
concentration of cytotoxic cells.


  We study the effects of interspecific correlations in a biological
coevolution model in which organisms are represented by genomes of bitstrings.
We present preliminary results for this model, indicating that these
correlations do not significantly affect the statistical behavior of the
system.


  Two ageing models were analysed whether they can confirm that the dauer
mutation of the nematode helps to preserve the species. As a result the Penna
model shows that populations with dauer larvae survive bad environmental
conditions, whereas populations without it die out. In the Stauffer model the
advantage of the dauer mutation for the survival is only given under certain
conditions.


  We present a method to quantify spatial fluctuations of the substitution rate
on different length scales throughout genomes of eukaryotes. The fluctuations
on large length scales are found to be predominantly a consequence of a
coarse-graining effect of fluctuations on shorter length scales. This is
verified for both the mouse and the human genome. We also found that both
species show similar standard deviation of fluctuations even though their mean
substitution rate differs by a factor of two. Our method furthermore allows to
determine time-resolved substitution rate maps from which we can compute
auto-correlation functions in order to quantify how fast the spatial
fluctuations in substitution rate change in time.


  Many real networks can be understood as two complementary networks with two
kind of nodes. This is the case of metabolic networks where the first network
has chemical compounds as nodes and the second one has nodes as reactions. The
second network can be related to the first one by a technique called line graph
transformation (i.e., edges in an initial network are transformed into nodes).
Recently, the main topological properties of the metabolic networks have been
properly described by means of a hierarchical model. In our work, we apply the
line graph transformation to a hierarchical network and the clustering
coefficient $C(k)$ is calculated for the transformed network, where $k$ is the
node degree. While $C(k)$ follows the scaling law $C(k)\sim k^{-1.1}$ for the
initial hierarchical network, $C(k)$ scales weakly as $k^{0.08}$ for the
transformed network. These results indicate that the reaction network can be
identified as a degree-independent clustering network.


  Biochemical networks are the analog computers of life. They allow living
cells to control a large number of biological processes, such as gene
expression and cell signalling. In biochemical networks, the concentrations of
the components are often low. This means that the discrete nature of the
reactants and the stochastic character of their interactions have to be taken
into account. Moreover, the spatial distribution of the components can be of
crucial importance. However, the current numerical techniques for simulating
biochemical networks either ignore the particulate nature of matter or treat
the spatial fluctuations in a mean-field manner. We have developed a new
technique, called Green's Function Reaction Dynamics (GFRD), that makes it
possible to simulate biochemical networks at the particle level and in both
time and space. In this scheme, a maximum time step is chosen such that only
single particles or pairs of particles have to be considered. For these
particles, the Smoluchowski equation can be solved analytically using Green's
functions. The main idea of GFRD is to exploit the exact solution of the
Smoluchoswki equation to set up an event-driven algorithm. This allows GFRD to
make large jumps in time when the particles are far apart from each other.
Here, we apply the technique to a simple model of gene expression. The
simulations reveal that the scheme is highly efficient. Under biologically
relevant conditions, GFRD is up to six orders of magnitude faster than
conventional particle-based techniques for simulating biochemical networks in
time and space.


  The evolution, regulation and sustenance of biological complexity is
determined by protein-protein interaction network that is filled with dynamic
events. Recent experimental evidences point out that clustering of proteins has
a vital role in many cellular processes. Upsurge in fluorescence imaging
methods has given a new spin to our ability to probe protein interactions in
cellular and sub-cellular compartments. Despite the increasing detection
sensitivity, quantitative information that can be obtained from these imaging
methods is limited. This is primarily due to (i) the difficulty in tracking the
problem analytically and (ii) limitations in spatio-temporal resolution that
can be achieved in interrogating living cells in real time. A novel point of
view based on diffusion-driven percolative clustering is proposed here that can
plausibly shed more light on the complex issues of protein-protein
interactions. Since this model is open to computational analysis, it is
quantitative in its premise. Besides being able to analyze the phenomenon, the
power of any model is gauged by its ability to predict interesting and novel
features of the phenomenon itself, which can subsequently be tested by
additional experiments. To this end, an experimental assay based on
fluorescence lifetime imaging is proposed to verify the validity of the
percolation model.


  Polymerization of actin proteins into dynamic structures is essential to
eukaryotic cell life. This has motivated a large body of in vitro experiments
measuring polymerization kinetics of individual filaments. Here we model these
kinetics, accounting for all relevant steps revealed by experiment:
polymerization, depolymerization, random ATP hydrolysis and release of
phosphate (Pi). We relate filament growth rates to the dynamics of ATP-actin
and ADP-Pi-actin caps which develop at filament ends. At the critical
concentration of the barbed end, c_crit, we find a short ATP cap and a long
fluctuation-stabilized ADP-Pi cap. We show that growth rates and the critical
concentration at the barbed end are intimately related to cap structure and
dynamics. Fluctuations in filament lengths are described by the length
diffusion coefficient, D. Recently Fujiwara et al. [Nature Cell Biol. (2002) 4,
666] and Kuhn and Pollard [Biohys. J. (2005) 88, 1387] observed large length
fluctuations slightly above c_crit, provoking speculation that growth may
proceed by oligomeric rather than monomeric on-off events. For the single
monomer growth process we find that $D$ exhibits a pronounced peak below
c_crit}, due to filaments alternating between capped and uncapped states, a
mild version of the dynamic instability of microtubules. Fluctuations just
above c_crit are enhanced but much smaller than those reported experimentally.
Future measurements of D as a function of concentration can help identify the
origin of the observed fluctuations.


  The three-dimensional structures of two common repeat motifs
Val$^1$-Pro$^2$-Gly$^3$-Val$^4$-Gly$^5$ and
Val$^1$-Gly$^2$-Val$^3$-Pro$^4$-Gly$^5$-Val$^6$-Gly$^7$-Val$^8$-Pro$^9$ of
tropoelastin are investigated by using the multicanonical simulation procedure.
By minimizing the energy structures along the trajectory the thermodynamically
most stable low-energy microstates of the molecule are determined. The
structural predictions are in good agreement with X-ray diffraction
experiments.


  We study the bifurcations of a set of nine nonlinear ordinary differential
equations that describe the regulation of the cyclin-dependent kinase that
triggers DNA synthesis and mitosis in the budding yeast, Saccharomyces
cerevisiae. We show that Clb2-dependent kinase exhibits bistability (stable
steady states of high or low kinase activity). The transition from low to high
Clb2-dependent kinase activity is driven by transient activation of
Cln2-dependent kinase, and the reverse transition is driven by transient
activation of the Clb2 degradation machinery. We show that a four-variable
model retains the main features of the nine-variable model. In a three-variable
model exhibiting birhythmicity (two stable oscillatory states), we explore
possible effects of extrinsic fluctuations on cell cycle progression.


  Models of the mind are based on the possibility of computing in brain
microtubules. From this point of view, information processing is the
fundamental issue for understanding the brain mechanisms that produce
consciousness. The cytoskeleton polymers could store and process information
through their dynamic coupling mediated by mechanical energy. We analyze the
problem of information transfer and storage in brain microtubules, considering
them as a communication channel. We discuss the implications of assuming that
consciousness is generated by the subneuronal process.


  We introduce a model for simulating mutation of prokaryote DNA sequences.
Using that model we can then evaluated traditional techniques like parsimony
and maximum likelihood methods for computing phylogenetic relationships. We
also use the model to mimic large scale genomic changes, and use this to
evaluate multifractal and related information theory techniques which take into
account these large changes in determining phylogenetic relationships.


  This paper explores fluctuations and noise in various facets of cancer
development. The three areas of particular focus are the stochastic progression
of cells to cancer, fluctuations of the tumor size during treatment, and noise
in cancer cell signalling. We explore the stochastic dynamics of tumor growth
and response to treatment using a Markov model, and fluctutions in tumor size
in response to treatment using partial differential equations. We also explore
noise within gene networks in cancer cells, and noise in inter-cell signalling.


  This paper examines two methods for finding whether long-range correlations
exist in DNA: a fractal measure and a mutual information technique. We evaluate
the performance and implications of these methods in detail. In particular we
explore their use comparing DNA sequences from a variety of sources. Using
software for performing in silico mutations, we also consider evolutionary
events leading to long range correlations and analyse these correlations using
the techniques presented. Comparisons are made between these virtual sequences,
randomly generated sequences, and real sequences. We also explore correlations
in chromosomes from different species.


  We address the controversial hot question concerning the validity of the
loose-coupling versus the lever-arm models in the actomyosin dynamics by
re-interpreting and extending the washboard potential model proposed by some of
us in a previous paper. In the new theory, a loose-coupling mechanism co-exists
with the deterministic lever-arm model. The synergetic action of a random
component, originating from the harnessed thermal energy, and of the
power-stroke generated by the lever-arm classical mechanism is seen to yield an
excellent fit of the set of data obtained in T. Yanagida's laboratory on the
sliding of Myosin II heads on actin filaments under various load conditions.
Our theoretical arguments are complemented by accurate numerical simulations,
and the robustness of theory is tested via different combination of parameters
and potential profiles.


  In simple models side chains are often represented implicitly (e.g., by
spin-states) or simplified as one atom. We study side chain effects using
square lattice and tetrahedral lattice models, with explicitly side chains of
two atoms. We distinguish effects due to chirality and effects due to side
chain flexibilities, since residues in proteins are L-residues, and their side
chains adopt different rotameric states. Short chains are enumerated
exhaustively. For long chains, we sample effectively rare events (eg, compact
conformations) and obtain complete pictures of ensemble properties of these
models at all compactness region. We find that both chirality and reduced side
chain flexibility lower the folding entropy significantly for globally compact
conformations, suggesting that they are important properties of residues to
ensure fast folding and stable native structure. This corresponds well with our
finding that natural amino acid residues have reduced effective flexibility, as
evidenced by analysis of rotamer libraries and side chain rotatable bonds. We
further develop a method calculating the exact side-chain entropy for a given
back bone structure. We show that simple rotamer counting often underestimates
side chain entropy significantly, and side chain entropy does not always
correlate well with main chain packing. Among compact backbones with maximum
side chain entropy, helical structures emerges as the dominating
configurations. Our results suggest that side chain entropy may be an important
factor contributing to the formation of alpha helices for compact
conformations.


  By taking into account base-base stacking interactions we improve the
Generalized Model of Polypeptide Chain (GMPC). Based on a one-dimensional
Potts-like model with many-particle interactions, the GMPC describes the
helix-coil transition in both polypeptides and polynucleotides. In the
framework of the GMPC we show that correctly introduced nearest-neighbor
stacking interactions against the background of hydrogen bonding lead to
increased stability (melting temperature) and, unexpectedly, to decreased
cooperativity (maximal correlation length). The increase in stability is
explained as due to an additional stabilizing interaction (stacking) and the
surprising decrease in cooperativity is seen as a result of mixing of
contributions of hydrogen bonding and stacking.


  We show that the contact map of the native structure of globular proteins can
be reconstructed starting from the sole knowledge of the contact map's
principal eigenvector, and present an exact algorithm for this purpose. Our
algorithm yields a unique contact map for all 221 globular structures of
PDBselect25 of length $N \le 120$. We also show that the reconstructed contact
maps allow in turn for the accurate reconstruction of the three-dimensional
structure. These results indicate that the reduced vectorial representation
provided by the principal eigenvector of the contact map is equivalent to the
protein structure itself. This representation is expected to provide a useful
tool in bioinformatics algorithms for protein structure comparison and
alignment, as well as a promising intermediate step towards protein structure
prediction.


  Networks have been used to model many real-world phenomena to better
understand the phenomena and to guide experiments in order to predict their
behavior. Since incorrect models lead to incorrect predictions, it is vital to
have a correct model. As a result, new techniques and models for analyzing and
modeling real-world networks have recently been introduced. One example of
large and complex networks involves protein-protein interaction (PPI) networks.
We demonstrate that the currently popular scale-free model of PPI networks
fails to fit the data in several respects. We show that a random geometric
model provides a much more accurate model of the PPI data.


  Active Brownian Particles are self-propelled particles that move in a
dissipative medium subject to random forces, or noise . Additionally, they can
be confined by an external field and/or they can interact with one another. The
external field may actually be an attractive marker, for example a light field
(as in the experiment) or an energy potential or a chemical gradient (as in the
theory). The potential energy can also be the result of interparticle
attractive and/or repulsive forces summed over all particles (a mean field
potential). Four, qualitatively different motions of the particles are
possible: at small particle density their motions are approximately independent
of one another subject only to the external field and the noise, which results
in moving randomly through or performing rotational motions about a central
point in space. At increasing densities interactions play an important role and
individuals form a swarm performing several types of self-organized collective
motion. We apply this model for the description of zooplankton Daphnia swarms.
In the case of the zooplankton Daphnia (and probably many other aquatic animals
that form similar motions as well) this vortex is hydrodynamical but motivated
by the self-propelled motion of the individuals. Similar vortex-type motions
have been observed for other creatures ranging in size from bacteria to flocks
of birds and schools of fish. However, our experiment with Daphnia is unique in
that all four motions can be observed in controlled laboratory conditions with
the same animal. Moreover, the theory, presented in both continuous
differential equation and random walk forms, offers a quantitative, physically
based explanation of the four motions.


  Peaks in the probabilities of loops or bubbles, helical segments, and
unzipping ends in melting DNA are found in this article using a peak finding
method that maps the hierarchical structure of certain energy landscapes. The
peaks indicate the alternative conformations that coexist in equilibrium and
the range of their fluctuations. This yields a representation of the
conformational ensemble at a given temperature, which is illustrated in a
single diagram called a stitch profile. This article describes the methodology
and discusses stitch profiles vs. the ordinary probability profiles using the
phage lambda genome as an example.


  We show how to incorporate fluctuations of the recombination rate along the
chromosome into standard gene-genealogical models for the decorrelation of gene
histories. This enables us to determine how small-scale fluctuations
(Poissonian hot-spot model) and large-scale variations [Kong et al. (2002)] of
the recombination rate influence this decorrelation. We find that the
empirically determined large-scale variations of the recombination rate give
rise to a significantly slower decay of correlations compared to the standard,
unstructured gene-genealogical model assuming constant recombination rate. A
model with long-range recombination-rate variations and with demographic
structure (divergent population) is found to be consistent with the empirically
observed slow decorrelation of gene histories. Conversely, we show that
small-scale recombination-rate fluctuations do not alter the large-scale
decorrelation of gene histories.


  The persistence of life requires populations to adapt at a rate commensurate
with the dynamics of their environment. Successful populations that inhabit
highly variable environments have evolved mechanisms to increase the likelihood
of successful adaptation. We introduce a $64 \times 64$ matrix to quantify
base-specific mutation potential, analyzing four different replicative systems,
error-prone PCR, mouse antibodies, a nematode, and Drosophila. Mutational
tendencies are correlated with the structural evolution of proteins. In systems
under strong selective pressure, mutational biases are shown to favor the
adaptive search of space, either by base mutation or by recombination. Such
adaptability is discussed within the context of the genetic code at the levels
of replication and codon usage.


  The vertebrate immune system is a wonder of modern evolution. Occasionally,
however, correlations within the immune system lead to inappropriate
recruitment of preexisting T cells against novel viral diseases. We present a
random energy theory for the correlations in the naive and memory T cell immune
responses. The non-linear susceptibility of the random energy model to
structural changes captures the correlations in the immune response to mutated
antigens. We show how the sequence-level diversity of the T cell repertoire
drives the dynamics of the immune response against mutated viral antigens.


  Much of the on-going statistical analysis of DNA sequences is focused on the
estimation of characteristics of coding and non-coding regions that would
possibly allow discrimination of these regions. In the current approach, we
concentrate specifically on genes and intergenic regions. To estimate the level
and type of correlation in these regions we apply various statistical methods
inspired from nonlinear time series analysis, namely the probability
distribution of tuplets, the Mutual Information and the Identical Neighbour
Fit. The methods are suitably modified to work on symbolic sequences and they
are first tested for validity on sequences obtained from well--known simple
deterministic and stochastic models. Then they are applied to the DNA sequence
of chromosome 1 of {\em arabidopsis thaliana}. The results suggest that
correlations do exist in the DNA sequence but they are weak and that intergenic
sequences tend to be more correlated than gene sequences. The use of
statistical tests with surrogate data establish these findings in a rigorous
statistical manner.


  The Gompertz model since 1825 has significantly contributed to interpretation
of ageing in biological and social sciences. However, in modern research
findings, it is clear that the Gompertz model is not successful to describe the
whole demographic trajectories. In this letter, a new demographic model is
introduced especially to describe human demographic trajectories, for example,
for Sweden (2002). The new model is derived from the Weibull model with an
age-dependent shape parameter, which seems to indicate the dynamical aspects of
biological systems for longevity. We will discuss the origin of the
age-dependent shape parameter. Finally, the new model presented here has
significant potential to change our understanding of definition of maximum
longevity.


  Are there limits to human longevity? We suggest a new demographic model to
describe human demographic trajectories. Specifically, the model mathematically
defines the limits of longevity. Through the demographic analysis of trends for
Sweden (between 1751 and 2002), Switzerland (between 1876 and 2002) and Japan
(between 1950 and 1999), which are the longest-lived countries, we would like
to demonstrate whether or not there is the ultimate limit to longevity. We
analyse the trends of new demographic indicators, the characteristic life and
the shape parameter, and calculate the mathematical limits of longevity. We
find out the surprising phenomenon that the mathematical limits of longevity
decrease as the longevity tendency increases in recent decades. These
paradoxical trends will be explained by the complementarity of longevity, which
is attributable to the nature of biological systems for longevity. According to
the regression analysis, the ultimate limit for humans is estimated to be
approximately 124 years, which may be considered to be the biological limit for
humans.


  Observations that rates of molecular evolution vary widely within and among
lineages have cast doubts upon the existence of a single molecular clock.
Differences in the timing of evolutionary events estimated from genetic and
fossil evidence have raised further questions about the existence of molecular
clocks and their use. Here we present a model of nucleotide substitution that
combines new theory on metabolic rate with the now classic neutral theory of
molecular evolution. The model quantitatively predicts rate heterogeneity, and
reconciles differences in molecular- and fossil-estimated dates of evolutionary
events. Model predictions are supported by extensive data from mitochondrial
and nuclear genomes. By accounting for the effects of body size and temperature
on metabolic rate, a single molecular clock explains heterogeneity in rates of
nucleotide substitution in different genes, taxa, and thermal environments.
This model suggests that there is indeed a general molecular clock, as
originally proposed by Zuckerkandl and Pauling, but that it ticks at a constant
substitution rate per unit mass-specific metabolic energy rather than per unit
time. More generally, the model suggests that body size and temperature combine
to control the overall rate of evolution through their effects on metabolism.


  Based on experimental observations in \textit{Daphnia}, we introduce an
agent-based model for the motion of single and swarms of animals. Each agent is
described by a stochastic equation that also considers the conditions for
active biological motion. An environmental potential further reflects local
conditions for \textit{Daphnia}, such as attraction to light sources. This
model is sufficient to describe the observed cycling behavior of single
\textit{Daphnia}. To simulate vortex swarming of many \textit{Daphnia}, i.e.
the collective rotation of the swarm in one direction, we extend the model by
considering avoidance of collisions. Two different ansatzes to model such a
behavior are developed and compared. By means of computer simulations of a
multi-agent system we show that local avoidance - as a special form of
asymmetric repulsion between animals - leads to the emergence of a vortex
swarm. The transition from uncorrelated rotation of single agents to the vortex
swarming as a function of the swarm size is investigated. Eventually, some
evidence of avoidance behavior in \textit{Daphnia} is provided by comparing
experimental and simulation results for two animals.


  A simple model of activatory-inhibitory interactions controlling the activity
of agents (substrates) through a "saturated response" dynamical rule in a
scale-free network is thoroughly studied. After discussing the most remarkable
dynamical features of the model, namely fragmentation and multistability, we
present a characterization of the temporal (periodic and chaotic) fluctuations
of the quasi-stasis asymptotic states of network activity. The double (both
structural and dynamical) source of entangled complexity of the system temporal
fluctuations, as an important partial aspect of the Correlation
Structure-Function problem, is further discussed to the light of the numerical
results, with a view on potential applications of these general results.


  The current picture of bacterial evolution is based largely on studies of 16S
rRNA. However, this is just one gene. It is known that horizontal gene transfer
can occur between bacterial species, although the frequency and implications of
this are not fully understood. If horizontal transfer were frequent, there
would be no single evolutionary tree for bacteria because each gene would
follow a different tree. We carried out phylogenetic analyses of rRNA and tRNA
genes from Proteobacteria (a diverse group for which many complete genome
sequences are available) using RNA-specific phylogenetic methods that account
for the conservation of the secondary structure. We compared trees for 16S rRNA
and 23S rRNA with those derived from concatenated alignments of 29 tRNA genes
that are found in all the genomes studied. The tRNA genes are scattered
throughout the genomes, and would not follow the same evolutionary history if
horizontal transfer were frequent. Nevertheless, the tRNA tree is consistent
with the rRNA tree in most respects. Minor differences can almost all be
attributed to uncertainty or unreliability of the phylogenetic method. We
therefore conclude that tRNA genes give a coherent picture of the phylogeny of
the organisms, and that horizontal transfer of tRNAs is too rare to obscure the
signal of the organismal tree. Some tRNA genes are not present in all the
genomes. We discuss possible explanations for the observed patterns of presence
and absence of genes: these involve gene deletion, gene duplication, and
mutations in the tRNA anticodons.


  The PHASE software package allows phylogenetic tree construction with a
number of evolutionary models designed specifically for use with RNA sequences
that have conserved secondary structure. Evolution in the paired regions of
RNAs occurs via compensatory substitutions, hence changes on either side of a
pair are correlated. Accounting for this correlation is important for
phylogenetic inference because it affects the likelihood calculation. In the
present study we use the complete set of tRNA and rRNA sequences from 69
complete mammalian mitochondrial genomes. The likelihood calculation uses two
evolutionary models simultaneously for different parts of the sequence: a
paired-site model for the paired sites and a single-site model for the unpaired
sites. We use Bayesian phylogenetic methods and a Markov chain Monte Carlo
algorithm is used to obtain the most probable trees and posterior probabilities
of clades. The results are well resolved for almost all the important branches
on the mammalian tree. They support the arrangement of mammalian orders within
the four supra-ordinal clades that have been identified by studies of much
larger data sets mainly comprising nuclear genes. Groups such as the hedgehogs
and the murid rodents, which have been problematic in previous studies with
mitochondrial proteins, appear in their expected position with the other
members of their order. Our choice of genes and evolutionary model appears to
be more reliable and less subject to biases caused by variation in base
composition than previous studies with mitochondrial genomes.


  A previously proved theorem gives sufficient conditions for an estimator of
the false discovery rate (FDR) to conservatively converge to the FDR with
probability 1 as the number of hypothesis tests increases, even for small
sample sizes. It does not follow that several thousand tests ensure that the
estimator has moderate variance under those conditions. In fact, they can hold
even if the test statistics have long-range correlations, which yield
unacceptably wide confidence intervals, as observed in genomic data when there
are 8 or 16 individuals (microarrays) per group. Thus, informative FDR
estimation will include some measure of its reliability.


  Motivation: The reconstruction of gene networks from gene expression
microarrays is gaining popularity as methods improve and as more data become
available. The reliability of such networks could be judged by the probability
that a connection between genes is spurious, resulting from chance fluctuations
rather than from a true biological relationship. Results: Unlike the false
discovery rate and positive false discovery rate, the decisive false discovery
rate (dFDR) is exactly equal to a conditional probability without assuming
independence or the randomness of hypothesis truth values. This property is
useful not only in the common application to the detection of differential gene
expression, but also in determining the probability of a spurious connection in
a reconstructed gene network. Estimators of the dFDR can estimate each of three
probabilities: 1. The probability that two genes that appear to be associated
with each other lack such association. 2. The probability that a time ordering
observed for two associated genes is misleading. 3. The probability that a time
ordering observed for two genes is misleading, either because they are not
associated or because they are associated without a lag in time. The first
probability applies to both static and dynamic gene networks, and the other two
only apply to dynamic gene networks. Availability: Cross-platform software for
network reconstruction, probability estimation, and plotting is free from
http://www.davidbickel.com as R functions and a Java application.


  We investigate the time-evolution and steady states of the stochastic
susceptible-infected-recovered-susceptible(SIRS) epidemic model on one- and
two- dimensional lattices. We compare the behavior of this system, obtained
from computer simulations, with those obtained from the mean-field
approximation(MFA) and pair-approximation(PA). The former(latter) approximates
higher order moments in terms of first(second) order ones. We find that the PA
gives consistently better results than the MFA. In one dimension the
improvement is even qualitative.


  Anaerobic glycolysis in yeast perturbed by the reduction of xenobiotic
ketones is studied numerically in two models which possess the same topology
but different levels of complexity. By comparing both models' predictions for
concentrations and fluxes as well as steady or oscillatory temporal behavior we
answer the question what phenomena require what kind of minimum model
abstraction. While mean concentrations and fluxes are predicted in agreement by
both models we observe different domains of oscillatory behavior in parameter
space. Generic properties of the glycolytic response to ketones are discussed.


  In the nuclei of eukaryotic cells, DNA is packaged through several levels of
compaction in an orderly retrievable way that enables the correct regulation of
gene expression. The functional dynamics of this assembly involves the
unwinding of the so-called 30 nm chromatin fiber and accordingly imposes strong
topological constraints. We present a general method for computing both the
twist and the writhe of any winding pattern. An explicit derivation is
implemented for the chromatin fiber which provides the linking number of DNA in
eukaryotic chromosomes. We show that there exists one and only one unwinding
path which satisfies both topological and mechanical constraints that DNA has
to deal with during condensation/decondensation processes.


  In this study, we explore nucleation and the transition state ensemble of the
ribosomal protein S6 using a Monte Carlo Go model in conjunction with
restraints from experiment. The results are analyzed in the context of
extensive experimental and evolutionary data. The roles of individual residues
in the folding nucleus are identified and the order of events in the S6 folding
mechanism is explored in detail. Interpretation of our results agrees with, and
extends the utility of, experiments that shift f-values by modulating
denaturant concentration and presents strong evidence for the realism of the
mechanistic details in our Monte Carlo Go model and the structural
interpretation of experimental f-values. We also observe plasticity in the
contacts of the hydrophobic core that support the specific nucleus. For S6,
which binds to RNA and protein after folding, this plasticity may result from
the conformational flexibility required to achieve biological function. These
results present a theoretical and conceptual picture that is relevant in
understanding the mechanism of nucleation in protein folding.


  Understanding the observed variability in the number of homologs of a gene is
a very important, unsolved problem that has broad implications for research
into co-evolution of structure and function, gene duplication, pseudogene
formation and possibly for emerging diseases. Here we attempt to define and
elucidate the reasons behind this observed unevenness in sequence space. We
present evidence that sequence variability and functional diversity of a gene
or fold family is influenced by certain quantitative characteristics of the
protein structure that reflect potential for sequence plasticity i.e. the
ability to accept mutation without losing thermodynamic stability.


  We report translocation experiments on double-strand DNA through a silicon
oxide nanopore. Samples containing DNA fragments with seven different lengths
between 2000 to 96000 basepairs have been electrophoretically driven through a
10 nm pore. We find a power-law scaling of the translocation time versus
length, with an exponent of 1.26 $\pm$ 0.07. This behavior is qualitatively
different from the linear behavior observed in similar experiments performed
with protein pores. We address the observed nonlinear scaling in a theoretical
model that describes experiments where hydrodynamic drag on the section of the
polymer outside the pore is the dominant force counteracting the driving. We
show that this is the case in our experiments and derive a power-law scaling
with an exponent of 1.18, in excellent agreement with our data.


  We study the spreading of families in two-dimensional multispecies
predator-prey systems, in which species cyclically dominate each other. In each
time step randomly chosen individuals invade one of the nearest sites of the
square lattice eliminating their prey. Initially all individuals get a
family-name which will be carried on by their descendants. Monte Carlo
simulations show that the systems with several species (N=3,4,5) are
asymptotically approaching the behavior of the voter model, i.e., the survival
probability of families, the mean-size of families and the mean-square distance
of descendants from their ancestor exhibit the same scaling behavior. The
scaling behavior of the survival probability of families has a logarithmic
correction. In case of the voter model this correction depends on the number of
species, while cyclic predator-prey models behave like the voter model with
infinite species. It is found that changing the rates of invasions does not
change this asymptotic behavior. As an application a three-species system with
a fourth species intruder is also discussed.


  This article introduces the physics of information in the context of
molecular biology and genomics. Entropy and information, the two central
concepts of Shannon's theory of information and communication, are often
confused with each other but play transparent roles when applied to statistical
ensembles (i.e., identically prepared sets) of symbolic sequences. Such an
approach can distinguish between entropy and information in genes, predict the
secondary structure of ribozymes, and detect the covariation between residues
in folded proteins. We also review applications to molecular sequence and
structure analysis, and introduce new tools in the characterization of
resistance mutations, and in drug design.


  Bacteria (e.g. E. Coli) are very sensitive to certain chemoattractants (e.g.
asparate) which they themselves produce. This leads to chemical instabilities
in a uniform population. We discuss here the different case of a single
bacterium, following the general scheme of Brenner, Levitov and Budrene. We
show that in one and two dimensions (in a capillary or in a thin film) the
bacterium can become self-trapped in its cloud of attractant. This should occur
if a certain coupling constant $g$ is larger than unity. We then estimate the
reduced diffusion D_eff of the bacterium in the strong coupling limit, and find
D_eff ~ 1/g.


  Function of proteins or a network of interacting proteins often involves
communication between residues that are well separated in sequence. The classic
example is the participation of distant residues in allosteric regulation.
Bioinformatic and structural analysis methods have been introduced to infer
residues that are correlated. Recently, increasing attention has been paid to
obtain the sequence properties that determine the tendency of disease related
proteins (Abeta peptides, prion proteins, transthyretin etc.) to aggregate and
form fibrils. Motivated in part by the need to identify sequence
characteristics that indicate a tendency to aggregate, we introduce a general
method that probes covariations in charged residues along the sequence in a
given protein family. The method, which involves computing the Sequence
Correlation Entropy (SCE) using the quenched probability Psk(i,j) of finding a
residue pair at a given sequence separation sk, allows us to classify protein
families in terms of their SCE. Our general approach may be a useful way in
obtaining evolutionary covariations of amino acid residues on a genome wide
level.


  We have exactly enumerated all sequences and conformations of HP proteins
with chains of up to 19 monomers on the simple cubic lattice. For two variants
of the hydrophobic-polar (HP) model, where only two types of monomers are
distinguished, we determined and statistically analyzed designing sequences,
i.e., sequences that have a non-degenerate ground state. Furthermore we were
interested in characteristic thermodynamic properties of HP proteins with
designing sequences. In order to be able to perform these exact studies, we
applied an efficient enumeration method based on contact sets.


  We present an analysis of the effects of global topology on the structural
stability of folded proteins in thermal equilibrium with a heat bath. For a
large class of single domain proteins, we computed the harmonic spectrum within
the Gaussian Network Model (GNM) and determined the spectral dimension, a
parameter describing the low frequency behaviour of the density of modes. We
find a surprisingly strong correlation between the spectral dimension and the
number of amino acids of the protein. Considering that larger spectral
dimension value relate to more topologically compact folded state, our results
indicate that for a given temperature and length of the protein, the folded
structure corresponds to the less compact folding compatible with thermodynamic
stability.


  Can complex engineered and biological networks be coarse-grained into smaller
and more understandable versions in which each node represents an entire
pattern in the original network? To address this, we define coarse-graining
units (CGU) as connectivity patterns which can serve as the nodes of a
coarse-grained network, and present algorithms to detect them. We use this
approach to systematically reverse-engineer electronic circuits, forming
understandable high-level maps from incomprehensible transistor wiring: first,
a coarse-grained version in which each node is a gate made of several
transistors is established. Then, the coarse-grained network is itself
coarse-grained, resulting in a high-level blueprint in which each node is a
circuit-module made of multiple gates. We apply our approach also to a
mammalian protein-signaling network, to find a simplified coarse-grained
network with three main signaling channels that correspond to cross-interacting
MAP-kinase cascades. We find that both biological and electronic networks are
'self-dissimilar', with different network motifs found at each level. The
present approach can be used to simplify a wide variety of directed and
nondirected, natural and designed networks.


  We present the results of an exact analysis of a model energy landscape of a
protein to clarify the notion of the transition state and the physical meaning
of the $\phi$ values determined in protein engineering experiments. We
benchmark our findings to various theoretical approaches proposed in the
literature for the identification and characterization of the transition state.


  In this paper we consider the problem of RNA folding with pseudoknots. We use
a graphical representation in which the secondary structures are described by
planar diagrams. Pseudoknots are identified as non-planar diagrams. We analyze
the non-planar topologies of RNA structures and propose a classification of RNA
pseudoknots according to the minimal genus of the surface on which the RNA
structure can be embedded. This classification provides a simple and natural
way to tackle the problem of RNA folding prediction in presence of pseudoknots.
Based on that approach, we describe a Monte Carlo algorithm for the prediction
of pseudoknots in an RNA molecule.


  Four reasons why you might wish to read this paper: 1. We have devised a new
statistical T test to determine differentially expressed genes (DEG) in the
context of microarray experiments. This statistical test adds a new member to
the traditional T-test family. 2. An exact formula for calculating the
detection power of this T test is presented, which can also be fairly easily
modified to cover the traditional T tests. 3. We have presented an accurate yet
computationally very simple method to estimate the fraction of non-DEGs in a
set of genes being tested. This method is superior to an existing one which is
computationally much involved. 4. We approach the multiple testing problem from
a fresh angle, and discuss its relation to the classical Bonferroni procedure
and to the FDR (false discovery rate) approach. This is most useful in the
analysis of microarray data, where typically several thousands of genes are
being tested simultaneously.


  Evolutionary models measure the probability of amino acid substitutions
occurring over different evolutionary distances. We examine various
evolutionary models based on empirically derived amino acid substitution
matrices. The models are constructed using the PAM and BLOSUM amino acid
substitution matrices. We rescale these matrices by raising them to powers to
model substitution patterns that account for different evolutionary distances.
We also examine models that account for the dissimilarity of substitution rates
along a protein sequence. We compare the models by computing the likelihood of
each model across different alignments. We also present a specific example to
illustrate the subtle differences in the estimation of evolutionary distance
computed using the different models.


  We consider a simple model for the unfolding of RNA tertiary structure under
dynamic loading. The opening of such a structure is regarded as a two step
process, each corresponding to the overcoming of a single energy barrier. The
resulting two-barrier energy landscape accounts for the dependence of the
unfolding kinetics on the pulling rate. Furthermore at intermediate force, the
two barriers cannot be distinguished by the analysis of the opening kinetic,
which turns out to be dominated by a single macro-barrier, whose properties
depend non-trivially on the two single barriers. Our results suggest that in
pulling experiments on RNA molecule containing tertiary structures, the details
of the single kinetic barriers can only be obtained using a low pulling rate
value, or in the high force regime.


  The Prisoner's Dilemma, a 2-person game in which the players can either
cooperate or defect, is a common paradigm for studying the evolution of
cooperation, when individuals exhibit variable degrees of cooperation. It is
known that in the presence of spatial structure, when individuals ``play
against'' their neighbours, and ``compare to'' them, cooperative investments
can evolve to considerable levels. Here we examine the effect of increasing the
neighbourhood size: we find that the mean-field limit of no cooperation is
reached for a critical neighbourhood size of about five neighbours. We also
find the related result that in a network of players, the critical average
degree (number of neighbours) of nodes for which defection is the final state
depends only on the network topology. This critical average degree is
considerably higher for clustered networks, than for distributed random
networks. This result strengthens the argument that clustering is the mechanism
which makes the development and maintenance of the cooperation possible.


  We present a simple physical model which demonstrates that the native state
folds of proteins can emerge on the basis of considerations of geometry and
symmetry. We show that the inherent anisotropy of a chain molecule, the
geometrical and energetic constraints placed by the hydrogen bonds and sterics,
and hydrophobicity are sufficient to yield a free energy landscape with broad
minima even for a homopolymer. These minima correspond to marginally compact
structures comprising the menu of folds that proteins choose from to house
their native-states in. Our results provide a general framework for
understanding the common characteristics of globular proteins.


  Signal processing (SP) techniques convert DNA and protein sequences into
information that lead to successful drug discovery. One must, however, be aware
about the difference between information and entropy1. Eight other physical
properties of DNA and protein segments are suggested for SP analysis other than
ones already used in the literature. QSAR formulations of these properties are
suggested for ranking the amino acids that maximize efficiency of the amino
acids in proteins. Multiobjective programs are suggested for constraining or
searching the components of such sequences. Geometric maps of the networks of
proteins are preferable to scale-free descriptions in most cases. The genetic
code is presented as graphlets which show interesting correspondence to each
other, leading to possible new revelations.


  Within a case study on the protein-protein interaction network (PIN) of
Drosophila melanogaster we investigate the relation between the network's
spectral properties and its structural features such as the prevalence of
specific subgraphs or duplicate nodes as a result of its evolutionary history.
The discrete part of the spectral density shows fingerprints of the PIN's
topological features including a preference for loop structures. Duplicate
nodes are another prominent feature of PINs and we discuss their representation
in the PIN's spectrum as well as their biological implications.


  With the number of fully-sequenced genomes now well over a hundred it has
become possible to start investigating if there are any quantitative
regularities in the genetic make-up of genomes. In (physics/0307001), I
originally showed that the numbers of genes in different functional categories
scale as power laws in the total number of genes in the genome. In this chapter
I revisit these results with more recent data and go into considerable more
depth regarding the implications of these scaling laws for our understanding of
the regulatory design of cells. In addition, I further develop the evolutionary
model first proposed in (physics/0307001), which suggests that the exponents of
the observed scaling laws correspond to fundamental constants of the
evolutionary process. In particular, I put forward an hypothesis for the
approximately quadratic scaling of regulatory and signal transducing genes with
genome size.


  Ecology and evolution are inseparable. Motivated by some recent experiments,
we have developed models of evolutionary ecology from the perspective of
dynamic networks. In these models, in addition to the intra-node dynamics,
which corresponds to an individual-based population dynamics of species, the
entire network itself changes slowly with time to capture evolutionary
processes. After a brief summary of our recent published works on these network
models of eco-systems, we extend the most recent version of the model
incorporating predators that wander into neighbouring spatial patches for food.


  We focused our attention on Ca$^{2+}$ release from the endoplasmic reticulum
through a cluster of inositol 1,4,5-trisphosphate (IP$_3$) receptor channels.
The random opening and closing of these receptors introduce stochastic effects
that have been observed experimentally. Here, we present a stochastic version
of Othmer-Tang model for IP$_3$ receptor clusters. We address the average
behavior of the channels in response to IP$_3$ stimuli. We found, by stochastic
simulation, that the shape of the receptor response to IP$_3$ (fraction of open
channels versus [IP$_3$]), is affected by the cytosolic Ca$^{2+}$ level. We
also study several aspects of the stochastic properties of Ca${2+}$ release and
we compare with experimental observations.


  Well protected human and laboratory animal populations with abundant
resources are evolutionary unprecedented, and their survival far beyond
reproductive age may be a byproduct rather than tool of evolution. Physical
approach, which takes advantage of their extensively quantified mortality,
establishes that its dominant fraction yields the exact law, and suggests its
unusual mechanism. The law is universal for all animals, from yeast to humans,
despite their drastically different biology and evolution. It predicts that the
universal mortality has short memory of the life history, at any age may be
reset to its value at a significantly younger age, and mean life expectancy
extended (by biologically unprecedented small changes) from its current maximal
value to immortality. Mortality change is rapid and stepwise. Demographic data
and recent experiments verify these predictions for humans, rats, flies,
nematodes and yeast. In particular, mean life expectancy increased 6-fold (to
"human" 430 years), with no apparent loss in health and vitality, in nematodes
with a small number of perturbed genes and tissues. Universality allows one to
study unusual mortality mechanism and the ways to immortality.


  Recent developments in complex networks have paved the way to a series of
important biological insights, such as the fact that many of the essential
proteins of S. cerevisae corresponds to the so-called hubs of the respective
protein-protein interaction networks. Despite the special importance of hubs,
other types of nodes such as those corresponding to the network border, as well
as the innermost nodes, also deserve special attention. This work reports on
how the application of the concept of distance transform to networks showed
that a great deal of the innermost nodes correspond to essential proteins, with
interesting biological implications.


  With the aim to study the relationship between protein sequences and their
native structures, we adopt vectorial representations for both sequence and
structure. The structural representation is based on the Principal Eigenvector
of the fold's contact matrix (PE). As recently shown, the latter encodes
sufficient information for reconstructing the whole contact matrix. The
sequence is represented through a Hydrophobicity Profile (HP), using a
generalized hydrophobicity scale that we obtain from the principal eigenvector
of a residue-residue interaction matrix and denote it as interactivity scale.
Using this novel scale, we define the optimal HP of a protein fold, and
predict, by means of stability arguments, that it is strongly correlated with
the PE of the fold's contact matrix. This prediction is confirmed through an
evolutionary analysis, which shows that the PE correlates with the HP of each
individual sequence adopting the same fold and, even more strongly, with the
average HP of this set of sequences. Thus, protein sequences evolve in such a
way that their average HP is close to the optimal one, implying that neutral
evolution can be viewed as a kind of motion in sequence space around the
optimal HP. Our results indicate that the correlation coefficient between
N-dimensional vectors constitutes a natural metric in the vectorial space in
which we represent both protein sequences and protein structures, which we call
Vectorial Protein Space. In this way, we define a unified framework for
sequence to sequence, sequence to structure, and structure to structure
alignments. We show that the interactivity scale is nearly optimal both for the
comparison of sequences with sequences and sequences with structures.


  Rebinding of dissociated ligands from cell surface proteins can confound
quantitative measurements of dissociation rates important for characterizing
the affinity of binding interactions. This can be true also for in vitro
techniques such as surface plasmon resonance (SPR). We present experimental
results using SPR for the interaction of insulin-like growth factor-I (IGF-I)
with one of its binding proteins, IGF binding protein-3 (IGFBP-3), and show
that rebinding, even with the addition of soluble heparin in the dissociation
phase, does not exhibit the expected exponential decay characteristic of a 1:1
binding reaction. We thus consider the effect of (multiple) rebinding events
and, within a self-consistent mean-field approximation, we derive the complete
mathematical form for the fraction of bound ligand as a function of time. We
show that, except for very low surface coverage/association rate, this function
is non-exponential at all times, indicating that multiple rebinding events
strongly influence dissociation even at early times. We compare the mean-field
results with numerical simulations and find good agreement, although deviations
are measurable in certain cases. Our analysis of the IGF-I-IGFBP-3 data
indicates that rebinding is prominent for this system and that the theoretical
predictions fit the experimental data well. Our results provide a means for
analyzing SPR biosensor data where rebinding is problematic and a methodology
to do so is presented.


  In the game-theoretic model war of attrition, players are subject to an
explicit cost proportional to the duration of contests. We construct a model
where the time cost is not explicitly given, but instead depends implicitly on
the strategies of the whole population. We identify and analyse the underlying
mechanisms responsible for the implicit time cost. Each player participates in
a series of games, where those prepared to wait longer win with higher
certainty but play less frequently. The model is characterised by the ratio of
the winner's score to the loser's score, in a single game. The fitness of a
player is determined by the accumulated score from the games played during a
generation. We derive the stationary distribution of strategies under the
replicator dynamics. When the score ratio is high, we find that the stationary
distribution is unstable, with respect to both evolutionary and dynamical
stability, and the dynamics converge to a limit cycle. When the ratio is low,
the dynamics converge to the stationary distribution. For an intermediate
interval of the ratio, the distribution is dynamically but not evolutionarily
stable. Finally, the implications of our results for previous models based on
the war of attrition are discussed.


  Bistable biochemical switches are ubiquitous in gene regulatory networks and
signal transduction pathways. Their switching dynamics, however, are difficult
to study directly in experiments or conventional computer simulations, because
switching events are rapid, yet infrequent. We present a simulation technique
that makes it possible to predict the rate and mechanism of flipping of
biochemical switches. The method uses a series of interfaces in phase space
between the two stable steady states of the switch to generate transition
trajectories in a ratchet-like manner. We demonstrate its use by calculating
the spontaneous flipping rate of a symmetric model of a genetic switch
consisting of two mutually repressing genes. The rate constant can be obtained
orders of magnitude more efficiently than using brute-force simulations. For
this model switch, we show that the switching mechanism, and consequently the
switching rate, depends crucially on whether the binding of one regulatory
protein to the DNA excludes the binding of the other one. Our technique could
also be used to study rare events and non-equilibrium processes in soft
condensed matter systems.


  We use large deviation methods to calculate rates of noise-induced
transitions between states in multistable genetic networks. We analyze a
synthetic biochemical circuit, the toggle switch, and compare the results to
those obtained from a numerical solution of the master equation.


  In this Letter, we extend the semiconservative quasispecies equations to
incorporate imperfect DNA lesion repair. We study the equilibrium behavior of
this model in the limit of infinite sequence length and population size, using
a single-fitness-peak landscape for which the master genome can sustain a
finite number of lesions and remain viable. We provide a full analytical
treatment of the problem, providing a general mathematical framework as well as
the full solution for a particular class of fitness landscapes. Stochastic
simulations using finite sequence lengths and populations agree well with the
analytical results. Applications to biological systems are briefly discussed.


  We describe the representation of the chemical affinity between the
antigen-combining site of the immunoglobulin molecule and the antigen molecule
as the probability of the two molecules existing in a bound state. Our model is
based on the identification of shape attractors in the configuration space for
the joint antibody / antigen combining site sequence. We parameterize
configuration space in terms of Ramachandran angles. The shape attractors allow
us to construct a Markov chain whose steady state distribution gives rise to
the desired attachment probability. As a result we are able to delineate the
enthalpic, entropic and kinetic components of affinity and study their
interactions.


  It is widely known that the spread of the HIV virus was slower than
exponential in several populations, even at the very beginning of the epidemic.
We show that this implies a significant reduction in the effective reproductive
rate of the epidemic, and describe a general mechanism, related to the
clustering properties of the disease transmission network, capable of
explaining this reduction. Our considerations provide a new angle on polynomial
epidemic processes, and may have implications for the choice of strategy
against such epidemics.


  We study the maturation of the antibody population following primary antigen
presentation as a global optimization problem. Emphasis is placed on the
trade-off between the safety of mutations that lead to local improvements to
the antibody's affinity and the necessity of eventual mutations that result in
global reconfigurations in the antibody's shape. The model described herein
gives evidence of the underlying optimization process from which the rapidity
and consistency of the biologic response could be derived.


  We define the concept of dependence among multiple variables using maximum
entropy techniques and introduce a graphical notation to denote the
dependencies. Direct inference of information theoretic quantities from data
uncovers dependencies even in undersampled regimes when the joint probability
distribution cannot be reliably estimated. The method is tested on synthetic
data. We anticipate it to be useful for inference of genetic circuits and other
biological signaling networks.


  We investigate by means of a simple theoretical model the emergence of prime
numbers as life cycles, as those seen for some species of cicadas. The cicadas,
more precisely, the Magicicadas spend most of their lives below the ground and
then emerge and die in a short period of time. The Magicicadas display an
uncommon behavior: their emergence is synchronized and these periods are
usually prime numbers. In the current work, we develop a spatially extended
model at which preys and predators coexist and can change their evolutionary
dynamics through the occurrence of mutations. We verified that prime numbers as
life cycles emerge as a result of the evolution of the population. Our results
seem to be a first step in order to prove that the development of such strategy
is selectively advantageous, especially for those organisms that are highly
vulnerable to attacks of predators.


  We present for the first time a complex network approach to the study of the
electrical properties of single protein devices. In particular, we consider an
electronic nanobiosensor based on a G-protein coupled receptor. By adopting a
coarse grain description, the protein is modeled as a complex network of
elementary impedances. The positions of the alpha-carbon atoms of each amino
acid are taken as the nodes of the network. The amino acids are assumed to
interact electrically among them. Consequently, a link is drawn between any
pair of nodes neighboring in space within a given distance and an elementary
impedance is associated with each link. The value of this impedance can be
related to the physical and chemical properties of the amino acid pair and to
their relative distance. Accordingly, the conformational changes of the
receptor induced by the capture of the ligand, are translated into a variation
of its electrical properties. Stochastic fluctuations in the value of the
elementary impedances of the network, which mimic different physical effects,
have also been considered. Preliminary results concerning the impedance
spectrum of the network and its fluctuations are presented and discussed for
different values of the model parameters.


  We fit the Fourier transforms of solvent accessibility and hydrophobicity
profiles of a representative set of proteins to a joint multi-variable
Gaussian. This allows us to separate the intrinsic tendencies of sequence and
structure profiles from the interactions that correlate them; for example, the
$\alpha$-helix periodicity in sequence hydrophobicity is dictated by the
solvent accessibility of structures. The distinct intrinsic tendencies of
sequence and structure profiles are most pronounced at long periods, where
sequence hydrophobicity fluctuates more, while solvent accessibility
fluctuations are less than average. Interestingly, correlations between the two
profiles can be interpreted as the Boltzmann weight of the solvation energy at
room temperature.


  It is shown that an ecosystem in equilibrium is generally organized in a
state which is poised in the vicinity of a novel phase transition.


  In this paper, we examine the mechanical role of the lipid bilayer in ion
channel conformation and function with specific reference to the case of the
mechanosensitive channel of large conductance (MscL). In a recent paper
(Wiggins and Phillips, 2004), we argued that mechanotransduction very naturally
arises from lipid-protein interactions by invoking a simple analytic model of
the MscL channel and the surrounding lipid bilayer. In this paper, we focus on
improving and expanding this analytic framework for studying lipid-protein
interactions with special attention to MscL. Our goal is to generate simple
scaling relations which can be used to provide qualitative understanding of the
role of membrane mechanics in protein function and to quantitatively interpret
experimental results. For the MscL channel, we find that the free energies
induced by lipid-protein interaction are of the same order as the free energy
differences between conductance states measured by Sukharev et al. (1999). We
therefore conclude that the mechanics of the bilayer plays an essential role in
determining the conformation and function of the channel. Finally, we compare
the predictions of our model to experimental results from the recent
investigations of the MscL channel by Perozo et al. (2002), Powl et al. (2003),
Yoshimura et al. (2004), and others and suggest a suite of new experiments.