  During batch growth on mixtures of two growth-limiting substrates, microbes
consume the substrates either sequentially or simultaneously. These growth
patterns are manifested in all types of bacteria and yeasts. The ubiquity of
these growth patterns suggests that they are driven by a universal mechanism
common to all microbial species. In previous work, we showed that a minimal
model accounting only for enzyme induction and dilution explains the phenotypes
observed in batch cultures of various wild-type and mutant/recombinant cells.
Here, we examine the extension of the minimal model to continuous cultures. We
show that: (1) Several enzymatic trends, usually attributed to specific
regulatory mechanisms such as catabolite repression, are completely accounted
for by dilution. (2) The bifurcation diagram of the minimal model for
continuous cultures, which classifies the substrate consumption pattern at any
given dilution rate and feed concentrations, provides a a precise explanation
for the empirically observed correlation between the growth patterns in batch
and continuous cultures. (3) Numerical simulations of the model are in
excellent agreement with the data. The model captures the variation of the
steady state substrate concentrations, cell densities, and enzyme levels during
the single- and mixed-substrate growth of bacteria and yeasts at various
dilution rates and feed concentrations. (4) This variation is well-approximated
by simple analytical expressions that furnish physical insights into the steady
states of continuous cultures. The minimal model provides a framework for
quantitating the effect of regulatory mechanisms. We illustrate this by
analyzing several data sets from the literature.


  We study the evolution of the probability density of an asexual, one locus
population under natural selection and random evolution. This evolution is
governed by a Fokker-Planck equation with degenerate coefficients on the
boundaries, supplemented by a pair of conservation laws. It is readily shown
that no classical or standard weak solution definition yields solvability of
the problem. We provide an appropriate definition of weak solution for the
problem, for which we show existence and uniqueness. The solution displays a
very distinctive structure and, for large time, we show convergence to a unique
stationary solution that turns out to be a singular measure supported at the
endpoints. An exponential rate of convergence to this steady state is also
proved.


  Gene activation in eukaryotes involves the concerted action of histone tail
modifiers, chromatin remodellers and transcription factors, whose precise
coordination is currently unknown. We demonstrate that the experimentally
observed interactions of the molecules are in accord with a kinetic
proofreading scheme. Our finding could provide a basis for the development of
quantitative models for gene regulation in eukaryotes based on the
combinatorical interactions of chromatin modifiers.


  Mechanochemical coupling was studied for two different types of myosin motors
in cells: myosin V, which carries cargo over long distances by as a single
molecule; and myosin II, which generates a contracting force in cooperation
with other myosin II molecules. Both mean and variance of myosin V velocity at
various [ATP] obeyed Michaelis-Menten mechanics, consistent with tight
mechanochemical coupling. Myosin II, working in an ensemble, however, was
explained by a loose coupling mechanism, generating variable step sizes
depending on the ATP concentration and realizing a much larger step (200 nm)
per ATP hydrolysis than myosin V through its cooperative nature at zero load.
These different mechanics are ideal for the respective myosin's physiological
functions.


  In this paper, we propose a thermodynamic mechanism for the formation of
transcriptional foci via the joint agglomeration of DNA-looping proteins and
protein-binding domains on DNA: The competition between the gain in protein-DNA
binding free energy and the entropy loss due to DNA looping is argued to result
in an effective attraction between loops. A mean-field approximation can be
described analytically via a mapping to a restricted random-graph ensemble
having local degree constraints and global constraints on the number of
connected components. It shows the emergence of protein clusters containing a
finite fraction of all looping proteins. If the entropy loss due to a single
DNA loop is high enough, this transition is found to be of first order.


  We argue that fractional order (FO) differential equations are more suitable
to model complex adaptive systems (CAS). Hence they are applied in replicator
equations for non-cooperative game. Rock-Scissors-Paper game is discussed. It
is known that its integer order model does not have a stable equilibrium. Its
fractional order model is shown to have a locally asymptotically stable
internal solution. A FO asymmetric game is shown to have a locally
asymptotically stable internal solution. This is not the case for its integer
order counterpart.


  Coevolutionary dynamics is investigated in chemical catalysis, biological
evolution, social and economic systems. The dynamics of these systems can be
analyzed within the unifying framework of evolutionary game theory. In this
Letter, we show that even in well-mixed finite populations, where the dynamics
is inherently stochastic, biodiversity is possible with three cyclic dominant
strategies. We show how the interplay of evolutionary dynamics, discreteness of
the population, and the nature of the interactions influences the coexistence
of strategies. We calculate a critical population size above which coexistence
is likely.


  The formation of out-of-equilibrium patterns is a characteristic feature of
spatially-extended, biodiverse, ecological systems. Intriguing examples are
provided by cyclic competition of species, as metaphorically described by the
`rock-paper-scissors' game. Both experimentally and theoretically, such
non-transitive interactions have been found to induce self-organization of
static individuals into noisy, irregular clusters. However, a profound
understanding and characterization of such patterns is still lacking. Here, we
theoretically investigate the influence of individuals' mobility on the spatial
structures emerging in rock-paper-scissors games. We devise a quantitative
approach to analyze the spatial patterns self-forming in the course of the
stochastic time evolution. For a paradigmatic model originally introduced by
May and Leonard, within an interacting particle approach, we demonstrate that
the system's behavior - in the proper continuum limit - is aptly captured by a
set of stochastic partial differential equations. The system's stochastic
dynamics is shown to lead to the emergence of entangled rotating spiral waves.
While the spirals' wavelength and spreading velocity is demonstrated to be
accurately predicted by a (deterministic) complex Ginzburg-Landau equation,
their entanglement results from the inherent stochastic nature of the
system.These findings and our methods have important applications for
understanding the formation of noisy patterns, e.g., in ecological and
evolutionary contexts, and are also of relevance for the kinetics of
(bio)-chemical reactions.


  We study metapopulation models for the spread of epidemics in which different
subpopulations (cities) are connected by fluxes of individuals (travelers).
This framework allows to describe the spread of a disease on a large scale and
we focus here on the computation of the arrival time of a disease as a function
of the properties of the seed of the epidemics and of the characteristics of
the network connecting the various subpopulations. Using analytical and
numerical arguments, we introduce an easily computable quantity which
approximates this average arrival time. We show on the example of a disease
spread on the world-wide airport network that this quantity predicts with a
good accuracy the order of arrival of the disease in the various subpopulations
in each realization of epidemic scenario, and not only for an average over
realizations. Finally, this quantity might be useful in the identification of
the dominant paths of the disease spread.


  Molecular motors are macromolecular complexes which use some form of input
energy to perform mechanical work. The filamentary tracks, on which these
motors move, are made of either proteins (e.g., microtubules) or nucleic acids
(DNA or RNA). Often, many such motors move simultaneously on the same track and
their collective properties have superficial similarities with vehicular
traffic on highways. The models we have developed provide ``unified''
description: in the low-density limit, a model captures the transport
properties of a single motor while, at higher densities the same model accounts
for the collective spatio-temporal organization of interacting motors. By
drawing analogy with vehicular traffic, we have introduced novel quantities for
characterizing the nature of the spatio-temporal organization of molecular
motors on their tracks. We show how the traffic-like intracellular collective
phenomena depend on the mechano-chemistry of the corresponding individual
motors.


  We prove that Nakhleh's latest dissimilarity measure for phylogenetic
networks separates distinguishable phylogenetic networks, and that a slight
modification of it provides a true distance on the class of all phylogenetic
networks.


  This Letter studies the quasispecies dynamics of a population capable of
genetic repair evolving on a time-dependent fitness landscape. We develop a
model that considers an asexual population of single-stranded, conservatively
replicating genomes, whose only source of genetic variation is due to copying
errors during replication. We consider a time-dependent, single-fitness-peak
landscape where the master sequence changes by a single point mutation every
time $ \tau $. We are able to analytically solve for the evolutionary dynamics
of the population in the point-mutation limit. In particular, our model
provides an analytical expression for the fraction of mutators in the dynamic
fitness landscape that agrees well with results from stochastic simulations.


  This paper deals with temporal enzyme distribution in the activation of
biochemical pathways. Pathway activation arises when production of a certain
biomolecule is required due to changing environmental conditions. Under the
premise that biological systems have been optimized through evolutionary
processes, a biologically meaningful optimal control problem is posed. In this
setup, the enzyme concentrations are assumed to be time dependent and
constrained by a limited overall enzyme production capacity, while the
optimization criterion accounts for both time and resource usage.
  Using geometric arguments we establish the bang-bang nature of the solution
and reveal that each reaction must be sequentially activated in the same order
as they appear in the pathway. The results hold for a broad range of enzyme
dynamics which includes, but is not limited to, Mass Action, Michaelis-Menten
and Hill Equation kinetics.


  In Soda'06, Chaudhuri, Chen, Mihaescu and Rao study algorithmic properties of
the tandem duplication - random loss model of genome rearrangement, well-known
in evolutionary biology. In their model, the cost of one step of
duplication-loss of width k is $\alpha^k$ for $\alpha =1$ or $\alpha >=2 $. In
this paper, we study a variant of this model, where the cost of one step of
width $k$ is 1 if $k <= K$ and $\infty$ if $k > K$, for any value of the
parameter $K in N$. We first show that permutations obtained after $p$ steps of
width $K$ define classes of pattern-avoiding permutations. We also compute the
numbers of duplication-loss steps of width $K$ necessary and sufficient to
obtain any permutation of $S_n$, in the worst case and on average. In this
second part, we may also consider the case $K=K(n)$, a function of the size $n$
of the permutation on which the duplication-loss operations are performed.


  We classify and predict the asymptotic dynamics of a class of swarming
models. The model consists of a conservation equation in one dimension
describing the movement of a population density field. The velocity is found by
convolving the density with a kernel describing attractive-repulsive social
interactions. The kernel's first moment and its limiting behavior at the origin
determine whether the population asymptotically spreads, contracts, or reaches
steady-state. For the spreading case, the dynamics approach those of the porous
medium equation. The widening, compactly-supported population has edges that
behave like traveling waves whose speed, density and slope we calculate. For
the contracting case, the dynamics of the cumulative density approach those of
Burgers' equation. We derive an analytical upper bound for the finite blow-up
time after which the solution forms one or more $\delta$-functions.


  The flexibility in gap cost enjoyed by Hidden Markov Models (HMMs) is
expected to afford them better retrieval accuracy than position-specific
scoring matrices (PSSMs). We attempt to quantify the effect of more general gap
parameters by separately examining the influence of position- and
composition-specific gap scores, as well as by comparing the retrieval accuracy
of the PSSMs constructed using an iterative procedure to that of the HMMs
provided by Pfam and SUPERFAMILY, curated ensembles of multiple alignments.
  We found that position-specific gap penalties have an advantage over uniform
gap costs. We did not explore optimizing distinct uniform gap costs for each
query. For Pfam, PSSMs iteratively constructed from seeds based on HMM
consensus sequences perform equivalently to HMMs that were adjusted to have
constant gap transition probabilities, albeit with much greater variance. We
observed no effect of composition-specific gap costs on retrieval performance.


  Most conspicuous organisms are multicellular and most multicellular organisms
develop somatic cells to perform specific, non-reproductive tasks. The ubiquity
of this division of labor suggests that it is highly advantageous. In this
paper, I present a model to study the evolution of specialized cells. The model
allows for unicellular and multicellular organisms that may contain somatic
(terminally differentiated) cells. Cells contribute additively to a
quantitative trait. The fitness of the organism depends on this quantitative
trait (via a benefit function), the size of the organism, and the number of
somatic cells. This model allows one to determine when somatic cells are
advantageous and to calculate the optimum number (or fraction) of reproductive
cells. I show that the fraction of reproductive cells is always surprisingly
high. If somatic cells are very small they can outnumber reproductive cells but
their biomass is still less than the biomass of reproductive cells. Only for
convex benefit functions can the biomass of somatic cell exceed the biomass of
reproductive cells. I discuss the biology of primitive multicellular organisms
with respect to the model predictions. I find good agreement and outline how
this work can be used to guide further quantitative studies of
multicellularity.


  We consider a simple-model population, whose individuals react with a certain
delay to temporal variations of their habitat. We investigate the impact of
such a delayed-answer on the survival chances of the population, both in a
periodically changing environment, and in the case of an abrupt change of it.
It is found that for population with low degree of mutation-induced
variability, being "slow-reacting" decreases the extinction risk face to
environmental changes. On the contrary, for populations with high mutation
amplitude, the delayed reaction reduces the survival chances.


  In mammals, dosage compensation of X linked genes in female cells is achieved
by inactivation of one of their two X chromosomes which is randomly chosen. The
earliest steps in X-inactivation (XCI), namely the mechanism whereby cells
count their X chromosomes and choose between two equivalent X, remain
mysterious. Starting from the recent discovery of X chromosome colocalization
at the onset of X-inactivation, we propose a Statistical Mechanics model of
XCI, which is investigated by computer simulations and checked against
experimental data. Our model describes how a `blocking factor' complex is
self-assembled and why only one is formed out of many diffusible molecules,
resulting in a spontaneous symmetry breaking (SB) in the binding to two
identical chromosomes. These results are used to derive a scenario of
biological implications describing all current experimental evidences, e.g.,
the importance of colocalization.


  Cyclic predator-prey models with four or six species are studied on a square
lattice when the invasion rates are varied. It is found that the cyclic
invasions maintain a self-organizing pattern as long as the deviation of the
invasion rate(s) from a uniform value does not exceed a threshold value. For
larger deviations the system exhibits a continuous phase transition into a
frozen distribution of odd (or even) label species.


  Covarion models of character evolution describe inhomogeneities in
substitution processes through time. In phylogenetics, such models are used to
describe changing functional constraints or selection regimes during the
evolution of biological sequences. In this work the identifiability of such
models for generic parameters on a known phylogenetic tree is established,
provided the number of covarion classes does not exceed the size of the
observable state space. `Generic parameters' as used here means all parameters
except possibly those in a set of measure zero within the parameter space.
Combined with earlier results, this implies both the tree and generic numerical
parameters are identifiable if the number of classes is strictly smaller than
the number of observable states.


  Motivation: Array-based comparative genomic hybridization (arrayCGH) has
recently become a popular tool to identify DNA copy number variations along the
genome. These profiles are starting to be used as markers to improve prognosis
or diagnosis of cancer, which implies that methods for automated supervised
classification of arrayCGH data are needed. Like gene expression profiles,
arrayCGH profiles are characterized by a large number of variables usually
measured on a limited number of samples. However, arrayCGH profiles have a
particular structure of correlations between variables, due to the spatial
organization of BACs along the genome. This suggests that classical
classification methods, often based on the selection of a small number of
discriminative features, may not be the most accurate methods and may not
produce easily interpretable prediction rules.
  Results: We propose a new method for supervised classification of arrayCGH
data. The method is a variant of support vector machine (SVM) that incorporates
the biological specificities of DNA copy number variations along the genome as
prior knowledge. The resulting classifier is a sparse linear classifier based
on a limited number of regions automatically selected on the chromosomes,
leading to easy interpretation and identification of discriminative regions of
the genome. We test this method on three classification problems for bladder
and uveal cancer, involving both diagnosis and prognosis. We demonstrate that
the introduction of the new prior on the classifier leads not only to more
accurate predictions, but also to the identification of known and new regions
of interest in the genome.
  Availability: All data and algorithms are publicly available.


  We consider an elastic rod model for twisted DNA in the plectonemic regime.
The molecule is treated as an impenetrable tube with an effective, adjustable
radius. The model is solved analytically and we derive formulas for the contact
pressure, twisting moment and geometrical parameters of the supercoiled region.
We apply our model to magnetic tweezer experiments of a DNA molecule subjected
to a tensile force and a torque, and extract mechanical and geometrical
quantities from the linear part of the experimental response curve. These
reconstructed values are derived in a self-contained manner, and are found to
be consistent with those available in the literature.


  We present a new analysis of relationships between disease incidence and the
prevalence of an experimentally defined state of `recent infection'. This leads
to a clean separation between biological parameters (properties of disease
progression as reflected in a test for recent infection), which need to be
calibrated, and epidemiological state variables, which are estimated in a
cross-sectional survey. The framework takes into account the possibility that
details of the assay and host/pathogen chemistry leave a (knowable) fraction of
the population in the recent category for all times. This systematically
addresses an issue which is the source of some controversy about the
appropriate use of the BED assay for defining recent HIV infection. Analysis of
relative contributions of error arising variously from statistical
considerations and simplifications of general expressions indicate that
statistical error dominates heavily over all sources of bias for realistic
epidemiological and biological scenarios. Numerical calculations validate the
approximations made in analytical relations.


  Nicodemi and Prisco recently proposed a model for X-chromosome inactivation
in mammals, explaining this phenomenon in terms of a spontaneous
symmetry-breaking mechanism [{\it Phys. Rev. Lett.} 99 (2007), 108104]. Here we
provide a mean-field version of their model.


  Background: The computation of the statistical properties of motif
occurrences has an obviously relevant practical application: for example,
patterns that are significantly over- or under-represented in the genome are
interesting candidates for biological roles. However, the problem is
computationally hard; as a result, virtually all the existing pipelines use
fast but approximate scoring functions, in spite of the fact that they have
been shown to systematically produce incorrect results. A few interesting exact
approaches are known, but they are very slow and hence not practical in the
case of realistic sequences. Results: We give an exact solution, solely based
on deterministic finite-state automata (DFAs), to the problem of finding not
only the p-value, but the whole relevant part of the Markovian probability
distribution function of a motif in a biological sequence. In particular, the
time complexity of the algorithm in the most interesting regimes is far better
than that of Nuel (2006), which was the fastest similar exact algorithm known
to date; in many cases, even approximate methods are outperformed. Conclusions:
DFAs are a standard tool of computer science for the study of patterns, but so
far they have been sparingly used in the study of biological motifs. Previous
works do propose algorithms involving automata, but there they are used
respectively as a first step to build a Finite Markov Chain Imbedding (FMCI),
or to write a generating function: whereas we only rely on the concept of DFA
to perform the calculations. This innovative approach can realistically be used
for exact statistical studies of very long genomes and protein sequences, as we
illustrate with some examples on the scale of the human genome.


  We present a theoretical investigation on possible selection of olfactory
receptors (ORs) as sensing components of nanobiosensors. Accordingly, we
generate the impedance spectra of the rat OR I7 in the native and activated
state and analyze their differences. In this way, we connect the protein
morphological transformation, caused by the sensing action, with its change of
electrical impedance. The results are compared with those obtained by studying
the best known protein of the GPCR family: bovine rhodopsin. Our investigations
indicate that a change in morphology goes with a change in impedance spectrum
mostly associated with a decrease of the static impedance up to about 60 % of
the initial value, in qualitative agreement with existing experiments on rat OR
I7. The predictiveness of the model is tested successfully for the case of
recent experiments on bacteriorhodopsin. The present results point to a
promising development of a new class of nanobiosensors based on the electrical
properties of GPCR and other sensing proteins.


  Using a combination of theory and computer simulations, we study the
translocation of an RNA molecule, pulled through a solid-state nanopore by an
optical tweezer, as a method to determine its secondary structure. The
resolution with which the elements of the secondary structure can be determined
is limited by thermal fluctuations. We present a detailed study of these
thermal fluctuations, including the frequency spectrum, and show that these
rule out single-nucleotide resolution under the experimental conditions which
we simulated. Two possible ways to improve this resolution are strong
stretching of the RNA with a back-pulling voltage across the membrane, and
stiffening of the translocated part of the RNA by biochemical means.


  This paper has been withdrawn by the author(s), due to the requirement of the
journal it currently submitted to


  We introduce a new phylogenetic reconstruction algorithm which, unlike most
previous rigorous inference techniques, does not rely on assumptions regarding
the branch lengths or the depth of the tree. The algorithm returns a forest
which is guaranteed to contain all edges that are: 1) sufficiently long and 2)
sufficiently close to the leaves. How much of the true tree is recovered
depends on the sequence length provided. The algorithm is distance-based and
runs in polynomial time.


  Interaction of a drug or chemical with a biological system can result in a
gene-expression profile or signature characteristic of the event. Using a
suitably robust algorithm these signatures can potentially be used to connect
molecules with similar pharmacological or toxicological properties. The
Connectivity Map was a novel concept and innovative tool first introduced by
Lamb et al to connect small molecules, genes, and diseases using genomic
signatures [Lamb et al (2006), Science 313, 1929-1935]. However, the
Connectivity Map had some limitations, particularly there was no effective
safeguard against false connections if the observed connections were considered
on an individual-by-individual basis. Further when several connections to the
same small-molecule compound were viewed as a set, the implicit null hypothesis
tested was not the most relevant one for the discovery of real connections.
Here we propose a simple and robust method for constructing the reference
gene-expression profiles and a new connection scoring scheme, which importantly
allows the valuation of statistical significance of all the connections
observed. We tested the new method with the two example gene-signatures (HDAC
inhibitors and Estrogens) used by Lamb et al and also a new gene signature of
immunosuppressive drugs. Our testing with this new method shows that it
achieves a higher level of specificity and sensitivity than the original
method. For example, our method successfully identified raloxifene and
tamoxifen as having significant anti-estrogen effects, while Lamb et al's
Connectivity Map failed to identify these. With these properties our new method
has potential use in drug development for the recognition of pharmacological
and toxicological properties in new drug candidates.


  In this paper, we focus on a spatial Holling-type IV predator-prey model
which contains some important factors, such as diffusion, noise (random
fluctuations) and external periodic forcing. By a brief stability and
bifurcation analysis, we arrive at the Hopf and Turing bifurcation surface and
derive the symbolic conditions for Hopf and Turing bifurcation in the spatial
domain. Based on the stability and bifurcation analysis, we obtain spiral
pattern formation via numerical simulation. Additionally, we study the model
with colored noise and external periodic forcing. From the numerical results,
we know that noise or external periodic forcing can induce instability and
enhance the oscillation of the species, and resonant response. Our results show
that modeling by reaction-diffusion equations is an appropriate tool for
investigating fundamental mechanisms of complex spatiotemporal dynamics.


  We investigate the delocalization transition appearing in an exclusion
process with two internal states resp. on two parallel lanes. At the
transition, delocalized domain walls form in the density profiles of both
internal states, in agreement with a mean-field approach. Remarkably, the
topology of the system's phase diagram allows for the delocalization of a
(localized) domain wall when approaching the transition. We quantify the domain
wall's delocalization close to the transition by analytic results obtained
within the framework of the domain wall picture. Power-law dependences of the
domain wall width on the distance to the delocalization transition as well as
on the system size are uncovered, they agree with numerical results.


  Formation and competition of associations are studied in a six-species
ecological model where each species has two predators and two prey. Each site
of a square lattice is occupied by an individual belonging to one of the six
species. The evolution of the spatial distribution of species is governed by
iterated invasions between the neighboring predator-prey pairs with species
specific rates and by site exchange between the neutral pairs with a
probability $X$. This dynamical rule yields the formation of five associations
composed of two or three species with proper spatiotemporal patterns. For large
$X$ a cyclic dominance can occur between the three two-species associations
whereas one of the two three-species associations prevails in the whole system
for low values of $X$ in the final state. Within an intermediate range of $X$
all the five associations coexist due to the fact that cyclic invasions between
the two-species associations reduce their resistance temporarily against the
invasion of three-species associations.


  Boolean networks have been the object of much attention, especially since S.
Kauffman proposed them in the 1960's as models for gene regulatory networks.
These systems are characterized by being defined on a Boolean state space and
by simultaneous updating at discrete time steps. Of particular importance for
biological applications are networks in which the indegree for each variable is
bounded by a fixed constant, as was stressed by Kauffman in his original
papers.
  An important question is which conditions on the network topology can rule
out exponentially long periodic orbits in the system. In this paper we consider
cooperative systems, i.e. systems with positive feedback interconnections among
all variables, which in a continuous setting guarantees a very stable dynamics.
In Part I of this paper we presented a construction that shows that for an
arbitrary constant 0<c<2 and sufficiently large n there exist n-dimensional
Boolean cooperative networks in which both the indegree and outdegree of each
for each variable is bounded by two (bi-quadratic networks) and which
nevertheless contain periodic orbits of length at least c^n.
  In this part, we prove an inverse result showing that for sufficiently large
n and for 0<c<2 sufficiently close to 2, any n-dimensional cooperative,
bi-quadratic Boolean network with a cycle of length at least c^n must have a
large proportion of variables with indegree 1. Such systems therefore share a
structural similarity to the systems constructed in Part I.


  A mutator is an allele that increases the mutation rate throughout the genome
by disrupting some aspect of DNA replication or repair. Mutators that increase
the mutation rate by the order of 100 fold have been observed to spontaneously
emerge and achieve high frequencies in natural populations and in long-term
laboratory evolution experiments with \textit{E. coli}. In principle, the
fixation of mutator alleles is limited by (i) competition with mutations in
wild-type backgrounds, (ii) additional deleterious mutational load, and (iii)
random genetic drift. Using a multiple locus model and employing both
simulation and analytic methods, we investigate the effects of these three
factors on the fixation probability $P_{fix}$ of an initially rare mutator as a
function of population size $N$, beneficial and deleterious mutation rates, and
the strength of mutations $s$. Our diffusion based approximation for $P_{fix}$
successfully captures effects (ii) and (iii) when selection is fast compared to
mutation ($\mu/s \ll 1$). This enables us to predict the conditions under which
mutators will be evolutionarily favored. Surprisingly, our simulations show
that effect (i) is typically small for strong-effect mutators. Our results
agree semi-quantitatively with existing laboratory evolution experiments and
suggest future experimental directions.


  We develop a theory of first passage processes in stochastic non-equilibrium
systems of birth-death type using two closely related epidemiological models as
examples. Our method employs the probability generating function technique in
conjunction with the eikonal approximation. In this way the problem is reduced
to finding the optimal path to extinction: a heteroclinic trajectory of an
effective multi-dimensional classical Hamiltonian system. We compute this
trajectory and mean extinction time of the disease numerically and uncover a
non-monotone, spiral path to extinction of a disease. We also obtain analytical
results close to a bifurcation point, where the problem is described by a
Hamiltonian previously identified in one-species population models.


  We investigate stochastic extinction in an epidemic model and the impact of
random vaccinations in large populations. We show that, in the absence of
vaccinations, the effective entropic barrier for extinction displays scaling
with the distance to the bifurcation point, with an unusual critical exponent.
Even a comparatively weak Poisson-distributed vaccination leads to an
exponential increase in the extinction rate, with the exponent that strongly
depends on the vaccination parameters.


  How molecular motors like Kinesin regulates the affinity to the rail protein
in the process of ATP hydrolysis remains to be uncovered. To understand the
regulation mechanism, we investigate the structural fluctuation of KIF1A in
different nucleotide states that are realized in the ATP hydrolysis process by
molecular dynamics simulations of Go-like model. We found that "alpha4 helix",
which is a part of the microtubule (MT) binding site, changes its fluctuation
systematically according to the nucleotide states. In particular, the frequency
of large fluctuations of alpha4 strongly correlates with the affinity of KIF1A
for microtubule. We also show how the strength of the thermal fluctuation and
the interaction with the nucleotide affect the dynamics of microtubule binding
site. These results suggest that KIF1A regulates the affinity to MT by changing
the flexibility of alpha4 helix according to the nucleotide states.


  A theoretical framework is developed to study the dynamics of protein
folding. The key insight is that the search for the native protein conformation
is influenced by the rate r at which external parameters, such as temperature,
chemical denaturant or pH, are adjusted to induce folding. A theory based on
this insight predicts that (1) proteins with non-funneled energy landscapes can
fold reliably to their native state, (2) reliable folding can occur as an
equilibrium or out-of-equilibrium process, and (3) reliable folding only occurs
when the rate r is below a limiting value, which can be calculated from
measurements of the free energy. We test these predictions against numerical
simulations of model proteins with a single energy scale.


  We present a Markov chain Monte Carlo algorithm for local alignments of
nucleotide sequences aiming to infer putative transcription factor binding
sites, referred to as the quadratic programming sampler. The new motif finder
incorporates detailed biophysical modeling of the transcription factor binding
site recognition which arises an intrinsic threshold discriminating putative
binding sites from other/background sequences. We validate the principal
functioning of the algorithm on a sample of four promoter regions from
Escherichia coli. The resulting description of the motif can be readily
evaluated on the whole genome to identify new putative binding sites.


  Bubbles in ion channel proteins have been proposed to be the bistable gates
that control current flow. Gating currents associated with channel gating would
then be an electrical signature of bubble breaking and formation, arising from
the change in dielectric coefficient as the bubble breaks or forms. A bubble
would have a dielectric coefficient of 1. A filled bubble would have a
dielectric coefficient (say) between 30 and 80. Transporters, pumps, and
channels would be expected to have gating currents.


  We predict various detectable mechanical responses to the presence of local
DNA defects which are defined as short DNA segments exhibiting mechanical
properties obviously different from the 50 nm persistence length based
semiflexible polymer model. The defects discussed are kinks and flexible hinges
either permanently fixed on DNA or thermally excited. Their effects on
extension shift, the effective persistence length, the end-to-end distance
distribution, and the cyclization probability are computed using a
transfer-matrix method. Our predictions will be useful in future experimental
designs to study DNA nicks or mismatch base pairs, mechanics of specific DNA
sequences, and specific DNA-protein interaction using magnetic tweezer,
fluorescence resonance energy transfer or plasmon resonance technique, and the
traditional biochemistry cyclization probability measurements.


  Despite the spontaneity of some in vitro protein folding reactions, native
folding in vivo often requires the participation of barrel-shaped multimeric
complexes known as chaperonins. Although it has long been known that chaperonin
substrates fold upon sequestration inside the chaperonin barrel, the precise
mechanism by which confinement within this space facilitates folding remains
unknown. In this study, we examine the possibility that the chaperonin mediates
a favorable reorganization of the solvent for the folding reaction. We begin by
discussing the effect of electrostatic charge on solvent-mediated hydrophobic
forces in an aqueous environment. Based on these initial physical arguments, we
construct a simple, phenomenological theory for the thermodynamics of density
and hydrogen bond order fluctuations in liquid water. Within the framework of
this model, we investigate the effect of confinement within a chaperonin-like
cavity on the configurational free energy of water by calculating solvent free
energies for cavities corresponding to the different conformational states in
the ATP- driven catalytic cycle of the prokaryotic chaperonin GroEL. Our
findings suggest that one function of chaperonins may be to trap unfolded
proteins and subsequently expose them to a micro-environment in which the
hydrophobic effect, a crucial thermodynamic driving force for folding, is
enhanced.


  Proteins are essential macromolecules of life and thus understanding their
function is of great importance. The number of functionally unclassified
proteins is large even for simple and well studied organisms such as baker's
yeast. Methods for determining protein function have shifted their focus from
targeting specific proteins based solely on sequence homology to analyses of
the entire proteome based on protein-protein interaction (PPI) networks. Since
proteins aggregate to perform a certain function, analyzing structural
properties of PPI networks may provide useful clues about the biological
function of individual proteins, protein complexes they participate in, and
even larger subcellular machines.
  We design a sensitive graph theoretic method for comparing local structures
of node neighborhoods that demonstrates that in PPI networks, biological
function of a node and its local network structure are closely related. The
method groups topologically similar proteins under this measure in a PPI
network and shows that these protein groups belong to the same protein
complexes, perform the same biological functions, are localized in the same
subcellular compartments, and have the same tissue expressions. Moreover, we
apply our technique on a proteome-scale network data and infer biological
function of yet unclassified proteins demonstrating that our method can provide
valuable guidelines for future experimental research.


  Ancestral maximum likelihood (AML) is a method that simultaneously
reconstructs a phylogenetic tree and ancestral sequences from extant data
(sequences at the leaves). The tree and ancestral sequences maximize the
probability of observing the given data under a Markov model of sequence
evolution, in which branch lengths are also optimized but constrained to take
the same value on any edge across all sequence sites. AML differs from the more
usual form of maximum likelihood (ML) in phylogenetics because ML averages over
all possible ancestral sequences. ML has long been known to be statistically
consistent -- that is, it converges on the correct tree with probability
approaching 1 as the sequence length grows. However, the statistical
consistency of AML has not been formally determined, despite informal remarks
in a literature that dates back 20 years. In this short note we prove a general
result that implies that AML is statistically inconsistent. In particular we
show that AML can `shrink' short edges in a tree, resulting in a tree that has
no internal resolution as the sequence length grows. Our results apply to any
number of taxa.


  In this work it is shown that 20 canonical amino acids (AAs) within genetic
code appear to be a whole system with strict distinction in Genetic Code Table
(GCT) into some different quantums: 20, 23, 61 amino acid molecules. These
molecules distinction is followed by specific balanced atom number and/or
nucleon number distinctions within those molecules. In this second version two
appendices are added; also a new version of Periodic system of numbers, whose
first verson is given in arXiv:1107.1998 [q-bio.OT].


  Background: It is of biological interest to make genome-wide predictions of
the locations of DNA melting bubbles using statistical mechanics models.
Computationally, this poses the challenge that a generic search through all
combinations of bubble starts and ends is quadratic.
  Results: An efficient algorithm is described, which shows that the time
complexity of the task is O(NlogN) rather than quadratic. The algorithm
exploits that bubble lengths may be limited, but without a prior assumption of
a maximal bubble length. No approximations, such as windowing, have been
introduced to reduce the time complexity. More than just finding the bubbles,
the algorithm produces a stitch profile, which is a probabilistic graphical
model of bubbles and helical regions. The algorithm applies a probability peak
finding method based on a hierarchical analysis of the energy barriers in the
Poland-Scheraga model.
  Conclusions: Exact and fast computation of genomic stitch profiles is thus
feasible. Sequences of several megabases have been computed, only limited by
computer memory. Possible applications are the genome-wide comparisons of
bubbles with promotors, TSS, viral integration sites, and other melting-related
regions.


  We have developed a mathematical model of regulation of expression of the
Escherichia coli lac operon, and have investigated bistability in its
steady-state induction behavior in the absence of external glucose. Numerical
analysis of equations describing regulation by artificial inducers revealed two
natural bistability parameters that can be used to control the range of inducer
concentrations over which the model exhibits bistability. By tuning these
bistability parameters, we found a family of biophysically reasonable systems
that are consistent with an experimentally determined bistable region for
induction by thio-methylgalactoside (Ozbudak et al. Nature 427:737, 2004). The
model predicts that bistability can be abolished when passive transport or
permease export becomes sufficiently large; the former case is especially
relevant to induction by isopropyl-beta, D-thiogalactopyranoside. To model
regulation by lactose, we developed similar equations in which allolactose, a
metabolic intermediate in lactose metabolism and a natural inducer of lac, is
the inducer. For biophysically reasonable parameter values, these equations
yield no bistability in response to induction by lactose; however, systems with
an unphysically small permease-dependent export effect can exhibit small
amounts of bistability for limited ranges of parameter values. These results
cast doubt on the relevance of bistability in the lac operon within the natural
context of E. coli, and help shed light on the controversy among existing
theoretical studies that address this issue. The results also suggest an
experimental approach to address the relevance of bistability in the lac operon
within the natural context of E. coli.


  We present a perturbation theory to find the response of an anisotropic DNA
to the external tension. It is shown that the anisotropy has a nonzero but
small contribution to the force-extension curve of the DNA. Thus an anisotropic
DNA behaves like an isotropic one with an effective bending constant equal to
the harmonic average of its soft and hard bending constants.


  In the presence of gratuitous inducers, the lac operon of Escherichia coli
exhibits bistability. Most models in the literature assume that the inducer
enters the cell via the carrier (permease), and exits by a diffusion-like
process. The diffusive influx and carrier efflux are neglected. However,
analysis of the data shows that in non-induced cells, the diffusive influx is
comparable to the carrier influx, and in induced cells, the carrier efflux is 7
times the diffusive efflux. Since bistability entails the coexistence of steady
states corresponding to both non-induced and induced cells, neither one of
these fluxes can be ignored. Here, we formulate a model accounting for both
fluxes. We show that: (a) The thresholds of bistability are profoundly affected
by both fluxes. The diffusive influx reduces the on threshold by enhancing
inducer accumulation in non-induced cells. The carrier efflux increases the off
threshold by decreasing inducer accumulation in induced cells. (b) Simulations
of the model with experimentally measured parameter values are in good
agreement with the data for IPTG. However, there are discrepancies with respect
to the data for TMG. They are most likely due to two features missing from the
model, namely, the variation of the inducer exclusion effect and the specific
growth rate with the lactose enzyme levels. (c) The steady states and
thresholds obtained in the presence of both fluxes are well approximated by
simple analytical expressions, which may be useful for the preliminary design
of the lac genetic switch in synthetic biology.


  Phenotypic heterogeneity in an isogenic, microbial population enables a
subset of the population to persist under stress. In mycobacteria, stresses
like nutrient and oxygen deprivation activate the stress response pathway
involving the two-component system MprAB and the sigma factor, SigE. SigE in
turn activates the expression of the stringent response regulator, rel. The
enzyme polyphosphate kinase 1 (PPK1) regulates this pathway by synthesizing
polyphosphate required for the activation of MprB. The precise manner in which
only a subpopulation of bacterial cells develops persistence, remains unknown.
Rel is required for mycobacterial persistence. Here we show that the
distribution of rel expression levels in a growing population of mycobacteria
is bimodal with two distinct peaks corresponding to low (L) and high (H)
expression states, and further establish that a positive feedback loop
involving the mprAB operon along with stochastic gene expression are
responsible for the phenotypic heterogeneity. Combining single cell analysis by
flow cytometry with theoretical modeling, we observe that during growth,
noise-driven transitions take a subpopulation of cells from the L to the H
state within a "window of opportunity" in time preceding the stationary phase.
We find evidence of hysteresis in the expression of rel in response to changing
concentrations of PPK1. Our results provide, for the first time, evidence that
bistability and stochastic gene expression could be important for the
development of "heterogeneity with an advantage" in mycobacteria.


  Bird flocking is a striking example of collective animal behaviour. A vivid
illustration of this phenomenon is provided by the aerial display of vast
flocks of starlings gathering at dusk over the roost and swirling with
extraordinary spatial coherence. Both the evolutionary justification and the
mechanistic laws of flocking are poorly understood, arguably because of a lack
of data on large flocks. Here, we report a quantitative study of aerial
display. We measured the individual three-dimensional positions in compact
flocks of up to 2700 birds. We investigated the main features of the flock as a
whole - shape, movement, density and structure - and discuss these as emergent
attributes of the grouping phenomenon. We find that flocks are relatively thin,
with variable sizes, but constant proportions. They tend to slide parallel to
the ground and, during turns, their orientation changes with respect to the
direction of motion. Individual birds keep a minimum distance from each other
that is comparable to their wingspan. The density within the aggregations is
non-homogeneous, as birds are packed more tightly at the border compared to the
centre of the flock. These results constitute the first set of large-scale data
on three-dimensional animal aggregations. Current models and theories of
collective animal behaviour can now be tested against these results.


  The most startling examples of collective animal behaviour are provided by
very large and cohesive groups moving in three dimensions. Paradigmatic
examples are bird flocks, fish schools and insect swarms. However, because of
the sheer technical difficulty of obtaining 3D data, empirical studies
conducted to date have only considered loose groups of a few tens of animals.
Moreover, these studies were very seldom conducted in the field. Recently the
STARFLAG project achieved the 3D reconstruction of thousands of birds under
field conditions, thus opening the way to a new generation of quantitative
studies of collective animal behaviour. Here, we review the main technical
problems in 3D data collection of large animal groups and we outline some of
the methodological solutions adopted by the STARFLAG project. In particular, we
explain how to solve the stereoscopic correspondence - or matching - problem,
which was the major bottleneck of all 3D studies in the past.


  The study of collective animal behaviour must progress through a comparison
between the theoretical predictions of numerical models and data coming from
empirical observations. To this aim it is important to develop methods of
three-dimensional (3D) analysis that are at the same time informative about the
structure of the group and suitable to empirical data. In fact, empirical data
are considerably noisier than numerical data, and they are subject to several
constraints. We review here the tools of analysis used by the STARFLAG project
to characterise the 3D structure of large flocks of starlings in the field. We
show how to avoid the most common pitfalls i the quantitative analysis of 3D
animal groups, with particular attention to the problem of the bias introduced
by the border of the group. By means of practical examples, we demonstrate that
neglecting border effects gives rise to artefacts when studying the 3D
structure of a group. Moreover, we show that mathematical rigour is essential
to distinguish important biological properties from trivial geometric features
of animal groups.


  The problem of RNA genomes packaged inside spherical viruses is studied. The
viral capsid is modeled as a hollowed sphere. The attraction between RNA
molecules and the inner viral capsid is assumed to be non-specific and occurs
at the inner capsid surface only. For small capsid attraction, it is found that
monomer concentration of RNA molecules is maximum at the center of the capsid
to maximize their configurational entropy. For stronger capsid attraction, RNA
concentration peaks at some distance near the capsid. In the latter case, the
competition between the branching of RNA secondary struture and its adsorption
to the inner capsid results in the formation of a dense layer of RNA near
capsid surface. The layer thickness is a slowly varying (logarithmic) function
of the capsid inner radius. Consequently, for immediate strength of RNA-capsid
interaction, the amount of RNA packaged inside a virus is proportional to the
capsid {\em area} (or the number of proteins) instead of its volume. The
numerical profiles describe reasonably well the experimentally observed RNA
nucleotide concentration profiles of various viruses.


  Edwards et al. [Nature 449, 1044-1048 (2007)] revisited well-known studies
reporting power-laws in the frequency distribution of flight duration of
wandering albatrosses, and concluded that no L\'evy process could model recent
observations with higher resolution. Here we show that this re-analysis suffers
from a conceptual misunderstanding, and that the new albatross data remain
consistent with a biological L\'evy flight.


  Protein aggregation in the form of amyloid fibrils has important biological
and technological implications. Although the self-assembly process is highly
efficient, aggregates not in the fibrillar form would also occur and it is
important to include these disordered species when discussing the thermodynamic
equilibrium behavior of the system. Here, we initiate such a task by
considering a mixture of monomeric proteins and the corresponding aggregates in
the disordered form (micelles) and in the fibrillar form (amyloid fibrils).
Starting with a model on the respective binding free energies for these
species, we calculate their concentrations at thermal equilibrium. We then
discuss how the incorporation of the disordered structure furthers our
understanding on the various amyloid promoting factors observed empirically,
and on the kinetics of fibrilization.


  The structure and dynamics of a typical biological system are complex due to
strong and inhomogeneous interactions between its constituents. The
investigation of such systems with classical mathematical tools, such as
differential equations for their dynamics, is not always suitable. The graph
theoretical models may serve as a rough but powerful tool in such cases. In
this thesis, I first consider the network modeling for the representation of
the biological systems. Both the topological and dynamical investigation tools
are developed and applied to the various model networks. In particular, the
attractor features' scaling with system size and distributions are explored for
model networks. Moreover, the theoretical robustness expressions are discussed
and computational studies are done for confirmation. The main biological
research in this thesis is to investigate the transcriptional regulation of
gene expression with synchronously and deterministically updated Boolean
network models. I explore the attractor structure and the robustness of the
known interaction network of the yeast, Saccharomyces Cerevisiae and compare
with the model networks. Furthermore, I discuss a recent model claiming a
possible root to the topology of the yeast's gene regulation network and
investigate this model dynamically. The thesis also included another study
which investigates a relation between folding kinetics with a new network
representation, namely, the incompatibility network of a protein's native
structure. I showed that the conventional topological aspects of these networks
are not statistically correlated with the phi-values, for the limited data that
is available.


  Heterogeneity is an important property of any population experiencing a
disease. Here we apply general methods of the theory of heterogeneous
populations to the simplest mathematical models in epidemiology. In particular,
an SIR (susceptible-infective-removed) model is formulated and analyzed for
different sources of heterogeneity. It is shown that a heterogeneous model can
be reduced to a homogeneous model with a nonlinear transmission function, which
is given in explicit form. The widely used power transmission function is
deduced from a heterogeneous model with the initial gamma-distribution of the
disease parameters. Therefore, a mechanistic derivation of the phenomenological
model, which mimics reality very well, is provided. The equation for the final
size of an epidemic for an arbitrary initial distribution is found. The
implications of population heterogeneity are discussed, in particular, it is
pointed out that usual moment-closure methods can lead to erroneous conclusions
if applied for the study of the long-term behavior of the model.


  The equilibrium statistical properties of DNA denaturation bubbles are
examined in detail within the framework of the Peyrard-Bishop-Dauxois model.
Bubble formation in homogeneous DNA is found to depend crucially on the
presence of nonlinear base-stacking interactions. Small bubbles extending over
less than 10 base pairs are associated with much larger free energies of
formation per site than larger bubbles. As the critical temperature is
approached, the free energy associated with further bubble growth becomes
vanishingly small. An analysis of average displacement profiles of bubbles of
varying sizes at different temperatures reveals almost identical scaled shapes
in the absence of nonlinear stacking; nonlinear stacking leads to distinct
scaled shapes of large and small bubbles.


  One of the important challenges facing high resolution molecular dynamics
calculations is to reproduce biological functions that occur in the macroscopic
world and involve macroscopic variables. Some of these functions can be well
described by reduced models, with few parameters, over a wide range of
conditions, even though the functions are complex. We are all eager to see how
molecular dynamics will accurately estimate the thermodynamic variables that
define and control so many macroscopic biological functions.


  A didactic introduction, dated by 1999, to the ideas of the papers
arXiv:q-bio/0701050 and arXiv:0704.0034


  The network paradigm is increasingly used to describe the dynamics of complex
systems. Here we review the current results and propose future development
areas in the assessment of perturbation waves, i.e. propagating structural
changes in amino acid networks building individual protein molecules and in
protein-protein interaction networks (interactomes). We assess the
possibilities and critically review the initial attempts for the application of
game theory to the often rather complicated process, when two protein molecules
approach each other, mutually adjust their conformations via multiple
communication steps and finally, bind to each other. We also summarize
available data on the application of percolation theory for the prediction of
amino acid network- and interactome-dynamics. Furthermore, we give an overview
of the dissection of signals and noise in the cellular context of various
perturbations. Finally, we propose possible applications of the reviewed
methodologies in drug design.


  Consensus methods provide a useful strategy for combining information from a
collection of gene trees. An important application of consensus methods is to
combine gene trees to estimate a species tree. To investigate the theoretical
properties of consensus trees that would be obtained from large numbers of loci
evolving according to a basic evolutionary model, we construct consensus trees
from independent gene trees that occur in proportion to gene tree probabilities
derived from coalescent theory. We consider majority-rule, rooted triple (R*),
and greedy consensus trees constructed from known gene trees, both in the
asymptotic case as numbers of gene trees approach infinity and for finite
numbers of genes. Our results show that for some combinations of species tree
branch lengths, increasing the number of independent loci can make the
majority-rule consensus tree more likely to be at least partially unresolved
and the greedy consensus tree less likely to match the species tree. However,
the probability that the R* consensus tree has the species tree topology
approaches 1 as the number of gene trees approaches infinity. Although the
greedy consensus algorithm can be the quickest to converge on the correct
species tree when increasing the number of gene trees, it can also be
positively misleading. The majority-rule consensus tree is not a misleading
estimator of the species tree topology, and the R* consensus tree is a
statistically consistent estimator of the species tree topology. Our results
therefore suggest a method for using multiple loci to infer the species tree
topology, even when it is discordant with the most likely gene tree.


  We construct and study an extended random matrix model of RNA (polymer)
folding. A perturbation which acts on all the nucleotides in the chain is added
to the action of the RNA partition function. The effect of this perturbation on
the partition function and the Genus Distributions is studied. This
perturbation distinguishes between the paired and unpaired bases. For example,
for $\alpha = 1$ (where $\alpha$ is the ratio of the strengths of the original
and perturbed term in the action) the partition function and genus distribution
for odd lengths vanish completely. This partition function and the genus
distribution is non-zero for even lengths where structures with fully paired
bases only remain. This implies that (i). the genus distributions are different
and (ii). there is a ``structural transition'' (from an ``unpaired-paired base
phase'' to a ``completely paired base phase'') as $\alpha$ approaches 1 in the
extended matrix models. We compare the results of the extended RNA model with
the results of G. Vernizzi, H. Orland and A. Zee in PRL 94, 168103(2005).


  An essential requirement for the representation of functional patterns in
complex neural networks, such as the mammalian cerebral cortex, is the
existence of stable network activations within a limited critical range. In
this range, the activity of neural populations in the network persists between
the extremes of quickly dying out, or activating the whole network. The nerve
fiber network of the mammalian cerebral cortex possesses a modular organization
extending across several levels of organization. Using a basic spreading model
without inhibition, we investigated how functional activations of nodes
propagate through such a hierarchically clustered network. The simulations
demonstrated that persistent and scalable activation could be produced in
clustered networks, but not in random networks of the same size. Moreover, the
parameter range yielding critical activations was substantially larger in
hierarchical cluster networks than in small-world networks of the same size.
These findings indicate that a hierarchical cluster architecture may provide
the structural basis for the stable and diverse functional patterns observed in
cortical networks.


  Many networks exhibit the small-world property of the neighborhood
connectivity being higher than in comparable random networks. However, the
standard measure of local neighborhood clustering is typically not defined if a
node has one or no neighbors. In such cases, local clustering has traditionally
been set to zero and this value influenced the global clustering coefficient.
Such a procedure leads to underestimation of the neighborhood clustering in
sparse networks. We propose to include $\theta$ as the proportion of leafs and
isolated nodes to estimate the contribution of these cases and provide a
formula for estimating a clustering coefficient excluding these cases from the
Watts and Strogatz (1998 Nature 393 440-2) definition of the clustering
coefficient. Excluding leafs and isolated nodes leads to values which are up to
140% higher than the traditional values for the observed networks indicating
that neighborhood connectivity is normally underestimated. We find that the
definition of the clustering coefficient has a major effect when comparing
different networks. For metabolic networks of 43 organisms, relations changed
for 58% of the comparisons when a different definition was applied. We also
show that the definition influences small-world features and that the
classification can change from non-small-world to small-world network. We
discuss the use of an alternative measure, disconnectedness D, which is less
influenced by leafs and isolated nodes.


  Bistability of MAP kinase (MAPK) activity has been suggested to contribute to
several cellular processes, including differentiation and long-term synaptic
potentiation. A recent model (48) predicts bistability due to interactions of
the kinases and phosphatases in the MAPK pathway, without feedback from MAPK to
earlier reactions. Using this model and enzyme concentrations appropriate for
neurons, we simulated bistable MAPK activity, but bistability only was present
within a relatively narrow range of activity of Raf, the first pathway kinase.
Stochastic fluctuations in molecule numbers eliminated bistability for small
molecule numbers, such as are expected in the volume of a dendritic spine.
However, positive feedback loops have been posited from MAPK up to Raf
activation. One proposed loop in which MAPK directly activates Raf was
incorporated into the model. We found that such feedback greatly enhanced the
robustness of both stable states of MAPK activity to stochastic fluctuations
and to parameter variations. Bistability was robust for molecule numbers
plausible for a dendritic spine volume. The upper state of MAPK activity was
resistant to inhibition of MEK activation for > 1 h, suggesting inhibitor
experiments have not sufficed to rule out a role for persistent MAPK activity
in LTP maintenance. These simulations suggest that persistent MAPK activity and
consequent upregulation of translation may contribute to LTP maintenance and to
long-term memory. Experiments using a fluorescent MAPK substrate may further
test this hypothesis.


  We present a top-down approach to the study of the dynamics of icosahedral
virus capsids, in which each protein is approximated by a point mass. Although
this represents a rather crude coarse-graining, we argue that it highlights
several generic features of vibrational spectra which have been overlooked so
far. We furthermore discuss the consequences of approximate inversion symmetry
as well as the role played by Viral Tiling Theory in the study of virus capsid
vibrations.


  Signal duration (e.g. the time scales over which an active signaling
intermediate persists) is a key regulator of biological decisions in myriad
contexts such as cell growth, proliferation, and developmental lineage
commitments. Accompanying differences in signal duration are numerous
downstream biological processes that require multiple steps of biochemical
regulation. Here, we present an analysis that investigates how simple
biochemical motifs that involve multiple stages of regulation can be
constructed to differentially process signals that persist at different time
scales. We compute the dynamic gain within these networks and resulting power
spectra to better understand how biochemical networks can integrate signals at
different time scales. We identify topological features of these networks that
allow for different frequency dependent signal processing properties. Our
studies suggest design principles for why signal duration in connection with
multiple steps of downstream regulation is a ubiquitous control motif in
biochemical systems.


  A scaling analysis is presented of the statistics of long DNA confined in
nanochannels and nanoslits. It is argued that there are several regimes in
between the de Gennes and Odijk limits introduced long ago. The DNA chain folds
back on itself giving rise to a global persistence length which may be very
large owing to entropic deflection. Moreover, there is an orientational
excluded-volume effect between the DNA segments imposed solely by the
nanoconfinement. These two effects cause the chain statistics to be intricate
leading to nontrivial power laws for the chain extension in the intermediate
regimes. It is stressed that DNA confinement within nanochannels differs from
that in nanoslits because the respective orientational excluded-volume effects
are not the same.


  In human societies the probability of strategy adoption from a given person
may be affected by the personal features. Now we investigate how an
artificially imposed restricted ability to reproduce, overruling ones fitness,
affects an evolutionary process. For this purpose we employ the evolutionary
prisoner's dilemma game on different complex graphs. Reproduction restrictions
can have a facilitative effect on the evolution of cooperation that sets in
irrespective of particularities of the interaction network. Indeed, an
appropriate fraction of less fertile individuals may lead to full supremacy of
cooperators where otherwise defection would be widespread. By studying
cooperation levels within the group of individuals having full reproduction
capabilities, we reveal that the recent mechanism for the promotion of
cooperation is conceptually similar to the one reported previously for
scale-free networks. Our results suggest that the diversity in the reproduction
capability, related to inherently different attitudes of individuals, can
enforce the emergence of cooperative behavior among selfish competitors.


  A generalized autocatalytic model for chiral polymerization is investigated
in detail. Apart from enantiomeric cross-inhibition, the model allows for the
autogenic (non-catalytic) formation of left and right-handed monomers from a
substrate with reaction rates $\epsilon_L$ and $\epsilon_R$, respectively. The
spatiotemporal evolution of the net chiral asymmetry is studied for models with
several values of the maximum polymer length, N. For N=2, we study the validity
of the adiabatic approximation often cited in the literature. We show that the
approximation obtains the correct equilibrium values of the net chirality, but
fails to reproduce the short time behavior. We show also that the autogenic
term in the full N=2 model behaves as a control parameter in a chiral symmetry-
breaking phase transition leading to full homochirality from racemic initial
conditions. We study the dynamics of the N -> infinity model with symmetric
($\epsilon_L = \epsilon_R$) autogenic formation, showing that it only achieves
homochirality for $\epsilon < \epsilon_c$, where $\epsilon_c$ is an N-dependent
critical value. For $\epsilon \leq \epsilon_c$ we investigate the behavior of
models with several values of N, showing that the net chiral asymmetry grows as
tanh(N). We show that for a given symmetric autogenic reaction rate, the net
chirality and the concentrations of chirally pure polymers increase with the
maximum polymer length in the model. We briefly discuss the consequences of our
results for the development of homochirality in prebiotic Earth and possible
experimental verification of our findings.


  Motivation: Algorithms for differential analysis of microarray data are vital
to modern biomedical research. Their accuracy strongly depends on effective
treatment of inter-gene correlation. Correlation is ordinarily accounted for in
terms of its effect on significance cut-offs. In this paper it is shown that
correlation can, in fact, be exploited {to share information across tests},
which, in turn, can increase statistical power.
  Results: Vastly and demonstrably improved differential analysis approaches
are the result of combining identifiability (the fact that in most microarray
data sets, a large proportion of genes can be identified a priori as
non-differential) with optimization criteria that incorporate correlation. As a
special case, we develop a method which builds upon the widely used two-sample
t-statistic based approach and uses the Mahalanobis distance as an optimality
criterion. Results on the prostate cancer data of Singh et al. (2002) suggest
that the proposed method outperforms all published approaches in terms of
statistical power.
  Availability: The proposed algorithm is implemented in MATLAB and in R. The
software, called Tellipsoid, and relevant data sets are available at
http://www.egr.msu.edu/~desaikey


  De novo prediction of protein structures, the prediction of structures from
amino-acid sequences which are not similar to those of hitherto resolved
structures, has been one of the major challenges in molecular biophysics. In
this paper, we develop a new method of de novo prediction, which combines the
fragment assembly method and the simulation of physical folding process:
Structures which have consistently assembled fragments are dynamically searched
by Langevin molecular dynamics of conformational change. The benchmarking test
shows that the prediction is improved when the candidate structures are
cross-checked by an empirically derived score function.


  A central model in theoretical ecology considers the competition of a range
of species for a broad spectrum of resources. Recent studies have shown that
essentially two different outcomes are possible. Either the species surviving
competition are more or less uniformly distributed over the resource spectrum,
or their distribution is 'lumped' (or 'clumped'), consisting of clusters of
species with similar resource use that are separated by gaps in resource space.
Which of these outcomes will occur crucially depends on the competition kernel,
which reflects the shape of the resource utilization pattern of the competing
species. Most models considered in the literature assume a Gaussian competition
kernel. This is unfortunate, since predictions based on such a Gaussian
assumption are not robust. In fact, Gaussian kernels are a border case
scenario, and slight deviations from this function can lead to either uniform
or lumped species distributions. Here we illustrate the non-robustness of the
Gaussian assumption by simulating different implementations of the standard
competition model with constant carrying capacity. In this scenario, lumped
species distributions can come about by secondary ecological or evolutionary
mechanisms or by details of the numerical implementation of the model. We
analyze the origin of this sensitivity and discuss it in the context of recent
applications of the model.


  A phenomenological model based on the three-dimensional theory of nonlinear
elasticity is developed to describe the phenomenon of overstretching in the
force-extension curve for dsDNA. By using the concept of a material with
multiple reference configurations a single formula is obtained to fit the
force-extension curve.


  Reaction-diffusion processes can be adopted to model a large number of
dynamics on complex networks, such as transport processes or epidemic
outbreaks. In most cases, however, they have been studied from a fermionic
perspective, in which each vertex can be occupied by at most one particle.
While still useful, this approach suffers from some drawbacks, the most
important probably being the difficulty to implement reactions involving more
than two particles simultaneously. Here we introduce a general framework for
the study of bosonic reaction-diffusion processes on complex networks, in which
there is no restriction on the number of interacting particles that a vertex
can host. We describe these processes theoretically by means of continuous time
heterogeneous mean-field theory and divide them into two main classes: steady
state and monotonously decaying processes. We analyze specific examples of both
behaviors within the class of one-species process, comparing the results
(whenever possible) with the corresponding fermionic counterparts. We find that
the time evolution and critical properties of the particle density are
independent of the fermionic or bosonic nature of the process, while
differences exist in the functional form of the density of occupied vertices in
a given degree class k. We implement a continuous time Monte Carlo algorithm,
well suited for general bosonic simulations, which allow us to confirm the
analytical predictions formulated within mean-field theory. Our results, both
at the theoretical and numerical level, can be easily generalized to tackle
more complex, multi-species, reaction-diffusion processes, and open a promising
path for a general study and classification of this kind of dynamical systems
on complex networks.


  We study the coupled dynamics of primary and secondary structure formation
(i.e. slow genetic sequence selection and fast folding) in the context of a
solvable microscopic model that includes both short-range steric forces and and
long-range polarity-driven forces. Our solution is based on the diagonalization
of replicated transfer matrices, and leads in the thermodynamic limit to
explicit predictions regarding phase transitions and phase diagrams at genetic
equilibrium. The predicted phenomenology allows for natural physical
interpretations, and finds satisfactory support in numerical simulations.


  Motivation: Recent studies of genomic-scale regulatory networks suggested
that a feed-forward loop (FFL) circuitry is a key component of many such
networks. This led to a study of the functional properties of different FFL
types, where the regulatory elements are transcription factors.
  Results: Here we investigate these properties when the mediating regulatory
element of the loop is a microRNA. We find that many of the FFL properties are
enhanced within this setup. We then proceed to identify all such FFLs in the D.
Melanogaster regulatory network. We observe that in FFLs rooted at the same
transcription factor there are significant correlations between the number of
predicted binding sites for the transcription factor and for the microRNA.
  Conclusions: Based on a modeling approach we suggest that these correlations
may be an outcome of the type of FFL preferred by the transcription factor.
This may help elucidate the type of regulation the TF confers.


  We consider a model of a population of fixed size N in which each individual
gets replaced at rate one and each individual experiences a mutation at rate
\mu. We calculate the asymptotic distribution of the time that it takes before
there is an individual in the population with m mutations. Several different
behaviors are possible, depending on how \mu changes with N. These results have
applications to the problem of determining the waiting time for regulatory
sequences to appear and to models of cancer development.


  Hydrodynamic forcing plays an important role in shaping the dynamics of
marine organisms, in particular of plankton. In this work we study the
planktonic biological activity in the wake of an island which is close to an
upwelling region. Our research is based on numerical analysis of a kinematic
flow mimicking the hydrodynamics in the wake, coupled to a three-component
plankton model. Depending on model parameters different phenomena are
described: a) The lack of transport of nutrients and plankton across the wake,
so that the influence of upwelling on primary production on the other side of
the wake is blocked. b) For sufficiently high vorticity, the role of the wake
in facilitating this transport and leading to an enhancement of primary
production. Finally c) we show that under certain conditions the interplay
between wake structures and biological growth leads to plankton blooms inside
mesoscale hydrodynamic vortices that act as incubators of primary production.


  Recently, using a numerical surface cooling approach, we have shown that
highly energetic discrete breathers (DB) can form in the stiffest parts of
nonlinear network models of large protein structures. In the present study,
using an analytical approach, we extend our previous results to low-energy
discrete breathers as well as to smaller proteins
  We confirm and further scrutinize the striking site selectiveness of energy
localisation in the presence of spatial disorder. In particular, we find that,
as a sheer consequence of disorder, a non-zero energy gap for exciting a DB at
a given site either exists or not. Remarkably, in the former case, the gaps
arise as result of the impossibility of exciting small-amplitude modes in the
first place. On the contrary, in the latter case, a small subset of linear edge
modes act as accumulation points, whereby DBs can be continued to arbitrary
small energies, while unavoidably approaching one of such normal modes. In
particular, the case of the edge mode seems peculiar, its dispersion relation
being simple and little system-dependent.
  Concerning the structure-dynamics relationship, we find that the regions of
protein structures where DBs form easily (zero or small gaps) are unfailingly
the most highly connected ones, also characterized by weak local clustering.
Remarkably, a systematic analysis on a large database of enzyme structures
reveals that amino-acid residues involved in enzymatic activity tend to be
located in such regions. This finding reinforces the idea that localised modes
of nonlinear origin may play an important biological role, e.g. by providing a
ready channel for energy storage and/or contributing to lower energy barriers
of chemical reactions.


  We study the dynamics of epidemic and reaction-diffusion processes in
metapopulation models with heterogeneous connectivity pattern. In SIR-like
processes, along with the standard local epidemic threshold, the system
exhibits a global invasion threshold. We provide an explicit expression of the
threshold that sets a critical value of the diffusion/mobility rate below which
the epidemic is not able to spread to a macroscopic fraction of subpopulations.
The invasion threshold is found to be affected by the topological fluctuations
of the metapopulation network. The presented results provide a general
framework for the understanding of the effect of travel restrictions in
epidemic containment.


  We propose a class of microstructurally informed models for the linear
elastic mechanical behavior of cross-linked polymer networks such as the actin
cytoskeleton. Salient features of the models include the possibility to
represent anisotropic mechanical behavior resulting from anisotropic filament
distributions, and a power-law scaling of the mechanical properties with the
filament density. Mechanical models within the class are parameterized by seven
different constants. We demonstrate a procedure for determining these constants
using finite element models of three-dimensional actin networks. Actin
filaments and cross-links were modeled as elastic rods, and the networks were
constructed at physiological volume fractions and at the scale of an image
voxel. We show the performance of the model in estimating the mechanical
behavior of the networks over a wide range of filament densities and degrees of
anisotropy.


  Background: The evolution of microRNA regulation in metazoans is a mysterious
process: MicroRNA sequences are highly conserved among distal organisms, but on
the other hand, there is no evident conservation of their targets. Results: We
study this extensive rewiring of the microRNA regulatory network by analyzing
the evolutionary trajectories of duplicated genes in D. melanogatser. We find
that in general microRNA-targeted genes tend to avoid gene duplication.
However, in cases where gene duplication is evident, we find that the gene that
displays high divergence from the ancestral gene at the sequence level is also
likely to be associated in an opposing manner with the microRNA regulatory
system - if the ancestral gene is a miRNA target then the divergent gene tends
not to be, and vice versa. Conclusions: This suggests that miRNAs not only have
a role in conferring expression robustness, as was suggested by previous works,
but are also an accessible tool in evolving expression divergence.


  Cancer invasion and metastasis depend on angiogenesis. The cellular processes
(growth, migration, and apoptosis) that occur during angiogenesis are tightly
regulated by signaling molecules. Thus, understanding how cells synthesize
multiple biochemical signals initiated by key external stimuli can lead to the
development of novel therapeutic strategies to combat cancer. In the face of
large amounts of disjoint experimental data generated from multitudes of
laboratories using various assays, theoretical signal transduction models
provide a framework to distill this vast amount of data. Such models offer an
opportunity to formulate and test new hypotheses, and can be used to make
experimentally verifiable predictions. This study is the first to propose a
network model that highlights the cross-talk between the key receptors involved
in angiogenesis, namely growth factor, integrin, and cadherin receptors. From
available experimental data, we construct a stochastic Boolean network model of
receptor cross-talk, and systematically analyze the dynamical stability of the
network under continuous-time Boolean dynamics with a noisy production
function. We find that the signal transduction network exhibits a robust and
fast response to external signals, independent of the internal cell state. We
derive an input-output table that maps external stimuli to cell phenotypes,
which is extraordinarily stable against molecular noise with one important
exception: an oscillatory feedback loop between the key signaling molecules
RhoA and Rac1 is unstable under arbitrarily low noise, leading to erratic,
dysfunctional cell motion. Finally, we show that the network exhibits an
apoptotic response rate that increases with noise, suggesting that the
probability of programmed cell death depends on cell health.


  We study the adaptation dynamics of an initially maladapted population
evolving via the elementary processes of mutation and selection. The evolution
occurs on rugged fitness landscapes which are defined on the multi-dimensional
genotypic space and have many local peaks separated by low fitness valleys. We
mainly focus on the Eigen's model that describes the deterministic dynamics of
an infinite number of self-replicating molecules. In the stationary state, for
small mutation rates such a population forms a {\it quasispecies} which
consists of the fittest genotype and its closely related mutants. The
quasispecies dynamics on rugged fitness landscape follow a punctuated (or
step-like) pattern in which a population jumps from a low fitness peak to a
higher one, stays there for a considerable time before shifting the peak again
and eventually reaches the global maximum of the fitness landscape. We
calculate exactly several properties of this dynamical process within a
simplified version of the quasispecies model.


  We study the interplay of hydrodynamic mesoscale structures and the growth of
plankton in the wake of an island, and its interaction with a coastal
upwelling. Our focus is on a mechanism for the emergence of localized plankton
blooms in vortices. Using a coupled system of a kinematic flow mimicking the
mesoscale structures behind the island and a simple three component model for
the marine ecosystem, we show that the long residence times of nutrients and
plankton in the vicinity of the island and the confinement of plankton within
vortices are key factors for the appearance of localized plankton blooms


  The paper analyses stochastic systems describing reacting molecular systems
with a combination of two types of state spaces, a finite-dimensional, and an
infinite dimenional part. As a typical situation consider the interaction of
larger macro-molecules, finite and small in numbers per cell (like protein
complexes), with smaller, very abundant molecules, for example metabolites. We
study the construction of the continuum approximation of the associated Master
Equation (ME) by using the Trotter approximation [27]. The continuum limit
shows regimes where the finite degrees of freedom evolve faster than the
infinite ones. Then we develop a rigourous asymptotic adiabatic theory upon the
condition that the jump process arising from the finite degrees of freedom of
the Markov Chain (MC, typically describing conformational changes of the
macro-molecules) occurs with large frequency. In a second part of this work,
the theory is applied to derive typical enzyme kinetics in an alternative way
and interpretation within this framework.


  In this paper we consider deterministic limits of molecular stochastic
systems with finite and infinite degrees of freedom. The method to obtain the
deterministic vector field is based on the continuum limit of such microscopic
systems which has been derived in [11]. With the aid of the theory we finally
develop a new approach for molecular systems that describe typical enzyme
kinetics or other interactions between molecular machines like genetic elements
and smaller communicating molecules. In contrast to the literature on enzyme
kinetics the resulting deterministic functional responses are not derived by
time-scale arguments on the macroscopic level, but are a result of time scaling
transition rates on the discrete microscopic level. We present several examples
of common functional responses found in the literature, like Michaelis-Menten
and Hill equation. We finally give examples of more complex but typical
macro-molecular machinery.


  In most natural sciences there is currently the insight that it is necessary
to bridge gaps between different processes which can be observed on different
scales. This is especially true in the field of chemical reactions where the
abilities to form bonds between different types of atoms and molecules create
much of the properties we experience in our everyday life, especially in all
biological activity. There are essentially two types of processes related to
biochemical reaction networks, the interactions among molecules and
interactions involving their conformational changes, so in a sense, their
internal state. The first type of processes can be conveniently approximated by
the so-called mass-action kinetics, but this is not necessarily so for the
second kind where molecular states do not define any kind of density or
concentration. In this paper we demonstrate the necessity to study reaction
networks in a stochastic formulation for which we can construct a coherent
approximation in terms of specific space-time scales and the number of
particles. The continuum limit procedure naturally creates equations of
Fokker-Planck type where the evolution of the concentration occurs on a slower
time scale when compared to the evolution of the conformational changes, for
example triggered by binding or unbinding events with other (typically smaller)
molecules. We apply the asymptotic theory to derive the effective, i.e.
macroscopic dynamics of the biochemical reaction system. The theory can also be
applied to other processes where entities can be described by finitely many
internal states, with changes of states occuring by arrival of other entities
described by a birth-death process.


  A simple flashing ratchet model in two dimensions is proposed to simulate the
hand-over-hand motion of two head molecular motors like kinesin. Extensive
Langevin simulations of the model are performed. Good qualitative agreement
with the expected behavior is observed. We discuss different regimes of motion
and efficiency depending of model parameters.


  A recent problem in community ecology lies in defining structures behind
matrices of species interactions. The interest in this area is to quantify the
nestedness degree of the matrix after its maximal packing. In this work we
evaluate nestedness using the sum of all distances of the occupied sites to the
vertex of the matrix. We calculate the distance for two artificial matrices
with the same size and occupancy: a random matrix and a perfect nested one.
Using these two benchmarks we develop a nestedness estimator. The estimator is
applied to a set of 23 real networks of insect-plant interactions.


  We show that the accuracy of reconstrucing an ancestral state is not an
increasing function of the size of taxon sampling.


  Filopodia are long, finger-like membrane tubes supported by cytoskeletal
filaments. Their shape is determined by the stiffness of the actin filament
bundles found inside them and by the interplay between the surface tension and
bending rigidity of the membrane. Although one might expect the Euler buckling
instability to limit the length of filopodia, we show through simple energetic
considerations that this is in general not the case. By further analyzing the
statics of filaments inside membrane tubes, and through computer simulations
that capture membrane and filament fluctuations, we show under which conditions
filopodia of arbitrary lengths are stable. We discuss several in vitro
experiments where this kind of stability has already been observed.
Furthermore, we predict that the filaments in long, stable filopodia adopt a
helical shape.


  A common feature of mycorrhizal observation is the growth of the infection on
the plant root as a percent of the infected root or root tip length. Often,
this is measured as a logistic curve with an eventual, though usually
transient, plateau. It is shown in this paper that the periods of stable
percent infection in the mycorrhizal growth cycle correspond to periods where
both the plant and mycorrhiza growth rates and likely metabolism are tightly
coupled.


  It is suggested that the degree distribution for networks of the
cell-metabolism for simple organisms reflects an ubiquitous randomness. This
implies that natural selection has exerted no or very little pressure on the
network degree distribution during evolution. The corresponding random network,
here termed the blind watchmaker network has a power-law degree distribution
with an exponent gamma >= 2. It is random with respect to a complete set of
network states characterized by a description of which links are attached to a
node as well as a time-ordering of these links. No a priory assumption of any
growth mechanism or evolution process is made. It is found that the degree
distribution of the blind watchmaker network agrees very precisely with that of
the metabolic networks. This implies that the evolutionary pathway of the
cell-metabolism, when projected onto a metabolic network representation, has
remained statistically random with respect to a complete set of network states.
This suggests that even a biological system, which due to natural selection has
developed an enormous specificity like the cellular metabolism, nevertheless
can, at the same time, display well defined characteristics emanating from the
ubiquitous inherent random element of Darwinian evolution. The fact that also
completely random networks may have scale-free node distributions gives a new
perspective on the origin of scale-free networks in general.


  A major obstacle in analyzing the evolution of information exchange and
processing is our insufficient understanding of the underlying signaling and
decision-making biological mechanisms. For instance, it is unclear why are
humans unique in developing such extensive communication abilities. To treat
this problem, a method based on the mutual information approach is developed
that evaluates the information content of communication between interacting
individuals through correlations of their behavior patterns (rather than
calculating the information load of exchanged discrete signals, e.g. Shannon
entropy). It predicts that correlated interactions of the indirect reciprocity
type together with affective behavior and selection rules changing with time
are necessary conditions for the emergence of significant information exchange.
Population size variations accelerate this development. These results are
supported by evidence of demographic bottlenecks, distinguishing human from
other species' (e.g. apes) evolution line. They indicate as well new pathways
for evolution of information based phenomena, such as intelligence and
complexity.


  We investigate how a catastrophic event (modeled as a temporary fall of the
reproduction rate) increases the extinction probability of an isolated
self-regulated stochastic population. Using a variant of the Verhulst logistic
model as an example, we combine the probability generating function technique
with an eikonal approximation to evaluate the exponentially large increase in
the extinction probability caused by the catastrophe. This quantity is given by
the eikonal action computed over "the optimal path" (instanton) of an effective
classical Hamiltonian system with a time-dependent Hamiltonian. For a general
catastrophe the eikonal equations can be solved numerically. For simple models
of catastrophic events analytic solutions can be obtained. One such solution
becomes quite simple close to the bifurcation point of the Verhulst model. The
eikonal results for the increase in the extinction probability caused by a
catastrophe agree well with numerical solutions of the master equation.


  Translocation of a single stranded DNA through genetically engineered
$\alpha$-hemolysin channels with positively charged walls is studied. It is
predicted that transport properties of such channels are dramatically different
from neutral wild type $\alpha$-hemolysin channel. We assume that the wall
charges compensate the fraction $x$ of the bare charge $q_{b}$ of the DNA piece
residing in the channel. Our prediction are as follows (i) At small
concentration of salt the blocked ion current decreases with $x$. (ii) The
effective charge $q$ of DNA piece, which is very small at $x = 0$ (neutral
channel) grows with $x$ and at $x=1$ reaches $q_{b}$. (iii) The rate of DNA
capture by the channel exponentially grows with $x$. Our theory is also
applicable to translocation of a double stranded DNA in narrow solid state
nanopores with positively charged walls.