  The induction of late long-term potentiation (L-LTP) involves complex
interactions among second messenger cascades. To gain insights into these
interactions, a mathematical model was developed for L-LTP induction in the CA1
region of the hippocampus. The differential equation-based model represents
actions of protein kinase A (PKA), MAP kinase (MAPK), and CaM kinase II
(CAMKII) in the vicinity of the synapse, and activation of transcription by CaM
kinase IV (CAMKIV) and MAPK. L-LTP is represented by increases in a synaptic
weight. Simulations suggest that steep, supralinear stimulus-response
relationships between stimuli (elevations in [Ca2+]) and kinase activation are
essential for translating brief stimuli into long-lasting gene activation and
synaptic weight increases. Convergence of multiple kinase activities to induce
L-LTP helps to generate a threshold whereby the amount of L-LTP varies steeply
with the number of tetanic electrical stimuli. The model simulates tetanic,
theta-burst, pairing-induced, and chemical L-LTP, as well as L-LTP due to
synaptic tagging. The model also simulates inhibition of L-LTP by inhibition of
MAPK, CAMKII, PKA, or CAMKIV. The model predicts results of experiments to
delineate mechanisms underlying L-LTP induction and expression. For example,
the cAMP antagonist RpcAMPs, which inhibits L-LTP induction, is predicted to
inhibit ERK activation. The model also appears useful to clarify similarities
and differences between hippocampal L-LTP and long-term synaptic strengthening
in other systems.


  Sex determination in the order Hymenoptera is based on arrhenotoky,
hymenopteran males are usually haploid and females diploid. Males of the
Ichneumonidae Diadromus pulchellus, solitary endoparasitoid of A. assectella
pupae are normally haploid, but diploid males are present in a natural
population and can be obtained in a experimental population. The future of an
ovocyte laid by mated females of the solitary endoparasitoid Diadromus
pulchellus was characterised by 6 probabilities related to the sex and the
development of the ovocyte. The probabilities of fertilisation of female
ovocyte (k1) or non-fertilisation (k3) showed that an inseminated female
functioned as a unmated female for half of the time (since k1= 0.492 and k3 =
0.455) with the probability of fertilisation of male ovocyte (k2 = 0.053).The
survival probabilities of each type of ovocyte showed that an ovocyte had a
high probability of developing up to the adult stage, although the difference
between the calculated sex ratios at laying (males / females = 1.032) and at
emergence (0.90) revealed a slight reduction in the number of haploid sons. The
probabilities of fertilisation and of viability of all the ovocytes laid by
each of the 33 mated females were analysed by an ascending hierarchical
classification of Euclidean distances and by an analysis of their principal
components. The 33 mothers were distributed into 4 distinct sub-groups
characterised by a sex ratio varying from an exclusive presence of females to
an exclusive presence of males. Our hypothesis was that this distribution in 4
sub-sets could not simply result from the random nature of the sample.


  Single-molecule tracking is a powerful way to look at the dynamic
organization of plasma membranes. However, there are some limitations to its
use. For example, it was recently observed, using numerical simulation, that
time-averaging effects inherent to the exposure time of detectors are likely to
bias the apparent motion of molecules confined in microdomains. Here, we solve
this apparently limiting issue analytically. We explore this phenomenon by
calculating its effects on the observed diffusion coefficients and domain
sizes. We demonstrate that the real parameters can be easily recovered from the
measured apparent ones. Interestingly, we find that single-molecule tracking
can be used to explore events occurring at a timescale smaller than the
exposure time.


This paper contains some immature results. It have been withdrawn.


  The influence of intrinsic channel noise on the spontaneous spiking activity
of poisoned excitable membrane patches is studied by use of a stochastic
generalization of the Hodgkin-Huxley model. Internal noise stemming from the
stochastic dynamics of individual ion channels is known to affect the
collective properties of the whole ion channel cluster. For example, there
exists an optimal size of the membrane patch for which the internal noise alone
causes a regular spontaneous generation of action potentials. In addition to
varying the size of ion channel clusters, living organisms may adapt the
densities of ion channels in order to optimally regulate the spontaneous
spiking activity. The influence of channel block on the excitability of a
membrane patch of certain size is twofold: First, a variation of ion channel
densities primarily yields a change of the conductance level. Second, a
down-regulation of working ion channels always increases the channel noise.
While the former effect dominates in the case of sodium channel block resulting
in a reduced spiking activity, the latter enhances the generation of
spontaneous action potentials in the case of a tailored potassium channel
blocking. Moreover, by blocking some portion of either potassium or sodium ion
channels, it is possible to either increase or to decrease the regularity of
the spike train.


  Motivation: Its well known the use of electric current in the cleaning of
viral diseases in plants instead a deep knowledge of this phenomenas
theoretical basis is not yet available. The description of the real causes of
nucleoprotein inactivation, will contribute to the optimization of further
experiments in order to obtain more and more efficient cleaning methodologies
for starting material supplies of vegetal species micropropagation process. The
dissipated energy as heat will depend on the specific treatment applied and on
the physical properties of the vegetal tissue and viral nucleoprotein. Results:
A hyperbolic dependence between the absorbance output values, obtained during a
diagnostic micro-ELISA experiment, and the electrical power (in watts) applied
to each explant was found. The former, demonstrate the cleaning process nature
is essentially the effect over viral nucleoprotein denaturation, by means of
the heat, to which they were exposed in the thermal bath that the vegetal
tissue constitutes. Our results are consistent with mathematical developed from
theoretical frame of harmonic oscillators in which molecular machine (viral
nucleoproteins) are redefined.


  Two complementary analyses of a cyclic negative feedback system with delay
are considered in this paper. The first analysis applies the work by Sontag,
Angeli, Enciso and others regarding monotone control systems under negative
feedback, and it implies the global attractiveness towards an equilibrium for
arbitrary delays. The second one concerns the existence of a Hopf bifurcation
on the delay parameter, and it implies the existence of nonconstant periodic
solutions for special delay values. A key idea is the use of the Schwarzian
derivative, and its application for the study of Michaelis-Menten
nonlinearities. The positive feedback case is also addressed.


  Recently experiments showed that some adhesive receptor-ligand complexes
increase their lifetimes when they are stretched by mechanical force, while the
force increase beyond some thresholds their lifetimes decrease. Several
specific chemical kinetic models have been developed to explain the intriguing
transitions from the "catch-bonds" to the "slip-bonds". In this work we suggest
that the counterintuitive forced dissociation of the complexes is a typical
rate process with dynamic disorder. An uniform one-dimension force modulating
Agmon-Hopfield model is used to quantitatively describe the transitions
observed in the single bond P-selctin glycoprotein ligand
1(PSGL-1)$-$P-selectin forced dissociation experiments, which were respectively
carried out on the constant force [Marshall, {\it et al.}, (2003) Nature {\bf
423}, 190-193] and the force steady- or jump-ramp [Evans {\it et al.}, (2004)
Proc. Natl. Acad. Sci. USA {\bf 98}, 11281-11286] modes. Our calculation shows
that the novel catch-slip bond transition arises from a competition of the two
components of external applied force along the dissociation reaction coordinate
and the complex conformational coordinate: the former accelerates the
dissociation by lowering the height of the energy barrier between the bound and
free states (slip), while the later stabilizes the complex by dragging the
system to the higher barrier height (catch).


  We study folding dynamics of protein-like sequences on square lattice using
physical move set that exhausts all possible conformational changes. By
analytically solving the master equation, we follow the time-dependent
probabilities of occupancy of all 802,075 conformations of 16-mers over
7-orders of time span. We find that (i) folding rates of these protein-like
sequences of same length can differ by 4-orders of magnitude, (ii) folding
rates of sequences of the same conformation can differ by a factor of 190, and
(iii) parameters of the native structures, designability, and thermodynamic
properties are weak predictors of the folding rates, rather, basin analysis of
the kinematic energy landscape defined by the moves can provide excellent
account of the observed folding rates.


  The amino acid sequences of proteins provide rich information for inferring
distant phylogenetic relationships and for predicting protein functions.
Estimating the rate matrix of residue substitutions from amino acid sequences
is also important because the rate matrix can be used to develop scoring
matrices for sequence alignment. Here we use a continuous time Markov process
to model the substitution rates of residues and develop a Bayesian Markov chain
Monte Carlo method for rate estimation. We validate our method using simulated
artificial protein sequences. Because different local regions such as binding
surfaces and the protein interior core experience different selection pressures
due to functional or stability constraints, we use our method to estimate the
substitution rates of local regions. Our results show that the substitution
rates are very different for residues in the buried core and residues on the
solvent exposed surfaces. In addition, the rest of the proteins on the binding
surfaces also have very different substitution rates from residues. Based on
these findings, we further develop a method for protein function prediction by
surface matching using scoring matrices derived from estimated substitution
rates for residues located on the binding surfaces. We show with examples that
our method is effective in identifying functionally related proteins that have
overall low sequence identity, a task known to be very challenging.


  This chapter discusses geometric models of biomolecules and geometric
constructs, including the union of ball model, the weigthed Voronoi diagram,
the weighted Delaunay triangulation, and the alpha shapes. These geometric
constructs enable fast and analytical computaton of shapes of biomoleculres
(including features such as voids and pockets) and metric properties (such as
area and volume). The algorithms of Delaunay triangulation, computation of
voids and pockets, as well volume/area computation are also described. In
addition, applications in packing analysis of protein structures and protein
function prediction are also discussed.


  Morphogens are proteins, often produced in a localised region, whose
concentrations spatially demarcate regions of differing gene expression in
developing embryos. The boundaries of expression must be set accurately and in
proportion to the size of the one-dimensional developing field; this cannot be
accomplished by a single gradient. Here, we show how a pair of morphogens
produced at opposite ends of a developing field can solve the pattern-scaling
problem. In the most promising scenario, the morphogens effectively interact
according to the annihilation reaction $A+B\to\emptyset$ and the switch occurs
according to the absolute concentration of $A$ or $B$. In this case embryonic
markers across the entire developing field scale approximately with system
size; this cannot be achieved with a pair of non-interacting gradients that
combinatorially regulate downstream genes. This scaling occurs in a window of
developing-field sizes centred at a few times the morphogen decay length.


  Pathogen-mediated competition, through which an invasive species carrying and
transmitting a pathogen can be a superior competitor to a more vulnerable
resident species, is one of the principle driving forces influencing
biodiversity in nature. Using an experimental system of bacteriophage-mediated
competition in bacterial populations and a deterministic model, we have shown
in [Joo et al 2005] that the competitive advantage conferred by the phage
depends only on the relative phage pathology and is independent of the initial
phage concentration and other phage and host parameters such as the
infection-causing contact rate, the spontaneous and infection-induced lysis
rates, and the phage burst size. Here we investigate the effects of stochastic
fluctuations on bacterial invasion facilitated by bacteriophage, and examine
the validity of the deterministic approach. We use both numerical and
analytical methods of stochastic processes to identify the source of noise and
assess its magnitude. We show that the conclusions obtained from the
deterministic model are robust against stochastic fluctuations, yet deviations
become prominently large when the phage are more pathological to the invading
bacterial strain.


  Scale-free foraging patterns are widespread among animals. These may be the
outcome of an optimal searching strategy to find scarce randomly distributed
resources, but a less explored alternative is that this behaviour may result
from the interaction of foraging animals with a particular distribution of
resources. We introduce a simple foraging model where individuals follow mental
maps and choose their displacements according to a maximum efficiency
criterion, in a spatially disordered environment containing many trees with a
heterogeneous size distribution. We show that a particular tree size frequency
distribution induces non-Gaussian movement patterns with multiple spatial
scales (L\'evy walks). These results are consistent with tree size variation
and Spider monkey (Ateles geoffroyi) foraging patterns. We discuss the
consequences that our results may have for the patterns of seed dispersal by
foraging primates.


  In addition to their self-renewal capabilities, hematopoietic stem cells
guarantee the continuous supply of fully differentiated, functional cells of
various types in the peripheral blood. The process which controls
differentiation into the different lineages of the hematopoietic system
(erythroid, myeloid, lymphoid) is referred to as lineage specification. It
requires a potentially multi-step decision sequence which determines the fate
of the cells and their successors. It is generally accepted that lineage
specification is regulated by a complex system of interacting transcription
factors. However, the underlying principles controlling this regulation are
currently unknown.
  Here, we propose a simple quantitative model describing the interaction of
two transcription factors. This model is motivated by experimental observations
on the transcription factors GATA-1 and PU.1, both known to act as key
regulators and potential antagonists in the erythroid vs. myeloid
differentiation processes of hematopoietic progenitor cells. We demonstrate the
ability of the model to account for the observed switching behavior of a
transition from a state of low expression of both factors (undifferentiated
state) to the dominance of one factor (differentiated state). Depending on the
parameter choice, the model predicts two different possibilities to explain the
experimentally suggested, stem cell characterizing priming state of low level
co-expression. Whereas increasing transcription rates are sufficient to induce
differentiation in one scenario, an additional system perturbation (by
stochastic fluctuations or directed impulses) of transcription factor levels is
required in the other case.


  This chapter discusses theoretical framework and methods for developing
knowledge-based potential functions essential for protein structure prediction,
protein-protein interaction, and protein sequence design. We discuss in some
details about the Miyazawa-Jernigan contact statistical potential,
distance-dependent statistical potentials, as well as geometric statistical
potentials. We also describe a geometric model for developing both linear and
non-linear potential functions by optimization. Applications of knowledge-based
potential functions in protein-decoy discrimination, in protein-protein
interactions, and in protein design are then described. Several issues of
knowledge-based potential functions are finally discussed.


  $\beta$-barrel membrane proteins are found in the outer membrane of
gram-negative bacteria, mitochondria, and chloroplasts. We have developed
probabilistic models to quantify propensities of residues for different spatial
locations and for interstrand pairwise contact interactions involving strong
H-bonds, side-chain interactions, and weak H-bonds. The propensity values and
p-values measuring statistical significance are calculated exactly by
analytical formulae we have developed. Contrary to the ``positive-inside'' rule
for helical membrane proteins, $\beta$-barrel membrane proteins follow a
significant albeit weaker ``positive-outside'' rule, in that the basic residues
Arg and Lys are disproportionately favored in the extracellular cap region and
disfavored in the periplasmic cap region. Different residue pairs prefer strong
backbone H-bonded interstrand pairings (e.g. Gly-Aromatic) or non-H-bonded
pairings (e.g. Aromatic-Aromatic). In addition, Tyr and Phe participate in
aromatic rescue by shielding Gly from polar environments. These propensities
can be used to predict the registration of strand pairs, an important task for
the structure prediction of $\beta$-barrel membrane proteins. Our accuracy of
44% is considerably better than random (7%) and other studies. Our results
imply several experiments that can help to elucidate the mechanisms of in vitro
and in vivo folding of $\beta$-barrel membrane proteins. See supplementary
material after the bibliography for detailed techniques.


  An effective potential function is critical for protein structure prediction
and folding simulation. Simplified protein models such as those requiring only
$C_\alpha$ or backbone atoms are attractive because they enable efficient
search of the conformational space. We show residue specific reduced discrete
state models can represent the backbone conformations of proteins with small
RMSD values. However, no potential functions exist that are designed for such
simplified protein models. In this study, we develop optimal potential
functions by combining contact interaction descriptors and local
sequence-structure descriptors. The form of the potential function is a
weighted linear sum of all descriptors, and the optimal weight coefficients are
obtained through optimization using both native and decoy structures. The
performance of the potential function in test of discriminating native protein
structures from decoys is evaluated using several benchmark decoy sets. Our
potential function requiring only backbone atoms or $C_\alpha$ atoms have
comparable or better performance than several residue-based potential functions
that require additional coordinates of side chain centers or coordinates of all
side chain atoms. By reducing the residue alphabets down to size 5 for local
structure-sequence relationship, the performance of the potential function can
be further improved. Our results also suggest that local sequence-structure
correlation may play important role in reducing the entropic cost of protein
folding.


  The evolution of two species with different fitness is investigated on
degree-heterogeneous graphs. The population evolves either by one individual
dying and being replaced by the offspring of a random neighbor (voter model
(VM) dynamics) or by an individual giving birth to an offspring that takes over
a random neighbor node (invasion process (IP) dynamics). The fixation
probability for one species to take over a population of N individuals depends
crucially on the dynamics and on the local environment. Starting with a single
fitter mutant at a node of degree k, the fixation probability is proportional
to k for VM dynamics and to 1/k for IP dynamics.


  Without invoking the Markov approximation, we derive formulas for vibrational
energy relaxation (VER) and dephasing for an anharmonic system oscillator using
a time-dependent perturbation theory. The system-bath Hamiltonian contains more
than the third order coupling terms since we take a normal mode picture as a
zeroth order approximation. When we invoke the Markov approximation, our theory
reduces to the Maradudin-Fein formula which is used to describe VER properties
of glass and proteins. When the system anharmonicity and the renormalization
effect due to the environment vanishes, our formulas reduce to those derived by
Mikami and Okazaki invoking the path-integral influence functional method [J.
Chem. Phys. 121 (2004) 10052]. We apply our formulas to VER of the amide I mode
of a small amino-acide like molecule, N-methylacetamide, in heavy water.


  The threat of avian influenza and the 2004-2005 influenza vaccine supply
shortage in the United States has sparked a debate about optimal vaccination
strategies to reduce the burden of morbidity and mortality caused by the
influenza virus. We present a comparative analysis of two classes of suggested
vaccination strategies: mortality-based strategies that target high risk
populations and morbidity-based that target high prevalence populations.
Applying the methods of contact network epidemiology to a model of disease
transmission in a large urban population, we evaluate the efficacy of these
strategies across a wide range of viral transmission rates and for two
different age-specific mortality distributions. We find that the optimal
strategy depends critically on the viral transmission level (reproductive rate)
of the virus: morbidity-based strategies outperform mortality-based strategies
for moderately transmissible strains, while the reverse is true for highly
transmissible strains. These results hold for a range of mortality rates
reported for prior influenza epidemics and pandemics. Furthermore, we show that
vaccination delays and multiple introductions of disease into the community
have a more detrimental impact on morbidity-based strategies than
mortality-based strategies. If public health officials have reasonable
estimates of the viral transmission rate and the frequency of new introductions
into the community prior to an outbreak, then these methods can guide the
design of optimal vaccination priorities. When such information is unreliable
or not available, as is often the case, this study recommends mortality-based
vaccination priorities.


  We have analyzed the relations between the mutational pressure, recombination
and selection pressure in the bit-string model with sexual reproduction. For
specific sets of these parameters we have found three phase transitions with
one phase where populations can survive. In this phase, recombination enhances
the survival probability. Even if recombination is associated, to some extent,
with additional mutations it could be advantageous to reproduction, indicating
that the frequencies of recombinations and recombination-associated mutations
can self-organize in Nature. Partitioning the diploid genome into pairs of
chromosomes independently assorted during gamete production enables
recombinations between groups of genes without the risk of mutations and is
also advantageous for the strategy of sexual reproduction.


  We have analysed an effect of the Bak-Sneppen predator-prey food-chain
self-organization on nucleotide content of evolving species. In our model,
genomes of the species under consideration have been represented by their
nucleotide genomic fraction and we have applied two-parameter Kimura model of
substitutions to include the changes of the fraction in time. The initial
nucleotide fraction and substitution rates were decided with the help of random
number generator. Deviation of the genomic nucleotide fraction from its
equilibrium value was playing the role of the fitness parameter, $B$, in
Bak-Sneppen model. Our finding is, that the higher is the value of the
threshold fitness, during the evolution course, the more frequent are large
fluctuations in number of species with strongly differentiated nucleotide
content; and it is more often the case that the oldest species, which survive
the food-chain competition, might have specific nucleotide fraction making
possible generating long genes


  Understanding the system level adaptive changes taking place in an organism
in response to variations in the environment is a key issue of contemporary
biology. Current modeling approaches such as the constraint-based flux balance
analyses (FBA) have proved highly successful in analyzing the capabilities of
cellular metabolism, including its capacity to predict deletion phenotypes, the
ability to calculate the relative flux values of metabolic reactions and the
properties of alternate optimal growth states. Here, we use FBA to thoroughly
assess the activity of the Escherichia coli, Helicobacter pylori, and
Saccharomyces cerevisiae metabolism in 30,000 diverse simulated environments.
We identify a set of metabolic reactions forming a connected metabolic core
that carry non-zero fluxes under all growth conditions, and whose flux
variations are highly correlated. Furthermore, we find that the enzymes
catalyzing the core reactions display a considerably higher fraction of
phenotypic essentiality and evolutionary conservation than those catalyzing
non-core reactions. Cellular metabolism is characterized by a large number of
species-specific conditionally-active reactions organized around an
evolutionary conserved always active metabolic core. Finally, we find that most
current antibiotics interfering with the bacterial metabolism target the core
enzymes, indicating that our findings may have important implications for
antimicrobial drug target discovery.


  If constraints are imposed on a macromolecule, two inequivalent classical
models may be used: the stiff and the rigid one. This work studies the effects
of such constraints on the Conformational Equilibrium Distribution (CED) of the
model dipeptide HCO-L-Ala-NH2 without any simplifying assumption. We use ab
initio Quantum Mechanics calculations including electron correlation at the MP2
level to describe the system, and we measure the conformational dependence of
all the correcting terms to the naive CED based in the Potential Energy Surface
(PES) that appear when the constraints are considered. These terms are related
to mass-metric tensors determinants and also occur in the Fixman's compensating
potential. We show that some of the corrections are non-negligible if one is
interested in the whole Ramachandran space. On the other hand, if only the
energetically lower region, containing the principal secondary structure
elements, is assumed to be relevant, then, all correcting terms may be
neglected up to peptides of considerable length. This is the first time, as far
as we know, that the analysis of the conformational dependence of these
correcting terms is performed in a relevant biomolecule with a realistic
potential energy function.


  Using coarse-grained model we have explored forced-unfolding of RNA hairpin
as a function of $f_S$ and the loading rate ($r_f$). The simulations and
theoretical analysis have been done without and with the handles that are
explicitly modeled by semiflexible polymer chains. The mechanisms and time
scales for denaturation by temperature jump and mechanical unfolding are vastly
different. The directed perturbation of the native state by $f_S$ results in a
sequential unfolding of the hairpin starting from their ends whereas thermal
denaturation occurs stochastically. From the dependence of the unfolding rates
on $r_f$ and $f_S$ we show that the position of the unfolding transition state
(TS) is not a constant but moves dramatically as either $r_f$ or $f_S$ is
changed. The TS movements are interpreted by adopting the Hammond postulate for
forced-unfolding. Forced-unfolding simulations of RNA, with handles attached to
the two ends, show that the value of the unfolding force increases (especially
at high pulling speeds) as the length of the handles increases. The pathways
for refolding of RNA from stretched initial conformation, upon quenching $f_S$
to the quench force $f_Q$, are highly heterogeneous. The refolding times, upon
force quench, are at least an order of magnitude greater than those obtained by
temperature quench. The long $f_Q$-dependent refolding times starting from
fully stretched states are analyzed using a model that accounts for the
microscopic steps in the rate limiting step which involves the trans to gauche
transitions of the dihedral angles in the GAAA tetraloop. The simulations with
explicit molecular model for the handles show that the dynamics of force-quench
refolding is strongly dependent on the interplay of their contour length and
the persistence length, and the RNA persistence length.


  Simplified Go models, where only native contacts interact favorably, have
proven useful to characterize some aspects of the folding of small proteins.
The success of these models is limited by the fact that all residues interact
in the same way, so that the folding features of a protein are determined only
by the geometry of its native conformation. We present an extended version of a
C-alpha based Go model where different residues interact with different
energies. The model is used to calculate the thermodynamics of three small
proteins (Protein G, SrcSH3 and CI2) and the effect of mutations on the
wildtype sequence. The model allows to investigate some of the most
controversial areas in protein folding such as its earliest stages, a subject
which has lately received particular attention. The picture which emerges for
the three proteins under study is that of a hierarchical process, where local
elementary structures (LES) (not necessarily coincident with elements of
secondary structure) are formed at the early stages of the folding and drive
the protein, through the transition state and the postcritical folding nucleus
(FN), resulting from the docking of the LES, to the native conformation.


  Post-transductional modifications tune the functions of proteins and regulate
the collective dynamics of biochemical networks that determine how cells
respond to environmental signals. For example, protein phosphorylation and
nitrosylation are well-known to play a pivotal role in the intracellular
transduction of activation and death signals. A protein can have multiple sites
where chemical groups can reversibly attach in processes such as
phosphorylation or nitrosylation. A microscopic description of these processes
must take into account the intrinsic probabilistic nature of the underlying
reactions. We apply combinatorial considerations to standard enzyme kinetics
and in this way we extend to the dynamic regime a simplified version of the
traditional models on the allosteric regulation of protein functions. We link a
generic modification chain to a downstream Michaelis-Menten enzymatic reaction
and we demonstrate numerically that this accounts both for thresholds and long
time delays in the conversion of the substrate by the enzyme. The proposed
mechanism is stable and robust and the higher the number of modification sites,
the greater the stability. We show that a high number of modification sites
converts a fast reaction into a slow process, and the slowing down depends on
the number of sites and may span many orders of magnitude; in this way
multisite modification of proteins stands out as a general mechanism that
allows the transfer of information from the very short time scales of enzyme
reactions (milliseconds) to the long time scale of cell response (hours).


  According to the 'old view', proteins fold along well-defined sequential
pathways, whereas the 'new view' sees protein folding as a highly parallel
stochastic process on funnel-shaped energy landscapes. We have analyzed
parallel and sequential processes on a large number of Molecular Dynamics
unfolding trajectories of the protein CI2 at high temperatures. Using rigorous
statistical measures, we quantify the degree of sequentiality on two structural
levels. The unfolding process is highly parallel on the microstructural level
of individual contacts. On a coarser, macrostructural level of contact
clusters, characteristic parallel and sequential events emerge. These
characteristic events can be understood from loop-closure dependencies between
the contact clusters. A correlation analysis of the unfolding times of the
contacts reveals a high degree of substructural cooperativity within the
contact clusters.


  The folding rates of two-state proteins have been found to correlate with
simple measures of native-state topology. The most prominent among these
measures is the relative contact order (CO), which is the average CO or
'localness' of all contacts in the native protein structure, divided by the
chain length. Here, we test whether such measures can be generalized to capture
the effect of chain crosslinks on the folding rate. Crosslinks change the chain
connectivity and therefore also the localness of some of the the native
contacts. These changes in localness can be taken into account by the
graph-theoretical concept of effective contact order (ECO). The relative ECO,
however, the natural extension of the relative CO for proteins with crosslinks,
overestimates the changes in the folding rates caused by crosslinks. We suggest
here a novel measure of native-state topology, the relative logCO, and its
natural extension, the relative logECO. The relative logCO is the average value
for the logarithm of the CO of all contacts, divided by the logarithm of the
chain length. The relative log(E)CO reproduces the folding rates of a set of 26
two-state proteins without crosslinks with essentially the same high
correlation coefficient as the relative CO. In addition, it also captures the
folding rates of 8 two-state proteins with crosslinks.


  Thermal unfolding of proteins is compared to folding and mechanical
stretching in a simple topology-based dynamical model. We define the unfolding
time and demonstrate its low-temperature divergence. Below a characteristic
temperature, contacts break at separate time scales and unfolding proceeds
approximately in a way reverse to folding. Features in these scenarios agree
with experiments and atomic simulations on titin.


  I address the question of the fluctuations in fishery landings. Using the
fishery statistics time-series collected by the Food and Agriculture
Organization of the United Nations since the early 1950s, I here analyze
fishing activities and find two scaling features of capture fisheries
production: (i) the standard deviation of growth rate of the domestically
landed catches decays as a power-law function of country landings with an
exponent of value 0.15; (ii) the average number of fishers in a country scales
to the 0.7 power of country landings. I show how these socio-ecological
patterns may be related, yielding a scaling relation between these exponents.
The predicted scaling relation implies that the width of the annual per capita
growth-rate distribution scales to the 0.2 power of country landings, i.e.
annual fluctuations in per capita landed catches increase with increased per
capita catches in highly producing countries. Beside the scaling behavior, I
report that fluctuations in the annual domestic landings have increased in the
last 30 years, while the mean of the annual growth rate declined significantly
after 1972.


  To develop and investigate detailed mathematical models of cellular metabolic
processes is one of the primary challenges in systems biology. However, despite
considerable advance in the topological analysis of metabolic networks,
explicit kinetic modeling based on differential equations is still often
severely hampered by inadequate knowledge of the enzyme-kinetic rate laws and
their associated parameter values. Here we propose a method that aims to give a
detailed and quantitative account of the dynamical capabilities of metabolic
systems, without requiring any explicit information about the particular
functional form of the rate equations. Our approach is based on constructing a
local linear model at each point in parameter space, such that each element of
the model is either directly experimentally accessible, or amenable to a
straightforward biochemical interpretation. This ensemble of local linear
models, encompassing all possible explicit kinetic models, then allows for a
systematic statistical exploration of the comprehensive parameter space. The
method is applied to two paradigmatic examples: The glycolytic pathway of yeast
and a realistic-scale representation of the photosynthetic Calvin cycle.


  Swimming bacteria detect chemical gradients by performing temporal
comparisons of recent measurements of chemical concentration. These comparisons
are described quantitatively by the chemotactic response function, which we
expect to optimize chemotactic behavioral performance. We identify two
independent chemotactic performance criteria: in the short run, a favorable
response function should move bacteria up chemoattractant gradients, while in
the long run, bacteria should aggregate at peaks of chemoattractant
concentration. Surprisingly, these two criteria conflict, so that when one
performance criterion is most favorable, the other is unfavorable. Since both
types of behavior are biologically relevant, we include both behaviors in a
composite optimization that yields a response function that closely resembles
experimental measurements. Our work suggests that the bacterial chemotactic
response function can be derived from simple behavioral considerations, and
sheds light on how the response function contributes to chemotactic
performance.


  Summary: Most cellular tasks are performed not by individual proteins, but by
groups of functionally associated proteins, often referred to as modules. In a
protein assocation network modules appear as groups of densely interconnected
nodes, also called communities or clusters. These modules often overlap with
each other and form a network of their own, in which nodes (links) represent
the modules (overlaps). We introduce CFinder, a fast program locating and
visualizing overlapping, densely interconnected groups of nodes in undirected
graphs, and allowing the user to easily navigate between the original graph and
the web of these groups. We show that in gene (protein) association networks
CFinder can be used to predict the function(s) of a single protein and to
discover novel modules. CFinder is also very efficient for locating the cliques
of large sparse graphs.
  Availability: CFinder (for Windows, Linux, and Macintosh) and its manual can
be downloaded from http://angel.elte.hu/clustering.
  Contact: cfinder@angel.elte.hu


  Among features of real immune responses which occur when antigens invade a
body,there are two remarkable features. One is that the amount of antibodies
produced in the secondary invasion by the same antigens is more than 10 times
larger than that in the primary invasion. The other is that more effective
antibodies which can neutralize the antigens more quickly are produced by
somatic hypermutation during the immune response. This phenomenon is named as
'affinity maturation'. In this paper, we try to reproduce these features by
dynamical system models and present possible factors to realize them. Further,
we present a model in which the memory of the invasion by antigens is realized
without immune memory cells.


  The chiral nature of DNA plays a crucial role in cellular processes. Here we
use magnetic tweezers to explore one of the signatures of this chirality, the
coupling between stretch and twist deformations. We show that the extension of
a stretched DNA molecule increases linearly by 0.42 nm per excess turn applied
to the double helix. This result contradicts the intuition that DNA should
lengthen as it is unwound and get shorter with overwinding. We then present
numerical results of energy minimizations of torsionally restrained DNA that
display a behaviour similar to the experimental data and shed light on the
molecular details of this surprising effect.


  Based on previously published experimental observations and mathematical
models for Hes1, p53 and NF-kB gene expression, we improve these models through
a distributed delay formulation of the time lag between transcription factor
binding and mRNA production. This description of natural variability for delays
introduces a transition from a stable steady state to limit cycle oscillations
and then a second transition back to a stable steady state which has not been
observed in previously published models. On the basis of our results and
following recent discussions about the role of delay-induced oscillations in
gene transcription we establish the hypothesis that the period of the
delay-induced cyclic changes should be characterized by an upper bound so that
it cannot be greater then the period of fundamental biological cycles. We
demonstrate our approach for two models. The first model describes Hes1
autorepression with equations for Hes1 mRNA production and Hes1 protein
translation. The second model describes Hes1 repression by the protein complex
Gro/TLE1/Hes1, where Gro/TLE1 is activated by Hes1 phosphorylation.


  An inchworm processive mechanism is proposed to explain the motion of dimeric
molecular motors such as kinesin. We present here preliminary results for this
mechanism focusing on observables like mean velocity, coupling ratio and
efficiency versus ATP concentration and the external load F.


  The in vivo free energy of looping double-stranded DNA by the lac repressor
has a remarkable behavior whose origins are not fully understood. In addition
to the intrinsic periodicity of the DNA double helix, the in vivo free energy
has an oscillatory component of about half the helical period and oscillates
asymmetrically with an amplitude significantly smaller than predicted by
current theories. Here, we show that the in vivo behavior is accurately
accounted for by the simultaneous presence of two distinct conformations of
looped DNA. Our analysis reveals that these two conformations have different
optimal free energies and phases and that they behave distinctly in the
presence of key architectural proteins.


  Protein folding is a problem of large interest since it concerns the
mechanism by which the genetic information is translated into proteins with
well defined three-dimensional (3D) structures and functions. Recently
theoretical models have been developed to predict the protein folding rate
considering the relationships of the process with tolopological parameters
derived from the native (atomic-solved) protein structures. Previous works
classified proteins in two different groups exhibiting either a
single-exponential or a multi-exponential folding kinetics. It is well known
that these two classes of proteins are related to different protein structural
features. The increasing number of available experimental kinetic data allows
the application to the problem of a machine learning approach, in order to
predict the kinetic order of the folding process starting from the experimental
data so far collected. This information can be used to improve the prediction
of the folding rate. In this work first we describe a support vector
machine-based method (SVM-KO) to predict for a given protein the kinetic order
of the folding process. Using this method we can classify correctly 78% of the
folding mechanisms over a set of 63 experimental data. Secondly we focus on the
prediction of the logarithm of the folding rate. This value can be obtained as
a linear regression task with a SVM-based method. In this paper we show that
linear correlation of the predicted with experimental data can improve when the
regression task is computed over two different sets, instead of one, each of
them composed by the proteins with a correctly predicted two state or
multistate kinetic order.


  We investigate a stochastic version of a simple enzymatic reaction which
follows the generic Michaelis-Menten kinetics. At sufficiently high
concentrations of reacting species, the molecular fluctuations can be
approximated as a realization of a Brownian dynamics for which the model
reaction kinetics takes on the form of a stochastic differential equation.
After eliminating a fast kinetics, the model can be rephrased into a form of a
one-dimensional overdamped Langevin equation. We discuss physical aspects of
environmental noises acting in such a reduced system, pointing out the
possibility of coexistence of dynamical regimes where noise-enhanced stability
and resonant activation phenomena can be observed together.


  We study a Langevin equation derived from the Michaelis-Menten (MM)
phenomenological scheme for catalysis accompanying a spontaneous replication of
molecules, which may serve as a simple model of cell-mediated immune
surveillance against cancer. We examine how two different and statistically
independent sources of noise - dichotomous multiplicative noise and additive
Gaussian white noise - influence the population's extinction time. This
quantity is identified as the mean first passage time of the system across the
zero population state. We observe the effects of resonant activation (RA) and
noise-enhanced stability (NES) and we report the evidence for competitive
co-occurrence of both phenomena in a given regime of noise parameters. We
discuss the statistics of first passage times in this regime and the role of
different pseudo-potential profiles on the RA and NES phenomena. The RA/NES
coexistence region brings an interesting interpretation for the growth kinetics
of cancer cells population, as the NES effect enhancing the stability of the
tumoral state becomes strongly reduced by the RA phenomenon.


  We study a generalised model of population growth in which the state variable
is population growth rate instead of population size. Stochastic parametric
perturbations, modelling phenotypic variability, lead to a Langevin system with
two sources of multiplicative noise. The stationary probability distributions
have two characteristic power-law scales. Numerical simulations show that noise
suppresses the explosion of the growth rate which occurs in the deterministic
counterpart. Instead, in different parameter regimes populations will grow with
``anomalous'' stochastic rates and (i) stabilise at ``random carrying
capacities'', or (ii) go extinct in random times. Using logistic fits to
reconstruct the simulated data, we find that even highly significant
estimations do not recover or reflect information about the deterministic part
of the process. Therefore, the logistic interpretation is not biologically
meaningful. These results have implications for distinct model-aided
calculations in biological situations because these kinds of estimations could
lead to spurious conclusions.


  We investigate the sequence-dependent behaviour of localised excitations in a
toy, nonlinear model of DNA base-pair opening originally proposed by Salerno.
Specifically we ask whether ``breather'' solitons could play a role in the
facilitated location of promoters by RNA polymerase. In an effective potential
formalism, we find excellent correlation between potential minima and {\em
Escherichia coli} promoter recognition sites in the T7 bacteriophage genome.
Evidence for a similar relationship between phage promoters and downstream
coding regions is found and alternative reasons for links between AT richness
and transcriptionally-significant sites are discussed. Consideration of the
soliton energy of translocation provides a novel dynamical picture of sliding:
steep potential gradients correspond to deterministic motion, while ``flat''
regions, corresponding to homogeneous AT or GC content, are governed by random,
thermal motion. Finally we demonstrate an interesting equivalence between
planar, breather solitons and the helical motion of a sliding protein
``particle'' about a bent DNA axis.


  Existing theories for the evolution of aging and death treat senescence as a
side-effect of strong selection for fertility. These theories are
well-developed mathematically, but fit poorly with emerging experimental data.
The data suggest that aging is an adaptation, selected for its own sake. But
aging contributes only negatively to fitness of the individual. What kind of
selection model would permit aging to emerge as a population-level adaptation?
I explore the thesis that population dynamics is inherently chaotic, and that
aging is selected for its role in smoothing demographic fluctuations. The
logistic equation provides a natural vehicle for this model because it has
played a central role in two sciences: Population growth in a resource-limited
niche has long been modeled by the differential LE; and, as a difference
equation, the LE is a canonical example of the emergence of chaos. Suppose that
feedback about depleted resources generally arrives too late to avoid a wave of
unsupportable population growth; then logistic population dynamics is subject
to chaotic fluctuations. It is my thesis that aging is an evolutionary
adaptation selected for its stabilizing effect on chaotic population dynamics.


  A general approach to the design of accurate classical potentials for protein
folding is described. It includes the introduction of a meaningful statistical
measure of the differences between approximations of the same potential energy,
the definition of a set of Systematic and Approximately Separable and Modular
Internal Coordinates (SASMIC), much convenient for the simulation of general
branched molecules, and the imposition of constraints on the most rapidly
oscillating degrees of freedom. All these tools are used to study the effects
of constraints in the Conformational Equilibrium Distribution (CED) of the
model dipeptide HCO-L-Ala-NH2. We use ab initio Quantum Mechanics calculations
including electron correlation at the MP2 level to describe the system, and we
measure the conformational dependence of the correcting terms to the naive CED
based in the Potential Energy Surface (PES) without any simplifying assumption.
These terms are related to mass-metric tensors determinants and also occur in
the Fixman's compensating potential. We show that some of the corrections are
non-negligible if one is interested in the whole Ramachandran space. On the
other hand, if only the energetically lower region, containing the principal
secondary structure elements, is assumed to be relevant, then, all correcting
terms may be neglected up to peptides of considerable length. This is the first
time, as far as we know, that the analysis of the conformational dependence of
these correcting terms is performed in a relevant biomolecule with a realistic
potential energy function.


  For most organisms with viscous population structure, spatially localized
growth drives the invasive advance of a favorable mutation. We model a
two-allele competition where recurrent mutation introduces a genotype with a
rate of local propagation exceeding the resident's rate. We capture
ecologically important properties of the rare invader's stochastic dynamics by
assuming discrete individuals and neighborhood interactions. To understand how
individual-level processes may govern population patterns, we invoke the
physical theory for nucleation of spatial systems. Nucleation theory
discriminates between single-cluster and multi-cluster dynamics. A sufficiently
low mutation rate, or a sufficiently small environment, generates
single-cluster dynamics, an inherently stochastic process; a favorable mutation
advances only if the invader cluster reaches a critical radius. For this mode
of invasion we identify the probability distribution of waiting times until the
favored allele advances to competitive dominance, and we ask how the critical
cluster size varies as propagation or mortality rates vary. Increasing the
mutation rate or system size generates multi-cluster invasion, where spatial
averaging produces nearly deterministic global dynamics. For this process, an
analytical approximation from nucleation theory, called Avrami's Law, describes
the time-dependent behavior of the genotype densities with remarkable accuracy.


  In the process of translation, ribosomes read the genetic code on an mRNA and
assemble the corresponding polypeptide chain. The ribosomes perform discrete
directed motion which is well modeled by a totally asymmetric simple exclusion
process (TASEP) with open boundaries. Using Monte Carlo simulations and a
simple mean-field theory, we discuss the effect of one or two ``bottlenecks''
(i.e., slow codons) on the production rate of the final protein. Confirming and
extending previous work by Chou and Lakatos, we find that the location and
spacing of the slow codons can affect the production rate quite dramatically.
In particular, we observe a novel ``edge'' effect, i.e., an interaction of a
single slow codon with the system boundary. We focus in detail on ribosome
density profiles and provide a simple explanation for the length scale which
controls the range of these interactions.


  In many models of genotypic evolution, the vector of genotype populations
satisfies a system of linear ordinary differential equations. This system of
equations models a competition between differential replication rates (fitness)
and mutation. Mutation operates as a generalized diffusion process on genotype
space. In the large time asymptotics, the replication term tends to produce a
single dominant quasispecies, unless the mutation rate is too high, in which
case the populations of different genotypes becomes de-localized. We introduce
a more macroscopic picture of genotypic evolution wherein a random replication
term in the linear model displays features analogous to Anderson localization.
When coupled with non-linearities that limit the population of any given
genotype, we obtain a model whose large time asymptotics display stable
genotypic diversity


  During the HIV infection several quasispecies of the virus arise, which are
able to use different coreceptors, in particular the CCR5 and CXCR4 coreceptors
(R5 and X4 phenotypes, respectively). The switch in coreceptor usage has been
correlated with a faster progression of the disease to the AIDS phase. As
several pharmaceutical companies are starting large phase III trials for R5 and
X4 drugs, models are needed to predict the co-evolutionary and competitive
dynamics of virus strains. We present a model of HIV early infection which
describes the dynamics of R5 quasispecies and a model of HIV late infection
which describes the R5 to X4 switch. We report the following findings: after
superinfection or coinfection, quasispecies dynamics has time scales of several
months and becomes even slower at low number of CD4+ T cells. The curve of CD4+
T cells decreases, during AIDS late stage, and can be described taking into
account the X4 related Tumor Necrosis Factor dynamics. Phylogenetic inference
of chemokine receptors suggests that viral mutational pathway may generate R5
variants able to interact with chemokine receptors different from CXCR4. This
may explain the massive signaling disruptions in the immune system observed
during AIDS late stages and may have relevance for vaccination and therapy.


  PCR (Polymerase Chain Reaction), a method which replicates a selected
sequence of DNA, has revolutionized the study of genomic material, but
mathematical study of the process has been limited to simple deterministic
models or descriptions relying on stochastic processes. In this paper we
develop a suite of deterministic models for the reactions of quantitative PCR
(Polymerase Chain Reaction) based on the law of mass action. Maps are created
from DNA copy number in one cycle to the next, with ordinary differential
equations describing the evolution of difference molecular species during each
cycle. Qualitative analysis is preformed at each stage and parameters are
estimated by fitting each model to data from Roche LightCycler (TM) runs.


  The problem of the rate and mechanisms of biological evolution was
considered. It was shown that species could not be formed due to undirected
mutations in characteristic times of about one million years. A mechanism of
deterministic molecular evolution assuming a directed change of the genome was
proposed.


  Experimental evidence suggests that epidermal growth factor receptor
(EGFR)-mediated activation of the signaling protein phospholipase C gamma plays
a critical role in a cancer cell's phenotypic decision to either proliferate or
to migrate at a given point in time. Here, we present a novel three-dimensional
multiscale agent-based model to simulate this cellular decision process in the
context of a virtual brain tumor. Each tumor cell is equipped with an EGFR
gene-protein interaction network module that also connects to a simplified cell
cycle description. The simulation results show that over time proliferative and
migratory cell populations not only oscillate but also directly impact the
spatio-temporal expansion patterns of the entire cancer system. The percentage
change in the concentration of the sub-cellular interaction network's molecular
components fluctuates, and, for the proliferation-to-migration switch we find
that the phenotype triggering molecular profile to some degree varies as the
tumor system grows and the microenvironment changes. We discuss potential
implications of these findings for experimental and clinical cancer research.


  Infectious diseases and computer malwares spread among humans and computers
through the network of contacts among them. These networks are characterized by
wide connectivity fluctuations, connectivity correlations and the small-world
property. In a previous work [A. Vazquez, Phys. Rev. Lett. 96, 038702 (2006)] I
have shown that the connectivity fluctuations together with the small-world
property lead to a novel spreading law, characterized by an initial power law
growth with an exponent determined by the average node distance on the network.
Here I extend these results to consider the influence of connectivity
correlations which are generally observed in real networks. I show that
assortative and disassortative connectivity correlations enhance and diminish,
respectively, the range of validity of this spreading law. As a corollary I
obtain the region of connectivity fluctuations and degree correlations
characterized by the absence of an epidemic threshold. These results are
relevant for the spreading of infectious diseases, rumors, and information
among humans and the spreading of computer viruses, email worms and hoaxes
among computer users.


  Protein Fragment Motif Finder (PFMFind) is a system that enables efficient
discovery of relationships between short fragments of protein sequences using
similarity search. It supports queries based on score matrices and PSSMs
obtained through an iterative procedure similar to PSI-BLAST. PSSM construction
is customisable through plugins written in Python. PFMFind consists of a GUI
client, an algorithm using an index for fast similarity search and a relational
database for storing search results and sequence annotations. It is written
mostly in Python. All components communicate between themselves using TCP/IP
sockets and can be located on different physical machines.


  A chemical kinetic model of the elongation dynamics of RNA polymerase along a
DNA sequence is introduced. The proposed model governs the discrete movement of
the RNA polymerase along a DNA template, with no consideration given to elastic
effects. The model's novel concept is a ``look-ahead'' feature, in which
nucleotides bind reversibly to the DNA prior to being incorporated covalently
into the nascent RNA chain. Results are presented for specific DNA sequences
that have been used in single-molecule experiments of the random walk of RNA
polymerase along DNA. By replicating the data analysis algorithm from the
experimental procedure, the model produces velocity histograms, enabling direct
comparison with these published results.


  We consider an individual-based two-dimensional spatial model with
nearest-neighbor preemptive competition to study front propagation between an
invader and a resident species. In particular, we investigate the asymptotic
front velocity and compare it with mean-field predictions.


  One of the main tasks of post-genomic informatics is to systematically
investigate all molecules and their interactions within a living cell so as to
understand how these molecules and the interactions between them relate to the
function of the organism, while networks are appropriate abstract description
of all kinds of interactions. In the past few years, great achievement has been
made in developing theory of complex networks for revealing the organizing
principles that govern the formation and evolution of various complex
biological, technological and social networks. This paper reviews the
accomplishments in constructing genome-based metabolic networks and describes
how the theory of complex networks is applied to analyze metabolic networks.


  Statistical analysis of protein-protein interactions shows anomalously high
frequency of homodimers [Ispolatov, I., et al. (2005) Nucleic Acids Res 33,
3629-35]. Furthermore, recent findings [Wright, C.F., et al. (2005) Nature 438,
878-81] demonstrate that maintaining low sequence identity is a key
evolutionary mechanism that inhibits protein aggregation. Here, we study
statistical properties of interacting protein-like surfaces and predict the
effect of universal, enhanced self-attraction of proteins. The effect
originates in the fact that a pattern self-match between two identical, even
randomly organized interacting protein surfaces is always stronger compared to
the pattern match between two different, promiscuous protein surfaces. This
finding implies an increased probability of homodimer selection in the course
of early evolution. Our simple model of early evolutionary selection of
interacting proteins accurately reproduces the experimental data on homodimer
interface aminoacid compositions. In addition, we predict that heterodimers
evolved from homodimers with the negative design evolutionary pressure applied
against promiscuous homodimer formation. We predict that the anti-homodimer
negative design evolutionary signal is conveyed through the enrichment of
heterodimeric interfaces in polar residues, and most profoundly in glutamic
acid and lysine, which is consistent with experimental findings. We predict
therefore that the negative design against homodimers is the


  Kinetics of folding of a protein held in a force-clamp are compared to an
unconstrained folding. The comparison is made within a simple topology-based
dynamical model of ubiquitin. We demonstrate that the experimentally observed
variations in the end-to-end distance reflect microscopic events during
folding. However, the folding scenarios in and out of the force-clamp are
distinct.


  There is a correspondence between the circulation of blood in all higher
animals and the circulation of sap in all higher plants - up to heights h of
140 m - through the xylem and phloem vessels. Plants suck in water from the
soil, osmotically through the roothair zone, and subsequently lift it
osmotically again, and by capillary suction (via their buds, leaves, and
fruits) into their crowns. In between happens a reverse osmosis - the
endodermis jump - realized by two layers of subcellular mechanical pumps in the
endodermis walls which are powered by ATP, or in addition by two analogous
layers of such pumps in the exodermis. The thus established root pressure helps
forcing the absorbed ground water upward, through the whole plant, and often
out again, in the form of guttation, or exudation.


  DNA looping plays a fundamental role in a wide variety of biological
processes, providing the backbone for long range interactions on DNA. Here we
develop the first model for DNA looping by an arbitrarily large number of
proteins and solve it analytically in the case of identical binding. We uncover
a switch-like transition between looped and unlooped phases and identify the
key parameters that control this transition. Our results establish the basis
for the quantitative understanding of fundamental cellular processes like DNA
recombination, gene silencing, and telomere maintenance.


  We discuss several qualitative properties of the solutions of
reaction-diffusion systems and equations of the form $u_t = \epsilon^2 D \Delta
u + f(u,x,\epsilon t)$, that are used in modeling pattern formation. We analyze
the diffusion neutral and the diffusion dependent situations that, in the time
autonomous case, are distinguished by considering the attractors of the shorted
equation $u_t = f(u,x)$. We discuss the consequences of being in one or in the
other of the two situations and present examples from developmental biology and
from fluid mechanics.


  We investigate aggregation mechanism of two proteins in a thermodynamically
unambiguous manner by considering the finite size effect of free energy
landscape of HP lattice protein model. Multi-Self-Overlap-Ensemble Monte Carlo
method is used for numerical calculations. We find that a dimer can be formed
spontaneously as a thermodynamically stable state when the system is small
enough. It implies the possibility that the aggregation of proteins in a cell
is triggered when they are confined in a small region by, for example, being
surrounded by other macromolecules.We also find that the dimer exhibits a
transition between unstable state and metastable state in the infinite system.


  Transforming Growth Factor-beta (TGF-beta) signalling is an important
regulator of cellular growth and differentiation. The principal intracellular
mediators of TGF-beta signalling are the Smad proteins, which upon TGF-beta
stimulation accumulate in the nucleus and regulate transcription of target
genes. To investigate the mechanisms of Smad nuclear accumulation, we developed
a simple mathematical model of canonical Smad signalling. The model was built
using both published data and our experimentally determined cellular Smad
concentrations (isoforms 2, 3, and 4). We found in mink lung epithelial cells
that Smad2 (8.5-12 x 10^4 molecules/cell) was present in similar amounts to
Smad4 (9.3-12 x 10^4 molecules/cell), while both were in excess of Smad3
(1.1-2.0 x 10^4 molecules/cell). Variation of the model parameters and
statistical analysis showed that Smad nuclear accumulation is most sensitive to
parameters affecting the rates of RSmad phosphorylation and dephosphorylation
and Smad complex formation/dissociation in the nucleus. Deleting Smad4 from the
model revealed that rate-limiting phospho-R-Smad dephosphorylation could be an
important mechanism for Smad nuclear accumulation. Furthermore, we observed
that binding factors constitutively localised to the nucleus do not efficiently
mediate Smad nuclear accumulation if dephosphorylation is rapid. We therefore
conclude that an imbalance in the rates of R-Smad phosphorylation and
dephosphorylation is likely an important mechanism of Smad nuclear accumulation
during TGF-beta signalling.


  Genetic switch systems with mutual repression of two transcription factors
are studied using deterministic and stochastic methods. Numerous studies have
concluded that cooperative binding is a necessary condition for the emergence
of bistability in these systems. Here we show that for a range of biologically
relevant conditions, a suitable combination of network structure and stochastic
effects gives rise to bistability even without cooperative binding.


  Although the importance of protein dynamics in protein function is generally
recognized, the role of protein fluctuations in allosteric effects scarcely has
been considered. To address this gap, the Kullback-Leibler divergence (Dx)
between protein conformational distributions before and after ligand binding
was proposed as a means of quantifying allosteric effects in proteins. Here,
previous applications of Dx to methods for analysis and simulation of proteins
are first reviewed, and their implications for understanding aspects of protein
function and protein evolution are discussed. Next, equations for Dx suggest
that k_{B}TDx should be interpreted as an allosteric free energy -- the free
energy associated with changing the ligand-free protein conformational
distribution to the ligand-bound conformational distribution. This
interpretation leads to a thermodynamic model of allosteric transitions that
unifies existing perspectives on the relation between ligand binding and
changes in protein conformational distributions. The definition of Dx is used
to explore some interesting mathematical relations among commonly recognized
thermodynamic and biophysical quantities, such as the total free energy change
upon ligand binding, and ligand-binding affinities for individual protein
conformations. These results represent the beginnings of a theoretical
framework for considering the full protein conformational distribution in
modeling allosteric transitions. Early applications of the framework have
produced results with implications both for methods for coarsed-grained
modeling of proteins, and for understanding the relation between ligand binding
and protein dynamics.


  Proteins, by virtue of their central role in most biological processes,
represent one of the key subjects of the study of molecular evolution. Inherent
to the indispensability of proteins for living cells is the fact that a given
protein can adopt a specific three-dimensional shape that is specified solely
by the proteins sequence of amino acids. Over the past several decades,
structural biologists have demonstrated that the array of structures that
proteins may adopt is quite astounding, and this has lead to a strong interest
in understanding how protein structures change and evolve over time. In this
review we consider a large body of recent work that attempts to illuminate this
structure-centric picture of protein evolution. Much of this work has focused
on the question of how completely new protein structures (i.e. new folds or
topologies) are discovered by protein sequences as they evolve. Pursuant to
this question of structural innovation has been a desire to describe and
understand the observation that certain types of protein structures are far
more abundant than others and how this uneven distribution of proteins
implicates on the process through which new shapes are discovered. We consider
a number of theoretical models that have been successful at explaining this
heterogeneity in protein populations and discuss the increasing amount of
evidence that indicates that the process of structural evolution involves the
divergence of protein sequences and structures from one another.


  We investigate the macroscopic effects of the ingredients that drive the
origin of species through sympatric speciation. In our model, sympatric
speciation is obtained as we tune up the strength of competition between
individuals with different phenotypes. As a function of this control parameter,
we can characterize, through the behavior of a macroscopic order parameter, a
phase transition from a non-speciation to a speciation state of the system. The
behavior of the first derivative of the order parameter with respect to the
control parameter is consistent with a phase transition and exhibits a sharp
peak at the transition point. For different resources distribution, the
transition point is shifted, an effect similar to pressure in PVT system. The
inverse of the parameter related to sexual selection strength behaves like an
external field in the system and, as thus, is also a control parameter. The
macroscopic effects of the biological parameters used in our model reveal thus
fingerprints typical of thermodynamic quantities in a phase transition of an
equilibrium physical system.


  RNA viruses provide prominent examples of measurably evolving populations. In
HIV infection, the development of drug resistance is of particular interest,
because precise predictions of the outcome of this evolutionary process are a
prerequisite for the rational design of antiretroviral treatment protocols. We
present a mutagenetic tree hidden Markov model for the analysis of longitudinal
clonal sequence data. Using HIV mutation data from clinical trials, we estimate
the order and rate of occurrence of seven amino acid changes that are
associated with resistance to the reverse transcriptase inhibitor efavirenz.


  In the template-assistance model, normal prion protein (PrPC), the pathogenic
cause of prion diseases such as Creutzfeldt-Jakob (CJD) in human, Bovine
Spongiform Encephalopathy (BSE) in cow, and scrapie in sheep, converts to
infectious prion (PrPSc) through an autocatalytic process triggered by a
transient interaction between PrPC and PrPSc. Conventional studies suggest the
S1-H1-S2 region in PrPC to be the template of S1-S2 $\beta$-sheet in PrPSc, and
the conformational conversion of PrPC into PrPSc may involve an unfolding of H1
in PrPC and its refolding into the $\beta$-sheet in PrPSc. Here we conduct a
series of simulation experiments to test the idea of transient interaction of
the template-assistance model. We find that the integrity of H1 in PrPC is
vulnerable to a transient interaction that alters the native dihedral angles at
residue Asn$^{143}$, which connects the S1 flank to H1, but not to interactions
that alter the internal structure of the S1 flank, nor to those that alter the
relative orientation between H1 and the S2 flank.


  The relationship between the shape of a fitness landscape and the underlying
gene interactions, or epistasis, has been extensively studied in the two-locus
case. Gene interactions among multiple loci are usually reduced to two-way
interactions. We present a geometric theory of shapes of fitness landscapes for
multiple loci. A central concept is the genotope, which is the convex hull of
all possible allele frequencies in populations. Triangulations of the genotope
correspond to different shapes of fitness landscapes and reveal all the gene
interactions. The theory is applied to fitness data from HIV and Drosophila
melanogaster. In both cases, our findings refine earlier analyses and reveal
previously undetected gene interactions.


  We analyze the role of the force-dependent kinetics of motor proteins in the
stability of antiparallel arrays of polar filaments, such as those in the
mitotic spindle. We determine the possible stable structures and show that
there exists an instability associated to the collective behavior of motors
that leads to the collapse of the spindle. Our analysis provides a general
framework to understand several experimental observations in eukaryotic cell
division.


  The large-scale shape and function of metabolic networks are intriguing
topics of systems biology. Such networks are on one hand commonly regarded as
modular (i.e. built by a number of relatively independent subsystems), but on
the other hand they are robust in a way not expected of a purely modular
system. To address this question we carefully discuss the partition of
metabolic networks into subnetworks. The practice of preprocessing such
networks by removing the most abundant substrates, "currency metabolites," is
formalized into a network-based algorithm. We study partitions for metabolic
networks of many organisms and find cores of currency metabolites and modular
peripheries of what we call "commodity metabolites." The networks are found to
be more modular than random networks but far from perfectly divisible into
modules. We argue that cross-modular edges are the key for the robustness of
metabolism.


  We describe a mechanism for pronounced biochemical oscillations, relevant to
microscopic systems, such as the intracellular environment. This mechanism
operates for reaction schemes which, when modeled using deterministic rate
equations, fail to exhibit oscillations for any values of rate constants. The
mechanism relies on amplification of the underlying stochasticity of reaction
kinetics within a narrow window of frequencies. This amplification allows
fluctuations to beat the central limit theorem, having a dominant effect even
though the number of molecules in the system is relatively large. The mechanism
is quantitatively studied within simple models of self-regulatory gene
expression, and glycolytic oscillations.


  The Yakushevich model of DNA torsion dynamics supports soliton solutions,
which are supposed to be of special interest for DNA transcription. In the
discussion of the model, one usually adopts the approximation $\ell_0 \to 0$,
where $\ell_0$ is a parameter related to the equilibrium distance between bases
in a Watson-Crick pair. Here we analyze the Yakushevich model without $\ell_0
\to 0$. The model still supports soliton solutions indexed by two winding
numbers $(n,m)$; we discuss in detail the fundamental solitons, corresponding
to winding numbers (1,0) and (0,1) respectively.


  The Yakushevich (Y) model provides a very simple pictures of DNA torsion
dynamics, yet yields remarkably correct predictions on certain physical
characteristics of the dynamics. In the standard Y model, the interaction
between bases of a pair is modelled by a harmonic potential, which becomes
anharmonic when described in terms of the rotation angles; here we substitute
to this different types of improved potentials, providing a more physical
description of the H-bond mediated interactions between the bases. We focus in
particular on soliton solutions; the Y model predicts the correct size of the
nonlinear excitations supposed to model the ``transcription bubbles'', and this
is essentially unchanged with the improved potential. Other features of soliton
dynamics, in particular curvature of soliton field configurations and the
Peierls-Nabarro barrier, are instead significantly changed.


  We study by simulation the effect of the diffusive motion of repressor
molecules on the noise in mRNA and protein levels in the case of a repressed
gene. We find that spatial fluctuations due to diffusion can drastically
enhance the noise in gene expression. For a fixed repressor strength, the noise
due to diffusion can be minimized by increasing the number of repressors or by
decreasing the rate of the open complex formation. We also show that the effect
of spatial fluctuations can be well described by a two-step kinetic scheme,
where formation of an encounter complex by diffusion and the subsequent
association reaction are treated separately. Our results also emphasize that
power spectra are a highly useful tool for studying the propagation of noise
through the different stages of gene expression.


  In this paper, we present a survey of the use of graph theoretical techniques
in Biology. In particular, we discuss recent work on identifying and modelling
the structure of bio-molecular networks, as well as the application of
centrality measures to interaction networks and research on the hierarchical
structure of such networks and network motifs. Work on the link between
structural network properties and dynamics is also described, with emphasis on
synchronization and disease propagation.


  Simple coarse-grained models, such as the Gaussian Network Model, have been
shown to capture some of the features of equilibrium protein dynamics. We
extend this model by using atomic contacts to define residue interactions and
introducing more than one interaction parameter between residues. We use
B-factors from 98 ultra-high resolution X-ray crystal structures to optimize
the interaction parameters. The average correlation between GNM fluctuation
predictions and the B-factors is 0.64 for the data set, consistent with a
previous large-scale study. By separating residue interactions into covalent
and noncovalent, we achieve an average correlation of 0.74, and addition of
ligands and cofactors further improves the correlation to 0.75. However,
further separating the noncovalent interactions into nonpolar, polar, and mixed
yields no significant improvement. The addition of simple chemical information
results in better prediction quality without increasing the size of the
coarse-grained model.


  The spread in time of a mutation through a population is studied analytically
and computationally in fully-connected networks and on spatial lattices. The
time, t_*, for a favourable mutation to dominate scales with population size N
as N^{(D+1)/D} in D-dimensional hypercubic lattices and as N ln N in
fully-connected graphs. It is shown that the surface of the interface between
mutants and non-mutants is crucial in predicting the dynamics of the system.
Network topology has a significant effect on the equilibrium fitness of a
simple population model incorporating multiple mutations and sexual
reproduction. Includes supplementary information.


  A forward diffusion equation describing the evolution of the allele frequency
spectrum is presented. The influx of mutations is accounted for by imposing a
suitable boundary condition. For a Wright-Fisher diffusion with or without
selection and varying population size, the boundary condition is $\lim_{x
\downarrow 0} x f(x,t)=\theta \rho(t)$, where $f(\cdot,t)$ is the frequency
spectrum of derived alleles at independent loci at time $t$ and $\rho(t)$ is
the relative population size at time $t$. When population size and selection
intensity are independent of time, the forward equation is equivalent to the
backwards diffusion usually used to derive the frequency spectrum, but the
forward equation allows computation of the time dependence of the spectrum both
before an equilibrium is attained and when population size and selection
intensity vary with time. From the diffusion equation, we derive a set of
ordinary differential equations for the moments of $f(\cdot,t)$ and express the
expected spectrum of a finite sample in terms of those moments. We illustrate
the use of the forward equation by considering neutral and selected alleles in
a highly simplified model of human history. For example, we show that
approximately 30% of the expected heterozygosity of neutral loci is
attributable to mutations that arose since the onset of population growth in
roughly the last $150,000$ years.


  An intense controversy has surrounded the mechanism of voltage-gating in ion
channels. We interpreted the two leading models of voltage-gating with respect
to the thermodynamic energetics of membrane insertion of the voltage-sensing
'module' from a comprehensive set of potassium channels. KvAP is an archaeal
voltage-gated potassium channel whose x-ray structure was the basis for
determining the general mechanism of voltage-gating. The free energy of
membrane insertion of the KvAP voltage sensor was revealed to be a single
outlier. This was due to its unusual sequence that facilitated large gating
movements in its native lipid membrane. This degree of free energy was the
least typical of the other voltage sensors, including the Shaker potassium
channel. We inferred that the two leading models of voltage-gating referred to
alternative mechanisms of voltage-gating: each is applicable to an independent
set of ion channels. The large motion of the voltage-sensor during gating
proposed by the KvAP-paddle model of gating is unlikely to be mirrored by the
majority of ion channels whose voltage sensors are not located at the
membrane-cytoplasm interface in the channel closed state.


  In the yeast Saccharomyces cerevisiae, the interplay between galactose,
Gal3p, Gal80p and Gal4p determines the transcriptional status of the genes
required for galactose utilization. After an increase in galactose
concentration, galactose molecules bind onto Gal3p. This event leads via Gal80p
to the activation of Gal4p, which then induces GAL3 and GAL80 gene
transcription. Here we propose a qualitative dynamic model, whereby these
molecular interaction events represent the first two stages of a functional
feedback loop that closes with the capture of activated Gal4p by newly
synthesized Gal3p and Gal80p, decreasing transcriptional activation and
creating again the protein complex that can bind incoming galactose molecules.
Based on the differential time scales of faster protein interactions versus
slower biosynthetic steps, this feedback loop functions as a derivative filter
where galactose is the input step signal, and released Gal4p is the output
derivative signal. One advantage of such a derivative filter is that GAL genes
are expressed in proportion to the cellular requirement. Furthermore, this
filter adaptively protects the cellular receptors from saturation by galactose,
allowing cells to remain sensitive to variations in galactose concentrations
rather than to absolute concentrations. Finally, this feedback loop, by
allowing phosphorylation of some active Gal4p, may be essential to initiate the
subsequent long-term response.


  Background: One-dimensional protein structures such as secondary structures
or contact numbers are useful for three-dimensional structure prediction and
helpful for intuitive understanding of the sequence-structure relationship.
Accurate prediction methods will serve as a basis for these and other purposes.
Results: We implemented a program CRNPRED which predicts secondary structures,
contact numbers and residue-wise contact orders. This program is based on a
novel machine learning scheme called critical random networks. Unlike most
conventional one-dimensional structure prediction methods which are based on
local windows of an amino acid sequence, CRNPRED takes into account the whole
sequence. CRNPRED achieves, on average per chain, Q3 = 81% for secondary
structure prediction, and correlation coefficients of 0.75 and 0.61 for contact
number and residue-wise contact order predictions, respectively. Conclusion:
CRNPRED will be a useful tool for computational as well as experimental
biologists who need accurate one-dimensional protein structure predictions.


  DNA torsion dynamics is essential in the transcription process; a simple
model for it, in reasonable agreement with experimental observations, has been
proposed by Yakushevich (Y) and developed by several authors; in this, the DNA
subunits made of a nucleoside and the attached nitrogen bases are described by
a single degree of freedom. In this paper we propose and investigate, both
analytically and numerically, a ``composite'' version of the Y model, in which
the nucleoside and the base are described by separate degrees of freedom. The
model proposed here contains as a particular case the Y model and shares with
it many features and results, but represents an improvement from both the
conceptual and the phenomenological point of view. It provides a more realistic
description of DNA and possibly a justification for the use of models which
consider the DNA chain as uniform. It shows that the existence of solitons is a
generic feature of the underlying nonlinear dynamics and is to a large extent
independent of the detailed modelling of DNA. The model we consider supports
solitonic solutions, qualitatively and quantitatively very similar to the Y
solitons, in a fully realistic range of all the physical parameters
characterizing the DNA.


  Affymetrix Genechip microarrays are used widely to determine the simultaneous
expression of genes in a given biological paradigm. Probes on the Genechip
array are atomic entities which by definition are randomly distributed across
the array and in turn govern the gene expression. In the present study, we make
several interesting observations. We show that there is considerable
correlation between the probe intensities across the array which defy the
independence assumption. While the mechanism behind such correlations is
unclear, we show that scaling behavior and the profiles of perfect match (PM)
as well as mismatch (MM) probes are similar and immune to background
subtraction. We believe that the observed correlations are possibly an outcome
of inherent non-stationarities or patchiness in the array devoid of biological
significance. This is demonstrated by inspecting their scaling behavior and
profiles of the PM and MM probe intensities obtained from publicly available
Genechip arrays from three eukaryotic genomes, namely: Drosophila Melanogaster,
Homo Sapiens and Mus musculus across distinct biological paradigms and across
laboratories, with and without background subtraction. The fluctuation
functions were estimated using detrended fluctuation analysis (DFA) with fourth
order polynomial detrending. The results presented in this study provide new
insights into correlation signatures of PM and MM probe intensities and
suggests the choice of DFA as a tool for qualitative assessment of Affymetrix
Genechip microarrays prior to their analysis. A more detailed investigation is
necessary in order to understand the source of these correlations.


  We provide a mean-field analysis of community structure of social and
biological networks assuming that actors are able to evaluate some tree-derived
distance to the other actors and tend to aggregate with the less distant. We
show that such networks have small components, and give exact descriptions for
the probability distribution of a typical community size and the number of
communities. In particular, we show that the probability distribution of the
community size is well-approximated by a power-law distribution with exponent
two. We illustrate the robustness of the mean-field analysis by comparing its
predictions on previously studied social networks and biological data.


  We advocate the use of qualitative models in the analysis of large biological
systems. We show how qualitative models are linked to theoretical differential
models and practical graphical models of biological networks. A new technique
for analyzing qualitative models is introduced, which is based on an efficient
representation of qualitative systems. As shown through several applications,
this representation is a relevant tool for the understanding and testing of
large and complex biological networks.


  To confer high specificity and affinity in binding, contacts at interfaces
between two interacting macromolecules are expected to exhibit pair preferences
for types of atoms or residues. Here we quantify these preferences by measuring
the mutual information of contacts for 895 protein-protein interfaces. The
information content is significant and is highest at the atomic resolution. A
simple phenomenological theory reveals a connection between information at
interfaces and the free energy spectrum of association. The connection is
presented in the form of a relation between mutual information and the energy
gap of the native bound state to off-target bound states. Measurement of
information content in designed lattice interfaces show the predicted scaling
behavior to the energy gap. Our theory also suggests that mutual information in
contacts emerges by a selection mechanism, and that strong selection, or high
conservation, of residues should lead to correspondingly high mutual
information. Amino acids which contribute more heavily to information content
are then expected to be more conserved. We verify this by showing a
statistically significant correlation between the conservation of each of the
twenty amino acids and their individual contribution to the information content
at protein-protein interfaces


  A novel theoretical model of dynamic instability of a system of linear (1D)
microtubules (MTs) in a bounded domain is introduced for studying the role of a
cell edge in vivo and analyzing the effect of competition for a limited amount
of tubulin. The model differs from earlier models in that the evolution of MTs
is based on the rates of single unit (e.g., a heterodimer per protofilament)
transformations, in contrast to postulating effective rates/frequencies of
larger-scale changes, extracted, e.g., from the length history plots of MTs.
Spontaneous GTP hydrolysis with finite rate after polymerization is assumed,
and theoretical estimates of an effective catastrophe frequency as well as
other parameters characterizing MT length distributions and cap size are
derived. We implement a simple cap model which does not include vectorial
hydrolysis. We demonstrate that our theoretical predictions, such as steady
state concentration of free tubulin, and parameters of MT length distributions,
are in agreement with the numerical simulations. The present model establishes
a quantitative link between microscopic parameters governing the dynamics of
MTs and macroscopic characteristics of MTs in a closed system. Lastly, we use a
computational Monte Carlo model to provide an explanation for non-exponential
MT length distributions observed in experiments. In particular, we show that
appearance of such non-exponential distributions in the experiments can occur
because the true steady state has not been reached, and/or due to the presence
of a cell edge.


  The search for similarity and dissimilarity measures on phylogenetic trees
has been motivated by the computation of consensus trees, the search by
similarity in phylogenetic databases, and the assessment of clustering results
in bioinformatics. The transposition distance for fully resolved phylogenetic
trees is a recent addition to the extensive collection of available metrics for
comparing phylogenetic trees. In this paper, we generalize the transposition
distance from fully resolved to arbitrary phylogenetic trees, through a
construction that involves an embedding of the set of phylogenetic trees with a
fixed number of labeled leaves into a symmetric group and a generalization of
Reidys-Stadler's involution metric for RNA contact structures. We also present
simple linear-time algorithms for computing it.


  Complex interactions between genes or proteins contribute a substantial part
to phenotypic evolution. Here we develop an evolutionarily grounded method for
the cross-species analysis of interaction networks by {\em alignment}, which
maps bona fide functional relationships between genes in different organisms.
Network alignment is based on a scoring function measuring mutual similarities
between networks taking into account their interaction patterns as well as
sequence similarities between their nodes. High-scoring alignments and optimal
alignment parameters are inferred by a systematic Bayesian analysis. We apply
this method to analyze the evolution of co-expression networks between human
and mouse. We find evidence for significant conservation of gene expression
clusters and give network-based predictions of gene function. We discuss
examples where cross-species functional relationships between genes do not
concur with sequence similarity.


  It was first suggested by Englander et al to model the nonlinear dynamics of
DNA relevant to the transcription process in terms of a chain of coupled
pendulums. In a related paper [q-bio.BM/0604014] we argued for the advantages
of an extension of this approach based on considering a chain of double
pendulums with certain characteristics. Here we study a simplified model of
this kind, focusing on its general features and nonlinear travelling wave
excitations; in particular, we show that some of the degrees of freedom are
actually slaved to others, allowing for an effective reduction of the relevant
equations.


  The molecular network in an organism consists of transcription/translation
regulation, protein-protein interactions/modifications and a metabolic network,
together forming a system that allows the cell to respond sensibly to the
multiple signal molecules that exist in its environment. A key part of this
overall system of molecular regulation is therefore the interface between the
genetic and the metabolic network. A motif that occurs very often at this
interface is a negative feedback loop used to regulate the level of the signal
molecules. In this work we use mathematical models to investigate the steady
state and dynamical behaviour of different negative feedback loops. We show, in
particular, that feedback loops where the signal molecule does not cause the
dissociation of the transcription factor from the DNA respond faster than loops
where the molecule acts by sequestering transcription factors off the DNA. We
use three examples, the bet, mer and lac systems in E. coli, to illustrate the
behaviour of such feedback loops.


  We present a stochastic approach to modeling the dynamics of coexistence of
prey and predator populations. It is assumed that the space of coexistence is
explicitly subdivided in a grid of cells. Each cell can be occupied by only one
individual of each species or can be empty. The system evolves in time
according to a probabilistic cellular automaton composed by a set of local
rules which describe interactions between species individuals and mimic the
process of birth, death and predation. By performing computational simulations,
we found that, depending on the values of the parameters of the model, the
following states can be reached: a prey absorbing state and active states of
two types. In one of them both species coexist in a stationary regime with
population densities constant in time. The other kind of active state is
characterized by local coupled time oscillations of prey and predator
populations. We focus on the self-organized structures arising from
spatio-temporal dynamics of the coexistence. We identify distinct spatial
patterns of prey and predators and verify that they are intimally connected to
the time coexistence behavior of the species. The occurrence of a prey
percolating cluster on the spatial patterns of the active states is also
examined.