  We study equilibrium configurations of swarming biological organisms subject
to exogenous and pairwise endogenous forces. Beginning with a discrete
dynamical model, we derive a variational description of the corresponding
continuum population density. Equilibrium solutions are extrema of an energy
functional, and satisfy a Fredholm integral equation. We find conditions for
the extrema to be local minimizers, global minimizers, and minimizers with
respect to infinitesimal Lagrangian displacements of mass. In one spatial
dimension, for a variety of exogenous forces, endogenous forces, and domain
configurations, we find exact analytical expressions for the equilibria. These
agree closely with numerical simulations of the underlying discrete model.The
exact solutions provide a sampling of the wide variety of equilibrium
configurations possible within our general swarm modeling framework. The
equilibria typically are compactly supported and may contain
$\delta$-concentrations or jump discontinuities at the edge of the support. We
apply our methods to a model of locust swarms, which are observed in nature to
consist of a concentrated population on the ground separated from an airborne
group. Our model can reproduce this configuration; quasi-two-dimensionality of
the model plays a critical role.


  We present a model for the actin contractile ring of adherent animal cells.
The model suggests that the actin concentration within the ring and
consequently the power that the ring exerts both increase during contraction.
We demonstrate the crucial role of actin polymerization and depolymerization
throughout cytokinesis, and the dominance of viscous dissipation in the
dynamics. The physical origin of two phases in cytokinesis dynamics ("biphasic
cytokinesis") follows from a limitation on the actin density. The model is
consistent with a wide range of measurements of the midzone of dividing animal
cells.


  Evolutionary $2 \times 2$ games are studied with players located on a square
lattice. During the evolution the randomly chosen neighboring players try to
maximize their collective income by adopting a random strategy pair with a
probability dependent on the difference of their summed payoffs between the
final and initial state assuming quenched strategies in their neighborhood. In
the case of the anti-coordination game this system behaves alike an
anti-ferromagnetic kinetic Ising model. Within a wide region of social dilemmas
this dynamical rule supports the formation of similar spatial arrangement of
the cooperators and defectors ensuring the optimum total payoff if the
temptation to choose defection exceeds a threshold value dependent on the
sucker's payoff. The comparison of the results with those achieved for pairwise
imitation and myopic strategy updates has indicated the relevant advantage of
pairwise collective strategy update in the maintenance of cooperation.


  Actin cytoskeletal protrusions in crawling cells, or lamellipodia, exhibit
various morphological properties such as two characteristic peaks in the
distribution of filament orientation with respect to the leading edge. To
understand these properties, using the dendritic nucleation model as a basis
for cytoskeletal restructuring, a kinetic-population model with
orientational-dependent branching (birth) and capping (death) is constructed
and analyzed. Optimizing for growth yields a relation between the branch angle
and filament orientation that explains the two characteristic peaks. The model
also exhibits a subdominant population that allows for more accurate modeling
of recent measurements of filamentous actin density along the leading edge of
lamellipodia in keratocytes. Finally, we explore the relationship between
orientational and spatial organization of filamentous actin in lamellipodia and
address recent observations of a prevalence of overlapping filaments to
branched filaments---a finding that is claimed to be in contradiction with the
dendritic nucleation model.


  We use computer simulation to study crystal-forming model proteins equipped
with interactions that are both orientationally specific and nonspecific.
Distinct dynamical pathways of crystal formation can be selected by tuning the
strengths of these interactions. When the nonspecific interaction is strong,
liquidlike clustering can precede crystallization; when it is weak, growth can
proceed via ordered nuclei. Crystal yields are in certain parameter regimes
enhanced by the nonspecific interaction, even though it promotes association
without local crystalline order. Our results suggest that equipping nanoscale
components with weak nonspecific interactions (such as depletion attractions)
can alter both their dynamical pathway of assembly and optimize the yield of
the resulting material.


  Theoretical ecologists have long sought to understand how the persistence of
populations depends on biotic and abiotic factors. Classical work showed that
demographic stochasticity causes the mean time to extinction to increase
exponentially with population size, whereas variation in environmental
conditions can lead to a power law scaling. Recent work has focused especially
on the influence of the autocorrelation structure ("color") of environmental
noise. In theoretical physics, there is a burst of research activity in
analyzing large fluctuations in stochastic population dynamics. This research
provides powerful tools for determining extinction times and characterizing the
pathway to extinction. It yields, therefore, sharp insights into extinction
processes and has great potential for further applications in theoretical
biology.


  Virus entry is a multistep process that triggers a variety of cellular
pathways interconnecting into a complex network, yet the molecular complexity
of this network remains largely unsolved. Here, by employing systems biology
approach, we reveal a systemic virus-entry network initiated by human
cytomegalovirus (HCMV), a widespread opportunistic pathogen. This network
contains all known interactions and functional modules (i.e. groups of
proteins) coordinately responding to HCMV entry. The number of both genes and
functional modules activated in this network dramatically declines shortly,
within 25 min post-infection. While modules annotated as receptor system, ion
transport, and immune response are continuously activated during the entire
process of HCMV entry, those for cell adhesion and skeletal movement are
specifically activated during viral early attachment, and those for immune
response during virus entry. HCMV entry requires a complex receptor network
involving different cellular components, comprising not only cell surface
receptors, but also pathway components in signal transduction, skeletal
development, immune response, endocytosis, ion transport, macromolecule
metabolism and chromatin remodeling. Interestingly, genes that function in
chromatin remodeling are the most abundant in this receptor system, suggesting
that global modulation of transcriptions is one of the most important events in
HCMV entry. Results of in silico knock out further reveal that this entire
receptor network is primarily controlled by multiple elements, such as EGFR
(Epidermal Growth Factor) and SLC10A1 (sodium/bile acid cotransporter family,
member 1). Thus, our results demonstrate that a complex systemic network, in
which components coordinating efficiently in time and space contributes to
virus entry.


  It is shown that the nucleotide sequences in DNA molecules have
cluster-scaling properties (discovered for the first time in turbulent
processes: Sreenivasan and Bershadskii, 2006, J. Stat. Phys., 125, 1141-1153.).
These properties are relevant to both types of nucleotide pair-bases
interactions: hydrogen bonds and stacking interactions. It is shown that taking
into account the cluster-scaling properties can help to improve heterogeneous
models of the DNA dynamics. Two human genes: BRCA2 and NRXN1, have been
considered as examples.


  This paper discusses the role of opportunistic punisher who may act selfishly
to free-ride cooperators or not to be exploited by defectors. To consider
opportunistic punisher, we make a change to the sequence of one-shot public
good game; instead of putting action choice first before punishment, the
commitment of punishment is declared first before choosing the action of each
participant. In this commitment-first setting, punisher may use information
about her team, and may defect to increase her fitness in the team. Reversing
sequence of public good game can induce different behavior of punisher, which
cannot be considered in standard setting where punisher always chooses
cooperation. Based on stochastic dynamics developed by evolutionary economists
and biologists, we show that opportunistic punisher can make cooperation evolve
where cooperative punisher fails. This alternative route for the evolution of
cooperation relies paradoxically on the players' selfishness to profit from
others' unconditional cooperation and defection.


  Stoichiometric constraints play a role in the dynamics of natural
populations, but are not explicitly considered in most mathematical models.
Recent theoretical works suggest that these constraints can have a significant
impact and should not be neglected. However, it is not yet resolved how
stoichiometry should be integrated in population dynamical models, as different
modeling approaches are found to yield qualitatively different results. Here we
investigate a unifying framework that reveals the differences and commonalities
between previously proposed models for producer-grazer systems. Our analysis
reveals that stoichiometric constraints affect the dynamics mainly by
increasing the intraspecific competition between producers and by introducing a
variable biomass conversion efficiency. The intraspecific competition has a
strongly stabilizing effect on the system, whereas the variable conversion
efficiency resulting from a variable food quality is the main determinant for
the nature of the instability once destabilization occurs. Only if the food
quality is high an oscillatory instability, as in the classical paradox of
enrichment, can occur. While the generalized model reveals that the generic
insights remain valid in a large class of models, we show that other details
such as the specific sequence of bifurcations encountered in enrichment
scenarios can depend sensitively on assumptions made in modeling stoichiometric
constraints.


  According to recent experimental evidence, the architecture of a promoter,
defined as the number, strength and regulatory role of the operators that
control the promoter, plays a major role in determining the level of
cell-to-cell variability in gene expression. These quantitative experiments
call for a corresponding modeling effort that addresses the question of how
changes in promoter architecture affect noise in gene expression in a
systematic rather than case-by-case fashion. In this article, we make such a
systematic investigation, based on a simple microscopic model of gene
regulation that incorporates stochastic effects. In particular, we show how
operator strength and operator multiplicity affect this variability. We examine
different modes of transcription factor binding to complex promoters
(cooperative, independent, simultaneous) and how each of these affects the
level of variability in transcription product from cell-to-cell. We propose
that direct comparison between in vivo single-cell experiments and theoretical
predictions for the moments of the probability distribution of mRNA number per
cell can discriminate between different kinetic models of gene regulation.


  Found in all eukaryotic cells, linker histones H1 are known to bind to and
rearrange nucleosomal linker DNA. In vitro, the fundamental nature of H1/DNA
interactions has attracted wide interest among research communities - for
biologists from a chromatin organization deciphering point of view, and for
physicists from the study of polyelectrolyte interactions point of view. Hence,
H1/DNA binding processes, structural and dynamical information about these
self-assemblies is of broad importance. Targeting a quantitative understanding
of H1 induced DNA compaction mechanisms our strategy is based on using small
angle X-ray microdiffraction in combination with microfluidics. The usage of
microfluidic hydrodynamic focusing devices facilitate a microscale control of
these self-assembly processes. In addition, the method enables time-resolved
access to structure formation in situ, in particular to transient intermediate
states. The observed time dependent structure evolution shows that the
interaction of H1 with DNA can be described as a two step process: an initial
unspecific binding of H1 to DNA is followed by a rearrangement of molecules
within the formed assemblies. The second step is most likely induced by
interactions between the charged side chains of the protein and DNA. This leads
to an increase in lattice spacing within the DNA/protein assembly and induces a
decrease in the correlation length of the mesophases, probably due to a local
bending of the DNA.


  We construct a nonstandard finite difference (NSFD) scheme for an SIRS
mathematical model of respiratory virus transmission. This discretization is in
full compliance with the NSFD methodology as formulated by R. E. Mickens. By
use of an exact conservation law satisfied by the SIRS differential equations,
we are able to determine the corresponding denominator function for the
discrete first-order time derivatives. Our scheme is dynamically consistent
with the SIRS differential equations since the conservation laws are preserved.
Further, the scheme is shown to satisfy a positivity condition for its
solutions for all values of the time step-size.


  Renewal-anomalous-heterogeneous files are solved. A simple file is made of
Brownian hard spheres that diffuse stochastically in an effective 1D channel.
Generally, Brownian files are heterogeneous: the spheres' diffusion
coefficients are distributed and the initial spheres' density is non-uniform.
In renewal-anomalous files, the distribution of waiting times for individual
jumps is exponential as in Brownian files, yet obeys: {\psi}_{\alpha}
(t)~t^(-1-{\alpha}), 0<{\alpha}<1. The file is renewal as all the particles
attempt to jump at the same time. It is shown that the mean square displacement
(MSD) in a renewal-anomalous-heterogeneous file, <r^2>, obeys,
<r^2>~[<r^2>_{nrml}]^{\alpha}, where <r^2 >_{nrml} is the MSD in the
corresponding Brownian file. This scaling is an outcome of an exact relation
(derived here) connecting probability density functions of Brownian files and
renewal-anomalous files. It is also shown that non-renewal-anomalous files are
slower than the corresponding renewal ones.


  Rotation of a single bacterial flagellar motor is powered by multiple stators
tethered to the cell wall. In a "power-stroke" model the observed independence
of the speed at low load on the number of stators is explained by a
torque-dependent stepping mechanism independent of the strength of the stator
tethering spring. On the other hand, in models that depend solely on the stator
spring to explain the observed behavior, exceedingly small stator spring
constants are required. To study the dynamics of the motor driven by external
forces (such as those exerted by an optical tweezer), back-stepping is
introduced when stators are driven far out of equilibrium. Our model with
back-stepping reproduces the observed absence of a barrier to backward
rotation, as well the behaviors in the high-speed negative-torque regime.
Recently measured temperature dependence of the motor speed near zero load
(Yuan & Berg 2010 Biophys J) is explained quantitatively by the thermally
activated stepping rates in our model. Finally, we suggest that the general
mechanical properties of all molecular motors (linear and rotary),
characterized by their force(torque)-speed curve, can be determined by their
power-stroke potentials and the dependence of the stepping rates on the
mechanical state of the motor (force or speed). The torque-speed curve for the
clockwise rotating flagellar motor has been observed for the first time
recently (Yuan et al. 2010 PNAS). Its quasi-linear behavior is quantitatively
reproduced by our model. In particular, we show that concave and convex shapes
of the torque-speed curve can be achieved by changing the interaction potential
from linear to quadratic form. We also show that reversing the stepping rate
dependence on force (torque) can lead to non-monotonicity in the speed-load
dependency.


  We understand the dynamics of the world around us as by associating pairs of
events, where one event has some influence on the other. These pairs of events
can be aggregated into a web of memories representing our understanding of an
episode of history. The events and the associations between them need not be
directly experienced-they can also be acquired by communication. In this paper
we take a network approach to study the dynamics of memories of history. First
we investigate the network structure of a data set consisting of reported
events by several individuals and how associations connect them. We focus our
measurement on degree distributions, degree correlations, cycles (which
represent inconsistencies as they would break the time ordering) and community
structure. We proceed to model effects of communication using an agent-based
model. We investigate the conditions for the memory webs of different
individuals to converge to collective memories, how groups where the
individuals have similar memories (but different from other groups) can form.
Our work outlines how the cognitive representation of memories and social
structure can co-evolve as a contagious process. We generate some testable
hypotheses including that the number of groups is limited as a function of the
total population size.


  We construct a coarse-grained model of parallel actin bundles crosslinked by
compact, globular bundling proteins, such as fascin and espin, necessary
components of filapodial and mechanosensory bundles. Consistent with structural
observations of bundles, we find that the optimal geometry for crosslinking is
overtwisted, requiring a coherent structural change of the helical geometry of
the filaments. We study the linker-dependent thermodynamic transition of
bundled actin filaments from their native state to the overtwisted state and
map out the "twist-state'' phase diagram in terms of the availability as well
as the flexibility of crosslinker proteins. We predict that the transition from
the uncrosslinked to fully-crosslinked state is highly sensitive to linker
flexibility: flexible crosslinking smoothly distorts the twist-state of bundled
filaments, while rigidly crosslinked bundles undergo a phase transition,
rapidly overtwisting filaments over a narrow range of free crosslinker
concentrations. Additionally, we predict a rich spectrum of intermediate
structures, composed of alternating domains of sparsely-bound (untwisted) and
strongly-bound (overtwisted) filaments. This model reveals that subtle
differences in crosslinking agents themselves modify not only the detailed
structure of parallel actin bundles, but also the thermodynamic pathway by
which they form.


  The ecological principle of limiting similarity dictates that species similar
in resource requirements will compete, with the superior eventually excluding
the inferior competitor from the community. The observation that nonetheless
apparently similar species comprise a significant proportion of the diversity
in any given community has led to suggestions that competition may not in fact
be an important regulator of community structure and assembly. Here we apply a
recently introduced metric of species interaction, fractional (relative)
abundance, to tree species of the tropical wet forest of Barro Colorado Island,
Panam\'a, the particular community that inspired the original model of
non-niche or 'neutral' community dynamics9. We show a distribution of
fractional abundances between pairs of most closely related congeneric tree
species differing from that expected of competitive exclusion, but also
inconsistent with expectations of simple similarity, whether such species
interchangeability (a fundamental requirement of neutrality) is inferred at the
community or the pair level. Similar evidence from a strikingly different dry
forest has been linked to the focused, stable competition of a temporal niche
dynamic. Taken together with these earlier findings, the results reported here
establish a potentially widespread and important role for species interaction
in the diversity and maintenance of natural communities that must be considered
when inferring process from pattern.


  We discuss explicit solutions of the logistic model with variable parameters.
Classical data on the sunflower seeds growth are revisited as a simple
application of the logistic model with periodic coefficients. Some applications
to related biological systems are briefly reviewed.


  Motivation: Second generation sequencing technology makes it feasible for
many researches to obtain enough sequence reads to attempt the de novo assembly
of higher eukaryotes (including mammals). De novo assembly not only provides a
tool for understanding wide scale biological variation, but within human
bio-medicine, it offers a direct way of observing both large scale structural
variation and fine scale sequence variation. Unfortunately, improvements in the
computational feasibility for de novo assembly have not matched the
improvements in the gathering of sequence data. This is for two reasons: the
inherent computational complexity of the problem, and the in-practice memory
requirements of tools.
  Results: In this paper we use entropy compressed or succinct data structures
to create a practical representation of the de Bruijn assembly graph, which
requires at least a factor of 10 less storage than the kinds of structures used
by deployed methods. In particular we show that when stored succinctly, the de
Bruijn assembly graph for homo sapiens requires only 23 gigabytes of storage.
Moreover, because our representation is entropy compressed, in the presence of
sequencing errors it has better scaling behaviour asymptotically than
conventional approaches.


  Two-dimensional gel electrophoresis has been instrumental in the birth and
developments of proteomics, although it is no longer the exclusive separation
tool used in the field of proteomics. In this review, a historical perspective
is made, starting from the days where two-dimensional gels were used and the
word proteomics did not even exist. The events that have led to the birth of
proteomics are also recalled, ending with a description of the now well-known
limitations of two-dimensional gels in proteomics. However, the
often-underestimated advantages of two-dimensional gels are also underlined,
leading to a description of how and when to use two-dimensional gels for the
best in a proteomics approach. Taking support of these advantages (robustness,
resolution, and ability to separate entire, intact proteins), possible future
applications of this technique in proteomics are also mentioned.


  Clustered regularly interspaced short palindromic repeats (CRISPR) in
bacterial and archaeal DNA have recently been shown to be a new type of
anti-viral immune system in these organisms. We here study the diversity of
spacers in CRISPR under selective pressure. We propose a population dynamics
model that explains the biological observation that the leader-proximal end of
CRISPR is more diversified and the leader-distal end of CRISPR is more
conserved. This result is shown to be in agreement with recent experiments. Our
results show thatthe CRISPR spacer structure is influenced by and provides a
record of the viral challenges that bacteria face.


  Hidden Markov Models (HMMs) are a commonly used tool for inference of
transcription factor (TF) binding sites from DNA sequence data. We exploit the
mathematical equivalence between HMMs for TF binding and the "inverse"
statistical mechanics of hard rods in a one-dimensional disordered potential to
investigate learning in HMMs. We derive analytic expressions for the Fisher
information, a commonly employed measure of confidence in learned parameters,
in the biologically relevant limit where the density of binding sites is low.
We then use techniques from statistical mechanics to derive a scaling principle
relating the specificity (binding energy) of a TF to the minimum amount of
training data necessary to learn it.


  Reliability on complex biological networks reconstructions remains a concern.
Although observations are getting more and more precise, the data collection
process is yet error prone and the proofs display uneven certitude. In the case
of metabolic networks, the currently employed confidence scoring system rates
reactions according to a discretized small set of labels denoting different
levels of experimental evidence or model-based likelihood. Here, we propose a
computational network-based system of reaction scoring that exploits the
complex hierarchical structure and the statistical regularities of the
metabolic network as a bipartite graph. We use the example of Escherichia coli
metabolism to illustrate our methodology. Our model is adjusted to the
observations in order to derive connection probabilities between individual
metabolite-reaction pairs and, after validation, we integrate individual link
information to assess the reliability of each reaction in probabilistic terms.
This network-based scoring system breaks the degeneracy of currently employed
scores, enables further confirmation of modeling results, uncovers very
specific reactions that could be functionally or evolutionary important, and
identifies prominent experimental targets for further verification. We foresee
a wide range of potential applications of our approach given the natural
network bipartivity of many biological interactions.


  The standard Penna ageing model with sexual reproduction is enlarged by
adding additional bit-strings for love: Marriage happens only if the male love
strings are sufficiently different from the female ones. We simulate at what
level of required difference the population dies out.


  We present a methodology for modelling population dynamics with formal means
of computer science. This allows unambiguous description of systems and
application of analysis tools such as simulators and model checkers. In
particular, the dynamics of a population of Aedes albopictus (a species of
mosquito) and its modelling with the Stochastic Calculus of Looping Sequences
(Stochastic CLS) are considered. The use of Stochastic CLS to model population
dynamics requires an extension which allows environmental events (such as
changes in the temperature and rainfalls) to be taken into account. A simulator
for the constructed model is developed via translation into the specification
language Maude, and used to compare the dynamics obtained from the model with
real data.


  We define an individual-based probabilistic model of a sole (Solea solea)
behaviour. The individual model is given in terms of an Extended Probabilistic
Discrete Timed Automaton (EPDTA), a new formalism that is introduced in the
paper and that is shown to be interpretable as a Markov decision process. A
given EPDTA model can be probabilistically model-checked by giving a suitable
translation into syntax accepted by existing model-checkers. In order to
simulate the dynamics of a given population of soles in different environmental
scenarios, an agent-based simulation environment is defined in which each agent
implements the behaviour of the given EPDTA model. By varying the probabilities
and the characteristic functions embedded in the EPDTA model it is possible to
represent different scenarios and to tune the model itself by comparing the
results of the simulations with real data about the sole stock in the North
Adriatic sea, available from the recent project SoleMon. The simulator is
presented and made available for its adaptation to other species.


  Metapopulations are models of ecological systems, describing the interactions
and the behavior of populations that live in fragmented habitats. In this
paper, we present a model of metapopulations based on the multivolume
simulation algorithm tau-DPP, a stochastic class of membrane systems, that we
utilize to investigate the influence that different habitat topologies can have
on the local and global dynamics of metapopulations. In particular, we focus
our analysis on the migration rate of individuals among adjacent patches, and
on their capability of colonizing the empty patches in the habitat. We compare
the simulation results obtained for each habitat topology, and conclude the
paper with some proposals for other research issues concerning metapopulations.


  Existing sequence alignment algorithms use heuristic scoring schemes which
cannot be used as objective distance metrics. Therefore one relies on measures
like the p- or log-det distances, or makes explicit, and often simplistic,
assumptions about sequence evolution. Information theory provides an
alternative, in the form of mutual information (MI) which is, in principle, an
objective and model independent similarity measure. MI can be estimated by
concatenating and zipping sequences, yielding thereby the "normalized
compression distance". So far this has produced promising results, but with
uncontrolled errors. We describe a simple approach to get robust estimates of
MI from global pairwise alignments. Using standard alignment algorithms, this
gives for animal mitochondrial DNA estimates that are strikingly close to
estimates obtained from the alignment free methods mentioned above. Our main
result uses algorithmic (Kolmogorov) information theory, but we show that
similar results can also be obtained from Shannon theory. Due to the fact that
it is not additive, normalized compression distance is not an optimal metric
for phylogenetics, but we propose a simple modification that overcomes the
issue of additivity. We test several versions of our MI based distance measures
on a large number of randomly chosen quartets and demonstrate that they all
perform better than traditional measures like the Kimura or log-det (resp.
paralinear) distances. Even a simplified version based on single letter Shannon
entropies, which can be easily incorporated in existing software packages, gave
superior results throughout the entire animal kingdom. But we see the main
virtue of our approach in a more general way. For example, it can also help to
judge the relative merits of different alignment algorithms, by estimating the
significance of specific alignments.


  We present a computational method for finding attractors (ergodic sets of
states) of Boolean networks under asynchronous update. The approach is based on
a systematic removal of state transitions to render the state transition graph
acyclic. In this reduced state transition graph, all attractors are fixed
points that can be enumerated with little effort in most instances. This
attractor set is then extended to the attractor set of the original dynamics.
Our numerical tests on standard Kauffman networks indicate that the method is
efficient in the sense that the total number of state vectors visited grows
moderately with the number of states contained in attractors.


  We extend a model of community isolation in the d-dimensional lattice onto
the case with an imposed imbalance between birth rates of competing
communities. We give analytical and numerical evidences that in the asymmetric
two-specie model there exists a well defined value of the asymmetry parameter
when the emergence of the isolated (blocked) subgroups is the fastest, i.e. the
characteristic time tc is minimal. This critical value of the parameter depends
only on the lattice dimensionality and is independent from the system size.
Similar phenomenon was observed in the multi-specie case with a geometric
distribution of the birth rates. We also show that blocked subgroups in the
multi-specie case are absent or very rare when either there is a strictly
dominant specie that outnumbers the others or when there is a large diversity
of species. The number of blocked species of different kinds decreases with the
dimension of the multi-specie system.


  Understanding protein structure is of crucial importance in science, medicine
and biotechnology. For about two decades, knowledge based potentials based on
pairwise distances -- so-called "potentials of mean force" (PMFs) -- have been
center stage in the prediction and design of protein structure and the
simulation of protein folding. However, the validity, scope and limitations of
these potentials are still vigorously debated and disputed, and the optimal
choice of the reference state -- a necessary component of these potentials --
is an unsolved problem. PMFs are loosely justified by analogy to the reversible
work theorem in statistical physics, or by a statistical argument based on a
likelihood function. Both justifications are insightful but leave many
questions unanswered. Here, we show for the first time that PMFs can be seen as
approximations to quantities that do have a rigorous probabilistic
justification: they naturally arise when probability distributions over
different features of proteins need to be combined. We call these quantities
reference ratio distributions deriving from the application of the reference
ratio method. This new view is not only of theoretical relevance, but leads to
many insights that are of direct practical use: the reference state is uniquely
defined and does not require external physical insights; the approach can be
generalized beyond pairwise distances to arbitrary features of protein
structure; and it becomes clear for which purposes the use of these quantities
is justified. We illustrate these insights with two applications, involving the
radius of gyration and hydrogen bonding. In the latter case, we also show how
the reference ratio method can be iteratively applied to sculpt an energy
funnel. Our results considerably increase the understanding and scope of energy
functions derived from known biomolecular structures.


  Motivation: Signaling pathways control a large variety of cellular processes.
However, currently, even within the same database signaling pathways are often
curated at different levels of detail. This makes comparative and cross-talk
analyses difficult. Results: We present SignaLink, a database containing 8
major signaling pathways from Caenorhabditis elegans, Drosophila melanogaster,
and humans. Based on 170 review and approx. 800 research articles, we have
compiled pathways with semi-automatic searches and uniform, well-documented
curation rules. We found that in humans any two of the 8 pathways can
cross-talk. We quantified the possible tissue- and cancer-specific activity of
cross-talks and found pathway-specific expression profiles. In addition, we
identified 327 proteins relevant for drug target discovery. Conclusions: We
provide a novel resource for comparative and cross-talk analyses of signaling
pathways. The identified multi-pathway and tissue-specific cross-talks
contribute to the understanding of the signaling complexity in health and
disease and underscore its importance in network-based drug target selection.
Availability: http://SignaLink.org


  Quantum effects are mainly used for the determination of molecular shapes in
molecular biology, but quantum information theory may be a more useful tool to
understand the physics of life. Organic molecules and quantum
circuits/protocols can be considered as hardware and software of living systems
that are co-optimized during evolution. We try to model DNA replication in this
sense as a multi-body entanglement swapping with a reliable qubit
representation of the nucleotides. In our model molecular recognition of a
nucleotide triggers an intrabase entanglement corresponding to a superposition
state of different tautomer forms. Then, base pairing occurs by swapping
intrabase entanglements with interbase entanglements. We examine possible
realizations of quantum circuits to be used to obtain intrabase entanglement
and swapping protocols to be employed to obtain interbase entanglement.
Finally, we discuss possible ways for computational and experimental
verification of the model.


  We derive an invertible transform linking two widely used measures of species
diversity: phylogenetic diversity and the expected proportions of segregating
(non-constant) sites. We assume a bi-allelic, symmetric, finite site model of
substitution. Like the Hadamard transform of Hendy and Penny, the transform can
be expressed completely independent of the underlying phylogeny. Our results
bridge work on diversity from two quite distinct scientific communities.


  Interactions between bacteria and protozoa is an increasing area of interest,
however there are a few systems that allow extensive observation of the
interactions. We examined a surface system consisting of non nutrient agar with
a uniform bacterial lawn that extended over the agar surface, and a spatially
localised central population of amoebae. The amoeba fed on bacteria and
migrated over the plate. Automated image analysis techniques were used to
locate and count amoebae, cysts and bacteria coverage in a series of spatial
images. Most algorithms were based on intensity thresholding, or a modification
of this idea with probabilistic models. Our strategy was two tiered, we
performed an automated analysis for object classification and bacteria counting
followed by user intervention/reclassification using custom written Graphical
User Interfaces.


  Using molecular dynamics simulations, we study the effect of pressure on the
binding propensity of small RNAs by calculating the free energy barrier
corresponding to the looped out conformations of unmatched base, which
presumably acts as the binding sites for ligands. We find that the free energy
associated with base looping out increases monotonically as the pressure is
increased. Furthermore, we calculate the mean first passage time of
conformational looping out of the base bulge using the diffusion of reaction
coordinate associated with the base flipping on the underlying free energy
surface. We find that the mean first passage time associated with bulge looping
out increases slowly upon increasing pressures $P$ upto $2$~kbar but changes
dramatically for $P>2$~kbar. Finally, we discuss our results in the light of
the role of hydration shell of water around RNA.


  This paper deals with extinction of an isolated population caused by
intrinsic noise. We model the population dynamics in a "refuge" as a Markov
process which involves births and deaths on discrete lattice sites and random
migrations between neighboring sites. In extinction scenario I the zero
population size is a repelling fixed point of the on-site deterministic
dynamics. In extinction scenario II the zero population size is an attracting
fixed point, corresponding to what is known in ecology as Allee effect.
Assuming a large population size, we develop WKB (Wentzel-Kramers-Brillouin)
approximation to the master equation. The resulting Hamilton's equations encode
the most probable path of the population toward extinction and the mean time to
extinction. In the fast-migration limit these equations coincide, up to a
canonical transformation, with those obtained, in a different way, by Elgart
and Kamenev (2004). We classify possible regimes of population extinction with
and without an Allee effect and for different types of refuge and solve several
examples analytically and numerically. For a very strong Allee effect the
extinction problem can be mapped into the over-damped limit of theory of
homogeneous nucleation due to Langer (1969). In this regime, and for very long
systems, we predict an optimal refuge size that maximizes the mean time to
extinction.


  Inspired by the ubiquity of composite filamentous networks in nature we
investigate models of biopolymer networks that consist of interconnected floppy
and stiff filaments. Numerical simulations carried out in three dimensions
allow us to explore the microscopic partitioning of stresses and strains
between the stiff and floppy fractions c_s and c_f, and reveal a non-trivial
relationship between the mechanical behavior and the relative fraction of stiff
polymer: when there are few stiff polymers, non-percolated stiff ``inclusions``
are protected from large deformations by an encompassing floppy matrix, while
at higher fractions of stiff material the stiff network is independently
percolated and dominates the mechanical response.


  Often questions arise about old or extinct networks. What proteins interacted
in a long-extinct ancestor species of yeast? Who were the central players in
the Last.fm social network 3 years ago? Our ability to answer such questions
has been limited by the unavailability of past versions of networks. To
overcome these limitations, we propose several algorithms for reconstructing a
network's history of growth given only the network as it exists today and a
generative model by which the network is believed to have evolved. Our
likelihood-based method finds a probable previous state of the network by
reversing the forward growth model. This approach retains node identities so
that the history of individual nodes can be tracked. We apply these algorithms
to uncover older, non-extant biological and social networks believed to have
grown via several models, including duplication-mutation with complementarity,
forest fire, and preferential attachment. Through experiments on both synthetic
and real-world data, we find that our algorithms can estimate node arrival
times, identify anchor nodes from which new nodes copy links, and can reveal
significant features of networks that have long since disappeared.


  The physical limit with which a cell senses external ligand concentration
corresponds to the perfect absorber, where all ligand particles are absorbed
and overcounting of same ligand particles does not occur. Here we analyze how
the lateral diffusion of receptors on the cell membrane affects the accuracy of
sensing ligand concentration. Specifically, we connect our modeling to
neurotransmission in neural synapses where the diffusion of glutamate receptors
is already known to refresh synaptic connections. We find that receptor
diffusion indeed increases the accuracy of sensing for both the glutamate AMPA
and NDMA receptors, although the NMDA receptor is overall much noiser. We
propose that the difference in accuracy of sensing of the two receptors can be
linked to their different roles in neurotransmission. Specifically, the high
accuracy in sensing glutamate is essential for the AMPA receptor to start
membrane depolarization, while the NMDA receptor is believed to work in a
second stage as a coincidence detector, involved in long-term potentiation and
memory.


  We develop a new framework for the study of complex continuous time dynamical
systems based on viewing them as collections of interacting control modules.
This framework is inspired by and builds upon the groupoid formalism of
Golubitsky, Stewart and their collaborators. Our approach uses the tools and
--- more importantly ---the stance of category theory. This enables us to put
the groupoid formalism in a coordinate-free setting and to extend it from
ordinary differential equations to vector fields on manifolds. In particular,
we construct combinatorial models for categories of modular continuous time
dynamical systems. Each such model, as a category, is a fibration over an
appropriate category of labeled directed graphs. This makes precise the
relation between dynamical systems living on networks and the combinatorial
structure of the underlying directed graphs, allowing us to exploit the
relation in new and interesting ways.


  Post-genomic research deals with challenging problems in screening genomes of
organisms for particular functions or potential for being the targets of
genetic engineering for desirable biological features. 'Phenotyping' of wild
type and mutants is a time-consuming and costly effort by many individuals.
This article is a preliminary progress report in research on large-scale
automation of phenotyping steps (imaging, informatics and data analysis) needed
to study plant gene-proteins networks that influence growth and development of
plants. Our results undermine the significance of phenotypic traits that are
implicit in patterns of dynamics in plant root response to sudden changes of
its environmental conditions, such as sudden re-orientation of the root tip
against the gravity vector. Including dynamic features besides the common
morphological ones has paid off in design of robust and accurate machine
learning methods to automate a typical phenotyping scenario, i.e. to
distinguish the wild type from the mutants.


  We consider a general, neutral, dynamical model of biodiversity. Individuals
have i.i.d. lifetime durations, which are not necessarily exponentially
distributed, and each individual gives birth independently at constant rate
\lambda. We assume that types are clonally inherited. We consider two classes
of speciation models in this setting. In the immigration model, new individuals
of an entirely new species singly enter the population at constant rate \mu
(e.g., from the mainland into the island). In the mutation model, each
individual independently experiences point mutations in its germ line, at
constant rate \theta. We are interested in the species abundance distribution,
i.e., in the numbers, denoted I_n(k) in the immigration model and A_n(k) in the
mutation model, of species represented by k individuals, k=1,2,...,n, when
there are n individuals in the total population. In the immigration model, we
prove that the numbers (I_t(k);k\ge 1) of species represented by k individuals
at time t, are independent Poisson variables with parameters as in Fisher's
log-series. When conditioning on the total size of the population to equal n,
this results in species abundance distributions given by Ewens' sampling
formula. In particular, I_n(k) converges as n\to\infty to a Poisson r.v. with
mean \gamma /k, where \gamma:=\mu/\lambda. In the mutation model, as
n\to\infty, we obtain the almost sure convergence of n^{-1}A_n(k) to a
nonrandom explicit constant. In the case of a critical, linear birth--death
process, this constant is given by Fisher's log-series, namely n^{-1}A_n(k)
converges to \alpha^{k}/k, where \alpha :=\lambda/(\lambda+\theta). In both
models, the abundances of the most abundant species are briefly discussed.


  The positions of nucleosomes in eukaryotic genomes determine which parts of
the DNA sequence are readily accessible for regulatory proteins and which are
not. Genome-wide maps of nucleosome positions have revealed a salient pattern
around transcription start sites, involving a nucleosome-free region (NFR)
flanked by a pronounced periodic pattern in the average nucleosome density.
While the periodic pattern clearly reflects well-positioned nucleosomes, the
positioning mechanism is less clear. A recent experimental study by Mavrich et
al. argued that the pattern observed in S. cerevisiae is qualitatively
consistent with a `barrier nucleosome model', in which the oscillatory pattern
is created by the statistical positioning mechanism of Kornberg and Stryer. On
the other hand, there is clear evidence for intrinsic sequence preferences of
nucleosomes, and it is unclear to what extent these sequence preferences affect
the observed pattern. To test the barrier nucleosome model, we quantitatively
analyze yeast nucleosome positioning data both up- and downstream from NFRs.
Our analysis is based on the Tonks model of statistical physics which
quantifies the interplay between the excluded-volume interaction of nucleosomes
and their positional entropy. We find that although the typical patterns on the
two sides of the NFR are different, they are both quantitatively described by
the same physical model, with the same parameters, but different boundary
conditions. The inferred boundary conditions suggest that the first nucleosome
downstream from the NFR (the +1 nucleosome) is typically directly positioned
while the first nucleosome upstream is statistically positioned via a
nucleosome-repelling DNA region. These boundary conditions, which can be
locally encoded into the genome sequence, significantly shape the statistical
distribution of nucleosomes over a range of up to ~1000 bp to each side.


  We have combined a custom implementation of the fast multiple-time-stepping
LN integrator with parallel tempering to explore folding properties of small
peptides in implicit solvent on the time scale of microseconds. We applied this
algorithm to the synthetic {\beta}-hairpin trpzip2 and one of its sequence
variants W2W9. Each simulation consisted of over 12 {\mu}s of aggregated
virtual time. Several measures of folding behavior showed convergence, allowing
comparison with experimental equilibrium properties. Our simulations suggest
that the electrostatic interaction of tryptophan sidechains is responsible for
much of the stability of the native fold. We conclude that the ff99 force field
combined with ff96 {\phi} and {\psi} dihedral energies and implicit solvent can
reproduce plausible folding behavior in both trpzip2 and W2W9.


  Galled networks, directed acyclic graphs that model evolutionary histories
with reticulation cycles containing only tree nodes, have become very popular
due to both their biological significance and the existence of polynomial time
algorithms for their reconstruction. In this paper we prove that Nakhleh's $m$
measure is a metric for this class of phylogenetic networks and hence it can be
safely used to evaluate galled network reconstruction methods.


  We consider a neutral dynamical model of biological diversity, where
individuals live and reproduce independently. They have i.i.d. lifetime
durations (which are not necessarily exponentially distributed) and give birth
(singly) at constant rate b. Such a genealogical tree is usually called a
splitting tree, and the population counting process (N_t;t\ge 0) is a
homogeneous, binary Crump--Mode--Jagers process. We assume that individuals
independently experience mutations at constant rate \theta during their
lifetimes, under the infinite-alleles assumption: each mutation instantaneously
confers a brand new type, called allele, to its carrier. We are interested in
the allele frequency spectrum at time t, i.e., the number A(t) of distinct
alleles represented in the population at time t, and more specifically, the
numbers A(k,t) of alleles represented by k individuals at time t,
k=1,2,...,N_t. We mainly use two classes of tools: coalescent point processes
and branching processes counted by random characteristics. We provide explicit
formulae for the expectation of A(k,t) in a coalescent point process
conditional on population size, which apply to the special case of splitting
trees. We separately derive the a.s. limits of A(k,t)/N_t and of A(t)/N_t
thanks to random characteristics. Last, we separately compute the expected
homozygosity by applying a method characterizing the dynamics of the tree
distribution as the origination time of the tree moves back in time, in the
spirit of backward Kolmogorov equations.


  We present an elementary model of random size varying population given by a
stationary continuous state branching process. For this model we compute the
joint distribution of: the time to the most recent common ancestor, the size of
the current population and the size of the population just before the most
recent common ancestor (MRCA). In particular we show a natural mild bottleneck
effect as the size of the population just before the MRCA is stochastically
smaller than the size of the current population. We also compute the number of
old families which corresponds to the number of individuals involved in the
last coalescent event of the genealogical tree. By studying more precisely the
genealogical structure of the population, we get asymptotics for the number of
ancestors just before the current time. We give explicit computations in the
case of the quadratic branching mechanism. In this case, the size of the
population at the MRCA is, in mean, less by 1/3 than size of the current
population size. We also provide in this case the fluctuations for the
renormalized number of ancestors.


  Continuous-time linear birth-death-immigration (BDI) processes are frequently
used in ecology and epidemiology to model stochastic dynamics of the population
of interest. In clinical settings, multiple birth-death processes can describe
disease trajectories of individual patients, allowing for estimation of the
effects of individual covariates on the birth and death rates of the process.
Such estimation is usually accomplished by analyzing patient data collected at
unevenly spaced time points, referred to as panel data in the biostatistics
literature. Fitting linear BDI processes to panel data is a nontrivial
optimization problem because birth and death rates can be functions of many
parameters related to the covariates of interest. We propose a novel
expectation--maximization (EM) algorithm for fitting linear BDI models with
covariates to panel data. We derive a closed-form expression for the joint
generating function of some of the BDI process statistics and use this
generating function to reduce the E-step of the EM algorithm, as well as
calculation of the Fisher information, to one-dimensional integration. This
analytical technique yields a computationally efficient and robust optimization
algorithm that we implemented in an open-source R package. We apply our method
to DNA fingerprinting of Mycobacterium tuberculosis, the causative agent of
tuberculosis, to study intrapatient time evolution of IS6110 copy number, a
genetic marker frequently used during estimation of epidemiological clusters of
Mycobacterium tuberculosis infections. Our analysis reveals previously
undocumented differences in IS6110 birth-death rates among three major lineages
of Mycobacterium tuberculosis, which has important implications for
epidemiologists that use IS6110 for DNA fingerprinting of Mycobacterium
tuberculosis.


  The computer simulations of Shennan (2001) are complemented by assuming the
environment to change randomly. For moderate change rates, fitness optimisation
through evolution is still possible.


  We study the thermal and mechanical behavior of DNA denaturation in the frame
of the mesoscopic Peyrard- Bishop-Dauxois model with the inclusion of solvent
interaction. By analyzing the melting transition of a homogeneous A-T sequence,
we are able to set suitable values of the parameters of the model and study the
formation and stability of bubbles in the system. Then, we focus on the case of
the P5 promoter sequence and use the Principal Component Analysis of the
trajectories to extract the main information on the dynamical behavior of the
system. We find that this analysis method gives an excellent agreement with
previous biological results.


  One-dimensional arrays of nucleosomes (DNA-bound histone octamers separated
by stretches of linker DNA) fold into higher-order chromatin structures which
ultimately make up eukaryotic chromosomes. Chromatin structure formation leads
to 10-11 base pair (bp) discretization of linker lengths caused by the smaller
free energy cost of packaging nucleosomes into a regular chromatin fiber if
their rotational setting (defined by DNA helical twist) is conserved. We
describe nucleosome positions along the fiber using a thermodynamic model of
finite-size particles with effective two-body interactions, subject to an
arbitrary external potential. We infer both one-body and two-body energies from
readily available large-scale maps of nucleosome positions. We show that
two-body forces play a leading role in establishing well-known 10-11 bp
genome-wide periodicity of nucleosome occupancies. They also explain nucleosome
ordering over transcribed regions observed in both in vitro and in vivo
high-throughput experiments.


  We propose the heterojunction effect in the analysis of the fluorescence
mechanism of the firefly chromophore. Following this analysis, and with respect
to the HOMO-LUMO gap alignment between the chromophore's functional fragments,
three main heterojunction types (I, II, and I*) are identified. Time-dependent
density-functional theory optical absorption calculations for the firefly
chromophore show that the strongest excitation appears in the deprotonated
anion state of the keto form. This can be explained by its high HOMO-LUMO
overlap due to strong bio-heterojunction confinement. It is also found that the
nitrogen atom in the thiazolyl rings, due to its larger electronegativity,
plays a key role in the emission process, its importance growing when HOMO and
LUMO overlap at its location. This principle is applied to enhance the
chromophore's fluorescence efficiency and to guide the functionalization of
molecular optoelectronic devices.


  We show the presence of both a minimum and clear oscillations in the
frequency dependence of the translocation time of a polymer described as a
unidimensional Rouse chain driven by a spatially localized oscillating linear
potential. The observed oscillations of the mean translocation time arise from
the synchronization between the very mean translocation time and the period of
the external force. We have checked the robustness of the frequency value for
the minimum translocation time by changing the damping parameter, finding a
very simple relationship between this frequency and the correspondent
translocation time. The translocation time as a function of the polymer length
has been also evaluated, finding a precise $L^2$ scaling. Furthermore, the role
played by the thermal fluctuations described as a Gaussian uncorrelated noise
has been also investigated, and the analogies with the resonant activation
phenomenon are commented.


  We introduce a multi-species lattice gas model for motor protein driven
collective cargo transport on cellular filaments. We use this model to describe
and analyze the collective motion of interacting vesicle cargoes being carried
by oppositely directed molecular motors, moving on a single biofilament.
Building on a totally asymmetric exclusion process (TASEP) to characterize the
motion of the interacting cargoes, we allow for mass exchange with the
environment, input and output at filament boundaries and focus on the role of
interconversion rates and how they affect the directionality of the net cargo
transport. We quantify the effect of the various different competing processes
in terms of non-equilibrium phase diagrams. The interplay of interconversion
rates, which allow for flux reversal and evaporation/deposition processes
introduce qualitatively new features in the phase diagrams. We observe regimes
of three-phase coexistence, the possibility of phase re-entrance and a
significant flexibility in how the different phase boundaries shift in response
to changes in control parameters. The moving steady state solutions of this
model allows for different possibilities for the spatial distribution of cargo
vesicles, ranging from homogeneous distribution of vesicles to polarized
distributions, characterized by inhomogeneities or {\it shocks}. Current
reversals due to internal regulation emerge naturally within the framework of
this model. We believe this minimal model will clarify the understanding of
many features of collective vesicle transport, apart from serving as the basis
for building more exact quantitative models for vesicle transport relevant to
various {\it in-vivo} situations.


  Ions in water are the liquid of life. Life occurs almost entirely in 'salt
water'. Water itself (without ions) is lethal to animal cells and damaging for
most proteins. Water must contain the right ions in the right amounts if it is
to sustain life. Physical chemistry is the language of electrolyte solutions.
Physical chemistry and biology are intertwined. Physical chemists and
biologists come from different traditions that separated for several decades as
biologists described the molecules of life. Communication is not easy between a
fundamentally descriptive tradition and a fundamentally analytical one.
Biologists have now learned to study well defined systems with physical
techniques, of considerable interest to physical chemists. Physical chemists
are increasingly interested in spatially inhomogeneous systems with structures
on the atomic scale so common in biology. Physical chemists will find it
productive to work on well defined systems built by evolution to be reasonably
robust, with input output relations insensitive to environmental insults. This
article deals with properties of ion channels that in my view can be dealt with
by 'physics as usual', with much the same tools that physical chemists apply to
other systems. Indeed, I introduce and use a tool of physicists-a field theory
(and boundary conditions) based on an energy variational approach developed by
Chun Liu-not too widely used among physical chemists. My goal is to provide the
knowledge base, and identify the assumptions, that biologists use in studying
ion channels, avoiding jargon. Rather simple models of selectivity and
permeation in ion channels work quite well in important cases. Those physical
models and cases are the main focus of this review because they demonstrate the
strong essential link between the traditional treatments of ions in chemical
physics, and the biological function of ion channels.


  Indirect effects are powerful influences in ecosystems that may maintain
species diversity and alter apparent relationships between species in
surprising ways. Here, we applied Network Environ Analysis to 50
empirically-based trophic ecosystem models to test the hypothesis that indirect
flows dominate direct flows in ecosystem networks. Further, we used Monte Carlo
based perturbations to investigate the robustness of these results to potential
error in the underlying data. To explain our findings, we further investigated
the importance of the microbial food web in recycling energy-matter using
components of the Finn Cycling Index and analysis of Environ Centrality. We
found that indirect flows dominate direct flows in 37/50 (74.0%) models. This
increases to 31/35 (88.5%) models when we consider only models that have
cycling structure and a representation of the microbial food web. The
uncertainty analysis reveals that there is less error in the I/D values than
the $\pm$ 5% error introduced into the models, suggesting the results are
robust to uncertainty. Our results show that the microbial food web mediates a
substantial percentage of cycling in some systems (median = 30.2%), but its
role is highly variable in these models, in agreement with the literature. Our
results, combined with previous work, strongly suggest that indirect effects
are dominant components of activity in ecosystems.


  We present a minimal model of plasma membrane heterogeneity that combines
criticality with connectivity to cortical cytoskeleton. Our model is motivated
by recent observations of micron-sized critical fluctuations in plasma membrane
vesicles that are detached from their cortical cytoskeleton. We incorporate
criticality using a conserved order parameter Ising model coupled to a simple
actin cytoskeleton interacting through point-like pinning sites. Using this
minimal model, we recapitulate several experimental observations of plasma
membrane raft heterogeneity. Small (r~20nm) and dynamic fluctuations at
physiological temperatures arise from criticality. Including connectivity to
cortical cytoskeleton disrupts large fluctuations, prevents macroscopic phase
separation at low temperatures (T<=22{\deg}C), and provides a template for long
lived fluctuations at physiological temperature (T=37{\deg}C).
Cytoskeleton-stabilized fluctuations produce significant barriers to the
diffusion of some membrane components in a manner that is weakly dependent on
the number of pinning sites and strongly dependent on criticality. More
generally, we demonstrate that critical fluctuations provide a physical
mechanism to organize and spatially segregate membrane components by providing
channels for interaction over large distances.


  The longitudinal response of single semiflexible polymers to sudden changes
in externally applied forces is known to be controlled by the propagation and
relaxation of backbone tension. Under many experimental circumstances,
realized, e.g., in nano-fluidic devices or in polymeric networks or solutions,
these polymers are effectively confined in a channel- or tube-like geometry. By
means of heuristic scaling laws and rigorous analytical theory, we analyze the
tension dynamics of confined semiflexible polymers for various generic
experimental setups. It turns out that in contrast to the well-known linear
response, the influence of confinement on the non-linear dynamics can largely
be described as that of an effective prestress. We also study the free
relaxation of an initially confined chain, finding a surprising superlinear
t^(9/8) growth law for the change in end-to-end distance at short times.


  The first realization of a designed, rather than natural, biochemical filter
process is reported and analyzed as a promising network component for
increasing the complexity of biomolecular logic systems. Key challenge in
biochemical logic research has been achieving scalability for complex network
designs. Various logic gates have been realized, but a "toolbox" of analog
elements for interconnectivity and signal processing has remained elusive.
Filters are important as network elements that allow control of noise in signal
transmission and conversion. We report a versatile biochemical filtering
mechanism designed to have sigmoidal response in combination with
signal-conversion process. Horseradish peroxidase-catalyzed oxidation of
chromogenic electron donor by hydrogen peroxide, was altered by adding
ascorbate, allowing to selectively suppress the output signal, modifying the
response from convex to sigmoidal. A kinetic model was developed for evaluation
of the quality of filtering. The results offer improved capabilities for design
of scalable biomolecular information processing systems.


  Signaling networks are designed to sense an environmental stimulus and adapt
to it. We propose and study a minimal model of signaling network that can sense
and respond to external stimuli of varying strength in an adaptive manner. The
structure of this minimal network is derived based on some simple assumptions
on its differential response to external stimuli. We employ stochastic
differential equations and probability distributions obtained from stochastic
simulations to characterize differential signaling response in our minimal
network model. We show that the proposed minimal signaling network displays two
distinct types of response as the strength of the stimulus is decreased. The
signaling network has a deterministic part that undergoes rapid activation by a
strong stimulus in which case cell-to-cell fluctuations can be ignored. As the
strength of the stimulus decreases, the stochastic part of the network begins
dominating the signaling response where slow activation is observed with
characteristic large cell-to-cell stochastic variability. Interestingly, this
proposed stochastic signaling network can capture some of the essential
signaling behaviors of a complex apoptotic cell death signaling network that
has been studied through experiments and large-scale computer simulations. Thus
we claim that the proposed signaling network is an appropriate minimal model of
apoptosis signaling. Elucidating the fundamental design principles of complex
cellular signaling pathways such as apoptosis signaling remains a challenging
task. We demonstrate how our proposed minimal model can help elucidate the
effect of a specific apoptotic inhibitor Bcl-2 on apoptotic signaling in a
cell-type independent manner. We also discuss the implications of our study in
elucidating the adaptive strategy of cell death signaling pathways.


  Cancer cells are widely known to be protected from apoptosis, which is a
major hurdle to successful anti-cancer therapy. Over-expression of several
anti-apoptotic proteins, or mutations in pro-apoptotic factors, has been
recognized to confer such resistance. Development of new experimental
strategies, such as in silico modeling of biological pathways, can increase our
understanding of how abnormal regulation of apoptotic pathway in cancer cells
can lead to tumour chemoresistance. Monte Carlo simulations are in particular
well suited to study inherent variability, such as spatial heterogeneity and
cell-to-cell variations in signaling reactions. Using this approach, often in
combination with experimental validation of the computational model, we
observed that large cell-to-cell variability could explain the kinetics of
apoptosis, which depends on the type of pathway and the strength of stress
stimuli. Most importantly, Monte Carlo simulations of apoptotic signaling
provides unexpected insights into the mechanisms of fractional cell killing
induced by apoptosis-inducing agents, showing that not only variation in
protein levels, but also inherent stochastic variability in signaling
reactions, can lead to survival of a fraction of treated cancer cells.


  This paper addresses the statistical significance of structures in random
data: Given a set of vectors and a measure of mutual similarity, how likely
does a subset of these vectors form a cluster with enhanced similarity among
its elements? The computation of this cluster p-value for randomly distributed
vectors is mapped onto a well-defined problem of statistical mechanics. We
solve this problem analytically, establishing a connection between the physics
of quenched disorder and multiple testing statistics in clustering and related
problems. In an application to gene expression data, we find a remarkable link
between the statistical significance of a cluster and the functional
relationships between its genes.


  This is an early but comprehensive review of the PNP Poisson Nernst Planck
theory of ion channels. Extensive reference is made to the earlier literature.
The starting place for this theory of open channels is a theory of
electrodiffusion rather like that used previously to describe membranes. The
theory uses Poisson's equation to describe how charge on ions and the channel
protein creates electrical potential; it uses the Nernst-Planck equations to
describe migration and diffusion of ions in gradients of concentration and
electrical potential. Combined, these are also the "drift-diffusion equations"
of solid state physics, which are widely, if not universally used to describe
the flow of current and the behavior of semiconductors.


  Equilibrium sampling of biomolecules remains an unmet challenge after more
than 30 years of atomistic simulation. Efforts to enhance sampling capability,
which are reviewed here, range from the development of new algorithms to
parallelization to novel uses of hardware. Special focus is placed on
classifying algorithms -- most of which are underpinned by a few key ideas --
in order to understand their fundamental strengths and limitations. Although
algorithms have proliferated, progress resulting from novel hardware use
appears to be more clear-cut than from algorithms alone, partly due to the lack
of widely used sampling measures.


  Since gene regulatory systems contain sometimes only a small number of
molecules, these systems are not described well by macroscopic rate equations;
a master equation approach is needed for such cases. We develop an
approximation scheme for dealing with the stochasticity of the gene regulatory
systems. Using an effective interaction concept, original master equations can
be reduced to simpler master equations, which can be solved analytically. We
apply the approximation scheme to self-regulating systems with monomer or dimer
interactions, and a two-gene system with an exclusive switch. The approximation
scheme can recover bistability of the exclusive switch adequately.


  A selective sweep describes the reduction of diversity due to strong positive
selection. If the mutation rate to a selectively beneficial allele is
sufficiently high, Pennings and Hermisson (2006a) have shown, that it becomes
likely, that a selective sweep is caused by several individuals. Such an event
is called a soft sweep and the complementary event of a single origin of the
beneficial allele, the classical case, a hard sweep. We give analytical
expressions for the linkage disequilibrium (LD) between two neutral loci linked
to the selected locus, depending on the recurrent mutation to the beneficial
allele, measured by $D$ and $\hat{\sigma_D^2}$, a quantity introduced by Ohta
and Kimura (1969), and conclude that the LD-pattern of a soft sweep differs
substantially from that of a hard sweep due to haplotype structure. We compare
our results with simulations.


  The analysis of stress response systems in microorganisms can reveal
molecular strategies for regulatory control and adaptation. Here, we focus on
the Cad module, a subsystem of E. coli's response to acidic stress, which is
conditionally activated at low pH only when lysine is available. When
expressed, the Cad system counteracts the elevated H+ concentration by
converting lysine to cadaverine under the consumption of H+, and exporting
cadaverine in exchange for external lysine. Surprisingly, the cad operon
displays a transient response, even when the conditions for its induction
persist. To quantitatively characterize the regulation of the Cad module, we
have experimentally recorded and theoretically modeled the dynamics of
important system variables. We establish a quantitative model that adequately
describes and predicts the transient expression behavior for various initial
conditions. Our quantitative analysis of the Cad system supports a negative
feedback by external cadaverine as the origin of the transient response.
Furthermore, the analysis puts causal constraints on the precise mechanism of
signal transduction via the regulatory protein CadC.


  The arabinose utilization system of E. coli displays a stochastic "all or
nothing" response at intermediate levels of arabinose, where the population
divides into a fraction catabolizing the sugar at a high rate (ON state) and a
fraction not utilizing arabinose (OFF state). Here we study this decision
process in individual cells, focusing on the dynamics of the transition from
the OFF to the ON state. Using quantitative time-lapse microscopy, we determine
the time delay between inducer addition and fluorescence onset of a GFP
reporter. Through independent characterization of the GFP maturation process,
we can separate the lag time caused by the reporter from the intrinsic
activation time of the arabinose system. The resulting distribution of
intrinsic time delays scales inversely with the external arabinose
concentration, and is compatible with a simple stochastic model for arabinose
uptake. Our findings support the idea that the heterogeneous timing of gene
induction is causally related to a broad distribution of uptake proteins at the
time of sugar addition.


  Annexins are a family of proteins that bind to anionic phospholipid membranes
in a Ca2+-dependent manner. Annexin A2 forms heterotetramers (Anx A2t) with the
S100A10 (p11) protein dimer. The tetramer is capable of bridging phospholipid
membranes and it has been suggested to play a role in Ca2+-dependent exocytosis
and cell-cell adhesion of metastatic cells. Here, we employ x-ray reflectivity
measurements to resolve the conformation of Anx A2t upon Ca2+-dependent binding
to single supported lipid bilayers (SLBs) composed of different mixtures of
anionic (POPS) and neutral (POPC) phospholipids. Based on our results we
propose that Anx A2t binds in a side-by-side configuration, i.e., both Anx A2
monomers bind to the bilayer with the p11 dimer positioned on top. Furthermore,
we observe a strong decrease of lipid mobility upon binding of Anx A2t to SLBs
with varying POPS content. X-ray reflectivity measurements indicate that
binding of Anx A2t also increases the density of the SLB. Interestingly, in the
protein-facing leaflet of the SLB the lipid density is higher than in the
substrate-facing leaflet. This asymmetric densification of the lipid bilayer by
Anx A2t and Ca2+ might have important implications for the biochemical
mechanism of Anx A2t-induced endo- and exocytosis.


  Patterns of word use both reflect and influence a myriad of human activities
and interactions. Like other entities that are reproduced and evolve, words
rise or decline depending upon a complex interplay between {their intrinsic
properties and the environments in which they function}. Using Internet
discussion communities as model systems, we define the concept of a word niche
as the relationship between the word and the characteristic features of the
environments in which it is used. We develop a method to quantify two important
aspects of the size of the word niche: the range of individuals using the word
and the range of topics it is used to discuss. Controlling for word frequency,
we show that these aspects of the word niche are strong determinants of changes
in word frequency. Previous studies have already indicated that word frequency
itself is a correlate of word success at historical time scales. Our analysis
of changes in word frequencies over time reveals that the relative sizes of
word niches are far more important than word frequencies in the dynamics of the
entire vocabulary at shorter time scales, as the language adapts to new
concepts and social groupings. We also distinguish endogenous versus exogenous
factors as additional contributors to the fates of words, and demonstrate the
force of this distinction in the rise of novel words. Our results indicate that
short-term nonstationarity in word statistics is strongly driven by individual
proclivities, including inclinations to provide novel information and to
project a distinctive social identity.


  The role of thermal pressure fluctuations in the ejection of tightly packaged
DNA from protein capsid shells is discussed in a model calculation. At
equilibrium before ejection we assume the DNA is folded many times into a
bundle of parallel segments that forms an equilibrium conformation at minimum
free energy, which presses tightly against internal capsid walls. Using a
canonical ensemble at temperature T we calculate internal pressure fluctuations
against a slowly moving or static capsid mantle for an elastic continuum model
of the folded DNA bundle. It is found that fluctuating pressure on the capsid
mantle from thermal excitation of longitudinal acoustic vibrations in the
bundle may have root-mean-square values which are several tens of atmospheres
for typically small phage dimensions.


  Certain short polycations, such as TAT and oligoarginine, rapidly pass
through the plasma membranes of mammalian cells by a mechanism called
transduction, as well as by endocytosis and macropinocytosis. These
cell-penetrating peptides (CPPs) can carry with them cargos of 30 amino acids,
more than the nominal limit of 500 Da and enough to be therapeutic. An analysis
of the electrostatics of a charge outside the cell membrane and some recent
experiments suggest that transduction may proceed by molecular electroporation.
Ways to target diseased cells, rather than all cells, are discussed.


  Cell-penetrating peptides (CPPs) such as HIV's trans-activating
transcriptional activator (TAT) and polyarginine rapidly pass through the
plasma membranes of mammalian cells by an unknown mechanism called
transduction. They may be medically useful when fused to well-chosen chains of
fewer than about 35 amino acids. I offer a simple model of transduction in
which phosphatidylserines and CPPs effectively form two plates of a capacitor
with a voltage sufficient to cause the formation of transient pores
(electroporation). The model is consistent with experimental data on the
transduction of oligoarginine into mouse C2-C12 myoblasts and makes three
testable predictions.


  In recent years, single molecule force techniques have opened a new avenue to
decipher the folding landscapes of biopolymers by allowing us to watch and
manipulate the dynamics of individual proteins and nucleic acids. In single
molecule force experiments, quantitative analyses of measurements employing
sound theoretical models and molecular simulations play central role more than
any other field. With a brief description of basic theories for force mechanics
and molecular simulation technique using self-organized polymer (SOP) model,
this chapter will discuss various issues in single molecule force spectroscopy
(SMFS) experiments, which include pulling speed dependent unfolding pathway,
measurement of energy landscape roughness, the in uence of molecular handles in
optical tweezers on measurement and molecular motion, and folding dynamics of
biopolymers under force quench condition.


  Phylogenomic approaches developed thus far are either too time-consuming or
lack a solid evolutionary basis. Moreover, no phylogenomic approach is capable
of constructing a tree directly from unassembled raw sequencing data. A new
phylogenomic method, CO-phylum, is developed to alleviate these flaws.
CO-phylum can generate a high-resolution and highly accurate tree using
complete genome or unassembled sequencing data of close related organisms, in
addition, CO-phylum distance is almost linear with p-distance.


  Self-organization of proteins in space and time is of crucial importance for
the functioning of cellular processes. Often, this organization takes place in
the presence of strong random fluctuations due to the small number of molecules
involved. We report on stochastic switching of the Min-protein distributions
between the two cell halves in short Escherichia coli cells. A computational
model provides strong evidence that the macroscopic switching is rooted in
microscopic noise on the molecular scale. In longer bacteria, the switching
turns into regular oscillations that are required for positioning of the
division plane. As the pattern becomes more regular, cell-to-cell variability
also lessens, indicating cell length-dependent regulation of Min-protein
activity.


  It has been reported that the number of transcription factors encoded in
prokaryotic genomes scales approximately quadratically with their total number
of genes. We propose a conceptual explanation of this finding and illustrate it
using a simple model in which metabolic and regulatory networks of prokaryotes
are shaped by horizontal gene transfer of coregulated metabolic pathways.
Adapting to a new environmental condition monitored by a new transcription
factor (e.g., learning to use another nutrient) involves both acquiring new
enzymes and reusing some of the enzymes already encoded in the genome. As the
repertoire of enzymes of an organism (its toolbox) grows larger, it can reuse
its enzyme tools more often and thus needs to get fewer new ones to master each
new task. From this observation, it logically follows that the number of
functional tasks and their regulators increases faster than linearly with the
total number of genes encoding enzymes. Genomes can also shrink, e.g., because
of a loss of a nutrient from the environment, followed by deletion of its
regulator and all enzymes that become redundant. We propose several simple
models of network evolution elaborating on this toolbox argument and
reproducing the empirically observed quadratic scaling. The distribution of
lengths of pathway branches in our model agrees with that of the real-life
metabolic network of Escherichia coli. Thus, our model provides a qualitative
explanation for broad distributions of regulon sizes in prokaryotes.


  In prokaryotic genomes the number of transcriptional regulators is known to
quadratically scale with the total number of protein-coding genes. Toolbox
model was recently proposed to explain this scaling for metabolic enzymes and
their regulators. According to its rules the metabolic network of an organism
evolves by horizontal transfer of pathways from other species. These pathways
are part of a larger "universal" network formed by the union of all
species-specific networks. It remained to be understood, however, how the
topological properties of this universal network influence the scaling law of
functional content of genomes. In this study we answer this question by first
analyzing the scaling properties of the toolbox model on arbitrary tree-like
universal networks. We mathematically prove that the critical branching
topology, in which the average number of upstream neighbors of a node is equal
to one, is both necessary and sufficient for the quadratic scaling. Conversely,
the toolbox model on trees with exponentially expanding, supercritical topology
is characterized by the linear scaling with logarithmic corrections. We further
generalize our model to include reactions with multiple substrates/products as
well as branched or cyclic metabolic pathways. Unlike the original model the
new version employs evolutionary optimized pathways with the smallest number of
reactions necessary to achieve their metabolic tasks. Numerical simulations of
this most realistic model on the universal network from the KEGG database again
produced approximately quadratic scaling. Our results demonstrate why, in spite
of their "small-world" topology, real-life metabolic networks are characterized
by a broad distribution of pathway lengths and sizes of metabolic regulons in
regulatory networks.


  We explore in detail the structural, mechanical and thermodynamic properties
of a coarse-grained model of DNA similar to that introduced in Thomas E.
Ouldridge, Ard A. Louis, Jonathan P.K. Doye, Phys. Rev. Lett. 104 178101
(2010). Effective interactions are used to represent chain connectivity,
excluded volume, base stacking and hydrogen bonding, naturally reproducing a
range of DNA behaviour. We quantify the relation to experiment of the
thermodynamics of single-stranded stacking, duplex hybridization and hairpin
formation, as well as structural properties such as the persistence length of
single strands and duplexes, and the torsional and stretching stiffness of
double helices. We also explore the model's representation of more complex
motifs involving dangling ends, bulged bases and internal loops, and the effect
of stacking and fraying on the thermodynamics of the duplex formation
transition.


  In this and the associated article 'BioBlender: Fast and Efficient All Atom
Morphing of Proteins Using Blender Game Engine', by Zini et al., we present
BioBlender, a complete instrument for the elaboration of motion (Zini et al.)
and the visualization (here) of proteins and other macromolecules, using
instruments of computer graphics. The availability of protein structures
enables the study of their surfaces and surface properties such as
electrostatic potential (EP) and hydropathy (MLP), based on atomic
contribution. Recent advances in 3D animation and rendering software have not
yet been exploited for the representation of proteins and other biological
molecules in an intuitive, animated form. Taking advantage of an open-source,
3D animation and rendering software, Blender, we developed BioBlender, a
package dedicated to biological work: elaboration of proteins' motions with the
simultaneous visualization of chemical and physical features. EP and MLP are
calculated using physico-chemical programs and custom programs and scripts,
organized and accessed within BioBlender interface. A new visual code is
introduced for MLP visualization: a range of optical features that permits a
photorealistic rendering of its spatial distribution on the surface of the
protein. EP is represented as animated line particles that flow along field
lines proportional to the total charge of the protein. Our system permits EP
and MLP visualization of molecules and, in the case of moving proteins, the
continuous perception of these features, calculated for each intermediate
conformation. Using real world tactile/sight feelings, the nanoscale world of
proteins becomes more understandable, familiar to our everyday life, making it
easier to introduce "un-seen" phenomena (concepts) such as hydropathy or
charges.


  The naming game (NG) describes the agreement dynamics of a population of
agents that interact locally in a pairwise fashion, and in recent years
statistical physics tools and techniques have greatly contributed to shed light
on its rich phenomenology. Here we investigate in details the role played by
the way in which the two agents update their states after an interaction. We
show that slightly modifying the NG rules in terms of which agent performs the
update in given circumstances (i.e. after a success) can either alter
dramatically the overall dynamics or leave it qualitatively unchanged. We
understand analytically the first case by casting the model in the broader
framework of a generalized NG. As for the second case, on the other hand, we
note that the modified rule reproducing the main features of the usual NG
corresponds in fact to a simplification of it consisting in the elimination of
feedback between the agents. This allows us to introduce and study a very
natural broadcasting scheme on networks that can be potentially relevant for
different applications, such as the design and implementation of autonomous
sensor networks, as pointed out in the recent literature.


  In this and the associated article 'BioBlender: A Software for Intuitive
Representation of Surface Properties of Biomolecules', (Andrei et al) we
present BioBlender as a complete instrument for the elaboration of motion
(here) and the visualization (Andrei et al) of proteins and other
macromolecules, using instruments of computer graphics. A vast number of
protein (if not most) exert their function through some extent of motion.
Despite recent advances in higly performant methods, it is very difficult to
obtain direct information on conformational changes of molecules. However,
several systems exist that can shed some light on the variability of
conformations of a single peptide chain; among them, NMR methods provide
collections of a number of static 'shots' of a moving protein. Starting from
this data, and assuming that if a protein exists in more than 1 conformation it
must be able to transit between the different states, we have elaborated a
system that makes ample use of the computational power of 3D computer graphics
technology. Considering information of all (heavy) atoms, we use animation and
game engine of Blender to obtain transition states. The model we chose to
elaborate our system is Calmodulin, a protein favorite among structural and
dynamic studies due to its (relative) simplicity of structure and small
dimension. Using Calmodulin we show a procedure that enables the building of a
'navigation map' of NMR models, that can help in the identification of
movements. In the process, a number of intermediate conformations is generated,
all of which respond to strict bio-physical and bio-chemical criteria. The
BioBlender system is available for download from the website
www.bioblender.net, together with examples, tutorial and other useful material.


  Motivation :Reconstructing the topology of a gene regulatory network is one
of the key tasks in systems biology. Despite of the wide variety of proposed
methods, very little work has been dedicated to the assessment of their
stability properties. Here we present a methodical comparison of the
performance of a novel method (RegnANN) for gene network inference based on
multilayer perceptrons with three reference algorithms (ARACNE, CLR, KELLER),
focussing our analysis on the prediction variability induced by both the
network intrinsic structure and the available data.
  Results: The extensive evaluation on both synthetic data and a selection of
gene modules of "Escherichia coli" indicates that all the algorithms suffer of
instability and variability issues with regards to the reconstruction of the
topology of the network. This instability makes objectively very hard the task
of establishing which method performs best. Nevertheless, RegnANN shows MCC
scores that compare very favorably with all the other inference methods tested.
  Availability: The software for the RegnANN inference algorithm is distributed
under GPL3 and it is available at the corresponding author home page
(http://mpba.fbk.eu/grimaldi/regnann-supmat)


  The catastrophic senescence of the Pacific salmon is among the initial tests
used to validate the Penna aging model. Based on the mutation accumulation
theory, the sudden decrease in fitness following reproduction may be solely
attributed to the semelparity of the species. In this work, we report other
consequences of mutation accumulation. Contrary to earlier findings, such
dramatic manifestation of aging depends not only on the choice of breeding
strategy but also on the value of the reproduction age, R, and the mutation
threshold, T. Senescence is catastrophic when $T \leq R$. As the organism's
tolerance for harmful genetic mutations increases, the aging process becomes
more gradual. We observe senescence that is threshold dependent whenever T>R.
That is, the sudden drop in survival rate occurs at age equal to the mutation
threshold value.


  A discrete chemotactic predator-prey model is proposed in which the prey
secrets a diffusing chemical which is sensed by the predator and vice versa.
Two dynamical states corresponding to catching and escaping are identified and
it is shown that steady hunting is unstable. For the escape process, the
predator-prey distance is diffusive for short times but exhibits a transient
subdiffusive behavior which scales as a power law $t^{1/3}$ with time $t$ and
ultimately crosses over to diffusion again. This allows to classify the
motility and dynamics of various predatory bacteria and phagocytes. In
particular, there is a distinct region in the parameter space where they prove
to be infallible predators.


  In synthetic biology, gene regulatory circuits are often constructed by
combining smaller circuit components. Connections between components are
achieved by transcription factors acting on promoters. If the individual
components behave as true modules and certain module interface conditions are
satisfied, the function of the composite circuits can in principle be
predicted. In this paper, we investigate one of the interface conditions:
fan-out. We quantify the fan-out, a concept widely used in electric
engineering, to indicate the maximum number of the downstream inputs that an
upstream output transcription factor can regulate. We show that the fan-out is
closely related to retroactivity studied by Del Vecchio, et al. We propose an
efficient operational method for measuring the fan-out that can be applied to
various types of module interfaces. We also show that the fan-out can be
enhanced by self-inhibitory regulation on the output. We discuss the potential
role of the inhibitory regulations found in gene regulatory networks and
protein signal pathways. The proposed estimation method for fanout not only
provides an experimentally efficient way for quantifying the level of
modularity in gene regulatory circuits but also helps characterize and design
module interfaces, enabling the modular construction of gene circuits.


  We show how one may analytically compute the stationary density of the
distribution of molecular constituents in populations of cells in the presence
of noise arising from either bursting transcription or translation, or noise in
degradation rates arising from low numbers of molecules. We have compared our
results with an analysis of the same model systems (either inducible or
repressible operons) in the absence of any stochastic effects, and shown the
correspondence between behaviour in the deterministic system and the stochastic
analogs. We have identified key dimensionless parameters that control the
appearance of one or two steady states in the deterministic case, or unimodal
and bimodal densities in the stochastic systems, and detailed the analytic
requirements for the occurrence of different behaviours. This approach
provides, in some situations, an alternative to computationally intensive
stochastic simulations. Our results indicate that, within the context of the
simple models we have examined, bursting and degradation noise cannot be
distinguished analytically when present alone.


  In this paper, we suggested eleven protein targets to be used as possible
vaccines against Plasmodium falciparum causative agent of almost two to three
million deaths per year. A comprehensive analysis of protein target have been
selected from the small experimental fragment of antigen in the P. falciparum
genome, all of them common to the four stages of the parasite life cycle (i.e.,
sporozoites, merozoites, trophozoites and gametocytes). The potential vaccine
candidates should be analyzed in silico technique using various bioinformatics
tools. Finally, the possible protein target according to PlasmoDB gene ID are
PFC0975c, PFE0660c, PF08_0071, PF10_0084, PFI0180w, MAL13P1.56, PF14_0192,
PF13_0141, PF14_0425, PF13_0322, y PF14_0598.