  The functioning of animal as well as human societies fundamentally relies on
cooperation. Yet, defection is often favorable for the selfish individual, and
social dilemmas arise. Selection by individuals' fitness, usually the basic
driving force of evolution, quickly eliminates cooperators. However, evolution
is also governed by fluctuations that can be of greater importance than fitness
differences, and can render evolution effectively neutral. Here, we investigate
the effects of selection versus fluctuations in social dilemmas. By studying
the mean extinction times of cooperators and defectors, a variable sensitive to
fluctuations, we are able to identify and quantify an emerging 'edge of neutral
evolution' that delineates regimes of neutral and Darwinian evolution. Our
results reveal that cooperation is significantly maintained in the neutral
regimes. In contrast, the classical predictions of evolutionary game theory,
where defectors beat cooperators, are recovered in the Darwinian regimes. Our
studies demonstrate that fluctuations can provide a surprisingly simple way to
partly resolve social dilemmas. Our methods are generally applicable to
estimate the role of random drift in evolutionary dynamics.


  Haplotyping is the bioinformatics problem of predicting likely haplotypes
based on given genotypes. It can be approached using Gusfield's perfect
phylogeny haplotyping (PPH) method for which polynomial and linear time
algorithms exist. These algorithm use sophisticated data structures or do a
stepwise transformation of the genotype data into haplotype data and,
therefore, need a linear amount of space. We are interested in the exact
computational complexity of PPH and show that it can be solved
space-efficiently by an algorithm that needs only a logarithmic amount of
space. Together with the recently proved L-hardness of PPH, we establish
L-completeness. Our algorithm relies on a new characterization for PPH in terms
of bipartite graphs, which can be used both to decide and construct perfect
phylogenies for genotypes efficiently.


  The promotion of cooperation on spatial lattices is an important issue in
evolutionary game theory. This effect clearly depends on the update rule: it
diminishes with stochastic imitative rules whereas it increases with
unconditional imitation. To study the transition between both regimes, we
propose a new evolutionary rule, which stochastically combines unconditional
imitation with another imitative rule. We find that, surprinsingly, in many
social dilemmas this rule yields higher cooperative levels than any of the two
original ones. This nontrivial effect occurs because the basic rules induce a
separation of timescales in the microscopic processes at cluster interfaces.
The result is robust in the space of 2x2 symmetric games, on regular lattices
and on scale-free networks.


  The transport by molecular motors along cytoskeletal filaments is studied
theoretically in the presence of static defects. The movements of single motors
are described as biased random walks along the filament as well as binding to
and unbinding from the filament. Three basic types of defects are
distinguished, which differ from normal filament sites only in one of the
motors' transition probabilities. Both stepping defects with a reduced
probability for forward steps and unbinding defects with an increased
probability for motor unbinding strongly reduce the velocities and the run
lengths of the motors with increasing defect density. For transport by single
motors, binding defects with a reduced probability for motor binding have a
relatively small effect on the transport properties. For cargo transport by
motors teams, binding defects also change the effective unbinding rate of the
cargo particles and are expected to have a stronger effect.


  We have compared a recently developed module-based algorithm LeMoNe for
reverse-engineering transcriptional regulatory networks to a mutual information
based direct algorithm CLR, using benchmark expression data and databases of
known transcriptional regulatory interactions for Escherichia coli and
Saccharomyces cerevisiae. A global comparison using recall versus precision
curves hides the topologically distinct nature of the inferred networks and is
not informative about the specific subtasks for which each method is most
suited. Analysis of the degree distributions and a regulator specific
comparison show that CLR is 'regulator-centric', making true predictions for a
higher number of regulators, while LeMoNe is 'target-centric', recovering a
higher number of known targets for fewer regulators, with limited overlap in
the predicted interactions between both methods. Detailed biological examples
in E. coli and S. cerevisiae are used to illustrate these differences and to
prove that each method is able to infer parts of the network where the other
fails. Biological validation of the inferred networks cautions against
over-interpreting recall and precision values computed using incomplete
reference networks.


  Complex diseases are multifactorial traits caused by both genetic and
environmental factors. They represent the most part of human diseases and
include those with largest prevalence and mortality (cancer, heart disease,
obesity, etc.). Despite of a large amount of information that have been
collected about both genetic and environmental risk factors, there are
relatively few examples of studies on their interactions in epidemiological
literature. One reason can be the incomplete knowledge of the power of
statistical methods designed to search for risk factors and their interactions
in this data sets. An improving in this direction would lead to a better
understanding and description of gene-environment interaction. To this aim, a
possible strategy is to challenge the different statistical methods against
data sets where the underlying phenomenon is completely known and fully
controllable, like for example simulated ones. We present a mathematical
approach that models gene-environment interactions. By this method it is
possible to generate simulated populations having gene-environment interactions
of any form. We implemented a simple version of this model in a
Gene-Environment iNteraction Simulator (GENS), a tool designed to simulate
case-control data sets where a one gene-one environment interaction influences
the disease risk. The main effort has been to allow user to describe
characteristics of population by using standard epidemiological measures and to
implement constraints to make the simulator behavior biologically meaningful.


  We study vibrational thermodynamic stability of small-world oscillator
networks, by relating the average mean-square displacement $S$ of oscillators
to the eigenvalue spectrum of the Laplacian matrix of networks. We show that
the cross-links suppress $S$ effectively and there exist two phases on the
small-world networks: 1) an unstable phase: when $p\ll1/N$, $S\sim N$; 2) a
stable phase: when $p\gg1/N$, $S\sim p^{-1}$, \emph{i.e.}, $S/N\sim
E_{cr}^{-1}$. Here, $p$ is the parameter of small-world, $N$ is the number of
oscillators, and $E_{cr}=pN$ is the number of cross-links. The results are
exemplified by various real protein structures that follow the same scaling
behavior $S/N\sim E_{cr}^{-1}$ of the stable phase. We also show that it is the
"small-world" property that plays the key role in the thermodynamic stability
and is responsible for the universal scaling $S/N\sim E_{cr}^{-1}$, regardless
of the model details.


  We investigate the effects of strong number fluctuations on traveling waves
in the Fisher-Kolmogorov reaction-diffusion system. Our findings are in stark
contrast to the commonly used deterministic and weak-noise approximations. We
compute the wave velocity in one and two spatial dimensions, for which we find
a linear and a square-root dependence of the speed on the particle density.
Instead of smooth sigmoidal wave profiles, we observe fronts composed of a few
rugged kinks that diffuse, annihilate, and rarely branch; this dynamics leads
to power-law tails in the distribution of the front sizes.


  Aligning protein-protein interaction (PPI) networks of different species has
drawn a considerable interest recently. This problem is important to
investigate evolutionary conserved pathways or protein complexes across
species, and to help in the identification of functional orthologs through the
detection of conserved interactions. It is however a difficult combinatorial
problem, for which only heuristic methods have been proposed so far. We
reformulate the PPI alignment as a graph matching problem, and investigate how
state-of-the-art graph matching algorithms can be used for that purpose. We
differentiate between two alignment problems, depending on whether strict
constraints on protein matches are given, based on sequence similarity, or
whether the goal is instead to find an optimal compromise between sequence
similarity and interaction conservation in the alignment. We propose new
methods for both cases, and assess their performance on the alignment of the
yeast and fly PPI networks. The new methods consistently outperform
state-of-the-art algorithms, retrieving in particular 78% more conserved
interactions than IsoRank for a given level of sequence similarity.
Availability:http://cbio.ensmp.fr/proj/graphm\_ppi/, additional data and codes
are available upon request. Contact: jean-philippe.vert@mines-paristech.fr


  Nanopores spanning synthetic membranes have been used as key components in
proof-of-principle nanofluidic applications, particularly those involving
manipulation of biomolecules or sequencing of DNA. The only practical way of
manipulating charged macromolecules near nanopores is through a voltage
difference applied across the nanopore-spanning membrane. However, recent
experiments have shown that salt concentration gradients applied across
nanopores can also dramatically enhance charged particle capture from a low
concentration reservoir of charged molecules at one end of the nanopore. This
puzzling effect has hitherto eluded a physically consistent theoretical
explanation. Here, we propose an electrokinetic mechanism of this enhanced
capture that relies on the electrostatic potential near the pore mouth. For
long pores with diameter much greater than the local screening length, we
obtain accurate analytic expressions showing how salt gradients control the
local conductivity which can lead to increased local electrostatic potentials
and charged analyte capture rates. We also find that the attractive
electrostatic potential may be balanced by an outward, repulsive electroosmotic
flow (EOF) that can in certain cases conspire with the salt gradient to further
enhance the analyte capture rate.


  Populations of heterogeneous cells play an important role in many biological
systems. In this paper we consider systems where each cell can be modelled by
an ordinary differential equation. To account for heterogeneity, parameter
values are different among individual cells, subject to a distribution function
which is part of the model specification.
  Experimental data for heterogeneous cell populations can be obtained from
flow cytometric fluorescence microscopy. We present a heuristic approach to use
such data for estimation of the parameter distribution in the population. The
approach is based on generating simulation data for samples in parameter space.
By convex optimisation, a suitable probability density function for these
samples is computed.
  To evaluate the proposed approach, we consider artificial data from a simple
model of the tumor necrosis factor (TNF) signalling pathway. Its main
characteristic is a bimodality in the TNF response: a certain percentage of
cells undergoes apoptosis upon stimulation, while the remaining part stays
alive. We show how our modelling approach allows to identify the reasons that
underly the differential response.


  New checkable criteria for persistence of chemical reaction networks are
proposed, which extend and complement those obtained by the authors in previous
work. The new results allow the consideration of reaction rates which are
time-varying, thus incorporating the effects of external signals, and also
relax the assumption of existence of global conservation laws, thus allowing
for inflows (production) and outflows (degradation). For time-invariant
networks parameter-dependent conditions for persistence of certain classes of
networks are provided. As an illustration, two networks arising in the systems
biology literature are analyzed, namely a hypoxia and an apoptosis network.


  Lipid membranes regulate the flow of materials and information between cells
and their organelles. Further, lipid composition and morphology can play a key
role in regulating a variety of biological processes. For example, viral
uptake, plasma membrane tension regulation, and the formation of caveolae all
require the creation and control of groups of lipids that adopt specific
morphologies. In this paper, we use a simplified model mixture of lipids and
cholesterol to examine the interplay between lipid phase-separation and bilayer
morphology. We observe and theoretically analyze three main features of
phase-separated giant unilamellar vesicles. First, by tracking the motion of
`dimpled' domains, we measure repulsive, elastic interactions that create
short--range translational and orientational order, leading to a stable
distribution of domain sizes, and hence maintaining lateral heterogeneity on
relatively short length scales and long time scales. Second, we examine the
transition to `budded' domain morphologies, showing that the transition is
size-selective, and has two kinetic regimes, as revealed by a calculated phase
diagram. Finally, using observations of the interactions between dimpled and
budded domains, we build a theoretical framework with an elastic model that
maps the free energies and allowed transitions in domain morphology upon
coalescence, to serve as an interpretive tool for understanding the algebra of
domain morphology. In all three cases, the two major factors that regulate
domain morphology and morphological transitions are the domain size and
membrane tension.


  The search for sample-variable associations is an important problem in the
exploratory analysis of high dimensional data. Biclustering methods search for
sample-variable associations in the form of distinguished submatrices of the
data matrix. (The rows and columns of a submatrix need not be contiguous.) In
this paper we propose and evaluate a statistically motivated biclustering
procedure (LAS) that finds large average submatrices within a given real-valued
data matrix. The procedure operates in an iterative-residual fashion, and is
driven by a Bonferroni-based significance score that effectively trades off
between submatrix size and average value. We examine the performance and
potential utility of LAS, and compare it with a number of existing methods,
through an extensive three-part validation study using two gene expression
datasets. The validation study examines quantitative properties of biclusters,
biological and clinical assessments using auxiliary information, and
classification of disease subtypes using bicluster membership. In addition, we
carry out a simulation study to assess the effectiveness and noise sensitivity
of the LAS search procedure. These results suggest that LAS is an effective
exploratory tool for the discovery of biologically relevant structures in high
dimensional data. Software is available at https://genome.unc.edu/las/.


  It has been difficult to measure efficacy of T cell-based vaccines and to
correlate efficacy of CD8$^+$T cell responses with protection against viral
infections. In part, this difficulty is due to our poor understanding of the in
vivo efficacy of CD8$^+$T cells. Using a recently developed experimental method
of in vivo cytotoxicity we investigated quantitative aspects of killing of
peptide-pulsed targets by effector and memory CD8$^+$T cells, specific to three
epitopes of lymphocytic choriomeningitis virus (LCMV), in the mouse spleen. By
analyzing data on killing of targets with varying number of epitope-specific
effector and memory CD8$^+$T cells, we find that killing of targets by
effectors follows the law of mass-action, that is the death rate of
peptide-pulsed targets is proportional to the frequency of CTLs in the spleen.
In contrast, killing of targets by memory CD8$^+$T cells does not follow the
mass action law because the death rate of targets saturates at high frequencies
of memory CD8$^+$T cells. For both effector and memory cells, we also find no
support for a killing term that includes the decrease of the death rate of
targets with increasing target cell density. Importantly, we find that at low
CD8$^+$T cell frequencies, effector and memory CD8$^+$T cells on the per capita
basis are equally efficient at killing peptide-pulsed targets. Our framework
provides the guideline for the calculation of the level of memory CD8$^+$T
cells required to provide sterilizing protection against viral infection. Our
results thus form a basis for quantitative understanding of the process of
killing of virus-infected cells by T cell responses in tissues and can be used
to correlate the phenotype of vaccine-induced memory CD8 T cells with their
killing efficacy in vivo.


  Viral infections can be broadly divided into infections that are cleared from
the host (acute) and those that persist (chronic). Why some viruses establish
chronic infections while other do not is poorly understood. One possibility is
that the host's immune response is impaired during chronic infections and is
unable to clear the virus from the host. In this report we use a recently
proposed framework to estimate the per capita killing efficacy of CD8$^+$ T
cells, specific for the MT389 epitope of polyoma virus (PyV), which establishes
a chronic infection in mice. Surprisingly, the estimated per cell killing
efficacy of MT389-specific effector CD8$^+$ T cells during the acute phase of
the infection was very similar to the previously estimated efficacy of effector
CD8$^+$ T cells specific to lymphocytic choriomeningitis virus
(LCMV-Armstrong), which is cleared from the host. We also find that during the
chronic phase of the infection the killing efficacy of PyV-specific CD8$^+$ T
cells was only half of that of cells in the acute phase. This decrease in the
killing efficacy is again surprisingly similar to the change in the killing
efficacy of LCMV-specific CD8$^+$ T cells from the peak of the response to the
memory phase. Interestingly, we also find that PyV-specific CD8$^+$ T cells in
the chronic phase of the infection require lower doses of antigen to kill a
target cell. In summary, we find little support for the hypothesis that
persistence of infections is caused by inability of the host to mount an
efficient immune response, and that even in the presence of an efficient
CD8$^+$ T cell response, some viruses can still establish a persistent
infection.


  The knowledge of the Free Energy Landscape topology is the essential key to
understand many biochemical processes. The determination of the conformers of a
protein and their basins of attraction takes a central role for studying
molecular isomerization reactions. In this work, we present a novel framework
to unveil the features of a Free Energy Landscape answering questions such as
how many meta-stable conformers are, how the hierarchical relationship among
them is, or what the structure and kinetics of the transition paths are.
Exploring the landscape by molecular dynamics simulations, the microscopic data
of the trajectory are encoded into a Conformational Markov Network. The
structure of this graph reveals the regions of the conformational space
corresponding to the basins of attraction. In addition, handling the
Conformational Markov Network, relevant kinetic magnitudes as dwell times or
rate constants, and the hierarchical relationship among basins, complete the
global picture of the landscape. We show the power of the analysis studying a
toy model of a funnel-like potential and computing efficiently the conformers
of a short peptide, the dialanine, paving the way to a systematic study of the
Free Energy Landscape in large peptides.


  The spatiotemporal oscillations of the Min proteins in the bacterium
Escherichia coli play an important role in cell division. A number of different
models have been proposed to explain the dynamics from the underlying
biochemistry. Here, we extend a previously described discrete polymer model
from a deterministic to a stochastic formulation. We express the stochastic
evolution of the oscillatory system as a map from the probability distribution
of maximum polymer length in one period of the oscillation to the probability
distribution of maximum polymer length half a period later and solve for the
fixed point of the map with a combined analytical and numerical technique. This
solution gives a theoretical prediction of the distributions of both lengths of
the polar MinD zones and periods of oscillations -- both of which are
experimentally measurable. The model provides an interesting example of a
stochastic hybrid system that is, in some limits, analytically tractable.


  We study the spread of susceptible-infected-recovered (SIR) infectious
diseases where an individual's infectiousness and probability of recovery
depend on his/her "age" of infection. We focus first on early outbreak stages
when stochastic effects dominate and show that epidemics tend to happen faster
than deterministic calculations predict. If an outbreak is sufficiently large,
stochastic effects are negligible and we modify the standard ordinary
differential equation (ODE) model to accommodate age-of-infection effects. We
avoid the use of partial differential equations which typically appear in
related models. We introduce a "memoryless" ODE system which approximates the
true solutions. Finally, we analyze the transition from the stochastic to the
deterministic phase.


  Biochemical computing attempts to process information with biomolecules and
biological objects. In this work we review our results on analysis and
optimization of single biochemical logic gates based on enzymatic reactions,
and a network of three gates, for reduction of the "analog" noise buildup. For
a single gate, optimization is achieved by analyzing the enzymatic reactions
within a framework of kinetic equations. We demonstrate that using
co-substrates with much smaller affinities than the primary substrate, a
negligible increase in the noise output from the logic gate is obtained as
compared to the input noise. A network of enzymatic gates is analyzed by
varying selective inputs and fitting standardized few-parameters response
functions assumed for each gate. This allows probing of the individual gate
quality but primarily yields information on the relative contribution of the
gates to noise amplification. The derived information is then used to modify
experimental single gate and network systems to operate them in a regime of
reduced analog noise amplification.


  An investigation into the properties of a two dimensional (2D+1) system of
self propelled particles (known as boids) in free space is conducted using a
Lagrangian Individual-Based Model. A potential, associated with each boid is
specified and a Lagrangian is subsequently derived in order to obtain the
equations of motion for each particle in the flock. The Morse potential and the
Lennard-Jones potential, both well understood in atomic and molecular physics,
are specified. In contrast to the original model proposed by Vicsek [Phys. Rev.
Lett. 75, 1226 (1995)] systems are considered with open boundary conditions.
These two models successfully replicate the phases observed in Vicsek's
original model, as well as several other significant phases, providing a
realistic model of a wide range of flocking phenomena.


  The bistable gene regulatory switch controlling the transition from lysogeny
to lysis in bacteriophage lambda presents a unique challenge to quantitative
modeling. Despite extensive characterization of this regulatory network, the
origin of the extreme stability of the lysogenic state remains unclear. We have
constructed a stochastic model for this switch. Using Forward Flux Sampling
simulations, we show that this model predicts an extremely low rate of
spontaneous prophage induction in a recA mutant, in agreement with experimental
observations. In our model, the DNA loop formed by octamerization of CI bound
to the O_L and O_R operator regions is crucial for stability, allowing the
lysogenic state to remain stable even when a large fraction of the total CI is
depleted by nonspecific binding to genomic DNA. DNA looping also ensures that
the switch is robust to mutations in the order of the O_R binding sites. Our
results suggest that DNA looping can provide a mechanism to maintain a stable
lysogenic state in the face of a range of challenges including noisy gene
expression, nonspecific DNA binding and operator site mutations.


  We propose the following simple stochastic model for phylogenetic trees. New
types are born and die according to a birth and death chain. At each birth we
associate a fitness to the new type sampled from a fixed distribution. At each
death the type with the smallest fitness is killed. We show that if the birth
(i.e. mutation) rate is subcritical we get a phylogenetic tree consistent with
an influenza tree (few types at any given time and one dominating type lasting
a long time). When the birth rate is supercritical we get a phylogenetic tree
consistent with an HIV tree (many types at any given time, none lasting very
long).


  We analyze populations of Kuramoto oscillators with a particular distribution
of natural frequencies. Inspired by networks where there are two groups of
nodes with opposite behaviors, as for instance in power-grids where energy is
either generated or consumed at different locations, we assume that the
frequencies can take only two different values. Correlations between the value
of the frequency of a given node and its topological localization are
considered in both regular and random topologies. Synchronization is enhanced
when nodes are surrounded by nodes of the opposite frequency. We find
analytical estimations for the minimum value of the coupling strength between
oscillators that guarantees the achievement of a globally synchronized state,
getting a very good agreement with the numerical simulations.


  The effect of ultralow-frequency or static magnetic and electric fields on
biological processes is of huge interest for researchers due to the resonant
change of the intensity of biochemical reactions although the energy in such
fields is small. A simplified model to study the effect of the weak magnetic
and electrical fields on fluctuation of the random ionic currents in blood and
to solve the $k_BT$ problem in magnetobiology is suggested. The analytic
expression for the kinetic energy of the molecules dissolved in certain liquid
media is obtained. The values of the magnetic field leading to resonant effects
in capillaries are estimated. The numerical estimates showed that the resonant
values of the energy of molecular in the capillaries and aorta are different:
under identical conditions a molecule of the aorta gets $10^{-9}$ times less
energy than the molecules in blood capillaries. So the capillaries are very
sensitive to the resonant effect, with an approach to the resonant value of the
magnetic field strength, the average energy of the molecule localized in the
capillary is increased by several orders of magnitude as compared to its
thermal energy, this value of the energy is sufficient for the deterioration of
the chemical bonds.


  Constraint-based flux balance analysis (FBA) has proven successful in
predicting the flux distribution of metabolic networks in diverse environmental
conditions. FBA finds one of the alternate optimal solutions that maximizes the
biomass production rate. Almaas et al have shown that the flux distribution
follows a power law, and it is possible to associate with most metabolites two
reactions which maximally produce and consume a give metabolite, respectively.
This observation led to the concept of high-flux backbone (HFB) in metabolic
networks. In previous work, the HFB has been computed using a particular optima
obtained using FBA. In this paper, we investigate the conservation of HFB of a
particular solution for a given medium across different alternate optima and
near-optima in metabolic networks of E. coli and S. cerevisiae. Using flux
variability analysis (FVA), we propose a method to determine reactions that are
guaranteed to be in HFB regardless of alternate solutions. We find that the HFB
of a particular optima is largely conserved across alternate optima in E. coli,
while it is only moderately conserved in S. cerevisiae. However, the HFB of a
particular near-optima shows a large variation across alternate near-optima in
both organisms. We show that the conserved set of reactions in HFB across
alternate near-optima has a large overlap with essential reactions and
reactions which are both uniquely consuming (UC) and uniquely producing (UP).
Our findings suggest that the structure of the metabolic network admits a high
degree of redundancy and plasticity in near-optimal flow patterns enhancing
system robustness for a given environmental condition.


  In this work we demonstrate both experimentally and theoretically that the
analog noise generation by a single enzymatic logic gate can be dramatically
reduced to yield gate operation with virtually no input noise amplification.
This is achieved by exploiting the enzyme's specificity when using a
co-substrate that has a much lower affinity than the primary substrate. Under
these conditions, we obtain a negligible increase in the noise output from the
logic gate as compared to the input noise level. Experimental realizations of
the AND logic gate with the enzyme horseradish peroxidase using hydrogen
peroxide and two different co-substrates,
2,2'-azino-bis(3-ethylbenzthiazoline-6-sulphonic acid) (ABTS) and ferrocyanide,
with vastly different rate constants confirmed our general theoretical
conclusions.


  The high computational cost of carrying out molecular dynamics simulations of
even small-size proteins is a major obstacle in the study, at atomic detail and
in explicit solvent, of the physical mechanism which is at the basis of the
folding of proteins. Making use of a biasing algorithm, based on the principle
of the ratchet-and-pawl, we have been able to calculate eight folding
trajectories (to an RMSD between 1.2A and 2.5A) of the B1 domain of protein G
in explicit solvent without the need of high-performance computing. The
simulations show that in the denatured state there is a complex network of
cause-effect relationships among contacts, which results in a rather
hierarchical folding mechanism. The network displays few local and nonlocal
native contacts which are cause of most of the others, in agreement with the
NOE signals obtained in mildly-denatured conditions. Also nonnative contacts
play an active role in the folding kinetics. The set of conformations
corresponding to the transition state display phi-values with a correlation
coefficient of 0.69 with the experimental ones. They are structurally quite
homogeneous and topologically native-like, although some of the side chains and
most of the hydrogen bonds are not in place.


  We calculate the optimality of a doublet precursor to the canonical genetic
code with respect to mitigating the effects of point mutations and compare our
results to corresponding ones for the canonical genetic code. We find that the
proposed precursor has much less optimality than that of the canonical code.
Our results render unlikely the notion that the doublet precursor was an
intermediate state in the evolution of the canonical genetic code. These
findings support the notion that code optimality reflects evolutionary
dynamics, and that if such a doublet code originally had a biochemical
significance, it arose before the emergence of translation.


  When someone mentions the name of a known person we immediately recall her
face and possibly many other traits. This is because we possess the so-called
associative memory, that is the ability to correlate different memories to the
same fact or event. Associative memory is such a fundamental and encompassing
human ability (and not just human) that the network of neurons in our brain
must perform it quite easily. The question is then whether electronic neural
networks (electronic schemes that act somewhat similarly to human brains) can
be built to perform this type of function. Although the field of neural
networks has developed for many years, a key element, namely the synapses
between adjacent neurons, has been lacking a satisfactory electronic
representation. The reason for this is that a passive circuit element able to
reproduce the synapse behaviour needs to remember its past dynamical history,
store a continuous set of states, and be "plastic" according to the
pre-synaptic and post-synaptic neuronal activity. Here we show that all this
can be accomplished by a memory-resistor (memristor for short). In particular,
by using simple and inexpensive off-the-shelf components we have built a
memristor emulator which realizes all required synaptic properties. Most
importantly, we have demonstrated experimentally the formation of associative
memory in a simple neural network consisting of three electronic neurons
connected by two memristor-emulator synapses. This experimental demonstration
opens up new possibilities in the understanding of neural processes using
memory devices, an important step forward to reproduce complex learning,
adaptive and spontaneous behaviour with electronic neural networks.


  Disease occurs due to aberrant expression of genes and modulation of the
biological pathways along which they lie. Inference of activated gene pathways,
using gene expression data during disease progression, is an important problem.
In this work, we have developed a generalizable framework for the
identification of interacting pathways while incorporating biological realism,
using functional data analysis and manifold embedding techniques. Additionally,
we have also developed a new method to query for the differential co-ordinated
activity of any desired pathway during disease progression. The methods
developed in this work can be generalized to any conditions of interest.


  A Charge Transport (CT) mechanism has been proposed in several papers (e.g.,
Yavin, et al. PNAS, v102 3546 (2005)) to explain the localization of Base
Excision Repair (BER) enzymes to lesions on DNA. The CT mechanism relies on
redox reactions of iron-sulfur cofactors that modify the enzyme's binding
affinity. These redox reactions are mediated by the DNA strand and involve the
exchange of electrons between BER enzymes along DNA. We propose a mathematical
model that incorporates enzyme binding/unbinding, electron transport, and
enzyme diffusion along DNA. Analysis of our model within a range of parameter
values suggests that the redox reactions can increase desorption of BER enzymes
not already bound to their targets, allowing the enzymes to be recycled, thus
accelerating the overall search process. This acceleration mechanism is most
effective when enzyme copy numbers and enzyme diffusivity along the DNA are
small. Under such conditions, we find that CT BER enzymes find their targets
more quickly than simple "passive" enzymes that simply attach to the DNA
without desorbing.


  We study the evolutionary dynamics of a haploid population of infinite size
recombining with a probability $r$ in a two locus model. Starting from a low
fitness locus, the population is evolved under mutation, selection and
recombination until a finite fraction of the population reaches the fittest
locus. An analytical method is developed to calculate the fixation time $T$ to
the fittest locus for various choices of epistasis. We find that (1) for
negative epistasis, $T$ decreases slowly for small $r$ but decays fast at
larger $r$ (2) for positive epistasis, $T$ increases linearly for small $r$ and
mildly for large $r$ (3) for compensatory mutation, $T$ diverges as a power law
with logarithmic corrections as the recombination fraction approaches a
critical value. Our calculations are seen to be in good agreement with the
exact numerical results.


  The plant microtubule cortical array is a striking feature of all growing
plant cells. It consists of a more or less homogeneously distributed array of
highly aligned microtubules connected to the inner side of the plasma membrane
and oriented transversely to the cell growth axis. Here we formulate a
continuum model to describe the origin of orientational order in such confined
arrays of dynamical microtubules. The model is based on recent experimental
observations that show that a growing cortical microtubule can interact through
angle dependent collisions with pre-existing microtubules that can lead either
to co-alignment of the growth, retraction through catastrophe induction or
crossing over the encountered microtubule. We identify a single control
parameter, which is fully determined by the nucleation rate and intrinsic
dynamics of individual microtubules. We solve the model analytically in the
stationary isotropic phase, discuss the limits of stability of this isotropic
phase, and explicitly solve for the ordered stationary states in a simplified
version of the model.


  Selection in a time-periodic environment is modeled via the two-player
replicator dynamics. For sufficiently fast environmental changes, this is
reduced to a multi-player replicator dynamics in a constant environment. The
two-player terms correspond to the time-averaged payoffs, while the three and
four-player terms arise from the adaptation of the morphs to their varying
environment. Such multi-player (adaptive) terms can induce a stable
polymorphism. The establishment of the polymorphism in partnership games
[genetic selection] is accompanied by decreasing mean fitness of the
population.


  The thermally induced denaturation of DNA in the presence of attractive solid
surface is studied. The two strands of DNA are modeled via two coupled flexible
chains without volume interactions. If the two strands are adsorbed on the
surface, the denaturation phase-transition disappears. Instead, there is a
smooth crossover to a weakly naturated state. Our second conclusion is that
even when the inter-strand attraction alone is too weak for creating a
naturated state at the given temperature, and also the surface-strand
attraction alone is too weak for creating an adsorbed state, the combined
effect of the two attractions can lead to a naturated and adsorbed state.


  We applied a periodic magnetic field of frequency 16 Hz and amplitude 16 nT
to a human I_Ks channel, expressed in a Xenopus oocyte and varied the membrane
depolarization between -100 mV and +100 mV. We observed a maximal increase of
about 9% in the potassium current at membrane depolarization between 0 mV and 8
mV (see Figure 3). A similar measurement of the potassium current in the KCNQ1
channel, expressed in an oocyte, gave a maximal increase of 16% at the same
applied magnetic field and membrane depolarization between -14 mV and -7 mV
(see Figure 4). We attribute this resonant behavior to stochastic resonance
between the thermal activations of the configuration of interacting ions in the
I_Ks channel over a low potential barrier inside the closed state of the
channel and the periodic electromotive force induced across the membrane by the
periodic magnetic field. The partial synchronization of the random jumps with
the periodic force changes the relative times spent on either side of the
barrier, thereby changing the open probability of the spontaneously gating open
channel. This, in turn, changes the conductance of the channel at the
particular depolarization level and frequency and is expressed in the
Hodgkin-Huxley equations as a bump at the given voltage in the
conductance-voltage relation. We integrate the modified Hodgkin-Huxley
equations for the current into the Luo-Rudy model of a Guinea pig ventricular
cardiac myocyte and obtain increased conductance during the plateau of the
action potential in the cell. This shortens both the action potential and the
cytosolic calcium concentration spike durations, lowers its amplitude,
increases cytosolic sodium, and lowers cytosolic potassium concentrations. The
shortening of the ventricular calcium signal shortens the QT period of the ECG.


  Estimation of division and death rates of lymphocytes in different conditions
is vital for quantitative understanding of the immune system. Deuterium, in the
form of deuterated glucose or heavy water, can be used to measure rates of
proliferation and death of lymphocytes in vivo. Inferring these rates from
labeling and delabeling curves has been subject to considerable debate with
different groups suggesting different mathematical models for that purpose. We
show that the three models that are most commonly used are in fact
mathematically identical and differ only in their interpretation of the
estimated parameters. By extending these previous models, we here propose a
more mechanistic approach for the analysis of data from deuterium labeling
experiments. We construct a model of "kinetic heterogeneity" in which the total
cell population consists of many sub-populations with different rates of cell
turnover. In this model, for a given distribution of the rates of turnover, the
predicted fraction of labeled DNA accumulated and lost can be calculated. Our
model reproduces several previously made experimental observations, such as a
negative correlation between the length of the labeling period and the rate at
which labeled DNA is lost after label cessation. We demonstrate the reliability
of the new explicit kinetic heterogeneity model by applying it to artificially
generated datasets, and illustrate its usefulness by fitting experimental data.
In contrast to previous models, the explicit kinetic heterogeneity model 1)
provides a mechanistic way of interpreting labeling data; 2) allows for a
non-exponential loss of labeled cells during delabeling, and 3) can be used to
describe data with variable labeling length.


  The study of biochemical pathways usually focuses on a small section of a
protein interactions network. Two distinct sources contribute to the noise in
such a system: intrinsic noise, inherent in the studied reactions, and
extrinsic noise generated in other parts of the network or in the environment.
We study the effect of extrinsic noise entering the system through a nonlinear
uptake reaction which acts as a nonlinear filter. Varying input noise intensity
varies the mean of the noise after the passage through the filter, which
changes the stability properties of the system. The steady-state displacement
due to small noise is independent on the kinetics of the system but it only
depends on the nonlinearity of the input function.
  For monotonically increasing and concave input functions such as the
Michaelis-Menten uptake rate, we give a simple argument based on the
small-noise expansion, which enables qualitative predictions of the
steady-state displacement only by inspection of experimental data: when weak
and rapid noise enters the system through a Michaelis-Menten reaction, then the
graph of the system's steady states vs. the mean of the input signal always
shifts to the right as noise intensity increases.
  We test the predictions on two models of lac operon, where TMG/lactose uptake
is driven by a Michaelis-Menten enzymatic process. We show that as a
consequence of the steady state displacement due to fluctuations in
extracellular TMG/lactose concentration the lac switch responds in an
asymmetric manner: as noise intensity increases, switching off lactose
metabolism becomes easier and switching it on becomes more difficult.


  Stochasticity and spatial heterogeneity are of great interest recently in
studying the spread of an infectious disease. The presented method solves an
inverse problem to discover the effectively decisive topology of a
heterogeneous network and reveal the transmission parameters which govern the
stochastic spreads over the network from a dataset on an infectious disease
outbreak in the early growth phase. Populations in a combination of
epidemiological compartment models and a meta-population network model are
described by stochastic differential equations. Probability density functions
are derived from the equations and used for the maximal likelihood estimation
of the topology and parameters. The method is tested with computationally
synthesized datasets and the WHO dataset on SARS outbreak.


  A repressilator consists of a loop made up of three repressively interacting
genes. We construct a hexagonal lattice with repressilators on each triangle,
and use this as a model system for multiple interacting feedback loops. Using
symmetry arguments and stability analysis we argue that the repressor-lattice
can be in a non-frustrated oscillating state with only three distinct phases.
If the system size is not commensurate with three, oscillating solutions of
several different phases are possible. As the strength of the interactions
between the nodes increases, the system undergoes many transitions, breaking
several symmetries. Eventually dynamical frustrated states appear, where the
temporal evolution is chaotic, even though there are no built-in frustrations.
Applications of the repressor-lattice to real biological systems, such as
tissues or biofilms, are discussed.


  There have been suggestions that heat caused by cerebral metabolic activity
may constrain mammalian brain evolution, architecture, and function. This
article investigates physical limits on brain wiring and corresponding changes
in brain temperature that are imposed by thermodynamics of heat balance
determined mainly by Na$^{+}$/K$^{+}$-ATPase, cerebral blood flow, and heat
conduction. It is found that even moderate firing rates cause significant
intracellular Na$^{+}$ build-up, and the ATP consumption rate associated with
pumping out these ions grows nonlinearly with frequency. Surprisingly, the
power dissipated by the Na$^{+}$/K$^{+}$ pump depends biphasically on
frequency, which can lead to the biphasic dependence of brain temperature on
frequency as well. Both the total power of sodium pumps and brain temperature
diverge for very small fiber diameters, indicating that too thin fibers are not
beneficial for thermal balance. For very small brains blood flow is not a
sufficient cooling mechanism deep in the brain. The theoretical lower bound on
fiber diameter above which brain temperature is in the operational regime is
strongly frequency dependent but finite due to synaptic depression. For normal
neurophysiological conditions this bound is at least an order of magnitude
smaller than average values of empirical fiber diameters, suggesting that
neuroanatomy of the mammalian brains operates in the thermodynamically safe
regime. Analytical formulas presented can be used to estimate average firing
rates in mammals, and relate their changes to changes in brain temperature,
which can have important practical applications. In general, activity in larger
brains is found to be slower than in smaller brains.


  Competitive birth-death processes often exhibit an oscillatory behavior. We
investigate a particular case where the oscillation cycles are marginally
stable on the mean-field level. An iconic example of such a system is the
Lotka-Volterra model of predator-prey competition. Fluctuation effects due to
discreteness of the populations destroy the mean-field stability and eventually
drive the system toward extinction of one or both species. We show that the
corresponding extinction time scales as a certain power-law of the population
sizes. This behavior should be contrasted with the extinction of models stable
in the mean-field approximation. In the latter case the extinction time scales
exponentially with size.


  Numerous studies over the past 30 years have suggested there is a causal
connection between the motion of the Sun through the Galaxy and terrestrial
mass extinctions or climate change. Proposed mechanisms include comet impacts
(via perturbation of the Oort cloud), cosmic rays and supernovae, the effects
of which are modulated by the passage of the Sun through the Galactic midplane
or spiral arms. Supposed periodicities in the fossil record, impact cratering
dates or climate proxies over the Phanerozoic (past 545 Myr) are frequently
cited as evidence in support of these hypotheses. This remains a controversial
subject, with many refutations and replies having been published. Here I review
both the mechanisms and the evidence for and against the relevance of
astronomical phenomena to climate change and evolution. This necessarily
includes a critical assessment of time series analysis techniques and
hypothesis testing. Some of the studies have suffered from flaws in
methodology, in particular drawing incorrect conclusions based on ruling out a
null hypothesis. I conclude that there is little evidence for intrinsic
periodicities in biodiversity, impact cratering or climate on timescales of
tens to hundreds of Myr. Furthermore, Galactic midplane and spiral arm
crossings seem to have little or no impact on biological or climate variation
above background level. (truncated)


  Bacteria communicate using secreted chemical signaling molecules called
autoinducers in a process known as quorum sensing. The quorum-sensing network
of the marine bacterium {\it Vibrio harveyi} employs three autoinducers, each
known to encode distinct ecological information. Yet how cells integrate and
interpret the information contained within the three autoinducer signals
remains a mystery. Here, we develop a new framework for analyzing signal
integration based on Information Theory and use it to analyze quorum sensing in
{\it V. harveyi}. We quantify how much the cells can learn about individual
autoinducers and explain the experimentally observed input-output relation of
the {\it V. harveyi} quorum-sensing circuit. Our results suggest that the need
to limit interference between input signals places strong constraints on the
architecture of bacterial signal-integration networks, and that bacteria likely
have evolved active strategies for minimizing this interference. Here we
analyze two such strategies: manipulation of autoinducer production and
feedback on receptor number ratios.


  Calcium activation of striated muscle is known to exhibit a strongly
cooperative dependency on calcium. Because the calcium receptor protein,
troponin (Tn) is known to bind calcium non-cooperatively and has yet to be
linked to a cooperative change in the myosin-blocking protein, tropomyosin
(Tm), we describe a model in which cooperativity is exclusively a
myosin-dependent mechanism. The model couples the energies of three
well-described reactions with actin, namely, actin-Tn, actin-Tm, and
actin-Tm-myosin, to the well-documented positions of Tm, B (blocking), C
(central), and M (myosin-dependent) respectively. Results of simulations with
and without data are consistent with a strand of Tm composed of ~20 subunits
being moved by the concerted action of 3-5 myosin heads resulting in an
all-or-none activation of the entire region of the thin filament overlapped by
myosin. Equations derived from the model fit both equilibrium myosin binding
data and steady-state calcium-dependent tension data and simulate
non-cooperative calcium binding both in the presence and absence of myosin. All
parameters of the model can be determined experimentally. The mechanism is
consistent with steric blocking being both necessary and sufficient for
regulation of striated muscle and can be applied to any actin-based contractile
system that includes Tm and filamentous myosin.


  Rigidity analysis using the "pebble game" has been applied to protein crystal
structures to obtain information on protein folding, assembly and t he
structure-function relationship. However, previous work using this technique
has not made clear how the set of hydrogen-bond constraints included in the
rigidity analysis should be chosen, nor how sensitive the results of rigidity
analysis are to small structural variations. We present a comparative study in
which "pebble game" rigidity analysis is applied to multiple protein crystal
structures, for each of six differen t protein families. We find that the
mainchain rigidity of a protein structure at a given hydrogen-bond energy
cutoff is quite sensitive to small structural variations, and conclude that the
hydrogen bond constraints in rigidity analysis should be chosen so as to form
and test specific hypotheses about the rigidity o f a particular protein. Our
comparative approach highlights two different characteristic patterns ("sudden"
or "gradual") for protein rigidity loss as constraints are re moved, in line
with recent results on the rigidity transitions of glassy networks.


  Tumor Necrosis Factor alpha (TNF) initiates a complex series of biochemical
events in the cell upon binding to its type R1 receptor (TNF-R1). Recent
experimental work has unravelled the molecular regulation of the recruitment of
initial signaling complexes that lead either to cell survival or death.
Survival signals are activated by direct binding of TNF to TNF-R1 at the cell
membrane whereas apoptotic signals by endocytosed TNF/TNF-R1 complexes. Here we
investigate these aspects by developing a quantitative mathematical model of
TNF binding, internalization and intracellular signaling. Model outputs compare
favorably with experimental data and allow to compute TNF-mediated cytotoxicity
as observed in different cell systems. We extensively study the space of
parameters to show that the model is structurally stable and robust over a
broad range of parameter values. Thus, our model is suitable for implementation
in multi-scale simulation programs that are presently under development to
study the behavior of large tumor cell populations.


  Based on the Zimm-Bragg model we study cooperative helix-coil transition
driven by a finite-speed change of temperature. There is an asymmetry between
the coil-to-helix and helix-to-coil transition: the latter is displayed already
for finite speeds, and takes shorter time than the former. This hysteresis
effect has been observed experimentally, and it is explained here via
quantifying system's stability in the vicinity of the critical temperature. A
finite-speed cooling induces a non-equilibrium helical phase with the
correlation length larger than in equilibrium. In this phase the characteristic
length of the coiled domain and the non-equilibrium specific heat can display
an anomalous response to temperature changes. Several pertinent experimental
results on the kinetics helical biopolymers are discussed in detail.


  As an alternative to parsimony analyses, stochastic models have been proposed
(Lewis, 2001), (Nylander, et al., 2004) for morphological characters, so that
maximum likelihood or Bayesian analyses may be used for phylogenetic inference.
A key feature of these models is that they account for ascertainment bias, in
that only varying, or parsimony-informative characters are observed. However,
statistical consistency of such model-based inference requires that the model
parameters be identifiable from the joint distribution they entail, and this
issue has not been addressed.
  Here we prove that parameters for several such models, with finite state
spaces of arbitrary size, are identifiable, provided the tree has at least 8
leaves. If the tree topology is already known, then 7 leaves suffice for
identifiability of the numerical parameters. The method of proof involves first
inferring a full distribution of both parsimony-informative and non-informative
pattern joint probabilities from the parsimony-informative ones, using
phylogenetic invariants. The failure of identifiability of the tree parameter
for 4-taxon trees is also investigated.


  The reaction-diffusion master equation (RDME) is commonly used to model
processes where both the spatial and stochastic nature of chemical reactions
need to be considered. We show that the RDME in many cases is inconsistent with
a microscopic description of diffusion limited chemical reactions and that this
will result in unphysical results. We describe how the inconsistency can be
reconciled if the association and dissociation rates used in the RDME are
derived from the underlying microscopic description. These rate constants will
however necessarily depend on the spatial discretization. At fine spatial
resolution the rates approach the microscopic rate constants defined at the
reaction radius. At low resolution the rates converge to the macroscopic
diffusion limited rate constants in 3D, whereas there is no limiting value in
2D. Our results make it possible to develop spatially discretized
reaction-diffusion models that correspond to a well-defined microscopic
description. We show that this is critical for a correct description of 2D
systems and systems that require high spatial resolution in 3D.


  Motivation: Many biochemical pathways are known, but the numerous parameters
required to correctly explore the dynamics of the pathways are not known. For
this reason, algorithms that can make inferences by looking at the topology of
a network are desirable. In this work, we are particular interested in the
question of whether a given pathway can potentially harbor multiple stable
steady states. In other words, the challenge is to find the set of parameters
such that the dynamical system defined by a set of ordinary differential
equations will contain multiple stable steady states. Being able to find
parameters that cause a network to be bistable may also be benefitial for
engineering synthetic bistable systems where the engineer needs to know a
working set of parameters.
  Result: We have developed an algorithm that optimizes the parameters of a
dynamical system so that the system will contain at least one saddle or
unstable point. The algorithm then looks at trajectories around this saddle or
unstable point to see whether the different trajectories converge to different
stable points. The algorithm returns the parameters that causes the system to
exhibit multiple stable points. Since this is an optimization algorithm, it is
not quaranteed to find a solution. Repeated runs are often required to find a
solution for systems where only a narrow set of parameters exhibit bistability.
  Availability: The C code for the algorithm is available at
http://tinkercell.googlecode.com


  Since proteins carry out biological processes by interacting with other
proteins, analyzing the structure of protein-protein interaction (PPI) networks
could explain complex biological mechanisms, evolution, and disease. Similarly,
studying protein structure networks, residue interaction graphs (RIGs), might
provide insights into protein folding, stability, and function. The first step
towards understanding these networks is finding an adequate network model that
closely replicates their structure. Evaluating the fit of a model to the data
requires comparing the model with real-world networks. Since network
comparisons are computationally infeasible, they rely on heuristics, or
"network properties." We show that it is difficult to assess the reliability of
the fit of a model with any individual network property. Thus, our approach
integrates a variety of network properties and further combines these with a
series of probabilistic methods to predict an appropriate network model for
biological networks. We find geometric random graphs, that model spatial
relationships between objects, to be the best-fitting model for RIGs. This
validates the correctness of our method, since RIGs have previously been shown
to be geometric. We apply our approach to noisy PPI networks and demonstrate
that their structure is also consistent with geometric random graphs.


  The species-area relationship is one of the central generalizations in
ecology however its origin has remained a puzzle. Since ecosystems are
understood as energy transduction systems, the regularities in species richness
are considered to result from ubiquitous imperatives in energy transduction.
  From a thermodynamic point of view, organisms are transduction mechanisms
that distribute an influx of energy down along the steepest gradients to the
ecosystem's diverse repositories of chemical energy, i.e., populations of
species. Transduction machineries, i.e. ecosystems assembled from numerous
species, may emerge and evolve toward high efficiency on large areas that hold
more matter than small ones. This results in the well-known logistic-like
relationship between the area and the number of species. The species-area
relationship is understood, in terms of thermodynamics, to be the skewed
cumulative curve of chemical energy distribution that is commonly known as the
species-abundance relationship.


  Biological macromolecules, proteins and nucleic acids are composed
exclusively of chirally pure monomers. The chirality consensus appears vital
for life and it has even been considered as a prerequisite of life. However the
primary cause for the ubiquitous handedness has remained obscure. We propose
that the chirality consensus is a kinetic consequence that follows from the
principle of increasing entropy, i.e. the 2nd law of thermodynamics. Entropy
increases when an open system evolves by decreasing gradients in free energy
with more and more efficient mechanisms of energy transduction. The rate of
entropy increase is the universal fitness criterion of natural selection that
favors diverse functional molecules and drives the system to the chirality
consensus to attain and maintain high-entropy non-equilibrium states.


  We study the effect of extrinsic noise in metabolic networks. We introduce
external random fluctuations at the kinetic level, and show how these lead to a
stochastic generalization of standard Metabolic Control Analysis. While
Summation and Connectivity Theorems hold true in presence of extrinsic noise,
control coefficients incorporate its effect through an explicit dependency on
the noise intensity. New elasticities and response coefficients are also
defined. Accordingly, the concept of control by noise is introduced as a way of
tuning the systemic behaviour of metabolisms. We argue that this framework
holds for intrinsic noise too, when time-scale separation is present in the
system.


  An important challenge in theoretical ecology is to find good, coarse-grained
representations of complex food webs. Here we use the approach of generalized
modeling to show that it may be possible to formulate a coarse-graining
algorithm that conserves the local dynamics of the model exactly. We show
examples of food webs with a different number of species that have exactly
identical local bifurcation diagrams. Based on these observations, we formulate
a conjecture governing which populations of complex food webs can be grouped
together into a single variable without changing the local dynamics. As an
illustration we use this conjecture to show that chaotic regions generically
exist in the parameter space of a class of food webs with more than three
trophic levels. While our conjecture is at present only applicable to
relatively special cases we believe that its applicability could be greatly
extended if a more sophisticated mapping of parameters were used in the model
reduction.


  A previously established multiscale population genetics model states that
fitness can be inferred from the physical properties of proteins under the
physiological assumption that a loss of stability by any protein confers the
lethal phenotype to an organism. Here we develop this model further by positing
that replication rate (fitness) of a bacterial or viral strain directly depends
on the copy number of folded proteins which determine its replication rate.
Using this model, and both numerical and analytical approaches, we studied the
adaptation process of bacteria and viruses at varied environmental
temperatures. We found that a broad distribution of protein stabilities
observed in the model and in experiment is the key determinant of thermal
response for viruses and bacteria. Our results explain most of the earlier
experimental observations: striking asymmetry of thermal response curves, the
absence of evolutionary trade-off which was expected but not found in
experiments, correlation between denaturation temperature for several protein
families and the Optimal Growth Temperature (OGT) of their host organisms, and
proximity of bacterial or viral OGTs to their evolutionary temperatures. Our
theory quantitatively and with high accuracy described thermal response curves
for 35 bacterial species using, for each species, only two adjustable
parameters, the number of replication rate determining genes and energy barrier
for metabolic reactions.


  The adhesion of cells is mediated by receptors and ligands anchored in
apposing membranes. A central question is how to characterize the binding
affinity of these membrane-anchored molecules. For soluble molecules, the
binding affinity is typically quantified by the binding equilibrium constant
K3D in the linear relation [RL] = K3D [R][L] between the volume concentration
[RL] of bound complexes and the volume concentrations [R] and [L] of unbound
molecules. For membrane-anchored molecules, it is often assumed by analogy that
the area concentration of bound complexes [RL] is proportional to the product
[R][L] of the area concentrations for the unbound receptor and ligand
molecules. We show here (i) that this analogy is only valid for two planar
membranes immobilized on rigid surfaces, and (ii) that the thermal roughness of
flexible membranes leads to cooperative binding of receptors and ligands. In
the case of flexible membranes, the area concentration [RL] of receptor-ligand
bonds is proportional to the square of [R][L] for typical lengths and
concentrations of receptors and ligands in cell adhesion zones. The cooperative
binding helps to understand why different experimental methods for measuring
the binding affinity of membrane-anchored molecules have led to values
differing by several orders of magnitude.


  We present a model study of gliding assays in which actin filaments are moved
by non-processive myosin motors. We show that even if the power stroke of the
motor protein has no lateral component, the filaments will rotate around their
axis while moving over the surface. Notably, the handedness of this twirling
motion is opposite from that of the actin filament structure. It stems from the
fact that the gliding actin filament has "target zones" where its subunits
point towards the surface and are therefore more accessible for myosin heads.
Each myosin head has a higher binding probability before it reaches the center
of the target zone than afterwards, which results in a left-handed twirling. We
present a stochastic simulation and an approximative analytical solution. The
calculated pitch of the twirling motion depends on the filament velocity (ATP
concentration). It reaches about 400nm for low speeds and increases with higher
speeds.


  In this paper, we analyze the temporal evolution of the age-dependent force
of infection and incidence of rubella, after the introduction of a very
specific vaccination programme in a previously nonvaccinated population where
rubella was in endemic steady state. We deduce an integral equation for the
age-dependent force of infection, which depends on a number of parameters that
can be estimated from the force of infection in steady state prior to the
vaccination program. We present the results of our simulations, which are
compared with observed data. We also examine the influence of contact patterns
among members of a community on the age-dependent intensity of transmission of
rubella and on the results of vaccination strategies. As an example of the
theory proposed, we calculate the effects of vaccination strategies for four
communities from Caieiras (Brazil), Huixquilucan (Mexico), Finland and the
United Kingdom. The results for each community differ considerably according to
the distinct intensity and pattern of transmission in the absence of
vaccination. We conclude that this simple vaccination program is not very
efficient (very slow) in the goal of eradicating the disease. This gives
support to a mixed strategy, proposed by Massad et al., accepted and
implemented by the government of the State of Sao Paulo, Brazil.


  The affect of demographic stochasticity of a system of globally coupled
chaotic maps is considered. A two-step model is studied, where the intra-patch
chaotic dynamics is followed by a migration step that coupled all patches; the
equilibrium number of agents on each site, $N$, controls the strength of the
discreteness-induced fluctuations. For small $N$ (large fluctuations) a
period-doubling cascade appears as the coupling (migration) increases. As $N$
grows an extremely slow dynamic emerges, leading to a flow along a
one-dimensional family of almost period 2 solutions. This manifold become a
true solutions in the deterministic limit. The degeneracy between different
attractors that characterizes the clustering phase of the deterministic system
is thus the $N \to \infty$ limit of the slow dynamics manifold.


  The amplified fragment length polymorphism (AFLP) method produces anonymous
genetic markers from throughout a genome. We extend the nucleotide substitution
model of AFLP evolution to additionally include insertion and deletion
processes. The new Sub-ID model relaxes the common assumption that markers are
independent and homologous. We build a Markov chain Monte Carlo methodology
tailored for the Sub-ID model to implement a Bayesian approach to infer AFLP
marker evolution. The method allows us to infer both the phylogenies and the
subset of markers that are possibly homologous. In addition, we can infer the
genome-wide relative rate of indels versus substitutions. In a case study with
AFLP markers from sedges, a grass-like plant common in North America, we find
that accounting for insertion and deletion makes a difference in phylogenetic
inference. The inference of topologies is not sensitive to the prior settings
and the Jukes--Cantor assumption for nucleotide substitution. The model for
insertion and deletion we introduce has potential value in other phylogenetic
applications.


  Galled trees, directed acyclic graphs that model evolutionary histories with
isolated hybridization events, have become very popular due to both their
biological significance and the existence of polynomial time algorithms for
their reconstruction. In this paper we establish to which extent several
distance measures for the comparison of evolutionary networks are metrics for
galled trees, and hence when they can be safely used to evaluate galled tree
reconstruction methods.


  Folding of RNA is subject to a competition between entropy, relevant at high
temperatures, and the random, or random looking, sequence, determining the low-
temperature phase. It is known from numerical simulations that for random as
well as biological sequences, high- and low-temperature phases are different,
e.g. the exponent rho describing the pairing probability between two bases is
rho = 3/2 in the high-temperature phase, and approximatively 4/3 in the
low-temperature (glass) phase. Here, we present, for random sequences, a field
theory of the phase transition separating high- and low-temperature phases. We
establish the existence of the latter by showing that the underlying theory is
renormalizable to all orders in perturbation theory. We test this result via an
explicit 2-loop calculation, which yields rho approximatively 1.36 at the
transition, as well as diverse other critical exponents, including the response
to an applied external force (denaturation transition).


  We study the mechanics of a reversible decohesion (unzipping) of an elastic
layer subjected to quasi-static end-point loading. At the micro level the
system is simulated by an elastic chain of particles interacting with a rigid
foundation through breakable springs. Such system can be viewed as prototypical
for the description of a wide range of phenomena from peeling of polymeric
tapes, to rolling of cells, working of gecko's fibrillar structures and
denaturation of DNA. We construct a rigorous continuum limit of the discrete
model which captures both stable and metastable configurations and present a
detailed parametric study of the interplay between elastic and cohesive
interactions. We show that the model reproduces the experimentally observed
abrupt transition from an incremental evolution of the adhesion front to a
sudden complete decohesion of a macroscopic segment of the adhesion layer. As
the microscopic parameters vary the macroscopic response changes from
quasi-ductile to quasi-brittle, with corresponding decrease in the size of the
adhesion hysteresis. At the micro-scale this corresponds to a transition from a
`localized' to a `diffuse' structure of the decohesion front (domain wall). We
obtain an explicit expression for the critical debonding threshold in the limit
when the internal length scales are much smaller than the size of the system.
The achieved parametric control of the microscopic mechanism can be used in the
design of new biological inspired adhesion devices and machines.


  Cell membranes interact via anchored receptor and ligand molecules. Central
questions on cell adhesion concern the binding affinity of these
membrane-anchored molecules, the mechanisms leading to the receptor-ligand
domains observed during adhesion, and the role of cytoskeletal and other active
processes. In this review, these questions are addressed from a theoretical
perspective. We focus on models in which the membranes are described as elastic
sheets, and the receptors and ligands as anchored molecules. In these models,
the thermal membrane roughness on the nanometer scale leads to a cooperative
binding of anchored receptor and ligand molecules, since the receptor-ligand
binding smoothens out the membranes and facilitates the formation of additional
bonds. Patterns of receptor domains observed in Monte Carlo simulations point
towards a joint role of spontaneous and active processes in cell adhesion. The
interactions mediated by the receptors and ligand molecules can be
characterized by effective membrane adhesion potentials that depend on the
concentrations and binding energies of the molecules.


  We introduce and investigate a series of models for an infection of a
diplodiploid host species by the bacterial endosymbiont \textit{Wolbachia}. The
continuous models are characterized by partial vertical transmission,
cytoplasmic incompatibility and fitness costs associated with the infection. A
particular aspect of interest is competitions between mutually incompatible
strains. We further introduce an age-structured model that takes into account
different fertility and mortality rates at different stages of the life cycle
of the individuals. With only a few parameters, the ordinary differential
equation models exhibit already interesting dynamics and can be used to predict
criteria under which a strain of bacteria is able to invade a population.
Interestingly, but not surprisingly, the age-structured model shows significant
differences concerning the existence and stability of equilibrium solutions
compared to the unstructured model.


  Modelling the process of recombination leads to a large coupled nonlinear
dynamical system. Here, we consider a particular case of recombination in {\em
discrete} time, allowing only for {\em single crossovers}. While the analogous
dynamics in {\em continuous} time admits a closed solution, this no longer
works for discrete time. A more general model (i.e. without the restriction to
single crossovers) has been studied before and was solved algorithmically by
means of Haldane linearisation. Using the special formalism introduced by Baake
and Baake (2003), we obtain further insight into the single-crossover dynamics
and the particular difficulties that arise in discrete time. We then transform
the equations to a solvable system in a two-step procedure: linearisation
followed by diagonalisation. Still, the coefficients of the second step must be
determined in a recursive manner, but once this is done for a given system,
they allow for an explicit solution valid for all times.


  We model the self-organization of the MinE ring that is observed during
subcellular oscillations of the proteins MinD and MinE within the rod-shaped
bacterium {\it Escherichia coli}. With a steady-state approximation, we can
study the MinE-ring generically -- apart from the other details of the Min
oscillation. Rebinding of MinE to depolymerizing MinD filament tips controls
MinE ring formation through a scaled cell shape parameter $\tilde{r}$. We find
two types of E-ring profiles near the filament tip: a strong plateau-like
E-ring controlled by 1D diffusion of MinE along the bacterial length, or a weak
cusp-like E-ring controlled by 3D diffusion near the filament tip. While the
width of a strong E-ring depends on $\tilde{r}$, the occupation fraction of
MinE at the MinD filament tip is saturated and hence the depolymerization speed
do not depend strongly on $\tilde{r}$. Conversely, for weak E-rings both
$\tilde{r}$ and the MinE to MinD stoichiometry strongly control the tip
occupation and hence the depolymerization speed. MinE rings {\em in vivo} are
close to the threshold between weak and strong, and so MinD-filament
depolymerization speed should be sensitive to cell shape, stoichiometry, and
the MinE-rebinding rate. We also find that the transient to MinE-ring formation
is quite long in the appropriate open geometry for assays of ATPase activity
{\it in vitro}, explaining the long delays of ATPase activity observed for
smaller MinE concentrations in those assays without the need to invoke
cooperative MinE activity.


  Receptors of bacterial chemotaxis form clusters at the cell poles, where
clusters act as "antennas" to amplify small changes in ligand concentration.
Interestingly, chemoreceptors cluster at multiple length scales. At the
smallest scale, receptors form dimers, which assemble into stable timers of
dimers. At a large scale, trimers form large polar clusters composed of
thousands of receptors. Although much is known about the signaling properties
emerging from receptor clusters, it is unknown how receptors localize at the
cell poles and what the cluster-size determining factors are. Here, we present
a model of polar receptor clustering based on coupled trimers of dimers, where
cluster size is determined as a minimum of the cluster-membrane free energy.
This energy has contributions from the cluster-membrane elastic energy,
penalizing large clusters due to their high intrinsic curvature, and
receptor-receptor coupling favoring large clusters. We find that the reduced
cluster-membrane curvature mismatch at the curved cell poles leads to large and
robust polar clusters in line with experimental observation, while lateral
clusters are efficiently suppressed.


  We consider a class of physiologically structured population models, a first
order nonlinear partial differential equation equipped with a nonlocal boundary
condition, with a constant external inflow of individuals. We prove that the
linearised system is governed by a quasicontraction semigroup. We also
establish that linear stability of equilibrium solutions is governed by a
generalized net reproduction function. In a special case of the model
ingredients we discuss the nonlinear dynamics of the system when the spectral
bound of the linearised operator equals zero, i.e. when linearisation does not
decide stability. This allows us to demonstrate, through a concrete example,
how immigration might be beneficial to the population. In particular, we show
that from a nonlinearly unstable positive equilibrium a linearly stable and
unstable pair of equilibria bifurcates. In fact, the linearised system exhibits
bistability, for a certain range of values of the external inflow, induced
potentially by All\'{e}e-effect.


  Ras GTPases are lipid-anchored G proteins which play a fundamental role in
cell signaling processes. Electron micrographs of immunogold-labeled Ras have
shown that membrane-bound Ras molecules segregate into nanocluster domains.
Several models have been developed in attempts to obtain quantitative
descriptions of nanocluster formation, but all have relied on assumptions such
as a constant, expression-level independent ratio of Ras in clusters to Ras
monomers (cluster/monomer ratio). However, this assumption is inconsistent with
the law of mass action. Here, we present a biophysical model of Ras clustering
based on short-range attraction and long-range repulsion between Ras molecules
in the membrane. To test this model, we performed Monte Carlo simulations and
compared statistical clustering properties with experimental data. We find that
we can recover the experimentally-observed clustering across a range of Ras
expression levels, without assuming a constant cluster/monomer ratio or the
existence of lipid rafts. In addition, our model makes predictions about the
signaling properties of Ras nanoclusters in support of the idea that Ras
nanoclusters act as an analog-digital-analog converter for high fidelity
signaling.


  Recent advances in genomics have underscored the surprising ubiquity of DNA
copy number variation (CNV). Fortunately, modern genotyping platforms also
detect CNVs with fairly high reliability. Hidden Markov models and algorithms
have played a dominant role in the interpretation of CNV data. Here we explore
CNV reconstruction via estimation with a fused-lasso penalty as suggested by
Tibshirani and Wang [Biostatistics 9 (2008) 18--29]. We mount a fresh attack on
this difficult optimization problem by the following: (a) changing the penalty
terms slightly by substituting a smooth approximation to the absolute value
function, (b) designing and implementing a new MM (majorization--minimization)
algorithm, and (c) applying a fast version of Newton's method to jointly update
all model parameters. Together these changes enable us to minimize the
fused-lasso criterion in a highly effective way. We also reframe the
reconstruction problem in terms of imputation via discrete optimization. This
approach is easier and more accurate than parameter estimation because it
relies on the fact that only a handful of possible copy number states exist at
each SNP. The dynamic programming framework has the added bonus of exploiting
information that the current fused-lasso approach ignores. The accuracy of our
imputations is comparable to that of hidden Markov models at a substantially
lower computational cost.


  Despite decades of extensive research, the large-scale analysis of membrane
proteins remains a difficult task. This is due to the fact that membrane
proteins require a carefully balanced hydrophilic and lipophilic environment,
which optimum varies with different proteins, while most protein chemistry
methods work mainly, if not only, in water-based media. Taking this review
[Santoni, Molloy and Rabilloud, Membrane proteins and proteomics: un amour
impossible? Electrophoresis 2000, 21, 1054-1070] as a pivotal paper, the
current paper analyzes how the field of membrane proteomics exacerbated the
trend in proteomics, i.e. developing alternate methods to the historical
two-dimensional electrophoresis, and thus putting more and more pressure on the
mass spectrometry side. However, in the case of membrane proteins, the
incentive in doing so is due to the poor solubility of membrane proteins. This
review also shows that in some situations, where this solubility problem is
less acute, two-dimensional electrophoresis remains a method of choice. Last
but not least, this review also critically examines the alternate approaches
that have been used for the proteomic analysis of membrane proteins.


  In this paper I will review twenty years of work on the question: When is
there coexistence in stochastic spatial models? The answer, announced in
Durrett and Levin [Theor. Pop. Biol. 46 (1994) 363--394], and that we explain
in this paper is that this can be determined by examining the mean-field ODE.
There are a number of rigorous results in support of this picture, but we will
state nine challenging and important open problems, most of which date from the
1990's.


  The radical pair mechanism is one of the two main hypotheses to explain the
navigability of animals in weak magnetic fields, enabling e.g. birds to see the
Earth's magnetic field. It also plays an essential role in the field of spin
chemistry. Here, we show how quantum control can be used to either enhance or
reduce the performance of such a chemical compass, providing a new route to
further study the radical pair mechanism and its applications. We study the
role of quantum entanglement in this mechanism, and demonstrate intriguing
connections between radical-pair entanglement and the magnetic field
sensitivity of the compass. Beyond their immediate application to the radical
pair mechanism, these results also demonstrate how state-of-the-art quantum
technologies could potentially be used to probe and control biological
functions.


  The importance of torsion vibration in the transmission of life information
is indicated. The localization of quantum torsion state is proved. Following
these analyses a formalism on the quantum theory of conformation-electron
system is proposed. The conformational-electronic transition is calculated by
non-adiabatic operator method. The protein folding is viewed from
conformational quantum transition and the folding rate is calculated. The
time-scale of microsecond to millisecond for the fundamental folding event
(nucleation, collapse, etc) is deduced. The dependence of transition rate W on
N inertial moments is given. It indicates how W increases with the number N of
torsion angles and decreases with the inertial moment I of atomic group in
cooperative transition. The temperature dependence is also deduced which is
different from chemical reaction in high-temperature region. It is demonstrated
that the conformational dynamics gives deep insights into the folding mechanism
and provides a useful tool for analyzing and explaining experimental facts on
the rate of protein folding.


  Many types of cells are able to accurately sense shallow gradients of
chemicals across their diameters, allowing the cells to move towards or away
from chemical sources. This chemotactic ability relies on the remarkable
capacity of cells to infer gradients from particles randomly arriving at
cell-surface receptors by diffusion. Whereas the physical limits of
concentration sensing by cells have been explored, there is no theory for the
physical limits of gradient sensing. Here, we derive such a theory, using as
models a perfectly absorbing sphere and a perfectly monitoring sphere, which,
respectively, infer gradients from the absorbed surface particle density or the
positions of freely diffusing particles inside a spherical volume. We find that
the perfectly absorbing sphere is superior to the perfectly monitoring sphere,
both for concentration and gradient sensing, since previously observed
particles are never remeasured. The superiority of the absorbing sphere helps
explain the presence at the surfaces of cells of signal degrading enzymes, such
as PDE for cAMP in Dictyostelium discoideum (Dicty) and BAR1 for mating factor
alpha in Saccharomyces cerevisiae (budding yeast). Quantitatively, our theory
compares favorably to recent measurements of Dicty moving up a cAMP gradient,
suggesting these cells operate near the physical limits of gradient detection.


  The path-difference metric is one of the oldest distances for the comparison
of fully resolved phylogenetic trees, but its statistical properties are still
quite unknown. In this paper we compute the mean value of the square of the
path-difference metric between two fully resolved rooted phylogenetic trees
with $n$ leaves, under the uniform distribution. This complements previous work
by Steel and Penny, who computed this mean value for fully resolved unrooted
phylogenetic trees.


  We carry out a theoretical study on the isotropic-nematic phase transition
and phase separation in amyloid fibril solutions. Borrowing the thermodynamic
model employed in the study of cylindrical micelles, we investigate the
variations in the fibril length distribution and phase behavior with respect to
changes in the protein concentration, fibril's rigidity, and binding energy. We
then relate our theoretical findings to the nematic ordering observed in Hen
Lysozyme fibril solution.


  We re-examine past suggestions of a close link between terrestrial climate
change and the Sun's transit of spiral arms in its path through the Milky Way
galaxy. These links produced concrete fits, deriving the unknown spiral pattern
speed from terrestrial climate correlations. We test these fits against new
data on spiral structure based on CO data that does not make simplifying
assumptions about symmetry and circular rotation. If we compare the times of
these transits to changes in the climate of Earth, not only do the claimed
correlations disappear, but also we find that they cannot be resurrected for
any reasonable pattern speed.


  Dendritic cells are known to be activated by a wide range of microbial
products, leading to cytokine production and increased levels of membrane
markers such as major histocompatibility complex class II molecules. Such
activated dendritic cells possess the capacity to activate na\"ive T cells. In
the present study we demonstrated that immature dendritic cells secrete both
the YM1 lectin and lipocalin-2. By testing the ligands of these two proteins,
chitosan and siderophores, respectively, we also demonstrated that chitosan, a
degradation product of various fungal and protozoal cell walls, induces an
activation of dendritic cells at the membrane level, as shown by the
up-regulation of membrane proteins such as class II molecules, CD80 and CD86
via a TLR4-dependent mechanism, but is not able to induce cytokine production.
This led to the production of activated dendritic cells unable to stimulate T
cells. However, costimulation with other microbial products overcame this
partial activation and restored the capacity of these activated dendritic cells
to stimulate T cells. In addition, successive stimulation with chitosan and
then by lipopolysaccharide induced a dose-dependent change in the cytokinic
IL-12/IL-10 balance produced by the dendritic cells.


  Recurrence formulas are presented for studying the accuracy of the Fitch
method for reconstructing the ancestral states in a given phylogenetic tree. As
their applications, we analyze the convergence of the accuracy of
reconstructing the root state in a complete binary tree of $2^n$ as $n$ goes to
infinity and also give a lower bound on the accuracy of reconstructing the root
state in an ultrametric tree.


  The cortical array is a structure consisting of highly aligned microtubules
which plays a crucial role in the characteristic uniaxial expansion of all
growing plant cells. Recent experiments have shown polymerization-driven
collisions between the membrane-bound cortical microtubules, suggesting a
possible mechanism for their alignment. We present both a coarse-grained
theoretical model and stochastic particle-based simulations of this mechanism,
and compare the results from these complementary approaches. Our results
indicate that collisions that induce depolymerization are sufficient to
generate the alignment of microtubules in the cortical array.


  Background: Community-acquired methicillin-resistant Staphylococcus aureus
(CA-MRSA), a novel strain of MRSA, has recently emerged and rapidly spread in
the community. Invasion into the hospital setting with replacement of the
hospital-acquired MRSA (HA-MRSA) has also been documented. Co-colonization with
both CA-MRSA and HA-MRSA would have important clinical implications given
differences in antimicrobial susceptibility profiles and the potential for
exchange of genetic information.
  Methods: A deterministic mathematical model was developed to characterize the
transmission dynamics of HA-MRSA and CA-MRSA in the hospital setting and to
quantify the emergence of co-colonization with both strains.
  Results: The model analysis shows that the state of co-colonization becomes
endemic over time and that there is no competitive exclusion of either strain.
Increasing the length of stay or rate of hospital entry among patients
colonized with CA-MRSA leads to a rapid increase in the co-colonized state.
Compared to MRSA decolonization strategy, improving hand hygiene compliance has
the greatest impact on decreasing the prevalence of HA-MRSA, CA-MRSA and the
co-colonized state.
  Conclusions: The model predicts that with the expanding community reservoir
of CA-MRSA, the majority of hospitalized patients will become colonized with
both CA-MRSA and HA-MRSA.


  Continuous and Boolean models for the Drosophila segment polarity network
have shown that the system is able to maintain the wild-type pattern when
subjected to sustained changes in the interaction parameters and initial
conditions. Embryo development is likely to occur under fluctuating
environmental conditions. We use a well-established Boolean model to explore
the ability of the segment polarity network to resist transient changes. We
identify paths along which alternate unviable states are reached, and hence
critical nodes whose state changes lead the system away from the wild-type
state. We find that the system appears to be more sensitive to changes that
involve activation of normally inactive nodes. Through a simulation of the heat
shock response, we show how a localized perturbation in one parasegment is more
deleterious than a global perturbation affecting all parasegments. We identify
the sequence of events involved in the recovery of the system from a global
transient heat shock condition. Finally we discuss these results in terms of
the robustness of the system response.


  The nonequilibrium effective potential is computed for the Frank model of
spontaneous mirror symmetry breaking (SMSB) in chemistry in which external
noise is introduced to account for random environmental effects. When these
fluctuations exceed a critical magnitude, mirror symmetry is restored. The
competition between ambient noise and the chiral bias due to physical fields
and polarized radiation can be explored with this potential.


  We typically observe large-scale outcomes that arise from the interactions of
many hidden, small-scale processes. Examples include age of disease onset,
rates of amino acid substitutions, and composition of ecological communities.
The macroscopic patterns in each problem often vary around a characteristic
shape that can be generated by neutral processes. A neutral generative model
assumes that each microscopic process follows unbiased stochastic fluctuations:
random connections of network nodes; amino acid substitutions with no effect on
fitness; species that arise or disappear from communities randomly. These
neutral generative models often match common patterns of nature. In this paper,
I present the theoretical background by which we can understand why these
neutral generative models are so successful. I show how the classic patterns
such as Poisson and Gaussian arise. Each classic pattern was often discovered
by a simple neutral generative model. The neutral patterns share a special
characteristic: they describe the patterns of nature that follow from simple
constraints on information. For example, any aggregation of processes that
preserves information only about the mean and variance attracts to the Gaussian
pattern; any aggregation that preserves information only about the mean
attracts to the exponential pattern; any aggregation that preserves information
only about the geometric mean attracts to the power law pattern. I present an
informational framework of the common patterns of nature based on the method of
maximum entropy. This framework shows that each neutral generative model is a
special case that helps to discover a particular set of informational
constraints; those informational constraints define a much wider domain of
non-neutral generative processes that attract to the same neutral pattern.


  Chemotactic cells of eukaryotic organisms are able to accurately sense
shallow chemical concentration gradients using cell-surface receptors. This
sensing ability is remarkable as cells must be able to spatially resolve small
fractional differences in the numbers of particles randomly arriving at
cell-surface receptors by diffusion. An additional challenge and source of
uncertainty is that particles, once bound and released, may rebind the same or
a different receptor, which adds to noise without providing any new information
about the environment. We recently derived the fundamental physical limits of
gradient sensing using a simple spherical-cell model, but not including
explicit particle-receptor kinetics. Here, we use a method based on the
fluctuation-dissipation theorem (FDT) to calculate the accuracy of gradient
sensing by realistic receptors. We derive analytical results for two receptors,
as well as two coaxial rings of receptors, e.g. one at each cell pole. For
realistic receptors, we find that particle rebinding lowers the accuracy of
gradient sensing, in line with our previous results.


  Evolutionary game dynamics describes the spreading of successful strategies
in a population of reproducing individuals. Typically, the microscopic
definition of strategy spreading is stochastic, such that the dynamics becomes
deterministic only in infinitely large populations. Here, we introduce a new
microscopic birth--death process that has a fully deterministic strong
selection limit in well--mixed populations of any size. Additionally, under
weak selection, from this new process the frequency dependent Moran process is
recovered. This makes it a natural extension of the usual evolutionary dynamics
under weak selection. We find simple expressions for the fixation probabilities
and average fixation times of the new process in evolutionary games with two
players and two strategies. For cyclic games with two players and three
strategies, we show that the resulting deterministic dynamics crucially depends
on the initial condition in a non--trivial way.