  The cell invasion mechanism of Trypanosoma cruzi has similarities with some
intracellular bacterial taxa especially regarding calcium mobilization. This
mechanism is not observed in other trypanosomatids, suggesting that the
molecules involved in this type of cell invasion were a product of (1) acquired
by horizontal gene transfer; (2) secondary loss in the other trypanosomatid
lineages of the mechanism inherited since the bifurcation Bacteria-Neomura (1.9
billion to 900 million years ago) or (3) de novo evolution from non-homologous
proteins via convergent evolution. Similar to T. cruzi, several bacterial
genera require increased host cell cytosolic calcium for intracellular
invasion. Among intracellular bacteria, the mechanism of host cell invasion of
genus Salmonella is the most similar to T. cruzi. The invasion of Salmonella
occurs by contact with the host's cell surface and is mediated by the type III
secretion system (T3SS) that promotes the contact-dependent translocation of
effector proteins directly into host's cell cytoplasm. Here we provide evidence
of distant sequence similarities and structurally conserved domains between T.
cruzi and Salmonella spp T3SS proteins. Exhaustive database searches were
directed to a wide range of intracellular bacteria and trypanosomatids,
exploring sequence patterns for comparison of structural similarities and
Bayesian phylogenies. Based on our data we hypothesize that T. cruzi acquired
genes for calcium mobilization mediated invasion by ancient horizontal gene
transfer from ancestral Salmonella lineages.


  Large biological datasets are being produced at a rapid pace and create
substantial storage challenges, particularly in the domain of high-throughput
sequencing (HTS). Most approaches currently used to store HTS data are either
unable to quickly adapt to the requirements of new sequencing or analysis
methods (because they do not support schema evolution), or fail to provide
state of the art compression of the datasets. We have devised new approaches to
store HTS data that support seamless data schema evolution and compress
datasets substantially better than existing approaches. Building on these new
approaches, we discuss and demonstrate how a multi-tier data organization can
dramatically reduce the storage, computational and network burden of
collecting, analyzing, and archiving large sequencing datasets. For instance,
we show that spliced RNA-Seq alignments can be stored in less than 4% the size
of a BAM file with perfect data fidelity. Compared to the previous compression
state of the art, these methods reduce dataset size more than 20% when storing
gene expression and epigenetic datasets. The approaches have been integrated in
a comprehensive suite of software tools (http://goby.campagnelab.org) that
support common analyses for a range of high-throughput sequencing assays.


  We present GobyWeb, a web-based system to facilitate the management and
analysis of high-throughput sequencing (HTS) projects. The software provides
integrated support for a broad set of HTS analyses and offers a simple plugin
extension mechanism. Analyses currently supported include quantification of
gene expression for messenger and small RNA sequencing, estimation of DNA
methylation (i.e., reduced bisulfite sequencing and whole genome methyl-seq),
or the detection of pathogens in sequenced data. In contrast to many analysis
pipelines developed for analysis of HTS data, GobyWeb requires significantly
less storage space, runs analyses efficiently on a parallel grid, scales
gracefully to process tens or hundreds of multi-gigabyte samples, yet can be
used effectively by researchers who are comfortable using a web browser.
GobyWeb can be obtained at http://gobyweb.campagnelab.org and is freely
available for non-commercial use.


  The reaction-diffusion master equation (RDME) is a lattice stochastic
reaction-diffusion model that has been used to study spatially distributed
cellular processes. The RDME is often interpreted as an approximation to
spatially-continuous models in which molecules move by Brownian motion and
react by one of several mechanisms when sufficiently close. In the limit that
the lattice spacing approaches zero, in two or more dimensions, the RDME has
been shown to lose bimolecular reactions. The RDME is therefore not a
convergent approximation to any spatially-continuous model that incorporates
bimolecular reactions. In this work we derive a new convergent RDME (CRDME) by
finite volume discretization of a spatially-continuous stochastic
reaction-diffusion model popularized by Doi. We demonstrate the numerical
convergence of reaction time statistics associated with the CRDME. For
sufficiently large lattice spacings or slow bimolecular reaction rates, we also
show the reaction time statistics of the CRDME may be approximated by those
from the RDME. The original RDME may therefore be interpreted as an
approximation to the CRDME in several asymptotic limits.


  Logistic functions are good models of biological population growth. They are
also popular in marketing in modelling demand-supply curves and in a different
context, to chart the sales of new products over time. Delays being inherent in
any biological system, we seek to analyse the effect of delays on the growth of
populations governed by the logistic equation. In this paper, the local
stability analysis, rate of convergence and local bifurcation analysis of the
logistic equation with one and two delays is carried out and it can be extended
to a system with multiple delays. Since fluctuating populations are susceptible
to extinction due to sudden and unforeseen environmental disturbances, a
knowledge of the conditions in which the population density is fluctuating or
stable is of great interest in planning and designing control as well as
management strategies.


  A major challenge in single particle reconstruction from cryo-electron
microscopy is to establish a reliable ab-initio three-dimensional model using
two-dimensional projection images with unknown orientations. Common-lines based
methods estimate the orientations without additional geometric information.
However, such methods fail when the detection rate of common-lines is too low
due to the high level of noise in the images. An approximation to the least
squares global self consistency error was obtained using convex relaxation by
semidefinite programming. In this paper we introduce a more robust global self
consistency error and show that the corresponding optimization problem can be
solved via semidefinite relaxation. In order to prevent artificial clustering
of the estimated viewing directions, we further introduce a spectral norm term
that is added as a constraint or as a regularization term to the relaxed
minimization problem. The resulted problems are solved by using either the
alternating direction method of multipliers or an iteratively reweighted least
squares procedure. Numerical experiments with both simulated and real images
demonstrate that the proposed methods significantly reduce the orientation
estimation error when the detection rate of common-lines is low.


  This paper investigates the computational geometry relevant to calculations
of the Frechet mean and variance for probability distributions on the
phylogenetic tree space of Billera, Holmes and Vogtmann, using the theory of
probability measures on spaces of nonpositive curvature developed by Sturm. We
show that the combinatorics of geodesics with a specified fixed endpoint in
tree space are determined by the location of the varying endpoint in a certain
polyhedral subdivision of tree space. The variance function associated to a
finite subset of tree space has a fixed $C^\infty$ algebraic formula within
each cell of the corresponding subdivision, and is continuously differentiable
in the interior of each orthant of tree space. We use this subdivision to
establish two iterative methods for producing sequences that converge to the
Frechet mean: one based on Sturm's Law of Large Numbers, and another based on
descent algorithms for finding optima of smooth functions on convex polyhedra.
We present properties and biological applications of Frechet means and extend
our main results to more general globally nonpositively curved spaces composed
of Euclidean orthants.


  The emergence of structure in cooperative relation is studied in a game
theoretical model. It is proved that specific types of reciprocity norm lead
individuals to split into two groups. The condition for the evolutionary
stability of the norms is also revealed. This result suggests a connection
between group formation and a specific type of reciprocity norm in our society.


  We introduce a method for calculating the extent to which chain non-crossing
is important in the most efficient, optimal trajectories or pathways for a
protein to fold. This involves recording all unphysical crossing events of a
ghost chain, and calculating the minimal uncrossing cost that would have been
required to avoid such events. A depth-first tree search algorithm is applied
to find minimal transformations to fold $\alpha$, $\beta$, $\alpha/\beta$, and
knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a
small fraction of the total distance travelled by a ghost chain. Different
structural classes may be distinguished by the amount of extra uncrossing
distance, and the effectiveness of such discrimination is compared with other
order parameters. It was seen that non-crossing distance over chain length
provided the best discrimination between structural and kinetic classes. The
scaling of non-crossing distance with chain length implies an inevitable
crossover to entanglement-dominated folding mechanisms for sufficiently long
chains. We further quantify the minimal folding pathways by collecting the
sequence of uncrossing moves, which generally involve leg, loop, and elbow-like
uncrossing moves, and rendering the collection of these moves over the unfolded
ensemble as a multiple-transformation "alignment". The consensus minimal
pathway is constructed and shown schematically for representative cases of an
$\alpha$, $\beta$, and knotted protein. An overlap parameter is defined between
pathways; we find that $\alpha$ proteins have minimal overlap indicating
diverse folding pathways, knotted proteins are highly constrained to follow a
dominant pathway, and $\beta$ proteins are somewhere in between. Thus we have
shown how topological chain constraints can induce dominant pathway mechanisms
in protein folding.


  A wide range of applications and research has been done with genome-scale
metabolic models. In this work we describe a methodology for comparing
metabolic networks constructed from genome-scale metabolic models and how to
apply this comparison in order to infer evolutionary distances between
different organisms. Our methodology allows a quantification of the metabolic
differences between different species from a broad range of families and even
kingdoms. This quantification is then applied in order to reconstruct
phylogenetic trees for sets of various organisms.


  Background: The study of genome-scale metabolic models and their underlying
networks is one of the most important fields in systems biology. The complexity
of these models and their description makes the use of computational tools an
essential element in their research. Therefore there is a strong need of
efficient and versatile computational tools for the research in this area.
  Results: In this manuscript we present PyNetMet, a Python library of tools to
work with networks and metabolic models. These are open-source free tools for
use in a Python platform, which adds considerably versatility to them when
compared with their desktop software similars. On the other hand these tools
allow one to work with different standards of metabolic models (OptGene and
SBML) and the fact that they are programmed in Python opens the possibility of
efficient integration with any other already existing Python tool.
  Conclusions: PyNetMet is, therefore, a collection of computational tools that
will facilitate the research work with metabolic models and networks.


  A major problem for inferring species trees from gene trees is that
evolutionary processes can sometimes favour gene tree topologies that conflict
with an underlying species tree. In the case of incomplete lineage sorting,
this phenomenon has recently been well-studied, and some elegant solutions for
species tree reconstruction have been proposed. One particularly simple and
statistically consistent estimator of the species tree under incomplete lineage
sorting is to combine three-taxon analyses, which are phylogenetically robust
to incomplete lineage sorting. In this paper, we consider whether such an
approach will also work under lateral gene transfer (LGT). By providing an
exact analysis of some cases of this model, we show that there is a zone of
inconsistency for triplet-based species tree reconstruction under LGT. However,
a triplet-based approach will consistently reconstruct a species tree under
models of LGT, provided that the expected number of LGT transfers is not too
high. Our analysis involves a novel connection between the LGT problem and
random walks on cyclic graphs. We have implemented a procedure for
reconstructing trees subject to LGT or lineage sorting in settings where taxon
coverage may be patchy and illustrate its use on two sample data sets.


  Transcription factor binding sites vary in their specificity, both within and
between species. Binding specificity has a strong impact on the evolution of
gene expression, because it determines how easily regulatory interactions are
gained and lost. Nevertheless, we have a relatively poor understanding of what
evolutionary forces determine the specificity of binding sites. Here we address
this question by studying regulatory modules composed of multiple binding
sites. Using a population-genetic model, we show that more complex regulatory
modules, composed of a greater number of binding sites, must employ binding
sites that are individually less specific, compared to less complex regulatory
modules. This effect is extremely general, and it hold regardless of the
regulatory logic of a module. We attribute this phenomenon to the inability of
stabilising selection to maintain highly specific sites in large regulatory
modules. Our analysis helps to explain broad empirical trends in the yeast
regulatory network: those genes with a greater number of transcriptional
regulators feature by less specific binding sites, and there is less variance
in their specificity, compared to genes with fewer regulators. Likewise, our
results also help to explain the well-known trend towards lower specificity in
the transcription factor binding sites of higher eukaryotes, which perform
complex regulatory tasks, compared to prokaryotes.


  We investigate the computational structure of a paradigmatic example of
distributed social interaction: that of the open-source Wikipedia community. We
examine the statistical properties of its cooperative behavior, and perform
model selection to determine whether this aspect of the system can be described
by a finite-state process, or whether reference to an effectively unbounded
resource allows for a more parsimonious description. We find strong evidence,
in a majority of the most-edited pages, in favor of a collective-state model,
where the probability of a "revert" action declines as the square root of the
number of non-revert actions seen since the last revert. We provide evidence
that the emergence of this social counter is driven by collective interaction
effects, rather than properties of individual users.


  Tree-size distribution is one of the most investigated subjects in plant
population biology. The forestry literature reports that tree-size distribution
trajectories vary across different stands and/or species, while the metabolic
scaling theory suggests that the tree number scales universally as -2 power of
diameter. Here, we propose a simple functional scaling model in which these two
opposing results are reconciled. Basic principles related to crown shape,
energy optimization and the finite size scaling approach were used to define a
set of relationships based on a single parameter, which allows us to predict
the slope of the tree-size distributions in a steady state condition. We tested
the model predictions on four temperate mountain forests. Plots (4 ha each,
fully mapped) were selected with different degrees of human disturbance
(semi-natural stands vs. formerly managed). Results showed that the size
distribution range successfully fitted by the model is related to the degree of
forest disturbance: in semi-natural forests the range is wide, while in
formerly managed forests, the agreement with the model is confined to a very
restricted range. We argue that simple allometric relationships, at individual
level, shape the structure of the whole forest community.


  Darwinism conceives evolution as a consequence of random variation and
natural selection, hence it is based on a materialistic, i.e. matter-based,
view of science inspired by classical physics. But matter in itself is
considered a very complex notion in modern physics. More specifically, at a
microscopic level, matter and energy are no longer retained within their simple
form, and quantum mechanical models are proposed wherein potential form is
considered in addition to actual form. In this paper we propose an alternative
to standard Neodarwinian evolution theory. We suggest that the starting point
of evolution theory cannot be limited to actual variation whereupon is
selected, but to variation in the potential of entities according to the
context. We therefore develop a formalism, referred to as Context driven
Actualization of Potential (CAP), which handles potentiality and describes the
evolution of entities as an actualization of potential through a reiterated
interaction with the context. As in quantum mechanics, lack of knowledge of the
entity, its context, or the interaction between context and entity leads to
different forms of indeterminism in relation to the state of the entity. This
indeterminism generates a non-Kolmogorovian distribution of probabilities that
is different from the classical distribution of chance described by Darwinian
evolution theory, which stems from a 'actuality focused', i.e. materialistic,
view of nature. We also present a quantum evolution game that highlights the
main differences arising from our new perspective and shows that it is more
fundamental to consider evolution in general, and biological evolution in
specific, as a process of actualization of potential induced by context, for
which its material reduction is only a special case.


  Almost two decades of research on applications of the mathematical formalism
of quantum theory as a modeling tool in domains different from the micro-world
has given rise to many successful applications in situations related to human
behavior and thought, more specifically in cognitive processes of
decision-making and the ways concepts are combined into sentences. In this
article, we extend this approach to animal behavior, showing that an analysis
of an interactive situation involving a mating competition between certain
lizard morphs allows to identify a quantum theoretic structure. More in
particular, we show that when this lizard competition is analyzed structurally
in the light of a compound entity consisting of subentities, the contextuality
provided by the presence of an underlying rock-paper-scissors cyclic dynamics
leads to a violation of Bell's inequality, which means it is of a non-classical
type. We work out an explicit quantum-mechanical representation in Hilbert
space for the lizard situation and show that it faithfully models a set of
experimental data collected on three throat-colored morphs of a specific lizard
species. Furthermore, we investigate the Hilbert space modeling, and show that
the states describing the lizard competitions contain entanglement for each one
of the considered confrontations of lizards with different competing
strategies, which renders it no longer possible to interpret these states of
the competing lizards as compositions of states of the individual lizards.


  Estimates of the fraction of nucleotide substitutions driven by positive
selection vary widely across different species. Accounting for different
estimates of positive selection has been difficult, in part because selection
on polymorphism within a species is known to obscure a signal of positive
selection between species. While methods have been developed to control for the
confounding effects of negative selection against deleterious polymorphism, the
impact of balancing selection on estimates of positive selection has not been
assessed. In Saccharomyces cerevisiae, there is no signal of positive selection
within protein coding sequences as the ratio of nonsynonymous to synonymous
polymorphism is higher than that of divergence. To investigate the impact of
balancing selection on estimates of positive selection we examined five genes
with high rates of nonsynonymous polymorphism in S. cerevisiae relative to
divergence from S. paradoxus. One of the genes, a high affinity zinc
transporter ZRT1, shows an elevated rate of synonymous polymorphism indicative
of balancing selection. The high rate of synonymous polymorphism coincides with
nonsynonymous divergence between three haplotype groups, which we find to be
functionally indistinguishable. We conclude that balancing selection is not
likely to be a common cause of genes harboring a large excess of nonsynonymous
polymorphism in yeast.


  Sequencing errors and biases in metagenomic datasets affect coverage-based
assemblies and are often ignored during analysis. Here, we analyze read
connectivity in metagenomes and identify the presence of problematic and likely
a-biological connectivity within metagenome assembly graphs. Specifically, we
identify highly connected sequences which join a large proportion of reads
within each real metagenome. These sequences show position-specific bias in
shotgun reads, suggestive of sequencing artifacts, and are only minimally
incorporated into contigs by assembly. The removal of these sequences prior to
assembly results in similar assembly content for most metagenomes and enables
the use of graph partitioning to decrease assembly memory and time
requirements.


  A biologically unavoidable sequence is an infinite gender sequence which
occurs in every gendered, infinite genealogical network satisfying certain tame
conditions. We show that every eventually periodic sequence is biologically
unavoidable (this generalizes Koenig's Lemma), and we exhibit some biologically
avoidable sequences. Finally we give an application of unavoidable sequences to
cellular automata.


  Configurational entropy is an important factor in the free energy change of
many macromolecular recognition and binding processes, and has been intensively
studied. Despite great progresses that have been made, the global sampling
remains to be a grand challenge in computational analysis of relevant
processes. Here we propose and demonstrate an entropy estimation method that is
based on physical partition of configurational space and can be readily
combined with currently available methodologies. Tests with two globular
proteins suggest that for flexible macromolecules with large and complex
configurational space, accurate configurational entropy estimation may be
achieved simply by considering the entropically most important subspace. This
conclusion effectively converts an exhaustive sampling problem into a local
sampling one, and defines entropically dominant state for proteins and other
complex macromolecules. The conceptional breakthrough is likely to positively
impact future theoretical analysis, computational algorithm development and
experimental design of diverse chemical and biological molecular systems.


  Invasion waves are a fundamental building block of theoretical ecology. In
this study we aim to take the first steps to link propagation failure and fast
acceleration of traveling waves to critical transitions (or tipping points).
The approach is based upon a detailed numerical study of various versions of
the Fisher-Kolmogorov-Petrovskii-Piscounov (FKPP) equation. The main motivation
of this work is to contribute to the following question: how much information
do statistics, collected by a stationary observer, contain about the speed and
bifurcations of traveling waves? We suggest warning signs based upon closeness
to carrying capacity, second-order moments and transients of localized initial
invasions.


  Viroids are small circular single-stranded infectious RNAs that are
characterized by a relatively high mutation level. Knowledge of their sequence
heterogeneity remains largely elusive, and, as yet, no strategy attempting to
address this question from a population dynamics point of view is in place. In
order to address these important questions, a GF305 indicator peach tree was
infected with a single variant of the Avsunviroidae family member Peach latent
mosaic viroid (PLMVd). Six months post-inoculation, full-length circular
conformers of PLMVd were isolated, deep-sequenced and the resulting sequences
analyzed using an original bioinformatics scheme specifically designed and
developed in order to evaluate the richness of a given the sequence's
population. Two distinct libraries were analyzed, and yielded 1125 and 1061
different PLMVd variants respectively, making this study the most productive to
date (by more than an order of magnitude) in terms of the reporting of novel
viroid sequences. Sequence variants exhibiting up to ~20% of mutations relative
to the inoculated viroid were retrieved, clearly illustrating the high
divergence dynamic inside a unique population. Using a novel hierarchical
clustering algorithm, the different variants obtained were grouped into either
7 or 8 clusters depending on the library being analyzed. Most of the sequences
contained, on average, between 4.6 and 6.3 mutations relative to the variant
used initially to inoculate the plant. Interestingly, it was possible to
reconstitute the sequence evolution between these clusters. On top of providing
a reliable pipeline for the treatment of viroid deep-sequencing, this study
sheds new light on the importance of the sequence variation that may take place
in a viroid population and which may result in the formation of a
quasi-species.


  Genomic data and biomedical imaging data are undergoing exponential growth.
However, our understanding of the phenotype-genotype connection linking the two
types of data is lagging behind. While there are many types of software that
enable the manipulation and analysis of image data and genomic data as separate
entities, there is no framework established for linking the two. We present a
generic set of software tools, BioDIG, that allows linking of image data to
genomic data. BioDIG tools can be applied to a wide range of research problems
that require linking images to genomes. BioDIG features the following: rapid
construction of web-based workbenches, community-based annotation, user
management, and web-services. By using BioDIG to create websites, researchers
and curators can rapidly annotate large number of images with genomic
information. Here we present the BioDIG software tools that include an image
module, a genome module and a user management module. We also introduce a
BioDIG-based website, MyDIG, which is being used to annotate images of
Mycoplasma.


  The United Nations (UN) Population Division is considering producing
probabilistic projections for the total fertility rate (TFR) using the Bayesian
hierarchical model of Alkema et al. (2011), which produces predictive
distributions of TFR for individual countries. The UN is interested in
publishing probabilistic projections for aggregates of countries, such as
regions and trading blocs. This requires joint probabilistic projections of
future country-specific TFRs, taking account of the correlations between them.
We propose an extension of the Bayesian hierarchical model that allows for
probabilistic projection of TFR for any set of countries. We model the
correlation between country forecast errors as a linear function of time
invariant covariates, namely whether the countries are contiguous, whether they
had a common colonizer after 1945, and whether they are in the same UN region.
The resulting correlation model is incorporated into the Bayesian hierarchical
model's error distribution. We produce predictive distributions of TFR for
1990-2010 for each of the UN's primary regions. We find that the proportions of
the observed values that fall within the prediction intervals from our method
are closer to their nominal levels than those produced by the current model.
Our results suggest that a significant proportion of the correlation between
forecast errors for TFR in different countries is due to countries' geographic
proximity to one another, and that if this correlation is accounted for, the
quality of probabilitistic projections of TFR for regions and other aggregates
is improved.


  Darwinian evolution can be illustrated as an uphill walk in a landscape,
where the surface consists of genotypes, the height coordinates represent
fitness, and each step corresponds to a point mutation. Epistasis, roughly
defined as the dependence between the fitness effects of mutations, is a key
concept in the theory of adaptation. Important recent approaches depend on
graphs and polytopes. Fitness graphs are useful for describing coarse
properties of a landscape, such as mutational trajectories and the number of
peaks. The graphs have been used for relating global and local properties of
fitness landscapes. The geometric theory of gene interaction, or the shape
theory, is the most fine-scaled approach to epistasis. Shapes, defined as
triangulations of polytopes for any number of loci, replace the well
established concepts of positive and negative epistasis for two mutations. From
the shape one can identify the fittest populations, i.e., populations where
allele shuffling (recombination) will not increase the mean fitness. Shapes and
graphs provide complementary information. The approaches make no structural
assumptions about the underlying fitness landscapes, which make them well
suited for empirical work.


  Predicting the response of a specific cancer to a therapy is a major goal in
modern oncology that should ultimately lead to a personalised treatment.
High-throughput screenings of potentially active compounds against a panel of
genomically heterogeneous cancer cell lines have unveiled multiple
relationships between genomic alterations and drug responses. Various
computational approaches have been proposed to predict sensitivity based on
genomic features, while others have used the chemical properties of the drugs
to ascertain their effect. In an effort to integrate these complementary
approaches, we developed machine learning models to predict the response of
cancer cell lines to drug treatment, quantified through IC50 values, based on
both the genomic features of the cell lines and the chemical properties of the
considered drugs. Models predicted IC50 values in a 8-fold cross-validation and
an independent blind test with coefficient of determination R2 of 0.72 and 0.64
respectively. Furthermore, models were able to predict with comparable accuracy
(R2 of 0.61) IC50s of cell lines from a tissue not used in the training stage.
Our in silico models can be used to optimise the experimental design of
drug-cell screenings by estimating a large proportion of missing IC50 values
rather than experimentally measure them. The implications of our results go
beyond virtual drug screening design: potentially thousands of drugs could be
probed in silico to systematically test their potential efficacy as anti-tumour
agents based on their structure, thus providing a computational framework to
identify new drug repositioning opportunities as well as ultimately be useful
for personalized medicine by linking the genomic traits of patients to drug
sensitivity.


  Large scale surveys in mammalian tissue culture cells suggest that the
protein expressed at the median abundance is present at 8,000 - 16,000
molecules per cell and that differences in mRNA expression between genes
explain only 10-40% of the differences in protein levels. We find, however,
that these surveys have significantly underestimated protein abundances and the
relative importance of transcription. Using individual measurements for 61
housekeeping proteins to rescale whole proteome data from Schwanhausser et al.,
we find that the median protein detected is expressed at 170,000 molecules per
cell and that our corrected protein abundance estimates show a higher
correlation with mRNA abundances than do the uncorrected protein data. In
addition, we estimated the impact of further errors in mRNA and protein
abundances, showing that mRNA levels explain at least 56% of the differences in
protein abundance for the genes detected by Schwanhausser et al., though
because one major source of error could not be estimated the true percent
contribution could be higher. We also employed a second, independent strategy
to determine the contribution of mRNA levels to protein expression. We show
that the variance in translation rates directly measured by ribosome profiling
is only 12% of that inferred by Schwanhausser et al. and that the measured and
inferred translation rates correlate only poorly (R2=0.13). Based on this, our
second strategy suggests that mRNA levels explain ~81% of the variance in
protein levels. We also determined the percent contributions of transcription,
RNA degradation, translation and protein degradation to the variance in protein
abundances using both of our strategies. While the magnitudes of the two
estimates vary, they both suggest that transcription plays a more important
role than the earlier studies implied and translation a much smaller role.


  Arabidopsis thaliana is an important model organism for understanding the
genetics and molecular biology of plants. Its highly selfing nature, together
with other important features, such as small size, short generation time, small
genome size, and wide geographic distribution, make it an ideal model organism
for understanding natural variation. Genome-wide association studies (GWAS)
have proven a useful technique for identifying genetic loci responsible for
natural variation in A. thaliana. Previously genotyped accessions (natural
inbred lines) can be grown in replicate under different conditions, and
phenotyped for different traits. These important features greatly simplify
association mapping of traits and allow for systematic dissection of the
genetics of natural variation by the entire Arabidopsis community. To
facilitate this, we present GWAPP, an interactive web-based application for
conducting GWAS in A. thaliana. Using an efficient Python implementation of a
linear mixed model, traits measured for a subset of 1386 publicly available
ecotypes can be uploaded and mapped with an efficient mixed model and other
methods in just a couple of minutes. GWAPP features an extensive, interactive,
and a user-friendly interface that includes interactive manhattan plots and
interactive local and genome-wide LD plots. It facilitates exploratory data
analysis by implementing features such as the inclusion of candidate SNPs in
the model as cofactors.


  We have developed a model for the interactions between keratin intermediate
filaments based on self-consistent field theory. The intermediate filaments are
represented by charged surfaces, and the disordered terminal domains of the
keratins are represented by charged heteropolymers grafted to these surfaces.
We estimate the system is close to a charge compensation point where the
heteropolymer grafting density is matched to the surface charge density. Using
a protein model with amino acid resolution for the terminal domains, we find
that the terminal chains can mediate a weak attraction between the keratin
surfaces. The origin of the attraction is a combination of bridging and
electrostatics. The attraction disappears when the system moves away from the
charge compensation point, or when excess small ions and/or NMF-representing
free amino acids are added. These results are in concordance with experimental
observations, and support the idea that the interaction between keratin
filaments, and ultimately in part the elastic properties of the
keratin-containing tissue, is controlled by a combination of the
physico-chemical properties of the disordered terminal domains and the
composition of the medium in the inter-filament region.


  The evolution of parasite virulence and host defences is affected by
population structure. This effect has been confirmed in studies focusing on
large spatial scales, whereas the importance of local structure is not well
understood. Slavemaking ants are social parasites that exploit workers of
another species to rear their offspring. Enslaved workers of the host species
Temnothorax longispinosus have been found to exhibit an effective
post-enslavement defence behaviour: enslaved workers were observed killing a
large proportion of the parasites' offspring. Since enslaved workers do not
reproduce, they gain no direct fitness benefit from this 'rebellion' behaviour.
However, there may be an indirect benefit: neighbouring host nests that are
related to 'rebel' nests can benefit from a reduced raiding pressure, as a
result of the reduction in parasite nest size due to the enslaved workers'
killing behaviour. We use a simple mathematical model to examine whether the
small-scale population structure of the host species could explain the
evolution of this potentially altruistic defence trait against slavemaking
ants. We find that this is the case if enslaved host workers are related to
nearby host nests. In a population genetic study we confirm that enslaved
workers are, indeed, more closely related to host nests within the raiding
range of their resident slavemaker nest, than to host nests outside the raiding
range. This small-scale population structure seems to be a result of polydomy
(e.g. the occupation of several nests in close proximity by a single colony)
and could have enabled the evolution of 'rebellion' by kin selection.


  Gene duplication is believed to play a major role in the evolution of genomic
complexity. The presence of a duplicate removes the constraint of natural
selection upon the gene, leading to its likely loss of function or,
occasionally, the gain of a novel one. Alternately, a pleiotropic gene might
partition its functions among its duplicates, thus preserving both copies.
Duplicate genes is not a novelty for diploid genotypes, but only for
haplotypes. In this paper, we study the consequences of regulatory interactions
in diploid genotypes and explore how the context of allelic interactions gives
rise to dynamical phenotypes that enable duplicate genes to spread in a
population. The regulatory network we study is that of a single autoregulatory
activator gene, and the two copies of the gene diverge either as alleles in a
diploid species or as duplicates in haploids. These differences are in their
transcriptional ability -- either via alterations to its activating domain, or
to its cis-regulatory binding repertoire. When cis-regulatory changes are
introduced that partition multiple regulatory triggers among the duplicates, it
is shown that mutually exclusive expression states of the duplicates that
emerge are accompanied by a back-up facility: when a highly expressed gene is
deleted, the previously unexpressed duplicate copy compensates for it. The
diploid version of the regulatory network model can account for allele-specific
expression variants, and a model of inheritance of the haplotype network
enables us to trace the evolutionary consequence of heterozygous phenotypes.
This is modelled for the variations in the activating domain of one copy,
whereby stable as well as transiently bursting oscillations ensue in single
cells. The evolutionary model shows that these phenotypic states accessible to
a diploid, heterozygous genotype enable the spread of a duplicated haplotype.


  In this work we apply a mathematical model of photosynthesis to quantify the
potential for photosynthetic life in the very Early Archean oceans. We assume
the presence of oceanic blockers of ultraviolet radiation, specifically ferrous
ions. For this scenario, our results suggest a potential for photosynthetic
life greater than or similar to that in later eras/eons, such as the Late
Archean and the current Phanerozoic eon.


  Availability of high-resolution crystal structures of ribosomal subunits of
different species opens a route to investigate about molecular interactions
between its constituents and stabilization strategy. Structural analysis of the
small ribosomal subunit shows that primary binder proteins are mainly employed
in stabilizing the folded ribosomal RNA by their high negative free energy of
association, where tertiary binders mainly help to stabilize protein-protein
interfaces. Secondary binders perform both the functions. Conformational
changes of prokaryotic and eukaryotic ribosomal proteins due to complexation
with 16S ribosomal RNA are linearly correlated with their RNA-interface area
and free energy of association. The proteins having long extensions buried
within ribosomal RNA have more flexible structures than those found on the
subunit surface. Thermus thermophilus ribosomal proteins undergo high
conformational changes compared to those of Escherichia coli, assuring
structural stability at high temperature environment. The general stabilization
strategy of ribosomal protein-RNA interfaces is shown, where high interface
polarity ensures high surface density of Hydrogen bonds even with low
base/backbone ratio. Polarity is regulated in evolutionary strategy of
ribosomal proteins. Thus, the habitat environmental conditions of the two
species sweet up their ribosomal protein-RNA interfaces to alter its physical
parameters in order of stabilization.


  Using Langevin dynamics simulations, we investigate the dynamics of polymer
translocation into a circular nanocontainer through a nanopore under a driving
force $F$. We observe that the translocation probability initially increases
and then saturates with increasing $F$, independent of $\phi$, which is the
average density of the whole chain in the nanocontainer. The translocation time
distribution undergoes a transition from a Gaussian distribution to an
asymmetric distribution with increasing $\phi$. Moreover, we find a
nonuniversal scaling exponent of the translocation time as chain length,
depending on $\phi$ and $F$. These results are interpreted by the conformation
of the translocated chain in the nanocontainer and the time of an individual
segment passing through the pore during translocation.


  We investigate the chain conformation of ring polymers confined to a
cylindrical nanochannel using both theoretical analysis and three dimensional
Langevin dynamics simulations. We predict that the longitudinal size of a ring
polymer scales with the chain length and the diameter of the channel in the
same manner as that for linear chains based on scaling analysis and Flory-type
theory. Moreover, Flory-type theory also gives the ratio of the longitudinal
sizes for a ring polymer and a linear chain with identical chain length. These
theoretical predictions are confirmed by numerical simulations. Finally, our
simulation results show that this ratio first decreases and then saturates with
increasing the chain stiffness, which has interpreted the discrepancy in
experiments. Our results have biological significance.


  Using Langevin dynamics simulations, we investigate the dynamics of a
flexible polymer translocation into a confined area under a driving force
through a nanopore. We choose an ellipsoidal shape for the confinement and
consider the dependence of the asymmetry of the ellipsoid measured by the
aspect ratio on the translocation time. Compared with an isotropic confinement
(sphere), an anisotropic confinement (ellipsoid) with the same volume slows
down the translocation, and the translocation time increases with increasing
the aspect ratio of the ellipsoid. We further find that it takes different time
for polymer translocation into the same ellipsoid through major-axis and
minor-axis directions, depending on the average density of the whole chain in
the ellipsoid, $\phi$. For $\phi$ lower than a critical value $\phi_c$, the
translocation through minor axis is faster, and vice versa. These complicated
behaviors are interpreted by the degree of the confinement and anisotropic
confinement induced folding of the translocated chain.


  We investigate the translocation of stiff polymers in the presence of binding
particles through a nanopore by two-dimensional Langevin dynamics simulations.
We find that the mean translocation time shows a minimum as a function of the
binding energy $\epsilon$ and the particle concentration $\phi$, due to the
interplay of the force from binding and the frictional force. Particularly, for
the strong binding the translocation proceeds with a decreasing translocation
velocity induced by a significant increase of the frictional force. In
addition, both $\epsilon$ and $\phi$ have an notable impact on the distribution
of the translocation time. With increasing $\epsilon$ and $\phi$, it undergoes
a transition from an asymmetric and broad distribution under the weak binding
to a nearly Gaussian one under the strong binding, and its width becomes
gradually narrower.


  Disease awareness in epidemiology can be modelled with adaptive contact
networks, where the interplay of disease dynamics and network alteration often
adds new phases to the standard models (Gross et al. 2006, Shaw et al. 2008)
and, in stochastic simulations, lets network topology settle down to a steady
state that can be static (in the frozen phase) or dynamic (in the endemic
phase). We show for the SIS model that, in the endemic phase, this steady state
does not depend on the initial network topology, only on the disease and
rewiring parameters and on the link density of the network, which is conserved.
We give an analytic description of the structure of this co-evolving network of
infection through its steady-state degree distribution.


  Iterated games are a fundamental component of economic and evolutionary game
theory. They describe situations where two players interact repeatedly and have
the possibility to use conditional strategies that depend on the outcome of
previous interactions. In the context of evolution of cooperation, repeated
games represent the mechanism of reciprocation. Recently a new class of
strategies has been proposed, so called 'zero determinant strategies'. These
strategies enforce a fixed linear relationship between one's own payoff and
that of the other player. A subset of those strategies are 'extortioners' which
ensure that any increase in the own payoff exceeds that of the other player by
a fixed percentage. Here we analyze the evolutionary performance of this new
class of strategies. We show that in reasonably large populations they can act
as catalysts for the evolution of cooperation, similar to tit-for-tat, but they
are not the stable outcome of natural selection. In very small populations,
however, relative payoff differences between two players in a contest matter,
and extortioners hold their ground. Extortion strategies do particularly well
in co-evolutionary arms races between two distinct populations: significantly,
they benefit the population which evolves at the slower rate - an instance of
the so-called Red King effect. This may affect the evolution of interactions
between host species and their endosymbionts.


  Regulatory networks have evolved to allow gene expression to rapidly track
changes in the environment as well as to buffer perturbations and maintain
cellular homeostasis in the absence of change. Theoretical work and empirical
investigation in Escherichia coli have shown that negative autoregulation
confers both rapid response times and reduced intrinsic noise, which is
reflected in the fact that almost half of Escherichia coli transcription
factors are negatively autoregulated. However, negative autoregulation is
exceedingly rare amongst the transcription factors of Saccharomyces cerevisiae.
This difference is all the more surprising because E. coli and S. cerevisiae
otherwise have remarkably similar profiles of network motifs. In this study we
first show that regulatory interactions amongst the transcription factors of
Drosophila melanogaster and humans have a similar dearth of negative
autoregulation to that seen in S. cerevisiae. We then present a model
demonstrating that this fundamental difference in the noise reduction
strategies used amongst species can be explained by constraints on the
evolution of negative autoregulation in diploids. We show that regulatory
interactions between pairs of homologous genes within the same cell can lead to
under-dominance - mutations which result in stronger autoregulation, and
decrease noise in homozygotes, paradoxically can cause increased noise in
heterozygotes. This severely limits a diploid's ability to evolve negative
autoregulation as a noise reduction mechanism. Our work offers a simple and
general explanation for a previously unexplained difference between the
regulatory architectures of E. coli and yeast, Drosophila and humans. It also
demonstrates that the effects of diploidy in gene networks can have
counter-intuitive consequences that may profoundly influence the course of
evolution.


  We revisit here the mathematical model for ATP production in mitochondria
introduced recently by Bertram, Pedersen, Luciani, and Sherman (BPLS) as a
simplification of the more complete but intricate Magnus and Keizer's model. We
identify some inaccuracies in the BPLS original approximations for two flux
rates, namely the adenine nucleotide translocator rate $J_{\rm ANT}$ and the
calcium uniporter rate $J_{\rm uni}$. We introduce new approximations for such
flux rates and then analyze some of the dynamical properties of the model. We
infer, from exhaustive numerical explorations, that the enhanced BPLS equations
have a unique attractor fixed point for physiologically acceptable ranges of
mitochondrial variables and respiration inputs, as one would indeed expect from
homeostasis. We determine, in the stationary regime, the dependence of the
mitochondrial variables on the respiration inputs, namely the cytosolic
concentration of calcium ${\rm Ca}_{\rm c}$ and the substrate fructose
1,6-bisphosphate FBP. The same dynamical effects of calcium and FBP saturations
reported for the original BPLS model are observed here. We find out, however, a
novel non-stationary effect which could be, in principle, physiologically
interesting: some response times of the model tend to increase considerably for
high concentrations of calcium and/or FBP. In particular, the larger the
concentrations of ${\rm Ca}_{\rm c}$ and/or FBP, the larger the necessary time
to attain homeostasis.


  Many of the stochastic models used in inference of phylogenetic trees from
biological sequence data have polynomial parameterization maps. The image of
such a map --- the collection of joint distributions for a model --- forms the
model space. Since the parameterization is polynomial, the Zariski closure of
the model space is an algebraic variety which is typically much larger than the
model space, but has been usefully studied with algebraic methods. Of ultimate
interest, however, is not the full variety, but only the model space. Here we
develop complete semialgebraic descriptions of the model space arising from the
k-state general Markov model on a tree, with slightly restricted parameters.
Our approach depends upon both recently-formulated analogs of Cayley's
hyperdeterminant, and the construction of certain quadratic forms from the
joint distribution whose positive (semi-)definiteness encodes information about
parameter values. We additionally investigate the use of Sturm sequences for
obtaining similar results.


  Cells often have tens of thousands of receptors, even though only a few
activated receptors can trigger full cellular responses. Reasons for the
overabundance of receptors remain unclear. We suggest that, in certain
conditions, the large number of receptors results in a competition among
receptors to be the first to activate the cell. The competition decreases the
variability of the time to cellular activation, and hence results in a more
synchronous activation of cells. We argue that, in simple models, this
variability reduction does not necessarily interfere with the receptor
specificity to ligands achieved by the kinetic proofreading mechanism. Thus
cells can be activated accurately in time and specifically to certain signals.
We predict the minimum number of receptors needed to reduce the coefficient of
variation for the time to activation following binding of a specific ligand.
Further, we predict the maximum number of receptors so that the kinetic
proofreading mechanism still can improve the specificity of the activation.
These predictions fall in line with experimentally reported receptor numbers
for multiple systems.


  Cycling tissues such as the intestinal epithelium, germ line, and hair
follicles, require a constant flux of differentiated cells. These tissues are
maintained by a population of stem cells, which generate differentiated
progenies and self-renew. Asymmetric division of each stem cell into one stem
cell and one differentiated cell can accomplish both tasks. However, in
mammalian cycling tissues, some stem cells divide symmetrically into two
differentiated cells and are replaced by a neighbor that divides symmetrically
into two stem cells. Besides this heterogeneity in fate (population asymmetry),
stem cells also exhibit heterogenous proliferation-rates; in the long run,
however, all stem cells proliferate at the same average rate (equipotency). We
construct and simulate a mathematical model based on these experimental
observations. We show that the complex steady-state dynamics of
population-asymmetric stem cells reduces the rate of replicative aging of the
tissue --potentially lowering the incidence of somatic mutations and genetics
diseases such as cancer. Essentially, slow-dividing stem cells proliferate and
purge the population of the fast-dividing --older-- cells which had undertaken
the majority of the tissue-generation burden. As the number of slow-dividing
cells grows, their proliferation-rate increases, eventually turning them into
fast-dividers, which are themselves replaced by newly emerging slow-dividers.
Going beyond current experiments, we propose a mechanism for equipotency that
can potentially halve the rate of replicative aging. Our results highlight the
importance of a population-level understanding of stem cells, and may explain
the prevalence of population asymmetry in a wide variety of cycling tissues.


  Using the Hamilton-Jacobi equation approach to study genomes of length $L$,
we obtain 1/L corrections for the steady state population distributions and
mean fitness functions for horizontal gene transfer model, as well as for the
diploid evolution model with general fitness landscapes. Our numerical
solutions confirm the obtained analytic equations. Our method could be applied
to the general case of nonlinear Markov models.


  We considered a {multi-block} molecular model of biological evolution, in
which fitness is a function of the mean types of alleles located at different
parts (blocks) of the genome. We formulated an infinite population model with
selection and mutation, and calculated the mean fitness. For the case of
recombination, we formulated a model with a multidimensional fitness landscape
(the dimension of the space is equal to the number of blocks) and derived a
theorem about the dynamics of initially narrow distribution. We also considered
the case of lethal mutations. We also formulated the finite population version
of the model in the case of lethal mutations. Our models, derived for the virus
evolution, are interesting also for the statistical mechanics and the
Hamilton-Jacobi equation as well.


  We consider finite population size effects for Crow-Kimura and Eigen
quasispecies models with single peak fitness landscape. We formulate accurately
the iteration procedure for the finite population models, then derive
Hamilton-Jacobi equation (HJE) to describe the dynamic of the probability
distribution. The steady state solution of HJE gives the variance of the mean
fitness. Our results are useful for understanding population sizes of virus in
which the infinite population models can give reliable results for the
biological evolution problems.


  We present atomistic molecular dynamics simulations of two Polyethylene
systems where all entanglements are trapped: a perfect network, and a melt with
grafted chain ends. We examine microscopically at what level topological
constraints can be considered as a collective entanglement effect, as in tube
model theories, or as certain pairwise uncrossability interactions, as in
slip-link models. A pairwise parameter, which varies between these limiting
cases, shows that, for the systems studied, the character of the entanglement
environment is more pairwise than collective.
  We employ a novel methodology, which analyzes entanglement constraints into a
complete set of pairwise interactions, similar to slip links. Entanglement
confinement is assembled by a plethora of links, with a spectrum of confinement
strengths, from strong to weak. The strength of interactions is quantified
through a link `persistence', which is the fraction of time for which the links
are active. By weighting links according to their strength, we show that
confinement is imposed mainly by the strong ones, and that the weak, trapped,
uncrossability interactions cannot contribute to the low frequency modulus of
an elastomer, or the plateau modulus of a melt.
  A self-consistent scheme for mapping topological constraints to specific,
strong binary links, according to a given entanglement density, is proposed and
validated. Our results demonstrate that slip links can be viewed as the
strongest pairwise interactions of a collective entanglement environment. The
methodology developed provides a basis for bridging the gap between atomistic
simulations and mesoscopic slip link models.


  Eigen's quasispecies system with explicit space and global regulation is
considered. Limit behavior and stability of the system in a functional space
under perturbations of a diffusion matrix with nonnegative spectrum are
investigated. It is proven that if the diffusion matrix has only positive
eigenvalues then the solutions of the distributed system converge to the
equilibrium solution of the corresponding local dynamical system. These results
imply that the error threshold does not change if the spatial interactions
under the principle of global regulation are taken into account.


  The genetic code underlying protein synthesis is a canonical example of a
degenerate biological system. Degeneracies in physical and biological systems
can be lifted by external perturbations thus allowing degenerate systems to
exhibit a wide range of behaviors. Here we show that the degeneracy of the
genetic code is lifted by environmental perturbations to regulate protein
levels in living cells. By measuring protein synthesis rates from a synthetic
reporter library in Escherichia coli, we find that environmental perturbations,
such as reduction of cognate amino acid supply, lift the degeneracy of the
genetic code by splitting codon families into a hierarchy of robust and
sensitive synonymous codons. Rates of protein synthesis associated with robust
codons are up to hundred-fold higher than those associated with sensitive
codons under these conditions. We find that the observed hierarchy between
synonymous codons is not determined by usual rules associated with tRNA
abundance and codon usage. Rather, competition among tRNA isoacceptors for
aminoacylation underlies the robustness of protein synthesis. Remarkably, the
hierarchy established using the synthetic library also explains the measured
robustness of synthesis for endogenous proteins in E. coli. We further found
that the same hierarchy is reflected in the fitness cost of synonymous
mutations in amino acid biosynthesis genes and in the transcriptional control
of sigma factor genes. Our study reveals that the degeneracy of the genetic
code can be lifted by environmental perturbations, and it suggests that
organisms can exploit degeneracy lifting as a general strategy to adapt protein
synthesis to their environment.


  We present Brownian dynamics simulations of the facilitated diffusion of a
protein, modelled as a sphere with a binding site on its surface, along DNA,
modelled as a semi-flexible polymer. We consider both the effect of DNA
organisation in 3D, and of sequence heterogeneity. We find that in a network of
DNA loops, as are thought to be present in bacterial DNA, the search process is
very sensitive to the spatial location of the target within such loops.
Therefore, specific genes might be repressed or promoted by changing the local
topology of the genome. On the other hand, sequence heterogeneity creates traps
which normally slow down facilitated diffusion. When suitably positioned,
though, these traps can, surprisingly, render the search process much more
efficient.


  The formation of fibrillar aggregates seems to be a common characteristic of
polypeptide chains, although the observation of these aggregates may depend on
appropriate experimental conditions. Partially folded intermediates seem to
have an important role in the generation of protein aggregates, and a mechanism
for this fibril formation considers that these intermediates also correspond to
metastable states with respect to the fibrillar ones. Here, using a
coarse-grained (CG) off-lattice model, we carry out a comparative analysis of
the thermodynamic aspects characterizing the folding transition with respect to
the propensity for aggregation of four different systems: two isoforms of the
amyloid $\beta$-protein, the Src SH3 domain, and the human prion proteins
(hPrP). Microcanonical analysis of the data obtained from replica exchange
method (REM) is conducted to evaluate the free-energy barrier and latent heat
in these models. The simulations of the amyloid $\beta$ isoforms and Src SH3
domain indicated that the folding process described by this CG model is related
to a negative specific heat, a phenomenon that can only be verified in the
microcanonical ensemble in first-order phase transitions. The CG simulation of
the hPrP heteropolymer yielded a continuous folding transition. The absence of
a free-energy barrier and latent heat favors the presence of partially unfolded
conformations, and in this context, this thermodynamic aspect could explain the
reason why the hPrP heteropolymer is more aggregation-prone than the other
heteropolymers considered in this study. We introduced the hydrophobic radius
of gyration as an order parameter and found that it can be used to obtain
reliable information about the hydrophobic packing and the transition
temperatures in the folding process.


  We present a model for host-parasite dynamics which incorporates both
vertical and horizontal transmission as well as spatial structure. Our model
consists of stacked contact processes (CP), where the dynamics of the host is a
simple CP on a lattice while the dynamics of the parasite is a secondary CP
which sits on top of the host-occupied sites. In the simplest case, where
infection does not incur any cost, we uncover a novel effect: a nonmonotonic
dependence of parasite prevalence on host turnover. Inspired by natural
examples of hyperparasitism, we extend our model to multiple levels of
parasites and identify a transition between the maintenance of a finite and
infinite number of levels, which we conjecture is connected to a roughening
transition in models of surface-growth.


  The Roma people, living throughout Europe, are a diverse population linked by
the Romani language and culture. Previous linguistic and genetic studies have
suggested that the Roma migrated into Europe from South Asia about 1000-1500
years ago. Genetic inferences about Roma history have mostly focused on the Y
chromosome and mitochondrial DNA. To explore what additional information can be
learned from genome-wide data, we analyzed data from six Roma groups that we
genotyped at hundreds of thousands of single nucleotide polymorphisms (SNPs).
We estimate that the Roma harbor about 80% West Eurasian ancestry-deriving from
a combination of European and South Asian sources- and that the date of
admixture of South Asian and European ancestry was about 850 years ago. We
provide evidence for Eastern Europe being a major source of European ancestry,
and North-west India being a major source of the South Asian ancestry in the
Roma. By computing allele sharing as a measure of linkage disequilibrium, we
estimate that the migration of Roma out of the Indian subcontinent was
accompanied by a severe founder event, which we hypothesize was followed by a
major demographic expansion once the population arrived in Europe.


  How do single cell fate decisions induced by activation of key signaling
proteins above threshold concentrations within a time interval are affected by
stochastic fluctuations in biochemical reactions? We address this question
using minimal models of stochastic chemical reactions commonly found in cell
signaling and gene regulatory systems. Employing exact solutions and
semi-analytical methods we calculate distributions of the maximum value ($N$)
of activated species concentrations ($P_{max}(N)$) and the time ($t$) taken to
reach the maximum value ($P_{max}(t)$) within a time window in the minimal
models. We find, the presence of positive feedback interactions make
$P_{max}(N)$ more spread out with a higher "peakedness" in $P_{max}(t)$. Thus
positive feedback interactions may help single cells to respond sensitively to
a stimulus when cell decision processes require upregulation of activated forms
of key proteins to a threshold number within a time window.


  Coalescent processes, including mutation, are derived from Moran type
population models admitting large offspring numbers. Including mutation in the
coalescent process allows for quantifying the turnover of alleles by computing
the distribution of the number of original alleles still segregating in the
population at a given time in the past. The turnover of alleles is considered
for specific classes of the Moran model admitting large offspring numbers.
Versions of the Kingman coalescent are also derived whose rates are functions
of the mean and variance of the offspring distribution. High variance in the
offspring distribution results in higher turnover and younger age of alleles
than predicted by the usual Kingman coalescent.


  Unraveling the structure of complex biological networks and relating it to
their functional role is an important task in systems biology. Here we attempt
to characterize the functional organization of the large-scale metabolic
networks of three microorganisms. We apply flux balance analysis to study the
optimal growth states of these organisms in different environments. By
investigating the differential usage of reactions across flux patterns for
different environments, we observe a striking bimodal distribution in the
activity of reactions. Motivated by this, we propose a simple algorithm to
decompose the metabolic network into three sub-networks. It turns out that our
reaction classifier which is blind to the biochemical role of pathways leads to
three functionally relevant sub-networks that correspond to input, output and
intermediate parts of the metabolic network with distinct structural
characteristics. Our decomposition method unveils a functional bow-tie
organization of metabolic networks that is different from the bow-tie structure
determined by graph-theoretic methods that do not incorporate functionality.


  Reconciling a gene tree with a species tree is an important task that reveals
much about the evolution of genes, genomes, and species, as well as about the
molecular function of genes. A wide array of computational tools have been
devised for this task under certain evolutionary events such as hybridization,
gene duplication/loss, or incomplete lineage sorting. Work on reconciling gene
tree with species phylogenies under two or more of these events have also begun
to emerge. Our group recently devised both parsimony and probabilistic
frameworks for reconciling a gene tree with a phylogenetic network, thus
allowing for the detection of hybridization in the presence of incomplete
lineage sorting. While the frameworks were general and could handle any
topology, they are computationally intensive, rendering their application to
large datasets infeasible. In this paper, we present two novel approaches to
address the computational challenges of the two frameworks that are based on
the concept of ancestral configurations. Our approaches still compute exact
solutions while improving the computational time by up to five orders of
magnitude. These substantial gains in speed scale the applicability of these
unified reconciliation frameworks to much larger data sets. We discuss how the
topological features of the gene tree and phylogenetic network may affect the
performance of the new algorithms. We have implemented the algorithms in our
PhyloNet software package, which is publicly available in open source.


  My analysis uses methods developed for data mining microarray experiments,
adapted for ageing research. Methods bridge knowledge of statistical mechanics
with data mining methods developed in statistical mathematics. Analyses can
reveal how the transcriptional regulation of genes might coincide, thereby
implicating proteins as parts of networks acting together towards a common
biological function. Such experiments are most useful for complex biological
traits that result from the presumed functioning of several molecular pathways.
Ageing is one such biological phenomenon that apparently incorporates numerous
molecular mechanisms underlying environmental stimulus sensing, metabolic
regulation, stress responses, reproductive signalling, hibernation, and
transcriptional regulation.


  Age-related muscle decline, a condition referred to as sarcopenia and defined
as loss in muscle mass and muscle strength over time, is one of the most
pervasive problems of the elderly, such that significant declines in strength
and mobility affects essentially every old person We have found that ageing C.
elegans body wall muscle undergoes a process remarkably reminiscent of human
sarcopenia. Both have mid-life onset and are characterized by progressive loss
of sarcomeres and cytoplasmic volume; both are associated with locomotory
decline. To extend understanding of this fundamental problem, I surveyed
expression of all known muscle related genes to describe a profile of
transcriptional changes in muscle that transpires during adult life and ageing.
Importantly, the intersection of this dataset with that from ageing flies and
some human studies can suggest conserved genes that might impact the process
most strongly. Hypotheses I formulate will be used to drive experiments at the
bench and perhaps to focus attention for human therapies.


  Motivation: Protein surface roughness is fractal in nature. Mass,
hydrophobicity, polarizability distributions of protein interior are fractal
too, as are the distributions of dipole moments, aromatic residues, and many
other structural determinants. The open-access server ProtFract, presents a
reliable way to obtain numerous fractal-dimension and correlation-dimension
based results to quantify the symmetry of self-similarity in distributions of
various properties of protein interior and exterior.
  Results: Fractal dimension based comparative analyses of various biophysical
properties of Ras superfamily proteins were conducted. Though the extent of
sequence and functional overlapping across Ras superfamily structures is
extremely high, results obtained from ProtFract investigation are found to be
sensitive to detect differences in the distribution of each of the properties.
For example, it was found that the RAN proteins are structurally most stable
amongst all Ras superfamily proteins, the ARFs possess maximum extent of unused
hydrophobicity in their structures, RAB protein interiors have
electrostatically least conducive environment, GEM/REM/RAD proteins possess
exceptionally high symmetry in the structural organization of their active
chiral centres, neither hydrophobicity nor columbic interactions play
significant part in stabilizing the RAS proteins but aromatic interactions do,
though cation-pi interactions are found to be more dominant in RAN than in RAS
proteins. Ras superfamily proteins can best be classified with respect to their
class-specific pi-pi and cation-pi interaction symmetries.
  Availability: ProtFract is freely available online at the URL:
http://bioinfo.net.in/protfract/index.html


  We consider Markovian susceptible-infectious-removed (SIR) dynamics on
time-invariant weighted contact networks where the infection and removal
processes are Poisson and where network links may be directed or undirected. We
prove that a particular pair-based moment closure representation generates the
expected infectious time series for networks with no cycles in the underlying
graph. Moreover, this ``deterministic'' representation of the expected
behaviour of a complex heterogeneous and finite Markovian system is
straightforward to evaluate numerically.


  We develop statistical enumeration methods for self-avoiding walks using a
powerful sampling technique called the multicanonical Monte Carlo method. Using
these methods, we estimate the numbers of the two dimensional N-step
self-avoiding walks up to N=256 with statistical errors. The developed methods
are based on statistical mechanical models of paths which include self-avoiding
walks. The criterion for selecting a suitable model for enumerating
self-avoiding walks is whether or not the configuration space of the model
includes a set for which the number of the elements can be exactly counted. We
call this set a scale fixing set. We selected the following two models which
satisfy the criterion: the Go model for lattice proteins and the Domb-Joyce
model for generalized random walks. There is a contrast between these two
models in the structures of the configuration space. The configuration space of
the Go model is defined as the universal set of self-avoiding walks, and the
set of the ground state conformation provides a scale fixing set. On the other
hand, the configuration space of the Domb-Joyce model is defined as the
universal set of random walks which can be used as a scale fixing set, and the
set of the ground state conformation is the same as the universal set of
self-avoiding walks. From the perspective of enumeration performance, we
conclude that the Domb-Joyce model is the better of the two. The reason for the
performance difference is partly explained by the existence of the first-order
phase transition of the Go model.


  In order to establish defined biomimetic systems, type I collagen was
functionalised with 1,3-Phenylenediacetic acid (Ph) as aromatic, bifunctional
segment. Following investigation on molecular organization and macroscopic
properties, material functionalities, i.e. degradability and bioactivity, were
addressed, aiming at elucidating the potential of this collagen system as
mineralization template. Functionalised collagen hydrogels demonstrated a
preserved triple helix conformation. Decreased swelling ratio and increased
thermo-mechanical properties were observed in comparison to state-of-the-art
carbodiimide (EDC)-crosslinked collagen controls. Ph-crosslinked samples
displayed no optical damage and only a slight mass decrease (~ 4 wt.-%)
following 1-week incubation in simulated body fluid (SBF), while nearly 50
wt.-% degradation was observed in EDC-crosslinked collagen. SEM/EDS revealed
amorphous mineral deposition, whereby increased calcium phosphate ratio was
suggested in hydrogels with increased Ph content. This investigation provides
valuable insights for the synthesis of triple helical collagen materials with
enhanced macroscopic properties and controlled degradation. In light of these
features, this system will be applied for the design of tissue-like scaffolds
for mineralized tissue formation.


  Stochastic resonance (SR) is a well known phenomenon in dynamical systems. It
consists of the amplification and optimization of the response of a system
assisted by stochastic noise. Here we carry out the first experimental study of
SR in single DNA hairpins which exhibit cooperatively folding/unfolding
transitions under the action of an applied oscillating mechanical force with
optical tweezers. By varying the frequency of the force oscillation, we
investigated the folding/unfolding kinetics of DNA hairpins in a periodically
driven bistable free-energy potential. We measured several SR quantifiers under
varied conditions of the experimental setup such as trap stiffness and length
of the molecular handles used for single-molecule manipulation. We find that
the signal-to-noise ratio (SNR) of the spectral density of measured
fluctuations in molecular extension of the DNA hairpins is a good quantifier of
the SR. The frequency dependence of the SNR exhibits a peak at a frequency
value given by the resonance matching condition. Finally, we carried out
experiments in short hairpins that show how SR might be useful to enhance the
detection of conformational molecular transitions of low SNR.


  The domestication and subsequent selection by humans to create breeds of
cattle undoubtedly altered the patterning of variation within their genomes.
Strong selection to fix advantageous large-effect mutations underlying
domesticability, breed characteristics or productivity created selective sweeps
in which variation was lost in the chromosomal region flanking the selected
allele. Selective sweeps have been identified in the genomes of many species
including humans, dogs, horses, and chickens. We attempt to identify regions of
the bovine genome that have been subjected to selective sweeps. Two datasets
were used for the discovery and validation of selective sweeps via the fixation
of alleles at a series of contiguous SNP loci. BovineSNP50 data were used to
identify 28 putative sweep regions among 14 cattle breeds. Affymetrix BOS 1
prescreening assay data for five breeds were used to identify 114 regions and
validate 5 regions identified using the BovineSNP50 data. Many genes are
located within these regions; however, phenotypes that we predict to have
historically been under strong selection include horned-polled, coat color,
stature, ear morphology, and behavior. The identified selective sweeps
represent recent events associated with breed formation rather than ancient
events associated with domestication. No sweep regions were shared between
indicine and taurine breeds reflecting their divergent selection histories. A
primary finding of this study is the sensitivity of results to assay
resolution. Despite the bias towards common SNPs in the BovineSNP50 design,
false positive sweep regions appear to be common due to the limited resolution
of the assay. This assay design bias leads to the detection of breed-specific
sweep regions, or regions shared by a small number of breeds, restricting the
suite of selected phenotypes detected to primarily those associated with breed
characteristics.


  In this paper we propose a model-order reduction method for chemical reaction
networks governed by general enzyme kinetics, including the mass-action and
Michaelis-Menten kinetics. The model-order reduction method is based on the
Kron reduction of the weighted Laplacian matrix which describes the graph
structure of complexes in the chemical reaction network. We apply our method to
a yeast glycolysis model, where the simulation result shows that the transient
behaviour of a number of key metabolites of the reduced-order model is in good
agreement with those of the full-order model.


  The recent explosion in available genetic data has led to significant
advances in understanding the demographic histories of and relationships among
human populations. It is still a challenge, however, to infer reliable
parameter values for complicated models involving many populations. Here we
present MixMapper, an efficient, interactive method for constructing
phylogenetic trees including admixture events using single nucleotide
polymorphism (SNP) genotype data. MixMapper implements a novel two-phase
approach to admixture inference using moment statistics, first building an
unadmixed scaffold tree and then adding admixed populations by solving systems
of equations that express allele frequency divergences in terms of mixture
parameters. Importantly, all features of the model, including topology, sources
of gene flow, branch lengths, and mixture proportions, are optimized
automatically from the data and include estimates of statistical uncertainty.
MixMapper also uses a new method to express branch lengths in easily
interpretable drift units. We apply MixMapper to recently published data for
HGDP individuals genotyped on a SNP array designed especially for use in
population genetics studies, obtaining confident results for 30 populations, 20
of them admixed. Notably, we confirm a signal of ancient admixture in European
populations---including previously undetected admixture in Sardinians and
Basques---involving a proportion of 20--40% ancient northern Eurasian ancestry.


  Cryptic genetic sequences have attenuated effects on phenotypes. In the
classic view, relaxed selection allows cryptic genetic diversity to build up
across individuals in a population, providing alleles that may later contribute
to adaptation when co-opted - e.g. following a mutation increasing expression
from a low, attenuated baseline. This view is described, for example, by the
metaphor of the spread of a population across a neutral network in genotype
space. As an alternative view, consider the fact that most phenotypic traits
are affected by multiple sequences, including cryptic ones. Even in a strictly
clonal population, the co-option of cryptic sequences at different loci may
have different phenotypic effects and offer the population multiple adaptive
possibilities. Here, we model the evolution of quantitative phenotypic
characters encoded by cryptic sequences, and compare the relative contributions
of genetic diversity and of variation across sites to the phenotypic potential
of a population. We show that most of the phenotypic variation accessible
through co-option would exist even in populations with no polymorphism. This is
made possible by a history of compensatory evolution, whereby the phenotypic
effect of a cryptic mutation at one site was balanced by mutations elsewhere in
the genome, leading to a diversity of cryptic effect sizes across sites rather
than across individuals. Cryptic sequences might accelerate adaptation and
facilitate large phenotypic changes even in the absence of genetic diversity,
as traditionally defined in terms of alternative alleles.


  Upon studying the B-Factors of all the atoms of all non-redundant proteins
belonging to 76 most commonly found structural domains of all four major
structural classes, it was found that the residue mobility has decreased during
the course of evolution. Though increased residue-flexibility was preferred in
the early stages of protein structure evolution, less flexibility is preferred
in the medieval and recent stages. GLU is found to be the most flexible residue
while VAL recorded to have the least flexibility. General trends in decrement
of B-Factors conformed to the general trend in the order of emergence of
protein structural domains. Decrement of B-Factor is observed to be most
decisive (monotonic and uniform) for VAL, while evolution of CYS and LYS
flexibility is found to be most skewed. Barring CYS, flexibility of all the
residues is found to have increased during evolution of alpha by beta folds,
however flexibility of all the residues (barring CYS) is found to have
decreased during evolution of all beta folds. Only in alpha by beta folds the
tendency of preferring higher residue mobility could be observed, neither alpha
plus beta, nor all alpha nor all beta folds were found to support higher
residue-mobility. In all the structural classes, the effect of evolutionary
constraint on polar residues is found to follow an exactly identical trend as
that on hydrophobic residues, only the extent of these effects are found to be
different. Though protein size is found to be decreasing during evolution,
residue mobility of proteins belonging to ancient and old structural domains
showed strong positive dependency upon protein size, however for medieval and
recent domains such dependency vanished. It is found that to optimize residue
fluctuations, alpha by beta class of proteins are subjected to more stringent
evolutionary constraints.


  The large volumes of sequencing data required to sample complex environments
deeply pose new challenges to sequence analysis approaches. De novo metagenomic
assembly effectively reduces the total amount of data to be analyzed but
requires significant computational resources. We apply two pre-assembly
filtering approaches, digital normalization and partitioning, to make large
metagenome assemblies more comput\ ationaly tractable. Using a human gut mock
community dataset, we demonstrate that these methods result in assemblies
nearly identical to assemblies from unprocessed data. We then assemble two
large soil metagenomes from matched Iowa corn and native prairie soils. The
predicted functional content and phylogenetic origin of the assembled contigs
indicate significant taxonomic differences despite similar function. The
assembly strategies presented are generic and can be extended to any
metagenome; full source code is freely available under a BSD license.


  With the aim of establishing a criterion for identifying when a lipid bilayer
has reached steady state using the molecular dynamics simulation technique,
lipid bilayers of different composition in their liquid crystalline phase were
simulated in aqueous solution in presence of CaCl$_2$ as electrolyte, at
different concentration levels. In this regard, we used two different lipid
bilayer systems: one composed by 288 DPPC (DiPalmitoylPhosphatidylCholine) and
another constituted by 288 DPPS (DiPalmitoylPhosphatidylSerine). In this sense,
for both type of lipid bilayers, we have studied the temporal evolution of some
lipids properties, such as the surface area per lipid, the deuterium order
parameter, the lipid hydration and the lipid-calcium coordination. From their
analysis, it became evident how each property has a different time to achieve
equilibrium. The following order was found, from faster property to slower
property: coordination of ions $\approx$ deuterium order parameter $>$ area per
lipid $\approx$ hydration. Consequently, when the hydration of lipids or the
mean area per lipid are stable, we can ensure that the lipid membrane has
reached the steady state.


  Ward and Kellis (Reports, September 5 2012) identify regulatory regions in
the human genome exhibiting lineage-specific constraint and estimate the extent
of purifying selection. There is no statistical rationale for the examples they
highlight, and their estimates of the fraction of the genome under constraint
are biased by arbitrary designations of completely constrained regions.


  Experimental data that prove the existence of the zygotic combinatorial
process occurring in an embryogenesis-entering zygote are presented in the
paper. The zygotic combinatorial process is found when analyzing F1 hybrid
plants obtained from crossing homozygous forms different, minimum, in two
marker enzymes, and it is found in that hybrid plant which, with one marker
enzyme heterozygous spectrum, has a homozygous spectrum of the other. The
zygotic combinatorial process leads to F1 hybrids uniformity aberration. The
zygotic combinatory process revealed in the study is supposed to be conditioned
by chromosome polyteny in mother plant cells and diminution of chromatin excess
from the embryogenesis-entering zygote. An obligatory condition for
combinatorial process is the presence of free exchange of cromatides among
homological chromosomes in an embryogenesis-entering cell, i.e. the presence of
crossing-over analogous to the one proceeding at meiosis. The found
combinatorial process and the earlier-obtained data confirm the hypothesis on
multi-dimensionality of inherited information coding. Differential polyteny of
certain chromosome regions can lead to differences among homozygous plants
having the same alleles in genes located in polytenized regions and controlling
morpho-physiological traits.


  Over the years numerous models of SIS (susceptible - infected - susceptible)
disease dynamics unfolding on networks have been proposed. Here, we discuss the
links between many of these models and how they can be viewed as more general
motif-based models. We illustrate how the different models can be derived from
one another and, where this is not possible, discuss extensions to established
models that enables this derivation. We also derive a general result for the
exact differential equations for the expected number of an arbitrary motif
directly from the Kolmogorov/master equations and conclude with a comparison of
the performance of the different closed systems of equations on networks of
varying structure.


  Studies of coevolution of amino acids within and between proteins have
revealed two types of coevolving units: coevolving contacts, which are pairs of
amino acids distant along the sequence but in contact in the three-dimensional
structure, and sectors, which are larger groups of structurally connected amino
acids that underlie the biochemical properties of proteins. By reconciling two
approaches for analyzing correlations in multiple sequence alignments, we link
these two findings together and with coevolving units of intermediate size,
called `sectons', which are shown to provide additional information. By
extending the analysis to the co-occurrence of orthologous genes in bacterial
genomes, we also show that the methods and results are general and relevant
beyond protein structures.


  Background: Predictive, stable and interpretable gene signatures are
generally seen as an important step towards a better personalized medicine.
During the last decade various methods have been proposed for that purpose.
However, one important obstacle for making gene signatures a standard tool in
clinics is the typical low reproducibility of these signatures combined with
the difficulty to achieve a clear biological interpretation. For that purpose
in the last years there has been a growing interest in approaches that try to
integrate information from molecular interaction networks. Results: We propose
a novel algorithm, called FrSVM, which integrates protein-protein interaction
network information into gene selection for prognostic biomarker discovery. Our
method is a simple filter based approach, which focuses on central genes with
large differences in their expression. Compared to several other competing
methods our algorithm reveals a significantly better prediction performance and
higher signature stability. More- over, obtained gene lists are highly enriched
with known disease genes and drug targets. We extendd our approach further by
integrating information on candidate disease genes and targets of disease
associated Transcript Factors (TFs).


  We introduce an extension to the Weighted Ensemble (WE) path sampling method
to restrict sampling to a one dimensional path through a high dimensional phase
space. Our method, which is based on the finite-temperature string method,
permits efficient sampling of both equilibrium and non-equilibrium systems.
Sampling obtained from the WE method guides the adaptive refinement of a
Voronoi tessellation of order parameter space, whose generating points, upon
convergence, coincide with the principle reaction pathway. We demonstrate the
application of this method to several simple, two-dimensional models of driven
Brownian motion and to the conformational change of the nitrogen regulatory
protein C receiver domain using an elastic network model. The simplicity of the
two-dimensional models allows us to directly compare the efficiency of the WE
method to conventional brute force simulations and other path sampling
algorithms, while the example of protein conformational change demonstrates how
the method can be used to efficiently study transitions in the space of many
collective variables.


  When an epidemic spreads in a population, individuals may adaptively change
the structure of their social contact network to reduce risk of infection. Here
we study the spread of an epidemic on an adaptive network with community
structure. We model the effect of two communities with different average
degrees. The disease model is susceptible-infected-susceptible (SIS), and
adaptation is rewiring of links between susceptibles and infectives. The
bifurcation structure is obtained, and a mean field model is developed that
accurately predicts the steady state behavior of the system. We show that an
epidemic can alter the community structure.


  Various approaches have explored the covariation of residues in
multiple-sequence alignments of homologous proteins to extract functional and
structural information. Among those are principal component analysis (PCA),
which identifies the most correlated groups of residues, and direct coupling
analysis (DCA), a global inference method based on the maximum entropy
principle, which aims at predicting residue-residue contacts. In this paper,
inspired by the statistical physics of disordered systems, we introduce the
Hopfield-Potts model to naturally interpolate between these two approaches. The
Hopfield-Potts model allows us to identify relevant 'patterns' of residues from
the knowledge of the eigenmodes and eigenvalues of the residue-residue
correlation matrix. We show how the computation of such statistical patterns
makes it possible to accurately predict residue-residue contacts with a much
smaller number of parameters than DCA. This dimensional reduction allows us to
avoid overfitting and to extract contact information from multiple-sequence
alignments of reduced size. In addition, we show that low-eigenvalue
correlation modes, discarded by PCA, are important to recover structural
information: the corresponding patterns are highly localized, that is, they are
concentrated in few sites, which we find to be in close contact in the
three-dimensional protein fold.


  How the very first nervous systems evolved remains a fundamental open
question. Molecular and genomic techniques have revolutionized our knowledge of
the molecular ingredients behind this transition but not yet provided a clear
picture of the morphological and tissue changes involved. Here we focus on a
behavioural perspective that centres on movement by muscle contraction.
Building on the finding that molecules for chemical neural signalling predate
multicellular animals, we investigate a gradual evolutionary scenario for
nervous systems that consists of two stages: A) Chemically transmission of
electrical activity between adjacent cells provided a primitive form of muscle
coordination in a contractile epithelial tissue. B) This primitive form of
coordination was subsequently improved upon by evolving the axodendritic
processes of modern neurons. We use computer simulations to investigate the
first stage. The simulations show that chemical transmission across a
contractile sheet can indeed produce useful body scale patterns, but only for
small-sized animals. For larger animals the noise in chemical neural signalling
interferes. Our results imply that a two-stage scenario is a viable approach to
nervous system evolution. The first stage could provide an initial behavioural
advantage, as well as a clear scaffold for subsequent improvements in
behavioural coordination.


  We investigate the microtubule polymerization dynamics with catastrophe and
rescue events for three different confinement scenarios, which mimic typical
cellular environments: (i) The microtubule is confined by rigid and fixed
walls, (ii) it grows under constant force, and (iii) it grows against an
elastic obstacle with a linearly increasing force. We use realistic catastrophe
models and analyze the microtubule dynamics, the resulting microtubule length
distributions, and force generation by stochastic and mean field calculations;
in addition, we perform stochastic simulations. We also investigate the force
dynamics if growth parameters are perturbed in dilution experiments. Finally,
we show the robustness of our results against changes of catastrophe models and
load distribution factors.


  Motivated by data-rich experiments in transcriptional regulation and sensory
neuroscience, we consider the following general problem in statistical
inference. When exposed to a high-dimensional signal S, a system of interest
computes a representation R of that signal which is then observed through a
noisy measurement M. From a large number of signals and measurements, we wish
to infer the "filter" that maps S to R. However, the standard method for
solving such problems, likelihood-based inference, requires perfect a priori
knowledge of the "noise function" mapping R to M. In practice such noise
functions are usually known only approximately, if at all, and using an
incorrect noise function will typically bias the inferred filter. Here we show
that, in the large data limit, this need for a pre-characterized noise function
can be circumvented by searching for filters that instead maximize the mutual
information I[M;R] between observed measurements and predicted representations.
Moreover, if the correct filter lies within the space of filters being
explored, maximizing mutual information becomes equivalent to simultaneously
maximizing every dependence measure that satisfies the Data Processing
Inequality. It is important to note that maximizing mutual information will
typically leave a small number of directions in parameter space unconstrained.
We term these directions "diffeomorphic modes" and present an equation that
allows these modes to be derived systematically. The presence of diffeomorphic
modes reflects a fundamental and nontrivial substructure within parameter
space, one that is obscured by standard likelihood-based inference.