  Based on a recent model of evolving viruses competing with an adapting immune
system [1], we study the conditions under which a viral quasispecies can
maximize its growth rate. The range of mutation rates that allows viruses to
thrive is limited from above due to genomic information deterioration, and from
below by insufficient sequence diversity, which leads to a quick eradication of
the virus by the immune system. The mutation rate that optimally balances these
two requirements depends to first order on the ratio of the inverse of the
virus' growth rate and the time the immune system needs to develop a specific
answer to an antigen. We find that a virus is most viable if it generates
exactly one mutation within the time it takes for the immune system to adapt to
a new viral epitope. Experimental viral mutation rates, in particular for HIV
(human immunodeficiency virus), seem to suggest that many viruses have achieved
their optimal mutation rate. [1] C.Kamp and S. Bornholdt, Phys. Rev. Lett., 88,
068104 (2002)


  We discuss the long-run behavior of stochastic dynamics of many interacting
players in spatial evolutionary games. In particular, we investigate the effect
of the number of players and the noise level on the stochastic stability of
Nash equilibria. We discuss similarities and differences between systems of
interacting players maximizing their individual payoffs and particles
minimizing their interaction energy. We use concepts and techniques of
statistical mechanics to study game-theoretic models. In order to obtain
results in the case of the so-called potential games, we analyze the
thermodynamic limit of the appropriate models of interacting particles.


  Diffusion processes in biological membranes are of interest to understand the
macromolecular organisation and function of several molecules. Fluorescence
Recovery After Photobleaching (FRAP) has been widely used as a method to
analyse this processes using classical Brownian diffusion model. In the first
part of this work, the analytical expression of the fluorescence recovery as a
function of time has been established for anomalous diffusion due to long
waiting times. Then, experimental fluorescence recoveries recorded in living
cells on a membrane-bound protein have been analysed using three different
models : normal Brownian diffusion, Brownian diffusion with an immobile
fraction and anomalous diffusion due to long waiting times.


  Chromatin is a complex of DNA and specific proteins forming an intermediary
level of organization of eukaryotic genomes, between double-stranded DNA and
chromosome. Within a generic modeling of the chromatin assembly, we investigate
the interplay between the mechanical properties of the chromatin fiber and its
biological functions. A quantitative step is to relate the mechanics at the DNA
level and the mechanics described at the chromatin fiber level. It allows to
calculate the complete set of chromatin elastic constants (twist and bend
persistence lengths, stretch modulus and twist-stretch coupling constant), in
terms of DNA elastic properties and geometric features of the fiber. These
elastic constants are strongly sensitive to the local architecture of the fiber
and we argue that this tunable elasticity might be a key feature in chromatin
functions, for instance in the initiation and regulation of transcription.
Moreover, this analysis provides a framework to interpret micromanipulations
studies of chromatin fiber and suggests further experiments involving
intercalators to scan the tunable elasticity of the fiber.


  The chromatin fiber is a complex of DNA and specific proteins called histones
forming the first structural level of organization of eukaryotic chromosomes.
In tightly organized chromatin fibers, the short segments of naked DNA linking
the nucleosomes are strongly end constrained. Longitudinal thermal fluctuations
in these linkers allow intercalative mode of protein binding. We show that
mechanical constraints generated in the first stage of the binding process
induce linker DNA buckling; buckling in turn modifies the binding energies and
activation barriers and creates a force of decondensation at the chromatin
fiber level. The unique structure and properties of DNA thus yield a novel
physical mechanism of buckling instability that might play a key role in the
regulation of gene expression.


  The effective free energy of globular protein chain is considered to be a
functional defined on smooth curves in three dimensional Euclidean space. From
the requirement of geometrical invariance, together with basic facts on
conformation of helical proteins and dynamical characteristics of the protein
chains, we are able to determine, in a unique way, the exact form of the free
energy functional. Namely, the free energy density should be a linear function
of the curvature of curves on which the free energy functional is defined. We
briefly discuss the possibility of using the model proposed in Monte Carlo
simulations of exhaustive searching the native stable state of the protein
chain. The relation of this model to the rigid relativistic particles and
strings is also considered.


  Vesicles self-assembled from amphiphilic diblock copolymers exhibit a wide
diversity of behavior upon electroporation, due to competitions between edge,
surface and bending energies that drive the system, while different viscous
dissipation mechanisms determine the timescales. These effects are
manifestations of the varying membrane thickness d from what are essentially
chemically identical molecules. For smaller d, we find large unstable pores and
the resulting membrane fragments reassemble into vesicles within minutes.
Vesicles with larger d form long-lived pores of nanometer dimension that could
be potentially useful for sieving applications and drug delivery.


  A model for soliton dynamics on a hydrogen-bond network in helical proteins
is proposed. It employs in three dimensions the formalism of fully integrable
Toda lattices which admits phonons as well as solitons along the hydrogen-bonds
of the helices. A simulation of the three dimensional Toda lattice system shows
that the solitons are spontaneously created and are stable and moving along the
helix axis. A perturbation on one of the three H-bond lines forms solitons on
the other H-bonds as well. The robust solitary wave may explain very long-lived
modes in the frequency range of 100 cm$^{-1}$ which are found in recent X-ray
laser experiments. The dynamics parameters of the Toda lattice are in
accordance with the usual Lennard-Jones parameters used for realistic H-bond
potentials in proteins.


  Random Boolean networks, originally invented as models of genetic regulatory
networks, are simple models for a broad class of complex systems that show rich
dynamical structures. From a biological perspective, the most interesting
networks lie at or near a critical point in parameter space that divides
``ordered'' from ``chaotic'' attractor dynamics. In the ordered regime, we show
rigorously that the average number of relevant nodes (the ones that determine
the attractor dynamics) remains constant with increasing system size N. For
critical networks, our analysis and numerical results show that the number of
relevant nodes scales like N^{1/3}. Numerical experiments also show that the
median number of attractors in critical networks grows faster than linearly
with N. The calculations explain why the correct asymptotic scaling is observed
only for very large N.


  We find exact solutions to a new one-dimensional (1D) interacting particle
theory and apply the results to the adsorption and wrapping of polymers (such
as DNA) around protein particles (such as histones). Each adsorbed protein is
represented by a Tonks gas particle. The length of each particle is a degree of
freedom that represents the degree of DNA wrapping around each histone.
Thermodynamic quantities are computed as functions of wrapping energy, adsorbed
histone density, and bulk histone concentration (or chemical potential); their
experimental signatures are also discussed. Histone density is found to undergo
a two-stage adsorption process as a function of chemical potential, while the
mean coverage by high affinity proteins exhibits a maximum as a function of the
chemical potential. However, {\it fluctuations} in the coverage are
concurrently maximal. Histone-histone correlation functions are also computed
and exhibit rich two length scale behavior.


  Myosin-V is a motor protein responsible for organelle and vesicle transport
in cells. Recent single-molecule experiments have shown that it is an efficient
processive motor that walks along actin filaments taking steps of mean size
close to 36 nm. A theoretical study of myosin-V motility is presented following
an approach used successfully to analyze the dynamics of conventional kinesin
but also taking some account of step-size variations. Much of the present
experimental data for myosin-V can be well described by a two-state chemical
kinetic model with three load-dependent rates. In addition, the analysis
predicts the variation of the mean velocity and of the randomness -- a
quantitative measure of the stochastic deviations from uniform, constant-speed
motion -- with ATP concentration under both resisting and assisting loads, and
indicates a {\it sub}step of size $d_{0} \simeq$ 13-14 nm (from the ATP-binding
site) that appears to accord with independent observations.


  Temporal gene expression data (Wen, et. al.) was analyzed using the recently
introduced inverse modeling technique of Embedded Complex Logistic Maps (ECLM).
Preliminary results indicate scale-free structure in the gene regulatory
network topology as well as in the network's energy dissipation
characteristics. General agreement was found with other recent studies
involving a dynamical model and analyses of other biochemical networks
exhibiting scale-free topology. The highlights of our study are briefly
outlined.


  We study selection acting on phenotype in a collection of agents playing
local games lacking Nash equilibria. After each cycle one of the agents losing
most games is replaced by a new agent with new random strategy and game
partner. The network generated can be considered critical in the sense that the
lifetimes of the agents is power law distributed. The longest surviving agents
are those with the lowest absolute score per time step. The emergent ecology is
characterized by a broad range of behaviors. Nevertheless, the agents tend to
be similar to their opponents in terms of performance.


  Cells receive a wide variety of cellular and environmental signals, which
must be processed combinatorially to generate specific and timely genetic
responses. We present here a theoretical study on the combinatorial control and
integration of transcription signals, with the finding that cis-regulatory
systems with specific protein-DNA interaction and glue-like protein-protein
interactions, supplemented by distal activation or repression mechanisms, have
the capability to execute a wide range of control functions encoded in the
regulatory DNA sequences. Using a quantitative model based on the
well-characterized bacterial transcription system, we show explicitly how
various regulatory logic functions can be implemented, by selecting the
strengths and relative positions of the relevant protein-binding DNA sequences
in the cis-regulatory region. The architecture that emerges is naturally
modular and highly evolvable. Our findings suggest that the transcription
regulatory apparatus is a "programmable" computing machine, belonging formally
to the class of Boltzmann machines. We also expose critical shortcomings of the
bacterial transcription systems, which limit the genome-wide adoption of
schemes for complex transcription control, and discuss how they may be overcome
by the eukaryotic transcription systems.


  Gene expression has a stochastic component owing to the single molecule
nature of the gene and the small number of copies of individual DNA binding
proteins in the cell. We show how the statistics of such systems can be mapped
on to quantum many body problems. The dynamics of a single gene switch
resembles the spin boson model of a two site polaron or an electron transfer
reaction. Networks of switches can be approximately described as quantum spin
systems by using an appropriate variational principle. In this way the concept
of frustration for magnetic systems can be taken over into gene networks. The
landscape of stable attractors depends on the degree and style of frustration,
much as for neural networks. We show the number of attractors, which may
represent cell types, is much smaller for appropriately designed, weakly
frustrated stochastic networks than for randomly connected networks.


  We study theoretically the in vitro evolution of a DNA sequence by binding to
a transcription factor. Using a simple model of protein-DNA binding and
available binding constants for the Mnt protein, we perform large-scale,
realistic simulations of evolution starting from a single DNA sequence. We
identify different parameter regimes characterized by distinct evolutionary
behaviors. For each regime we find analytical estimates which agree well with
simulation results. For small population sizes, the DNA evolutional path is a
random walk on a smooth landscape. While for large population sizes, the
evolution dynamics can be well described by a mean-field theory. We also study
how the details of the DNA-protein interaction affect the evolution.


  Spatial models for spread of an epidemic may be mapped onto bond percolation.
We point out that with disorder in the strength of contacts between individuals
patchiness in the spread of the epidemic is very likely, and the criterion for
epidemic outbreak depends strongly on the disorder because the critical region
of the corresponding percolation model is broadened. In some networks the
percolation threshold is zero if another kind of disorder is present, namely
divergent fluctuations in the number of contacts. We give an example, a network
with a well defined geography, where this is not necessarily so, and discuss
whether real infection networks are likely to have this property.


  We investigate the equilibrium state of the model of Peng, \textit{et al.}
for molecular breeding. In the model, a population of DNA sequences is
successively culled by removing the sequences with the lowest binding affinity
to a particular target sequence. The remaining sequences are then amplified to
restore the original population size, undergoing some degree of
point-substitution of nucleotides in the process. Working in the infinite
population size limit, we derive an equation for the equilibrium distribution
of binding affinity, here modeled by the number of matches to the target
sequence. The equation is then solved approximately in the limit of large
sequence length, in the three regimes of strong, intermediate and weak
selection. The approximate solutions are verified via comparison to exact
numerical results.


  Differences in the regional substitution patterns in the human genome created
patterns of large-scale variation of base composition known as genomic
isochores. To gain insight into the origin of the genomic isochores we develop
a maximum likelihood approach to determine the history of substitution patterns
in the human genome. This approach utilizes the vast amount of repetitive
sequence deposited in the human genome over the past ~250 MYR. Using this
approach we estimate the frequencies of seven types of substitutions: the four
transversions, two transitions, and the methyl-assisted transition of cytosine
in CpG. Comparing substitutional patterns in repetitive elements of various
ages, we reconstruct the history of the base-substitutional process in the
different isochores for the past 250 Myr. At around 90 Myr ago (around the time
of the mammalian radiation), we find an abrupt 4- to 8-fold increase of the
cytosine transition rate in CpG pairs compared to that of the reptilian
ancestor. Further analysis of nucleotide substitutions in regions with
different GC-content reveals concurrent changes in the substitutional patterns.
While the substitutional pattern was dependent on the regional GC-content in
such ways that it preserved the regional GC-content before the mammalian
radiation, it lost this dependence afterwards. The substitutional pattern
changed from an isochore-preserving to an isochore-degrading one. We conclude
that isochores have been established before the radiation of the eutherian
mammals and have been subject to the process of homogenization since then.


  The regulation of a gene depends on the binding of transcription factors to
specific sites located in the regulatory region of the gene. The generation of
these binding sites and of cooperativity between them are essential building
blocks in the evolution of complex regulatory networks. We study a theoretical
model for the sequence evolution of binding sites by point mutations. The
approach is based on biophysical models for the binding of transcription
factors to DNA. Hence we derive empirically grounded fitness landscapes, which
enter a population genetics model including mutations, genetic drift, and
selection. We show that the selection for factor binding generically leads to
specific correlations between nucleotide frequencies at different positions of
a binding site. We demonstrate the possibility of rapid adaptive evolution
generating a new binding site for a given transcription factor by point
mutations. The evolutionary time required is estimated in terms of the neutral
(background) mutation rate, the selection coefficient, and the effective
population size. The efficiency of binding site formation is seen to depend on
two joint conditions: the binding site motif must be short enough and the
promoter region must be long enough. These constraints on promoter architecture
are indeed seen in eukaryotic systems. Furthermore, we analyse the adaptive
evolution of genetic switches and of signal integration through binding
cooperativity between different sites. Experimental tests of this picture
involving the statistics of polymorphisms and phylogenies of sites are
discussed.


  A protocell model consisting of mutually catalyzing molecules is studied in
order to investigate how chemical compositions are transferred recursively
through cell divisions under replication errors. Depending on the path rate,
the numbers of molecules and species, three phases are found: fast switching
state without recursive production, recursive production, and itinerancy
between the above two states. The number distributions of the molecules in the
recursive states are shown to be log-normal except for those species that form
a core hypercycle, and are explained with the help of a heuristic argument.


  Understanding the subgraph distribution in random networks is important for
modelling complex systems. In classic Erdos networks, which exhibit a
Poissonian degree distribution, the number of appearances of a subgraph G with
n nodes and g edges scales with network size as \mean{G} ~ N^{n-g}. However,
many natural networks have a non-Poissonian degree distribution. Here we
present approximate equations for the average number of subgraphs in an
ensemble of random sparse directed networks, characterized by an arbitrary
degree sequence. We find new scaling rules for the commonly occurring case of
directed scale-free networks, in which the outgoing degree distribution scales
as P(k) ~ k^{-\gamma}. Considering the power exponent of the degree
distribution, \gamma, as a control parameter, we show that random networks
exhibit transitions between three regimes. In each regime the subgraph number
of appearances follows a different scaling law, \mean{G} ~ N^{\alpha}, where
\alpha=n-g+s-1 for \gamma<2, \alpha=n-g+s+1-\gamma for 2<\gamma<\gamma_c, and
\alpha=n-g for \gamma>\gamma_c, s is the maximal outdegree in the subgraph, and
\gamma_c=s+1. We find that certain subgraphs appear much more frequently than
in Erdos networks. These results are in very good agreement with numerical
simulations. This has implications for detecting network motifs, subgraphs that
occur in natural networks significantly more than in their randomized
counterparts.


  We examine fluctuation effects due to the low copy number of proteins
involved in pattern-forming dynamics within a bacterium. We focus on a
stochastic model of the oscillating MinCDE protein system regulating accurate
cell division in E. coli. We find that, for some parameter regions, the protein
concentrations are low enough that fluctuations are essential for the
generation of patterns. We also examine the role of fluctuations in
constraining protein concentration levels.


  We study spatial networks constructed by randomly placing nodes on a manifold
and joining two nodes with an edge whenever their distance is less than a
certain cutoff. We derive the general expression for the connectivity
distribution of such networks as a functional of the distribution of the nodes.
We show that for regular spatial densities, the corresponding spatial network
has a connectivity distribution decreasing faster than an exponential. In
contrast, we also show that scale-free networks with a power law decreasing
connectivity distribution are obtained when a certain information measure of
the node distribution (integral of higher powers of the distribution) diverges.
We illustrate our results on a simple example for which we present simulation
results. Finally, we speculate on the role played by the limiting case P(k)=1/k
which appears empirically to be relevant to spatial networks of biological
origin such as the ones constructed from gene expression data.


  A theoretical model for the folding of proteins containing disulfide bonds is
introduced. The model exploits the knowledge of the native state to favour the
progressive establishment of native interactions. At variance with traditional
approaches based on native topology, not all native bonds are treated in the
same way; in particular, a suitable energy term is introduced to account for
the special strength of disulfide bonds (irrespective of whether they are
native or not) as well as their ability to undergo intra-molecular reshuffling.
The model thus possesses the minimal ingredients necessary to investigated the
much debated issue of whether the re-folding process occurs through partially
structured intermediates with native or non-native disulfide bonds. This
strategy is applied to a context of particular interest, the re-folding process
of Hirudin, a thrombin-specific protease inhibitor, for which conflicting
folding pathways have been proposed. We show that the only two parameters in
the model (temperature and disulfide strength) can be tuned to reproduce well a
set of experimental transitions between species with different number of formed
disulfide. This model is then used to provide a characterisation of the folding
process and a detailed description of the species involved in the rate-limiting
step of Hirudin refolding.


  The response of a semiflexible polyelectrolyte chain to stretching in the
regimes of moderate and weak screening is studied theoretically, with a special
focus on DNA experiments. By using the nonlinear Poisson--Boltzmann description
of electrostatic self--interactions of the chain, we explicitly demonstrate the
applicability of the concept of effective charge to certain aspects of the
problem. This charge can be extracted from the far--field asymptotic behavior
of electrostatic potential of the fully aligned chain. Surprisingly, in terms
of the effective charge, the electrostatically renormalized persistence length
can be formally described by the classical Odijk-Skolnick-Fixman (OSF) formula,
whose domain of applicability is limited to linearized Debye Huckel (DH)
approximation. However, the short--length behavior of the chain in the
nonlinear regime deviates from the of DH--based result, even upon the charge
renormalization. This difference is revealed in calculated stretching curves of
strongly charged DNA, and our results are in good agreement with the recent
experiments. In the limit of weak screening we predict the elastic response to
have a distinctive two-stage character, with a peculiar intermediate
"unstretchable" regime.


  We propose a lattice model for RNA based on a self-interacting two-tolerant
trail. Self-avoidance and elements of tertiary structure are taken into
account. We investigate a simple version of the model in which the native state
of RNA consists of just one hairpin. Using exact arguments and Monte Carlo
simulations we determine the phase diagram for this case. We show that the
denaturation transition is first order and can either occur directly or through
an intermediate molten phase.


  Recent progress has been made in the understanding of the physical properties
of chromatin -- the dense complex of DNA and histone proteins that occupies the
nuclei of plant and animal cells. Here I will focus on the two lowest levels of
the hierarchy of DNA folding into the chromatin complex: (i) the nucleosome,
the chromatin repeating unit consisting of a globular aggregate of eight
histone proteins with the DNA wrapped around: its overcharging, the DNA
unwrapping transition, the ''sliding'' of the octamer along the DNA. (ii) The
30nm chromatin fiber, the necklace-like structure of nucleosomes connected via
linker DNA: its geometry, its mechanical properties under stretching and its
response to changing ionic conditions. I will stress that chromatin combines
two seemingly contradictory features: (1) high compaction of DNA within the
nuclear envelope and at the same time (2) accessibility to genes, promoter
regions and gene regulatory sequences.


  In this paper, we study the large-scale protein interaction network of yeast
uti lizing a stochastic method based upon percolation of random graphs. In
order to find the global features of connectivities in the network, we
introduce numeric al measures that quantify (1) how strongly a protein ties
with the other parts o f the network and (2) how significantly an interaction
contributes to the integr ity of the network. Our study shows that the
distribution of essential proteins is distinct from the background in terms of
global connectivities. This observ ation highlights a fundamental difference
between the essential and the non-esse ntial proteins in the network.
Furthermore, we find that the interaction data o btained from different
experimental methods such as immunoprecipitation and two- hybrid techniques
possess different characteristics. We discuss the biological implications of
these observations.


  The Saccharomyces cerevisiae protein-protein interaction map, as well as many
natural and man-made networks, shares the scale-free topology. The preferential
attachment model was suggested as a generic network evolution model that yields
this universal topology. However, it is not clear that the model assumptions
hold for the protein interaction network. Using a cross genome comparison we
show that (a) the older a protein, the better connected it is, and (b) The
number of interactions a protein gains during its evolution is proportional to
its connectivity. Therefore, preferential attachment governs the protein
network evolution. The evolutionary mechanism leading to such preference and
some implications are discussed.


  Networks of filamentous proteins play a crucial role in cell mechanics. These
cytoskeletal networks, together with various crosslinking and other associated
proteins largely determine the (visco)elastic response of cells. In this letter
we study a model system of crosslinked, stiff filaments in order to explore the
connection between the microstructure under strain and the macroscopic response
of cytoskeletal networks. We find two distinct regimes as a function primarily
of crosslink density and filament rigidity: one characterized by affine
deformation and one by non-affine deformation. We characterize the crossover
between these two.


  Only about 1,000 qualitatively different protein folds are believed to exist
in nature. Here, we review theoretical studies which suggest that some folds
are intrinsically more designable than others, {\it i.e.} are lowest energy
states of an unusually large number of sequences. The sequences associated with
these folds are also found to be unusually thermally stable. The connection
between highly designable structures and highly stable sequences is generally
known as the "designability principle". The designability principle may help
explain the small number of natural folds, and may also guide the design of new
folds.


  The dynamical response of a lipid membrane to a local perturbation of its
molecular symmetry is investigated theoretically. A density asymmetry between
the two membrane leaflets is predominantly released by in-plane lipid diffusion
or membrane curvature, depending upon the spatial extent of the perturbation.
It may result in the formation of non-equilibrium structures (buds), for which
a dynamical size selection is observed. A preferred size in the micrometer
range is predicted, as a signature of the crossover between membrane and
solvent dominated dynamical membrane response.


  The budding yeast {\it Saccharomyces cerevisiae} is the first eukaryote whose
genome has been completely sequenced. It is also the first eukaryotic cell
whose proteome (the set of all proteins) and interactome (the network of all
mutual interactions between proteins) has been analyzed. In this paper we study
the structure of the yeast protein complex network in which weighted edges
between complexes represent the number of shared proteins. It is found that the
network of protein complexes is a small world network with scale free behavior
for many of its distributions. However we find that there are no strong
correlations between the weights and degrees of neighboring complexes. To
reveal non-random features of the network we also compare it with a null model
in which the complexes randomly select their proteins. Finally we propose a
simple evolutionary model based on duplication and divergence of proteins.


  We present a theoretical study of single-stranded DNA under stretching.
Within the proposed framework, the effects of basepairing on the mechanical
response of the molecule can be studied in combination with an arbitrary
underlying model of chain elasticity. In a generic case, we show that the
stretching curve of ssDNA exhibits two distinct features: the second-order
"unfolding" phase transition, and a sharp crossover, reminiscent of the
first-order "unzipping" transition in dsDNA. We apply the theory to the
particular cases of Worm-like Chain (WLC) and Freely-Joint Chain (FJC) models,
and discuss the universal and model--dependent features of the mechanical
response of ssDNA. In particular, we show that variation of the width of the
unzipping crossover with interaction strength is very sensitive to the
energetics of hairpin loops. This opens a new way of testing the elastic
properties of ssDNA.


  We consider a stochastic SIS infection model for a population partitioned
into $m$ households assuming random mixing. We solve the model in the limit $m
\to \infty$ by using the self-consistent field method of statistical physics.
We derive a number of explicit results, and give numerical illustrations. We
then do numerical simulations of the model for finite $m$ and without random
mixing. We find in many of these cases that the self-consistent field method is
a very good approximation.


  As proteins typically have charges of around 10, they will interact strongly
with charged surfaces. We calculate the electrostatic contribution to the
interaction of crystals of protein with charged surfaces. The surfaces repel
like-charged crystals and attract oppositely-charged crystals, with free
energies which can be easily several kT per protein molecule brought into
contact with the surface. This means that oppositely charged surfaces can act
as a nucleant, they can induce nucleation of a protein crystal by lowering the
free energy barrier to heterogeneous nucleation of the crystal from a dilute
solution.


  For many years it had been believed that steric compatibility of helix
interfaces could be the source of the observed preference for particular angles
between neighbouring helices as emerging from statistical analysis of protein
databanks. Several elegant models describing how side chains on helices can
interdigitate without steric clashes were able to account quite reasonably for
the observed distributions. However, it was later recognized (Bowie, 1997 and
Walther, 1998) that the ``bare'' measured angle distribution should be
corrected to avoid statistical bias. Disappointingly, the rescaled
distributions dramatically lost their similarity with theoretical predictions
casting many doubts on the validity of the geometrical assumptions and models.
In this report we elucidate a few points concerning the proper choice of the
random reference distribution. In particular we show the existence of crucial
corrections due to the correct implementation of the approach used to
discriminate whether two helices are in contact or not and to measure their
relative orientations. By using this new rescaling, the ``true'' packing angle
preferences are well described, even more than with the original ``bare''
distribution, by regular packing models.


  We have measured the temperature driven denaturing, or melting transition in
poly d(A)-poly d(T) DNA oligomers of various lengths in different buffer
conditions. Our findings are in clear disagreement with two state, reaction
kinetics model, and we find that the so-called zipper model, where denaturing
proceeds through opening of the duplex at the ends describes well the
temperature dependence of the average number of open base pairs. Analysis of
the length dependence of the transition parameters however suggest that bubble
formation is important and that the transition, in the thermodynamic limit, is
continuous, albeit close to first order.


  Different aspects of protein folding are illustrated by simplified polymer
models. Stressing the diversity of side chains (residues) leads one to view
folding as the freezing transition of an heteropolymer. Technically, the most
common approach to diversity is randomness, which is usually implemented in two
body interactions (charges, polar character,..). On the other hand, the
(almost) universal character of the protein backbone suggests that folding may
also be viewed as the crystallization transition of an homopolymeric chain, the
main ingredients of which are the peptide bond and chirality (proline and
glycine notwithstanding). The model of a chiral dipolar chain leads to a
unified picture of secondary structures, and to a possible connection of
protein structures with ferroelectric domain theory.


  The electrostatic contribution to twist rigidity of DNA is studied, and it is
shown that the Coulomb self-energy of the double-helical sugar-phosphate
backbone contributes considerably to twist rigidity of DNA--the electrostatic
twist rigidity of DNA is found as $C_{\rm elec}\approx 5$ nm, which makes up
about 7% of its total twist rigidity ($C_{\rm DNA}\approx 75$ nm). The
electrostatic twist rigidity is found, however, to only weakly depend on the
salt concentration, because of a competition between two different screening
mechanisms: (1) Debye screening by the salt ions in the bulk, and (2)
structural screening by the periodic charge distribution along the backbone of
the helical polyelectrolyte. It is found that depending on the parameters, the
electrostatic contribution could stabilize or destabilize the structure of a
helical polyelectrolyte.


  There has been some confusion concerning the animal group-size: an
exponential distribution was deduced by maximizing the entropy; lognormal
distributions were practically used; a power-law decay with exponent {3/2} was
proposed in physical analogy to aerosol condensation. Here I show that the
animal group-size distribution follows a power-law decay with exponent 1, and
is truncated at a cut-off size which is the expected size of the groups an
arbitrary individual engages in. An elementary model of animal aggregation
based on binary splitting and coalescing on contingent encounter is presented.
The model predicted size distribution holds for various data from pelagic
fishes and mammalian herbivores in the wild.


  We presented a predator-prey lattice model containing moveable wolves and
sheep, which are characterized by Penna double bit strings. Sexual reproduction
and child-care strategies are considered. To implement this model in an
efficient way, we build a fully parallel Cellular Automata based on a new
definition of the neighborhood. We show the roles played by the initial
densities of the populations, the mutation rate and the linear size of the
lattice in the evolution of this model.


  We have simulated demographic changes in the human population using the Penna
microscopic model, based on the simple Monte Carlo method. The results of
simulations have shown that during a few generations changes in the genetic
pool of a population are negligible, while improving the methods of
compensation of genetic defects or genetically determined proneness to many
disorders drastically affects the average life span of organisms. Age
distribution and mortality of the simulated populations correspond very well to
real demographic data available from different countries. Basing on the
comparison of structures of real human populations and the results of
simulations it is possible to predict changes in the age structure of
populations in the future.


  Gene regulation in Eukaryotes is mainly effected through transcription
factors binding to rather short recognition motifs generally located upstream
of the coding region. We present a novel computational method to identify
regulatory elements in the upstream region of Eukaryotic genes. The genes are
grouped in sets sharing an overrepresented short motif in their upstream
sequence. For each set, the average expression level from a microarray
experiment is determined: if this level is significantly higher or lower than
the average taken over the whole genome, then the overrepresented motif shared
by the genes in the set is likely to play a role in their regulation. We
illustrate the method by applying it to the genome of {\it S. cerevisiae}, for
which many datasets of microarray experiments are publicly available. Several
known binding motifs are correctly recognized by our algorithm, and a new
candidate is suggested for experimental verification.


  Recently we have used a cellular automata model which describes the dynamics
of a multi-connected network to reproduce the refractory behavior and aging
effects obtained in immunization experiments performed with mice when subjected
to multiple perturbations. In this paper we investigate the similarities
between the aging dynamics observed in this multi-connected network and the one
observed in glassy systems, by using the usual tools applied to analyze the
latter. An interesting feature we show here is that the model reproduces the
biological aspects observed in the experiments during the long transient time
it takes to reach the stationary state. Depending on the initial conditions,
and without any perturbation, the system may reach one of a family of
long-period attractors. The pertrubations may drive the system from its natural
attractor to other attractors of the same family. We discuss the different
roles played by the small random perturbations (noise) and by the large
periodic perturbations (immunizations).


  We study the thermodynamical properties of a topology-based model proposed by
Galzitskaya and Finkelstein for the description of protein folding. We devise
and test three different mean-field approaches for the model, that simplify the
treatment without spoiling the description. The validity of the model and its
mean-field approximations is checked by applying them to the $\beta$-hairpin
fragment of the immunoglobulin-binding protein (GB1) and making a comparison
with available experimental data and simulation results. Our results indicate
that this model is a rather simple and reasonably good tool for interpreting
folding experimental data, provided the parameters of the model are carefully
chosen. The mean-field approaches substantially recover all the relevant exact
results and represent reliable alternatives to the Monte Carlo simulations.


  In several recent papers new gene-detection algorithms were proposed for
detecting protein-coding regions without requiring learning dataset of already
known genes. The fact that unsupervised gene-detection is possible closely
connected to existence of a cluster structure in oligomer frequency
distributions. In this paper we study cluster structure of several genomes in
the space of their triplet frequencies, using pure data exploration strategy.
Several complete genomic sequences were analyzed, using visualization of tables
of triplet frequencies in a sliding window. The distribution of 64-dimensional
vectors of triplet frequencies displays a well-detectable cluster structure.
The structure was found to consist of seven clusters, corresponding to
protein-coding information in three possible phases in one of the two
complementary strands and in the non-coding regions with high accuracy (higher
than 90% on the nucleotide level). Visualizing and understanding the structure
allows to analyze effectively performance of different gene-prediction tools.
Since the method does not require extraction of ORFs, it can be applied even
for unassembled genomes. The information content of the triplet distributions
and the validity of the mean-field models are analysed.


  We characterize the equilibrium thermodynamics of a thick polymer confined in
a spherical region of space. This is used to gain insight into the DNA
packaging process. The experimental reference system for the present study is
the recent characterization of the loading process of the genome inside the
$\phi$29 bacteriophage capsid. Our emphasis is on the modelling of
double-stranded DNA as a flexible thick polymer (tube) instead of a
beads-and-springs chain. By using finite-size scaling to extrapolate our
results to genome lengths appropriate for $\phi$29, we find that the
thickness-induced force may account for up to half the one measured
experimentally at high packing densities. An analogous agreement is found for
the total work that has to be spent in the packaging process. Remarkably, such
agreement can be obtained in the absence of any tunable parameters and is a
mere consequence of the DNA thickness. Furthermore, we provide a quantitative
estimate of how the persistence length of a polymer depends on its thickness.
The expression accounts for the significant difference in the persistence
lengths of single- and double-stranded DNA (again with the sole input of their
respective sections and natural nucleotide/base-pair spacing).


  In this paper we review the trajectory of a model proposed by Stauffer and
Weisbuch in 1992 to describe the evolution of the immune repertoire and present
new results about its dynamical behavior. Ten years later this model, which is
based on the ideas of the immune network as proposed by Jerne, has been able to
describe a multi-connected network and could be used to reproduce immunization
and aging experiments performed with mice. Besides its biological implications,
the physical aspects of the complex dynamics of this network is very
interesting {\it per se}. The immunization protocol is simulated by introducing
small and large perturbations (damages), and in this work we discuss the role
of both. In a very recent paper we studied the aging effects by using
auto-correlation functions, and the results obtained apparently indicated that
the small perturbations would be more important than the large ones, since
their cumulative effects may change the attractor of the dynamics. However our
new results indicate that both types of perturbations are important. It is the
cooperative effects between both that lead to the complex behavior which allows
to reproduce experimental results.


  We investigate theoretically the translocation of structured RNA/DNA
molecules through narrow pores which allow single but not double strands to
pass. The unzipping of basepaired regions within the molecules presents
significant kinetic barriers for the translocation process. We show that this
circumstance may be exploited to determine the full basepairing pattern of
polynucleotides, including RNA pseudoknots. The crucial requirement is that the
translocation dynamics (i.e., the length of the translocated molecular segment)
needs to be recorded as a function of time with a spatial resolution of a few
nucleotides. This could be achieved, for instance, by applying a mechanical
driving force for translocation and recording force-extension curves (FEC's)
with a device such as an atomic force microscope or optical tweezers. Our
analysis suggests that with this added spatial resolution, nanopores could be
transformed into a powerful experimental tool to study the folding of nucleic
acids.


  The role of the rigidity of a peptide chain in its equilibrium dynamics is
investigated within a realistic model with stringent microscopically derived
coupling interaction potential and effective on-site potential. The coupling
interaction characterizing the chain rigidity and the effective on-site
potentials are calculated for three main types of protein secondary structure.
The coupling interaction is found to be surprisingly weak for all of them but
different in character: repulsive for alpha-helix and anti-parallel beta-sheet
structures and attractive for parallel beta-sheet structure. The effective
on-site potential is found to be a hard one for alpha-helix and anti-parallel
beta-sheet and a soft one for parallel beta-sheet. In all three types of
protein secondary structures a stable zig-zag shape discrete breather (DB)
associated with the oscillations of torsional (dihedral) angles can exist due
to weakness of the coupling interaction. However, since the absorption of far
infrared radiation (IR) by proteins is known to require the existence of rather
long chains of hydrogen bonds that takes place only in alpha-helicies, then one
can conclude that the excitation of a DB in such a way is possible in
alpha-helix and seems to be hardly possible in beta-sheet structures. The
interpretation of the recent experiments of Xie et al. on far IR laser pulse
spectroscopy of proteins is suggested. The frequency of a DB in an the
alpha-helix is obtained in the region of 115 cm^{-1} in accordance with these
experiments.


  It is well known amongst molecular biologists that proteins with a common
ancestor and that perform the same function in similar organisms, can have
rather different amino-acid sequences. Mutations have altered the amino-acid
sequences without affecting the function; this is called neutral evolution. A
simple model of a protein in which the interactions are encoded by sequences of
bits is introduced, and used to study how mutations can change these bits, and
hence the interactions, while maintaining the stability of the protein
solution. This stability is a simple minimal requirement on our model proteins
which mimics part of the requirement on a real protein to be functional. The
properties of our model protein, such as its second virial coefficient, are
found to vary significantly from one model protein to another. It is suggested
that this may also be the case for real proteins in vivo.


  We show that Kramers rate theory gives a straightforward, accurate estimate
of the closing time $\tau_c$ of a semiflexible polymer that is valid in cases
of physical interest. The calculation also reveals how the time scales of chain
relaxation and closing are intertwined, illuminating an apparent conflict
between two ways of calculating $\tau_c$ in the flexible limit.


  Using a general Green function formulation, we re-derive, both, (i) Spitzer
and his followers results for the winding angle distribution of the planar
Brownian motion, and (ii) Edwards-Prager-Frisch results on the statistical
mechanics of a ring polymer entangled with a straight bar. In the statistical
mechanics part, we consider both cases of quenched and annealed topology. Among
new results, we compute exactly the (expectation value of) the surface area of
the locus of points such that each of them has linking number $n$ with a given
closed random walk trajectory ($=$ ring polymer). We also consider the
generalizations of the problem for the finite diameter (disc-like) obstacle and
winding within a cavity.


  A global optimization method, conformational space annealing (CSA), is
applied to study a 46-residue protein with the sequence
B_9N_3(LB)_4N_3B_9N_3(LB)_5L, where B, L and N designate hydrophobic,
hydrophilic, and neutral residues, respectively. The 46-residue BLN protein is
folded into the native state of a four-stranded beta-barrel. It has been a
challenging problem to locate the global minimum of the 46-residue BLN protein
since the system is highly frustrated and consequently its energy landscape is
quite rugged. The CSA successfully located the global minimum of the 46-mer for
all 100 independent runs. The CPU time for CSA is about seventy times less than
that for simulated annealing (SA), and its success rate (100 %) to find the
global minimum is about eleven times higher. The amount of computational
efforts used for CSA is also about ten times less than that of the best global
optimization method yet applied to the 46-residue BLN protein, the quantum
thermal annealing with renormalization. The 100 separate CSA runs produce the
global minimum 100 times as well as other 5950 final conformations
corresponding to a total of 2361 distinct local minima of the protein. Most of
the final conformations have relatively small RMSD values from the global
minimum, independent of their diverse energy values. Very close to the global
minimum, there exist quasi-global-minima which are frequently obtained as one
of the final answers from SA runs. We find that there exist two largest energy
gaps between the quasi-global-minima and the other local minima. Once a SA run
is trapped in one of these quasi-global-minima, it cannot be folded into the
global minimum before crossing over the two large energy barriers, clearly
demonstrating the reason for the poor success rate of SA.


  The determination of protein functions is one of the most challenging
problems of the post-genomic era. The sequencing of entire genomes and the
possibility to access gene's co-expression patterns has moved the attention
from the study of single proteins or small complexes to that of the entire
proteome. In this context, the search for reliable methods for proteins'
function assignment is of uttermost importance. Previous approaches to deduce
the unknown function of a class of proteins have exploited sequence
similarities or clustering of co-regulated genes, phylogenetic profiles,
protein-protein interactions, and protein complexes. We propose to assign
functional classes to proteins from their network of physical interactions, by
minimizing the number of interacting proteins with different categories. The
function assignment is made on a global scale and depends on the entire
connectivity pattern of the protein network. Multiple functional assignments
are made possible as a consequence of the existence of multiple equivalent
solutions. The method is applied to the yeast Saccharomices Cerevisiae
protein-protein interaction network. Robustness is tested in presence of a high
percentage of unclassified proteins and under deletion/insertion of
interactions.


  We develop simple models for the global spread of infectious diseases,
emphasizing human mobility via air travel and the variation of public health
infrastructure from region to region. We derive formulas relating the total and
peak number of infections in two countries to the rate of travel between them
and their respective epidemiological parameters.


  We show that the icosahedral packings of protein capsomeres proposed by
Caspar and Klug for spherical viruses become unstable to faceting for
sufficiently large virus size, in analogy with the buckling instability of
disclinations in two-dimensional crystals. Our model, based on the nonlinear
physics of thin elastic shells, produces excellent one parameter fits in real
space to the full three-dimensional shape of large spherical viruses. The
faceted shape depends only on the dimensionless Foppl-von Karman number
\gamma=YR^2/\kappa, where Y is the two-dimensional Young's modulus of the
protein shell, \kappa is its bending rigidity and R is the mean virus radius.
The shape can be parameterized more quantitatively in terms of a spherical
harmonic expansion. We also investigate elastic shell theory for extremely
large \gamma, 10^3 < \gamma < 10^8, and find results applicable to icosahedral
shapes of large vesicles studied with freeze fracture and electron microscopy.


  Very recently we have developed a dynamic network model for eco-systems that
achieved ``unification'' of ``micro'' and ``macro''-evolution. We now propose
an extension of our model so as to stabilize the eco-system and describe {\it
speciation} in a more realistic manner.


  We study voltage driven translocation of a single stranded (ss) DNA through a
membrane channel. Our model, based on a master equation (ME) approach,
investigates the probability density function (pdf) of the translocation times,
and shows that it can be either double or mono-peaked, depending on the system
parameters. We show that the most probable translocation time is proportional
to the polymer length, and inversely proportional to the first or second power
of the voltage, depending on the initial conditions. The model recovers
experimental observations on hetro-polymers when using their properties inside
the pore, such as stiffness and polymer-pore interaction.


  We study the evolution of the network properties of a populated network
embedded in a genotype space characterised by either a low or a high number of
potential links, with particular emphasis on the connectivity and clustering.
Evolution produces two distinct types of network. When a specific genotype is
only able to influence a few other genotypes, the ecology consists of separate
non-interacting clusters in genotype space. When different types may influence
a large number of other sites, the network becomes one large interconnected
cluster. The distribution of interaction strengths -- but not the number of
connections -- changes significantly with time. We find that the species
abundance is only realistic for a high level of species connectivity. This
suggests that real ecosystems form one interconnected whole in which selection
leads to stronger interactions between the different types. Analogies with
niche and neutral theory are also considered.


  Modes of speciation have been the subject of a century's debate.
Traditionally, most speciations are believed to be caused by spatial separation
of populations (allopatry). Recent observations (Meyer 1990, Schliewen 1994,
Schliewen 2001, Rico 2002) and models (MaynardSmith 1966, Antonovics 1971,
Dickinson 1973, Rosenzweig 1978, T urner 1995, Noest 1997, Geritz 1998,
Kondrashov 1999, Dieckmann 1999, Doebeli 2000, Slatkin 1980), show that
speciation can also take place in sympatry. We discuss a comprehensive model of
coupled differentiation in phenotype, mating, and space, showing that spatial
segregation can be an induced process following a sympatric differentiation.
This is found to be a generic mechanism of adaptation to heterogeneous
environments, for which we propose the term diapatric speciation (Greek). It
explains the ubiquitous spatial patching of newly formed species, despite their
sympatric origin (Schliewen 1994, Schliewen 2001, Rico 2002).


  Molecular dynamics simulation methods are used to study the folding of
polymer chains into packed cubic states. The polymer model, based on a chain of
linked sites moving in the continuum, includes both excluded volume and
torsional interactions. Different native-state packing arrangements and chain
lengths are explored; the organization of the native state is found to affect
both the ability of the chain to fold successfully and the nature of the
folding pathway as the system is gradually cooled. An order parameter based on
contact counts is used to provide information about the folding process, with
contacts additionally classified according to criteria such as core and surface
sites or local and distant site pairs. Fully detailed contact maps and their
evolution are also examined.


  Many complex networks, including human societies, the Internet, the World
Wide Web and power grids, have surprising properties that allow vertices
(individuals, nodes, Web pages, etc.) to be in close contact and information to
be transferred quickly between them. Nothing is known of the emerging
properties of animal societies, but it would be expected that similar trends
would emerge from the topology of animal social networks. Despite its small
size (64 individuals), the Doubtful Sound community of bottlenose dolphins has
the same characteristics. The connectivity of individuals follows a complex
distribution that has a scale-free power-law distribution for large k. In
addition, the ability for two individuals to be in contact is unaffected by the
random removal of individuals. The removal of individuals with many links to
others does affect the length of the information path between two individuals,
but, unlike other scale-free networks, it does not fragment the cohesion of the
social network. These self-organizing phenomena allow the network to remain
united, even in the case of catastrophic death events.


  Jumps between neighboring minima in the energy landscape of both
homopolymeric and heteropolymeric chains are numerically investigated by
determining the average escape time from different valleys. The numerical
results are compared to the theoretical expression derived by Langer [J.S.
Langer, Ann. Phys. 54 (1969) 258] with reference to a 2N-dimensional space. Our
simulations indicate that the dynamics within the native valley is well
described by a sequence of thermally activated process up to temperatures well
above the folding temperature. At larger temperatures, systematic deviations
from the Langer's estimate are instead observed. Several sources for such
discrepancies are thoroughly discussed.


  Duplication graphs are graphs that grow by duplication of existing vertices,
and are important models of biological networks, including protein-protein
interaction networks and gene regulatory networks. Three models of graph growth
are studied: pure duplication growth, and two two-parameter models in which
duplication forms one element of the growth dynamics. A power-law degree
distribution is found to emerge in all three models. However, the parameter
space of the latter two models is characterized by a range of parameter values
for which duplication is the predominant mechanism of graph growth. For
parameter values that lie in this ``duplication-dominated'' regime, it is shown
that the degree distribution either approaches zero asymptotically, or
approaches a non-zero power-law degree distribution very slowly. In either
case, the approach to the true asymptotic degree distribution is characterized
by a dependence of the scaling exponent on properties of the initial degree
distribution. It is therefore conjectured that duplication-dominated,
scale-free networks may contain identifiable remnants of their early structure.
This feature is inherited from the idealized model of pure duplication growth,
for which the exact finite-size degree distribution is found and its asymptotic
properties studied.


  We propose and study a hydrodynamic model for pattern formation in mixtures
of molecular motors and microtubules. The steady state patterns we obtain in
different regimes of parameter space include arrangements of vortices and
asters separately as well as aster-vortex mixtures and fully disordered states.
Such stable steady states are observed in experiments in vitro. The sequence of
patterns obtained in the experiments can be associated with smooth trajectories
in a non-equilibrium phase diagram for our model.


  In nonlinear voter models the transitions between two states depend in a
nonlinear manner on the frequencies of these states in the neighborhood. We
investigate the role of these nonlinearities on the global outcome of the
dynamics for a homogeneous network where each node is connected to $m=4$
neighbors. The paper unfolds in two directions. We first develop a general
stochastic framework for frequency dependent processes from which we derive the
macroscopic dynamics for key variables, such as global frequencies and
correlations. Explicite expressions for both the mean-field limit and the pair
approximation are obtained. We then apply these equations to determine a phase
diagram in the parameter space that distinguishes between different dynamic
regimes. The pair approximation allows us to identify three regimes for
nonlinear voter models: (i) complete invasion, (ii) random coexistence, and --
most interestingly -- (iii) correlated coexistence. These findings are
contrasted with predictions from the mean-field phase diagram and are confirmed
by extensive computer simulations of the microscopic dynamics.


  It is not merely the position of residues that are of utmost importance in
protein function and stability, but the interactions between them. We
illustrate, by using a network construction on a set of 595 non-homologous
proteins, that regular packing is preserved in short-range interactions, but
short average path lengths are achieved through some long-range contacts. Thus,
lying between the two extremes of regularity and randomness, residues in folded
proteins are distributed according to a "small-world" topology. Using this
topology, we show that the core residues have the same local packing
arrangements irrespective of protein size. Furthermore, we find that the
average shortest path lengths are highly correlated with residue fluctuations,
providing a link between the spatial arrangement of the residues and protein
dynamics.


  Non-reversing relaxation enthalpies (DHnr) at glass transitions Tg(x) in the
PxGexSe1-2x ternary display a wide, sharp and deep global minimum (~0) in the
0.09 < x < 0.145 range, within which Tg becomes thermally reversing. In the
reversibility window these glasses are found not to age, in contrast to aging
observed for fragile glass compositions outside the window. Thermal
reversibility and lack of aging are paradigms that molecular glasses in the
window share with proteins in transition states, which result from structural
self-organization in both systems. In proteins the self-organized structures
appear to be at places where life sustaining repeating foldings and unfoldings
occur.


  In complex systems such as spin systems and protein systems, conventional
simulations in the canonical ensemble will get trapped in states of energy
local minima. We employ the generalized-ensemble algorithms in order to
overcome this multiple-minima problem. Two well-known generalized-ensemble
algorithms, namely, multicanonical algorithm and replica-exchange method, are
described. We then present four new generalized-ensemble algorithms as further
extensions of the two methods. Effectiveness of the new methods are illustrated
with a Potts model, Lennard-Jones fluid system, and protein system.


  A new family of compound Poisson distribution functions from statistical
linguistic is used to study the n-tuples and nucleotide composition features of
DNA sequences. The relative frequency distribution of the 6-tuples and 7-
tuples occurrence studies suggest that most of the DNA sequences follow the
general shape of the compound Poisson distribution. It is also noted that the
$\chi$-square test indicated that some of the sequences follow this
distribution with a reasonable level of goodness of fit. The compositional
segmentation study fits quite well using this new family of distribution
functions. Furthermore, the absolute values of the relative frequency come out
naturally from the linguistic model without ambiguity. It is suggesting that
DNA sequences are not random sequences and they could possibly have subsequence
structures.


  We investigate the translocation of a single stranded DNA through a pore
which fluctuates between two conformations, using coupled master equations. The
probability density function of the first passage times (FPT) of the
translocation process is calculated, displaying a triple, double or mono peaked
behavior, depending on the interconversion rates between the conformations, the
applied electric field, and the initial conditions. The cumulative probability
function of the FPT, in a field-free environment, is shown to have two regimes,
characterized by fast and slow timescales. An analytical expression for the
mean first passage time of the translocation process is derived, and provides,
in addition to the interconversion rates, an extensive characterization of the
translocation process. Relationships to experimental observations are
discussed.


  Increasing every year the retirement age by the same amount as the increase
of the life expectancy gives roughly stable ratios of the number of retired to
working-age people in industrialized countries. Continuous influx of
immigrants, below one percent per year of the total population, is needed for
this stabilization.


  Interaction networks are of central importance in post-genomic molecular
biology, with increasing amounts of data becoming available by high-throughput
methods. Examples are gene regulatory networks or protein interaction maps. The
main challenge in the analysis of these data is to read off biological
functions from the topology of the network. Topological motifs, i.e., patterns
occurring repeatedly at different positions in the network have recently been
identified as basic modules of molecular information processing. In this paper,
we discuss motifs derived from families of mutually similar but not necessarily
identical patterns. We establish a statistical model for the occurrence of such
motifs, from which we derive a scoring function for their statistical
significance. Based on this scoring function, we develop a search algorithm for
topological motifs called graph alignment, a procedure with some analogies to
sequence alignment. The algorithm is applied to the gene regulation network of
E. coli.


  Using each node's degree as a proxy for its importance, the topological
hierarchy of a complex network is introduced and quantified. We propose a
simple dynamical process used to construct networks which are either maximally
or minimally hierarchical. Comparison with these extremal cases as well as with
random scale-free networks allows us to better understand hierarchical versus
modular features in several real-life complex networks. For random scale-free
topologies the extent of topological hierarchy is shown to smoothly decline
with $\gamma$ -- the exponent of a degree distribution -- reaching its highest
possible value for $\gamma \leq 2$ and quickly approaching zero for $\gamma>3$.


  The binding of basic fibroblast growth factor (FGF-2) to its cell surface
receptor (CSR) and subsequent signal transduction is known to be enhanced by
Heparan Sulfate Proteoglycans (HSPGs). HSPGs bind FGF-2 with low affinity and
likely impact CSR-mediated signaling via stabilization of FGF-2-CSR complexes
via association with both the ligand and the receptor. What is unknown is
whether HSPG associates with CSR in the absence of FGF-2. In this paper, we
determine conditions by which pre-association would impact CSR-FGF-2-HSPG triad
formation assuming diffusion-limited surface reactions. Using mean-field rate
equations, we show that (i) when [HSPG] is much higher than [CSR], the presence
of pre-formed complexes does not affect the steady state of FGF-2 binding, and
(ii) when the concentrations are comparable, the presence of pre-formed
complexes substantially increases the steady state concentration of FGF-2 bound
to CSR. These findings are supported by explicit cellular automaton
simulations, which justify the mean-field treatment. We discuss the advantages
of such a two-receptor system compared to a single receptor model, when the
parameters are comparable. Further, we speculate that the observed high
concentration of HSPG in intact cells ([HSPG] ~ 100[CSR]) provides a way to
ensure that the binding levels of FGF-2 to its signaling receptor remains high,
irrespective of the presence of pre-formed CSR-HSPG complexes on the cell
surface, while allowing the cell to finely tune the response to FGF-2 via
down-regulation of the signaling receptor.


  Segments with the amino acid sequence EKAYLRT appear in natural occurring
proteins both in $\alpha$-helices and $\beta$-sheets. For this reason, we have
use this peptide to study how secondary structure formation in proteins depends
on the local environment. Our data rely on multicanonical Monte Carlo
simulations where the interactions among all atoms are taken into account.
Results in gas phase are compared with that in an implicit solvent. We find
that both in gas phase and solvated EKAYLRT forms an $\alpha$-helix when not
interacting with other molecules. However, in the vicinity of a $\beta$-strand,
the peptide forms a $\beta$-strand. Because of this change in secondary
structure our peptide may provide a simple model for the $\alpha \to \beta$
transition that is supposedly related to the outbreak of Prion diseases and
similar illnesses.


  Epidemiological processes are studied within a recently proposed hierarchical
network model using the susceptible-infected-refractory dynamics of an
epidemic. Within the network model, a population may be characterized by $H$
independent hierarchies or dimensions, each of which consists of groupings of
individuals into layers of subgroups. Detailed numerical simulations reveal
that for $H>1$, global spreading results regardless of the degree of homophily
of the individuals forming a social circle. For H=1, a transition from global
to local spread occurs as the population becomes decomposed into increasingly
homophilous groups. Multiple dimensions in classifying individuals (nodes) thus
make a society (computer network) highly susceptible to large scale outbreaks
of infectious diseases (viruses).


  We study the interplay of activity, order and flow through a set of
coarse-grained equations governing the hydrodynamic velocity, concentration and
stress fields in a suspension of active, energy-dissipating particles. We make
several predictions for the rheology of such systems, which can be tested on
bacterial suspensions, cell extracts with motors and filaments, or artificial
machines in a fluid. The phenomena of cytoplasmic streaming, elastotaxis and
active mechanosensing find natural explanations within our model.


  We twice modify the Penna model for biological ageing. First we introduce
back (good) mutations and a memory for them into the model. It allows us to
observe an improvement of the species fitness over long time scales as well as
punctuated equilibrium. Second we adopt a food/space competition factor that
depends on the number of accumulated mutations in the individuals genomes, and
get rid of the fixed limiting number of allowed mutations. Besides reproducing
the main results of the standard model, we also observe a mortality maximum for
the oldest old.


  We removed from the Penna model for biological ageing any random killing
Verhulst factor. Deaths are due only to genetic diseases and the population
size is fixed, instead of fluctuating around some constant value. We show that
these modifications give qualitatively the same results obtained in an earlier
paper, where the random killings (used to avoid an exponential increase of the
population) were applied only to newborns.


  We introduce a model of protein evolution to explain limitations in the
immune system response to vaccination and disease. The phenomenon of original
antigenic sin, wherein vaccination creates memory sequences that can
\emph{increase} susceptibility to future exposures to the same disease, is
explained as stemming from localization of the immune system response in
antibody sequence space. This localization is a result of the roughness in
sequence space of the evolved antibody affinity constant for antigen and is
observed for diseases with high year-to-year mutation rates, such as influenza.


  We introduce a population dynamics model, where individual genomes are
represented by bit-strings. Selection is described by death probabilities which
depend on these genomes, and new individuals continuously replace the ones that
die, keeping the population constant. An offspring has the same genome as its
(randomly chosen) parent, except for a small amount of (also random) mutations.
Chance may thus generate a newborn with a genome that is better than that of
its parent, and the newborn will have a smaller death probability. When this
happens, this individual is a would-be founder of a new lineage. A new lineage
is considered created if its alive descendence grows above a certain previously
defined threshold. The time evolution of populations evolving under these rules
is followed by computer simulations and the probability densities of lineage
duration and size, among others, are computed. These densities show a
scale-free behaviour, in accordance with some conjectures in paleoevolution,
and suggesting a simple mechanism as explanation for the ubiquity of these
power-laws.


  We propose a novel method to optimize existing force-field parameters for
protein systems. The method consists of minimizing the summation of the square
of the force acting on each atom in the proteins with the structures from the
Protein Data Bank. We performed this optimization to the partial-charge and
torsion-energy parameters of the AMBER parm96 force field, using 100 molecules
from the Protein Data Bank. We then performed folding simulations of
$\alpha$-helical and $\beta$-hairpin peptides. The optimized force-field
parameters gave structures more similar to the experimental implications than
the original AMBER force field.


  Protein folds are built primarily from the packing together of two types of
structures: alpha-helices and beta-sheets. Neither structure is rigid, and the
flexibility of helices and sheets is often important in determining the final
fold ({\it e.g.}, coiled coils and beta-barrels). Recent work has quantified
the flexibility of alpha-helices using a principal-component analysis (PCA) of
database helical structures (Emberly, 2003). Here, we extend the analysis to
beta-sheet flexibility using PCA on a database of beta-sheet structures. For
sheets of varying dimension and geometry, we find two dominant modes of
flexibility: twist and bend. The distributions of amplitudes for these modes
are found to be Gaussian and independent, suggesting that the PCA twist and
bend modes can be identified as the soft elastic normal modes of sheets. We
consider the scaling of mode eigenvalues with sheet size and find that parallel
beta-sheets are more rigid than anti-parallel sheets over the entire range
studied. Lastly, we discuss the application of our PCA results to modeling and
design of beta-sheet proteins.


  Generalized-ensemble algorithm and diffusion theory have been combined in
order to compute the dynamical properties monitored by nuclear magnetic
resonance experiments from efficient and reliable evaluation of statistical
averages. Replica-exchange Monte Carlo simulations have been performed with a
C-peptide analogue of ribonuclease A, and Smoluchowski diffusion equations have
been applied. A fairly good agreement between the calculated and measured
$^1$H-NOESY NMR cross peaks has been obtained. The combination of these
advanced and continuously improving statistical tools allows the calculation of
a wide variety of dynamical properties routinely obtained by experiments.


  We report results showing that thermally-induced openings of double stranded
DNA coincide with the location of functionally relevant sites for
transcription. Investigating both viral and bacterial DNA gene promoter
segments, we found that the most probable opening occurs at the transcription
start site. Minor openings appear to be related to other regulatory sites. Our
results suggest that coherent thermal fluctuations play an important role in
the initiation of transcription. Essential elements of the dynamics, in
addition to sequence specificity, are nonlinearity and entropy, provided by
local base-pair constraints.


  Using a lattice model based on Monte Carlo simulations, we study the role of
the reproduction pattern on the fate of an evolving population. Each individual
is under the selection pressure from the environment and random mutations. The
habitat ("climate") is changing periodically. Evolutions of populations
following two reproduction patterns are compared, asexual and sexual. We show,
via Monte Carlo simulations, that sexual reproduction by keeping more
diversified populations gives them better chances to adapt themselves to the
changing environment. However, in order to obtain a greater chance to mate, the
birth rate should be high. In the case of low birth rate and high mutation
probability there is a preference for the asexual reproduction.


  A statistical mechanical distance constraint model (DCM) is presented that
explicitly accounts for network rigidity among constraints present within a
system. Constraints are characterized by local microscopic free energy
functions. Topological re-arrangements of thermally fluctuating constraints are
permitted. The partition function is obtained by combining microscopic free
energies of individual constraints using network rigidity as an underlying
long-range mechanical interaction -- giving a quantitative explanation for the
non-additivity in component entropies exhibited in molecular systems. Two
exactly solved 2-dimensional toy models representing flexible molecules that
can undergo conformational change are presented to elucidate concepts, and to
outline a DCM calculation scheme applicable to many types of physical systems.
It is proposed that network rigidity plays a central role in balancing the
energetic and entropic contributions to the free energy of bio-polymers, such
as proteins. As a demonstration, the distance constraint model is solved
exactly for the alpha-helix to coil transition in homogeneous peptides.
Temperature and size independent model parameters are fitted to Monte Carlo
simulation data, which includes peptides of length 10 for gas phase, and
lengths 10, 15, 20 and 30 in water. The DCM is compared to the Lifson-Roig
model. It is found that network rigidity provides a mechanism for cooperativity
in molecular structures including their ability to spontaneously self-organize.
In particular, the formation of a characteristic topological arrangement of
constraints is associated with the most probable microstates changing under
different thermodynamic conditions.


  We compare the speed with which a sexual, respectively an asexual, population
is able to respond to a biased selective pressure. Our model focuses on the
Weismann hypothesis that the extra variation caused by crossing-over and
recombination during sexual reproduction allows a sexual population to adapt
faster. We find, however, that the extra variation amongst the progeny produced
during sexual reproduction for most model parameters is unable to overcome the
effect that parents with a high individual fitness in general must mate with
individuals of lower individual fitness resulting in a moderate reproductive
fitness for the pair.


  Biological and biomimetic membranes often contain aggregates of embedded or
adsorbed macromolecules. In this article, the indirect interactions of
cylindrical objects adhering to a planar membrane are considered theoretically.
The adhesion of the cylinders causes a local perturbation of the equilibrium
membrane shape, which leads to membrane-mediated interactions. For a planar
membrane under lateral tension, the interaction is repulsive for a pair of
cylinders adhering to the same side of the membrane, and attractive for
cylinders adhering at opposite membrane sides. For a membrane in an external
harmonic potential, the interaction of adsorbed cylinders is always attractive
and increases if forces perpendicular to the membrane act on the cylinders.


  We propose a method for predicting helical membrane protein structures by
computer simulations. Our method consists of two parts. In the first part,
amino-acid sequences of the transmembrane helix regions are obtained from one
of existing WWW servers. In the second part, we perform a replica-exchange
simulation of these transmembrane helices with some constraints and identify
the predicted structure as the global-minimum-energy state. We have tested the
method with the dimeric transmembrane domain of glycophorin A. The structure
obtained from the prediction was in close agreement with the experimental data.


  We present a short review of various experiments that measure charge transfer
and charge transport in DNA. Some general comments are made on the possible
connection between 'chemistry-style' charge transfer experiments that probe
fluorescence quenching and remote oxidative damage and 'physics-style'
measurements that measure transport properties as defined typically in the
solid-state. We then describe measurements performed by our group on the
millimeter wave response of DNA. By measuring over a wide range of humidity
conditions and comparing the response of single strand DNA and double strand
DNA, we show that the appreciable AC conductivity of DNA is not due to photon
assisted hopping between localized states, but instead due to dissipation from
dipole motion in the surrounding water helix.


  We studied the crystallization of a globular protein, lysozyme, in the cubic
phase of the lipid monoolein. The solubility of lysozyme in salt solution
decreased by a factor of $\sim 4$ when confined in cubic phase. Calculations
and Monte Carlo simulations show that this can be explained by the {\it
confinement} of lysozyme molecules to the narrow water cells in the cubic
phase.


  We explore and quantify the physical and biochemical mechanisms that may be
relevant in the regulation of translation. After elongation and detachment from
the 3' termination site of mRNA, parts of the ribosome machinery can diffuse
back to the initiation site, especially if it is held nearby, enhancing overall
translation rates. The elongation steps of the mRNA-bound ribosomes are modeled
using exact and asymptotic results of the totally asymmetric exclusion process
(TASEP).Since the ribosome injection rates of the TASEP depend on the local
concentrations at the initiation site, a source of ribosomes emanating from the
termination end can feed back to the initiation site, leading to a
self-consistent set of equations for the steady-state ribosome throughput.
Additional mRNA binding factors can also promote loop formation, or
cyclization, bringing the initiation and termination sites into close
proximity. The probability distribution of the distance between the initiation
and termination sites is described using simple noninteracting polymer models.
We find that the initiation, or initial ribosome adsorption binding required
for maximal throughput can vary dramatically depending on certain values of the
bulk ribosome concentration and diffusion constant. If cooperative interactions
among the loop-promoting proteins and the initiation/termination sites are
considered, the throughput can be further regulated in a nonmonotonic manner.
Potential experiments to test the hypothesized physical mechanisms are
discussed.


  We provide evidence that the energy landscapes of folded proteins do not
shift with temperature, but the onset of functional dynamics is associated with
its effective sampling. The motion of the backbone is described by three
distinct regimes. One is associated with slow time scales of the activity along
the envelope of the energy surface defining the folded protein. Another, with
fast time scales, is due to activity along the pockets decorating the
folded-state envelope. The intermediate regime emerges at temperatures where
jumps between the pockets become possible, leading to an active protein.


  We study the static and dynamical properties of DNA in the vicinity of its
melting transition, i.e. the separation of the two strands upon heating. The
investigation is based on a simple mechanical model which includes the
helicoidal geometry of the molecule and allows an exact numerical evaluation of
its thermodynamical properties. Dynamical simulations of long-enough molecular
segments allow the study of the structure factors and of the properties of the
denaturated regions. Simulations of finite chains display the hallmarks of a
first order transition for sufficiently long-ranged stacking forces although a
study of the model's ``universality class'' strongly suggests the presence of
an ``underlying'' continuous transition.


  We address the relation between long range correlations and charge transfer
efficiency in aperiodic artificial or genomic DNA sequences. Coherent charge
transfer through the HOMO states of the guanine nucleotide is studied using the
transmission approach, and focus is made on how the sequence-dependent
backscattering profile can be inferred from correlations between base pairs.


  We apply the general protocol of parameter optimization (Lee, J. et al. Phys.
Chem. B 2001, 105, 7291) to the UNRES potential. In contrast to the earlier
works where only the relative weights of various interaction terms were
optimized, we optimize all linear parameters of the potential. The method
exploits the high efficiency of the conformal space annealing method in finding
distinct low energy conformations. For a given training set of proteins, the
parameters are modified to make the native-like conformations energetically
more favorable than the non-native ones. Linear approximation is used to
estimate the energy change due to the parameter modification. The parameter
change is followed by local energy reminimization and new conformational
searches to find the energies of native-like and non-native local minima of the
energy function with new parameters. These steps are repeated until the
potential predicts a native-like conformation as one of the low energy
conformations for each protein in the training set. We consider a training set
of crambin (PDB ID 1ejg), 1fsd, and the 10-55 residue fragment of
staphylococcal protein A (PDB ID 1bdd). As the first check for the feasibility
of our protocol, we optimize the parameters separately for these proteins and
find an optimal set of parameters for each of them. Next we apply the method
simultaneously to these three proteins. By refining all linear parameters, we
obtain an optimal set of parameters from which the native-like conformations of
the all three proteins are retrieved as the global minima, without introducing
additional multi-body energy terms.


  We propose a concept of autonomous energy transducer at a molecular scale,
where output is produced with small input energy, of the same order of the
thermal energy, without restriction of magnitude or timing of input, and
without any control after the input. As an example that satisfies these
requisites, a dynamical systems model with several degrees of freedom is
proposed, which transduces input energy to output motion on the average. It is
shown that this transduction is robust and the coupling between the input and
output is generally loose. How this transducer works is analyzed in terms of
dynamical systems theory, where chaotic dynamics of the internal degrees of
freedom, as well as duration of active state which is self-organized with the
energy flow, is essential. We also discuss possible relationships to enzyme
dynamics or protein motors.


  A new general algorithm for optimization of potential functions for protein
folding is introduced. It is based upon gradient optimization of the
thermodynamic stability of native folds of a training set of proteins with
known structure. The iterative update rule contains two thermodynamic averages
which are estimated by (generalized ensemble) Monte Carlo. We test the learning
algorithm on a Lennard-Jones (LJ) force field with a torsional angle
degrees-of-freedom and a single-atom side-chain. In a test with 24 peptides of
known structure, none folded correctly with the initial potential functions,
but two-thirds came within 3{\AA} to their native fold after optimizing the
potential functions.


  The denaturation of double-stranded DNA as function of force and temperature
is discussed. At room temperature, sequence heterogeneity dominates the physics
of single molecule force-extension curves starting about 7 piconewtons of below
a ~15 pN unzipping transition. The dynamics of the unzipping fork exhibits
anomalous drift and diffusion in a similar range above this transition. Energy
barriers near the transition scale as the square root of the genome size.
Recent observations of jumps and plateaus in the unzipping of lambda phage DNA
at constant force are consistent with these predictions.


  We study a generic reaction-diffusion model for single-species population
dynamics that includes reproduction, death, and competition. The population is
assumed to be confined in a refuge beyond which conditions are so harsh that
they lead to certain extinction. Standard continuum mean field models in one
dimension yield a critical refuge length L_c such that a population in a refuge
larger than this is assured survival. Herein we extend the model to take into
account the discreteness and finiteness of the population, which leads us to a
stochastic description. We present a particular critical criterion for likely
extinction, namely, that the standard deviation of the population be equal to
the mean. According to this criterion, we find that while survival can no
longer be guaranteed for any refuge size, for sufficiently weak competition one
can make the refuge large enough (certainly larger than L_c) to cause
extinction to be unlikely. However, beyond a certain value of the competition
rate parameter it is no longer possible to escape a likelihood of extinction
even in an infinite refuge. These unavoidable fluctuations therefore have a
severe impact on refuge design issues.


  We consider the general properties of the quasispecies dynamical system from
the standpoint of its evolution and stability. Vector field analysis as well as
spectral properties of such system has been studied. Mathematical modelling of
the system under consideration has been performed.


  Three results are presented: First, we solve the problem of persistence of
dissipation for reduction of kinetic models. Kinetic equations with
thermodynamic Lyapunov functions are studied. Uniqueness of the thermodynamic
projector is proven: There exists only one projector which transforms the
arbitrary vector field equipped with the given Lyapunov function into a vector
field with the same Lyapunov function for a given anzatz manifold which is not
tangent to the Lyapunov function levels. Moreover, from the requirement of
persistence of the {\it sign} of dissipation follows that the {\it value} of
dissipation (the entropy production) persists too. The explicit construction of
this {\it thermodynamic projector} is described. Second, we use this projector
for developing the short memory approximation and coarse-graining for general
nonlinear dynamic systems. We prove that in this approximation the entropy
production grows. ({\it The theorem about entropy overproduction.}) In example
we apply the thermodynamic projector to derivation the equations of reduced
kinetics for the Fokker-Planck equation. The new class of closures is
developed: The kinetic multipeak polyhedrons. Distributions of this type are
expected to appear in each kinetic model with multidimensional instability as
universally, as Gaussian distribution appears for stable systems. The number of
possible relatively stable states of nonequilibrium system grows as $2^m$, and
the number of macroscopic parameters is in order $mn$, where $n$ is the
dimension of configuration space, and $m$ is the number of independent unstable
directions in this space. The elaborated class of closures and equations
pretends to describe the effects of "molecular individualism". This is the
third result.


  This paper extends Eigen's quasispecies equations to account for the
semiconservative nature of DNA replication. We solve the equations in the limit
of infinite sequence length for the simplest case of a static, sharply peaked
fitness landscape. We show that the error catastrophe occurs when $ \mu $, the
product of sequence length and per base pair mismatch probability, exceeds $ 2
\ln \frac{2}{1 + 1/k} $, where $ k > 1 $ is the first order growth rate
constant of the viable ``master'' sequence (with all other sequences having a
first-order growth rate constant of $ 1 $). This is in contrast to the result
of $ \ln k $ for conservative replication. In particular, as $ k \to \infty $,
the error catastrophe is never reached for conservative replication, while for
semiconservative replication the critical $ \mu $ approaches $ 2 \ln 2 $.
Semiconservative replication is therefore considerably less robust than
conservative replication to the effect of replication errors. We also show that
the mean equilibrium fitness of a semiconservatively replicating system is
given by $ k (2 e^{-\mu/2} - 1) $ below the error catastrophe, in contrast to
the standard result of $ k e^{-\mu} $ for conservative replication (derived by
Kimura and Maruyama in 1966).


  We propose an automated protocol for designing the energy landscape of a
protein energy function by optimizing its parameters. The parameters are
optimized so that not only the global minimum energy conformation becomes
native-like, but also the conformations distinct from the native structure have
higher energies than those close to the native one.
  We successfully apply our protocol to the parameter optimization of the UNRES
potential energy, using the training set of betanova, 1fsd, the 36-residue
subdomain of chicken villin headpiece (PDB ID 1vii), and the 10-55 residue
fragment of staphylococcal protein A (PDB ID 1bdd). The new protocol of the
parameter optimization shows better performance than earlier methods where only
the difference between the lowest energies of native-like and non-native
conformations was adjusted without considering various degrees of
native-likeness of the conformations. We also perform jackknife tests on other
proteins not included in the training set and obtain promising results. The
results suggest that the parameters we obtained using the training set of the
four proteins are transferable to other proteins to some extent.


  The inert polymer polyethylene glycol (PEG) induces a "bundling" phenomenon
in F-actin solutions when its concentration exceeds a critical onset value C_o.
Over a limited range of PEG molecular weight and ionic strength, C_o can be
expressed as a function of these two variables. The process is reversible, but
hysteresis is also observed in the dissolution of the bundles, with ionic
strength having a large influence. Additional actin filaments are able to join
previously formed bundles. Little, if any, polymer is associated with the
bundle structure.
  Continuum estimates of the Asakura-Oosawa depletion force, Coulomb repulsion,
and van der Waals potential are combined for a partial explanation of the
bundling effect and hysteresis. Conjectures are presented concerning the
apparent limit in bundle size.


  The free energy of inserting a protein into a membrane is determined by
considering the variation in the spectrum of thermal fluctuations in response
to the presence of a rigid inclusion. Both numerically and through a simple
analytical approximation, we find that the primary effect of fluctuations is to
reduce the effective surface tension, hampering the insertion at low surface
tension. Our results, which should also be relevant for membrane pores, suggest
(in contrast to classical nucleation theory) that a finite surface tension is
necessary to facilitate the opening of a pore.


  The present review deals with the computer simulation of biological ageing as
well as its demographic consequences for industrialized societies.


  The free energy of globular protein chain is considered to be a functional
defined on smooth curves in three dimensional Euclidean space. From the
requirement of geometrical invariance, together with basic facts on
conformation of helical proteins and dynamical characteristics of the protein
chains, we are able to determine, in a unique way, the exact form of the free
energy functional. Namely, the free energy density should be a linear function
of the curvature of curves on which the free energy functional is defined. This
model can be used, for example, in Monte Carlo simulations of exhaustive
searching the native stable state of the protein chain.


  The theoretical description of the evolution of cooperation presented by
Bergstrom based on assortative matching with partner choice allows to model the
population dynamics in a game of Nonrepetitive Prisoners Dilemma. In this paper
we present a short analysis of asymmetric effects brought into the game by self
knowledge of the participants, that is the knowledge of one's own strategy.
Within the same conceptual framework as introduced by Bergstrom inclusion of
selfknowledge leads to different behaviour of the assortativity index and
changes in payoffs for different strategies.


  We studied the electrical conductivity of DNA samples as function of the
number of DNA molecules. We showed that the insulating gap (no current at low
voltage) increases from ~1-2 V for bundles and large ropes to ~4-7 V for few
DNA molecules. From the distance dependent variation of the current, a unique
hopping distance of ~3 nm is calculated (polaron-hopping model) independently
of the number of DNA in the sample. The highly resistive behavior of the single
DNA is correlated with its flattened conformation on the surface (reduced
thickness, ~0.5-1.5 nm, compared to its nominal value, ~2 nm).