  Upon growth on dinitrogen, the filamentous cyanobacterium Nostoc PCC 7120
initiates metabolic and morphological changes. We analyzed the expression of
1249 genes from major metabolic categories under nitrogen fixing and
non-nitrogen fixing growth. The expression data were correlated with potential
target secondary structures, probe GC-content, predicted operon structures, and
nitrogen content of gene products. Of the selected genes, 494 show a more than
2-fold difference in the two conditions analyzed. Under nitrogen-fixing
conditions 465 genes, mainly involved in energy metabolism, photosynthesis,
respiration and nitrogen-fixation, were found to be stronger expressed, whereas
29 genes showed a stronger expression under non-nitrogen fixing conditions.
Analysis of the nitrogen content of regulated genes shows that Nostoc PCC 7120
growing on dinitrogen is freed from any constraints to save nitrogen. For the
first time the expression of high light-induced stress proteins (HLIP-family)
is shown to be linked to the nitrogen availability.


  We present a novel method to investigate the dynamics of a single
semiflexible polymer, subject to anisotropic friction in a viscous fluid. In
contrast to previous approaches, we do not rely on a discrete bead-rod model,
but introduce a suitable normal mode decomposition of a continuous space curve.
By means of a perturbation expansion for stiff filaments we derive a closed set
of coupled Langevin equations in mode space for the nonlinear dynamics in two
dimensions, taking into account exactly the local constraint of
inextensibility. The stochastic differential equations obtained this way are
solved numerically, with parameters adjusted to describe the motion of actin
filaments. As an example, we show results for the tumbling motion in shear
flow.


  Reputed intractable, the question of the origin of viruses has long been
neglected. In the modern literature 'Virus evolution' has come to refer to
study more akin to population genetics, such as the world-wide scrutiny on new
polymorphisms appearing daily in the H5N1 avian flu virus [1], than to the
fundamental interrogation: where do viruses come from? This situation is now
rapidly changing, due to the coincidence of bold new ideas (and sometimes the
revival of old ones), the unexpected features exhibited by recently isolated
spectacular viruses [2] (see at URL: www.giantvirus.org), as well as the steady
increase of genomic sequences for 'regular' viruses and cellular organisms
enhancing the power of comparative genomics [3]. After being considered
non-living and relegated in the wings by a majority of biologists, viruses are
now pushed back on the center stage: they might have been at the origin of DNA,
of the eukaryotic cell, and even of today's partition of biological organisms
into 3 domains of life: bacteria, archaea and eukarya. Here, I quickly survey
some of the recent discoveries and the new evolutionary thoughts they have
prompted, before adding to the confusion with one interrogation of my own: what
if we totally missed the true nature of (at least some) viruses?


  We develop the qualitative theory of the solutions of the McKendrick partial
differential equation of population dynamics. We calculate explicitly the weak
solutions of the McKendrick equation and of the Lotka renewal integral equation
with time and age dependent birth rate. Mortality modulus is considered age
dependent. We show the existence of demography cycles. For a population with
only one reproductive age class, independently of the stability of the weak
solutions and after a transient time, the temporal evolution of the number of
individuals of a population is always modulated by a time periodic function.
The periodicity of the cycles is equal to the age of the reproductive age
class, and a population retains the memory from the initial data through the
amplitude of oscillations. For a population with a continuous distribution of
reproductive age classes, the amplitude of oscillation is damped. The
periodicity of the damped cycles is associated with the age of the first
reproductive age class. Damping increases as the dispersion of the fertility
function around the age class with maximal fertility increases. In general, the
period of the demography cycles is associated with the time that a species
takes to reach the reproductive maturity.


  I study the spreading of infectious diseases on heterogeneous populations. I
represent the population structure by a contact-graph where vertices represent
agents and edges represent disease transmission channels among them. The
population heterogeneity is taken into account by the agent's subdivision in
types and the mixing matrix among them. I introduce a type-network
representation for the mixing matrix allowing an intuitive understanding of the
mixing patterns and the analytical calculations. Using an iterative approach I
obtain recursive equations for the probability distribution of the outbreak
size as a function of time. I demonstrate that the expected outbreak size and
its progression in time are determined by the largest eigenvalue of the
reproductive number matrix and the characteristic distance between agents on
the contact-graph. Finally, I discuss the impact of intervention strategies to
halt epidemic outbreaks. This work provides both a qualitative understanding
and tools to obtain quantitative predictions for the spreading dynamics on
heterogeneous populations.


  We have developed an alignment-free method that calculates phylogenetic
distances using a maximum likelihood approach for a model of sequence change on
patterns that are discovered in unaligned sequences. To evaluate the
phylogenetic accuracy of our method, and to conduct a comprehensive comparison
of existing alignment-free methods (freely available as Python package decaf+py
at http://www.bioinformatics.org.au), we have created a dataset of reference
trees covering a wide range of phylogenetic distances. Amino acid sequences
were evolved along the trees and input to the tested methods; from their
calculated distances we infered trees whose topologies we compared to the
reference trees.
  We find our pattern-based method statistically superior to all other tested
alignment-free methods on this dataset. We also demonstrate the general
advantage of alignment-free methods over an approach based on automated
alignments when sequences violate the assumption of collinearity. Similarly, we
compare methods on empirical data from an existing alignment benchmark set that
we used to derive reference distances and trees. Our pattern-based approach
yields distances that show a linear relationship to reference distances over a
substantially longer range than other alignment-free methods. The pattern-based
approach outperforms alignment-free methods and its phylogenetic accuracy is
statistically indistinguishable from alignment-based distances.


  The exploration of the structural topology and the organizing principles of
genome-based large-scale metabolic networks is essential for studying possible
relations between structure and functionality of metabolic networks.
Topological analysis of graph models has often been applied to study the
structural characteristics of complex metabolic networks.In this work,
metabolic networks of 75 organisms were investigated from a topological point
of view. Network decomposition of three microbes (Escherichia coli, Aeropyrum
pernix and Saccharomyces cerevisiae) shows that almost all of the sub-networks
exhibit a highly modularized bow-tie topological pattern similar to that of the
global metabolic networks. Moreover, these small bow-ties are hierarchically
nested into larger ones and collectively integrated into a large metabolic
network, and important features of this modularity are not observed in the
random shuffled network. In addition, such a bow-tie pattern appears to be
present in certain chemically isolated functional modules and spatially
separated modules including carbohydrate metabolism, cytosol and mitochondrion
respectively. The highly modularized bow-tie pattern is present at different
levels and scales, and in different chemical and spatial modules of metabolic
networks, which is likely the result of the evolutionary process rather than a
random accident. Identification and analysis of such a pattern is helpful for
understanding the design principles and facilitate the modelling of metabolic
networks.


  As a discrete approach to genetic regulatory networks, Boolean models provide
an essential qualitative description of the structure of interactions among
genes and proteins. Boolean models generally assume only two possible states
(expressed or not expressed) for each gene or protein in the network as well as
a high level of synchronization among the various regulatory processes. In this
paper, we discuss and compare two possible methods of adapting qualitative
models to incorporate the continuous-time character of regulatory networks. The
first method consists of introducing asynchronous updates in the Boolean model.
In the second method, we adopt the approach introduced by L. Glass to obtain a
set of piecewise linear differential equations which continuously describe the
states of each gene or protein in the network. We apply both methods to a
particular example: a Boolean model of the segment polarity gene network of
Drosophila melanogaster. We analyze the dynamics of the model, and provide a
theoretical characterization of the model's gene pattern prediction as a
function of the timescales of the various processes.


  The extraction of information form high-throughput experiments is a key
aspect of modern biology. Early in the development of microarray technology,
researchers recognized that the size of the datasets and the limitations of
both computational and visualization techniques restricted their ability to
find the biological meaning hidden in the data. In addition, most researchers
wanted to make their datasets accessible to others. This resulted in the
development of new and advanced data storage, analysis, and visualization tools
enabling the cross-platform validation of the experiments and the
identification of previously undetected patterns. In order to reap the benefits
of this microarray data, researchers have needed to implement database
management systems providing integration of different experiments and data
types. Moreover, it was necessary to standardize the basic data structure and
experimental techniques for the standardization of microarray platforms. In
this chapter, we introduce the reader to the major concepts related to the use
of controlled vocabularies (ontologies), the definition of Minimum Information
About a Microarray Experiment (MIAME) and provide an overview of different
microarray data management strategies in use today. We summarize the main
characteristics of microarray data storage and sharing strategies including
warehouses, datamarts, and federations. The fundamental challenges involved in
the distribution, and retrieval of microarray data are presented, along with an
overview of some emerging technologies.


  In the past couple of years several studies have shown that hybridization in
Affymetrix DNA microarrays can be rather well understood on the basis of simple
models of physical chemistry. In the majority of the cases a Langmuir isotherm
was used to fit experimental data. Although there is a general consensus about
this approach, some discrepancies between different studies are evident. For
instance, some authors have fitted the hybridization affinities from the
microarray fluorescent intensities, while others used affinities obtained from
melting experiments in solution. The former approach yields fitted affinities
that at first sight are only partially consistent with solution values. In this
paper we show that this discrepancy exists only superficially: a sufficiently
complete model provides effective affinities which are fully consistent with
those fitted to experimental data. This link provides new insight on the
relevant processes underlying the functioning of DNA microarrays.


  We analyze publicly available data on Affymetrix microarrays spike-in
experiments on the human HGU133 chipset in which sequences are added in
solution at known concentrations. The spike-in set contains sequences of
bacterial, human and artificial origin. Our analysis is based on a recently
introduced molecular-based model [E. Carlon and T. Heim, Physica A 362, 433
(2006)] which takes into account both probe-target hybridization and
target-target partial hybridization in solution. The hybridization free
energies are obtained from the nearest-neighbor model with experimentally
determined parameters. The molecular-based model suggests a rescaling that
should result in a "collapse" of the data at different concentrations into a
single universal curve. We indeed find such a collapse, with the same
parameters as obtained before for the older HGU95 chip set. The quality of the
collapse varies according to the probe set considered. Artificial sequences,
chosen by Affymetrix to be as different as possible from any other human genome
sequence, generally show a much better collapse and thus a better agreement
with the model than all other sequences. This suggests that the observed
deviations from the predicted collapse are related to the choice of probes or
have a biological origin, rather than being a problem with the proposed model.


  Gene regulatory dynamics is governed by molecular processes and therefore
exhibits an inherent stochasticity. However, for the survival of an organism it
is a strict necessity that this intrinsic noise does not prevent robust
functioning of the system. It is still an open question how dynamical stability
is achieved in biological systems despite the omnipresent fluctuations. In this
paper we investigate the cell-cycle of the budding yeast Saccharomyces
cerevisiae as an example of a well-studied organism.
  We study a genetic network model of eleven genes that coordinate the
cell-cycle dynamics using a modeling framework which generalizes the concept of
discrete threshold dynamics. By allowing for fluctuations in the
transcription/translation times, we introduce noise in the model, accounting
for the effects of biochemical stochasticity. We study the dynamical attractor
of the cell cycle and find a remarkable robustness against fluctuations of this
kind. We identify mechanisms that ensure reliability in spite of fluctuations:
'Catcher' states and persistence of activity levels contribute significantly to
the stability of the yeast cell cycle despite the inherent stochasticity.


  Biological functions in living cells are controlled by protein interaction
and genetic networks. These molecular networks should be dynamically stable
against various fluctuations which are inevitable in the living world. In this
paper, we propose and study a stochastic model for the network regulating the
cell cycle of the budding yeast. The stochasticity in the model is controlled
by a temperature-like parameter $\beta$. Our simulation results show that both
the biological stationary state and the biological pathway are stable for a
wide range of "temperature". There is, however, a sharp transition-like
behavior at $\beta_c$, below which the dynamics is dominated by noise. We also
define a pseudo energy landscape for the system in which the biological pathway
can be seen as a deep valley.


  We propose a simple model for mass transport within a fungal hypha and its
subsequent growth. Inspired by the role of microtubule-transported vesicles, we
embody the internal dynamics of mass inside a hypha with mutually excluding
particles progressing stochastically along a growing one-dimensional lattice.
The connection between long range transport of materials for growth, and the
resulting extension of the hyphal tip has not previously been addressed in the
modelling literature. We derive and analyse mean-field equations for the model
and present a phase diagram of its steady state behaviour, which we compare to
simulations. We discuss our results in the context of the filamentous fungus,
Neurospora crassa.


  We have derived exact Langevin equations for a model of quasispecies
dynamics. The inherent multiplicative reaction noise is complex and its
statistical properties are specified completely. The numerical simulation of
the complex Langevin equations is carried out using the Cholesky decomposition
for the noise covariance matrix. This internal noise, which is due to
diffusion-limited reactions, produces unavoidable spatio-temporal density
fluctuations about the mean field value. In two dimensions, this noise strictly
vanishes only in the perfectly mixed limit, a situation difficult to attain in
practice.


  We investigate the application of hierarchical classification schemes to the
annotation of gene function based on several characteristics of protein
sequences including phylogenic descriptors, sequence based attributes, and
predicted secondary structure. We discuss three Bayesian models and compare
their performance in terms of predictive accuracy. These models are the
ordinary multinomial logit (MNL) model, a hierarchical model based on a set of
nested MNL models, and a MNL model with a prior that introduces correlations
between the parameters for classes that are nearby in the hierarchy. We also
provide a new scheme for combining different sources of information. We use
these models to predict the functional class of Open Reading Frames (ORFs) from
the E. coli genome. The results from all three models show substantial
improvement over previous methods, which were based on the C5 algorithm. The
MNL model using a prior based on the hierarchy outperforms both the
non-hierarchical MNL model and the nested MNL model. In contrast to previous
attempts at combining these sources of information, our approach results in a
higher accuracy rate when compared to models that use each data source alone.
Together, these results show that gene function can be predicted with higher
accuracy than previously achieved, using Bayesian models that incorporate
suitable prior information.


  A stochastic model for the activation of T-cells is analysed. T-cells are
part of the immune system and recognize foreign antigens against a background
of the body's own molecules. The model under consideration is a slight
generalization of a model introduced by Van den Berg, Rand and Burroughs in
2001, and is capable of explaining how this recognition works on the basis of
rare stochastic events. With the help of a refined large deviation theorem and
numerical evaluation it is shown that, for a wide range of parameters, T-cells
can distinguish reliably between foreign antigens and self-antigens.


  Using a theoretical model for spontaneous partial DNA unwrapping from
histones, we study the transient exposure of protein-binding DNA sites within
nucleosomes. We focus on the functional dependence of the rates for site
exposure and reburial on the site position, which is measurable experimentally
and pertinent to gene regulation. We find the dependence to be roughly
described by a random walker model. Close inspection reveals a surprising
physical effect of flexibility-assisted barrier crossing, which we characterize
within a toy model, the "semiflexible Brownian rotor."


  The Fast Fourier Transform (FFT) correlation approach to protein-protein
docking can evaluate the energies of billions of docked conformations on a grid
if the energy is described in the form of a correlation function. Here, this
restriction is removed, and the approach is efficiently used with pairwise
interactions potentials that substantially improve the docking results. The
basic idea is approximating the interaction matrix by its eigenvectors
corresponding to the few dominant eigenvalues, resulting in an energy
expression written as the sum of a few correlation functions, and solving the
problem by repeated FFT calculations. In addition to describing how the method
is implemented, we present a novel class of structure based pairwise
intermolecular potentials. The DARS (Decoys As the Reference State) potentials
are extracted from structures of protein-protein complexes and use large sets
of docked conformations as decoys to derive atom pair distributions in the
reference state. The current version of the DARS potential works well for
enzyme-inhibitor complexes. With the new FFT-based program, DARS provides much
better docking results than the earlier approaches, in many cases generating
50\% more near-native docked conformations. Although the potential is far from
optimal for antibody-antigen pairs, the results are still slightly better than
those given by an earlier FFT method. The docking program PIPER is freely
available for non-commercial applications.


  In the present article the recent works to formulate laws in Darwinian
evolutionary dynamics are discussed. Although there is a strong consensus that
general laws in biology may exist, opinions opposing such suggestion are
abundant. Based on recent progress in both mathematics and biology, another
attempt to address this issue is made in the present article. Specifically,
three laws which form a mathematical framework for the evolutionary dynamics in
biology are postulated. The second law is most quantitative and is explicitly
expressed in the unique form of a stochastic differential equation. Salient
features of Darwinian evolutionary dynamics are captured by this law: the
probabilistic nature of evolution, ascendancy, and the adaptive landscape. Four
dynamical elements are introduced in this formulation: the ascendant matrix,
the transverse matrix, the Wright evolutionary potential, and the stochastic
drive. The first law may be regarded as a special case of the second law. It
gives the reference point to discuss the evolutionary dynamics. The third law
describes the relationship between the focused level of description to its
lower and higher ones, and defines the dichotomy of deterministic and
stochastic drives. It is an acknowledgement of the hierarchical structure in
biology. A new interpretation of Fisher's fundamental theorem of natural
selection is provided in terms of the F-Theorem. The proposed laws are based on
continuous representation in both time and population. Their generic nature is
demonstrated through their equivalence to classical formulations. The present
three laws appear to provide a coherent framework for the further development
of the subject.


  The problem Parsimony Haplotyping (PH) asks for the smallest set of
haplotypes which can explain a given set of genotypes, and the problem Minimum
Perfect Phylogeny Haplotyping (MPPH) asks for the smallest such set which also
allows the haplotypes to be embedded in a perfect phylogeny, an evolutionary
tree with biologically-motivated restrictions. For PH, we extend recent work by
further mapping the interface between ``easy'' and ``hard'' instances, within
the framework of (k,l)-bounded instances where the number of 2's per column and
row of the input matrix is restricted. By exploring, in the same way, the
tractability frontier of MPPH we provide the first concrete, positive results
for this problem, and the algorithms underpinning these results offer new
insights about how MPPH might be further tackled in the future. In addition, we
construct for both PH and MPPH polynomial time approximation algorithms, based
on properties of the columns of the input matrix. We conclude with an overview
of intriguing open problems in PH and MPPH.


  Scientists investigate the dynamics of complex systems with quantitative
models, employing them to synthesize knowledge, to explain observations, and to
forecast future system behavior. Complete specification of systems is
impossible, so models must be simplified abstractions. Thus, the art of
modeling involves deciding which system elements to include and determining how
they should be represented. We view modeling as search through a space of
candidate models that is guided by model objectives, theoretical knowledge, and
empirical data. In this contribution, we introduce a method for representing
process-based models that facilitates the discovery of models that explain
observed behavior. This representation casts dynamic systems as interacting
sets of processes that act on entities. Using this approach, a modeler first
encodes relevant ecological knowledge into a library of generic entities and
processes, then instantiates these theoretical components, and finally
assembles candidate models from these elements. We illustrate this methodology
with a model of the Ross Sea ecosystem.


  Replacement of papers q-bio/0605026 and q-bio/0607043 by a single manuscript
(q-bio/0701029) follows suggestions from a magazine's referee


  We increasingly rely on the network approach to understand the complexity of
cellular functions. Chaperones (heat shock proteins) are key "networkers",
which have among their functions to sequester and repair damaged protein. In
order to link the network approach and chaperones with the aging process, we
first summarize the properties of aging networks suggesting a "weak link theory
of aging". This theory suggests that age-related random damage primarily
affects the overwhelming majority of the low affinity, transient interactions
(weak links) in cellular networks leading to increased noise, destabilization
and diversity. These processes may be further amplified by age-specific network
remodelling and by the sequestration of weakly linked cellular proteins to
protein aggregates of aging cells. Chaperones are weakly linked hubs [i.e.,
network elements with a large number of connections] and inter-modular bridge
elements of protein-protein interaction, signalling and mitochondrial networks.
As aging proceeds, the increased overload of damaged proteins is an especially
important element contributing to cellular disintegration and destabilization.
Additionally, chaperone overload may contribute to the increase of "noise" in
aging cells, which leads to an increased stochastic resonance resulting in a
deficient discrimination between signals and noise. Chaperone- and other
multi-target therapies, which restore the missing weak links in aging cellular
networks, may emerge as important anti-aging interventions.


  We investigated the structural relaxation of myosin motor domain from the
pre-power stroke state to the near-rigor state using molecular dynamics
simulation of a coarse-grained protein model. To describe the structural
change, we propose a "dual Go-model," a variant of the Go-like model that has
two reference structures. The nucleotide dissociation process is also studied
by introducing a coarse-grained nucleotide in the simulation. We found that the
myosin structural relaxation toward the near-rigor conformation cannot be
completed before the nucleotide dissociation. Moreover, the relaxation and the
dissociation occurred cooperatively when the nucleotide was tightly bound to
the myosin head. The result suggested that the primary role of the nucleotide
is to suppress the structural relaxation.


  Multivariate polynomial dynamical systems over finite fields have been
studied in several contexts, including engineering and mathematical biology. An
important problem is to construct models of such systems from a partial
specification of dynamic properties, e.g., from a collection of state
transition measurements. Here, we consider static models, which are directed
graphs that represent the causal relationships between system variables,
so-called wiring diagrams. This paper contains an algorithm which computes all
possible minimal wiring diagrams for a given set of state transition
measurements. The paper also contains several statistical measures for model
selection. The algorithm uses primary decomposition of monomial ideals as the
principal tool. An application to the reverse-engineering of a gene regulatory
network is included. The algorithm and the statistical measures are implemented
in Macaulay2 and are available from the authors.


  The observed structure of protein interaction networks is corrupted by many
false positive/negative links. This observational incompleteness is abstracted
as random link removal and a specific, experimentally motivated (spoke) link
rearrangement. Their impact on the structural properties of
gene-duplication-and-mutation network models is studied. For the degree
distribution a curve collapse is found, showing no sensitive dependence on the
link removal/rearrangement strengths and disallowing a quantitative extraction
of model parameters. The spoke link rearrangement process moves other
structural observables, like degree correlations, cluster coefficient and motif
frequencies, closer to their counterparts extracted from the yeast data. This
underlines the importance to take a precise modeling of the observational
incompleteness into account when network structure models are to be
quantitatively compared to data.


  Tree shape statistics quantify some aspect of the shape of a phylogenetic
tree. They are commonly used to compare reconstructed trees to evolutionary
models and to find evidence of tree reconstruction bias. Historically, to find
a useful tree shape statistic, formulas have been invented by hand and then
evaluated for utility. This article presents the first method which is capable
of optimizing over a class of tree shape statistics, called Binary Recursive
Tree Shape Statistics (BRTSS). After defining the BRTSS class, a set of
algebraic expressions is defined which can be used in the recursions. The tree
shape statistics definable using these expressions in the BRTSS is very
general, and includes many of the statistics with which phylogenetic
researchers are already familiar. We then present a practical genetic algorithm
which is capable of performing optimization over BRTSS given any objective
function. The chapter concludes with a successful application of the methods to
find a new statistic which indicates a significant difference between two
distributions on trees which were previously postulated to have similar
properties.


  A dynamical theory for the evolution of the genetic code is presented, which
accounts for its universality and optimality. The central concept is that a
variety of collective, but non-Darwinian, mechanisms likely to be present in
early communal life generically lead to refinement and selection of
innovation-sharing protocols, such as the genetic code. Our proposal is
illustrated using a simplified computer model, and placed within the context of
a sequence of transitions that early life may have made, prior to the emergence
of vertical descent.


  Self-similar properties of the ribosome in terms of the mass fractal
dimension are investigated. We find that both the 30S subunit and the 16S rRNA
have fractal dimensions of 2.58 and 2.82, respectively; while the 50S subunit
as well as the 23S rRNA has the mass fractal dimension close to 3, implying a
compact three dimensional macromolecule. This finding supports the dynamic and
active role of the 30S subunit in the protein synthesis, in contrast to the
pass role of the 50S subunit.


  We present an extremely simplified model of multiple-domains polymer
stretching in an atomic force microscopy experiment. We portray each module as
a binary set of contacts and decompose the system energy into a harmonic term
(the cantilever) and long-range interactions terms inside each domain. Exact
equilibrium computations and Monte Carlo simulations qualitatively reproduce
the experimental saw-tooth pattern of force-extension profiles, corresponding
(in our model) to first-order phase transitions. We study the influence of the
coupling induced by the cantilever and the pulling speed on the relative
heights of the force peaks. The results suggest that the increasing height of
the critical force for subsequent unfolding events is an out-of-equilibrium
effect due to a finite pulling speed. The dependence of the average unfolding
force on the pulling speed is shown to reproduce the experimental logarithmic
law.


  Our chances to halt epidemic outbreaks rely on how accurately we represent
the population structure underlying the disease spread. When analyzing global
epidemics this force us to consider metapopulation models taking into account
intra- and inter-community interactions. Recently Watts et al introduced a
metapopulation model which accounts for several features observed in real
outbreaks [Watts et al, PNAS 102, 11157 (2005)]. In this work I provide an
analytical solution to this model, enhancing our understanding of the model and
the epidemic outbreaks it represents. First, I demonstrate that depending on
the intra-community expected outbreak size and the fraction of social bridges
the epidemic outbreaks die out or there is a finite probability to observe a
global epidemics. Second, I show that the global scenario is characterized by
resurgent epidemics, their number increasing with increasing the
intra-community average distance between individuals. Finally, I present
empirical data for the AIDS epidemics supporting the model predictions.


  Naturally evolving proteins gradually accumulate mutations while continuing
to fold to thermodynamically stable native structures. This process of neutral
protein evolution is an important mode of genetic change, and forms the basis
for the molecular clock. Here we present a mathematical theory that predicts
the number of accumulated mutations, the index of dispersion, and the
distribution of stabilities in an evolving protein population from knowledge of
the stability effects (ddG values) for single mutations. Our theory
quantitatively describes how neutral evolution leads to marginally stable
proteins, and provides formulae for calculating how fluctuations in stability
cause an overdispersion of the molecular clock. It also shows that the
structural influences on the rate of sequence evolution that have been observed
in earlier simulations can be calculated using only the single-mutation ddG
values. We consider both the case when the product of the population size and
mutation rate is small and the case when this product is large, and show that
in the latter case proteins evolve excess mutational robustness that is
manifested by extra stability and increases the rate of sequence evolution. Our
basic method is to treat protein evolution as a Markov process constrained by a
minimal requirement for stable folding, enabling an evolutionary description of
the proteins solely in terms of the experimentally measureable ddG values. All
of our theoretical predictions are confirmed by simulations with model lattice
proteins. Our work provides a mathematical foundation for understanding how
protein biophysics helps shape the process of evolution.


  Cyclic dominance of species has been identified as a potential mechanism to
maintain biodiversity, see e.g. B. Kerr, M. A. Riley, M. W. Feldman and B. J.
M. Bohannan [Nature {\bf 418}, 171 (2002)] and B. Kirkup and M. A. Riley
[Nature {\bf 428}, 412 (2004)]. Through analytical methods supported by
numerical simulations, we address this issue by studying the properties of a
paradigmatic non-spatial three-species stochastic system, namely the
`rock-paper-scissors' or cyclic Lotka-Volterra model. While the deterministic
approach (rate equations) predicts the coexistence of the species resulting in
regular (yet neutrally stable) oscillations of the population densities, we
demonstrate that fluctuations arising in the system with a \emph{finite number
of agents} drastically alter this picture and are responsible for extinction:
After long enough time, two of the three species die out. As main findings we
provide analytic estimates and numerical computation of the extinction
probability at a given time. We also discuss the implications of our results
for a broad class of competing population systems.


  We investigate a social system of agents faced with a binary choice. We
assume there is a correct, or beneficial, outcome of this choice. Furthermore,
we assume agents are influenced by others in making their decision, and that
the agents can obtain information that may guide them towards making a correct
decision. The dynamic model we propose is of nonequilibrium type, converging to
a final decision. We run it on random graphs and scale-free networks. On random
graphs, we find two distinct regions in terms of the "finalizing time" -- the
time until all agents have finalized their decisions. On scale-free networks on
the other hand, there does not seem to be any such distinct scaling regions.


  In this work we propose a physical model of organismal evolution, where
phenotype, organism life expectancy, is directly related to genotype i.e. the
stability of its proteins which can be determined exactly in the model.
Simulating the model on a computer, we consistently observe the Big Bang
scenario whereby exponential population growth ensues as favorable
sequence-structure combinations (precursors of stable proteins) are discovered.
After that, random diversity of the structural space abruptly collapses into a
small set of preferred structural motifs. We observe that protein folds remain
stable and abundant in the population at time scales much greater than mutation
or organism lifetime, and the distribution of the lifetimes of dominant folds
in a population approximately follows a power law. The separation of
evolutionary time scales between discovery of new folds and generation of new
sequences gives rise to emergence of protein families and superfamilies whose
sizes are power-law distributed, closely matching the same distributions for
real proteins. The network of structural similarities of the universe of
evolved proteins has the same scale-free like character as the actual protein
domain universe graph (PDUG). Further, the model predicts that ancient protein
domains represent a highly connected and clustered subset of all protein
domains, in complete agreement with reality. Together, these results provide a
microscopic first principles picture of how protein structures and gene
families evolved in the course of evolution.


  We model the transcription factor based regulation network of yeast using a
content-based network model that mimicks the recognition of binding motifs on
the regulatory regions of the genes. We are thereby able to faithfully
reproduce many of the topological features of the gene regulatory network of
yeast once the parameters of the yeast genome, in particular the distribution
of information coded by the "binding sequences" within the promoter regions is
provided as input. The length distribution for the promoter regions is fixed by
comparing the k-core analysis of the model network with that of yeast. Our
results strongly point to the possibility that the observed topological
features are generic to networks formed via sequence-matching between random
strings obeying certain length distributions.


  This paper develops a set of simplified dynamical models with which to
explore the conditions under which division of labor leads to optimized system
output, as measured by the rate of production of a given product. We consider
two models: In the first model, we consider the flow of some resource into a
compartment, and the conversion of this resource into some product. In the
second model, we consider the resource-limited growth of autoreplicating
systems. In this case, we divide the replication and metabolic tasks among
different agents. The general features that emerge from our models is that
division of labor is favored when the resource to agent ratio is at
intermediate values, and when the time cost associated with transporting
intermediate products is small compared to characteristic process times. We
discuss the results of this paper in the context of simulations with digital
life. We also argue that division of labor in the context of our replication
model suggests an evolutionary basis for the emergence of the stem-cell-based
tissue architecture in complex organisms.


  Phi-values are experimental measures of the effects of mutations on the
folding kinetics of a protein. A central question is which structural
information Phi-values contain about the transition state of folding.
Traditionally, a Phi-value is interpreted as the 'nativeness' of a mutated
residue in the transition state. However, this interpretation is often
problematic because it assumes a linear relation between the nativeness of the
residue and its free-energy contribution. We present here a better structural
interpretation of Phi-values for mutations within a given helix. Our
interpretation is based on a simple physical model that distinguishes between
secondary and tertiary free-energy contributions of helical residues. From a
linear fit of our model to the experimental data, we obtain two structural
parameters: the extent of helix formation in the transition state, and the
nativeness of tertiary interactions in the transition state. We apply our model
to all proteins with well-characterized helices for which more than 10
Phi-values are available: protein A, CI2, and protein L. The model captures
nonclassical Phi-values <0 or >1 in these helices, and explains how different
mutations at a given site can lead to different Phi-values.


  In human cells, estrogenic signals induce cyclical association and
dissociation of specific proteins with the DNA in order to activate
transcription of estrogen-responsive genes. These oscillations can be modeled
by assuming a large number of sequential reactions represented by linear
kinetics with random kinetic rates. Application of the model to experimental
data predicts robust binding sequences in which proteins associate with the DNA
at several different phases of the oscillation. Our methods circumvent the need
to derive detailed kinetic graphs, and are applicable to other oscillatory
biological processes involving a large number of sequential steps.


  Clonal organisms present a particular challenge in population genetics
because, in addition to the possible existence of replicates of the same
genotype in a given sample, some of the hypotheses and concepts underlying
classical population genetics models are irreconcilable with clonality. The
genetic structure and diversity of clonal populations was examined using a
combination of new tools to analyze microsatellite data in the marine
angiosperm Posidonia oceanica. These tools were based on examination of the
frequency distribution of the genetic distance among ramets, termed the
spectrum of genetic diversity (GDS), and of networks built on the basis of
pairwise genetic distances among genets. The properties and topology of
networks based on genetic distances showed a "small-world" topology,
characterized by a high degree of connectivity among nodes, and a substantial
amount of substructure, revealing organization in sub-families of closely
related individuals.
  Keywords: genetic networks; small-world networks; genetic diversity; clonal
organisms


  Mathematical theory of selection systems is developed for a wide class of
dynamical models of inhomogeneous populations with discrete time. The Price
equation and its particular case, the Fisher Fundamental theorem of natural
selection (FTNS), are well known general results of the theory. It is known
that the Price equation being a mathematical identity is not dynamically
sufficient, i.e., it does not allow one to predict changes in the mean of a
trait beyond the immediate response if only the value of covariance of the
trait and fitness at this moment is known. We show that the problem of
dynamically insufficiency for the Price equations and for the FTNS can be
correctly overcome if to study these equations subject to given initial
distribution in framework of exact models of the population dynamics. The
knowledge of the entire distribution at any given instant allows making the
exact prediction for indefinite time and this prediction dramatically depends
on the initial distribution. For these models, the current trait distribution
and hence all statistical characteristics of interest, such as mean values of
the fitness or any trait could be computed effectively at any time moment.


  The hierarchical system of forest ecosystem models based on the theory of
individual-based (structured) models of populations and communities is briefly
described. New self-thinning models are integrated with tree stand models
within a hierarchical system of models, aimed at assessing (quasi-) stable
states of a forest ecosystem under given environment and management. Within the
structural model of succession, the ergodic theorems in biology are proved.
They assert that the areas of biocenoses in a climax state of association
should be proportional to the proper times of their development in a succession
line; the other characteristics of biocenoses should be also proportional to
their own proper regeneration times in a climax state. A simple theoretical
approach to estimating the deviation of the current state of forest ecosystem
from its steady state is derived and applied to data from the Prioksko-Terrasny
Reserve. The results obtained can serve as a theoretical basis in numerical
processing of data on ecological monitoring of undisturbed forest lands.


  Non-linear maps can possess various dynamical behaviors varying from stable
steady states and cycles to chaotic oscillations. Most models assume that
individuals within a given population are identical ignoring the fundamental
role of variation. Here we develop a theory of inhomogeneous maps and apply the
general approach to modeling heterogeneous populations with discrete
evolutionary time step. We show that the behavior of the inhomogeneous maps may
possess complex transition regimes, which depends both on the mean and the
variance of the initial parameter distribution. The examples of inhomogeneous
models are discussed.


  A model for the unidirectional movement of dynein is presented based on
structural observations and biochemical experimental results available. In this
model, the binding affinity of dynein for microtubule is independent of its
nucleotide state and the change between strong and weak microtubule-binding is
determined naturally by the variation of relative orientation between the stalk
and microtubule as the stalk rotates following nucleotide-state transition.
Thus the enigmatic communication from the ATP binding site in the globular
domain to the far MT-binding site in the tip of the stalk, which is
prerequisite in conventional models, is not required. Using the present model,
the previous experimental results such as the effect of ATP and ADP bindings on
dissociation of dynein from microtubule, the processive movement of
single-headed axonemal dyneins at saturating ATP concentration, the load
dependence of step size for the processive movement of two-headed cytoplasmic
dyneins and the dependence of stall force on ATP concentration can be well
explained.


  Microtubules (MTs) are important cytoskeletal structures, engaged in a number
of specific cellular activities, including vesicular traffic, cell
cyto-architecture and motility, cell division, and information processing
within neuronal processes. MTs have also been implicated in higher neuronal
functions, including memory, and the emergence of "consciousness". How MTs
handle and process electrical information, however, is heretofore unknown. Here
we show new electrodynamic properties of MTs. Isolated, taxol-stabilized
microtubules behave as bio-molecular transistors capable of amplifying
electrical information. Electrical amplification by MTs can lead to the
enhancement of dynamic information, and processivity in neurons can be
conceptualized as an "ionic-based" transistor, which may impact among other
known functions, neuronal computational capabilities.


  Alignment algorithms usually rely on simplified models of gaps for
computational efficiency. Based on an isomorphism between alignments and
physical helix-coil models, we show in statistical mechanics that alignments
with realistic laws for gaps can be computed with fast algorithms. Improved
performances of probabilistic alignments with realistic models of gaps are
illustrated. Probabilistic and optimization formulations are compared, with
potential implications in many fields and perspectives for computationally
efficient extensions to Markov models with realistic long-range interactions.


  An invitation to simulation is made into a realizable challenge: the reader
is taken from scratch to games with words, to mathematical models, to
statistical analysis, to simulation of evolution. No special requirements are
needed: we work with the inbuilt Visual Basic language that is attached to
Microsoft(R) Excel and Word. The Tex source must be downloaded to copy the code
that must be run.


  Over the last 10-15 years a general understanding of the chemical reaction of
protein folding has emerged from statistical mechanics. The lessons learned
from protein folding kinetics based on energy landscape ideas have benefited
protein structure prediction, in particular the development of coarse grained
models. We survey results from blind structure prediction. We explore how
second generation prediction energy functions can be developed by introducing
information from an ensemble of previously simulated structures. This procedure
relies on the assumption of a funnelled energy landscape keeping with the
principle of minimal frustration. First generation simulated structures provide
an improved input for associative memory energy functions in comparison to the
experimental protein structures chosen on the basis of sequence alignment.


  The precise details of how myosin-V coordinates the biochemical reactions and
mechanical motions of its two head elements to engineer effective processive
molecular motion along actin filaments remain unresolved. We compare a
quantitative kinetic model of the myosin-V walk, consisting of five basic
states augmented by two further states to allow for futile hydrolysis and
detachments, with experimental results for run lengths, velocities, and dwell
times and their dependence on bulk nucleotide concentrations and external loads
in both directions. The model reveals how myosin-V can use the internal strain
in the molecule to synchronise the motion of the head elements. Estimates for
the rate constants in the reaction cycle and the internal strain energy are
obtained by a computational comparison scheme involving an extensive
exploration of the large parameter space. This scheme exploits the fact that we
have obtained analytic results for our reaction network, e.g. for the velocity
but also the run length, diffusion constant and fraction of backward steps. The
agreement with experiment is often reasonable but some open problems are
highlighted, in particular the inability of such a general model to reproduce
the reported dependence of run length on ADP. The novel way that our approach
explores parameter space means that any confirmed discrepancies should give new
insights into the reaction network model.


  The prion protein (PrP) binds Cu2+ ions in the octarepeat domain of the
N-terminal tail up to full occupancy at pH=7.4. Recent experiments show that
the HGGG octarepeat subdomain is responsible for holding the metal bound in a
square planar coordination. By using first principle ab initio molecular
dynamics simulations of the Car-Parrinello type, the Cu coordination mode to
the binding sites of the PrP octarepeat region is investigated. Simulations are
carried out for a number of structured binding sites. Results for the complexes
Cu(HGGGW)+(wat), Cu(HGGG) and the 2[Cu(HGGG)] dimer are presented. While the
presence of a Trp residue and a H2O molecule does not seem to affect the nature
of the Cu coordination, high stability of the bond between Cu and the amide
Nitrogens of deprotonated Gly's is confirmed in the case of the Cu(HGGG)
system. For the more interesting 2[Cu(HGGG)] dimer a dynamically entangled
arrangement of the two monomers, with intertwined N-Cu bonds, emerges. This
observation is consistent with the highly packed structure seen in experiments
at full Cu occupancy.


  Human chromosome 21 is the only chromosome in human genome that exhibits
oscillation of (G+C)-content of cycle length of hundreds kilobases (500 kb near
the right telomere). We aim at establishing the existence of similar
periodicity in structure-related sequence features in order to relate this
(G+C)% oscillation to other biological phenomena. The following quantities are
shown to oscillate with the same 500kb periodicity in human chromosome 21:
binding energy calculated by two sets of dinucleotide-based thermodynamic
parameters, AA/TT and AAA/TTT bi-/tri-nucleotide density, 5'-TA-3' dinucleotide
density, and signal for 10/11-base periodicity of AA/TT or AAA/TTT. These
intrinsic quantities are related to structural features of the double helix of
DNA molecules, such as base-pair binding, untwisting/unwinding, stiffness, and
a putative tendency for nucleosome formation.


  We study coarse-grained (group-level) alignment dynamics of individual-based
animal group models for {\it heterogeneous} populations consisting of informed
(on preferred directions) and uninformed individuals. The orientation of each
individual is characterized by an angle, whose dynamics are nonlinearly coupled
with those of all the other individuals, with an explicit dependence on the
difference between the individual's orientation and the instantaneous average
direction. Choosing convenient coarse-grained variables (suggested by
uncertainty quantification methods) that account for rapidly developing
correlations during initial transients, we perform efficient computations of
coarse-grained steady states and their bifurcation analysis. We circumvent the
derivation of coarse-grained governing equations, following an equation-free
computational approach.


  In this paper, we propose a mean-field model which attempts to bridge the gap
between random Boolean networks and more realistic stochastic modeling of
genetic regulatory networks. The main idea of the model is to replace all
regulatory interactions to any one gene with an average or effective
interaction, which takes into account the repression and activation mechanisms.
We find that depending on the set of regulatory parameters, the model exhibits
rich nonlinear dynamics. The model also provides quantitative support to the
earlier qualitative results obtained for random Boolean networks.


  A new heuristic based on vertex invariants is developed to rapidly
distinguish non-isomorphic graphs to a desired level of accuracy. The method is
applied to sample subgraphs from an E.coli protein interaction network, and as
a probe for discovery of extended motifs. The network's structure is described
using statistical properties of its $N$-node subgraphs for $N\leq 14$. The Zipf
plots for subgraph occurrences are robust power laws that do not change when
rewiring the network while fixing the degree sequence -- although the specific
subgraphs may exchange ranks. However the exponent depends on $N$. The study of
larger subgraphs highlights some striking patterns for various $N$. Motifs, or
connected pieces that are over-abundant in the ensemble of subgraphs, have more
edges, for a given number of nodes, than antimotifs and generally display a
bipartite structure or tend towards a complete graph. In contrast, antimotifs,
which are under-abundant connected pieces, are mostly trees or contain at most
a single, small loop. The extension to directed graphs is straightforward.


  We study the adaptation dynamics of an initially maladapted asexual
population with genotypes represented by binary sequences of length $L$. The
population evolves in a maximally rugged fitness landscape with a large number
of local optima. We find that whether the evolutionary trajectory is
deterministic or stochastic depends on the effective mutational distance
$d_{\mathrm{eff}}$ upto which the population can spread in genotype space. For
$d_{\mathrm{eff}}=L$, the deterministic quasispecies theory operates while for
$d_{\mathrm{eff}} < 1$, the evolution is completely stochastic. Between these
two limiting cases, the dynamics are described by a local quasispecies theory
below a crossover time $T_{\times}$ while above $T_{\times}$, the population
gets trapped at a local fitness peak and manages to find a better peak either
via stochastic tunneling or double mutations. In the stochastic regime
$d_\mathrm{eff} < 1$, we identify two subregimes associated with clonal
interference and uphill adaptive walks, respectively. We argue that our
findings are relevant to the interepretation of evolution experiments with
microbial populations.


  The hybrid form is a combination of the Rydberg potential and the London
inverse-sixth-power energy. It is accurate at all relevant distance scales and
simple enough for use in all-atom simulations of biomolecules. One may compute
the parameters of the hybrid potential for the ground state of a pair of
neutral atoms from their internuclear separation, the depth and curvature of
their potential at its minimum, and from their van der Waals coefficient of
dispersion.


  Undulatory locomotion is common to nematodes as well as to limbless
vertebrates, but its control is not understood in spite of the identification
of hundred of genes involved in Caenorhabditis elegans locomotion. To reveal
the mechanisms of nematode undulatory locomotion, we quantitatively analyzed
the movement of C. elegans with genetic perturbations to neurons, muscles, and
skeleton (cuticle). We also compared locomotion of different Caenorhabditis
species. We constructed a theoretical model that combines mechanics and
biophysics, and that is constrained by the observations of propulsion and
muscular velocities, as well as wavelength and amplitude of undulations. We
find that normalized wavelength is a conserved quantity among wild-type C.
elegans individuals, across mutants, and across different species. The velocity
of forward propulsion scales linearly with the velocity of the muscular wave
and the corresponding slope is also a conserved quantity and almost optimal;
the exceptions are in some mutants affecting cuticle structure. In theoretical
terms, the optimality of the slope is equivalent to the exact balance between
muscular and visco-elastic body reaction bending moments. We find that the
amplitude and frequency of undulations are inversely correlated and provide a
theoretical explanation for this fact. These experimental results are valid
both for young adults and for all larval stages of wild type C. elegans. In
particular, during development, the amplitude scales linearly with the
wavelength, consistent with our theory. We also investigated the influence of
substrate firmness on motion parameters, and found that it does not affect the
above invariants. In general, our biomechanical model can explain the observed
robustness of the mechanisms controlling nematode undulatory locomotion.


  In this work we analyze the possibility that soliton dynamics in a simple
nonlinear model allows functionally relevant predictions of the behaviour of
DNA. This suggestion was first put forward by Salerno [Phys. Rev. A, vol. 44,
p. 5292 (1991)] by showing results indicating that sine-Gordon kinks were set
in motion at certain regions of a DNA sequence that include promoters. We
revisit that system and show that the observed behaviour has nothing to do with
promoters; on the contrary, it originates from the bases at the boundary, which
are not part of the studied genome. We explain this phenomenology in terms of
an effective potential for the kink center. This is further extended to
disprove recent claims that the dynamics of kinks [Lenholm and H\"ornquist,
Physica D, vol. 177, p. 233 (2003)] or breathers [Bashford, J. Biol. Phys.,
vol. 32, p. 27 (2006)] has functional significance. We conclude that no such
information can be extracted from this simple nonlinear model or its associated
effective potential.


  In protein-protein interaction networks certain topological properties appear
to be recurrent: networks maps are considered scale-free. It is possible that
this topology is reflected in the protein structure. In this paper we
investigate the role of protein disorder in the network topology. We find that
the disorder of a protein (or of its neighbors) is independent of its number of
protein-protein interactions. This result suggests that protein disorder does
not play a role in the scale-free architecture of protein networks.


  The biochemical changes on normal and degenerated tissues of rotator cuff
supraspinatus tendons were probed by FT-Raman spectroscopy. The Raman spectra
showed differences on the spectral regions of cysteine, amino acids, nucleic
acids, carbohydrates, and lipids. These spectral differences were assigned to
pathological biochemical alterations due to the degenerative process of the
tendon. Principal Components Analysis was performed on the spectral data and
enabled the correct classification of the spectra as normal (grade 1) and
degenerated (grades 2 and 3). These findings indicate that Raman spectroscopy
could be a very promising tool for the rotator cuff supraspinatus tendon
diagnosis and for quantification of their degenerative degree.


  We formulate a simple solvation potential based on a coarsed-grain
representation of amino acids with two spheres modeling the $C_\alpha$ atom and
an effective side-chain centroid. The potential relies on a new method for
estimating the buried area of residues, based on counting the effective number
of burying neighbours in a suitable way. This latter quantity shows a good
correlation with the buried area of residues computed from all atom
crystallographic structures. We check the discriminatory power of the solvation
potential alone to identify the native fold of a protein from a set of decoys
and show the potential to be considerably selective.


  We describe here the new concept of $\epsilon$-Homomorphisms of Probabilistic
Regulatory Gene Networks(PRN). The $\epsilon$-homomorphisms are special
mappings between two probabilistic networks, that consider the algebraic action
of the iteration of functions and the probabilistic dynamic of the two
networks. It is proved here that the class of PRN, together with the
homomorphisms, form a category with products and coproducts. Projections are
special homomorphisms, induced by invariant subnetworks. Here, it is proved
that an $\epsilon$-homomorphism for 0 <$\epsilon$< 1 produces simultaneous
Markov Chains in both networks, that permit to introduce the concepts of
$\epsilon$-isomorphism of Markov Chains, and similar networks.


  Evolutionary game theory has traditionally assumed that all individuals in a
population interact with each other between reproduction events. We show that
eliminating this restriction by explicitly considering the time scales of
interaction and selection leads to dramatic changes in the outcome of
evolution. Examples include the selection of the inefficient strategy in the
Harmony and Stag-Hunt games, and the disappearance of the coexistence state in
the Snowdrift game. Our results hold for any population size and in the
presence of a background of fitness.


  We have used the sexual Penna ageing model to show that the relation between
dominance and recessiveness could be a force which optimizes the genome size.
While the possibility of complementation of the damaged allele by its
functional counterparts (recessiveness) leads to the redundancy of genetic
information, the dominant effect of defective genes tends to diminish the
number of alleles fulfilling the same function. Playing with the fraction of
dominant loci in the genome it is possible to obtain the condition where the
diploid state of the genome is optimal. If the status of each bit position as
dominant or recessive mutations is changed for each individual randomly and
rarely, then after a long time a stationary equilibrium of many recessive and
few dominant loci is established in the sexual Penna model. This effect
vanishes if the same changing distribution of dominant loci applies to all
individuals.


  Background: Radiotherapy outcomes are usually predicted using the Linear
Quadratic model. However, this model does not integrate complex features of
tumor growth, in particular cell cycle regulation.
  Methods: In this paper, we propose a multiscale model of cancer growth based
on the genetic and molecular features of the evolution of colorectal cancer.
The model includes key genes, cellular kinetics, tissue dynamics, macroscopic
tumor evolution and radiosensitivity dependence on the cell cycle phase. We
investigate the role of gene-dependent cell cycle regulation in the response of
tumors to therapeutic irradiation protocols.
  Results: Simulation results emphasize the importance of tumor tissue features
and the need to consider regulating factors such as hypoxia, as well as tumor
geometry and tissue dynamics, in predicting and improving radiotherapeutic
efficacy.
  Conclusion: This model provides insight into the coupling of complex
biological processes, which leads to a better understanding of oncogenesis.
This will hopefully lead to improved irradiation therapy.


  Successive whole genome duplications have recently been firmly established in
all major eukaryote kingdoms. It is not clear, however, how such dramatic
evolutionary process has contributed to shape the large scale topology of
protein-protein interaction (PPI) networks. We propose and analytically solve a
generic model of PPI network evolution under successive whole genome
duplications. This demonstrates that the observed scale-free degree
distributions and conserved multi-protein complexes may have concomitantly
arised from i) intrinsic exponential dynamics of PPI network evolution and ii)
asymmetric divergence of gene duplicates. This requirement of asymmetric
divergence is in fact "spontaneously" fulfilled at the level of protein-binding
domains. In addition, domain shuffling of multi-domain proteins is shown to
provide a powerful combinatorial source of PPI network innovation, while
preserving essential structures of the underlying single-domain interaction
network. Finally, large scale features of PPI networks reflecting the
"combinatorial logic" behind direct and indirect protein interactions are well
reproduced numerically with only two adjusted parameters of clear biological
significance.


  We study fixation probabilities and times as a consequence of neutral genetic
drift in subdivided populations, motivated by a model of the cultural
evolutionary process of language change that is described by the same
mathematics as the biological process. We focus on the growth of fixation times
with the number of subpopulations, and variation of fixation probabilities and
times with initial distributions of mutants. A general formula for the fixation
probability for arbitrary initial condition is derived by extending a duality
relation between forwards- and backwards-time properties of the model from a
panmictic to a subdivided population. From this we obtain new formulae,
formally exact in the limit of extremely weak migration, for the mean fixation
time from an arbitrary initial condition for Wright's island model, presenting
two cases as examples. For more general models of population subdivision,
formulae are introduced for an arbitrary number of mutants that are randomly
located, and a single mutant whose position is known. These formulae contain
parameters that typically have to be obtained numerically, a procedure we
follow for two contrasting clustered models. These data suggest that variation
of fixation time with the initial condition is slight, but depends strongly on
the nature of subdivision. In particular, we demonstrate conditions under which
the fixation time remains finite even in the limit of an infinite number of
demes. In many cases - except this last where fixation in a finite time is seen
- the time to fixation is shown to be in precise agreement with predictions
from formulae for the asymptotic effective population size.


  Having in mind the large-scale analysis of gene regulatory networks, we
review a graph decimation algorithm, called "leaf-removal", which can be used
to evaluate the feedback in a random graph ensemble. In doing this, we consider
the possibility of analyzing networks where the diagonal of the adjacency
matrix is structured, that is, has a fixed number of nonzero entries. We test
these ideas on a network model with fixed degree, using both numerical and
analytical calculations. Our results are the following. First, the leaf-removal
behavior for large system size enables to distinguish between different regimes
of feedback. We show their relations and the connection with the onset of
complexity in the graph. Second, the influence of the diagonal structure on
this behavior can be relevant.


  Interactions in ecological communities are inherently nonlinear and can lead
to complex population dynamics including irregular fluctuations induced by
chaos. Chaotic population dynamics can exhibit violent oscillations with
extremely small or large population abundances that might cause extinction and
recurrent outbreaks, respectively. We present a simple method that can guide
management efforts to prevent crashes, peaks, or any other undesirable state.
At the same time, the irregularity of the dynamics can be preserved when chaos
is desirable for the population. The control scheme is easy to implement
because it relies on time series information only. The method is illustrated by
two examples: control of crashes in the Ricker map and control of outbreaks in
a stage-structured model of the flour beetle Tribolium. It turns out to be
effective even with few available data and in the presence of noise, as is
typical for ecological settings.


  FT-Raman spectroscopy was employed to study (in-vitro) non-altered human
colorectal tissues aiming evaluate their spectral differences. The samples were
collected from 39 patients, adding 144 spectra. The results enable one
estabilish 3 well defined spectroscopic groups of what was consistently checked
by statistical (clustering) and biological (histopathology) analysis. The
similarity within each group was better than 90%. Groups 1 and 2 had 80% of
similarity while presenting both 70% of similarity to group 3. All groups
presented connective and epithelial tissues as histopathological
characteristic. Group 1 differs from others by the presence of soft muscle
while group 3 by presence of fatty tissue. Tissues from group 2 are composed by
samples with connective and epithelial tissues. This study is very relevant in
order to establish the non-altered or normal spectroscopic standard for future
studies and applications on optical biopsy and diagnosis of colorectal cancer.


  The aim of this work is to elucidate how physical principles of protein
design are reflected in natural sequences that evolved in response to the
thermal conditions of the environment. Using an exactly solvable lattice model,
we design sequences with selected thermal properties. Compositional analysis of
designed model sequences and natural proteomes reveals a specific trend in
amino acid compositions in response to the requirement of stability at elevated
environmental temperature, i.e. the increase of fractions of hydrophobic and
charged amino acid residues at the expense of polar ones. We show that this
from both ends of hydrophobicity scale trend is due to positive (to stabilize
the native state) and negative (to destabilize misfolded states) components of
protein design. Negative design strengthens specific repulsive nonnative
interactions that appear in misfolded structures. A pressure to preserve
specific repulsive interactions in non-native conformations may result in
correlated mutations between amino acids which are far apart in the native
state but may be in contact in misfolded conformations. Such correlated
mutations are indeed found in TIM barrel and other proteins.


  Prokaryotes living at extreme environmental temperatures exhibit pronounced
signatures in the amino acid composition of their proteins and nucleotide
compositions of their genomes reflective of adaptation to their thermal
environments. However, despite significant efforts, the definitive answer of
what are the genomic and proteomic compositional determinants of Optimal Growth
Temperature of prokaryotic organisms remained elusive. Here the authors
performed a comprehensive analysis of amino acid and nucleotide compositional
signatures of thermophylic adaptation by exhaustively evaluating all
combinations of amino acids and nucleotides as possible determinants of Optimal
Growth Temperature for all prokaryotic organisms with fully sequences genomes..
The authors discovered that total concentration of seven amino acids in
proteomes, IVYWREL, serves as a universal proteomic predictor of Optimal Growth
Temperature in prokaryotes. Resolving the old-standing controversy the authors
determined that the variation in nucleotide composition (increase of purine
load, or A+G content with temperature) is largely a consequence of thermal
adaptation of proteins. However, the frequency with which A and G nucleotides
appear as nearest neighbors in genome sequences is strongly and independently
correlated with Optimal Growth Temperature. as a result of codon bias in
corresponding genomes. Together these results provide a complete picture of
proteomic and genomic determinants of thermophilic adaptation.


  In this paper we consider two populations whose generations are not
overlapping and whose size is large. The number of males and females in both
populations is constant. Any generation is replaced by a new one and any
individual has two parents for what concerns nuclear DNA and a single one (the
mother) for what concerns mtDNA. Moreover, at any generation some individuals
migrate from the first population to the second.
  In a finite random time $T$, the mtDNA of the second population is completely
replaced by the mtDNA of the first. In the same time, the nuclear DNA is not
completely replaced and a fraction $F$ of the ancient nuclear DNA persists. We
compute both $T$ and $F$. Since this study shows that complete replacement of
mtDNA in a population is compatible with the persistence of a large fraction of
nuclear DNA, it may have some relevance for the Out of Africa/Multiregional
debate in Paleoanthropology.


  This paper develops a set of simplified dynamical models with which to
explore the conditions under which temporal differentiation leads to optimized
system output. By temporal differentiation, we mean a division of labor whereby
different subtasks associated with performing a given task are done at
different times. The idea is that, by focusing on one particular set of
subtasks at a time, it is possible to increase the efficiency with which each
subtask is performed, thereby allowing for faster completion of the overall
task. For this paper, we consider a process whereby some resource is converted
into some final product in a series of three agent-mediated steps. Temporal
differentiation is incorporated by allowing the agents to oscillate between
performing the first two steps and performing the last step. We find that
temporal differentiation is favored when the number of agents is small, and
when the process intermediates have a much longer lifetime than the original
resource. Within the framework of biological systems, we argue that these
results provide a possible evolutionary basis for the emergence of sleep, and
also of distinct REM and non-REM sleep states. We also discuss our use of a
three-step model. Briefly, in order for temporal differentiation to increase
product output in a mean-field description of resource metabolism, it is
necessary for temporal differentiation to have a nonlinear effect on individual
process rates. For stochastic models, we argue that temporal differentiation
can increase product output even in fundamentally linear systems.


  The distribution of most genes is not random, and functionally linked genes
are often found in clusters. Several theories have been put forward to explain
the emergence and persistence of operons in bacteria. Careful analysis of
genomic data favours the co-regulation model, where gene organization into
operons is driven by the benefits of coordinated gene expression and
regulation. Direct evidence that co-expression increases the individual's
fitness enough to ensure operon formation and maintenance is, however, still
lacking. Here, a previously described quantitative model of the network that
controls the transcription factor sigmaF during sporulation in Bacillus
subtilis is employed to quantify the benefits arising from both organisation of
the sporulation genes into the spoIIA operon and from translational coupling.
The analysis shows that operon organization, together with translational
coupling, is important because of the inherent stochastic nature of gene
expression which skews the ratios between protein concentrations in the absence
of co- regulation. The predicted impact of different forms of gene regulation
on fitness and survival agrees quantitatively with published sporulation
efficiencies.


  The distribution of bubble lengths in double-stranded DNA is presented for
segments of varying guanine-cytosine (GC) content, obtained with Monte Carlo
simulations using the Peyrard-Bishop-Dauxois model at 310 K. An analytical
description of the obtained distribution in the whole regime investigated,
i.e., up to bubble widths of the order of tens of nanometers, is available. We
find that the decay lengths and characteristic exponents of this distribution
show two distinct regimes as a function of GC content. The observed
distribution is attributed to the anharmonic interactions within base pairs.
The results are discussed in the framework of the Poland-Scheraga and the
Peyrard-Bishop (with linear instead of nonlinear stacking interaction) models.


  Our work stems from the consideration that the spreading of a disease is
modulated by the individual's perception of the infected neighborhood and
his/her strategy to avoid being infected as well. We introduced a general
``cellular agent'' model that accounts for a hetereogeneous and variable
network of connections. The probability of infection is assumed to depend on
the perception that an individual has about the spreading of the disease in her
local neighborhood and on broadcasting media. In the one-dimensional
homogeneous case the model reduces to the DK one, while for long-range coupling
the dynamics exhibits large fluctuations that may lead to the complete
extinction of the disease.


  Statistical mechanics of relative species abundance (RSA) patterns in
biological networks is presented. The theory is based on multispecies
replicator dynamics equivalent to the Lotka-Volterra equation, with diverse
interspecies interactions. Various RSA patterns observed in nature are derived
from a single parameter related to productivity or maturity of a community. The
abundance distribution is formed like a widely observed left-skewed lognormal
distribution. It is also found that the ``canonical hypothesis'' is supported
in some parameter region where the typical RSA patterns are observed. As the
model has a general form, the result can be applied to similar patterns in
other complex biological networks, e.g. gene expression.


  We express the Crow-Kimura and Eigen models of quasispecies theory in a
functional integral representation. We formulate the spin coherent state
functional integrals using the Schwinger Boson method. In this formulation, we
are able to deduce the long-time behavior of these models for arbitrary
replication and degradation functions.
  We discuss the phase transitions that occur in these models as a function of
mutation rate. We derive for these models the leading order corrections to the
infinite genome length limit.


  No abstract given. Confirms earlier simulatiobns of the self-organization of
dominance in the sexual Penna model, and the advantage of hermaphroditic over
sexual reproduction.


  Toxic or allelopathic compounds liberated by toxin-producing phytoplankton
(TPP) acts as a strong mediator in plankton dynamics. On an analysis of a set
of phytoplankton biomass-data that have been collected by our group in the
North-West part of the Bay of Bengal, and by analysis of a three-component
mathematical model under a constant as well as a stochastic environment, we
explore the role of toxin-allelopathy in determining the dynamic behaviour of
the competing-phytoplankton species. The overall results, based on analytical
and numerical wings, demonstrate that toxin-allelopathy due to the
toxin-producing phytoplankton (TPP) promotes a stable coexistence of those
competitive phytoplankton that would otherwise exhibit competitive exclusion of
the weak species. Our study suggests that TPP might be a potential candidate
for maintaining the coexistence and diversity of competing phytoplankton
species.


  The manner in which groups of neurons represent events in the external world
is fundamental to neuroscience. Here, we analyze the population code of the
retina during naturalistic stimulation and show that the information conveyed
by specific multi-neuronal firing patterns can be very different from the sum
of the pattern's parts. Synchronous spikes convey more information than either
of the participating cells, but almost always less than their sum, making them
redundant coding symbols. Surprisingly, patterns of spiking and silence are
mostly synergistic - carrying information that is accessible only by observing
the whole pattern of activity, rather than its components, and signifying
unique features in the stimulus. These results demonstrate that the retina uses
a combinatorial code and that the brain can benefit significantly from
recognizing multi-neuronal firing patterns.


  Using basic properties of p-adic numbers, we consider a simple new approach
to describe main aspects of DNA sequence and genetic code. Central role in our
investigation plays an ultrametric p-adic information space which basic
elements are nucleotides, codons and genes. We show that a 5-adic model is
appropriate for DNA sequence. This 5-adic model, combined with 2-adic distance,
is also suitable for genetic code and for a more advanced employment in
genomics. We find that genetic code degeneracy is related to the p-adic
distance between codons.


  Chemical reaction systems operating in nonequilibrium open-system states
arise in a great number of contexts, including the study of living organisms,
in which chemical reactions, in general, are far from equilibrium. Here we
introduce a theorem that relates forward and re-verse fluxes and free energy
for any chemical process operating in a steady state. This rela-tionship, which
is a generalization of equilibrium conditions to the case of a chemical process
occurring in a nonequilibrium steady state, provides a novel equivalent
definition for chemical reaction free energy. In addition, it is shown that
previously unrelated theories introduced by Ussing and Hodgkin and Huxley for
transport of ions across membranes, Hill for catalytic cycle fluxes, and Crooks
for entropy production in microscopically reversible systems, are united in a
common framework based on this relationship.


  In a recent paper [Phys. Rev. E 68, 011906 (2003)], Naef and Magnasco
suggested that the "bright" mismatches observed in Affymetrix microarray
experiments are caused by the fluorescent molecules used to label RNA target
sequences, which would impede target-probe hybridization. Their conclusion is
based on the observation of "unexpected" asymmetries in the affinities obtained
by fitting microarray data from publicly available experiments. We point out
here that the observed asymmetry is due to the inequivalence of RNA and DNA,
and that the reported affinities are consistent with stacking free energies
obtained from melting experiments of unlabeled nucleic acids in solution. The
conclusion of Naef and Magnasco is therefore based on an unjustified
assumption.


  We present the methods and results of a two-stage modeling process that
generates candidate gene-regulatory networks of the bacterium B. subtilis from
experimentally obtained, yet mathematically underdetermined microchip array
data. By employing a computational, linear correlative procedure to generate
these networks, and by analyzing the networks from a graph theoretical
perspective, we are able to verify the biological viability of our inferred
networks, and we demonstrate that our networks' graph theoretical properties
are remarkably similar to those of other biological systems. In addition, by
comparing our inferred networks to those of a previous, noisier implementation
of the linear inference process [17], we are able to identify trends in graph
theoretical behavior that occur both in our networks as well as in their
perturbed counterparts. These commonalities in behavior at multiple levels of
complexity allow us to ascertain the level of complexity to which our process
is robust to noise.


  Cellular signaling networks have evolved to cope with intrinsic fluctuations,
coming from the small numbers of constituents, and the environmental noise.
Stochastic chemical kinetics equations govern the way biochemical networks
process noisy signals. The essential difficulty associated with the master
equation approach to solving the stochastic chemical kinetics problem is the
enormous number of ordinary differential equations involved. In this work, we
show how to achieve tremendous reduction in the dimensionality of specific
reaction cascade dynamics by solving variationally an equivalent quantum field
theoretic formulation of stochastic chemical kinetics. The present formulation
avoids cumbersome commutator computations in the derivation of evolution
equations, making more transparent the physical significance of the variational
method. We propose novel time-dependent basis functions which work well over a
wide range of rate parameters. We apply the new basis functions to describe
stochastic signaling in several enzymatic cascades and compare the results so
obtained with those from alternative solution techniques. The variational
ansatz gives probability distributions that agree well with the exact ones,
even when fluctuations are large and discreteness and nonlinearity are
important. A numerical implementation of our technique is many orders of
magnitude more efficient computationally compared with the traditional Monte
Carlo simulation algorithms or the Langevin simulations.


  The increasing availability of high throughput data arising from gene
expression studies leads to the necessity of methods for summarizing the
available information. As annotation quality improves it is becoming common to
rely on the Gene Ontology (GO) to build functional profiles that characterize a
set of genes using the frequency of use of each GO term or group of terms in
the array. In this work we describe a statistical model for such profiles,
provide methods to compare profiles and develop inferential procedures to
assess this comparison. An R-package implementing the methods is available.


  The building of minimal self-reproducing systems with a physical embodiment
(generically called protocells) is a great challenge, with implications for
both theory and applied sciences. Although the classical view of a living
protocell assumes that it includes information-carrying molecules as an
essential ingredient, a dividing cell-like structure can be built from a
metabolism-container coupled system, only. An example of such a system, modeled
with dissipative particle dynamics, is presented here. This article
demonstrates how a simple coupling between a precursor molecule and surfactant
molecules forming micelles can experience a growth-division cycle in a
predictable manner, and analyzes the influence of crucial parameters on this
replication cycle. Implications of these results for origins of cellular life
and living technology are outlined.


  We used various analytical and numerical techniques to elucidate signal
propagation in a small enzymatic cascade which is subjected to external and
internal noise. The nonlinear character of catalytic reactions, which underlie
protein signal transduction cascades, renders stochastic signaling dynamics in
cytosol biochemical networks distinct from the usual description of stochastic
dynamics in gene regulatory networks. For a simple 2-step enzymatic cascade
which underlies many important protein signaling pathways, we demonstrated that
the commonly used techniques such as the linear noise approximation and the
Langevin equation become inadequate when the number of proteins becomes too
low. Consequently, we developed a new analytical approximation, based on mixing
the generating function and distribution function approaches, to the solution
of the master equation that describes nonlinear chemical signaling kinetics for
this important class of biochemical reactions. Our techniques work in a much
wider range of protein number fluctuations than the methods used previously. We
found that under certain conditions the burst-phase noise may be injected into
the downstream signaling network dynamics, resulting possibly in unusually
large macroscopic fluctuations. In addition to computing first and second
moments, which is the goal of commonly used analytical techniques, our new
approach provides the full time-dependent probability distributions of the
colored non-Gaussian processes in a nonlinear signal transduction cascade.


  We show, that the specific distribution of gene's length, which is observed
in natural genomes, might be a result of a growth process, in which a single
length scale $L(t)$ develops that grows with time as $t^{1/3}$. This length
scale could be associated with the length of the longest gene in an evolving
genome. The growth kinetics of the genes resembles the one observed in physical
systems with conserved ordered parameter. We show, that in genome this
conservation is guaranteed by compositional compensation along DNA strands of
the purine-like trends introduced by genes. The presented mathematical model is
the modified Bak-Sneppen model of critical self-organization applied to the
one-dimensional system of $N$ spins. The spins take discrete values, which
represent gene's length.


  We present a novel topological classification of RNA secondary structures
with pseudoknots. It is based on the topological genus of the circular diagram
associated to the RNA base-pair structure. The genus is a positive integer
number, whose value quantifies the topological complexity of the folded RNA
structure. In such a representation, planar diagrams correspond to pure RNA
secondary structures and have zero genus, whereas non planar diagrams
correspond to pseudoknotted structures and have higher genus. We analyze real
RNA structures from the databases wwPDB and Pseudobase, and classify them
according to their topological genus. We compare the results of our statistical
survey with existing theoretical and numerical models. We also discuss possible
applications of this classification and show how it can be used for identifying
new RNA structural motifs.


  In this work we suggest and consider an original, simple mathematical model
of a "quasi-rapid" extinction population dynamics. It describes a decrease and
final extinction of the population of one prey species by a "quasi-rapid"
interaction with one predator species with increasing population. This
"quasi-rapid" interaction means ecologically that prey species behaves
practically quite passively (since there is no time for any reaction, i.e.
defense), like an appropriate environment, in respect to "quasi-rapid" activity
of the predator species that can have different "quasi-rapid" hunting
abilities. Mathematically, our model is based on a non-Lotka-Volterraian system
of two differential equations of the first order, first of which is linear
while second, depending of a parameter that characterizes hunting ability is
nonlinear. We compare suggested "quasi-rapid" extinction population dynamics
and the global model of the overkill of the prehistoric megafauna (mammoths).
We demonstrate that our "quasi-rapid" extinction population dynamics is able to
restitute successfully correlations between empirical (archeological) data and
overkill theory in North America as well as Australia. For this reason, we
conclude that global overkill theory, completely mathematically modelable by
"quasi-rapid" extinction population dynamics can consistently explain the
Pleistocene extinctions of the megafauna.


  Although atomistic simulations of proteins and other biological systems are
approaching microsecond timescales, the quality of trajectories has remained
difficult to assess. Such assessment is critical not only for establishing the
relevance of any individual simulation but also in the extremely active field
of developing computational methods. Here we map the trajectory assessment
problem onto a simple statistical calculation of the ``effective sample size''
- i.e., the number of statistically independent configurations. The mapping is
achieved by asking the question, ``How much time must elapse between snapshots
included in a sample for that sample to exhibit the statistical properties
expected for independent and identically distributed configurations?'' The
resulting ``structural de-correlation time'' is robustly calculated using exact
properties deduced from our previously developed ``structural histograms,''
without any fitting parameters. We show the method is equally and directly
applicable to toy models, peptides, and a 72-residue protein model. Variants of
our approach can readily be applied to a wide range of physical and chemical
systems.


  Many people are familiar with the physico-chemical properties of gene
sequences. In this paper I present a mathematical perspective: how do
mathematical principles such as information theory, coding theory, and
combinatorics influence the beginnings of life and the formation of the genetic
codes we observe today? What constraints on possible life forms are imposed by
information-theoretical concepts? Further, I detail how mathematical principles
can help us to analyse the genetic sequences we observe in the world today.


  F-actin bundles constitute principal components of a multitude of
cytoskeletal processes including stereocilia, filopodia, microvilli,
neurosensory bristles, cytoskeletal stress fibers, and the sperm acrosome. The
bending, buckling, and stretching behaviors of these processes play key roles
in cellular functions ranging from locomotion to mechanotransduction and
fertilization. Despite their central importance to cellular function, F-actin
bundle mechanics remain poorly understood. Here, we demonstrate that bundle
bending stiffness is a state-dependent quantity with three distinct regimes
that are mediated by bundle dimensions in addition to crosslink properties. We
calculate the complete state-dependence of the bending stiffness and elucidate
the mechanical origin of each. A generic set of design parameters delineating
the regimes in state-space is derived and used to predict the bending stiffness
of a variety of F-actin bundles found in cells. Finally, the broad and direct
implications that the isolated state-dependence of F-actin bundle stiffness has
on the interpretation of the bending, buckling, and stretching behavior of
cytoskeletal bundles is addressed.


  Replacement of papers q-bio/0605026 and q-bio/0607043 by a single manuscript
(q-bio/0701029) follows suggestions from a magazine's referee


  We prove that the Laessig-Wiese (LW) field theory for the freezing transition
of the secondary structure of random RNA is renormalizable to all orders in
perturbation theory. The proof relies on a formulation of the model in terms of
random walks and on the use of the multilocal operator product expansion.
Renormalizability allows us to work in the simpler scheme of open polymers, and
to obtain the critical exponents at 2-loop order. It also allows to prove some
exact exponent identities, conjectured in LW.


  The phase diagrams survival-extinction for the Penna model with parameters:
(mutations rate)-(birth rate), (mutation rate)-(harmful mutations threshold),
(harmful mutation threshold)-(minimal reproduction age) are presented. The
extinction phase may be caused by either mutational meltdown or overpopulation.
When the Verhulst factor is responsible for removing only newly born babies and
does not act on adults the overpopulation is avoided and only genetic factors
may lead to species extinction.


  Major advances in large-scale yeast two hybrid (Y2H) screening have provided
a global view of binary protein-protein interactions across species as
dissimilar as human, yeast, and bacteria. Remarkably, these analyses have
revealed that all species studied have a degree distribution of protein-protein
binding that is approximately scale-free (varies as a power law) even though
their evolutionary divergence times differ by billions of years. The universal
power-law shows only the surface of the rich information harbored by these
high-throughput data. We develop a detailed mathematical model of the
protein-protein interaction network based on association free energy, the
biochemical quantity that determines protein-protein interaction strength. This
model reproduces the degree distribution of all of the large-scale Y2H data
sets available and allows us to extract the distribution of free energy, the
likelihood that a pair of proteins of a given species will bind. We find that
across-species interactomes have significant differences that reflect the
strengths of the protein-protein interaction. Our results identify a global
evolutionary shift: more evolved organisms have weaker binary protein-protein
binding. This result is consistent with the evolution of increased protein
unfoldedness and challenges the dogma that only specific protein-protein
interactions can be biologically functional..