  We consider a simple working hypothesis that all permeation properties of
open ionic channels can be predicted by understanding electrodiffusion in fixed
structures, without invoking conformation changes, or changes in chemical
bonds. We know, of course, that ions can bind to specific protein structures,
and that this binding is not easily described by the traditional electrostatic
equations of physics textbooks, that describe average electric fields, the
so-called `mean field'. The question is which specific properties can be
explained just by mean field electrostatics and which cannot. I believe the
best way to uncover the specific chemical properties of channels is to invoke
them as little as possible, seeking to explain with mean field electrostatics
first. Then, when phenomena appear that cannot be described that way, by the
mean field alone, we turn to chemically specific explanations, seeking the
appropriate tools (of electrochemistry, Langevin, or molecular dynamics, for
example) to understand them. In this spirit, we turn now to the structure of
open ionic channels, apply the laws of electrodiffusion to them, and see how
many of their properties we can predict just that way.


  We realize a biochemical filtering process by introducing a buffer in a
biocatalytic signal-transduction logic system based on the function of an
enzyme, esterase. The input, ethyl butyrate, is converted into butyric acid-the
output signal, which in turn is measured by the drop in the pH value. The
developed approach offers a versatile "network element" for increasing the
complexity of biochemical information processing systems. Evaluation of an
optimal regime for quality filtering is accomplished in the framework of a
kinetic rate-equation model.


  A comparative classification scheme provides a good basis for several
approaches to understand proteins, including prediction of relations between
their structure and biological function. But it remains a challenge to combine
a classification scheme that describes a protein starting from its well
organized secondary structures and often involves direct human involvement,
with an atomary level Physics based approach where a protein is fundamentally
nothing more than an ensemble of mutually interacting carbon, hydrogen, oxygen
and nitrogen atoms. In order to bridge these two complementary approaches to
proteins, conceptually novel tools need to be introduced. Here we explain how
the geometrical shape of entire folded proteins can be described analytically
in terms of a single explicit elementary function that is familiar from
nonlinear physical systems where it is known as the kink-soliton. Our approach
enables the conversion of hierarchical structural information into a
quantitative form that allows for a folded protein to be characterized in terms
of a small number of global parameters that are in principle computable from
atomary level considerations. As an example we describe in detail how the
native fold of the myoglobin 1M6C emerges from a combination of kink-solitons
with a very high atomary level accuracy. We also verify that our approach
describes longer loops and loops connecting $\alpha$-helices with
$\beta$-strands, with same overall accuracy.


  The extinction transition in the presence of a localized quenched defect is
studied numerically. When the bulk is at criticality, the correlation length
diverges and even an infinite system cannot "decouple" from the defect. The
results presented here suggest that, in 1+1 dimensions, the critical exponent
$\delta$ that controls the asymptotic power-law decay depends on the strength
of the local perturbation. On the other hand, the exponent was found to be
independent of the local arrangement of the defect. In higher dimensions the
defect seems to induce a transient behavior that decays algebraically in time.


  The last few million years on planet Earth have witnessed two remarkable
phases of hominid development, starting with a phase of biological evolution
characterised by rather rapid increase of the size of the brain. This has been
followed by a phase of even more rapid technological evolution and concomitant
expansion of the size of the population, that began when our own particular
`sapiens' species emerged, just a few hundred thousand years ago. The present
investigation exploits the analogy between the neo-Darwinian genetic evolution
mechanism governing the first phase, and the memetic evolution mechanism
governing the second phase. From the outset of the latter until very recently
-- about the year 2000 -- the growth of the global population N was roughly
governed by an equation of the form dN/Ndt= N/T*, in which T* is a coefficient
introduced (in 1960) by von Foerster, who evaluated it empirically as about two
hundred thousand million years. It is shown here how the value of this hitherto
mysterious timescale governing the memetic phase is explicable in terms of what
happenned in the preceding genetic phase. The outcome is that the order of
magnitude of the Foerster timescale can be accounted for as the product of the
relevant (human) generation timescale, about 20 years, with the number of bits
of information in the genome, of the order of ten thousand million. Whereas the
origin of our `homo' genus may well have involved an evolutionary hard step, it
transpires that the emergence of our particular `sapiens' species was rather an
automatic process.


  We determine the Renyi entropies K_q of symbol sequences generated by human
chromosomes. These exhibit nontrivial behaviour as a function of the scanning
parameter q. In the thermodynamic formalism, there are phase transition-like
phenomena close to the q=1 region. We develop a theoretical model for this
based on the superposition of two multifractal sets, which can be associated
with the different statistical properties of coding and non-coding DNA
sequences. This model is in good agreement with the human chromosome data.


  The influence of the salts KCl, NaCl, and NaI at molar concentrations on the
{\alpha}-helical folding kinetics of the alanine-based oligopeptide
Ace-AEAAAKEAAAKA-Nme is investigated by means of (explicit-water) molecular
dynamics simulations and a diffusional analysis. The mean first passage times
for folding and unfolding are found to be highly salt-specific. In particular,
the folding times increase about one order of magnitude for the sodium salts.
The drastic slowing down can be traced back to long-lived, compact
configurations of the partially folded peptide, in which sodium ions are
tightly bound by several carbonyl and carboxylate groups. This multiple
trapping is found to lead to a non-exponential residence time distribution of
the cations in the first solvation shell of the peptide. The analysis of
{\alpha}-helical folding in the framework of diffusion in a reduced
(one-dimensional) free energy landscape further shows that the salt not only
specifically modifies equilibrium properties, but also induces kinetic barriers
due to individual ion binding. In the sodium salts, for instance, the peptide's
configurational mobility (or "diffusivity") can decrease about one order of
magnitude. This study demonstrates the highly specific action of ions and
highlights the intimate coupling of intramolecular friction and solvent effects
in protein folding.


  We study a complementarity game with multiple populations whose members'
offered contributions are put together towards some common aim. When the sum of
the players' offers reaches or exceeds some threshold K, they each receive K
minus their own offers. Else, they all receive nothing. Each player tries to
offer as little as possible, hoping that the sum of the contributions still
reaches K, however. The game is symmetric at the individual level, but has many
equilibria that are more or less favorable to the members of certain
populations. In particular, it is possible that the members of one or several
populations do not contribute anything, a behavior called defecting, while the
others still contribute enough to reach the threshold. Which of these
equilibria then is attained is decided by the dynamics at the population level
that in turn depends on the strategic options the players possess. We find that
defecting occurs when more than 3 populations participate in the game, even
when the strategy scheme employed is very simple, if certain conditions for the
system parameters are satisfied. The results are obtained through systematic
simulations.


  We study a complementarity game as a systematic tool for the investigation of
the interplay between individual optimization and population effects and for
the comparison of different strategy and learning schemes. The game randomly
pairs players from opposite populations. The game is symmetric at the
individual level, but has many equilibria that are more or less favorable to
the members of the two populations. Which of these equilibria then is attained
is decided by the dynamics at the population level. Players play repeatedly,
but in each round with a new opponent. They can learn from their previous
encounters and translate this into their actions in the present round on the
basis of strategic schemes. The schemes can be quite simple, or very elaborate.
We can then break the symmetry in the game and give the members of the two
populations access to different strategy spaces. Typically, simpler strategy
types have an advantage because they tend to go more quickly towards a
favorable equilibrium which, once reached, the other population is forced to
accept. Also, populations with bolder individuals that may not fare so well at
the level of individual performance may obtain an advantage towards ones with
more timid players. By checking the effects of parameters such as the
generation length or the mutation rate, we are able to compare the relative
contributions of individual learning and evolutionary adaptations.


  We present a technique to characterize differentially expressed genes in
terms of their position in a high-dimensional co-expression network. The set-up
of Gaussian graphical models is used to construct representations of the
co-expression network in such a way that redundancy and the propagation of
spurious information along the network are avoided. The proposed inference
procedure is based on the minimization of the Bayesian Information Criterion
(BIC) in the class of decomposable graphical models. This class of models can
be used to represent complex relationships and has suitable properties that
allow to make effective inference in problems with high degree of complexity
(e.g. several thousands of genes) and small number of observations (e.g.
10-100) as typically occurs in high throughput gene expression studies. Taking
advantage of the internal structure of decomposable graphical models, we
construct a compact representation of the co-expression network that allows to
identify the regions with high concentration of differentially expressed genes.
It is argued that differentially expressed genes located in highly
interconnected regions of the co-expression network are less informative than
differentially expressed genes located in less interconnected regions. Based on
that idea, a measure of uncertainty that resembles the notion of relative
entropy is proposed. Our methods are illustrated with three publically
available data sets on microarray experiments (the larger involving more than
50,000 genes and 64 patients) and a short simulation study.


  Boolean threshold networks have recently been proposed as useful tools to
model the dynamics of genetic regulatory networks, and have been successfully
applied to describe the cell cycles of \textit{S. cerevisiae} and \textit{S.
pombe}. Threshold networks assume that gene regulation processes are additive.
This, however, contrasts with the mechanism proposed by S. Kauffman in which
each of the logic functions must be carefully constructed to accurately take
into account the combinatorial nature of gene regulation. While Kauffman
Boolean networks have been extensively studied and proved to have the necessary
properties required for modeling the fundamental characteristics of genetic
regulatory networks, not much is known about the essential properties of
threshold networks. Here we study the dynamical properties of these networks
with different connectivities, activator-repressor proportions,
activator-repressor strengths and different thresholds. Special attention is
paid to the way in which the threshold value affects the dynamical regime in
which the network operates and the structure of the attractor landscape. We
find that only for a very restricted set of parameters, these networks show
dynamical properties consistent with what is observed in biological systems.
The virtues of these properties and the possible problems related with the
restrictions are discussed and related to earlier work that uses these kind of
models.


  If we can lower the number of people needed to vaccinate for a community to
be immune against contagious diseases, we can save resources and life. A key to
reach such a lower threshold of immunization is to find and vaccinate people
who, through their behavior, are more likely to become infected and effective
to spread the disease than the average. Fortunately, the very behavior that
makes these people important to vaccinate can help us finding them. People you
have met recently are more likely to be socially active and thus central in the
contact pattern, and important to vaccinate. We propose two immunization
schemes exploiting temporal contact patterns. Both of these rely only on
obtainable, local information and could implemented in practice. We show that
these schemes outperform benchmark protocols in four real data sets under
various epidemic scenarios. The data sets are dynamic, which enables us to make
more realistic evaluations than other studies - we use information only about
the past to perform the vaccination and the future to simulate disease
outbreaks. We also use models to elucidate the mechanisms behind how the
temporal structures make our immunization protocols efficient.


  Food webs, networks of feeding relationships among organisms, provide
fundamental insights into mechanisms that determine ecosystem stability and
persistence. Despite long-standing interest in the compartmental structure of
food webs, past network analyses of food webs have been constrained by a
standard definition of compartments, or modules, that requires many links
within compartments and few links between them. Empirical analyses have been
further limited by low-resolution data for primary producers. In this paper, we
present a Bayesian computational method for identifying group structure in food
webs using a flexible definition of a group that can describe both functional
roles and standard compartments. The Serengeti ecosystem provides an
opportunity to examine structure in a newly compiled food web that includes
species-level resolution among plants, allowing us to address whether groups in
the food web correspond to tightly-connected compartments or functional groups,
and whether network structure reflects spatial or trophic organization, or a
combination of the two. We have compiled the major mammalian and plant
components of the Serengeti food web from published literature, and we infer
its group structure using our method. We find that network structure
corresponds to spatially distinct plant groups coupled at higher trophic levels
by groups of herbivores, which are in turn coupled by carnivore groups. Thus
the group structure of the Serengeti web represents a mixture of trophic guild
structure and spatial patterns, in contrast to the standard compartments
typically identified in ecological networks. From data consisting only of nodes
and links, the group structure that emerges supports recent ideas on spatial
coupling and energy channels in ecosystems that have been proposed as important
for persistence.


  Fat tailed statistics and power-laws are ubiquitous in many complex systems.
Usually the appearance of of a few anomalously successful individuals
(bio-species, investors, websites) is interpreted as reflecting some inherent
"quality" (fitness, talent, giftedness) as in Darwin's theory of natural
selection. Here we adopt the opposite, "neutral", outlook, suggesting that the
main factor explaining success is merely luck. The statistics emerging from the
neutral birth-death-mutation (BDM) process is shown to fit marvelously many
empirical distributions. While previous neutral theories have focused on the
power-law tail, our theory economically and accurately explains the entire
distribution. We thus suggest the BDM distribution as a standard neutral model:
effects of fitness and selection are to be identified by substantial deviations
from it.


  Evolution is the fundamental physical process that gives rise to biological
phenomena. Yet it is widely treated as a subset of population genetics, and
thus its scope is artificially limited. As a result, the key issues of how
rapidly evolution occurs, and its coupling to ecology have not been
satisfactorily addressed and formulated. The lack of widespread appreciation
for, and understanding of, the evolutionary process has arguably retarded the
development of biology as a science, with disastrous consequences for its
applications to medicine, ecology and the global environment. This review
focuses on evolution as a problem in non-equilibrium statistical mechanics,
where the key dynamical modes are collective, as evidenced by the plethora of
mobile genetic elements whose role in shaping evolution has been revealed by
modern genomic surveys. We discuss how condensed matter physics concepts might
provide a useful perspective in evolutionary biology, the conceptual failings
of the modern evolutionary synthesis, the open-ended growth of complexity, and
the quintessentially self-referential nature of evolutionary dynamics.


  Phylogenetic mixture models are statistical models of character evolution
allowing for heterogeneity. Each of the classes in some unknown partition of
the characters may evolve by different processes, or even along different
trees. The fundamental question of whether parameters of such a model are
identifiable is difficult to address, due to the complexity of the
parameterization. We analyze mixture models on large trees, with many mixture
components, showing that both numerical and tree parameters are indeed
identifiable in these models when all trees are the same. We also explore the
extent to which our algebraic techniques can be employed to extend the result
to mixtures on different trees.


  In this paper we consider a model for the spread of a stochastic SIR
(Susceptible $\to$ Infectious $\to$ Recovered) epidemic on a network of
individuals described by a random intersection graph. Individuals belong to a
random number of cliques, each of random size, and infection can be transmitted
between two individuals if and only if there is a clique they both belong to.
Both the clique sizes and the number of cliques an individual belongs to follow
mixed Poisson distributions. An infinite-type branching process approximation
(with type being given by the length of an individual's infectious period) for
the early stages of an epidemic is developed and made fully rigorous by proving
an associated limit theorem as the population size tends to infinity. This
leads to a threshold parameter $R_*$, so that in a large population an epidemic
with few initial infectives can give rise to a large outbreak if and only if
$R_*>1$. A functional equation for the survival probability of the
approximating infinite-type branching process is determined; if $R_*\le1$, this
equation has no nonzero solution, while if $R_*>1$, it is shown to have
precisely one nonzero solution. A law of large numbers for the size of such a
large outbreak is proved by exploiting a single-type branching process that
approximates the size of the susceptibility set of a typical individual.


  Homologous recombination facilitates the exchange of genetic material between
homologous DNA molecules. This crucial process requires detecting a specific
homologous DNA sequence within a huge variety of heterologous sequences. The
detection is mediated by RecA in E. coli, or members of its superfamily in
other organisms. Here we examine how well is the RecA-DNA interaction adjusted
to its task. By formulating the DNA recognition process as a signal detection
problem, we find the optimal value of binding energy that maximizes the ability
to detect homologous sequences. We show that the experimentally observed
binding energy is nearly optimal. This implies that the RecA-induced
deformation and the binding energetics are fine-tuned to ensure optimal
sequence detection. Our analysis suggests a possible role for DNA extension by
RecA, in which deformation enhances detection. The present signal detection
approach provides a general recipe for testing the optimality of other
molecular recognition systems.


  We investigate evolutionary dynamics related to periodicity fossil
biodiversity. Coherent periodic fluctuation in origination/extinction of marine
genera that survive <45 million years is the source of an observed ~62 million
year periodicity analyzed in Paper I. We also show that the evolutionary
dynamics of "long-lived" genera (those that survive >45 million years) do not
participate in the periodic fluctuation in diversity and differ from those of
"short-lived" genera. The difference between the evolutionary dynamics of these
2 genera classes indicates that the periodic pattern is not an artifact of
variation in quality of the geologic record. The interplay of these two
previously undifferentiated systems, together with the secular increase in
abundance of "long-lived" genera, is probably the source of heretofore
unexplained differences in evolutionary dynamics between the Paleozoic and
post-Paleozoic as reported by others. Testing for cycles similar to the 62 Myr
cycle in fossil biodiversity superimposed on the long-term trends of the
Phanerozoic as described in Paper I, we find a significant (but weaker) signal
in sedimentary rock packages, particularly carbonates, which suggests a
connection. The presence of a periodic pattern in evolutionary dynamics of the
vulnerable "short-lived" component of marine fauna demonstrates that a
long-term periodic fluctuation in environmental conditions capable of affecting
evolution in the marine realm characterizes our planet. Coincidence in timing
is more consistent with a common cause than sampling bias. A previously
identified set of mass extinctions preferentially occur during the declining
phase of the 62 Myr periodicity, supporting the idea that the periodicity
relates to variation in biotically important stresses. Further work should
focus on finding links to physical phenomena that might reveal the causal
system or systems.


  Generalizing the cyclically competing three-species model (often referred to
as the rock-paper-scissors game), we consider a simple system of population
dynamics without spatial structures that involves four species. Unlike the
previous model, the four form alliance pairs which resemble partnership in the
game of Bridge. In a finite system with discrete stochastic dynamics, all but 4
of the absorbing states consist of coexistence of a partner-pair. From a master
equation, we derive a set of mean field equations of evolution. This approach
predicts complex time dependence of the system and that the surviving
partner-pair is the one with the larger product of their strengths (rates of
consumption). Simulations typically confirm these scenarios. Beyond that, much
richer behavior is revealed, including complicated extinction probabilities and
non-trivial distributions of the population ratio in the surviving pair. These
discoveries naturally raise a number of intriguing questions, which in turn
suggests a variety of future avenues of research, especially for more realistic
models of multispecies competition in nature.


  Different numerical mappings of the DNA sequences have been studied using a
new cluster-scaling method and the well known spectral methods. It is shown, in
particular, that the nucleotide sequences in DNA molecules have robust
cluster-scaling properties. These properties are relevant to both types of
nucleotide pair-bases interactions: hydrogen bonds and stacking interactions.
It is shown that taking into account the cluster-scaling properties can help to
improve heterogeneous models of the DNA dynamics. It is also shown that a
chaotic (deterministic) order, rather than a stochastic randomness, controls
the energy minima positions of the stacking interactions in the DNA sequences
on large scales. The chaotic order results in a large-scale chaotic coherence
between the two complimentary DNA-duplex's sequences. A competition between
this broad-band chaotic coherence and the resonance coherence produced by
genetic code has been briefly discussed. The Arabidopsis plant genome (which is
a model plant for genome analysis) and two human genes: BRCA2 and NRXN1, have
been considered as examples.


  We describe an approach to model genetic regulatory networks at the level of
promotion-inhibition circuitry through a class of stochastic spin models that
includes spatial and temporal density fluctuations in a natural way. The
formalism can be viewed as an agent-based model formalism with agent behavior
ruled by a classical spin-like pseudo-Hamiltonian playing the role of a local,
individual objective function. A particular but otherwise generally applicable
choice for the microscopic transition rates of the models also makes them of
independent interest. To illustrate the formalism, we investigate (by Monte
Carlo simulations) some stationary state properties of the repressilator, a
synthetic three-gene network of transcriptional regulators that possesses
oscillatory behavior.


  Phylogenetic tree reconstruction is traditionally based on multiple sequence
alignments (MSAs) and heavily depends on the validity of this information
bottleneck. With increasing sequence divergence, the quality of MSAs decays
quickly. Alignment-free methods, on the other hand, are based on abstract
string comparisons and avoid potential alignment problems. However, in general
they are not biologically motivated and ignore our knowledge about the
evolution of sequences. Thus, it is still a major open question how to define
an evolutionary distance metric between divergent sequences that makes use of
indel information and known substitution models without the need for a multiple
alignment. Here we propose a new evolutionary distance metric to close this
gap. It uses finite-state transducers to create a biologically motivated
similarity score which models substitutions and indels, and does not depend on
a multiple sequence alignment. The sequence similarity score is defined in
analogy to pairwise alignments and additionally has the positive semi-definite
property. We describe its derivation and show in simulation studies and
real-world examples that it is more accurate in reconstructing phylogenies than
competing methods. The result is a new and accurate way of determining
evolutionary distances in and beyond the twilight zone of sequence alignments
that is suitable for large datasets.


  We investigate the error threshold for the emergence of quasispecies in the
Eigen model. By mapping to to an effective Hamiltonian ruled by the
"imaginary-time" Schr\"odinger equation, a variational ansatz is proposed and
applied to calculate various quantities associated with the quasispecies. The
variational ansatz gives correct predictions for the survival population of the
wild-type sequence and also reveals an unexpected universal scaling behaviors
near the error threshold. We check the validity of the variational ansatz by
numerical methods and find excellent agreement. Though the emergence of the
scaling behaviors is not yet fully understood, it is remarkable that the
universal scaling function reigns even for relatively short genome length such
as L=16. Further investigations may reveal the mechanism of the universal
scaling and extract the essential ingredients for the emergence of the
quasispecies in molecular evolution.


  To be the fittest is central to proliferation in evolutionary games.
Individuals thus adopt the strategies of better performing players in the hope
of successful reproduction. In structured populations the array of those that
are eligible to act as strategy sources is bounded to the immediate neighbors
of each individual. But which one of these strategy sources should potentially
be copied? Previous research dealt with this question either by selecting the
fittest or by selecting one player uniformly at random. Here we introduce a
parameter $u$ that interpolates between these two extreme options. Setting $u$
equal to zero returns the random selection of the opponent, while positive $u$
favor the fitter players. In addition, we divide the population into two
groups. Players from group $A$ select their opponents as dictated by the
parameter $u$, while players from group $B$ do so randomly irrespective of $u$.
We denote the fraction of players contained in groups $A$ and $B$ by $v$ and
$1-v$, respectively. The two parameters $u$ and $v$ allow us to analyze in
detail how aspirations in the context of the prisoner's dilemma game influence
the evolution of cooperation. We find that for sufficiently positive values of
$u$ there exist a robust intermediate $v \approx 0.5$ for which cooperation
thrives best. The robustness of this observation is tested against different
levels of uncertainty in the strategy adoption process $K$ and for different
interaction networks. We also provide complete phase diagrams depicting the
dependence of the impact of $u$ and $v$ for different values of $K$, and
contrast the validity of our conclusions by means of an alternative model where
individual aspiration levels are subject to evolution as well. Our study
indicates that heterogeneity in aspirations may be key for the sustainability
of cooperation in structured populations.


  The consistency of the species abundance distribution across diverse
communities has attracted widespread attention. In this paper, I argue that the
consistency of pattern arises because diverse ecological mechanisms share a
common symmetry with regard to measurement scale. By symmetry, I mean that
different ecological processes preserve the same measure of information and
lose all other information in the aggregation of various perturbations. I frame
these explanations of symmetry, measurement, and aggregation in terms of a
recently developed extension to the theory of maximum entropy. I show that the
natural measurement scale for the species abundance distribution is log-linear:
the information in observations at small population sizes scales
logarithmically and, as population size increases, the scaling of information
grades from logarithmic to linear. Such log-linear scaling leads naturally to a
gamma distribution for species abundance, which matches well with the observed
patterns. Much of the variation between samples can be explained by the
magnitude at which the measurement scale grades from logarithmic to linear.
This measurement approach can be applied to the similar problem of allelic
diversity in population genetics and to a wide variety of other patterns in
biology.


  In this work we study, at the single molecular level, the thermodynamic and
dynamic characteristics of an enzymatic reaction comprising a rate limiting
step. We investigate how the stability of the enzyme-state stationary
probability distribution, the reaction velocity, and its efficiency of energy
conversion depend on the system parameters. We employ in this study a recently
introduced formalism for performing a multiscale thermodynamic analysis in
continuous-time discrete-state stochastic systems.


  Diversification is nested, and early models suggested this could lead to a
great deal of evolutionary redundancy in the Tree of Life. This result is based
on a particular set of branch lengths produced by the common coalescent, where
pendant branches leading to tips can be very short compared to branches deeper
in the tree. Here, we analyze alternative and more realistic Yule and
birth-death models. We show how censoring at the present both makes average
branches one half what we might expect and makes pendant and interior branches
roughly equal in length. Although dependent on whether we condition on the size
of the tree, its age, or both, these results hold both for the Yule model and
for birth-death models with moderate extinction. Importantly, the rough
equivalency in interior and exterior branch lengths means the loss of
evolutionary history with loss of species can be roughly linear. Under these
models, the Tree of Life may offer limited redundancy in the face of ongoing
species loss.


  We answer, in the affirmative, the following question proposed by Mike Steel
as a $100 challenge: "Is the following problem NP-hard? Given a ternary
phylogenetic X-tree T and a collection Q of quartet subtrees on X, is T the
only tree that displays Q ?"


  The patterns of motion of mobile agents has received recently wide attention
in the literature. There is a number of recent studies centered around the
motion behavior of many agents ranging from albatrosses to human beings.
Special attention has been given to the covered distances statistical
distributions. In some cases, due to the lack of accurate data about the motion
of the agents it has been necessary to plan very clever experiments to obtain
them. These experiments try to infer the statistical properties of the agents'
real motion from the observed positions in consecutive time intervals. The
length of the time intervals is a random variable taking values from a
previously known statistical distribution or from a distribution deduced from
empirical data. The aim of this work is to demonstrate that for a Gaussian
Random Walker it is, in general, impossible to recover the real motion patterns
distribution from the stroboscopic observation of the agents. Moreover, it is
also shown that the distances distribution strongly depends on the agents'
observation time intervals. These claims are sustained by numerical
experiments.


  The science of networks has revolutionised research into the dynamics of
interacting elements. It could be argued that epidemiology in particular has
embraced the potential of network theory more than any other discipline. Here
we review the growing body of research concerning the spread of infectious
diseases on networks, focusing on the interplay between network theory and
epidemiology. The review is split into four main sections, which examine: the
types of network relevant to epidemiology; the multitude of ways these networks
can be characterised; the statistical methods that can be applied to infer the
epidemiological parameters on a realised network; and finally simulation and
analytical methods to determine epidemic dynamics on a given network. Given the
breadth of areas covered and the ever-expanding number of publications, a
comprehensive review of all work is impossible. Instead, we provide a
personalised overview into the areas of network epidemiology that have seen the
greatest progress in recent years or have the greatest potential to provide
novel insights. As such, considerable importance is placed on analytical
approaches and statistical methods which are both rapidly expanding fields.
Throughout this review we restrict our attention to epidemiological issues.


  Stress fibers are contractile actomyosin bundles commonly observed in the
cytoskeleton of metazoan cells. The spatial profile of the polarity of actin
filaments inside contractile actomyosin bundles is either monotonic (graded) or
periodic (alternating). In the framework of linear irreversible thermodynamics,
we write the constitutive equations for a polar, active, elastic
one-dimensional medium. An analysis of the resulting equations for the dynamics
of polarity shows that the transition from graded to alternating polarity
patterns is a nonequilibrium Lifshitz point. Active contractility is a
necessary condition for the emergence of sarcomeric, alternating polarity
patterns.


  Inferring dynamic biochemical networks is one of the main challenges in
systems biology. Given experimental data, the objective is to identify the
rules of interaction among the different entities of the network. However, the
number of possible models fitting the available data is huge and identifying a
biologically relevant model is of great interest. Nested canalyzing functions,
where variables in a given order dominate the function, have recently been
proposed as a framework for modeling gene regulatory networks. Previously we
described this class of functions as an algebraic toric variety. In this paper,
we present an algorithm that identifies all nested canalyzing models that fit
the given data. We demonstrate our methods using a well-known Boolean model of
the cell cycle in budding yeast.


  The concept of a carrying capacity is essential in most models to prevent
unlimited growth. Despite the large amount of deaths it introduces, the actual
influence of the Verhulst term in simulations is often times not accounted for.
Generally, it is treated merely as a scaling parameter that functions to keep
simulated populations within computer limits. Here, we compare two different
implementations of the concept in the Penna model - Vehulst applied to all
individuals (VA) and to newborns only (VB). We observe variations in certain
model features when random Verhulst deaths are restricted to a single age
group.


  We present two models for electron-driven uphill proton transport across
lipid membranes, with the electron energy converted to the proton gradient via
the electrostatic interaction. In the first model, associated with the
cytochrome c oxidase complex in the inner mitochondria membranes, the
electrostatic coupling to the site occupied by an electron lowers the energy
level of the proton-binding site, making the proton transfer possible. In the
second model, roughly describing the redox loop in a nitrate respiration of E.
coli bacteria, an electron displaces a proton from the negative side of the
membrane to a shuttle, which subsequently diffuses across the membrane and
unloads the proton to its positive side. We show that both models can be
described by the same approach, which can be significantly simplified if the
system is separated into several clusters, with strong Coulomb interaction
inside each cluster and weak transfer couplings between them. We derive and
solve the equations of motion for the electron and proton creation/annihilation
operators, taking into account the appropriate Coulomb terms, tunnel couplings,
and the interaction with the environment. For the second model, these equations
of motion are solved jointly with a Langevin-type equation for the shuttle
position. We obtain expressions for the electron and proton currents and
determine their dependence on the electron and proton voltage build-ups,
on-site charging energies, reorganization energies, temperature, and other
system parameters. We show that the quantum yield in our models can be up to
100% and the power-conversion efficiency can reach 35%.


  Over the past decade, a number of researchers in systems biology have sought
to relate the function of biological systems to their network-level
descriptions -- lists of the most important players and the pairwise
interactions between them. Both for large networks (in which statistical
analysis is often framed in terms of the abundance of repeated small subgraphs)
and for small networks which can be analyzed in greater detail (or even
synthesized in vivo and subjected to experiment), revealing the relationship
between the topology of small subgraphs and their biological function has been
a central goal. We here seek to pose this revelation as a statistical task,
illustrated using a particular setup which has been constructed experimentally
and for which parameterized models of transcriptional regulation have been
studied extensively. The question "how does function follow form" is here
mathematized by identifying which topological attributes correlate with the
diverse possible information-processing tasks which a transcriptional
regulatory network can realize. The resulting method reveals one form-function
relationship which had earlier been predicted based on analytic results, and
reveals a second for which we can provide an analytic interpretation. Resulting
source code is distributed via http://formfunction.sourceforge.net.


  Methods of phylogenetic inference use more and more complex models to
generate trees from data. However, even simple models and their implications
are not fully understood.
  Here, we investigate the two-state Markov model on a tripod tree, inferring
conditions under which a given set of observations gives rise to such a model.
This type of investigation has been undertaken before by several scientists
from different fields of research.
  In contrast to other work we fully analyse the model, presenting conditions
under which one can infer a model from the observation or at least get support
for the tree-shaped interdependence of the leaves considered.
  We also present all conditions under which the results can be extended from
tripod trees to quartet trees, a step necessary to reconstruct at least a
topology. Apart from finding conditions under which such an extension works we
discuss example cases for which such an extension does not work.


  Humans and other animals can adapt their social behavior in response to
environmental cues including the feedback obtained through experience.
Nevertheless, the effects of the experience-based learning of players in
evolution and maintenance of cooperation in social dilemma games remain
relatively unclear. Some previous literature showed that mutual cooperation of
learning players is difficult or requires a sophisticated learning model. In
the context of the iterated Prisoner's Dilemma, we numerically examine the
performance of a reinforcement learning model. Our model modifies those of
Karandikar et al. (1998), Posch et al. (1999), and Macy and Flache (2002) in
which players satisfice if the obtained payoff is larger than a dynamic
threshold. We show that players obeying the modified learning mutually
cooperate with high probability if the dynamics of threshold is not too fast
and the association between the reinforcement signal and the action in the next
round is sufficiently strong. The learning players also perform efficiently
against the reactive strategy. In evolutionary dynamics, they can invade a
population of players adopting simpler but competitive strategies. Our version
of the reinforcement learning model does not complicate the previous model and
is sufficiently simple yet flexible. It may serve to explore the relationships
between learning and evolution in social dilemma situations.


  Bacteriophages, phages for short, are viruses of bacteria. The majority of
phages contain a double-stranded DNA genome packaged in a capsid at a density
of ~500 mg/ml. This high density requires substantial compression of the normal
B form helix, leading to the conjecture that DNA in mature phage virions is
under significant pressure, and that pressure is used to eject the DNA during
infection. A large number of theoretical, computer simulation and in vitro
experimental studies surrounding this conjecture has revealed many --- though
often isolated and/or contradictory --- aspects of packaged DNA. This prompts
us to present a unified view of the statistical physics and thermodynamics of
DNA packaged in phage capsids. We argue that the DNA in a mature phage is in a
(meta)stable state, wherein electrostatic self-repulsion is balanced by
curvature stress due to confinement in the capsid. We show that in addition to
the osmotic pressure associated with the packaged DNA and its counterions,
there are four different pressures within the capsid: pressure on the DNA,
hydrostatic pressure, the pressure experienced by the capsid, and the pressure
associated with the chemical potential of DNA ejection. Significantly, we
analyze the mechanism of force transmission in the packaged DNA, and
demonstrate that the pressure on DNA is not important for ejection. We derive
equations showing a strong hydrostatic pressure difference across the capsid
shell. We propose that when a phage is triggered to eject by interaction with
its receptor in vitro, the (thermodynamic) incentive of water molecules to
enter the phage capsid flushes the DNA out of the capsid. In vivo, the
difference between the osmotic pressures in the bacterial cell cytoplasm and
the culture medium similarly results in a water flow that drags the DNA out of
the capsid and into the bacterial cell.


  When microbes compete for limited resources, they often engage in chemical
warfare using bacterial toxins. This competition can be understood in terms of
evolutionary game theory (EGT). We study the predictions of EGT for the
bacterial "suicide bomber" game in terms of the phase portraits of population
dynamics, for parameter combinations that cover all interesting games for
two-players, and seven of the 38 possible phase portraits of the three-player
game. We compare these predictions to simulations of these competitions in
finite well-mixed populations, but also allowing for probabilistic rather than
pure strategies, as well as Darwinian adaptation over tens of thousands of
generations. We find that Darwinian evolution of probabilistic strategies
stabilizes games of the rock-paper-scissors type that emerge for parameters
describing realistic bacterial populations, and point to ways in which the
population fixed point can be selected by changing those parameters.


  Crosstalk is defined as the set of unwanted interactions among the different
entities of a network. Crosstalk is present in various degrees in every system
where information is transmitted through a means that is accessible by all the
individual units of the network. Using concepts from graph theory, we introduce
a quantifiable measure for sensitivity to crosstalk, and analytically derive
the structure of the networks in which it is minimized. It is shown that
networks with an inhomogeneous degree distribution are more robust to crosstalk
than corresponding homogeneous networks. We provide a method to construct the
graph with the minimum possible sensitivity to crosstalk, given its order and
size. Finally, for networks with a fixed degree sequence, we present an
algorithm to find the optimal interconnection structure among their vertices.


  RSC is an essential chromatin remodeling factor that is required for the
control of several processes including transcription, repair and replication.
The ability of RSC to relocate centrally positioned mononucleosomes at the end
of nucleosomal DNA is firmly established, but the data on RSC action on
oligo-nucleosomal templates remains still scarce. By using Atomic Force
Microscopy (AFM) imaging, we have quantitatively studied the RSC- induced
mobilization of positioned di- and trinucleosomes as well as the directionality
of mobilization on mononucleosomal template labeled at one end with
streptavidin. AFM imaging showed only a limited set of distinct configurational
states for the remodeling products. No stepwise or preferred directionality of
the nucleosome motion was observed. Analysis of the corresponding reaction
pathways allows deciphering the mechanistic features of RSC-induced nucleosome
relocation. The final outcome of RSC remodeling of oligosome templates is the
packing of the nucleosomes at the edge of the template, providing large
stretches of DNA depleted of nucleosomes. This feature of RSC may be used by
the cell to overcome the barrier imposed by the presence of nucleosomes.


  Metabolic networks have two properties that are generally regarded as
unrelated: One, they have metabolic reactions whose single knockout is lethal
for the organism, and two, they have correlated sets of reactions forming
functional modules. In this review we argue that both essentiality and
modularity seem to arise as a consequence of the same structural property: the
existence of low degree metabolites. This observation allows a prediction of
(a) essential metabolic reactions which are potential drug targets in
pathogenic microorganisms and (b) regulatory modules within biological
networks, from purely structural information about the metabolic network.


  In the small phylogeny problem we, are given a phylogenetic tree and gene
orders of the extant species and our goal is to reconstruct all of the
ancestral genomes so that the number of evolutionary operations is minimized.
Algorithms for reconstructing evolutionary history from gene orders are usually
based on repeatedly computing medians of genomes at neighbouring vertices of
the tree. We propose a new, more general approach, based on an iterative local
optimization procedure. In each step, we propose candidates for ancestral
genomes and choose the best ones by dynamic programming. We have implemented
our method and used it to reconstruct evolutionary history of 16 yeast mtDNAs
and 13 Campanulaceae cpDNAs.


  The neocortex is widely believed to be the seat of intelligence and "mind".
However, it's unclear what "mind" is, or how the special features of neocortex
enable it, though likely "connectionist" principles are involved *A. The key to
intelligence1 is learning relationships between large numbers of signals (such
as pixel values), rather than memorizing explicit patterns. Causes (such as
objects) can then be inferred from a learned internal model. These
relationships fall into 2 classes: simple pairwise or second-order correlations
(socs), and complex, and vastly more numerous, higher-order correlations
(hocsB), such as the product of 3 or more pixels averaged over a set of images.
Thus if 3 pixels correlate, they may give an "edge". Neurons with "Hebbian"
synapses (changing strength in response to input-output spike-coincidences) are
sensitive to such correlations, and it's likely that learned internal models
use such neurons. Because output firing depends on input firing via the
relevant connection strengths, Hebbian learning provides, in a feedback manner,
sensitivity to input correlations. Hocs are vital, since they express
"interesting" structure2 (e.g. edges), but their detection requires nonlinear
rules operating at synapses of individual neurons. Here we report that in
single model neurons learning from hocs fails, and defaults to socs, if
nonlinear Hebbian rules are not sufficiently connection-specific. Such failure
would inevitably occur if a neuron's input synapses were too crowded, and would
undermine biological connectionism. Since the cortex must be hoc-sensitive to
achieve the type of learning enabling mind, we propose it uses known, detailed
but poorly understood circuitry and physiology to "proofread" Hebbian
connections. Analogous DNA proofreading allows evolution of complex genomes
(i.e. "life").


  The electrochemical behaviour of the biomedical and metallic alloys,
especially in the orthopaedic implants fields, raises many questions. This
study is dedicated for studying the Ti-6Al-4V alloy, by electrochemical
impedance spectroscopy, EIS, in various physiological media,: Ringer solution,
phosphate buffered solution (PBS), PBS solution and albumin, PBS solution with
calf serum and PBS solution with calf serum and an antioxidant (sodium azide).
Moreover, the desionised water was considered as the reference solution. The
tests reproducibility was investigated. The time-frequency-Module graphs
highlighted that the desionised water is the most protective for the Ti-6Al-4V
alloy. This biomedical alloy is the less protected in the solution constituted
by PBS and albumin. The time-frequency graph allows pointing out the graphic
signatures of adsorption for organic and inorganic species (differences between
the modules means in studied solution and the modules mean in the reference
solution). --- Le comportement \'electrochimique des alliages m\'etalliques
biom\'edicaux, notamment dans le domaine des implants orthop\'ediques, pose
encore de nombreuses questions. Ce travail propose d'\'etudier l'alliage de
titane Ti-6Al-4V, par spectroscopie d'imp\'edance \'electrochimique, SIE, dans
diff\'erents milieux physiologiques : solution de Ringer, solution \`a base
d'un tampon phosphate (PBS), solution PBS avec de l'albumine, solution PBS avec
du s\'erum bovin et une solution PBS avec du s\'erum bovin et un antioxydant
(azoture de sodium). De plus, une solution d'eau ultra-pure servira de
r\'ef\'erence. La reproductibilit\'e des tests a \'et\'e \'etudi\'ee. Les
repr\'esentations temps-fr\'equence des modules ont mis en \'evidence que l'eau
d\'esionis\'ee est la solution qui pr\'esente le caract\`ere le plus protecteur
pour le Ti-6Al-4V. Cet alliage de titane est le moins prot\'eg\'e dans la
solution de PBS contenant de l'albumine. Cette repr\'esentation permet de
mettre en \'evidence des signatures graphiques d'adsorption des esp\`eces
inorganiques et organiques (diff\'erences entre les moyennes des modules dans
les solutions \'etudi\'ees et la moyenne des modules dans la solution de
r\'ef\'erence).


  The dynamics of well-mixed biological populations is usually studied by
mean-field methods and weak-noise expansions. Similar methods have been applied
also in spatially extended problems, relying on the fact that these populations
are organized in colonies with a large local density of individuals. We provide
a counterexample discussing a one-dimensional neutral population with negative
frequency-dependent selection. The system exhibits a continuous phase
transition between genetic fixation and coexistence unexpected from weak-noise
arguments. We show that the behavior is a non-perturbative effect of the
internal noise that is amplified by presence of spatial correlations
(strong-noise regime).


  We introduce a dengue model (SEIR) where the human individuals are treated on
an individual basis (IBM) while the mosquito population, produced by an
independent model, is treated by compartments (SEI). We study the spread of
epidemics by the sole action of the mosquito. Exponential, deterministic and
experimental distributions for the (human) exposed period are considered in two
weather scenarios, one corresponding to temperate climate and the other to
tropical climate. Virus circulation, final epidemic size and duration of
outbreaks are considered showing that the results present little sensitivity to
the statistics followed by the exposed period provided the median of the
distributions are in coincidence. Only the time between an introduced
(imported) case and the appearance of the first symptomatic secondary case is
sensitive to this distribution. We finally show that the IBM model introduced
is precisely a realization of a compartmental model, and that at least in this
case, the choice between compartmental models or IBM is only a matter of
convenience.


  Networks coming from protein-protein interactions, transcriptional
regulation, signaling, or metabolism may appear to have "unusual" properties.
To quantify this, it is appropriate to randomize the network and test the
hypothesis that the network is not statistically different from expected in a
motivated ensemble. However, when dealing with metabolic networks, the
randomization of the network using edge exchange generates fictitious reactions
that are biochemically meaningless. Here we provide several natural ensembles
of randomized metabolic networks. A first constraint is to use valid
biochemical reactions. Further constraints correspond to imposing appropriate
functional constraints. We explain how to perform these randomizations with the
help of Markov Chain Monte Carlo (MCMC) and show that they allow one to
approach the properties of biological metabolic networks. The implication of
the present work is that the observed global structural properties of real
metabolic networks are likely to be the consequence of simple biochemical and
functional constraints.


  We present a stochastic approach for ion transport at the mesoscopic level.
The description takes into account the self-consistent electric field generated
by the fixed and mobile charges as well as the discrete nature of these latter.
As an application we study the noise in the ion transport process, including
the effect of the displacement current generated by the fluctuating electric
field. The fluctuation theorem is shown to hold for the electric current with
and without the displacement current.


  We study the evolution of a population in a two-locus genotype space, in
which the negative effects of two single mutations are overcompensated in a
high fitness double mutant. We discuss how the interplay of finite population
size, $N$, and sexual recombination at rate $r$ affects the escape times
$t_\mathrm{esc}$ to the double mutant. For small populations demographic noise
generates massive fluctuations in $t_\mathrm{esc}$. The mean escape time varies
non-monotonically with $r$, and grows exponentially as $\ln t_{\mathrm{esc}}
\sim N(r - r^\ast)^{3/2}$ beyond a critical value $r^\ast$.


  There has been much recent interest in modelling epidemics on networks,
particularly in the presence of substantial clustering. Here, we develop
pairwise methods to answer questions that are often addressed using epidemic
models, in particular: on the basis of potential observations early in an
outbreak, what can be predicted about the epidemic outcomes and the levels of
intervention necessary to control the epidemic? We find that while some results
are independent of the level of clustering (early growth predicts the level of
`leaky' vaccine needed for control and peak time, while the basic reproductive
ratio predicts the random vaccination threshold) the relationship between other
quantities is very sensitive to clustering.


  The biological world, especially its majority microbial component, is
strongly interacting and may be dominated by collective effects. In this
review, we provide a brief introduction for statistical physicists of the way
in which living cells communicate genetically through transferred genes, as
well as the ways in which they can reorganize their genomes in response to
environmental pressure. We discuss how genome evolution can be thought of as
related to the physical phenomenon of annealing, and describe the sense in
which genomes can be said to exhibit an analogue of information entropy. As a
direct application of these ideas, we analyze the variation with ocean depth of
transposons in marine microbial genomes, predicting trends that are consistent
with recent observations using metagenomic surveys.


  Our aim is to detect mechanistic interaction between the effects of two
causal factors on a binary response, as an aid to identifying situations where
the effects are mediated by a common mechanism. We propose a formalization of
mechanistic interaction which acknowledges asymmetries of the kind "factor A
interferes with factor B, but not viceversa". A class of tests for mechanistic
interaction is proposed, which works on discrete or continuous causal
variables, in any combination. Conditions under which these tests can be
applied under a generic regime of data collection, be it interventional or
observational, are discussed in terms of conditional independence assumptions
within the framework of Augmented Directed Graphs. The scientific relevance of
the method and the practicality of the graphical framework are illustrated with
the aid of two studies in coronary artery disease. Our analysis relies on the
"deep determinism" assumption that there exists some relevant set V - possibly
unobserved - of "context variables", such that the response Y is a
deterministic function of the values of V and of the causal factors of
interest. Caveats regarding this assumption in real studies are discussed.


  Genomic distance between two genomes, i.e., the smallest number of genome
rearrangements required to transform one genome into the other, is often used
as a measure of evolutionary closeness of the genomes in comparative genomics
studies. However, in models that include rearrangements of significantly
different "power" such as reversals (that are "weak" and most frequent
rearrangements) and transpositions (that are more "powerful" but rare), the
genomic distance typically corresponds to a transformation with a large
proportion of transpositions, which is not biologically adequate.
  Weighted genomic distance is a traditional approach to bounding the
proportion of transpositions by assigning them a relative weight {\alpha} > 1.
A number of previous studies addressed the problem of computing weighted
genomic distance with {\alpha} \leq 2.
  Employing the model of multi-break rearrangements on circular genomes, that
captures both reversals (modelled as 2-breaks) and transpositions (modelled as
3-breaks), we prove that for {\alpha} \in (1,2], a minimum-weight
transformation may entirely consist of transpositions, implying that the
corresponding weighted genomic distance does not actually achieve its purpose
of bounding the proportion of transpositions. We further prove that for
{\alpha} \in (1,2), the minimum-weight transformations do not depend on a
particular choice of {\alpha} from this interval. We give a complete
characterization of such transformations and show that they coincide with the
transformations that at the same time have the shortest length and make the
smallest number of breakages in the genomes.
  Our results also provide a theoretical foundation for the empirical
observation that for {\alpha} < 2, transpositions are favored over reversals in
the minimum-weight transformations.


  The internal motions of DNA immersed in bio-fluid are investigated. The
interactions between the fragments of DNA and the surrounding bio-fluid are
modeled using the gauge fluid lagrangian. In the model, the bio-fluid is
coupled to the standard gauge invariant bosonic lagrangian describing the DNA.
It is shown that at non-relativistic limit various equation of motions, from
the well-known Sine-Gordon equation to the simultaneous nonlinear equations,
can be constructed within a single framework. The effects of bio-fluid are
investigated for two cases : single and double stranded DNA. It is argued that
the small and large amplitudes of a single stranded DNA motion immersed in
bio-fluid can be explained in a natural way within the model as a solitonic
wave regardless with the fluid velocity. In contrary the double stranded DNA
behaves as regular or damped harmonic oscillator and is highly depending on the
fluid velocity.


  X-ray crystallography is a powerful tool to determine the protein 3D
structure. However, it is time-consuming and expensive, and not all proteins
can be successfully crystallized, particularly for membrane proteins. Although
nuclear magnetic resonance (NMR) spectroscopy is indeed a very powerful tool in
determining the 3D structures of membrane proteins, it is also time-consuming
and costly. To the best of the authors' knowledge, there is little structural
data available on the AGAAAAGA palindrome in the hydrophobic region (113-120)
of prion proteins due to the noncrystalline and insoluble nature of the amyloid
fibril, although many experimental studies have shown that this region has
amyloid fibril forming properties and plays an important role in prion
diseases. In view of this, the present study is devoted to address this problem
from computational approaches such as global energy optimization, simulated
annealing, and structural bioinformatics. The optimal atomic-resolution
structures of prion AGAAAAGA amyloid fibils reported in this paper have a value
to the scientific community in its drive to find treatments for prion diseases.


  The widespread relevance of increasingly complex networks requires methods to
extract meaningful coarse-grained representations of such systems. For
undirected graphs, standard community detection methods use criteria largely
based on density of connections to provide such representations. We propose a
method for grouping nodes in directed networks based on the role of the nodes
in the network, understood in terms of patterns of incoming and outgoing flows.
The role groupings are obtained through the clustering of a similarity matrix,
formed by the distances between feature vectors that contain the number of in
and out paths of all lengths for each node. Hence nodes operating in a similar
flow environment are grouped together although they may not themselves be
densely connected. Our method, which includes a scale factor that reveals
robust groupings based on increasingly global structure, provides an
alternative criterion to uncover structure in networks where there is an
implicit flow transfer in the system. We illustrate its application to a
variety of data from ecology, world trade and cellular metabolism.


  We analyze theoretically the effects of excluded-volume interactions between
motors on the dynamics of a cargo driven by multiple motors. The model
considered shares many commons with other recently proposed in the literature,
with the addition of direct interaction between motors and motor back steps.
The cargo is assumed to follow a continuum Langevin dynamics, while individual
motors evolve following a Monte Carlo algorithm based on experimentally
accessible probabilities for discrete forward and backward jumps, and
attachment and detachment rates. The links between cargo and motors are
considered as non linear springs. By means of numerical simulations we compute
the relevant quantities characterizing the dynamical properties of the system,
and we compare the results to those for non interacting motors. We find that
interactions lead to quite relevant changes in the force-velocity relation for
cargo, with a considerable reduction of the stall force, and cause also a
notable decrease of the run length. These effects are mainly due to
traffic-like phenomena in the microtubule. The consideration of several
parallel tracks for motors reduces such effects. However, we find that for
realistic values of the number of motors and the number of tracks, the
influence of interactions on the global parameters of transport of cargo are
far from being negligible. Our studies provide also an analysis of the
relevance of motor back steps on the modeling, and of the influence of
different assumptions for the detachment rates. In particular, we discuss these
two aspects in connection with the possibility of observing processive back
motion of cargo at large load forces.


  The interplay between mutation and selection plays a fundamental role in the
behaviour of evolutionary algorithms (EAs). However, this interplay is still
not completely understood. This paper presents a rigorous runtime analysis of a
non-elitist population-based EA that uses the linear ranking selection
mechanism. The analysis focuses on how the balance between parameter $\eta$,
controlling the selection pressure in linear ranking, and parameter $\chi$
controlling the bit-wise mutation rate, impacts the runtime of the algorithm.
The results point out situations where a correct balance between selection
pressure and mutation rate is essential for finding the optimal solution in
polynomial time. In particular, it is shown that there exist fitness functions
which can only be solved in polynomial time if the ratio between parameters
$\eta$ and $\chi$ is within a narrow critical interval, and where a small
change in this ratio can increase the runtime exponentially. Furthermore, it is
shown quantitatively how the appropriate parameter choice depends on the
characteristics of the fitness function. In addition to the original results on
the runtime of EAs, this paper also introduces a very useful analytical tool,
i.e., multi-type branching processes, to the runtime analysis of non-elitist
population-based EAs.


  For a single enzyme or molecular motor operating in an aqueous solution of
non-equilibrated solute concentrations, a thermodynamic description is
developed on the level of an individual trajectory of transitions between
states. The concept of internal energy, intrinsic entropy and free energy for
states follows from a microscopic description using one assumption on
time-scale separation. A first law energy balance then allows the unique
identification of the heat dissipated in one transition. Consistency with the
second law on the ensemble level enforces both stochastic entropy as third
contribution to the entropy change involved in one transition and the local
detailed balance condition for the ratio between forward and backward rates for
any transition. These results follow without assuming weak coupling between the
enzyme and the solutes, ideal solution behavior or mass action law kinetics.
The present approach highlights both the crucial role of the intrinsic entropy
of each state and the physically questionable role of chemiostats for deriving
the first law for molecular motors subject to an external force under realistic
conditions.


  We describe a simple conductance-based model neuron that includes intra- and
extra-cellular ion concentration dynamics and show that this model exhibits
periodic bursting. The bursting arises as the fast spiking behavior of the
neuron is modulated by the slow oscillatory behavior in the ion concentration
variables, and vice versa. By separating these time scales and studying the
bifurcation structure of the neuron, we catalog several qualitatively different
bursting profiles that are strikingly similar to those seen in experimental
preparations. Our work suggests that ion concentration dynamics may play an
important role in modulating neuronal excitability in real biological systems.


  The quasi-species equation describes the evolution of the probability that a
random individual in a population carries a given genome. Here we map the
quasi-species equation for individuals of a self-reproducing population to an
ensemble of scalar field elementary units undergoing a creation and
annihilation process. In this mapping, the individuals of the population are
mapped to field units and their genome to the field value. The selective
pressure is mapped to an inverse temperature \beta of the system regulating the
evolutionary dynamics of the fields. We show that the quasi-species equation if
applied to an ensemble of field units gives in the small \beta limit can be put
in relation with existing stochastic quantization approaches. The ensemble of
field units described by the quasi-species equation relaxes to the fundamental
state, describing an intrinsically dissipative dynamics. For a quadratic
dispersion relation the mean energy <U> of the system changes as a function of
the inverse temperature \beta. For small values of \beta the average energy <U>
takes a relativistic form, for large values of \beta, the average energy <U>
takes a classical form.


  Spin models of neural networks and genetic networks are considered elegant as
they are accessible to statistical mechanics tools for spin glasses and
magnetic systems. However, the conventional choice of variables in spin systems
may cause problems in some models when parameter choices are unrealistic from a
biological perspective. Obviously, this may limit the role of a model as a
template model for biological systems. Perhaps less obviously, also ensembles
of random networks are affected and may exhibit different critical properties.
We consider here a prototypical network model that is biologically plausible in
its local mechanisms. We study a discrete dynamical network with two
characteristic properties: Nodes with binary states 0 and 1, and a modified
threshold function with $\Theta_0(0)=0$. We explore the critical properties of
random networks of such nodes and find a critical connectivity $K_c=2.0$ with
activity vanishing at the critical point. Finally, we observe that the present
model allows a more natural implementation of recent models of budding yeast
and fission yeast cell-cycle control networks.


  We introduce a one-parametric family of tree growth models, in which
branching probabilities decrease with branch age $\tau$ as $\tau^{-\alpha}$.
Depending on the exponent $\alpha$, the scaling of tree depth with tree size
$n$ displays a transition between the logarithmic scaling of random trees and
an algebraic growth. At the transition ($\alpha=1$) tree depth grows as $(\log
n)^2$. This anomalous scaling is in good agreement with the trend observed in
evolution of biological species, thus providing a theoretical support for
age-dependent speciation and associating it to the occurrence of a critical
point.


  We study the role of quantum fluctuations of atomic nuclei in the real-time
dynamics of non-equilibrium macro-molecular transitions. To this goal we
introduce an extension of the Dominant Reaction Pathways (DRP) formalism, in
which the quantum corrections to the classical overdamped Langevin dynamics are
rigorously taken into account to order h^2 . We first illustrate our approach
in simple cases, and compare with the results of the instanton theory. Then we
apply our method to study the C7_eq to C7_ax transition of alanine dipeptide.
We find that the inclusion of quantum fluctuations can significantly modify the
reaction mechanism for peptides. For example, the energy difference which is
overcome along the most probable pathway is reduced by as much as 50%.


  Gene regulation involves a hierarchy of events that extend from specific
protein-DNA interactions to the combinatorial assembly of nucleoprotein
complexes. The effects of DNA sequence on these processes have typically been
studied based either on its quantitative connection with single-domain binding
free energies or on empirical rules that combine different DNA motifs to
predict gene expression trends on a genomic scale. The middle-point approach
that quantitatively bridges these two extremes, however, remains largely
unexplored. Here, we provide an integrated approach to accurately predict gene
expression from statistical sequence information in combination with detailed
biophysical modeling of transcription regulation by multidomain binding on
multiple DNA sites. For the regulation of the prototypical lac operon, this
approach predicts within 0.3-fold accuracy transcriptional activity over a
10,000-fold range from DNA sequence statistics for different intracellular
conditions.


  Recently, we have proposed a redox molecular hypothesis about the natural
biophysical substrate of visual perception and imagery (B\'okkon, 2009.
BioSystems; B\'okkon and D'Angiulli, 2009. Bioscience Hypotheses). Namely, the
retina transforms external photon signals into electrical signals that are
carried to the V1 (striate cortex). Then, V1 retinotopic electrical signals
(spike-related electrical signals along classical axonal-dendritic pathways)
can be converted into regulated ultraweak bioluminescent photons (biophotons)
through redox processes within retinotopic visual neurons that make it possible
to create intrinsic biophysical pictures during visual perception and imagery.
However, the consensus opinion is to consider biophotons as by-products of
cellular metabolism. This paper argues that biophotons are not by-products,
other than originating from regulated cellular radical/redox processes. It also
shows that the biophoton intensity can be considerably higher inside cells than
outside. Our simple calculations, within a level of accuracy, suggest that the
real biophoton intensity in retinotopic neurons may be sufficient for creating
intrinsic biophysical picture representation of a single-object image during
visual perception.


  Complex networks can often be decomposed into less complex sub-networks whose
structures can give hints about the functional organization of the network as a
whole. However, these structural motifs can only tell one part of the
functional story because in this analysis each node and edge is treated on an
equal footing. In real networks, two motifs that are topologically identical
but whose nodes perform very different functions will play very different roles
in the network. Here, we combine structural information derived from the
topology of the neuronal network of the nematode C. elegans with information
about the biological function of these nodes, thus coloring nodes by function.
We discover that particular colorations of motifs are significantly more
abundant in the worm brain than expected by chance, and have particular
computational functions that emphasize the feed-forward structure of
information processing in the network, while evading feedback loops.
Interneurons are strongly over-represented among the common motifs, supporting
the notion that these motifs process and transduce the information from the
sensor neurons towards the muscles. Some of the most common motifs identified
in the search for significant colored motifs play a crucial role in the system
of neurons controlling the worm's locomotion. The analysis of complex networks
in terms of colored motifs combines two independent data sets to generate
insight about these networks that cannot be obtained with either data set
alone. The method is general and should allow a decomposition of any complex
network into its functional (rather than topological) motifs as long as both
wiring and functional information is available.


  The design of protocols to suppress the propagation of viral infections is an
enduring enterprise, especially hindered by limited knowledge of the mechanisms
through which extinction of infection propagation comes about. We here report
on a mechanism causing extinction of a propagating infection due to
intraspecific competition to infect susceptible hosts. Beneficial mutations
allow the pathogen to increase the production of progeny, while the host cell
is allowed to develop defenses against infection. When the number of
susceptible cells is unlimited, a feedback runaway co-evolution between host
resistance and progeny production occurs. However, physical space limits the
advantage that the virus can obtain from increasing offspring numbers, thus
infection clearance may result from an increase in host defenses beyond a
finite threshold. Our results might be relevant to better understand
propagation of viral infections in tissues with mobility constraints, and the
implications that environments with different geometrical properties might have
in devising control strategies.


  Scaling mobility patterns have been widely observed for animals. In this
paper, we propose a deterministic walk model to understand the scaling mobility
patterns, where walkers take the least-action walks on a lattice landscape and
prey. Scaling laws in the displacement distribution emerge when the amount of
prey resource approaches the critical point. Around the critical point, our
model generates ordered collective movements of walkers with a quasi-periodic
synchronization of walkers' directions. These results indicate that the
co-evolution of walkers' least-action behavior and the landscape could be a
potential origin of not only the individual scaling mobility patterns, but also
the flocks of animals. Our findings provide a bridge to connect the individual
scaling mobility patterns and the ordered collective movements.


  In genetics the Moran model describes the neutral evolution of a bi-allelic
gene in a population of haploid individuals subjected to mutations. We show in
this paper that this model can be mapped into an influence dynamical process on
networks subjected to external influences. The panmictic case considered by
Moran corresponds to fully connected networks and can be completely solved in
terms of hypergeometric functions. Other types of networks correspond to
structured populations, for which approximate solutions are also available.
This new approach to the classic Moran model leads to a relation between
regular networks based on spatial grids and the mechanism of isolation by
distance. We discuss the consequences of this connection for topopatric
speciation and the theory of neutral speciation and biodiversity. We show that
the effect of mutations in structured populations, where individuals can mate
only with neighbors, is greatly enhanced with respect to the panmictic case. If
mating is further constrained by genetic proximity between individuals, a
balance of opposing tendencies take place: increasing diversity promoted by
enhanced effective mutations versus decreasing diversity promoted by similarity
between mates. Stabilization occurs with speciation via pattern formation. We
derive an explicit relation involving the parameters characterizing the
population that indicates when speciation is possible.


  We study the spatio-temporal dynamics of the generalized repressilator, a
system of coupled repressing genes arranged in a directed ring topology, and
give analytical conditions for the emergence of a cascade of unstable periodic
orbits (UPOs) that lead to reachable long-lived oscillating transients. Such
transients dominate the finite time horizon dynamics that is relevant in
confined, noisy environments such as bacterial cells (see our previous work
[Strelkowa and Barahona, 2010]) and are therefore of interest for
bioengineering and synthetic biology. We show that the family of unstable
orbits possesses spatial symmetries and can also be understood in terms of
traveling wave solutions of kink-like topological defects. The long-lived
oscillatory transients correspond to the propagation of quasistable two-kink
configurations that unravel over a long time. We also assess the similarities
between the generalized repressilator model and other unidirectionally coupled
electronic systems, such as magnetic flux gates, which have been implemented
experimentally.


  Comparison of human and chimpanzee genomes has received much attention,
because of paramount role for understanding evolutionary step distinguishing us
from our closest living relative. In order to contribute to insight into Y
chromosome evolutionary history, we study and compare tandems, higher order
repeats (HORs), and regularly dispersed repeats in human and chimpanzee Y
chromosome contigs, using robust Global Repeat Map algorithm. We find a new
type of long-range acceleration, human-accelerated HOR regions. In peripheral
domains of 35mer human alphoid HORs, we find riddled features with ten
additional repeat monomers. In chimpanzee, we identify 30mer alphoid HOR. We
construct alphoid HOR schemes showing significant human-chimpanzee difference,
revealing rapid evolution after human-chimpanzee separation. We identify and
analyze over 20 large repeat units, most of them reported here for the first
time as: chimpanzee and human ~1.6 kb 3mer secondary repeat unit (SRU) and
~23.5 kb tertiary repeat unit (~0.55 kb primary repeat unit, PRU); human 10848,
15775, 20309, 60910, and 72140 bp PRUs; human 3mer SRU (~2.4 kb PRU); 715mer
and 1123mer SRUs (5mer PRU); chimpanzee 5096, 10762, 10853, 60523 bp PRUs; and
chimpanzee 64624 bp SRU (10853 bp PRU). We show that substantial
human-chimpanzee differences are concentrated in large repeat structures, at
the level of as much as ~70% divergence, sizably exceeding previous numerical
estimates for some selected noncoding sequences. Smeared over the whole
sequenced assembly (25 Mb) this gives ~14% human--chimpanzee divergence. This
is significantly higher estimate of divergence between human and chimpanzee
than previous estimates.


  Many cell functions are accomplished thanks to intracellular transport
mechanisms of macromolecules along filaments. Molecular motors such as dynein
or kinesin are proteins playing a primary role in these processes. The behavior
of such proteins is quite well understood when there is only one of them moving
a cargo particle. Indeed, numerous in vitro experiments have been performed to
derive accurate models for a single molecular motor. However, in vivo
macromolecules are often carried by multiple motors. The main focus of this
paper is to provide an analysis of the behavior of more molecular motors
interacting together in order to improve the understanding of their actual
physiological behavior. Previous studies provide analyses based on results
obtained from Monte Carlo simulations. Different from these studies, we derive
an equipollent probabilistic model to describe the dynamics of multiple
proteins coupled together and provide an exact theoretical analysis. We are
capable of obtaining the probability density function of the motor protein
configurations, thus enabling a deeper understanding of their behavior.


  [Background] Several studies have mentioned network modularity -- that a
network can easily be decomposed into subgraphs that are densely connected
within and weakly connected between each other -- as a factor affecting
metabolic robustness. In this paper we measure the relation between network
modularity and several aspects of robustness directly in a model system of
metabolism. [Methodology/Principal Findings] By using a model for generating
chemical reaction systems where one can tune the network modularity, we find
that robustness increases with modularity for changes in the concentrations of
metabolites, whereas it decreases with changes in the expression of enzymes.
The same modularity scaling is true for the speed of relaxation after the
perturbations. [Conclusions/Significance] Modularity is not a general principle
for making metabolism either more or less robust; this question needs to be
addressed specifically for different types of perturbations of the system.


  Protein footprinting is a new methodology that is based on probing, typically
with the use of mass spectrometry, of reactivity of different aminoacid
residues to a modifying reagent. Data thus obtained allow one to make
inferences about protein conformations and their intermolecular interactions.
Most of the protein footprinting studies so far have been performed on
individual proteins in vitro. We explore whether a similar approach is possible
with the proteins inside of living cells, employing dimethylsulfate (DMS), a
reagent widely used for the in vivo footprinting of nucleic acids. DMS can
induce methylation of the lysine, histidine and glutamate residues on proteins.
Using models of the histone H2B/H2AZ heterodimer assembled in vitro and from
chromatin treated in vivo, we show that the methylation by deuterated DMS
allows one to distinguish the accessibility of a particular residue in and out
of the protein's environmental/structural context. The detection of changes in
protein conformations or their interactions in vivo can provide a new approach
to the identification of proteins involved in various intracellular pathways
and help in the search for perspective drug targets and biomarkers of diseases.


  Neurotransmitter receptor molecules, concentrated in postsynaptic domains
along with scaffold and a number of other molecules, are key regulators of
signal transmission across synapses. Employing experiment and theory, we
develop a quantitative description of synaptic receptor domains in terms of a
reaction-diffusion model. We show that interactions between only receptor and
scaffold molecules, together with the rapid diffusion of receptors on the cell
membrane, are sufficient for the formation and stable characteristic size of
synaptic receptor domains. Our work reconciles long-term stability of synaptic
receptor domains with rapid turnover and diffusion of individual receptors.


  Not all nodes in a network are created equal. Differences and similarities
exist at both individual node and group levels. Disentangling single node from
group properties is crucial for network modeling and structural inference.
Based on unbiased generative probabilistic exponential random graph models and
employing distributive message passing techniques, we present an efficient
algorithm that allows one to separate the contributions of individual nodes and
groups of nodes to the network structure. This leads to improved detection
accuracy of latent class structure in real world data sets compared to models
that focus on group structure alone. Furthermore, the inclusion of hitherto
neglected group specific effects in models used to assess the statistical
significance of small subgraph (motif) distributions in networks may be
sufficient to explain most of the observed statistics. We show the predictive
power of such generative models in forecasting putative gene-disease
associations in the Online Mendelian Inheritance in Man (OMIM) database. The
approach is suitable for both directed and undirected uni-partite as well as
for bipartite networks.


  The catalytic conversion ATP + AMP -> 2ADP by the enzyme adenylate kinase
(ADK) involves the binding of one ATP molecule to the LID domain and one AMP
molecule to the NMP domain. The latter is followed by a phosphate transfer, and
then the release of two ADP molecules. We have computed a novel two dimensional
configurational free energy surface (2DCFES), with one reaction coordinate each
for the LID and the NMP domain motions, with explicit interactions with water.
Our computed 2DCFES clearly reveals the existence of a stable
half-open-half-closed (HOHC) intermediate state of the enzyme. Cycling of the
enzyme through the HOHC state reduces the conformational free energy barrier
for the reaction by about 20 kJ/mol. We find that the stability of the
half-open-half-closed state (missed in all earlier studies with implicit
solvent model) is largely because of the increase of specific interactions of
the polar amino acid side chains with water, particularly with the arginine and
the histidine residues. Free energy surface of the LID domain is rather rugged,
which can conveniently slow down LID's conformational motion, thus facilitating
a new substrate capture after the product release in the catalytic cycle.


  Genome-wide association studies have become increasingly common due to
advances in technology and have permitted the identification of differences in
single nucleotide polymorphism (SNP) alleles that are associated with diseases.
However, while typical GWAS analysis techniques treat markers individually,
complex diseases are unlikely to have a single causative gene. There is thus a
pressing need for multi-SNP analysis methods that can reveal system-level
differences in cases and controls. Here, we present a novel multi-SNP GWAS
analysis method called Pathways of Distinction Analysis (PoDA). The method uses
GWAS data and known pathway-gene and gene-SNP associations to identify pathways
that permit, ideally, the distinction of cases from controls. The technique is
based upon the hypothesis that if a pathway is related to disease risk, cases
will appear more similar to other cases than to controls for the SNPs
associated with that pathway. By systematically applying the method to all
pathways of potential interest, we can identify those for which the hypothesis
holds true, i.e., pathways containing SNPs for which the samples exhibit
greater within-class similarity than across classes. Importantly, PoDA improves
on existing single-SNP and SNP-set enrichment analyses in that it does not
require the SNPs in a pathway to exhibit independent main effects. This permits
PoDA to reveal pathways in which epistatic interactions drives risk. In this
paper, we detail the PoDA method and apply it to two GWA studies: one of breast
cancer, and the other of liver cancer. The results obtained strongly suggest
that there exist pathway-wide genomic differences that contribute to disease
susceptibility. PoDA thus provides an analytical tool that is complementary to
existing techniques and has the power to enrich our understanding of disease
genomics at the systems-level.


  The role of thermal pressure fluctuation excited within tightly packaged DNA
prior to ejection from protein capsid shells is discussed in a model
calculation. At equilibrium before ejection we assume the DNA is folded many
times into a bundle of parallel segments that forms an equilibrium conformation
at minimum free energy, which presses tightly against internal capsid walls.
Using a canonical ensemble at temperature T we calculate internal pressure
fluctuations against a slowly moving or static capsid mantle for an elastic
continuum model of the folded DNA bundle. It is found that fluctuating pressure
on the capsid internal wall from thermal excitation of longitudinal acoustic
vibrations in the bundle may have root-mean-square values which are several
tens of atmospheres for typically small phage dimensions. Comparisons are given
with measured data on three mutants of lambda phage with different base pair
lengths and total genome ejection pressures.


  Detection of false-positive motifs is one of the main causes of low
performance in motif finding methods. It is generally assumed that
false-positives are mostly due to algorithmic weakness of motif-finders. Here,
however, we derive the theoretical dependence of false positives on dataset
size and find that false positives can arise as a result of large dataset size,
irrespective of the algorithm used. Interestingly, the false-positive strength
depends more on the number of sequences in the dataset than it does on the
sequence length. As expected, false-positives can be reduced by decreasing the
sequence length or by adding more sequences to the dataset. The dependence on
number of sequences, however, diminishes and reaches a plateau after which
adding more sequences to the dataset does not reduce the false-positive rate
significantly. Based on the theoretical results presented here, we provide a
number of intuitive rules of thumb that may be used to enhance motif-finding
results in practice.


  In this work, we modify the superparamagnetic clustering algorithm (SPC) by
adding an extra weight to the interaction formula that considers which genes
are regulated by the same transcription factor. With this modified algorithm
that we call SPCTF, we analyze Spellman et al. microarray data for cell cycle
genes in yeast, and find clusters with a higher number of elements compared
with those obtained with the SPC algorithm. Some of the incorporated genes by
using SPCFT were not detected at first by Spellman et al. but were later
identified by other studies, whereas several genes still remain unclassified.
The clusters composed by unidentified genes were analyzed with MUSA, the motif
finding using an unsupervised approach algorithm, and this allow us to select
the clusters whose elements contain cell cycle transcription factor binding
sites as clusters worth of further experimental studies because they would
probably lead to new cell cycle genes. Finally, our idea of introducing
available information about transcription factors to optimize the gene
classification could be implemented for other distance-based clustering
algorithms


  It is known that the Kimura 3ST model of sequence evolution on phylogenetic
trees can be extended quite naturally to arbitrary split systems. However, this
extension relies heavily on mathematical peculiarities of the K3ST model, and
providing an analogous augmentation of the general Markov model has thus far
been elusive. In this paper we rectify this shortcoming by showing how to
extend the general Markov model on trees to to include arbitrary splits; and
even further to more general network models. This is achieved by exploring the
algebra of the generators of the continuous-time Markov chain together with the
"splitting" operator that generates the branching process on phylogenetic
trees. For simplicity we proceed by discussing the two state case and note that
our results are easily extended to more states with little complication.
Intriguingly, upon restriction of the two state general Markov model to the
parameter space of the binary symmetric model, our extension is
indistinguishable from the previous approach only on trees; as soon as any
incompatible splits are introduced the two approaches give rise to differing
probability distributions with disparate structure. Through exploration of a
simple example, we give a tentative argument that our approach to extending to
more general networks has desirable properties that the previous approaches do
not share. In particular, our construction allows for the possibility of
convergent evolution of previously divergent lineages; a property that is of
significant interest for biological applications.


  X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy are
two powerful tools to determine the protein 3D structure. However, not all
proteins can be successfully crystallized, particularly for membrane proteins.
Although NMR spectroscopy is indeed very powerful in determining the 3D
structures of membrane proteins, same as X-ray crystallography, it is still
very time-consuming and expensive. Under many circumstances, due to the
noncrystalline and insoluble nature of some proteins, X-ray and NMR cannot be
used at all. Computational approaches, however, allow us to obtain a
description of the protein 3D structure at a submicroscopic level.
  To the best of the author's knowledge, there is little structural data
available to date on the AGAAAAGA palindrome in the hydrophobic region
(113--120) of prion proteins, which falls just within the N-terminal
unstructured region (1--123) of prion proteins. Many experimental studies have
shown that the AGAAAAGA region has amyloid fibril forming properties and plays
an important role in prion diseases. However, due to the noncrystalline and
insoluble nature of the amyloid fibril, little structural data on the AGAAAAGA
is available. This paper introduces a simple optimization strategy approach to
address the 3D atomic-resolution structure of prion AGAAAAGA amyloid fibrils.
Atomic-resolution structures of prion AGAAAAGA amyloid fibrils got in this
paper are useful for the drive to find treatments for prion diseases in the
field of medicinal chemistry.


  We consider a fixed size population that undergoes an evolutionary adaptation
in the weak mutuation rate limit, which we model as a biased Langevin process
in the genotype space. We show analytically and numerically that, if the
fitness landscape has a small highly epistatic (rough) and time-varying
component, then the population genotype exhibits a high effective diffusion in
the genotype space and is able to escape local fitness minima with a large
probability. We argue that our principal finding that even very small
time-dependent fluctuations of fitness can substantially speed up evolution is
valid for a wide class of models.