  In this paper, we introduce an ASEP-like transport model for bidirectional
motion of particles on a multi-lane lattice. The model is motivated by {\em in
vivo} experiments on organelle motility along a microtubule (MT), consisting of
thirteen protofilaments, where particles are propelled by molecular motors
(dynein and kinesin). In the model, organelles (particles) can switch
directions of motion due to "tug-of-war" events between counteracting motors.
Collisions of particles on the same lane can be cleared by switching to
adjacent protofilaments (lane changes).
  We analyze transport properties of the model with no-flux boundary conditions
at one end of a MT ("plus-end" or tip). We show that the ability of lane
changes can affect the transport efficiency and the particle-direction change
rate obtained from experiments is close to optimal in order to achieve
efficient motor and organelle transport in a living cell. In particular, we
find a nonlinear scaling of the mean {\em tip size} (the number of particles
accumulated at the tip) with injection rate and an associated phase transition
leading to {\em pulsing states} characterized by periodic filling and emptying
of the system.


  The dynamics of strongly immunizing childhood infections is still not well
understood. Although reports of successful modeling of several incidence data
records can be found in the literature, the key determinants of the observed
temporal patterns have not been clearly identified. In particular, different
models of immunity waning and degree of protection applied to disease and
vaccine induced immunity have been debated in the literature on pertussis. Here
we study the effect of disease acquired immunity on the long term patterns of
pertussis prevalence. We compare five minimal models, all of which are
stochastic, seasonally forced, well-mixed models of infection based on
susceptible-infective-recovered dynamics in a closed population. These models
reflect different assumptions about the immune response of naive hosts, namely
total permanent immunity, immunity waning, immunity waning together with
immunity boosting, reinfection of recovered, and repeat infection after partial
immunity waning. The power spectra of the output prevalence time series
characterize the long term dynamics of the models. For epidemiological
parameters consistent with published data for pertussis, the power spectra show
quantitative and even qualitative differences that can be used to test their
assumptions by comparison with ensembles of several decades long
pre-vaccination data records. We illustrate this strategy on two publicly
available historical data sets.


  Motivated by the dynamics of cultural change and diversity, we generalize the
three-species constrained voter model on a complete graph introduced in [J.
Phys. A 37, 8479 (2004)]. In this opinion dynamics model, a population of size
N is composed of "leftists" and "rightists" that interact with "centrists": a
leftist and centrist can both become leftists with rate (1+q)/2 or centrists
with rate (1-q)/2 (and similarly for rightists and centrists), where q denotes
the bias towards extremism (q>0) or centrism (q<0). This system admits three
absorbing fixed points and a "polarization" line along which a frozen mixture
of leftists and rightists coexist. In the realm of Fokker-Planck equation, and
using a mapping onto a population genetics model, we compute the fixation
probability of ending in every absorbing state and the mean times for these
events. We therefore show, especially in the limit of weak bias and large
population size when |q|~1/N and N>>1, how fluctuations alter the mean field
predictions: polarization is likely when q>0, but there is always a finite
probability to reach a consensus; the opposite happens when q<0. Our findings
are corroborated by stochastic simulations.


  Nonequilibrium experiments of single biomolecules such as force-induced
unfolding reveal details about a few degrees of freedom of a complex system.
Molecular dynamics simulations can provide complementary information, but
exploration of the space of possible configurations is often hindered by large
barriers in phase space that separate metastable regions. To solve this
problem, enhanced sampling methods have been developed that divide a phase
space into regions and integrate trajectory segments in each region. These
methods boost the probability of passage over barriers, and facilitate
parallelization since integration of the trajectory segments does not require
communication, aside from their initialization and termination. Here we present
a parallel version of an enhanced sampling method suitable for systems driven
far from equilibrium: nonequilibrium umbrella sampling (NEUS). We apply this
method to a coarse-grained model of a 262-nucleotide RNA molecule that unfolds
and refolds in an explicit flow field modeled with stochastic rotation
dynamics. Using NEUS we are able to observe extremely rare unfolding events
that have mean first passage times as long as 1.4 s (3.4 E13 dynamics steps).
We examine the unfolding process for a range of flow rates of the medium, and
we describe two competing pathways in which different intramolecular contacts
are broken.


  We study the role that parallelism plays in time complexity of Winfree's
abstract Tile Assembly Model (aTAM), a model of molecular algorithmic
self-assembly. In the "hierarchical" aTAM, two assemblies, both consisting of
multiple tiles, are allowed to aggregate together, whereas in the "seeded"
aTAM, tiles attach one at a time to a growing assembly. Adleman, Cheng, Goel,
and Huang ("Running Time and Program Size for Self-Assembled Squares", STOC
2001) showed how to assemble an n x n square in O(n) time in the seeded aTAM
using O(log n / log log n) unique tile types, where both of these parameters
are optimal. They asked whether the hierarchical aTAM could allow a tile system
to use the ability to form large assemblies in parallel before they attach to
break the Omega(n) lower bound for assembly time. We show that there is a tile
system with the optimal O(log n / log log n) tile types that assembles an n x n
square using O(log^2 n) parallel "stages", which is close to the optimal
Omega(log n) stages, forming the final n x n square from four n/2 x n/2
squares, which are themselves recursively formed from n/4 x n/4 squares, etc.
However, despite this nearly maximal parallelism, the system requires
superlinear time to assemble the square. We extend the definition of *partial
order tile systems* studied by Adleman et al. in a natural way to hierarchical
assembly and show that no hierarchical partial order tile system can build any
shape with diameter N in less than time Omega(N), demonstrating that in this
case the hierarchical model affords no speedup whatsoever over the seeded
model. We strengthen the Omega(N) time lower bound for deterministic seeded
systems of Adleman et al. to nondeterministic seeded systems. Finally, we show
that for infinitely many n, a tile system can assemble an n x n' rectangle,
with n > n', in time O(n^{4/5} log n), breaking the linear-time lower bound.


  Atomic-accuracy structure prediction of macromolecules is a long-sought goal
of computational biophysics. Accurate modeling should be achievable by
optimizing a physically realistic energy function but is presently precluded by
incomplete sampling of a biopolymer's many degrees of freedom. We present
herein a working hypothesis, called the "stepwise ansatz", for recursively
constructing well-packed atomic-detail models in small steps, enumerating
several million conformations for each monomer and covering all build-up paths.
By implementing the strategy in Rosetta and making use of high-performance
computing, we provide first tests of this hypothesis on a benchmark of fifteen
RNA loop modeling problems drawn from riboswitches, ribozymes, and the
ribosome, including ten cases that were not solvable by prior knowledge based
modeling approaches. For each loop problem, this deterministic stepwise
assembly (SWA) method either reaches atomic accuracy or exposes flaws in
Rosetta's all-atom energy function, indicating the resolution of the
conformational sampling bottleneck. To our knowledge, SWA is the first
enumerative, ab initio build-up method to systematically outperform existing
Monte Carlo and knowledge-based methods for 3D structure prediction. As a
rigorous experimental test, we have applied SWA to a small RNA motif of
previously unknown structure, the C7.2 tetraloop/tetraloop-receptor, and
stringently tested this blind prediction with nucleotide-resolution structure
mapping data.


  Tiling arrays make possible a large scale exploration of the genome thanks to
probes which cover the whole genome with very high density until 2 000 000
probes. Biological questions usually addressed are either the expression
difference between two conditions or the detection of transcribed regions. In
this work we propose to consider simultaneously both questions as an
unsupervised classification problem by modeling the joint distribution of the
two conditions. In contrast to previous methods, we account for all available
information on the probes as well as biological knowledge like annotation and
spatial dependence between probes. Since probes are not biologically relevant
units we propose a classification rule for non-connected regions covered by
several probes. Applications to transcriptomic and ChIP-chip data of
Arabidopsis thaliana obtained with a NimbleGen tiling array highlight the
importance of a precise modeling and the region classification.


  Large-scale proteomic analysis is emerging as a powerful technique in biology
and relies heavily on data acquired by state-of-the-art mass spectrometers. As
with any other field in Systems Biology, computational tools are required to
deal with this ocean of data. iTRAQ (isobaric Tags for Relative and Absolute
quantification) is a technique that allows simultaneous quantification of
proteins from multiple samples. Although iTRAQ data gives useful insights to
the biologist, it is more complex to perform analysis and draw biological
conclusions because of its multi-plexed design. One such problem is to find
proteins that behave in a similar way (i.e. change in abundance) among various
time points since the temporal variations in the proteomics data reveal
important biological information. Distance based methods such as Euclidian
distance or Pearson coefficient, and clustering techniques such as k-mean etc,
are not able to take into account the temporal information of the series. In
this paper, we present an linear-time algorithm for clustering similar patterns
among various iTRAQ time course data irrespective of their absolute values. The
algorithm, referred to as Temporal Pattern Mining(TPM), maps the data from a
Cartesian plane to a discrete binary plane. After the mapping a dynamic
programming technique allows mining of similar data elements that are
temporally closer to each other. The proposed algorithm accurately clusters
iTRAQ data that are temporally closer to each other with more than 99%
accuracy. Experimental results for different problem sizes are analyzed in
terms of quality of clusters, execution time and scalability for large data
sets. An example from our proteomics data is provided at the end to demonstrate
the performance of the algorithm and its ability to cluster temporal series
irrespective of their distance from each other.


  Loops are essential secondary structure elements in folded DNA and RNA
molecules and proliferate close to the melting transition. Using a theory for
nucleic acid secondary structures that accounts for the logarithmic entropy c
ln m for a loop of length m, we study homopolymeric single-stranded nucleic
acid chains under external force and varying temperature. In the thermodynamic
limit of a long strand, the chain displays a phase transition between a low
temperature / low force compact (folded) structure and a high temperature /
high force molten (unfolded) structure. The influence of c on phase diagrams,
critical exponents, melting, and force extension curves is derived
analytically. For vanishing pulling force, only for the limited range of loop
exponents 2 < c < 2.479 a melting transition is possible; for c <= 2 the chain
is always in the folded phase and for 2.479 < c always in the unfolded phase. A
force induced melting transition with singular behavior is possible for all
loop exponents c < 2.479 and can be observed experimentally by single molecule
force spectroscopy. These findings have implications for the hybridization or
denaturation of double stranded nucleic acids. The Poland-Scheraga model for
nucleic acid duplex melting does not allow base pairing between nucleotides on
the same strand in denatured regions of the double strand. If the sequence
allows these intra-strand base pairs, we show that for a realistic loop
exponent c ~ 2.1 pronounced secondary structures appear inside the single
strands. This leads to a lower melting temperature of the duplex than predicted
by the Poland-Scheraga model. Further, these secondary structures renormalize
the effective loop exponent c^, which characterizes the weight of a denatured
region of the double strand, and thus affect universal aspects of the duplex
melting transition.


  We describe the structure of connected graphs with the minimum and maximum
average distance, radius, diameter, betweenness centrality, efficiency and
resistance distance, given their order and size. We find tight bounds on these
graph qualities for any arbitrary number of nodes and edges and analytically
derive the form and properties of such networks.


  We consider an asexual population under strong selection-weak mutation
conditions evolving on rugged fitness landscapes with many local fitness peaks.
Unlike the previous studies in which the initial fitness of the population is
assumed to be high, here we start the adaptation process with a low fitness
corresponding to a population in a stressful novel environment. For generic
fitness distributions, using an analytic argument we find that the average
number of steps to a local optimum varies logarithmically with the genotype
sequence length and increases as the correlations amongst genotypic fitnesses
increase. When the fitnesses are exponentially or uniformly distributed, using
an evolution equation for the distribution of population fitness, we
analytically calculate the fitness distribution of fixed beneficial mutations
and the walk length distribution.


  The surface accessibility of {\alpha}-bungarotoxin has been investigated by
using Gd2L7, a newly designed paramagnetic NMR probe. Signal attenuations
induced by Gd2L7 on {\alpha}-bungarotoxin C{\alpha}H peaks of 1H-13C HSQC
spectra have been analyzed and compared with the ones previously obtained in
the presence of GdDTPA-BMA. In spite of the different molecular size and shape,
for the two probes a common pathway of approach to the {\alpha}-bungarotoxin
surface can be observed with an equally enhanced access of both GdDTPA-BMA and
Gd2L7 towards the protein surface side where the binding site is located.
Molecular dynamics simulations suggest that protein backbone flexibility and
surface hydration contribute to the observed preferential approach of both
gadolinium complexes specifically to the part of the {\alpha}-bungarotoxin
surface which is involved in the interaction with its physiological target, the
nicotinic acetylcholine receptor.


  Life occurs in ionic solutions, not pure water. The ionic mixtures of these
solutions are very different from water and have dramatic effects on the cells
and molecules of biological systems, yet theories and simulations cannot
calculate their properties. I suggest the reason is that existing theories stem
from the classical theory of ideal or simple gases in which (to a first
approximation) atoms do not interact. Even the law of mass action describes
reactants as if they were ideal. I propose that theories of ionic solutions
should start with the theory of complex fluids because that theory is designed
to deal with interactions from the beginning. The variational theory of complex
fluids is particularly well suited to describe mixtures like the solutions in
and outside biological cells. When a component or force is added to a solution,
the theory derives - by mathematics alone - a set of partial differential
equations that captures the resulting interactions self-consistently. Such a
theory has been implemented and shown to be computable in biologically relevant
systems but it has not yet been thoroughly tested in equilibrium or flow.


  Approximate Bayesian computation (ABC) is a class of algorithmic methods in
Bayesian inference using statistical summaries and computer simulations. ABC
has become popular in evolutionary genetics and in other branches of biology.
However model selection under ABC algorithms has been a subject of intense
debate during the recent years. Here we propose novel approaches to model
selection based on posterior predictive distributions and approximations of the
deviance. We argue that this framework can settle some contradictions between
the computation of model probabilities and posterior predictive checks using
ABC posterior distributions. A simulation study and an analysis of a
resequencing data set of human DNA show that the deviance criteria lead to
sensible results in a number of model choice problems of interest to population
geneticists.


  The nonequilibrium stationary state of an exclusive genetic switch is
considered. The model comprises two competing species and a single binding site
which, when bound to by a protein of one species, causes the other species to
be repressed. The model may be thought of as a minimal model of the power
struggle between two competing parties. Exact solutions are given for the
limits of vanishing binding/unbinding rates and infinite binding/unbinding
rates. A mean field theory is introduced which is exact in the limit of
vanishing binding/unbinding rates. The mean field theory and numerical
simulations reveal that generically bistability occurs and the system is in a
symmetry broken state. An exact perturbative solution which in principle allows
the nonequilibrium stationary state to be computed is also developed and
computed to first and second order.


  Fitness effects of mutations fall on a continuum ranging from lethal to
deleterious to beneficial. The distribution of fitness effects (DFE) among
random mutations is an essential component of every evolutionary model and a
mathematical portrait of robustness. Recent experiments on five viral species
all revealed a characteristic bimodal shaped DFE, featuring peaks at neutrality
and lethality. However, the phenotypic causes underlying observed fitness
effects are still unknown, and presumably thought to vary unpredictably from
one mutation to another. By combining population genetics simulations with a
simple biophysical protein folding model, we show that protein thermodynamic
stability accounts for a large fraction of observed mutational effects. We
assume that moderately destabilizing mutations inflict a fitness penalty
proportional to the reduction in folded protein, which depends continuously on
folding free energy (\Delta G). Most mutations in our model affect fitness by
altering \Delta G, while, based on simple estimates, \approx10% abolish
activity and are unconditionally lethal. Mutations pushing \Delta G>0 are also
considered lethal. Contrary to neutral network theory, we find that, in
mutation/selection/drift steady-state, high mutation rates (m) lead to less
stable proteins and a more dispersed DFE, i.e. less mutational robustness.
Small population size (N) also decreases stability and robustness. In our
model, a continuum of non-lethal mutations reduces fitness by \approx2% on
average, while \approx10-35% of mutations are lethal, depending on N and m.
Compensatory mutations are common in small populations with high mutation
rates. More broadly, we conclude that interplay between biophysical and
population genetic forces shapes the DFE.


  We investigate a model where strong noise in a sub-population creates a
metastable state in an otherwise unstable two-population system. The induced
metastable state is vortex-like, and its persistence time grows exponentially
with the noise strength. A variety of distinct scaling relations are observed
depending on the relative strength of the sub-population noises.


  Dense cooperative networks are an essential element of social capital for a
prosperous society. These networks enable individuals to overcome collective
action dilemmas by enhancing trust. In many biological and social settings,
network structures evolve endogenously as agents exit relationships and build
new ones. However, the process by which evolutionary dynamics lead to
self-organization of dense cooperative networks has not been explored. Our
large group prisoner's dilemma experiments with exit and partner choice options
show that core-periphery segregation of cooperators and defectors drives the
emergence of cooperation. Cooperators' Quit-for-Tat and defectors' Roving
strategy lead to a highly asymmetric core and periphery structure. Densely
connected to each other, cooperators successfully isolate defectors and earn
larger payoffs than defectors. Our analysis of the topological characteristics
of evolving networks illuminates how social capital is generated.


  We study biological evolution in a high-dimensional genotype space in the
regime of rare mutations and strong selection. The population performs an
uphill walk which terminates at local fitness maxima. Assigning fitness
randomly to genotypes, we show that the mean walk length is logarithmic in the
number of initially available beneficial mutations, with a prefactor determined
by the tail of the fitness distribution. This result is derived analytically in
a simplified setting where the mutational neighborhood is fixed during the
adaptive process, and confirmed by numerical simulations.


  Robustness of biochemical systems has become one of the central questions in
Systems Biology, although it is notoriously difficult to formally capture its
multifaceted nature. Maintenance of normal system function depends not only on
the stoichiometry of the underlying interrelated components, but also on a
multitude of kinetic parameters. For given parameter values, recent findings
have aimed at characterizing the property of the system components to exhibit
same concentrations in the resulting steady states, termed absolute
concentration robustness (ACR). However, the existing method for determining
system components exhibiting ACR is applicable only to one class of mass-action
networks for which this property can be confirmed, but not discarded. Here we
design a new method which relies on biochemical network decompositions into
subnetworks, called elementary flux modes, to identify ACR in a broader class
of mass-action networks by using only the given stoichiometry. This approach
reduces the problem of determining ACR to that of solving parameterized systems
of linear equations, rendering it amenable to networks of larger sizes. Our
unified framework will be helpful in analyzing this biologically important type
of robustness as well as detection of novel systemic properties independent of
the kinetic parameters for more complex biochemical networks.


  Understanding patterns of selectively neutral genetic variation is essential
in order to model deviations from neutrality, caused for example by different
forms of selection. Best understood is neutral genetic variation at a single
locus, but additional insights can be gained by investigating genetic variation
at multiple loci. The corresponding patterns of variation reflect linkage
disequilibrium and provide information about the underlying multi-locus gene
genealogies. The statistical properties of two-locus genealogies have been
intensively studied for populations of constant census size, as well as for
simple demographic histories such as exponential population growth, and single
bottlenecks. By contrast, the combined effect of recombination and sustained
demographic fluctuations is poorly understood. Addressing this issue, we study
a two-locus Wright-Fisher model of a population subject to recurrent
bottlenecks. We derive coalescent approximations for the covariance of the
times to the most recent common ancestor at two loci. We find, first, that an
effective population-size approximation describes the numerically observed
linkage disequilibrium provided that recombination occurs either much faster or
much more slowly than the population size changes. Second, when recombination
occurs frequently between bottlenecks but rarely within bottlenecks, we observe
long-range linkage disequilibrium. Third, we show that in the latter case, a
commonly used measure of linkage disequilibrium, sigma_d^2 (closely related to
r^2), fails to capture long-range linkage disequilibrium because constituent
terms, each reflecting long-range linkage disequilibrium, cancel. Fourth, we
analyse a limiting case in which long-range linkage disequilibrium can be
described in terms of a Xi-coalescent process allowing for simultaneous
multiple mergers of ancestral lines.


  To explain the observed dynamics in equilibrium single-molecule measurements
of biomolecules, the experimental observable is often chosen as a putative
reaction coordinate along which kinetic behavior is presumed to be governed by
diffusive dynamics. Here, we invoke the splitting probability as a test of the
suitability of such a proposed reaction coordinate. Comparison of the observed
splitting probability with that computed from the kinetic model provides a
simple test to reject poor reaction coordinates. We demonstrate this test for a
force spectroscopy measurement of a DNA hairpin.


  Resonance energy transfer methods are in wide use for evaluating
protein-protein interactions and protein conformational changes in living
cells. Fluorescence resonance energy transfer (FRET) measures energy transfer
as a function of the acceptor:donor ratio, generating FRET saturation curves.
Modeling these curves by Michaelis-Menten kinetics allows characterization by
two parameters, which serve to evaluate apparent affinity between two proteins
and to compare this affinity in different experimental conditions. To reduce
the effect of sampling variability, several statistical samples of the
saturation curve are generated in the same biological conditions. Here we study
three procedures to determine whether statistical samples in a collection are
homogeneous, in the sense that they are extracted from the same regression
model. From the hypothesis testing viewpoint, we considered an F test and a
procedure based on bootstrap resampling. The third method analyzed the problem
from the model selection viewpoint, and used the Akaike information criterion
(AIC). Although we only considered the Michaelis-Menten model, all statistical
procedures would be applicable to any other nonlinear regression model. We
compared the performance of the homogeneity testing methods in a Monte Carlo
study and through analysis in living cells of FRET saturation curves for
dimeric complexes of CXCR4, a seven-transmembrane receptor of the G
protein-coupled receptor family. We show that the F test, the bootstrap
procedure and the model selection method lead in general to similar
conclusions, although AIC gave the best results when sample sizes were small,
whereas the F test and the bootstrap method were more appropriate for large
samples. In practice, all three methods are easy to use simultaneously and show
consistency, facilitating conclusions on sample homogeneity.


  We extend the Moran model with single-crossover recombination to include
general recombination and mutation. We show that, in the case without
resampling, the expectations of products of marginal processes defined via
partitions of sites form a closed hierarchy, which is exhaustively described by
a finite system of differential equations. One thus has the exceptional
situation of moment closure in a nonlinear system. Surprisingly, this property
is lost when resampling (i.e., genetic drift) is included.


  When a population inhabits an inhomogeneous environment, the fitness value of
traits can vary with the position in the environment. Gene flow caused by
random mating can nevertheless prevent that a sexually reproducing population
splits into different species under such circumstances. This is the problem of
sympatric speciation. However, mating need not be entirely random. Here, we
present a model where the individually advantageous preference for partners of
high fitness can lead to genetic clustering as a precondition for speciation.
In simulations, in appropriate parameter regimes, our model leads to the rapid
fixation of the corresponding alleles.


  In this paper we present a biologically detailed mathematical model of
tripartite synapses, where astrocytes modulate short-term synaptic plasticity.
The model consists of a pre-synaptic bouton, a post-synaptic dendritic
spine-head, a synaptic cleft and a peri-synaptic astrocyte controlling Ca2+
dynamics inside the synaptic bouton. This in turn controls glutamate release
dynamics in the cleft. As a consequence of this, glutamate concentration in the
cleft has been modeled, in which glutamate reuptake by astrocytes has also been
incorporated. Finally, dendritic spine-head dynamics has been modeled. As an
application, this model clearly shows synaptic potentiation in the hippocampal
region, i.e., astrocyte Ca2+ mediates synaptic plasticity, which is in
conformity with the majority of the recent findings (Perea & Araque, 2007;
Henneberger et al., 2010; Navarrete et al., 2012).


  In eukaryotic cell nuclei, a variety of DNA interactions with nuclear
elements occur, which, in combination with intra- and inter- chromosomal
cross-talks, shape a functional 3D architecture. In some cases they are
organized by active, i.e. actin/myosin, motors. More often, however, they have
been related to passive diffusion mechanisms. Yet, the crucial questions on how
DNA loci recognize their target and are reliably shuttled to their destination
by Brownian diffusion are still open. Here, we complement the current
experimental scenario by considering a physics model, in which the interaction
between distant loci is mediated by diffusing bridging molecules. We show that,
in such a system, the mechanism underlying target recognition and
colocalization is a thermodynamic switch-like process (a phase transition) that
only occurs if the concentration and affinity of binding molecules is above a
threshold, or else stable contacts are not possible. We also briefly discuss
the kinetics of this "passive-shuttling" process, as produced by random
diffusion of DNA loci and their binders, and derive predictions based on the
effects of genomic modifications and deletions.


  Genomic copy number variation (CNV) is a large source of variation between
organisms, and its consequences include phenotypic differences and genetic
disorders. CNVs are commonly detected by hybridizing genomic DNA to microarrays
of nucleic acid probes. System noise caused by operational and probe
performance variability complicates the interpretation of these data. To
minimize the distortion of genetic signal by system noise, we have explored the
latter in an archive of hybridizations in which no genetic signal is expected.
This archive is obtained by comparative genomic hybridization (CGH) of a sample
in one channel to the same sample in the other channel, or 'self-self' data.
These self-self hybridizations trap a variety of system noise inherent in
sample-reference (test) data. Through singular value decomposition (SVD) of
self-self data, we have determined the principal components of system noise.
Assuming simple linear models of noise generation, the linear correction of
test data with self-self data -or 'system normalization'- reduces local and
long-range correlations and improves signal-to-noise metrics, yet does not
introduce detectable spurious signal. Using this method, 90% of hybridizations
displayed improved signal-to-noise ratios with an average increase of 7.0%, due
mainly to a reduced median average deviation (MAD). In addition, we have found
that principal component loadings correlate with specific probe variables
including array coordinates, base composition, and proximity to the 5' ends of
genes. The correlation of the principal component loadings with the test data
depends on operational variables, such as the temporal order of processing and
the localization of individual samples within 96-well plates.


  Statistical thermodynamics basis of energy and residue position fluctuations
is explained for native proteins. The protein and its surroundings are treated
as a canonical system with emphasis on the effects of energy exchange between
the two. Fluctuations of the energy are related to fluctuations of residue
positions, which in turn are related to the connectivity matrix of the protein,
thus establishing a connection between energy fluctuation pathways and protein
architecture. The model gives the locations of hotspots for ligand binding, and
identifies the pathways of energy conduction within the protein. Results are
discussed in terms of two sets of models, the BPTI and twelve proteins that
contain the PDZ domain. Possible use of the model for determining functionally
similar domains in a diverse set of proteins is pointed out.


  The replicator equation is ubiquitous for many areas of mathematical biology.
One of major shortcomings of this equation is that it does not allow for an
explicit spatial structure. Here we review analytical approaches to include
spatial variables to the system. We also provide a concise exposition of the
results concerning the appearance of spatial patterns in replicator
reaction-diffusion systems.


  In living cells, molecular motors convert chemical energy into mechanical
work. Its thermodynamic energy efficiency, i.e. the ratio of output mechanical
work to input chemical energy, is usually high. However, using two-state
models, we found the motion of molecular motors is loosely coupled to the
chemical cycle. Only part of the input energy can be converted into mechanical
work. Others is dissipated into environment during substeps without
contributions to the macro scale unidirectional movement.


  In this study, a two-state mechanochemical model is presented to describe the
dynamic instability of microtubules (MTs) in cells. The MTs switches between
two states, assembly state and disassembly state. In assembly state, the growth
of MTs includes two processes: free GTP-tubulin binding to the tip of
protofilament (PF) and conformation change of PF, during which the first
tubulin unit which curls outwards is rearranged into MT surface using the
energy released from the hydrolysis of GTP in the penultimate tubulin unit. In
disassembly state, the shortening of MTs includes also two processes, the
release of GDP-tibulin from the tip of PF and one new tubulin unit curls out of
the MT surface. Switches between these two states, which are usually called
rescue and catastrophe, happen stochastically with external force dependent
rates. Using this two-state model with parameters obtained by fitting the
recent experimental data, detailed properties of MT growth are obtained, we
find that MT is mainly in assembly state, its mean growth velocity increases
with external force and GTP-tubulin concentration, MT will shorten in average
without external force. To know more about the external force and GTP-tubulin
concentration dependent properties of MT growth, and for the sake of the future
experimental verification of this two-state model, eleven {\it critical forces}
are defined and numerically discussed.


  The position of propagating population fronts fluctuates because of the
discreteness of the individuals and stochastic character of processes of birth,
death and migration. Here we consider a Markov model of a population front
propagating into a metastable state, and focus on the weak noise limit. For
typical, small fluctuations the front motion is diffusive, and we calculate the
front diffusion coefficient. We also determine the probability distribution of
rare, large fluctuations of the front position and, for a given average front
velocity, find the most likely population density profile of the front.
Implications of the theory for population extinction risk are briefly
considered.


  The quasispecies model of biological evolution for asexual organisms such as
bacteria and viruses has attracted considerable attention of biological
physicists. Many variants of the model have been proposed and subsequently
solved using the methods of statistical physics. In this paper I will put
forward important but largely overlooked relations between localization theory,
random matrices, and the quasispecies model. These relations will help me to
study the dynamics of this model. In particular, I will show that the
distribution of times between evolutionary jumps in the genotype space follows
a power law, in agreement with recent findings in the shell model - a
simplified version of the quasispecies model.


  Understanding the influence of macromolecular crowding and nanoparticles on
the formation of in-register $\beta$-sheets, the primary structural component
of amyloid fibrils, is a first step towards describing \emph{in vivo} protein
aggregation and interactions between synthetic materials and proteins. Using
all atom molecular simulations in implicit solvent we illustrate the effects of
nanoparticle size, shape, and volume fraction on oligomer formation of an
amyloidogenic peptide from the transthyretin protein. Surprisingly, we find
that inert spherical crowding particles destabilize in-register $\beta$-sheets
formed by dimers while stabilizing $\beta$-sheets comprised of trimers and
tetramers. As the radius of the nanoparticle increases crowding effects
decrease, implying smaller crowding particles have the largest influence on the
earliest amyloid species. We explain these results using a theory based on the
depletion effect. Finally, we show that spherocylindrical crowders destabilize
the ordered $\beta$-sheet dimer to a greater extent than spherical crowders,
which underscores the influence of nanoparticle shape on protein aggregation.


  Perfect phylogeny consisting of determining the compatibility of a set of
characters is known to be NP-complete. We propose in this article a conjecture
on the necessary and sufficient conditions of compatibility: Given a set
$\mathcal{C}$ of $r$-states full characters, there exists a function $f(r)$
such that $\mathcal{C}$ is compatible iff every set of $f(r)$ characters of
$\mathcal{C}$ is compatible. Some previous work showed that $f(2)=2$, $f(3)=3$
and $f(r) \ge r-1$. Gusfield et al. 09 conjectured that $f(r) = r$ for any $r
\ge 2$. In this paper, we present an example showing that $f(4) \ge 5$ and then
a closure operation for chordal sandwich graphs. The later problem is a common
approach of perfect phylogeny. This operation can be the first step to simplify
the problem before solving some particular cases $f(4), f(5), ... $, and
determining the function $f(r)$.


  An explicit solution for a general two-type birth-death branching process
with one way mutation is presented. This continuous time process mimics the
evolution of resistance to treatment, or the onset of an extra driver mutation
during tumor progression. We obtain the exact generating function of the
process at arbitrary times, and derive various large time scaling limits. In
the simultaneous small mutation rate and large time scaling limit, the
distribution of the mutant cells develops some atypical properties, including a
power law tail and diverging average.


  Water and water-mediated interactions determine thermodynamic and kinetics of
protein folding, protein aggregation and self-assembly in confined spaces. To
obtain insights into the role of water in the context of folding problems, we
describe computer simulations of a few related model systems. The dynamics of
collapse of eicosane shows that upon expulsion of water the linear hydrocarbon
chain adopts an ordered helical hairpin structure with 1.5 turns. The structure
of dimer of eicosane molecules has two well ordered helical hairpins that are
stacked perpendicular to each other. As a prelude to studying folding in
confined spaces we used simulations to understand changes in hydrophobic and
ionic interactions in nano droplets. Solvation of hydrophobic and charged
species change drastically in nano water droplets. Hydrophobic species are
localized at the boundary. The tendency of ions to be at the boundary where
water density is low increases as the charge density decreases. Interaction
between hydrophobic, polar, and charged residue are also profoundly altered in
confined spaces. Using the results of computer simulations and accounting for
loss of chain entropy upon confinement we argue and then demonstrate, using
simulations in explicit water, that ordered states of generic amphiphilic
peptide sequences should be stabilized in cylindrical nanopores.


  Given an endogenous timescale set by invasion in a constant environment, we
introduced periodic temporal variation in competitive superiority by
alternating the species' propagation rates. By manipulating habitat size and
introduction rate, we simulated environments where successful invasion proceeds
through growth of many spatial clusters, and where invasion can occur only as a
single-cluster process. In the multi-cluster invasion regime, rapid
environmental variation produced spatial mixing of the species and
non-equilibrium coexistence. The dynamics' dominant response effectively
averaged environmental fluctuation, so that each species could avoid
competitive exclusion. Increasing the environment's half-period to match the
population-dynamic timescale let the (initially) more abundant resident
repeatedly repel the invader. Periodic transition in propagation-rate advantage
rarely interrupted the exclusion process when the more abundant species had
competitive advantage. However, at infrequent and randomly occurring times, the
rare species could invade and reverse the density pattern by rapidly eroding
the resident's preemption of space. In the single-cluster invasion regime,
environmental variation occurring faster than the population-dynamic timescale
prohibited successful invasion; the first species to reach its stationary
density (calculated for a constant environment) continued to repel the other
during long simulations. When the endogenous and exogenous timescales matched,
the species randomly reversed roles of resident and invader; the waiting times
for reversal of abundances indicate stochastic resonance. For both invasion
regimes, environmental fluctuation occurring much slower than the endogenous
dynamics produced symmetric limit cycles, alternations of the
constant-environment pattern.


  We derive a convex optimization problem on a steady-state nonequilibrium
network of biochemical reactions, with the property that energy conservation
and the second law of thermodynamics both hold at the problem solution. This
suggests a new variational principle for biochemical networks that can be
implemented in a computationally tractable manner. We derive the Lagrange dual
of the optimization problem and use strong duality to demonstrate that a
biochemical analogue of Tellegen's theorem holds at optimality. Each optimal
flux is dependent on a free parameter that we relate to an elementary kinetic
parameter when mass action kinetics is assumed.


  In analogon to successful artificial evolution simulations as Tierra or
avida, this text presents a way to perform artificial evolution in a native x86
system. The implementation of the artificial chemistry and first results of
statistical experiments are presented.


  Quantum simulations constructing probability tensors of biological multi-taxa
in phylogenetic trees are proposed, in terms of positive trace preserving maps,
describing evolving systems of quantum walks with multiple walkers. Basic
phylogenetic models applying on trees of various topologies are simulated
following appropriate decoherent quantum circuits. Quantum simulations of
statistical inference for aligned sequences of biological characters are
provided in terms of a quantum pruning map operating on likelihood operator
observables, utilizing state-observable duality and measurement theory.


  Recent experimental studies have shown that amyloid fibril formed by
aggregation of {\beta} peptide exhibits excellent mechanical properties
comparable to other protein materials such as actin filaments and microtubules.
These excellent mechanical properties of amyloid fibrils are related to their
functional role in disease expression. This indicates the necessity to
understand how an amyloid fibril achieves the remarkable mechanical properties
through self-aggregation with structural hierarchy, highlighting the
structure-property-function relationship for amyloids, whereas such
relationship still remains elusive. In this work, we have studied the
mechanical properties of human islet amyloid polypeptide (hIAPP) with respect
to its structural hierarchies and structural shapes by coarse-grained normal
mode analysis. Our simulation shows that hIAPP fibril can achieve the excellent
bending rigidity via specific aggregation pattern such as antiparallel stacking
of {\beta} peptides. Moreover, we have found the length-dependent mechanical
properties of amyloids. This length-dependent property has been elucidated from
Timoshenko beam model that takes into account the shear effect on the bending
of amyloids. In summary, our study sheds light on the importance of not only
the molecular architecture, which encodes the mechanical properties of the
fibril, but also the shear effect on the mechanical (bending) behavior of the
fibril.


  We present and discuss new importance sampling schemes for the approximate
computation of the sample probability of observed genetic types in the
infinitely many sites model from population genetics. More specifically, we
extend the 'classical framework', where genealogies are assumed to be governed
by Kingman's coalescent, to the more general class of Lambda-coalescents and
develop further Hobolth et. al.'s (2008) idea of deriving importance sampling
schemes based on 'compressed genetrees'. The resulting schemes extend earlier
work by Griffiths and Tavar\'e (1994), Stephens and Donnelly (2000), Birkner
and Blath (2008) and Hobolth et. al. (2008). We conclude with a performance
comparison of classical and new schemes for Beta- and Kingman coalescents.


  Forty years ago, Robert May questioned a central belief in ecology by proving
that sufficiently large or complex ecological networks have probability of
persisting close to zero. To prove this point, he analyzed large networks in
which species interact at random. However, in natural systems pairs of species
have well-defined interactions (e.g., predator-prey, mutualistic or
competitive). Here we extend May's results to these relationships and find
remarkable differences between predator-prey interactions, which increase
stability, and mutualistic and competitive, which are destabilizing. We provide
analytic stability criteria for all cases. These results have broad
applicability in ecology. For example, we show that, surprisingly, the
probability of stability for predator-prey networks is decreased when we impose
realistic food web structure or we introduce a large preponderance of weak
interactions. Similarly, stability is negatively impacted by nestedness in
bipartite mutualistic networks.


  Classical ecological theory predicts that environmental stochasticity
increases extinction risk by reducing the average per-capita growth rate of
populations. To understand the interactive effects of environmental
stochasticity, spatial heterogeneity, and dispersal on population growth, we
study the following model for population abundances in $n$ patches: the
conditional law of $X_{t+dt}$ given $X_t=x$ is such that when $dt$ is small the
conditional mean of $X_{t+dt}^i-X_t^i$ is approximately $[x^i\mu_i+\sum_j(x^j
D_{ji}-x^i D_{ij})]dt$, where $X_t^i$ and $\mu_i$ are the abundance and per
capita growth rate in the $i$-th patch respectivly, and $D_{ij}$ is the
dispersal rate from the $i$-th to the $j$-th patch, and the conditional
covariance of $X_{t+dt}^i-X_t^i$ and $X_{t+dt}^j-X_t^j$ is approximately $x^i
x^j \sigma_{ij}dt$. We show for such a spatially extended population that if
$S_t=(X_t^1+...+X_t^n)$ is the total population abundance, then $Y_t=X_t/S_t$,
the vector of patch proportions, converges in law to a random vector $Y_\infty$
as $t\to\infty$, and the stochastic growth rate $\lim_{t\to\infty}t^{-1}\log
S_t$ equals the space-time average per-capita growth rate
$\sum_i\mu_i\E[Y_\infty^i]$ experienced by the population minus half of the
space-time average temporal variation $\E[\sum_{i,j}\sigma_{ij}Y_\infty^i
Y_\infty^j]$ experienced by the population. We derive analytic results for the
law of $Y_\infty$, find which choice of the dispersal mechanism $D$ produces an
optimal stochastic growth rate for a freely dispersing population, and
investigate the effect on the stochastic growth rate of constraints on
dispersal rates. Our results provide fundamental insights into "ideal free"
movement in the face of uncertainty, the persistence of coupled sink
populations, the evolution of dispersal rates, and the single large or several
small (SLOSS) debate in conservation biology.


  We establish that mass conserving single terminal-linkage networks of
chemical reactions admit positive steady states regardless of network
deficiency and the choice of reaction rate constants. This result holds for
closed systems without material exchange across the boundary, as well as for
open systems with material exchange at rates that satisfy a simple sufficient
and necessary condition. Our proof uses a fixed point of a novel convex
optimization formulation to find the steady state behavior of chemical reaction
networks that satisfy the law of mass-action kinetics. A fixed point iteration
can be used to compute these steady states, and we show that it converges for
weakly reversible homogeneous systems. We report the results of our algorithm
on numerical experiments.


  The rates of protein folding with photon absorption or emission and the cross
section of photon -protein inelastic scattering are calculated from the quantum
folding theory by use of standard field-theoretical method. All these protein
photo-folding processes are compared with common protein folding without
interaction of photons (nonradiative folding). It is demonstrated that there
exists a common factor (thermo-averaged overlap integral of vibration wave
function, TAOI) for protein folding and protein photo-folding. Based on this
finding it is predicted that: 1) the stimulated photo-folding rates show the
same temperature dependence as protein folding; 2) the spectral line of
electronic transition is broadened to a band which includes abundant vibration
spectrum without and with conformational transition and the width of the
vibration spectral line is largely reduced; 3) the resonance fluorescence cross
section changes with temperature obeying the same law (Luo-Lu's law). The
particular form of the folding rate - temperature relation and the abundant
spectral structure imply the existence of a set of quantum oscillators in the
transition process and these oscillators are mainly of torsion type of low
frequency, imply the quantum tunneling between protein conformations does exist
in folding and photo-folding processes and the tunneling is rooted deeply in
the coherent motion of the conformational-electronic system.


  The present article is an invited contribution to the Encyclopedia of
Complexity and System Science, Robert A. Meyers Ed., Springer New York (2009).
It is a review of the biophysical mechanisms that underly cell motility. It
mainly focuses on the eukaryotic cytoskeleton and cell-motility mechanisms.
Bacterial motility as well as the composition of the prokaryotic cytoskeleton
is only briefly mentioned. The article is organized as follows. In Section III,
I first present an overview of the diversity of cellular motility mechanisms,
which might at first glance be categorized into two different types of
behaviors, namely "swimming" and "crawling". Intracellular transport, mitosis -
or cell division - as well as other extensions of cell motility that rely on
the same essential machinery are briefly sketched. In Section IV, I introduce
the molecular machinery that underlies cell motility - the cytoskeleton - as
well as its interactions with the external environment of the cell and its main
regulatory pathways. Sections IV D to IV F are more detailed in their
biochemical presentations; readers primarily interested in the theoretical
modeling of cell motility might want to skip these sections in a first reading.
I then describe the motility mechanisms that rely essentially on
polymerization-depolymerization dynamics of cytoskeleton filaments in Section
V, and the ones that rely essentially on the activity of motor proteins in
Section VI. Finally, Section VII is devoted to the description of the
integrated approaches that have been developed recently to try to understand
the cooperative phenomena that underly self-organization of the cell
cytoskeleton as a whole.


  The analysis of complex networks permeates all sciences, from biology to
sociology. A fundamental, unsolved problem is how to characterize the community
structure of a network. Here, using both standard and novel benchmarks, we show
that maximization of a simple global parameter, which we call Surprise (S),
leads to a very efficient characterization of the community structure of
complex synthetic networks. Particularly, S qualitatively outperforms the most
commonly used criterion to define communities, Newman and Girvan's modularity
(Q). Applying S maximization to real networks often provides natural,
well-supported partitions, but also sometimes counterintuitive solutions that
expose the limitations of our previous knowledge. These results indicate that
it is possible to define an effective global criterion for community structure
and open new routes for the understanding of complex networks.


  Apoptosis is a complex pathway regulated by the concerted action of multiple
pro- and anti-apoptotic molecules. The intrinsic (mitochondrial) pathway of
apoptosis is governed up-stream of mitochondria, by the family of Bcl-2
proteins, and down-stream of mitochondria, by low-probability events, such as
apoptosome formation, and by feedback circuits involving caspases and inhibitor
of apoptosis proteins (IAPs), such as XIAP. All these regulatory mechanisms
ensure that cells only commit to death once a threshold of damage has been
reached and the anti-apoptotic reserve of the cell is overcome. As cancer cells
are invariably exposed to strong intracellular and extracellular stress
stimuli, they are particularly reliant on the expression of anti-apoptotic
proteins. Hence, many cancer cells undergo apoptosis when exposed to agents
that inhibit anti-apoptotic Bcl-2 molecules, such as BH3 mimetics, while normal
cells remain relatively insensitive to single agent treatments with the same
class of molecules. Targeting different proteins within the apoptotic network
with combinatorial treatment approaches often achieves even greater
specificity. This led us to investigate the sensitivity of leukemia and
lymphoma cells to a pro-apoptotic action of a BH3 mimetic combined with a small
molecule inhibitor of XIAP. Using computational probabilistic model of
apoptotic pathway, verified by experimental results from human leukemia and
lymphoma cell lines, we show that inhibition of XIAP has a non-linear effect on
sensitization towards apoptosis induced by the BH3 mimetic HA14-1. This study
justifies further ex vivo and animal studies on the potential of the treatment
of leukemia and lymphoma with a combination of BH3 mimetics and XIAP
inhibitors.


  In this technical note we calculate the dynamics of a linear feedback model
of progression in the cell cycle in the case that the cells are organized into
k=3 clusters. We examine the dynamics in detail for a specific subset of
parameters with non-empty interior.
  There is an interior fixed point of the Poincare' map defined by the system.
This fixed point corresponds to a periodic solution with period $T$ in which
the three cluster exchange positions after time $T/3$. We call this solution
3=cyclic. In all the parameters studied, the fixed point is either:
  * isolated and locally unstable, or,
  * contained in a neutrally stable set of period 3 points.
  In the later case the edges of the neutrally stable set are unstable. This
case exists if either the three clusters are isolated from each other, or, if
they interact in a non-essential way. In both cases the orbits of all other
interior points are asymptotic to the boundary. Thus 3-cyclic solutions are
practically unstable in the sense the arbitrarily small perturbations may lead
to loss of stability and eventual merger of clusters. Since the single cluster
solution (synchronization) is the only solution that is asymptotically stable,
it would seem to be the most likely to be observed in application if the
feedback is similar to the form we propose and is positive.


  In self-assembly processes, kinetic trapping effects often hinder the
formation of thermodynamically stable ordered states. In a model of viral
capsid assembly and in the phase transformation of a lattice gas, we show how
simulations in a self-assembling steady state can be used to identify two
distinct mechanisms of kinetic trapping. We argue that one of these mechanisms
can be adequately captured by kinetic rate equations, while the other involves
a breakdown of theories that rely on cluster size as a reaction coordinate. We
discuss how these observations might be useful in designing and optimising
self-assembly reactions.


  In the present paper we propose to describe gene networks in biological
systems using probabilistic algorithms. We describe gene duplication in the
process of biological evolution using introduction of the replica procedure for
probabilistic algorithms. We construct the examples of such a replica procedure
for hidden Markov models. We introduce the family of hidden Markov models where
the set of hidden states is a finite additive group with a p-adic metric and
build the replica procedure for this family of markovian models.


  We study the totally asymmetric simple exclusion process (TASEP) on complex
networks, as a paradigmatic model for transport subject to excluded volume
interactions. Building on TASEP phenomenology on a single segment and borrowing
ideas from random networks we investigate the effect of connectivity on
transport. In particular, we argue that the presence of disorder in the
topology of vertices crucially modifies the transport features of a network:
irregular networks involve homogeneous segments and have a bimodal distribution
of edge densities, whereas regular networks are dominated by shocks leading to
a unimodal density distribution. The proposed numerical approach of solving for
mean-field transport on networks provides a general framework for studying
TASEP on large networks, and is expected to generalize to other transport
processes.


  Network analysis became a powerful tool in recent years. Heat shock is a
well-characterized model of cellular dynamics. S. cerevisiae is an appropriate
model organism, since both its protein-protein interaction network
(interactome) and stress response at the gene expression level have been well
characterized. However, the analysis of the reorganization of the yeast
interactome during stress has not been investigated yet. We calculated the
changes of the interaction-weights of the yeast interactome from the changes of
mRNA expression levels upon heat shock. The major finding of our study is that
heat shock induced a significant decrease in both the overlaps and connections
of yeast interactome modules. In agreement with this the weighted diameter of
the yeast interactome had a 4.9-fold increase in heat shock. Several key
proteins of the heat shock response became centers of heat shock-induced local
communities, as well as bridges providing a residual connection of modules
after heat shock. The observed changes resemble to a "stratus-cumulus" type
transition of the interactome structure, since the unstressed yeast interactome
had a globally connected organization, similar to that of stratus clouds,
whereas the heat shocked interactome had a multifocal organization, similar to
that of cumulus clouds. Our results showed that heat shock induces a partial
disintegration of the global organization of the yeast interactome. This change
may be rather general occurring in many types of stresses. Moreover, other
complex systems, such as single proteins, social networks and ecosystems may
also decrease their inter-modular links, thus develop more compact modules, and
display a partial disintegration of their global structure in the initial phase
of crisis. Thus, our work may provide a model of a general, system-level
adaptation mechanism to environmental changes.


  Recently, a three-state model is presented to describe the intracellular
traffic of unconventional (single-headed) kinesin KIF1A [Phys. Rev. Lett. {\bf
95}, 118101 (2005)], in which each motor can bind strongly or weakly to its
microtubule track, and each binding site of the track might be empty or
occupied by one motor. As the usual two-state model, i.e. the totally
asymmetric simple exclusion process (TASEP) with motor detachment and
attachment, in steady state of the system, this three-state model also exhibits
shock (or domain wall separating the high-density and low density phases) and
boundary layers. In this study, using mean-field analysis, the conditions of
existence of shock and boundary layers are obtained theoretically. Combined
with numerical calculations, the properties of shock are also studied. This
study will be helpful to understand the biophysical properties of the
collective transport of kinesin KIF1A.


  We discuss the distribution of various estimators for extracting the
diffusion constant of single Brownian trajectories obtained by fitting the
squared displacement of the trajectory. The analysis of the problem can be
framed in terms of quadratic functionals of Brownian motion that correspond to
the Euclidean path integral for simple Harmonic oscillators with time dependent
frequencies. Explicit analytical results are given for the distribution of the
diffusion constant estimator in a number of cases and our results are confirmed
by numerical simulations.


  Predictive understanding of the myriads of signal transduction pathways in a
cell is an outstanding challenge of systems biology. Such pathways are
primarily mediated by specific but transient protein-protein interactions,
which are difficult to study experimentally. In this study, we dissect the
specificity of protein-protein interactions governing two-component signaling
(TCS) systems ubiquitously used in bacteria. Exploiting the large number of
sequenced bacterial genomes and an operon structure which packages many pairs
of interacting TCS proteins together, we developed a computational approach to
extract a molecular interaction code capturing the preferences of a small but
critical number of directly interacting residue pairs. This code is found to
reflect physical interaction mechanisms, with the strongest signal coming from
charged amino acids. It is used to predict the specificity of TCS interaction:
Our results compare favorably to most available experimental results, including
the prediction of 7 (out of 8 known) interaction partners of orphan signaling
proteins in Caulobacter crescentus. Surveying among the available bacterial
genomes, our results suggest 15~25% of the TCS proteins could participate in
out-of-operon "crosstalks". Additionally, we predict clusters of crosstalking
candidates, expanding from the anecdotally known examples in model organisms.
The tools and results presented here can be used to guide experimental studies
towards a system-level understanding of two-component signaling.


  On the example of exceptional families of viruses we i) show the existence of
a completely new type of matter organization in nanoparticles, in which the
regions with a chiral pentagonal quasicrystalline order of protein positions
are arranged in a structure commensurate with the spherical topology and
dodecahedral geometry, ii) generalize the classical theory of quasicrystals
(QCs) to explain this organization, and iii) establish the relation between
local chiral QC order and nonzero curvature of the dodecahedral capsid faces.


  The method of generalized modeling has been applied successfully in many
different contexts, particularly in ecology and systems biology. It can be used
to analyze the stability and bifurcations of steady-state solutions. Although
many dynamical systems in mathematical biology exhibit steady-state behaviour
one also wants to understand nonlocal dynamics beyond equilibrium points. In
this paper we analyze predator-prey dynamical systems and extend the method of
generalized models to periodic solutions. First, we adapt the equilibrium
generalized modeling approach and compute the unique Floquet multiplier of the
periodic solution which depends upon so-called generalized elasticity and scale
functions. We prove that these functions also have to satisfy a flow on
parameter (or moduli) space. Then we use Fourier analysis to provide computable
conditions for stability and the moduli space flow. The final stability
analysis reduces to two discrete convolutions which can be interpreted to
understand when the predator-prey system is stable and what factors enhance or
prohibit stable oscillatory behaviour. Finally, we provide a sampling algorithm
for parameter space based on nonlinear optimization and the Fast Fourier
Transform which enables us to gain a statistical understanding of the stability
properties of periodic predator-prey dynamics.


  The response of complex networks to perturbations is of utmost importance in
areas as diverse as ecosystem management, emergency response, and cell
reprogramming. A fundamental property of networks is that the perturbation of
one node can affect other nodes, in a process that may cause the entire or
substantial part of the system to change behavior and possibly collapse. Recent
research in metabolic and food-web networks has demonstrated the concept that
network damage caused by external perturbations can often be mitigated or
reversed by the application of compensatory perturbations. Compensatory
perturbations are constrained to be physically admissible and amenable to
implementation on the network. However, the systematic identification of
compensatory perturbations that conform to these constraints remains an open
problem. Here, we present a method to construct compensatory perturbations that
can control the fate of general networks under such constraints. Our approach
accounts for the full nonlinear behavior of real complex networks and can bring
the system to a desirable target state even when this state is not directly
accessible. Applications to genetic networks show that compensatory
perturbations are effective even when limited to a small fraction of all nodes
in the network and that they are far more effective when limited to the
highest-degree nodes. The approach is conceptually simple and computationally
efficient, making it suitable for the rescue, control, and reprogramming of
large complex networks in various domains.


  A poorly understood step in the transition from a chemical to a biological
world is the emergence of self-replicating molecular systems. We study how a
precursor for such a replicator might arise in a hydrothermal RNA reactor,
which accumulates longer sequences from unbiased monomer influx and random
ligation. In the reactor, intra- and inter-molecular basepairing locally
protects from random cleavage. By analyzing stochastic simulations, we find
temporal sequence correlations that constitute a signature of information
transmission, weaker but of the same form as in a true replicator.


  The introduction and emergence of agriculture into Eastern North America
(ENA) and Europe proceeded very differently in both subcontinents: it varied in
timing, speed, and mechanism. Common to both regions, agricultural subsistence
profited from the introduction of major staple crops which had been
domesticated elsewhere; in both regions, the temperate climate and originally
predominant forest vegetation provided an environmental context conducive to
agriculture. To understand the different paths to agricultural subsistence, an
integrated view of technological innovation, domestication and exchange of
domesticates, migration and trade is required within the constraints imposed by
the environmental context and geography. This study makes use of a numerical
model of regional socio-technological evolution which builds on adaptation of
important characteristics of prehistoric societies such as technology and
subsistence diversity. With this Global Land Use and Technological Evolution
Simulator (GLUES), regional transitions from foraging to agropastoralism, and
associated land use and demographic changes are realistically hindcasted. I
show that the model is capable of explaining the different timing, speed, and
transition mechanisms: Europe received a large package of foreign domesticates
and converted rapidly to agriculture, this fast transition can be best
explained with demic and cultural diffusion followed by fast adoption by
resident foragers. In contrast, ENA trajectories show a more gradual
transition. Hunting-gathering and agropastoral life style coexisted for a long
time and agriculture was adopted slowly into the existing subsistence scheme.


  Despite their importance for the spread of zoonotic diseases, our
understanding of the dynamical aspects characterizing the movements of farmed
animal populations remains limited as these systems are traditionally studied
as static objects and through simplified approximations. By leveraging on the
network science approach, here we are able for the first time to fully analyze
the longitudinal dataset of Italian cattle movements that reports the mobility
of individual animals among farms on a daily basis. The complexity and
inter-relations between topology, function and dynamical nature of the system
are characterized at different spatial and time resolutions, in order to
uncover patterns and vulnerabilities fundamental for the definition of targeted
prevention and control measures for zoonotic diseases. Results show how the
stationarity of statistical distributions coexists with a strong and
non-trivial evolutionary dynamics at the node and link levels, on all
timescales. Traditional static views of the displacement network hide important
patterns of structural changes affecting nodes' centrality and farms' spreading
potential, thus limiting the efficiency of interventions based on partial
longitudinal information. By fully taking into account the longitudinal
dimension, we propose a novel definition of dynamical motifs that is able to
uncover the presence of a temporal arrow describing the evolution of the system
and the causality patterns of its displacements, shedding light on mechanisms
that may play a crucial role in the definition of preventive actions.


  We present Molecular Dynamics simulations of a single stranded unprotonated
DNA i-motif in explicit solvent. Our results indicate that the native structure
in non-acidic solution at 300 K is unstable and completely vanishes on a time
scale up to 10 ns. Two unfolding mechanisms with decreasing connectivity
between the initially interacting nucleobases can be identified where one
pathway is characterized as entropically more favorable. The entropic
preference can be mainly explained by strong water ordering effects due to
hydrogen bonds for several occurring structures along the pathways. Finally we
are able to indicate via free energy calculations the most stable
configurations belonging to distinct hairpin structures in good agreement to
experimental results.


  Three, interrelated biologically-relevant examples of biased random walks are
presented: (1) A model for DNA melting, modelled as DNA unzipping, which
provides a way to illustrate the role of the Boltzmann factor in a venue
well-known to biology and pre-medical students; (2) the activity of helicase
motor proteins in unzipping double-stranded DNA, for example, at the
replication fork, which is an example of a Brownian ratchet; (3) force
generation by actin polymerization, which is another Brownian ratchet, and for
which the force and actin-concentration dependence of the velocity of actin
polymerization is determined.


  It is well-established that including spatial structure and stochastic noise
in models for predator-prey interactions invalidates the classical
deterministic Lotka-Volterra picture of neutral population cycles. In contrast,
stochastic models yield long-lived, but ultimately decaying erratic population
oscillations, which can be understood through a resonant amplification
mechanism for density fluctuations. In Monte Carlo simulations of spatial
stochastic predator-prey systems, one observes striking complex spatio-temporal
structures. These spreading activity fronts induce persistent correlations
between predators and prey. In the presence of local particle density
restrictions (finite prey carrying capacity), there exists an extinction
threshold for the predator population. The accompanying continuous
non-equilibrium phase transition is governed by the directed-percolation
universality class. We employ field-theoretic methods based on the Doi-Peliti
representation of the master equation for stochastic particle interaction
models to (i) map the ensuing action in the vicinity of the absorbing state
phase transition to Reggeon field theory, and (ii) to quantitatively address
fluctuation-induced renormalizations of the population oscillation frequency,
damping, and diffusion coefficients in the species coexistence phase.


  Cellular behavior is governed by gene regulatory processes that are
intrinsically dynamic and nonlinear, and are subject to non-negligible amounts
of random fluctuations. Such conditions are ubiquitous in physical systems,
where they have been studied for decades using the tools of statistical and
nonlinear physics. The goal of this review is to show how approaches
traditionally used in physics can help in reaching a systems-level
understanding of living cells. To that end, we present an overview of the
dynamical phenomena exhibited by genetic circuits and their functional
significance. We also describe the theoretical and experimental approaches that
are being used to unravel the relationship between circuit structure and
function in dynamical cellular processes under the influence of noise, both at
the single-cell level and in cellular populations, where intercellular coupling
plays an important role.


  Nested structure, which is non-random, controls cooperation dynamics and
biodiversity in plant-animal mutualistic networks. This structural pattern has
been explained in a static (non-growth) network models. However, evolutionary
processes might also influence the formation of such a structural pattern. We
thereby propose an evolving network model for plant-animal interactions and
show that non-random patterns such as nested structure and heterogeneous
connectivity are both qualitatively and quantitatively predicted through simple
evolutionary processes. This finding implies that network models can be
simplified by considering evolutionary processes, and also that another
explanation exists for the emergence of non-random patterns and might provide
more comprehensible insights into the formation of plant-animal mutualistic
networks from the evolutionary perspective.


  Living organisms produce metabolites of many types via their metabolisms.
Especially, flavonoids, a kind of secondary metabolites, of plant species are
interesting examples. Since plant species are believed to have specific
flavonoids with respect to diverse environment, elucidation of design
principles of metabolite distributions across plant species is important to
understand metabolite diversity and plant evolution. In the previous work, we
found heterogeneous connectivity in metabolite distributions, and proposed a
simple model to explain a possible origin of heterogeneous connectivity. In
this paper, we show further structural properties in the metabolite
distribution among families inspired by analogy with plant-animal mutualistic
networks: nested structure and modular structure. An earlier model represents
that these structural properties in bipartite relationships are determined
based on traits of elements and external factors. However, we find that the
architecture of metabolite distributions is described by simple evolution
processes without trait-based mechanisms by comparison between our model and
the earlier model. Our model can better predict nested structure and modular
structure in addition to heterogeneous connectivity both qualitatively and
quantitatively. This finding implies an alternative possible origin of these
structural properties, and suggests simpler formation mechanisms of metabolite
distributions across plant species than expected.


  We examine temperature dependent picosecond dynamics as a function of
structure and function for lysozyme and cytochrome c using temperature
dependent terahertz permittivity measurements. A double Arrhenius temperature
dependence with activation energies E1 ~ 0.1 kJ/mol and E2 ~10 kJ/mol fits the
native state response. The higher activation energy is consistent with the
so-called protein dynamical transition associated with beta relaxations at the
solvent-protein interface. The lower activation energy is consistent with
correlated structural motions. When the structure is removed by denaturing the
lower activation energy process is no longer present. Additionally the lower
activation energy process is diminished with ligand binding, but not for
changes in internal oxidation state. We suggest that the lower energy
activation process is associated with collective structural motions that are no
longer accessible with denaturing or binding.


  Self organization mechanisms are essential for the cytoskeleton to adapt to
the requirements of living cells. They rely on the intricate interplay of
cytoskeletal filaments, crosslinking proteins and molecular motors. Here we
present an in vitro minimal model system consisting of actin filaments, fascin
and myosin-II filaments exhibiting pulsative collective long range dynamics.
The reorganizations in the highly dynamic steady state of the active gel are
characterized by alternating periods of runs and stalls resulting in a
superdiffusive dynamics of the network's constituents. They are dominated by
the complex competition of crosslinking molecules and motor filaments in the
network: Collective dynamics are only observed if the relative strength of the
binding of myosin-II filaments to the actin network allows exerting high enough
forces to unbind actin/fascin crosslinks. The feedback between structure
formation and dynamics can be resolved by combining these experiments with
phenomenological simulations based on simple interaction rules.


  Motivation: Identifying the molecular pathways more prone to disruption
during a pathological process is a key task in network medicine and, more in
general, in systems biology.
  Results: In this work we propose a pipeline that couples a machine learning
solution for molecular profiling with a recent network comparison method. The
pipeline can identify changes occurring between specific sub-modules of
networks built in a case-control biomarker study, discriminating key groups of
genes whose interactions are modified by an underlying condition. The proposal
is independent from the classification algorithm used. Three applications on
genomewide data are presented regarding children susceptibility to air
pollution and two neurodegenerative diseases: Parkinson's and Alzheimer's.
  Availability: Details about the software used for the experiments discussed
in this paper are provided in the Appendix.


  There is great interest in the dynamics of health behaviors in social
networks and how they affect collective public health outcomes, but measuring
population health behaviors over time and space requires substantial resources.
Here, we use publicly available data from 101,853 users of online social media
collected over a time period of almost six months to measure the
spatio-temporal sentiment towards a new vaccine. We validated our approach by
identifying a strong correlation between sentiments expressed online and CDC-
estimated vaccination rates by region. Analysis of the network of opinionated
users showed that information flows more often between users who share the same
sentiments - and less often between users who do not share the same sentiments
- than expected by chance alone. We also found that most communities are
dominated by either positive or negative sentiments towards the novel vaccine.
Simulations of infectious disease transmission show that if clusters of
negative vaccine sentiments lead to clusters of unprotected individuals, the
likelihood of disease outbreaks are greatly increased. Online social media
provide unprecedented access to data allowing for inexpensive and efficient
tools to identify target areas for intervention efforts and to evaluate their
effectiveness.


  We explore the connection between a stochastic simulation model and an
ordinary differential equations (ODEs) model of the dynamics of an excitable
gene circuit that exhibits noise-induced oscillations. Near a bifurcation point
in the ODE model, the stochastic simulation model yields behavior dramatically
different from that predicted by the ODE model. We analyze how that behavior
depends on the gene copy number and find very slow convergence to the large
number limit near the bifurcation point. The implications for understanding the
dynamics of gene circuits and other birth-death dynamical systems with small
numbers of constituents are discussed.


  Recent work has discussed the importance of multiplicative closure for the
Markov models used in phylogenetics. For continuous-time Markov chains, a
sufficient condition for multiplicative closure of a model class is ensured by
demanding that the set of rate-matrices belonging to the model class form a Lie
algebra. It is the case that some well-known Markov models do form Lie algebras
and we refer to such models as "Lie Markov models". However it is also the case
that some other well-known Markov models unequivocally do not form Lie
algebras. In this paper, we will discuss how to generate Lie Markov models by
demanding that the models have certain symmetries under nucleotide
permutations. We show that the Lie Markov models include, and hence provide a
unifying concept for, "group-based" and "equivariant" models. For each of two,
three and four character states, the full list of Lie Markov models with
maximal symmetry is presented and shown to include interesting examples that
are neither group-based nor equivariant. We also argue that our scheme is
pleasing in the context of applied phylogenetics, as, for a given symmetry of
nucleotide substitution, it provides a natural hierarchy of models with
increasing number of parameters.


  Proteins with multiple binding sites play important roles in cell signaling
systems by nucleating protein complexes in which, for example, enzymes and
substrates are co-localized. Proteins that specialize in this function are
called by a variety names, including adapter, linker and scaffold.
Scaffold-mediated nucleation of protein complexes can be either constitutive or
induced. Induced nucleation is commonly mediated by a docking site on a
scaffold that is activated by phosphorylation. Here, by considering minimalist
mathematical models, which recapitulate scaffold effects seen in more
mechanistically detailed models, we obtain analytical and numerical results
that provide insights into scaffold function. These results elucidate how
recruitment of a pair of ligands to a scaffold depends on the concentrations of
the ligands, on the binding constants for ligand-scaffold interactions, on
binding cooperativity, and on the milieu of the scaffold, as ligand recruitment
is affected by competitive ligands and decoy receptors. For the case of a
bivalent scaffold, we obtain an expression for the unique scaffold
concentration that maximally recruits a pair of monovalent ligands. Through
simulations, we demonstrate that a bivalent scaffold can nucleate distinct sets
of ligands to equivalent extents when the scaffold is present at different
concentrations. Thus, the function of a scaffold can potentially change
qualitatively with a change in copy number. We also demonstrate how a scaffold
can change the catalytic efficiency of an enzyme and the sensitivity of the
rate of reaction to substrate concentration. The results presented here should
be useful for understanding scaffold function and for engineering scaffolds to
have desired properties.


  We investigate unfolding/folding force kinetics in DNA hairpins exhibiting
two and three states with newly designed short dsDNA handles (29 bp) using
optical tweezers. We show how the higher stiffness of the molecular setup
moderately enhances the signal-to-noise ratio (SNR) in hopping experiments as
compared to conventional long handles constructs (approximately 700 bp). The
shorter construct results in a signal of higher SNR and slower
folding/unfolding kinetics, thereby facilitating the detection of otherwise
fast structural transitions. A novel analysis of the elastic properties of the
molecular setup, based on high-bandwidth measurements of force fluctuations
along the folded branch, reveals that the highest SNR that can be achieved with
short handles is potentially limited by the marked reduction of the effective
persistence length and stretch modulus of the short linker complex.


  Turnover of regulatory sequence and function is an important part of
molecular evolution. But what are the modes of sequence evolution leading to
rapid formation and loss of regulatory sites? Here, we show that a large
fraction of neighboring transcription factor binding sites in the fly genome
have formed from a common sequence origin by local duplications. This mode of
evolution is found to produce regulatory information: duplications can seed new
sites in the neighborhood of existing sites. Duplicate seeds evolve
subsequently by point mutations, often towards binding a different factor than
their ancestral neighbor sites. These results are based on a statistical
analysis of 346 cis-regulatory modules in the Drosophila melanogaster genome,
and a comparison set of intergenic regulatory sequence in Saccharomyces
cerevisiae. In fly regulatory modules, pairs of binding sites show
significantly enhanced sequence similarity up to distances of about 50 bp. We
analyze these data in terms of an evolutionary model with two distinct modes of
site formation: (i) evolution from independent sequence origin and (ii)
divergent evolution following duplication of a common ancestor sequence. Our
results suggest that pervasive formation of binding sites by local sequence
duplications distinguishes the complex regulatory architecture of higher
eukaryotes from the simpler architecture of unicellular organisms.


  We consider a realization of the XOR logic gate in a process biocatalyzed by
an enzyme (here horseradish peroxidase: HRP), the function of which can be
inhibited by a substrate (hydrogen peroxide for HRP), when the latter is
inputted at large enough concentrations. A model is developed for describing
such systems in an approach suitable for evaluation of the analog noise
amplification properties of the gate. The obtained data are fitted for gate
quality evaluation within the developed model, and we discuss aspects of
devising XOR gates for functioning in "biocomputing" systems utilizing
biomolecules for information processing.


  Background: Zipf's law and Heaps' law are two representatives of the scaling
concepts, which play a significant role in the study of complexity science. The
coexistence of the Zipf's law and the Heaps' law motivates different
understandings on the dependence between these two scalings, which is still
hardly been clarified.
  Methodology/Principal Findings: In this article, we observe an evolution
process of the scalings: the Zipf's law and the Heaps' law are naturally shaped
to coexist at the initial time, while the crossover comes with the emergence of
their inconsistency at the larger time before reaching a stable state, where
the Heaps' law still exists with the disappearance of strict Zipf's law. Such
findings are illustrated with a scenario of large-scale spatial epidemic
spreading, and the empirical results of pandemic disease support a universal
analysis of the relation between the two laws regardless of the biological
details of disease. Employing the United States(U.S.) domestic air
transportation and demographic data to construct a metapopulation model for
simulating the pandemic spread at the U.S. country level, we uncover that the
broad heterogeneity of the infrastructure plays a key role in the evolution of
scaling emergence.
  Conclusions/Significance: The analyses of large-scale spatial epidemic
spreading help understand the temporal evolution of scalings, indicating the
coexistence of the Zipf's law and the Heaps' law depends on the collective
dynamics of epidemic processes, and the heterogeneity of epidemic spread
indicates the significance of performing targeted containment strategies at the
early time of a pandemic disease.


  Unraveling the evolutionary forces shaping bacterial diversity can today be
tackled using a growing amount of genomic data. While the genome of eukaryotes
is highly stable, bacterial genomes from cells of the same species highly vary
in gene content. This huge variation in gene content led to the concepts of the
distributed genome of bacteria and their pangenome (Tettelin et al.,2005;
Ehrlich et al.,2005). We present a population genetic model for gene content
evolution which accounts for several mechanisms. Gene uptake from the
environment is modeled by events of gene gain along the genealogical tree
relating the population. Pseudogenization may lead to deletion of genes and is
incoporated by gene loss. These two mechanisms were studied by Huson and Steel
(2004) using a fixed phylogenetic tree. Taking the random genealogy given by
the coalescent (Kingman, 1982; Hudson, 1983), we studied the resulting genomic
diversity already in Baumdicker et al. (2010). In the present paper, we extend
the model in order to incorporate events of interspecies horizontal gene
transfer. Within this model, we derive expectations for the gene frequency
spectrum and other quantities of interest.


  Can we objectively distinguish chemical systems that are able to process
meaningful information from those that are not suitable for information
processing? Here, we present a formal method to assess the semantic capacity of
a chemical reaction network. The semantic capacity of a network can be measured
by analyzing the capability of the network to implement molecular codes. We
analyzed models of real chemical systems (Martian atmosphere chemistry and
various combustion chemistries), bio-chemical systems (gene expression, gene
translation, and phosphorylation signaling cascades), as well as an artificial
chemistry and random networks. Our study suggests that different chemical
systems posses different semantic capacities. Basically no semantic capacity
was found in the atmosphere chemistry of Mars and all studied combustion
chemistries, as well as in highly connected random networks, i.e., with these
chemistries molecular codes cannot be implemented. High semantic capacity was
found in the bio-chemical systems, as well as in random networks where the
number of second order reactions is at the number of species. Hypotheses
concern the origin and evolution of life. We conclude that our approach can be
applied to evaluate the information processing capabilities of a chemical
system and may thus be a useful tool to understand the origin and evolution of
meaningful information, e.g., at the origin of life.


  Numerous transcription factors self-assemble into different order oligomeric
species in a way that is actively regulated by the cell. Until now, no general
functional role has been identified for this widespread process. Here we
capture the effects of modulated self-assembly in gene expression with a novel
quantitative framework. We show that this mechanism provides precision and
flexibility, two seemingly antagonistic properties, to the sensing of diverse
cellular signals by systems that share common elements present in transcription
factors like p53, NF-kappaB, STATs, Oct, and RXR. Applied to the nuclear
hormone receptor RXR, this framework accurately reproduces a broad range of
classical, previously unexplained, sets of gene expression data and
corroborates the existence of a precise functional regime with flexible
properties that can be controlled both at a genome-wide scale and at the
individual promoter level.


  We develop a theory of enzyme catalysis within biological cells where the
substrate concentration [S](t) is time dependent, in contrast to the
Michaelis-Menten theory that assumes a steady state. We find that the time
varying concentration can combine, in a non-linear way, with the ruggedness of
the free energy landscape of enzymes (discovered both in single molecule
studies and in simulations) to provide a highly efficient switch (or,
bifurcation) between two catalytically active states, at a critical substrate
concentration. This allows a dynamic control of product synthesis in cell.


  The flexibility in the structure of calmodulin (CaM) allows its binding to
over 300 target proteins in the cell. To investigate the structure-function
relationship of CaM, we combined methods of computer simulation and experiments
based on circular dichroism (CD) to investigate the structural characteristics
of CaM that influence its target recognition in crowded cell-like conditions.
We developed a unique multiscale solution of charges computed from quantum
chemistry, together with protein reconstruction, coarse-grained molecular
simulations, and statistical physics, to represent the charge distribution in
the transition from apoCaM to holoCaM upon calcium binding. Computationally, we
found that increased levels of macromolecular crowding, in addition to calcium
binding and ionic strength typical of that found inside cells, can impact the
conformation, helicity and the EF hand orientation of CaM. Because EF hand
orientation impacts the affinity of calcium binding and the specificity of
CaM's target selection, our results may provide unique insight into
understanding the promiscuous behavior of calmodulin in target selection inside
cells.


  Traveling waves are ubiquitous in nature and control the speed of many
important dynamical processes, including chemical reactions, epidemic
outbreaks, and biological evolution. Despite their fundamental role in complex
systems, traveling waves remain elusive because they are often dominated by
rare fluctuations in the wave tip, which have defied any rigorous analysis so
far. Here, we show that by adjusting nonlinear model details, noisy traveling
waves can be solved exactly. The moment equations of these tuned models are
closed and have a simple analytical structure resembling the deterministic
approximation supplemented by a nonlocal cutoff term. The peculiar form of the
cutoff shapes the noisy edge of traveling waves and is critical for the correct
prediction of the wave speed and its fluctuations. Our approach is illustrated
and benchmarked using the example of fitness waves arising in simple models of
microbial evolution, which are highly sensitive to number fluctuations. We
demonstrate explicitly how these models can be tuned to account for finite
population sizes and determine how quickly populations adapt as a function of
population size and mutation rates. More generally, our method is shown to
apply to a broad class of models, in which number fluctuations are generated by
branching processes. Because of this versatility, the method of model tuning
may serve as a promising route toward unraveling universal properties of
complex discrete particle systems.


  Human micro RNAs (miRNAs) target about 90% of the coding genes and form a
complex regulatory network. We study the community structure of the miRNA
co-target network considering miRNAs as the nodes which are connected by
weighted links. The weight of link that connects a pair of miRNAs denote the
total number of common transcripts targeted by that pair. We argue that the
network consists of about 74 modules, quite similar to the components (or
clusters) obtained earlier [Online J Bioinformatics, 10,280 ], indicating that
the components of the miRNA co-target network are self organized in a way to
maximize the modularity.


  In this note, we focus on a selection model problem: a mono-exponential model
versus a bi-exponential one. This is done in the biological context of living
cells, where small data are available. Classical statistics are revisited to
improve existing results. Some unavoidable limits are also pointed out.


  The widespread popularity of replica exchange and expanded ensemble
algorithms for simulating complex molecular systems in chemistry and biophysics
has generated much interest in enhancing phase space mixing of these protocols,
thus improving their efficiency. Here, we demonstrate how both of these classes
of algorithms can be considered a form of Gibbs sampling within a Markov chain
Monte Carlo (MCMC) framework. While the update of the conformational degrees of
freedom by Metropolis Monte Carlo or molecular dynamics unavoidably generates
correlated samples, we show how judicious updating of the thermodynamic state
indices---corresponding to thermodynamic parameters such as temperature or
alchemical coupling variables---associated with these configurations can
substantially increase mixing while still sampling from the desired
distributions. We show how state update methods in common use lead to
suboptimal mixing, and present some simple, inexpensive alternatives that can
increase mixing of the overall Markov chain, reducing simulation times
necessary to obtain estimates of the desired precision. These improved schemes
are demonstrated for several common applications, including an alchemical
expanded ensemble simulation, parallel tempering, and multidimensional replica
exchange umbrella sampling.