  Motivated by a biased diffusion of molecular motors with the bias dependent
on the state of the substrate, we investigate a random walk on a
one-dimensional lattice that contains weak links (called "bridges'') which are
affected by the walker. Namely, a bridge is destroyed with probability p when
the walker crosses it; the walker is not allowed to cross it again and this
leads to a directed motion. The velocity of the walker is determined
analytically for equidistant bridges. The special case of p=1 is more tractable
-- both the velocity and the diffusion constant are calculated for uncorrelated
locations of bridges, including periodic and random distributions.


  We study theoretically the entropic elasticity of a semi-flexible polymer,
such as DNA, confined to two dimensions. Using the worm-like-chain model we
obtain an exact analytical expression for the partition function of the polymer
pulled at one end with a constant force. The force-extension relation for the
polymer is computed in the long chain limit in terms of Mathieu characteristic
functions. We also present applications to the interaction between a
semi-flexible polymer and a nematic field, and derive the nematic order
parameter and average extension of the polymer in a strong field.


  The phase behavior of charged rods in the presence of inter-rod linkers is
studied theoretically as a model for the equilibrium behavior underlying the
organization of actin filaments by linker proteins in the cytoskeleton. The
presence of linkers in the solution modifies the effective inter-rod
interaction and can lead to inter-filament attraction. Depending on the
system's composition and physical properties such as linker binding energies,
filaments will either orient perpendicular or parallel to each other, leading
to network-like or bundled structures. We show that such a system can have one
of three generic phase diagrams, one dominated by bundles, another by networks,
and the third containing both bundle and network-like phases. The first two
diagrams can be found over a wide range of interaction energies, while the
third occurs only for a narrow range. These results provide theoretical
understanding of the classification of linker proteins as bundling proteins or
crosslinking proteins. In addition, they suggest possible mechanisms by which
the cell may control cytoskeletal morphology.


  Mapping a complex network of $N$coupled identical oscillators to a quantum
system, the nearest neighbor level spacing (NNLS) distribution is used to
identify collective chaos in the corresponding classical dynamics on the
complex network. The classical dynamics on an Erdos-Renyi network with the
wiring probability $p_{ER} \le \frac{1}{N}$ is in the state of collective
order, while that on an Erdos-Renyi network with $p_{ER} > \frac{1}{N}$ in the
state of collective chaos. The dynamics on a WS Small-world complex network
evolves from collective order to collective chaos rapidly in the region of the
rewiring probability $p_r \in [0.0,0.1]$, and then keeps chaotic up to $p_r =
1.0$. The dynamics on a Growing Random Network (GRN) is in a special state
deviates from order significantly in a way opposite to that on WS small-world
networks. Each network can be measured by a couple values of two parameters
$(\beta ,\eta)$.


  We consider ion transport through protein ion channels in lipid membranes and
water-filled nanopores in silicon films. It is known that, due to the large
ratio of dielectric constants of water and the surrounding material, an ion
placed inside the channel faces a large electrostatic self-energy barrier. The
barrier leads to an exponentially large resistance of the channel. We study
reduction of the electrostatic barrier by immobile charges located on the
internal walls of the channel. We show that the barrier practically vanishes
already at relatively small concentration of wall charges.


  We study a discrete spatial model for invasive allele spread in which two
alleles compete preemptively, initially only the "residents" (weaker
competitors) being present. We find that the spread of the advantageous
mutation is well described by homogeneous nucleation; in particular, in large
systems the time-dependent global density of the resident allele is well
approximated by Avrami's law.


  We model the dynamics of a processive or rotary molecular motor using a
renewal processes, in line with the work initiated by Svoboda, Mitra and Block.
We apply a functional technique to compute different types of multiple-time
correlation functions of the renewal process, which have applications to
bead-assay experiments performed both with processive molecular motors, such as
myosin V and kinesin, and rotary motors, such as F1-ATPase.


  We propose a new statistical mechanics model for the melting transition of
DNA. Base pairing and stacking are treated as separate degrees of freedom, and
the interplay between pairing and stacking is described by a set of local rules
which mimic the geometrical constraints in the real molecule. This microscopic
mechanism intrinsically accounts for the cooperativity related to the free
energy penalty of bubble nucleation. The model describes both the unpairing and
unstacking parts of the spectroscopically determined experimental melting
curves. Furthermore, the model explains the observed temperature dependence of
the effective thermodynamic parameters used in models of the nearest neighbor
(NN) type. We compute the partition function for the model through the transfer
matrix formalism, which we also generalize to include non local chain entropy
terms. This part introduces a new parametrization of the Yeramian-like transfer
matrix approach to the Poland-Scheraga description of DNA melting. The model is
exactly solvable in the homogeneous thermodynamic limit, and we calculate all
observables without use of the grand partition function. As is well known,
models of this class have a first order or continuous phase transition at the
temperature of complete strand separation depending on the value of the
exponent of the bubble entropy.


  The effect of sequence heterogeneity on the dynamics of molecular motors is
reviewed and analyzed using a set of recently introduced lattice models. First,
we review results for the influence of heterogenous tracks such as a
single-strand of DNA or RNA on the dynamics of the motors. We stress how the
predicted behavior might be observed experimentally in anomalous drift and
diffusion of motors over a wide range of parameters near the stall force and
discuss the extreme limit of strongly biased motors with one-way hopping. We
then consider the dynamics in an environment containing a variety of different
fuels which supply chemical energy for the motor motion, either on a
heterogeneous or on a periodic track. The results for motion along a periodic
track are relevant to kinesin motors in a solution with a mixture of different
nucleotide triphosphate fuel sources.


  Individual-based models of chemical or biological dynamics usually consider
individual entities diffusing in space and performing a birth-death type
dynamics. In this work we study the properties of a model in this class where
the birth dynamics is mediated by the local, within a given distance, density
of particles. Groups of individuals are formed in the system and in this paper
we concentrate on the study of the properties of these clusters (lifetime,
size, and collective diffusion). In particular, in the limit of the interaction
distance approaching the system size, a unique cluster appears which helps to
understand and characterize the clustering dynamics of the model.


  We introduce a novel model for DNA condensation (whip-toroid transition)
using the path integral method in the framework of the non-linear sigma model
on a line segment. We show that some of its classical configurations exhibit
toroidal forms, and the system has phase transitions from a whip to toroidal
phases with a parameter $c=W/(2l) (L/(2\pi))^2$. We also discuss stability and
finite size effect on these states.


  The finely tuned structures of membrane channels allow selective passage of
ions through the available aqueous pores. In order to understand channel
function, it is crucial to locate the pore and study its physical and chemical
properties. Recently obtained X-ray crystal structures of bacterial chloride
channel homologues reveal a complicated topology with curvilinear pores. The
commonly used HOLE program encounters difficulties in studying such pores. Here
we propose a new pore-searching algorithm (TransPath) which uses the
Configurational Bias Monte Carlo (CBMC) method to generate transmembrane
trajectories driven by both geometric and electrostatic features. The
trajectories are binned into groups determined by a vector distance criterion.
  From each group, a representative trajectory is selected based on the
Rosenbluth weight, and the geometrically optimal path is obtained by simulated
annealing. Candidate ion pathways can then be determined by analysis of the
radius and potential profiles. The proposed method and its implementation are
illustrated using the bacterial KcsA potassium channel as an example. The
procedure is then applied to the more complex structures of the bacterial E.
coli ClC channel homologues.


  Polymerization dynamics of single actin filaments is investigated
theoretically using a stochastic model that takes into account the hydrolysis
of ATP-actin subunits, the geometry of actin filament tips, the lateral
interactions between the monomers as well as the processes at both ends of the
polymer. Exact analytical expressions are obtained for a mean growth velocity
and for dispersion in length fluctuations. It is found that the ATP hydrolysis
has a strong effect on dynamic properties of single actin filaments. At high
concentrations of free actin monomers the mean size of unhydrolyzed ATP-cap is
very large, and the dynamics is governed by association/dissociation of
ATP-actin subunits. However, at low concentrations the size of the cap becomes
finite, and the dissociation of ADP-actin subunits makes a significant
contribution to overall dynamics. Actin filament length fluctuations reach the
maximum at the boundary between two dynamic regimes, and this boundary is
always larger than the critical concentration. Random and vectorial mechanisms
of hydrolysis are compared, and it is found that they predict qualitatively
similar dynamic properties. The possibility of attachment and detachment of
oligomers is also discussed. Our theoretical approach is successfully applied
to analyze the latest experiments on the growth and length fluctuations of
individual actin filaments.


  Membrane tubes (spicules) arise in cells, or artificial membranes, in the
nonlinear deformation regime due to, e.g. the growth of microtubules, actin
filaments or sickle hemoglobin fibers towards a membrane. We calculate the
axial force exerted by the cylindrical membrane tube, and its average radius,
by taking into account steric interactions between the fluctuating membrane and
the enclosed rod. The force required to confine a fluctuating membrane near the
surface of the enclosed rod diverges as the separation approaches zero. This
results in a smooth crossover of the axial force between a square root and a
linear dependence on the membrane tension as the tension increases and the tube
radius shrinks. This crossover can occur at the most physiologically relevant
membrane tensions. Our work may be important in (i) interpreting experiments in
which axial force is related to the tube radius or membrane tension (ii)
dynamical theories for biopolymer growth in narrow tubes where these
fluctuation effects control the tube radius.


  Enyzme kinetics are cyclic. We study a Markov renewal process model of
single-enzyme turnover in nonequilibrium steady-state (NESS) with sustained
concentrations for substrates and products. We show that the forward and
backward cycle times have idential non-exponential distributions:
$\QQ_+(t)=\QQ_-(t)$. This equation generalizes the Haldane relation in
reversible enzyme kinetics. In terms of the probabilities for the forward
($p_+$) and backward ($p_-$) cycles, $k_BT\ln(p_+/p_-)$ is shown to be the
chemical driving force of the NESS, $\Delta\mu$. More interestingly, the moment
generating function of the stochastic number of substrate cycle $\nu(t)$,
$<e^{-\lambda\nu(t)}>$ follows the fluctuation theorem in the form of
Kurchan-Lebowitz-Spohn-type symmetry. When $\lambda$ = $\Delta\mu/k_BT$, we
obtain the Jarzynski-Hatano-Sasa-type equality: $<e^{-\nu(t)\Delta\mu/k_BT}>$
$\equiv$ 1 for all $t$, where $\nu\Delta\mu$ is the fluctuating chemical work
done for sustaining the NESS. This theory suggests possible methods to
experimentally determine the nonequilibrium driving force {\it in situ} from
turnover data via single-molecule enzymology.


  We present the results of analytic calculations and numerical simulations of
the behaviour of a new class of chain molecules which we call thick polymers.
The concept of the thickness of such a polymer, viewed as a tube, is
encapsulated by a special three body interaction and impacts on the behaviour
both locally and non-locally. When thick polymers undergo compaction due to an
attractive self-interaction, we find a new type of phase transition between a
compact phase and a swollen phase at zero temperature on increasing the
thickness. In the vicinity of this transition, short tubes form space filling
helices and sheets as observed in protein native state structures. Upon
increasing the chain length, or the number of chains, we numerically find a
crossover from secondary structure motifs to a quite distinct class of
structures akin to the semi-crystalline phase of polymers or amyloid fibers in
polypeptides.


  A method to reconstruct the energy landscape of small peptides is presented
with reference to a 2d off--lattice model. The starting point is a statistical
analysis of the configurational distances between generic minima and directly
connected pairs (DCP). As the mutual distance of DCP is typically much smaller
than that of generic pairs, a metric criterion can be established to identify
the great majority of DCP. Advantages and limits of this approach are
thoroughly analyzed for three different heteropolymeric chains. A funnel--like
structure of the energy landscape is found in all of the three cases, but the
escape rates clearly reveal that the native configuration is more easily
accessible (and is significantly more stable) for the sequence that is expected
to behave as a real protein.


  The Bak-Sneppen model displaying punctuated equilibria in biological
evolution is studied on random complex networks. By using the rate equation and
the random walk approaches, we obtain the analytic solution of the fitness
threshold $x_c$ to be 1/(<k>_f+1), where <k>_f=<k^2>/<k> (=<k>) in the quenched
(annealed) updating case, where <k^n> is the n-th moment of the degree
distribution. Thus, the threshold is zero (finite) for the degree exponent
\gamma <3 (\gamma > 3) for the quenched case in the thermodynamic limit. The
theoretical value x_c fits well to the numerical simulation data in the
annealed case only. Avalanche size, defined as the duration of successive
mutations below the threshold, exhibits a critical behavior as its distribution
follows a power law, P_a(s) ~ s^{-3/2}.


  We study electronic transport in long DNA chains using the tight-binding
approach for a ladder-like model of DNA. We find insulating behavior with
localizaton lengths xi ~ 25 in units of average base-pair seperation.
Furthermore, we observe small, but significant differences between lambda-DNA,
centromeric DNA, promoter sequences as well as random-ATGC DNA.


  The effect of interactions on dynamics of coupled motor proteins is
investigated theoretically. A simple stochastic discrete model, that allows to
calculate explicitly the dynamic properties of the system, is developed. It is
shown that there are two dynamic regimes, depending on the interaction between
the particles. For strong interactions the motor proteins move as one tight
cluster, while for weak interactions there is no correlation in the motion of
the proteins, and the particle separation increases steadily with time. The
boundary between two dynamic phases is specified by a critical interaction that
has a non-zero value only for the coupling of the asymmetric motor proteins,
and it depends on the temperature and the transitions rates. At the critical
interaction there is a change in a slope for the mean velocities and a
discontinuity in the dispersions of the motor proteins as a function of the
interaction energy.


  We study the translocation process of a polymer in the absence of external
fields for various pore diameters $b$ and membrane thickness $L$. The polymer
performs Rouse and reptation dynamics. The mean translocation time $<\tau_t>$
that the polymer needs to escape from a cell, and the mean dwell time
$<\tau_d>$ that the polymer spends in the pore during the translocation
process, obey scaling relations in terms of the polymer length $N$, $L$ and
$b/R_g$, where $R_g$ is the radius of gyration for the polymer. We explain
these relations using simple arguments based on polymer dynamics and the
equilibrium properties of polymers.


  We propose a model for evolving networks by merging building blocks
represented as complete graphs, reminiscent of modules in biological system or
communities in sociology. The model shows power-law degree distributions,
power-law clustering spectra and high average clustering coefficients
independent of network size. The analytical solutions indicate that a degree
exponent is determined by the ratio of the number of merging nodes to that of
all nodes in the blocks, demonstrating that the exponent is tunable, and are
also applicable when the blocks are classical networks such as
Erd\H{o}s-R\'enyi or regular graphs. Our model becomes the same model as the
Barab\'asi-Albert model under a specific condition.


  The Penna bit-string model successfully encompasses many phenomena of
population evolution, including inheritance, mutation, evolution and ageing. If
we consider social interactions among individuals in the Penna model, the
population will form a complex network. In this paper, we first modify the
Verhulst factor to control only the birth rate, and introduce activity-based
preferential reproduction of offspring in the Penna model. The social
interactions among individuals are generated by both inheritance and
activity-based preferential increase. Then we study the properties of the
complex network generated by the modified Penna model. We find that the
resulting complex network has a small-world effect and the assortative mixing
property.


  We study the growth of a directed network, in which the growth is constrained
by the cost of adding links to the existing nodes. We propose a new
preferential-attachment scheme, in which a new node attaches to an existing
node i with probability proportional to 1/k_i, where k_i is the number of
outgoing links at i. We calculate the degree distribution for the outgoing
links in the asymptotic regime (t->infinity), both analytically and by Monte
Carlo simulations. The distribution decays like k c^k/Gamma(k) for large k,
where c is a constant. We investigate the effect of this
preferential-attachment scheme, by comparing the results to an equivalent
growth model with a degree-independent probability of attachment, which gives
an exponential outdegree distribution. Also, we relate this mechanism to simple
food-web models by implementing it in the cascade model. We show that the
low-degree preferential-attachment mechanism breaks the symmetry between in-
and outdegree distributions in the cascade model. It also causes a faster decay
in the tails of the outdegree distributions for both our network growth model
and the cascade model.


  The dielectrophoresis method for trapping and attaching nanoscale
double-stranded DNA between nanoelectrodes was developed. The method gives a
high yield of trapping single or a few molecules only which enables transport
measurements at the single molecule level. Electrical conductivity of
individual 140-nm-long DNA molecules was measured, showing insulating behavior
in dry conditions. In contrast, clear enhancement of conductivity was observed
in moist conditions, relating to the interplay between the conformation of DNA
molecules and their conductivity.


  It appears that thermally activated DNA bubbles of different sizes play
central roles in important genetic processes. Here we show that the probability
for the formation of such bubbles is regulated by the number of soft AT pairs
in specific regions with lengths which at physiological temperatures are of the
order of (but not equal to) the size of the bubble. The analysis is based on
the Peyrard- Bishop-Dauxois model, whose equilibrium statistical properties
have been accurately calculated here with a transfer integral approach.


  We present a thermodynamically consistent mesoscopic model of protein
adsorption at liquid-solid interfaces. First describing the equilibrium state
under varying protein concentration of the solution and binding conditions, we
predict a non-trivial (non- monotonic) dependence of the experimentally
observable properties of the adsorbed layer (such as the surface density and
surface coverage) on these parameters. We subsequently proceed to develop a
dynamical model consistent with the equilibrium description, which
qualitatively reproduces known experimental phenomena and offers a promising
way of studying the exchange of the adsorbed proteins by the proteins of the
solution.


  We investigate the dynamics of polymer translocation through a nanopore under
an externally applied field using the 2D fluctuating bond model with
single-segment Monte Carlo moves. We concentrate on the influence of the field
strength $E$, length of the chain $N$, and length of the pore $L$ on forced
translocation. As our main result, we find a crossover scaling for the
translocation time $\tau $ with the chain length from $\tau \sim N^{2\nu}$ for
relatively short polymers to $\tau \sim N^{1 + \nu}$ for longer chains, where
$\nu $ is the Flory exponent. We demonstrate that this crossover is due to the
change in the dependence of the translocation velocity v on the chain length.
For relatively short chains $v \sim N^{- \nu}$, which crosses over to $v \sim
N^{- 1}$ for long polymers. The reason for this is that with increasing $N$
there is a high density of segments near the exit of the pore, which slows down
the translocation process due to slow relaxation of the chain. For the case of
a long nanopore for which $R_\parallel $, the radius of gyration $R_{g}$ along
the pore, is smaller than the pore length, we find no clear scaling of the
translocation time with the chain length. For large $N$, however, the
asymptotic scaling $\tau \sim N^{1 + \nu}$ is recovered. In this regime, $\tau
$ is almost independent of $L$. We have previously found that for a polymer,
which is initially placed in the middle of the pore, there is a minimum in the
escape time for $R_\parallel \approx L$. We show here that this minimum
persists for a weak fields $E$ such that $EL$ is less than some critical value,
but vanishes for large values of $EL$.


  A noise source model, consisting of a pulse sequence at random times with
memory, is presented. By varying the memory we can obtain variable randomness
of the stochastic process. The delay time between pulses, i. e. the noise
memory, produces different kinds of correlated noise ranging from white noise,
without delay, to quasi-periodical process, with delay close to the average
period of the pulses. The spectral density is calculated. This type of noise
could be useful to describe physical and biological systems where some delay is
present. In particular it could be useful in population dynamics. A simple
dynamical model for epidemiological infection with this noise source is
presented. We find that the time behavior of the illness depends on the noise
parameters. Specifically the amplitude and the memory of the noise affect the
number of infected people.


  We study the effect of the composition of the genetic sequence on the melting
temperature of double stranded DNA, using some simple analytically solvable
models proposed in the framework of the wetting problem. We review previous
work on disordered versions of these models and solve them when there were not
preexistent solutions. We check the solutions with Monte Carlo simulations and
transfer matrix numerical calculations. We present numerical evidence that
suggests that the logarithmic corrections to the critical temperature due to
disorder, previously found in RSOS models, apply more generally to ASOS and
continuous models. The agreement between the theoretical models and
experimental data shows that, in this context, disorder should be the crucial
ingredient of any model while other aspects may be kept very simple, an
approach that can be useful for a wider class of problems. Our work has also
implications for the existence of correlations in DNA sequences.


  The flow of motor proteins on a filamental track is modelled within the the
framework of lattice driven diffusive systems. Motors, considered as hopping
particles, perform a highly biased asymmetric exclusion process when bound to
the filament. With a certain rate, they detach from the filament and execute
unbiased random walk in the bulk which is considered as a closed cubic
compartment. Motors are injected (extracted) from the leftmost (rightmost) site
of the filament located along the symmetry axis of the compartment. We explore
the transport properties of this system and investigate the bulk-boundary
interplay on the system stationary states. It is shown that the detachment rate
notably affects the system properties. In particular and in contrast to ASEP,
it is shown that the density profile of bound particles exhibit different types
of non monotonic behaviours when the detachment rate varies. It is shown that
in certain situations, the density profile of the filament consists of
coexisting high and low regions.


  In this paper, we analyse large random Boolean networks in terms of a
constraint satisfaction problem. We first develop an algorithmic scheme which
allows to prune simple logical cascades and under-determined variables,
returning thereby the computational core of the network. Second we apply the
cavity method to analyse number and organisation of fixed points. We find in
particular a phase transition between an easy and a complex regulatory phase,
the latter one being characterised by the existence of an exponential number of
macroscopically separated fixed-point clusters. The different techniques
developed are reinterpreted as algorithms for the analysis of single Boolean
networks, and they are applied to analysis and in silico experiments on the
gene-regulatory networks of baker's yeast (saccaromices cerevisiae) and the
segment-polarity genes of the fruit-fly drosophila melanogaster.


  Based on our previous paper [cond-mat/0507477], we establish a general model
for the whip-toroid transitions of a semiflexible homopolymer chain using the
path integral method and the O(3) nonlinear sigma model on a line segment with
the local inextensibility constraint. We exactly solve the energy levels of
classical solutions, and show that some of its classical configurations exhibit
toroidal forms, and the system has phase transitions from a whip to toroidal
states with a conformation parameter $c=\frac{W}{2 l} (\frac{L}{2\pi})^2$. We
also discuss the stability of the toroid states and propose the low-energy
effective Green function. Finally, with the finite size effect on the toroid
states, predicted toroidal properties are successfully compared to experimental
results of DNA condensation.


  Molecular motors perform active movements along cytoskeletal filaments and
drive the traffic of organelles and other cargo particles in cells. In contrast
to the macroscopic traffic of cars, however, the traffic of molecular motors is
characterized by a finite walking distance (or run length) after which a motor
unbinds from the filament along which it moves. Unbound motors perform Brownian
motion in the surrounding aqueous solution until they rebind to a filament. We
use variants of driven lattice gas models to describe the interplay of their
active movements, the unbound diffusion, and the binding/unbinding dynamics. If
the motor concentration is large, motor-motor interactions become important and
lead to a variety of cooperative traffic phenomena such as traffic jams on the
filaments, boundary-induced phase transitions, and spontaneous symmetry
breaking in systems with two species of motors. If the filament is surrounded
by a large reservoir of motors, the jam length, i.e., the extension of the
traffic jams is of the order of the walking distance. Much longer jams can be
found in confined geometries such as tube-like compartments.


  We consider a theoretical model for membranes with adhesive receptors, or
stickers, that are actively switched between two conformational states. In
their 'on'-state, the stickers bind to ligands in an apposing membrane, whereas
they do not interact with the ligands in their 'off'-state. We show that the
adhesiveness of the membranes depends sensitively on the rates of the
conformational switching process. This dependence is reflected in a resonance
at intermediate switching rates, which can lead to large membrane separations
and unbinding. Our results may provide insights into novel mechanisms for the
controlled adhesion of biological or biomimetic membranes.


  We develop a phenomenological vector model of polar liquids capable to
describe aqueous interactions of macroscopic bodies. It is shown that a strong,
long-range and orientationally dependent interaction between macroscopic
objects appears as a result of competition between short-range (hydrogen
bonding) and the long-range dipole-dipole interactions of the solvent
molecules. Spontaneous polarization of molecular dipoles next to a hydrophobic
boundaries leads to formation of globally ordered network of hydrogen-bonded
molecules with ferroelectric properties. The proposed vector model naturally
describes topological excitations on the solute boundaries and can be used to
explain the hydrogen bonds networks and order-disorder phase transitions in the
hydration water layer.


  A realistic physical model for the so called rate promoting vibration (RPV)
at enzyme action is constructed. The origin of the RPV is assumed to be an
oscillating electric field produced by long-lived localized vibrational modes
in protein dynamics, namely, by the so called discrete breather (DB) in
secondary structure. The strength of interaction of the RPV with the reaction
coordinate is evaluated and its effect on the reaction acceleration is assessed
within the framework of modern theory for thermally activated escape rate at
periodic driving. We reveal the phenomenon of resonant activation in our model
elucidating why the frequency of the RPV in the range $100\div200 cm^{-1}$ was
chosen by the evolution of enzymes as an optimal one. The effect of the RPV on
the reaction acceleration is shown to vary from moderate one (up to
$10^3\div10^4$) in the case of three-site DB to enormous (up to $10^6\div10^8$)
in the case of five-site DB and thus can significantly contribute into enzyme
catalytic efficiency. Also the model is shown to be compatible with the known
functional dependence of enzymatic reaction rates on solvent viscosity.


  We consider the evolution of a population of fixed size with no selection.
The number of generations $G$ to reach the first common ancestor evolves in
time. This evolution can be described by a simple Markov process which allows
one to calculate several characteristics of the time dependence of $G$. We also
study how $G$ is correlated to the genetic diversity.


  By resorting to the thick-chain model we discuss how the stretching response
of a polymer is influenced by the self-avoidance entailed by its finite
thickness. The characterization of the force versus extension curve for a thick
chain is carried out through extensive stochastic simulations. The
computational results are captured by an analytic expression that is used to
fit experimental stretching measurements carried out on DNA and single-stranded
RNA (poly-U) in various solutions. This strategy allows us to infer the
apparent diameter of two biologically-relevant polyelectrolytes, namely DNA and
poly-U, for different ionic strengths. Due to the very different degree of
flexibility of the two molecules, the results provide insight into how the
apparent diameter is influenced by the interplay between the
(solution-dependent) Debye screening length and the polymers' ``bare''
thickness. For DNA, the electrostatic contribution to the effective radius,
$\Delta$, is found to be about 5 times larger than the Debye screening length,
consistently with previous theoretical predictions for highly-charged stiff
rods. For the more flexible poly-U chains the electrostatic contribution to
$\Delta$ is found to be significantly smaller than the Debye screening length.


  Loop formation between monomers in the interior of semiflexible chains
describes elementary events in biomolecular folding and DNA bending. We
calculate analytically the interior distance distribution function for
semiflexible chains using a mean-field approach. Using the potential of mean
force derived from the distance distribution function we present a simple
expression for the kinetics of interior looping by adopting Kramers theory. For
the parameters, that are appropriate for DNA, the theoretical predictions in
comparison to the case are in excellent agreement with explicit Brownian
dynamics simulations of worm-like chain (WLC) model. The interior looping times
($\tau_{IC}$) can be greatly altered in cases when the stiffness of the loop
differs from that of the dangling ends. If the dangling end is stiffer than the
loop then $\tau_{IC}$ increases for the case of the WLC with uniform
persistence length. In contrast, attachment of flexible dangling ends enhances
rate of interior loop formation. The theory also shows that if the monomers are
charged and interact via screened Coulomb potential then both the cyclization
($\tau_c$) and interior looping ($\tau_{IC}$) times greatly increase at low
ionic concentration. Because both $\tau_c$ and $\tau_{IC}$ are determined
essentially by the effective persistence length ($l_p^{(R)}$) we computed
$l_p^{(R)}$ by varying the range of the repulsive interaction between the
monomers. For short range interactions $l_p^{(R)}$ nearly coincides with the
bare persistence length which is determined largely by the backbone chain
connectivity. This finding rationalizes the efficacy of describing a number of
experimental observations (response of biopolymers to force and cyclization
kinetics) in biomolecules using WLC model with an effective persistence length.


  Differential microcalorimetric measurements have been performed in aqueous
dispersions of dimyristoylphosphatidylcholine (DMPC) multilamellar vesicles in
a wide range of temperatures covering the whole transition between the lamellar
gel and lamellar fluid states (the chain-melting/ordering transition). The
shape of calorimetric curves was analyzed in a temperature range some distance
away but close to the chain-ordering transition point where small nuclei of gel
(solid state) are formed. In this range, where the fraction of the "new-state"
lipid is small, the nucleation process can be considered independent of the
interlayer interactions and determined mainly by the lateral interactions. The
small-scale nucleation was analyzed in terms of a two-dimensional Ising-like
lattice model. The gel-fluid contact energy related to the critical temperature
(in terms of Ising model) was estimated for DMPC studied here and
dipalmitoylphosphatidylcholine studied earlier. The contact energy was found to
be not high enough to provide the discontinuous, first-order character of
transition. Therefore, the signs of first-order character observed in the
vesicles made of these lipids are not due to lateral but rather interlayer
3D-interactions. The extrapolation to longer lipids shows that the transition
discontinuity inherent to a lipid layer (i.e., determined by the lateral
interactions themselves) is expected for the saturated phosphatidylcholines
having more than 16-18 C-atoms per chain. Interestingly, it is the saturated
chains of this length that are the most abundant in biological membranes.
Probably, the biological membrane "prefers" to be near the critical state where
the system is the most responsive to physical actions.


  We have proposed an efficient parameterization method for a recent variant of
the Gay-Berne potential for dissimilar and biaxial particles and demonstrated
it for a set of small organic molecules. Compared to the previously proposed
coarse-grained models, the new potential exhibits a superior performance in
close contact and large distant interactions. The repercussions of thermal
vibrations and elasticity has been studied through a statistical method. The
study justifies that the potential of mean force is representable with the same
functional form, extending the application of this coarse-grained description
to a broader range of molecules. Moreover, the advantage of employing
coarse-grained models over truncated atomistic summations with large distance
cutoffs has been briefly studied.


  Proteins can locate their specific targets on DNA up to two orders of
magnitude faster than the Smoluchowski three-dimensional diffusion rate. This
happens due to non-specific adsorption of proteins to DNA and subsequent
one-dimensional sliding along DNA. We call such one-dimensional route towards
the target "antenna". We studied the role of the dispersion of nonspecific
binding energies within the antenna due to quasi random sequence of natural
DNA. Random energy profile for sliding proteins slows the searching rate for
the target. We show that this slowdown is different for the macroscopic and
mesoscopic antennas.


  We study the geometric properties of the energy landscape of coarse-grained,
off-lattice models of polymers by endowing the configuration space with a
suitable metric, depending on the potential energy function, such that the
dynamical trajectories are the geodesics of the metric. Using numerical
simulations, we show that the fluctuations of the curvature clearly mark the
folding transition, and that this quantity allows to distinguish between
polymers having a protein-like behavior (i.e., that fold to a unique
configuration) and polymers which undergo a hydrophobic collapse but do not
have a folding transition. These geometrical properties are defined by the
potential energy without requiring any prior knowledge of the native
configuration.


  We analyze the nonequilibrium dynamics of single inextensible semiflexible
biopolymers as stretching forces are applied at the ends. Based on different
(contradicting) heuristic arguments, various scaling laws have been proposed
for the propagation speed of the backbone tension which is induced in response
to stretching. Here, we employ a newly developed unified theory to
systematically substantiate, restrict, and extend these approaches. Introducing
the practically relevant scenario of a chain equilibrated under some
prestretching force $f_\text{pre}$ that is suddenly exposed to a different
external force $f_\text{ext}$ at the ends, we give a concise physical
explanation of the underlying relaxation processes by means of an intuitive
blob picture. We discuss the corresponding intermediate asymptotics, derive
results for experimentally relevant observables, and support our conclusions by
numerical solutions of the coarse-grained equations of motion for the tension.


  We study the elasticity of random fiber networks. Starting from a microscopic
picture of the non-affine deformation fields we calculate the macroscopic
elastic moduli both in a scaling theory and a self-consistent effective medium
theory. By relating non-affinity to the low-energy excitations of the network
(``floppy-modes'') we achieve a detailed characterization of the non-affine
deformations present in fibrous networks.


  Biological macromolecules have complex and non-trivial energy landscapes,
endowing them a unique conformational adaptability and diversity in function.
Hence, understanding the processes of elasticity and dissipation at the
nanoscale is important to molecular biology and also emerging fields such as
nanotechnology. Here we analyse single molecule fluctuations in an atomic force
microscope (AFM) experiment using a generic model of biopolymer viscoelasticity
that importantly includes sources of local `internal' conformational
dissipation. Comparing two biopolymers, dextran and cellulose, polysaccharides
with and without the well-known `chair-to-boat' transition, reveals a signature
of this simple conformational change as minima in both the elasticity and
internal friction around a characteristic force. A calculation of two-state
populations dynamics offers a simple explanation in terms of an elasticity
driven by the entropy, and friction by barrier-controlled hopping, of
populations on a landscape. The microscopic model, allows quantitative mapping
of features of the energy landscape, revealing unexpectedly slow dynamics,
suggestive of an underlying roughness to the free energy.


  We propose a model to explain finite-size effects in intracellular
microrheology observed in experiments. The constrained dynamics of the
particles in the intracellular medium, treated as a viscoelastic medium, is
described by means of a diffusion equation in which interactions of the
particles with the cytoskeleton are modelled by a harmonic force. The model
reproduces the observed power-law behavior of the mean-square displacement in
which the exponent depends on the ratio between
particle-to-cytoskeleton-network sizes.


  The current-voltage characteristics of the alpha-Hemolysin protein pore
during the passage of single-stranded DNA under varying ionic strength, C, are
studied experimentally. We observe strong blockage of the current, weak
super-linear growth of the current as a function of voltage, and a minimum of
the current as a function of C. These observations are interpreted as the
result of the ion electrostatic self-energy barrier originating from the large
difference in the dielectric constants of water and the lipid bilayer. The
dependence of DNA capture rate on C also agrees with our model.


  We investigate theoretically the fluid-crystal coexistence of solutions of
globular charged nanoparticles like proteins and inorganic colloids. The
thermodynamic properties of the fluid phase are computed via the optimized
Baxter model. This is done specifically for lysozyme and silicotungstates for
which the bare adhesion parameters are evaluated via the experimental second
virial coefficients. The electrostatic free energy of the crystal is
approximated by supposing the cavities in the interstitial phase between the
particles are spherical in form. In the salt-free case a Poisson-Boltzmann
equation is solved to calculate the effective charge on a particle and a Donnan
approximation is used to derive the chemical potential and osmotic pressure in
the presence of salt. The coexistence data of lysozyme and silicotungstates are
analyzed within this scheme, especially with regard to the ionic-strength
dependence of the chemical potentials. The latter agree within the two phases
provided some upward adjustment of the effective charge is allowed for.


  The co-evolution of network topology and dynamics is studied in an
evolutionary Boolean network model that is a simple model of gene regulatory
network. We find that a critical state emerges spontaneously resulting from
interplay between topology and dynamics during the evolution. The final evolved
state is shown to be independent of initial conditions. The network appears to
be driven to a random Boolean network with uniform in-degree of two in the
large network limit. However, for biologically realized network sizes,
significant finite-size effects are observed including a broad in-degree
distribution and an average in-degree connection between two and three. These
results may be important for explaining properties of gene regulatory networks.


  Equilibrium phase transitions are associated with rearrangements of minima of
a (Lagrangian) potential. Treatment of non-equilibrium systems requires
doubling of degrees of freedom, which may be often interpreted as a transition
from the ``coordinate'' to the ``phase'' space representation. As a result, one
has to deal with the Hamiltonian formulation of the field theory instead of the
Lagrangian one. We suggest a classification scheme of phase transitions in
reaction-diffusion models based on the topology of the phase portraits of
corresponding Hamiltonians. In models with an absorbing state such a topology
is fully determined by intersecting curves of zero ``energy''. We identify four
families of topologically distinct classes of phase portraits stable upon RG
transformations.


  Effect of bending anisotropy on a planar DNA loop, using energy minimization
and neglecting entropic effects, is studied. We show that the anisotropy
results in polygonal shape of the loop and increasing the anisotropy makes the
edges sharper. Calculating the energy of such a loop lets us to find effective
persistence length as the geometrical mean of hard and soft rigidities, which
is quite different from harmonic mean for an unconstrained long DNA.


  We derive exact expressions for the forces and torques between biaxial
molecules interacting via the RE-squared potential, a recent variant of the
Gay-Berne potential. Moreover, efficient routines have been provided for rigid
body MD simulations, resulting in 1.6 times speedup compared to the two-point
finite difference approach. It has also been shown that the time cost of a MD
simulation will be almost equal to a similar MC simulation, making use of the
provided routines.


  We present Monte Carlo simulation results for square-well homopolymers at a
series of bond lengths. Although the model contains only isotropic pairwise
interactions, under appropriate conditions this system shows spontaneous chiral
symmetry breaking, where the chain exists in either a left- or a right-handed
helical structure. We investigate how this behavior depends upon the ratio
between bond length and monomer radius.


  We study the threshold control protocol for a collective flashing ratchet. In
particular, we analyze the dependence of the current on the values of the
thresholds. We have found analytical expressions for the small threshold
dependence both for the few and for the many particle case. For few particles
the current is a decreasing function of the thresholds, thus, the maximum
current is reached for zero thresholds. In contrast, for many particles the
optimal thresholds have a nonzero finite value. We have numerically checked the
relation that allows to obtain the optimal thresholds for an infinite number of
particles from the optimal period of the periodic protocol. These optimal
thresholds for an infinite number of particles give good results for many
particles. In addition, they also give good results for few particles due to
the smooth dependence of the current up to these threshold values.


  In a variety of contexts, physicists study complex, nonlinear models with
many unknown or tunable parameters to explain experimental data. We explain why
such systems so often are sloppy; the system behavior depends only on a few
`stiff' combinations of the parameters and is unchanged as other `sloppy'
parameter combinations vary by orders of magnitude. We contrast examples of
sloppy models (from systems biology, variational quantum Monte Carlo, and
common data fitting) with systems which are not sloppy (multidimensional linear
regression, random matrix ensembles). We observe that the eigenvalue spectra
for the sensitivity of sloppy models have a striking, characteristic form, with
a density of logarithms of eigenvalues which is roughly constant over a large
range. We suggest that the common features of sloppy models indicate that they
may belong to a common universality class. In particular, we motivate focusing
on a Vandermonde ensemble of multiparameter nonlinear models and show in one
limit that they exhibit the universal features of sloppy models.


  From the underlying Master equations we derive one-dimensional stochastic
processes that describe generalized ensemble simulations as well as tempering
(simulated and parallel) simulations. The representations obtained are either
in the form of a one-dimensional Fokker-Planck equation or a hopping process on
a one-dimensional chain. In particular, we discuss the conditions under which
these representations are valid approximate Markovian descriptions of the
random walk in order parameter or control parameter space. They allow a unified
discussion of the stationary distribution on, as well as of the stationary flow
across each space. We demonstrate that optimizing the flow is equivalent to
minimizing the first passage time for crossing the space, and discuss the
consequences of our results for optimizing simulations. Finally, we point out
the limitations of these representations under conditions of broken ergodicity.


  When mixed together, DNA and polyaminoamide (PAMAM) dendrimers form fibers
that condense into a compact structure. We use optical tweezers to pull
condensed fibers and investigate the decondensation transition by measuring
force-extension curves (FECs). A characteristic plateau force (around 10 pN)
and hysteresis between the pulling and relaxation cycles are observed for
different dendrimer sizes, indicating the existence of a first-order transition
between two phases (condensed and extended) of the fiber. The fact that we can
reproduce the same FECs in the absence of additional dendrimers in the buffer
medium indicates that dendrimers remain irreversibly bound to the DNA backbone.
Upon salt variation FECs change noticeably confirming that electrostatic forces
drive the condensation transition. Finally, we propose a simple model for the
decondensing transition that qualitatively reproduces the FECs and which is
confirmed by AFM images.


  A simplified form of the time dependent evolutionary dynamics of a
quasispecies model with a rugged fitness landscape is solved via a mapping onto
a random flux model whose asymptotic behavior can be described in terms of a
random walk. The statistics of the number of changes of the dominant genotype
from a finite set of genotypes are exactly obtained confirming existing
conjectures based on numerics.


  We develop a spectral formulation and a stationary WKB approximation for
calculating the probabilities of rare events (large deviations from the mean)
in systems of reacting particles with infinite-range interaction, describable
by a master equation. We compare the stationary WKB approximation with a recent
time-dependent semiclassical approximation developed, for the same class of
problems, by Elgart and Kamenev. As a benchmark we use an exactly solvable
problem of the binary annihilation reaction 2A -> 0.


  We apply Kramers theory to investigate the dissociation of multiple bonds
under mechanical force and interpret experimental results for the
unfolding/refolding force distributions of an RNA hairpin pulled at different
loading rates using laser tweezers. We identify two different kinetic regimes
depending on the range of forces explored during the unfolding and refolding
process. The present approach extends the range of validity of the two-states
approximation by providing a theoretical framework to reconstruct free-energy
landscapes and identify force-induced structural changes in molecular
transition states using single molecule pulling experiments. The method should
be applicable to RNA hairpins with multiple kinetic barriers.


  We derive exact Langevin-type equations governing quasispecies dynamics. The
inherent multiplicative noise has both real and imaginary parts. The numerical
simulation of the underlying complex stochastic partial differential equations
is carried out employing the Cholesky decomposition for the noise covariance
matrix. This noise produces unavoidable spatio-temporal density fluctuations
about the mean field value. In two dimensions, the fluctuations are suppressed
only when the diffusion time scale is much smaller than the amplification time
scale for the master species.


  We systematically study the design of simple patchy sphere models that
reversibly self-assemble into monodisperse icosahedral clusters. We find that
the optimal patch width is a compromise between structural specificity (the
patches must be narrow enough to energetically select the desired clusters) and
kinetic accessibility (they must be sufficiently wide to avoid kinetic traps).
Similarly, for good yields the temperature must be low enough for the clusters
to be thermodynamically stable, but the clusters must also have enough thermal
energy to allow incorrectly formed bonds to be broken. Ordered clusters can
form through a number of different dynamic pathways, including direct
nucleation and indirect pathways involving large disordered intermediates. The
latter pathway is related to a reentrant liquid-to-gas transition that occurs
for intermediate patch widths upon lowering the temperature. We also find that
the assembly process is robust to inaccurate patch placement up to a certain
threshold, and that it is possible to replace the five discrete patches with a
single ring patch with no significant loss in yield.


  We consider a simplified model of protein folding, with binary degrees of
freedom, whose equilibrium thermodynamics is exactly solvable. Based on this
exact solution, the kinetics is studied in the framework of a local equilibrium
approach, for which we prove that (i) the free energy decreases with time, (ii)
the exact equilibrium is recovered in the infinite time limit, and (iii) the
folding rate is an upper bound of the exact one. The kinetics is compared to
the exact one for a small peptide and to Monte Carlo simulations for a longer
protein, then rates are studied for a real protein and a model structure.


  We study a stochastic lattice predator-prey system by means of Monte Carlo
simulations that do not impose any restrictions on the number of particles per
site, and discuss the similarities and differences of our results with those
obtained for site-restricted model variants. In accord with the classic
Lotka-Volterra mean-field description, both species always coexist in two
dimensions. Yet competing activity fronts generate complex, correlated
spatio-temporal structures. As a consequence, finite systems display transient
erratic population oscillations with characteristic frequencies that are
renormalized by fluctuations. For large reaction rates, when the processes are
rendered more local, these oscillations are suppressed. In contrast with
site-restricted predator-prey model, we observe species coexistence also in one
dimension. In addition, we report results on the steady-state prey age
distribution.


  Game theory is one of the key paradigms behind many scientific disciplines
from biology to behavioral sciences to economics. In its evolutionary form and
especially when the interacting agents are linked in a specific social network
the underlying solution concepts and methods are very similar to those applied
in non-equilibrium statistical physics. This review gives a tutorial-type
overview of the field for physicists. The first three sections introduce the
necessary background in classical and evolutionary game theory from the basic
definitions to the most important results. The fourth section surveys the
topological complications implied by non-mean-field-type social network
structures in general. The last three sections discuss in detail the dynamic
behavior of three prominent classes of models: the Prisoner's Dilemma, the
Rock-Scissors-Paper game, and Competing Associations. The major theme of the
review is in what sense and how the graph structure of interactions can modify
and enrich the picture of long term behavioral patterns emerging in
evolutionary games.


  Dielectrophoresis can potentially be used as an efficient trapping tool in
the fabrication of molecular devices. For nanoscale objects, however, the
Brownian motion poses a challenge. We show that the use of carbon nanotube
electrodes makes it possible to apply relatively low trapping voltages and
still achieve high enough field gradients for trapping nanoscale objects, e.g.,
single molecules. We compare the efficiency and other characteristics of
dielectrophoresis between carbon nanotube electrodes and lithographically
fabricated metallic electrodes, in the case of trapping nanoscale DNA
molecules. The results are analyzed using finite element method simulations and
reveal information about the frequency dependent polarizability of DNA.


  We investigate the effect of randomness in both relationships and decisions
on the evolution of cooperation. Simulation results show, in such randomness'
presence, the system evolves to a more frequency cooperation state than in its
absence. Specifically, there is an optimal amount of randomness, which can
induce the highest level of cooperation. The mechanism of randomness promoting
cooperation resembles a coherence-resonance-like fashion, which could be of
particular interest in evolutionary game dynamics in economic, biological and
social systems.


  We study the effect of random forces on a double stranded DNA in unzipping
the two strands, analogous to the problem of an adsorbed polymer under a random
force. The ground state develops bubbles of various lengths as the random force
fluctuation is increased. The unzipping phase diagram is shown to be
drastically different from the pure case.


  Translocation of a single stranded DNA (ssDNA) through an alpha-hemolysin
channel in a lipid membrane driven by applied transmembrane voltage V was
extensively studied recently. While the bare charge of the ssDNA piece inside
the channel is approximately 12 (in units of electron charge) measurements of
different effective charges resulted in values between one and two. We explain
these challenging observations by a large self-energy of a charge in the narrow
water filled gap between ssDNA and channel walls, related to large difference
between dielectric constants of water and lipid, and calculate effective
charges of ssDNA. We start from the most fundamental stall charge $q_s$, which
determines the force $F_s= q_s V/L$ stalling DNA against the voltage V (L is
the length of the channel). We show that the stall charge $q_s$ is proportional
to the ion current blocked by DNA, which is small due to the self-energy
barrier. Large voltage V reduces the capture barrier which DNA molecule should
overcome in order to enter the channel by $|q_c|V$, where $q_c$ is the
effective capture charge. We expressed it through the stall charge $q_s$. We
also relate the stall charge $q_s$ to two other effective charges measured for
ssDNA with a hairpin in the back end: the charge $q_u$ responsible for
reduction of the barrier for unzipping of the hairpin and the charge $q_e$
responsible for DNA escape in the direction of hairpin against the voltage. At
small V we explain reduction of the capture barrier with the salt
concentration.


  We study a reaction model that presents stochastic resonance purely due to
internal noise. This means that the only source of fluctuations comes from the
discrete character of the reactants, and no more noises enter into the system.
Our analysis reveals that the phenomenon is highly complex, and that is
generated by the interplay of different stochasticity at the three fixed points
of a bistable system.


  When pulled along its axis, double-stranded DNA elongates abruptly at a force
of about 65 pN. Two physical pictures have been developed to describe this
overstretched state. The first proposes that strong forces induce a phase
transition to a molten state consisting of unhybridized single strands. The
second picture instead introduces an elongated hybridized phase, called S-DNA,
structurally and thermodynamically distinct from standard B-DNA. Little
thermodynamic evidence exists to discriminate directly between these competing
pictures. Here we show that within a microscopic model of DNA we can
distinguish between the dynamics associated with each. In experiment,
considerable hysteresis in a cycle of stretching and shortening develops as
temperature is increased. Since there are few possible causes of hysteresis in
a system whose extent is appreciable in only one dimension, such behavior
offers a discriminating test of the two pictures of overstretching. Most
experiments are performed upon nicked DNA, permitting the detachment
(`unpeeling') of strands. We show that the long-wavelength progression of the
unpeeled front generates hysteresis, the character of which agrees with
experiment only if we assume the existence of S-DNA. We also show that internal
melting (distinct from unpeeling) can generate hysteresis, the degree of which
is strongly dependent upon the nonextensive loop entropy of single-stranded
DNA.


  The nonequilibrium phase transition in the triplet-creation model is
investigated using critical spreading and the conservative diffusive contact
process. The results support the claim that at high enough diffusion the phase
transition becomes discontinuous. As the diffusion probability increases the
critical exponents change continuously from the ordinary directed percolation
(DP) class to the compact directed percolation (CDP). The fractal dimension of
the critical cluster, however, switches abruptly between those two universality
classes. Strong crossover effects in both methods make it difficult, if not
impossible, to establish the exact location of the tricritical point.


  The enigmatic stability of population oscillations within ecological systems
is analyzed. The underlying mechanism is presented in the framework of two
interacting species free to migrate between two spatial patches. It is shown
that that the combined effects of migration and noise cannot account for the
stabilization. The missing ingredient is the dependence of the oscillations'
frequency upon their amplitude; with that, noise-induced differences between
patches are amplified due to the frequency gradient. Migration among
desynchronized regions then stabilizes a "soft" limit cycle in the vicinity of
the homogenous manifold. A simple model of diffusively coupled oscillators
allows the derivation of quantitative results, like the functional dependence
of the desynchronization upon diffusion strength and frequency differences. The
oscillations' amplitude is shown to be (almost) noise independent. The results
are compared with a numerical integration of the marginally stable
Lotka-Volterra equations. An unstable system is extinction-prone for small
noise, but stabilizes at larger noise intensity.


  Stochastic dynamics of chemical reactions in a mutually repressing two-gene
circuit is numerically simulated. The circuit has a rich variety of different
states when the kinetic change of DNA status is slow. The stochastic switching
transition between those states are compared with the theoretical estimation of
the switching rate derived from the idea similar to the transition state
theory. Even though the circuit is kept far from equilibrium, the method gives
a consistent explanation of the switching kinetics for a wide range of
parameters. The transition state theory-like estimation, however, fails to
describe transitions involving the state which has the extremely small numbers
of protein molecules.


  I review single-molecule experiments (SME) in biological physics. Recent
technological developments have provided the tools to design and build
scientific instruments of high enough sensitivity and precision to manipulate
and visualize individual molecules and measure microscopic forces. Using SME it
is possible to: manipulate molecules one at a time and measure distributions
describing molecular properties; characterize the kinetics of biomolecular
reactions and; detect molecular intermediates. SME provide the additional
information about thermodynamics and kinetics of biomolecular processes. This
complements information obtained in traditional bulk assays. In SME it is also
possible to measure small energies and detect large Brownian deviations in
biomolecular reactions, thereby offering new methods and systems to scrutinize
the basic foundations of statistical mechanics. This review is written at a
very introductory level emphasizing the importance of SME to scientists
interested in knowing the common playground of ideas and the interdisciplinary
topics accessible by these techniques. The review discusses SME from an
experimental perspective, first exposing the most common experimental
methodologies and later presenting various molecular systems where such
techniques have been applied. I briefly discuss experimental techniques such as
atomic-force microscopy (AFM), laser optical tweezers (LOT), magnetic tweezers
(MT), biomembrane force probe (BFP) and single-molecule fluorescence (SMF). I
then present several applications of SME to the study of nucleic acids (DNA,
RNA and DNA condensation), proteins (protein-protein interactions, protein
folding and molecular motors). Finally, I discuss applications of SME to the
study of the nonequilibrium thermodynamics of small systems and the
experimental verification of fluctuation theorems. I conclude with a discussion
of open questions and future perspectives.


  Dielectrophoretic trapping of six different DNA fragments, sizes varying from
the 27 to 8416 bp, has been studied using confocal microscopy. The effect of
the DNA length and the size of the constriction between nanoscale fingertip
electrodes on the trapping efficiency have been investigated. Using finite
element method simulations in conjunction with the analysis of the experimental
data, the polarizabilities of the different size DNA fragments have been
calculated for different frequencies. Also the immobilization of trapped
hexanethiol- and DTPA-modified 140 nm long DNA to the end of gold
nanoelectrodes was experimentally quantified and the observations were
supported by density functional theory calculations.


  We apply a theorem of Wick to rewrite certain classes of exponential measures
on random graphs as integrals of Feynman-Gibbs type, on the real line. The
analytic properties of these measures can then be studied in terms of phase
transitions; spaces of scale-free trees are a particularly interesting example.


  We suggest a general spectral method for calculating statistics of multi-step
birth-death processes and chemical reactions of the type mA->nA (m and n are
positive integers) which possess an absorbing state. The method yields accurate
results for the extinction statistics, and for the quasi-stationary probability
distribution, including large deviations, of the metastable state. The power of
the method is demonstrated on the example of binary annihilation and triple
branching 2A->0 and A->3A, representative of the rather general class of
dissociation-recombination reactions.


  In this paper, we investigate the microscopic dynamics of a polymer of length
$N$ translocating through a narrow pore. Characterization of its purportedly
anomalous dynamics has so far remained incomplete. We show that the polymer
dynamics is anomalous until the Rouse time $\tau_{R}\sim N^{1+2\nu}$, with a
mean square displacement through the pore consistent with
$t^{(1+\nu)/(1+2\nu)}$, with $\nu\approx0.588$ the Flory exponent. This is
shown to be directly related to a decay in time of the excess monomer density
near the pore as $t^{-(1+\nu)/(1+2\nu)}\exp(-t/\tau_{R})$. Beyond the Rouse
time translocation becomes diffusive. In consequence of this, the dwell-time
$\tau_{d}$, the time a translocating polymer typically spends within the pore,
scales as $N^{2+\nu}$, in contrast to previous claims.


  The mechanical unfolding of proteins is investigated by extending the
Wako-Saito-Munoz-Eaton model, a simplified protein model with binary degrees of
freedom, which has proved successful in describing the kinetics of protein
folding. Such a model is generalized by including the effect of an external
force, and its thermodynamics turns out to be exactly solvable. We consider two
molecules, the 27th immunoglobulin domain of titin and protein PIN1. In the
case of titin we determine equilibrium force-extension curves and study
nonequilibrium phenomena in the frameworks of dynamic loading and force clamp
protocols, verifying theoretical laws and finding the position of the kinetic
barrier which hinders the unfolding of the molecule. The PIN1 molecule is used
to check the possibility of computing the free energy landscape as a function
of the molecule length by means of an extended form of the Jarzynski equality.


  Since their introduction, Boolean networks have been traditionally studied in
view of their rich dynamical behavior under different update protocols and for
their qualitative analogy with cell regulatory networks. More recently, tools
borrowed from statistical physics of disordered systems and from computer
science have provided a more complete characterization of their equilibrium
behavior. However, the largest part of the results have been obtained in the
thermodynamic limit, which is often far from being reached when dealing with
realistic instances of the problem. The numerical analysis presented here aims
at comparing - for a specific family of models - the outcomes given by the
heuristic belief propagation algorithm with those given by exhaustive
enumeration. In the second part of the paper some analytical considerations on
the validity of the annealed approximation are discussed.


  We introduce a model of negotiation dynamics whose aim is that of mimicking
the mechanisms leading to opinion and convention formation in a population of
individuals. The negotiation process, as opposed to ``herding-like'' or
``bounded confidence'' driven processes, is based on a microscopic dynamics
where memory and feedback play a central role. Our model displays a
non-equilibrium phase transition from an absorbing state in which all agents
reach a consensus to an active stationary state characterized either by
polarization or fragmentation in clusters of agents with different opinions. We
show the exystence of at least two different universality classes, one for the
case with two possible opinions and one for the case with an unlimited number
of opinions. The phase transition is studied analytically and numerically for
various topologies of the agents' interaction network. In both cases the
universality classes do not seem to depend on the specific interaction
topology, the only relevant feature being the total number of different
opinions ever present in the system.


  We apply the spectral method, recently developed by the authors, to calculate
the statistics of a reaction-limited multi-step birth-death process, or
chemical reaction, that includes as elementary steps branching A->2A and
annihilation 2A->0. The spectral method employs the generating function
technique in conjunction with the Sturm-Liouville theory of linear differential
operators. We focus on the limit when the branching rate is much higher than
the annihilation rate, and obtain accurate analytical results for the complete
probability distribution (including large deviations) of the metastable
long-lived state, and for the extinction time statistics. The analytical
results are in very good agreement with numerical calculations. Furthermore, we
use this example to settle the issue of the "lacking" boundary condition in the
spectral formulation.


  We study a six-species Lotka-Volterra type system on different
two-dimensional lattices when each species has two superior and two inferior
partners. The invasion rates from predator sites to a randomly chosen
neighboring prey's site depend on the predator-prey pair, whereby cyclic
symmetries within the two three-species defensive alliances are conserved.
Monte Carlo simulations reveal an unexpected non-monotonous dependence of
alliance survival on the difference of alliance-specific invasion rates. This
behavior is qualitatively reproduced by a four-point mean-field approximation.
The study addresses fundamental problems of stability for the competition of
two defensive alliances and thus has important implications in natural and
social sciences.


  Theoretical models and rate equations relevant to the Soai reaction are
reviewed. It is found that in a production of chiral molecules from an achiral
substrate autocatalytic processes can induce either enantiomeric excess (ee)
amplification or chiral symmetry breaking. Former terminology means that the
final ee value is larger than the initial value but depends on this, whereas
the latter means the selection of a unique value of the final ee, independent
of the initial value. The ee amplification takes place in an irreversible
reaction such that all the substrate molecules are converted to chiral products
and the reaction comes to a halt. The chiral symmetry breaking is possible when
recycling processes are incorporated. Reactions become reversible and the
system relaxes slowly to a unique final state. The difference between the two
behaviors is apparent in the flow diagram in the phase space of chiral molecule
concentrations. The ee amplification takes place when the flow terminates on a
line of fixed points (or a fixed line), whereas symmetry breaking corresponds
to the dissolution of the fixed line accompanied by the appearance of fixed
points. Relevance of the Soai reaction to the homochirality in life is also
discussed.


  Stochastic aspects of chemical reaction models related to the Soai reactions
as well as to the homochirality in life are studied analytically and
numerically by the use of the master equation and random walk model. For
systems with a recycling process, a unique final probability distribution is
obtained by means of detailed balance conditions. With a nonlinear
autocatalysis the distribution has a double-peak structure, indicating the
chiral symmetry breaking. This problem is further analyzed by examining
eigenvalues and eigenfunctions of the master equation. In the case without
recycling process, final probability distributions depend on the initial
conditions. In the nonlinear autocatalytic case, time-evolution starting from a
complete achiral state leads to a final distribution which differs from that
deduced from the nonzero recycling result. This is due to the absence of the
detailed balance, and a directed random walk model is shown to give the correct
final profile. When the nonlinear autocatalysis is sufficiently strong and the
initial state is achiral, the final probability distribution has a double-peak
structure, related to the enantiomeric excess amplification. It is argued that
with autocatalyses and a very small but nonzero spontaneous production, a
single mother scenario could be a main mechanism to produce the homochirality.


  Cells in multicellular organisms adhere to the extracellular matrix through
two-dimensional clusters spanning a size range from very few to thousands of
adhesion bonds. For many common receptor-ligand systems, the ligands are
tethered to a surface via polymeric spacers with finite binding range, thus
adhesion cluster stability crucially depends on receptor-ligand distance. We
introduce a one-step master equation which incorporates the effect of
cooperative binding through a finite number of polymeric ligand tethers. We
also derive Fokker-Planck and mean field equations as continuum limits of the
master equation. Polymers are modeled either as harmonic springs or as
worm-like chains. In both cases, we find bistability between bound and unbound
states for intermediate values of receptor-ligand distance and calculate the
corresponding switching times. For small cluster sizes, stochastic effects
destabilize the clusters at large separation, as shown by a detailed analysis
of the stochastic potential resulting from the Fokker-Planck equation.


  Nonequilibrium collective motion is ubiquitous in nature and often results in
a rich collection of intringuing phenomena, such as the formation of shocks or
patterns, subdiffusive kinetics, traffic jams, and nonequilibrium phase
transitions. These stochastic many-body features characterize transport
processes in biology, soft condensed matter and, possibly, also in nanoscience.
Inspired by these applications, a wide class of lattice-gas models has recently
been considered. Building on the celebrated {\it totally asymmetric simple
exclusion process} (TASEP) and a generalization accounting for the exchanges
with a reservoir, we discuss the qualitative and quantitative nonequilibrium
properties of these model systems. We specifically analyze the case of a
dimeric lattice gas, the transport in the presence of pointwise disorder and
along coupled tracks.


  We study the extremal dynamics emerging in an out-of-equilibrium
one-dimensional Jepsen gas of $(N+1)$ hard-point particles. The particles
undergo binary elastic collisions, but move ballistically in-between
collisions. The gas is initally uniformly distributed in a box $[-L,0]$ with
the "leader" (or the rightmost particle) at X=0, and a random positive
velocity, independently drawn from a distribution $\phi(V)$, is assigned to
each particle. The gas expands freely at subsequent times. We compute
analytically the distribution of the leader's velocity at time $t$, and also
the mean and the variance of the number of collisions that are undergone by the
leader up to time $t$. We show that in the thermodynamic limit and at fixed
time $t\gg 1$ (the so-called "growing regime"), when interactions are strongly
manifest, the velocity distribution exhibits universal scaling behavior of only
three possible varieties, depending on the tail of $\phi(V)$. The associated
scaling functions are novel and different from the usual extreme-value
distributions of uncorrelated random variables. In this growing regime the mean
and the variance of the number of collisions of the leader up to time $t$
increase logarithmically with $t$, with universal prefactors that are computed
exactly. The implications of our results in the context of biological evolution
modeling are pointed out.


  We systematically study and compare damage spreading at the sparse
percolation (SP) limit for random boolean and threshold networks with
perturbations that are independent of the network size $N$. This limit is
relevant to information and damage propagation in many technological and
natural networks. Using finite size scaling, we identify a new characteristic
connectivity $K_s$, at which the average number of damaged nodes $\bar d$,
after a large number of dynamical updates, is independent of $N$. Based on
marginal damage spreading, we determine the critical connectivity
$K_c^{sparse}(N)$ for finite $N$ at the SP limit and show that it
systematically deviates from $K_c$, established by the annealed approximation,
even for large system sizes. Our findings can potentially explain the results
recently obtained for gene regulatory networks and have important implications
for the evolution of dynamical networks that solve specific computational or
functional tasks.


  It now becomes apparent that condensed DNA toroid which emerges in a poor
solvent condition can be realised in the framework of the non-linear sigma
model on a line segment. In fact, the classical solutions of the model, i.e.,
of the bending potential exhibit toroidal forms, and fit well with the
delta-function term of the attractive interaction. The proposed theory is in
good agreement with experimental observations that the toroid is the ground
state. In this paper, we give a rigorous proof that the solutions are indeed
the exact solutions of the equations of motion with the first derivatives of
the attractive interaction term. We also show a refined mapping to experimental
data, considering the finite size effects of the cross section and of the chain
length.


  We investigate how different microscopic interactions between semiflexible
chain segments can qualitatively alter the physical properties of the condensed
toroid. We propose a general form of the Hamiltonian of the toroid and discuss
its analytic properties. For different interactions, the theory predicts
different scaling behaviours of the mean toroidal and cross sectional radii,
$r_c$ and $r_{cross}$, as functions of the contour length L: $(r_c, r_{cross})
\sim L^{\nu(N_c)}$ with $\nu=(1/5, 2/5)$ for the van der Waals type,
$\nu=(-1/3, 2/3)$ for the Coulomb type, $\nu=(-1, 1)$ for the delta function
type attractions in the asymptotic limit. For the toroids with finite winding
number $N_c=100 \sim 400$, we find $\nu \simeq 0$ for the Yukawa interaction
with screening parameter $\kappa=0.5 \sim 1.0$, and $\nu=0.1 \sim 0.13$ for the
van der Waals type interactions. These findings could provide possible
explanation for the experimentally well known observation $\nu \simeq 0$ of the
condensed DNA toroids. Conformational transitions are also discussed.


  We have studied the competition between helix formation and aggregation for a
simple polymer model. We present simulation results for a system of two such
polymers, examining the potential of mean force, the balance between inter and
intramolecular interactions, and the promotion or disruption of secondary
structure brought on by the proximity of the two molecules. In particular, we
demonstrate that proximity between two such molecules can stabilize secondary
structure. However, for this model, observed secondary structure is not stable
enough to prevent collapse of the system into an unstructured globule.


  We offer a theoretical design of new systems that show promise for digital
biochemical computing, including realizations of error correction by utilizing
redundancy, as well as signal rectification. The approach includes information
processing using encoded DNA sequences, DNAzyme biocatalyzed reactions and the
use of DNA-functionalized magnetic nanoparticles. Digital XOR and NAND logic
gates and copying (fanout) are designed using the same components.


  We study the elastic properties and mechanical stability of viral capsids
under external force-loading with computer simulations. Our approach allows the
implementation of specific geometries corresponding to specific phages such as
$\phi$29 and CCMV. We demonstrate how in a combined numerical and experimental
approach the elastic parameters can be determined with high precision. The
experimentally observed bimodality of elastic spring constants is shown to be
of geometrical origin, namely the presence of pentavalent units in the viral
shell. A criterion for capsid breakage is defined, which explains well the
experimentally observed rupture. From our numerics we find for the dependence
of the rupture force on the F\"oppl-von K\'arm\'an (FvK) number a crossover
from $\gamma^{2/3}$ to $\gamma^{1/2}$. For filled capsids high internal
pressures lead to a stronger destabilization of viruses with a buckled ground
state than unbuckled ones. Finally, we show how our numerically calculated
energy maps can be used to extract information about the strength of
protein-protein interactions from rupture experiments.


  A novel combination of discontinuous molecular dynamics and the Langevin
equation, together with an intermediate-resolution model, are used to carry out
long (several $\mu$s) simulation and study folding transition and transport of
proteins in slit nanopores. Both attractive ($U^+$) and repulsive ($U^-$)
interaction potentials between the proteins and the pore walls are considered.
Near the folding temperature $T_f$ and in the presence of $U^+$ the proteins
undergo a repeating sequence of folding/partially-folding/ unfolding
transitions, while $T_f$ decreases with decreasing pore sizes. The opposite is
true when $U^-$ is present. The proteins' effective diffusivity $D$ is computed
as a function of their length (number of the amino acid groups), temperature
$T$, the pore size, and the interaction potentials $U^\pm$. Far from $T_f$, $D$
increases (roughly) linearly with $T$, but due to the thermal fluctuations and
their effect on the proteins' structure near $T_f$, the dependence of $D$ on
$T$ in this region is nonlinear. Under certain conditions, transport of
proteins in smaller pores can be {\it faster} than that in larger pores.


  We consider a polymer of length $N$ translocating through a narrow pore in
the absence of external fields. Characterization of its purportedly anomalous
dynamics has so far remained incomplete. We show that the polymer dynamics is
anomalous until the Rouse time $\tau_{R}\sim N^{1+2\nu}$, with a mean square
displacement through the pore consistent with $t^{(1+\nu)/(1+2\nu)}$, with
$\nu\approx0.588$ the Flory exponent. This is shown to be directly related to a
decay in time of the excess monomer density near the pore as
$t^{-(1+\nu)/(1+2\nu)}\exp(-t/\tau_{R})$. Beyond the Rouse time translocation
becomes diffusive. In consequence of this, the dwell-time $\tau_{d}$, the time
a translocating polymer typically spends within the pore, scales as
$N^{2+\nu}$, in contrast to previous claims.


  We give a overview of stochastic models of evolution that have found
applications in genetics, ecology and linguistics for an audience of
nonspecialists, especially statistical physicists. In particular, we focus
mostly on neutral models in which no intrinsic advantage is ascribed to a
particular type of the variable unit, for example a gene, appearing in the
theory. In many cases these models are exactly solvable and furthermore go some
way to describing observed features of genetic, ecological and linguistic
systems.


  Motor proteins are active enzyme molecules that play a crucial role in many
biological processes. They transform the chemical energy into the mechanical
work and move unidirectionally along rigid cytoskeleton filaments.
Single-molecule experiments suggest that motor proteins, consisting of two
motor domains, move in a hand-over-hand mechanism when each subunit changes
between trailing and leading positions in alternating steps, and these subunits
do not interact with each other. However, recent experiments on heterodimeric
kinesins suggest that the motion of motor domains is not independent, but
rather strongly coupled and coordinated, although the mechanism of these
interactions are not known. We propose a simple discrete stochastic model to
describe the dynamics of homodimeric and heterodimeric two-headed motor
proteins. It is argued that interactions between motor domains modify free
energy landscapes of each motor subunit, and motor proteins still move via the
hand-over-hand mechanism but with different transitions rates. Our calculations
of biophysical properties agree with experimental observations. Several ways to
test the theoretical model are proposed.


  We solve the problem of discrete translocation of a polymer through a pore,
driven by the irreversible, random sequential adsorption of particles on one
side of the pore. Although the kinetics of the wall motion and the deposition
are coupled, we find the exact steady-state distribution for the gap between
the wall and the nearest deposited particle. This result enables us to
construct the mean translocation velocity demonstrating that translocation is
faster when the adsorbing particles are smaller. Monte-Carlo simulations also
show that smaller particles gives less dispersion in the ratcheted motion. We
also define and compare the relative efficiencies of ratcheting by deposition
of particles with different sizes and we describe an associated
"zone-refinement" process.


  Proteins contain a large fraction of regular, repeating conformations, called
secondary structure. A simple, generic definition of secondary structure is
presented which consists of measuring local correlations along the protein
chain. Using this definition and a simple model for proteins, the forces
driving the formation of secondary structure are explored. The relative role of
energy and entropy are examined. Recent work has indicated that compaction is
sufficient to create secondary structure. We test this hypothesis, using simple
non-lattice protein models.


  Using a simple three-dimensional lattice copolymer model and Monte Carlo
dynamics, we study the collapse and folding of protein-like heteropolymers. The
polymers are 27 monomers long and consist of two monomer types. Although these
chains are too long for exhaustive enumeration of all conformations, it is
possible to enumerate all the maximally compact conformations, which are 3x3x3
cubes. This allows us to select sequences that have a unique global minimum. We
then explore the kinetics of collapse and folding and examine what features
determine the various rates. The folding time has a plateau over a broad range
of temperatures and diverges at both high and low temperatures. The folding
time depends on sequence and is related to the amount of energetic frustration
in the native state. The collapse times of the chains are sequence independent
and are a few orders of magnitude faster than the folding times, indicating a
two-phase folding process. Below a certain temperature the chains exhibit
glass-like behavior, characterized by a slowing down of time scales and loss of
self-averaging behavior. We explicitly define the glass transition temperature
(Tg), and by comparing it to the folding temperature (Tf), we find two classes
of sequences: good folders with Tf > Tg and non-folders with Tf < Tg.


  We consider the Derrida-Higgs (DH) statistical model of species formation in
the case where the population is geographically distributed in discrete
locations, and mating only takes place within one location. Keeping the rate of
migration between neighbouring locations at a fixed value, we change the
mutation rate, changing therefore the average overlap between genotypes. When
the overlap between individuals living in different locations falls below a
fecundity threshold, speciation occurs. When more species coexist, the genetic
structure of the population (as described by the overlap distribution $P(q)$)
fluctuates. However, the average overlap, both within one location and among
neighbouring locations, appears to vary according to the same laws as in the
absence of speciation. The model provides a reasonable estimate of the
parameter values necessary to observe geographic speciation, which is found to
be much more likely than the sympatric speciation of the original DH model.
Applications to the case of circular invasion, where the concept of biological
species appears to run into difficulties, are sketched.


  We will review some of the theoretical progresses that have been in the study
of complex systems in physics and of their applications to biology.


  Contents:
  1. Introduction
  2. Amphiphilic molecules and the phases they form
  3. Isolated membranes: the Helfrich hamiltonian
  4. Vesicle shapes
  5. Shape fluctuations in vesicles
  6. Interacting fluid membranes
  7. Conclusions
  A. Differential equations for vesicle shapes
  B. The Faddeev-Popov determinant
  C. One-loop calculation of the renormalization group
  D. The Liouville model


  The concept of fitness is introduced, and a simple derivation of the
Fundamental Theorem of Natural Selection (which states that the average fitness
of a population increases if its variance is nonzero) is given. After a short
discussion of the adaptative walk model, a short review is given of the
quasispecies approach to molecular evolution and to the error threshold. The
relevance of flat fitness landscapes to molecular evolution is stressed.
Finally a few examples which involve wider concepts of fitness, and in
particular two-level selection, are shortly reviewed.


  In a statistical approach to protein structure analysis, Miyazawa and
Jernigan (MJ) derived a $20\times 20$ matrix of inter-residue contact energies
between different types of amino acids. Using the method of eigenvalue
decomposition, we find that the MJ matrix can be accurately reconstructed from
its first two principal component vectors as $M_{ij}=C_0+C_1(q_i+q_j)+C_2 q_i
q_j$, with constant $C$'s, and 20 $q$ values associated with the 20 amino
acids. This regularity is due to hydrophobic interactions and a force of
demixing, the latter obeying Hildebrand's solubility theory of simple liquids.


  A simple growth model for DNA evolution is introduced which is analytically
solvable and reproduces the observed statistical behavior of real sequences.


  The quantitative description of model protein folding kinetics using a
diffusive collective reaction coordinate is examined. Direct folding kinetics,
diffusional coefficients and free energy profiles are determined from Monte
Carlo simulations of a 27-mer, 3 letter code lattice model, which corresponds
roughly to a small helical protein. Analytic folding calculations, using simple
diffusive rate theory, agree extremely well with the full simulation results.
Folding in this system is best seen as a diffusive, funnel-like process.


  Protein structures in nature often exhibit a high degree of regularity
(secondary structures, tertiary symmetries, etc.) absent in random compact
conformations. We demonstrate in a simple lattice model of protein folding that
structural regularities are related to high designability and evolutionary
stability. We measure the designability of each compact structure by the number
of sequences which can design the structure, i.e., which possess the structure
as their nondegenerate ground state. We find that compact structures are
drastically different in terms of their designability; highly designable
structures emerge with a number of associated sequences much larger than the
average. These structures are found to have ``protein like'' secondary
structure and even tertiary symmetries. In addition, they are also
thermodynamically more stable than ordinary structures. These results suggest
that protein structures are selected because they are easy to design and stable
against mutations, and that such a selection simutaneously leads to
thermodynamic stability.


  We study the evolution of a system of $N$ interacting species which mimics
the dynamics of a cyclic food chain. On a one-dimensional lattice with N<5
species, spatial inhomogeneities develop spontaneously in initially homogeneous
systems. The arising spatial patterns form a mosaic of single-species domains
with algebraically growing size, $\ell(t)\sim t^\alpha$, where $\alpha=3/4$
(1/2) and 1/3 for N=3 with sequential (parallel) dynamics and N=4,
respectively. The domain distribution also exhibits a self-similar spatial
structure which is characterized by an additional length scale, ${\cal
L}(t)\sim t^\beta$, with $\beta=1$ and 2/3 for N=3 and 4, respectively. For
$N\geq 5$, the system quickly reaches a frozen state with non interacting
neighboring species. We investigate the time distribution of the number of
mutations of a site using scaling arguments as well as an exact solution for
N=3. Some possible extensions of the system are analyzed.


  We consider the quasispecies description of a population evolving in both the
"master sequence" landscape (where a single sequence is evolutionarily
preferred over all others) and the REM landscape (where the fitness of
different sequences is an independent, identically distributed, random
variable). We show that, in both cases, the error threshold is analogous to a
first order thermodynamical transition, where the overlap between the average
genotype and the optimal one drops discontinuously to zero.


  We study a homopolymer model of a protein chain, where each monomer carries a
dipole moment. To mimic the geometry of the peptidic bond, these dipoles are
constrained to be locally perpendicular to the chain. The tensorial character
of the dipolar interaction leads naturally to a (tensorial) liquid crystal-like
order parameter. For non chiral chains, a mean field study of this model shows
that a classical $\theta$ collapse transition occurs first; at lower
temperature, nematic order sets in. For chiral chains, an anisotropic
(tensorial) collapse transition may occur before the $\theta$ temperature is
reached: the ordered phase can be described as a ``compact phase of secondary
structures'', and possesses great similarities with the liquid crystal blue
phases.


  We study a simple heteropolymer model containing sequence-independent local
interactions on both square and triangular lattices. Sticking to a two-letter
code, we investigate the model for varying strength $\kappa$ of the local
interactions; $\kappa=0$ corresponds to the well-known HP model [K.F. Lau and
K.A. Dill, Macromolecules 22, 3986 (1989)]. By exhaustive enumerations for
short chains, we obtain all structures which act as a unique and pronounced
energy minimum for at least one sequence. We find that the number of such
designable structures depends strongly on $\kappa$. Also, we find that the
number of designable structures can differ widely for the two lattices at a
given $\kappa$. This is the case, for example, at $\kappa=0$, which implies
that the HP model exhibits different behavior on the two lattices. Our findings
clearly show that sequence-independent local properties of the chains can play
an important role in the formation of unique minimum energy structures.


  DNA and other biopolymers differ from classical polymers due to their
torsional stiffness. This property changes the statistical character of their
conformations under tension from a classical random walk to a problem we call
the `torsional directed walk'. Motivated by a recent experiment on single
lambda-DNA molecules [Strick et al., Science 271 (1996) 1835], we formulate the
torsional directed walk problem and solve it analytically in the appropriate
force regime. Our technique affords a direct physical determination of the
microscopic twist stiffness C and twist-stretch coupling D relevant for DNA
functionality. The theory quantitatively fits existing experimental data for
relative extension as a function of overtwist over a wide range of applied
force; fitting to the experimental data yields the numerical values C=120nm and
D=50nm. Future experiments will refine these values. We also predict that the
phenomenon of reduction of effective twist stiffness by bend fluctuations
should be testable in future single-molecule experiments, and we give its
analytic form.