  We use a classical combinatorial inequality to establish a Markov inequality
for multivariate binary Markov processes on trees. We then apply this result,
alongside with the FKG inequality, to compare the expected loss of biodiversity
under two models of species extinction. One of these models is the generalized
version of an earlier model in which extinction is influenced by some trait
that can be classified into two states and which evolves on a tree according to
a Markov process. Since more than one trait can affect the rates of species
extinction, it is reasonable to allow, in the generalized model, k binary
states that influence extinction rates. We compare this model to one that has
matching marginal extinction probabilities for each species but for which the
species extinction events are stochastically independent.


  Network or graph structures are ubiquitous in the study of complex systems.
Often, we are interested in complexity trends of these system as it evolves
under some dynamic. An example might be looking at the complexity of a food web
as species enter an ecosystem via migration or speciation, and leave via
extinction.
  In a previous paper, a complexity measure of networks was proposed based on
the {\em complexity is information content} paradigm. To apply this paradigm to
any object, one must fix two things: a representation language, in which
strings of symbols from some alphabet describe, or stand for the objects being
considered; and a means of determining when two such descriptions refer to the
same object. With these two things set, the information content of an object
can be computed in principle from the number of equivalent descriptions
describing a particular object.
  The previously proposed representation language had the deficiency that the
fully connected and empty networks were the most complex for a given number of
nodes. A variation of this measure, called zcomplexity, applied a compression
algorithm to the resulting bitstring representation, to solve this problem.
Unfortunately, zcomplexity proved too computationally expensive to be
practical.
  In this paper, I propose a new representation language that encodes the
number of links along with the number of nodes and a representation of the
linklist. This, like zcomplexity, exhibits minimal complexity for fully
connected and empty networks, but is as tractable as the original measure.
  ...


  The quantum-mechanical mechanisms by which the enzymes catalyze the hydrogen
transfer in biochemical reactions are considered. Up to date it was established
both experimentally and theoretically that in many cases the proton tunnelling
through the intermolecular potential barrier is essential. We argue that in
this case the enzyme excitation and internal motion facilitate proton transfer
between reactants by squeezing the potential barrier which otherwise is
practically impenetrable. In the similar fashion, the enzymes can facilitate
the formation of hydrogen (H) bonds between the molecules. By means of barrier
squeezing, the enzymes not only facilitate such reactions but also can control
their rate and their final outcome, depending of enzyme excitation. In
particular, such effects can play the major role in DNA polymerization
reactions where preliminary DNTP selection is quite important.


  We have developed a generalized semi-analytic approach for efficiently
computing cyclization and looping $J$ factors of DNA under arbitrary binding
constraints. Many biological systems involving DNA-protein interactions impose
precise boundary conditions on DNA, which necessitates a treatment beyond the
Shimada-Yamakawa model for ring cyclization. Our model allows for DNA to be
treated as a heteropolymer with sequence-dependent intrinsic curvature and
stiffness. In this framework, we independently compute enthlapic and entropic
contributions to the $J$ factor and show that even at small length scales
$(\sim \ell_{p})$ entropic effects are significant. We propose a simple
analytic formula to describe our numerical results for a homogenous DNA in
planar loops, which can be used to predict experimental cyclization and loop
formation rates as a function of loop size and binding geometry. We also
introduce an effective torsional persistence length that describes the coupling
between twist and bending of DNA when looped.


  A linear external perturbation is introduced in the action of the partition
function of the random matrix model of RNA [G. Vernizzi, H. Orland and A. Zee,
Phys. Rev. Lett. 94, 168103 (2005)]. It is seen that (i). the perturbation
distinguishes between paired and unpaired bases in that there are structural
changes, from unpaired and paired base structures ($0 \leq \alpha < 1$) to
completely paired base structures ($\alpha=1$), as the perturbation parameter
$\alpha$ approaches 1 ($\alpha$ is the ratio of interaction strengths of
original and perturbed terms in the action of the partition function), (ii).
the genus distributions exhibit small differences for small even and odd
lengths $L$, (iii). the partition function of the linear interacting matrix
model is related via a scaling formula to the re-scaled partition function of
the random matrix model of RNA, (iv). the free energy and specific heat are
plotted as functions of $L$, $\alpha$ and temperature $T$ and their first
derivative with respect to $\alpha$ is plotted as a function of $\alpha$. The
free energy shows a phase transition at $\alpha=1$ for odd (both small and
large) lengths and for even lengths the transition at $\alpha=1$ gets sharper
and sharper as more pseudoknots are included (that is for large lengths).


  We demonstrate with a thought experiment that fitness-based population
dynamical approaches to evolution are not able to make quantitative,
falsifiable predictions about the long-term behavior of evolutionary systems. A
key characteristic of evolutionary systems is the ongoing endogenous production
of new species. These novel entities change the conditions for already existing
species. Even {\em Darwin's Demon}, a hypothetical entity with exact knowledge
of the abundance of all species and their fitness functions at a given time,
could not pre-state the impact of these novelties on established populations.
We argue that fitness is always {\it a posteriori} knowledge -- it measures but
does not explain why a species has reproductive success or not. To overcome
these conceptual limitations, a variational principle is proposed in a
spin-model-like setup of evolutionary systems. We derive a functional which is
minimized under the most general evolutionary formulation of a dynamical
system, i.e. evolutionary trajectories causally emerge as a minimization of a
functional. This functional allows the derivation of analytic solutions of the
asymptotic diversity for stochastic evolutionary systems within a mean-field
approximation. We test these approximations by numerical simulations of the
corresponding model and find good agreement in the position of phase
transitions in diversity curves. The model is further able to reproduce
stylized facts of timeseries from several man-made and natural evolutionary
systems. Light will be thrown on how species and their fitness landscapes
dynamically co-evolve.


  The system of partial differential equations describing tumor-immune dynamics
with angiogenesis taken into account is presented. For spatially homogeneous
case, the steady state analysis of the model is carried out. The effects of
single photodynamic impact are numerically simulated. In the case of strong
immune response we found that the photodynamic therapy (PDT) gives rise to the
substantial shrinkage of tumor size which is accompanied by the increase of
interleukin-2 concentration. On the contrary, the photodynamic stimulation of
weak immune response is shown to be insufficient to reduce the tumor. These
findings indicate the important role of anti-cancer immune response in the
long-term tumor control after PDT.


  We study the explosive character of the percolation transition in a
real-world network. We show that the emergence of a spanning cluster in the
Human Protein Homology Network (H-PHN) exhibits similar features to an
Achlioptas-type process and is markedly different from regular random
percolation. The underlying mechanism of this transition can be described by
slow-growing clusters that remain isolated until the later stages of the
process, when the addition of a small number of links leads to the rapid
interconnection of these modules into a giant cluster. Our results indicate
that the evolutionary-based process that shapes the topology of the H-PHN
through duplication-divergence events may occur in sudden steps, similarly to
what is seen in first-order phase transitions.


  We introduce a coarse-grained model of DNA with bases modeled as rigid-body
ellipsoids to capture their anisotropic stereochemistry. Interaction potentials
are all physicochemical and generated from all-atom simulation/parameterization
with minimal phenomenology. Persistence length, degree of stacking, and twist
are studied by molecular dynamics simulation as functions of temperature, salt
concentration, sequence, interaction potential strength, and local position
along the chain, for both single- and double-stranded DNA where appropriate.
The model of DNA shows several phase transitions and crossover regimes in
addition to dehybridization, including unstacking, untwisting, and collapse
which affect mechanical properties such as rigidity and persistence length. The
model also exhibits chirality with a stable right-handed and metastable
left-handed helix.


  From bird flocks to fish schools, animal groups often seem to react to
environmental perturbations as if of one mind. Most studies in collective
animal behaviour have aimed to understand how a globally ordered state may
emerge from simple behavioural rules. Less effort has been devoted to
understanding the origin of collective response, namely the way the group as a
whole reacts to its environment. Yet collective response is the adaptive key to
survivor, especially when strong predatory pressure is present. Here we argue
that collective response in animal groups is achieved through scale-free
behavioural correlations. By reconstructing the three-dimensional position and
velocity of individual birds in large flocks of starlings, we measured to what
extent the velocity fluctuations of different birds are correlated to each
other. We found that the range of such spatial correlation does not have a
constant value, but it scales with the linear size of the flock. This result
indicates that behavioural correlations are scale-free: the change in the
behavioural state of one animal affects and is affected by that of all other
animals in the group, no matter how large the group is. Scale-free correlations
extend maximally the effective perception range of the individuals, thus
compensating for the short-range nature of the direct inter-individual
interaction and enhancing global response to perturbations. Our results suggest
that flocks behave as critical systems, poised to respond maximally to
environmental perturbations.


  On the basis of the physico-chemical principles underlying silver-staining of
proteins, which are recalled in this paper, several methods of silver-staining
of proteins after SDS electrophoresis in polyacrylamide gels or isoelectric
focusing were tested. The most valuable protocols are presented in this report,
including standard methods for unsupported gels and new methods devised for
thin (0.5 mm) supported gels for SDS electrophoresis or isoelectric focusing
and for staining of small peptides. Generally speaking, the most rapid methods
were found to be less sensitive and less reproducible than more time-consuming
ones. Among the long methods, those using silver-diammine complex gave the most
uniform sensitivity. They require however special home-made gels and cannot be
applied to several electrophoretic systems (e.g. systems using tricine or
bicine as the trailing ion, or isoelectric focusing in immobilized pH
gradients). For these reasons, protocols based on silver nitrate are of a more
general use and might be favored. Future trends for silver-staining will also
be discussed.


  We present a bigraphical framework suited for modeling biological systems
both at protein level and at membrane level. We characterize formally bigraphs
corresponding to biologically meaningful systems, and bigraphic rewriting rules
representing biologically admissible interactions. At the protein level, these
bigraphic reactive systems correspond exactly to systems of kappa-calculus.
Membrane-level interactions are represented by just two general rules, whose
application can be triggered by protein-level interactions in a well-de\"ined
and precise way. This framework can be used to compare and merge models at
different abstraction levels; in particular, higher-level (e.g. mobility)
activities can be given a formal biological justification in terms of low-level
(i.e., protein) interactions. As examples, we formalize in our framework the
vesiculation and the phagocytosis processes.


  We introduce the BioBeta Framework, a meta-model for both protein-level and
membrane-level interactions of living cells. This formalism aims to provide a
formal setting where to encode, compare and merge models at different
abstraction levels; in particular, higher-level (e.g. membrane) activities can
be given a formal biological justification in terms of low-level (i.e.,
protein) interactions. A BioBeta specification provides a protein signature
together a set of protein reactions, in the spirit of the kappa-calculus.
Moreover, the specification describes when a protein configuration triggers one
of the only two membrane interaction allowed, that is "pinch" and "fuse". In
this paper we define the syntax and semantics of BioBeta, analyse its
properties, give it an interpretation as biobigraphical reactive systems, and
discuss its expressivity by comparing with kappa-calculus and modelling
significant examples. Notably, BioBeta has been designed after a bigraphical
metamodel for the same purposes. Hence, each instance of the calculus
corresponds to a bigraphical reactive system, and vice versa (almost).
Therefore, we can inherith the rich theory of bigraphs, such as the automatic
construction of labelled transition systems and behavioural congruences.


  The jasmonate (JA) signaling pathway in plants is activated as defense
response to a number of stresses like attacks by pests or pathogens and
wounding by animals. Some recent experiments provide significant new knowledge
on the molecular detail and connectivity of the pathway. The pathway has two
major components in the form of feedback loops, one negative and the other
positive. We construct a minimal mathematical model, incorporating the feedback
loops, to study the dynamics of the JA signaling pathway. The model exhibits
transient gene expression activity in the form of JA pulses in agreement with
experimental observations. The dependence of the pulse amplitude, duration and
peak time on the key parameters of the model is determined computationally. The
deterministic and stochastic aspects of the pathway dynamics are investigated
using both the full mathematical model as well as a reduced version of it. We
also compare the mechanism of pulse formation with the known mechanisms of
pulse generation in some bacterial and viral systems.


  We present the MIM calculus, a modeling formalism with a strong biological
basis, which provides biologically-meaningful operators for representing the
interaction capabilities of molecular species. The operators of the calculus
are inspired by the reaction symbols used in Molecular Interaction Maps (MIMs),
a diagrammatic notation used by biologists. Models of the calculus can be
easily derived from MIM diagrams, for which an unambiguous and executable
interpretation is thus obtained. We give a formal definition of the syntax and
semantics of the MIM calculus, and we study properties of the formalism. A case
study is also presented to show the use of the calculus for modeling
biomolecular networks.


  Considering all available non-redundant protein structures across different
structural classes, present study identified the probabilistic characteristics
that describe several facets of the occurrence of 3(10) and Pi-helices in
proteins. Occurrence profile of 3(10) and Pi-helices revealed that, their
presence follows Poisson flow on the primary structure; implying that, their
occurrence profile is rare, random and accidental. Structural class-specific
statistical analyses of sequence intervals between consecutive occurrences of
3(10) and Pi-helices revealed that these could be best described by gamma and
exponential distributions, across structural classes. Comparative study of
normalized percentage of non-glycine and non-proline residues in 3(10), Pi and
alpha-helices revealed a considerably higher proportion of 3(10) and Pi-helix
residues in disallowed, generous and allowed regions of Ramachandran map. Probe
into these findings in the light of evolution suggested clearly that 3(10) and
Pi-helices should appropriately be viewed as evolutionary intermediates on long
time scale, for not only the {\alpha}-helical conformation but also for the
'turns', equiprobably. Hence, accidental and random nature of occurrences of
3(10) and Pi-helices, and their evolutionary non-conservation, could be
described and explained from an invariant quantitative framework. Extent of
correctness of two previously proposed hypotheses on 3(10) and Pi-helices, have
been investigated too. Alongside these, a new algorithm to differentiate
between related sequences is proposed, which reliably studies evolutionary
distance with respect to protein secondary structures.


  The vast majority of biochemical systems involve the exchange of information
between different compartments, either in the form of transportation or via the
intervention of membrane proteins which are able to transmit stimuli between
bordering compartments. The correct quantitative handling of compartments is,
therefore, extremely important when modelling real biochemical systems. The
Bio-PEPA process algebra is equipped with the capability of explicitly defining
quantitative information such as compartment volumes and membrane surface
areas. Furthermore, the recent development of the Bio-PEPA Eclipse Plug-in
allows us to perform a correct stochastic simulation of multi-compartmental
models.
  Here we present a Bio-PEPA compartmental model of the cAMP/PKA/MAPK pathway.
We analyse the system using the Bio-PEPA Eclipse Plug-in and we show the
correctness of our model by comparison with an existing ODE model. Furthermore,
we perform computational experiments in order to investigate certain properties
of the pathway. Specifically, we focus on the system response to the inhibition
and strengthening of feedback loops and to the variation in the activity of key
pathway reactions and we observe how these modifications affect the behaviour
of the pathway. These experiments are useful to understand the control and
regulatory mechanisms of the system.


  The Stochastic Calculus of Looping Sequences is suitable to describe the
evolution of microbiological systems, taking into account the speed of the
described activities. We propose a type system for this calculus that models
how the presence of positive and negative catalysers can modify these speeds.
We claim that types are the right abstraction in order to represent the
interaction between elements without specifying exactly the element positions.
Our claim is supported through an example modelling the lactose operon.


  This paper explores the reduced form of horse cytochrome c confined in
reverse micelles (RM) of so-dium bis-(2-ethylhexyl) sulfosuccinate (AOT) in
isooctane by molecular dynamics simulation. RMs of two sizes were constructed
at a water content of Wo = [H2O]/[AOT] = 5.5 and 9.1. Our results show that the
protein secondary structure and the heme conformation both depend on micellar
hydration. At low hydration, the protein structure and the heme moiety remain
stable, whereas at high water content the protein becomes unstable and starts
to unfold. At Wo = 9.1, according to the X-ray structure, conforma-tional
changes are mainly localized on protein loops and around the heme moiety, where
we observe a partial opening of the heme crevice. These findings suggest that
within our time window (10 ns), the structural changes observed at the heme
level are the first steps of the protein denaturation process, pre-viously
described experimentally in micellar solutions. In addition, a specific binding
of AOT molecules to a few lysine residues of the protein was found only in the
small-sized RM.


  Gene regulatory networks constitute the first layer of the cellular
computation for cell adaptation and surveillance. In these webs, a set of
causal relations is built up from thousands of interactions between
transcription factors and their target genes. The large size of these webs and
their entangled nature make difficult to achieve a global view of their
internal organisation. Here, this problem has been addressed through a
comparative study for {\em Escherichia coli}, {\em Bacillus subtilis} and {\em
Saccharomyces cerevisiae} gene regulatory networks. We extract the minimal core
of causal relations, uncovering the hierarchical and modular organisation from
a novel dynamical/causal perspective. Our results reveal a marked top-down
hierarchy containing several small dynamical modules for \textit{E. coli} and
\textit{B. subtilis}. Conversely, the yeast network displays a single but large
dynamical module in the middle of a bow-tie structure. We found that these
dynamical modules capture the relevant wiring among both common and
organism-specific biological functions such as transcription initiation,
metabolic control, signal transduction, response to stress, sporulation and
cell cycle. Functional and topological results suggest that two fundamentally
different forms of logic organisation may have evolved in bacteria and yeast.


  During development, morphogen gradients precisely determine the position of
gene expression boundaries despite the inevitable presence of fluctuations.
Recent experiments suggest that some morphogen gradients may be interpreted
prior to reaching steady-state. Theoretical work has predicted that such
systems will be more robust to embryo-to-embryo fluctuations. By analysing two
experimentally motivated models of morphogen gradient formation, we investigate
the positional precision of gene expression boundaries determined by
pre-steady-state morphogen gradients in the presence of embryo-to-embryo
fluctuations, internal biochemical noise and variations in the timing of
morphogen measurement. Morphogens that are direct transcription factors are
found to be particularly sensitive to internal noise when interpreted prior to
steady-state, disadvantaging early measurement, even in the presence of large
embryo-to-embryo fluctuations. Morphogens interpreted by cell-surface receptors
can be measured prior to steady-state without significant decrease in
positional precision provided fluctuations in the timing of measurement are
small. Applying our results to experiment, we predict that Bicoid, a
transcription factor morphogen in Drosophila, is unlikely to be interpreted
prior to reaching steady-state. We also predict that Activin in Xenopus and
Nodal in zebrafish, morphogens interpreted by cell-surface receptors, can be
decoded in pre-steady-state.


  We re-examine the evolutionary dynamics of RNA secondary structures under
directional selection towards an optimum RNA structure. We find that the
punctuated equilibria lead to a very slow approach to the optimum, following on
average an inverse power of the evolutionary time. In addition, our study of
the trajectories shows that the out-of-equilibrium effects due to the
evolutionary process are very weak. In particular, the distribution of
genotypes is close to that arising during equilibrium stabilizing selection. As
a consequence, the evolutionary dynamics leave almost no measurable
out-of-equilibrium trace, only the transition genotypes (close to the border
between different periods of stasis) have atypical mutational properties.


  Recent work on synthetic rescues has shown that the targeted deletion of
specific metabolic genes can often be used to rescue otherwise non-viable
mutants. This raises a fundamental biophysical question: to what extent can the
whole-cell behavior of a large metabolic network be controlled by constraining
the flux of one or more reactions in the network? This touches upon the issue
of the number of degrees of freedom contained by one such network. Using the
metabolic network of Escherichia coli as a model system, here we address this
question theoretically by exploring not only reaction deletions, but also a
continuum of all possible reaction expression levels. We show that the behavior
of the metabolic network can be largely manipulated by the pinned expression of
a single reaction. In particular, a relevant fraction of the metabolic
reactions exhibits canalizing interactions, in that the specification of one
reaction flux determines cellular growth as well as the fluxes of most other
reactions in optimal steady states. The activity of individual reactions can
thus be used as surrogates to monitor and possibly control cellular growth and
other whole-cell behaviors. In addition to its implications for the study of
control processes, our methodology provides a new approach to study how the
integrated dynamics of the entire metabolic network emerges from the
coordinated behavior of its component parts.


  We consider neutral evolution of a large population subject to changes in its
population size. For a population with a time-variable carrying capacity we
have computed the distributions of the total branch lengths of its sample
genealogies. Within the coalescent approximation we have obtained a general
expression, Eq. (27), for the moments of these distributions for an arbitrary
smooth dependence of the population size on time. We investigate how the
frequency of population-size variations alters the distributions. This allows
us to discuss their influence on the distribution of the number of mutations,
and on the population homozygosity in populations with variable size.


  We use the Dominant Reaction Pathway (DRP) approach to study the dynamics of
the folding of a beta-hairpin, within a model which accounts for both native
and non-native interactions. We compare the most probable folding pathways
calculated with the DRP method with those obtained directly from molecular
dynamics (MD) simulations. We find that the two approaches give completely
consistent results. We investigate the effects of the non-native hydrophobic
interactions on the folding dynamics found them to be small.


  In 'no common mechanism' (NCM) models of character evolution, each character
can evolve on a phylogenetic tree under a partially or totally separate process
(e.g. with its own branch lengths).
  In such cases, the usual conditions that suffice to establish the statistical
consistency of tree reconstruction by methods such as maximum likelihood (ML)
break down, suggesting that such methods may be prone to statistical
inconsistency (SIN). In this paper we ask whether we can avoid SIN for tree
topology reconstruction when adopting such models, either by using ML or any
other method that could be devised. We prove that it is possible to avoid SIN
for certain NCM models, but not for others, and the results depend delicately
on the tree reconstruction method employed. We also describe the biological
relevance of some recent mathematical results for the more usual 'common
mechanism' setting. Our results are not intended to justify NCM, rather to set
in place a framework within which such questions can be formally addressed.


  We study the evolution of quantum entanglement during exciton energy transfer
(EET) in a network model of the Fenna-Matthews-Olson (FMO) complex, a
biological pigment-protein complex involved in the early steps of
photosynthesis in sulphur bacteria. The influence of Markovian, as well as
spatially and temporally correlated (non-Markovian) noise on the generation of
entanglement across distinct chromophores (site entanglement) and different
excitonic eigenstates (mode entanglement) is studied for different injection
mechanisms, including thermal and coherent laser excitation. Additionally, we
study the entangling power of the FMO complex under natural operating
conditions. While quantum information processing tends to favor maximal
entanglement, near unit EET is achieved as the result of an intricate interplay
between coherent and noisy processes where the initial part of the evolution
displays intermediate values of both forms of entanglement.


  One of the most striking effect of fluctuations in evolutionary game theory
is the possibility for mutants to fixate (take over) an entire population.
Here, we generalize a recent WKB-based theory to study fixation in evolutionary
games under non-vanishing selection, and investigate the relation between
selection intensity w and demographic (random) fluctuations. This allows the
accurate treatment of large fluctuations and yields the probability and mean
times of fixation beyond the weak selection limit. The power of the theory is
demonstrated on prototypical models of cooperation dilemmas with multiple
absorbing states. Our predictions compare excellently with numerical
simulations and, for finite w, significantly improve over those of the
Fokker-Planck approximation.


  Background: Network communities help the functional organization and
evolution of complex networks. However, the development of a method, which is
both fast and accurate, provides modular overlaps and partitions of a
heterogeneous network, has proven to be rather difficult. Methodology/Principal
Findings: Here we introduce the novel concept of ModuLand, an integrative
method family determining overlapping network modules as hills of an influence
function-based, centrality-type community landscape, and including several
widely used modularization methods as special cases. As various adaptations of
the method family, we developed several algorithms, which provide an efficient
analysis of weighted and directed networks, and (1) determine pervasively
overlapping modules with high resolution; (2) uncover a detailed hierarchical
network structure allowing an efficient, zoom-in analysis of large networks;
(3) allow the determination of key network nodes and (4) help to predict
network dynamics. Conclusions/Significance: The concept opens a wide range of
possibilities to develop new approaches and applications including network
routing, classification, comparison and prediction.


  Hybrid systems, and especially piecewise affine (PWA) systems, are often used
to model gene regulatory networks. In this paper we elaborate on previous work
about control problems for this class of models, using also some recent results
guaranteeing the existence and uniqueness of limit cycles, based solely on a
discrete abstraction of the system and its interaction structure. Our aim is to
control the transition graph of the PWA system to obtain an oscillatory
behaviour, which is indeed of primary functional importance in numerous
biological networks; we show how it is possible to control the appearance or
disappearance of a unique stable limit cycle by hybrid qualitative action on
the degradation rates of the PWA system, both by static and dynamic feedback,
i.e. the adequate coupling of a controlling subnetwork. This is illustrated on
two classical gene network modules, having the structure of mixed feedback
loops.


  Development of computer technology in chemistry, bring many application of
chemistry. Not only the application to visualize the structure of molecule but
also to molecular dynamics simulation. One of them is Gromacs. Gromacs is an
example of molecular dynamics application developed by Groningen University.
This application is a non-commercial and able to work in the operating system
Linux. The main ability of Gromacs is to perform molecular dynamics simulation
and minimization energy. In this paper, the author discusses about how to work
Gromacs in molecular dynamics simulation of some protein. In the molecular
dynamics simulation, Gromacs does not work alone. Gromacs interact with pymol
and Grace. Pymol is an application to visualize molecule structure and Grace is
an application in Linux to display graphs. Both applications will support
analysis of molecular dynamics simulation.


  Mechanical unfolding of the fourth domain of Distyostelium discoideum filamin
(DDFLN4) was studied in detail using the C$_{\alpha}$-Go model. We show that
unfolding pathways of this protein depend on the pulling speed. The agreement
between theoretical and experimental results on the sequencing of unfolding
events is achieved at low loading rates. The unfolding free energy landscape is
also constructed using dependencies of unfolding forces on pulling speeds.


  Mechanical unfolding of the fourth domain of Distyostelium discoideum filamin
(DDFLN4) was studied by all-atom molecular dynamics simulations, using the
GROMOS96 force field 43a1 and the simple point charge explicit water solvent.
Our study reveals an important role of non-native interactions in the unfolding
process. Namely, the existence of a peak centered at the end-to-end extension
22 nm in the force-extension curve, is associated with breaking of non-native
hydrogen bonds. Such a peak has been observed in experiments but not in Go
models, where non-native interactions are neglected. We predict that an
additional peak occurs at 2 nm using not only GROMOS96 force field 43a1 but
also Amber 94 and OPLS force fields. This result would stimulate further
experimental studies on elastic properties of DDFLN4.


  Urea titration of RNA by urea is an effective approach to investigate the
forces stabilizing this biologically important molecule. We used all atom
molecular dynamics simulations using two urea force fields and two RNA
constructs to elucidate in atomic detail the destabilization mechanism of
folded RNA in aqueous urea solutions. Urea denatures RNA by forming multiple
hydrogen bonds with the RNA bases and has little influence on the
phosphodiester backbone. Most significantly we discovered that urea engages in
stacking interactions with the bases. We also estimate, for the first time,
m-value for RNA, which is a measure of the strength of urea-RNA interactions.
Our work provides a conceptual understanding of the mechanism by which urea
enhances RNA folding rates.


  The ubiquity of oscillations in epidemics presents a long standing challenge
for the formulation of epidemic models. Whether they are external and
seasonally driven, or arise from the intrinsic dynamics is an open problem. It
is known that fixed time delays destabilize the steady state solution of the
standard SIRS model, giving rise to stable oscillations for certain parameters
values. In this contribution, starting from the classical SIRS model, we make a
general treatment of the recovery and loss of immunity terms. We present
oscillation diagrams (amplitude and period) in terms of the parameters of the
model, showing how oscillations can be destabilized by the shape of the
distributions of the two characteristic (infectious and immune) times. The
formulation is made in terms of delay equation which are both numerical
integrated and linearized. Results from simulation are included showing where
they support the linear analysis and explaining why not where they do not.
Considerations and comparison with real diseases are presented along.


  Background: The pandemic of influenza A (H1N1) is a serious on-going global
public crisis. Understanding its spreading dynamics is of fundamental
importance for both public health and scientific researches. Recent studies
have focused mainly on evaluation and prediction of on-going spreading, which
strongly depends on detailed information about the structure of social
contacts, human traveling patterns and biological activity of virus, etc.
  Methodology/Principal Findings: In this work we analyzed the distributions of
confirmed cases of influenza A (H1N1) in different levels and find the Zipf's
law and Heaps' law. Similar scaling properties were also observed for severe
acute respiratory syndrome (SARS) and bird cases of H5N1. We also found a
hierarchical spreading pattern from countries with larger population and GDP to
countries with smaller ones. We proposed a model that considers generic control
effects on both the local growth and transregional transmission, without the
need of the above mentioned detailed information. We studied in detail the
impact of control effects and heterogeneity on the spreading dynamics in the
model and showed that they are responsible for the scaling and hierarchical
spreading properties observed in empirical data.
  Conclusions/Significance: Our analysis and modeling showed that although
strict control measures for interregional travelers could delay the outbreak in
the regions without local cases, the focus should be turned to local prevention
after the outbreak of local cases. Target control on a few regions with the
largest number of active interregional travelers can efficiently prevent the
spreading. This work provided not only a deeper understanding of the generic
mechanisms underlying the spread of infectious diseases, but also some
practical guidelines for decision makers to adopt suitable control strategies.


  Apoptosis is a highly regulated cell death mechanism involved in many
physiological processes. A key component of extrinsically activated apoptosis
is the death receptor Fas, which, on binding to its cognate ligand FasL,
oligomerize to form the death-inducing signaling complex. Motivated by recent
experimental data, we propose a mathematical model of death ligand-receptor
dynamics where FasL acts as a clustering agent for Fas, which form locally
stable signaling platforms through proximity-induced receptor interactions.
Significantly, the model exhibits hysteresis, providing an upstream mechanism
for bistability and robustness. At low receptor concentrations, the bistability
is contingent on the trimerism of FasL. Moreover, irreversible bistability,
representing a committed cell death decision, emerges at high concentrations,
which may be achieved through receptor pre-association or localization onto
membrane lipid rafts. Thus, our model provides a novel theory for these
observed biological phenomena within the unified context of bistability.
Importantly, as Fas interactions initiate the extrinsic apoptotic pathway, our
model also suggests a mechanism by which cells may function as bistable
life/death switches independently of any such dynamics in their downstream
components. Our results highlight the role of death receptors in deciding cell
fate and add to the signal processing capabilities attributed to receptor
clustering.


  Major efforts to sequence cancer genomes are now occurring throughout the
world. Though the emerging data from these studies are illuminating, their
reconciliation with epidemiologic and clinical observations poses a major
challenge. In the current study, we provide a novel mathematical model that
begins to address this challenge. We model tumors as a discrete time branching
process that starts with a single driver mutation and proceeds as each new
driver mutation leads to a slightly increased rate of clonal expansion. Using
the model, we observe tremendous variation in the rate of tumor development -
providing an understanding of the heterogeneity in tumor sizes and development
times that have been observed by epidemiologists and clinicians. Furthermore,
the model provides a simple formula for the number of driver mutations as a
function of the total number of mutations in the tumor. Finally, when applied
to recent experimental data, the model allows us to calculate, for the first
time, the actual selective advantage provided by typical somatic mutations in
human tumors in situ. This selective advantage is surprisingly small, 0.005 +-
0.0005, and has major implications for experimental cancer research.


  Recently there have been several attempts to provide a whole set of
generators of the ideal of the algebraic variety associated to a phylogenetic
tree evolving under an algebraic model. These algebraic varieties have been
proven to be useful in phylogenetics. In this paper we prove that, for
phylogenetic reconstruction purposes, it is enough to consider generators
coming from the edges of the tree, the so-called edge invariants. This is the
algebraic analogous to Buneman's Splits Equivalence Theorem. The interest of
this result relies on its potential applications in phylogenetics for the
widely used evolutionary models such as Jukes-Cantor, Kimura 2 and 3
parameters, and General Markov models.


  This paper presents an algorithm for approximating certain types of dynamical
systems given by a system of ordinary delay differential equations by a Boolean
network model. Often Boolean models are much simpler to understand than complex
differential equations models. The motivation for this work comes from
mathematical systems biology. While Boolean mechanisms do not provide
information about exact concentration rates or time scales, they are often
sufficient to capture steady states and other key dynamics. Due to their
intuitive nature, such models are very appealing to researchers in the life
sciences. This paper is focused on dynamical systems that exhibit bistability
and are desc ribedby delay equations. It is shown that if a certain motif
including a feedback loop is present in the wiring diagram of the system, the
Boolean model captures the bistability of molecular switches. The method is
appl ied to two examples from biology, the lac operon and the phage lambda
lysis/lysogeny switch.


  With the aid of quantum mechanical calculations we investigate the electronic
structure of the full length (FL) potassium channel protein, FL-KcsA, in its
closed conformation, and the electronic structure of the ClC chloride channel.
The results indicate that both ion channels are strongly polarized towards the
extracellular region with respect to the membrane mean plane. FL-KcsA possesses
an electric dipole moment of magnitude 403 Debye while ClC has a macrodipole
whose magnitude is about five times larger, 1983 Debye, thereby contributing to
differentiate their membrane electric barriers. The dipole vectors of both
proteins are aligned along the corresponding selectivity filters. This result
suggests that potassium and chloride ion channels are not passive with respect
to the movement of ions across the membrane and the ionic motion might be
partially driven by the electric field of the protein in conjunction with the
electrochemical potential of the membrane.


  We use stochastic computer simulations to study the transport of a spherical
cargo particle along a microtubule-like track on a planar substrate by several
kinesin-like processive motors. Our newly developed adhesive motor dynamics
algorithm combines the numerical integration of a Langevin equation for the
motion of a sphere with kinetic rules for the molecular motors. The Langevin
part includes diffusive motion, the action of the pulling motors, and
hydrodynamic interactions between sphere and wall. The kinetic rules for the
motors include binding to and unbinding from the filament as well as active
motor steps. We find that the simulated mean transport length increases
exponentially with the number of bound motors, in good agreement with earlier
results. The number of motors in binding range to the motor track fluctuates in
time with a Poissonian distribution, both for springs and cables being used as
models for the linker mechanics. Cooperativity in the sense of equal load
sharing only occurs for high values for viscosity and attachment time.


  We investigate the nature of genetic drift acting at the leading edge of
range expansions, building on recent results in [Hallatschek et al., Proc.\
Natl.\ Acad.\ Sci., \textbf{104}(50): 19926 - 19930 (2007)]. A well mixed
population of two fluorescently labeled microbial species is grown in a
circular geometry. As the population expands, a coarsening process driven by
genetic drift gives rise to sectoring patterns with fractal boundaries, which
show a non-trivial asymptotic distribution. Using simplified lattice based
Monte Carlo simulations as a generic caricature of the above experiment, we
present detailed numerical results to establish a model for sector boundaries
as time changed Brownian motions. This is used to derive a general one-to-one
mapping of sector statistics between circular and linear geometries, which
leads to a full understanding of the sectoring patterns in terms of
annihilating diffusions.


  Simulated nucleotide sequences are widely used in theoretical and empirical
molecular evolution studies. Conventional simulators generally use fixed
parameter time-homogeneous Markov model for sequence evolution. In this work,
we use the folding free energy of the secondary structure of an RNA as a proxy
for its phenotypic fitness, and simulate RNA macroevolution by a
mutation-selection population genetics model. Because the two-step process is
conditioned on an RNA and its mutant ensemble, we no longer have a global
substitution matrix, nor do we explicitly assume any for this inhomogeneous
stochastic process. After introducing the base model of RNA evolution, we
outline the heuristic implementation algorithm and several model improvements.
We then discuss the calibration of the model parameters and demonstrate that in
phylogeny reconstruction with both the parsimony method and the likelihood
method, the sequences generated by our simulator, rnasim, have greater
statistical complexity than those by two standard simulators, ROSE and Seq-Gen,
and are close to empirical sequences.


  We present a simple mechanism for intra-molecular refrigeration, where parts
of a molecule are actively cooled below the environmental temperature. We
discuss the potential role and applications of such a mechanism in biology, in
particular in enzymatic reactions.


  Biochemical and regulatory interactions central to biological networks are
expected to cause extensive genetic interactions or epistasis affecting the
heritability of complex traits and the distribution of genotypes in
populations. However, the inference of epistasis from the observed
phenotype-genotype correlation is impeded by statistical difficulties, while
the theoretical understanding of the effects of epistasis remains limited, in
turn limiting our ability to interpret data. Of particular interest is the
biologically relevant situation of numerous interacting genetic loci with small
individual contributions to fitness. Here, we present a computational model of
selection dynamics involving many epistatic loci in a recombining population.
We demonstrate that a large number of polymorphic interacting loci can, despite
frequent recombination, exhibit cooperative behavior that locks alleles into
favorable genotypes leading to a population consisting of a set of competing
clones. When the recombination rate exceeds a certain critical value that
depends on the strength of epistasis, this "genotype selection" regime
disappears in an abrupt transition, giving way to "allele selection"-the regime
where different loci are only weakly correlated as expected in sexually
reproducing populations. We show that large populations attain highest fitness
at a recombination rate just below critical. Clustering of interacting sets of
genes on a chromosome leads to the emergence of an intermediate regime, where
blocks of cooperating alleles lock into genetic modules. These haplotype blocks
disappear in a second transition to pure allele selection. Our results
demonstrate that the collective effect of many weak epistatic interactions can
have dramatic effects on the population structure.


  The evolutionary dynamics of HIV during the chronic phase of infection is
driven by the host immune response and by selective pressures exerted through
drug treatment. To understand and model the evolution of HIV quantitatively,
the parameters governing genetic diversification and the strength of selection
need to be known. While mutation rates can be measured in single replication
cycles, the relevant effective recombination rate depends on the probability of
coinfection of a cell with more than one virus and can only be inferred from
population data. However, most population genetic estimators for recombination
rates assume absence of selection and are hence of limited applicability to
HIV, since positive and purifying selection are important in HIV evolution.
Here, we estimate the rate of recombination and the distribution of selection
coefficients from time-resolved sequence data tracking the evolution of HIV
within single patients. By examining temporal changes in the genetic
composition of the population, we estimate the effective recombination to be
r=1.4e-5 recombinations per site and generation. Furthermore, we provide
evidence that selection coefficients of at least 15% of the observed
non-synonymous polymorphisms exceed 0.8% per generation. These results provide
a basis for a more detailed understanding of the evolution of HIV. A
particularly interesting case is evolution in response to drug treatment, where
recombination can facilitate the rapid acquisition of multiple resistance
mutations. With the methods developed here, more precise and more detailed
studies will be possible, as soon as data with higher time resolution and
greater sample sizes is available.


  Understanding predator-prey relationships among insects is a challenging task
in the domain of insect-colony research. This is due to several factors
involved, such as determining whether a particular behavior is the result of a
predator-prey interaction, a friend-foe interaction or another kind of
interaction. In this paper, we analyze a series of predator-prey and friend-foe
interactions in two colonies of carpenter ants to better understand and predict
such behavior. Using the data gathered, we have also come up with a preliminary
model for predicting such behavior under the specific conditions the experiment
was conducted in. In this paper, we present the results of our data analysis as
well as an overview of the processes involved.


  Molecular phylogenetic techniques do not generally account for such common
evolutionary events as site insertions and deletions (known as indels). Instead
tree building algorithms and ancestral state inference procedures typically
rely on substitution-only models of sequence evolution. In practice these
methods are extended beyond this simplified setting with the use of heuristics
that produce global alignments of the input sequences--an important problem
which has no rigorous model-based solution. In this paper we consider a new
version of the multiple sequence alignment in the context of stochastic indel
models. More precisely, we introduce the following {\em trace reconstruction
problem on a tree} (TRPT): a binary sequence is broadcast through a tree
channel where we allow substitutions, deletions, and insertions; we seek to
reconstruct the original sequence from the sequences received at the leaves of
the tree. We give a recursive procedure for this problem with strong
reconstruction guarantees at low mutation rates, providing also an alignment of
the sequences at the leaves of the tree. The TRPT problem without indels has
been studied in previous work (Mossel 2004, Daskalakis et al. 2006) as a
bootstrapping step towards obtaining optimal phylogenetic reconstruction
methods. The present work sets up a framework for extending these works to
evolutionary models with indels.


  Increasing evidence of the effects of changing climate on physical ocean
conditions and long-term changes in fish populations adds to the need to
understand the effects of stochastic forcing on marine populations. Cohort
resonance is of particular interest because it involves selective sensitivity
to specific time scales of environmental variability, including that of mean
age of reproduction, and, more importantly, very low frequencies (i.e.,
trends). We present an age-structured model for two Pacific salmon species with
environmental variability in survival rate and in individual growth rate, hence
spawning age distribution. We use computed frequency response curves and
analysis of the linearized dynamics to obtain two main results. First, the
frequency response of the population is affected by the life history stage at
which variability affects the population; varying growth rate tends to excite
periodic resonance in age structure, while varying survival tends to excite
low-frequency fluctuation with more effect on total population size. Second,
decreasing adult survival strengthens the cohort resonance effect at all
frequencies, a finding that addresses the question of how fishing and climate
change will interact.


  Driven by advancements in high-throughput biological technologies and the
growing number of sequenced genomes, the construction of in silico models at
the genome scale has provided powerful tools to investigate a vast array of
biological systems and applications. Here, we review comprehensively the uses
of such models in industrial and medical biotechnology, including biofuel
generation, food production, and drug development. While the use of in silico
models is still in its early stages for delivering to industry, significant
initial successes have been achieved. For the cases presented here,
genome-scale models predict engineering strategies to enhance properties of
interest in an organism or to inhibit harmful mechanisms of pathogens. Going
forward, genome-scale in silico models promise to extend their application and
analysis scope to become a transformative tool in biotechnology.


  Recent years have witnessed an increasing interest in neuron-glia
communication. This interest stems from the realization that glia participates
in cognitive functions and information processing and is involved in many brain
disorders and neurodegenerative diseases. An important process in neuron-glia
communications is astrocyte encoding of synaptic information transfer: the
modulation of intracellular calcium dynamics in astrocytes in response to
synaptic activity. Here, we derive and investigate a concise mathematical model
for glutamate-induced astrocytic intracellular Ca2+ dynamics that captures the
essential biochemical features of the regulatory pathway of inositol
1,4,5-trisphosphate (IP3). Starting from the well-known two-state Li-Rinzel
model for calcium-induced-calcium release, we incorporate the regulation of the
IP3 production and phosphorylation. Doing so we extended it to a three-state
model (referred as the G-ChI model), that could account for Ca2+ oscillations
triggered by endogenous IP3 metabolism as well as by IP3 production by external
glutamate signals. Compared to previous similar models, our three-state models
include a more realistic description of the IP3 production and degradation
pathways, lumping together their essential nonlinearities within a concise
formulation. Using bifurcation analysis and time simulations, we demonstrate
the existence of new putative dynamical features. The cross-couplings between
IP3 and Ca2+ pathways endows the system with self-consistent oscillator
properties and favor mixed frequency-amplitude encoding modes over pure
amplitude modulation ones. These and additional results of our model are in
general agreement with available experimental data and may have important
implications on the role of astrocytes in the synaptic transfer of information.


  Many cells use calcium signalling to carry information from the extracellular
side of the plasma membrane to targets in their interior. Since virtually all
cells employ a network of biochemical reactions for Ca2+ signalling, much
effort has been devoted to understand the functional role of Ca2+ responses and
to decipher how their complex dynamics is regulated by the biochemical network
of Ca2+-related signal transduction pathways. Experimental observations show
that Ca2+ signals in response to external stimuli encode information via
frequency modulation or alternatively via amplitude modulation. Although
minimal models can capture separately both types of dynamics, they fail to
exhibit different and more advanced encoding modes. By arguments of bifurcation
theory, we propose instead that under some biophysical conditions more complex
modes of information encoding can also be manifested by minimal models. We
consider the minimal model of Li and Rinzel and show that information encoding
can occur by amplitude modulation (AM) of Ca2+ oscillations, by frequency
modulation (FM) or by both modes (AFM). Our work is motivated by calcium
signalling in astrocytes, the predominant type of cortical glial cells that is
nowadays recognized to play a crucial role in the regulation of neuronal
activity and information processing of the brain. We explain that our results
can be crucial for a better understanding of synaptic information transfer.
Furthermore, our results might also be important for better insight on other
examples of physiological processes regulated by Ca2+ signalling.


  With the development of high throughput sequencing technology, it becomes
possible to directly analyze mutation distribution in a genome-wide fashion,
dissociating mutation rate measurements from the traditional underlying
assumptions. Here, we sequenced several genomes of Escherichia coli from
colonies obtained after chemical mutagenesis and observed a strikingly
nonrandom distribution of the induced mutations. These include long stretches
of exclusively G to A or C to T transitions along the genome and orders of
magnitude intra- and inter-genomic differences in mutation density. Whereas
most of these observations can be explained by the known features of enzymatic
processes, the others could reflect stochasticity in the molecular processes at
the single-cell level. Our results demonstrate how analysis of the molecular
records left in the genomes of the descendants of an individual mutagenized
cell allows for genome-scale observations of fixation and segregation of
mutations, as well as recombination events, in the single genome of their
progenitor.


  We present a derivation of the classical SIR model through a mean-field
approximation from a discrete version of SIR. We then obtain a hyperbolic
forward Kolmogorov equation, and show that its projected characteristics
recover the standard SIR model. Moreover, we show that the long time limit of
the evolution will be a Dirac measure. The exact position will depend on the
well-know $R_0$ parameter, and it will be supported on the corresponding stable
SIR equilibrium.


  Extinction of a long-lived isolated stochastic population can be described as
an exponentially slow decay of quasi-stationary probability distribution of the
population size. We address extinction of a population in a two-population
system in the case when the population turnover -- renewal and removal -- is
much slower than all other processes. In this case there is a time scale
separation in the system which enables one to introduce a short-time
quasi-stationary extinction rate W_1 and a long-time quasi-stationary
extinction rate W_2, and develop a time-dependent theory of the transition
between the two rates. It is shown that W_1 and W_2 coincide with the
extinction rates when the population turnover is absent, and present but very
slow, respectively. The exponentially large disparity between the two rates
reflects fragility of the extinction rate in the population dynamics without
turnover.


  We present a quantitative measure of physical complexity, based on the amount
of information required to build a given physical structure through
self-assembly. Our procedure can be adapted to any given geometry, and thus to
any given type of physical system. We illustrate our approach using
self-assembling polyominoes, and demonstrate the breadth of its potential
applications by quantifying the physical complexity of molecules and protein
complexes. This measure is particularly well suited for the detection of
symmetry and modularity in the underlying structure, and allows for a
quantitative definition of structural modularity. Furthermore we use our
approach to show that symmetric and modular structures are favoured in
biological self-assembly, for example of protein complexes. Lastly, we also
introduce the notions of joint, mutual and conditional complexity, which
provide a useful distance measure between physical structures.


  Messenger RNA translation is often studied by means of statistical-mechanical
models based on the Asymmetric Simple Exclusion Process (ASEP), which considers
hopping particles (the ribosomes) on a lattice (the polynucleotide chain). In
this work we extend this class of models and consider the two fundamental steps
of the ribosome's biochemical cycle following a coarse-grained perspective. In
order to achieve a better understanding of the underlying biological processes
and compare the theoretical predictions with experimental results, we provide a
description lying between the minimal ASEP-like models and the more detailed
models, which are analytically hard to treat. We use a mean-field approach to
study the dynamics of particles associated with an internal stepping cycle. In
this framework it is possible to characterize analytically different phases of
the system (high density, low density or maximal current phase). Crucially, we
show that the transitions between these different phases occur at different
parameter values than the equivalent transitions in a standard ASEP, indicating
the importance of including the two fundamental steps of the ribosome's
biochemical cycle into the model.


  The problem of DNA-DNA interaction mediated by divalent counterions is
studied using computer simulation. Although divalent counterions cannot
condense free DNA molecules in solution, we show that if DNA configurational
entropy is restricted, divalent counterions can cause DNA reentrant
condensation similar to that caused by tri- or tetra-valent counterions.
DNA-DNA interaction is strongly repulsive at small or large counterion
concentration and is negligible or slightly attractive for a concentration in
between. Implications of our results to experiments of DNA ejection from
bacteriophages are discussed. The quantitative result serves to understand
electrostatic effects in other experiments involving DNA and divalent
counterions.


  For regulated exocytosis, we model the morphology and dynamics of the making
of the fusion pore or porosome as a cup-shaped lipoprotein structure (a dimple
or pit) on the cytosol side of the plasma membrane. We describe the forming of
the dimple by a free boundary problem. We discuss the various forces acting and
analyze the magnetic character of the wandering electromagnetic field wave
produced by intracellular spatially distributed pulsating (and well observed)
release and binding of calcium ions anteceding the bilayer membrane vesicle
fusion of exocytosis. Our approach explains the energy efficiency of the
observed dimple forming prior to hemifusion and fusion pore, and the observed
flickering in secretion. It provides a frame to relate characteristic time
length of exocytosis to the frequency, amplitude and direction of propagation
of the underlying electromagnetic field wave.


  Genetically identical cells under the same environmental conditions can show
strong variations in protein copy numbers due to inherently stochastic events
in individual cells. We here develop a theoretical framework to address how
variations in enzyme abundance affect the collective kinetics of metabolic
reactions observed within a population of cells. Kinetic parameters measured at
the cell population level are shown to be systematically deviated from those of
single cells, even within populations of homogeneous parameters. Because of
these considerations, Michaelis-Menten kinetics can even be inappropriate to
apply at the population level. Our findings elucidate a novel origin of
discrepancy between in vivo and in vitro kinetics, and offer potential utility
for analysis of single-cell metabolomic data.


  Chromatin is a complex of DNA, RNA and proteins whose primary function is to
package genomic DNA into the tight confines of a cell nucleus. A fundamental
repeating unit of chromatin is the nucleosome, an octamer of histone proteins
around which 147 base pairs of DNA are wound in almost two turns of a
left-handed superhelix. Chromatin is a dynamic structure which exerts profound
influence on regulation of gene expression and other cellular functions. These
chromatin-directed processes are facilitated by optimizing nucleosome positions
throughout the genome and by remodeling nucleosomes in response to various
external and internal signals such as environmental perturbations. Here we
discuss large-scale maps of nucleosome positions made available through recent
advances in parallel high-throughput sequencing and microarray technologies. We
show that these maps reveal common features of nucleosome organization in
eukaryotic genomes. We also survey computational models designed to predict
nucleosome formation scores or energies, and demonstrate how these predictions
can be used to position multiple nucleosome on the genome without steric
overlap.


  Recent experiments have demonstrated that the timescale of adaptation of
single neurons and ion channel populations to stimuli slows down as the length
of stimulation increases; in fact, no upper bound on temporal time-scales seems
to exist in such systems. Furthermore, patch clamp experiments on single ion
channels have hinted at the existence of large, mostly unobservable,
inactivation state spaces within a single ion channel. This raises the question
of the relation between this multitude of inactivation states and the observed
behavior. In this work we propose a minimal model for ion channel dynamics
which does not assume any specific structure of the inactivation state space.
The model is simple enough to render an analytical study possible. This leads
to a clear and concise explanation of the experimentally observed exponential
history-dependent relaxation in sodium channels in a voltage clamp setting, and
shows that their recovery rate from slow inactivation must be voltage
dependent. Furthermore, we predict that history-dependent relaxation cannot be
created by overly sparse spiking activity. While the model was created with ion
channel populations in mind, its simplicity and genericalness render it a good
starting point for modeling similar effects in other systems, and for scaling
up to higher levels such as single neurons which are also known to exhibit
multiple time scales.


  Epidemic spreading of infectious diseases is ubiquitous and has often
considerable impact on public health and economic wealth. The large variability
in spatio-temporal patterns of epidemics prohibits simple interventions and
demands for a detailed analysis of each epidemic with respect to its infectious
agent and the corresponding routes of transmission. To facilitate this
analysis, a mathematical framework is introduced which links epidemic patterns
to the topology and dynamics of the underlying transmission network. The
evolution both in disease prevalence and transmission network topology are
derived from a closed set of partial differential equations for infections
without recovery which are in excellent agreement with complementarily
conducted agent based simulations. The mutual influence between the epidemic
process and its transmission network is shown by several case studies on HIV
epidemics in synthetic populations. They reveal context dependent key processes
which drive the epidemic and which in turn can be exploited for targeted
intervention strategies. The mathematical framework provides a capable toolbox
to analyze epidemics from first principles. This allows for fast in silico
modeling - and manipulation - of epidemics which is especially powerful if
complemented with adequate empirical data for parametrization.


  A binary matrix has the Consecutive Ones Property (C1P) if its columns can be
ordered in such a way that all 1's on each row are consecutive. A Minimal
Conflicting Set is a set of rows that does not have the C1P, but every proper
subset has the C1P. Such submatrices have been considered in comparative
genomics applications, but very little is known about their combinatorial
structure and efficient algorithms to compute them. We first describe an
algorithm that detects rows that belong to Minimal Conflicting Sets. This
algorithm has a polynomial time complexity when the number of 1's in each row
of the considered matrix is bounded by a constant. Next, we show that the
problem of computing all Minimal Conflicting Sets can be reduced to the joint
generation of all minimal true clauses and maximal false clauses for some
monotone boolean function. We use these methods on simulated data related to
ancestral genome reconstruction to show that computing Minimal Conflicting Set
is useful in discriminating between true positive and false positive ancestral
syntenies. We also study a dataset of yeast genomes and address the reliability
of an ancestral genome proposal of the Saccahromycetaceae yeasts.


  We fit the parameters of a differential equations model describing the
production of gap gene proteins Hunchback and Knirps along the antero-posterior
axis of the embryo of \emph{Drosophila}. As initial data for the differential
equations model, we take the antero-posterior distribution of the proteins
Bicoid, Hunchback and Tailless at the beginning of cleavage cycle 14. We
calibrate and validate the model with experimental data using single- and
multi-objective evolutionary optimization techniques. In the multi-objective
optimization technique, we compute the associated Pareto fronts. We analyze the
cross regulation mechanism between the gap-genes protein pair Hunchback-Knirps
and we show that the posterior distribution of Hunchback follow the
experimental data if Hunchback is negatively regulated by the Huckebein
protein. This approach enables to predict the posterior localization on the
embryo of the protein Huckebein, and we validate with the experimental data the
genetic regulatory network responsible for the antero-posterior distribution of
the gap gene protein Hunchback. We discuss the importance of Pareto
multi-objective optimization techniques in the calibration and validation of
biological models.


  Post-translational modifications of histone proteins are an important factor
in epigenetic control that serve to regulate transcription, depending on the
particular modification states of the histone proteins. We study the stochastic
dynamics of histone protein states, taking into account a feedback mechanism
where modified nucleosomes recruit enzymes that diffuse to adjacent
nucleosomes. We map the system onto a quantum spin system whose dynamics is
generated by a non-Hermitian Hamiltonian. Making an ansatz for the solution as
a tensor product state leads to nonlinear partial differential equations that
describe the dynamics of the system. Multiple stable histone states appear in a
parameter regime whose size increases with increasing number of modification
sites. We discuss the role of the spatial dependance, and we consider the
effects of spatially heterogeneous enzymatic activity. Finally, we consider
multistability in a model of several types of correlated post-translational
modifications.


  Gene trees are evolutionary trees representing the ancestry of genes sampled
from multiple populations. Species trees represent populations of individuals
-- each with many genes -- splitting into new populations or species. The
coalescent process, which models ancestry of gene copies within populations, is
often used to model the probability distribution of gene trees given a fixed
species tree. This multispecies coalescent model provides a framework for
phylogeneticists to infer species trees from gene trees using maximum
likelihood or Bayesian approaches. Because the coalescent models a branching
process over time, all trees are typically assumed to be rooted in this
setting. Often, however, gene trees inferred by traditional phylogenetic
methods are unrooted.
  We investigate probabilities of unrooted gene trees under the multispecies
coalescent model. We show that when there are 4 species with one gene sampled
per species, the distribution of unrooted gene tree topologies identifies the
unrooted species tree topology and some, but not all, information in the
species tree edges (branch lengths). The location of the root on the species
tree is not identifiable in this situation. However, for 5 or more species with
one gene sampled per species, we show that the distribution of unrooted gene
tree topologies identifies the rooted species tree topology and all its
internal branch lengths. The length of any pendent branch leading to a leaf of
the species tree is also identifiable for any species from which more than one
gene is sampled.


  In [2] it was proven that the Cass algorithm is a polynomial-time algorithm
for constructing level<=2 networks from clusters. Here we demonstrate, for each
k>=0, a polynomial-time algorithm for constructing level-k phylogenetic
networks from clusters. Unlike Cass the algorithm scheme given here is only of
theoretical interest. It does, however, strengthen the hope that efficient
polynomial-time algorithms (and perhaps fixed parameter tractable algorithms)
exist for this problem.


  Finding optimal policies to reduce the morbidity and mortality of the ongoing
pandemic is a top public health priority. Using a compartmental model with age
structure and vaccination status, we examined the effect of age specific
scheduling of vaccination during a pandemic influenza outbreak, when there is a
race between the vaccination campaign and the dynamics of the pandemic. Our
results agree with some recent studies on that age specificity is paramount to
vaccination planning. However, little is known about the effectiveness of such
control measures when they are applied during the outbreak. Comparing five
possible strategies, we found that age specific scheduling can have a huge
impact on the outcome of the epidemic. For the best scheme, the attack rates
were up to 10% lower than for other strategies. We demonstrate the importance
of early start of the vaccination campaign, since ten days delay may increase
the attack rate by up to 6%. Taking into account the delay between developing
immunity and vaccination is a key factor in evaluating the impact of
vaccination campaigns. We provide a general framework which will be useful for
the next pandemic waves as well.


  Biochemical computing is an emerging field of unconventional computing that
attempts to process information with biomolecules and biological objects using
digital logic. In this work we survey filtering in general, in biochemical
computing, and summarize the experimental realization of an AND logic gate with
sigmoid response in one of the inputs. The logic gate is realized with
electrode-immobilized glucose-6-phosphate dehydrogenase enzyme that catalyzes a
reaction corresponding to the Boolean AND functions. A kinetic model is also
developed and used to evaluate the extent to which the performance of the
experimentally realized logic gate is close to optimal.


  This paper develops a formulation of the quasispecies equations appropriate
for polysomic, semiconservatively replicating genomes. This paper is an
extension of previous work on the subject, which considered the case of haploid
genomes. Here, we develop a more general formulation of the quasispecies
equations that is applicable to diploid and even polyploid genomes.
Interestingly, with an appropriate classification of population fractions, we
obtain a system of equations that is formally identical to the haploid case. As
with the work for haploid genomes, we consider both random and immortal DNA
strand chromosome segregation mechanisms. However, in contrast to the haploid
case, we have found that an analytical solution for the mean fitness is
considerably more difficult to obtain for the polyploid case. Accordingly,
whereas for the haploid case we obtained expressions for the mean fitness for
the case of an analogue of the single-fitness-peak landscape for arbitrary
lesion repair probabilities (thereby allowing for non-complementary genomes),
here we solve for the mean fitness for the restricted case of perfect lesion
repair.


  Dramatic rise of mutators has been found to accompany adaptation of bacteria
in response to many kinds of stress. Two views on the evolutionary origin of
this phenomenon emerged: the pleiotropic hypothesis positing that it is a
byproduct of environmental stress or other specific stress response mechanisms
and the second order selection which states that mutators hitchhike to fixation
with unrelated beneficial alleles. Conventional population genetics models
could not fully resolve this controversy because they are based on certain
assumptions about fitness landscape. Here we address this problem using a
microscopic multiscale model, which couples physically realistic molecular
descriptions of proteins and their interactions with population genetics of
carrier organisms without assuming any a priori fitness landscape. We found
that both pleiotropy and second order selection play a crucial role at
different stages of adaptation: the supply of mutators is provided through
destabilization of error correction complexes or fluctuations of production
levels of prototypic mismatch repair proteins (pleiotropic effects), while rise
and fixation of mutators occur when there is a sufficient supply of beneficial
mutations in replication-controlling genes. This general mechanism assures a
robust and reliable adaptation of organisms to unforeseen challenges. This
study highlights physical principles underlying physical biological mechanisms
of stress response and adaptation.


  The quasispecies model describes processes related to the origin of life and
viral evolutionary dynamics. We discuss how the error catastrophe that reflects
the transition from localized to delocalized quasispecies population is
affected by catalytic replication of different reaction orders. Specifically,
we find that 2nd order mechanisms lead to 1st order discontinuous phase
transitions in the viable population fraction, and conclude that the "higher"
the interaction the "lower" the transition. We discuss potential implications
for understanding the replication of highly mutating RNA viruses.


  This work illustrates potentials for recognition within {\em ad hoc} sensor
networks if their nodes possess individual inter-related biologically inspired
genetic codes. The work takes ideas from natural immune systems protecting
organisms from infection. Nodes in the present proposal have individual gene
sets fitting into a self organised phylogenetic tree. Members of this
population are genetically ''relatives''. Outsiders cannot easily copy or
introduce a new node in the network without going through a process of
conception between two nodes in the population. Related nodes can locally
decide to check each other for their genetic relation without directly
revealing their gene sets. A copy/clone of a gene sequence or a random gene set
will appear as alien. Nodes go through a cycle of introduction (conception or
''birth'') with parents in the network and later exit from it (''death'').
Hence the phylogenetic tree is dynamic or possesses a genetic drift. Typical
lifetimes of gene sets and number of offspring from different parents affect
this genetic drift and the level of correlation between gene sets. The
frequency of mutations similarly affects the gene pool. Correlation between
genes of the nodes implies a common secret for cryptographic material for
communication and consistency check facilitating intrusion detection and
tracing of events. A node can by itself (non-specifically) recognise an
adversary if it does not respond properly according to its genes. Nodes can
also collaborate to recognise adversaries by communicating response from
intruders to check for consistency with the whole gene pool (phylogenetic
tree).


  Recognition of pathogens relies on families of proteins showing great
diversity. Here we construct maximum entropy models of the sequence repertoire,
building on recent experiments that provide a nearly exhaustive sampling of the
IgM sequences in zebrafish. These models are based solely on pairwise
correlations between residue positions, but correctly capture the higher order
statistical properties of the repertoire. Exploiting the interpretation of
these models as statistical physics problems, we make several predictions for
the collective properties of the sequence ensemble: the distribution of
sequences obeys Zipf's law, the repertoire decomposes into several clusters,
and there is a massive restriction of diversity due to the correlations. These
predictions are completely inconsistent with models in which amino acid
substitutions are made independently at each site, and are in good agreement
with the data. Our results suggest that antibody diversity is not limited by
the sequences encoded in the genome, and may reflect rapid adaptation to
antigenic challenges. This approach should be applicable to the study of the
global properties of other protein families.


  We study the oscillatory dynamics in the generic three-species
rock-paper-scissors games with mutations. In the mean-field limit, different
behaviors are found: (a) for high mutation rate, there is a stable interior
fixed point with coexistence of all species; (b) for low mutation rates, there
is a region of the parameter space characterized by a limit cycle resulting
from a Hopf bifurcation; (c) in the absence of mutations, there is a region
where heteroclinic cycles yield oscillations of large amplitude (not robust
against noise). After a discussion on the main properties of the mean-field
dynamics, we investigate the stochastic version of the model within an
individual-based formulation. Demographic fluctuations are therefore naturally
accounted and their effects are studied using a diffusion theory complemented
by numerical simulations. It is thus shown that persistent erratic oscillations
(quasi-cycles) of large amplitude emerge from a noise-induced resonance
phenomenon. We also analytically and numerically compute the average escape
time necessary to reach a (quasi-)cycle on which the system oscillates at a
given amplitude.


  Phylogenomics, even more so than traditional phylogenetics, needs to
represent the uncertainty in evolutionary trees due to systematic error. Here
we illustrate the analysis of genome-scale alignments of yeast, using robust
measures of the additivity of the fit of distances to tree when using flexi
Weighted Least Squares. A variety of DNA and protein distances are used. We
explore the nature of the residuals, standardize them, and then create
replicate data sets by resampling these residuals. Under the model, the results
are shown to be very similar to the conventional sequence bootstrap. With real
data they show up uncertainty in the tree that is either due to underestimating
the stochastic error (hence massively overestimating the effective sequence
length) and/or systematic error. The methods are extended to the very fast BME
criterion with similarly promising results.


  Huge numbers of short reads are being generated for mapping back to the
genome to discover the frequency of transcripts, miRNAs, DNAase hypersensitive
sites, FAIRE regions, nucleosome occupancy, etc. Since these reads are
typically short (e.g., 36 base pairs) and since many eukaryotic genomes,
including humans, have highly repetitive sequences then many of these reads map
to two or more locations in the genome. Current mapping of these reads, grading
them according to 0, 1 or 2 mismatches wastes a great deal of information.
These short sequences are typically mapped with no account of the accuracy of
the sequence, even in company software when per base error rates are being
reported by another part of the machine. Further, multiply mapping locations
are frequently discarded altogether or allocated with no regard to where other
reads are accumulating. Here we show how to combine probabilistic mapping of
reads with an EM algorithm to iteratively improve the empirical likelihood of
the allocation of short reads. Mapping using LAST takes into account the per
base accuracy of the read, plus insertions and deletions, plus anticipated
occasional errors or SNPs with respect to the parent genome. The probabilistic
EM algorithm iteratively allocates reads based on the proportion of reads
mapping within windows on the previous cycle, along with any prior information
on where the read best maps. The methods are illustrated with FAIRE ENCODE data
looking at the very important head-to-head gene combination of TP53 and WRAP
53.


  Gene regulatory networks (GRN) are being studied with increasingly precise
quantitative tools and can provide a testing ground for ideas regarding the
emergence and evolution of complex biological networks. We analyze the global
statistical properties of the transcriptional regulatory network of the
prokaryote Escherichia coli, identifying each operon with a node of the
network. We propose a null model for this network using the content-based
approach applied earlier to the eukaryote Saccharomyces cerevisiae. (Balcan et
al., 2007) Random sequences that represent promoter regions and binding
sequences are associated with the nodes. The length distributions of these
sequences are extracted from the relevant databases. The network is constructed
by testing for the occurrence of binding sequences within the promoter regions.
The ensemble of emergent networks yields an exponentially decaying in-degree
distribution and a putative power law dependence for the out-degree
distribution with a flat tail, in agreement with the data. The clustering
coefficient, degree-degree correlation, rich club coefficient and k-core
visualization all agree qualitatively with the empirical network to an extent
not yet achieved by any other computational model, to our knowledge. The
significant statistical differences can point the way to further research into
non-adaptive and adaptive processes in the evolution of the E. coli GRN.


  Geometric and structural constraints greatly restrict the selection of folds
adapted by protein backbones, and yet, folded proteins show an astounding
diversity in functionality. For structure to have any bearing on function, it
is thus imperative that, apart from the protein backbone, other tunable degrees
of freedom be accountable. Here, we focus on side-chain interactions, which
non-covalently link amino acids in folded proteins to form a network structure.
At a coarse-grained level, we show that the network conforms remarkably well to
realizations of random graphs and displays associated percolation behavior.
Thus, within the rigid framework of the protein backbone that restricts the
structure space, the side-chain interactions exhibit an element of randomness,
which account for the functional flexibility and diversity shown by proteins.
However, at a finer level, the network exhibits deviations from these random
graphs which, as we demonstrate for a few specific examples, reflect the
intrinsic uniqueness in the structure and stability, and perhaps specificity in
the functioning of biological proteins.


  Networks are ubiquitous in science and have become a focal point for
discussion in everyday life. Formal statistical models for the analysis of
network data have emerged as a major topic of interest in diverse areas of
study, and most of these involve a form of graphical representation.
Probability models on graphs date back to 1959. Along with empirical studies in
social psychology and sociology from the 1960s, these early works generated an
active network community and a substantial literature in the 1970s. This effort
moved into the statistical literature in the late 1970s and 1980s, and the past
decade has seen a burgeoning network literature in statistical physics and
computer science. The growth of the World Wide Web and the emergence of online
networking communities such as Facebook, MySpace, and LinkedIn, and a host of
more specialized professional network communities has intensified interest in
the study of networks and network data. Our goal in this review is to provide
the reader with an entry point to this burgeoning literature. We begin with an
overview of the historical development of statistical network modeling and then
we introduce a number of examples that have been studied in the network
literature. Our subsequent discussion focuses on a number of prominent static
and dynamic network models and their interconnections. We emphasize formal
model descriptions, and pay special attention to the interpretation of
parameters and their estimation. We end with a description of some open
problems and challenges for machine learning and statistics.


  Proteins with nontrivial topology, containing knots and slipknots, have the
ability to fold to their native states without any additional external forces
invoked. A mechanism is suggested for folding of these proteins, such as YibK
and YbeA, which involves an intermediate configuration with a slipknot. It
elucidates the role of topological barriers and backtracking during the folding
event. It also illustrates that native contacts are sufficient to guarantee
folding in around 1-2% of the simulations, and how slipknot intermediates are
needed to reduce the topological bottlenecks. As expected, simulations of
proteins with similar structure but with knot removed fold much more
efficiently, clearly demonstrating the origin of these topological barriers.
Although these studies are based on a simple coarse-grained model, they are
already able to extract some of the underlying principles governing folding in
such complex topologies.


  Central to the functioning of a living cell is its ability to control the
readout or expression of information encoded in the genome. In many cases, a
single transcription factor protein activates or represses the expression of
many genes. As the concentration of the transcription factor varies, the target
genes thus undergo correlated changes, and this redundancy limits the ability
of the cell to transmit information about input signals. We explore how
interactions among the target genes can reduce this redundancy and optimize
information transmission. Our discussion builds on recent work [Tkacik et al,
Phys Rev E 80, 031920 (2009)], and there are connections to much earlier work
on the role of lateral inhibition in enhancing the efficiency of information
transmission in neural circuits; for simplicity we consider here the case where
the interactions have a feed forward structure, with no loops. Even with this
limitation, the networks that optimize information transmission have a
structure reminiscent of the networks found in real biological systems.


  Theoretical studies of stretching proteins with slipknots reveal a surprising
growth of their unfolding times when the stretching force crosses an
intermediate threshold. This behavior arises as a consequence of the existence
of alternative unfolding routes that are dominant at different force ranges.
Responsible for longer unfolding times at higher forces is the existence of an
intermediate, metastable configuration where the slipknot is jammed.
Simulations are performed with a coarsed grained model with further
quantification using a refined description of the geometry of the slipknots.
The simulation data is used to determine the free energy landscape (FEL) of the
protein, which supports recent analytical predictions.


  We present a generic mechanism by which reproducing microorganisms, with a
diffusivity that depends on the local population density, can form stable
patterns. It is known that a decrease of swimming speed with density can
promote separation into bulk phases of two coexisting densities; this is
opposed by the logistic law for birth and death which allows only a single
uniform density to be stable. The result of this contest is an arrested
nonequilibrium phase separation in which dense droplets or rings become
separated by less dense regions, with a characteristic steady-state length
scale. Cell division mainly occurs in the dilute regions and cell death in the
dense ones, with a continuous flux between these sustained by the diffusivity
gradient. We formulate a mathematical model of this in a case involving
run-and-tumble bacteria, and make connections with a wider class of mechanisms
for density-dependent motility. No chemotaxis is assumed in the model, yet it
predicts the formation of patterns strikingly similar to those believed to
result from chemotactic behavior.


  A new method for determining whether or not a mitrochondrial DNA (mtDNA)
sequence belongs to a vertebrate is described and tested. This method only
needs the mtDNA sequence of the organism in question, and unlike alignment
based methods, it does not require it to be compared with anything else. The
method is tested on all 1877 mtDNA sequences that were on NCBI's nucleotide
database on August 12, 2009, and works in 94.57% of the cases. Furthermore, all
organisms on which this method failed are closely related phylogenetically in
comparison to all other organisms included in the study. A list of potential
extensions to this method and open problems that emerge out of this study is
presented at the end.


  We consider optimal vaccination protocol where the vaccine is in short
supply. In this case, disease extinction results from a large and rare
fluctuation. We show that the probability of such fluctuation can be
exponentially increased by vaccination. For periodic vaccination with fixed
average rate, the optimal vaccination protocol is model independent and
presents a sequence of short pulses. The effect of vaccination can be
resonantly enhanced if the pulse period coincides with the characteristic
period of the disease dynamics or its multiples. This resonant effect is
illustrated using a simple epidemic model. If the system is periodically
modulated, the pulses must be synchronized with the modulation, whereas in the
case of a wrong phase the vaccination can lead to a negative result. The
analysis is based on the theory of fluctuation-induced population extinction in
periodically modulated systems that we develop.


  Dispersal is an important strategy that allows organisms to locate and
exploit favorable habitats. The question arises: given competition in a
spatially heterogeneous landscape, what is the optimal rate of dispersal?
Continuous population models predict that a species with a lower dispersal rate
always drives a competing species to extinction in the presence of spatial
variation of resources. However, the introduction of intrinsic demographic
stochasticity can reverse this conclusion. We present a simple model in which
competition between the exploitation of resources and stochastic fluctuations
leads to victory by either the faster or slower of two species depending on the
environmental parameters. A simplified limiting case of the model, analyzed by
closing the moment and correlation hierarchy, quantitatively predicts which
species will win in the complete model under given parameters of spatial
variation and average carrying capacity.


  The problem of motif detection can be formulated as the construction of a
discriminant function to separate sequences of a specific pattern from
background. In computational biology, motif detection is used to predict DNA
binding sites of a transcription factor (TF), mostly based on the weight matrix
(WM) model or the Gibbs free energy (FE) model. However, despite the wide
applications, theoretical analysis of these two models and their predictions is
still lacking. We derive asymptotic error rates of prediction procedures based
on these models under different data generation assumptions. This allows a
theoretical comparison between the WM-based and the FE-based predictions in
terms of asymptotic efficiency. Applications of the theoretical results are
demonstrated with empirical studies on ChIP-seq data and protein binding
microarray data. We find that, irrespective of underlying data generation
mechanisms, the FE approach shows higher or comparable predictive power
relative to the WM approach when the number of observed binding sites used for
constructing a discriminant decision is not too small.


  We show that minuscule entropic forces, on the order of 100 fN, can prevent
the formation of DNA loops--a ubiquitous means of regulating the expression of
genes. We observe a tenfold decrease in the rate of LacI-mediated DNA loop
formation when a tension of 200 fN is applied to the substrate DNA, biasing the
thermal fluctuations that drive loop formation and breakdown events.
Conversely, once looped, the DNA-protein complex is insensitive to applied
force. Our measurements are in excellent agreement with a simple polymer model
of loop formation in DNA, and show that an anti-parallel topology is the
preferred LacI-DNA loop conformation for a generic loop-forming construct.


  The celebrated Kermack-McKendric model of epidemics studies the transmission
of a disease in a population where each individual is initially susceptible
(S), may become infective (I) and then removed or recovered (R) and plays no
further epidemiological role. This ODE model arises as the limiting case of a
network model where each individual has an equal chance of infecting every
other. More recent work gives explicit consideration to the network of social
interaction and attendant probability of transmission for each interacting
pair. The state of such a network is an assignment of the values {S,I,R} to its
members. Given such a network, an initial state and a particular susceptible
individual, we would like to compute their probability of becoming infected in
the course of an epidemic. It turns out that this problem is NP-hard. In
particular, it belongs in a class of problems all of whose known solutions
require an exponential amount of computation and for which it is unlikely that
there will be more efficient solutions.