  The nuclear pore supports molecular communication between cytoplasm and
nucleus in eukaryotic cells. Selective transport of proteins is mediated by
soluble receptors, whose regulation by the small GTPase Ran leads to cargo
accumulation in, or depletion from the nucleus, i.e., nuclear import or nuclear
export. We consider the operation of this transport system by a combined
analytical and experimental approach. Provocative predictions of a simple model
were tested using cell-free nuclei reconstituted in Xenopus egg extract, a
system well suited to quantitative studies. We found that accumulation capacity
is limited, so that introduction of one import cargo leads to egress of
another. Clearly, the pore per se does not determine transport directionality.
Moreover, different cargo reach a similar ratio of nuclear to cytoplasmic
concentration in steady-state. The model shows that this ratio should in fact
be independent of the receptor-cargo affinity, though kinetics may be strongly
influenced. Numerical conservation of the system components highlights a
conflict between the observations and the popular concept of transport cycles.
We suggest that chemical partitioning provides a framework to understand the
capacity to generate concentration gradients by equilibration of the
receptor-cargo intermediary.


  We introduce a toy model of ecosystem assembly for which we are able to map
out all assembly pathways generated by external invasions. The model allows to
display the whole phase space in the form of an assembly graph whose nodes are
communities of species and whose directed links are transitions between them
induced by invasions. We characterize the process as a finite Markov chain and
prove that it exhibits a unique set of recurrent states (the endstate of the
process), which is therefore resistant to invasions. This also shows that the
endstate is independent on the assembly history. The model shares all features
with standard assembly models reported in the literature, with the advantage
that all observables can be computed in an exact manner.


  We investigate the distribution of flavonoids, a major category of plant
secondary metabolites, across species. Flavonoids are known to show high
species specificity, and were once considered as chemical markers for
understanding adaptive evolution and characterization of living organisms. We
investigate the distribution among species using bipartite networks, and find
that two heterogeneous distributions are conserved among several families: the
power-law distributions of the number of flavonoids in a species and the number
of shared species of a particular flavonoid. In order to explain the possible
origin of the heterogeneity, we propose a simple model with, essentially, a
single parameter. As a result, we show that two respective power-law statistics
emerge from simple evolutionary mechanisms based on a multiplicative process.
These findings provide insights into the evolution of metabolite diversity and
characterization of living organisms that defy genome sequence analysis for
different reasons.


  Background: Hybridogenesis is a very interesting example of reproduction
which seems to integrate the sexual and clonal processes in one system. In a
case of frogs, described in the paper, two parental species - Rana lessonae and
Rana ridibunda can form fertile hybrid individuals - Rana esculenta. Hybrid
individuals eliminate one parental haplotype from their germ line cells before
meiosis (end before recombination) which implicates clonal reproduction of the
haplotype transferred to the gametes. All three "species" are called "complex
species". To study the evolution of genomes in the hybridogenetic fraction of
this complex species we have used the Monte Carlo based model rendering the age
structured populations. The model enables the analysis of distribution of
defective alleles in the individual genomes as well as in the genetic pool of
the whole populations.
  Results: We have shown that longer isolation of hybrids' populations leads to
the speciation through emerging the specific sets of complementing haplotypes
in their genetic pool. The fraction of defective alleles increases but the
defects are complemented in the heterozygous loci. Nevertheless, even small
supply of new hybrids generated by the two parental species or crossbreeding
between hybrids and one of the parental species prevents the speciation and
changes the strategy of the genome evolution from the complementing to the
purifying Darwinian selection.


  Epidemic models with inhomogeneous populations have been used to study major
outbreaks and recently Britton and Lindenstrand \cite{BL} described the case
when latency and infectivity have independent gamma distributions. They found
that variability in these random variables had opposite effects on the epidemic
growth rate. That rate increased with greater variability in latency but
decreased with greater variability in infectivity. Here we extend their result
by using the McKay bivariate gamma distribution for the joint distribution of
latency and infectivity, recovering the above effects of variability but
allowing possible correlation. We use methods of stochastic rate processes to
obtain explicit solutions for the growth of the epidemic and the evolution of
the inhomogeneity and information entropy. We obtain a closed analytic solution
to the evolution of the distribution of the number of uninfected individuals as
the epidemic proceeds, and a concomitant expression for the decay of entropy.
The family of McKay bivariate gamma distributions has a tractable information
geometry which provides a framework in which the evolution of distributions can
be studied as the outbreak grows, with a natural distance structure for
quantitative tracking of progress.


  Computer simulations of the Penna ageing model suggest that already a small
fraction of births with enhanced number of new mutations can negatively
influence the whole population.


  The AGBNP2 implicit solvent model, an evolution of the Analytical Generalized
Born plus Non-Polar (AGBNP) model we have previously reported, is presented
with the aim of modeling hydration effects beyond those described by
conventional continuum dielectric representations. A new empirical hydration
free energy component based on a procedure to locate and score hydration sites
on the solute surface is introduced to model first solvation shell effects,
such as hydrogen bonding, which are poorly described by continuum dielectric
models. This new component is added to the Generalized Born and non-polar AGBNP
models which have been improved with respect to the description of the solute
volume description. We have introduced an analytical Solvent Excluding Volume
(SEV) model which reduces the effect of spurious high-dielectric interstitial
spaces present in conventional van der Waals representations of the solute
volume. The AGBNP2 model is parametrized and tested with respect to
experimental hydration free energies of small molecules and the results of
explicit solvent simulations. Modeling the granularity of water is one of the
main design principles employed for the the first shell solvation function and
the SEV model, by requiring that water locations have a minimum available
volume based on the size of a water molecule. We show that the new volumetric
model produces Born radii and surface areas in good agreement with accurate
numerical evaluations. The results of Molecular Dynamics simulations of a
series of mini-proteins show that the new model produces conformational
ensembles in substantially better agreement with reference explicit solvent
ensembles than the original AGBNP model with respect to both structural and
energetics measures.


  ChIP-seq, which combines chromatin immunoprecipitation with massively
parallel short-read sequencing, can profile in vivo genome-wide transcription
factor-DNA association with higher sensitivity, specificity and spatial
resolution than ChIP-chip. While it presents new opportunities for research,
ChIP-seq poses new challenges for statistical analysis that derive from the
complexity of the biological systems characterized and the variability and
biases in its digital sequence data. We propose a method called PICS
(Probabilistic Inference for ChIP-seq) for extracting information from ChIP-seq
aligned-read data in order to identify regions bound by transcription factors.
PICS identifies enriched regions by modeling local concentrations of
directional reads, and uses DNA fragment length prior information to
discriminate closely adjacent binding events via a Bayesian hierarchical
t-mixture model. Its per-event fragment length estimates also allow it to
remove from analysis regions that have atypical lengths. PICS uses
pre-calculated, whole-genome read mappability profiles and a truncated
t-distribution to adjust binding event models for reads that are missing due to
local genome repetitiveness. It estimates uncertainties in model parameters
that can be used to define confidence regions on binding event locations and to
filter estimates. Finally, PICS calculates a per-event enrichment score
relative to a control sample, and can use a control sample to estimate a false
discovery rate. We compared PICS to the alternative methods MACS, QuEST, and
CisGenome, using published GABP and FOXA1 data sets from human cell lines, and
found that PICS' predicted binding sites were more consistent with
computationally predicted binding motifs.


  We improve the lower bound on the extremal version of the Maximum Agreement
Subtree problem. Namely we prove that two binary trees on the same $n$ leaves
have subtrees with the same $\geq c\log\log n$ leaves which are homeomorphic,
such that homeomorphism is identity on the leaves.


  How intracellular and extracellular signals are integrated by transcription
factors is essential for understanding complex cellular patterns at the
population level. In this Letter, by using a synthetic genetic oscillator
coupled to a quorum-sensing apparatus, we propose an experimentally feasible
cis-regulatory module (CRM) which performs four possible logic operations
(ANDN, ORN, NOR and NAND) of input signals. We show both numerically and
theoretically that these different CRMs drive fundamentally different dynamic
patterns, such as synchronization, clustering and splay state.


  We show the existence of a high interrelation between the different loops
that may appear in a DNA segment. Conformational changes in a chain segment
caused by the formation of a particular loop may either promote or prevent the
appearance of another. The underlying loop selection mechanism is analyzed by
means of a Hamiltonian model from which the looping free energy and the
corresponding repression level can be computed. We show significant differences
between the probability of single and multiple loop formation. The consequences
that these collective effects might have on gene regulation processes are
outlined.


  We examine the problem of family size statistics (the number of individuals
carrying the same surname, or the same DNA sequence) in a given size subsample
of an exponentially growing population. We approach the problem from two
directions. In the first, we construct the family size distribution for the
subsample from the stable distribution for the full population. This latter
distribution is calculated for an arbitrary growth process in the limit of slow
growth, and is seen to depend only on the average and variance of the number of
children per individual, as well as the mutation rate. The distribution for the
subsample is shifted left with respect to the original distribution, tending to
eliminate the part of the original distribution reflecting the small families,
and thus increasing the mean family size. From the subsample distribution,
various bulk quantities such as the average family size and the percentage of
singleton families are calculated. In the second approach, we study the past
time development of these bulk quantities, deriving the statistics of the
genealogical tree of the subsample. This approach reproduces that of the first
when the current statistics of the subsample is considered. The surname
distribution from th e 2000 U.S. Census is examined in light of these findings,
and found to misrepresent the population growth rate by a factor of 1000.


  Developmental transcriptional networks in plants and animals operate in both
space and time. To understand these transcriptional networks it is essential to
obtain whole-genome expression data at high spatiotemporal resolution.
Substantial amounts of spatial and temporal microarray expression data
previously have been obtained for the Arabidopsis root; however, these two
dimensions of data have not been integrated thoroughly. Complicating this
integration is the fact that these data are heterogeneous and incomplete, with
observed expression levels representing complex spatial or temporal mixtures.
Given these partial observations, we present a novel method for reconstructing
integrated high resolution spatiotemporal data. Our method is based on a new
iterative algorithm for finding approximate roots to systems of bilinear
equations.


  This paper proposes an agent-based model which reproduces different
structures of animal groups. The shape and structure of the group is the effect
of simple interaction rules among individuals: each animal deploys itself
depending on the position of a limited number of close group mates. The
proposed model is shown to produce clustered formations, as well as lines and
V-like formations. The key factors which trigger the onset of different
patterns are argued to be the relative strength of attraction and repulsion
forces and, most important, the anisotropy in their application.


  This note represents the further progress in understanding the determination
of the genetic code by Golden mean (Rakocevic, 1998). Three classes of amino
acids that follow from this determination (the 7 "golden" amino acids, 7 of
their complements, and 6 non-complements) are observed now together with two
further possible splittings into 4 x 5 and 5 x 4 amino acids.


  Multipotent stem or progenitor cells undergo a sequential series of binary
fate decisions, which ultimately generate the diversity of differentiated
cells. Efforts to understand cell fate control have focused on simple gene
regulatory circuits that predict the presence of multiple stable states,
bifurcations and switch-like transitions. However, existing gene network models
do not explain more complex properties of cell fate dynamics such as the
hierarchical branching of developmental paths. Here, we construct a generic
minimal model of the genetic regulatory network controlling cell fate
determination, which exhibits five elementary characteristics of cell
differentiation: stability, directionality, branching, exclusivity, and
promiscuous expression. We argue that a modular architecture comprising
repeated network elements reproduces these features of differentiation by
sequentially repressing selected modules and hence restricting the dynamics to
lower dimensional subspaces of the high-dimensional state space. We implement
our model both with ordinary differential equations (ODEs), to explore the role
of bifurcations in producing the one-way character of differentiation, and with
stochastic differential equations (SDEs), to demonstrate the effect of noise on
the system. We further argue that binary cell fate decisions are prevalent in
cell differentiation due to general features of the underlying dynamical
system. This minimal model makes testable predictions about the structural
basis for directional, discrete and diversifying cell phenotype development and
thus can guide the evaluation of real gene regulatory networks that govern
differentiation.


  The gram-negative bacterium Yersinia pestis, the aetiological agent of
bubonic plague, is one the deadliest pathogens known to man. Despite its
historical reputation, plague is a modern disease which annually afflicts
thousands of people. Public safety considerations greatly limit clinical
experimentation on this organism and thus development of theoretical tools to
analyze the capabilities of this pathogen is of utmost importance. Here, we
report the first genome-scale metabolic model of Yersinia pestis biovar
Mediaevalis based both on its recently annotated genome, and physiological and
biochemical data from literature. Our model demonstrates excellent agreement
with Y. pestis known metabolic needs and capabilities. Since Y. pestis is a
meiotrophic organism, we have developed CryptFind, a systematic approach to
identify all candidate cryptic genes responsible for known and theoretical
meiotrophic phenomena. In addition to uncovering every known cryptic gene for
Y. pestis, our analysis of the rhamnose fermentation pathway suggests that betB
is the responsible cryptic gene. Despite all of our medical advances, we still
do not have a vaccine for bubonic plague. Recent discoveries of antibiotic
resistant strains of Yersinia pestis coupled with the threat of plague being
used as a bioterrorism weapon compel us to develop new tools for studying the
physiology of this deadly pathogen. Using our theoretical model, we can study
the cells phenotypic behavior under different circumstances and identify
metabolic weaknesses which may be harnessed for the development of
therapeutics. Additionally, the automatic identification of cryptic genes
expands the usage of genomic data for pharmaceutical purposes.


  Ratio-dependent predator-prey models have been increasingly favored by field
ecologists where predator-prey interactions have to be taken into account the
process of predation search. In this paper we study the conditions of the
existence and stability properties of the equilibrium solutions in a
reaction-diffusion model in which predator mortality is neither a constant nor
an unbounded function, but it is increasing with the predator abundance. We
show that analytically at a certain critical value a diffusion driven (Turing
type) instability occurs, i.e. the stationary solution stays stable with
respect to the kinetic system (the system without diffusion). We also show that
the stationary solution becomes unstable with respect to the system with
diffusion and that Turing bifurcation takes place: a spatially non-homogenous
(non-constant) solution (structure or pattern) arises. A numerical scheme that
preserve the positivity of the numerical solutions and the boundedness of prey
solution will be presented. Numerical examples are also included.


  In order to survive, reproduce and (in multicellular organisms)
differentiate, cells must control the concentrations of the myriad different
proteins that are encoded in the genome. The precision of this control is
limited by the inevitable randomness of individual molecular events. Here we
explore how cells can maximize their control power in the presence of these
physical limits; formally, we solve the theoretical problem of maximizing the
information transferred from inputs to outputs when the number of available
molecules is held fixed. We start with the simplest version of the problem, in
which a single transcription factor protein controls the readout of one or more
genes by binding to DNA. We further simplify by assuming that this regulatory
network operates in steady state, that the noise is small relative to the
available dynamic range, and that the target genes do not interact. Even in
this simple limit, we find a surprisingly rich set of optimal solutions.
Importantly, for each locally optimal regulatory network, all parameters are
determined once the physical constraints on the number of available molecules
are specified. Although we are solving an over--simplified version of the
problem facing real cells, we see parallels between the structure of these
optimal solutions and the behavior of actual genetic regulatory networks.
Subsequent papers will discuss more complete versions of the problem.


  Skeptics of biological evolution often claim that evolution requires a
decrease in entropy, giving rise to a conflict with the second law of
thermodynamics. This argument is fallacious because it neglects the large
increase in entropy provided by sunlight striking the Earth. A recent article
provided a quantitative assessment of the entropies involved and showed
explicitly that there is no conflict. That article rests on an unjustified
assumption about the amount of entropy reduction involved in evolution. I
present a refinement of the argument that does not rely on this assumption.


  Biological systems possess negative entropy. In them, one form of order
produces another, more organized form of order. We propose a formal scheme to
calculate robustness of an entire biological system by quantifying the negative
entropy present in it. Our Methodology is based upon a computational
implementation of two-person non-cooperative finite zero-sum game between
positive (physico-chemical) and negative (biological) entropy, present in the
system(TCA cycle, for this work). Biochemical analogue of Nash equilibrium,
proposed here, could measure the robustness in TCA cycle in exact numeric
terms, whereas the mixed strategy game between these entropies could quantitate
the progression of stages of biological adaptation. Synchronization profile
amongst macromolecular concentrations (even under environmental perturbations)
is found to account for negative entropy and biological robustness. Emergence
of synchronization profile was investigated with dynamically varying metabolite
concentrations. Obtained results were verified with that from the deterministic
simulation methods. Categorical plans to apply this algorithm in Cancer studies
and anti-viral therapies are proposed alongside. From theoretical perspective,
this work proposes a general, rigorous and alternative view of immunology.


  We develop a numerical method to simulate mechanical objects in a viscous
medium at a scale where inertia is negligible. Fibers, spheres and other
voluminous objects are represented with points. Different types of connections
are used to link the points together and in this way create composite
mechanical structures. The motion of such structures in a Brownian environment
is described by a first-order multivariate Langevin equation. We propose a
computationally efficient method to integrate the equation, and illustrate the
applicability of the method to cytoskeletal modeling with several examples.


  The mitotic spindle is an important intermediate structure in eukaryotic cell
division, in which each of a pair of duplicated chromosomes is attached through
microtubules to centrosomal bodies located close to the two poles of the
dividing cell. Several mechanisms are at work towards the formation of the
spindle, one of which is the `capture' of chromosome pairs, held together by
kinetochores, by randomly searching microtubules. Although the entire cell
cycle can be up to 24 hours long, the mitotic phase typically takes only less
than an hour. How does the cell keep the duration of mitosis within this limit?
Previous theoretical studies have suggested that the chromosome search and
capture is optimized by tuning the microtubule dynamic parameters to minimize
the search time. In this paper, we examine this conjecture. We compute the mean
search time for a single target by microtubules from a single nucleating site,
using a systematic and rigorous theoretical approach, for arbitrary kinetic
parameters. The result is extended to multiple targets and nucleating sites by
physical arguments. Estimates of mitotic time scales are then obtained for
different cells using experimental data. In yeast and mammalian cells, the
observed changes in microtubule kinetics between interphase and mitosis are
beneficial in reducing the search time. In {\it Xenopus} extracts, by contrast,
the opposite effect is observed, in agreement with the current understanding
that large cells use additional mechanisms to regulate the duration of the
mitotic phase.


  We report on numerical investigations concerning the propagation of solitons
in a real DNA chain (the Human Adenovirus 2) using a realistic model of DNA
torsional dynamics; this takes fully into account the inhomogeneities in the
real chain. We find that twist solitons propagate for considerable distances
(2-10 times their diameters) before stopping due to phonon emission. Our
results show that twist solitons may exist in real DNA chains; and on a more
general level that solitonic propagation can take place in highly inhomogeneous
media.


  Studying the topology of so-called {\em real networks}, that is networks
obtained from sociological or biological data for instance, has become a major
field of interest in the last decade. One way to deal with it is to consider
that networks are built from small functional units called {\em motifs}, which
can be found by looking for small subgraphs whose numbers of occurrences in the
whole network of interest are surprisingly high. In this paper, we propose to
define motifs through a local over-representation in the network and develop a
statistic which allows us to detect them limiting the number of false positives
and without time-consuming simulations. We apply it to the Yeast gene
interaction data and show that the known biologically relevant motifs are found
again and that our method gives some more information than the existing ones.


  Of the twenty amino acids used in proteins, ten were formed in Miller's
atmospheric discharge experiments. The two other major proposed sources of
prebiotic amino acid synthesis include formation in hydrothermal vents and
delivery to Earth via meteorites. We combine observational and experimental
data of amino acid frequencies formed by these diverse mechanisms and show
that, regardless of the source, these ten early amino acids can be ranked in
order of decreasing abundance in prebiotic contexts. This order can be
predicted by thermodynamics. The relative abundances of the early amino acids
were most likely reflected in the composition of the first proteins at the time
the genetic code originated. The remaining amino acids were incorporated into
proteins after pathways for their biochemical synthesis evolved. This is
consistent with theories of the evolution of the genetic code by stepwise
addition of new amino acids. These are hints that key aspects of early
biochemistry may be universal.


  The bacterial flagellar motor is a highly efficient rotary machine used by
many bacteria to propel themselves. It has recently been shown that at low
speeds its rotation proceeds in steps [Sowa et al. (2005) Nature 437,
916--919]. Here we propose a simple physical model that accounts for this
stepping behavior as a random walk in a tilted corrugated potential that
combines torque and contact forces. We argue that the absolute angular position
of the rotor is crucial for understanding step properties, and show this
hypothesis to be consistent with the available data, in particular the
observation that backward steps are smaller on average than forward steps. Our
model also predicts a sublinear torque-speed relationship at low torque, and a
peak in rotor diffusion as a function of torque.


  The effect of signal integration through cis-regulatory modules (CRMs) on
synchronization and clustering of populations of two-component genetic
oscillators coupled by quorum sensing is in detail investigated. We find that
the CRMs play an important role in achieving synchronization and clustering.
For this, we investigate 6 possible cis-regulatory input functions (CRIFs) with
AND, OR, ANDN, ORN, XOR, and EQU types of responses in two possible kinds of
cell-to-cell communications: activator-regulated communication (i.e., the
autoinducer regulates the activator) and repressor-regulated communication
(i.e., the autoinducer regulates the repressor). Both theoretical analysis and
numerical simulation show that different CRMs drive fundamentally different
cellular patterns, such as complete synchronization, various cluster-balanced
states and several cluster-nonbalanced states.


  Glycosylation is a highly complex process to produce a diverse repertoire of
cellular glycans that are attached to proteins and lipids. Glycans are involved
in fundamental biological processes, including protein folding and clearance,
cell proliferation and apoptosis, development, immune responses, and
pathogenesis. One of the major types of glycans, N-linked glycans, is formed by
sequential attachments of monosaccharides to proteins by a limited number of
enzymes. Many of these enzymes can accept multiple N-linked glycans as
substrates, thereby generating a large number of glycan intermediates and their
intermingled pathways. Motivated by the quantitative methods developed in
complex network research, we investigated the large-scale organization of such
N-linked glycosylation pathways in mammalian cells. The N-linked glycosylation
pathways are extremely modular, and are composed of cohesive topological
modules that directly branch from a common upstream pathway of glycan
synthesis. This unique structural property allows the glycan production between
modules to be controlled by the upstream region. Although the enzymes act on
multiple glycan substrates, indicating cross-talk between modules, the impact
of the cross-talk on the module-specific enhancement of glycan synthesis may be
confined within a moderate range by transcription-level control. The findings
of the present study provide experimentally-testable predictions for
glycosylation processes, and may be applicable to therapeutic glycoprotein
engineering.


  This unit describes strategies for studying the proteomes of organelles,
which is one example of targeted proteomics. It relies heavily on previously
published units dealing with organelle preparation, protein solubilization, and
proteomics techniques. A specific commentary for organelle proteomics is
provided. Specific protocols for the isolation of nuclei from various sources
(cell cultures, tissues) are also provided.


  In recent times whole-genome gene expression analysis has turned out to be a
highly important tool to study the coordinated function of a very large number
of genes within their corresponding cellular environment, especially in
relation to phenotypic diversity and disease. A wide variety of methods of
quantitative analysis have been developed to cope with high throughput data
sets generated by gene expression profiling experiments. Due to the complexity
associated with transcriptomics, specially in the case of gene regulation
phenomena, most of these methods are of a probabilistic or statistical nature.
Even if these methods have reached a central status in the development of an
integrative, systematic understanding of the associated biological processes,
they very rarely constitute a concrete guide to the actual physicochemical
mechanisms behind biological function and the role of these methods is more on
a hypotheses generating line. An important improvement could be done with the
development of a thermodynamic theory for gene expression and transcriptional
regulation that will build the foundations for a proper integration of the vast
amount of molecular biophysical data and could lead, in the future, to a
systemic view of genetic transcription and regulation.


  Because of the outstanding ability of two-dimensional electrophoresis to
separate complex mixtures of intact proteins, it would be advantageous to apply
it to all types of proteins, including hydrophobic and membrane proteins.
Unfortunately, poor solubility hampers the analysis of these molecules. As
these problems arise mainly in the extraction and isoelectric focusing steps,
the solution is to improve protein solubility under the conditions prevailing
during isoelectric focusing. This chapter describes the use of chaotropes and
novel detergents to enhance protein solubility during sample extraction and
isoelectric focussing, and discusses the contribution of these compounds to
improving proteomic analysis of membrane proteins.


  the MIAPE Gel Electrophoresis (MIAPE-GE) guidelines specify the minimum
information that should be provided when reporting the use of n-dimensional gel
electrophoresis in a proteomics experiment. Developed through a joint effort
between the gel-based analysis working group of the Human Proteome
Organisation's Proteomics Standards Initiative (HUPO-PSI;
http://www.psidev.info/) and the wider proteomics community, they constitute
one part of the overall Minimum Information about a Proteomics Experiment
(MIAPE) documentation system published last August in Nature Biotechnology


  In the last ten years, the field of proteomics has expanded at a rapid rate.
A range of exciting new technology has been developed and enthusiastically
applied to an enormous variety of biological questions. However, the degree of
stringency required in proteomic data generation and analysis appears to have
been underestimated. As a result, there are likely to be numerous published
findings that are of questionable quality, requiring further confirmation
and/or validation. This manuscript outlines a number of key issues in proteomic
research, including those associated with experimental design, differential
display and biomarker discovery, protein identification and analytical
incompleteness. In an effort to set a standard that reflects current thinking
on the necessary and desirable characteristics of publishable manuscripts in
the field, a minimal set of guidelines for proteomics research is then
described. These guidelines will serve as a set of criteria which editors of
PROTEOMICS will use for assessment of future submissions to the Journal.


  Biochemical reaction networks in living cells usually involve reversible
covalent modification of signaling molecules, such as protein phosphorylation.
Under conditions of small molecule numbers, as is frequently the case in living
cells, mass action theory fails to describe the dynamics of such systems.
Instead, the biochemical reactions must be treated as stochastic processes that
intrinsically generate concentration fluctuations of the chemicals. We
investigate the stochastic reaction kinetics of covalent modification cycles
(CMCs) by analytical modeling and numerically exact Monte-Carlo simulation of
the temporally fluctuating concentration. Depending on the parameter regime, we
find for the probability density of the concentration qualitatively distinct
classes of distribution functions, including power law distributions with a
fractional and tunable exponent. These findings challenge the traditional view
of biochemical control networks as deterministic computational systems and
suggest that CMCs in cells can function as versatile and tunable noise
generators.


  Background: If biology is modular then clusters, or communities, of proteins
derived using only protein interaction network structure should define protein
modules with similar biological roles. We investigate the link between
biological modules and network communities in yeast and its relationship to the
scale at which we probe the network.
  Results: Our results demonstrate that the functional homogeneity of
communities depends on the scale selected, and that almost all proteins lie in
a functionally homogeneous community at some scale. We judge functional
homogeneity using a novel test and three independent characterizations of
protein function, and find a high degree of overlap between these measures. We
show that a high mean clustering coefficient of a community can be used to
identify those that are functionally homogeneous. By tracing the community
membership of a protein through multiple scales we demonstrate how our approach
could be useful to biologists focusing on a particular protein.
  Conclusions: We show that there is no one scale of interest in the community
structure of the yeast protein interaction network, but we can identify the
range of resolution parameters that yield the most functionally coherent
communities, and predict which communities are most likely to be functionally
homogeneous.


  Chloroplast microsatellites have been widely used in population genetic
studies of conifers in recent years. However, their haplotype configurations
suggest that they could have high levels of homoplasy, thus limiting the power
of these molecular markers. A coalescent-based computer simulation was used to
explore the influence of homoplasy on measures of genetic diversity based on
chloroplast microsatellites. The conditions of the simulation were defined to
fit isolated populations originating from the colonization of one single
haplotype into an area left available after a glacial retreat. Simulated data
were compared with empirical data available from the literature for a species
of Pinus that has expanded north after the Last Glacial Maximum. In the
evaluation of genetic diversity, homoplasy was found to have little influence
on Nei's unbiased haplotype diversity (H(E)) while Goldstein's genetic distance
estimates (D2sh) were much more affected. The effect of the number of
chloroplast microsatellite loci for evaluation of genetic diversity is also
discussed.


  It is presented that the positions of amino acids within Genetic Code Table
follow from strict their physical and chemical properties as well as from a
pure formal determination by the Golden mean.


  An analysis of a variety of existing experimental data leads to the
conclusion on the existence of a resonance mechanism allowing weak magnetic
fields to affect biological processes. These fields may either be static
magnetic fields comparable in magnitude with the magnetic field of the earth or
weak ultra-low frequency time-dependent fields. So far, a generally accepted
theoretical model allowing one to understand the effect of magnetic and
electric fields on biological processes is not available. By this reason, it is
not clear which characteristics of the fields, like magnetic and electric field
strength, frequency of change of the field, shape of the electromagnetic wave,
the duration of the magnetic or electric influence or some particular
combination of them, are responsible for the biological effect. In the present
analysis it is shown that external time-independent magnetic fields may cause a
resonance amplification of ionic electric currents in biological tissues and,
in particular, in the vasculature system due to a Brownian motion of charges.
These resonance electric currents may cause necrotic changes in the tissues or
blood circulation and in this way significantly affect the biological organism.
The magnitude of the magnetic fields leading to resonance effects is estimated,
it is shown that it depends significantly on the radius of the blood
capillaries.


  In order to empirically assess the effect of actual seed sampling strategy on
genetic diversity of holm oak (Quercus ilex) forestations in Sicily, we have
analysed the genetic composition of two seedling lots (nursery stock and
plantation) and their known natural seed origin stand by means of six nuclear
microsatellite loci. Significant reduction in genetic diversity and significant
difference in genetic composition of the seedling lots compared to the seed
origin stand were detected. The female and the total effective number of
parents were quantified by means of maternity assignment of seedlings and
temporal changes in allele frequencies. Extremely low effective maternity
numbers were estimated (Nfe $\approx$ 2-4) and estimates accounting for both
seed and pollen donors gave also low values (Ne $\approx$ 35-50). These values
can be explained by an inappropriate forestry seed harvest strategy limited to
a small number of spatially close trees.


  "Module networks" are a framework to learn gene regulatory networks from
expression data using a probabilistic model in which coregulated genes share
the same parameters and conditional distributions. We present a method to infer
ensembles of such networks and an averaging procedure to extract the
statistically most significant modules and their regulators. We show that the
inferred probabilistic models extend beyond the data set used to learn the
models.


  Cellular membranes are a heterogeneous mix of lipids, proteins and small
molecules. Special groupings of saturated lipids and cholesterol form a
liquid-ordered phase, known as `lipid rafts,' serving as platforms for
signaling, trafficking and material transport throughout the secretory pathway.
Questions remain as to how the cell maintains heterogeneity of a fluid membrane
with multiple phases, through time, on a length-scale consistent with the fact
that no large-scale phase separation is observed. We have utilized a
combination of mechanical modeling and in vitro experiments to show that
membrane morphology can be a key player in maintaining this heterogeneity and
organizing such domains in the membrane. We demonstrate that lipid domains can
adopt a flat or dimpled morphology, where the latter facilitates a repulsive
interaction that slows coalescence and tends to organize domains. These forces,
that depend on domain morphology, play an important role in regulating lipid
domain size and in the lateral organization of lipids in the membrane.


  We describe and test an implicit solvent all-atom potential for simulations
of protein folding and aggregation. The potential is developed through studies
of structural and thermodynamic properties of 17 peptides with diverse
secondary structure. Results obtained using the final form of the potential are
presented for all these peptides. The same model, with unchanged parameters, is
furthermore applied to a heterodimeric coiled-coil system, a mixed alpha/beta
protein and a three-helix-bundle protein, with very good results. The
computational efficiency of the potential makes it possible to investigate the
free-energy landscape of these 49--67-residue systems with high statistical
accuracy, using only modest computational resources by today's standards.


  We investigate the structural and dynamical properties of the transcriptional
regulatory network of the yeast {\it Saccharomyces cerevisiae} and compare it
with two unbiased ensembles: one obtained by reshuffling the edges and the
other generated by mimicking the transcriptional regulation mechanism within
the cell. Both ensembles reproduce the degree distributions (the first -by
construction- exactly and the second approximately), degree-degree correlations
and the $k$-core structure observed in Yeast. An exceptionally large
dynamically relevant core network found in Yeast in comparison with the second
ensemble points to a strong bias towards a collective organization which is
achieved by subtle modifications in the network's degree distributions. We use
a Boolean model of regulatory dynamics with various classes of update functions
to represent in vivo regulatory interactions. We find that the Yeast's core
network has a qualitatively different behaviour, accommodating on average
multiple attractors unlike typical members of both reference ensembles which
converge to a single dominant attractor. Finally, we investigate the robustness
of the networks and find that the stability depends strongly on the used
function class. The robustness measure is squeezed into a narrower band around
the order-chaos boundary when Boolean inputs are required to be nonredundant on
each node. However, the difference between the reference models and the Yeast's
core is marginal, suggesting that the dynamically stable network elements are
located mostly on the peripherals of the regulatory network. Consistently, the
statistically significant three-node motifs in the dynamical core of Yeast turn
out to be different from and less stable than those found in the full
transcriptional regulatory network.


  Biochemical processes typically involve huge numbers of individual reversible
steps, each with its own dynamical rate constants. For example, kinetic
proofreading processes rely upon numerous sequential reactions in order to
guarantee the precise construction of specific macromolecules. In this work, we
study the transient properties of such systems and fully characterize their
first passage (completion) time distributions. In particular, we provide
explicit expressions for the mean and the variance of the completion time for a
kinetic proofreading process and computational analyses for more complicated
biochemical systems. We find that, for a wide range of parameters, as the
system size grows, the completion time behavior simplifies: it becomes either
deterministic or exponentially distributed, with a very narrow transition
between the two regimes. In both regimes, the dynamical complexity of the full
system is trivial compared to its apparent structural complexity. Similar
simplicity is likely to arise in the dynamics of many complex multi-step
biochemical processes. In particular, these findings suggest not only that one
may not be able to understand individual elementary reactions from macroscopic
observations, but also that such understanding may be unnecessary.


  Morphogen profiles play a vital role in biology by specifying position in
embryonic development. However, the factors that influence the shape of a
morphogen profile remain poorly understood. Since morphogens should provide
precise positional information, one significant factor is the robustness of the
profile to noise. We compare three classes of morphogen profiles (linear,
exponential, algebraic) to see which is most precise when subject to both
external embryo-to-embryo fluctuations and internal fluctuations due to
intrinsically random processes such as diffusion. We find that both the kinetic
parameters and the overall gradient shape (e.g. exponential versus algebraic)
can be optimised to generate maximally precise positional information.


  Population extinction is a rare event which requires overcoming an effective
barrier. We show that the extinction rate can be fragile: a small change in the
system parameters leads to an exponentially strong change of the rate, with the
barrier height depending on the parameters nonanalytically. General conditions
of the fragility are established. The fragility is found in one of the
best-known models of epidemiology, the SIS model. The analytical expressions
are compared with simulations.


  Does bearing children shorten a woman's life expectancy? Several demographic
studies, historic and current, have found no such effect. But the Caerphilly
cohort study is far the most prominent and frequently-cited, and it answers in
the affirmative. Why has this study found an effect that others fail to see?
Their analysis is based on Poisson regression, a statistical technique that is
accurate only if the underlying data are Poisson distributed. But the
distribution of the number of children born to women in the Caerphilly database
departs strongly from Poisson at the high end. This makes the result overly
sensitive to a handful of women with 15 children or more who lived before 1700.
When these 5 women are removed from a database of more than 2,900, the Poisson
regression no longer shows a significant result. Bi-linear regression relating
life span to fertility and date of birth results in a positive coefficient for
fertility.


  Although DNA is often bent in vivo, it is unclear how DNA-bending forces
modulate DNA-protein binding affinity. Here, we report how a range of
DNA-bending forces modulates the binding of the Integration Host Factor (IHF)
protein to various DNAs. Using solution fluorimetry and electrophoretic
mobility shift assays, we measured the affinity of IHF for DNAs with different
bending forces and sequence mutations. Bending force was adjusted by varying
the fraction of double-stranded DNA in a circular substrate, or by changing the
overall size of the circle (1). DNA constructs contained a pair of Forster
Resonance Energy Transfer dyes that served as probes for affinity assays, and
read out bending forces measured by optical force sensors (2). Small bending
forces significantly increased binding affinity; this effect saturated beyond
~3 pN. Surprisingly, when DNA sequences that bound IHF only weakly were
mechanically bent by circularization, they bound IHF more tightly than the
linear "high-affinity" binding sequence. These findings demonstrate that small
bending forces can greatly augment binding at sites that deviate from a
protein's consensus binding sequence. Since cellular DNA is subject to
mechanical deformation and condensation, affinities of architectural proteins
determined in vitro using short linear DNAs may not reflect in vivo affinities.


  Response time-delay is an ubiquitous phenomenon in biological systems. Here
we use a simple stochastic population model with time-delayed switching-rate
conversion to quantitatively study the biological influence of the response
time-delay on the survival fitness of cells living in periodically fluctuating
and stochastic environment, respectively. Our calculation and simulation show
that for the cells having a slow rate transition into fit phenotype and a fast
rate transition into unfit phenotype, the response time-delay can always
enhance their fitness during the environment change. Particularly, in the
periodic or stochastic environment with small variance, the optimal fitness
achieved by these cells is superior to that of cells with reverse switching
rates even if the latter exhibits rapid response. These results suggest that
the response time delay may be utilized by cells to enhance their adaptation to
the fluctuating environment.


  Subspace clustering has gained increasing popularity in the analysis of gene
expression data. Among subspace cluster models, the recently introduced
order-preserving sub-matrix (OPSM) has demonstrated high promise. An OPSM,
essentially a pattern-based subspace cluster, is a subset of rows and columns
in a data matrix for which all the rows induce the same linear ordering of
columns. Existing OPSM discovery methods do not scale well to increasingly
large expression datasets. In particular, twig clusters having few genes and
many experiments incur explosive computational costs and are completely pruned
off by existing methods. However, it is of particular interest to determine
small groups of genes that are tightly coregulated across many conditions. In
this paper, we present KiWi, an OPSM subspace clustering algorithm that is
scalable to massive datasets, capable of discovering twig clusters and
identifying negative as well as positive correlations. We extensively validate
KiWi using relevant biological datasets and show that KiWi correctly assigns
redundant probes to the same cluster, groups experiments with common clinical
annotations, differentiates real promoter sequences from negative control
sequences, and shows good association with cis-regulatory motif predictions.


  Meso- and submesoscales (fronts, eddies, filaments) in surface ocean flow
have a crucial influence on marine ecosystems. Their dynamics partly control
the foraging behaviour and the displacement of marine top predators (tuna,
birds, turtles, and cetaceans). In this work we focus on the role of
submesoscale structures in the Mozambique Channel on the distribution of a
marine predator, the Great Frigatebird. Using a newly developed dynamical
concept, namely the Finite-Size Lyapunov Exponent (FSLE), we have identified
Lagrangian coherent structures (LCSs) present in the surface flow in the
Channel over a 2-month observation period (August and September 2003). By
comparing seabirds' satellite positions with LCSs locations, we demonstrate
that frigatebirds track precisely these structures in the Mozambique Channel,
providing the first evidence that a top predator is able to track these FSLE
ridges to locate food patches. After comparing bird positions during long and
short trips, and different parts of these trips, we propose several hypotheses
to understand how frigatebirds can follow these LCSs. The birds might use
visual and/or olfactory cues and/or atmospheric current changes over the
structures to move along these biological corridors. The birds being often
associated to tuna schools around foraging areas, a thorough comprehension of
their foraging behaviour and movement during the breeding season is crucial not
only to seabirds' ecology but also to an appropriate ecosystemic approach of
fisheries in the Channel.


  The first passage probability (FPP), of trafficked intracellular particles
reaching a displacement L, in a given time t or inverse velocity S = t/L, can
be calculated robustly from measured particle tracks, and gives a measure of
particle movement in which different types of motion, e.g. diffusion, ballistic
motion, and transient run-rest motion, can readily be distinguished in a single
graph, and compared with mathematical models. The FPP is attractive in that it
offers a means of reducing the data in the measured tracks, without making
assumptions about the mechanism of motion: for example, it does not employ
smoothing, segementation or arbitrary thresholds to discriminate between
different types of motion in a particle track. Taking experimental data from
tracked endocytic vesicles, and calculating the FPP, we see how three molecular
treatments affect the trafficking. We show the FPP can quantify complicated
movement which is neither completely random nor completely deterministic,
making it highly applicable to trafficked particles in cell biology.


  We propose a stochastic model for intracellular transport processes
associated with the activity of molecular motors. This out-of-equilibrium
model, based on a generalized Langevin equation, considers a particle immersed
in a viscoelastic environment and simultaneously driven by an external random
force that models the motors activity. An analytical expression for the mean
square displacement is derived, which exhibits a subdiffusive to superdiffusive
transition. We show that the experimentally accessible statistical properties
of the diffusive particle motion can be reproduced by this model.


  Intracellular transport of large cargoes, such as organelles, vesicles or
large proteins, is a complex dynamical process that involves the interplay of
ATP-consuming molecular motors, cytoskeleton filaments and the viscoelastic
cytoplasm. The displacements of particles or probes in the cell cytoplasm as a
function of time are characterized by different (anomalous) diffusion regimes.
We investigate here the motion of pigment organelles (melanosomes) driven by
myosin-V motors in \emph{Xenopus laevis} melanocytes using a high
spatio-temporal resolution tracking technique. By analyzing the mean square
displacement (MSD) of the obtained trajectories as a function of the time lag,
we show that the melanosomes display a transition between subdiffusive to
superdiffusive behavior. A stochastic theoretical model is introduced to
generalize the interpretation of our data. Starting from a generalized Langevin
equation that explicitly considers the collective action of the molecular
motors we derive an analytical expression for the MSD as a function of the time
lag, which also takes into account the experimental noise. By fitting our model
to the experimental data we were able to discriminate the exponents that
characterize the passive and active contributions to melanosome dynamics. The
model also estimates the "global" motor forces correctly. In this sense, our
model gives a quantitative description of active transport in living cells with
a reduced number of parameters.


  Protein-surface interactions cause the desirable effect of controlled protein
adsorption onto biodevices as well as the undesirable effect of protein
fouling. The key to controlling protein-surface adsorptions is to identify and
quantify the main adsorption mechanisms: adsorptions that occur (1) while
depositing a protein solution onto dry surfaces and (2) after the deposition
onto wet surfaces. Bulk measurements cannot reveal the dynamic protein
adsorption pathways and thus cannot differentiate between the two adsorption
mechanisms. We imaged the interactions of single streptavidin molecules with
hydrophobic fused-silica surfaces in real-time. We observed both adsorbed
proteins on surfaces and diffusing proteins near surfaces and analyzed their
adsorption kinetics. Our analysis shows that the protein solution deposition
process is the primary mechanism of streptavidin adsorption onto surfaces at
the sub-nanomolar to nanomolar protein concentrations. Furthermore, we found
that hydrophilic fused-silica surfaces can prevent the adsorption of
streptavidin molecules.


  Deterministic and stochastic Boolean network models are build for the
dynamics of negative feedback loops of the p53 pathways. It is shown that the
main function of the negative feedback in the p53 pathways is to keep p53 at a
low steady state level, and each sequence of protein states in the negative
feedback loops, is globally attracted to a closed cycle of the p53 dynamics
after being perturbed by outside signal (e.g. DNA damage). Our theoretical and
numerical studies show that both the biological stationary state and the
biological oscillation after being perturbed are stable for a wide range of
noise level. Applying the mathematical circulation theory of Markov chains, we
investigate their stochastic synchronized dynamics and by comparing the network
dynamics of the stochastic model with its corresponding deterministic network
counterpart, a dominant circulation in the stochastic model is the natural
generalization of the deterministic limit cycle in the deterministic system.
Moreover, the period of the main peak in the power spectrum, which is in common
use to characterize the synchronized dynamics, perfectly corresponds to the
number of states in the main cycle with dominant circulation. Such a large
separation in the magnitude of the circulations, between a dominant, main cycle
and the rest, gives rise to the stochastic synchronization phenomenon.


  Applying the mathematical circulation theory of Markov chains, we investigate
the synchronized stochastic dynamics of a discrete network model of yeast
cell-cycle regulation where stochasticity has been kept rather than being
averaged out. By comparing the network dynamics of the stochastic model with
its corresponding deterministic network counterpart, we show that the
synchronized dynamics can be soundly characterized by a dominant circulation in
the stochastic model, which is the natural generalization of the deterministic
limit cycle in the deterministic system. Moreover, the period of the main peak
in the power spectrum, which is in common use to characterize the synchronized
dynamics, perfectly corresponds to the number of states in the main cycle with
dominant circulation. Such a large separation in the magnitude of the
circulations, between a dominant, main cycle and the rest, gives rise to the
stochastic synchronization phenomenon.


  Galled trees, evolutionary networks with isolated reticulation cycles, have
appeared under several slightly different definitions in the literature. In
this paper we establish the actual relationships between the main four such
alternative definitions: namely, the original galled trees, level-1 networks,
nested networks with nesting depth 1, and evolutionary networks with
arc-disjoint reticulation cycles.


  In this paper, we extend a previously developed recursive entropic
segmentation scheme for applications to biological sequences. Instead of
Bernoulli chains, we model the statistically stationary segments in a
biological sequence as Markov chains, and define a generalized Jensen-Shannon
divergence for distinguishing between two Markov chains. We then undertake a
mean-field analysis, based on which we identify pitfalls associated with the
recursive Jensen-Shannon segmentation scheme. Following this, we explain the
need for segmentation optimization, and describe two local optimization schemes
for improving the positions of domain walls discovered at each recursion stage.
We also develop a new termination criterion for recursive Jensen-Shannon
segmentation based on the strength of statistical fluctuations up to a minimum
statistically reliable segment length, avoiding the need for unrealistic null
and alternative segment models of the target sequence. Finally, we compare the
extended scheme against the original scheme by recursively segmenting the
Escherichia coli K-12 MG1655 genome.


  Experiments show that for many two state folders the free energy of the
native state DG_ND([C]) changes linearly as the denaturant concentration [C] is
varied. The slope, m = d DG_ND([C])/d[C], is nearly constant. The m-value is
associated with the difference in the surface area between the native (N) and
the denatured (D) state, which should be a function of DR_g^2, the difference
in the square of the radius of gyration between the D and N states. Single
molecule experiments show that the denatured state undergoes an equilibrium
collapse transition as [C] decreases, which implies m also should be
[C]-dependent. We resolve the conundrum between constant m-values and
[C]-dependent changes in Rg using molecular simulations of a coarse-grained
representation of protein L, and the Molecular Transfer Model, for which the
equilibrium folding can be accurately calculated as a function of denaturant
concentration. We find that over a large range of denaturant concentration (> 3
M) the m-value is a constant, whereas under strongly renaturing conditions (< 3
M) it depends on [C]. The m-value is a constant above [C]> 3 M because the
[C]-dependent changes in the surface area of the backbone groups, which make
the largest contribution to m, is relatively narrow in the denatured state. The
burial of the backbone gives rise to substantial surface area changes below
[C]< 3 M, leading to collapse in the denatured state. The midpoint of
transition of individual residues vary significantly even though global folding
can be described as an all-or-none transition. Collapse driven by the loss of
favorable residue-solvent interactions and a concomitant increase in the
strength of intrapeptide interactions with decreasing [C]. These interactions
are non-uniformly distributed throughout the native structure of protein L.


  Protein amyloid fibrils are a form of linear protein aggregates that are
implicated in many neurodegenerative diseases. Here, we study the dynamics of
amyloid fibril elongation by performing Langevin dynamic simulations on a
coarse-grained model of peptides. Our simulation results suggest that the
elongation process is dominated by a series of local minimum due to frustration
in monomer-fibril interactions. This rugged energy landscape picture indicates
that the amount of recycling of monomers at the fibrils' ends before being
fibrilized is substantially reduced in comparison to the conventional two-step
elongation model. This picture, along with other predictions discussed, can be
tested with current experimental techniques.


  Polymerization of RNA from a template DNA is carried out by a molecular
machine called RNA polymerase (RNAP). It also uses the template as a track on
which it moves as a motor utilizing chemical energy input. The time it spends
at each successive monomer of DNA is random; we derive the exact distribution
of these "dwell times" in our model. The inverse of the mean dwell time
satisfies a Michaelis-Menten-like equation and is also consistent with a
general formula derived earlier by Fisher and Kolomeisky for molecular motors
with unbranched mechano-chemical cycles. Often many RNAP motors move
simultaneously on the same track. Incorporating the steric interactions among
the RNAPs in our model, we also plot the three-dimensional phase diagram of our
model for RNAP traffic using an extremum current hypothesis.


  A statistical model of protein families, called profile conditional random
fields (CRFs), is proposed. This model may be regarded as an integration of the
profile hidden Markov model (HMM) and the Finkelstein-Reva (FR) theory of
protein folding. While the model structure of the profile CRF is almost
identical to the profile HMM, it can incorporate arbitrary correlations in the
sequences to be aligned to the model. In addition, like in the FR theory, the
profile CRF can incorporate long-range pairwise interactions between model
states via mean-field-like approximations. We give the detailed formulation of
the model, self-consistent approximations for treating long-range interactions,
and algorithms for computing partition functions and marginal probabilities. We
also outline the methods for the global optimization of model parameters as
well as a Bayesian framework for parameter learning and selection of optimal
alignments.


  In this paper, we describe the context sensitivity problem encountered in
partitioning a heterogeneous biological sequence into statistically homogeneous
segments. After showing signatures of the problem in the bacterial genomes of
Escherichia coli K-12 MG1655 and Pseudomonas syringae DC3000, when these are
segmented using two entropic segmentation schemes, we clarify the contextual
origins of these signatures through mean-field analyses of the segmentation
schemes. Finally, we explain why we believe all sequence segmentation schems
are plagued by the context sensitivity problem.


  Mathew-Fenn et al. (Science (2008) 322, 446-9) measured end-to-end distances
of short DNA and concluded that stretching fluctuations in several consecutive
turns of the double helix should be strongly correlated. I argue that this
conclusion is based on incorrect assumptions, notably, on a simplistic
treatment of the excluded volume effect of reporter labels. Contrary to the
author's claim, their conclusion is not supported by other data.


  Summary: Founded upon diffusion with damping, ITM Probe is an application for
modeling information flow in protein interaction networks without prior
restriction to the sub-network of interest. Given a context consisting of
desired origins and destinations of information, ITM Probe returns the set of
most relevant proteins with weights and a graphical representation of the
corresponding sub-network. With a click, the user may send the resulting
protein list for enrichment analysis to facilitate hypothesis formation or
confirmation.
  Availability: ITM Probe web service and documentation can be found at
www.ncbi.nlm.nih.gov/CBBresearch/qmbp/mn/itm_probe


  In this paper, we investigate how stochastic reaction processes are affected
by external perturbations. We describe an extension of the deterministic
metabolic control analysis (MCA) to the stochastic regime. We introduce
stochastic sensitivities for mean and covariance values of reactant
concentrations and reaction fluxes and show that there exist MCA-like summation
theorems among these sensitivities. The summation theorems for flux variances
are shown to depend on the size of the measurement time window ($\epsilon$),
within which reaction events are counted for measuring a single flux. The
degree of the $\epsilon$-dependency can become significant for processes
involving multi-time-scale dynamics and is estimated by introducing a new
measure of time scale separation. This $\epsilon$-dependency is shown to be
closely related to the power-law scaling observed in flux fluctuations in
various complex networks. We propose a systematic way to control fluctuations
of reactant concentrations while minimizing changes in mean concentration
levels. Such orthogonal control is obtained by introducing a control vector
indicating the strength and direction of parameter perturbations leading to a
sensitive control. We also propose a possible implication in the control of
flux fluctuation: The control distribution for flux fluctuations changes with
the measurement time window size, $\epsilon$. When a control engineer applies a
specific control operation on a reaction system, the system can respond
contrary to what is expected, depending on the time window size $\epsilon$.


  We present two examples of how single-molecule experimental techniques
applied to biological systems can give insight into problems within the scope
of equilibrium and nonequilibrium mesoscopic thermodynamics. The first example
is the mapping of the free energy landscape of a macromolecule, the second the
experimental verification of Crooks' fluctuation theorem. In both cases the
experimental setup comprises optical tweezers and DNA molecules.


  The folding dynamics of small single-domain proteins is a current focus of
simulations and experiments. Many of these proteins are 'two-state folders',
i.e. proteins that fold rather directly from the denatured state to the native
state, without populating metastable intermediate states. A central question is
how to characterize the instable, partially folded conformations of two-state
proteins, in particular the rate-limiting transition-state conformations
between the denatured and the native state. These partially folded
conformations are short-lived and cannot be observed directly in experiments.
However, experimental data from detailed mutational analyses of the folding
dynamics provide indirect access to transition states. The interpretation of
these data, in particular the reconstruction of transition-state conformations,
requires simulation and modeling. The traditional interpretation of the
mutational data aims to reconstruct the degree of structure formation of
individual residues in the transition state, while a novel interpretation aims
at degrees of structure formation of cooperative substructures such as
alpha-helices and beta-hairpins. By splitting up mutation-induced free energy
changes into secondary and tertiary structural components, the novel
interpretation resolves some of the inconsistencies of the traditional
interpretation.


  In Sphingomonas CHY-1, a single ring-hydroxylating dioxygenase is responsible
for the initial attack of a range of polycyclic aromatic hydrocarbons (PAHs)
composed of up to five rings. The components of this enzyme were separately
purified and characterized. The oxygenase component (ht-PhnI) was shown to
contain one Rieske-type [2Fe-2S] cluster and one mononuclear Fe center per
alpha subunit, based on EPR measurements and iron assay. Steady-state kinetic
measurements revealed that the enzyme had a relatively low apparent Michaelis
constant for naphthalene (Km= 0.92 $\pm$ 0.15 $\mu$M), and an apparent
specificity constant of 2.0 $\pm$ 0.3 $\mu$M-1 s-1. Naphthalene was converted
to the corresponding 1,2-dihydrodiol with stoichiometric oxidation of NADH. On
the other hand, the oxidation of eight other PAHs occurred at slower rates, and
with coupling efficiencies that decreased with the enzyme reaction rate.
Uncoupling was associated with hydrogen peroxide formation, which is
potentially deleterious to cells and might inhibit PAH degradation. In single
turnover reactions, ht-PhnI alone catalyzed PAH hydroxylation at a faster rate
in the presence of organic solvent, suggesting that the transfer of substrate
to the active site is a limiting factor. The four-ring PAHs chrysene and
benz[a]anthracene were subjected to a double ring-dihydroxylation, giving rise
to the formation of a significant proportion of bis-cis-dihydrodiols. In
addition, the dihydroxylation of benz[a]anthracene yielded three dihydrodiols,
the enzyme showing a preference for carbons in positions 1,2 and 10,11. This is
the first characterization of a dioxygenase able to dihydroxylate PAHs made up
of four and five rings.


  We develop a model for speciation due to postzygotic incompatibility
generated by autoimmune reactions. The model is based on predator-prey
interactions between a host plants and their pathogens. Such interactions are
often frequency-dependent, so that pathogen attack is focused on the most
abundant plant phenotype, while rare plant types may escape pathogen attack.
Thus, frequency dependence can generate disruptive selection, which can give
rise to speciation if distant phenotypes become reproductively isolated. Based
on recent experimental evidence from {\it Arabidopsis}, we assume that at the
molecular level, incompatibility between strains is caused by epistatic
interactions between two proteins in the plant immune system, the guard and the
guardee. Within each plant strain, immune reactions occur when the guardee
protein is modified by a pathogen effector, and the guard subsequently binds to
the guardee, thus precipitating an immune response. However, when guard and
guardee proteins come from phenotypically distant parents, a hybrid's immune
system can be triggered by erroneous interactions between these proteins even
in the absence of pathogen attack, leading to severe autoimmune reactions in
hybrids. Our model shows how phenotypic variation generated by
frequency-dependent host-pathogen interactions can lead to postzygotic
incompatibility between extremal types, and hence to speciation.


  Protein solubilization for two-dimensional electrophoresis (2DE) has to break
molecular interactions to separate the biological contents of the material of
interest into isolated and intact polypeptides. This must be carried out in
conditions compatible with the first dimension of 2DE, namely isoelectric
focusing. In addition, the extraction process must enable easy removal of any
nonprotein component interfering with the isoelectric focusing. The constraints
brought in this process by the peculiar features of isoelectric focusing are
discussed, as well as their consequences in terms of possible solutions and
limits for the solubilization process.


  Silver staining detects proteins after electrophoretic separation on
polyacrylamide gels. Its main positive features are its excellent sensitivity
(in the low nanogram range) and the use of very simple and cheap equipment and
chemicals. The sequential phases of silver staining are protein fixation, then
sensitization, then silver impregnation, and finally image development. Several
variants of silver staining are described here, which can be completed in a
time range from 2 h to 1 day after the end of the electrophoretic separation.
Once completed, the stain is stable for several weeks.


  Circular double stranded DNA has different topological states which are
defined by their linking numbers. Equilibrium distribution of linking numbers
can be obtained by closing a linear DNA into a circle by ligase. Using Monte
Carlo simulation, we predict the temperature dependence of the linking number
distribution of small circular DNAs. Our predictions are based on flexible
defect excitations resulted from local melting or unstacking of DNA base pairs.
We found that the reduced bending rigidity alone can lead to measurable changes
of the variance of linking number distribution of short circular DNAs. If the
defect is accompanied by local unwinding, the effect becomes much more
prominent. The predictions can be easily investigated in experiments, providing
a new method to study the micromechanics of sharply bent DNAs and the thermal
stability of specific DNA sequences. Furthermore, the predictions are directly
applicable to the studies of binding of DNA distorting proteins that can
locally reduce DNA rigidity, form DNA kinks, or introduce local unwinding.


  Initial reactions involved in the bacterial degradation of polycyclic
aromatic hydrocarbons (PAHs) include a ring-dihydroxylation catalyzed by a
dioxygenase and a subsequent oxidation of the dihydrodiol products by a
dehydrogenase. In this study, the dihydrodiol dehydrogenase from the
PAH-degrading Sphingomonas strain CHY-1 has been characterized. The bphB gene
encoding PAH dihydrodiol dehydrogenase (PDDH) was cloned and overexpressed as a
His-tagged protein. The recombinant protein was purified as a homotetramer with
an apparent Mr of 110,000. PDDH oxidized the cis-dihydrodiols derived from
biphenyl and eight polycyclic hydrocarbons, including chrysene,
benz[a]anthracene, and benzo[a]pyene, to corresponding catechols. Remarkably,
the enzyme oxidized pyrene 4,5-dihydrodiol, whereas pyrene is not metabolized
by strain CHY-1. The PAH catechols produced by PDDH rapidly auto-oxidized in
air but were regenerated upon reaction of the o-quinones formed with NADH.
Kinetic analyses performed under anoxic conditions revealed that the enzyme
efficiently utilized two- to four-ring dihydrodiols, with Km values in the
range of 1.4 to 7.1 $\mu$M, and exhibited a much higher Michaelis constant for
NAD+ (Km of 160 $\mu$M). At pH 7.0, the specificity constant ranged from (1.3
$\pm$ 0.1) x 106 M?1 s?1 with benz[a]anthracene 1,2-dihydrodiol to (20.0 $\pm$
0.8) x 106 M?1 s?1 with naphthalene 1,2-dihydrodiol. The catalytic activity of
the enzyme was 13-fold higher at pH 9.5. PDDH was subjected to inhibition by
NADH and by 3,4-dihydroxyphenanthrene, and the inhibition patterns suggested
that the mechanism of the reaction was ordered Bi Bi. The regulation of PDDH
activity appears as a means to prevent the accumulation of PAH catechols in
bacterial cells.


  An individual-based model (IBM) of a spatiotemporal terrestrial ecological
population is proposed. This model is spatially explicit and features the
position of each individual together with another characteristic, such as the
size of the individual, which evolves according to a given stochastic model.
The population is locally regulated through an explicit competition kernel. The
IBM is represented as a measure-valued branching/diffusing stochastic process.
The approach allows (i) to describe the associated Monte Carlo simulation and
(ii) to analyze the limit process under large initial population size
asymptotic. The limit macroscopic model is a deterministic integro-differential
equation.


  2-Ethyhexyl nitrate (2-EHN) is a major additive of fuel which is used to
comply with the cetane number of diesel. Because of its wide use and possible
accidental release, 2-EHN is a potential pollutant of the environment. In this
study, Mycobacterium austroafricanum IFP 2173 was selected among several
strains as the best 2-EHN degrader. The 2-EHN biodegradation rate was increased
in biphasic cultures where the hydrocarbon was dissolved in an inert
non-aqueous phase liquid (NAPL), suggesting that the transfer of the
hydrophobic substrate to the cells was a growth-limiting factor. Carbon balance
calculation as well as organic carbon measurement indicated a release of
metabolites in the culture medium. Further analysis by gas chromatography
revealed that a single metabolite accumulated during growth. This metabolite
had a molecular mass of 114 Da as determined by GC/MS and was provisionally
identified as 4-ethyldihydrofuran-2(3H)-one by LC-MS/MS analysis.
Identification was confirmed by analysis of the chemically synthesized lactone.
Based on these results, a plausible catabolic pathway is proposed whereby 2-EHN
is converted to 4-ethyldihydrofuran-2(3H)-one, which cannot be metabolised
further by strain IFP 2173. This putative pathway provides an explanation for
the low energetic efficiency of 2-EHN degradation and its poor
biodegradability.


  In this study, the genes involved in the initial attack on fluorene by
Sphingomonas sp. LB126 were investigated. The ? and ? subunits of a dioxygenase
complex (FlnA1A2), showing 63% and 51% sequence identity respectively, with the
subunits of an angular dioxygenase from Gram-positive Terrabacter sp. DBF63,
were identified. When overexpressed in E. coli, FlnA1A2 was responsible for the
angular oxidation of fluorene, fluorenol, fluorenone, dibenzofuran and
dibenzo-p-dioxin. Moreover, FlnA1A2 was able to oxidize polycyclic aromatic
hydrocarbons and heteroaromatics, some of which were not oxidized by the
dioxygenase from Terrabacter sp. DBF63. Quantification of resulting oxidation
products showed that fluorene and phenanthrene were preferred substrates.


  We extend recent methods for parametric sequence alignment to the parameter
space for scoring RNA folds. This involves the construction of an RNA polytope.
A vertex of this polytope corresponds to RNA secondary structures with common
branching. We use this polytope and its normal fan to study the effect of
varying three parameters in the free energy model that are not determined
experimentally. Our results indicate that variation of these specific
parameters does not have a dramatic effect on the structures predicted by the
free energy model. We additionally map a collection of known RNA secondary
structures to the RNA polytope.


  We study a simplified model of the RNA molecule proposed by G. Vernizzi, H.
Orland and A. Zee in the regime of strong concentration of positive ions in
solution. The model considers a flexible chain of equal bases that can pairwise
interact with any other one along the chain, while preserving the property of
saturation of the interactions. In the regime considered, we observe the
emergence of a critical temperature T_c separating two phases that can be
characterized by the topology of the predominant configurations: in the large
temperature regime, the dominant configurations of the molecule have very large
genera (of the order of the size of the molecule), corresponding to a complex
topology, whereas in the opposite regime of low temperatures, the dominant
configurations are simple and have the topology of a sphere. We determine that
this topological phase transition is of first order and provide an analytic
expression for T_c. The regime studied for this model exhibits analogies with
that for the dense polymer systems studied by de Gennes


  To avoid pathological platelet aggregation by von Willebrand factor (VWF),
VWF multimers are regulated in size and reactivity for adhesion by
ADAMTS13-mediated proteolysis in a shear flow dependent manner. We examined if
tensile stress in VWF under shear flow activates the VWF A2 domain for cleavage
by ADAMTS13 using molecular dynamics simulations. We indeed observed stepwise
unfolding of A2 and exposure of its deeply buried ADAMTS13 cleavage site.
Interestingly, disulfide bonds in the adjacent and highly homologous VWF A1 and
A3 domains obstruct their mechanical unfolding. We generated a full length
mutant VWF featuring a homologous disulfide bond in A2 (N1493C and C1670S), in
an attempt to lock A2 against unfolding. We find this mutant to feature
ADAMTS13-resistant behavior in vitro. Our results yield molecular-detail
evidence for the force-sensoring function of VWF A2, by revealing how tension
in VWF due to shear flow selectively exposes the A2 proteolysis site to
ADAMTS13 for cleavage while keeping the folded remainder of A2 intact and
functional. We find the unconventional knotted Rossman fold of A2 to be the key
to this mechanical response, tailored for regulating VWF size and activity.
Based on our model we can explain the pathomechanism of some natural mutations
in the VWF A2 domain that significantly increase the cleavage by ADAMTS13
without shearing or chemical denaturation, and provide with the
cleavage-activated A2 conformation a structural basis for the design of
inhibitors for VWF type 2 diseases.


  Previous studies have described the concept of peptide strings in qualitative
terms and illustrated it by applying its corrolaries in order to elucidate
basic questions in oncology and rheumatology. The present investigation is the
first to quantify these potential sub- and transcellular phenomena.
Accordingly, the propagation of peptide strings is proposed here to occur by
way of waves that in turn are subject to the energy equation established by
Planck. As a result of these insights, widespread future applications can now
be envisaged for peptide strings both in molecular medicine and quantum optics.


  The study of the response of complex dynamical social, biological, or
technological networks to external perturbations has numerous applications.
Random Boolean Networks (RBNs) are commonly used a simple generic model for
certain dynamics of complex systems. Traditionally, RBNs are interconnected
randomly and without considering any spatial extension and arrangement of the
links and nodes. However, most real-world networks are spatially extended and
arranged with regular, power-law, small-world, or other non-random connections.
Here we explore the RBN network topology between extreme local connections,
random small-world, and pure random networks, and study the damage spreading
with small perturbations. We find that spatially local connections change the
scaling of the relevant component at very low connectivities ($\bar{K} \ll 1$)
and that the critical connectivity of stability $K_s$ changes compared to
random networks. At higher $\bar{K}$, this scaling remains unchanged. We also
show that the relevant component of spatially local networks scales with a
power-law as the system size N increases, but with a different exponent for
local and small-world networks. The scaling behaviors are obtained by
finite-size scaling. We further investigate the wiring cost of the networks.
  From an engineering perspective, our new findings provide the key design
trade-offs between damage spreading (robustness), the network's wiring cost,
and the network's communication characteristics.


  A discrete-state model of the F1-ATPase molecular motor is developed which
describes not only the dependences of the rotation and ATP consumption rates on
the chemical concentrations of ATP, ADP, and inorganic phosphate, but also on
mechanical control parameters such as the friction coefficient and the external
torque. The dependence on these mechanical parameters is given to the
discrete-state model by fitting its transition rates to the continuous-angle
model of P. Gaspard and E. Gerritsma [J. Theor. Biol. 247 (2007) 672-686]. This
discrete-state model describes the behavior of the F1 motor in the regime of
tight coupling between mechanical motion and chemical reaction. In this way,
kinetic and thermodynamic properties of the F1 motor are obtained such as the
Michaelis-Menten dependence of the rotation and ATP consumption rates on ATP
concentration and its extension in the presence of ADP and Pi, their
dependences on friction and external torque, as well as the chemical and
mechanical thermodynamic efficiencies.


  Game Theory has been frequently applied in biological research since 1970s.
While the key idea of Game Theory is Nash Equilibrium, it is critical to
understand and figure out the payoff matrix in order to calculate Nash
Equilibrium. In this paper we present a dynamic programming implemented method
to compute 2x2 non-cooperative finite resource allocation game's payoff matrix.
We assume in one population there exists two types of individuals, aggressive
and non-aggressive and each individual has equal and finite resource. The
strength of individual could be described by a function of resource consumption
in discrete development stages. Each individual undergoes logistic growth hence
we divide the development into three stages: initialization, quasilinear growth
and termination. We first discuss the theoretical frame of how to dynamic
programming to calculate payoff matrix then give three numerical examples
representing three different types of aggressive individuals and calculate the
payoff matrix for each of them respectively. Based on the numerical payoff
matrix we further investigate the evolutionary stable strategies (ESS) of the
games.


  One of the characteristic features of genetic networks is their inherent
robustness, that is, their ability to retain functionality in spite of the
introduction of random errors. In this paper, we seek to better understand how
robustness is achieved and what functionalities can be maintained robustly. Our
goal is to formalize some of the language used in biological discussions in a
reasonable mathematical framework, where questions can be answered in a
rigorous fashion. These results provide basic conceptual understanding of
robust regulatory networks that should be valuable independent of the details
of the formalism.
  We model the gene regulatory network as a boolean network, a general and
well-established model introduced by Stuart Kauffman. A boolean network is said
to be in a viable configuration if the node states of the network at its
fixpoint satisfy some given constraint. We specify how mutations affect the
behavior of the boolean network. A network is then said to be robust if most
random mutations to the network reach a viable configuration. The main question
investigated in our study is: given a constraint on the fixpoint configuration,
does there exist a network that is robust with respect to it and, if so, what
is its structure? We demonstrate both explicit constructions of robust networks
as well as negative results disproving their existence.


  Siphons in a chemical reaction system are subsets of the species that have
the potential of being absent in a steady state. We present a characterization
of minimal siphons in terms of primary decomposition of binomial ideals, we
explore the underlying geometry, and we demonstrate the effective computation
of siphons using computer algebra software. This leads to a new method for
determining whether given initial concentrations allow for various boundary
steady states.


  Chromatin organization and dynamics is studied in this work at scales ranging
from single nucleosome to nucleosomal array by using a unique combination of
biochemical assays, single molecule imaging technique and numerical modeling.
We demonstrate that a subtle modification in the nucleosome structure induced
by the histone variant H2A.Bbd drastically modifies the higher order
organization of the nucleosomal arrays. Importantly, as directly visualized by
AFM, conventional H2A nucleosomal arrays exhibit specific local organization,
in contrast to H2A.Bbd arrays, which show "beads on a string" structure. The
combination of systematic image analysis and theoretical modeling allows a
quantitative description relating the observed gross structural changes of the
arrays to their local organization. Our results strongly suggest that
higher-order organization of H1-free nucleosomal arrays is mainly determined by
the fluctuation properties of individual nucleosomes. Moreover, numerical
simulations suggest the existence of attractive interactions between
nucleosomes to provide the degree of compaction observed for conventional
chromatin fibers.


  We study the problem of computing outer bounds for the region of steady
states of biochemical reaction networks modelled by ordinary differential
equations, with respect to parameters that are allowed to vary within a
predefined region. Using a relaxed version of the corresponding feasibility
problem and its Lagrangian dual, we show how to compute certificates for
regions in state space not containing any steady states. Based on these
results, we develop an algorithm to compute outer bounds for the region of all
feasible steady states. We apply our algorithm to the sensitivity analysis of a
Goldbeter--Koshland enzymatic cycle, which is a frequent motif in reaction
networks for regulation of metabolism and signal transduction.


  We review and extend results for mutation, selection, genetic drift, and
migration in a one-dimensional continuous population. The population is
described by a continuous limit of the stepping stone model, which leads to the
stochastic Fisher-Kolmogorov-Petrovsky-Piscounov equation with additional terms
describing mutations. Although the stepping stone model was first proposed for
population genetics, it is closely related to "voter models" of interest in
nonequilibrium statistical mechanics. The stepping stone model can also be
regarded as an approximation to the dynamics of a thin layer of actively
growing pioneers at the frontier of a colony of microorganisms undergoing a
range expansion on a Petri dish. We find that the population tends to segregate
into monoallelic domains. This segregation slows down genetic drift and
selection because these two evolutionary forces can only act at the boundaries
between the domains; the effects of mutation, however, are not significantly
affected by the segregation. Although fixation in the neutral well-mixed (or
"zero dimensional") model occurs exponentially in time, it occurs only
algebraically fast in the one-dimensional model. If selection is weak,
selective sweeps occur exponentially fast in both well-mixed and
one-dimensional populations, but the time constants are different. We also find
an unusual sublinear increase in the variance of the spatially averaged allele
frequency with time. Although we focus on two alleles or variants, q-allele
Potts-like models of gene segregation are considered as well. We also
investigate the effects of geometry at the frontier by considering growth of
circular colonies. Our analytical results are checked with simulations, and
could be tested against recent spatial experiments on range expansions off
linear inoculations of Escherichia coli and Saccharomyces cerevisiae.


  We define a new parameter to quantify the antigenic distance between two H3N2
influenza strains: we use this parameter to measure antigenic distance between
circulating H3N2 strains and the closest vaccine component of the influenza
vaccine. For the data between 1971 and 2004, the measure of antigenic distance
correlates better with efficacy in humans of the H3N2 influenza A annual
vaccine than do current state of the art measures of antigenic distance such as
phylogenetic sequence analysis or ferret antisera inhibition assays. We suggest
that this measure of antigenic distance could be used to guide the design of
the annual flu vaccine. We combine the measure of antigenic distance with a
multiple-strain avian influenza transmission model to study the threat of
simultaneous introduction of multiple avian influenza strains. For H3N2
influenza, the model is validated against observed viral fixation rates and
epidemic progression rates from the World Health Organization FluNet - Global
Influenza Surveillance Network. We find that a multiple-component avian
influenza vaccine is helpful to control a simultaneous multiple introduction of
bird-flu strains. We introduce Population at Risk (PaR) to quantify the risk of
a flu pandemic, and calculate by this metric the improvement that a multiple
vaccine offers.


  We add here compared to our former arXiv version an explicit expression for
the descendant ratio along the generations The equation that is added here
appeared in the Hebrew version published in BDD, Bar Ilan University Press, 23,
71, 2010, titled The distribution route from ancestors to descendants in
Equation 11. Otherwise, we don not change nor show here the paper but leave and
direct the reader to the former arXiv version with the above addition or to the
published paper in Hebrew.


  The functioning of a living cell is largely determined by the structure of
its regulatory network, comprising non-linear interactions between regulatory
genes. An important factor for the stability and evolvability of such
regulatory systems is neutrality - typically a large number of alternative
network structures give rise to the necessary dynamics. Here we study the
discretized regulatory dynamics of the yeast cell cycle [Li et al., PNAS, 2004]
and the set of networks capable of reproducing it, which we call functional.
Among these, the empirical yeast wildtype network is close to optimal with
respect to sparse wiring. Under point mutations, which establish or delete
single interactions, the neutral space of functional networks is fragmented
into 4.7 * 10^8 components. One of the smaller ones contains the wildtype
network. On average, functional networks reachable from the wildtype by
mutations are sparser, have higher noise resilience and fewer fixed point
attractors as compared with networks outside of this wildtype component.


  The kinetics for the assembly of viral proteins into a population of capsids
can be measured in vitro with size exclusion chromatography or dynamic light
scattering, but extracting mechanistic information from these studies is
challenging. For example, it is not straightforward to determine the critical
nucleus size or the elongation time (the time required for a nucleated partial
capsid to grow completion). We show that, for two theoretical models of capsid
assembly, the critical nucleus size can be determined from the concentration
dependence of the assembly reaction half-life and the elongation time is
revealed by the length of the lag phase. Furthermore, we find that the system
becomes kinetically trapped when nucleation becomes fast compared to
elongation. Implications of this constraint for determining elongation
mechanisms from experimental assembly data are discussed.


  Single molecule Forster resonance energy transfer (FRET) experiments are used
to infer the properties of the denatured state ensemble (DSE) of proteins. From
the measured average FRET efficiency, <E>, the distance distribution P(R) is
inferred by assuming that the DSE can be described as a polymer. The single
parameter in the appropriate polymer model (Gaussian chain, Wormlike chain, or
Self-avoiding walk) for P(R) is determined by equating the calculated and
measured <E>. In order to assess the accuracy of this "standard procedure," we
consider the generalized Rouse model (GRM), whose properties [<E> and P(R)] can
be analytically computed, and the Molecular Transfer Model for protein L for
which accurate simulations can be carried out as a function of guanadinium
hydrochloride (GdmCl) concentration. Using the precisely computed <E> for the
GRM and protein L, we infer P(R) using the standard procedure. We find that the
mean end-to-end distance can be accurately inferred (less than 10% relative
error) using <E> and polymer models for P(R). However, the value extracted for
the radius of gyration (Rg) and the persistence length (lp) are less accurate.
The relative error in the inferred R-g and lp, with respect to the exact
values, can be as large as 25% at the highest GdmCl concentration. We propose a
self-consistency test, requiring measurements of <E> by attaching dyes to
different residues in the protein, to assess the validity of describing DSE
using the Gaussian model. Application of the self-consistency test to the GRM
shows that even for this simple model the Gaussian P(R) is inadequate. Analysis
of experimental data of FRET efficiencies for the cold shock protein shows that
at there are significant deviations in the DSE P(R) from the Gaussian model.


  A mathematical model to describe period-memorizing behavior in Physarum
plasmodium are reported. In constructing the model, we first examine the basic
characteristics required for the class of models, then create a minimal linear
model to fulfill these requirements. We also propose two modifications of the
minimal model, nonlinearization and noise addition, which improve the
reproducibility of experimental evidences. Differences in the mechanisms and in
the reproducibility of experiments between our models and the previous models
are discussed.