  In unicellular organisms such as bacteria the same acquired mutations
beneficial in one environment can be restrictive in another. However, evolving
Escherichia coli populations demonstrate remarkable flexibility in adaptation.
The mechanisms sustaining genetic flexibility remain unclear. In E. coli the
transcriptional regulation of gene expression involves both dedicated
regulators binding specific DNA sites with high affinity and also global
regulators - abundant DNA architectural proteins of the bacterial chromoid
binding multiple low affinity sites and thus modulating the superhelical
density of DNA. The first form of transcriptional regulation is dominantly
pairwise and specific, representing digitial control, while the second form is
(in strength and distribution) continuous, representing analog control. Here we
look at the properties of effective networks derived from significant gene
expression changes under variation of the two forms of control and find that
upon limitations of one type of control (caused e.g. by mutation of a global
DNA architectural factor) the other type can compensate for compromised
regulation. Mutations of global regulators significantly enhance the digital
control; in the presence of global DNA architectural proteins regulation is
mostly of the analog type, coupling spatially neighboring genomic loci;
together our data suggest that two logically distinct types of control are
balancing each other. By revealing two distinct logical types of control, our
approach provides basic insights into both the organizational principles of
transcriptional regulation and the mechanisms buffering genetic flexibility. We
anticipate that the general concept of distinguishing logical types of control
will apply to many complex biological networks.


  We study the dynamics of parallel tempering simulations, also known as the
replica exchange technique, which has become the method of choice for
simulation of proteins and other complex systems. Recent results for the
optimal choice of the control parameter discretization allow a treatment
independent of the system in question. Analyzing mean first passage times
across control parameter space, we find an expression for the optimal number of
replicas in simulations covering a given temperature range. Our results suggest
a particular protocol to optimize the number of replicas in actual simulations.


This paper has been withdrawn.


  Small single-domain proteins often exhibit only a single free-energy barrier,
or transition state, between the denatured and the native state. The folding
kinetics of these proteins is usually explored via mutational analysis. A
central question is which structural information on the transition state can be
derived from the mutational data. In this article, we model and structurally
interpret mutational Phi-values for two small beta-sheet proteins, the PIN and
the FBP WW domain. The native structure of these WW domains comprises two
beta-hairpins that form a three-stranded beta-sheet. In our model, we assume
that the transition state consists of two conformations in which either one of
the hairpins is formed. Such a transition state has been recently observed in
Molecular Dynamics folding-unfolding simulations of a small designed
three-stranded beta-sheet protein. We obtain good agreement with the
experimental data (i) by splitting up the mutation-induced free-energy changes
into terms for the two hairpins and for the small hydrophobic core of the
proteins, and (ii) by fitting a single parameter, the relative degree to which
hairpin 1 and 2 are formed in the transition state. The model helps to
understand how mutations affect the folding kinetics of WW domains, and
captures also negative Phi-values that have been difficult to interpret.


  Simple theoretical concepts and models have been helpful to understand the
folding rates and routes of single-domain proteins. As reviewed in this
article, a physical principle that appears to underly these models is loop
closure.


  A model of autocrine signaling in cultures of suspended cells is developed on
the basis of the effective medium approximation. The fraction of autocrine
ligands, the mean and distribution of distances traveled by paracrine ligands
before binding, as well as the mean and distribution of the ligand lifetime are
derived. Interferon signaling by dendritic immune cells is considered as an
illustration.


  We find that discrete noise of inhibiting (signal) molecules can greatly
delay the extinction of plasmids in a plasmid replication system: a
prototypical biochemical regulatory network. We calculate the probability
distribution of the metastable state of the plasmids and show on this example
that the reaction rate equations may fail in predicting the average number of
regulated molecules even when this number is large, and the time is much
shorter than the mean extinction time.


  In many intracellular processes, the length distribution of microtubules is
controlled by depolymerizing motor proteins. Experiments have shown that,
following non-specific binding to the surface of a microtubule, depolymerizers
are transported to the microtubule tip(s) by diffusion or directed walk and,
then, depolymerize the microtubule from the tip(s) after accumulating there. We
develop a quantitative model to study the depolymerizing action of such a
generic motor protein, and its possible effects on the length distribution of
microtubules. We show that, when the motor protein concentration in solution
exceeds a critical value, a steady state is reached where the length
distribution is, in general, non-monotonic with a single peak. However, for
highly processive motors and large motor densities, this distribution
effectively becomes an exponential decay. Our findings suggest that such motor
proteins may be selectively used by the cell to ensure precise control of MT
lengths. The model is also used to analyze experimental observations of
motor-induced depolymerization.


  We review theoretical results for the adhesion-induced phase behavior of
biomembranes. The focus is on models in which the membranes are represented as
discretized elastic sheets with embedded adhesion molecules. We present several
mechanism that lead to the formation of domains during adhesion, and discuss
the time-dependent evolution of domain patterns obtained in Monte-Carlo
simulations. The simulated pattern dynamics has striking similarities to the
pattern evolution observed during T cell adhesion.


  The diversity of a community that cannot be fully counted must be inferred.
The two preeminent inference methods are the MaxEnt method, which uses
information in the form of constraints and Bayes' rule which uses information
in the form of data. It has been shown that these two methods are special cases
of the method of Maximum (relative) Entropy (ME). We demonstrate how this
method can be used as a measure of diversity that not only reproduces the
features of Shannon's index but exceeds them by allowing more types of
information to be included in the inference. A specific example is solved in
detail. Additionally, the entropy that is found is the same form as the
thermodynamic entropy.


  The isotopic composition, for example, 14C/12C, 13C/12C, 2H/1H, 15N/14N and
18O/16O, of the elements of matter is heterogeneous. It is ruled by physical,
chemical and biological mechanisms. Isotopes can be employed to follow the fate
of mineral and organic compounds during biogeochemical transformations. The
determination of the isotopic composition of organic substances occurring at
trace level in very complex mixtures such as sediments, soils and blood, has
been made possible during the last 20 years due to the rapid development of
molecular level isotopic techniques. After a brief glance at pioneering studies
revealing isotopic breakthroughs at the molecular and intramolecular levels,
this paper reviews selected applications of compound-specific isotope analysis
in various scientific fields.


  In addition to variation in terms of single nucleotide polymorphisms (SNPs),
whole regions ranging from several kilobases up to a megabase in length differ
in copy number among individuals. These differences are referred to as Copy
Number Variants (CNVs) and extensive mapping of these is underway. Recent
studies have highlighted their great prevalence in the human genome. Segmental
Duplications (SDs) are long (>1kb) stretches of duplicated DNA with high
sequence identity. First, we analyzed the co-localization of SDs and find that
SDs are significantly co-localized with each other, resulting in a power-law
distribution, which suggests a preferential attachment mechanism, i.e. existing
SDs are likely to be involved in creating new ones nearby. Second, we look at
the relationship of CNVs/SDs with various types of repeats. We we find that the
previously recognized association of SDs with Alu elements is significantly
stronger for older SDs and is sharply decreasing for younger ones. While it
might be expected that the patterns should be similar for SDs and CNVs, we
find, surprisingly, no association of CNVs with Alu elements. This trend is
consistent with the decreasing correlation between Alu elements and younger
SDs, the activity of Alu elements has been decreasing and by now it they seem
no longer active. Furthermore, we find a striking association of SDs with
processed pseudogenes suggesting that they may also have mediated SD formation.
Moreover, find strong association with microsatellites for both SDs and CNVs
that suggests a role for satellites in the formation of both.


  The identification of genes essential for survival is important for the
understanding of the minimal requirements for cellular life and for drug
design. As experimental studies with the purpose of building a catalog of
essential genes for a given organism are time-consuming and laborious, a
computational approach which could predict gene essentiality with high accuracy
would be of great value. We present here a novel computational approach, called
NTPGE (Network Topology-based Prediction of Gene Essentiality), that relies on
network topology features of a gene to estimate its essentiality. The first
step of NTPGE is to construct the integrated molecular network for a given
organism comprising protein physical, metabolic and transcriptional regulation
interactions. The second step consists in training a decision tree-based
machine learning algorithm on known essential and non-essential genes of the
organism of interest, considering as learning attributes the network topology
information for each of these genes. Finally, the decision tree classifier
generated is applied to the set of genes of this organism to estimate
essentiality for each gene. We applied the NTPGE approach for discovering
essential genes in Escherichia coli and then assessed its performance.


  Changes in a cell's external or internal conditions are usually reflected in
the concentrations of the relevant transcription factors. These proteins in
turn modulate the expression levels of the genes under their control and
sometimes need to perform non-trivial computations that integrate several
inputs and affect multiple genes. At the same time, the activities of the
regulated genes would fluctuate even if the inputs were held fixed, as a
consequence of the intrinsic noise in the system, and such noise must
fundamentally limit the reliability of any genetic computation. Here we use
information theory to formalize the notion of information transmission in
simple genetic regulatory elements in the presence of physically realistic
noise sources. The dependence of this "channel capacity" on noise parameters,
cooperativity and cost of making signaling molecules is explored
systematically. We find that, at least in principle, capacities higher than one
bit should be achievable and that consequently genetic regulation is not
limited the use of binary, or "on-off", components.


  We examine an evolutionary naming-game model where communicating agents are
equipped with an evolutionarily selected learning ability. Such a coupling of
biological and linguistic ingredients results in an abrupt transition: upon a
small change of a model control parameter a poorly communicating group of
linguistically unskilled agents transforms into almost perfectly communicating
group with large learning abilities. When learning ability is kept fixed, the
transition appears to be continuous. Genetic imprinting of the learning
abilities proceeds via Baldwin effect: initially unskilled communicating agents
learn a language and that creates a niche in which there is an evolutionary
pressure for the increase of learning ability.Our model suggests that when
linguistic (or cultural) processes became intensive enough, a transition took
place where both linguistic performance and biological endowment of our species
experienced an abrupt change that perhaps triggered the rapid expansion of
human civilization.


  We introduce a simple algorithm for reconstructing phylogenies from multiple
gene trees in the presence of incomplete lineage sorting, that is, when the
topology of the gene trees may differ from that of the species tree. We show
that our technique is statistically consistent under standard stochastic
assumptions, that is, it returns the correct tree given sufficiently many
unlinked loci. We also show that it can tolerate moderate estimation errors.


  Noise and spatial degrees of freedom characterize most ecosystems. Some
aspects of their influence on the coevolution of populations with cyclic
interspecies competition have been demonstrated in recent experiments [e.g. B.
Kerr et al., Nature {\bf 418}, 171 (2002)]. To reach a better theoretical
understanding of these phenomena, we consider a paradigmatic spatial model
where three species exhibit cyclic dominance. Using an individual-based
description, as well as stochastic partial differential and deterministic
reaction-diffusion equations, we account for stochastic fluctuations and
spatial diffusion at different levels, and show how fascinating patterns of
entangled spirals emerge. We rationalize our analysis by computing the
spatio-temporal correlation functions and provide analytical expressions for
the front velocity and the wavelength of the propagating spiral waves.


  A hypothesis of the evolution of the genetic code is proposed, the leading
mechanism of which is the nucleotide spontaneous damage leading to
AT-enrichment of the genome. The hypothesis accounts for stability of the
genetic code towards point mutations, the presence of code dialects, and the
symmetry of the genetic code table.


  Protein detection on SDS gels or on 2-D gels must combine several features,
such as sensitivity, homogeneity from one protein to another, speed, low cost,
and user-friendliness. For some applications, it is also interesting to have a
nonfixing stain, so that proteins can be mobilized from the gel for further use
(electroelution, blotting). We show here that coelectrophoretic staining by
fluorophores of the oxacarbocyanine family, and especially
diheptyloxacarbocyanine, offers several positive features. The sensitivity is
intermediate between the one of colloidal CBB and the one of fluroescent
ruthenium complexes. Detection is achieved within 1 h after the end of the
electrophoretic process and does not use any fixing or toxic agent. The
fluorescent SDS-carbocyanine-protein complexes can be detected either with a
laser scanner with an excitation wavelength of 488 nm or with a UV table
operating at 302 nm. Excellent sequence coverage in subsequent MS analysis of
proteolytic peptides is also achieved with this detection method.


  Metadynamics is a powerful computational tool to obtain the free energy
landscape of complex systems. The Monte Carlo algorithm has proven useful to
calculate thermodynamic quantities associated with simplified models of
proteins, and thus to gain an ever-increasing understanding on the general
principles underlying the mechanism of protein folding. We show that it is
possible to couple metadynamics and Monte Carlo algorithms to obtain the free
energy of model proteins in a way which is computationally very economical.


  We develop a matrix-based approach to predict and verify indirect
interactions in gene and protein regulatory networks. It is based on the
approximate transitivity of indirect regulations (e.g. A regulates B and B
regulates C often implies that A regulates C) and optimally takes into account
the length of a cascade and signs of intermediate interactions. Our method is
at its most powerful when applied to large and densely interconnected networks.
It successfully predicts both the yet unknown indirect regulations, as well as
the sign (activation or repression) of already known ones. The reliability of
sign predictions was calibrated using the gold-standard sets of positive and
negative interactions. We fine-tuned the parameters of our algorithm by
maximizing the area under the Receiver Operating Characteristic (ROC) curve. We
then applied the optimized algorithm to large literature-derived networks of
all direct and indirect regulatory interactions in several model organisms
(Homo sapiens, Saccharomyces cerevisiae, Arabidopsis thaliana and Drosophila
melanogaster).


  Extracting the intrinsic kinetic information of biological molecule from its
single-molecule kinetic data is of considerable biophysical interest. In this
work, we theoretically investigate the feasibility of inferring single RNA's
intrinsic kinetic parameters from the time series obtained by forced
folding/unfolding experiment done in the light tweezer, where the molecule is
flanked by long double-stranded DNA/RNA handles and tethered between two big
beads. We first construct a coarse-grain physical model of the experimental
system. The model has captured the major physical factors: the Brownian motion
of the bead, the molecular structural transition, and the elasticity of the
handles and RNA. Then based on an analytic solution of the model, a Bayesian
method using Monte Carlo Markov Chain is proposed to infer the intrinsic
kinetic parameters of the RNA from the noisy time series of the distance or
force. Because the force fluctuation induced by the Brownian motion of the bead
and the structural transition can significantly modulate the transition rates
of the RNA, we prove that, this statistic method is more accurate and efficient
than the conventional histogram fitting method in inferring the molecule's
intrinsic parameters.


  In eukaryotic genes the protein coding sequence is split into several
fragments, the exons, separated by non-coding DNA stretches, the introns.
Prokaryotes do not have introns in their genome. We report the calculations of
stability domains of actin genes for various organisms in the animal, plant and
fungi kingdoms. Actin genes have been chosen because they have been highly
conserved during evolution. In these genes all introns were removed so as to
mimic ancient genes at the time of the early eukaryotic development, i.e.
before introns insertion. Common stability boundaries are found in evolutionary
distant organisms, which implies that these boundaries date from the early
origin of eukaryotes. In general boundaries correspond with introns positions
of vertebrates and other animals actins, but not much for plants and fungi. The
sharpest boundary is found in a locus where fungi, algae and animals have
introns in positions separated by one nucleotide only, which identifies a
hot-spot for insertion. These results suggest that some introns may have been
incorporated into the genomes through a thermodynamic driven mechanism, in
agreement with previous observations on human genes. They also suggest a
different mechanism for introns insertion in plants and animals.


  Tethered particle motion is an experimental technique to monitor
conformational changes in single molecules of DNA in real time, by observing
the position fluctuations of a micrometer-size particle attached to the DNA.
This article reviews some recent work on theoretical problems inherent in the
interpretation of TPM experiments, both in equilibrium and dynamical aspects.


  We discovered a dynamic phase transition induced by sexual reproduction. The
dynamics is a pure Darwinian rule with both fundamental ingredients to drive
evolution: 1) random mutations and crossings which act in the sense of
increasing the entropy (or diversity); and 2) selection which acts in the
opposite sense by limiting the entropy explosion. Selection wins this
competition if mutations performed at birth are few enough. By slowly
increasing the average number m of mutations, however, the population suddenly
undergoes a mutational degradation precisely at a transition point mc. Above
this point, the "bad" alleles spread over the genetic pool of the population,
overcoming the selection pressure. Individuals become selectively alike, and
evolution stops. Only below this point, m < mc, evolutionary life is possible.
  The finite-size-scaling behaviour of this transition is exhibited for large
enough "chromosome" lengths L. One important and surprising observation is the
L-independence of the transition curves, for large L. They are also independent
on the population size. Another is that mc is near unity, i.e. life cannot be
stable with much more than one mutation per diploid genome, independent of the
chromosome length, in agreement with reality. One possible consequence is that
an eventual evolutionary jump towards larger L enabling the storage of more
genetic information would demand an improved DNA copying machinery in order to
keep the same total number of mutations per offspring.


  What determines species diversity is dramatic concern in science. Here we
report the effect of doping on diversity in spatiotemporal rock-paper-scissors
(RPS) games, which can be observed directly in ecological, biological and
social systems in nature. Doping means that there exists some buffer patches
which do not involve the main procession of the conflicts but occupied the game
space. Quantitative lattices simulation finds that (1) decrease of extinction
possibility is exponential dependent on the increase of doping rate, (2) the
possibility of the conflict is independent of doping rate at well mix evolution
beginning, and is buffered by doping in long time coexistence procession.
Practical meaning of doping are discussed. To demonstrate the importance of
doping, we present one practical example for microbial laboratory efficient
operation and one theoretical example for human-environment co-existence system
better understanding. It suggests that, for diversity, doping can not be
neglected.


  Intracellular transport along microtubules or actin filaments, powered by
molecular motors such as kinesins, dyneins or myosins, has been recently
modeled using one-dimensional driven lattice gases. We discuss some
generalizations of these models, that include extended particles and defects.
We investigate the feasibility of single molecule experiments aiming to measure
the average motor density and to locate the position of traffic jams by mean of
a tracer particle. Finally, we comment on preliminary single molecule
experiments performed in living cells.


  The genetic repressilator circuit consists of three transcription factors, or
repressors, which negatively regulate each other in a cyclic manner. This
circuit was synthetically constructed on plasmids in {\it Escherichia coli} and
was found to exhibit oscillations in the concentrations of the three
repressors. Since the repressors and their binding sites often appear in low
copy numbers, the oscillations are noisy and irregular. Therefore, the
repressilator circuit cannot be fully analyzed using deterministic methods such
as rate-equations. Here we perform stochastic analysis of the repressilator
circuit using the master equation and Monte Carlo simulations. It is found that
fluctuations modify the range of conditions in which oscillations appear as
well as their amplitude and period, compared to the deterministic equations.
The deterministic and stochastic approaches coincide only in the limit in which
all the relevant components, including free proteins, plasmids and bound
proteins, appear in high copy numbers. We also find that subtle features such
as cooperative binding and bound-repressor degradation strongly affect the
existence and properties of the oscillations.


  Comparative genomics methods are widely used to aid the functional annotation
of non coding DNA regions. However, aligning non coding sequences requires new
algorithms and strategies, in order to take into account extensive
rearrangements and turnover during evolution. Here we present a novel large
scale alignment strategy which aims at drawing a precise map of conserved non
coding regions between genomes, even when these regions have undergone small
scale rearrangments events and a certain degree of sequence variability. We
applied our alignment approach to obtain a genome-wide catalogue of conserved
non coding blocks (CNBs) between Drosophila melanogaster and 11 other
Drosophila species. Interestingly, we observe numerous small scale
rearrangement events, such as local inversions, duplications and
translocations, which are not observable in the whole genome alignments
currently available. The high rate of observed low scale reshuffling show that
this database of CNBs can constitute the starting point for several
investigations, related to the evolution of regulatory DNA in Drosophila and
the in silico identification of unannotated functional elements.


  Cellular networks undergo rearrangements during stress and diseases. In
un-stressed state the yeast protein-protein interaction network (interactome)
is highly compact, and the centrally organized modules have a large overlap.
During stress several original modules became more separated, and a number of
novel modules also appear. A few basic functions, such as the proteasome
preserve their central position. However, several functions with high energy
demand, such the cell-cycle regulation loose their original centrality during
stress. A number of key stress-dependent protein complexes, such as the
disaggregation-specific chaperone, Hsp104, gain centrality in the stressed
yeast interactome. Molecular chaperones, heat shock, or stress proteins form
complex interaction networks (the chaperome) with each other and their
partners. Here we show that the human chaperome recovers the segregation of
protein synthesis-coupled and stress-related chaperones observed in yeast
recently. Examination of yeast and human interactomes shows that (1) chaperones
are inter-modular integrators of protein-protein interaction networks, which
(2) often bridge hubs and (3) are favorite candidates for extensive
phosphorylation. Moreover, chaperones (4) become more central in the
organization of the isolated modules of the stressed yeast protein-protein
interaction network, which highlights their importance in the de-coupling and
re-coupling of network modules during and after stress. Chaperone-mediated
evolvability of cellular networks may play a key role in cellular adaptation
during stress and various polygenic and chronic diseases, such as cancer,
diabetes or neurodegeneration.


  In this work, charge transport (CT) properties of the p53 gene are
numerically studied by the transfer matrix method, and using either single or
double strand effective tight-binding models. A statistical analysis of the
consequences of known p53 point mutations on CT features is performed. It is
found that in contrast to other kind of mutation defects, cancerous mutations
result in much weaker changes of CT efficiency. Given the envisioned role
played by CT in the DNA-repairing mechanism, our theoretical results suggest an
underlying physical explanation at the origin of carcinogenesis.


  Gene regulation is a complex process involving the role of several genomic
elements which work in concert to drive spatio-temporal expression. The
experimental characterization of gene regulatory elements is a very complex and
resource-intensive process. One of the major goals in computational biology is
the \textit{in-silico} annotation of previously uncharacterized elements using
results from the subset of known, previously annotated, regulatory elements.
  The recent results of the ENCODE project (\emph{http://encode.nih.gov})
presented in-depth analysis of such functional (regulatory) non-coding elements
for 1% of the human genome. It is hoped that the results obtained on this
subset can be scaled to the rest of the genome. This is an extremely important
effort which will enable faster dissection of other functional elements in key
biological processes such as disease progression and organ development
(\cite{Kleinjan2005},\cite{Lieb2006}. The computational annotation of these
hitherto uncharacterized regions would require an identification of features
that have good predictive value.
  In this work, we study transcriptional regulation as a problem in
heterogeneous data integration, across sequence, expression and interactome
level attributes. Using the example of the \textit{Gata2} gene and its recently
discovered urogenital enhancers \cite{Khandekar2004} as a case study, we
examine the predictive value of various high throughput functional genomic
assays (from projects like ENCODE and SymAtlas) in characterizing these
enhancers and their regulatory role. Observing results from the application of
modern statistical learning methodologies for each of these data modalities, we
propose a set of features that are most discriminatory to find these enhancers.


  The conformational dynamics of a single protein molecule in a shear flow is
investigated using Brownian dynamics simulations. A structure-based coarse
grained model of a protein is used. We consider two proteins, ubiquitin and
integrin, and find that at moderate shear rates they unfold through a sequence
of metastable states - a pattern which is distinct from a smooth unraveling
found in homopolymers. Full unfolding occurs only at very large shear rates.
Furthermore, the hydrodynamic interactions between the amino acids are shown to
hinder the shear flow unfolding. The characteristics of the unfolding process
depend on whether a protein is anchored or not, and if it is, on the choice of
an anchoring point.


  We demonstrate that a common-line method can assemble a 3D oversampled
diffracted intensity distribution suitable for high-resolution structure
solution from a set of measured 2D diffraction patterns, as proposed in
experiments with an X-ray free electron laser (XFEL) (Neutze {\it et al.},
2000). Even for a flat Ewald sphere, we show how the ambiguities due to
Friedel's Law may be overcome. The method breaks down for photon counts below
about 10 per detector pixel, almost 3 orders of magnitude higher than expected
for scattering by a 500 kDa protein with an XFEL beam focused to a 0.1 micron
diameter spot. Even if 10**3 orientationally similar diffraction patterns could
be identified and added to reach the requisite photon count per pixel, the need
for about 10**6 orientational classes for high-resolution structure
determination suggests that about ~ 10**9 diffraction patterns must be
recorded. Assuming pulse and read-out rates of 100 Hz, such measurements would
require ~ 10**7 seconds, i.e. several months of continuous beam time.


  The evolution of complex molecular traits such as disulphide bridges often
requires multiple mutations. The intermediate steps in such evolutionary
trajectories are likely to be selectively neutral or deleterious. Therefore,
large populations and long times may be required to evolve such traits. We
propose that errors in transcription and translation may allow selection for
the intermediate mutations if the final trait provides a large enough selective
advantage. We test this hypothesis using a population based model of protein
evolution. If an individual acquires one of two mutations needed for a novel
trait, the second mutation can be introduced into the phenotype due to
transcription and translation errors. If the novel trait is advantageous
enough, the allele with only one mutation will spread through the population,
even though the gene sequence does not yet code for the complete trait. The
first mutation then has a higher frequency than expected without phenotypic
mutations giving the second mutation a higher probability of fixation. Thus,
errors allow protein sequences to ''look-ahead'' for a more direct path to a
complex trait.


  This thesis develops and expands upon known techniques of mathematical
physics relevant to the analysis of the popular Markov model of phylogenetic
trees required in biology to reconstruct the evolutionary relationships of
taxonomic units from biomolecular sequence data. The techniques of mathematical
physics are plethora and have been developed for some time. The Markov model of
phylogenetics and its analysis is a relatively new technique where most
progress to date has been achieved by using discrete mathematics. This thesis
takes a group theoretical approach to the problem by beginning with a
remarkable mathematical parallel to the process of scattering in particle
physics. This is shown to equate to branching events in the evolutionary
history of molecular units. The major technical result of this thesis is the
derivation of existence proofs and computational techniques for calculating
polynomial group invariant functions on a multi-linear space where the group
action is that relevant to a Markovian time evolution. The practical results of
this thesis are an extended analysis of the use of invariant functions in
distance based methods and the presentation of a new reconstruction technique
for quartet trees which is consistent with the most general Markov model of
sequence evolution.


  This article concerns the following question arising in computational
evolutionary biology. For a given subclass of phylogenetic networks, what is
the maximum value of 0 <= p <= 1 such that for every input set T of rooted
triplets, there exists some network N(T) from the subclass such that at least
p|T| of the triplets are consistent with N(T)? Here we prove that the set
containing all triplets (the full triplet set) in some sense defines p, and
moreover that any network N achieving fraction p' for the full triplet set can
be converted in polynomial time into an isomorphic network N'(T) achieving >=
p' for an arbitrary triplet set T. We demonstrate the power of this result for
the field of phylogenetics by giving worst-case optimal algorithms for level-1
phylogenetic networks (a much-studied extension of phylogenetic trees),
improving considerably upon the 5/12 fraction obtained recently by Jansson,
Nguyen and Sung. For level-2 phylogenetic networks we show that p >= 0.61. We
note that all the results in this article also apply to weighted triplet sets.


  The coexistence between different types of templates has been the choice
solution to the information crisis of prebiotic evolution, triggered by the
finding that a single RNA-like template cannot carry enough information to code
for any useful replicase. In principle, confining $d$ distinct templates of
length $L$ in a package or protocell, whose survival depends on the coexistence
of the templates it holds in, could resolve this crisis provided that $d$ is
made sufficiently large. Here we review the prototypical package model of
Niesert et al. 1981 which guarantees the greatest possible region of viability
of the protocell population, and show that this model, and hence the entire
package approach, does not resolve the information crisis. This is so because
to secure survival the total information content of the protocell, $Ld$, must
tend to a constant value that depends only on the spontaneous error rate per
nucleotide of the template replication mechanism. As a result, an increase of
$d$ must be followed by a decrease of $L$ to ensure the protocell viability, so
that the net information gain is null.


  We study a spatial cyclic predator-prey model with an even number of species
(for n=4, 6, and 8) that allows the formation of two defective alliances
consisting of the even and odd label species. The species are distributed on
the sites of a square lattice. The evolution of spatial distribution is
governed by iteration of two elementary processes on neighboring sites chosen
randomly: if the sites are occupied by a predator-prey pair then the predator
invades the prey's site; otherwise the species exchange their site with a
probability $X$. For low $X$ values a self-organizing pattern is maintained by
cyclic invasions. If $X$ exceeds a threshold value then two types of domains
grow up that formed by the odd and even label species, respectively. Monte
Carlo simulations indicate the blocking of this segregation process within a
range of X for n=8.


  Two crucial elements facilitate the understanding and control of communicable
disease spread within a social setting. These components are, the underlying
contact structure among individuals that determines the pattern of disease
transmission; and the evolution of this pattern over time. Mathematical models
of infectious diseases, which are in principle analytically tractable, use two
general approaches to incorporate these elements. The first approach, generally
known as compartmental modeling, addresses the time evolution of disease spread
at the expense of simplifying the pattern of transmission. On the other hand,
the second approach uses network theory to incorporate detailed information
pertaining to the underlying contact structure among individuals. However,
while providing accurate estimates on the final size of outbreaks/epidemics,
this approach, in its current formalism, disregards the progression of time
during outbreaks. So far, the only alternative that enables the integration of
both aspects of disease spread simultaneously has been to abandon the
analytical approach and rely on computer simulations. We offer a new analytical
framework based on percolation theory, which incorporates both the complexity
of contact network structure and the time progression of disease spread.
Furthermore, we demonstrate that this framework is equally effective on finite-
and "infinite"-size networks. Application of this formalism is not limited to
disease spread; it can be equally applied to similar percolation phenomena on
networks in other areas in science and technology.


  Mathematical models are increasingly used in both academia and the
pharmaceutical industry to understand how phenotypes emerge from systems of
molecular interactions. However, their current construction as monolithic sets
of equations presents a fundamental barrier to progress. Overcoming this
requires modularity, enabling sub-systems to be specified independently and
combined incrementally, and abstraction, enabling general properties to be
specified independently of specific instances. These in turn require models to
be represented as programs rather than as datatypes. Programmable modularity
and abstraction enables libraries of modules to be created for generic
biological processes, which can be instantiated and re-used repeatedly in
different contexts with different components. We have developed a computational
infrastructure to support this. We show here why these capabilities are needed,
what is required to implement them and what can be accomplished with them that
could not be done previously.


  We investigate the selective forces that promote the emergence of modularity
in nature. We demonstrate the spontaneous emergence of modularity in a
population of individuals that evolve in a changing environment. We show that
the level of modularity correlates with the rapidity and severity of
environmental change. The modularity arises as a synergistic response to the
noise in the environment in the presence of horizontal gene transfer. We
suggest that the hierarchical structure observed in the natural world may be a
broken symmetry state, which generically results from evolution in a changing
environment.


  A branching random walk in presence of an absorbing wall moving at a constant
velocity $v$ undergoes a phase transition as the velocity $v$ of the wall
varies. Below the critical velocity $v_c$, the population has a non-zero
survival probability and when the population survives its size grows
exponentially. We investigate the histories of the population conditioned on
having a single survivor at some final time $T$. We study the quasi-stationary
regime for $v<v_c$ when $T$ is large. To do so, one can construct a modified
stochastic process which is equivalent to the original process conditioned on
having a single survivor at final time $T$. We then use this construction to
show that the properties of the quasi-stationary regime are universal when
$v\to v_c$. We also solve exactly a simple version of the problem, the
exponential model, for which the study of the quasi-stationary regime can be
reduced to the analysis of a single one-dimensional map.


  Branching process approximation to the initial stages of an epidemic process
has been used since the 1950's as a technique for providing stochastic
counterparts to deterministic epidemic threshold theorems. One way of
describing the approximation is to construct both branching and epidemic
processes on the same probability space, in such a way that their paths
coincide for as long as possible. In this paper, it is shown, in the context of
a Markovian model of parasitic infection, that coincidence can be achieved with
asymptotically high probability until o(N^{2/3}) infections have occurred,
where N denotes the total number of hosts.


  The evolution of complexity has been a central theme for Biology [2] and
Artificial Life research [1]. It is generally agreed that complexity has
increased in our universe, giving way to life, multi-cellularity, societies,
and systems of higher complexities. However, the mechanisms behind the
complexification and its relation to evolution are not well understood.
Moreover complexification can be used to mean different things in different
contexts. For example, complexification has been interpreted as a process of
diversification between evolving units [2] or as a scaling process related to
the idea of transitions between different levels of complexity [7].
Understanding the difference or overlap between the mechanisms involved in both
situations is mandatory to create acceptable synthetic models of the process,
as is required in Artificial Life research. (...)


  A description for regulatory genetic network based on generalized potential
energy is constructed. The potential energy is derived from the steady state
solution of linearized Fokker-Plank equation, and the result is shown to be
equivalent to the system of coupled oscillators. The correspondence between the
quantities from the mechanical picture and the steady-state fluctuations is
established. Explicit calculation is given for auto-regulatory networks in
which, the force constant associated with the degree of protein is very weak.
Negative feedback not only suppresses the fluctuations but also increases the
steepness of the potential. The results for the fluctuations agree completely
with those obtained from linear noise Fokker-Planck equation.


  The role of post-translational modification of histones in eukaryotic gene
regulation is well recognized. Epigenetic silencing of genes via heritable
chromatin modifications plays a major role in cell fate specification in higher
organisms. We formulate a coarse-grained model of chromatin silencing in yeast
and study the conditions under which the system becomes bistable, allowing for
different epigenetic states. We also study the dynamics of the boundary between
the two locally stable states of chromatin: silenced and unsilenced. The model
could be of use in guiding the discussion on chromatin silencing in general. In
the context of silencing in budding yeast, it helps us understand the phenotype
of various mutants, some of which may be non-trivial to see without the help of
a mathematical model. One such example is a mutation that reduces the rate of
background acetylation of particular histone side-chains that competes with the
deacetylation by Sir2p. The resulting negative feedback due to a Sir protein
depletion effect gives rise to interesting counter-intuitive consequences. Our
mathematical analysis brings forth the different dynamical behaviors possible
within the same molecular model and guides the formulation of more refined
hypotheses that could be addressed experimentally.


  Background. All-atom crystallographic refinement of proteins is a laborious
manually driven procedure, as a result of which, alternative and multiconformer
interpretations are not routinely investigated.
  Results. We describe efficient loop sampling procedures in Rappertk and
demonstrate that single loops in proteins can be automatically and accurately
modelled with few positional restraints. Loops constructed with a composite
CNS/Rappertk protocol consistently have better Rfree than those with CNS alone.
This approach is extended to a more realistic scenario where there are often
large positional uncertainties in loops along with small imperfections in the
secondary structural framework. Both ensemble and collection methods are used
to estimate the structural heterogeneity of loop regions.
  Conclusion. Apart from benchmarking Rappertk for the all-atom protein
refinement task, this work also demonstrates its utility in both aspects of
loop modelling - building a single conformer and estimating structural
heterogeneity the loops can exhibit.


  Background. Dramatic increases in RNA structural data have made it possible
to recognize its conformational preferences much better than a decade ago. This
has created an opportunity to use discrete restraint-based conformational
sampling for modelling RNA and automating its crystallographic refinement.
Results. All-atom sampling of entire RNA chains, termini and loops is achieved
using the Richardson RNA backbone rotamer library and an unbiased distribution
for glycosidic dihedral angle. Sampling behaviour of Rappertk on a diverse
dataset of RNA chains under varying spatial restraints is benchmarked. The
iterative composite crystallographic refinement protocol developed here is
demonstrated to outperform CNS-only refinement on parts of tRNA(Asp) structure.
Conclusion. This work opens exciting possibilities for further work in RNA
modelling and crystallography.


  To counterbalance the views presented here by Suzana Moss de Oliveira, we
explain here the truth: How men are oppressed by Mother Nature, who may have
made an error inventing us, and by living women, who could get rid of most of
us. Why do women live longer than us? Why is the Y chromosome for men so small?
What are the dangers of marital fidelity? In an appendix we mention the
demographic challenges of the future with many old and few young people.


  The impact of environmental fluctuation on species diversity is studied with
a model of the evolutionary ecology of microorganisms. We show that
environmental fluctuation induces evolutionary branching and assures the
consequential coexistence of multiple species. Pairwise invasibility analysis
is applied to illustrate the speciation process. We also discuss how
fluctuation affects species diversity.


  We explore the large-scale behavior of nucleotide compositional strand
asymmetries along human chromosomes. As we observe for 7 of 9 origins of
replication experimentally identified so far, the (TA+GC) skew displays rather
sharp upward jumps, with a linear decreasing profile in between two successive
jumps. We present a model of replication with well positioned replication
origins and random terminations that accounts for the observed characteristic
serrated skew profiles. We succeed in identifying 287 pairs of putative
adjacent replication origins with an origin spacing approximately 1-2 Mbp that
are likely to correspond to replication foci observed in interphase nuclei and
recognized as stable structures that persist throughout subsequent cell
generations.


  Motivation: Measurements of gene expression over time enable the
reconstruction of transcriptional networks. However, Bayesian networks and many
other current reconstruction methods rely on assumptions that conflict with the
differential equations that describe transcriptional kinetics. Practical
approximations of kinetic models would enable inferring causal relationships
between genes from expression data of microarray, tag-based and conventional
platforms, but conclusions are sensitive to the assumptions made.
  Results: The representation of a sufficiently large portion of genome enables
computation of an upper bound on how much confidence one may place in
influences between genes on the basis of expression data. Information about
which genes encode transcription factors is not necessary but may be
incorporated if available. The methodology is generalized to cover cases in
which expression measurements are missing for many of the genes that might
control the transcription of the genes of interest. The assumption that the
gene expression level is roughly proportional to the rate of translation led to
better empirical performance than did either the assumption that the gene
expression level is roughly proportional to the protein level or the Bayesian
model average of both assumptions.
  Availability: http://www.oisb.ca points to R code implementing the methods (R
Development Core Team 2004).
  Supplementary information: http://www.davidbickel.com


  Although conservation of thermodynamics is much less studied than of
sequences and structures, thermodynamic details are biophysical features
different from but as important as structural features. As a succession of
previous research which revealed the important relationships between
thermodynamic features and miRNA maturization, this article applies
interspecies conservation of miRNA thermodynamics and structures to study miRNA
maturization. According to a thermodynamic and structural consensus principle,
miRBase is categorized by conservation subgroups, which imply various
functions. These subgroups are divided without the introduction of functional
information. This suggests the consistency between the two processes of miRNA
maturization and functioning. Different from prevailing methods which predict
extended miRNA precursors, a learning-based algorithm is proposed to predict
~22bp mature parts of 2780 test miRNA genes in 44 species with a rate of 79.4%.
This is the first attempt of a general interspecies prediction of mature
miRNAs. Suboptimal structures that most fit the consensus thermodynamic and
structural profiles are chosen to improve structure prediction. Distribution of
miRNA locations on corresponding pri-miRNA stem-loop structures is then
studied. Existing research on Drosha cleavage site is not generally true across
species. Instead, the distance between mature miRNA and center loop normalized
by stem length is a more conserved structural feature in animals, and the
normalized distance between mature miRNA and ss-RNA tail is the counterpart in
plants. This suggests two possibly-updating mechanisms of miRNA maturization in
animals and plants. All in all, conservations of thermodynamics together with
other features are shown closely related to miRNA maturization.


  The concept of robustness of regulatory networks has been closely related to
the nature of the interactions among genes, and the capability of pattern
maintenance or reproducibility. Defining this robustness property is a
challenging task, but mathematical models have often associated it to the
volume of the space of admissible parameters. Not only the volume of the space
but also its topology and geometry contain information on essential aspects of
the network, including feasible pathways, switching between two parallel
pathways or distinct/disconnected active regions of parameters. A general
method is presented here to characterize the space of admissible parameters, by
writing it as a semi-algebraic set, and then theoretically analyzing its
topology and geometry, as well as volume. This method provides a more objective
and complete measure of the robustness of a developmental module. As an
illustration, the segment polarity gene network is analyzed.


  We study the influence of a structural disorder on the thermodynamical
properties of 2D-elastic chains submitted to mechanical/topological constraint
as loops. The disorder is introduced via a spontaneous curvature whose
distribution along the chain presents either no correlation or long-range
correlations (LRC). The equilibrium properties of the one-loop system are
derived numerically and analytically for weak disorder. LRC are shown to favor
the formation of small loop, larger the LRC, smaller the loop size. We use the
mean first passage time formalism to show that the typical short time loop
dynamics is superdiffusive in the presence of LRC. Potential biological
implications on nucleosome positioning and dynamics in eukaryotic chromatin are
discussed.


  DNA torsion dynamics is essential in the transcription process; simple models
for it have been proposed by several authors, in particular Yakushevich (Y
model). These are strongly related to models of DNA separation dynamics such as
the one first proposed by Peyrard and Bishop (and developed by Dauxois, Barbi,
Cocco and Monasson among others), but support topological solitons. We recently
developed a ``composite'' version of the Y model, in which the sugar-phosphate
group and the base are described by separate degrees of freedom. This at the
same time fits experimental data better than the simple Y model, and shows
dynamical phenomena, which are of interest beyond DNA dynamics. Of particular
relevance are the mechanism for selecting the speed of solitons by tuning the
physical parameters of the non linear medium and the hierarchal separation of
the relevant degrees of freedom in ``master'' and ``slave''. These mechanisms
apply not only do DNA, but also to more general macromolecules, as we show
concretely by considering polyethylene.


  Nested Effects Models (NEMs) are a class of graphical models introduced to
analyze the results of gene perturbation screens. NEMs explore noisy subset
relations between the high-dimensional outputs of phenotyping studies, e.g. the
effects showing in gene expression profiles or as morphological features of the
perturbed cell.
  In this paper we expand the statistical basis of NEMs in four directions:
First, we derive a new formula for the likelihood function of a NEM, which
generalizes previous results for binary data. Second, we prove model
identifiability under mild assumptions. Third, we show that the new formulation
of the likelihood allows to efficiently traverse model space. Fourth, we
incorporate prior knowledge and an automated variable selection criterion to
decrease the influence of noise in the data.


  Among the multiple steps constituting the kinesin's mechanochemical cycle,
one of the most interesting events is observed when kinesins move an 8-nm step
from one microtubule (MT)-binding site to another. The stepping motion that
occurs within a relatively short time scale (~100 microsec) is, however, beyond
the resolution of current experiments, therefore a basic understanding to the
real-time dynamics within the 8-nm step is still lacking. For instance, the
rate of power stroke (or conformational change), that leads to the
undocked-to-docked transition of neck-linker, is not known, and the existence
of a substep during the 8-nm step still remains a controversial issue in the
kinesin community. By using explicit structures of the kinesin dimer and the MT
consisting of 13 protofilaments (PFs), we study the stepping dynamics with
varying rates of power stroke (kp). We estimate that 1/kp <~ 20 microsec to
avoid a substep in an averaged time trace. For a slow power stroke with 1/kp>20
microsec, the averaged time trace shows a substep that implies the existence of
a transient intermediate, which is reminiscent of a recent single molecule
experiment at high resolution. We identify the intermediate as a conformation
in which the tethered head is trapped in the sideway binding site of the
neighboring PF. We also find a partial unfolding (cracking) of the binding
motifs occurring at the transition state ensemble along the pathways prior to
binding between the kinesin and MT.


  Within the frame of an effective, coarse-grained hydrophobic-polar protein
model, we employ multicanonical Monte Carlo simulations to investigate
free-energy landscapes and folding channels of exemplified heteropolymer
sequences, which are permutations of each other. Despite the simplicity of the
model, the knowledge of the free-energy landscape in dependence of a suitable
system order parameter enables us to reveal complex folding characteristics
known from real bioproteins and synthetic peptides, such as two-state folding,
folding through weakly stable intermediates, and glassy metastability.


  We investigate stochastic extinction in an epidemic model and the impact of
random vaccinations in large populations formulated in terms of an optimal
escape path. We find that different random vaccination strategies can have
widely different results in decreasing expected time till extinction, for the
same total amount of vaccines used. Vaccination strategies are considered in
terms of two parameters: average frequency of vaccinations, given by $\gamma$,
and the amplitude of the vaccinations, $\epsilon$, where $\epsilon \ll 1$
refers to the proportion of the population being vaccinated at some particular
instant. It is found that while the average number of individuals vaccinated
per unit time, $\gamma \epsilon$, is kept constant, the particular values of
$\gamma$ and $\epsilon$ can play a highly significant role in increasing the
chance of epidemic extinction. The findings suggest that expected time till
extinction can be significantly shortened if less frequent vaccinations occur
in larger groups, corresponding to low $\gamma$, high $\epsilon$ strategy.


  Local protein structure analysis is informative to protein structure analysis
and has been used successfully in protein structure prediction and others.
Proteins have recurring structural features, such as helix caps and beta turns,
which often have strong amino acid sequence preferences. And the challenges for
local structure analysis have been identification and assignment of such common
short structural motifs.
  This paper proposes a new mathematical framework that can be applied to
analysis of the local structure of proteins, where local conformations of
protein backbones are described using differential geometry of folded
tetrahedron sequences. Using the framework, we could capture the recurring
structural features without any structural templates, which makes local
structure analysis not only simpler, but also more objective.
  Programs and examples are available from http://www.genocript.com .


  We study pore blockade times for a translocating polymer of length $N$,
driven by a field $E$ across the pore in three dimensions. The polymer performs
Rouse dynamics, i.e., we consider polymer dynamics in the absence of
hydrodynamical interactions. We find that the typical time the pore remains
blocked during a translocation event scales as $\sim N^{(1+2\nu)/(1+\nu)}/E$,
where $\nu\simeq0.588$ is the Flory exponent for the polymer. In line with our
previous work, we show that this scaling behaviour stems from the polymer
dynamics at the immediate vicinity of the pore -- in particular, the memory
effects in the polymer chain tension imbalance across the pore. This result,
along with the numerical results by several other groups, violates the lower
bound $\sim N^{1+\nu}/E$ suggested earlier in the literature. We discuss why
this lower bound is incorrect and show, based on conservation of energy, that
the correct lower bound for the pore-blockade time for field-driven
translocation is given by $\eta N^{2\nu}/E$, where $\eta$ is the viscosity of
the medium surrounding the polymer.


  We are interested in a stochastic model of trait and age-structured
population undergoing mutation and selection. We start with a continuous time,
discrete individual-centered population process. Taking the large population
and rare mutations limits under a well-chosen time-scale separation condition,
we obtain a jump process that generalizes the Trait Substitution Sequence
process describing Adaptive Dynamics for populations without age structure.
Under the additional assumption of small mutations, we derive an age-dependent
ordinary differential equation that extends the Canonical Equation. These
evolutionary approximations have never been introduced to our knowledge. They
are based on ecological phenomena represented by PDEs that generalize the
Gurtin-McCamy equation in Demography. Another particularity is that they
involve a fitness function, describing the probability of invasion of the
resident population by the mutant one, that can not always be computed
explicitly. Examples illustrate how adding an age-structure enrich the
modelling of structured population by including life history features such as
senescence. In the cases considered, we establish the evolutionary
approximations and study their long time behavior and the nature of their
evolutionary singularities when computation is tractable. Numerical procedures
and simulations are carried.


  The popular neighbor-joining (NJ) algorithm used in phylogenetics is a greedy
algorithm for finding the balanced minimum evolution (BME) tree associated to a
dissimilarity map. From this point of view, NJ is ``optimal'' when the
algorithm outputs the tree which minimizes the balanced minimum evolution
criterion. We use the fact that the NJ tree topology and the BME tree topology
are determined by polyhedral subdivisions of the spaces of dissimilarity maps
${\R}_{+}^{n \choose 2}$ to study the optimality of the neighbor-joining
algorithm. In particular, we investigate and compare the polyhedral
subdivisions for $n \leq 8$. A key requirement is the measurement of volumes of
spherical polytopes in high dimension, which we obtain using a combination of
Monte Carlo methods and polyhedral algorithms. We show that highly unrelated
trees can be co-optimal in BME reconstruction, and that NJ regions are not
convex. We obtain the $l_2$ radius for neighbor-joining for $n=5$ and we
conjecture that the ability of the neighbor-joining algorithm to recover the
BME tree depends on the diameter of the BME tree.


  Questions of understanding and quantifying the representation and amount of
information in organisms have become a central part of biological research, as
they potentially hold the key to fundamental advances. In this paper, we
demonstrate the use of information-theoretic tools for the task of identifying
segments of biomolecules (DNA or RNA) that are statistically correlated. We
develop a precise and reliable methodology, based on the notion of mutual
information, for finding and extracting statistical as well as structural
dependencies. A simple threshold function is defined, and its use in
quantifying the level of significance of dependencies between biological
segments is explored. These tools are used in two specific applications. First,
for the identification of correlations between different parts of the maize
zmSRp32 gene. There, we find significant dependencies between the 5'
untranslated region in zmSRp32 and its alternatively spliced exons. This
observation may indicate the presence of as-yet unknown alternative splicing
mechanisms or structural scaffolds. Second, using data from the FBI's Combined
DNA Index System (CODIS), we demonstrate that our approach is particularly well
suited for the problem of discovering short tandem repeats, an application of
importance in genetic profiling.


  Interconvertible enzyme cascades, exemplified by the mitogen activated
protein kinase (MAPK) cascade, are a frequent mechanism in signal transduction
pathways. There has been much speculation as to the role of these pathways, and
how their structure is related to their function. A common conclusion is that
the cascades serve to amplify biochemical signals so that a single bound ligand
molecule might produce a multitude of second messengers. Some recent work has
focused on a particular feature present in some MAPK pathways -- a negative
feedback loop which spans the length of the cascade. This is a feature that is
shared by a man-made engineering device, the feedback amplifier. We propose a
novel interpretation: that by wrapping a feedback loop around an amplifier,
these cascades may be acting as biochemical feedback amplifiers which imparts
i) increased robustness with respect to internal perturbations; ii) a linear
graded response over an extended operating range; iii) insulation from external
perturbation, resulting in functional modularization. We also report on the
growing list of experimental evidence which supports a graded response of MAPK
with respect to Epidermal Growth Factor. This evidence supports our hypothesis
that in these circumstances MAPK cascade, may be acting as a feedback
amplifier.


  We show how accurate kinetic information, such as the rates of protein
folding and unfolding, can be extracted from replica-exchange molecular
dynamics (REMD) simulations. From the brief and continuous trajectory segments
between replica exchanges, we estimate short-time propagators in conformation
space and use them to construct a master equation. For a helical peptide in
explicit water, we determine the rates of transitions both locally between
microscopic conformational states and globally for folding and unfolding. We
show that accurate rates in the ~1/(100 ns) to ~1/(1 ns) range can be obtained
from REMD with exchange times of 5 ps, in excellent agreement with results from
long equilibrium molecular dynamics.


  We present an approach, called the "Shadow Method," for the identification of
disease loci from dense genetic marker maps in complex, potentially incomplete
pedigrees. "Shadow" is a simple method based on an analysis of the patterns of
obligate meiotic recombination events in genotypic data. This method can be
applied to any high density marker map and was specifically designed to exploit
the fact that extremely dense marker maps are becoming more readily available.
We also describe how to interpret and associate meaningful P-Values to the
results. Shadow has significant advantages over traditional parametric linkage
analysis methods in that it can be readily applied even in cases in which the
topology of a pedigree or pedigrees can only be partially determined. In
addition, Shadow is robust to variability in a range of parameters and in
particular does not require prior knowledge of mode of inheritance, penetrance
or clinical misdiagnosis rate. Shadow can be used for any SNP data, but is
especially effective when applied to dense samplings. Our primary example uses
data from Affymetrix 100k SNPChip samples in which we illustrate our approach
by analyzing simulated data as well as genome-wide SNP data from two pedigrees
with inherited forms of kidney failure, one of which is compared with a typical
LOD score analysis.


  Models for RNA secondary structures (the topology of folded RNA) without
pseudo knots are disordered systems with a complex state-space below a critical
temperature. Hence, a complex dynamical (glassy) behavior can be expected, when
performing Monte Carlo simulation. Interestingly, in contrast to most other
complex systems, the ground states and the density of states can be computed in
polynomial time exactly using transfer matrix methods. Hence, RNA secondary
structure is an ideal model to study the relation between static/thermodynamic
properties and dynamics of algorithms. Also they constitute an ideal benchmark
system for new Monte Carlo methods.
  Here we considered three different recent Monte Carlo approaches: entropic
sampling using flat histograms, optimized-weights ensembles, and ParQ, which
estimates the density of states from transition matrices.
  These methods were examined by comparing the obtained density of states with
the exact results. We relate the complexity seen in the dynamics of the Monte
Carlo algorithms to static properties of the phase space by studying the
correlations between tunneling times, sampling errors, amount of meta-stable
states and degree of ultrametricity at finite temperature.


  What could cause the emergence of non-encoding codons in the course of
evolution of the genetic code? Hypothesis of evolution of the genetic code from
GC to the AGUC-alphabet account for existence of stop codons.


  Models for species-area relationships up to now have focused on the mean
richness as a function of area. We present MFp1p2, a self-similar multifractal.
It explicitly models both trend and variation in richness as a function of
area, and is a generalisation of the model of scaling of mean species richness
due to Harte et al (1999). The construction is based on a cascade of bisections
of a rectangle. The two parameters of the model are p1, the proportion of
species that occur in the richer half, and p2, the proportion of species that
occur in the poorer half. Equivalent parameterisations are a = (p1 + p2)/2 and
b = p1/p2. These parameters are interpreted as follows: a gives the scaling of
mean density, b gives the scaling of spatial variability. Several properties of
MFp1p2 are derived, a generalisation is noted and some applications are
suggested.


  We study the coupled dynamics of two populations of random replicators by
means of statistical mechanics methods, and focus on the effects of relative
population size, strategy correlations and heterogeneities in the respective
co-operation pressures. To this end we generalise existing path-integral
approaches to replicator systems with random asymmetric couplings. This
technique allows one to formulate an effective dynamical theory, which is exact
in the thermodynamic limit and which can be solve for persistent order
parameters in a fixed-point regime regardless of the symmetry of the
interactions. The onset of instability can be determined self-consistently. We
calculate quantities such as the diversity of the respective populations and
their fitnesses in the stationary state, and compare results with data from a
numerical integration of the replicator equations


  Time has come for a survey of our knowledge on the physical interaction
between the growth-promoting insulin molecule and retinoblastoma tumor
suppressor protein (RB). Theoretical and experimental observations over the
past 15 years reviewed here indicate that the insulin-RB dimer may represent an
essential molecular crossroads involved in major physiological and pathological
conditions. Within this system, the putative tumor suppressor insulin-degrading
enzyme (IDE) should be an important modulator. Perhaps most remarkably, the
abstraction of this encounter between insulin and RB, two growth-regulatory
giants acting either in concert or against each other depending on the
respective cellular requirements, reveals that Nature may compute in
controlling cell fate and we could follow in its footsteps towards developing
more efficient therapeutics as well as novel technical devices.


  Until presently, most approaches for understanding physiological or
pathological phenomena have been based on the assumption that these processes
start from individual cells, a concept introduced primarily by Virchow in the
19th century. Yet, it has also been increasingly recognized that this
perception is insufficient, at least when it comes to grasping the mechanisms
underlying largely incurable diseases such as metastatic cancer or rheumatoid
arthritis. Despite this insight, even recently founded disciplines such as
systems biology are still locked in this century-old mind-set of cellular
building blocks and thus predictably of limited usefulness. Other studies
conducted over the past years, however, suggest that there is something more
fundamental to life and its various conditions than the cell: peptide strings.
Here, I review the origin and nature of these sub- and trans-cellular elements
as well as their potential to provide the long-sought answers hitherto
inaccessible to cell biology.


  In systems biology new ways are required to analyze the large amount of
existing data on regulation of cellular processes. Recent work can be roughly
classified into either dynamical models of well-described subsystems, or
coarse-grained descriptions of the topology of the molecular networks at the
scale of the whole organism. In order to bridge these two disparate approaches
one needs to develop simplified descriptions of dynamics and topological
measures which address the propagation of signals in molecular networks. Here,
we consider the directed network of protein regulation in E. coli,
characterizing its modularity in terms of its potential to transmit signals. We
demonstrate that the simplest measure based on identifying sub-networks of
strong components, within which each node could send a signal to every other
node, indeed partitions the network into functional modules. We then suggest
measures to quantify the cost and spread associated with sending a signal
between any particular pair of proteins. Thereby, we address the signalling
specificity within and between modules, and show that in the regulation of
E.coli there is a systematic reduction of the cost and spread for signals
traveling over more than two intermediate reactions.


  The Caspar-Klug classification of viruses whose protein shell, called viral
capsid, exhibits icosahedral symmetry, has recently been extended to
incorporate viruses whose capsid proteins are exclusively organised in
pentamers. The approach, named `Viral Tiling Theory', is inspired by the theory
of quasicrystals, where aperiodic Penrose tilings enjoy 5-fold and 10-fold
local symmetries. This paper analyzes the extent to which this classification
approach informs dynamical properties of the viral capsids, in particular the
pattern of Raman active modes of vibrations, which can be observed
experimentally.


  Genetic competence is a phenotypic state of a bacterial cell in which it is
capable of importing DNA, presumably to hasten its exploration of alternate
genes in its quest for survival under stress. Recently, it was proposed that
this transition is uncorrelated among different cells in the colony. Motivated
by several discovered signaling mechanisms which create colony-level responses,
we present a model for the influence of quorum-sensing signals on a colony of
B. Subtilis cells during the transition to genetic competence. Coupling to the
external signal creates an effective inhibitory mechanism, which results in
anti-correlation between the cycles of adjacent cells. We show that this
scenario is consistent with the specific experimental measurement, which fails
to detect some underlying collective signaling mechanisms. Rather, we suggest
other parameters that should be used to verify the role of a quorum-sensing
signal. We also study the conditions under which phenotypic spatial patterns
may emerge.


  A primordial genetic code is proposed, having only four codons assigned, GGC
meaning glycine, GAC meaning aspartate/glutamate, GCC meaning alanine-like and
GUC meaning valine-like. Pathways of ambiguity reduction enlarged the codon
repertoire with CUC meaning leucine, AUC meaning isoleucine, ACC meaning
threonine-like and GAG meaning glutamate. Introduction of UNN anticodons, in a
next episode of code evolution in which nonsense elimination was the leading
theme, introduced a family box structure superposed on the original mirror
structure. Finally, growth rate was the leading theme during the remaining
repertoire expansion, explaining the ordered phylogenetic pattern of
aminoacyl-tRNA synthetases. The special role of natural aptamers in the process
is high-lighted, and the error robustness characteristics of the code are shown
to have evolved by way of a stepwise, restricted enlargement of the tRNA
repertoire, instead of by an exhaustive selection process testing myriads of
codes.


  The theory of transition between $\alpha$-helix, $\beta$-sheet and random
coil conformation of a protein is discussed through a simple model, that
includes both short and long-range interactions. Besides the bonding parameter
and helical initiation factor in Zimm-Bragg model, three new parameters are
introduced to describe beta structure: the local constraint factor for a single
residue to be contained in a $\beta$-strand, the long-range bonding parameter
that accounts for the interaction between a pair of bonded $\beta$-strands, and
a correction factor for the initiation of a $\beta$-sheet. Either increasing
local constraint factor or long-range bonding parameter can cause a transition
from $\alpha$-helix or random coil conformation to $\beta$-sheet conformation.
The sharpness of transition depends on the competition between short and
long-range interactions. Other effective factors, such as the chain length and
temperature, are also discussed. In this model, the entropy due to different
ways to group $\beta$-strands into different $\beta$-sheets gives rise to
significant contribution to partition function, and makes major differences
between beta structure and helical structure.


  We study the spreading of an infection within an SIS epidemiological model on
a network. Susceptible agents are given the opportunity of breaking their links
with infected agents. Broken links are either permanently removed or
reconnected with the rest of the population. Thus, the network coevolves with
the population as the infection progresses. We show that a moderate
reconnection frequency is enough to completely suppress the infection. A
partial, rather weak isolation of infected agents suffices to eliminate the
endemic state.


  We present a theoretical study of the folding of small proteins inside
confining potentials. Proteins are described in the framework of an effective
potential model which contains the Ramachandran angles as degrees of freedom
and does not need any {\it a priori} information about the native state.
Hydrogen bonds, dipole-dipole- and hydrophobic interactions are taken
explicitly into account. An interesting feature displayed by this potential is
the presence of some intermediates between the unfolded and native states. We
consider different types of confining potentials in order to study the
structural properties of proteins folding inside cages with repulsive or
attractive walls. Using the Wang-Landau algorithm we determine the density of
states (DOS) and analyze in detail the thermodynamical properties of the
confined proteins for different sizes of the cages. We show that confinement
dramatically reduces the phase space available to the protein and that the
presence of intermediate states can be controlled by varying the properties of
the confining potential. Cages with strongly attractive walls lead to the
disappearance of the intermediate states and to a two-state folding into a less
stable configuration. However, cages with slightly attractive walls make the
native structure more stable than in the case of pure repulsive potentials, and
the folding process occurs through intermediate configurations. In order to
test the metastable states we analyze the free energy landscapes as a function
of the configurational energy and of the end-to-end distance as an order
parameter.


  Position-specific scoring matrices (PSSMs) are useful for detecting weak
homology in protein sequence analysis, and they are thought to contain some
essential signatures of the protein families. In order to elucidate what kind
of ingredients constitute such family-specific signatures, we apply singular
value decomposition to a set of PSSMs and examine the properties of dominant
right and left singular vectors. The first right singular vectors were
correlated with various amino acid indices including relative mutability, amino
acid composition in protein interior, hydropathy, or turn propensity, depending
on proteins. A significant correlation between the first left singular vector
and a measure of site conservation was observed. It is shown that the
contribution of the first singular component to the PSSMs act to disfavor
potentially but falsely functionally important residues at conserved sites. The
second right singular vectors were highly correlated with hydrophobicity
scales, and the corresponding left singular vectors with contact numbers of
protein structures. It is suggested that sequence alignment with a PSSM is
essentially equivalent to threading supplemented with functional information.
The presented method may be used to separate functionally important sites from
structurally important ones, and thus it may be a useful tool for predicting
protein functions.


  A method to search for local structural similarities in proteins at atomic
resolution is presented. It is demonstrated that a huge amount of structural
data can be handled within a reasonable CPU time by using a conventional
relational database management system with appropriate indexing of geometric
data. This method, which we call geometric indexing, can enumerate ligand
binding sites that are structurally similar to sub-structures of a query
protein among more than 160,000 possible candidates within a few hours of CPU
time on an ordinary desktop computer. After detecting a set of high scoring
ligand binding sites by the geometric indexing search, structural alignments at
atomic resolution are constructed by iteratively applying the Hungarian
algorithm, and the statistical significance of the final score is estimated
from an empirical model based on a gamma distribution. Applications of this
method to several protein structures clearly shows that significant
similarities can be detected between local structures of non-homologous as well
as homologous proteins.


  We present the results of our numerical analysis of a "composite" model of
DNA which generalizes a well-known elementary torsional model of Yakushevich by
allowing bases to move independently from the backbone. The model shares with
the Yakushevich model many features and results but it represents an
improvement from both the conceptual and the phenomenological point of view. It
provides a more realistic description of DNA and possibly a justification for
the use of models which consider the DNA chain as uniform. It shows that the
existence of solitons is a generic feature of the underlying nonlinear dynamics
and is to a large extent independent of the detailed modelling of DNA. As
opposite to the Yakushevich model, where it is needed to use an unphysical
value for the torsion in order to induce the correct velocity of sound, the
model we consider supports solitonic solutions, qualitatively and
quantitatively very similar to the Yakushevich solitons, in a fully realistic
range of all the physical parameters characterizing the DNA.


  Transcriptional pulsing has been observed in both prokaryotes and eukaryotes
and plays a crucial role in cell to cell variability of protein and mRNA
numbers. The issue is how the time constants associated with episodes of
transcriptional bursting impact cellular mRNA and protein distributions and
reciprocally, to what extent experimentally observed distributions can be
attributed to transcriptional pulsing. We address these questions by
investigating the exact time-dependent solution of the Master equation for a
transcriptional pulsing model of mRNA distributions. We find a plethora of
results: we show that, among others, bimodal and long-tailed (power law)
distributions occur in the steady state as the rate constants are varied over
biologically significant time scales. Since steady state distributions may not
be reached experimentally we present results for the time evolution of the
distributions. Because cellular behavior is essentially determined by proteins,
we investigate the effect of the different mRNA distributions on the
corresponding protein distributions. We delineate the regimes of rate constants
for which the protein distribution mimics the mRNA distribution and those for
which the protein distribution deviates significantly from the mRNA
distribution.


  We examine a simple model of proton pumping through the inner membrane of
mitochondria in the living cell. We demonstrate that the pumping process can be
described using approaches of condensed matter physics. In the framework of
this model, we show that the resonant F\"orster-type energy exchange due to
electron-proton Coulomb interaction can provide an unidirectional flow of
protons against an electrochemical proton gradient, thereby accomplishing
proton pumping. The dependence of this effect on temperature as well as
electron and proton voltage build-ups are obtained taking into account
electrostatic forces and noise in the environment. We find that the proton pump
works with maximum efficiency in the range of temperatures and transmembrane
electrochemical potentials which correspond to the parameters of living cells.


  The proliferation and migration dichotomy of the tumor cell invasion is
examined within a two-component continuous time random walk (CTRW) model. The
balance equations for the cancer cells of two phenotypes with random switching
between cell proliferation and migration are derived. The transport of tumor
cells is formulated in terms of the CTRW with an arbitrary waiting time
distribution law, while proliferation is modelled by a logistic growth. The
overall rate of tumor cell invasion for normal diffusion and subdiffusion is
determined.


  In this paper we enumerate $k$-noncrossing RNA pseudoknot structures with
given minimum stack-length. We show that the numbers of $k$-noncrossing
structures without isolated base pairs are significantly smaller than the
number of all $k$-noncrossing structures. In particular we prove that the
number of 3- and 4-noncrossing RNA structures with stack-length $\ge 2$ is for
large $n$ given by $311.2470 \frac{4!}{n(n-1)...(n-4)}2.5881^n$ and $1.217\cdot
10^{7} n^{-{21/2}} 3.0382^n$, respectively. We furthermore show that for
$k$-noncrossing RNA structures the drop in exponential growth rates between the
number of all structures and the number of all structures with stack-size $\ge
2$ increases significantly. Our results are of importance for prediction
algorithms for pseudoknot-RNA and provide evidence that there exist neutral
networks of RNA pseudoknot structures.


  Robustness to mutations and noise has been shown to evolve through
stabilizing selection for optimal phenotypes in model gene regulatory networks.
The ability to evolve robust mutants is known to depend on the network
architecture. How do the dynamical properties and state-space structures of
networks with high and low robustness differ? Does selection operate on the
global dynamical behavior of the networks? What kind of state-space structures
are favored by selection? We provide damage propagation analysis and an
extensive statistical analysis of state spaces of these model networks to show
that the change in their dynamical properties due to stabilizing selection for
optimal phenotypes is minor. Most notably, the networks that are most robust to
both mutations and noise are highly chaotic. Certain properties of chaotic
networks, such as being able to produce large attractor basins, can be useful
for maintaining a stable gene-expression pattern. Our findings indicate that
conventional measures of stability, such as the damage-propagation rate, do not
provide much information about robustness to mutations or noise in model gene
regulatory networks.


  Network science is already making an impact on the study of complex systems
and offers a promising variety of tools to understand their formation and
evolution (1-4) in many disparate fields from large communication networks
(5,6), transportation infrastructures (7) and social communities (8,9) to
biological systems (1,10,11). Even though new highthroughput technologies have
rapidly been generating large amounts of genomic data, drug design has not
followed the same development, and it is still complicated and expensive to
develop new single-target drugs. Nevertheless, recent approaches suggest that
multi-target drug design combined with a network-dependent approach and
large-scale systems-oriented strategies (12-14) create a promising framework to
combat complex multigenetic disorders like cancer or diabetes. Here, we
investigate the human network corresponding to the interactions between all US
approved drugs and human therapies, defined by known drug-therapy
relationships. Our results show that the key paths in this network are shorter
than three steps, indicating that distant therapies are separated by a
surprisingly low number of chemical compounds. We also identify a sub-network
composed by drugs with high centrality measures (15), which represent the
structural back-bone of the drug-therapy system and act as hubs routing
information between distant parts of the network. These findings provide for
the first time a global map of the largescale organization of all known drugs
and associated therapies, bringing new insights on possible strategies for
future drug development. Special attention should be given to drugs which
combine the two properties of (a) having a high centrality value and (b) acting
on multiple targets.


  We critically examine the evidence for the idea that encephalization
quotients increase with time. We find that human-like intelligence is not a
convergent feature of evolution. Implications for the search for
extraterrestrial intelligence are discussed.


  We consider the evolution of an asexually reproducing population in an
uncorrelated random fitness landscape in the limit of infinite genome size,
which implies that each mutation generates a new fitness value drawn from a
probability distribution $g(w)$. This is the finite population version of
Kingman's house of cards model [J.F.C. Kingman, \textit{J. Appl. Probab.}
\textbf{15}, 1 (1978)]. In contrast to Kingman's work, the focus here is on
unbounded distributions $g(w)$ which lead to an indefinite growth of the
population fitness. The model is solved analytically in the limit of infinite
population size $N \to \infty$ and simulated numerically for finite $N$. When
the genome-wide mutation probability $U$ is small, the long time behavior of
the model reduces to a point process of fixation events, which is referred to
as a \textit{diluted record process} (DRP). The DRP is similar to the standard
record process except that a new record candidate (a number that exceeds all
previous entries in the sequence) is accepted only with a certain probability
that depends on the values of the current record and the candidate. We develop
a systematic analytic approximation scheme for the DRP. At finite $U$ the
fitness frequency distribution of the population decomposes into a stationary
part due to mutations and a traveling wave component due to selection, which is
shown to imply a reduction of the mean fitness by a factor of $1-U$ compared to
the $U \to 0$ limit.


  Increased permeability of the blood-brain barrier (BBB) is important in
neurological disorders. Neuroinflammation is associated with increased BBB
breakdown and brain injury. Tumor necrosis factor-a (TNF-a) is involved in BBB
injury and edema formation through a mechanism involving matrix
metalloproteinase (MMP) upregulation. There is emerging evidence indicating
that cyclooxygenase (COX) inhibition limits BBB disruption following ischemic
stroke and bacterial meningitis, but the mechanisms involved are not known. We
used intracerebral injection of TNF-a to study the effect of COX inhibition on
TNF-a-induced BBB breakdown, MMP expression/activity and oxidative stress. BBB
disruption was evaluated by the uptake of 14C-sucrose into the brain and by
magnetic resonance imaging (MRI) utilizing Gd-DTPA as a paramagnetic contrast
agent. Using selective inhibitors of each COX isoform, we found that COX-1
activity is more important than COX-2 in BBB opening. TNF-a induced a
significant upregulation of gelatinase B (MMP-9), stromelysin-1 (MMP-3) and
COX-2. In addition, TNF-a significantly depleted glutathione as compared to
saline. Indomethacin (10 mg/kg; i.p.), an inhibitor of COX-1 and COX-2, reduced
BBB damage at 24 h. Indomethacin significantly attenuated MMP-9 and MMP-3
expression and activation, and prevented the loss of endogenous radical
scavenging capacity following intracerebral injection of TNF-a. Our results
show for the first time that BBB disruption during neuroinflammation can be
significantly reduced by administration of COX inhibitors. Modulation of COX in
brain injury by COX inhibitors or agents modulating prostaglandin E2
formation/signaling may be useful in clinical settings associated with BBB
disruption.


  The concept of (auto)catalytic systems has become a cornerstone in
understanding evolutionary processes in various fields. The common ground is
the observation that for the production of new species/goods/ideas/elements
etc. the pre-existence of specific other elements is a necessary condition. In
previous work some of us showed that the dynamics of the catalytic network
equation can be understood in terms of topological recurrence relations paving
a path towards the analytic tractability of notoriously high dimensional
evolution equations. We apply this philosophy to studies in socio-physics,
bio-diversity and massive events of creation and destruction in technological
and biological networks. Cascading events, triggered by small exogenous
fluctuations, lead to dynamics strongly resembling the qualitative picture of
Schumpeterian economic evolution. Further we show that this new methodology
allows to mathematically treat a variant of the threshold voter-model of
opinion formation on networks. For fixed topology we find distinct phases of
mixed opinions and consensus.


  The relationship between the regulatory design and the functionality of
molecular networks is a key issue in biology. Modules and motifs have been
associated to various cellular processes, thereby providing anecdotal evidence
for performance based localization on molecular networks. To quantify
structure-function relationship we investigate similarities of proteins which
are close in the regulatory network of the yeast Saccharomyces Cerevisiae. We
find that the topology of the regulatory network show weak remnants of its
history of network reorganizations, but strong features of co-regulated
proteins associated to similar tasks. This suggests that local topological
features of regulatory networks, including broad degree distributions, emerge
as an implicit result of matching a number of needed processes to a finite
toolbox of proteins.


  Motivated by rolling adhesion of white blood cells in the vasculature, we
study how cells move in linear shear flow above a wall to which they can adhere
via specific receptor-ligand bonds. Our computer simulations are based on a
Langevin equation accounting for hydrodynamic interactions, thermal
fluctuations and adhesive interactions. In contrast to earlier approaches, our
model not only includes stochastic rules for the formation and rupture of
bonds, but also fully resolves both receptor and ligand positions. We identify
five different dynamic states of motion in regard to the translational and
angular velocities of the cell. The transitions between the different states
are mapped out in a dynamic state diagram as a function of the rates for bond
formation and rupture. For example, as the cell starts to adhere under the
action of bonds, its translational and angular velocities become synchronized
and the dynamic state changes from slipping to rolling. We also investigate the
effect of non-molecular parameters. In particular, we find that an increase in
viscosity of the medium leads to a characteristic expansion of the region of
stable rolling to the expense of the region of firm adhesion, but not to the
expense of the regions of free or transient motion. Our results can be used in
an inverse approach to determine single bond parameters from flow chamber data
on rolling adhesion.


  The effects of the finite size of the network on the evolutionary dynamics of
a Boolean network are analyzed. In the model considered, Boolean networks
evolve via a competition between nodes that punishes those in the majority. It
is found that finite size networks evolve in a fundamentally different way than
infinitely large networks do. The symmetry of the evolutionary dynamics of
infinitely large networks that selects for canalizing Boolean functions is
broken in the evolutionary dynamics of finite size networks. In finite size
networks there is an additional selection for input inverting Boolean functions
that output a value opposite to the majority of input values. These results are
revealed through an empirical study of the model that calculates the frequency
of occurrence of the different possible Boolean functions. Classes of functions
are found to occur with the same frequency. Those classes depend on the
symmetry of the evolutionary dynamics and correspond to orbits of the relevant
symmetry group. The empirical results match analytic results, determined by
utilizing Polya's theorem, for the number of orbits expected in both finite
size and infinitely large networks. The reason for the symmetry breaking in the
evolutionary dynamics is found to be due to the need for nodes in finite size
networks to behave differently in order to cooperate so that the system
collectively performs as well as possible. The results suggest that both finite
size effects and symmetry are important for understanding the evolution of
real-world complex networks, including genetic regulatory networks.


  In this paper we enumerate $k$-noncrossing RNA pseudoknot structures with
given minimum arc- and stack-length. That is, we study the numbers of RNA
pseudoknot structures with arc-length $\ge 3$, stack-length $\ge \sigma$ and in
which there are at most $k-1$ mutually crossing bonds, denoted by ${\sf
T}_{k,\sigma}^{[3]}(n)$. In particular we prove that the numbers of 3, 4 and
5-noncrossing RNA structures with arc-length $\ge 3$ and stack-length $\ge 2$
satisfy ${\sf T}_{3,2}^{[3]}(n)^{}\sim K_3 n^{-5} 2.5723^n$, ${\sf
T}^{[3]}_{4,2}(n)\sim K_4 n^{-{21/2}} 3.0306^n$, and ${\sf
T}^{[3]}_{5,2}(n)\sim K_5 n^{-18} 3.4092^n$, respectively, where $K_3,K_4,K_5$
are constants. Our results are of importance for prediction algorithms for RNA
pseudoknot structures.