  The inference of the interactions between organisms in an ecosystem from
observational data is an important problem in ecology. This paper presents a
mathematical inference method, originally developed for the inference of
biochemical networks in molecular biology, adapted for the inference of
networks of ecological interactions. The method is applied to a network of
invertebrate families (taxa) in a rice field.


  We show that biological networks with serial regulation (each node regulated
by at most one other node) are constrained to {\it direct functionality}, in
which the sign of the effect of an environmental input on a target species
depends only on the direct path from the input to the target, even when there
is a feedback loop allowing for multiple interaction pathways. Using a
stochastic model for a set of small transcriptional regulatory networks that
have been studied experimentally, we further find that all networks can achieve
all functions permitted by this constraint under reasonable settings of
biochemical parameters. This underscores the functional versatility of the
networks.


  A phylogenetic network is a directed acyclic graph that visualises an
evolutionary history containing so-called reticulations such as recombinations,
hybridisations or lateral gene transfers. Here we consider the construction of
a simplest possible phylogenetic network consistent with an input set T, where
T contains at least one phylogenetic tree on three leaves (a triplet) for each
combination of three taxa. To quantify the complexity of a network we consider
both the total number of reticulations and the number of reticulations per
biconnected component, called the level of the network. We give polynomial-time
algorithms for constructing a level-1 respectively a level-2 network that
contains a minimum number of reticulations and is consistent with T (if such a
network exists). In addition, we show that if T is precisely equal to the set
of triplets consistent with some network, then we can construct such a network
with smallest possible level in time O(|T|^(k+1)), if k is a fixed upper bound
on the level of the network.


  In the evolutionary Prisoner's Dilemma (PD) game, agents play with each other
and update their strategies in every generation according to some microscopic
dynamical rule. In its spatial version, agents do not play with every other
but, instead, interact only with their neighbors, thus mimicking the existing
of a social or contact network that defines who interacts with whom. In this
work, we explore evolutionary, spatial PD systems consisting of two types of
agents, each with a certain update (reproduction, learning) rule. We
investigate two different scenarios: in the first case, update rules remain
fixed for the entire evolution of the system; in the second case, agents update
both strategy and update rule in every generation. We show that in a well-mixed
population the evolutionary outcome is always full defection. We subsequently
focus on two-strategy competition with nearest-neighbor interactions on the
contact network and synchronized update of strategies. Our results show that,
for an important range of the parameters of the game, the final state of the
system is largely different from that arising from the usual setup of a single,
fixed dynamical rule. Furthermore, the results are also very different if
update rules are fixed or evolve with the strategies. In these respect, we have
studied representative update rules, finding that some of them may become
extinct while others prevail. We describe the new and rich variety of final
outcomes that arise from this co-evolutionary dynamics. We include examples of
other neighborhoods and asynchronous updating that confirm the robustness of
our conclusions. Our results pave the way to an evolutionary rationale for
modelling social interactions through game theory with a preferred set of
update rules.


  We present a comparative analysis of large-scale topological and evolutionary
properties of transcription networks in three species, the two distant bacteria
E. coli and B. subtilis, and the yeast S. cerevisiae. The study focuses on the
global aspects of feedback and hierarchy in transcriptional regulatory
pathways. While confirming that gene duplication has a significant impact on
the shaping of all the analyzed transcription networks, our results point to
distinct trends between the bacteria, where time constraints in the
transcription of downstream genes might be important in shaping the
hierarchical structure of the network, and yeast, which seems able to sustain a
higher wiring complexity, that includes the more feedback, intricate hierarchy,
and the combinatorial use of heterodimers made of duplicate transcription
factors.


  We establish a log-supermodularity property for probability distributions on
binary patterns observed at the tips of a tree that are generated under any
2--state Markov process. We illustrate the applicability of this result in
phylogenetics by deriving an inequality relevant to estimating expected future
phylogenetic diversity under a model of species extinction. In a further
application of the log-supermodularity property, we derive a purely
combinatorial inequality for the parsimony score of a binary character. The
proofs of our results exploit two classical theorems in the combinatorics of
finite sets.


  A new paradigm emerging in the description of magnetic-sensitive
radical-ion-pair recombination reactions is presented. This paradigm is founded
on the realization that the recombination process of radical-ion pairs is a
continuous quantum measurement. In the regime of low magnetic fields we
describe the appearance of purely quantum phenomena, that fade away as the
magnetic field increases. We will analyze the magnetic sensitivity of these
reactions under this perspective and bridge the new full quantum theory with
the existing classical reaction theory applicable at high magnetic fields.
Based on the quantum theory of recombination we will then explain experimental
observations incompatible with classical reaction theory, in particular the
effect of deuteration on the magnetic sensitivity of radical-ion pair
recombination yields.


  The Ca-sensitive regulatory switch of cardiac muscle is a paradigmatic
example of protein assemblies that communicate ligand binding through
allosteric change. The switch is a dimeric complex of troponin C (TnC), an
allosteric sensor for Ca, and troponin I (TnI), an allosteric reporter.
Time-resolved equilibrium FRET measurements suggest that the switch activates
in two steps: a TnI-independent Ca-priming step followed by TnI-dependent
opening. To resolve the mechanistic role of TnI in activation we performed
stopped-flow FRET measurements of activation following rapid addition of a
lacking component (Ca or TnI) and deactivation following rapid chelation of Ca.
The time-resolved measurements, stopped-flow measurements, and Ca-titration
measurements were globally analyzed in terms of a new quantitative dynamic
model of TnC-TnI allostery. The analysis provided a mesoscopic parameterization
of distance changes, free energy changes, and transition rates among the
accessible coarse-grained states of the system. The results reveal (i) the
Ca-induced priming step, which precedes opening, is the rate limiting step in
activation, (ii) closing is the rate limiting step in deactivation, (iii) TnI
induces opening, (iv) an incompletely deactivated population when regulatory Ca
is not bound, which generates an accessory pathway of activation, and (v)
incomplete activation by Ca--when regulatory Ca is bound, a 3:2 mixture of
dynamically inter-converting open (active) and primed-closed (partially active)
conformers is observed (15 C). Temperature-dependent stopped-flow FRET
experiments provide a near complete thermo-kinetic parametrization of opening.
<Abstract Truncated>


  The unique Hamiltonian description of neuro- and psycho-dynamics at the
macroscopic, classical, inter-neuronal level of brain's neural networks, and
microscopic, quantum, intra-neuronal level of brain's microtubules, is
presented in the form of open Liouville equations. This implies the arrow of
time in both neuro- and psycho-dynamic processes and shows the existence of the
formal neuro-biological space-time self-similarity.
  Keywords: Neuro- and psycho-dynamics, Brain microtubules, Hamiltonian and
Liouville dynamics, Neuro-biological self-similarity


  We study the effect of permeabilizing electric fields applied to two
different types of giant unilamellar vesicles, the first formed from EggPC
lipids and the second formed from DOPC lipids. Experiments on vesicles of both
lipid types show a decrease in vesicle radius which is interpreted as being due
to lipid loss during the permeabilization process. We show that the decrease in
size can be qualitatively explained as a loss of lipid area which is
proportional to the area of the vesicle which is permeabilized. Three possible
mechanisms responsible for lipid loss were directly observed: pore formation,
vesicle formation and tubule formation.


  We describe a novel approach to the analysis of toxicogenomics data and
elucidation of biological networks affected by drug treatments. In this method
approximately 15,000 linear pathway modules were generated from manually
assembled pathway maps from MetaCore (GeneGo, Inc.). Microarray expression data
from livers of rat exposed to phenobarbital, mestranol and tamoxifen were
mapped onto these modules. Using different analytical techniques we have
identified sets of "differential" pathways featuring highly correlated
expression among multiple repeats of the same treatment while showing strong
anti-correlation across different treatments. Network modules distinguishing
chemical treatments were re-assembled based on these pathways. Unlike
traditional statistical and clustering procedures in expression profiling, our
method takes into account both network connectivity and gene expression in the
course of the analysis. We demonstrate that it enables identification of
important cellular mechanisms involved in drug response that would have been
missed by the analysis based on individual gene expression profiles.


  How does immune system evolve functional proteins - potent antibodies - in
such a short time? We address this question using a microscopic, protein-level,
sequence-based model of humoral immune response with explicitly defined
interactions between Immunoglobulins, host and pathogen proteins. Potent
Immunoglobulins are discovered in this model via clonal selection and affinity
maturation. Possible outcomes of an infection (extinction of cells, survival
with complete elimination of viruses, or persistent infection) crucially depend
on mutation rates of viral and Immunoglobulin proteins. The model predicts that
there is an optimal Somatic Hypermutation (SHM) rate close to experimentally
observed 10-3 per nucleotide per replication. Further, we developed an
analytical theory which explains the physical reason for an optimal SHM program
as a compromise between deleterious effects of random mutations on nascent
maturing Immunoglobulins (adversity) and the need to generate diverse pool of
mutated antibodies from which highly potent ones can be drawn (diversity). The
theory explains such effects as dependence of B cell fate on affinity for an
incoming antigen, ceiling in affinity of mature antibodies, Germinal Center
sizes and maturation times. The theory reveals the molecular factors which
determine the efficiency of affinity maturation, providing insight into
variability of immune response to cytopathic (direct response by germline
antibodies) and poorly cytopathic viruses (crucial role of SHM in response).
These results demonstrate the feasibility and promise of microscopic
sequence-based models of immune system, where population dynamics of evolving
Immunoglobulins is explicitly tied to their molecular properties.


  Rapidly expanding insect populations, deforestation, and global climate
change threaten to destabilize key planetary carbon pools, especially the
Earth's forests which link the micro-ecology of insect infestation to climate.
To the extent mean temperature increases, insect populations accelerate
deforestation. This alters climate via the loss of active carbon sequestration
by live trees and increased carbon release from decomposing dead trees. A
positive feedback loop can emerge that is self-sustaining--no longer requiring
independent climate-change drivers. Current research regimes and insect control
strategies are insufficient at present to cope with the present regional scale
of insect-caused deforestation, let alone its likely future global scale.
Extensive field recordings demonstrate that bioacoustic communication plays a
role in infestation dynamics and is likely to be a critical link in the
feedback loop. These results open the way to novel detection and monitoring
strategies and nontoxic control interventions.


  Closeness is a global measure of centrality in networks, and a proxy for how
influential actors are in social networks. In most network models, and many
empirical networks, closeness is strongly correlated with degree. However, in
social networks there is a cost of maintaining social ties. This leads to a
situation (that can occur in the professional social networks of executives,
lobbyists, diplomats and so on) where agents have the conflicting objectives of
aiming for centrality while simultaneously keeping the degree low. We
investigate this situation in an adaptive network-evolution model where agents
optimize their positions in the network following individual strategies, and
using only local information. The strategies are also optimized, based on the
success of the agent and its neighbors. We measure and describe the time
evolution of the network and the agents' strategies.


  In eukaryotic genomes, nucleosomes function to compact DNA and to regulate
access to it both by simple physical occlusion and by providing the substrate
for numerous covalent epigenetic tags. While nucleosome positions in vitro are
determined by sequence alone, in vivo competition with other DNA-binding
factors and action of chromatin remodeling enzymes play a role that needs to be
quantified. We developed a biophysical model for the sequence dependence of DNA
bending energies, and validated it against a collection of in vitro free
energies of nucleosome formation and a nucleosome crystal structure; we also
successfully designed both strong and poor histone binding sequences ab initio.
For in vivo data from S.cerevisiae, the strongest positioning signal came from
the competition with other factors. Based on sequence alone, our model predicts
that functional transcription factor binding sites have a tendency to be
covered by nucleosomes, but are uncovered in vivo because functional sites
cluster within a single nucleosome footprint, making transcription factors bind
cooperatively. Similarly a weak enhancement of nucleosome binding in the TATA
region for naked DNA becomes a strong depletion when the TATA-binding protein
is included, in quantitative agreement with experiment. Predictions at specific
loci were also greatly enhanced by including competing factors. Our physically
grounded model distinguishes multiple ways in which genomic sequence can
influence nucleosome positions and thus provides an alternative explanation for
several important experimental findings.


  The problem of reliability of the dynamics in biological regulatory networks
is studied in the framework of a generalized Boolean network model with
continuous timing and noise. Using well-known artificial genetic networks such
as the repressilator, we discuss concepts of reliability of rhythmic
attractors. In a simple evolution process we investigate how overall network
structure affects the reliability of the dynamics. In the course of the
evolution, networks are selected for reliable dynamics. We find that most
networks can be easily evolved towards reliable functioning while preserving
the original function.


  The length of coding sequence series in microbial genomes were regarded as a
fluctuating system and characterized by the methods of statistical physics. The
distribution and the correlatin properties of 50 genomes including bacteria and
several archaea were investigated. The distribution was investigated by
rank-size analysis (Zipf's law. We found that coding sequence lengths series do
not obey Zipf's law contrary to natural languages. The distribution was found
to be more closely to an exponential distribution. The correlation appeared to
be similar to natural languages. Segmentation analysis of the series showed to
be short range memory systems.


  I compare two quantum-theoretical approaches to the phenomenon of adaptive
mutations, termed here Q-cell and Q-genome. I use 'fluctuation trapping' model
as a general framework. I introduce notions of R-error and D-error and argue
that the 'fluctuation trapping' model has to employ a correlation between the
R- and D- errors. Further, I compare how the two approaches can justify the
R-D-error correlation, focusing on the advantages of the Q-cell approach. The
positive role of environmentally induced decoherence (EID) on both steps of the
adaptation process is emphasized. A starving bacterial cell is proposed to be
in an einselected state. The intracellular dynamics in this state has a unitary
character and I propose to interpret it as 'exponential growth in imaginary
time', analogously to the commonly considered 'diffusion' interpretation of the
Schroedinger equation. Addition of a substrate leads to Wick rotation and a
switch from 'imaginary time' reproduction to a 'real time' reproduction regime.
Due to the variations at the genomic level (such as base tautomery), the
starving cell has to be represented as a superposition of different components,
all 'reproducing in imaginary time'. Adidtion of a selective substrate,
allowing only one of these components to amplify, will cause Wick rotation and
amplification of this component, thus justifying the occurence of the R-D-error
correlation. Further ramifications of the proposed ideas for evolutionary
theory are discussed.


  We prove the existence of infinitely many periodic solutions and complicated
dynamics, due to the presence of a topological horseshoe, for the classical
Volterra predator--prey model with a periodic harvesting. The proof relies on
some recent results about chaotic planar maps combined with the study of
geometric features which are typical of linked twist maps.


  In this work we focus on how noise propagates in biochemical reaction
networks and affects sensitivities of the system. We discover that the
stochastic fluctuations can enhance sensitivities in one region of the value of
control parameters by reducing sensitivities in another region. Based on this
compensation principle, we designed a concentration detector in which enhanced
amplification is achieved by an incoherent feedforward reaction network.


  We analyze the dynamics of rotary biomotors within a simple
nano-electromechanical model, consisting of a stator part and a ring-shaped
rotor having twelve proton-binding sites. This model is closely related to the
membrane-embedded F$_0$ motor of adenosine triphosphate (ATP) synthase, which
converts the energy of the transmembrane electrochemical gradient of protons
into mechanical motion of the rotor. It is shown that the Coulomb coupling
between the negative charge of the empty rotor site and the positive stator
charge, located near the periplasmic proton-conducting channel (proton source),
plays a dominant role in the torque-generating process. When approaching the
source outlet, the rotor site has a proton energy level higher than the energy
level of the site, located near the cytoplasmic channel (proton drain). In the
first stage of this torque-generating process, the energy of the
electrochemical potential is converted into potential energy of the
proton-binding sites on the rotor. Afterwards, the tangential component of the
Coulomb force produces a mechanical torque. We demonstrate that, at low
temperatures, the loaded motor works in the shuttling regime where the energy
of the electrochemical potential is consumed without producing any
unidirectional rotation. The motor switches to the torque-generating regime at
high temperatures, when the Brownian ratchet mechanism turns on. In the
presence of a significant external torque, created by ATP hydrolysis, the
system operates as a proton pump, which translocates protons against the
transmembrane potential gradient. Here we focus on the F$_0$ motor, even though
our analysis is applicable to the bacterial flagellar motor.


  The set of known dialects of the genetic code (GC) is analyzed from the
viewpoint of the genetic octave Yin-Yang-algebra. This algebra was described in
the previous author's publications. The algebra was discovered on the basis of
structural features of the GC in the matrix form of its presentation ("matrix
genetics"). The octave Yin-Yang-algebra is considered as the pre-code or as the
model of the GC. From the viewpoint of this algebraic model, for example, the
sets of 20 amino acids and of 64 triplets consist of sub-sets of "male",
"female" and "androgynous" molecules, etc. This algebra permits to reveal
hidden peculiarities of the structure and evolution of the GC and to propose
the conception of "sexual" relationships among genetic molecules. The first
results of the analysis of the GC systems from such algebraic viewpoint say
about the close connection between evolution of the GC and this algebra. They
include 8 evolutionary rules of the dialects of the GC. The evolution of the GC
is appeared as the struggle between male and female beginnings. The hypothesis
about new biophysical factor of "sexual" interactions among genetic molecules
is put forward. The matrix forms of presentation of elements of the genetic
octave Yin-Yang-algebra are connected with Hadamard matrices by means of the
simple U-algorithm. Hadamard matrices play a significant role in the theory of
quantum computers, in particular. It gives new opportunities for possible
understanding the GC systems as quantum computer systems. Revealed algebraic
properties of the GC permit to put forward the problem of algebraization of
bioinformatics on the basis of the algebras of the GC. Our investigation is
connected with the question: what is life from the viewpoint of algebra? The
algebraic version of the origin of the GC is discussed.


  The Cambrian explosion is a grand challenge to science today and involves
multidisciplinary study. This event is generally believed as a result of
genetic innovations, environmental factors and ecological interactions, even
though there are many conflicts on nature and timing of metazoan origins. The
crux of the matter is that an entire roadmap of the evolution is missing to
discern the biological complexity transition and to evaluate the critical role
of the Cambrian explosion in the overall evolutionary context. Here we
calculate the time of the Cambrian explosion by an innovative and accurate
"C-value clock"; our result (560 million years ago) quite fits the fossil
records. We clarify that the intrinsic reason of genome evolution determined
the Cambrian explosion. A general formula for evaluating genome size of
different species has been found, by which major questions of the C-value
enigma can be solved and the genome size evolution can be illustrated. The
Cambrian explosion is essentially a major transition of biological complexity,
which corresponds to a turning point in genome size evolution. The observed
maximum prokaryotic complexity is just a relic of the Cambrian explosion and it
is supervised by the maximum information storage capability in the observed
universe. Our results open a new prospect of studying metazoan origins and
molecular evolution.


  This is the fourth paper, the last one, on solution to the problem of absence
of detailed balance in nonequilibrium processes. It is an approach based on
another known universal dynamics: The evolutionary dynamics first conceived by
Darwin and Wallace, referring to as Darwinian dynamics in the present paper,
has been found to be universally valid in biology; The statistical mechanics
and thermodynamics, while enormously successful in physics, have been in an
awkward situation of wanting a consistent dynamical understanding; Here we
present from a formal point of view an exploration of the connection between
thermodynamics and Darwinian dynamics and a few related topics. We first show
that the stochasticity in Darwinian dynamics implies the existence temperature,
hence the canonical distribution of Boltzmann-Gibbs type. In term of relative
entropy the Second Law of thermodynamics is dynamically demonstrated without
detailed balance condition, and is valid regardless of size of the system. In
particular, the dynamical component responsible for breaking detailed balance
condition does not contribute to the change of the relative entropy. Two types
of stochastic dynamical equalities of current interest are explicitly discussed
in the present approach: One is based on Feynman-Kac formula and another is a
generalization of Einstein relation. Both are directly accessible to
experimental tests. Our demonstration indicates that Darwinian dynamics
represents logically a simple and straightforward starting point for
statistical mechanics and thermodynamics and is complementary to and consistent
with conservative dynamics that dominates the physical sciences. Present
exploration suggests the existence of a unified stochastic dynamical framework
both near and far from equilibrium.


  Much evolutionary information is stored in the fluctuations of protein length
distributions. The genome size and non-coding DNA content can be calculated
based only on the protein length distributions. So there is intrinsic
relationship between the coding DNA size and non-coding DNA size. According to
the correlations and quasi-periodicity of protein length distributions, we can
classify life into three domains. Strong evidences are found to support the
order in the structures of protein length distributions.


  We show that the Adomian decomposition method, the time--series expansion,
the homotopy--perturbation method, and the variational--iteration method
completely fail to provide a reasonable description of the dynamics of the
simplest prey--predator system.


  Double-stranded DNA `overstretches' at a pulling force of about 65 pN,
increasing in length by a factor of 1.7. The nature of the overstretched state
is unknown, despite its considerable importance for DNA's biological function
and technological application. Overstretching is thought by some to be a
force-induced denaturation, and by others to consist of a transition to an
elongated, hybridized state called S-DNA. Within a statistical mechanical model
we consider the effect upon overstretching of extreme sequence heterogeneity.
`Chimeric' sequences possessing halves of markedly different AT composition
elongate under fixed external conditions via distinct, spatially segregated
transitions. The corresponding force-extension data display two plateaux at
forces whose difference varies with pulling rate in a manner that depends
qualitatively upon whether the hybridized S-form is accessible. This
observation implies a test for S-DNA that could be performed in experiment. Our
results suggest that qualitatively different, spatially segregated
conformational transitions can occur at a single thermodynamic state within
single molecules of DNA.


  The conformation of circular DNA molecules of various lengths adsorbed in a
2D conformation on a mica surface is studied. The results confirm the
conjecture that the critical exponent $\nu$ is topologically invariant and
equal to the SAW value (in the present case $\nu=3/4$), and that the topology
and dimensionality of the system strongly influences the cross-over between the
rigid regime and the self-avoiding regime at a scale $L\approx 8 \ell_p$.
Additionally, the bond correlation function scales with the molecular length
$L$ as predicted. For molecular lengths $L\leq5 \ell_p$, circular DNA behaves
like a stiff molecule with approximately elliptic shape.


  Prediction and control of cancer invasion is a vital problem in medical
science. This paper proposes a modern geometric Ricci-flow and entropy based
model for control of avascular multicellular tumor spheroid growth and decay.
As a tumor growth/decay control tool, a monoclonal antibody therapy is
proposed.
  Keywords: avascular tumor growth and decay, multicellular tumor spheroid,
Ricci flow and entropy, nonlinear heat equation, monoclonal antibody cancer
therapy


  Chemical reactions involving radical-ion pairs are ubiquitous in biology,
since not only are they at the basis of the photosynthetic reaction chain, but
are also assumed to underlie the biochemical magnetic compass used by avian
species for navigation. Recent experiments with magnetic-sensitive radical-ion
pair reactions provided strong evidence for the radical-ion-pair
magnetoreception mechanism, verifying the expected magnetic sensitivities and
chemical product yield changes. It is here shown that the theoretical
description of radical-ion-pair reactions used since the 70's cannot explain
the observed data, because it is based on phenomenological equations masking
quantum coherence effects. The fundamental density matrix equation derived here
from basic quantum measurement theory considerations naturally incorporates the
quantum Zeno effect and readily explains recent experimental observations on
low- and high-magnetic-field radical-ion-pair reactions.


  We present the first numerical simulation of actin-driven propulsion by
elastic filaments. Specifically, we use a Brownian dynamics formulation of the
dendritic nucleation model of actin-driven propulsion. We show that the model
leads to a self-assembled network that exerts forces on a disk and pushes it
with an average speed. This simulation approach is the first to observe a speed
that varies non-monotonically with the concentration of branching proteins
(Arp2/3), capping protein and depolymerization rate (ADF), in accord with
experimental observations. Our results suggest a new interpretation of the
origin of motility that can be tested readily by experiment.


  We propose a model for evolution aiming to reproduce statistical features of
fossil data, in particular the distributions of extinction events, the
distribution of species per genus and the distribution of lifetimes, all of
which are known to be of power law type. The model incorporates both
species-species interactions and ancestral relationships. The main novelty of
this work is to show the feasibility of k-core percolation as a selection
mechanism. We give theoretical predictions for the observable distributions,
confirm their validity by computer simulation and report good agreement with
the fossil data. A key feature of the proposed model is a co-evolving fitness
landscape determined by the topology of the underlying species interactions,
ecological niches emerge naturally. The predicted distributions are independent
of the rate of speciation, i.e. whether one adopts an gradualist or punctuated
view of evolution.


  We explore the use of a top-down approach to analyse the dynamics of
icosahedral virus capsids and complement the information obtained from
bottom-up studies of viral vibrations available in the literature. A normal
mode analysis based on protein association energies is used to study the
frequency spectrum, in which we reveal a universal plateau of low-frequency
modes shared by a large class of Caspar-Klug capsids. These modes break
icosahedral symmetry and are potentially relevant to the genome release
mechanism. We comment on the role of viral tiling theory in such dynamical
considerations.


  We study the dynamics of a deterministic walk confined in a narrow
two-dimensional space randomly filled with point-like targets. At each step,
the walker visits the nearest target not previously visited. Complex dynamics
is observed at some intermediate values of the domain width, when, while
drifting, the walk performs long intermittent backward excursions. As the width
is increased, evidence of a transition from ballistic motion to a weakly
non-ergodic regime is shown, characterized by sudden inversions of the drift
velocity with a probability slowly decaying with time, as $1/t$ at leading
order. Excursion durations, first-passage times and the dynamics of unvisited
targets follow power-law distributions. For parameter values below this scaling
regime, precursory patterns in the form of "wild" outliers are observed, in
close relation with the presence of log-oscillations in the probability
distributions. We discuss the connections between this model and several
evolving biological systems.


  Ecdysteroid hormones are major regulators in reproduction and development of
insects, including larval molts and metamorphosis. The functional ecdysone
receptor is a heterodimer of ECR (NR1H1) and USP-RXR (NR2B4), which is the
orthologue of vertebrate retinoid X receptors (RXR alpha, beta, gamma). Both
proteins belong to the superfamily of nuclear hormone receptors,
ligand-dependent transcription factors that share two conserved domains: the
DNA-binding domain (DBD) and the ligand-binding domain (LBD). In order to gain
further insight into the evolution of metamorphosis and gene regulation by
ecdysone in arthropods, we performed a phylogenetic analysis of both partners
of the heterodimer ECR/USP-RXR. Overall, 38 USP-RXR and 19 ECR protein
sequences, from 33 species, have been used for this analysis. Interestingly,
sequence alignments and structural comparisons reveal high divergence rates,
for both ECR and USP-RXR, specifically among Diptera and Lepidoptera. The most
impressive differences affect the ligand-binding domain of USP-RXR. In
addition, ECR sequences show variability in other domains, namely the
DNA-binding and the carboxy-terminal F domains. Our data provide the first
evidence that ECR and USP-RXR may have coevolved during holometabolous insect
diversification, leading to a functional divergence of the ecdysone receptor.
These results have general implications on fundamental aspects of insect
development, evolution of nuclear receptors, and the design of specific
insecticides.


  We show that spatial models of simple predator-prey interactions predict that
predator and prey numbers oscillate in time and space. These oscillations are
not seen in the deterministic versions of the models, but are due to stochastic
fluctuations about the time-independent solutions of the deterministic
equations which are amplified due to the existence of a resonance. We calculate
the power spectra of the fluctuations analytically and show that they agree
well with results obtained from stochastic simulations. This work extends the
analysis of these quasi-cycles from that previously developed for well-mixed
systems to spatial systems, and shows that the ideas and methods used for
non-spatial models naturally generalize to the spatial case.


  Isochores are long genome segments relatively homogeneous in G+C. A heuristic
algorithm based on entropic segmentation has been developed by our group, and a
web server implementing all the required components is available. However, a
researcher may want to perform batch processing of many sequences
simultaneously in its local machine, instead of analyzing them on one by one
basis through the web. To this end, standalone versions are required. We report
here the implementation of two standalone programs, able to predict isochores
at the sequence level: 1) a command-line version (IsoFinder) for Windows and
Linux systems; and 2) a user-friendly version (IsoFinderWin) running under
Windows.


  The Hadamard transform of \cite{Hendy89a, Hendy89} provides a way to work
with stochastic models for sequence evolution without having to deal with the
complications of tree space and the graphical structure of trees. Here we
demonstrate that the transform can be expressed in terms of the familiar
$\bP[\bt] = e^{\bQ[\bt]}$ formula for Markov chains. The key idea is to study
the evolution of vectors of states, one vector entry for each taxa; we call
this the $n$-taxon process. We derive transition probabilities for the process.
Significantly, the findings show that tree-based models are indeed in the
family of (multi-variate) exponential distributions.


  We calculate the probability of DNA loop formation mediated by regulatory
proteins such as Lac repressor (LacI), using a mathematical model of DNA
elasticity. Our model is adapted to calculating quantities directly observable
in Tethered Particle Motion (TPM) experiments, and it accounts for all the
entropic forces present in such experiments. Our model has no free parameters;
it characterizes DNA elasticity using information obtained in other kinds of
experiments. [...] We show how to compute both the "looping J factor" (or
equivalently, the looping free energy) for various DNA construct geometries and
LacI concentrations, as well as the detailed probability density function of
bead excursions. We also show how to extract the same quantities from recent
experimental data on tethered particle motion, and then compare to our model's
predictions. [...] Our model successfully reproduces the detailed distributions
of bead excursion, including their surprising three-peak structure, without any
fit parameters and without invoking any alternative conformation of the LacI
tetramer. Indeed, the model qualitatively reproduces the observed dependence of
these distributions on tether length (e.g., phasing) and on LacI concentration
(titration). However, for short DNA loops (around 95 basepairs) the experiments
show more looping than is predicted by the harmonic-elasticity model, echoing
other recent experimental results. Because the experiments we study are done in
vitro, this anomalously high looping cannot be rationalized as resulting from
the presence of DNA-bending proteins or other cellular machinery. We also show
that it is unlikely to be the result of a hypothetical "open" conformation of
the LacI tetramer.


  Spatial structure is known to have an impact on the evolution of cooperation,
and so it has been intensively studied during recent years. Previous work has
shown the relevance of some features, such as the synchronicity of the
updating, the clustering of the network or the influence of the update rule.
This has been done, however, for concrete settings with particular games,
networks and update rules, with the consequence that some contradictions have
arisen and a general understanding of these topics is missing in the broader
context of the space of 2x2 games. To address this issue, we have performed a
systematic and exhaustive simulation in the different degrees of freedom of the
problem. In some cases, we generalize previous knowledge to the broader context
of our study and explain the apparent contradictions. In other cases, however,
our conclusions refute what seems to be established opinions in the field, as
for example the robustness of the effect of spatial structure against changes
in the update rule, or offer new insights into the subject, e.g. the relation
between the intensity of selection and the asymmetry between the effects on
games with mixed equilibria.


  We use a biophysical model of a local neuronal circuit to study the
implications of synaptic plasticity for the detection of weak sensory stimuli.
Networks with fast plastic coupling show behavior consistent with stochastic
resonance. Addition of an additional slow coupling that accounts for the
asynchronous release of neurotransmitter results in qualitatively different
properties of signal detection, and also leads to the appearance of transient
post-stimulus bistability. Our results suggest testable hypothesis with regard
to the self-organization and dynamics of local neuronal circuits.


  A dynamic logic method was developed to analyze molecular networks of cells
by combining Kauffman and Thomas's logic operations with molecular interaction
parameters. The logic operations characterize the discrete interactions between
biological components. The interaction parameters (e.g. response times)
describe the quantitative kinetics. The combination of the two quantitatively
characterizes the discrete biological interactions. A number of simple networks
were analyzed. The main results include: we proved the theorems to determine
bistable states and oscillation behaviors of networks, we showed that time
delays are essential for oscillation structures, we proved that single variable
networks do not have chaotic behaviors, and we explained why one signal can
have multiply responses. In addition, we applied the present method to the
analysis of the MAPK cascade, feed-forward loops, and mitosis cycle of budding
yeast cells.


  In many cases, transcriptional regulation involves the binding of
transcription factors at sites on the DNA that are not immediately adjacent to
the promoter of interest. This action at a distance is often mediated by the
formation of DNA loops: Binding at two or more sites on the DNA results in the
formation of a loop, which can bring the transcription factor into the
immediate neighborhood of the relevant promoter. Though there have been a
variety of insights into the combinatorial aspects of transcriptional control,
the mechanism of DNA looping as an agent of combinatorial control in both
prokaryotes and eukaryotes remains unclear. We use single-molecule techniques
to dissect DNA looping in the lac operon. In particular, we measure the
propensity for DNA looping by the Lac repressor as a function of the
concentration of repressor protein and as a function of the distance between
repressor binding sites. As with earlier single-molecule studies, we find (at
least) two distinct looped states and demonstrate that the presence of these
two states depends both upon the concentration of repressor protein and the
distance between the two repressor binding sites. We find that loops form even
at interoperator spacings considerably shorter than the DNA persistence length,
without the intervention of any other proteins to prebend the DNA. The
concentration measurements also permit us to use a simple statistical
mechanical model of DNA loop formation to determine the free energy of DNA
looping, or equivalently, the J-factor for looping.


  We study the effects of switching social contacts as a strategy to control
epidemic outbreaks. Connections between susceptible and infective individuals
can be broken by either individual, and then reconnected to a randomly chosen
member of the population. It is assumed that the reconnecting individual has no
previous information on the epidemiological condition of the new contact. We
show that reconnection can completely suppress the disease, both by continuous
and discontinuous transitions between the endemic and the infection-free
states. For diseases with an asymptomatic phase, we analyze the conditions for
the suppression of the disease, and show that, even when these conditions are
not met, the increase of the endemic infection level is usually rather small.
We conclude that, within some simple epidemiological models, contact switching
is a quite robust and effective control strategy. This suggests that it may
also be an efficient method in more complex situations.


  The Raman spectrum of microscopic droplets of lipid material on the lumen of
the aorta of an apolipoprotein E knock-out mouse were reported in the
proceedings of an SPIE conference. Based on the absence of the carbonyl band in
the spectrum, at that time it was determined that the spectrum represents free
fatty acids rather than triglycerides. More recent examination of the spectrum
indicates that these droplets contain dissolved cholesterol, and thus can be
used as an early indicator of the atherosclerosis process in animal models
during drug development.


  Dissimilarity measures for (possibly weighted) phylogenetic trees based on
the comparison of their vectors of path lengths between pairs of taxa, have
been present in the systematics literature since the early seventies. But, as
far as rooted phylogenetic trees goes, these vectors can only separate
non-weighted binary trees, and therefore these dissimilarity measures are
metrics only on this class. In this paper we overcome this problem, by
splitting in a suitable way each path length between two taxa into two lengths.
We prove that the resulting splitted path lengths matrices single out arbitrary
rooted phylogenetic trees with nested taxa and arcs weighted in the set of
positive real numbers. This allows the definition of metrics on this general
class by comparing these matrices by means of metrics in spaces of real-valued
$n\times n$ matrices. We conclude this paper by establishing some basic facts
about the metrics for non-weighted phylogenetic trees defined in this way using
$L^p$ metrics on these spaces of matrices.


  The accurate prediction and characterization of DNA melting domains by
computational tools could facilitate a broad range of biological applications.
However, no algorithm for melting domain prediction has been available until
now. The main challenges include the difficulty of mathematically mapping a
qualitative description of DNA melting domains to quantitative statistical
mechanics models, as well as the absence of 'gold standards' and a need for
generality. In this paper, we introduce a new approach to identify the twostate
regions and melting fork regions along a given DNA sequence. Compared with an
ad hoc segmentation used in one of our previous studies, the new algorithm is
based on boundary probability profiles, rather than standard melting maps. We
demonstrate that a more detailed characterization of the DNA melting domain map
can be obtained using our new method, and this approach is independent of the
choice of DNA melting model. We expect this work to drive our understanding of
DNA melting domains one step further.


  Hypercycles are information integration systems which are thought to overcome
the information crisis of prebiotic evolution by ensuring the coexistence of
several short templates. For imperfect template replication, we derive a simple
expression for the maximum number of distinct templates $n_m$ that can coexist
in a hypercycle and show that it is a decreasing function of the length $L$ of
the templates. In the case of high replication accuracy we find that the
product $n_m L$ tends to a constant value, limiting thus the information
content of the hypercycle. Template coexistence is achieved either as a
stationary equilibrium (stable fixed point) or a stable periodic orbit in which
the total concentration of functional templates is nonzero. For the hypercycle
system studied here we find numerical evidence that the existence of an
unstable fixed point is a necessary condition for the presence of periodic
orbits.


  The evolution of genetic systems has been analyzed through the use of
modifier gene models, in which a neutral gene is posited to control the
transmission of other genes under selection. Analysis of modifier gene models
has found the manifestations of an "evolutionary reduction principle": in a
population near equilibrium, a new modifier allele that scales equally all
transition probabilities between different genotypes under selection can invade
if and only if it reduces the transition probabilities. Analytical results on
the reduction principle have always required some set of constraints for
tractability: limitations to one or two selected loci, two alleles per locus,
specific selection regimes or weak selection, specific genetic processes being
modified, extreme or infinitesimal effects of the modifier allele, or tight
linkage between modifier and selected loci. Here, I prove the reduction
principle in the absence of any of these constraints, confirming a twenty-year
old conjecture. The proof is obtained by a wider application of Karlin's
Theorem 5.2 (1982) and its extension to ML-matrices, substochastic matrices,
and reducible matrices.


  We study the spreading of excitations in 2D systems of mobile agents where
the excitation is transmitted when a quiescent agent keeps contact with an
excited one during a non-vanishing time. We show that the steady states
strongly depend on the spatial agent dynamics. Moreover, the coupling between
exposition time ($\omega$) and agent-agent contact rate (CR) becomes crucial to
understand the excitation dynamics, which exhibits three regimes with CR: no
excitation for low CR, an excited regime in which the number of quiescent
agents (S) is inversely proportional to CR, and for high CR, a novel third
regime, model dependent, here S scales with an exponent $\xi -1$, with $\xi$
being the scaling exponent of $\omega$ with CR.


  Elucidating the fitness measures optimized during the evolution of complex
biological systems is a major challenge in evolutionary theory. We present
experimental evidence and an analytical framework demonstrating how biochemical
networks exploit optimal control strategies in their evolutionary dynamics.
Optimal control theory explains a striking pattern of extremization in the
redox potentials of electron transport proteins, assuming only that their
fitness measure is a control objective functional with bounded controls.


  The complex dynamics of intracellular calcium regulates cellular responses to
information encoded in extracellular signals. Here, we study the encoding of
these external signals in the context of the Li-Rinzel model. We show that by
control of biophysical parameters the information can be encoded in amplitude
modulation, frequency modulation or mixed (AM and FM) modulation. We briefly
discuss the possible implications of this new role of information encoding for
astrocytes.


  The sequence of amino acids in a protein is believed to determine its native
state structure, which in turn is related to the functionality of the protein.
In addition, information pertaining to evolutionary relationships is contained
in homologous sequences. One powerful method for inferring these sequence
attributes is through comparison of a query sequence with reference sequences
that contain significant homology and whose structure, function, and/or
evolutionary relationships are already known. In spite of decades of concerted
work, there is no simple framework for deducing structure, function, and
evolutionary (SF&E) relationships directly from sequence information alone,
especially when the pair-wise identity is less than a threshold figure ~25%
[1,2]. However, recent research has shown that sequence identity as low as 8%
is sufficient to yield common structure/function relationships and sequence
identities as large as 88% may yet result in distinct structure and function
[3,4]. Starting with a basic premise that protein sequence encodes information
about SF&E, one might ask how one could tease out these measures in an unbiased
manner. Here we present a unified framework for inferring SF&E from sequence
information using a knowledge-based approach which generates phylogenetic
profiles in an unbiased manner. We illustrate the power of phylogenetic
profiles generated using the Gestalt Domain Detection Algorithm Basic Local
Alignment Tool (GDDA-BLAST) to derive structural domains, functional
annotation, and evolutionary relationships for a host of ion-channels and human
proteins of unknown function. These data are in excellent accord with published
data and new experiments. Our results suggest that there is a wealth of
previously unexplored information in protein sequence.


  In evolutionary biology, genetic sequences carry with them a trace of the
underlying tree that describes their evolution from a common ancestral
sequence. The question of how many sequence sites are required to recover this
evolutionary relationship accurately depends on the model of sequence
evolution, the substitution rate, divergence times and the method used to infer
phylogenetic history. A particularly challenging problem for phylogenetic
methods arises when a rapid divergence event occurred in the distant past. We
analyse an idealised form of this problem in which the terminal edges of a
symmetric four--taxon tree are some factor ($p$) times the length of the
interior edge. We determine an order $p^2$ lower bound on the growth rate for
the sequence length required to resolve the tree (independent of any particular
branch length). We also show that this rate of sequence length growth can be
achieved by existing methods (including the simple `maximum parsimony' method),
and compare these order $p^2$ bounds with an order $p$ growth rate for a model
that describes low-homoplasy evolution. In the final section, we provide a
generic bound on the sequence length requirement for a more general class of
Markov processes.


  We discuss the problem of proteasomal degradation of proteins. Though
proteasomes are important for all aspects of the cellular metabolism, some
details of the physical mechanism of the process remain unknown. We introduce a
stochastic model of the proteasomal degradation of proteins, which accounts for
the protein translocation and the topology of the positioning of cleavage
centers of a proteasome from first principles. For this model we develop the
mathematical description based on a master-equation and techniques for
reconstruction of the cleavage specificity inherent to proteins and the
proteasomal translocation rates, which are a property of the proteasome specie,
from mass spectroscopy data on digestion patterns. With these properties
determined, one can quantitatively predict digestion patterns for new
experimental set-ups. Additionally we design an experimental set-up for a
synthetic polypeptide with a periodic sequence of amino acids, which enables
especially reliable determination of translocation rates.


  Evolutionary dynamics can be studied in well-mixed or structured populations.
Population structure typically arises from the heterogeneous distribution of
individuals in physical space or on social networks. Here we introduce a new
type of space to evolutionary game dynamics: phenotype space. The population is
well-mixed in the sense that everyone is equally likely to interact with
everyone else, but the behavioral strategies depend on distance in phenotype
space. Individuals might behave differently towards those who look similar or
dissimilar. Individuals mutate to nearby phenotypes. We study the `phenotypic
space walk' of populations. We present analytic calculations that bring
together ideas from coalescence theory and evolutionary game dynamics. As a
particular example, we investigate the evolution of cooperation in phenotype
space. We obtain a precise condition for natural selection to favor cooperators
over defectors: for a one-dimensional phenotype space and large population size
the critical benefit-to-cost ratio is given by b/c=1+2/sqrt{3}. We derive the
fundamental condition for any evolutionary game and explore higher dimensional
phenotype spaces.


  One approach to studying the system-wide organization of biochemistry is to
use statistical graph theory. Even in such a heavily simplified method, which
disregards most of the dynamic aspects of biochemistry, one is faced with
fundamental questions, such as how the chemical reaction systems should be
reduced to a graph retaining as much functional information as possible from
the original reaction system. In such graph representations, should the edges
go between substrates and products, or substrates and substrates, or both?
Should vertices represent substances or reactions? Different definitions encode
different information about the reaction system. In this paper we evaluate four
different graph representations of metabolism, applied to data from different
organisms and databases. The graph representations are evaluated by comparing
the overlap between clusters (network modules) and annotated functions, and
also by comparing the set of identified currency metabolites with those that
other authors have identified using qualitative biological arguments. We find
that a "substance network," where all metabolites participating in a reaction
are connected, is relatively better than others, evaluated both with respect to
the functional overlap between modules and functions and to the number and
identity of identified currency metabolites.


  Reverse engineering of gene regulatory networks presents one of the big
challenges in systems biology. Gene regulatory networks are usually inferred
from a set of single-gene over-expressions and/or knockout experiments.
Functional relationships between genes are retrieved either from the steady
state gene expressions or from respective time series. We present a novel
algorithm for gene network reconstruction on the basis of steady-state
gene-chip data from over-expression experiments. The algorithm is based on a
straight forward solution of a linear gene-dynamics equation, where
experimental data is fed in as a first predictor for the solution. We compare
the algorithm's performance with the NIR algorithm, both on the well known
E.Coli experimental data and on in-silico experiments. We show superiority of
the proposed algorithm in the number of correctly reconstructed links and
discuss computational time and robustness. The proposed algorithm is not
limited by combinatorial explosion problems and can be used in principle for
large networks of thousands of genes.


  We apply a simulational proxy of the phi-value analysis and perform extensive
mutagenesis experiments to identify the nucleating residues in the folding
reactions of two small lattice Go polymers with different native geometries.
These results are compared with those obtained from an accurate analysis based
on the reaction coordinate folding probability Pfold, and on structural
clustering methods. For both protein models, the transition state ensemble is
rather heterogeneous and splits-up into structurally different populations. For
the more complex geometry the identified subpopulations are actually
structurally disjoint. For the less complex native geometry we found a broad
transition state with microscopic heterogeneity. For both geometries, the
identification of the folding nucleus via the Pfold analysis agrees with the
identification of the folding nucleus carried out with the phi-value analysis.
For the most complex geometry, however, the apllied methodologies give more
consistent results than for the more local geometry. The study of the
transition state' structure reveals that the nucleus residues are not
necessarily fully native in the transition state. Indeed, it is only for the
more complex geometry that two of the five critical residues show a
considerably high probability of having all its native bonds formed in the
transition state. Therefore, one concludes that in general the phi-value
correlates with the acceleration/deceleration of folding induced by mutation,
rather than with the degree of nativeness of the transition state, and that the
traditional interpretation of phi-values may provide a more realistic picture
of the structure of the transition state only for more complex native
geometries.


  The classical electrostatic interaction between DNA molecules in water in the
presence of counterions is reconsidered and we propose it is governed by a
modified Poisson-Boltzmann equation. Quantum fluctuations are then studied and
shown to lead to a vacuum interaction that is numerically computed for several
configurations of many DNA strands and found to be strongly many-body. This
Casimir vacuum interaction can be the ``glue'' holding together DNA molecules
into aggregates.


  Shuttle-assisted charge transfer is pivotal for the efficient energy
transduction from the food-stuff electrons to protons in the respiratory chain
of animal cells and bacteria. The respiratory chain consists of four
metalloprotein Complexes (I-IV) embedded in the inner membrane of a
mitochondrion. Three of these complexes pump protons across the membrane,
fuelled by the energy of food-stuff electrons. Despite extensive biochemical
and biophysical studies, the physical mechanism of this proton pumping is still
not well understood. Here we present a nanoelectromechanical model of the
electron-driven proton pump related to the second loop of the respiratory
chain, where a lipid-soluble ubiquinone molecule shuttles between the Complex I
and Complex III, carrying two electrons and two protons. We show that the
energy of electrons can be converted to the transmembrane proton potential
gradient via the electrostatic interaction between electrons and protons on the
shuttle. We find that the system can operate either as a proton pump, or, in
the reverse regime, as an electron pump. For membranes with various
viscosities, we demonstrate that the uphill proton current peaks near the body
temperature $T \approx 37 ^{\circ}$C.


  Most of the statistical tests currently used to detect differentially
expressed genes are based on asymptotic results, and perform poorly for low
expression tags. Another problem is the common use of a single canonical cutoff
for the significance level (p-value) of all the tags, without taking into
consideration the type II error and the highly variable character of the sample
size of the tags.
  This work reports the development of two significance tests for the
comparison of digital expression profiles, based on frequentist and Bayesian
points of view, respectively. Both tests are exact, and do not use any
asymptotic considerations, thus producing more correct results for low
frequency tags than the chi-square test. The frequentist test uses a
tag-customized critical level which minimizes a linear combination of type I
and type II errors. A comparison of the Bayesian and the frequentist tests
revealed that they are linked by a Beta distribution function. These tests can
be used alone or in conjunction, and represent an improvement over the
currently available methods for comparing digital profiles.


  We study an S--I type epidemic model in an age-structured population, with
mortality due to the disease. A threshold quantity is found that controls the
stability of the disease-free equilibrium and guarantees the existence of an
endemic equilibrium. We obtain conditions on the age-dependence of the
susceptibility to infection that imply the uniqueness of the endemic
equilibrium. An example with two endemic equilibria is shown. Finally, we
analyse numerically how the stability of the endemic equilibrium is affected by
extra-mortality and by the possible periodicities induced by the demographic
age-structure.


  In a recent study, (Grigorov, 2006) analyzed temporal gene expression
profiles (Arbeitman et al., 2002) generated in a Drosophila experiment using
SSA in conjunction with Monte-Carlo SSA. The author (Grigorov, 2006) makes
three important claims in his article, namely:
  Claim1: A new method based on the theory of nonlinear time series analysis is
used to capture the global dynamics of the fruit-fly cycle temporal gene
expression profiles.
  Claim 2: Flattening of a significant part of the eigen-spectrum confirms the
hypothesis about an underly-ing high-dimensional chaotic generating process.
  Claim 3: Monte-Carlo SSA can be used to establish whether a given time series
is distinguishable from any well-defined process including deterministic chaos.
  In this report we present fundamental concerns with respect to the above
claims (Grigorov, 2006) in a systematic manner with simple examples. The
discussion provided especially discourages the choice of SSA for inferring
nonlinear dynamical structure form time series obtained in any biological
paradigm.


  In this paper, the authors develop a method of detecting correlations between
epidemic patterns in different regions that are due to human movement and
introduce a null model in which the travel-induced correlations are cancelled.
They apply this method to the well-documented cases of seasonal influenza
outbreaks in the United States and France. In the United States (using data for
1972-2002), the authors observed strong short-range correlations between
several states and their immediate neighbors, as well as robust long-range
spreading patterns resulting from large domestic air-traffic flows. The
stability of these results over time allowed the authors to draw conclusions
about the possible impact of travel restrictions on epidemic spread. The
authors also applied this method to the case of France (1984-2004) and found
that on the regional scale, there was no transportation mode that clearly
dominated disease spread. The simplicity and robustness of this method suggest
that it could be a useful tool for detecting transmission channels in the
spread of epidemics.


  Associative memory Hamiltonian structure prediction potentials are not overly
rugged, thereby suggesting their landscapes are like those of actual proteins.
In the present contribution we show how basin-hopping global optimization can
identify low-lying minima for the corresponding mildly frustrated energy
landscapes. For small systems the basin-hopping algorithm succeeds in locating
both lower minima and conformations closer to the experimental structure than
does molecular dynamics with simulated annealing. For large systems the
efficiency of basin-hopping decreases for our initial implementation, where the
steps consist of random perturbations to the Cartesian coordinates. We
implemented umbrella sampling using basin-hopping to further confirm when the
global minima are reached. We have also improved the energy surface by
employing bioinformatic techniques for reducing the roughness or variance of
the energy surface. Finally, the basin-hopping calculations have guided
improvements in the excluded volume of the Hamiltonian, producing better
structures. These results suggest a novel and transferable optimization scheme
for future energy function development.


  We assess the reliability of the recently developed approach denominated
Dominant Reaction Pathways (DRP) by studying the folding of a 16-residue
beta-hairpin, within a coarse-grained Go-type model. We show that the DRP
predictions are in quantitative agreement with the results of Molecular
Dynamics simulations, performed in the same model. On the other hand, in the
DRP approach, the computational difficulties associated to the decoupling of
time scales are rigorously bypassed. The analysis of the important transition
pathways supports a picture of the beta-hairpin folding in which the reaction
is initiated by the collapse of the hydrophobic cluster.


  The ability to absorb mutations while retaining structure and function, or
mutational robustness, is a remarkable property of natural proteins. In this
Letter, we use a computational model of organismic evolution [Zeldovich et al,
PLOS Comp Biol 3(7):e139 (2007)], which explicitly couples protein physics and
population dynamics, to study mutational robustness of evolved model proteins.
We find that dominant protein structures which evolved in the simulations are
highly designable ones, in accord with some of the earlier observations. Next,
we compare evolved sequences with the ones designed to fold into the same
dominant structures and having the same thermodynamic stability, and find that
evolved sequences are more robust against point mutations, being less likely to
be destabilized upon them. These results point to sequence evolution as an
important method of protein engineering if mutational robustness of the
artificially developed proteins is desired. On the biological side, mutational
robustness of proteins appears to be a natural consequence of the
mutation-selection evolutionary process.


  Recent in vitro experiments have shown that DNA ejection from bacteriophage
can be partially stopped by surrounding osmotic pressure when ejected DNA is
digested by DNase I on the course of ejection. We argue in this work by
combination of experimental techniques (osmotic suppression without DNaseI
monitored by UV absorbance, pulse-field electrophoresis, and cryo-EM
visualization) and simple scaling modeling that intact genome (i.e. undigested)
ejection in a crowded environment is, on the contrary, enhanced or eventually
complete with the help of a pulling force resulting from DNA condensation
induced by the osmotic stress itself. This demonstrates that in vivo, the
osmotically stressed cell cytoplasm will promote phage DNA ejection rather than
resisting it. The further addition of DNA-binding proteins under crowding
conditions is shown to enhance the extent of ejection. We also found some
optimal crowding conditions for which DNA content remaining in the capsid upon
ejection is maximum, which correlates well with the optimal conditions of
maximum DNA packaging efficiency into viral capsids observed almost 20 years
ago. Biological consequences of this finding are discussed.


  Motivated by experimental results on the interplay between molecular motors
and tau proteins, we extend lattice-based models of intracellular transport to
include a second species of particle which locally influences the
motor-filament attachment rate. We consider various exactly solvable limits of
a stochastic multi-particle model before focusing on the low-motor-density
regime. Here, an approximate treatment based on the random walk behaviour of
single motors gives good quantitative agreement with simulation results for the
tau-dependence of the motor current. Finally, we discuss the possible
physiological implications of our results.


  The plasma membrane of human red blood cells consists of a lipid bilayer
attached to a regular network of underlying cytoskeletal polymers. We model
this system at a dynamic coarse-grained level, treating the bilayer as an
elastic sheet and the cytoskeletal network as a series of phantom entropic
springs. In contrast to prior simulation efforts, we explicitly account for
dynamics of the cytoskeletal network, both via motion of the protein anchors
that attach the cytoskeleton to the bilayer and through breaking and
reconnection of individual cytoskeletal filaments. Simulation results are
explained in the context of a simple mean-field percolation model and
comparison is made to experimental measurements of red blood cell fluctuation
amplitudes.


  In advancing discrete-based computational cancer models towards clinical
applications, one faces the dilemma of how to deal with an ever growing amount
of biomedical data that ought to be incorporated eventually in one form or
another. Model scalability becomes of paramount interest. In an effort to start
addressing this critical issue, here, we present a novel multi-scale and
multi-resolution agent-based in silico glioma model. While "multi-scale" refers
to employing an epidermal growth factor receptor (EGFR)-driven molecular
network to process cellular phenotypic decisions within the micro-macroscopic
environment, "multi-resolution" is achieved through algorithms that classify
cells to either active or inactive spatial clusters, which determine the
resolution they are simulated at. The aim is to assign computational resources
where and when they matter most for maintaining or improving the predictive
power of the algorithm, onto specific tumor areas and at particular times.
Using a previously described 2D brain tumor model, we have developed four
different computational methods for achieving the multi-resolution scheme,
three of which are designed to dynamically train on the high-resolution
simulation that serves as control. To quantify the algorithms' performance, we
rank them by weighing the distinct computational time savings of the simulation
runs versus the methods' ability to accurately reproduce the high-resolution
results of the control. Finally, to demonstrate the flexibility of the
underlying concept, we show the added value of combining the two highest-ranked
methods.


  Aberrantly regulated cell motility is a hallmark of cancer cells. A hybrid
agent-based model has been developed to investigate the synergistic and
antagonistic cell motility-impacting effects of three microenvironment
variables simultaneously: chemoattraction, haptotactic permission, and
biomechanical constraint or resistance. Reflecting distinct cell-specific
intracellular machinery, the cancer cells are modelled as processing a variety
of spatial search strategies that respond to these three influencing factors
with differential weights attached to each. While responding exclusively to
chemoattraction optimizes cell displacement effectiveness, incorporating
permission and resistance components becomes increasingly important with
greater distance to the chemoattractant source and/or after reducing the
ligand's effective diffusion coefficient. Extending this to a heterogeneous
population of cells shows that displacement effectiveness increases with clonal
diversity as characterized by the Shannon index. However, the resulting data
can be fit best to an exponential function, suggesting that there is a level of
population heterogeneity beyond which its added value to the cancer system
becomes minimal as directionality ceases to increase. Possible experimental
extensions and potential clinical implications are discussed.


  We have recently introduced an efficient method for the detection and
identification of modules in complex networks, based on the de-synchronization
properties (dynamical clustering) of phase oscillators. In this paper we apply
the dynamical clustering tecnique to the identification of communities of
marine organisms living in the Chesapeake Bay food web. We show that our
algorithm is able to perform a very reliable classification of the real
communities existing in this ecosystem by using different kinds of dynamical
oscillators. We compare also our results with those of other methods for the
detection of community structures in complex networks.


  For a cell moving in hydrodynamic flow above a wall, translational and
rotational degrees of freedom are coupled by the Stokes equation. In addition,
there is a close coupling of convection and diffusion due to the
position-dependent mobility. These couplings render calculation of the mean
encounter time between cell surface receptors and ligands on the substrate very
difficult. Here we show for a two-dimensional model system how analytical
progress can be achieved by treating motion in the vertical direction by an
effective reaction term in the mean first passage time equation for the
rotational degree of freedom. The strength of this reaction term can either be
estimated from equilibrium considerations or used as a fit parameter. Our
analytical results are confirmed by computer simulations and allow to assess
the relative roles of convection and diffusion for different scaling regimes of
interest.


  It is now well established that amphibians are suffering widespread decline
and extinctions. Among other causes, urbanization is responsible for habitat
reduction, habitat fragmentation and massive road kills. In this context, it is
urgent to develop and assess appropriate conservation measures. Using yearly
censuses of migrating adults of two anuran species at one location in
Switzerland, we examined the impact of a road mitigation measure - permanent
under-road tunnels with guiding trenches - along a road separating wintering
forests from breeding wetlands. We observe that the adult migrating populations
do not exhibit any long-term trend but undergo a transient increase a few years
after the installation of the road mitigation measure. Using additional
datasets like climatic data and censuses obtained in a control area, we show
that the observed pattern of migrating populations cannot be explained by any
other data at our disposal. We then checked as a working hypothesis whether the
installation of under-road tunnels could explain the observed transient or not.
To this end, we use a simple population model and show that the road mitigation
measure together with competition for resources can successfully explain the
experimental observations. We conclude by discussing the requirements for
further assessment of this hypothesis as well as consequences for conservation
planners.


  We discuss the possibility of existence of entanglement in biological
systems. Our arguments centre on the fact that biological systems are
thermodynamic open driven systems far from equilibrium. In such systems error
correction can occur which may maintain entanglement despite high levels of
de-coherence. We also discuss the possibility of cooling (classical or quantum)
at molecular level.


  We analyse an epidemiological model of competing strains of pathogens and
hence differences in transmission for first versus secondary infection due to
interaction of the strains with previously aquired immunities, as has been
described for dengue fever (in dengue known as antibody dependent enhancement,
ADE). Such models show a rich variety of dynamics through bifurcations up to
deterministic chaos. Including temporary cross-immunity even enlarges the
parameter range of such chaotic attractors, and also gives rise to various
coexisting attractors, which are difficult to identify by standard numerical
bifurcation programs using continuation methods. A combination of techniques,
including classical bifurcation plots and Lyapunov exponent spectra has to be
applied in comparison to get further insight into such dynamical structures.
Here we present for the first time multi-parameter studies in a range of
biologically plausible values for dengue. The multi-strain interaction with the
immune system is expected to also have implications for the epidemiology of
other diseases.


  The performance of single folding predictors and combination scores is
critically evaluated. We test mean packing, mean pairwise energy and the new
index gVSL2 on a dataset of 743 folded proteins and 81 natively unfolded
proteins. These predictors have an individual performance comparable or even
better than other proposed methods. We introduce here a strictly unanimous
score S_{SU} that combines them but leaves undecided those sequences
differently classified by two single predictors. The performance of the single
predictors on a dataset purged from the proteins left unclassified by S_{SU},
significantly increases, indicating that unclassified proteins are mainly false
predictions. Amino acid composition is the main determinant considered by these
predictors, therefore unclassified proteins have a composition compatible with
both folded and unfolded status. This is why purging a dataset from these
ambiguous proteins increases the performance of single predictors. The
percentage of proteins predicted as natively unfolded by S_{SU} in the three
kingdoms are: 4.1% for Bacteria, 1.0% for Archaea and 20.0% for Eukarya;
compatible with previous determinations. Evidence is given of a scaling law
relating the number of natively unfolded proteins with the total number of
proteins in a genome; a first estimate of the critical exponent is 1.95 +- 0.21


  In parametric sequence alignment, optimal alignments of two sequences are
computed as a function of the penalties for mismatches and spaces, producing
many different optimal alignments. Here we give a 3/(2^{7/3}\pi^{2/3})n^{2/3}
+O(n^{1/3} \log n) lower bound on the maximum number of distinct optimal
alignment summaries of length-n binary sequences. This shows that the upper
bound given by Gusfield et. al. is tight over all alphabets, thereby disproving
the "square root of n conjecture". Thus the maximum number of distinct optimal
alignment summaries (i.e. vertices of the alignment polytope) over all pairs of
length-n sequences is Theta(n^{2/3}).


  In a recent paper, Michael J. Saxton proposes to interpret as anomalous
diffusion the occurrence of apparent transient sub-diffusive regimes in
mean-squared displacements (MSD) plots, calculated from experimental
trajectories of molecules diffusing in living cells, acquired by Single
Particle (or Molecule) Tracking techniques (SPT or SMT). In this Comment,
without questioning the existence of sub-diffusive behaviors, which certainly
play a key role in numbers of mechanisms in living systems, we point out that
the data used by J.M. Saxton can as well be fitted by a simple law, resulting
from confined diffusion at short times, with a slower free diffusion
superimposed at larger times. When visualizing MSD plots, the transition from
short-term diffusion confined in domains of size L, to slower, longer-term free
diffusion, can be confused with anomalous diffusion over several orders of
magnitude of time.


  Hybridization networks are representations of evolutionary histories that
allow for the inclusion of reticulate events like recombinations,
hybridizations, or lateral gene transfers. The recent growth in the number of
hybridization network reconstruction algorithms has led to an increasing
interest in the definition of metrics for their comparison that can be used to
assess the accuracy or robustness of these methods. In this paper we establish
some basic results that make it possible the generalization to tree-child time
consistent (TCTC) hybridization networks of some of the oldest known metrics
for phylogenetic trees: those based on the comparison of the vectors of path
lengths between leaves. More specifically, we associate to each hybridization
network a suitably defined vector of `splitted' path lengths between its
leaves, and we prove that if two TCTC hybridization networks have the same such
vectors, then they must be isomorphic. Thus, comparing these vectors by means
of a metric for real-valued vectors defines a metric for TCTC hybridization
networks. We also consider the case of fully resolved hybridization networks,
where we prove that simpler, `non-splitted' vectors can be used.


  Self-arrangement of individuals into spatial patterns often accompanies and
promotes species diversity in ecological systems. Here, we investigate pattern
formation arising from cyclic dominance of three species, operating near a
bifurcation point. In its vicinity, an Eckhaus instability occurs, leading to
convectively unstable "blurred" patterns. At the bifurcation point, stochastic
effects dominate and induce counterintuitive effects on diversity: Large
patterns, emerging for medium values of individuals' mobility, lead to rapid
species extinction, while small patterns (low mobility) promote diversity, and
high mobilities render spatial structures irrelevant. We provide a quantitative
analysis of these phenomena, employing a complex Ginzburg-Landau equation.


  Active centres and hot spots of proteins have a paramount importance in
enzyme action, protein complex formation and drug design. Recently a number of
publications successfully applied the analysis of residue networks to predict
active centres in proteins. Most real-world networks show a number of
properties, such as small-worldness or scale-free degree distribution, which
are rather general features of networks, from molecules to society at large.
Using analogy I propose that existing findings and methodology already enable
us to detect active centres in cells, and can be expanded to social networks
and ecosystems. Members of these active centres are termed here as creative
elements of their respective networks, which may help them to survive
unprecedented, novel challenges, and play a key role in the development,
survival and evolvability of complex systems.


  W. D. Hamilton's celebrated formula for the age-specific force of natural
selection furnishes predictions for senescent mortality due to mutation
accumulation, at the price of reliance on a linear approximation. Applying to
Hamilton's setting the full non-linear demographic model for mutation
accumulation of Evans et al. (2007), we find surprising differences. Non-linear
interactions cause the collapse of Hamilton-style predictions in the most
commonly studied case, refine predictions in other cases, and allow Walls of
Death at ages before the end of reproduction. Haldane's Principle for genetic
load has an exact but unfamiliar generalization.


  Intracellular transport is based on molecular motors that pull cargos along
cytoskeletal filaments. One motor species always moves in one direction, e.g.
conventional kinesin moves to the microtubule plus end, while cytoplasmic
dynein moves to the microtubule minus end. However, many cellular cargos are
observed to move bidirectionally, involving both plus-end and minus-end
directed motors. The presumably simplest mechanism for such bidirectional
transport is provided by a tug-of-war between the two motor species. This
mechanism is studied theoretically using the load-dependent transport
properties of individual motors as measured in single-molecule experiments. In
contrast to previous expectations, such a tug-of-war is found to be highly
cooperative and to exhibit seven different motility regimes depending on the
precise values of the single motor parameters. The sensitivity of the transport
process to small parameter changes can be used by the cell to regulate its
cargo traffic.


  Perturbation experiments are carried out by contact process and its
mean-field version. Here, the mortality rate is increased or decreased
suddenly. It is known that the fluctuation enhancement (FE) occurs after the
perturbation, where FE means a population uncertainty. In the present paper, we
develop a new theory of stochastic differential equation. The agreement between
the theory and the mean-field simulation is almost perfect. This theory enables
us to find much stronger FE than reported previously. We discuss the population
uncertainty in the recovering process of endangered species.


  I believe an atomic biology is needed to supplement present day molecular
biology, if we are to design and understand proteins, as well as define, make,
and use them. Topics in the paper are molecular biology and atomic biology.
Electrodiffusion in the open channel. Electrodiffusion in mixed electrolytes.
Models of permeation. State Models of Permeation are Inconsistent with the
Electric Field. Making models in atomic biology. Molecular dynamics. Temporal
Limitations; Spatial Limitations; Periodic boundary conditions. Hierarchy of
models of the open channel. Stochastic Motion of the Channel. Langevin
Dynamics. Simulations of the Reaction Path: the Permion. Chemical reactions.
What was wrong? Back to the hierarchy: Occam's razor can slit your throat.
Poisson-Nernst-Planck PNP Models Flux Ratios; Pumping by Field Coupling. Gating
in channels of one conformation. Gating by Field Switching; Gating Current;
Gating in Branched Channels; Blocking. Back to the hierarchy: Linking levels.
Is there a theory? At what level will the adaptation be found? Simplicity,
evolution, and natural function.


  The paper shows how the diffusive movement of ions through a channel protein
can be described as a chemical reaction over an arbitrary shaped potential
barrier. The result is simple and intuitive but without approximation beyond
the electrodiffusion description of ion movement.


  Signaling events in eukaryotic cells are often guided by a scaffolding
protein. Scaffold proteins assemble multiple proteins in a spatially localized
signaling complex and exert numerous physical effects on signaling pathways. To
study these effects, we consider a minimal, three-state kinetic model of
scaffold mediated kinase activation. We first introduce and apply a path
summation technique to obtain approximate solutions to a single molecule master
equation that governs protein kinase activation. We then consider exact
numerical solutions. We comment on when this approximation is appropriate and
then use this analysis to illustrate the competition of processes occurring at
many time scales involved in signal transduction in the presence of a scaffold
protein. The findings are consistent with recent experiments and simulation
data. Our results provide a framework and offer a mechanism for understanding
how scaffold proteins can influence the shape of the waiting time distribution
of kinase activation and effectively broaden the times over which protein
kinases are activated in the course of cell signaling.


  Eukaryote genomes contain excessively introns, inter-genic and other
non-genic sequences that appear to have no vital functional role or phenotype
manifestation. Their existence, a long-standing puzzle, is viewed from the
principle of increasing entropy. According to thermodynamics of open systems,
genomes evolve toward diversity by various mechanisms that increase, decrease
and distribute genomic material in response to thermodynamic driving forces.
Evolution results in an excessive genome, a high-entropy ecosystem of its own,
where copious non-coding segments associate with low-level functions and
conserved sequences code coordinated activities. The rate of entropy increase,
equivalent to the rate of free energy decrease, is identified with the
universal fitness criterion of natural selection that governs populations of
genomic entities as well as other species.


  Molecular motors transport various cargos along cytoskeletal filaments,
analogous to trucks on roads. In contrast to vehicles, however, molecular
motors do not work alone but in small teams. We describe a simple model for the
transport of a cargo by one team of motors and by two teams of motors, which
walk into opposite directions. The cooperation of one team of motors generates
long-range transport, which we observed experimentally in vitro. Transport by
two teams of motors leads to a variety of bidirectional motility behaviour and
to dynamic instabilities reminiscent of spontaneous symmetry breaking. We also
discuss how cargo transport by teams of motors allows the cell to generate
robust long-range bidirectional transport.


  Methods of modeling cellular regulatory networks as diverse as differential
equations and Boolean networks co-exist, however, without any closer
correspondence to each other. With the example system of the fission yeast cell
cycle control network, we here set the two approaches in relation to each
other. We find that the Boolean network can be formulated as a specific
coarse-grained limit of the more detailed differential network model for this
system. This lays the mathematical foundation on which Boolean networks can be
applied to biological regulatory networks in a controlled way.


  This review article deals with the specificities of the proteomics analysis
of membrane proteins.


  The oxidative stress response is characterized by various effects on a range
of biologic molecules. When examined at the protein level, both expression
levels and protein modifications are altered by oxidative stress. While these
effects have been studied in the past by classic biochemical methods, the
recent onset of proteomics methods has allowed the oxidative stress response to
be studied on a much wider scale. The input of proteomics in the study of
oxidative stress response and in the evidence of an oxidative stress component
in biologic phenomena is reviewed in this paper.


  Biomarker discovery for clinical purposes is one of the major areas in which
proteomics is used. However, despite considerable effort, the successes have
been relatively scarce. In this perspective paper, we try to highlight and
analyze the main causes for this limited success, and to suggest alternate
strategies, which will avoid them, without eluding the foreseeable weak points
of these strategies. Two major strategies are analyzed, namely, the switch from
body fluids to cell and tissues for the initial biomarker discovery step or, if
body fluids must be analyzed, the implementation of highly selective protein
selection strategies.


  Mitochondria are complex organelles, and their proteomics analysis requires a
combination of techniques. The emphasis in this chapter is made first on
mitochondria preparation from cultured mammalian cells, then on the separation
of the mitochondrial proteins with two-dimensional electrophoresis (2DE),
showing some adjustment over the classical techniques to improve resolution of
the mitochondrial proteins. This covers both the protein solubilization, the
electrophoretic part per se, and the protein detection on the gels, which makes
the interface with the protein identification part relying on mass
spectrometry.


  Biological functions are generated as a result of developmental dynamics that
form phenotypes governed by genotypes. The dynamical system for development is
shaped through genetic evolution following natural selection based on the
fitness of the phenotype. Here we study how this dynamical system is robust to
noise during development and to genetic change by mutation. We adopt a
simplified transcription regulation network model to govern gene expression,
which gives a fitness function. Through simulations of the network that
undergoes mutation and selection, we show that a certain level of noise in gene
expression is required for the network to acquire both types of robustness. The
results reveal how the noise that cells encounter during development shapes any
network's robustness, not only to noise but also to mutations. We also
establish a relationship between developmental and mutational robustness
through phenotypic variances caused by genetic variation and epigenetic noise.
A universal relationship between the two variances is derived, akin to the
fluctuation-dissipation relationship known in physics.


  The sequence-dependent structural variability and conformational dynamics of
DNA play pivotal roles in many biological milieus, such as in the site-specific
binding of transcription factors to target regulatory elements. To better
understand DNA structure, function, and dynamics in general, and protein-DNA
recognition in the 'kB' family of genetic regulatory elements in particular, we
performed molecular dynamics simulations of a 20-base pair DNA encompassing a
cognate kB site recognized by the proto-oncogenic 'c-Rel' subfamily of NF-kB
transcription factors. Simulations of the kB DNA in explicit water were
extended to microsecond duration, providing a broad, atomically-detailed
glimpse into the structural and dynamical behavior of double helical DNA over
many timescales. Of particular note, novel (and structurally plausible)
conformations of DNA developed only at the long times sampled in this
simulation -- including a peculiar state arising at ~ 0.7 us and characterized
by cross-strand intercalative stacking of nucleotides within a
longitudinally-sheared base pair, followed (at ~ 1 us) by spontaneous base
flipping of a neighboring thymine within the A-rich duplex. Results and
predictions from the us-scale simulation include implications for a dynamical
NF-kB recognition motif, and are amenable to testing and further exploration
via specific experimental approaches that are suggested herein.


  Neuroinflammation is one of the key pathological events involved in the
progression of brain damage caused by cerebral ischemia. Metabolism of
arachidonic acid through cyclooxygenase (COX) enzymes is known to be actively
involved in the neuroinflammatory events leading to neuronal death after
ischemia. Two isoforms of COX, termed COX-1 and COX-2, have been identified.
Unlike COX-1, COX-2 expression is dramatically induced by ischemia and appears
to be an effector of tissue damage. This review article will focus specifically
on the involvement of COX isozymes in brain ischemia. We will discuss issues
related to the biochemistry and selective pharmacological inhibition of COX
enzymes, and further refer to their expression in the brain under normal
conditions and following excitotoxicity and ischemic cerebral injury. We will
review present knowledge of the relative contribution of each COX isoform to
the brain ischemic pathology, based on data from investigations utilizing
selective COX-1/COX-2 inhibitors and genetic knockout mouse models. The
mechanisms of neurotoxicity associated with increased COX activity after
ischemia will also be examined. Finally, we will provide a critical evaluation
of the therapeutic potential of COX inhibitors in cerebral ischemia and discuss
new targets downstream of COX with potential neuroprotective ability.


  An analysis of traveling wave solutions of pure cross-diffusion systems,
i.e., systems that lack reaction and self-diffusion terms, is presented. Using
the qualitative theory of phase plane analysis the conditions for existence of
different types of wave solutions are formulated. In particular, it is shown
that family of wave trains is a generic phenomenon in pure cross-diffusion
systems. The results can be used for construction and analysis of different
mathematical models describing systems with directional movement.


  We explore how inherent flexibility of a protein molecule influences the
mechanism controlling the kinetics of allosteric transitions using a
variational model inspired from work in protein folding. The striking
differences in the predicted transition mechanism for the opening of the two
domains of calmodulin (CaM) emphasizes that inherent flexibility is key to
understanding the complex conformational changes that occur in proteins. In
particular, the C-terminal domain of CaM (cCaM) which is inherently less
flexible than its N-terminal domain (nCaM) reveals "cracking" or local partial
unfolding during the open/closed transition. This result is in harmony with the
picture that cracking relieves local stresses due to conformational
deformations of a sufficiently rigid protein. We also compare the
conformational transition in a recently studied "even-odd" paired fragment of
CaM. Our results rationalize the different relative binding affinities of the
EF-hands in the engineered fragment compared to the intact "odd-even" paired
EF-hands (nCaM and cCaM) in terms of changes in flexibility along the
transition route. Aside from elucidating general theoretical ideas about the
cracking mechanism, these studies also emphasize how the remarkable intrinsic
plasticity of CaM underlies conformational dynamics essential for its diverse
functions.


  Phylogenetic trees describe the evolutionary history of a group of
present-day species from a common ancestor. These trees are typically
reconstructed from aligned DNA sequence data. In this paper we analytically
address the following question: is the amount of sequence data required to
accurately reconstruct a tree significantly more than the amount required to
test whether or not a candidate tree was the `true' tree? By `significantly',
we mean that the two quantities behave the same way as a function of the number
of species being considered. We prove that, for a certain type of model, the
amount of information required is not significantly different; while for
another type of model, the information required to test a tree is independent
of the number of leaves, while that required to reconstruct it grows with this
number. Our results combine probabilistic and combinatorial arguments.