  The ubiquity of modules in biological networks may result from an
evolutionary benefit of a modular organization. For instance, modularity may
increase the rate of adaptive evolution, because modules can be easily combined
into new arrangements that may benefit their carrier. Conversely, modularity
may emerge as a by-product of some trait. We here ask whether this last
scenario may play a role in genome-scale metabolic networks that need to
sustain life in one or more chemical environments. For such networks, we define
a network module as a maximal set of reactions that are fully coupled, i.e.,
whose fluxes can only vary in fixed proportions. This definition overcomes
limitations of purely graph based analyses of metabolism by exploiting the
functional links between reactions. We call a metabolic network viable in a
given chemical environment if it can synthesize all of an organism's biomass
compounds from nutrients in this environment. An organism's metabolism is
highly versatile if it can sustain life in many different chemical
environments. We here ask whether versatility affects the modularity of
metabolic networks.


  The internal dynamics of strongly interacting systems and that of
biomolecules such as proteins display several important analogies, despite the
huge difference in their characteristic energy and length scales. For example,
in all such systems, collective excitations, cooperative transitions and phase
transitions emerge as the result of the interplay of strong correlations with
quantum or thermal fluctuations. In view of such an observation, some
theoretical methods initially developed in the context of theoretical nuclear
physics have been adapted to investigate the dynamics of biomolecules. In this
talk, we review some of our recent studies performed along this direction. In
particular, we discuss how the path integral formulation of the molecular
dynamics allows to overcome some of the long-standing problems and limitations
which emerge when simulating the protein folding dynamics at the atomistic
level of detail.


  Heat-bath cooling is a component of practicable algorithmic cooling of spins,
an approach which might be useful for in vivo 13C spectroscopy, in particular
for prolonged metabolic processes where substrates that are hyperpolarized
ex-vivo are not effective. We applied heat-bath cooling to 1,2-13C2-amino
acids, using the alpha protons to shift entropy from selected carbons to the
environment. For glutamate and glycine, both carbons were cooled by about
2.5-fold, and in other experiments the polarization of C1 nearly doubled while
all other spins had equilibrium polarization, indicating reduction in total
entropy. The effect of adding Magnevist, a gadolinium contrast agent, on
heat-bath cooling of glutamate was investigated.


  In this paper we find and classify all patterns for a single locus three- and
four-allele population genetics models in continuous time. A pattern for a
$k$-allele model means all coexisting locally stable equilibria with respect to
the flow defined by the equations $\dot{p}_i = p_i(r_i-r), i=1,...,k,$ where
$p_i, r_i$ are the frequency and marginal fitness of allele $A_i$,
respectively, and $r$ is the mean fitness of the population. It is well known
that for the two-allele model there are only three patterns depending on the
relative fitness between the homozygotes and the heterozygote. It turns out
that for the three-allele model there are 14 patterns and for the four-allele
model there are 117 patterns. With the help of computer simulations, we find
2351 patterns for the five-allele model. For the six-allele model, there are
more than 60,000 patterns. In addition, for each pattern of the three-allele
model, we also determine the asymptotic behavior of solutions of the above
system of equations as $t \to \infty$. The problem of finding patterns has been
studied in the past and it is an important problem because the results can be
used to predict the long-term genetic makeup of a population.


  In a recent paper on 'Estimating Species Trees from Unrooted Gene Trees' Liu
and Yu observe that the distance matrix on the underlying taxon set, which is
built up from expected internode distances on gene trees under the multispecies
coalescent, is tree-like, and that the underlying additive tree has the same
topology as the true species tree. Hence they suggest to use (observed) average
internode distances on gene trees as an input for the neighbor joining
algorithm to estimate the underlying species tree in a statistically consistent
way. In this note we give a rigorous proof of their above mentioned
observation.


  Extracting meaningful patterns from voluminous amount of biological data is a
very big challenge. Motifs are biological patterns of great interest to
biologists. Many different versions of the motif finding problem have been
identified by researchers. Examples include the Planted $(l, d)$ Motif version,
those based on position-specific score matrices, etc. A comparative study of
the various motif search algorithms is very important for several reasons. For
example, we could identify the strengths and weaknesses of each. As a result,
we might be able to devise hybrids that will perform better than the individual
components. In this paper we (either directly or indirectly) compare the
performance of PMSprune (an algorithm based on the $(l, d)$ motif model) and
several other algorithms in terms of seven measures and using well established
benchmarks
  In this paper, we (directly or indirectly) compare the quality of motifs
predicted by PMSprune and 14 other algorithms. We have employed several
benchmark datasets including the one used by Tompa, et.al. These comparisons
show that the performance of PMSprune is competitive when compared to the other
14 algorithms tested.
  We have compared (directly or indirectly) the performance of PMSprune and 14
other algorithms using the Benchmark dataset provided by Tompa, et.al. It is
observed that both PMSprune and DME (an algorithm based on position-specific
score matrices) in general perform better than the 13 algorithms reported in
Tompa et. al.. Subsequently we have compared PMSprune and DME on other
benchmark data sets including ChIP-Chip, ChIP-seq, and ABS. Between PMSprune
and DME, PMSprune performs better than DME on six measures. DME performs better
than PMSprune on one measure (namely, specificity).


  Chemical reaction systems are dynamical systems that arise in chemical
engineering and systems biology. In this work, we consider the question of
whether the minimal (in a precise sense) multistationary chemical reaction
networks, which we propose to call `atoms of multistationarity,' characterize
the entire set of multistationary networks. Our main result states that the
answer to this question is `yes' in the context of fully open continuous-flow
stirred-tank reactors (CFSTRs), which are networks in which all chemical
species take part in the inflow and outflow. In order to prove this result, we
show that if a subnetwork admits multiple steady states, then these steady
states can be lifted to a larger network, provided that the two networks share
the same stoichiometric subspace. We also prove an analogous result when a
smaller network is obtained from a larger network by `removing species.' Our
results provide the mathematical foundation for a technique used by
Siegal-Gaskins et al. of establishing bistability by way of `network ancestry.'
Additionally, our work provides sufficient conditions for establishing
multistationarity by way of atoms and moreover reduces the problem of
classifying multistationary CFSTRs to that of cataloging atoms of
multistationarity. As an application, we enumerate and classify all 386
bimolecular and reversible two-reaction networks. Of these, exactly 35 admit
multiple positive steady states. Moreover, each admits a unique minimal
multistationary subnetwork, and these subnetworks form a poset (with respect to
the relation of `removing species') which has 11 minimal elements (the atoms of
multistationarity).


  We introduce an implicit solvent Molecular Dynamics approach for calculating
ionic fluxes through narrow nano-pores and transmembrane channels. The method
relies on a dual-control- volume grand-canonical molecular dynamics (DCV-GCMD)
simulation and the analytical solution for the electrostatic potential inside a
cylindrical nano-pore recently obtained by Levin [Europhys. Lett., 76, 163
(2006)]. The theory is used to calculate the ionic fluxes through an artificial
trans-membrane c hannel which mimics the antibacterial gramicidin A channel.
Both current-voltage and current-concentration relations are calculated under
various experimental conditions. We show that our results are comparable to the
characteristics associated to the gramicidin A pore, specially the existence of
two binding sites inside the pore and the observed saturation in the
current-concentration profiles.


  In this paper, we show that the coarse grain model for DNA, which has been
proposed recently by Knotts, Rathore, Schwartz and de Pablo (J. Chem. Phys.
126, 084901 (2007)), can be adapted to describe the thermal and mechanical
denaturation of long DNA sequences by adjusting slightly the base pairing
contribution. The adjusted model leads to (i) critical temperatures for long
homogeneous sequences that are in good agreement with both experimental ones
and those obtained from statistical models, (ii) a realistic step-like
denaturation behaviour for long inhomogeneous sequences, and (iii) critical
forces at ambient temperature of the order of 10 pN, close to measured values.
The adjusted model furthermore supports the conclusion that the thermal
denaturation of long homogeneous sequences corresponds to a first-order phase
transition and yields a critical exponent for the critical force equal to
sigma=0.70. This model is both geometrically and energetically realistic, in
the sense that the helical structure and the grooves, where most proteins bind,
are satisfactorily reproduced, while the energy and force required to break a
base pair lie in the expected range. It therefore represents a promising tool
for studying the dynamics of DNA-protein specific interactions at an
unprecedented detail level.


  In this paper we suggest that, under suitable conditions, supervised learning
can provide the basis to formulate at the microscopic level quantitative
questions on the phenotype structure of multicellular organisms. The problem of
explaining the robustness of the phenotype structure is rephrased as a real
geometrical problem on a fixed domain. We further suggest a generalization of
path integrals that reduces the problem of deciding whether a given molecular
network can generate specific phenotypes to a numerical property of a
robustness function with complex output, for which we give heuristic
justification. Finally, we use our formalism to interpret a pointedly
quantitative developmental biology problem on the allowed number of pairs of
legs in centipedes.


  Human language is still an embarrassment for evolutionary theory, as the
speaker's benefit remains unclear. The willingness to communicate information
is shown here to be an evolutionary stable strategy (ESS), even if acquiring
original information from the environment involves significant cost and
communicating it provides no material benefit to addressees. In this study,
communication is used to advertise the emitter's ability to obtain novel
information. We found that communication strategies can take two forms,
competitive and uniform, that these two strategies are stable and that they
necessarily coexist.


  The common techniques to study protein-protein proximity in vivo are not
well-adapted to the capabilities and the expertise of a standard proteomics
laboratory, typically based on the use of mass spectrometry. With the aim of
closing this gap, we have developed PUB-MS (for Proximity Utilizing
Biotinylation and Mass Spectrometry), an approach to monitor protein-protein
proximity, based on biotinylation of a protein fused to a biotin-acceptor
peptide (BAP) by a biotin-ligase, BirA, fused to its interaction partner. The
biotinylation status of the BAP can be further detected by either Western
analysis or mass spectrometry. The BAP sequence was redesigned for easy
monitoring of the biotinylation status by LC-MS/MS. In several experimental
models, we demonstrate that the biotinylation in vivo is specifically enhanced
when the BAP- and BirA- fused proteins are in proximity to each other. The
advantage of mass spectrometry is demonstrated by using BAPs with different
sequences in a single experiment (allowing multiplex analysis) and by the use
of stable isotopes. Finally, we show that our methodology can be also used to
study a specific subfraction of a protein of interest that was in proximity
with another protein at a predefined time before the analysis.


  We show that in experimental atomic force microscopy studies of the lifetime
distribution of mechanically stressed folded proteins the effects of externally
applied fluctuations can not be distinguished from those of internally present
fluctuations. In certain circumstances this leads to artificially
non-exponential lifetime distributions which can be misinterpreted as a
signature of protein complexity. This work highlights the importance of fully
characterizing and controlling external sources of fluctuation in mechanical
studies of proteins before drawing conclusions on the physics at play on the
molecular level.


  We introduce a new distance-based phylogeny reconstruction technique which
provably achieves, at sufficiently short branch lengths, a logarithmic
sequence-length requirement---improving significantly over previous polynomial
bounds for distance-based methods and matching existing results for general
methods. The technique is based on an averaging procedure that implicitly
reconstructs ancestral sequences.
  In the same token, we extend previous results on phase transitions in
phylogeny reconstruction to general time-reversible models. More precisely, we
show that in the so-called Kesten-Stigum zone (roughly, a region of the
parameter space where ancestral sequences are well approximated by "linear
combinations" of the observed sequences) sequences of length $O(\log n)$
suffice for reconstruction when branch lengths are discretized. Here $n$ is the
number of extant species.
  Our results challenge, to some extent, the conventional wisdom that estimates
of evolutionary distances alone carry significantly less information about
phylogenies than full sequence datasets.


  We summarize some aspects of electrostatic interactions in the context of
viruses. A simplified but, within well defined limitations, reliable approach
is used to derive expressions for electrostatic energies and the corresponding
osmotic pressures in single-stranded RNA viruses and double-stranded DNA
bacteriophages. The two types of viruses differ crucially in the spatial
distribution of their genome charge which leads to essential differences in
their free energies, depending on the capsid size and total charge in a quite
different fashion. Differences in the free energies are trailed by the
corresponding characteristics and variations in the osmotic pressure between
the inside of the virus and the external bathing solution.


  Certain membrane proteins, peptides, nanoparticles and nanotubes have rigid
structure and fixed shape. They are often viewed as spheres and cylinders with
certain surface properties. Single Chain Mean Field theory is used to model the
equilibrium insertion of nanoscale spheres and rods into the phospholipid
bilayer. The equilibrium structures and the resulting free energies of the
nano-objects in the bilayer allow to distinguish different orientations in the
bilayer and estimate the energy barrier of insertion.


  We study analytically and numerically a generic continuum model of an
isotropic active solid with internal stresses generated by non-equilibrium
`active' mechano-chemical reactions. Our analysis shows that the gel can be
tuned through three classes of dynamical states by increasing motor activity: a
constant unstrained state of homogeneous density, a state where the local
density exhibits sustained oscillations, and a steady-state which is
spontaneously contracted, with a uniform mean density.


  The complex interactions involved in regulation of a cell's function are
captured by its interaction graph. More often than not, detailed knowledge
about enhancing or suppressive regulatory influences and cooperative effects is
lacking and merely the presence or absence of directed interactions is known.
Here we investigate to which extent such reduced information allows to forecast
the effect of a knock-out or a combination of knock-outs. Specifically we ask
in how far the lethality of eliminating nodes may be predicted by their network
centrality, such as degree and betweenness, without knowing the function of the
system. The function is taken as the ability to reproduce a fixed point under a
discrete Boolean dynamics. We investigate two types of stochastically generated
networks: fully random networks and structures grown with a mechanism of node
duplication and subsequent divergence of interactions. On all networks we find
that the out-degree is a good predictor of the lethality of a single node
knock-out. For knock-outs of node pairs, the fraction of successors shared
between the two knocked-out nodes (out-overlap) is a good predictor of
synthetic lethality. Out-degree and out-overlap are locally defined and
computationally simple centrality measures that provide a predictive power
close to the optimal predictor.


  It is well-known that biological and social interaction networks have a
varying degree of redundancy, though a consensus of the precise cause of this
is so far lacking. In this paper, we introduce a topological redundancy measure
for labeled directed networks that is formal, computationally efficient and
applicable to a variety of directed networks such as cellular signaling,
metabolic and social interaction networks. We demonstrate the computational
efficiency of our measure by computing its value and statistical significance
on a number of biological and social networks with up to several thousands of
nodes and edges. Our results suggest a number of interesting observations: (1)
social networks are more redundant that their biological counterparts, (2)
transcriptional networks are less redundant than signaling networks, (3) the
topological redundancy of the C. elegans metabolic network is largely due to
its inclusion of currency metabolites, and (4) the redundancy of signaling
networks is highly (negatively) correlated with the monotonicity of their
dynamics.


  We present the first application of the hypothesis-rich mathematical theory
to genome-wide association data. The Hamza et al. late-onset sporadic
Parkinson's disease genome-wide association study dataset was analyzed. We
found a rare, coding, non-synonymous SNP variant in the gene DZIP1 that confers
increased susceptibility to Parkinson's disease. The association of DZIP1 with
Parkinson's disease is consistent with a Parkinson's disease stem-cell ageing
theory.


  We use a multitype continuous time Markov branching process model to describe
the dynamics of the spread of parasites of two types that can mutate into each
other in a common host population. Instead of using a single virulence
characteristic which is typical of most mathematical models for infectious
diseases, our model uses a combination of two characteristics: lethality and
transmissibility. This makes the model capable of reproducing the empirically
observed fact that the increase in the host density can lead to the prevalence
of the more virulent pathogen type. We provide some numerical illustrations and
discuss the effects of the size of the enclosure containing the host population
on the encounter rate in our model that plays the key role in determining what
pathogen type will eventually prevail. We also present a multistage extension
of the model to situations where there are several populations and parasites
can be transmitted from one of them to another.


  Highlighting similarities and differences between networks is an informative
task in investigating many biological processes. Typical examples are detecting
differences between an inferred network and the corresponding gold standard, or
evaluating changes in a dynamic network along time. Although fruitful insights
can be drawn by qualitative or feature-based methods, a distance must be used
whenever a quantitative assessment is required. Here we introduce the
Ipsen-Mikhailov metric for biological network comparison, based on the
difference of the distributions of the Laplacian eigenvalues of the compared
graphs. Being a spectral measure, its focus is on the general structure of the
net so it can overcome the issues affecting local metrics such as the edit
distances. Relation with the classical Matthews Correlation Coefficient (MCC)
is discussed, showing the finer discriminant resolution achieved by the
Ipsen-Mikhailov metric. We conclude with three examples of application in
functional genomic tasks, including stability of network reconstruction as
robustness to data subsampling, variability in dynamical networks and
differences in networks associated to a classification task.


  We construct a descriptive toy model that considers quantum effects on
biological evolution starting from Chaitin's classical framework. There are
smart evolution scenarios in which a quantum world is as favorable as classical
worlds for evolution to take place. However, in more natural scenarios, the
rate of evolution depends on the degree of entanglement present in quantum
organisms with respect to classical organisms. If the entanglement is maximal,
classical evolution turns out to be more favorable.


  Both symmetry and organized breaking of symmetry have a pivotal r\^ole in our
understanding of structure and pattern formation in physical systems, including
the origin of mass in the Universe and the chiral structure of biological
macromolecules. Here we report on a new symmetry breaking phenomenon that takes
place in all biologically active proteins, thus this symmetry breaking relates
to the inception of life. The unbroken symmetry determines the covalent bond
geometry of a sp3 hybridized carbon atom. It dictates the tetrahedral
architecture of atoms around the central carbon of an amino acid. Here we show
that in a biologically active protein this symmetry becomes broken. Moreover,
we show that the pattern of symmetry breaking is in a direct correspondence
with the local secondary structure of the folded protein.


  We analyze the statistical properties of Poincar\'e recurrences of Homo
sapiens, mammalian and other DNA sequences taken from Ensembl Genome data base
with up to fifteen billions base pairs. We show that the probability of
Poincar\'e recurrences decays in an algebraic way with the Poincar\'e exponent
$\beta \approx 4$ even if oscillatory dependence is well pronounced. The
correlations between recurrences decay with an exponent $\nu \approx 0.6$ that
leads to an anomalous super-diffusive walk. However, for Homo sapiens
sequences, with the largest available statistics, the diffusion coefficient
converges to a finite value on distances larger than million base pairs. We
argue that the approach based on Poncar\'e recurrences determines new proximity
features between different species and shed a new light on their evolution
history.


  We develop and study the evolution of a model of radiation induced apoptosis
in cells using stochastic simulations, and identified key protein targets for
effective mitigation of radiation damage. We identified several key proteins
associated with cellular apoptosis using an extensive literature survey. In
particular, we focus on the p53 transcription dependent and p53 transcription
independent pathways for mitochondrial apoptosis. Our model reproduces known
p53 oscillations following radiation damage. The key, experimentally testable
hypotheses that we generate are - inhibition of PUMA is an effective strategy
for mitigation of radiation damage if the treatment is administered
immediately, at later stages following radiation damage, inhibition of tBid is
more effective.


  Consider a large Boolean network with a feed forward structure. Given a
probability distribution on the inputs, can one find, possibly small,
collections of input nodes that determine the states of most other nodes in the
network? To answer this question, a notion that quantifies the determinative
power of an input over the states of the nodes in the network is needed. We
argue that the mutual information (MI) between a given subset of the inputs X =
{X_1, ..., X_n} of some node i and its associated function f_i(X) quantifies
the determinative power of this set of inputs over node i. We compare the
determinative power of a set of inputs to the sensitivity to perturbations to
these inputs, and find that, maybe surprisingly, an input that has large
sensitivity to perturbations does not necessarily have large determinative
power. However, for unate functions, which play an important role in genetic
regulatory networks, we find a direct relation between MI and sensitivity to
perturbations. As an application of our results, we analyze the large-scale
regulatory network of Escherichia coli. We identify the most determinative
nodes and show that a small subset of those reduces the overall uncertainty of
the network state significantly. Furthermore, the network is found to be
tolerant to perturbations of its inputs.


  Motivation: High-throughput sequencing enables expression analysis at the
level of individual transcripts. The analysis of transcriptome expression
levels and differential expression estimation requires a probabilistic approach
to properly account for ambiguity caused by shared exons and finite read
sampling as well as the intrinsic biological variance of transcript expression.
  Results: We present BitSeq (Bayesian Inference of Transcripts from Sequencing
data), a Bayesian approach for estimation of transcript expression level from
RNA-seq experiments. Inferred relative expression is represented by Markov
chain Monte Carlo (MCMC) samples from the posterior probability distribution of
a generative model of the read data. We propose a novel method for differential
expression analysis across replicates which propagates uncertainty from the
sample-level model while modelling biological variance using an
expression-level-dependent prior. We demonstrate the advantages of our method
using simulated data as well as an RNA-seq dataset with technical and
biological replication for both studied conditions.
  Availability: The implementation of the transcriptome expression estimation
and differential expression analysis, BitSeq, has been written in C++.


  Assembly and disassembly dynamics of microtubules (MTs) is tightly controlled
by MT associated proteins. Here, we investigate how plus-end-directed
depolymerases of the kinesin-8 family regulate MT depolymerization dynamics.
Employing an individual-based model, we reproduce experimental findings.
Moreover, crowding is identified as the key regulatory mechanism of
depolymerization dynamics. Our analysis gives two qualitatively distinct
regimes. For motor densities above a particular threshold, a macroscopic
traffic jam emerges at the plus-end and the MT dynamics become independent of
the motor concentration. Below this threshold, microscopic traffic jams at the
tip arise which cancel out the effect of the depolymerization kinetics such
that the depolymerization speed is solely determined by the motor density.
Because this density changes over the MT length, length-dependent regulation is
possible. Remarkably, motor cooperativity does not affect the depolymerization
speed but only the end-residence time of depolymerases.


  Alzheimer's disease causes severe neurodegeneration in the brain that leads
to a certain death. The defining factor is the formation of extracellular
senile amyloid plaques in the brain. However, therapeutic approaches to remove
them have not been effective in humans, and so our understanding of the cause
of Alzheimer's disease remains incomplete. Here we investigate physical
processes that might relate to its onset. Instead of the extracellular amyloid,
we scrutinize the intracellular domain of its precursor protein. We argue for a
phenomenon that has never before been discussed in the context of polymer
physics: Like ice and water together, the intracellular domain of the amyloid
precursor protein forms a state of phase coexistence with another protein. This
leads to an inherent instability that could well be among the missing pieces in
the puzzle of Alzheimer's disease.


  We consider a population of genotype sequences evolving on a rugged fitness
landscape with many local fitness peaks. The population walks uphill until it
encounters a local fitness maximum. We find that the statistical properties of
the walk length depend on whether the underlying fitness distribution has a
finite mean. If the mean is finite, all the walk length cumulants grow with the
sequence length but approach a constant otherwise. Experimental implications of
our analytical results are also discussed.


  Verhulst logistic curve either grows OR decays, depending on the {\it growth
rate} parameter value. A similar situation is found in the Gompertz law about
human mortality. However, growth can neither be infinite nor reach a finite
steady state at an infinite asymptotic time. Moreover before some decay, some
growth must have occurred. For example, some levelling-off could occur at
finite time, followed either by some growth again or then by some decay.
Numerous examples show that the Verhulst and Gompertz modelisations are too
reductive (or restrictive) descriptions of their original purpose. It is aimed,
in the present note, to encompass into ONE simple differential equation the
growth AND decay features of, e.g., population sizes, or numbers, but also of
many other measured characteristics found in social and physical science
systems. Previous generalisations of Verhulst or Gompertz functions are
recalled. It is shown that drastic growth or decay jumps or turnovers can be
readily described through drastic changes in values of the growth or decay
rate. However smoother descriptions can be found if the growth or decay rate is
modified in order to take into account some time or size dependence. Similar
arguments can be carried through, but not so easily, for the so called carrying
capacity, indeed leading to more elaborate algebraic work.


  The design of genetic networks with specific functions is one of the major
goals of synthetic biology. However, constructing biological devices that work
"as required" remains challenging, while the cost of uncovering flawed designs
experimentally is large. To address this issue, we propose a fully automated
framework that allows the correctness of synthetic gene networks to be formally
verified in silico from rich, high level functional specifications.
  Given a device, we automatically construct a mathematical model from
experimental data characterizing the parts it is composed of. The specific
model structure guarantees that all experimental observations are captured and
allows us to construct finite abstractions through polyhedral operations. The
correctness of the model with respect to temporal logic specifications can then
be verified automatically using methods inspired by model checking.
  Overall, our procedure is conservative but it can filter through a large
number of potential device designs and select few that satisfy the
specification to be implemented and tested further experimentally. Illustrative
examples of the application of our methods to the design of simple synthetic
gene networks are included.


  The rapid growth of cancer genome structural information provides an
opportunity for a better understanding of the mutational mechanisms of genomic
alterations in cancer and the forces of selection that act upon them. Here we
test the evidence for two major forces, spatial chromosome structure and
purifying (or negative) selection, that shape the landscape of somatic
copy-number alterations (SCNAs) in cancer1. Using a maximum likelihood
framework we compare SCNA maps and three-dimensional genome architecture as
determined by genome-wide chromosome conformation capture (HiC) and described
by the proposed fractal-globule (FG) model2. This analysis provides evidence
that the distribution of chromosomal alterations in cancer is spatially related
to three-dimensional genomic architecture and additionally suggests that
purifying selection as well as positive selection shapes the landscape of SCNAs
during somatic evolution of cancer cells.


  Gene regulatory networks are collections of genes that interact with one
other and with other substances in the cell. By measuring gene expression over
time using high-throughput technologies, it may be possible to reverse
engineer, or infer, the structure of the gene network involved in a particular
cellular process. These gene expression data typically have a high
dimensionality and a limited number of biological replicates and time points.
Due to these issues and the complexity of biological systems, the problem of
reverse engineering networks from gene expression data demands a specialized
suite of statistical tools and methodologies. We propose a non-standard
adaptation of a simulation-based approach known as Approximate Bayesian
Computing based on Markov chain Monte Carlo sampling. This approach is
particularly well suited for the inference of gene regulatory networks from
longitudinal data. The performance of this approach is investigated via
simulations and using longitudinal expression data from a genetic repair system
in Escherichia coli.


  Natural ecosystems are characterized by striking diversity of form and
functions and yet exhibit deep symmetries emerging across scales of space, time
and organizational complexity. Species-area relationships and species-abundance
distributions are examples of emerging patterns irrespective of the details of
the underlying ecosystem functions. Here we present empirical and theoretical
evidence for a new macroecological pattern related to the distributions of
local species persistence times, defined as the timespans between local
colonizations and extinctions in a given geographic region. Empirical
distributions pertaining to two different taxa, breeding birds and herbaceous
plants, analyzed in a new framework that accounts for the finiteness of the
observational period, exhibit power-law scaling limited by a cut-off determined
by the rate of emergence of new species. In spite of the differences between
taxa and spatial scales of analysis, the scaling exponents are statistically
indistinguishable from each other and significantly different from those
predicted by existing models. We theoretically investigate how the scaling
features depend on the structure of the spatial interaction network and show
that the empirical scaling exponents are reproduced once a two-dimensional
isotropic texture is used, regardless of the details of the ecological
interactions. The framework developed here also allows to link the cut-off
timescale with the spatial scale of analysis, and the persistence-time
distribution to the species-area relationship. We conclude that the inherent
coherence obtained between spatial and temporal macroecological patterns points
at a seemingly general feature of the dynamical evolution of ecosystems.


  We present a model of an ATP-fueled molecular machine which push a polymer
through a pore channel. The machine acts between two levels (working-waiting),
and the working one remains active for a fixed time giving a constant force.
The machine activation rate can be put in relationship with the available ATP
concentration in the solution, which gives the necessary energy supply. The
translocation time shows a monotonic behavior as a function of the activation
frequency and the velocity follows a Michaelis-Menten law that arises naturally
in this description. The estimation of the stall force of the motor follows a
corrected Michaelis-Menten law which still is to be checked in experimental
investigation. The results presented agree with recent biological experimental
findings.


  Gene delivery of nucleic acid to the cell nucleus is a fundamental step in
gene therapy. In this review of modeling drug and gene delivery, we focus on
the particular stage of plasmid DNA or virus cytoplasmic trafficking. A
challenging problem is to quantify the success of this limiting stage. We
present some models and simulations of plasmid trafficking and of the limiting
phase of DNA-polycation escape from an endosome and discuss virus cytoplasmic
trafficking. The models can be used to assess the success of viral escape from
endosomes, to quantify the early step of viral-cell infection, and to propose
new simulation tools for designing new hybrid-viruses as synthetic vectors.


  We consider chemical reaction networks taken with mass action kinetics. The
steady states of such a system are solutions to a system of polynomial
equations. Even for small systems the task of finding the solutions is
daunting. We develop an algebraic framework and procedure for linear
elimination of variables. The procedure reduces the variables in the system to
a set of "core" variables by eliminating variables corresponding to a set of
non-interacting species. The steady states are parameterized algebraically by
the core variables, and a graphical condition is given for when a steady state
with positive core variables necessarily have all variables positive. Further,
we characterize graphically the sets of eliminated variables that are
constrained by a conservation law and show that this conservation law takes a
specific form.


  We consider the dynamics of Q-learning in two-player two-action games with a
Boltzmann exploration mechanism. For any non-zero exploration rate the dynamics
is dissipative, which guarantees that agent strategies converge to rest points
that are generally different from the game's Nash Equlibria (NE). We provide a
comprehensive characterization of the rest point structure for different games,
and examine the sensitivity of this structure with respect to the noise due to
exploration. Our results indicate that for a class of games with multiple NE
the asymptotic behavior of learning dynamics can undergo drastic changes at
critical exploration rates. Furthermore, we demonstrate that for certain games
with a single NE, it is possible to have additional rest points (not
corresponding to any NE) that persist for a finite range of the exploration
rates and disappear when the exploration rates of both players tend to zero.


  The success of imitation as an evolutionary driving force in spatial games
has often been questioned, especially for social dilemmas such as the snowdrift
game, where the most profitable may be the mixed phase sustaining both the
cooperative as well as the defective strategy. Here we reexamine this
assumption by investigating the evolution of cooperation in spatial social
dilemma games, where instead of pure strategies players can adopt emotional
profiles of their neighbors. For simplicity, the emotional profile of each
player is determined by two pivotal factors only, namely how it behaves towards
less and how towards more successful neighbors. We find that imitating emotions
such as goodwill and envy instead of pure strategies from the more successful
players reestablishes imitation as a tour de force for resolving social
dilemmas on structured populations without any additional assumptions or
strategic complexity.


  We introduce a model of biological evolution where species evolve in response
to biotic interactions and a fluctuating environmental stress. The species may
either become extinct or mutate to acquire a new fitness value when the
effective stress level is greater than their individual fitness. The model
exhibits a phase transition to a completely extinct phase as the environmental
stress or the mutation rate is varied. We discuss the generic conditions for
which this transition is continuous. The model is exactly solvable and the
critical behavior is characterized by an unusual dynamic exponent z=1/3. Apart
from predicting large scale evolution, the model can be applied to understand
the trends in the available fossil data.


  Topological network motifs represent functional relationships within and
between regulatory and protein-protein interaction networks. Enriched motifs
often aggregate into self-contained units forming functional modules.
Theoretical models for network evolution by duplication-divergence mechanisms
and for network topology by hierarchical scale-free networks have suggested a
one-to-one relation between network motif enrichment and aggregation, but this
relation has never been tested quantitatively in real biological interaction
networks. Here we introduce a novel method for assessing the statistical
significance of network motif aggregation and for identifying clusters of
overlapping network motifs. Using an integrated network of transcriptional,
posttranslational and protein-protein interactions in yeast we show that
network motif aggregation reflects a local modularity property which is
independent of network motif enrichment. In particular our method identified
novel functional network themes for a set of motifs which are not enriched yet
aggregate significantly and challenges the conventional view that network motif
enrichment is the most basic organizational principle of complex networks.


  Metabolism is a fascinating cell machinery underlying life and disease and
genome-scale reconstructions provide us with a captivating view of its
complexity. However, deciphering the relationship between metabolic structure
and function remains a major challenge. In particular, turning observed
structural regularities into organizing principles underlying systemic
functions is a crucial task that can be significantly addressed after endowing
complex network representations of metabolism with the notion of geometric
distance. Here, we design a cartographic map of metabolic networks by embedding
them into a simple geometry that provides a natural explanation for their
observed network topology and that codifies node proximity as a measure of
hidden structural similarities. We assume a simple and general connectivity law
that gives more probability of interaction to metabolite/reaction pairs which
are closer in the hidden space. Remarkably, we find an astonishing congruency
between the architecture of E. coli and human cell metabolisms and the
underlying geometry. In addition, the formalism unveils a backbone-like
structure of connected biochemical pathways on the basis of a quantitative
cross-talk. Pathways thus acquire a new perspective which challenges their
classical view as self-contained functional units.


  We analytically study the input-output properties of a neuron whose active
dendritic tree, modeled as a Cayley tree of excitable elements, is subjected to
Poisson stimulus. Both single-site and two-site mean-field approximations
incorrectly predict a non-equilibrium phase transition which is not allowed in
the model. We propose an excitable-wave mean-field approximation which shows
good agreement with previously published simulation results [Gollo et al., PLoS
Comput. Biol. 5(6) e1000402 (2009)] and accounts for finite-size effects. We
also discuss the relevance of our results to experiments in neuroscience,
emphasizing the role of active dendrites in the enhancement of dynamic range
and in gain control modulation.


  For decades, the hemagglutination-inhibition (HI) assay has been used in
epidemiological and basic biological studies of influenza viruses. The
mechanistic basis of the assay results (called titers) is not well understood.
The first part of this document describes a biophysical model of HI that
illuminates the mechanistic basis of and provides the theoretical motivation
for new ways of interpreting titers. The biophysical model is applicable to
other serological assays; this fact is illustrated using the neutralization
assay. The second part of the document describes some new ways of interpreting
titers, which involve, among other methods, singular value decomposition and
probabilistic multidimensional scaling. The third part of the document
discusses biological and mathematical issues related to the determination of
the effective dimensionality of titers, and describes an algorithm for
recovering unavailable titers.


  Darwin introduced the concept of the "living fossil" to describe species
belonging to lineages that have experienced little evolutionary change, and
suggested that species in more slowly evolving lineages are more prone to
extinction (1). Recent studies revealed that some living fossils such as the
lungfish are indeed evolving more slowly than other vertebrates (2, 3). The
reason for the slower rate of evolution in these lineages remains unclear, but
the same observations suggest a possible genome size effect on rates of
evolution. Genome size (C-value) in vertebrates varies over 200 fold ranging
from pufferfish (0.4 pg) to lungfish (132.8 pg) (4). Variation in genome size
and architecture is a fundamental cellular adaptation that remains poorly
understood (5). C-value is correlated with several allometric traits such as
body size and developmental rates in many, but not all, organisms (6, 7). To
date, no consensus exists concerning the mechanisms driving genome size
evolution or the effect that genome size has on species traits such as
evolutionary rates (8-12). In the following we show that: 1) within the same
range of divergence times, genetic diversity decreases as genome size increases
and 2) average rates of molecular evolution decline with increasing genome size
in vertebrates. Together, these observations indicate that genome size is an
important factor influencing rates of speciation and extinction.


  Prion diseases cover a large range of neurodegenerative diseases in humans
and animals, which are invariably fatal and highly infectious. By now there
have not been some effective therapeutic approaches or medications to treat all
prion diseases. Fortunately, numerous experimental experiences have showed that
rabbits are resistant to infection from prion diseases isolated from other
species, and recently the molecular structures of rabbit prion protein and its
mutants were released into protein data bank. Prion diseases are "protein
structural conformational" diseases. Thus, in order to reveal some secrets of
prion diseases, it is amenable to study rabbits by techniques of the molecular
structure and its dynamics. Wen et al. (PLoS One 5(10) e13273 (2010), Journal
of Biological Chemistry 285(41) 31682-31693 (2010)) reported the surface of NMR
RaPrPC(124-228) molecular snapshot has a large land of continuous positive
charge distribution, which contributes to the structural stability of rabbit
prion protein. This is just the property at one snapshot. This paper confirms a
large land of positive charge distribution has still been reserved during the
long molecular dynamics of 30 nanoseconds in many environments, but the
continuous of the land has not been always reserved yet. These results may be
useful for the medicinal treatment of prion diseases.


  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 authors' 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. Due to the noncrystalline and insoluble
nature of the amyloid fibril, little structural data on the AGAAAAGA is
available. This paper introduces a simple molecular modelling strategy 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.


  Many applications in genetic analyses utilize sampling distributions, which
describe the probability of observing a sample of DNA sequences randomly drawn
from a population. In the one-locus case with special models of mutation such
as the infinite-alleles model or the finite-alleles parent-independent mutation
model, closed-form sampling distributions under the coalescent have been known
for many decades. However, no exact formula is currently known for more general
models of mutation that are of biological interest. In this paper, models with
finitely-many alleles are considered, and an urn construction related to the
coalescent is used to derive approximate closed-form sampling formulas for an
arbitrary irreducible recurrent mutation model or for a reversible recurrent
mutation model, depending on whether the number of distinct observed allele
types is at most three or four, respectively. It is demonstrated empirically
that the formulas derived here are highly accurate when the per-base mutation
rate is low, which holds for many biological organisms.


  Certain invasive plants may rely on interference mechanisms (allelopathy,
e.g.) to gain competitive superiority over native species. But expending
resources on interference presumably exacts a cost in another life-history
trait, so that the significance of interference competition for invasion
ecology remains uncertain. We model ecological invasion when combined effects
of preemptive and interference competition govern interactions at the
neighborhood scale. We consider three cases. Under "novel weapons," only the
initially rare invader exercises interference. For "resistance zones" only the
resident species interferes, and finally we take both species as interference
competitors. Interference increases the other species' mortality, opening space
for colonization. However, a species exercising greater interference has
reduced propagation, which can hinder its colonization of open sites.
Interference never enhances a rare invader's growth in the homogeneously mixing
approximation to our model. But interference can significantly increase an
invader's competitiveness, and its growth when rare, if interactions are
structured spatially. That is, interference can increase an invader's success
when colonization of open sites depends on local, rather than global, species
densities. In contrast, interference enhances the common, resident species'
resistance to invasion independently of spatial structure, unless the
propagation-cost is too great. Increases in background mortality (i.e.,
mortality not due to interference) always reduce the effectiveness of
interference competition.


  A common feature of biological networks is the geometric property of
self-similarity. Molecular regulatory networks through to circulatory systems,
nervous systems, social systems and ecological trophic networks, show
self-similar connectivity at multiple scales. We analyze the relationship
between topology and signaling in contrasting classes of such topologies. We
find that networks differ in their ability to contain or propagate signals
between arbitrary nodes in a network depending on whether they possess
branching or loop-like features. Networks also differ in how they respond to
noise, such that one allows for greater integration at high noise, and this
performance is reversed at low noise. Surprisingly, small-world topologies,
with diameters logarithmic in system size, have slower dynamical timescales,
and may be less integrated (more modular) than networks with longer path
lengths. All of these phenomena are essentially mesoscopic, vanishing in the
infinite limit but producing strong effects at sizes and timescales relevant to
biology.


  We studied how to obtain a distribution for the number of ancestors in
species of sexual reproduction. Present models concentrate on the estimation of
distributions repetitions of ancestors in genealogical trees. It has been shown
that is not possible to reconstruct the genealogical history of each species
along all its generations by means of a geometric progression. This analysis
demonstrates that it is possible to rebuild the tree of progenitors by modeling
the problem with a Markov chain. For each generation, the maximum number of
possible ancestors is different. This brings huge problems for the resolution.
We found a solution through a dilation of the sample space, although the
distribution defined there takes smaller values respect to the initial problem.
In order to correct the distribution for each generation, we introduced the
invariance under gauge (local) group of dilations. These ideas can used to
study the interaction of several processes and provide a new approach on the
problem of the common ancestor. In the same direction, this model also provide
some elements that can be used to improve models of animal reproduction.


  Lattice-model simulations and experiments of some small proteins suggest that
folding is essentially controlled by a few conserved contacts. Residues of
these conserved contacts form the minimum set of native contacts needed to
ensure foldability. Keeping such conserved specific contacts in mind, we
examine contacts made by two secondary structure elements of different helices
or sheets and look for possible 'packing units' of the protein structure. Two
short backbone fragments of width five centred at the C? atoms in contact is
called an H-form, which serves as a candidate for the packing units. The
structural alignment of protein family members or even across families
indicates that there are conservative H-forms which are similar both in their
sequences and local geometry, and consistent with the structural alignment.
Carrying strong sequence signals, such packing units would provide 3D
constraints as a complement of the potential functions for the structure
prediction.


  We describe a large class of chemical reaction networks, those endowed with a
subtle structural property called concordance. We show that the class of
concordant networks coincides precisely with the class of networks which, when
taken with any weakly monotonic kinetics, invariably give rise to kinetic
systems that are injective --- a quality that, among other things, precludes
the possibility of switch-like transitions between distinct positive steady
states. We also provide persistence characteristics of concordant networks,
instability implications of discordance, and consequences of stronger variants
of concordance. Some of our results are in the spirit of recent ones by Banaji
and Craciun, but here we do not require that every species suffer a degradation
reaction. This is especially important in studying biochemical networks, for
which it is rare to have all species degrade.


  The availability of high-throughput parallel methods for sequencing microbial
communities is increasing our knowledge of the microbial world at an
unprecedented rate. Though most attention has focused on determining
lower-bounds on the alpha-diversity i.e. the total number of different species
present in the environment, tight bounds on this quantity may be highly
uncertain because a small fraction of the environment could be composed of a
vast number of different species. To better assess what remains unknown, we
propose instead to predict the fraction of the environment that belongs to
unsampled classes. Modeling samples as draws with replacement of colored balls
from an urn with an unknown composition, and under the sole assumption that
there are still undiscovered species, we show that conditionally unbiased
predictors and exact prediction intervals (of constant length in logarithmic
scale) are possible for the fraction of the environment that belongs to
unsampled classes. Our predictions are based on a Poissonization argument,
which we have implemented in what we call the Embedding algorithm. In fixed
i.e. non-randomized sample sizes, the algorithm leads to very accurate
predictions on a sub-sample of the original sample. We quantify the effect of
fixed sample sizes on our prediction intervals and test our methods and others
found in the literature against simulated environments, which we devise taking
into account datasets from a human-gut and -hand microbiota. Our methodology
applies to any dataset that can be conceptualized as a sample with replacement
from an urn. In particular, it could be applied, for example, to quantify the
proportion of all the unseen solutions to a binding site problem in a random
RNA pool, or to reassess the surveillance of a certain terrorist group,
predicting the conditional probability that it deploys a new tactic in a next
attack.


  Our work is motivated by and illustrated with application of association
networks in computational biology, specifically in the context of gene/protein
regulatory networks. Association networks represent systems of interacting
elements, where a link between two different elements indicates a sufficient
level of similarity between element attributes. While in reality relational
ties between elements can be expected to be based on similarity across multiple
attributes, the vast majority of work to date on association networks involves
ties defined with respect to only a single attribute. We propose an approach
for the inference of multi-attribute association networks from measurements on
continuous attribute variables, using canonical correlation and a
hypothesis-testing strategy. Within this context, we then study the impact of
partial information on multi-attribute network inference and characterization,
when only a subset of attributes is available. We consider in detail the case
of two attributes, wherein we examine through a combination of analytical and
numerical techniques the implications of the choice and number of node
attributes on the ability to detect network links and, more generally, to
estimate higher-level network summary statistics, such as node degree,
clustering coefficients, and measures of centrality. Illustration and
applications throughout the paper are developed using gene and protein
expression measurements on human cancer cell lines from the NCI-60 database.


  We investigate how the pattern of contacts between species in mutualistic
ecosystems is affected by the phylogenetic proximity between the species of
each guild. We develop a dynamical model geared to establish the role of such
proximity in the emergence of a nested pattern of contacts. We also define a
parameter that provides a direct measure of the influence of phylogenetic
proximity in a given pattern of contacts. We conclude that although
phylogenetic proximity is compatible with nestedness it can not be claimed to
be a cause of it. We find that nestedness can instead be attributed to a
general rule by which species tend to hold contacts with counterparts that
already have a large number of contacts. If the phylogenetic structure of both
guilds is brought into the analysis, this rule is equivalent to maximize the
phylogenetic diversity of the mutualistic counterparts of species of either
guild.


  Two-dimensional electrophoresis of proteins has preceded, and accompanied,
the birth of proteomics. Although it is no longer the only experimental scheme
used in modern proteomics, it still has distinct features and advantages. The
purpose of this tutorial paper is to guide the reader through the history of
the field, then through the main steps of the process, from sample preparation
to in-gel detection of proteins, commenting the constraints and caveats of the
technique. Then the limitations and positive features of two-dimensional
electrophoresis are discussed (e.g. its unique ability to separate complete
proteins and its easy interfacing with immunoblotting techniques), so that the
optimal type of applications of this technique in current and future proteomics
can be perceived. This is illustrated by a detailed example taken from the
literature and commented in detail. This Tutorial is part of the International
Proteomics Tutorial Programme (IPTP 2).


  We study a protein-DNA target search model with explicit DNA dynamics
applicable to in vitro experiments. We show that the DNA dynamics plays a
crucial role for the effectiveness of protein "jumps" between sites distant
along the DNA contour but close in 3D space. A strongly binding protein that
searches by 1D sliding and jumping alone, explores the search space less
redundantly when the DNA dynamics is fast on the timescale of protein jumps
than in the opposite "frozen DNA" limit. We characterize the crossover between
these limits using simulations and scaling theory. We also rationalize the slow
exploration in the frozen limit as a subtle interplay between long jumps and
long trapping times of the protein in "islands" within random DNA
configurations in solution.


  Recently, considerable effort has been put into developing fast algorithms to
reconstruct a rooted phylogenetic network that explains two rooted phylogenetic
trees and has a minimum number of hybridization vertices. With the standard
approach to tackle this problem being combinatorial, the reconstructed network
is rarely unique. From a biological point of view, it is therefore of
importance to not only compute one network, but all possible networks. In this
paper, we make a first step towards approaching this goal by presenting the
first algorithm---called allMAAFs---that calculates all
maximum-acyclic-agreement forests for two rooted binary phylogenetic trees on
the same set of taxa.


  Much of the complexity observed in gene regulation originates from
cooperative protein-DNA binding. While studies of the target search of proteins
for their specific binding sites on the DNA have revealed design principles for
the quantitative characteristics of protein-DNA interactions, no such
principles are known for the cooperative interactions between DNA-binding
proteins. We consider a simple theoretical model for two interacting
transcription factor (TF) species, searching for and binding to two adjacent
target sites hidden in the genomic background. We study the kinetic competition
of a dimer search pathway and a monomer search pathway, as well as the
steady-state regulation function mediated by the two TFs over a broad range of
TF-TF interaction strengths. Using a transcriptional AND-logic as exemplary
functional context, we identify the functionally desirable regime for the
interaction. We find that both weak and very strong TF-TF interactions are
favorable, albeit with different characteristics. However, there is also an
unfavorable regime of intermediate interactions where the genetic response is
prohibitively slow.


  We carry out an analytic investigation of stochastic oscillations in a
susceptible-infected-recovered model of disease spread on a network of $n$
cities. In the model a fraction $f_{jk}$ of individuals from city $k$ commute
to city $j$, where they may infect, or be infected by, others. Starting from a
continuous time Markov description of the model the deterministic equations,
which are valid in the limit when the population of each city is infinite, are
recovered. The stochastic fluctuations about the fixed point of these equations
are derived by use of the van Kampen system-size expansion. The fixed point
structure of the deterministic equations is remarkably simple: a unique
non-trivial fixed point always exists and has the feature that the fraction of
susceptible, infected and recovered individuals is the same for each city
irrespective of its size. We find that the stochastic fluctuations have an
analogously simple dynamics: all oscillations have a single frequency, equal to
that found in the one city case. We interpret this phenomenon in terms of the
properties of the spectrum of the matrix of the linear approximation of the
deterministic equations at the fixed point.


  We derive exact formulae for the allele frequency spectrum under the
coalescent with mutation, conditioned on allele counts at some fixed time in
the past. We consider unlinked biallelic markers mutating according to a finite
sites, or infinite sites, model. This work extends the coalescent theory of
unlinked biallelic markers, enabling fast computations of allele frequency
spectra in multiple populations. Our results have applications to demographic
inference, species tree inference, and the analysis of genetic variation in
closely related species more generally.


  We analyse the dynamics of overlapping antiparallel treadmilling microtubules
in the presence of crosslinking processive motor proteins that counterbalance
an external force. We show that coupling the force-dependent velocity of motors
and the kinetics of motor exchange with a bath in the presence of treadmilling
leads generically to oscillatory behavior. In addition we show that coupling
the polymerization kinetics to the external force through the kinetics of the
crosslinking motors can stabilize the oscillatory instability into
finite-amplitude nonlinear oscillations and may lead to other scenarios,
including bistability.


  During the past few years the development of experimental techniques has
allowed the quantitative analysis of biological systems ranging from
neurobiology and molecular biology. This work focuses on the quantitative
description of these systems by means of theoretical and numerical tools
ranging from statistical physics to probability theory.
  This dissertation is divided in three parts, each of which has a different
biological system as its focus.
  The first such system is Infotaxis, an olfactory search algorithm proposed by
Vergassola et al. in 2007: we give a continuous formulation and we characterize
its performances.
  Secondly we will focus on single-molecule experiments, especially unzipping
of DNA molecules, whose experimental traces depend strongly on the DNA
sequence: we develop a detailed model of the dynamics for this kind of
experiments and then we propose several inference algorithm aiming at the
characterization of the genetic sequence.
  The last section is devoted to the description of an algorithm that allows
the inference of interactions between neurons given the recording of neural
activity from multi-electrode experiments; we propose an integrated software
that will allow the analysis of these data.


  In the recent paper of Valleriani {\it et al} [Phys. Rev. E {\bf 83}, 042903
(2011)], a simple model for describing the translation of messenger RNA (mRNA)
by ribosomes is presented, and an expression of the translational ratio $r$,
defined as the ratio of translation rate $\omega_{\rm tl}$ of protein from mRNA
to degradation rate $\omega_p$ of protein, is obtained. The key point to get
this ratio $r$ is to get the translation rate $\omega_{\rm tl}$. In the study
of Valleriani {\it et al}, $\omega_{\rm tl}$ is assumed to be the mean value of
measured translation rate, i.e. the mean value of ratio of the translation
number of protein to the lifetime of mRNA. However, in experiments different
methods might be used to get $\omega_{\rm tl}$. Therefore, for the sake of
future application of their model to more experimental data analysis, in this
comment three methods to get the translation rate $\omega_{\rm tl}$, and
consequently the translational ratio $r$, are provided. Based on one of the
methods which might be employed in most of the experiments, we find that the
translational ratio $r$ decays exponentially with the length of mRNA in
prokaryotic cells, and decays reciprocally with the length of mRNA in
eukaryotic cells. This result is slight different from that obtained in
Valleriani's study.


  We develop an agent-based model of the motion and pattern formation of
vesicles. These intracellular particles can be found in four different modes of
(undirected and directed) motion and can fuse with other vesicles. While the
size of vesicles follows a log-normal distribution that changes over time due
to fusion processes, their spatial distribution gives rise to distinct
patterns. Their occurrence depends on the concentration of proteins which are
synthesized based on the transcriptional activities of some genes. Hence,
differences in these spatio-temporal vesicle patterns allow indirect
conclusions about the (unknown) impact of these genes.
  By means of agent-based computer simulations we are able to reproduce such
patterns on real temporal and spatial scales. Our modeling approach is based on
Brownian agents with an internal degree of freedom, $\theta$, that represents
the different modes of motion. Conditions inside the cell are modeled by an
effective potential that differs for agents dependent on their value $\theta$.
Agent's motion in this effective potential is modeled by an overdampted
Langevin equation, changes of $\theta$ are modeled as stochastic transitions
with values obtained from experiments, and fusion events are modeled as
space-dependent stochastic transitions. Our results for the spatio-temporal
vesicle patterns can be used for a statistical comparison with experiments. We
also derive hypotheses of how the silencing of some genes may affect the
intracellular transport, and point to generalizations of the model.


  We investigate the ejection dynamics of a ring polymer out of a cylindrical
nanochannel using both theoretical analysis and three dimensional Langevin
dynamics simulations. The ejection dynamics for ring polymers shows two regimes
like for linear polymers, depending on the relative length of the chain
compared with the channel. For long chains with length $N$ larger than the
critical chain length $N_{c}$, at which the chain just fully occupies the
nanochannel, the ejection for ring polymers is faster compared with linear
chains of identical length due to a larger entropic pulling force; while for
short chains ($N<N_c$), it takes longer time for ring polymers to eject out of
the channel due to a longer distance to be diffused to reach the exit of the
channel before experiencing the entropic pulling force. These results can help
understand many biological processes, such as bacterial chromosome segregation.


  Structural classification shows that the number of different protein folds is
surprisingly small. It also appears that proteins are built in a modular
fashion, from a relatively small number of components. Here we propose to
identify the modular building blocks of proteins with the dark soliton solution
of a generalized discrete nonlinear Schrodinger equation. For this we show that
practically all protein loops can be obtained simply by scaling the size and by
joining together a number of copies of the soliton, one after another. The
soliton has only two loop specific parameters and we identify their possible
values in Protein Data Bank. We show that with a collection of 200 sets of
parameters, each determining a soliton profile that describes a different short
loop, we cover over 90% of all proteins with experimental accuracy. We also
present two examples that describe how the loop library can be employed both to
model and to analyze the structure of folded proteins.


  Background: The accumulation of deleterious mutations of a population
directly contributes to the fate as to how long the population would exist.
Muller's ratchet provides a quantitative framework to study the effect of
accumulation. Adaptive landscape as a powerful concept in system biology
provides a handle to describe complex and rare biological events. In this
article we study the evolutionary process of a population exposed to Muller's
ratchet from the new viewpoint of adaptive landscape which allows us estimate
the single click of the ratchet starting with an intuitive understanding.
  Methods: We describe how Wright-Fisher process maps to Muller's ratchet. We
analytically construct adaptive landscape from general diffusion equation. It
shows that the construction is dynamical and the adaptive landscape is
independent of the existence and normalization of the stationary distribution.
We generalize the application of diffusion model from adaptive landscape
viewpoint.
  Results: We develop a novel method to describe the dynamical behavior of the
population exposed to Muller's ratchet, and analytically derive the decaying
time of the fittest class of populations as a mean first passage time. Most
importantly, we describe the absorption phenomenon by adaptive landscape, where
the stationary distribution is non-normalizable. These results suggest the
method may be used to understand the mechanism of populations evolution and
describe the biological processes quantitatively.


  We use the Litvinov-Maslov correspondence principle to reduce and hybridize
networks of biochemical reactions. We apply this method to a cell cycle
oscillator model. The reduced and hybridized model can be used as a hybrid
model for the cell cycle. We also propose a practical recipe for detecting
quasi-equilibrium QE reactions and quasi-steady state QSS species in
biochemical models with rational rate functions and use this recipe for model
reduction. Interestingly, the QE/QSS invariant manifold of the smooth model and
the reduced dynamics along this manifold can be put into correspondence to the
tropical variety of the hybridization and to sliding modes along this variety,
respectively


  Cooperative interactions among sensory receptors provide a general mechanism
to increase the sensitivity of signal transduction. In particular, bacterial
chemotaxis receptors interact cooperatively to produce an ultrasensitive
response to chemoeffector concentrations. However, cooperativity between
receptors in large macromolecular complexes is necessarily based on local
interactions and consequently is fundamentally connected to slowing of receptor
conformational dynamics, which increases intrinsic noise. Therefore, it is not
clear whether or under what conditions cooperativity actually increases the
precision of the concentration measurement. We explictly calculate the
signal-to-noise ratio (SNR) for sensing a concentration change using a simple,
Ising-type model of receptor-receptor interactions, generalized via scaling
arguments, and find that the optimal SNR is always achieved by independent
receptors.


  Despite the recognized importance of the multi-scale spatio-temporal
organization of proteins, most computational tools can only access a limited
spectrum of time and spatial scales, thereby ignoring the effects on protein
behavior of the intricate coupling between the different scales. Starting from
a physico-chemical atomistic network of interactions that encodes the structure
of the protein, we introduce a methodology based on multi-scale graph
partitioning that can uncover partitions and levels of organization of proteins
that span the whole range of scales, revealing biological features occurring at
different levels of organization and tracking their effect across scales.
Additionally, we introduce a measure of robustness to quantify the relevance of
the partitions through the generation of biochemically-motivated surrogate
random graph models. We apply the method to four distinct conformations of
myosin tail interacting protein, a protein from the molecular motor of the
malaria parasite, and study properties that have been experimentally addressed
such as the closing mechanism, the presence of conserved clusters, and the
identification through computational mutational analysis of key residues for
binding.


  Mass action type deterministic kinetic models of ion channels are usually
constructed in such a way as to obey the principle of detailed balance (or,
microscopic reversibility) for two reasons: first, the authors aspire to have
models harmonizing with thermodynamics, second, the conditions to ensure
detailed balance reduce the number of reaction rate coefficients to be
measured. We investigate a series of ion channel models which are asserted to
obey detailed balance, however, these models violate mass conservation and in
their case only the necessary conditions (the so-called circuit conditions) are
taken into account. We show that ion channel models have a very specific
structure which makes the consequences true in spite of the imprecise
arguments. First, we transform the models into mass conserving ones, second, we
show that the full set of conditions ensuring detailed balance (formulated by
Feinberg) leads to the same relations for the reaction rate constants in these
special cases, both for the original models and the transformed ones.


  Water near hydrophobic surfaces is like that at a liquid-vapor interface,
where fluctuations in water density are substantially enhanced compared to that
in bulk water. Here we use molecular simulations with specialized sampling
techniques to show that water density fluctuations are similarly enhanced, even
near hydrophobic surfaces of complex biomolecules, situating them at the edge
of a dewetting transition. Consequently, water near these surfaces is sensitive
to subtle changes in surface conformation, topology, and chemistry, any of
which can tip the balance towards or away from the wet state, and thus
significantly alter biomolecular interactions and function. Our work also
resolves the long-standing puzzle of why some biological surfaces dewet and
other seemingly similar surfaces do not.


  Living systems are forced away from thermodynamic equilibrium by exchange of
mass and energy with their environment. In order to model a biochemical
reaction network in a non-equilibrium state one requires a mathematical
formulation to mimic this forcing. We provide a general formulation to force an
arbitrary large kinetic model in a manner that is still consistent with the
existence of a non-equilibrium steady state. We can guarantee the existence of
a non-equilibrium steady state assuming only two conditions; that every
reaction is mass balanced and that continuous kinetic reaction rate laws never
lead to a negative molecule concentration. These conditions can be verified in
polynomial time and are flexible enough to permit one to force a system away
from equilibrium. In an expository biochemical example we show how a
reversible, mass balanced perpetual reaction, with thermodynamically infeasible
kinetic parameters, can be used to perpetually force a kinetic model of
anaerobic glycolysis in a manner consistent with the existence of a steady
state. Easily testable existence conditions are foundational for efforts to
reliably compute non-equilibrium steady states in genome-scale biochemical
kinetic models.


  The emergence of biochemical activities in a protein seem to commence with
the onset of atomic mean-square displacements along the protein lattice. The
ensuing protein dynamical transition has been discussed extensively in the
literature, and often with conflicting conclusions. Here we clarify the
phenomenon by establishing a deep connection between the dynamical transition
and the pseudogap state where high-temperature superconductivity comes to its
end.
  For this we first show how to endow proteins with an order parameter akin the
quasiparticle wave function in superconductors. We then present universality
arguments to claim that the protein dynamical transition takes place in tandem
with a pseudogap transmutation. We confirm that available experimental data
fully supports our proposal.


  Latent tree graphical models are widely used in computational biology, signal
and image processing, and network tomography. Here we design a new efficient,
estimation procedure for latent tree models, including Gaussian and discrete,
reversible models, that significantly improves on previous sample requirement
bounds. Our techniques are based on a new hidden state estimator which is
robust to inaccuracies in estimated parameters. More precisely, we prove that
latent tree models can be estimated with high probability in the so-called
Kesten-Stigum regime with $O(log^2 n)$ samples where $n$ is the number of
nodes.


  A brief informal overview on the BioPAX and SBML standards for formal
presentation of complex biological knowledge.


  In complex systems, feedback loops can build intricate emergent phenomena, so
that a description of the whole system cannot be easily derived from the
properties of the individual parts. Here we propose that inter-molecular
frustration mechanisms can provide non trivial feedback loops which can develop
nontrivial specificity amplification. We show that this mechanism can be seen
as a more general form of a kinetic proofreading mechanism, with an interesting
new property, namely the ability to tune the specificity amplification by
changing the reactants concentrations. This contrasts with the classical
kinetic proofreading mechanism in which specificity is a function of only the
reaction rate constants involved in a chemical pathway. These results are also
interesting because they show that a wide class of frustration models exists
that share the same underlining kinetic proofreading mechanisms, with even
richer properties. These models can find applications in different areas such
as evolutionary biology, immunology and biochemistry.


  We present an efficient phylogenetic reconstruction algorithm allowing
insertions and deletions which provably achieves a sequence-length requirement
(or sample complexity) growing polynomially in the number of taxa. Our
algorithm is distance-based, that is, it relies on pairwise sequence
comparisons. More importantly, our approach largely bypasses the difficult
problem of multiple sequence alignment.


  The impact of various types of external potentials to the Peyrard-Bishop DNA
denaturation is investigated through statistical mechanics approach. The
partition function is obtained using transfer integral method, and further the
stretching of hydrogen bond is calculated using time independent perturbation
method. It is shown that all types of external potentials accelerate the
denaturation processes at lower temperature. In particular, it is argued that
the Gaussian potential with infinitesimal width reproduces a constant force at
one end of DNA sequence as already done in some previous works.


  Quasispecies theory predicts that there is a critical mutation probability
above which a viral population will go extinct. Above this threshold the virus
loses the ability to replicate the best adapted genotype, leading to a
population composed of low replicating mutants that is eventually doomed. We
propose a new branching model that shows that this is not necessarily so. That
is, a population composed of ever changing mutants may survive.


  It has long been known that antibiotic treatment will not completely kill off
a bacteria population. For many species a small fraction of bacteria is not
sensitive to antibiotics. These bacteria are said to persist. Recently it has
been shown that persistence is not a permanent state and that in fact a
bacterium can switch back and forth between persistent and non persistent
states. We introduce two stochastic models for bacteria persistence. In both
models there are mass killings of non persistent bacteria at certain times. The
first model has deterministic killing times and the second one has random
killing times. Both models suggest that persistence may be a successful
strategy for a wide range of parameter values.


  We provide a Jacobian criterion that applies to arbitrary chemical reaction
networks taken with mass-action kinetics to preclude the existence of multiple
positive steady states within any stoichiometric class for any choice of rate
constants. We are concerned with the characterization of injective networks,
that is, networks for which the species formation rate function is injective in
the interior of the positive orthant within each stoichiometric class. We show
that a network is injective if and only if the determinant of the Jacobian of a
certain function does not vanish. The function consists of components of the
species formation rate function and a maximal set of independent conservation
laws. The determinant of the function is a polynomial in the species
concentrations and the rate constants (linear in the latter) and its
coefficients are fully determined. The criterion also precludes the existence
of degenerate steady states. Further, we relate injectivity of a chemical
reaction network to that of the chemical reaction network obtained by adding
outflow, or degradation, reactions for all species.


  Phosphorelays are a class of signaling mechanisms used by cells to respond to
changes in their environment. Phosphorelays (of which two-component systems
constitute a special case) are particularly abundant in prokaryotes and have
been shown to be involved in many fundamental processes such as stress
response, osmotic regulation, virulence, and chemotaxis. We develop a general
model of phosphorelays extending existing models of phosphorelays and
two-component systems. We analyze the model analytically under the assumption
of mass-action kinetics and prove that a phosphorelay has a unique stable
steady-state. Furthermore, we derive explicit functions relating stimulus to
the response in any layer of a phosphorelay and show that a limited degree of
ultrasensitivity (the ability to respond to changes in stimulus in a
switch-like manner) in the bottom layer of a phosphorelay is an intrinsic
feature which does not depend on any reaction rates or substrate amounts. On
the other hand, we show how adjusting reaction rates and substrate amounts may
lead to higher degrees of ultrasensitivity in intermediate layers. The explicit
formulas also enable us to prove how the response changes with alterations in
stimulus, kinetic parameters, and substrate amounts. Aside from providing
biological insight, the formulas may also be used to avoid time-consuming
simulations in numerical analyses and simulations.


  Allosteric drugs are increasingly used because they produce fewer side
effects. Allosteric signal propagation does not stop at the 'end' of a protein,
but may be dynamically transmitted across the cell. Here, we propose that the
concept of allosteric drugs can be broadened to allo-network drugs, whose
effects can propagate either within a protein, or across several proteins, to
enhance or inhibit specific interactions along a pathway. We posit that current
allosteric drugs are a special case of allo-network drugs, and suggest that
allo-network drugs can achieve specific, limited changes at the systems level,
and in this way can achieve fewer side effects and lower toxicity. Finally, we
propose steps and methods to identify allo-network drug targets and sites
outlining a new paradigm in systems-based drug design.


  In this paper we provide theoretical results that relate steady states of
continuous and discrete models arising from biology.


  RNA thermometers are mRNA strands with a temperature dependent secondary
structure: depending on the spatial conformation, the mRNA strand can get
translated (on-state) or can be inaccessible for ribosomes binding (off-state).
These have been found in a number of microorganisms (mainly pathogens), where
they are used to adaptively regulate the gene expression, in response to
changes in the environmental temperature. Besides naturally occurring RNA
thermometers, synthetic RNA thermometers have been recently designed by
modifying their natural counterparts (Hofacker et al., 2003). The newly
designed RNA thermometers are simpler, and exhibit a sharper switching between
off- and on-states. However, the proposed trial-and-error design procedure has
not been algorithmically formalized, and the switching temperature is rigidly
determined by the natural RNA thermometer used as template for the design.
  We developed a general algorithmic procedure (consensus distribution) for the
design of RNA thermo-switches with a tunable switching temperature that can be
decided in advance by the designer. A software tool with a user friendly GUI
has been written to automate the design of RNA thermo-switches with a desired
threshold temperature. Starting from a natural template, a new RNA thermometer
has been designed by our method for a new desired threshold temperature of 32C.
The designed RNA thermo-switch has been experimentally validated by using it to
control the expression of lucifarase. A 9.2 fold increase of luminescence has
been observed between 30C and 37C, whereas between 20C and 30C the luminescence
increase is less than 3-fold.
  This work represents a first step towards the design of flexible and tunable
RNA thermometers that can be used for a precise control of gene expression
without the need of external chemicals and possibly for temperature
measurements at a nano-scale resolution.


  We study the fraction $f$ of nucleotides involved in the formation of a
cactus--like secondary structure of random heteropolymer RNA--like molecules.
In the low--temperature limit we study this fraction as a function of the
number $c$ of different nucleotide species. We show, that with changing $c$,
the secondary structures of random RNAs undergo a morphological transition:
$f(c)\to 1$ for $c \le c_{\rm cr}$ as the chain length $n$ goes to infinity,
signaling the formation of a virtually "perfect" gapless secondary structure;
while $f(c)<1$ for $c>c_{\rm cr}$, what means that a non-perfect structure with
gaps is formed. The strict upper and lower bounds $2 \le c_{\rm cr} \le 4$ are
proven, and the numerical evidence for $c_{\rm cr}$ is presented. The relevance
of the transition from the evolutional point of view is discussed.


  Although taxonomy is often used informally to evaluate the results of
phylogenetic inference and find the root of phylogenetic trees, algorithmic
methods to do so are lacking. In this paper we formalize these procedures and
develop algorithms to solve the relevant problems. In particular, we introduce
a new algorithm that solves a "subcoloring" problem for expressing the
difference between the taxonomy and phylogeny at a given rank. This algorithm
improves upon the current best algorithm in terms of asymptotic complexity for
the parameter regime of interest; we also describe a branch-and-bound algorithm
that saves orders of magnitude in computation on real data sets. We also
develop a formalism and an algorithm for rooting phylogenetic trees according
to a taxonomy. All of these algorithms are implemented in freely-available
software.


  We study the standard SIS model of epidemic spreading on networks where
individuals have a fluctuating number of connections around a preferred degree
$\kappa $. Using very simple rules for forming such preferred degree networks,
we find some unusual statistical properties not found in familiar
Erd\H{o}s-R\'{e}nyi or scale free networks. By letting $\kappa $ depend on the
fraction of infected individuals, we model the behavioral changes in response
to how the extent of the epidemic is perceived. In our models, the behavioral
adaptations can be either `blind' or `selective' -- depending on whether a node
adapts by cutting or adding links to randomly chosen partners or selectively,
based on the state of the partner. For a frozen preferred network, we find that
the infection threshold follows the heterogeneous mean field result
$\lambda_{c}/\mu =<k>/<k^{2}>$ and the phase diagram matches the predictions of
the annealed adjacency matrix (AAM) approach. With `blind' adaptations,
although the epidemic threshold remains unchanged, the infection level is
substantially affected, depending on the details of the adaptation. The
`selective' adaptive SIS models are most interesting. Both the threshold and
the level of infection changes, controlled not only by how the adaptations are
implemented but also how often the nodes cut/add links (compared to the time
scales of the epidemic spreading). A simple mean field theory is presented for
the selective adaptations which capture the qualitative and some of the
quantitative features of the infection phase diagram.


  We investigate the unfolding pathway of a cytosine-rich DNA hairpin structure
via Molecular Dynamics simulations. Our results indicate a hidden complexity
present in the unfolding process. W e show that this complex behavior arises
due to non-stochastic contributions which have an important impact on the
description of the dynamics in terms of collective variables.


  Many studies have analyzed how variability in reproductive success affects
fitness. However, each study tends to focus on a particular problem, leaving
unclear the overall structure of variability in populations. This fractured
conceptual framework often causes particular applications to be incomplete or
improperly analyzed. In this paper, I present a concise introduction to the two
key aspects of the theory. First, all measures of fitness ultimately arise from
the relative comparison of the reproductive success of individuals or genotypes
with the average reproductive success in the population. That relative measure
creates a diminishing relation between reproductive success and fitness.
Diminishing returns reduce fitness in proportion to variability in reproductive
success. The relative measurement of success also induces a frequency
dependence that favors rare types. Second, variability in populations has a
hierarchical structure. Variable success in different traits of an individual
affects that individual's variation in reproduction. Correlation between
different individuals' reproduction affects variation in the aggregate success
of particular alleles across the population. One must consider the hierarchical
structure of variability in relation to different consequences of temporal,
spatial, and developmental variability. Although a complete analysis of
variability has many separate parts, this simple framework allows one to see
the structure of the whole and to place particular problems in their proper
relation to the general theory. The biological understanding of relative
success and the hierarchical structure of variability in populations may also
contribute to a deeper economic theory of returns under uncertainty.


  In classical evolutionary theory, genetic variation provides the source of
heritable phenotypic variation on which natural selection acts. Against this
classical view, several theories have emphasized that developmental variability
and learning enhance nonheritable phenotypic variation, which in turn can
accelerate evolutionary response. In this paper, I show how developmental
variability alters evolutionary dynamics by smoothing the landscape that
relates genotype to fitness. In a fitness landscape with multiple peaks and
valleys, developmental variability can smooth the landscape to provide a
directly increasing path of fitness to the highest peak. Developmental
variability also allows initial survival of a genotype in response to novel or
extreme environmental challenge, providing an opportunity for subsequent
adaptation. This initial survival advantage arises from the way in which
developmental variability smooths and broadens the fitness landscape.
Ultimately, the synergism between developmental processes and genetic variation
sets evolutionary rate.


  Summary: We describe a tool for quantifying the uniformity of mapped reads in
high-throughput sequencing experiments. Our statistic directly measures the
uniformity of both read position and fragment length, and we explain how to
compute a p-value that can be used to quantify biases arising from experimental
protocols and mapping procedures. Our method is useful for comparing different
protocols in experiments such as RNA-Seq.
  Availability and Implementation: We provide a freely available and open
source python script that can be used to analyze raw read data or reads mapped
to transcripts in BAM format at http://www.math.miami.edu/~vhower/ReadSpy.html .
  Contact: lpachter@math.berkeley.edu


  It has been generally recognized that stochasticity can play an important
role in the information processing accomplished by reaction networks in
biological cells. Most treatments of that stochasticity employ Gaussian noise
even though it is a priori obvious that this approximation can violate physical
constraints, such as the positivity of chemical concentrations. Here, we show
that even when such nonphysical fluctuations are rare, an exact solution of the
Gaussian model shows that the model can yield unphysical results. This is done
in the context of a simple incoherent-feedforward model which exhibits perfect
adaptation in the deterministic limit. We show how one can use the natural
separation of time scales in this model to yield an approximate model, that is
analytically solvable, including its dynamical response to an environmental
change. Alternatively, one can employ a cutoff procedure to regularize the
Gaussian result.


  It is widely recognized nowadays that complex diseases are caused by, amongst
the others, multiple genetic factors. The recent advent of genome-wide
association study (GWA) has triggered a wave of research aimed at discovering
genetic factors underlying common complex diseases. While the number of
reported susceptible genetic variants is increasing steadily, the application
of such findings into diseases prognosis for the general population is still
unclear, and there are doubts about whether the size of the contribution by
such factors is significant. In this respect, some recent simulation-based
studies have shed more light to the prospect of genetic tests. In this report,
we discuss several aspects of simulation-based studies: their parameters, their
assumptions, and the information they provide.