  We evolve network topology of an asymmetrically connected threshold network
by a simple local rewiring rule: quiet nodes grow links, active nodes lose
links. This leads to convergence of the average connectivity of the network
towards the critical value $K_c =2$ in the limit of large system size $N$. How
this principle could generate self-organization in natural complex systems is
discussed for two examples: neural networks and regulatory networks in the
genome.


  We suggest to simulate evolution of complex organisms constrained by the sole
requirement of robustness in their expression patterns. This scenario is
illustrated by evolving discrete logical networks with epigenetic properties.
Evidence for dynamical features in the evolved networks is found that can be
related to biological observables.


  The idea of this review is to connect the different models of evolution to
those of biological ageing through Darwin's theory. We start with the Eigen
model of quasispecies for microevolution, then introduce the Bak-Sneppen model
for macroevolution and, finally, present the Penna model for biological ageing
and some of its most important results. We also explore the concept of
coevolution using this model.


  A method for generating angular forces around $\sigma$-bonded transition
metal ions is generalized to treat $\pi$-bonded configurations. The theoretical
approach is based on an analysis of a ligand-field Hamiltonian based on the
moments of the electron state distribution. The functional forms that are
obtained involve a modification of the usual expression of the binding energy
as a sum of ligand-ligand interactions, which however requires very little
increased in CPU time. The angular interactions have simple forms involving sin
and cos functions, whose relative weights depend on whether the ligands are
$\sigma$- or $\pi$-bonded. They describe the ligand-field stabilization energy
to an accuracy of about 10%. The resulting force field is used to model the
structure of small clusters, including fragments of the copper blue protein
structure. Large deviations from the typical square copper coordination are
found when $\pi$-bonded ligands are present.


  We study a class of micromanipulation experiments, exemplified by the pulling
apart of the two strands of double-stranded DNA (dsDNA). When the pulling force
is increased to a critical value, an ``unzipping'' transition occurs. For
random DNA sequences with short-ranged correlations, we obtain exact results
for the number of monomers liberated and the specific heat, including the
critical behavior at the transition. Related systems include a random
heteropolymer pulled away from an adsorbing surface and a vortex line in a type
II superconductor tilted away from a fragmented columnar defect.


  We examine the behavior of a model which describes the melting of
double-stranded DNA chains. The model, with displacement-dependent stiffness
constants and a Morse on-site potential, is analyzed numerically; depending on
the stiffness parameter, it is shown to have either (i) a second-order
transition with "nu_perpendicular" = - beta = 1, "nu_parallel" = gamma/2 = 2
(characteristic of short range attractive part of the Morse potential) or (ii)
a first-order transition with finite melting entropy, discontinuous fraction of
bound pairs, divergent correlation lengths, and critical exponents
"nu_perpendicular" = - beta = 1/2, "nu_parallel" = gamma/2 = 1.


  Many processes in biology, from the regulation of gene expression in bacteria
to memory in the brain, involve switches constructed from networks of
biochemical reactions. Crucial molecules are present in small numbers, raising
questions about noise and stability. Analysis of noise in simple reaction
schemes indicates that switches stable for years and switchable in milliseconds
can be built from fewer than one hundred molecules. Prospects for direct tests
of this prediction, as well as implications, are discussed.


  Native protein folds often have a high degree of symmetry. We study the
relationship between the symmetries of native proteins, and their
designabilities -- how many different sequences encode a given native
structure. Using a two-dimensional lattice protein model based on
hydrophobicity, we find that those native structures that are encoded by the
largest number of different sequences have high symmetry. However only certain
symmetries are enhanced, e.g. x/y-mirror symmetry and $180^o$ rotation, while
others are suppressed. If it takes a large number of mutations to destabilize
the native state of a protein, then, by definition, the state is highly
designable. Hence, our findings imply that insensitivity to mutation implies
high symmetry. It appears that the relationship between designability and
symmetry results because protein substructures are also designable. Native
protein folds may therefore be symmetric because they are composed of repeated
designable substructures.


  An asexual set of primitive bacteria is simulated with a bit-string Penna
model with a Fermi function for survival. A recent hypothesis by Jan, Stauffer
and Moseley on the evolution of sex from asexual cells as a strategy for trying
to escape the effects of deleterious mutations is checked. This strategy is
found to provide a successful scenario for the evolution of a stable
macroscopic sexual population.


  For asexual organisms point mutations correspond to local displacements in
the genotypic space, while other genotypic rearrangements represent long-range
jumps. We investigate the spreading properties of an initially homogeneous
population in a flat fitness landscape, and the equilibrium properties on a
smooth fitness landscape. We show that a small-world effect is present: even a
small fraction of quenched long-range jumps makes the results indistinguishable
from those obtained by assuming all mutations equiprobable. Moreover, we find
that the equilibrium distribution is a Boltzmann one, in which the fitness
plays the role of an energy, and mutations that of a temperature.


  A stochastic genetic model for biological aging is introduced bridging the
gap between the bit-string Penna model and the Pletcher-Neuhauser approach. The
phenomenon of exponentially increasing mortality function at intermediate ages
and its deceleration at advanced ages is reproduced for both the evolutionary
steady-state population and the genetically homogeneous individuals.


  We have studied folding mechanisms of three small globular proteins: crambin
(CRN), chymotrypsin inhibitor 2 (CI2) and the fyn Src Homology 3 domain (SH3)
which are modelled by a Go-like Hamiltonian with the Lennard-Jones
interactions. It is shown that folding is dominated by a well-defined
sequencing of events as determined by establishment of particular contacts. The
order of events depends primarily on the geometry of the native state.
Variations in temperature, coupling strengths and viscosity affect the
sequencing scenarios to a rather small extent. The sequencing is strongly
correlated with the distance of the contacting aminoacids along the sequence.
Thus $\alpha$-helices get established first. Crambin is found to behave like a
single-route folder, whereas in CI2 and SH3 the folding trajectories are more
diversified. The folding scenarios for CI2 and SH3 are consistent with
experimental studies of their transition states.


  This work addresses the question of whether it is possible to define simple
pair-wise interaction terms to approximate free energies of proteins or
polymers. Rather than ask how reliable a potential of mean force is, one can
ask how reliable it could possibly be. In a two-dimensional, infinite lattice
model system one can calculate exact free energies by exhaustive enumeration. A
series of approximations were fitted to exact results to assess the feasibility
and utility of pair-wise free energy terms. Approximating the true free energy
with pair-wise interactions gives a poor fit with little transferability
between systems of different size. Adding extra artificial terms to the
approximation yields better fits, but does not improve the ability to
generalise from one system size to another. Further, one cannot distinguish
folding from non-folding sequences via the approximated free energies. Most
usefully, the methodology shows how one can assess the utility of various terms
in lattice protein/polymer models.


  We characterize the distributions of short cycles in a large metabolic
network previously shown to have small world characteristics and a power law
degree distribution. Compared with three classes of random networks, including
Erdoes-Renyi random graphs and synthetic small world networks of the same
connectivity, the metabolic network has a particularly large number of
triangles and a deficit in large cycles. Short cycles reduce the length of
detours when a connection is clipped, so we propose that long cycles in
metabolism may have been selected against in order to shorten transition times
and reduce the likelihood of oscillations in response to external
perturbations.


  A novel mechanism for sympatric speciation that takes into account complex
bio-processes within each individual organism is proposed. According to
dynamical systems theory, organisms with identical genotypes can possess
differentiated physiological states and may coexist `symbiotically' through
appropriate mutual interaction. With mutations, the phenotypically
differentiated organisms gradually come to possess distinct genotypes, while
maintaining their symbiotic relationship. This symbiotic speciation is robust
against sexual recombination, because offspring of mixed parentage, with
intermediate genotypes, are less fit than their parents. This leads to
sterility of the hybrid. Accordingly, a basis for mating preference also
arises.


  A coarse grained description of a two-dimensional prey-predator system is
given in terms of a 3-state lattice model containing two control parameters:
the spreading rates of preys and predators. The properties of the model are
investigated by dynamical mean-field approximations and extensive numerical
simulations. It is shown that the stationary state phase diagram is divided
into two phases: a pure prey phase and a coexistence phase of preys and
predators in which temporal and spatial oscillations can be present. The
different type of phase transitions occuring at the boundary of the prey
absorbing phase, as well as the crossover phenomena occuring between the
oscillatory and non-oscillatory domains of the coexistence phase are studied.
The importance of finite size effects are discussed and scaling relations
between different quantities are established. Finally, physical arguments,
based on the spatial structure of the model, are given to explain the
underlying mechanism leading to oscillations.


  Optimal geometrical arrangements, such as the stacking of atoms, are of
relevance in diverse disciplines. A classic problem is the determination of the
optimal arrangement of spheres in three dimensions in order to achieve the
highest packing fraction; only recently has it been proved that the answer for
infinite systems is a face-centred-cubic lattice. This simply stated problem
has had a profound impact in many areas, ranging from the crystallization and
melting of atomic systems, to optimal packing of objects and subdivision of
space. Here we study an analogous problem--that of determining the optimal
shapes of closely packed compact strings. This problem is a mathematical
idealization of situations commonly encountered in biology, chemistry and
physics, involving the optimal structure of folded polymeric chains. We find
that, in cases where boundary effects are not dominant, helices with a
particular pitch-radius ratio are selected. Interestingly, the same geometry is
observed in helices in naturally-occurring proteins.


  A novel scheme is introduced to capture the spatial correlations of
consecutive amino acids in naturally occurring proteins. This knowledge-based
strategy is able to carry out optimally automated subdivisions of protein
fragments into classes of similarity. The goal is to provide the minimal set of
protein oligomers (termed ``oligons'' for brevity) that is able to represent
any other fragment. At variance with previous studies where recurrent local
motifs were classified, our concern is to provide simplified protein
representations that have been optimised for use in automated folding and/or
design attempts. In such contexts it is paramount to limit the number of
degrees of freedom per amino acid without incurring in loss of accuracy of
structural representations. The suggested method finds, by construction, the
optimal compromise between these needs. Several possible oligon lengths are
considered. It is shown that meaningful classifications cannot be done for
lengths greater than 6 or smaller than 4. Different contexts are considered
were oligons of length 5 or 6 are recommendable. With only a few dozen of
oligons of such length, virtually any protein can be reproduced within typical
experimental uncertainties. Structural data for the oligons is made publicly
available.


  We study the statistical properties of hydrophobic/polar model sequences with
unique native states on the square lattice. It is shown that this ensemble of
sequences differs from random sequences in significant ways in terms of both
the distribution of hydrophobicity along the chains and total hydrophobicity.
Whenever statistically feasible, the analogous calculations are performed for a
set of real enzymes, too.


  We study the thermodynamic behavior of a model protein with 54 amino acids
that forms a three-helix bundle in its native state. The model contains three
types of amino acids and five to six atoms per amino acid and has the
Ramachandran torsional angles $\phi_i$, $\psi_i$ as its degrees of freedom. The
force field is based on hydrogen bonds and effective hydrophobicity forces. For
a suitable choice of the relative strength of these interactions, we find that
the three-helix-bundle protein undergoes an abrupt folding transition from an
expanded state to the native state. Also shown is that the corresponding one-
and two-helix segments are less stable than the three-helix sequence.


  We present a quantitative analysis of recent data on the kinetics of ATP
hydrolysis, which has presented a puzzle regarding the load dependence of the
Michaelis constant. Within the framework of coarse grained two-state ratchet
models, our analysis not only explains the puzzling data, but provides a
modified Michaelis law, which could be useful as a guide for future
experiments.


  We reply to the Comment by Hartmann (cond-mat/9908132) on our paper Phys.
Rev. Lett. 84 (2000) 2026 (also cond-mat/9907125).


  The prediction of the three-dimensional structures of the native state of
proteins from the sequences of their amino acids is one of the most important
challenges in molecular biology. An essential ingredient to solve this problem
within coarse-grained models is the task of deducing effective interaction
potentials between the amino acids. Over the years several techniques have been
developed to extract potentials that are able to discriminate satisfactorily
between the native and non-native folds of a pre-assigned protein sequence. In
general, when these potentials are used in actual dynamical folding
simulations, they lead to a drift of the native structure outside the
quasi-native basin. In this study, we present and validate an approach to
overcome this difficulty. By exploiting several numerical and analytical tools
we set up a rigorous iterative scheme to extract potentials satisfying a
pre-requisite of any viable potential: the stabilization of proteins within
their native basin (less than 3-4 \AA$ $ cRMS). The scheme is flexible and is
demonstrated to be applicable to a variety of parametrizations of the energy
function and provides, in each case, the optimal potentials.


  Noise induced changes in the critical and oscillatory behavior of a
Prey-Predator system are studied using power spectrum density and Spectral
Amplification Factor (SAF) analysis. In the absence of external noise, the
population densities exhibit three kinds of asymptotic behavior, namely:
Absorbing State, Fixed Point (FP) and an Oscillatory Regime (OR) with a well
defined proper (natural) frequency. The addition of noise destabilizes the FP
phase inducing a transition to a new OR. Surprisingly, it is found that when a
periodic signal is added to the control parameter, the system responds
robustly, without relevant changes in its behavior. Nevertheless, the "Coherent
Stochastic Resonance" phenomenon is found only at the proper frequency. Also, a
method based on SAF allows us to locate very accurately the transition points
between the different regimes.


  The long-term behaviour of dynamic systems can be classified in two different
regimes, regular or chaotic, depending on the values of the control parameters,
which are kept constant during the time evolution. Starting from slightly
different initial conditions, a regular system converges to the same final
trajectory, whereas a chaotic system follows two distinct trajectories
exponentially diverging from each other. In spite of these differences, regular
and chaotic systems share a common property: both arrive exponentially fast to
their final destiny, becoming trapped there. In both cases one has finite
transient times. This is not a profitable property in what concerns
evolutionary strategies, where the eternal search for new forms, better than
the current one, is imperative. That is why evolutionary dynamic systems tend
to tune themselves in very particular situations in between regular and chaotic
regimes. These particular situations present eternal transients, and the system
actually never reaches its final destiny, preserving diversity. This feature
allows the system to visit other regions of the space of possibilities, not
only the tiny region covered by its final attractor.


  Drug resistance to HIV-1 Protease involves accumulation of multiple mutations
in the protein. Here we investigate the role of these mutations by using
molecular dynamics simulations which exploit the influence of the native-state
topology in the folding process. Our calculations show that sites contributing
to phenotypic resistance of FDA-approved drugs are among the most sensitive
positions for the stability of partially folded states and should play a
relevant role in the folding process. Furthermore, associations between amino
acid sites mutating under drug treatment are shown to be statistically
correlated. The striking correlation between clinical data and our calculations
suggest a novel approach to the design of drugs tailored to bind regions
crucial not only for protein function but also for folding.


  We describe quantitatively a RNA molecule under the influence of an external
force exerted at its two ends as in a typical single-molecule experiment. Our
calculation incorporates the interactions between nucleotides by using the
experimentally-determined free energy rules for RNA secondary structure and
models the polymeric properties of the exterior single-stranded regions
explicitly as elastic freely-jointed chains. We find that in spite of
complicated secondary structures, force-extension curves are typically smooth
in quasi-equilibrium. We identify and characterize two
sequence/structure-dependent mechanisms that, in addition to the
sequence-independent entropic elasticity of the exterior single-stranded
regions, are responsible for the smoothness. These involve compensation between
different structural elements on which the external force acts simultaneously,
and contribution of suboptimal structures, respectively. We estimate how many
features a force-extension curve recorded in non-equilibrium, where the pulling
proceeds faster than rearrangements in the secondary structure of the molecule,
could show in principle. Our software is available to the public through a
`RNA-pulling server'.


  Force-induced structural transitions both in relatively random and in
designed single-stranded DNA (ssDNA) chains are studied theoretically. At high
salt conditions, ssDNA forms compacted hairpin patterns stabilized by
base-pairing and base-pair stacking interactions, and a threshold external
force is needed to pull the hairpinned structure into a random coiled one. The
base-pair stacking interaction in the ssDNA chain makes this hairpin-coil
conversion a discontinuous (first-order) phase transition process characterized
by a force plateau in the force-extension curve, while lowering this potential
below some critical level turns this transition into continuous (second-order)
type, no matter how strong the base-pairing interaction is. The phase diagram
(including hairpin-I, -II, and random coil) is discussed as a function of
stacking potential and external force. These results are in quantitative
agreement with recent experimental observations of different ssDNA sequences,
and they reveal the necessity to consider the base-pair stacking interactions
in order to understand the structural formation of RNA, a polymer designed by
nature itself. The theoretical method used may be extended to study the
long-range interaction along double-stranded DNA caused by the topological
constraint of fixed linking number.


  We find that the hypothesis made by Jan, Stauffer and Moseley [Theory in
Biosc., 119, 166 (2000)] for the evolution of sex, namely a strategy devised to
escape extinction due to too many deleterious mutations, is sufficient but not
necessary for the successful evolution of a steady state population of sexual
individuals within a finite population. Simply allowing for a finite
probability for conversion to sex in each generation also gives rise to a
stable sexual population, in the presence of an upper limit on the number of
deleterious mutations per individual. For large values of this probability, we
find a phase transition to an intermittent, multi-stable regime. On the other
hand, in the limit of extremely slow drive, another transition takes place to a
different steady state distribution, with fewer deleterious mutations within
the asexual population.


  We have carried out experiments to probe the static and dynamic interfacial
properties of $\beta$--casein monolayers spread at the air-buffer interface,
and analysed these results in the context of models of weak polyampholytes.
Measurements have been made systematically over a wide range of ionic strength
and pH. In the semi-dilute regime of surface concentration a scaling exponent,
which can be linked to the degree of chain swelling, is found. This shows that
at pH close to the isoelectric point, the protein is compact. At pH away from
the isoelectric pH the protein is extended. The transition between compact and
extended states is continuous. As a function of increasing ionic strength, we
observe swelling of the protein at the isoelectric pH but contraction of the
protein at pH values away from it. These behaviours are typical of a those
predicted theoretically for a weak polyampholyte. Dilational moduli
measurements, made as a function of surface concentration exhibit maxima that
are linked to the collapse of hydrophilic regions of the protein into the
subphase. Based on this data we present a configuration map of the protein
configuration in the monolayer. These findings are supported by strain (surface
pressure) relaxation measurements and surface quasi-elastic light scattering
(SQELS) measurements which suggest the existence of loops and tails in the
subphase at higher surface concentrations.


  We present a theoretical analysis of the structural and mechanical properties
of the 30-nm chromatin fiber. Our study is based on the two-angle model
introduced by Woodcock et al. (Woodcock, C. L., S. A. Grigoryev, R. A.
Horowitz, and N. Whitaker. 1993. PNAS 90:9021-9025) that describes the
chromatin fiber geometry in terms of the entry-exit angle of the nucleosomal
DNA and the rotational setting of the neighboring nucleosomes with respect to
each other. We explore analytically the different structures that arise from
this building principle, and demonstrate that the geometry with the highest
density is close to the one found in native chromatin fibers under
physiological conditions. On the basis of this model we calculate mechanical
properties of the fiber under stretching. We obtain expressions for the
stress-strain characteristics which show good agreement with the results of
recent stretching experiments (Cui, Y., and C. Bustamante. 2000. PNAS
97:127-132) and computer simulations (Katritch, V., C. Bustamante, and W. K.
Olson. 2000. J. Mol. Biol. 295:29-40), and which provide simple physical
insights into correlations between the structural and elastic properties of
chromatin.


  We consider how beads can diffuse along a chain that wraps them, without
becoming displaced from the chain; our proposed mechanism is analogous to the
reptation of "stored length" in more familiar situations of polymer dynamics.
The problem arises in the case of globular aggregates of proteins (histones)
that are wound by DNA in the chromosomes of plants and animals; these beads
(nucleosomes) are multiply wrapped and yet are able to reposition themselves
over long distances, while remaining bound by the DNA chain.


  We consider the statistical properties of primary sequences of two-letter HP
copolymers (H for hydrophobic and P for polar) designed to have water soluble
globular conformations with H monomers shielded from water inside the shell of
P monomers. We show, both by computer simulations and by exact analytical
calculation, that for large globules and flexible polymers such sequences
exhibit long-range correlations which can be described by Levy-flight
statistics.


  Using a bit-string model of evolution, we find a successful route to diploidy
and sex in simple organisms. Allowing the sexually reproducing diploid
individuals to also perform mitosis, as they do in a haploid-diploid cycle,
leads to the complete takeover of the population by sexual diploids. This
mechanism is so robust, that even the accidental conversion and pairing of only
two diploids give rise to a sexual population.


  We consider DNA-cationic lipid complexes that form lamellar stacks of lipid
bilayers with parallel DNA strands intercalated in between. We calculate the
electrostatically induced elastic deformations of the lipid bilayers. It is
found that the membranes undulate with a periodicity that is set by the DNA
interaxial distance. As a consequence the lamellar repeat distance changes
resulting in a swelling or compression of the lamellar stack. Such undulations
may be responsible for the intermembrane coupling between DNA strands in
different layers as it is observed experimentally.


  We investigate the mechanical stability of a segment of DNA wrapped around a
histone in the nucleosome configuration. The assumption underlying this
investigation is that the proper model for this packaging arrangement is that
of an elastic rod that is free to twist and that writhes subject to mechanical
constraints. We find that the number of constraints required to stabilize the
nuclesome configuration is determined by the length of the segment, the number
of times the DNA wraps around the histone spool, and the specific constraints
utilized. While it can be shown that four constraints suffice, in principle, to
insure stability of the nucleosome, a proper choice must be made to guarantee
the effectiveness of this minimal number. The optimal choice of constraints
appears to bear a relation to the existence of a spiral ridge on the surface of
the histone octamer. The particular configuration that we investigate is
related to the 30 nanometer spiral, a higher-order organization of DNA in
chromatin.


  The critical properties of a simple prey-predator model are revisited. For
some values of the control parameters, the model exhibits a line of directed
percolation like transitions to a single absorbing state. For other values of
the control parameters one finds a second line of continuous transitions toward
infinite number of absorbing states, and the corresponding steady-state
exponents are mean-field like. The critical behavior of the special point T
(bicritical point), where the two transition lines meet, belongs to a different
universality class. The use of dynamical Monte-Carlo method shows that a
particular strategy for preparing the initial state should be devised to
correctly describe the physics of the system near the second transition line.
Relationships with a forest fire model with immunization are also discussed.


  This paper focuses on mechanical aspects of chromatin biological functioning.
Within a basic geometric modeling of the chromatin assembly, we give for the
first time the complete set of elastic constants (twist and bend persistence
lengths, stretch modulus and twist-stretch coupling constant) of the so-called
30-nm chromatin fiber, in terms of DNA elastic properties and geometric
properties of the fiber assembly. The computation naturally embeds the fiber
within a current analytical model known as the ``extensible worm-like rope'',
allowing a straightforward prediction of the force-extension curves. We show
that these elastic constants are strongly sensitive to the linker length, up to
1 bp, or equivalently to its twist, and might locally reach very low values,
yielding a highly flexible and extensible domain in the fiber. In particular,
the twist-stretch coupling constant, reflecting the chirality of the chromatin
fiber, exhibits steep variations and sign changes when the linker length is
varied.
  We argue that this tunable elasticity might be a key feature for chromatin
function, for instance in the initiation and regulation of transcription.


  As sections of a strand duplexed DNA denature when exposed to high
temperature, the excess linking number is taken up by the undenatured portions
of the molecule. The mechanical energy that arises because of the overwinding
of the undenatured sections can, in principle, alter the nature of the thermal
denaturation process. Assuming that the strains associated with this
overwinding are not relieved, we find that a simple model of strain-altered
melting leads to a suppression of the melting transition when the unaltered
transition is continuous. When the melting transition is first order in the
absence of strain associated with overwinding, the modification is to a third
order phase transition.


  We review recent advances in the physics of strongly interacting charged
systems functioning in water at room temperature. We concentrate on the
phenomena which go beyond the framework of mean field theories, whether linear
Debye-Huckel or non-linear Poisson-Boltzmann. We place major emphasis on charge
inversion - a counterintuitive phenomenon in which a strongly charged particle,
called macroion, binds so many counterions that its net charge changes sign. We
discuss the universal theory of charge inversion based on the idea of a
strongly correlated liquid of adsorbed counterions, similar to a Wigner
crystal. This theory has a vast array of applications, particularly in biology
and chemistry; for example, the DNA double helix in the presence of positive
multivalent ions (e.g., polycations) acquires a net positive charge and drifts
as a positive particle in electric field. This simplifies DNA uptake by the
cell as needed for gene therapy, because the cell membrane is negatively
charged. We discuss also the analogies of charge inversion in other fields of
physics.


  We study the dynamics of protein folding via statistical energy-landscape
theory. In particular, we concentrate on the local-connectivity case with the
folding progress described by the fraction of native conformations. We obtain
information for the first passage-time (FPT) distribution and its moments. The
results show a dynamic transition temperature below which the FPT distribution
develops a power-law tail, a signature of the intermittency phenomena of the
folding dynamics. We also discuss the possible application of the results to
single-molecule dynamics experiments.


  I consider a class of fitness landscapes, in which the fitness is a function
of a finite number of phenotypic "traits", which are themselves linear
functions of the genotype. I show that the stationary trait distribution in
such a landscape can be explicitly evaluated in a suitably defined
"thermodynamic limit", which is a combination of infinite-genome and strong
selection limits. These considerations can be applied in particular to identify
relevant features of the evolution of promoter binding sites, in spite of the
shortness of the corresponding sequences.


  We investigate the orientational order of transverse polarization vectors of
long, stiff polymer molecules and their coupling to bond orientational and
positional order in high density mesophases. Homogeneous ordering of transverse
polarization vector promotes distortions in the hexatic phase, whereas
inhomogeneous ordering precipitates crystalization of the 2D sections with
different orientations of the transverse polarization vector on each molecule
in the unit cell. We propose possible scenarios for going from the hexatic
phase, through the distorted hexatic phase to the crystalline phase with an
orthorhombic unit cell observed experimentally for the case of DNA.


  Formation of RNA secondary structures is an example of the sequence-structure
problem omnipresent in biopolymers. A theoretical question of recent interest
is whether a random RNA sequence undergoes a finite temperature glass
transition. We answer this question affirmatively by first establishing the
perturbative stability of the high temperature phase via a two replica
calculation. Subsequently, we show that this phase cannot persist down to zero
temperature by considering energetic contributions due to rare regions of
complementary subsequences.


  We consider a spatial model related to bond percolation for the spread of a
disease that includes variation in the susceptibility to infection. We work on
a lattice with random bond strengths and show that with strong disorder, i.e. a
wide range of variation of susceptibility, patchiness in the spread of the
epidemic is very likely, and the criterion for epidemic outbreak depends
strongly on the disorder. These results are qualitatively different from those
of standard models in epidemiology, but correspond to real effects.


  We study the thermodynamic behavior of a model protein with 54 amino acids
that is designed to form a three-helix bundle in its native state. The model
contains three types of amino acids and five to six atoms per amino acid, and
has the Ramachandran torsion angles as its only degrees of freedom. The force
field is based on hydrogen bonds and effective hydrophobicity forces. We study
how the character of the collapse transition depends on the strengths of these
forces. For a suitable choice of these two parameters, it is found that the
collapse transition is first-order-like and coincides with the folding
transition. Also shown is that the corresponding one- and two-helix segments
make less stable secondary structure than the three-helix sequence.


  The formation of secondary structures by a random RNA sequence is studied as
a model system for the sequence-structure problem omnipresent in biopolymers.
Several toy energy models are introduced to allow detailed analytical and
numerical studies. First, a two-replica calculation is performed. By mapping
the two-replica problem to the denaturation of a single homogeneous RNA in
6-dimensional embedding space, we show that sequence disorder is perturbatively
irrelevant, i.e., an RNA molecule with weak sequence disorder is in a molten
phase where many secondary structures with comparable total energy coexist. A
numerical study of various models at high temperature reproduces behaviors
characteristic of the molten phase. On the other hand, a scaling argument based
on the extremal statistics of rare regions can be constructed to show that the
low temperature phase is unstable to sequence disorder. We performed a detailed
numerical study of the low temperature phase using the droplet theory as a
guide, and characterized the statistics of large-scale, low-energy excitations
of the secondary structures from the ground state structure. We find the
excitation energy to grow very slowly (i.e., logarithmically) with the length
scale of the excitation, suggesting the existence of a marginal glass phase.
The transition between the low temperature glass phase and the high temperature
molten phase is also characterized numerically. It is revealed by a change in
the coefficient of the logarithmic excitation energy, from being disorder
dominated to entropy dominated.


  We consider statistical patterns of variation in growth rates for over 400
species of breeding birds across North America surveyed from 1966 to 1998. We
report two results. First, the standard deviation of population growth rates
decays as a power-law function of total population size with an exponent $\beta
= 0.36\pm0.02$. Second, the number of subpopulations, measured as the number of
survey locations with non-zero counts, scales to the 3/4-power of total number
of birds counted in a given species. We show how these patterns may be related,
and discuss a simple stochastic growth model for a geographically subdivided
population that formalizes the relationship. We also examine reasons that may
explain why some species deviate from these scaling-laws.


  We propose a scenario for the prebiotic co-evolution of RNA and of fast
folding proteins with large entropy gaps as observed today. We show from very
general principles that the folding and unfolding of the proteins synthesized
by RNA can function as a heat pump. Rock surfaces can facilitate the folding of
amino acid chains having polar and hydrophobic residues, with an accompanying
heat loss to the surrounding rock. These chains then absorb heat from the soup
as they unfold. This opens the way to the enhancement of RNA replication rates,
by the enzymatic action of folded proteins present in greater numbers at
reduced temperatures. This gives an evolutionary advantage to those RNA coding
amino acid sequences with non-degenerate folded states which would provide the
most efficient refrigeration.


  We present an extensive theoretical investigation of the mechanical unzipping
of double-stranded DNA under the influence of an applied force. In the limit of
long polymers, there is a thermodynamic unzipping transition at a critical
force value of order 10 pN, with different critical behavior for homopolymers
and for random heteropolymers. We extend results on the disorder-averaged
behavior of DNA's with random sequences to the more experimentally accessible
problem of unzipping a single DNA molecule. As the applied force approaches the
critical value, the double-stranded DNA unravels in a series of discrete,
sequence-dependent steps that allow it to reach successively deeper energy
minima. Plots of extension versus force thus take the striking form of a series
of plateaus separated by sharp jumps. Similar qualitative features should
reappear in micromanipulation experiments on proteins and on folded RNA
molecules. Despite their unusual form, the extension versus force curves for
single molecules still reveal remnants of the disorder-averaged critical
behavior. Above the transition, the dynamics of the unzipping fork is related
to that of a particle diffusing in a random force field; anomalous,
disorder-dominated behavior is expected until the applied force exceeds the
critical value for unzipping by roughly 5 pN.


  Using a fast tree-searching algorithm and a Pentium cluster, we enumerated
all the sequences and compact conformations (structures) for a protein folding
model on a cubic lattice of size $4\times3\times3$. We used two types of amino
acids -- hydrophobic (H) and polar (P) -- to make up the sequences, so there
were $2^{36} \approx 6.87 \times 10^{10}$ different sequences. The total number
of distinct structures was 84,731,192. We made use of a simple solvation model
in which the energy of a sequence folded into a structure is minus the number
of hydrophobic amino acids in the ``core'' of the structure. For every
sequence, we found its ground state or ground states, i.e., the structure or
structures for which its energy is lowest. About 0.3% of the sequences have a
unique ground state. The number of structures that are unique ground states of
at least one sequence is 2,662,050, about 3% of the total number of structures.
However, these ``designable'' structures differ drastically in their
designability, defined as the number of sequences whose unique ground state is
that structure. To understand this variation in designability, we studied the
distribution of structures in a high dimensional space in which each structure
is represented by a string of 1's and 0's, denoting core and surface sites,
respectively.


  The Prisoner's dilemma is the main game theoretical framework in which the
onset and maintainance of cooperation in biological populations is studied. In
the spatial version of the model, we study the robustness of cooperation in
heterogeneous ecosystems in spatial evolutionary games by considering site
diluted lattices. The main result is that due to disorder, the fraction of
cooperators in the population is enhanced. Moreover, the system presents a
dynamical transition at $\rho^*$, separating a region with spatial chaos from
one with localized, stable groups of cooperators.


  We introduce a simple theoretical approach for an equilibrium study of
proteins with known native state structures. We test our approach with results
on well-studied globular proteins, Chymotrypsin Inhibitor (2ci2), Barnase and
the alpha spectrin SH3 domain and present evidence for a hierarchical onset of
order on lowering the temperature with significant organization at the local
level even at high temperatures. A further application to the folding process
of HIV-1 protease shows that the model can be reliably used to identify key
folding sites that are responsible for the development of drug resistance .


  We introduce a graph generating model aimed at representing the evolution of
protein interaction networks. The model is based on the hypotesis of evolution
by duplications and divergence of the genes which produce proteins. The
obtained graphs shows multifractal properties recovering the absence of a
characteristic connectivity as found in real data of protein interaction
networks. The error tolerance of the model to random or targeted damage is in
very good agreement with the behavior obtained in real protein networks
analysis. The proposed model is a first step in the identification of the
evolutionary dynamics leading to the development of protein functions and
interactions.


  We propose an algorithmic strategy for improving the efficiency of Monte
Carlo searches for the low-energy states of proteins. Our strategy is motivated
by a model of how proteins alter their shapes. In our model when proteins fold
under physiological conditions, their backbone dihedral angles change
synchronously in groups of four or more so as to avoid steric clashes and
respect the kinematic conservation laws. They wriggle; they do not thrash. We
describe a simple algorithm that can be used to incorporate wriggling in Monte
Carlo simulations of protein folding. We have tested this wriggling algorithm
against a code in which the dihedral angles are varied independently
(thrashing). Our standard of success is the average root-mean-square distance
(rmsd) between the alpha-carbons of the folding protein and those of its native
structure. After 100,000 Monte Carlo sweeps, the relative decrease in the mean
rmsd, as one switches from thrashing to wriggling, rises from 11% for the
protein 3LZM with 164 amino acids (aa) to 40% for the protein 1A1S with 313 aa
and 47% for the protein 16PK with 415 aa. These results suggest that wriggling
is useful and that its utility increases with the size of the protein. One may
implement wriggling on a parallel computer or a computer farm.


  Computer modelling for evolutionary systems consists in: 1) to store in the
memory the individual features of each member of a large population; and 2) to
update the whole system repeatedly, as time goes by, according to some
prescribed rules (reproduction, death, ageing, etc) where some degree of
randomness is included through pseudo-random number sequences. Compared to
direct observation of Nature, this approach presents two distinguishing
features. First, one can follow the characteristics of the system in real time,
instead of only observing the current, static situation which is a long-term
consequence of a remote past completely unknown except for some available
fossil snapshots. In particular, one can repeat the whole dynamical process,
starting from the same initial population, using the same randomness, changing
only some minor contingency during the process, in order to study its long-term
consequences. Second, evolution necessarily follows a critical dynamics with
long-term memory characteristics, equivalent to the long-range correlations
responsible for the well known universality properties of static critical
phenomena. Accordingly, some strong simplifications can be applied, allowing
one to obtain many characteristics of real populations from toy models easily
implementable on the computer.


  Co-evolution of two coupled quasispecies is studied, motivated by the
competition between viral evolution and adapting immune response. In this
co-adaptive model, besides the classical error catastrophe for high virus
mutation rates, a second ``adaptation-'' catastrophe occurs, when virus
mutation rates are too small to escape immune attack. Maximizing both regimes
of viral error catastrophes is a possible strategy for an optimal immune
response, reducing the range of allowed viral mutation rates to a minimum. From
this requirement one obtains constraints on B-cell mutation rates and receptor
lengths, yielding an estimate of somatic hypermutation rates in the germinal
center in accordance with observation.


  We characterize the low temperature phase of a simple model for RNA secondary
structures by determining the typical energy scale E(l) of excitations
involving l bases. At zero temperature, we find a scaling law E(l) \sim
l^\theta with \theta \approx 0.23, and this same scaling holds at low enough
temperatures. Above a critical temperature, there is a different phase
characterized by a relatively flat free energy landscape resembling that of a
homopolymer with a scaling exponent \theta=1. These results strengthen the
evidence in favour of the existence of a glass phase at low temperatures.


  In the diffusion-collision model, the unfolding rates are given by the
likelihood of secondary structural cluster dissociation. In this work, we
introduce an unfolding rate calculation for proteins whose secondary structural
elements are $\alpha$-helices, modeled from thermal escape over a barrier which
arises from the free energy in buried hydrophobic residues. Our results are in
good agreement with currently accepted values for the attempt rate.


  We investigate by numerical simulations and analytical calculations the
Bak-Sneppen model for biological evolution in scale-free networks. By using
large scale numerical simulations, we study the avalanche size distribution and
the activity time behavior at nodes with different connectivities. We argue the
absence of a critical barrier and its associated critical behavior for infinite
size systems. These findings are supported by a single site mean-field analytic
treatment of the model.


  Chain-like macromolecules in solution, whether biological or synthetic,
transform from a spatially extended conformation to a compact one upon change
of temperature or solvent qualities. This sharp transition plays a key role in
various phenomena, including DNA condensation, protein folding, and the
behaviour of polymer solutions. In biological processes such as DNA
condensation the collapse is sensitively induced by a small amount of added
molecules. Here we derive a general criterion for the effect of such agents on
conformational transitions. We find two different scenarios depending on chain
stiffness. If the persistence length --the characteristic distance along which
the chain retains its direction-- is smaller than the range of attractive
correlations induced by the agent (typically up to several nanometres), the
chain contracts gradually. Stiffer chains undergo sharp collapse. We thereby
suggest that the enhanced rigidity of double-stranded DNA as compared to the
single strand is a prerequisite for sharp, controllable conformational
transitions.


  The propagation of model epidemics on a small-world network under the action
of immunization is studied. Although the connectivity in this kind of networks
is rather uniform, a vaccination strategy focused on the best connected
individuals yields a considerable improvement of disease control. The model
exhibits a transition from disease localization to propagation as the disorder
of the underlying network grows. As a consequence, for fixed disorder, a
threshold immunization level exists above which the disease remains localized.


  We propose that the degree of sweeling of the 30nm chromatin fiber (a
"measure" of its transcriptional activity) is mainly determined by the
short-ranged electrostatical interaction between different sections of the
"folded" DNA chain. These sections constitute only a small fraction of the
chain and they are located close to the entry-exit points of the DNA chain at
the nucleosome core particles. We present a model that allows to estimate the
degree of swelling of chromatin fibers as a function of salt concentration,
charge density of the strands etc. Different mechanisms by which the state of
chromatin can be controlled in vitro and in vivo are discussed.


  We report a method to control the conformation of a weak polyampholyte (the
protein $\beta$-casein) in Langmuir monolayers by light, even though the
protein is not photosensitive. Our approach is to couple the monolayer state to
a photochemical reaction excited in the liquid subphase. The conformational
transition of the protein molecule is triggered through its sensitivity to a
subphase bulk field (pH in this study), changing in the course of the
photochemical process. Thus, reaction of photoaquation of the ferrocyanide ion,
which increases the subphase pH from 7.0 to about 8.3, produces a change in the
surface monolayer pressure, $\Delta\Pi$, between -0.5 and +1.5 ${\rm mN/m}$
(depending on the surface concentration), signalling a conformational switch.
The approach proposed here can be used to selectively target and influence
different interfacial properties by light, without embedding photosensitisers
in the matrix.


  Fungal colonies are able to exhibit different morphologies depending on the
enviromental conditions. This allows them to cope with and adapt to external
changes. When grown in solid or semi-solid media the bulk of the colony is
compact and several morphological transitions have been reported to occur as
the external conditions are varied. Here we show how a unified simple
mathematical model, which includes the effect of the accumulation of toxic
metabolites, can account for the morphological changes observed. Our numerical
results are in excellent agreement with experiments carried out with the fungus
Aspergillus oryzae on solid agar.


  A reduced protein model with five to six atoms per amino acid and five amino
acid types is developed and tested on a three-helix-bundle protein, a 46-amino
acid fragment from staphylococcal protein A. The model does not rely on the
widely used Go approximation where non-native interactions are ignored. We find
that the collapse transition is considerably more abrupt for the protein A
sequence than for random sequences with the same composition. The chain
collapse is found to be at least as fast as helix formation. Energy
minimization restricted to the thermodynamically favored topology gives a
structure that has a root-mean-square deviation of 1.8 A from the native
structure. The sequence-dependent part of our potential is pairwise additive.
Our calculations suggest that fine-tuning this potential by parameter
optimization is of limited use.


  A simple model to encapsulate the essential growth properties of \emph{the
web of human sexual contacts} is presented. In the model only heterosexual
connection is considered and represented by a random growing bipartite graph
where both male-female contact networks grow simultaneously. The time evolution
of the model is analysed by a rate equation approach leading to confirm that
male and female sexual contact distributions decay as power laws with exponents
depending on influx and charisma of the sexes.


  The pair interaction between two stiff parallel linear DNA molecules depends
not only on the distance between their axes but on their azimuthal orientation.
The positional and orientational order in columnar B-DNA assemblies in solution
is investigated, based on the DNA-DNA electrostatic pair potential that takes
into account DNA helical symmetry and the amount and distribution of adsorbed
counterions. A phase diagram obtained by lattice sum calculations predicts a
variety of positionally and azimuthally ordered phases and bundling transitions
strongly depending on the counterion adsorption patterns.


  Positioning of the midcell division plane within the bacterium E. coli is
controlled by the min system of proteins: MinC, MinD and MinE. These proteins
coherently oscillate from end to end of the bacterium. We present a
reaction--diffusion model describing the diffusion of min proteins along the
bacterium and their transfer between the cytoplasmic membrane and cytoplasm.
Our model spontaneously generates protein oscillations in good agreement with
experiments. We explore the oscillation stability, frequency and wavelength as
a function of protein concentration and bacterial length.


  Mechanical stretching of secondary structures is studied through molecular
dynamics simulations of a Go-like model. Force vs. displacement curves are
studied as a function of the stiffness and velocity of the pulling device. The
succession of stretching events, as measured by the order in which contacts are
ruptured, is compared to the sequencing of events during thermal folding and
unfolding. Opposite cross-correlations are found for an $\alpha$-helix and a
$\beta$-hairpin structure. In a tandem of two $\alpha$-helices, the two
constituent helices unravel nearly simultaneously. A simple condition for
simultaneous vs. sequential unraveling of repeat units is presented.


  Mechanical stretching of the I27 domain of titin and of its double and triple
repeats are studied through molecular dynamics simulations of a Go-like model
with Lennard-Jones contact interactions. We provide a thorough characterization
of the system and correlate the sequencing of the folding and unraveling events
with each other and with the contact order. The roles of cantilever stiffness
and pulling rate are studied. Unraveling of tandem titin structures has a
serial nature. The force-displacement curves in this coarse-grained model are
similar to those obtained through all atom calculations.


  The hydrophobic/polar HP model on the square lattice has been widely used to
investigate basics of protein folding. In the cases where all designing
sequences (sequences with unique ground states) were enumerated without
restrictions on the number of contacts, the upper limit on the chain length N
has been 18-20 because of the rapid exponential growth of the numbers of
conformations and sequences. We show how a few optimizations push this limit by
about 5 units. Based on these calculations, we study the statistical
distribution of hydrophobicity along designing sequences. We find that the
average number of hydrophobic and polar clumps along the chains is larger for
designing sequences than for random ones, which is in agreement with earlier
findings for N up to 18 and with results for real enzymes. We also show that
this deviation from randomness disappears if the calculations are restricted to
maximally compact structures.


  The relevance of various residue positions for the stability and the folding
characteristics of the prion protein are investigated by using molecular
dynamics simulations of models exploiting the topology of the native state.
Highly significant correlations are found between the most relevant sites in
our analysis and the single point mutations known to be associated with the
arousal of the genetic forms of prion disease (caused by the conformational
change from the cellular to the scrapie isoform). Considerable insight into the
conformational change is provided by comparing the folding process of prion and
doppel (a newly discovered protein) sharing very similar native state topology:
the folding pathways of the former can be grouped in two main classes according
to which tertiary structure contacts are formed first enroute to the native
state. For the latter a single class of pathways leads to the native state. Our
results are consistent and supportive of the recent experimental findings that
doppel lacks the scrapie isoform and that such remarkably different behavior
results from differences in the region containing the two $\beta-$strands and
the intervening helix.


  We use LINUS, a procedure developed by Srinivasan and Rose, to provide a
physical interpretation of and to predict the secondary structures of proteins.
The secondary structure type at a given site is identified by the largest
conformational bias during short time simulations. We examine the rate of
successful prediction as a function of temperature and the interaction window.
At high temperatures, there is a large propensity for the establishment of
$\beta$-strands whereas $\alpha$-helices appear only when the temperature is
lower than a certain threshold value. It is found that there exists an optimal
temperature at which the correct secondary structures are predicted most
accurately. We find that this temperature is close to the peak temperature of
the specific heat. Changing the interaction window or carrying out longer
simulations approaching equilibrium lead to little change in the optimal
success rate. Our findings are in accord with the observation by Srinivasan and
Rose that the secondary structures are mainly determined by local interactions
and they appear in the early stage of folding.


  A simple extension of existing models for protein crystallisation is
described, in which salt ions and charge neutrality are explicitly
incorporated. This provides a straightforward explanation for the shape of
protein crystallisation boundaries, the associated scaling properties seen for
lysozyme, and can also explain much of the salt dependence of the second virial
coefficient. The analysis has wider implications for the use of pair potentials
to understand protein crystallisation.


  The origin of the unusual incubation period distribution in the development
of AIDS is largely unresolved. A key factor in understanding the observed
distribution of latency periods, as well as the occurrence of infected
individuals not developing AIDS at all, is the dynamics of the long lasting
struggle between HIV and the immune system. Using a computer simulation, we
study the diversification of viral genomes under mutation and the selective
pressure of the immune system.In common infections vast spreading of viral
genomes usually does not takes place. In the case of an HIV infection this may
occur, as the virus successively weakens the immune system by depletion of CD4+
cells.In a sequence space framework, this leads to a dynamically induced
percolation transition, corresponding to the onset of AIDS. As a result, we
obtain the prolongated shape of the incubation period distribution, as well as
a finite fraction of non-progressors that do not develop AIDS, comparing well
with results from recent clinical research.


  Surnames and nonrecombining alleles are inherited from a single parent in a
highly similar way. A simple birth-death model with mutations can accurately
describe this process. Exponentially growing and constant populations are
investigated, and we study how different compositions of the founder
populations can be observed in present-day diversity distributions. We analyse
different quantities in the statistically stationary state, both through
analytic and numerical methods. Our results compare favourably to field data
for family sizes in several countries. We discuss the relationship between the
distribution of surnames and the genetic diversity of a population.


  We present an individual based model of evolutionary ecology.
  The reproduction rate of individuals characterized by their genome depends on
the composition of the population in genotype space. Ecological features such
as the taxonomy and the macro-evolutionary mode of the dynamics are emergent
properties.
  The macro-dynamics exhibit intermittent two mode switching with a gradually
decreasing extinction rate. The generated ecologies become gradually better
adapted as well as more complex in a collective sense.
  The form of the species abundance curve compares well with observed
functional forms. The models error threshold can be understood in terms of the
characteristics of the two dynamical modes of the system.


  A variety of experimental and theoretical studies have established that the
folding process of monomeric proteins is strongly influenced by the topology of
the native state. In particular, folding times have been shown to correlate
well with the contact order, a measure of contact locality. Our investigation
focuses on identifying additional topologic properties that correlate with
experimentally measurable quantities, such as folding rates and transition
state placement, for both two- and three-state folders. The validation against
data from forty experiments shows that a particular topologic property which
measures the interdepedence of contacts, termed cliquishness or clustering
coefficient, can account with significant accuracy both for the transition
state placement and especially for folding rates, the linear correlation
coefficient being $r=0.71$. This result can be further improved to $r=0.74$, by
optimally combining the distinct topologic information captured by cliquishness
and contact order.


  In this review, we give an introduction to the structural and functional
properties of the biological networks. We focus on three major themes: topology
of complex biological networks like the metabolic and protein-protein
interaction networks, nonlinear dynamics in gene regulatory networks and in
particular the design of synthetic genetic networks using the concepts and
techniques of nonlinear physics and lastly the effect of stochasticity on the
dynamics. The examples chosen illustrate the usefulness of interdisciplinary
approaches in the study of biological networks.


  We present a general, physically motivated non-linear and non-local advection
equation in which the diffusion of interacting random walkers competes with a
local drift arising from a kind of peer pressure. We show, using a mapping to
an integrable dynamical system, that on varying a parameter, the steady state
behaviour undergoes a transition from the standard diffusive behavior to a
localized stationary state characterized by a tailed distribution. Finally, we
show that recent empirical laws on economic growth can be explained as a
collective phenomenon due to peer pressure interaction.


  Molecular dynamics simulations are used to study structure formation in
simple model polymer chains that are subject to excluded volume and torsional
interactions. The changing conformations exhibited by chains of different
lengths under gradual cooling are followed until each reaches a state from
which no further change is possible. The interactions are chosen so that the
true ground state is a helix, and a high proportion of simulation runs succeed
in reaching this state; the fraction that manage to form defect-free helices is
a function of both chain length and cooling rate. In order to demonstrate
behavior analogous to the formation of protein tertiary structure, additional
attractive interactions are introduced into the model, leading to the
appearance of aligned, antiparallel helix pairs. The simulations employ a
computational approach that deals directly with the internal coordinates in a
recursive manner; this representation is able to maintain constant bond lengths
and angles without the necessity of treating them as an algebraic constraint
problem supplementary to the equations of motion.


  We show how concepts from statistical physics, such as order parameter,
thermodynamic limit, and quantum phase transition, translate into biological
concepts in mutation-selection models for sequence evolution and can be used
there. The article takes a biological point of view within a population
genetics framework, but contains an appendix for physicists, which makes this
correspondence clear. We analyze the equilibrium behavior of deterministic
haploid mutation-selection models. Both the forward and the time-reversed
evolution processes are considered. The stationary state of the latter is
called the ancestral distribution, which turns out as a key for the study of
mutation-selection balance. We find that it determines the sensitivity of the
equilibrium mean fitness to changes in the fitness values and discuss
implications for the evolution of mutational robustness. We further show that
the difference between the ancestral and the population mean fitness, termed
mutational loss, provides a measure for the sensitivity of the equilibrium mean
fitness to changes in the mutation rate. For a class of models in which the
number of mutations in an individual is taken as the trait value, and fitness
is a function of the trait, we use the ancestor formulation to derive a simple
maximum principle, from which the mean and variance of fitness and the trait
may be derived; the results are exact for a number of limiting cases, and
otherwise yield approximations which are accurate for a wide range of
parameters. These results are applied to (error) threshold phenomena caused by
the interplay of selection and mutation. They lead to a clarification of
concepts, as well as criteria for the existence of thresholds.


  We investigate a model protein interaction network whose links represent
interactions between individual proteins. This network evolves by the
functional duplication of proteins, supplemented by random link addition to
account for mutations. When link addition is dominant, an infinite-order
percolation transition arises as a function of the addition rate. In the
opposite limit of high duplication rate, the network exhibits giant structural
fluctuations in different realizations. For biologically-relevant growth rates,
the node degree distribution has an algebraic tail with a peculiar rate
dependence for the associated exponent.


  Molecular motors interacting with cytoskeletal filaments undergo peculiar
random walks consisting of alternating sequences of directed movements along
the filaments and diffusive motion in the surrounding solution. An ensemble of
motors is studied which interacts with a single filament in two and three
dimensions. The time evolution of the probability distribution for the bound
and unbound motors is determined analytically. The diffusion of the motors is
strongly enhanced parallel to the filament. The analytical expressions are in
excellent agreement with the results of Monte Carlo simulations.


  We study the designability of all compact 3x3x3 and 6x6 lattice-protein
structures using the Miyazawa-Jernigan (MJ) matrix. The designability of a
structure is the number of sequences that design the structure, i.e. sequences
that have that structure as their unique lowest-energy state. Previous studies
of hydrophobic-polar (HP) models showed a wide distribution of structure
designabilities. Recently, questions were raised concerning the use of a
2-letter (HP) code in such studies. Here we calculate designabilities using all
20 amino acids, with empirically determined interaction potentials (MJ matrix),
and compare with HP model results. We find good qualitative agreement between
the two models. In particular, highly designable structures in the HP model are
also highly designable in the MJ model--and vice versa--with the associated
sequences having enhanced thermodynamic stability.


  We study the free energy of the worm-like-chain model, in the
constant-extension ensemble, as a function of the stiffness for finite chains
of length L. We find that the polymer properties obtained in this ensemble are
"qualitatively" different from those obtained using constant-force ensembles.
In particular we find that as we change the stiffness parameter, the polymer
makes a transition from the flexible to the rigid phase and there is an
intermediate regime of parameter values where the free energy has three minima
and both phases are stable. This leads to interesting features in the
force-extension curves.


  We present a numerical solution of the Worm-Like Chain (WLC) model for
semi-flexible polymers. We display graphs for the end-to-end distance
distribution and the force-extension relation expected from the model. We
predict the expected level of fluctuations around the mean value in
force-extension curves. Our treatment analyses the entire range of polymer
lengths and reproduces interesting qualitative features seen in recent computer
simulations for polymers of intermediate length. These results can be tested
against experiments on single molecules. This study is relevant to mechanical
properties of biological molecules.


  We present a method for solving the Worm Like Chain (WLC) model for twisting
semiflexible polymers to any desired accuracy. We show that the WLC free energy
is a periodic function of the applied twist with period 4 pi. We develop an
analogy between WLC elasticity and the geometric phase of a spin half system.
These analogies are used to predict elastic properties of twist-storing
polymers. We graphically display the elastic response of a single molecule to
an applied torque. This study is relevant to mechanical properties of
biopolymers like DNA.


  A representation of the genetic code as a six-dimensional Boolean hypercube
is proposed. It is assumed here that this structure is the result of the
hierarchical order of the interaction energies of the bases in codon-anticodon
recognition. The proposed structure demonstrates that in the genetic code there
is a balance between conservatism and innovation. Comparing aligned positions
in homologous protein sequences two different behaviors are found: a) There are
sites in which the different amino acids present may be explained by one or two
"attractor nodes" (coding for the dominating amino acid(s)) and their one-bit
neighbors in the codon hypercube, and b) There are sites in which the amino
acids present correspond to codons located in closed paths in the hypercube.
The structure of the code facilitates evolution: the variation found at the
variable positions of proteins do not corresponds to random jumps at the codon
level, but to well defined regions of the hypercube.


  We investigate correlations in information carriers, e.g. texts and pieces of
music, which are represented by strings of letters. For information carrying
strings generated by one source (i.e. a novel or a piece of music) we find
correlations on many length scales. The word distribution, the higher order
entropies and the transinformation are calculated. The analogy to strings
generated through symbolic dynamics by nonlinear systems in critical states is
discussed.


  We investigate symbolic sequences and in particular information carriers as
e.g. books and DNA-strings. First the higher order Shannon entropies are
calculated, a characteristic root law is detected. Then the algorithmic entropy
is estimated by using Lempel-Ziv compression algorithms. In the third section
the correlation function for distant letters, the low frequency Fourier
spectrum and the characteristic scaling exponents are calculated. We show that
all these measures are able to detect long-range correlations. However, as
demonstrated by shuffling experiments, different measures operate on different
length scales. The longest correlations found in our analysis comprise a few
hundreds or thousands of letters and may be understood as long-wave
fluctuations of the composition.


  A representation of the genetic code as a six-dimensional Boolean hypercube
is described. This structure is the result of the hierarchical order of the
interaction energies of the bases in codon-anticodon recognition. In this paper
it is applied to study molecular evolution in vivo and in vitro. In the first
case we compared aligned positions in homologous protein sequences and found
two different behaviors: a) There are sites in which the different amino acids
may be explained by one or two "attractor nodes" (coding for the dominating
amino acid(s)) and their one-bit neighbors in the codon hypercube, and b) There
are sites in which the amino acids correspond to codons located in closed paths
in the hypercube. In the second case we studied the "Sexual PCR" experiment
described by Stemmer and found that the success of this combination of usual
PCR and recombination is in part due to the Gray code structure of the genetic
code.


  We study competitive DNA sequence evolution directed by {\it in vitro}
protein binding. The steady-state dynamics of this process is well described by
a shape-preserving pulse which decelerates and eventually reaches equilibrium.
We explain this dynamical behavior within a continuum mean-field framework.
Analytical results obtained on the motion of the pulse agree with simulations.
Furthermore, finite population correction to the mean-field results are found
to be insignificant.


  We carry out a theoretical study of the vibrational and relaxation properties
of naturally-occurring proteins with the purpose of characterizing both the
folding and equilibrium thermodynamics. By means of a suitable model we provide
a full characterization of the spectrum and eigenmodes of vibration at various
temperatures by merely exploiting the knowledge of the protein native
structure. It is shown that the rate at which perturbations decay at the
folding transition correlates well with experimental folding rates. This
validation is carried out on a list of about 30 two-state folders. Furthermore,
the qualitative analysis of residues mean square displacements (shown to
accurately reproduce crystallographic data) provides a reliable and
statistically accurate method to identify crucial folding sites/contacts. This
novel strategy is validated against clinical data for HIV-1 Protease. Finally,
we compare the spectra and eigenmodes of vibration of natural proteins against
randomly-generated compact structures and regular random graphs. The comparison
reveals a distinctive enhanced flexibility of natural structures accompanied by
slow relaxation times at the folding temperature. The fact that these
properties are intimately connected to the presence and assembly of secondary
motifs hints at the special criteria adopted by evolution in the selection of
viable folds.


  We address the problem of the statistical analysis of a time series generated
by complex dynamics with a new method: the Diffusion Entropy Analysis (DEA)
(Fractals, {\bf 9}, 193 (2001)). This method is based on the evaluation of the
Shannon entropy of the diffusion process generated by the time series imagined
as a physical source of fluctuations, rather than on the measurement of the
variance of this diffusion process, as done with the traditional methods. We
compare the DEA to the traditional methods of scaling detection and we prove
that the DEA is the only method that always yields the correct scaling value,
if the scaling condition applies. Furthermore, DEA detects the real scaling of
a time series without requiring any form of de-trending. We show that the joint
use of DEA and variance method allows to assess whether a time series is
characterized by L\'{e}vy or Gauss statistics. We apply the DEA to the study of
DNA sequences, and we prove that their large-time scales are characterized by
L\'{e}vy statistics, regardless of whether they are coding or non-coding
sequences. We show that the DEA is a reliable technique and, at the same time,
we use it to confirm the validity of the dynamic approach to the DNA sequences,
proposed in earlier work.


  We develop a statistical theory of networks. A network is a set of vertices
and links given by its adjacency matrix $\c$, and the relevant statistical
ensembles are defined in terms of a partition function $Z=\sum_{\c} \exp
{[}-\beta \H(\c) {]}$. The simplest cases are uncorrelated random networks such
as the well-known Erd\"os-R\'eny graphs. Here we study more general
interactions $\H(\c)$ which lead to {\em correlations}, for example, between
the connectivities of adjacent vertices. In particular, such correlations occur
in {\em optimized} networks described by partition functions in the limit
$\beta \to \infty$. They are argued to be a crucial signature of evolutionary
design in biological networks.


  We study the electrical double layer at the interface between a protein
crystal and a salt solution or a dilute solution of protein, and estimate the
double layer's contribution to the interfacial tension of this interface. This
contribution is negative and decreases in magnitude with increasing salt
concentration. We also consider briefly the interaction between a pair of
protein surfaces.


  We discuss the physics of RNA as described by its secondary structure. We
examine the static properties of a homogeneous RNA-model that includes pairing
and base stacking energies as well as entropic costs for internal loops. For
large enough costs the model exhibits a thermal denaturation transition which
we analyze in terms of the radius of gyration. We point out an inconsistency in
the standard approach to RNA secondary structure prediction for large
molecules. Under an external force a second order phase transition between a
globular and an extended phase takes place. A Harris-type criterion shows that
sequence disorder does not affect the correlation length exponent while the
other critical exponents are modified in the glass phase. However, at high
temperatures, on a coarse-grained level, disordered RNA is well described by a
homogeneous model. The characteristics of force-extension curves are discussed
as a function of the energy parameters. We show that the force transition is
always second order. A re-entrance phenomenon relevant for real disordered RNA
is predicted.


  We present a method to decompose biochemical networks into subnetworks based
on the global geometry of the network. This method enables us to analyse the
full hierarchical organisation of biochemical networks and is applied to 43
organisms from the WIT database. Two types of biochemical networks are
considered: metabolic networks and whole-cellular networks (also including e.g.
information processes). Conceptual and quantitative ways of describing the
hierarchical ordering are discussed. The general picture of the metabolic
networks arising from our study is that of a few core-clusters centred around
the most highly connected substances enclosed by other substances in outer
shells, and a few other well-defined subnetworks.


  We study the force-induced unfolding of random disordered RNA or
single-stranded DNA polymers. The system undergoes a second order phase
transition from a collapsed globular phase at low forces to an extensive
necklace phase with a macroscopic end-to-end distance at high forces. At low
temperatures, the sequence inhomogeneities modify the critical behaviour. We
provide numerical evidence for the universality of the critical exponents
which, by extrapolation of the scaling laws to zero force, contain useful
information on the ground state (f=0) properties. This provides a good method
for quantitative studies of scaling exponents characterizing the collapsed
globule. In order to get rid of the blurring effect of thermal fluctuations we
restrict ourselves to the groundstate at fixed external force. We analyze the
statistics of rearrangements, in particular below the critical force, and point
out its implications for force-extension experiments on single molecules.


  A typical protein structure is a compact packing of connected alpha-helices
and/or beta-strands. We have developed a method for generating the ensemble of
compact structures a given set of helices and strands can form. The method is
tested on structures composed of four alpha-helices connected by short turns.
All such natural four-helix bundles that are connected by short turns seen in
nature are reproduced to closer than 3.6 Angstroms per residue within the
ensemble. Since structures with no natural counterpart may be targets for ab
initio structure design, the designability of each structure in the ensemble --
defined as the number of sequences with that structure as their lowest energy
state -- is evaluated using a hydrophobic energy. For the case of four
alpha-helices, a small set of highly designable structures emerges, most of
which have an analog among the known four-helix fold families, however several
novel packings and topologies are identified.


  We study the kinetics of protein folding via statistical energy landscape
theory. We concentrate on the local-connectivity case, where the
configurational changes can only occur among neighboring states, with the
folding progress described in terms of an order parameter given by the fraction
of native conformations. The non-Markovian diffusion dynamics is analyzed in
detail and an expression for the mean first-passage time (MFPT) from non-native
unfolded states to native folded state is obtained. It was found that the MFPT
has a V-shaped dependence on the temperature. We also find that the MFPT is
shortened as one increases the gap between the energy of the native and average
non-native folded states relative to the fluctuations of the energy landscape.
The second- and higher-order moments are studied to infer the first-passage
time (FPT) distribution. At high temperature, the distribution becomes close to
a Poisson distribution, while at low temperatures the distribution becomes a
L\'evy-like distribution with power-law tails, indicating a non-self-averaging
intermittent behavior of folding dynamics. We note the likely relevance of this
result to single-molecule dynamics experiments, where a power law (L\'evy)
distribution of the relaxation time of the underlined protein energy landscape
is observed.


  George and Wilson [Acta. Cryst. D 50, 361 (1994)] looked at the distribution
of values of the second virial coefficient of globular proteins, under the
conditions at which they crystallise. They found the values to lie within a
fairly narrow range. We have defined a simple model of a generic globular
protein. We then generate a set of proteins by picking values for the
parameters of the model from a probability distribution. At fixed solubility,
this set of proteins is found to have values of the second virial coefficient
that fall within a fairly narrow range. The shape of the probability
distribution of the second virial coefficient is Gaussian because the second
virial coefficient is a sum of contributions from different patches on the
protein surface.


  Motivated by the problem of the definition of a distance between two
sequences of characters, we investigate the so-called learning process of
typical sequential data compression schemes. We focus on the problem of how a
compression algorithm optimizes its features at the interface between two
different sequences $A$ and $B$ while zipping the sequence $A+B$ obtained by
simply appending $B$ after $A$. We show the existence of a universal scaling
function (the so-called learning function) which rules the way in which the
compression algorithm learns a sequence $B$ after having compressed a sequence
$A$. In particular it turns out that it exists a crossover length for the
sequence $B$, which depends on the relative entropy between $A$ and $B$, below
which the compression algorithm does not learn the sequence $B$ (measuring in
this way the relative entropy between $A$ and $B$) and above which it starts
learning $B$, i.e. optimizing the compression using the specific features of
$B$. We check the scaling function on three main classes of systems: Bernoulli
schemes, Markovian sequences and the symbolic dynamic generated by a non
trivial chaotic system (the Lozi map). As a last application of the method we
present the results of a recognition experiment, namely recognize which
dynamical systems produced a given time sequence. We finally point out the
potentiality of these results for segmentation purposes, i.e. the
identification of homogeneous sub-sequences in heterogeneous sequences (with
applications in various fields from genetic to time-series analysis).


  We construct a gene network based on expression data from DNA microarray
experiments, by establishing a link between two genes whenever the Pearson's
correlation coefficient between their expression profiles is higher than a
certain cutoff. The resulting connectivity distribution is compatible with a
power-law decay with exponent ~1, corrected by an exponential cutoff at large
connectivity. The biological relevance of such network is demonstrated by
showing that there is a strong statistical correlation between high
connectivity number and lethality: in close analogy to what happens for protein
interaction networks, essential genes are strongly overerpresented among the
hubs of the network, that is the genes connected to many other genes.


  We present an effective immunization strategy for computer networks and
populations with broad and, in particular, scale-free degree distributions. The
proposed strategy, acquaintance immunization, calls for the immunization of
random acquaintances of random nodes (individuals). The strategy requires no
knowledge of the node degrees or any other global knowledge, as do targeted
immunization strategies. We study analytically the critical threshold for
complete immunization. We also study the strategy with respect to the
susceptible-infected-removed epidemiological model. We show that the
immunization threshold is dramatically reduced with the suggested strategy, for
all studied cases.


  The response of a short DNA segment to bending is studied, taking into
account the anisotropy in the bending rigidities caused by the double-helical
structure. It is shown that the anisotropy introduces an effective nonlinear
twist-bend coupling that can lead to the formation of kinks and modulations in
the curvature and/or in the twist, depending on the values of the elastic
constants and the imposed deflection angle. The typical wavelength for the
modulations, or the distance between the neighboring kinks is found to be set
by half of the DNA pitch.


  Molecular dynamics simulations in simplified models allow one to study the
scaling properties of folding times for many proteins together under a
controlled setting. We consider three variants of the Go models with different
contact potentials and demonstrate scaling described by power laws and no
correlation with the relative contact order parameter. We demonstrate existence
of at least three kinetic universality classes which are correlated with the
types of structure: the alpha-, alpha--beta-, and beta- proteins have the
scaling exponents of about 1.7, 2.5, and 3.2 respectively. The three classes
merge into one when the contact range is truncated at a 'reasonable' value. We
elucidate the role of the potential associated with the chirality of a protein.


  A lattice of uniformly charged, infinitesimally thin, rods decorated with an
ordered array of counterions exhibits anomalous behavior as the spacing between
the rods is varied. In particular, the counterion lattice undergoes a sequence
of structural shearing, or "tilting," phase transformations as the spacing
between the rods decreases. The potential implications of this behavior with
respect to the packaging of biologically relevant polyelectrolytic molecules
are commented upon.


  Charge transfer along the base-pair stack in DNA is modeled in terms of
thermally-assisted tunneling between adjacent base pairs. Central to our
approach is the notion that tunneling between fluctuating pairs is rate-limited
by the requirement of their optimal alignment. We focus on this aspect of the
process by modeling two adjacent base pairs in terms of a classical damped
oscillator subject to thermal fluctuations as described by a Fokker-Planck
equation. We find that the process is characterized by two time scales, a
result that is in accord with experimental findings.


  The structure of molecular networks derives from dynamical processes on
evolutionary time scales. For protein interaction networks, global statistical
features of their structure can now be inferred consistently from several
large-throughput datasets. Understanding the underlying evolutionary dynamics
is crucial for discerning random parts of the network from biologically
important properties shaped by natural selection. We present a detailed
statistical analysis of the protein interactions in Saccharomyces cerevisiae
based on several large-throughput datasets. Protein pairs resulting from gene
duplications are used as tracers into the evolutionary past of the network.
  From this analysis, we infer rate estimates for two key evolutionary
processes shaping the network: (i) gene duplications and (ii) gain and loss of
interactions through mutations in existing proteins, which are referred to as
link dynamics. Importantly, the link dynamics is asymmetric, i.e., the
evolutionary steps are mutations in just one of the binding parters. The link
turnover is shown to be much faster than gene duplications. According to this
model, the link dynamics is the dominant evolutionary force shaping the
statistical structure of the network, while the slower gene duplication
dynamics mainly affects its size. Specifically, the model predicts (i) a broad
distribution of the connectivities (i.e., the number of binding partners of a
protein) and (ii) correlations between the connectivities of interacting
proteins.


  We study theoretically the denaturation of single RNA molecules by mechanical
stretching, focusing on signatures of the (un)folding pathway in molecular
fluctuations. Our model describes the interactions between nucleotides by
incorporating the experimentally determined free energy rules for RNA secondary
structure, while exterior single stranded regions are modeled as freely jointed
chains. For exemplary RNA sequences (hairpins and the Tetrahymena thermophila
group I intron), we compute the quasi-equilibrium fluctuations in the
end-to-end distance as the molecule is unfolded by pulling on opposite ends.
Unlike the average quasi-equilibrium force-extension curves, these fluctuations
reveal clear signatures from the unfolding of individual structural elements.
We find that the resolution of these signatures depends on the spring constant
of the force-measuring device, with an optimal value intermediate between very
rigid and very soft. We compare and relate our results to recent experiments by
Liphardt et al. [Science 292, 733-737 (2001)].


  Neutral evolution is the simplest model of molecular evolution and thus it is
most amenable to a comprehensive theoretical investigation. In this paper, we
characterize the statistical properties of neutral evolution of proteins under
the requirement that the native state remains thermodynamically stable, and
compare them to the ones of Kimura's model of neutral evolution. Our study is
based on the Structurally Constrained Neutral (SCN) model which we recently
proposed. We show that, in the SCN model, the substitution rate decreases as
longer time intervals are considered, and fluctuates strongly from one branch
of the evolutionary tree to another, leading to a non-Poissonian statistics for
the substitution process. Such strong fluctuations are also due to the fact
that neutral substitution rates for individual residues are strongly correlated
for most residue pairs. Interestingly, structurally conserved residues,
characterized by a much below average substitution rate, are also much less
correlated to other residues and evolve in a much more regular way. Our results
could improve methods aimed at distinguishing between neutral and adaptive
substitutions as well as methods for computing the expected number of
substitutions occurred since the divergence of two protein sequences.


  We discuss a population of sequences subject to mutations and
frequency-dependent selection, where the fitness of a sequence depends on the
composition of the entire population. This type of dynamics is crucial to
understand the evolution of genomic regulation. Mathematically, it takes the
form of a reaction-diffusion problem that is nonlinear in the population state.
In our model system, the fitness is determined by a simple mathematical game,
the hawk-dove game. The stationary population distribution is found to be a
quasispecies with properties different from those which hold in fixed fitness
landscapes.


  Bacteria typically have a few thousand different proteins. The number of
proteins with a given charge is a roughly Gaussian function of charge - centred
near zero, and with a width around ten (in units of the charge on the proton).
We have used the charges on E. coli's proteins to estimate the changes in the
second virial coefficients of all its proteins as the concentration of a 1:1
salt is increased. The second virial coefficient has dimensions of volume and
we find that on average it decreases by about twice the average volume of a
protein when the salt concentration is increased from 0.2 to 1 Molar. The
standard deviation of the decrease is of the same order. The consequences of
this for the complex mixture of proteins inside an E. coli cell, are briefly
discussed.


  An exactly solvable model based on the topology of a protein native state is
applied to identify bottlenecks and key-sites for the folding of HIV-1
Protease. The predicted sites are found to correlate well with clinical data on
resistance to FDA-approved drugs. It has been observed that the effects of drug
therapy are to induce multiple mutations on the protease. The sites where such
mutations occur correlate well with those involved in folding bottlenecks
identified through the deterministic procedure proposed in this study. The high
statistical significance of the observed correlations suggests that the
approach may be promisingly used in conjunction with traditional techniques to
identify candidate locations for drug attacks.


  Following recent X-ray diffraction experiments by Wong, Li, and Safinya on
biopolymer gels, we apply Onsager excluded volume theory to a nematic mixture
of rigid rods and strong ``$\pi/2$'' cross-linkers obtaining a long-ranged,
highly anisotropic depletion attraction between the linkers. This attraction
leads to breakdown of the percolation theory for this class of gels, to
breakdown of Onsager's second-order virial method, and to formation of
heterogeneities in the form of raft-like ribbons.


  We simulate neutral evolution of proteins imposing conservation of the
thermodynamic stability of the native state in the framework of an effective
model of folding thermodynamics. This procedure generates evolutionary
trajectories in sequence space which share two universal features for all of
the examined proteins. First, the number of neutral mutations fluctuates
broadly from one sequence to another, leading to a non-Poissonian substitution
process. Second, the number of neutral mutations displays strong correlations
along the trajectory, thus causing the breakdown of self-averaging of the
resulting evolutionary substitution process.


  The electronic structure of an infinite poly-guanine - poly-cytosine DNA
molecule in its dry A-helix structure is studied by means of density-functional
calculations. An extensive study of 30 nucleic base pairs is performed to
validate the method. The electronic energy bands of DNA close to the Fermi
level are then analyzed in order to clarify the electron transport properties
in this particularly simple DNA realization, probably the best suited candidate
for conduction. The energy scale found for the relevant band widths, as
compared with the energy fluctuations of vibrational or genetic-sequence
origin, makes highly implausible the coherent transport of electrons in this
system. The possibility of diffusive transport with sub-nanometer mean free
paths is, however, still open. Information for model Hamiltonians for
conduction is provided.


  $\alpha$-helices stand out as common and relatively invariant secondary
structural elements of proteins. However, $\alpha$-helices are not rigid bodies
and their deformations can be significant in protein function ({\it e.g.}
coiled coils). To quantify the flexibility of $\alpha$-helices we have
performed a structural principal-component analysis of helices of different
lengths from a representative set of protein folds in the Protein Data Bank. We
find three dominant modes of flexibility: two degenerate bend modes and one
twist mode. The data are consistent with independent Gaussian distributions for
each mode. The mode eigenvalues, which measure flexibility, follow simple
scaling forms as a function of helix length. The dominant bend and twist modes
and their harmonics are reproduced by a simple spring model, which incorporates
hydrogen-bonding and excluded volume. As an application, we examine the amount
of bend and twist in helices making up several coiled-coil proteins.
Incorporation of $\alpha$-helix flexibility into structure refinement and
design is discussed.