  Expression quantitative trait loci (eQTL) mapping aims to determine genomic
regions that regulate gene transcription. Expression QTL is used to study the
regulatory structure of normal tissues and to search for genetic factors in
complex diseases such as cancer, diabetes, and cystic fibrosis. A modern eQTL
dataset contains millions of SNPs and thousands of transcripts measured for
hundreds of samples. This makes the analysis computationally complex as it
involves independent testing for association for every transcript-SNP pair. The
heavy computational burden makes eQTL analysis less popular, often forces
analysts to restrict their attention to just a subset of transcripts and SNPs.
As larger genotype and gene expression datasets become available, the demand
for fast tools for eQTL analysis increases. We present a new method for fast
eQTL analysis via linear models, called Matrix eQTL. Matrix eQTL can model and
test for association using both linear regression and ANOVA models. The models
can include covariates to account for such factors as population structure,
gender, and clinical variables. It also supports testing of heteroscedastic
models and models with correlated errors. In our experiment on large datasets
Matrix eQTL was thousands of times faster than the existing popular software
for QTL/eQTL analysis. Matrix eQTL is implemented as both Matlab and R packages
and thus can easily be run on Windows, Mac OS, and Linux systems. The software
is freely available at the following address:
http://www.bios.unc.edu/research/genomic_software/Matrix_eQTL


  We consider a model for the evolution of cooperation in a population where
individuals may have one of a number of different heritable and distinguishable
markers or tags. Individuals interact with each of their neighbours on a square
lattice by either cooperating by donating some benefit at a cost to themselves
or defecting by doing nothing. The decision to cooperate or defect is
contingent on each individual's perception of its interacting partner's tag.
Unlike in other tag-based models individuals do not compare their own tag to
that of their interaction partner. That is, there is no {\em self-matching}.
When perception is perfect the cooperation rate is substantially higher than in
the usual spatial prisoner's dilemma game when the cost of cooperation is high.
The enhancement in cooperation is positively correlated with the number of
different tags. The more diverse a population is the more cooperative it
becomes. When individuals start with an inability to perceive tags the
population evolves to a state where individuals gain at least partial
perception. With some reproduction mechanisms perfect perception evolves, but
with others the ability to perceive tags is imperfect. We find that perception
of tags evolves to lower levels when the cost of cooperation is higher.


  We consider here an extension of a previous work by Neves and Serva, still
unpublished, which estimates the amount of interbreeding between anatomically
modern Africans and Neandertals necessary for explaining the experimental fact
that 1 to 4% of the DNA in non-African living humans is of Neandertal origin.
In that work we considered that Africans and Neandertals had the same fitness
(neutral hypothesis) and Neandertal extinction was thus an event of fortune. In
this work we consider that Africans had larger fitnesses. We show results for
four values for the fitness difference: 1%, 5%, 10% and 20% and compare them
with the corresponding neutral results. Some technical differences with respect
to the neutral case appear. We conclude that even with 1% fitness difference
Neandertals extinction comes up in too small a time, so the neutral model looks
more suitable for explaining the known data on occupation of some caves in
Israel for a very long time, alternately by Africans and Neandertals.


  In this contribution, a design of a synthetic calibration genetic circuit to
characterize the relative strength of different sensing promoters is proposed
and its specifications and performance are analyzed via an effective
mathematical model. Our calibrator device possesses certain novel and useful
features like modularity (and thus the possibility of being used in many
different biological contexts), simplicity, being based on a single cell, high
sensitivity and fast response. To uncover the critical model parameters and the
corresponding parameter domain at which the calibrator performance will be
optimal, a sensitivity analysis of the model parameters was carried out over a
given range of sensing protein concentrations (acting as input). Our analysis
suggests that the half saturation constants for repression, sensing and
difference in binding cooperativity (Hill coefficients) for repression are the
key to the performance of the proposed device. They furthermore are determinant
for the sensing speed of the device, showing that it is possible to produce
detectable differences in the repression protein concentrations and in turn in
the corresponding fluorescence in less than two hours. This analysis paves the
way for the design, experimental construction and validation of a new family of
functional genetic circuits for the purpose of calibrating promoters.


  Phenotypic evolution implies sequential fixations of new genomic sequences.
The speed at which these mutations fixate depends, in part, on the relative
fitness (selection coefficient) of the mutant vs. the ancestor. Using a simple
population dynamics model we show that the relative fitness in dynamical
environments is not equal to the fitness averaged over individual environments.
Instead it includes a term that explicitly depends on the sequence of the
environments. This term is geometric in nature and depends only on the oriented
area enclosed by the trajectory taken by the system in the environment state
space. It is related to the well-studied geometric phases in classical and
quantum physical systems. We discuss possible biological implications of these
observations, focusing on evolution of novel metabolic or stress-resistant
functions.


  Coexpression of genes or, more generally, similarity in the expression
profiles poses an unsurmountable obstacle to inferring the gene regulatory
network (GRN) based solely on data from DNA microarray time series. Clustering
of genes with similar expression profiles allows for a course-grained view of
the GRN and a probabilistic determination of the connectivity among the
clusters. We present a model for the temporal evolution of a gene cluster
network which takes into account interactions of gene products with genes and,
through a non-constant degradation rate, with other gene products. The number
of model parameters is reduced by using polynomial functions to interpolate
temporal data points. In this manner, the task of parameter estimation is
reduced to a system of linear algebraic equations, thus making the computation
time shorter by orders of magnitude. To eliminate irrelevant networks, we test
each GRN for stability with respect to parameter variations, and impose
restrictions on its behavior near the steady state. We apply our model and
methods to DNA microarray time series' data collected on Escherichia coli
during glucose-lactose diauxie and infer the most probable cluster network for
different phases of the experiment.


  The bacterial flagellar motor plays a crucial role in both bacterial
locomotion and chemotaxis. Recent experiments reveal that the switching
dynamics of the motor depends on the motor rotation speed, and thus the motor
torque, non-monotonically. Here we present a unified mathematical model which
models motor torque generation based on experimental torque-speed curves and
torque-dependent switching based on the conformational spread model. The model
successfully reproduces the observed switching rate as a function of the
rotation speed, and provides a generic physical explanation independent of most
details. A stator affects the switching dynamics through two mechanisms:
accelerating the conformation flipping rates of individual rotor switching
units, which favours slower motor speed and thus increasing torque; and
affecting more switching units within unit time, which favours faster speed.
Consequently, the switching rate shows a maximum at intermediate speed. Our
model predicts that a motor switches more often with more stators. The
load-switching relation may serve as a mechanism for sensing the physical
environment, similar to the chemotaxis system for sensing the chemical
environment. It may also coordinate the switch dynamics of motors within a
cell.


  Based on the phosphorelay kinetics operative within BvgAS two component
system we propose a mathematical framework for signal transduction and gene
regulation of phenotypic phases in Bordetella pertussis. The proposed model
identifies a novel mechanism of transcriptional interference between two
promoters present in the bvg locus. To understand the system behavior under
elevated temperature, the developed model has been studied in two different
ways. First, a quasi-steady state analysis has been carried out for the two
component system, comprising of sensor BvgS and response regulator BvgA. The
quasi-steady state analysis reveals temperature induced sharp molecular switch,
leading to amplification in the output of BvgA. Accumulation of a large pool of
BvgA thus results into differential regulation of the downstream genes,
including the gene encoding toxin. Numerical integration of the full network
kinetics is then carried out to explore time dependent behavior of different
system components, that qualitatively capture the essential features of
experimental results performed in vivo. Furthermore, the developed model has
been utilized to study mutants that are impaired in their ability to
phosphorylate the transcription factor, BvgA, of the signaling network.


  In this paper, a co-infection of malaria and visceral helminthiasis with
immune stimulation and impairment is studied. We have assumed that an
individual gets malaria infection during invasion by helminths larvae. In
absence of immune response, our results show that the antigens invade the blood
system if the rate of red cell rupture per invading merozoite is greater than
one. If fewer merozoites are released, the initial state is globally
asymptotically stable. If more merozoites are released, there exists
malaria-only endemic point. However, both antigens coexist if the mean
infection burden is greater than one. In this case, there is a threshold value
for drug action below which no recovery of host is expected. In presence of
immune response, three equilibrium states exist. The initial invasion state,
and secondly, the unstable state when the immune population has been triggered
and the antigens have been eliminated. The third state represents the endemic
state which is stable if the infection-induced rupture rate per erythrocyte
with immune activation is less than the total mortality of the erythrocytes. A
model for severity of the co-infection shows that an immune response will be
delayed until immunological barrier values for malaria and helminths are
exceeded.


  We have studied the structure of salt-free lysozyme at 293 K and pH 7.8 using
molecular simulations and experimental SAXS effective potentials between
proteins at three volume fractions, 0.012, 0.033, and 0.12. We found that the
structure of lysozyme near physiological conditions strongly depends on the
volume fraction of proteins. The studied lysozyme solutions are dominated by
monomers only for <0.012; for the strong dilution 70% of proteins are in a form
of monomers. For 0.033 only 20% of proteins do not belong to a cluster. The
clusters are mainly elongated. For 0.12 almost no individual particles exits,
and branched, irregular clusters of large extent appear. Our simulation study
provides new insight into the formation of equilibrium clusters in charged
protein solutions near physiological conditions.


  The process of polymerizing a protein by a ribosome, using a messenger RNA
(mRNA) as the corresponding template, is called {\it translation}. Ribosome may
be regarded as a molecular motor for which the mRNA template serves also as the
track. Often several ribosomes may translate the same (mRNA) simultaneously.
The ribosomes bound simultaneously to a single mRNA transcript are the members
of a polyribosome (or, simply, {\it polysome}). Experimentally measured {\it
polysome profile} gives the distribution of polysome {\it sizes}. Recently a
breakthrough in determining the instantaneous {\it positions} of the ribosomes
on a given mRNA track has been achieved and the technique is called {\it
ribosome profiling} \cite{ingolia10,guo10}. Motivated by the success of these
techniques, we have studied the spatio-temporal organization of ribosomes by
extending a theoretical model that we have reported elsewhere \cite{sharma11}.
This extended version of our model incorporates not only (i) mechano-chemical
cycle of individual ribomes, and (ii) their steric interactions, but also (iii)
the effects of (a) kinetic proofreading, (b) translational infidelity, (c)
ribosome recycling, and (d) sequence inhomogeneities. The theoretical framework
developed here will serve in guiding further experiments and in analyzing the
data to gain deep insight into various kinetic processes involved in
translation.


  An epidemiological model which incorporates synergistic effects that allow
the infectivity and/or susceptibility of hosts to be dependent on the number of
infected neighbours is proposed. Constructive synergy induces an exploitative
behaviour which results in a rapid invasion that infects a large number of
hosts. Interfering synergy leads to a slower and sparser explorative foraging
strategy that traverses larger distances by infecting fewer hosts. The model
can be mapped to a dynamical bond-percolation with spatial correlations that
affect the mechanism of spread but do not influence the critical behaviour of
epidemics.


  Many growing phenomena in both nature and society can be predicted with
sigmoid function. The growth curve of the level of urbanization is a typical
S-shaped one, and can be described by using logistic function. The logistic
model implies a replacement process, and the logistic substitution suggests
non-linear dynamical behaviors such as bifurcation and chaos. Using
mathematical transform and numerical computation, we can demonstrate that the
1-dimensional map comes from a 2-dimensional two-group interaction map. By
analogy with urbanization, a general theory of replacement dynamics is
presented in this paper, and the replacement process can be simulated with the
2-dimansional map. If the rate of replacement is too high, periodic
oscillations and chaos will arise, and the system maybe breaks down. The
replacement theory can be used to interpret various complex interaction and
conversion in physical and human systems. The replacement dynamics provides a
new way of looking at Volterra-Lotka's predator-prey interaction, man-land
relation, and dynastic changes resulting from peasant uprising, and so on.


  It has recently been proposed that proteins embedded in lipidic bio-membranes
can spontaneously self-organize into stable small clusters, or membrane
nano-domains, due to the competition between short-range attractive and
longer-range repulsive forces between proteins, specific to these systems. In
this paper, we carry on our investigation, by Monte Carlo simulations, of
different aspects of cluster phases of proteins in bio-membranes. First, we
compare different long-range potentials (including notably three-body terms) to
demonstrate that the existence of cluster phases should be quite generic.
Furthermore, a real membrane contains hundreds of different protein species
that are far from being randomly distributed in these nano-domains. We take
this protein diversity into account by modulating protein-protein interaction
potentials both at short and longer range. We confirm theoretical predictions
in terms of biological cluster specialization by deciphering how clusters
recruit only a few protein species. In this respect, we highlight that cluster
phases can turn out to be an advantage at the biological level, for example by
enhancing the cell response to external stimuli.


  Chemical reaction networks taken with mass-action kinetics are dynamical
systems that arise in chemical engineering and systems biology. In general,
determining whether a chemical reaction network admits multiple steady states
is difficult, as this requires determining existence of multiple positive
solutions to a large system of polynomials with unknown coefficients. However,
in certain cases, various easy criteria can be applied. One such test is the
Jacobian Criterion, due to Craciun and Feinberg, which gives sufficient
conditions for ruling out the possibility of multiple steady states. A chemical
reaction network is said to pass the Jacobian Criterion if all terms in the
determinant expansion of its parametrized Jacobian matrix have the same sign.
In this article, we present a procedure which simplifies the application of the
Jacobian Criterion, and as a result, we identify a new class of networks for
which multiple steady states is precluded: those in which all chemical species
have total molecularity of at most two. The total molecularity of a species
refers to the sum of all of its stoichiometric coefficients in the network. We
illustrate our results by examining enzyme catalysis networks.


  We analyze how temporal variability in local demography and dispersal combine
to affect the rate of spread of an invading species. Our model combines
state-structured local demography (specified by an integral or matrix
projection model) with general dispersal distributions that may depend on the
state of the individual or its parent, and it allows very general patterns of
stationary temporal variation in both local demography and in the frequency and
distribution of dispersal distances. We show that expressions for the
asymptotic spread rate and its sensitivity to parameters, that have been
derived previously for less general models, continue to hold. Using these
results, we show that random temporal variability in dispersal can accelerate
population spread. Demographic variability can further accelerate spread if it
is positively correlated with dispersal variability, for example if
high-fecundity years are also years in which juveniles tend to settle further
away from their parents. A simple model for the growth and spread of patches of
an invasive plant (perennial pepperweed, Lepidium latifolium) illustrates these
effects and shows that they can have substantial impacts on the predicted speed
of an invasion wave. Temporal variability in dispersal has gotten very little
attention in both the theoretical and empirical literatures on invasive species
spread. Our results suggest that this needs to change.


  A systematic expansion of the many-body master equation for active matter, in
which motors power configurational changes as in the cytoskeleton, is shown to
yield a description of the steady state and responses in terms of an effective
temperature. The effective temperature depends on the susceptibility of the
motors and a Peclet number which measures their strength relative to thermal
Brownian diffusion. The analytic prediction is shown to agree with previous
numerical simulations and experiments. The mapping also establishes a
description of aging in active matter that is also kinetically jammed.


  This study presents a method to predict the growth fluctuation of firms
interdependent in a network economy. The risk of downward growth fluctuation of
firms is calculated from the statistics on Japanese industry.


  We say that a collection $\Cc$ of subsets of $X$ is {\em bureaucratic} if
every maximal hierarchy on $X$ contained in $\Cc$ is also maximum. We
characterise bureaucratic set systems and show how they arise in phylogenetics.
This framework has several useful algorithmic consequences: we generalize some
earlier results and derive a polynomial-time algorithm for a parsimony problem
arising in phylogenetic networks.


  Motivated by recent experiments on bilayer polyhedra composed of amphiphilic
molecules, we study the elastic bending energies of bilayer vesicles forming
polyhedral shapes. Allowing for segregation of excess amphiphiles along the
ridges of polyhedra, we find that bilayer polyhedra can indeed have lower
bending energies than spherical bilayer vesicles. However, our analysis also
implies that, contrary to what has been suggested on the basis of experiments,
the snub dodecahedron, rather than the icosahedron, generally represents the
energetically favorable shape of bilayer polyhedra.


  We investigated the ability of diamond nanoparticles (nanodiamonds, NDs) to
deliver small interfering RNA (siRNA) in Ewing sarcoma cells, in the
perspective of in vivo anti-cancer nucleic acid drug delivery. siRNA was
adsorbed onto NDs previously coated with cationic polymer. Cell uptake of NDs
has been demonstrated by taking advantage of NDs intrinsic fluorescence coming
from embedded color center defects. Cell toxicity of these coated NDs was shown
to be low. Consistent with the internalization efficacy, we have shown a
specific inhibition of EWS/Fli-1 gene expression at the mRNA and protein level
by the ND vectorized siRNA in a serum containing medium.


  Negative and positive transcriptional feedback loops are present in natural
and synthetic genetic oscillators. A single gene with negative transcriptional
feedback needs a time delay and sufficiently strong nonlinearity in the
transmission of the feedback signal in order to produce biochemical rhythms. A
single gene with only positive transcriptional feedback does not produce
oscillations. Here, we demonstrate that this single-gene network in conjunction
with a simple negative interaction can also easily produce rhythms. We examine
a model comprised of two well-differentiated parts. The first is a positive
feedback created by a protein that binds to the promoter of its own gene and
activates the transcription. The second is a negative interaction in which a
repressor molecule prevents this protein from binding to its promoter. A
stochastic study shows that the system is robust to noise. A deterministic
study identifies that the dynamics of the oscillator are mainly driven by two
types of biomolecules: the protein, and the complex formed by the repressor and
this protein. The main conclusion of this paper is that a simple and usual
negative interaction, such as degradation, sequestration or inhibition, acting
on the positive transcriptional feedback of a single gene is a sufficient
condition to produce reliable oscillations. One gene is enough and the positive
transcriptional feedback signal does not need to activate a second repressor
gene. This means that at the genetic level an explicit negative feedback loop
is not necessary. The model needs neither cooperative binding reactions nor the
formation of protein multimers. Therefore, our findings could help to clarify
the design principles of cellular clocks and constitute a new efficient tool
for engineering synthetic genetic oscillators.


  Motivation: Protein interactions are fundamental building blocks of
biochemical reaction systems underlying cellular functions. The complexity and
functionality of such systems emerge not from the protein interactions
themselves but from the dependencies between these interactions. Therefore, a
comprehensive approach for integrating and using information about such
dependencies is required. Results: We present an approach for endowing protein
networks with interaction dependencies using propositional logic, thereby
obtaining protein hypernetworks. First we demonstrate how this framework
straightforwardly improves the prediction of protein complexes. Next we show
that modeling protein perturbations in hypernetworks, rather than in networks,
allows to better infer the functional necessity of proteins for yeast.
Furthermore, hypernetworks improve the prediction of synthetic lethal
interactions in yeast, indicating their capability to capture high-order
functional relations between proteins. Conclusion: Protein hypernetworks are a
consistent formal framework for modeling dependencies between protein
interactions within protein networks. First applications of protein
hypernetworks on the yeast interactome indicate their value for inferring
functional features of complex biochemical systems.


  Canalization is a classic concept in Developmental Biology that is thought to
be an important feature of evolving systems. In a Boolean network it is a form
of network robustness in which a subset of the input signals control the
behavior of a node regardless of the remaining input. It has been shown that
Boolean networks can become canalized if they evolve through a frustrated
competition between nodes. This was demonstrated for large networks in which
each node had K=3 inputs. Those networks evolve to a critical steady-state at
the boarder of two phases of dynamical behavior. Moreover, the evolution of
these networks was shown to be associated with the symmetry of the evolutionary
dynamics. We extend these results to the more highly connected K>3 cases and
show that similar canalized critical steady states emerge with the same
associated dynamical symmetry, but only if the evolutionary dynamics is biased
toward homogeneous Boolean functions.


  de Bruijn graph-based algorithms are one of the two most widely used
approaches for de novo genome assembly. A major limitation of this approach is
the large computational memory space requirement to construct the de Bruijn
graph, which scales with k-mer length and total diversity (N) of unique k-mers
in the genome expressed in base pairs or roughly (2k+8)N bits. This limitation
is particularly important with large-scale genome analysis and for sequencing
centers that simultaneously process multiple genomes. We present a sparse de
Bruijn graph structure, based on which we developed SparseAssembler that
greatly reduces memory space requirements. The structure also allows us to
introduce a novel method for the removal of substitution errors introduced
during sequencing. The sparse de Bruijn graph structure skips g intermediate
k-mers, therefore reducing the theoretical memory space requirement to
~(2k/g+8)N. We have found that a practical value of g=16 consumes approximately
10% of the memory required by standard de Bruijn graph-based algorithms but
yields comparable results. A high error rate could potentially derail the
SparseAssembler. Therefore, we developed a sparse de Bruijn graph-based
denoising algorithm that can remove more than 99% of substitution errors from
datasets with a \leq 2% error rate. Given that substitution error rates for the
current generation of sequencers is lower than 1%, our denoising procedure is
sufficiently effective to safeguard the performance of our algorithm. Finally,
we also introduce a novel Dijkstra-like breadth-first search algorithm for the
sparse de Bruijn graph structure to circumvent residual errors and resolve
polymorphisms.


  This paper examines the objective of optimally harvesting a single species in
a stochastic environment. This problem has previously been analyzed in Alvarez
(2000) using dynamic programming techniques and, due to the natural payoff
structure of the price rate function (the price decreases as the population
increases), no optimal harvesting policy exists. This paper establishes a
relaxed formulation of the harvesting model in such a manner that existence of
an optimal relaxed harvesting policy can not only be proven but also
identified. The analysis embeds the harvesting problem in an
infinite-dimensional linear program over a space of occupation measures in
which the initial position enters as a parameter and then analyzes an auxiliary
problem having fewer constraints. In this manner upper bounds are determined
for the optimal value (with the given initial position); these bounds depend on
the relation of the initial population size to a specific target size. The more
interesting case occurs when the initial population exceeds this target size; a
new argument is required to obtain a sharp upper bound. Though the initial
population size only enters as a parameter, the value is determined in a
closed-form functional expression of this parameter.


  Many recent statistical applications involve inference under complex models,
where it is computationally prohibitive to calculate likelihoods but possible
to simulate data. Approximate Bayesian Computation (ABC) is devoted to these
complex models because it bypasses evaluations of the likelihood function using
comparisons between observed and simulated summary statistics. We introduce the
R abc package that implements several ABC algorithms for performing parameter
estimation and model selection. In particular, the recently developed
non-linear heteroscedastic regression methods for ABC are implemented. The abc
package also includes a cross-validation tool for measuring the accuracy of ABC
estimates, and to calculate the misclassification probabilities when performing
model selection. The main functions are accompanied by appropriate summary and
plotting tools. Considering an example of demographic inference with population
genetics data, we show the potential of the R package.
  R is already widely used in bioinformatics and several fields of biology. The
R abc package will make the ABC algorithms available to the large number of R
users. abc is a freely available R package under the GPL license, and it can be
downloaded at http://cran.r-project.org/web/packages/abc/index.html.


  Cancer results from a sequence of genetic and epigenetic changes which lead
to a variety of abnormal phenotypes including increased proliferation and
survival of somatic cells, and thus, to a selective advantage of pre-cancerous
cells. The notion of cancer progression as an evolutionary process has been
experiencing increasing interest in recent years. Many efforts have been made
to better understand and predict the progression to cancer using mathematical
models; these mostly consider the evolution of a well-mixed cell population,
even though pre-cancerous cells often evolve in highly structured epithelial
tissues. We propose a novel model of cancer progression that considers a
spatially structured cell population where clones expand via adaptive waves.
This model is used to asses two different paradigms of asexual evolution that
have been suggested to delineate the process of cancer progression. The
standard scenario of periodic selection assumes that driver mutations are
accumulated strictly sequentially over time. However, when the mutation supply
is sufficiently high, clones may arise simultaneously on distinct genetic
backgrounds, and clonal adaptation waves interfere with each other. We find
that in the presence of clonal interference, spatial structure increases the
waiting time for cancer, leads to a patchwork structure of non-uniformly sized
clones, decreases the survival probability of virtually neutral (passenger)
mutations, and that genetic distance begins to increase over a characteristic
length scale, determined here. These characteristic features of clonal
interference may help to predict the onset of cancers with pronounced spatial
structure and to interpret spatially-sampled genetic data obtained from
biopsies. Our estimates suggest that clonal interference likely occurs in the
progressing colon cancer, and possibly other cancers where spatial structure
matters.


  The enteric bacterium Proteus mirabilis, which is a pathogen that forms
biofilms in vivo, can swarm over hard surfaces and form concentric ring
patterns in colonies. Colony formation involves two distinct cell types:
swarmer cells that dominate near the surface and the leading edge, and swimmer
cells that prefer a less viscous medium, but the mechanisms underlying pattern
formation are not understood. New experimental investigations reported here
show that swimmer cells in the center of the colony stream inward toward the
inoculation site and in the process form many complex patterns, including
radial and spiral streams, in addition to concentric rings. These new
observations suggest that swimmers are motile and that indirect interactions
between them are essential in the pattern formation. To explain these
observations we develop a hybrid cell-based model that incorporates a
chemotactic response of swimmers to a chemical they produce. The model predicts
that formation of radial streams can be explained as the modulation of the
local attractant concentration by the cells, and that the chirality of the
spiral streams can be predicted by incorporating a swimming bias of the cells
near the surface of the substrate. The spatial patterns generated from the
model are in qualitative agreement with the experimental observations.


  Recently, ultra high-throughput sequencing of RNA (RNA-Seq) has been
developed as an approach for analysis of gene expression. By obtaining tens or
even hundreds of millions of reads of transcribed sequences, an RNA-Seq
experiment can offer a comprehensive survey of the population of genes
(transcripts) in any sample of interest. This paper introduces a statistical
model for estimating isoform abundance from RNA-Seq data and is flexible enough
to accommodate both single end and paired end RNA-Seq data and sampling bias
along the length of the transcript. Based on the derivation of minimal
sufficient statistics for the model, a computationally feasible implementation
of the maximum likelihood estimator of the model is provided. Further, it is
shown that using paired end RNA-Seq provides more accurate isoform abundance
estimates than single end sequencing at fixed sequencing depth. Simulation
studies are also given.


  The accidental oil spill in the Gulf of Mexico in 2010 has caused perceptible
damage to marine and freshwater ecosystems. The large quantity of oil leaking
at a constant rate and the long duration of the event caused an exponentially
increasing mortality of vertebrates. Using data provided by NOAA and USFWS, we
assessed the effects of this event on birds, sea turtles, and mammals.
Mortality rates (measured as the number of carcasses recorded per day) were
exponential for all three groups. Birds were the most affected group, as
indicated by the steepest increase of mortality rates over time. For sea
turtles and mammals, an exponential increase in mortality was observed after an
initial delay. These exponential behaviors are consistent with a unified
scenario for the mortality rate for tetrapod vertebrates. However, at least for
mammals, pre-spill data seem to indicate that the growth in the mortality rate
is not entirely a consequence of the spill.


  Many species engage in acts that could be called creative. However, human
creativity is unique in that it has transformed our planet. Given that the
anatomy of the human brain is not so different from that of the great apes,
what enables us to be so creative? Recent collaborations at the frontier of
anthropology, archaeology, psychology, and cognitive science are culminating in
speculative but increasingly sophisticated efforts to answer to this question.
Examining the skeletons of our ancestors gives cues as to anatomical
constraints that hindered or made possible various kinds of creative
expression. Relics of the past have much to tell us about the thoughts,
beliefs, and creative abilities of the people who invented and used them. How
the spectacular creativity of humans came about is the first topic addressed in
this chapter. Studies at the intersection of creativity and evolution are not
limited to investigations into the biological evolution of a highly creative
species. Creative ideas themselves might be said to evolve through culture.
Human creativity is distinctive because of the adaptive and open-ended manner
in which change accumulates. Inventions build on previous ones in ways that
enhance their utility or aesthetic appeal, or make them applicable in different
situations. There is no a priori limit to how a creative idea might unfold. It
is this proclivity to take an idea and make it our own, or 'put our own spin on
it', that makes creative ideas evolve. The next section of this chapter
investigates in what sense creative ideas evolve through culture. Finally, we
address what forces supported the evolution of creativity. Does being creative
help us live longer, or attract mates? Perhaps creative projects can sometimes
interfere with survival and reproductive fitness; are there non-biological
factors that compel us to create? This is a third topic addressed in this
chapter.


  Human mobility is a key factor in spatial disease dynamics and related
phenomena. In computational models host mobility is typically modelled by
diffusion in space or on metapolulation networks. Alternatively, an effective
force of infection across distance has been introduced to capture spatial
dispersal implicitly. Both approaches do not account for important aspects of
natural human mobility, diffusion does not capture the high degree of
predictability in natural human mobility patters, e.g. the high percentage of
return movements to individuals' base location, the effective force of
infection approach assumes immediate equilibrium with respect to dispersal.
These conditions are typically not met in natural scenarios. We investigate an
epidemiological model that explicitly captures natural individual mobility
patterns. We systematically investigate generic dynamical features of the model
on regular lattices as well as metapopulation networks and show that generally
the model exhibits significant dynamical differences in comparison to ordinary
diffusion and effective force of infection models. For instance, the natural
human mobility model exhibits a saturation of wave front speeds and a novel
type of invasion threshold that is a function of the return rate in mobility
patterns. In the light of these new findings and with the availability of
precise and pervasive data on human mobility our approach provides a framework
for a more sophisticated modeling of spatial disease dynamics.


  We study competition between two biological species advected by a
compressible velocity field. Individuals are treated as discrete Lagrangian
particles that reproduce or die in a density-dependent fashion. In the absence
of a velocity field and fitness advantage, number fluctuations lead to a
coarsening dynamics typical of the stochastic Fisher equation. We then study
three examples of compressible advecting fields: a shell model of turbulence, a
sinusoidal velocity field and a linear velocity sink. In all cases, advection
leads to a striking drop in the fixation time, as well as a large reduction in
the global carrying capacity. Despite localization on convergence zones, one
species goes extinct much more rapidly than in well-mixed populations. For a
weak harmonic potential, one finds a bimodal distribution of fixation times.
The long-lived states in this case are demixed configurations with a single
boundary, whose location depends on the fitness advantage.


  In this paper we will review various aspects of the biology of prions and
focus on what is currently known about the mammalian PrP prion. Also we briefly
describe the prions of yeast and other fungi. Prions are infectious proteins
behaving like genes, i.e. proteins that not only contain genetic information in
its tertiary structure, i.e. its shape, but are also able to transmit and
replicate in a manner analogous to genes but through very different mechanisms.
The term prion is derived from "proteinaceous infectious particle" and arose
from the Prusiner hypothesis that the infectious agent of certain
neurodegenerative diseases was only in a protein, without the participation of
nucleic acids. Currently there are several known types of prion, in addition to
the originally described, which are pathogens of mammals, yeast and other
fungi. Prion proteins are ubiquitous and not always detrimental to their hosts.
This vision of the prion as a causative agent of disease is changing, finding
more and more evidence that they could have important roles in cells and
contribute to the phenotypic plasticity of organisms through the mechanisms of
evolution.


  The position of a reaction front, propagating into an unstable state,
fluctuates because of the shot noise. What is the probability that the
fluctuating front moves considerably slower than its deterministic counterpart?
Can the noise arrest the front motion for some time, or even make it move in
the wrong direction? We present a WKB theory that assumes many particles in the
front region and answers these questions for the microscopic model A->2A, 2A->A
and random walk.


  We introduce several methods of decomposition for two player normal form
games. Viewing the set of all games as a vector space, we exhibit explicit
orthonormal bases for the subspaces of potential games, zero-sum games, and
their orthogonal complements which we call anti-potential games and
anti-zero-sum games, respectively. Perhaps surprisingly, every anti-potential
game comes either from the Rock-Paper-Scissors type games (in the case of
symmetric games) or from the Matching Pennies type games (in the case of
asymmetric games). Using these decompositions, we prove old (and some new)
cycle criteria for potential and zero-sum games (as orthogonality relations
between subspaces). We illustrate the usefulness of our decomposition by (a)
analyzing the generalized Rock-Paper-Scissors game, (b) completely
characterizing the set of all null-stable games, (c) providing a large class of
strict stable games, (d) relating the game decomposition to the decomposition
of vector fields for the replicator equations, (e) constructing Lyapunov
functions for some replicator dynamics, and (f) constructing Zeeman games
-games with an interior asymptotically stable Nash equilibrium and a pure
strategy ESS.


  Evolutionary game theory is an abstract and simple, but very powerful way to
model evolutionary dynamics. Even complex biological phenomena can sometimes be
abstracted to simple two-player games. But often, the interaction between
several parties determines evolutionary success. Rather than pair-wise
interactions, in this case we must take into account the interactions between
many players, which are inherently more complicated than the usual two-player
games, but can still yield simple results. In this manuscript we derive the
composition of a many-player multiple strategy system in the mutation-selection
equilibrium. This results in a simple expression which can be obtained by
recursions using coalescence theory. This approach can be modified to suit a
variety of contexts, e.g. to find the equilibrium frequencies of a finite
number of alleles in a polymorphism or that of different strategies in a social
dilemma in a cultural context.


  In this paper we study some aspects of search for an immobile target by a
swarm of N non-communicating, randomly moving searchers (numbered by the index
k, k = 1, 2,..., N), which all start their random motion simultaneously at the
same point in space. For each realization of the search process, we record the
unordered set of time moments \{\tau_k\}, where \tau_k is the time of the first
passage of the k-th searcher to the location of the target. Clearly, \tau_k's
are independent, identically distributed random variables with the same
distribution function \Psi(\tau). We evaluate then the distribution P(\omega)
of the random variable \omega \sim \tau_1/bar{\tau}, where bar{\tau} = N^{-1}
\sum_{k=1}^N \tau_k is the ensemble-averaged realization-dependent first
passage time. We show that P(\omega) exhibits quite a non-trivial and sometimes
a counterintuitive behaviour. We demonstrate that in some well-studied cases
e.g., Brownian motion in finite d-dimensional domains) the \textit{mean} first
passage time is not a robust measure of the search efficiency, despite the fact
that \Psi(\tau) has moments of arbitrary order. This implies, in particular,
that even in this simplest case (not saying about complex systems and/or
anomalous diffusion) first passage data extracted from a single particle
tracking should be regarded with an appropriate caution because of the
significant sample-to-sample fluctuations.


  The human genome sequence database contains DNA sequences very like those of
mycoplasma molds. It appears such moulds infect not only molecular Biology
laboratories but were picked up by experimenters from contaminated samples and
inserted into GenBank as if they were human. At least one mouldy EST (Expressed
Sequence Tag) has transferred from public databases to commercial tools
(Affymetrix HG-U133 plus 2.0 microarrays). We report a second example
(DA466599) and suggest there is a need to clean up genomic databases but fear
current tools will be inadequate to catch genes which have jumped the silicon
barrier.


  A mesoscopic model for heterogeneous DNA denaturation is developed in the
framework of the path integral formalism. The base pair stretchings are treated
as one-dimensional, time dependent paths contributing to the partition
function. The size of the paths ensemble, which measures the degree of
cooperativity of the system, is computed versus temperature consistently with
the model potential physical requirements. It is shown that the ensemble size
strongly varies with the molecule backbone stiffness providing a quantitative
relation between stacking and features of the melting transition. The latter is
an overall smooth crossover which begins from the \emph{adenine-thymine} rich
portions of the fragment. The harmonic stacking coupling shifts, along the
$T$-axis, the occurrence of the multistep denaturation but it does not change
the character of the crossover. The methods to compute the fractions of open
base pairs versus temperature are discussed: by averaging the base pair
displacements over the path ensemble we find that such fractions signal the
multisteps of the transition in good agreement with the indications provided by
the specific heat plots.


  McNamara and Dall (2011) identified novel relationships between the abundance
of a species in different environments, the temporal properties of
environmental change, and selection for or against dispersal. Here, the
mathematics underlying these relationships in their two-environment model are
investigated for arbitrary numbers of environments. The effect they described
is quantified as the fitness-abundance covariance. The phase in the life cycle
where the population is censused is crucial for the implications of the
fitness-abundance covariance. These relationships are shown to connect to the
population genetics literature on the Reduction Principle for the evolution of
genetic systems and migration. Conditions that produce selection for increased
unconditional dispersal are found to be new instances of departures from
reduction described by the "Principle of Partial Control" proposed for the
evolution of modifier genes. According to this principle, variation that only
partially controls the processes that transform the transmitted information of
organisms may be selected to increase these processes. Mathematical methods of
Karlin, Friedland, and Elsner, Johnson, and Neumann, are central in
generalizing the analysis. Analysis of the adaptive landscape of the model
shows that the evolution of conditional dispersal is very sensitive to the
spectrum of genetic variation the population is capable of producing, and
suggests that empirical study of particular species will require an evaluation
of its variational properties.


  Within systems biology there is an increasing interest in the stochastic
behavior of genetic and biochemical reaction networks. An appropriate
stochastic description is provided by the chemical master equation, which
represents a continuous time Markov chain (CTMC). In this paper we consider the
stochastic properties of a biochemical circuit, known to control eukaryotic
cell cycle and possibly involved in oncogenesis, recently proposed in the
literature within a deterministic framework. Due to the inherent stochasticity
of biochemical processes and the small number of molecules involved, the
stochastic approach should be more correct in describing the real system: we
study the agreement between the two approaches by exploring the system
parameter space. We address the problem by proposing a simplified version of
the model that allows analytical treatment, and by performing numerical
simulations for the full model. We observed optimal agreement between the
stochastic and the deterministic description of the circuit in a large range of
parameters, but some substantial differences arise in at least two cases: 1)
when the deterministic system is in the proximity of a transition from a
monostable to a bistable configuration, and 2) when bistability (in the
deterministic system) is "masked" in the stochastic system by the distribution
tails. The approach provides interesting estimates of the optimal number of
molecules involved in the toggle. Our discussion of the points of strengths,
potentiality and weakness of the chemical master equation in systems biology
and the differences with respect to deterministic modeling are leveraged in
order to provide useful advice for both the bioinformatician practitioner and
the theoretical scientist.


  Dual phospho/dephosphorylation cycles, as well as covalent
enzymatic-catalyzed modifications of substrates, are widely diffused within
cellular systems and are crucial for the control of complex responses such as
learning, memory and cellular fate determination. Despite the large body of
deterministic studies and the increasing work aimed to elucidate the effect of
noise in such systems, some aspects remain unclear. Here we study the
stationary distribution provided by the two-dimensional Chemical Master
Equation for a well known model of a two step phospho/dephosphorylation cycle
using the quasi steady state approximation of the enzymatic kinetics. Our aim
is to analyze the role of fluctuations and the molecules distribution
properties in the transition to a bistable regime. When detailed balance
conditions are satisfied it is possible to compute equilibrium distributions in
a closed and explicit form. When detailed balance is not satisfied, the
stationary non-equilibrium state is strongly influenced by the chemical fluxes.
In the last case, we show how the Helmholtz decomposition of the external field
associated to the generation and recombination transition rates, into a
conservative and a rotational (irreversible) part, allows to compute the
stationary distribution via a perturbative approach . For a finite number of
molecules, there exists diffusion dynamics in a macroscopic region of the state
space, where a relevant transition rate between the two critical points is
observed. Further, the stationary distribution function can be approximated by
the solution of a Fokker-Planck equation. We illustrate the theoretical results
using several numerical simulations.


  The present work investigates environmental effects on the growth of fungal
colonies of P. vitreus by using a lattice-free discrete modelling approach
called FGM (Fuhr et al. (2010), arXiv:1101.1747), in which hyphae and nutrients
are considered as discrete structures. A discrete modelling approach allows
studying the underlying mechanistic rule concerning the basic architecture and
dynamic of fungal networks on the scale of a single colony. By comparing
simulations of the FGM with laboratory experiments of growing fungal colonies
on malt extract agar we show that combined effect of temperature, pH and water
activity on the radial growth rate of a fungal colony on a macroscopic scale
may be explained by a power law for the growth costs of hyphal expansion on a
microscopic scale. The information about the response of the fungal mycelium on
a microscopic scale to environmental conditions is essential to simulate its
behavior in complex structure substrates such as wood, where the impact of the
fungus to the wood (i.e. the degradation of bordered pits or the creation of
bore holes and cavities) changes the local environmental condition, e.g. the
permeability of the substrate and therefore the water activity level in a
colonized wood cell lumen. A combination of diffusion and moisture processes
into wood with the FGM may brighten the knowledge about the colonization
strategy of P. vitreus and helps to optimize its growth behavior for
biotechnological application such as bioincising.


  Biological systems (among others) may respond to a large variety of distinct
external stimuli, or signals. These perturbations will generally be presented
to the system not singly, but in various combinations, so that a proper
understanding of the system response requires assessment of the degree to which
the effects of one signal modulate the effects of another. This paper develops
a pair of structural metrics for sparse differential equation models of complex
dynamic systems and demonstrates that said metrics correlate with proxies of
the susceptibility of one signal-response to be altered in the context of a
second signal. One of these metrics may be interpreted as a normalized arc
density in the neighborhood of certain influential nodes; this metric appears
to correlate with increased independence of signal response.


  Prion diseases are invariably fatal and highly infectious neurodegenerative
diseases affecting humans and animals. By now there have not been some
effective therapeutic approaches to treat all these prion diseases. In 2008,
canine mammals including dogs (canis familials) were the first time
academically reported to be resistant to prion diseases (Vaccine 26: 2601--2614
(2008)). Rabbits are the mammalian species known to be resistant to infection
from prion diseases from other species (Journal of Virology 77: 2003--2009
(2003)). Horses were reported to be resistant to prion diseases too
(Proceedings of the National Academy of Sciences USA 107: 19808--19813 (2010)).
By now all the NMR structures of dog, rabbit and horse prion proteins had been
released into protein data bank respectively in 2005, 2007 and 2010
(Proceedings of the National Academy of Sciences USA 102: 640--645 (2005),
Journal of Biomolecular NMR 38:181 (2007), Journal of Molecular Biology 400:
121--128 (2010)). Thus, at this moment it is very worth studying the NMR
molecular structures of horse, dog and rabbit prion proteins to obtain insights
into their immunity prion diseases.
  This article reports the findings of the molecular structural dynamics of
wild-type horse, dog, and rabbit prion proteins. The dog and horse prion
proteins have stable molecular structures whether under neutral or low pH
environments. Rabbit prion protein has been found having stable molecular
structures under neutral pH environment, but without structural stability under
low pH environment. Under low pH environment, the salt bridges such as
D177--R163 were broken and caused the collapse of the stable $\alpha$-helical
molecular structures.


  The cascade model generates a food web at random. In it the species are
labeled from 0 to $m$, and arcs are given at random between pairs of the
species. For an arc with endpoints $i$ and $j$ ($i<j$), the species $i$ is
eaten by the species labeled $j$. The chain length (height), generated at
random, models the length of food chain in ecological data. The aim of this
note is to introduce the random sequential generation of intervals as a Poisson
model which gives naturally an analogous behavior to the cascade model.


  A framework for the mathematical modeling of evolution in group structured
populations is introduced. The population is divided into a fixed large number
of groups of fixed size. From generation to generation, new groups are formed
that descend from previous groups, through a two-level Fisher-Wright process,
with selection between groups and within groups and with migration between
groups at rate $m$. When $m=1$, the framework reduces to the often used
trait-group framework, so that our setting can be seen as an extension of that
approach. Our framework allows the analysis of previously introduced models in
which altruists and non-altruists compete, and provides new insights into these
models. We focus on the situation in which initially there is a single
altruistic allele in the population, and no further mutations occur. The main
questions are conditions for the viability of that altruistic allele to spread,
and the fashion in which it spreads when it does. Because our results and
methods are rigorous, we see them as shedding light on various controversial
issues in this field, including the role of Hamilton's rule, and of the Price
equation, the relevance of linearity in fitness functions and the need to only
consider pairwise interactions, or weak selection. In this paper we analyze the
early stages of the evolution, during which the number of altruists is small
compared to the size of the population. We show that during this stage the
evolution is well described by a multitype branching process. The driving
matrix for this process can be obtained, reducing the problem of determining
when the altruistic gene is viable to a comparison between the leading
eigenvalue of that matrix, and the fitness of the non-altruists before the
altruistic gene appeared. This leads to a generalization of Hamilton's
condition for the viability of a mutant gene.


  The generating functional method is employed to investigate the synchronous
dynamics of Boolean networks, providing an exact result for the system dynamics
via a set of macroscopic order parameters. The topology of the networks studied
and its constituent Boolean functions represent the system's quenched disorder
and are sampled from a given distribution. The framework accommodates a variety
of topologies and Boolean function distributions and can be used to study both
the noisy and noiseless regimes; it enables one to calculate correlation
functions at different times that are inaccessible via commonly used
approximations. It is also used to determine conditions for the annealed
approximation to be valid, explore phases of the system under different levels
of noise and obtain results for models with strong memory effects, where
existing approximations break down. Links between BN and general Boolean
formulas are identified and common results to both system types are
highlighted.


  Biological molecular machines are proteins that operate under isothermal
conditions hence are referred to as free energy transducers. They can be
formally considered as enzymes that simultaneously catalyze two chemical
reactions: the free energy-donating reaction and the free energy-accepting one.
Most if not all biologically active proteins display a slow stochastic dynamics
of transitions between a variety of conformational substates composing their
native state. In the steady state, this dynamics is characterized by mean
first-passage times between transition substates of the catalyzed reactions. On
taking advantage of the assumption that each reaction proceeds through a single
pair (the gate) of conformational transition substates of the enzyme-substrates
complex, analytical formulas were derived for the flux-force dependence of the
both reactions, the respective stalling forces and the degree of coupling
between the free energy-accepting (output) reaction flux and the free
energy-donating (input) one. The theory is confronted with the results of
random walk simulations on the 5-dimensional hypercube. The formal proof is
given that in the case of reactions proceeding through single gates, the degree
of coupling cannot exceed unity. As some experiments suggest such exceeding,
looking for conditions of increasing the degree of coupling over unity
challenges theory. Though no analytical formulas for models involving more
transition substates are available, study simulations of random walks on
several model networks indicate that the case of the degree of coupling value
higher than one occurs in a natural way for scale-free tree-like networks. This
supports a hypothesis that the protein conformational transition networks, like
higher level biological networks: the proteome and the metabolome, have evolved
in a process of self-organized criticality.


  Evolutionary game theory is employed to study topological conditions of
scale-free networks for the evolution of cooperation. We show that Apollonian
Networks (ANs) are perfect scale-free networks, on which cooperation can spread
to all individuals, even though there are initially only 3 or 4 hubs occupied
by cooperators and all the others by defectors. Local topological features such
as degree, clustering coefficient, gradient as well as topology potential are
adopted to analyze the advantages of ANs in cooperation enhancement.
Furthermore, a degree-skeleton underlying ANs is uncovered for understanding
the cooperation diffusion. Constructing this kind degree-skeleton for random
scale-free networks promotes cooperation level close to that of
Barab\'asi-Albert networks, which gives deeper insights into the origin of the
latter on organization and further promotion of cooperation.


  Drug delivery devices based on nanocomposite membranes containing
thermoresponsive nanogels and superparamagnetic nanoparticles have been
demonstrated to provide reversible, on-off drug release upon application (and
removal) of an oscillating magnetic field. The dose of drug delivered can be
tuned by engineering the phase transition temperature of the nanogel, the
loading of nanogels in the membrane, and the membrane thickness, allowing for
the delivery of drugs over several orders of magnitude of release rates. The
zero-order kinetics of drug release through the membranes permit drug doses
from a specific device to be tuned according to the duration of the magnetic
field. Drugs over a broad range of molecular weights (500-40,000 Da) can be
delivered by the same membrane device. Membrane-to-membrane and cycle-to-cycle
reproducibility is demonstrated, suggesting the general utility of these
membranes for drug delivery.


  The coarsening instability typically disrupts steady-state cluster-size
distributions. We show that degradation coupled to the cluster size, such as
arising from biological proteolysis, leads to a novel fixed-point cluster size.
Stochastic evaporative and condensative fluxes determine the width of the
fixed-point size distribution. At the fixed-point, we show how the peak size
and width depend on number, interactions, and proteolytic rate. This
proteolytic size-control mechanism is consistent with the phenomenology of
pseudo-pilus length control in the general secretion pathway of bacteria.


  We find a statistical mechanism that can adjust orientations of intracellular
filaments to cell geometry in absence of organizing centers. The effect is
based on random and isotropic filament (de-)polymerization dynamics and is
independent of filament interactions and explicit regulation. It can be
understood by an analogy to electrostatics and appears to be induced by the
confining boundaries; for periodic boundary conditions no orientational bias
emerges. Including active transport of particles, the model reproduces
experimental observations of vesicle accumulations in transected axons.


  The Q-cycle mechanism plays an important role in the conversion of the redox
energy into the energy of the proton electrochemical gradient across the
biomembrane. The bifurcated electron transfer reaction, which is built into
this mechanism, recycles one electron, thus, allowing to translocate two
protons per one electron moving to the high-potential redox chain. We study a
kinetic model of the Q-cycle mechanism in an artificial system which mimics the
bf complex of plants and cyanobacteria in the regime of ferredoxin-dependent
cyclic electron flow. Using methods of condensed matter physics, we derive a
set of master equations and describe a time sequence of electron and proton
transfer reactions in the complex. We find energetic conditions when the
bifurcation of the electron pathways at the positive side of the membrane
occurs naturally, without any additional gates. For reasonable parameter
values, we show that this system is able to translocate more than 1.8 protons,
on average, per one electron, with a thermodynamic efficiency of the order of
32% or higher.


  Repeated interaction between individuals is the main mechanism for
maintaining cooperation in social dilemma situations. Variants of tit-for-tat
(repeating the previous action of the opponent) and the win-stay lose-shift
strategy are known as strong competitors in iterated social dilemma games. On
the other hand, real repeated interaction generally allows plasticity (i.e.,
learning) of individuals based on the experience of the past. Although
plasticity is relevant to various biological phenomena, its role in repeated
social dilemma games is relatively unexplored. In particular, if
experience-based learning plays a key role in promotion and maintenance of
cooperation, learners should evolve in the contest with nonlearners under
selection pressure. By modeling players using a simple reinforcement learning
model, we numerically show that learning enables the evolution of cooperation.
We also show that numerically estimated adaptive dynamics appositely predict
the outcome of evolutionary simulations. The analysis of the adaptive dynamics
enables us to capture the obtained results as an affirmative example of the
Baldwin effect, where learning accelerates the evolution to optimality.


  We consider the edge-based compartmental models for epidemic spread developed
in Part I. We show conditions under which simpler models may be substituted for
more detailed models, and in so doing we define a hierarchy of epidemic models.
In particular we provide conditions under which it is appropriate to use the
standard mass action SIR model, and we show what happens when these conditions
fail. Using our hierarchy, we provide a procedure leading to the choice of the
appropriate model for a given population. Our result about the convergence of
models to the Mass Action model gives clear, rigorous conditions under which
the Mass Action model is accurate.


  The primary tool for predicting infectious disease spread and intervention
effectiveness is the mass action Susceptible-Infected-Recovered model of
Kermack and McKendrick. Its usefulness derives largely from its conceptual and
mathematical simplicity; however, it incorrectly assumes all individuals have
the same contact rate and contacts are fleeting. This paper is the first of
three investigating edge-based compartmental modeling, a technique eliminating
these assumptions. In this paper, we derive simple ordinary differential
equation models capturing social heterogeneity (heterogeneous contact rates)
while explicitly considering the impact of contact duration. We introduce a
graphical interpretation allowing for easy derivation and communication of the
model. This paper focuses on the technique and how to apply it in different
contexts. The companion papers investigate choosing the appropriate level of
complexity for a model and how to apply edge-based compartmental modeling to
populations with various sub-structures.


  We consider the edge-based compartmental models for infectious disease spread
introduced in Part I. These models allow us to consider standard SIR diseases
spreading in random populations. In this paper we show how to handle deviations
of the disease or population from the simplistic assumptions of Part I. We
allow the population to have structure due to effects such as demographic
detail or multiple types of risk behavior the disease to have more complicated
natural history. We introduce these modifications in the static network
context, though it is straightforward to incorporate them into dynamic
networks. We also consider serosorting, which requires using the dynamic
network models. The basic methods we use to derive these generalizations are
widely applicable, and so it is straightforward to introduce many other
generalizations not considered here.


  The study of how photosynthetic organisms convert light offers insight not
only into nature's evolutionary process, but may also give clues as to how best
to design and manipulate artificial photosynthetic systems -- and also how far
we can drive natural photosynthetic systems beyond normal operating conditions,
so that they can harvest energy for us under otherwise extreme conditions. In
addition to its interest from a basic scientific perspective, therefore, the
goal to develop a deep quantitative understanding of photosynthesis offers the
potential payoff of enhancing our current arsenal of alternative energy sources
for the future.
  In the following Chapter, we consider the trade-off between dynamics,
structure and function of light harvesting membranes in Rps. Photometricum
purple bacteria, as a model to highlight the priorities that arise when
photosynthetic organisms adapt to deal with the ever-changing natural
environment conditions.


  The inference of networks of dependencies by Gaussian Graphical models on
high-throughput data is an open issue in modern molecular biology. In this
paper we provide a comparative study of three methods to obtain small sample
and high dimension estimates of partial correlation coefficients: the
Moore-Penrose pseudoinverse (PINV), residual correlation (RCM) and
covariance-regularized method $(\ell_{2C})$. We first compare them on simulated
datasets and we find that PINV is less stable in terms of AUC performance when
the number of variables changes. The two regularized methods have comparable
performances but $\ell_{2C}$ is much faster than RCM. Finally, we present the
results of an application of $\ell_{2C}$ for the inference of a gene network
for isoprenoid biosynthesis pathways in Arabidopsis thaliana.


  We discuss general thermodynamic properties of molecular structure formation
processes like protein folding by means of simplified, coarse-grained models.
The conformational transitions accompanying these processes exhibit
similarities to thermodynamic phase transitions, but also significant
differences as the systems that we investigate here are very small. The
usefulness of a microcanonical statistical analysis of these transitions in
comparison with a canonical interpretation is emphasized. The results are
obtained by employing sophisticated generalized-ensemble Markov-chain Monte
Carlo methodologies.


  By means of multicanonical computer simulations, we investigate thermodynamic
properties of the aggregation of interacting semiflexible polymers. We analyze
a mesoscopic bead-stick model, where nonbonded monomers interact via
Lennard-Jones forces. Aggregation turns out to be a process, in which the
constituents experience strong structural fluctuations, similar to peptides in
coupled folding-binding cluster formation processes. In contrast to a recently
studied related proteinlike hydrophobic-polar heteropolymer model, aggregation
and crystallization are separate processes for a homopolymer with the same
small bending rigidity. Rather stiff semiflexible polymers form a
liquid-crystal-like phase, as expected. In analogy to the heteropolymer study,
we find that the first-order-like aggregation transition of the complexes is
accompanied by strong system-size dependent hierarchical surface effects. In
consequence, the polymer aggregation is a phase-separation process with entropy
reduction.


  By means of contact-density chain-growth simulations, we investigate a simple
lattice model of a flexible polymer interacting with an attractive substrate.
The contact density is a function of the numbers of monomer-substrate and
monomer-monomer contacts. These contact numbers represent natural order
parameters and allow for a comprising statistical study of the conformational
space accessible to the polymer in dependence of external parameters such as
the attraction strength of the substrate and the temperature. Since the contact
density is independent of the energy scales associated to the interactions, its
logarithm is an unbiased measure for the entropy of the conformational space.
By setting explicit energy scales, the thus defined, highly general
microcontact entropy can easily be related to the microcanonical entropy of the
corresponding hybrid polymer-substrate system.


  Protein folding cooperativity is defined by the nature of the finite-size
thermodynamic transition exhibited upon folding: two-state transitions show a
free energy barrier between the folded and unfolded ensembles, while downhill
folding is barrierless. A microcanonical analysis, where the energy is the
natural variable, has shown better suited to unambiguously characterize the
nature of the transition compared to its canonical counterpart. Replica
exchange molecular dynamics simulations of a high resolution coarse-grained
model allow for the accurate evaluation of the density of states, in order to
extract precise thermodynamic information, and measure its impact on structural
features. The method is applied to three helical peptides: a short helix shows
sharp features of a two-state folder, while a longer helix and a three-helix
bundle exhibit downhill and two-state transitions, respectively. Extending the
results of lattice simulations and theoretical models, we find that it is the
interplay between secondary structure and the loss of non-native tertiary
contacts which determines the nature of the transition.


  We discuss the hierarchy of subphase transitions in first-order-like
nucleation processes for an exemplified aggregation transition of
heteropolymers. We perform an analysis of the microcanonical entropy, i.e., the
density of states is considered as the central statistical system quantity
since it connects system-specific entropic and energetic information in a
natural and unique way.


  Two-state cooperativity is an important characteristic in protein folding. It
is defined by a depletion of states lying energetically between folded and
unfolded conformations. While there are different ways to test for two-state
cooperativity, most of them probe indirect proxies of this depletion. Yet,
generalized-ensemble computer simulations allow to unambiguously identify this
transition by a microcanonical analysis on the basis of the density of states.
Here we perform a detailed characterization of several helical peptides using
coarse-grained simulations. The level of resolution of the coarse-grained model
allows to study realistic structures ranging from small alpha-helices to a de
novo three-helix bundle - without biasing the force field toward the native
state of the protein. Linking thermodynamic and structural features shows that
while short alpha-helices exhibit two-state cooperativity, the type of
transition changes for longer chain lengths because the chain forms multiple
helix nucleation sites, stabilizing a significant population of intermediate
states. The helix bundle exhibits the signs of two-state cooperativity owing to
favorable helix-helix interactions, as predicted from theoretical models. The
detailed analysis of secondary and tertiary structure formation fits well into
the framework of several folding mechanisms and confirms features observed so
far only in lattice models.


  Interactions among neighboring birds in a flock cause an alignment of their
flight directions. We show that the minimally structured (maximum entropy)
model consistent with these local correlations correctly predicts the
propagation of order throughout entire flocks of starlings, with no free
parameters. These models are mathematically equivalent to the Heisenberg model
of magnetism, and define an "energy" for each configuration of flight
directions in the flock. Comparing flocks of different densities, the range of
interactions that contribute to the energy involves a fixed number of
(topological) neighbors, rather than a fixed (metric) spatial range. Comparing
flocks of different sizes, the model correctly accounts for the observed scale
invariance of long ranged correlations among the fluctuations in flight
direction.


  Deleterious genetic variants can be evaluated as quantitative traits using
information theory-based sequence analysis of recognition sites. To assess the
effect of such variants, fitness and genetic load of SNPs which alter binding
site affinity are derived from changes in individual information and allele
frequencies. Human SNPs that alter mRNA splicing are partitioned according to
their genetic load. SNPs with high genetic loads (>0.5) are common in the
genome and, in many instances, predicted effects are supported by gene
expression studies.


  We extend recent analyses of stochastic effects in game dynamical learning to
cases of multi-player games, and to games defined on networked structures. By
means of an expansion in the noise strength we consider the weak-noise limit,
and present an analytical computation of spectral properties of fluctuations in
multi-player public good games. This extends existing work on two-player games.
In particular we show that coherent cycles may emerge driven by noise in the
adaptation dynamics. These phenomena are not too dissimilar from cyclic
strategy switching observed in experiments of behavioural game theory.


  The last decade saw the advent of increasingly realistic epidemic models that
leverage on the availability of highly detailed census and human mobility data.
Data-driven models aim at a granularity down to the level of households or
single individuals. However, relatively little systematic work has been done to
provide coupled behavior-disease models able to close the feedback loop between
behavioral changes triggered in the population by an individual's perception of
the disease spread and the actual disease spread itself. While models lacking
this coupling can be extremely successful in mild epidemics, they obviously
will be of limited use in situations where social disruption or behavioral
alterations are induced in the population by knowledge of the disease. Here we
propose a characterization of a set of prototypical mechanisms for
self-initiated social distancing induced by local and non-local
prevalence-based information available to individuals in the population. We
characterize the effects of these mechanisms in the framework of a
compartmental scheme that enlarges the basic SIR model by considering separate
behavioral classes within the population. The transition of individuals in/out
of behavioral classes is coupled with the spreading of the disease and provides
a rich phase space with multiple epidemic peaks and tipping points. The class
of models presented here can be used in the case of data-driven computational
approaches to analyze scenarios of social adaptation and behavioral change.


  Community ecologists are principally occupied with the proposition that
natural assemblages of species exhibit orderliness and with identifying its
causes. Plant-pollinator networks exhibit a variety of orderly properties, one
of which is 'nestedness'. Nestedness has been attributed to various causes, but
we propose a further influence arising from the phylogenetic structure of the
biochemical constraints on the pollen diets of bees. We use an artificial
assemblage as an opportunity to isolate the action of this mechanism. The
properties of the network that we studied are consistent with the proposition
that nestedness is caused by the phylogeny of diet range in bees, but the claim
is preliminary and we propose that valuable progress in understanding
plant-pollinator systems may be made through applying the techniques of
chemical ecology at the community scale.


  The design of hybrid peptide-solid interfaces for nanotechnological
applications such as biomolecular nanoarrays requires a deep understanding of
the basic mechanisms of peptide binding and assembly at solid substrates. Here
we show by means of experimental and computational analyses that the adsorption
properties of mutated synthetic peptides at semiconductors exhibit a clear
sequence-dependent adhesion specificity. Our simulations of a novel hybrid
peptide-substrate model reveal the correspondence between proline mutation and
binding affinity to a clean silicon substrate. After synthesizing theoretically
suggested amino-acid sequences with different binding behavior, we confirm the
relevance of the selective mutations upon adhesion in our subsequent atomic
force microscopy experiments.


  ATP-dependent chromatin remodeling enzymes (CRE) are bio-molecular motors in
eukaryotic cells. These are driven by a chemical fuel, namely, adenosine
triphosphate (ATP). CREs actively participate in many cellular processes that
require accessibility of specific segments of DNA which are packaged as
chromatin. The basic unit of chromatin is a nucleosome where 146 bp $\sim$ 50
nm of a double stranded DNA (dsDNA) is wrapped around a spool formed by histone
proteins. The helical path of histone-DNA contact on a nucleosome is also
called "footprint". We investigate the mechanism of footprint traversal by a
CRE that translocates along the dsDNA. Our two-state model of a CRE captures
effectively two distinct chemical (or conformational) states in the
mechano-chemical cycle of each ATP-dependent CRE. We calculate the mean time of
traversal. Our predictions on the ATP-dependence of the mean traversal time can
be tested by carrying out {\it in-vitro} experiments on mono-nucleosomes.


  We revisit the classical problem of nucleated polymerisation and derive a
range of exact results describing polymerisation in systems intermediate
between the well-known limiting cases of a reaction starting from purely
soluble material and for a reaction where no new growth nuclei are formed.


  The hopping model for cargo transport by molecular motors introduced in Refs.
goldman1, goldman2, is extended here in order to incorporate the movement of
cargo-motor complexes. In this context, hopping process expresses the
possibility for cargo to be exchanged between neighbor motors at a microtubule
where the transport takes place. Jamming of motors is essential for cargos to
execute long-range movement in this way. Results from computer simulations
performed using the extended model indicate that cargo may execute
bidirectional movement in the presence of motors of a single polarity,
confirming previous analytical results. Moreover, these results suggest the
existence of a balance between cargo hopping and the movement of the complex
that may control the efficiency of cargo transfer and cargo delivering.
Considerations about the energy involved in the transport process show that the
model presented here offers a considerable advantage over other models in the
literature for which cargo movement is restricted to the movement of
cargo-motor complexes.


  We study a class of processes that are akin to the Wright-Fisher model, with
transition probabilities weighted in terms of the frequency-dependent fitness
of the population types. By considering an approximate weak formulation of the
discrete problem, we are able to derive a corresponding continuous weak
formulation for the probability density. Therefore, we obtain a family of
partial differential equations (PDE) for the evolution of the probability
density, and which will be an approximation of the discrete process in the
joint large population, small time-steps and weak selection limit. If the
fitness functions are sufficiently regular, we can recast the weak formulation
in a more standard formulation, without any boundary conditions, but
supplemented by a number of conservation laws. The equations in this family can
be purely diffusive, purely hyperbolic or of convection-diffusion type, with
frequency dependent convection. The particular outcome will depend on the
assumed scalings. The diffusive equations are of the degenerate type; using a
duality approach, we also obtain a frequency dependent version of the Kimura
equation without any further assumptions. We also show that the convective
approximation is related to the replicator dynamics and provide some estimate
of how accurate is the convective approximation, with respect to the
convective-diffusion approximation. In particular, we show that the mode, but
not the expected value, of the probability distribution is modelled by the
replicator dynamics. Some numerical simulations that illustrate the results are
also presented.


  Metallic nanoparticles are more and more widely used in an increasing number
of applications. Consequently, they are more and more present in the
environment, and the risk that they may represent for human health must be
evaluated. This requires to increase our knowledge of the cellular responses to
nanoparticles. In this context, macrophages appear as an attractive system.
They play a major role in eliminating foreign matter, e.g. pathogens or
infectious agents, by phagocytosis and inflammatory responses, and are thus
highly likely to react to nanoparticles. We have decided to study their
responses to nanoparticles by a combination of classical and wide-scope
approaches such as proteomics. The long term goal of this study is the better
understanding of the responses of macrophages to nanoparticles, and thus to
help to assess their possible impact on human health. We chose as a model
system bone marrow-derived macrophages and studied the effect of commonly used
nanoparticles such as TiO2 and Cu. Classical responses of macrophage were
characterized and proteomic approaches based on 2D gels of whole cell extracts
were used. Preliminary proteomic data resulting from whole cell extracts showed
different effects for TiO2-NPs and Cu-NPs. Modifications of the expression of
several proteins involved in different pathways such as, for example, signal
transduction, endosome-lysosome pathway, Krebs cycle, oxidative stress response
have been underscored. These first results validate our proteomics approach and
open a new wide field of investigation for NPs impact on macrophages


  It has recently been shown that in some DNA microarrays the time needed to
reach thermal equilibrium may largely exceed the typical experimental time,
which is about 15h in standard protocols (Hooyberghs et al. Phys. Rev. E 81,
012901 (2010)). In this paper we discuss how this breakdown of thermodynamic
equilibrium could be detected in microarray experiments without resorting to
real time hybridization data, which are difficult to implement in standard
experimental conditions. The method is based on the analysis of the
distribution of fluorescence intensities I from different spots for probes
carrying base mismatches. In thermal equilibrium and at sufficiently low
concentrations, log I is expected to be linearly related to the hybridization
free energy $\Delta G$ with a slope equal to $1/RT_{exp}$, where $T_{exp}$ is
the experimental temperature and R is the gas constant. The breakdown of
equilibrium results in the deviation from this law. A model for hybridization
kinetics explaining the observed experimental behavior is discussed, the
so-called 3-state model. It predicts that deviations from equilibrium yield a
proportionality of $\log I$ to $\Delta G/RT_{eff}$. Here, $T_{eff}$ is an
effective temperature, higher than the experimental one. This behavior is
indeed observed in some experiments on Agilent arrays. We analyze experimental
data from two other microarray platforms and discuss, on the basis of the
results, the attainment of equilibrium in these cases. Interestingly, the same
3-state model predicts a (dynamical) saturation of the signal at values below
the expected one at equilibrium.


  This paper is devoted to the analysis of the complex damage of DNA irradiated
by ions. The analysis and assessment of complex damage is important because
cells in which it occurs are less likely to survive because the DNA repair
mechanisms may not be sufficiently effective. We studied the flux of secondary
electrons through the surface of nucleosomes and calculated the radial dose and
the distribution of clustered damage around the ion's track. The calculated
radial dose distribution is compared to simulations. The radial distribution of
the complex damage is found to be different from that of the dose. Comparison
with experiments may solve the question of what is more lethal for the cell,
damage complexity or absorbed energy. We suggest a way to calculate the
probability of cell death based on the complexity of the damage. This work is
done within the framework of the phenomenon-based multiscale approach to
radiation damage by ions.


  Diffusion coefficient measurements are important for many biological and
material investigations, such as particle dynamics, kinetics, and size
determinations. Amongst current measurement methods, single particle tracking
(SPT) offers the unique capability of providing location and diffusion
information of a molecule simultaneously while using only femptomoles of
sample. However, the temporal resolution of SPT is limited to seconds for
single-color labeled samples. By directly imaging three dimensional (3D)
diffusing fluorescent proteins and studying the widths of their intensity
profiles, we determine the proteins' diffusion coefficients using single
protein images of sub-millisecond exposure times. This simple method improves
the temporal resolution of diffusion coefficient measurements to
sub-millisecond, and can be readily applied to a range of particle sizes in SPT
investigations and applications where diffusion coefficient measurements are
needed, such as reaction kinetics and particle size determinations.


  This study reveals the essence of ligand recognition mechanisms by which
calmodulin (CaM) controls a variety of Ca2+ signaling processes. We study eight
forms of calcium-loaded CaM each with distinct conformational states. Reducing
the structure to two degrees of freedom conveniently describes main features of
conformational changes of CaM via simultaneous twist-bend motions of the two
lobes. We utilize perturbation-response scanning (PRS) technique, coupled with
molecular dynamics simulations to analyze conformational preferences of
calcium-loaded CaM, initially in extended form. PRS is comprised of sequential
application of directed forces on residues followed by recording the resulting
coordinates. We show that manipulation of a single residue, E31 located in one
of the EF hand motifs, reproduces structural changes to compact forms, and the
flexible linker acts as a transducer of binding information to distant parts of
the protein. Independently, using four different pKa calculation strategies, we
find E31 to be the charged residue (out of 52), whose ionization state is most
sensitive to subtle pH variations in the physiological range. It is proposed
that at relatively low pH, CaM structure is less flexible. By gaining charged
states at specific sites at a pH value around 7, local conformational changes
in the protein will lead to shifts in the energy landscape, paving the way to
other conformational states. These findings are in accordance with FRET
measured shifts in conformational distributions towards more compact forms with
decreased pH. They also corroborate mutational studies and proteolysis results
which point to the significant role of E31 in CaM dynamics.


  Of the major deuterostome groups, the echinoderms with their multiple forms
and complex development are arguably the most mysterious. Although larval
echinoderms are bilaterally symmetric, the adult body seems to abandon the
larval body plan and to develop independently a new structure with different
symmetries. The prevalent pentamer structure, the asymmetry of Loven's rule and
the variable location of the periproct and madrepore present enormous
difficulties in homologizing structures across the major clades, despite the
excellent fossil record. This irregularity in body forms seems to place
echinoderms outside the other deuterostomes. Here I propose that the
predominant five-ray structure is derived from a hexamer structure that is
grounded directly in the structure of the bilaterally symmetric larva. This
hypothesis implies that the adult echinoderm body can be derived directly from
the larval bilateral symmetry and thus firmly ranks even the adult echinoderms
among the bilaterians. In order to test the hypothesis rigorously, a model is
developed in which one ray is missing between rays IV-V (Loven's schema) or
rays C-D (Carpenter's schema). The model is used to make predictions, which are
tested and verified for the process of metamorphosis and for the morphology of
recent and fossil forms. The theory provides fundamental insight into the
M-plane and the Ubisch', Loven's and Carpenter's planes and generalizes them
for all echinoderms. The theory also makes robust predictions about the
evolution of the pentamer structure and its developmental basis. *** including
corrections (see footnotes) ***


  The analysis of non-equilibrium steady states of biochemical reaction
networks relies on finding the configurations of fluxes and chemical potentials
satisfying stoichiometric (mass balance) and thermodynamic (energy balance)
constraints. Efficient methods to explore such states are crucial to predict
reaction directionality, calculate physiologic ranges of variability, estimate
correlations, and reconstruct the overall energy balance of the network from
the underlying molecular processes. While different techniques for sampling the
space generated by mass balance constraints are currently available,
thermodynamics is generically harder to incorporate. Here we introduce a method
to sample the free energy landscape of a reaction network at steady state. In
its most general form, it allows to calculate distributions of fluxes and
concentrations starting from trial functions that may contain prior biochemical
information. We apply our method to the human red blood cell's metabolic
network, whose space of mass-balanced flux states has been sampled extensively
in recent years. Specifically, we profile its thermodynamically feasible flux
configurations, characterizing in detail how fluctuations of fluxes and
potentials are correlated. Based on this, we derive the cell's energy balance
in terms of entropy production, chemical work done and thermodynamic
efficiency.


  Time-resolved single-molecule biophysical experiments yield data that contain
a wealth of dynamic information, in addition to the equilibrium distributions
derived from histograms of the time series. In typical force spectroscopic
setups the molecule is connected via linkers to a read-out device, forming a
mechanically coupled dynamic network. Deconvolution of equilibrium
distributions, filtering out the influence of the linkers, is a straightforward
and common practice. We have developed an analogous dynamic deconvolution
theory for the more challenging task of extracting kinetic properties of
individual components in networks of arbitrary complexity and topology. Our
method determines the intrinsic linear response functions of a given molecule
in the network, describing the power spectrum of conformational fluctuations.
The practicality of our approach is demonstrated for the particular case of a
protein linked via DNA handles to two optically trapped beads at constant
stretching force, which we mimic through Brownian dynamics simulations. Each
well in the protein free energy landscape (corresponding to folded, unfolded,
or possibly intermediate states) will have its own characteristic equilibrium
fluctuations. The associated linear response function is rich in physical
content, since it depends both on the shape of the well and its diffusivity---a
measure of the internal friction arising from such processes like the transient
breaking and reformation of bonds in the protein structure. Starting from the
autocorrelation functions of the equilibrium bead fluctuations measured in this
force clamp setup, we show how an experimentalist can accurately extract the
state-dependent protein diffusivity using a straightforward two-step procedure.


  We provide a logarithmic upper bound for the disentangling number on
unordered lists of leaf labeled trees. This results is useful for analyzing
phylogenetic mixture models. The proof depends on interpreting multisets of
trees as high dimensional contingency tables.


  Genomic DNA is constantly subjected to various mechanical stresses arising
from its biological functions and cell packaging. If the local mechanical
properties of DNA change under torsional and tensional stress, the activity of
DNA-modifying proteins and transcription factors can be affected and regulated
allosterically. To check this possibility, appropriate steady forces and
torques were applied in the course of all-atom molecular dynamics simulations
of DNA with AT- and GC-alternating sequences. It is found that the stretching
rigidity grows with tension as well as twisting. The torsional rigidity is not
affected by stretching, but it varies with twisting very strongly, and
differently for the two sequences. Surprisingly, for AT-alternating DNA it
passes through a minimum with the average twist close to the experimental value
in solution. For this fragment, but not for the GC-alternating sequence, the
bending rigidity noticeably changes with both twisting and stretching. The
results have important biological implications and shed light upon earlier
experimental observations.


  Most biological processes are described as a series of interactions between
proteins and other molecules, and interactions are in turn described in terms
of atomic structures. To annotate protein functions as sets of interaction
states at atomic resolution, and thereby to better understand the relation
between protein interactions and biological functions, we conducted exhaustive
all-against-all atomic structure comparisons of all known binding sites for
ligands including small molecules, proteins and nucleic acids, and identified
recurring elementary motifs. By integrating the elementary motifs associated
with each subunit, we defined composite motifs which represent
context-dependent combinations of elementary motifs. It is demonstrated that
function similarity can be better inferred from composite motif similarity
compared to the similarity of protein sequences or of individual binding sites.
By integrating the composite motifs associated with each protein function, we
define meta-composite motifs each of which is regarded as a time-independent
diagrammatic representation of a biological process. It is shown that
meta-composite motifs provide richer annotations of biological processes than
sequence clusters. The present results serve as a basis for bridging atomic
structures to higher-order biological phenomena by classification and
integration of binding site structures.