  Identifying the complete set of functional elements within the human genome
would be a windfall for multiple areas of biological research including
medicine, molecular biology, and evolution. Complete knowledge of function
would aid in the prioritization of loci when searching for the genetic bases of
disease or adaptive phenotypes. Because mutations that disrupt function are
disfavored by natural selection, purifying selection leaves a detectable
signature within functional elements; accordingly this signal has been
exploited through the use of genomic comparisons of distantly related species.
However, the functional complement of the genome changes extensively across
time and between lineages, therefore, evidence of the current action of
purifying selection is essential. Because the removal of deleterious mutations
by natural selection also reduces within-species genetic diversity within
functional loci, dense population genetic data have the potential to reveal
genomic elements that are currently functional. Here we assess the potential of
this approach using 16,411 human mitochondrial genomes. We show that the high
density of polymorphism in this dataset precisely delineates regions
experiencing purifying selection. Further, we show that the number of
segregating alleles at a site is strongly correlated with its divergence across
species after accounting for known mutational biases in human mtDNA. These two
measures track one another at a remarkably fine scale across many loci--a
correlation that is purely the result of natural selection. Our results
demonstrate that genetic variation has the potential to reveal exactly which
nucleotides in the genome are currently performing important functions and
likely to have deleterious fitness effects when mutated. As more complete
genomes are sequenced, similar power to reveal purifying selection may be
achievable in the human nuclear genome.


  Markov state models (MSMs) have been demonstrated to be a powerful method for
computationally studying intramolecular processes such as protein folding and
macromolecular conformational changes. In this article, we present a new
approach to construct MSMs that is applicable to modeling a broad class of
multi-molecular assembly reactions. Distinct structures formed during assembly
are distinguished by their undirected graphs, which are defined by strong
subunit interactions. Spatial inhomogeneities of free subunits are accounted
for using a recently developed Gaussian-based signature. Simplifications to
this state identification are also investigated. The feasibility of this
approach is demonstrated on two different coarse-grained models for virus
self-assembly. We find good agreement between the dynamics predicted by the
MSMs and long, unbiased simulations, and that the MSMs can reduce overall
simulation time by orders of magnitude.


  In the post-genome era, biologists have sought to measure the complete
complement of proteins, termed proteomics. Currently, the most effective method
to measure the proteome is with shotgun, or bottom-up, proteomics, in which the
proteome is digested into peptides that are identified followed by protein
inference. Despite continuous improvements to all steps of the shotgun
proteomics workflow, observed proteome coverage is often low; some proteins are
identified by a single peptide sequence. Complete proteome sequence coverage
would allow comprehensive characterization of RNA splicing variants and all
post translational modifications, which would drastically improve the accuracy
of biological models. There are many reasons for the sequence coverage deficit,
but ultimately peptide length determines sequence observability. Peptides that
are too short are lost because they match many protein sequences and their true
origin is ambiguous. The maximum observable peptide length is determined by
several analytical challenges. This paper explores computationally how peptide
lengths produced from several common proteome digestion methods limit
observable proteome coverage. Iterative proteome cleavage strategies are also
explored. These simulations reveal that maximized proteome coverage can be
achieved by use of an iterative digestion protocol involving multiple proteases
and chemical cleavages that theoretically allow 91.1% proteome coverage.


  Collective phenomena, whereby agent-agent interactions determine spatial
patterns, are ubiquitous in the animal kingdom. On the other hand, movement and
space use are also greatly influenced by the interactions between animals and
their environment. Despite both types of interaction fundamentally influencing
animal behaviour, there has hitherto been no unifying framework for the models
proposed in both areas. Here, we construct a general method for inferring
population-level spatial patterns from underlying individual movement and
interaction processes, a key ingredient in building a statistical mechanics for
ecological systems. We show that resource selection functions, as well as
several examples of collective motion models, arise as special cases of our
framework, thus bringing together resource selection analysis and collective
animal behaviour into a single theory. In particular, we focus on combining the
various mechanistic models of territorial interactions in the literature with
step selection functions, by incorporate interactions into the step selection
framework and demonstrating how to derive territorial patterns from the
resulting models. We demonstrate the efficacy of our model by application to a
population of insectivore birds in the Amazon rainforest.


  F1-ATPase is the soluble portion of the membrane-embedded enzyme FoF1-ATP
synthase that catalyzes the production of adenosine triphosphate in eukaryotic
and eubacterial cells. In reverse, the F1 part can also hydrolyze ATP quickly
at three catalytic binding sites. Therefore, catalysis of 'non-productive' ATP
hydrolysis by F1 (or FoF1) must be minimized in the cell. In bacteria, the
epsilon subunit is thought to control and block ATP hydrolysis by mechanically
inserting its C-terminus into the rotary motor region of F1. We investigate
this proposed mechanism by labeling F1 specifically with two fluorophores to
monitor the C-terminus of the epsilon subunit by F\"orster resonance energy
transfer. Single F1 molecules are trapped in solution by an Anti-Brownian
electrokinetic trap which keeps the FRET-labeled F1 in place for extended
observation times of several hundreds of milliseconds, limited by
photobleaching. FRET changes in single F1 and FRET histograms for different
biochemical conditions are compared to evaluate the proposed regulatory
mechanism.


  In this paper we make the first steps to bridge the gap between classic
control theory and modern, network-based epidemic models. In particular, we
apply nonlinear model predictive control (NMPC) to a pairwise ODE model which
we use to model a susceptible-infectious-susceptible (SIS) epidemic on
non-trivial contact structures. While classic control of epidemics concentrates
on aspects such as vaccination, quarantine and fast diagnosis, our novel setup
allows us to deliver control by altering the contact network within the
population. Moreover, the ideal outcome of control is to eradicate the disease
while keeping the network well connected. The paper gives a thorough and
detailed numerical investigation of the impact and interaction of system and
control parameters on the controllability of the system. The analysis reveals,
that for certain set parameters it is possible to identify critical control
bounds above which the system is controllable. We foresee, that our approach
can be extended to even more realistic or simulation-based models with the aim
to apply these to real-world situations.


  Mutations in proteins can have deleterious effects on a protein's stability
and function, which ultimately causes particular diseases. Genetically
inherited muscular dystrophies (MDs) include several genetic diseases, which
cause increasing weakness in muscles and disability to perform muscular
functions progressively. Different types of mutations in the gene coding
translates into defunct proteins cause different neuro-muscular diseases.
Defunct protein interactions in human proteome may cause a stress to its
neighboring proteins and its modules. We therefore aimed to understand the
effects of mutated proteins on interacting partners in different muscular
dystrophies utilizing network biology to understand system properties of these
MDs subnetworks .We investigated rigidity and flexibility of protein-protein
interaction subnetworks associated with causative mutated genes showing high
mean interference values in muscular dystrophy. Rigid component related to
EEF1A1 subnetwork and members of 14.3.3 protein family formed the core of
network showed involvement in molecular function related to protein domain
specific binding. CACNA1S and CALM1 showing functionality related to
Voltage-dependent calcium channel demonstrated highest flexibility. The
interconnected subnets of proteins corresponding to known causative genes
having large genetic variants are shared in different muscular dystrophies
inferred towards comorbidity in diseases. The studies demonstrates core network
of MDs as highly rigid, constituting of large intermodular edges and
interconnected hub nodes suggesting high information transfer flow. The core
skeleton of the network is organized in protein specific domain binding. This
suggests neuro-muscular disorders may initiate due to interruption in molecular
function related with the core and its aggression may depend on the tolerance
level of the networks.


  One of the central questions of metacommunity theory is how dispersal of
organisms affects species diversity. Here we show that the diversity-dispersal
relationship should not be studied in isolation of other abiotic and biotic
flows in the metacommunity. We study a mechanistic metacommunity model in which
consumer species compete for an abiotic or biotic resource. We consider both
consumer species specialized to a habitat patch, and generalist species capable
of using the resource throughout the metacommunity. We present analytical
results for different limiting values of consumer dispersal and resource
dispersal, and complement these results with simulations for intermediate
dispersal values. Our analysis reveals generic patterns for the combined
effects of consumer and resource dispersal on the metacommunity diversity of
consumer species, and shows that hump-shaped relationships between local
diversity and dispersal are not universal. Diversity-dispersal relationships
can also be monotonically increasing or multimodal. Our work is a new step
towards a general theory of metacommunity diversity integrating dispersal at
multiple trophic levels.


  Reconstructing past population size from present day genetic data is a major
goal of population genetics. Recent empirical studies infer population size
history using coalescent-based models applied to a small number of individuals.
Here we provide tight bounds on the amount of exact coalescence time data
needed to recover the population size history of a single, panmictic population
at a certain level of accuracy. In practice, coalescence times are estimated
from sequence data and so our lower bounds should be taken as rather
conservative.


  We analyse a model consisting of a population of individuals which is
subdivided into a finite set of demes, each of which has a fixed but differing
number of individuals. The individuals can reproduce, die and migrate between
the demes according to an arbitrary migration network. They are haploid, with
two alleles present in the population; frequency independent selection is also
incorporated, where the strength and direction of selection can vary from deme
to deme. The system is formulated as an individual-based model, and the
diffusion approximation systematically applied to express it as a set of
nonlinear coupled stochastic differential equations. These can be made amenable
to analysis through the elimination of fast-time variables. The resulting
reduced model is analysed in a number of situations, including
migration-selection balance leading to a polymorphic equilibrium of the two
alleles, and an illustration of how the subdivision of the population can lead
to non-trivial behaviour in the case where the network is a simple hub. The
method we develop is systematic, may be applied to any network, and agrees well
with the results of simulations in all cases studied and across a wide range of
parameter values.


  We show that evolutionary computation can be implemented as standard
Markov-chain Monte-Carlo (MCMC) sampling. With some care, `genetic algorithms'
can be constructed that are reversible Markov chains that satisfy detailed
balance; it follows that the stationary distribution of populations is a Gibbs
distribution in a simple factorised form. For some standard and popular
nonparametric probability models, we exhibit Gibbs-sampling procedures that are
plausible genetic algorithms. At mutation-selection equilibrium, a population
of genomes is analogous to a sample from a Bayesian posterior, and the genomes
are analogous to latent variables. We suggest this is a general, tractable, and
insightful formulation of evolutionary computation in terms of standard machine
learning concepts and techniques.
  In addition, we show that evolutionary processes in which selection acts by
differences in fecundity are not reversible, and also that it is not possible
to construct reversible evolutionary models in which each child is produced by
only two parents.


  A looming challenge for agriculture is sustainable intensification of food
production to feed the growing human population. Current chemical and genetic
technologies used to manage plant diseases are highly vulnerable to pathogen
evolution and are not sustainable. Pathogen evolution is facilitated by the
genetic uniformity underlying modern agroecosystems, suggesting that one path
to sustainable disease control lies through increasing genetic diversity at the
field scale by using genetically diverse host mixtures. We investigate how host
mixtures can improve disease control using a population dynamics model. We find
that when a population of crop plants is exposed to host-specialized pathogen
species or strains, the overall disease severity is smaller in the mixture of
two host varieties than in each of the corresponding pure stands. The disease
severity can be minimized over a range of mixing ratios. These findings may
help in designing host mixtures that efficiently control diseases of crops. We
then generalize the model to describe host mixtures with many components. We
find that when pathogens exhibit host specialization, the overall disease
severity decreases with the number of components in the mixture. As the degree
of specialization increases, the decrease in disease severity becomes larger.
Using these model outcomes, we propose ways to optimize the use of host
mixtures to decrease disease in agroecosystems.


  Models for genome-wide prediction and association studies usually target a
single phenotypic trait. However, in animal and plant genetics it is common to
record information on multiple phenotypes for each individual that will be
genotyped. Modeling traits individually disregards the fact that they are most
likely associated due to pleiotropy and shared biological basis, thus providing
only a partial, confounded view of genetic effects and phenotypic interactions.
In this paper we use data from a Multiparent Advanced Generation Inter-Cross
(MAGIC) winter wheat population to explore Bayesian networks as a convenient
and interpretable framework for the simultaneous modeling of multiple
quantitative traits. We show that they are equivalent to multivariate genetic
best linear unbiased prediction (GBLUP), and that they are competitive with
single-trait elastic net and single-trait GBLUP in predictive performance.
Finally, we discuss their relationship with other additive-effects models and
their advantages in inference and interpretation. MAGIC populations provide an
ideal setting for this kind of investigation because the very low population
structure and large sample size result in predictive models with good power and
limited confounding due to relatedness.


  Much of the current theory of adaptation is based on Gillespie's mutational
landscape model (MLM), which assumes that the fitness values of genotypes
linked by single mutational steps are independent random variables. On the
other hand, a growing body of empirical evidence shows that real fitness
landscapes, while possessing a considerable amount of ruggedness, are smoother
than predicted by the MLM. In the present article we propose and analyse a
simple fitness landscape model with tunable ruggedness based on the Rough Mount
Fuji (RMF) model originally introduced by Aita et al. [Biopolymers 54:64-79
(2000)] in the context of protein evolution. We provide a comprehensive
collection of results pertaining to the topographical structure of RMF
landscapes, including explicit formulae for the expected number of local
fitness maxima, the location of the global peak, and the fitness correlation
function. The statistics of single and multiple adaptive steps on the RMF
landscape are explored mainly through simulations, and the results are compared
to the known behavior in the MLM model. Finally, we show that the RMF model can
explain the large number of second-step mutations observed on a highly-fit
first step backgound in a recent evolution experiment with a microvirid
bacteriophage [Miller et al., Genetics 187:185-202 (2011)].


  The outcomes of evolution are determined by which mutations occur and fix. In
rapidly adapting microbial populations, this process is particularly hard to
predict because lineages with different beneficial mutations often spread
simultaneously and interfere with one another's fixation. Hence to predict the
fate of any individual variant, we must know the rate at which new mutations
create competing lineages of higher fitness. Here, we directly measured the
effect of this interference on the fates of specific adaptive variants in
laboratory Saccharomyces cerevisiae populations and used these measurements to
infer the distribution of fitness effects of new beneficial mutations. To do
so, we seeded marked lineages with different fitness advantages into replicate
populations and tracked their subsequent frequencies for hundreds of
generations. Our results illustrate the transition between strongly
advantageous lineages which decisively sweep to fixation and more moderately
advantageous lineages that are often outcompeted by new mutations arising
during the course of the experiment. We developed an approximate likelihood
framework to compare our data to simulations and found that the effects of
these competing beneficial mutations were best approximated by an exponential
distribution, rather than one with a single effect size. We then used this
inferred distribution of fitness effects to predict the rate of adaptation in a
set of independent control populations. Finally, we discuss how our
experimental design can serve as a screen for rare, large-effect beneficial
mutations.


  The effectiveness and validity of applying variation partitioning methods in
community ecology has been questioned. Here, using mathematical deduction and
numerical simulation, we made an attempt to uncover the underlying mechanisms
determining the effectiveness of variation partitioning techniques. The
covariance among independent variables determines the under-fitting and
over-fitting problem with the variation partitioning process. Ideally, it is
assumed that the covariance among independent variables will be zero (no
correlation at all), however, typically there will be some colinearities.
Therefore, we analyzed the role of slight covariance on influencing species
variation partitioning. We concluded that when the covariance between spatial
and environmental predictors is positive, all the three components-pure
environmental, spatial variations and mixed covariation were over-fitted, with
the sign of the true covariation being negative. In contrast, when the
covariance is negative, all the three components were under-fitted with the
sign of true covariation being positive. Other factors, including extra noise
levels, the strengths of variable coefficients and the patterns of landscape
gradients, could reduce the fitting problems caused by the covariance of
variables. The conventional calculation of mixed covariation is incorrect and
misleading, as the true and estimated covariations are always sign-opposite. In
conclusion, I challenge the conventional three-step procedure of variation
partitioning, suggesting that a full regression model with all variables
together is robust enough to correctly partition variations.


  Polymers exposed to shear flow exhibit a rich tumbling dynamics. While rigid
rods rotate on Jeffery orbits, flexible polymers stretch and coil up during
tumbling. Theoretical results show that in both of these asymptotic regimes the
tumbling frequency f_c in a linear shear flow of strength \gamma scales as a
power law Wi^(2/3) in the Weissenberg number Wi=\gamma \tau, where \tau is a
characteristic time of the polymer's relaxational dynamics. For flexible
polymers these theoretical results are well confirmed by experimental single
molecule studies. However, for the intermediate semiflexible regime the
situation is less clear. Here we perform extensive Brownian dynamics
simulations to explore the tumbling dynamics of semiflexible polymers over a
broad range of shear strength and the polymer's persistence length l_p. We find
that the Weissenberg number alone does not suffice to fully characterize the
tumbling dynamics, and the classical scaling law breaks down. Instead, both the
polymer's stiffness and the shear rate are relevant control parameters. Based
on our Brownian dynamics simulations we postulate that in the parameter range
most relevant for cytoskeletal filaments there is a distinct scaling behavior
with f_c \tau*=Wi^(3/4) f_c (x) with Wi=\gamma \tau* and the scaling variable
x=(l_p/L)(Wi)^(-1/3); here \tau* is the time the polymer's center of mass
requires to diffuse its own contour length L. Comparing these results with
experimental data on F-actin we find that the Wi^(3/4) scaling law agrees
quantitatively significantly better with the data than the classical Wi^(2/3)
law. Finally, we extend our results to single ring polymers in shear flow, and
find similar results as for linear polymers with slightly different power laws.


  The probability distribution describing the state of a Stochastic Reaction
Network evolves according to the Chemical Master Equation (CME). It is common
to estimated its solution using Monte Carlo methods such as the Stochastic
Simulation Algorithm (SSA). In many cases these simulations can take an
impractical amount of computational time. Therefore many methods have been
developed that approximate the Stochastic Process underlying the Chemical
Master Equation. Prominent strategies are Hybrid Models that regard the firing
of some reaction channels as being continuous and applying the quasi-stationary
assumption to approximate the dynamics of fast subnetworks. However as the
dynamics of a Stochastic Reaction Network changes with time these
approximations might have to be adapted during the simulation. We develop a
method that approximates the solution of a CME by automatically partitioning
the reaction dynamics into discrete/continuous components and applying the
quasi-stationary assumption on identifiable fast subnetworks. Our method does
not require user intervention and it adapts to exploit the changing timescale
separation between reactions and/or changing magnitudes of copy numbers of
constituent species. We demonstrate the efficiency of the proposed method by
considering examples from Systems Biology and showing that very good
approximations to the exact probability distributions can be achieved in
significantly less computational time.


  Recently, Press and Dyson have proposed a new class of probabilistic and
conditional strategies for the two-player iterated Prisoner's Dilemma,
so-called zero-determinant strategies. A player adopting zero-determinant
strategies is able to pin the expected payoff of the opponents or to enforce a
linear relationship between his own payoff and the opponents' payoff, in a
unilateral way. This paper considers zero-determinant strategies in the
iterated public goods game, a representative multi-player evolutionary game
where in each round each player will choose whether or not put his tokens into
a public pot, and the tokens in this pot are multiplied by a factor larger than
one and then evenly divided among all players. The analytical and numerical
results exhibit a similar yet different scenario to the case of two-player
games: (i) with small number of players or a small multiplication factor, a
player is able to unilaterally pin the expected total payoff of all other
players; (ii) a player is able to set the ratio between his payoff and the
total payoff of all other players, but this ratio is limited by an upper bound
if the multiplication factor exceeds a threshold that depends on the number of
players.


  Motivation: Molecular signatures for predicting breast cancer prognosis could
greatly improve care through personalization of treatment. Computational
analyses of genome-wide expression datasets have identified such signatures,
but these signatures leave much to be desired in terms of accuracy,
reproducibility and biological interpretability. Methods that take advantage of
structured prior knowledge (e.g. protein interaction networks) show promise in
helping to define better signatures but most knowledge remains unstructured.
  Crowdsourcing via scientific discovery games is an emerging methodology that
has the potential to tap into human intelligence at scales and in modes
previously unheard of. Here, we developed and evaluated a game called The Cure
on the task of gene selection for breast cancer survival prediction. Our
central hypothesis was that knowledge linking expression patterns of specific
genes to breast cancer outcomes could be captured from game players. We
envisioned capturing knowledge both from the players prior experience and from
their ability to interpret text related to candidate genes presented to them in
the context of the game.
  Results: Between its launch in Sept. 2012 and Sept. 2013, The Cure attracted
more than 1,000 registered players who collectively played nearly 10,000 games.
Gene sets assembled through aggregation of the collected data clearly
demonstrated the accumulation of relevant expert knowledge. In terms of
predictive accuracy, these gene sets provided comparable performance to gene
sets generated using other methods including those used in commercial tests.
The Cure is available at http://genegames.org/cure/


  This paper proposes a non-phenomenological model of population growth that is
based on the interactions between the individuals that compose the system. It
is assumed that the individuals interact cooperatively and competitively. As a
consequence of this interaction, it is shown that some well-known
phenomenological population growth models (such as the Malthus, Verhulst,
Gompertz, Richards, Von Foerster, and power-law growth models) are special
cases of the model presented herein. Moreover, other ecological behaviors can
be seen as the emergent behavior of such interactions. For instance, the Allee
effect, which is the characteristic of some populations to increase the
population growth rate at a small population size, is observed. Whereas the
models presented in the literature explain the Allee effect with
phenomenological ideas, the model presented here explains this effect by the
interactions between the individuals. The model is tested with empirical data
to justify its formulation. Other interesting macroscopic emergent behavior
from the model proposed here is the observation of a regime of population
divergence at a finite time. It is interesting that this characteristic is
observed in humanity's global population growth. It is shown that in a regime
of cooperation, the model fits very well to the human population growth data
since 1000 A.D.


  The compartmentalization of distinct templates in protocells and the exchange
of templates between them (migration) are key elements of a modern scenario for
prebiotic evolution. Here we use the diffusion approximation of population
genetics to study analytically the steady-state properties of such prebiotic
scenario. The coexistence of distinct template types inside a protocell is
achieved by a selective pressure at the protocell level (group selection)
favoring protocells with a mixed template composition. In the degenerate case,
where the templates have the same replication rate, we find that a vanishingly
small migration rate suffices to eliminate the segregation effect of random
drift and so to promote coexistence. In the non-degenerate case, a small
migration rate greatly boosts coexistence as compared with the situation where
there is no migration. However, increase of the migration rate beyond a
critical value leads to the complete dominance of the more efficient template
type (homogeneous regime). In this case, we find a continuous phase transition
separating the homogeneous and the coexistence regimes, with the order
parameter vanishing linearly with the distance to the transition point.


  The idea that all life on earth traces back to a common beginning dates back
at least to Charles Darwin's {\em Origin of Species}. Ever since, biologists
have tried to piece together parts of this `tree of life' based on what we can
observe today: fossils, and the evolutionary signal that is present in the
genomes and phenotypes of different organisms. Mathematics has played a key
role in helping transform genetic data into phylogenetic (evolutionary) trees
and networks. Here, I will explain some of the central concepts and basic
results in phylogenetics, which benefit from several branches of mathematics,
including combinatorics, probability and algebra.


  Subunit epsilon is an intrinsic regulator of the bacterial FoF1-ATP synthase,
the ubiquitous membrane-embedded enzyme that utilizes a proton motive force in
most organisms to synthesize adenosine triphosphate (ATP). The C-terminal
domain of epsilon can extend into the central cavity formed by the alpha and
beta subunits, as revealed by the recent X-ray structure of the F1 portion of
the Escherichia coli enzyme. This insertion blocks the rotation of the central
gamma subunit and, thereby, prevents wasteful ATP hydrolysis. Here we aim to
develop an experimental system that can reveal conditions under which epsilon
inhibits the holoenzyme FoF1-ATP synthase in vitro. Labeling the C-terminal
domain of epsilon and the gamma subunit specifically with two different
fluorophores for single-molecule Foerster resonance energy transfer (smFRET)
allowed monitoring of the conformation of epsilon in the reconstituted enzyme
in real time. New mutants were made for future three-color smFRET experiments
to unravel the details of regulatory conformational changes in epsilon.


  This paper presents and tests a previously unrecognised mechanism for driving
a replicating molecular system on the prebiotic earth. It is proposed that
cell-free RNA replication in the primordial soup may have been driven by
self-sustained oscillatory thermochemical reactions. To test this hypothesis a
well-characterised hydrogen peroxide oscillator was chosen as the driver and
complementary RNA strands with known association and melting kinetics were used
as the substrate. An open flow system model for the self-consistent, coupled
evolution of the temperature and concentrations in a simple autocatalytic
scheme is solved numerically, and it is shown that thermochemical cycling
drives replication of the RNA strands. For the (justifiably realistic) values
of parameters chosen for the simulated example system, the mean amount of
replicant produced at steady state is 6.56 times the input amount, given a
constant supply of substrate species. The spontaneous onset of sustained
thermochemical oscillations via slowly drifting parameters is demonstrated, and
a scheme is given for prebiotic production of complementary RNA strands on rock
surfaces.


  We study theoretically in the present work the self-assembly of molecules in
an open system, which is fed by monomers and depleted in partial or complete
clusters. Such a scenario is likely to occur for example in the context of
viral self-assembly. We provide a general formula for the mean-field size
distribution which is valid both at equilibrium in a closed system, and in the
stationary state in an open system. This allows us to explore in a simple way
out-of-equilibrium features for self-assembly and compare them to equilibrium
properties. In particular, we identify a region of parameter space for which
the out-of-equilibrium size distribution in the presence of external fluxes is
equal to the equilibrium size distribution in the absence of external fluxes,
up to a constant renormalization factor. The range of validity of this result
and its consequences are discussed.


  HIV-1, an enveloped RNA virus, produces viral particles that are known to be
much more heterogeneous in size than is typical of non-enveloped viruses. We
present here a novel strategy to study HIV-1 Viral Like Particles (VLP)
assembly by measuring the size distribution of these purified VLPs and
subsequent viral cores thanks to Atomic Force Microscopy imaging and
statistical analysis. This strategy allowed us to identify whether the presence
of viral RNA acts as a modulator for VLPs and cores size heterogeneity in a
large population of particles. These results are analyzed in the light of a
recently proposed statistical physics model for the self-assembly process. In
particular, our results reveal that the modulation of size distribution by the
presence of viral RNA is qualitatively reproduced, suggesting therefore an
entropic origin for the modulation of RNA uptake by the nascent VLP.


  In small volumes, the kinetics of filamentous protein self-assembly is
expected to show significant variability, arising from intrinsic molecular
noise. This is not accounted for in existing deterministic models. We introduce
a simple stochastic model including nucleation and autocatalytic growth via
elongation and fragmentation, which allows us to predict the effects of
molecular noise on the kinetics of autocatalytic self-assembly. We derive an
analytic expression for the lag-time distribution, which agrees well with
experimental results for the fibrillation of bovine insulin. Our expression
decomposes the lag time variability into contributions from primary nucleation
and autocatalytic growth and reveals how each of these scales with the key
kinetic parameters. Our analysis shows that significant lag-time variability
can arise from both primary nucleation and from autocatalytic growth, and
should provide a way to extract mechanistic information on early-stage
aggregation from small-volume experiments.


  Microtubules (MTs) nucleated by centrosomes form star-shaped structures
referred to as asters. Aster motility and dynamics is vital for genome
stability, cell division, polarization and differentiation. Asters move either
towards the cell center or away from it. Here, we focus on the centering
mechanism in a membrane independent system of Xenopus cytoplasmic egg extracts.
Using live microscopy and single particle tracking, we find that asters move
towards chromatinized DNA structures. The velocity and directionality profiles
suggest a random walk with drift directed towards DNA. We have developed a
theoretical model that can explain this movement as a result of a gradient of
MT length dynamics and MT gliding on immobilized dynein motors. In simulations,
the antagonistic action of the motor species on the radial array of MTs leads
to a tug-of-war purely due to geometric considerations and aster motility
resembles a directed random-walk. Additionally our model predicts that aster
velocities do not change greatly with varying initial distance from DNA. The
movement of asymmetric asters becomes increasingly super-diffusive with
increasing motor density, but for symmetric asters it becomes less
super-diffusive. The transition of symmetric asters from superdiffusive to
diffusive mobility is the result of number fluctuations in bound motors in the
tug-of-war. Overall, our model is in good agreement with experimental data in
Xenopus cytoplasmic extracts and predicts novel features of the collective
effects of motor-MT interactions.


  During M phase, Endosulfine (Endos) family proteins are phosphorylated by
Greatwall kinase (Gwl), and the resultant pEndos inhibits the phosphatase
PP2A-B55, which would otherwise prematurely reverse many CDK-driven
phosphorylations. We show here that PP2A-B55 is the enzyme responsible for
dephosphorylating pEndos during M phase exit. The kinetic parameters for
PP2A-B55's action on pEndos are orders of magnitude lower than those for
CDK-phosphorylated substrates, suggesting a simple model for PP2A-B55
regulation that we call inhibition by unfair competition. As the name suggests,
during M phase PP2A-B55's attention is diverted to pEndos, which binds much
more avidly and is dephosphorylated more slowly than other substrates. When Gwl
is inactivated during the M phase-to-interphase transition, the dynamic balance
changes: pEndos dephosphorylated by PP2A-B55 cannot be replaced, so the
phosphatase can refocus its attention on CDK-phosphorylated substrates. This
mechanism explains simultaneously how PP2A-B55 and Gwl together regulate
pEndos, and how pEndos controls PP2A-B55.


  The nature of the quantitative distribution of the 64 DNA codons in the human
genome has been an issue of debate for over a decade. Some groups have proposed
that the quantitative distribution of the DNA codons ordered as a
rank-frequency plot follows a well-known power law called Zipf's law. Others
have shown that the DNA codon distribution is best fitted to an exponential
function. However, the reason for such scaling behavior has not yet been
addressed. In the present study, we demonstrate that the nonlinearity of the
DNA codon distribution is a direct consequence of the frequency recurrence of
the codon usage (i.e., the repetitiveness of codon usage frequencies at the
whole genome level). We discover that if frequency recurrence is absent from
the human genome, the frequency of occurrence of codons scales linearly with
the codon rank. We also show that DNA codons of both low and high frequency of
occurrence in the genome are best fitted by an exponential function and provide
strong evidence to suggest that the coding region of the human genome does not
follow Zipf's law. Information-theoretic methods and entropy calculations are
applied to the DNA codon distribution and a new approach, called the lariat
method, is proposed to quantitatively analyze the DNA codon distribution in
Homo sapiens.


  We study the existence and nonexistence of traveling waves of general
diffusive Kermack-McKendrick SIR models with standard incidence where the total
population is not constant. The three classes, susceptible $S$, infected $I$
and removed $R$, are all involved in the traveling wave solutions. We show that
the minimum speed for the existence of traveling waves for this
three-dimensional non-monotonic system can be derived from its linearizaion at
the initial disease-free equilibrium. The proof in this paper is based on
Schauder fixed point theorem and Laplace transform and provides a promising
method to deal with high dimensional epidemic models.


  We have studied the apo (Fe3+ free) form of periplasmic ferric binding
protein (FbpA) under different conditions and we have monitored the changes in
the binding and release dynamics of H2PO4- that acts as a synergistic anion in
the presence of Fe3+. Our simulations predict a dissociation constant of
2.2$\pm$0.2 mM which is in remarkable agreement with the experimentally
measured value of 2.3$\pm$0.3 mM under the same ionization strength and pH
conditions. We apply perturbations relevant for changes in environmental
conditions as (i) different values of ionic strength (IS), and (ii) protonation
of a group of residues to mimic a different pH environment. Local perturbations
are also studied by protonation or mutation of a site distal to the binding
region that is known to mechanically manipulate the hinge-like motions of FbpA.
We find that while the average conformation of the protein is intact in all
simulations, the H2PO4- dynamics may be substantially altered by the changing
conditions. In particular, the bound fraction which is 20$\%$ for the wild type
system is increased to 50$\%$ with a D52A mutation/protonation and further to
over 90$\%$ at the protonation conditions mimicking those at pH 5.5. The change
in the dynamics is traced to the altered electrostatic distribution on the
surface of the protein which in turn affects hydrogen bonding patterns at the
active site. The observations are quantified by rigorous free energy
calculations. Our results lend clues as to how the environment versus single
residue perturbations may be utilized for regulation of binding modes in hFbpA
systems in the absence of conformational changes.


  Living at high-altitude is one of the most difficult challenges that humans
had to cope with during their evolution. Whereas several genomic studies have
revealed some of the genetic bases of adaptations in Tibetan, Andean and
Ethiopian populations, relatively little evidence of convergent evolution to
altitude in different continents has accumulated. This lack of evidence can be
due to truly different evolutionary responses, but it can be also due to the
low power of former studies that have mainly focused on populations from a
single geographical region or performed separate analyses on multiple pairs of
populations to avoid problems linked to shared histories between some
populations. We introduce here a hierarchical Bayesian method to detect local
adaptation that can deal with complex demographic histories. Our method can
identify selection occurring at different scales, as well as convergent
adaptation in different regions. We apply our approach to the analysis of a
large SNP dataset from low- and high-altitude human populations from America
and Asia. The simultaneous analysis of these two geographic areas allows us to
identify several candidate genome regions for altitudinal selection, and we
show that convergent evolution among continents has been quite common. In
addition to identifying several genes and biological processes involved in high
altitude adaptation, we identify two specific biological pathways that could
have evolved in both continents to counter toxic effects induced by hypoxia.


  Current modeling of infectious diseases allows for the study of complex and
realistic scenarios that go from the population to the individual level of
description. However, most epidemic models assume that the spreading process
takes place on a single level (be it a single population, a meta-population
system or a network of contacts). In particular, interdependent contagion
phenomena can only be addressed if we go beyond the scheme one pathogen-one
network. In this paper, we propose a framework that allows describing the
spreading dynamics of two concurrent diseases. Specifically, we characterize
analytically the epidemic thresholds of the two diseases for different
scenarios and also compute the temporal evolution characterizing the unfolding
dynamics. Results show that there are regions of the parameter space in which
the onset of a disease's outbreak is conditioned to the prevalence levels of
the other disease. Moreover, we show, for the SIS scheme, that under certain
circumstances, finite and not vanishing epidemic thresholds are found even at
the thermodynamic limit for scale-free networks. For the SIR scenario, the
phenomenology is richer and additional interdependencies show up. We also find
that the secondary thresholds for the SIS and SIR models are different, which
results directly from the interaction between both diseases. Our work thus
solve an important problem and pave the way towards a more comprehensive
description of the dynamics of interacting diseases.


  The role of skill (fitness) and luck (randomness) as driving forces on the
dynamics of resource accumulation in a myriad of systems have long puzzled
scientists. Fueled by undisputed inequalities that emerge from actual
competitions, there is a pressing need for better understanding the effects of
skill and luck in resource accumulation. When such competitions are driven by
externalities such as cumulative advantage (CA), the rich-get-richer effect,
little is known with respect to fundamental properties such as their duration
and intensity. In this work we provide a mathematical understanding of how CA
exacerbates the role of luck in detriment of skill in simple and well-studied
competition models. We show, for instance, that if two agents are competing for
resources that arrive sequentially at each time unit, an early stroke of luck
can place the less skilled in the lead for an extremely long period of time, a
phenomenon we call "struggle of the fittest". In the absence of CA, the more
skilled quickly prevails despite any early stroke of luck that the less skilled
may have. We prove that duration of a simple skill and luck competition model
exhibit power law tails when CA is present, regardless of skill difference,
which is in sharp contrast to exponential tails when CA is absent. Our findings
have important implications to competitions not only in complex social systems
but also in contexts that leverage such models.


  Expression quantitative trait loci (eQTL) mapping constitutes a challenging
problem due to, among other reasons, the high-dimensional multivariate nature
of gene-expression traits. Next to the expression heterogeneity produced by
confounding factors and other sources of unwanted variation, indirect effects
spread throughout genes as a result of genetic, molecular and environmental
perturbations. From a multivariate perspective one would like to adjust for the
effect of all of these factors to end up with a network of direct associations
connecting the path from genotype to phenotype. In this paper we approach this
challenge with mixed graphical Markov models, higher-order conditional
independences and q-order correlation graphs. These models show that additive
genetic effects propagate through the network as function of gene-gene
correlations. Our estimation of the eQTL network underlying a well-studied
yeast data set leads to a sparse structure with more direct genetic and
regulatory associations that enable a straightforward comparison of the genetic
control of gene expression across chromosomes. Interestingly, it also reveals
that eQTLs explain most of the expression variability of network hub genes.


  The uniform sampling of convex regions in high dimension is an important
computational issue, from both theoretical and applied point of view. The
hit-and-run montecarlo algorithms are the most efficient methods known to
perform it and one of their bottlenecks relies in the difficulty of escaping
from tight corners in high dimension. Inspired by optimized montecarlo methods
used in statistical mechanics we define a new algorithm by overelaxing the
hit-and-run dynamics. We made numerical simulations on high dimensional
simplices in order to test its performances, pointing out its improved ability
to escape from angles and finally apply it to an inference problem in the
steady state dynamics of metabolic networks.


  The Ultimatum Game (UG) is an economic game where two players (proposer and
responder) decide how to split a certain amount of money. While traditional
economic theories based on rational decision making predict that the proposer
should make a minimal offer and the responder should accept it, human subjects
tend to behave more fairly in UG. Previous studies suggested that extra
information such as reputation, empathy, or spatial structure is needed for
fairness to evolve in UG. Here we show that fairness can evolve without
additional information if players make decisions probabilistically and may
continue interactions when the offer is rejected, which we call the Not Quite
Ultimatum Game (NQUG). Evolutionary simulations of NQUG showed that the
probabilistic decision making contributes to the increase of proposers' offer
amounts to avoid rejection, while the repetition of the game works to
responders' advantage because they can wait until a good offer comes. These
simple extensions greatly promote evolution of fairness in both proposers'
offers and responders' acceptance thresholds.


  Recent research has identified the heterogeneity as crucial for the evolution
of cooperation in spatial population. However, the influence of heterogeneous
noise is still lack. Inspired by this interesting question, in this work, we
try to incorporate heterogeneous noise into the evaluation of utility, where
only a proportion of population possesses noise, whose range can also be tuned.
We find that increasing heterogeneous noise monotonously promotes cooperation
and even translates the full defection phase (of the homogeneous version) into
the complete cooperation phase. Moreover, the promotion effect of this
mechanism can be attributed to the leading role of cooperators who have the
heterogeneous noise. These type of cooperators can attract more agents
penetrating into the robust cooperator clusters, which is beyond the text of
traditional spatial reciprocity. We hope that our work may shed light on the
understanding of the cooperative behavior in the society.


  Background: Short-read aligners have recently gained a lot of speed by
exploiting the massive parallelism of GPU. An uprising alternative to GPU is
Intel MIC; supercomputers like Tianhe-2, currently top of TOP500, is built with
48,000 MIC boards to offer ~55 PFLOPS. The CPU-like architecture of MIC allows
CPU-based software to be parallelized easily; however, the performance is often
inferior to GPU counterparts as an MIC board contains only ~60 cores (while a
GPU board typically has over a thousand cores). Results: To better utilize
MIC-enabled computers for NGS data analysis, we developed a new short-read
aligner MICA that is optimized in view of MICs limitation and the extra
parallelism inside each MIC core. Experiments on aligning 150bp paired-end
reads show that MICA using one MIC board is 4.9 times faster than the BWA-MEM
(using 6-core of a top-end CPU), and slightly faster than SOAP3-dp (using a
GPU). Furthermore, MICAs simplicity allows very efficient scale-up when
multiple MIC boards are used in a node (3 cards give a 14.1-fold speedup over
BWA-MEM). Summary: MICA can be readily used by MIC-enabled supercomputers for
production purpose. We have tested MICA on Tianhe-2 with 90 WGS samples (17.47
Tera-bases), which can be aligned in an hour less than 400 nodes. MICA has
impressive performance even though the current MIC is at its initial stage of
development (the next generation of MIC has been announced to release in late
2014).


  Non-African populations have experienced major bottlenecks in the time since
their split from West Africans, which has led to the hypothesis that natural
selection to remove weakly deleterious mutations may have been less effective
in non-Africans. To directly test this hypothesis, we measure the per-genome
accumulation of deleterious mutations across diverse humans. We fail to detect
any significant differences, but find that archaic Denisovans accumulated
non-synonymous mutations at a higher rate than modern humans, consistent with
the longer separation time of modern and archaic humans. We also revisit the
empirical patterns that have been interpreted as evidence for less effective
removal of deleterious mutations in non-Africans than in West Africans, and
show they are not driven by differences in selection after population
separation, but by neutral evolution.


  Gliding of microtubule filaments on surfaces coated with the motor protein
kinesin has potential applications for nano-scale devices. The ability to guide
the gliding direction in three dimensions allows the fabrication of tracks of
arbitrary geometry in space. Here, we achieve this by using kinesin-coated
glass wires of micrometer diameter range. Unlike previous methods in which the
guiding tracks are fixed on flat two-dimensional surfaces, the flexibility of
glass wires in shape and size facilitates building in-vitro devices that have
deformable tracks.


  Evolution of cooperation is a widely studied problem in biology, social
science, economics, and artificial intelligence. Most of the existing
approaches that explain cooperation rely on some notion of direct or indirect
reciprocity. These reciprocity based models assume agents recognize their
partner and know their previous interactions, which requires advanced cognitive
abilities. In this paper we are interested in developing a model that produces
cooperation without requiring any explicit memory of previous game plays. Our
model is based on the notion of, a concept introduced within behavioral
economics, whereby individuals care about payoff equality in outcomes. Here we
explore the effect of using income inequality to guide partner selection and
interaction. We study our model by considering both the well-mixed and the
spatially structured population and present the conditions under which
cooperation becomes dominant. Our results support the hypothesis that inequity
aversion promotes cooperative relationship among nonkin.


  In many key applications of metabolomics, such as toxicology or
nutrigenomics, it is of interest to profile and detect changes in metabolic
processes, usually represented in the form of pathways. As an alternative, a
broader point of view would enable investigators to better understand the
relations between entities that exist in different processes. Therefore,
relating a possible perturbation to several known processes represents a new
approach to this field of study. We propose to use a network representation of
metabolism in terms of reactants, enzyme and metabolite. To model these systems
it is possible to describe both reactions and relations among enzymes and
metabolites. In this way, analysis of the impact of changes in some metabolites
or enzymes on different processes are easier to understand, detect and predict.
Results: We release the MetaboX library, an open source PHP framework for
developing metabolic networks from a set of compounds. This library uses data
stored in Kyoto Encyclopedia for Genes and Genomes (KEGG) database using its
RESTful Application Programming Interfaces (APIs), and methods to enhance
manipulation of the information returned from KEGG webservice. The MetaboX
library includes methods to extract information about a resource of interest
(e.g. metabolite, reaction and enzyme) and to build reactants networks,
bipartite enzyme-metabolite and unipartite enzyme networks. These networks can
be exported in different formats for data visualization with standard tools. As
a case study, the networks built from a subset of the Glycolysis pathway are
described and discussed. Conclusions: The advantages of using such a library
imply the ability to model complex systems with few starting information
represented by a collection of metabolites KEGG IDs.


  We explore the evolutionary dynamics of two games - the Prisoner's Dilemma
and the Snowdrift Game - played within distinct networks (layers) of
interdependent networks. In these networks imitation and interaction between
individuals of opposite layers is established through interlinks. We explore an
update rule in which revision of strategies is a biased imitation process:
individuals imitate neighbors from the same layer with probability p, and
neighbors from the second layer with complementary probability 1 - p. We
demonstrate that a small decrease of p from p = 1 (which corresponds to
forbidding strategy transfer between layers) is sufficient to promote
cooperation in the Prisoner's Dilemma subpopulation. This, on the other hand,
is detrimental for cooperation in the Snowdrift Game subpopulation. We provide
results of extensive computer simulations for the case in which layers are
modelled as regular random networks, and support this study with analytical
results for coupled well-mixed populations.


  Search for possible relationships between phylogeny and ontogeny is one of
the most important issues in the field of evolutionary developmental biology.
By representing developmental dynamics of spatially located cells with gene
expression dynamics with cell-to-cell interaction under external morphogen
gradient, evolved are gene regulation networks under mutation and selection
with the fitness to approach a prescribed spatial pattern of expressed genes.
For most of thousands of numerical evolution experiments, evolution of pattern
over generations and development of pattern by an evolved network exhibit
remarkable congruence. Here, both the pattern dynamics consist of several
epochs to form successive stripe formations between quasi-stationary regimes.
In evolution, the regimes are generations needed to hit relevant mutations,
while in development, they are due to the emergence of slowly varying
expression that controls the pattern change. Successive pattern changes are
thus generated, which are regulated by successive combinations of feedback or
feedforward regulations under the upstream feedforward network that reads the
morphogen gradient. By using a pattern generated by the upstream feedforward
network as a boundary condition, downstream networks form later stripe
patterns. These epochal changes in development and evolution are represented as
same bifurcations in dynamical-systems theory, and this agreement of
bifurcations lead to the evolution-development congruences. Violation of the
evolution-development congruence, observed exceptionally, is shown to be
originated in alteration of the boundary due to mutation at the upstream
feedforward network. Our results provide a new look on developmental stages,
punctuated equilibrium, developmental bottlenecks, and evolutionary acquisition
of novelty in morphogenesis.


  Studying strategy update rules in the framework of evolutionary game theory,
one can differentiate between imitation processes and aspiration-driven
dynamics. In the former case, individuals imitate the strategy of a more
successful peer. In the latter case, individuals adjust their strategies based
on a comparison of their payoffs from the evolutionary game to a value they
aspire, called the level of aspiration. Unlike imitation processes of pairwise
comparison, aspiration-driven updates do not require additional information
about the strategic environment and can thus be interpreted as being more
spontaneous. Recent work has mainly focused on understanding how aspiration
dynamics alter the evolutionary outcome in structured populations. However, the
baseline case for understanding strategy selection is the well-mixed population
case, which is still lacking sufficient understanding. We explore how
aspiration-driven strategy-update dynamics under imperfect rationality
influence the average abundance of a strategy in multi-player evolutionary
games with two strategies. We analytically derive a condition under which a
strategy is more abundant than the other in the weak selection limiting case.
This approach has a long standing history in evolutionary game and is mostly
applied for its mathematical approachability. Hence, we also explore strong
selection numerically, which shows that our weak selection condition is a
robust predictor of the average abundance of a strategy. The condition turns
out to differ from that of a wide class of imitation dynamics, as long as the
game is not dyadic. Therefore a strategy favored under imitation dynamics can
be disfavored under aspiration dynamics. This does not require any population
structure thus highlights the intrinsic difference between imitation and
aspiration dynamics.


  A central part of population genomics consists of finding genomic regions
implicated in local adaptation. Population genomic analyses are based on
genotyping numerous molecular markers and looking for outlier loci in terms of
patterns of genetic differentiation. One of the most common approach for
selection scan is based on statistics that measure population differentiation
such as $F_{ST}$. However they are important caveats with approaches related to
$F_{ST}$ because they require grouping individuals into populations and they
additionally assume a particular model of population structure. Here we
implement a more flexible individual-based approach based on Bayesian factor
models. Factor models capture population structure with latent variables called
factors, which can describe clustering of individuals into populations or
isolation-by-distance patterns. Using hierarchical Bayesian modeling, we both
infer population structure and identify outlier loci that are candidates for
local adaptation. As outlier loci, the hierarchical factor model searches for
loci that are atypically related to population structure as measured by the
latent factors. In a model of population divergence, we show that the factor
model can achieve a 2-fold or more reduction of false discovery rate compared
to the software BayeScan or compared to a $F_{ST}$ approach. We analyze the
data of the Human Genome Diversity Panel to provide an example of how factor
models can be used to detect local adaptation with a large number of SNPs. The
Bayesian factor model is implemented in the open-source PCAdapt software.


  Interactive molecular graphics applications facilitate analysis of three
dimensional protein structures. Naturally, non-interactive 2-D snapshots of the
protein structures do not convey the same level of geometric detail. Several
2-D visualization methods have been in use to summarize structural information,
including contact maps and 2-D cartoon views. We present a new approach for 2-D
visualization of protein structures where amino acid residues are displayed on
a circle and spatially close residues are depicted by links. Furthermore,
residue-specific properties, such as conservation, accessibility, temperature
factor, can be displayed as plots on the same circular view.


  The way for performing multiple sequence alignment is based on the criterion
of the maximum scored information content computed from a weight matrix, but it
is possible to have two or more alignments to have the same highest score
leading to ambiguities in selecting the best alignment. This paper addresses
this issue by introducing the concept of joint weight matrix to eliminate the
randomness in selecting the best multiple sequence alignment. Alignments with
equal scores are iteratively rescored with the joint weight matrix of
increasing level (nucleotide pairs, triplets and so on) until one single best
alignment is eventually found. This method for resolving ambiguity in multiple
sequence alignment can be easily implemented by use of the improved scoring
matrix.


  In this paper, a branching process model of the Northern Spotted Owl is
simulated. We focus on the time until extinction. It is shown how an
approximation of the model with a multivariate autoregressive process works
well near the equilibrium, but does not give a good estimate of the time until
extinction. We also show that introduction of randomness in some of the
parameters previously assumed to be constants shortens the time until
extinction considerably.


  A statistical study of cis-regulatory modules (CRMs) is presented based on
the estimation of similar-word set distribution. It is observed that CRMs tend
to have a fat-tail distribution. A new statistical fat-tail test with two
kurtosis-based fatness coefficients is proposed to distinguish CRMs from
non-CRMs. As compared with the existing fluffy-tail test, the first fatness
coefficient is designed to reduce computational time, making the novel fat-tail
test very suitable for long sequences and large database analysis in the
post-genome time and the second one to improve separation accuracy between CRMs
and non-CRMs. These two fatness coefficients may be served as valuable
filtering indexes to predict CRMs experimentally.


  The conventional way of identifying DNA motifs, solely based on match
alignment information, is susceptible to a high number of spurious sites. A
novel scoring system has been introduced by taking both match and mismatch
alignment information into account. The mismatch alignment information is
useful to remove spurious sites encountered in DNA motif searching. As an
example, a correct TATA box site in Homo sapiens H4/g gene has successfully
been identified based on match and mismatch alignment information.


  We study the impact on shape parameters of an underlying
Bienaym\'e-Galton-Watson branching process (height, width and first hitting
time), of having a non-spatial branching mechanism with infinite variance.
Aiming at providing a comparative study of the spread of an epidemics whose
dynamics is given by the modulus of a branching Brownian motion (BBM) we then
consider spatial branching processes in dimension d, not necessarily integer.
The underlying branching mechanism is then either a binary branching model or
one presenting infinite variance. In particular we evaluate the chance p(x) of
being hit if the epidemics started away at distance x. We compute the large x
tail probabilities of this event, both when the branching mechanism is regular
and when it exhibits very large fluctuations.


  Linear mixed models (LMMs) are a powerful and established tool for studying
genotype-phenotype relationships. A limiting assumption of LMMs is that the
residuals are Gaussian distributed, a requirement that rarely holds in
practice. Violations of this assumption can lead to false conclusions and
losses in power, and hence it is common practice to pre-process the phenotypic
values to make them Gaussian, for instance by applying logarithmic or other
non-linear transformations. Unfortunately, different phenotypes require
different specific transformations, and choosing a "good" transformation is in
general challenging and subjective. Here, we present an extension of the LMM
that estimates an optimal transformation from the observed data. In extensive
simulations and applications to real data from human, mouse and yeast we show
that using such optimal transformations lead to increased power in genome-wide
association studies and higher accuracy in heritability estimates and phenotype
predictions.


  We propose and numerically analyze a PDE model of ant foraging behavior. Ant
foraging is a prime example of individuals following simple behavioral rules
based on local information producing complex, organized and ``intelligent''
strategies at the population level. One of its main aspects is the widespread
use of pheromones, which are chemical compounds laid by the ants used to
attract other ants to a food source. In this work, we consider a continuous
description of a population of ants and simulate numerically the foraging
behavior using a system of PDEs of chemotaxis type. We show that, numerically,
this system accurately reproduces observed foraging behavior, such as trail
formation and efficient removal of food sources.


  The unique inheritance pattern of X chromosomes makes them preferential
targets of adaptive evolution. We here investigate natural selection on the X
chromosome in all species of great apes. We find that diversity is more
strongly reduced around genes on the X compared with autosomes, and that a
higher proportion of substitutions results from positive selection. Strikingly,
the X exhibits several megabase long regions where diversity is reduced more
than five fold. These regions overlap significantly among species, and have a
higher singleton proportion, population differentiation, and nonsynonymous to
synonymous substitution ratio. We rule out background selection and soft
selective sweeps as explanations for these observations, and conclude that
several strong selective sweeps have occurred independently in similar regions
in several species. Since these regions are strongly associated with ampliconic
sequences we propose that intra-genomic conflict between the X and the Y
chromosomes is a major driver of X chromosome evolution.


  Gene expression data, or transcription data, are surrogates for actual
protein concentrations in the cells. In addition protein-protein interactions
are static diagrams of all the protein-protein interactions in the cell. These
interactions may consist of covalent bonding or maybe just secondary bonding
such as hydrogen bonding. Given these two surrogate data types we show a
technique to compute the Gibbs free energy of a cell. We apply this to yeast
cell cycle and to cancer.


  Mutualistic interactions benefit both partners, promoting coexistence and
genetic diversity. Spatial structure can promote cooperation, but spatial
expansions may also make it hard for mutualistic partners to stay together,
since genetic drift at the expansion front creates regions of low genetic and
species diversity. To explore the antagonism between mutualism and genetic
drift, we grew cross-feeding strains of the budding yeast S. cerevisiae on agar
surfaces as a model for mutualists undergoing spatial expansions. By supplying
varying amounts of the exchanged nutrients, we tuned strength and symmetry of
the mutualistic interaction. Strong mutualism suppresses genetic demixing
during spatial expansions and thereby maintains diversity, but weak or
asymmetric mutualism is overwhelmed by genetic drift even when mutualism is
still beneficial, slowing growth and reducing diversity. Theoretical modeling
using experimentally measured parameters predicts the size of demixed regions
and how strong mutualism must be to survive a spatial expansion.


  Conjugation is the primary mechanism of horizontal gene transfer that spreads
antibiotic resistance among bacteria. Although conjugation normally occurs in
surface-associated growth (e.g., biofilms), it has been traditionally studied
in well-mixed liquid cultures lacking spatial structure, which is known to
affect many evolutionary and ecological processes. Here we visualize spatial
patterns of gene transfer mediated by F plasmid conjugation in a colony of
Escherichia coli growing on solid agar, and we develop a quantitative
understanding by spatial extension of traditional mass-action models. We found
that spatial structure suppresses conjugation in surface-associated growth
because strong genetic drift leads to spatial isolation of donor and recipient
cells, restricting conjugation to rare boundaries between donor and recipient
strains. These results suggest that ecological strategies, such as enforcement
of spatial structure and enhancement of genetic drift, could complement
molecular strategies in slowing the spread of antibiotic resistance genes.


  In contrast to proteins much less attention has been focused on development
of computational models for describing RNA molecules, which are being
recognized as playing key roles in many cellular functions. Current atomically
detailed force fields are not accurate enough to capture the properties of even
simple nucleic acid constructs. In this article, we review our efforts to
develop coarse-grained (CG) models that capture the underlying physics for the
particular length scale of interest. Two models are discussed. One of them is
the Three Interaction Site (TIS) model in which each nucleotide is represented
by three beads corresponding to sugar, phosphate, and base. The other is the
Self-Organized Polymer (SOP) model in which each nucleotide is represented as a
single interaction center. Applications of the TIS model to study the
complexity of hairpin formation and the effects of crowding in shifting
equilibrium between two conformations in human telomerase pseudoknot are
described. The work on crowding illustrates a direct link to the activity of
telomerase. We use the SOP model to describe the response of {\it Tetrahymena}
ribozyme to force. The simulated unfolding pathways agree well with single
molecule pulling experiments. We also review predictions for the unfolding
pathways for {\it Azoarcus} ribozyme. The success of the CG applications to
describe dynamics in RNA gives hope that more complex processes involving
RNA-protein interactions can be tackled using variants of the proposed models.


  We study the conditions under which species interaction, as described by
continuous versions of the competitive Lotka-Volterra model (namely the
nonlocal Kolmogorov-Fisher model, and its differential approximation), can
support the existence of localized states, i.e. patches of species with
enhanced population surrounded in niche space by species at smaller densities.
These states would arise from species interaction, and not by any preferred
niche location or better fitness. In contrast to previous works we include only
quadratic nonlinearities, so that the localized patches appear on a background
of homogeneously distributed species coexistence, instead than on top of the
no-species empty state. For the differential model we find and describe in
detail the stable localized states. For the full nonlocal model, however
competitive interactions alone do not allow the conditions for the observation
of self-localized states, and we show how the inclusion of additional
facilitative interactions lead to the appearance of them.


  In this paper we have proposed a model for the distribution of allelic
probabilities for generating populations as reliably as possible. Our objective
was to develop such a model which would allow simulating allelic probabilities
with different observed truncation and de- gree of noise. In addition, we have
also introduced here a complete new approach to analyze a genome-wide
association study (GWAS) dataset, starting from a new test of association with
a statistical distribution and two effect sizes of each genotype. The new
methodologi- cal approach was applied to a real data set together with a Monte
Carlo experiment which showed the power performance of our new method. Finally,
we compared the new method based on beta distribution with the conventional
method (based on Chi-Squared distribu- tion) using the agreement Kappa index
and a principal component analysis (PCA). Both the analyses show found
differences existed between both the approaches while selecting the single
nucleotide polymorphisms (SNPs) in association.


  There is now a certain consensus that Transcription Factors (TFs) reach their
target sites, where they regulate gene transcription, via a mechanism dubbed
facilitated diffusion (FD). In FD, the TF cycles between events of 3D-diffusion
in solution (jumps), 1D-diffusion along DNA (sliding), and small jumps
(hopping) achieving association rates higher than for 3D-diffusion only. We
investigate the FD phenomenology through Molecular Dynamics simulations in the
framework of coarse-grained modeling. Despite the crude approximations, we show
that the model generates, upon varying equilibrium distance of the DNA-TF
interaction, a phenomenology matching a number of experimental and numerical
results obtained with more refined models. In particular, by focusing on the
kinematics of the process, we characterize the geometrical properties of TF
trajectories during sliding. We find that sliding occurs via helical paths
around the DNA helix leading to a coupling of translation along the DNA-axis
with rotation around it. The 1D-diffusion constant measured in simulations is
found to be interwoven with the geometrical properties of sliding and we
develop a simple argument able to quantitatively reproduce the measured values.


  We study the question of existence of positive steady states of nonlinear
evolution equations. We recast the steady state equation in the form of
eigenvalue problems for a parametrised family of unbounded linear operators,
which are generators of strongly continuous semigroups; and a fixed point
problem. In case of irreducible governing semigroups we consider evolution
equations with non-monotone nonlinearities of dimension two, and we establish a
new fixed point theorem for set-valued maps. In case of reducible governing
semigroups we establish results for monotone nonlinearities of any finite
dimension $n$. In addition, we establish a non-quasinilpotency result for a
class of strictly positive operators, which are neither irreducible nor
compact, in general. We illustrate our theoretical results with examples of
partial differential equations arising in structured population dynamics. In
particular, we establish existence of positive steady states of a
size-structured juvenile-adult and a structured consumer-resource population
model, as well as for a selection-mutation model with distributed recruitment
process.


  To understand biological diversification, it is important to account for
large-scale processes that affect the evolutionary history of groups of
co-distributed populations of organisms. Such events predict temporally
clustered divergences times, a pattern that can be estimated using genetic data
from co-distributed species. I introduce a new approximate-Bayesian method for
comparative phylogeographical model-choice that estimates the temporal
distribution of divergences across taxa from multi-locus DNA sequence data. The
model is an extension of that implemented in msBayes. By reparameterizing the
model, introducing more flexible priors on demographic and divergence-time
parameters, and implementing a non-parametric Dirichlet-process prior over
divergence models, I improved the robustness, accuracy, and power of the method
for estimating shared evolutionary history across taxa. The results demonstrate
the improved performance of the new method is due to (1) more appropriate
priors on divergence-time and demographic parameters that avoid prohibitively
small marginal likelihoods for models with more divergence events, and (2) the
Dirichlet-process providing a flexible prior on divergence histories that does
not strongly disfavor models with intermediate numbers of divergence events.
The new method yields more robust estimates of posterior uncertainty, and thus
greatly reduces the tendency to incorrectly estimate models of shared
evolutionary history with strong support.


  Cancer progression is an example of a rapid adaptive process where evolving
new traits is essential for survival and requires a high mutation rate.
Precancerous cells acquire a few key mutations that drive rapid population
growth and carcinogenesis. Cancer genomics demonstrates that these few 'driver'
mutations occur alongside thousands of random 'passenger' mutations-a natural
consequence of cancer's elevated mutation rate. Some passengers can be
deleterious to cancer cells, yet have been largely ignored in cancer research.
In population genetics, however, the accumulation of mildly deleterious
mutations has been shown to cause population meltdown. Here we develop a
stochastic population model where beneficial drivers engage in a tug-of-war
with frequent mildly deleterious passengers. These passengers present a barrier
to cancer progression that is described by a critical population size, below
which most lesions fail to progress, and a critical mutation rate, above which
cancers meltdown. We find support for the model in cancer age-incidence and
cancer genomics data that also allow us to estimate the fitness advantage of
drivers and fitness costs of passengers. We identify two regimes of adaptive
evolutionary dynamics and use these regimes to rationalize successes and
failures of different treatment strategies. We find that a tumor's load of
deleterious passengers can explain previously paradoxical treatment outcomes
and suggest that it could potentially serve as a biomarker of response to
mutagenic therapies. Collective deleterious effect of passengers is currently
an unexploited therapeutic target. We discuss how their effects might be
exacerbated by both current and future therapies.


  Microbes require several complex organic molecules for growth. A species may
obtain a required factor by taking up molecules released by other species or by
synthesizing the molecule. The patterns of uptake and synthesis set a flow of
resources through the multiple species that create a microbial community. This
article analyzes a simple mathematical model of the tradeoff between uptake and
synthesis. Key factors include the influx rate from external sources relative
to the outflux rate, the rate of internal decay within cells, and the cost of
synthesis. Aspects of demography also matter, such as cellular birth and death
rates, the expected time course of a local resource flow, and the associated
lifespan of the local population. Spatial patterns of genetic variability and
differentiation between populations may also strongly influence the evolution
of metabolic regulatory controls of individual species and thus the structuring
of microbial communities. The widespread use of optimality approaches in recent
work on microbial metabolism has ignored demography and genetic structure.


  It has not been known whether preferential dispersal is adaptive in
fluctuating environments. We investigate the effect of preferential and random
dispersals in bet-hedging systems by using a discrete stochastic metapopulation
model, where each site fluctuates between good and bad environments with
temporal correlation. To explore the optimal migration pattern, an analytical
estimation of the total growth is derived by mean field approximation. We found
that the preference for fertile sites is disadvantageous when transportation
among sites has a cost or the sensitivity of preference is high.


  Establishing that a set of population-splitting events occurred at the same
time can be a potentially persuasive argument that a common process affected
the populations. Oaks et al. (2013) assessed the ability of an
approximate-Bayesian method (msBayes) to estimate such a pattern of
simultaneous divergence across taxa, to which Hickerson et al. (2014)
responded. Both papers agree the method is sensitive to prior assumptions and
often erroneously supports shared divergences; the papers differ about the
explanation and solution. Oaks et al. (2013) suggested the method's behavior is
caused by the strong weight of uniform priors on divergence times leading to
smaller marginal likelihoods of models with more divergence-time parameters
(Hypothesis 1); they proposed alternative priors to avoid strongly weighted
posteriors. Hickerson et al. (2014) suggested numerical approximation error
causes msBayes analyses to be biased toward models of clustered divergences
(Hypothesis 2); they proposed using narrow, empirical uniform priors. Here, we
demonstrate that the approach of Hickerson et al. (2014) does not mitigate the
method's tendency to erroneously support models of clustered divergences, and
often excludes the true parameter values. Our results also show that the
tendency of msBayes analyses to support models of shared divergences is
primarily due to Hypothesis 1. This series of papers demonstrate that if our
prior assumptions place too much weight in unlikely regions of parameter space
such that the exact posterior supports the wrong model of evolutionary history,
no amount of computation can rescue our inference. Fortunately, more flexible
distributions that accommodate prior uncertainty about parameters without
placing excessive weight in vast regions of parameter space with low likelihood
increase the method's robustness and power to detect temporal variation in
divergences.


  We briefly review some of the models used to describe DNA denaturation
dynamics, focusing on the value of the dynamical exponent $z$, which governs
the scaling of the characteristic time $\tau\sim L^z$ as a function of the
sequence length $L$. The models contain different degrees of simplifications,
in particular sometimes they do not include a description for helical
entanglement: we discuss how this aspect influences the value of $z$, which
ranges from $z=0$ to $z \approx 3.3$. Connections with experiments are also
mentioned.


  Literature data on the properties of DNA open states are reviewed and
analyzed. These states are formed as a result of strong DNA fluctuations and
have a great impact on a number of biochemical processes; among them is charge
transfer in DNA, for example. A comparative analysis of experimental data on
the kinetics and thermodynamics of DNA open states for a wide temperature range
was carried out. Discrepancies between the results of various experiments have
been explained. Three types of DNA open states are recognized based on their
differences in thermodynamic properties and other characteristics. Besides, an
up-to-date definition of the term "open state" is given. A review is carried
out for simple mathematical models of DNA in most of which the state of one
pair is described by one or two variables. The main problems arising in
theoretical investigations of heterogeneous DNA in the framework of models of
this level are considered. The role of each group of models in interpretation
of experimental data is discussed. Special consideration is given to the
studies of the transfer and localization of the nucleotide pairs oscillations'
energy by mechanical models. These processes are shown to play a key role in
the dynamics of a heterogeneous duplex. Their theoretical interpretation is
proven to be very important for the development of modern molecular biology and
biophysics. The main features of the theoretical approaches are considered
which enabled describing various experimental data. Prospects of the models'
development are described, particular details of their optimization are
suggested, and possible ways of modernization of some experimental techniques
are discussed.


  The formation of dynamical clusters of proteins is ubiquitous in cellular
membranes and is in part regulated by the recycling of membrane components.
Mean-field models of out-of-equilibrium cluster formation with recycling
predict a broad cluster size distribution for infinite systems and must be
corrected for finite-size effects for small systems such as cellular
organelles. We show, using stochastic simulations and analytic modelling, that
tuning the system size is an efficient way to control the size of lateral
membrane heterogeneities. We apply these findings to a chain of enzymatic
reaction sensitive to membrane protein clustering. The reaction efficiency is
found to be a non-monotonic function of the system size, and can be optimal for
sizes comparable to those of cellular organelles.


  Game theory provides a quantitative framework for analyzing the behavior of
rational agents. The Iterated Prisoner's Dilemma in particular has become a
standard model for studying cooperation and cheating, with cooperation often
emerging as a robust outcome in evolving populations. Here we extend
evolutionary game theory by allowing players' strategies as well as their
payoffs to evolve in response to selection on heritable mutations. In nature,
many organisms engage in mutually beneficial interactions, and individuals may
seek to change the ratio of risk to reward for cooperation by altering the
resources they commit to cooperative interactions. To study this, we construct
a general framework for the co-evolution of strategies and payoffs in arbitrary
iterated games. We show that, as payoffs evolve, a trade-off between the
benefits and costs of cooperation precipitates a dramatic loss of cooperation
under the Iterated Prisoner's Dilemma; and eventually to evolution away from
the Prisoner's Dilemma altogether. The collapse of cooperation is so extreme
that the average payoff in a population may decline, even as the potential
payoff for mutual cooperation increases. Our work offers a new perspective on
the Prisoner's Dilemma and its predictions for cooperation in natural
populations; and it provides a general framework to understand the co-evolution
of strategies and payoffs in iterated interactions.


  Heterogeneity in biological molecules, resulting in molecule-to-molecule
variations in their dynamics and function, is an emerging theme. To elucidate
the consequences of heterogeneous behavior at the single molecule level, we
propose an exactly solvable model in which the unfolding rate due to mechanical
force depends parametrically on an auxiliary variable representing an entropy
barrier arising from fluctuations in internal dynamics. When the rate of
fluctuations, a measure of dynamical disorder, is comparable to or smaller than
the rate of force-induced unbinding, we show that there are two experimentally
observable consequences: non-exponential survival probability at constant
force, and a heavy-tailed rupture force distribution at constant loading rate.
By fitting our analytical expressions to data from single molecule pulling
experiments on proteins and DNA, we quantify the extent of disorder. We show
that only by analyzing data over a wide range of forces and loading rates can
the role of disorder due to internal dynamics be quantitatively assessed.


  Motivated by phylogenetics, our aim is to obtain a system of equations that
define a phylogenetic variety on an open set containing the biologically
meaningful points. In this paper we consider phylogenetic varieties defined via
group-based models. For any finite abelian group $G$, we provide an explicit
construction of $codim X$ phylogenetic invariants (polynomial equations) of
degree at most $|G|$ that define the variety $X$ on a Zariski open set $U$. The
set $U$ contains all biologically meaningful points when $G$ is the group of
the Kimura 3-parameter model. In particular, our main result confirms a
conjecture by the third author and, on the set $U$, a couple of conjectures by
Bernd Sturmfels and Seth Sullivant.


  Recent work identified the fundamental limits on the information requirements
in terms of read length and coverage depth required for successful de novo
genome reconstruction from shotgun sequencing data, based on the idealistic
assumption of no errors in the reads (noiseless reads). In this work, we show
that even when there is noise in the reads, one can successfully reconstruct
with information requirements close to the noiseless fundamental limit. A new
assembly algorithm, X-phased Multibridging, is designed based on a
probabilistic model of the genome. It is shown through analysis to perform well
on the model, and through simulations to perform well on real genomes.


  In this paper we explore the eco-evolutionary dynamics of a predator-prey
model, where the prey population is structured according to a certain life
history trait. The trait distribution within the prey population is the result
of interplay between genetic inheritance and mutation, as well as selectivity
in the consumption of prey by the predator. The evolutionary processes are
considered to take place on the same time scale as ecological dynamics, i.e. we
consider the evolution to be rapid. Previously published results show that
population structuring and rapid evolution in such predator-prey system can
stabilise an otherwise globally unstable dynamics even with an unlimited
carrying capacity of prey. However, those findings were only based on direct
numerical simulation of equations and obtained for particular parametrisations
of model functions, which obviously calls into question the correctness and
generality of the previous results. The main objective of the current study is
to treat the model analytically and consider various parametrisations of
predator selectivity and inheritance kernel. We investigate the existence of a
coexistence stationary state in the model and carry out stability analysis of
this state. We derive expressions for the Hopf bifurcation curve which can be
used for constructing bifurcation diagrams in the parameter space without the
need for a direct numerical simulation of the underlying integro-differential
equations. We analytically show the possibility of stabilisation of a globally
unstable predator-prey system with prey structuring. We prove that the
coexistence stationary state is stable when the saturation in the predation
term is low. Finally, for a class of kernels describing genetic inheritance and
mutation we show that stability of the predator-prey interaction will require a
selectivity of predation according to the life trait.


  Topological features of gene regulatory networks can be successfully
reproduced by a model population evolving under selection for short dynamical
attractors. The evolved population of networks exhibit motif statistics,
summarized by significance profiles, which closely match those of {\it E. coli,
S. cerevsiae} and {\it B. subtilis}, in such features as the excess of linear
motifs and feed-forward loops, and deficiency of feedback loops. The slow
relaxation to stasis is a hallmark of a rugged fitness landscape, with
independently evolving populations exploring distinct valleys strongly
differing in network properties.


  We present a model for obligated mutualistic associations, in which two
transmissible diseases are allowed to infect just one population. As the
general model proves too hard to be fully analytically investigated, some
special cases are analysed. Among our findings, the coexistence of the two
strains does not appear possible, under the model assumptions. Furthermore, in
particularly unfavorable circumstances the ecosystem may disappear. In this
respect, an accurate computation of the basin of attraction of the origin is
provided using novel techniques. For this obligated mutualistic system the
presence of the diseases appears to be less relevant than in many other
circumstances in ecoepidemiology, including also the case of facultative
symbiotic associations.


  One of the challenges in Synthetic Biology is to design circuits with
increasing levels of complexity. While circuits in Biology are complex and
subject to natural tradeoffs, most synthetic circuits are simple in terms of
the number of regulatory regions, and have been designed to meet a single
design criterion. In this contribution we introduce a multiobjective
formulation for the design of biocircuits. We set up the basis for an advanced
optimization tool for the modular and systematic design of biocircuits capable
of handling high levels of complexity and multiple design criteria. Our
methodology combines the efficiency of global Mixed Integer Nonlinear
Programming solvers with multiobjective optimization techniques. Through a
number of examples we show the capability of the method to generate non
intuitive designs with a desired functionality setting up a priori the desired
level of complexity. The presence of more than one competing objective provides
a realistic design setting where every design solution represents a trade-off
between different criteria. The tool can be useful to explore and identify
different design principles for synthetic gene circuits.


  The recent increase in Bordetella pertussis incidence (whooping cough)
presents a challenge to global health. Recent studies have called into question
the effectiveness of acellular B. pertussis vaccination in reducing
transmission. Here we examine the epidemiological consequences of an
ineffective B. pertussis vaccine. Using a dynamic transmission model, we find
that: 1) an ineffective vaccine can account for the observed increase in B.
pertussis incidence; 2) asymptomatic infections can bias surveillance and upset
situational awareness of B. pertussis; and 3) vaccinating individuals in close
contact with infants too young to receive vaccine (so called "cocooning"
unvaccinated children) may be ineffective. Our results have important
implications for B. pertussis vaccination policy and paint a complicated
picture for achieving herd immunity and possible B. pertussis eradication.


  Atomistic simulations were performed on hydrated model lipid multilayers that
are representative of the lipid matrix in the outer skin (stratum corneum). We
find that cholesterol transfers easily between adjacent leaflets belonging to
the same bilayer via fast orientational diffusion (tumbling) in the
inter-leaflet disordered region, while at the same time there is a large free
energy cost against swelling. This fast flip-flop may play an important role in
accommodating the variety of curvatures that would be required in the three
dimensional arrangement of the lipid multilayers in skin, and for enabling
mechanical or hydration induced strains without large curvature elastic costs.


  Gregarious animals need to make collective decisions in order to keep their
cohesiveness. Several species of them live in multilevel societies, and form
herds composed of smaller communities. We present a model for the development
of a leadership hierarchy in a herd consisting of loosely connected sub-groups
(e.g. harems) by combining self organization and social dynamics. It starts
from unfamiliar individuals without relationships and reproduces the emergence
of a hierarchical and modular leadership network that promotes an effective
spreading of the decisions from more capable individuals to the others, and
thus gives rise to a beneficial collective decision. Our results stemming from
the model are in a good agreement with our observations of a Przewalski horse
herd (Hortob\'agy, Hungary). We find that the harem-leader to harem-member
ratio observed in Przewalski horses corresponds to an optimal network in this
approach regarding common success, and that the observed and modeled harem size
distributions are close to a lognormal.


  Correlation patterns in multiple sequence alignments of homologous proteins
can be exploited to infer information on the three-dimensional structure of
their members. The typical pipeline to address this task, which we in this
paper refer to as the three dimensions of contact prediction, is to: (i) filter
and align the raw sequence data representing the evolutionarily related
proteins; (ii) choose a predictive model to describe a sequence alignment;
(iii) infer the model parameters and interpret them in terms of structural
properties, such as an accurate contact map. We show here that all three
dimensions are important for overall prediction success. In particular, we show
that it is possible to improve significantly along the second dimension by
going beyond the pair-wise Potts models from statistical physics, which have
hitherto been the focus of the field. These (simple) extensions are motivated
by multiple sequence alignments often containing long stretches of gaps which,
as a data feature, would be rather untypical for independent samples drawn from
a Potts model. Using a large test set of proteins we show that the combined
improvements along the three dimensions are as large as any reported to date.