  A rigorous formalism for estimating noncovalent binding free energies and
thermodynamic expectations from calculations in which receptor configurations
are sampled independently from the ligand is derived. Due to this separation,
receptor configurations only need to be sampled once, facilitating the use of
binding free energy calculations in virtual screening. Demonstrative
calculations on a host-guest system yield good agreement with previous free
energy calculations and isothermal titration calorimetry measurements. Implicit
ligand theory provides guidance on how to improve existing molecular docking
algorithms and insight into the concepts of induced fit and conformational
selection in noncovalent macromolecular recognition.


  We study an individual based model describing competition in space between
two different alleles. Although the model is similar in spirit to classic
models of spatial population genetics such as the stepping stone model, here
however space is continuous and the total density of competing individuals
fluctuates due to demographic stochasticity. By means of analytics and
numerical simulations, we study the behavior of fixation probabilities,
fixation times, and heterozygosity, in a neutral setting and in cases where the
two species can compete or cooperate. By concluding with examples in which
individuals are transported by fluid flows, we argue that this model is a
natural choice to describe competition in marine environments.


  Conditional sampling distributions (CSDs), sometimes referred to as copying
models, underlie numerous practical tools in population genomic analyses.
Though an important application that has received much attention is the
inference of population structure, the explicit exchange of migrants at
specified rates has not hitherto been incorporated into the CSD in a principled
framework. Recently, in the case of a single panmictic population, a
sequentially Markov CSD has been developed as an accurate, efficient
approximation to a principled CSD derived from the diffusion process dual to
the coalescent with recombination. In this paper, the sequentially Markov CSD
framework is extended to incorporate subdivided population structure, thus
providing an efficiently computable CSD that admits a genealogical
interpretation related to the structured coalescent with migration and
recombination. As a concrete application, it is demonstrated empirically that
the CSD developed here can be employed to yield accurate estimation of a wide
range of migration rates.


  Characterizing time-evolution of allele frequencies in a population is a
fundamental problem in population genetics. In the Wright-Fisher diffusion,
such dynamics is captured by the transition density function, which satisfies
well-known partial differential equations. For a multi-allelic model with
general diploid selection, various theoretical results exist on representations
of the transition density, but finding an explicit formula has remained a
difficult problem. In this paper, a technique recently developed for a
diallelic model is extended to find an explicit transition density for an
arbitrary number of alleles, under a general diploid selection model with
recurrent parent-independent mutation. Specifically, the method finds the
eigenvalues and eigenfunctions of the generator associated with the
multi-allelic diffusion, thus yielding an accurate spectral representation of
the transition density. Furthermore, this approach allows for efficient,
accurate computation of various other quantities of interest, including the
normalizing constant of the stationary distribution and the rate of convergence
to this distribution.


  Creel surveys were performed over a three week period in late spring, 2011,
in the Lafourche and Calcasieu area estuaries of the Louisiana Gulf Coast.
Weights and lengths were measured for black drum (Pogonias cromis), red drum
(Sciaenops ocellatus), and spotted seatrout (Cynoscion nebulosus), and relative
condition factors were calculated relative to expected weights from the long
term (5 year) Louisiana data for each species. A normal relative condition
factor is 1.00. The mean relative condition factors in the Lafourche area were
black drum, 0.955 (0.020); red drum, 0.955 (0.011); spotted seatrout, 0.994
(0.009). In the Calcasieu area, the mean relative condition factors were black
drum, 0.934 (0.017); red drum, 0.965 (0.014); spotted seatrout, 0.971 (0.010).
Uncertainties are in parentheses. Results suggest that the abundance of primary
food sources for black drum and red drum in Lafourche, including oysters and
crab, were likely affected by the oil spill and continued to be reduced one
year later. Increased harvest of oysters and blue crab in the Calcasieu area
(in part to make up for the ban in most of Louisiana) resulted in less food for
the black drum and red drum there, also. Spotted seatrout eat mainly shrimp and
small fish and showed no significant reduction in relative condition factor in
Lafourche and a slight reduction in Calcasieu one year after the oil spill.


  A simple model of DNA based on two interacting polymers has been used to
study the unzipping of a double stranded DNA subjected to a periodic force. We
propose a dynamical transition, where without changing the physiological
condition, it is possible to bring DNA from the zipped/unzipped state to a new
dynamic (hysteretic) state by varying the frequency of the applied force. Our
studies reveal that the area of the hystersis loop grows with the same
exponents as of the isotropic spin systems. These exponents are amenable to
verification in the force spectroscopic experiments.


  We consider a single-species diffusion-limited annihilation reaction with
reactants confined to a two-dimensional surface with one arbitrarily large
dimension and the other comparable in size to interparticle distances. This
situation could describe reactants which undergo both longitudinal and
transverse diffusion on long filamentous molecules (such as microtubules), or
molecules that undergo truly one-dimensional translational diffusion (e.g. a
transcription factor on DNA) but simultaneously exhibit diffusive behavior in a
second dimension corresponding to a rotational or conformational degree of
freedom. We combine simple analytical arguments and Monte Carlo simulations to
show that the reaction rate law exhibits a crossover from one-dimensional to
two-dimensional diffusion as a function of particle concentration and the size
of the smaller dimension. In the case of a reversible binding reaction, the
diffusion-limited reaction rate is given by the Smoluchowski expression, but
the crossover is revealed in the statistics of particle collision histories.
The results can also be applied to a particle-antiparticle annihilation
reaction.


  Noise in transduction of chemotactic stimuli to the flagellar motor of E.
coli will affect the random run-and-tumble motion of the cell and the ability
to perform chemotaxis. Here we use numerical simulations to show that an
intermediate level of noise in the slow methylation dynamics enhances drift
while not compromising localisation near concentration peaks. A minimal model
shows how such an optimal noise level arises from the interplay of noise and
the dependence of the motor response on the network output. Our results suggest
that cells can exploit noise to improve chemotactic performance.


  In this work we examine the properties of a recently described ordinary
differential equation that relates the age-specific prevalence of a chronic
disease with the incidence and mortalities of the diseased and healthy persons.
The equation has been used to estimate the incidence from prevalence data,
which is an inverse problem. The ill-posedness of this problem is proven, too.


  Epidemics occur in all shapes and forms: infections propagating in our sparse
sexual networks, rumours and diseases spreading through our much denser social
interactions, or viruses circulating on the Internet. With the advent of large
databases and efficient analysis algorithms, these processes can be better
predicted and controlled. In this study, we use different characteristics of
network organization to identify the influential spreaders in 17 empirical
networks of diverse nature using 2 epidemic models. We find that a judicious
choice of local measures, based either on the network's connectivity at a
microscopic scale or on its community structure at a mesoscopic scale, compares
favorably to global measures, such as betweenness centrality, in terms of
efficiency, practicality and robustness. We also develop an analytical
framework that highlights a transition in the characteristic scale of different
epidemic regimes. This allows to decide which local measure should govern
immunization in a given scenario.


  We discuss the conditions under which a population of anomalously diffusing
individuals can be characterized by demographic fluctuations that are
anomalously scaling themselves. Two examples are provided in the case of
individuals migrating by Gaussian diffusion, and by a sequence of L\'evy
flights.


  Mounting evidence suggests that natural populations can harbor extensive
fitness diversity with numerous genomic loci under selection. It is also known
that genealogical trees for populations under selection are quantifiably
different from those expected under neutral evolution and described
statistically by Kingman's coalescent. While differences in the statistical
structure of genealogies have long been used as a test for the presence of
selection, the full extent of the information that they contain has not been
exploited. Here we shall demonstrate that the shape of the reconstructed
genealogical tree for a moderately large number of random genomic samples taken
from a fitness diverse, but otherwise unstructured asexual population can be
used to predict the relative fitness of individuals within the sample. To
achieve this we define a heuristic algorithm, which we test in silico using
simulations of a Wright-Fisher model for a realistic range of mutation rates
and selection strength. Our inferred fitness ranking is based on a linear
discriminator which identifies rapidly coalescing lineages in the reconstructed
tree. Inferred fitness ranking correlates strongly with actual fitness, with a
genome in the top 10% ranked being in the top 20% fittest with false discovery
rate of 0.1-0.3 depending on the mutation/selection parameters. The ranking
also enables to predict the genotypes that future populations inherit from the
present one. While the inference accuracy increases monotonically with sample
size, samples of 200 nearly saturate the performance. We propose that our
approach can be used for inferring relative fitness of genomes obtained in
single-cell sequencing of tumors and in monitoring viral outbreaks.


  In this paper, we study the $SIS$ (susceptible-infected-susceptible) and
$SIR$ (susceptible-infected-removed) epidemic models on undirected, weighted
networks by deriving pairwise-type approximate models coupled with
individual-based network simulation. Two different types of
theoretical/synthetic weighted network models are considered. Both models start
from non-weighted networks with fixed topology followed by the allocation of
link weights in either (i) random or (ii) fixed/deterministic way. The pairwise
models are formulated for a general discrete distribution of weights, and these
models are then used in conjunction with network simulation to evaluate the
impact of different weight distributions on epidemic threshold and dynamics in
general. For the $SIR$ dynamics, the basic reproductive ratio $R_0$ is
computed, and we show that (i) for both network models $R_{0}$ is maximised if
all weights are equal, and (ii) when the two models are equally matched, the
networks with a random weight distribution give rise to a higher $R_0$ value.
The models are also used to explore the agreement between the pairwise and
simulation models for different parameter combinations.


  Multiple interlinked positive feedback loops shape the stimulus responses of
various biochemical systems, such as the cell cycle or intracellular calcium
release. Recent studies with simplified models have identified two advantages
of coupling fast and slow feedback loops. Namely, this dual-time structure
enables a fast response while enhancing resistances of responses and
bistability to stimulus noise. We now find that in addition: 1) the dual-time
structure confers resistance to internal noise due to molecule number
fluctuations, and 2) model variants with altered coupling, which better
represent some specific systems, share all the above advantages. We develop a
similar bistable model with a fast autoactivation loop coupled to a slow loop,
which minimally represents positive feedback that may be essential for
long-term synaptic potentiation (LTP). The advantages of fast response and
noise resistance carry over to this model. Empirically, LTP develops resistance
to reversal over ~1 h. The model suggests this resistance may result from
increased amounts of synaptic kinases involved in positive feedback.


  Cancer progression is driven by a small number of genetic alterations
accumulating in a neoplasm. These few driver alterations reside in a cancer
genome alongside tens of thousands of other mutations that are widely believed
to have no role in cancer and termed passengers. Many passengers, however, fall
within protein coding genes and other functional elements and can possibly have
deleterious effects on cancer cells. Here we investigate a potential of mildly
deleterious passengers to accumulate and alter the course of neoplastic
progression. Our approach combines evolutionary simulations of cancer
progression with the analysis of cancer sequencing data. In our simulations,
individual cells stochastically divide, acquire advantageous driver and
deleterious passenger mutations, or die. Surprisingly, despite selection
against them, passengers accumulate and largely evade selection during
progression. Although individually weak, the collective burden of passengers
alters the course of progression leading to several phenomena observed in
oncology that cannot be explained by a traditional driver-centric view. We
tested predictions of the model using cancer genomic data. We find that many
passenger mutations are likely to be damaging and that, in agreement with the
model, they have largely evaded purifying selection. Finally, we used our model
to explore cancer treatments that exploit the load of passengers by either 1)
increasing the mutation rate; or 2) exacerbating their deleterious effects.
While both approaches lead to cancer regression, the later leads to less
frequent relapse. Our results suggest a new framework for understanding cancer
progression as a balance of driver and passenger mutations.


  Helices are a key folding motif in protein structure. The question which
factors determine helix stability for a given polypeptide or protein is an
ongoing challenge. Here we use van der Waals corrected density-functional
theory to address a part of this question in a bottom-up approach. We show how
intrinsic helical structure is stabilized with length and temperature for a
series of experimentally well studied unsolvated alanine based polypeptides,
Ac-Alan-LysH+. By exploring extensively the conformational space of these
molecules, we find that helices emerge as the preferred structure in the length
range n=4-8 not just due to enthalpic factors (hydrogen bonds and their
cooperativity, van der Waals dispersion interactions, electrostatics), but
importantly also by a vibrational entropic stabilization over competing
conformers at room temperature. The stabilization is shown to be due to softer
low-frequency vibrational modes in helical conformers than in more compact
ones. This observation is corroborated by including anharmonic effects
explicitly through \emph{ab initio} molecular dynamics, and generalized by
testing different terminations and considering larger helical peptide models.


  In this Chapter we review our works on force fields for molecular simulations
of protein systems. We first discuss the functional forms of the force fields
and present some extensions of the conventional ones. We then present various
methods for force-field parameter optimizations. Finally, some examples of our
applications of these parameter optimization methods are given and they are
compared with the results from the existing force-fields.


  Summary: The Smith Waterman (SW) algorithm, which produces the optimal
pairwise alignment between two sequences, is frequently used as a key component
of fast heuristic read mapping and variation detection tools, but current
implementations are either designed as monolithic protein database searching
tools or are embedded into other tools. To facilitate easy integration of the
fast Single Instruction Multiple Data (SIMD) SW algorithm into third party
software, we wrote a C/C++ library, which extends Farrars Striped SW (SSW) to
return alignment information in addition to the optimal SW score. Availability:
SSW is available both as a C/C++ software library, as well as a stand alone
alignment tool wrapping the librarys functionality at
https://github.com/mengyao/Complete- Striped-Smith-Waterman-Library Contact:
marth@bc.edu


  The contact structure between hosts has a critical influence on disease
spread. However, most networkbased models used in epidemiology tend to ignore
heterogeneity in the weighting of contacts. This assumption is known to be at
odds with the data for many contact networks (e.g. sexual contact networks) and
to have a strong effect on the predictions of epidemiological models. One of
the reasons why models usually ignore heterogeneity in transmission is that we
currently lack tools to analyze weighted networks, such that most studies rely
on numerical simulations. Here, we present a novel framework to estimate key
epidemiological variables, such as the rate of early epidemic expansion and the
basic reproductive ratio, from joint probability distributions of number of
partners (contacts) and number of interaction events through which contacts are
weighted. This framework also allows for a derivation of the full time course
of epidemic prevalence and contact behaviour which is validated using numerical
simulations. Our framework allows for the incorporation of more realistic
contact networks into epidemiological models, thus improving predictions on the
spread of emerging infectious diseases.


  Phylogenetic diversity (PD) depends on sampling intensity, which complicates
the comparison of PD between samples of different depth. One approach to
dealing with differing sample depth for a given diversity statistic is to
rarefy, which means to take a random subset of a given size of the original
sample. Exact analytical formulae for the mean and variance of species richness
under rarefaction have existed for some time but no such solution exists for
PD. We have derived exact formulae for the mean and variance of PD under
rarefaction. We show that these formulae are correct by comparing exact
solution mean and variance to that calculated by repeated random (Monte Carlo)
subsampling of a dataset of stem counts of woody shrubs of Toohey Forest,
Queensland, Australia. We also demonstrate the application of the method using
two examples: identifying hotspots of mammalian diversity in Australasian
ecoregions, and characterising the human vaginal microbiome. There is a very
high degree of correspondence between the analytical and random subsampling
methods for calculating mean and variance of PD under rarefaction, although the
Monte Carlo method requires a large number of random draws to converge on the
exact solution for the variance. Rarefaction of mammalian PD of ecoregions in
Australasia to a common standard of 25 species reveals very different rank
orderings of ecoregions, indicating quite different hotspots of diversity than
those obtained for unrarefied PD. The application of these methods to the
vaginal microbiome shows that a classical score used to quantify bacterial
vaginosis is correlated with the shape of the rarefaction curve. The analytical
formulae for the mean and variance of PD under rarefaction are both exact and
more efficient than repeated subsampling. Rarefaction of PD allows for many
applications where comparisons of samples of different depth is required.


  Motivated by claims about the nature of the observed timescales in protein
systems said to fold downhill, we have studied the finite, linear master
equation which is a model of the downhill process. By solving for the system
eigenvalues, we prove the often stated claim that in situations where there is
no free energy barrier, a transition between single and multi-exponential
kinetics occurs at sufficient bias (towards the native state). Consequences for
protein folding, especially the downhill folding scenario, are briefly
discussed.


  Saccharomyces cerevisiae is one of the premier model systems for studying the
genomics and evolution of transposable elements. The availability of the S.
cerevisiae genome led to many insights into its five known transposable element
families (Ty1-Ty5) in the years shortly after its completion. However,
subsequent advances in bioinformatics tools for analysing transposable elements
and the recent availability of genome sequences for multiple strains and
species of yeast motivates new investigations into Ty evolution in S.
cerevisiae. Here we provide a comprehensive phylogenetic and population genetic
analysis of Ty families in S. cerevisiae based on a reannotation of Ty elements
in the S288c reference genome. We show that previous annotation efforts have
underestimated the total copy number of Ty elements for all known families. In
addition, we identify a new family of Ty3-like elements related to the S.
paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic
analyses of LTR sequences identified three families with short-branch, recently
active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4),
one family with essentially all recently active elements (Ty2) and two families
with only inactive elements (Ty3p and Ty5). Population genomic data from 38
additional strains of S. cerevisiae show that elements present in active clades
are predominantly polymorphic, whereas most of the inactive elements are fixed.
Finally, we use comparative genomic data to provide evidence that the Ty2 and
Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer.
Our results demonstrate that the genome of a single individual contains
important information about the state of TE population dynamics within a
species and suggest that horizontal transfer may play an important role in
shaping the diversity of transposable elements in unicellular eukaryotes.


  Background: Transposable elements (TEs) in eukaryote genomes are
quantitatively the main components affecting genome size, structure and
expression. The dynamics of their insertion and deletion depend on diverse
factors varying in strength and nature along the genome. We address here how TE
sequence evolution is affected by neighboring genes and the chromatin status
(euchromatin or heterochromatin) at their insertion site. Results: We estimated
ages of TE sequences in Arabidopsis thaliana, and found that they depend on the
distance to the nearest genes: TEs located close to genes are older than those
that are more distant. Consequently, TE sequences in heterochromatic regions,
which are gene-poor regions, are surprisingly younger and longer than that
elsewhere. Conclusions: We provide evidence for biased TE age distribution
close or near to genes. Interestingly, TE sequences in euchromatin and those in
heterochromatin evolve at different rates, and as a result, could explain that
TE sequences in heterochromatin tend to be younger and longer. Then, we revisit
models of TE sequence dynamics and point out differences for TE-rich genomes,
such as maize and wheat, compared to TE-poor genomes such as fly and A.
thaliana.


  We investigate the fine structure of the sp3 hybridized covalent bond
geometry that governs the tetrahedral architecture around the central
C$_\alpha$ carbon of a protein backbone, and for this we develop new
visualization techniques to analyze high resolution X-ray structures in Protein
Data Bank. We observe that there is a correlation between the deformations of
the ideal tetrahedral symmetry and the local secondary structure of the
protein. We propose a universal coarse grained energy function to describe the
ensuing side-chain geometry in terms of the C$_\beta$ carbon orientations. The
energy function can model the side-chain geometry with a sub-atomic precision.
As an example we construct the C$_\alpha$-C$_\beta$ structure of HP35 chicken
villin headpiece. We obtain a configuration that deviates less than 0.4 \.A in
root-mean-square distance from the experimental X-ray structure.


  Folded proteins have a modular assembly. They are constructed from regular
secondary structures like alpha-helices and beta-strands that are joined
together by loops. Here we develop a visualization technique that is adapted to
describe this modular structure. In complement to the widely employed
Ramachandran plot that is based on toroidal geometry, our approach utilizes the
geometry of a two-sphere. Unlike the more conventional approaches that only
describea given peptide unit, ours is capable of describing the entire backbone
environment including the neighboring peptide units. It maps the positions of
each atom to the surface of the two-sphere exactly how these atoms are seen by
an observer who is located at the position of the central C-alpha atom. At each
level of side-chain atoms we observe a strong correlation between the
positioning of the atom and the underlying local secondary structure with very
little if any variation between the different amino acids. As a concrete
example we analyze the left-handed helix region of non-glycyl amino acids. This
region corresponds to an isolated and highly localized residue independent
sector in the direction of the C-beta carbons on the two-sphere. We show that
the residue independent localization extends to C-gamma and C-delta carbons,
and to side-chain oxygen and nitrogen atoms in the case of asparagine and
aspartic acid. When we extend the analysis to the side-chain atoms of the
neighboring residues, we observe that left-handed beta-turns display a regular
and largely amino acid independent structure that can extend to seven
consecutive residues. This collective pattern is duu to the presence of a
backbone soliton. We show how one can use our visualization techniques to
analyze and classify the different solitons in terms of selection rules that we
describe in detail.


  The Hox gene collinearity enigma has often been approached using models based
on biomolecular mechanisms. The biophysical model, is an alternative approach,
speculating that collinearity is caused by physical forces pulling the Hox
clusters from a territory where they are inactive to a distinct spatial domain
where they are activated in a step by step manner.
  Hox gene translocations have recently been observed in support of the
biophysical model. Furthermore, genetic engineering experiments, performed in
embryonic mice, gave rise to some unexpected mutant expressions that
biomolecular models could not predict. In several cases when anterior Hoxd
genes are deleted, the expression of the genes whose expression is probed in
the mutants are impossible to anticipate. On the contrary, the biophysical
model offers convincing explanation.
  All these experimental results support the idea of physical forces being
responsible for Hox gene collinearity. In order to test the validity of the
various models further, certain experiment involving gene deletions are
proposed. The biophysical and biomolecular models predict different results for
these experiments, hence the expected outcome will confirm or question the
validity of these models.


  Anti-retroviral drugs can reduce the infectiousness of people living with HIV
by about 96%--treatment as prevention or TasP--and can reduce the risk of being
infected by an HIV positive person by about 70%--pre-exposure prophylaxis or
PrEP--raising the prospect of using anti-retroviral drugs to stop the epidemic
of HIV. The question as to which is more effective, more affordable and more
cost effective, and under what conditions, continues to be debated in the
scientific literature. Here we compare TasP and PreP in order to determine the
conditions under which each strategy is favourable. This analysis suggests that
where the incidence of HIV is less than 5% or the risk-reduction under PrEP is
less than 50%, TasP is favoured over PrEP; otherwise PrEP is favoured over
TasP. The potential for using PreP should therefore be restricted to those
among whom the annual incidence of HIV is greater than 5% and TasP reduces
transmission by more than 50%. PreP should be considered for commercial sex
workers, young women aged about 20 to 25 years, men-who-have-sex with men, or
intravenous drug users, but only where the incidence of HIV is high.


  Colonies of bacteria grown on thin agar plate exhibit fractal patterns as a
result of adaptation to their environments. The bacterial colony pattern
formation is regulated crucially by chemotaxis, the movement of cells along a
chemical concentration gradient. Here, the dynamics of pattern formation in
bacterial colony is investigated theoretically through a continuum model that
considers chemotaxis. In the case of the gradient sensed by the bacterium is
nearly uniform, the bacterial colony patterns are self-similar, which they look
the same at every scale. The scaling law of the bacterial colony growth has
been revealed explicitly. Chemotaxis biases the movement of bacterial
population in colony trend toward the chemical attractant. Moreover, the
bacterial colonies evolve long time as the traveling wave with sharp front.


  A common view in evolutionary biology is that mutation rates are minimised.
However, studies in combinatorial optimisation and search have shown a clear
advantage of using variable mutation rates as a control parameter to optimise
the performance of evolutionary algorithms. Much biological theory in this area
is based on Ronald Fisher's work, who used Euclidean geometry to study the
relation between mutation size and expected fitness of the offspring in
infinite phenotypic spaces. Here we reconsider this theory based on the
alternative geometry of discrete and finite spaces of DNA sequences. First, we
consider the geometric case of fitness being isomorphic to distance from an
optimum, and show how problems of optimal mutation rate control can be solved
exactly or approximately depending on additional constraints of the problem.
Then we consider the general case of fitness communicating only partial
information about the distance. We define weak monotonicity of fitness
landscapes and prove that this property holds in all landscapes that are
continuous and open at the optimum. This theoretical result motivates our
hypothesis that optimal mutation rate functions in such landscapes will
increase when fitness decreases in some neighbourhood of an optimum, resembling
the control functions derived in the geometric case. We test this hypothesis
experimentally by analysing approximately optimal mutation rate control
functions in 115 complete landscapes of binding scores between DNA sequences
and transcription factors. Our findings support the hypothesis and find that
the increase of mutation rate is more rapid in landscapes that are less
monotonic (more rugged). We discuss the relevance of these findings to living
organisms.


  From a coarse-grained perspective the motif of a self-activating species,
activating a second species which acts as its own repressor, is widely found in
biological systems, in particular in genetic systems with inherent oscillatory
behavior. Here we consider a specific realization of this motif as a genetic
circuit, in which genes are described as directly producing proteins, leaving
out the intermediate step of mRNA production. We focus on the effect that
inherent time scales on the underlying fine-grained scale can have on the
bifurcation patterns on a coarser scale in time. Time scales are set by the
binding and unbinding rates of the transcription factors to the promoter
regions of the genes. Depending on the ratio of these rates to the decay times
of the proteins, the appropriate averaging procedure for obtaining a
coarse-grained description changes and leads to sets of deterministic
equations, which differ in their bifurcation structure. In particular the
desired intermediate range of regular limit cycles fades away when the binding
rates of genes are of the same order or less than the decay time of at least
one of the proteins. Our analysis illustrates that the common topology of the
widely found motif alone does not necessarily imply universal features in the
dynamics.


  To better understand and manage complex systems like ecosystems it is
critical to know the relative contribution of system components to system
functioning. Ecologists and social scientists have described many ways that
individuals can be important; This paper makes two key contributions to this
research area. First, it shows that throughflow, the total energy-matter
entering or exiting a system component, is a global indicator of the relative
contribution of the component to the whole system activity. It is global
because it includes the direct and indirect exchanges among community members.
Further, throughflow is a special case of Hubbell status as defined in social
science. This recognition effectively joins the concepts, enabling ecologists
to use and build on the broader centrality research in network science. Second,
I characterize the distribution of throughflow in 45 empirically-based trophic
ecosystem models. Consistent with expectations, this analysis shows that a
small fraction of the system components are responsible for the majority of the
system activity. In 73% of the ecosystem models, 20% or less of the nodes
generate 80% or more of the total system throughflow. Four or fewer dominant
nodes are required to account for 50% of the total system activity. 121 of the
130 dominant nodes in the 45 ecosystem models could be classified as primary
producers, dead organic matter, or bacteria. Thus, throughflow centrality
indicates the rank power of the ecosystems components and shows the power
concentration in the primary production and decomposition cycle. Although these
results are specific to ecosystems, these techniques build on flow analysis
based on economic input-output analysis. Therefore these results should be
useful for ecosystem ecology, industrial ecology, the study of urban
metabolism, as well as other domains using input-output analysis.


  We apply recently developed inference methods based on general coalescent
processes to DNA sequence data obtained from various marine species. Several of
these species are believed to exhibit so-called shallow gene genealogies,
potentially due to extreme reproductive behaviour, e.g. via Hedgecock's
"reproduction sweepstakes". Besides the data analysis, in particular the
inference of mutation rates and the estimation of the (real) time to the most
recent common ancestor, we briefly address the question whether the genealogies
might be adequately described by so-called Beta coalescents (as opposed to
Kingman's coalescent), allowing multiple mergers of genealogies.
  The choice of the underlying coalescent model for the genealogy has drastic
implications for the estimation of the above quantities, in particular the
real-time embedding of the genealogy.


  Metabolic networks consist of linked functional components, or modules. The
mechanism underlying metabolic network modularity is of great interest not only
to researchers of basic science but also to those in fields of engineering.
Previous studies have suggested a theoretical model, which proposes that a
change in the evolutionary goal (system-specific purpose) increases network
modularity, and this hypothesis was supported by statistical data analysis.
Nevertheless, further investigation has uncovered additional possibilities that
might explain the origin of network modularity. In this work, we propose an
evolving network model without tuning parameters to describe metabolic
networks. We demonstrate, quantitatively, that metabolic network modularity can
arise from simple growth processes, independent of the change in the
evolutionary goal. Our model is applicable to a wide range of organisms, and
appears to suggest that metabolic network modularity can be more simply
determined than previously thought. Nonetheless, our proposition does not serve
to contradict the previous model; it strives to provide an insight from a
different angle in the ongoing efforts to understand metabolic evolution, with
the hope of eventually achieving the synthetic engineering of metabolic
networks.


  The graph isomorphism problem has a long history in mathematics and computer
science, with applications in computational chemistry and biology, and it is
believed to be neither solvable in polynomial time nor NP-complete. E. Luks
proposed in 1982 the best algorithm so far for the solution of this problem,
which moreover runs in polynomial time if an upper bound for the degrees of the
nodes in the graphs is taken as a constant. Unfortunately, Luks' algorithm is
purely theoretical, very difficult to use in practice, and, in particular, we
have not been able to find any implementation of it in the literature. The main
goal of this paper is to present an efficient implementation of this algorithm
for ternary graphs in the SAGE system, as well as an adaptation to fully
resolved rooted phylogenetic networks on a given set of taxa.


  The X chromosome is present as a single copy in the heterogametic sex, and
this hemizygosity is expected to drive unusual patterns of evolution on the X
relative to the autosomes. For example, the hemizgosity of the X may lead to a
lower chromosomal effective population size compared to the autosomes
suggesting that the X might be more strongly affected by genetic drift.
However, the X may also experience stronger positive selection than the
autosomes because recessive beneficial mutations will be more visible to
selection on the X where they will spend less time being masked by the
dominant, less beneficial allele - a proposal known as the faster-X hypothesis.
Thus, empirical studies demonstrating increased genetic divergence on the X
chromosome could be indicative of either adaptive or non-adaptive evolution. We
measured gene expression in Drosophila species and in D. melanogaster inbred
strains for both embryos and adults. In the embryos we found that expression
divergence is on average more than 20% higher for genes on the X chromosome
relative to the autosomes, but in contrast, in the inbred strains gene
expression variation is significantly lower on the X chromosome. Furthermore,
expression divergence of genes on Muller's D element is significantly greater
along the branch leading to the obscura sub-group, in which this element
segregates as a neo-X chromosome. In the adults, divergence is greatest on the
X chromosome for males, but not for females, yet in both sexes inbred strains
harbour the lowest level of gene expression variation on the X chromosome. We
consider different explanations for our results and conclude that they are most
consistent within the framework of the faster-X hypothesis.


  Self-associates of nucleic acid components (stacking trimers and tetramers of
the base pairs of nucleic acids) and short fragments of nucleic acids are
nanoparticles (linear sizes of these particles are more than 10 A. Modern
quantum-mechanical methods and softwares allow one to perform ab initio
calculations of the systems consisting of 150-200 atoms with enough large basis
sets (for example, 6-31G*). The aim of this work is to reveal the peculiarities
of molecular and electronic structures, as well as the energy features of
nanoparticles of nucleic acid components.


  A conservation law and stability, recovering phenomena and characteristic
patterns of a nonlinear dynamical system have been studied and applied to
biological and ecological systems. In our previous study, we proposed a system
of symmetric 2n-dimensional conserved nonlinear differential equations with
external perturbations. In this paper, competitive systems described by
2-dimensional nonlinear dynamical (ND) model with external perturbations are
applied to population cycles and recovering phenomena of systems from microbes
to mammals. The famous 10-year cycle of population density of Canadian lynx and
snowshoe hare is numerically analyzed. We find that a nonlinear dynamical
system with a conservation law is stable and generates a characteristic rhythm
(cycle) of population density, which we call the {\it standard rhythm} of a
nonlinear dynamical system. The stability and restoration phenomena are
strongly related to a conservation law and balance of a system. The {\it
standard rhythm} of population density is a manifestation of the survival of
the fittest to the balance of a nonlinear dynamical system.


  Self-replication is a capacity common to every species of living thing, and
simple physical intuition dictates that such a process must invariably be
fueled by the production of entropy. Here, we undertake to make this intuition
rigorous and quantitative by deriving a lower bound for the amount of heat that
is produced during a process of self-replication in a system coupled to a
thermal bath. We find that the minimum value for the physically allowed rate of
heat production is determined by the growth rate, internal entropy, and
durability of the replicator, and we discuss the implications of this finding
for bacterial cell division, as well as for the pre-biotic emergence of
self-replicating nucleic acids.


  Both linear mixed models (LMMs) and sparse regression models are widely used
in genetics applications, including, recently, polygenic modeling in
genome-wide association studies. These two approaches make very different
assumptions, so are expected to perform well in different situations. However,
in practice, for a given data set one typically does not know which assumptions
will be more accurate. Motivated by this, we consider a hybrid of the two,
which we refer to as a "Bayesian sparse linear mixed model" (BSLMM) that
includes both these models as special cases. We address several key
computational and statistical issues that arise when applying BSLMM, including
appropriate prior specification for the hyper-parameters, and a novel Markov
chain Monte Carlo algorithm for posterior inference. We apply BSLMM and compare
it with other methods for two polygenic modeling applications: estimating the
proportion of variance in phenotypes explained (PVE) by available genotypes,
and phenotype (or breeding value) prediction. For PVE estimation, we
demonstrate that BSLMM combines the advantages of both standard LMMs and sparse
regression modeling. For phenotype prediction it considerably outperforms
either of the other two methods, as well as several other large-scale
regression methods previously suggested for this problem. Software implementing
our method is freely available from
http://stephenslab.uchicago.edu/software.html


  This paper generalizes a previously published differential equation that
describes the relation between the age-specific incidence, remission, and
mortality of a disease with its prevalence. The underlying model is a simple
compartment model with three states (illness-death model). In contrast to the
former work, migration- and calendar time-effects are included. As an
application of the theoretical findings, a hypothetical example of an
irreversible disease is treated.


  The Taylor-Frank method for making kin selection models when fitness is a
nonlinear function of a continuous phenotype requires this function to be
differentiable. This assumption sometimes fails for biologically important
fitness functions, for instance in microbial data and the theory of repeated
n-person games, even when fitness functions are smooth and continuous. In these
cases, the Taylor-Frank methodology cannot be used, and a more general form of
direct fitness must replace the standard one to account for kin selection, even
under weak selection.


  The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes
and genetically-encoded protein catalysts, has been central to models for the
early evolution of life on Earth. A key part of such models is continuity
between the earliest stages in the evolution of life and the RNA repertoires of
extant lineages. Some assessments seem consistent with a diverse RNA world, yet
direct continuity between modern RNAs and an RNA world has not been
demonstrated for the majority of RNA families, and, anecdotally, many RNA
functions appear restricted in their distribution. Despite much discussion of
the possible antiquity of RNA families, no systematic analyses of RNA family
distribution have been performed. To chart the broad evolutionary history of
known RNA families, we performed comparative genomic analysis of over 3 million
RNA annotations spanning 1446 families from the Rfam 10 database. We report
that 99% of known RNA families are restricted to a single domain of life,
revealing discrete repertoires for each domain. For the 1% of RNA
families/clans present in more than one domain, over half show evidence of
horizontal gene transfer, and the rest show a vertical trace, indicating the
presence of a complex protein synthesis machinery in the Last Universal Common
Ancestor (LUCA) and consistent with the evolutionary history of the most
ancient protein-coding genes. However, with limited interdomain transfer and
few RNA families exhibiting demonstrable antiquity as predicted under RNA world
continuity, our results indicate that the majority of modern cellular RNA
repertoires have primarily evolved in a domain-specific manner.


  In addition to natural selection, adaptive evolution requires genetic
variation to proceed. Yet the G-matrix may have limited 'genetic degrees of
freedom', with certain combinations of trait values unavailable to evolution.
Such limitations are often referred to as genetic constraints. Unfortunately,
clear predictions about when to expect constraints are rarely available.
Therefore, we developed an experimental system that provides specific
predictions regarding constraints. Such tests are important as disagreements
persist regarding the evidence for genetic constraints, possibly due to
differences in methodology, study system or both. Numerous measures of genetic
constraints have been suggested, and generally focus on whether some axes of G
have eigenvalues=~0, indicating a lack of genetic variance.The mutation
Ultrabithorax1 causes a mild homeotic transformation of segmental identity. We
predicted that this mutation would induce a genetic constraint due to this
homeosis. We measured genetic co-variation for a set of traits in a panel of
strains with and without Ubx1. As expected, Ubx1 induced homeotic
transformations, and altered patterns of allometry. Yet, no changes in
correlational structure nor in the distribution of eigenvalues of G were
observed. We discuss the role of using genetic manipulations to refine
hypotheses of constraints in natural systems.


  We consider a stochastic process for the generation of species which combines
a Yule process with a simple model for hybridization between pairs of
co-existent species. We assume that the origin of the process, when there was
one species, occurred at an unknown time in the past, and we condition the
process on producing n species via the Yule process and a single hybridization
event. We prove results about the distribution of the time of the hybridization
event. In particular we calculate a formula for all moments, and show that
under various conditions, the distribution tends to an exponential with rate
twice that of the birth rate for the Yule process.


  Fungal symbionts are often overlooked in studies of plant invasion.
Nevertheless, their role could be essential to the competitive success of the
invader. We studied fungal endophytes in the widespread invasive Centaurea
stoebe (common knapweed). A preliminary experiment showed that endophytes in
roots of C. stoebe significantly reduced the biomass of evolutionarily na\"ive
neighbours (Festuca idahoensis), compared to endophyte-free C. stoebe. In the
main experiment non-clavicipitaceous endophytes belonging to six phylotypes,
were employed as root inoculants. Each of these endophytes again reduced the
growth of na\"ive neighbours (F. idahoensis); and remarkably, each also
increased the growth of adapted neighbours (F. ovina) that were tested for the
first time. Four of the six endophytes caused C. stoebe to gain a competitive
advantage over its na\"ive neighbour that was significantly greater than the
endophyte-free C. stoebe over that same neighbour. However, endophyte-free C.
stoebe had no greater competitive advantage over F. idahoensis than it had over
F. ovina. Therefore, plant-plant interactions were dramatically affected by the
presence of endophytes in a way that would favor invasion.


  In a matching paper (arXiv:1207.5027), I proved that Conservation of Size and
Information in a discrete token based system is overwhelmingly likely to lead
to a power-law component size distribution with respect to the size of its
unique alphabet. This was substantiated to a very high level of significance
using some 55 million lines of source code of mixed provenance. The principle
was also applied to show that average gene length should be constant in an
animal kingdom where the same constraints appear to hold, the implication being
that Conservation of Information plays a similar role in discrete token-based
systems as the Conservation of Energy does in physical systems.
  In this part 2, the role of defect will be explored and a functional
behaviour for defect derived to be consistent with the power-law behaviour
substantiated above.
  This will be supported by further experimental data and the implications
explored.


  The number of algorithms available to reconstruct a biological network from a
dataset of high-throughput measurements is nowadays overwhelming, but
evaluating their performance when the gold standard is unknown is a difficult
task. Here we propose to use a few reconstruction stability tools as a
quantitative solution to this problem. We introduce four indicators to
quantitatively assess the stability of a reconstructed network in terms of
variability with respect to data subsampling. In particular, we give a measure
of the mutual distances among the set of networks generated by a collection of
data subsets (and from the network generated on the whole dataset) and we rank
nodes and edges according to their decreasing variability within the same set
of networks. As a key ingredient, we employ a global/local network distance
combined with a bootstrap procedure. We demonstrate the use of the indicators
in a controlled situation on a toy dataset, and we show their application on a
miRNA microarray dataset with paired tumoral and non-tumoral tissues extracted
from a cohort of 241 hepatocellular carcinoma patients.


  The enterobacteria lambda phage is a paradigm temperate bacteriophage. Its
lysogenic and lytic life cycles echo competition between the DNA binding
$\lambda$-repressor (CI) and CRO proteins. Here we scrutinize the structure,
stability and folding pathways of the $\lambda$-repressor protein, that
controls the transition from the lysogenic to the lytic state. We first
investigate the super-secondary helix-loop-helix composition of its backbone.
We use a discrete Frenet framing to resolve the backbone spectrum in terms of
bond and torsion angles. Instead of four, there appears to be seven individual
loops. We model the putative loops using an explicit soliton Ansatz. It is
based on the standard soliton profile of the continuum nonlinear Schr\"odinger
equation. The accuracy of the Ansatz far exceeds the B-factor fluctuation
distance accuracy of the experimentally determined protein configuration. We
then investigate the folding pathways and dynamics of the $\lambda$-repressor
protein. We introduce a coarse-grained energy function to model the backbone in
terms of the C$_\alpha$ atoms and the side-chains in terms of the relative
orientation of the C$_\beta$ atoms. We describe the folding dynamics in terms
of relaxation dynamics, and find that the folded configuration can be reached
from a very generic initial configuration. We conclude that folding is
dominated by the temporal ordering of soliton formation. In particular, the
third soliton should appear before the first and second. Otherwise, the DNA
binding turn does not acquire its correct structure. We confirm the stability
of the folded configuration by repeated heating and cooling simulations.


  Information processing and decision making is based upon logic operations,
which in cellular networks has been well characterized at the level of
transcription. In recent years however, both experimentalists and theorists
have begun to appreciate that cellular decision making can also be performed at
the level of a single protein, giving rise to the notion of protein logic. Here
we systematically explore protein logic using a well known statistical
mechanical model. As an example system, we focus on receptors which bind either
one or two ligands, and their associated dimers. Notably, we find that a single
heterodimer can realize any of the 16 possible logic gates, including the XOR
gate, by variation of biochemical parameters. We then introduce the novel idea
that a set of receptors with fixed parameters can encode functionally unique
logic gates simply by forming different dimeric combinations. An exhaustive
search reveals that the simplest set of receptors (two single-ligand receptors
and one double-ligand receptor) can realize several different groups of three
unique gates, a result for which the parametric analysis of single receptors
and dimers provides a clear interpretation. Both results underscore the
surprising functional freedom readily available to cells at the single-protein
level.


  In the commons, communities whose growth depends on public goods, individuals
often rely on surprisingly simple strategies, or heuristics, to decide whether
to contribute to the common good (at risk of exploitation by free-riders).
Although this appears a limitation, here we show how four heuristics lead to
sustainable growth by exploiting specific environmental constraints. The two
simplest ones --contribute permanently or switch stochastically between
contributing or not-- are first shown to bring sustainability when the public
good efficiently promotes growth. If efficiency declines and the commons is
structured in small groups, the most effective strategy resides in contributing
only when a majority of individuals are also contributors. In contrast, when
group size becomes large, the most effective behavior follows a minimal-effort
rule: contribute only when it is strictly necessary. Both plastic strategies
are observed in natural systems what presents them as fundamental social motifs
to successfully manage sustainability.


  The fitness of a biological strategy is typically measured by its expected
reproductive rate, the first moment of its offspring distribution. However,
strategies with high expected rates can also have high probabilities of
extinction. A similar situation is found in gambling and investment, where
strategies with a high expected payoff can also have a high risk of ruin. We
take inspiration from the gambler's ruin problem to examine how extinction is
related to population growth. Using moment theory we demonstrate how higher
moments can impact the probability of extinction. We discuss how moments can be
used to find bounds on the extinction probability, focusing on s-convex
ordering of random variables, a method developed in actuarial science. This
approach generates "best case" and "worst case" scenarios to provide upper and
lower bounds on the probability of extinction. Our results demonstrate that
even the most fit strategies can have high probabilities of extinction.


  Myopia, or nearsightedness, is the most common eye disorder, resulting
primarily from excess elongation of the eye. The etiology of myopia, although
known to be complex, is poorly understood. Here we report the largest ever
genome-wide association study (43,360 participants) on myopia in Europeans. We
performed a survival analysis on age of myopia onset and identified 19
significant associations (p < 5e-8), two of which are replications of earlier
associations with refractive error. These 19 associations in total explain 2.7%
of the variance in myopia age of onset, and point towards a number of different
mechanisms behind the development of myopia. One association is in the gene
PRSS56, which has previously been linked to abnormally small eyes; one is in a
gene that forms part of the extracellular matrix (LAMA2); two are in or near
genes involved in the regeneration of 11-cis-retinal (RGR and RDH5); two are
near genes known to be involved in the growth and guidance of retinal ganglion
cells (ZIC2, SFRP1); and five are in or near genes involved in neuronal
signaling or development. These novel findings point towards multiple genetic
factors involved in the development of myopia and suggest that complex
interactions between extracellular matrix remodeling, neuronal development, and
visual signals from the retina may underlie the development of myopia in
humans.


  In 1922, JBS Haldane discovered an intriguing bias of postzygotic isolation
during early speciation: the heterogametic sex of F1 hybrids between closely
related species or subspecies is more susceptible to sterility or inviability
than the homogametic sex. This phenomenon, now known as Haldane's rule, has
been repeatedly confirmed across broad taxa in diecious animals and plants.
Currently, the dominant view in the field of speciation genetics believes that
Haldane's rule for sterility, inviability, male heterogamety and female
heterogametic belongs to different entities; and Haldane's rule in these
subdivisions has different causes, which operate coincidentally and/or
collectively resulting in this striking bias against the heterogametic sex in
hybridization. This view, known as the composite theory, was developed after
many unsuccessful quests in searching for a unitary genetic mechanism. The
composite theory has multiple sub-theories. The dominance theory and the faster
male theory are the major ones. In this note, I challenge the composite theory
and its scientific validity. By declaring Haldane's rule as a composite
phenomenon caused by multiple mechanisms coincidentally/collectively, the
composite theory becomes a self-fulfilling prophecy and untestable. I believe
that the composite theory is an ad hoc hypothesis that lacks falsifiability,
refutability and testability that a scientific theory requires. It is my belief
that the composite theory does not provide meaningful insights for the study of
speciation and should be abandoned.


  The leaves of the Coriandrum sativum plant, known as cilantro or coriander,
are widely used in many cuisines around the world. However, far from being a
benign culinary herb, cilantro can be polarizing---many people love it while
others claim that it tastes or smells foul, often like soap or dirt. This soapy
or pungent aroma is largely attributed to several aldehydes present in
cilantro. Cilantro preference is suspected to have a genetic component, yet to
date nothing is known about specific mechanisms. Here we present the results of
a genome-wide association study among 14,604 participants of European ancestry
who reported whether cilantro tasted soapy, with replication in a distinct set
of 11,851 participants who declared whether they liked cilantro. We find a
single nucleotide polymorphism (SNP) significantly associated with soapy-taste
detection that is confirmed in the cilantro preference group. This SNP,
rs72921001, (p=6.4e-9, odds ratio 0.81 per A allele) lies within a cluster of
olfactory receptor genes on chromosome 11. Among these olfactory receptor genes
is OR6A2, which has a high binding specificity for several of the aldehydes
that give cilantro its characteristic odor. We also estimate the heritability
of cilantro soapy-taste detection in our cohort, showing that the heritability
tagged by common SNPs is low, about 0.087. These results confirm that there is
a genetic component to cilantro taste perception and suggest that cilantro
dislike may stem from genetic variants in olfactory receptors. We propose that
OR6A2 may be the olfactory receptor that contributes to the detection of a
soapy smell from cilantro in European populations.


  Coupling the control of expression stochasticity (noise) to the ability of
expression change (plasticity) can alter gene function and influence
adaptation. A number of factors, such as transcription re-initiation, strong
chromatin regulation or genome neighboring organization, underlie this
coupling. However, these factors do not necessarily combine in equivalent ways
and strengths in all genes. Can we identify then alternative architectures that
modulate in distinct ways the linkage of noise and plasticity? Here we first
show that strong chromatin regulation, commonly viewed as a source of coupling,
can lead to plasticity without noise. The nature of this regulation is relevant
too, with plastic but noiseless genes being subjected to general activators
whereas plastic and noisy genes experience more specific repression.
Contrarily, in genes exhibiting poor transcriptional control, it is
translational efficiency what separates noise from plasticity, a pattern
related to transcript length. This additionally implies that genome neighboring
organization -as modifier- appears only effective in highly plastic genes. In
this class, we confirm bidirectional promoters (bipromoters) as a configuration
capable to reduce coupling by abating noise but also reveal an important
trade-off, since bipromoters also decrease plasticity. This presents ultimately
a paradox between intergenic distances and modulation, with short intergenic
distances both associated and disassociated to noise at different plasticity
levels. Balancing the coupling among different types of expression variability
appears as a potential shaping force of genome regulation and organization.
This is reflected in the use of different control strategies at genes with
different sets of functional constraints.


  We report self-consistent Brownian dynamics simulations of a simple
electrostatic model of the selectivity filters (SF) of calcium ion channels.
They reveal regular structure in the conductance and selectivity as functions
of the fixed negative charge Qf at the SF. This structure consists of distinct
regions of high conductance (conduction bands) M0, M1, M2 separated by regions
of zero-conductance (stop-bands). Two of these conduction bands, M1 and M2,
demonstrate high calcium selectivity and prominent anomalous mole fraction
effects and can be identified with the L-type and RyR calcium channels.


  It is a classical result that any finite tree with positively weighted edges,
and without vertices of degree 2, is uniquely determined by the weighted path
distance between each pair of leaves. Moreover, it is possible for a (small)
strict subset $\cl$ of leaf pairs to suffice for reconstructing the tree and
its edge weights, given just the distances between the leaf pairs in $\cl$. It
is known that any set $\cl$ with this property for a tree in which all interior
vertices have degree 3 must form a {\em cover} for $T$ -- that is, for each
interior vertex $v$ of $T$, $\cl$ must contain a pair of leaves from each pair
of the three components of $T-v$. Here we provide a partial converse of this
result by showing that if a set $\cl$ of leaf pairs forms a cover of a certain
type for such a tree $T$ then $T$ and its edge weights can be uniquely
determined from the distances between the pairs of leaves in $\cl$. Moreover,
there is a polynomial-time algorithm for achieving this reconstruction. The
result establishes a special case of a recent question concerning `triplet
covers', and is relevant to a problem arising in evolutionary genomics.


  One of the key challenges in modeling the dynamics of contagion phenomena is
to understand how the structure of social interactions shapes the time course
of a disease. Complex network theory has provided significant advances in this
context. However, awareness of an epidemic in a population typically yields
behavioral changes that correspond to changes in the network structure on which
the disease evolves. This feedback mechanism has not been investigated in
depth. For example, one would intuitively expect susceptible individuals to
avoid other infecteds. However, doctors treating patients or parents tending
sick children may also increase the amount of contact made with an infecteds,
in an effort to speed up recovery but also exposing themselves to higher risks
of infection. We study the role of these caretaker links in an adaptive network
models where individuals react to a disease by increasing or decreasing the
amount of contact they make with infected individuals. We find that pure
avoidance, with only few caretaker links, is the best strategy for curtailing
an SIS disease in networks that possess a large topological variability. In
more homogeneous networks, disease prevalence is decreased for low
concentrations of caretakers whereas a high prevalence emerges if caretaker
concentration passes a well defined critical value.


  The haloarchaea are aerobic, heterotrophic, photophosphorylating prokaryotes,
whose supposed closest relatives and ancestors, the methanogens, are
CO2-reducing, anaerobic chemolithotrophs. Using two available haloarchaeal
genomes we firstly confirmed the methanogenic ancestry of the group and then
investigated those individual genes in the haloarchaea that differ in their
phylogenetic signal to this relationship. We found that almost half the genes,
about which we can make strong statements, have bacterial ancestry and are
likely a result of multiple horizontal transfer events. Futhermore their
functions specifically relate to the phenotypic changes required for a
chemolithotroph to become a heterotroph. If this phylogenetic relationship is
correct, it implies the development of the haloarchaeal phenotype was among the
most extreme changes in cellular physiology fuelled by horizontal gene
transfer.


  Heat shock response is a stress response to temperature changes and a
consecutive increase in amounts of unfolded proteins. To restore homeostasis,
cells upregulate chaperones facilitating protein folding by means of
transcription factors (TF). We here investigate two heat shock systems: one
characteristic to gram negative bacteria, mediated by transcriptional activator
sigma32 in E. coli, and another characteristic to gram positive bacteria,
mediated by transcriptional repressor HrcA in L. lactis. We construct simple
mathematical model of the two systems focusing on the negative feedbacks, where
free chaperons suppress sigma32 activation in the former, while they activate
HrcA repression in the latter. We demonstrate that both systems, in spite of
the difference at the TF regulation level, are capable of showing very similar
heat shock dynamics. We find that differences in regulation impose distinct
constrains on chaperone-TF binding affinities: the binding constant of free
sigma32 to chaperon DnaK, known to be in 100 nM range, set the lower limit of
amount of free chaperon that the system can sense the change at the heat shock,
while the binding affinity of HrcA to chaperon GroE set the upper limit and
have to be rather large extending into the micromolar range.


  Hox proteins are one of the best studied sets of transcription factors in
developmental biology. They are major determinants for establishing
morphological differences along the anterior-posterior axis of animals and are
generally regarded as highly conserved in function. This view is based on
experiments comparing a few (anterior) Hox proteins, however, the extent to
which central or abdominal Hox proteins share sequence features or functions
remains largely unexplored. To shed light on the origin and functional
divergence of the central Hox proteins, we combine a powerful bioinformatics
tool (CLANS) with a large-scale phylogeny of species. CLANS is used to
differentiate between the various types of central Hox protein sequences, while
the phylogeny provides an evolutionary context to the analysis. The combination
of both enables us to infer the relative timepoint at which a given central Hox
protein type arose. We identify seven distinct central Hox protein sequence
types, only one of which is common to all protostome and deuterostome clades
(Antp/Hox7). Together, these results lead us to suggest reevaluating the
usefulness of the increasingly depicted synteny-based classification scheme
that assumes a one-to-one orthology between protostome and deuterostome central
Hox proteins. Instead, we propose that the use of sequence-based classification
schemes able to resolve the central and posterior Hox proteins provides a more
promising and biologically meaningful alternative to resolving these groups.
This analysis, which provides a unique overview of the Hox protein sequence
types present across protostomes and deuterostomes as well as a relative dating
for the emergence of the central Hox protein types, provides a crucial clue to
illuminate how and when the distinct developmental blueprints for organisms
evolved within the evolutionarily immensely successful bilaterian lineage.


  The scarcity of water characterising drylands forces vegetation to adopt
appropriate survival strategies. Some of these generate water-vegetation
feedback mechanisms that can lead to spatial self-organisation of vegetation,
as it has been shown with models representing plants by a density of biomass,
varying continuously in time and space. However, although plants are usually
quite plastic they also display discrete qualities and stochastic behaviour.
These features may give rise to demographic noise, which in certain cases can
influence the qualitative dynamics of ecosystem models. In the present work we
explore the effects of demographic noise on the resilience of a model semi-arid
ecosystem. We introduce a spatial stochastic eco-hydrological hybrid model in
which plants are modelled as discrete entities subject to stochastic dynamical
rules, while the dynamics of surface and soil water are described by continuous
variables. The model has a deterministic approximation very similar to previous
continuous models of arid and semi-arid ecosystems. By means of numerical
simulations we show that demographic noise can have important effects on the
extinction and recovery dynamics of the system. In particular we find that the
stochastic model escapes extinction under a wide range of conditions for which
the corresponding deterministic approximation predicts absorption into desert
states.


  Influenza virus contains two highly variable envelope glycoproteins,
hemagglutinin (HA) and neuraminidase (NA). The structure and properties of HA,
which is responsible for binding the virus to the cell that is being infected,
change significantly when the virus is transmitted from avian or swine species
to humans. Here we focus on much smaller human individual evolutionary amino
acid mutational changes in NA, which cleaves sialic acid groups and is required
for influenza virus replication. We show that very small amino acid changes can
be monitored very accurately across many Uniprot and NCBI strains using
hydropathicity scales to quantify the roughness of water film packages.
Quantitative sequential analysis is most effective with the differential
hydropathicity scale based on protein self-organized criticality (SOC). NA
exhibits punctuated evolution at the molecular scale, millions of times smaller
than the more familiar species scale, and thousands of times smaller than the
genomic scale. Our analysis shows that large-scale vaccination programs have
been responsible for a very large convergent reduction in influenza severity in
the last century, a reduction which is hidden from short-term studies of
vaccine effectiveness. Hydropathic analysis is capable of interpreting and even
predicting trends of functional changes in mutation prolific viruses.


  Evolution is simultaneously driven by a number of processes such as mutation,
competition and random sampling. Understanding which of these processes is
dominating the collective evolutionary dynamics in dependence on system
properties is a fundamental aim of theoretical research. Recent works
quantitatively studied coevolutionary dynamics of competing species with a
focus on linearly frequency-dependent interactions, derived from a
game-theoretic viewpoint. However, several aspects of evolutionary dynamics,
e.g. limited resources, may induce effectively nonlinear frequency
dependencies. Here we study the impact of nonlinear frequency dependence on
evolutionary dynamics in a model class that covers linear frequency dependence
as a special case. We focus on the simplest non-trivial setting of two
genotypes and analyze the co-action of nonlinear frequency dependence with
asymmetric mutation rates. We find that their co-action may induce novel
metastable states as well as stochastic switching dynamics between them. Our
results reveal how the different mechanisms of mutation, selection and genetic
drift contribute to the dynamics and the emergence of metastable states,
suggesting that multistability is a generic feature in systems with
frequency-dependent fitness.


  We have developed a mathematical model for in-host virus dynamics that
includes spatial chemotaxis and diffusion across a two dimensional surface
representing the vaginal or rectal epithelium at primary HIV infection. A
linear stability analysis of the steady state solutions identified conditions
for Turing instability pattern formation. We have solved the model equations
numerically using parameter values obtained from previous experimental results
for HIV infections. Simulations of the model for this surface show hot spots of
infection. Understanding this localization is an important step in the ability
to correctly model early HIV infection. These spatial variations also have
implications for the development and effectiveness of microbicides against HIV.


  Histone deacetylase (HDAC) and Histone acetyl-transferase (HAT) are enzymes
that influence transcription by selectively deacetylating or acetylating the
(epsilon)-amino groups of lysine located near the amino termini of core histone
proteins. Over expression of HDACs noted in many forms of cancers including
leukemia and breast cancer. HDAC inhibitors have been shown to be potent
inducers of growth arrest, differentiation, and/or apoptotic cell death. There
is a growing interest in the development of histone deacetylase inhibitors as
anti cancer agents. Three known ligands of HDAC-8 were taken and docked. The
best scores were analyzed and structures similar to these ligands were
downloaded using carol and corina databases and docked. Also large databases of
small molecules were computationally screened using molecular docking for hits
that can conformationally and chemically fit to the active site. Molecules
which got high scores for both GoldScore and ChemScore were selected and
compared with the previous results. Those with best results were then taken for
calculating H-bond interactions and close contacts. Bioactivity prediction of
the best ranked ligands was done. Their physicochemical properties were also
analyzed. Four new molecules were identified and suggested for further testing
in the wet lab.


  Currently, we are overwhelmed by a deluge of experimental data, and network
physics has the potential to become an invaluable method to increase our
understanding of large interacting datasets. However, this potential is often
unrealized for two reasons: uncovering the hidden community structure of a
network, known as community detection, is difficult, and further, even if one
has an idea of this community structure, it is not a priori obvious how to
efficiently use this information. Here, to address both of these issues, we,
first, identify optimal community structure of given networks in terms of
modularity by utilizing a recently introduced community detection method.
Second, we develop an approach to use this community information to extract
hidden information from a network. When applied to a protein-protein
interaction network, the proposed method outperforms current state-of-the-art
methods that use only the local information of a network. The method is
generally applicable to networks from many areas.


  Most common SNPs are popularly assumed to be neutral. We here developed novel
methods to examine in animal models and humans whether extreme amount of minor
alleles (MAs) carried by an individual may represent extreme trait values and
common diseases. We analyzed panels of genetic reference populations and
identified the MAs in each panel and the MA content (MAC) that each strain
carried. We also analyzed 21 published GWAS datasets of human diseases and
identified the MAC of each case or control. MAC was nearly linearly linked to
quantitative variations in numerous traits in model organisms, including life
span, tumor susceptibility, learning and memory, sensitivity to alcohol and
anti-psychotic drugs, and two correlated traits poor reproductive fitness and
strong immunity. Similarly, in Europeans or European Americans, enrichment of
MAs of fast but not slow evolutionary rate was linked to autoimmune and
numerous other diseases, including type 2 diabetes, Parkinson's disease,
psychiatric disorders, alcohol and cocaine addictions, cancer, and less life
span. Therefore, both high and low MAC correlated with extreme values in many
traits, indicating stabilizing selection on most MAs. The methods here are
broadly applicable and may help solve the missing heritability problem in
complex traits and diseases.


  Ion channels are specific proteins present in the membranes of living cells.
They control the flow of specific ions through a cell, initiated by an ion
channel's electrochemical gradient. In doing so, they control important
physiological processes such as muscle contraction and neuronal connectivity,
which cannot be properly activated if these channels go haywire, leading to
life-threatening diseases and psychological disorders. Here, we will develop a
generalized microscopic theory of ion selectivity applicable to KcsA, Na$_{\rm
v}$Rh and Ca$_{\rm v}$ (L-type) ion channels. We unambiguously expose why and
how a given ion-channel can be highly selective, and yet has a conductance of
the order of one million ions per second, or higher. We will identify and prove
the correct physico-biochemical mechanisms that are responsible for the high
selectivity of a particular ion in a given ion channel. The above mechanisms
consist of five conditions, which can be directly associated to these
parameters - (i) dehydration energy, (ii) concentration of the "correct" ions
(iii) Coulomb-van-der-Waals attraction, (iv) pore and ionic sizes, and
indirectly to (v) the thermodynamic stability and (vi) the "knock-on" assisted
permeation.


  One of the most salient spatio-temporal patterns in population ecology is the
synchronization of fluctuating local populations across vast spatial extent.
Synchronization of abundance has been widely observed across a range of spatial
scales in relation to rate of dispersal among discrete populations. However,
the dependence of synchrony on patterns of among-patch movement across
heterogeneous landscapes has been largely ignored. Here we consider the
duration of movement between two predator-prey communities connected by weak
dispersal, and its effect on population synchrony. More specifically, we
introduce time delayed dispersal to incorporate the finite transmission time
between discrete populations across a continuous landscape. Reducing the system
to a phase model using weakly connected network theory, it is found that the
time delay is an important factor determining the nature and stability of
phase-locked states. Our analysis predicts enhanced convergence to stable
synchronous fluctuations in general, and a decreased ability of systems to
produce in-phase synchronization dynamics in the presence of delayed dispersal.
These results introduce delayed dispersal as a tool for understanding the
importance of dispersal time across a landscape matrix in affecting
metacommunity dynamics. They further highlight the importance of landscape and
dispersal patterns for predicting the onset of synchrony between weakly-coupled
populations.


  Comparing allele frequencies among populations that differ in environment has
long been a tool for detecting loci involved in local adaptation. However, such
analyses are complicated by an imperfect knowledge of population allele
frequencies and neutral correlations of allele frequencies among populations
due to shared population history and gene flow. Here we develop a set of
methods to robustly test for unusual allele frequency patterns, and
correlations between environmental variables and allele frequencies while
accounting for these complications based on a Bayesian model previously
implemented in the software Bayenv. Using this model, we calculate a set of
`standardized allele frequencies' that allows investigators to apply tests of
their choice to multiple populations, while accounting for sampling and
covariance due to population history. We illustrate this first by showing that
these standardized frequencies can be used to calculate powerful tests to
detect non-parametric correlations with environmental variables, which are also
less prone to spurious results due to outlier populations. We then demonstrate
how these standardized allele frequencies can be used to construct a test to
detect SNPs that deviate strongly from neutral population structure. This test
is conceptually related to FST but should be more powerful as we account for
population history. We also extend the model to next-generation sequencing of
population pools, which is a cost-efficient way to estimate population allele
frequencies, but it implies an additional level of sampling noise. The utility
of these methods is demonstrated in simulations and by re-analyzing human SNP
data from the HGDP populations. An implementation of our method will be
available from http://gcbias.org.


  While modern structural biology has provided us with a rich and diverse
picture of membrane proteins, the biological function of membrane proteins is
often influenced by the mechanical properties of the surrounding lipid bilayer.
Here we explore the relation between the shape of membrane proteins and the
cooperative function of membrane proteins induced by membrane-mediated elastic
interactions. For the experimental model system of mechanosensitive ion
channels we find that the sign and strength of elastic interactions depend on
the protein shape, yielding distinct cooperative gating curves for distinct
protein orientations. Our approach predicts how directional elastic
interactions affect the molecular structure, organization, and biological
function of proteins in crowded membranes.


  Among aquaglyceroporins that transport both water and glycerol across the
cell membrane, Escherichia coli glycerol uptake facilitator (GlpF) is the most
thoroughly studied. However, one question remains: Does glycerol modulate water
permeation? This study answers this fundamental question by determining the
chemical-potential profile of glycerol along the permeation path through GlpF's
conducting pore. There is a deep well near the Asn-Pro-Ala (NPA) motifs
(dissociation constant 14 microM) and a barrier near the selectivity filter
(10.1 kcal/mol above the well bottom). This profile owes its existence to
GlpF's perfect steric arrangement: The glycerol-protein van der Waals
interactions are attractive near the NPA but repulsive elsewhere in the
conducting pore. In light of the single-file nature of waters and glycerols
lining up in GlpF's amphipathic pore, it leads to the following conclusion:
Glycerol modulates water permeation in the microM range. At mM concentrations,
GlpF is glycerol-saturated and a glycerol dwelling in the well occludes the
conducting pore. Therefore, water permeation is fully correlated to glycerol
dissociation that has an Arrhenius activation barrier of 6.5 kcal/mol.
Validation of this theory is based on the existent in vitro data, some of which
have not been given the proper attention they deserved: The Arrhenius
activation barriers were found to be 7 kcal/mol for water permeation and 9.6
kcal/mol for glycerol permeation; The presence of up to 100 mM glycerol did not
affect the kinetics of water transport with very low permeability, in apparent
contradiction with the existent theories that predicted high permeability (0 M
glycerol).


  Swarming behaviors in animals have been extensively studied due to their
implications for the evolution of cooperation, social cognition, and
predator-prey dynamics. An important goal of these studies is discerning which
evolutionary pressures favor the formation of swarms. One hypothesis is that
swarms arise because the presence of multiple moving prey in swarms causes
confusion for attacking predators, but it remains unclear how important this
selective force is. Using an evolutionary model of a predator-prey system, we
show that predator confusion provides a sufficient selection pressure to evolve
swarming behavior in prey. Furthermore, we demonstrate that the evolutionary
effect of predator confusion on prey could in turn exert pressure on the
structure of the predator's visual field, favoring the frontally oriented,
high-resolution visual systems commonly observed in predators that feed on
swarming animals. Finally, we provide evidence that when prey evolve swarming
in response to predator confusion, there is a change in the shape of the
functional response curve describing the predator's consumption rate as prey
density increases. Thus, we show that a relatively simple perceptual
constraint--predator confusion--could have pervasive evolutionary effects on
prey behavior, predator sensory mechanisms, and the ecological interactions
between predators and prey.


  A recursion for the joint moments of the external branch lengths for
coalescents with multiple collisions ($\Lambda$-coalescents) is provided. This
recursion is used to derive asymptotic results as the sample size $n$ tends to
infinity for the joint moments of the external branch lengths and for the
moments of the total external branch length of the Bolthausen-Sznitman
coalescent. These asymptotic results are based on a differential equation
approach, which is as well useful to obtain exact solutions for the joint
moments of the external branch lengths for the Bolthausen-Sznitman coalescent.
The results for example show that the lengths of two randomly chosen external
branches are positively correlated for the Bolthausen--Sznitman coalescent,
whereas they are negatively correlated for the Kingman coalescent provided that
$n\ge 4$.


  How natural selection acts to limit the proliferation of transposable
elements (TEs) in genomes has been of interest to evolutionary biologists for
many years. To describe TE dynamics in populations, many previous studies have
used models of transposition-selection equilibrium that rely on the assumption
of a constant rate of transposition. However, since TE invasions are known to
happen in bursts through time, this assumption may not be reasonable in natural
populations. Here we propose a test of neutrality for TE insertions that does
not rely on the assumption of a constant transposition rate. We consider the
case of TE insertions that have been ascertained from a single haploid
reference genome sequence and have subsequently had their allele frequency
estimated in a population sample. By conditioning on the age of an individual
TE insertion (using information contained in the number of substitutions that
have occurred within the TE sequence since insertion), we determine the
probability distribution for the insertion allele frequency in a population
sample under neutrality. Taking models of varying population size into account,
we then evaluate predictions of our model against allele frequency data from
190 retrotransposon insertions sampled from North American and African
populations of Drosophila melanogaster. Using this non-equilibrium model, we
are able to explain about 80% of the variance in TE insertion allele
frequencies based on age alone. Controlling both for nonequilibrium dynamics of
transposition and host demography, we provide evidence for negative selection
acting against most TEs as well as for positive selection acting on a small
subset of TEs. Our work establishes a new framework for the analysis of the
evolutionary forces governing large insertion mutations like TEs, gene
duplications or other copy number variants.


  There has been a considerable effort to understand and quantify the spatial
distribution of species across different ecosystems. Relative species abundance
(RSA), beta diversity and species area relationship (SAR) are among the most
used macroecological measures to characterize plants communities in forests. In
this article we introduce a simple phenomenological model based on Poisson
cluster processes which allows us to exactly link RSA and beta diversity to
SAR. The framework is spatially explicit and accounts for the spatial
aggregation of conspecific individuals. Under the simplifying assumption of
neutral theory, we derive an analytical expression for the SAR which reproduces
tri-phasic behavior as sample area increases from local to continental scales,
explaining how the tri-phasic behavior can be understood in terms of simple
geometric arguments. We also find an expression for the endemic area
relationship (EAR) and for the scaling of the RSA.


  In a single molecule assay, the motion of a molecular motor is often inferred
from measuring the stochastic trajectory of a large probe particle attached to
it. We discuss a simple model for this generic set-up taking into account
explicitly the elastic coupling between probe and motor. The combined dynamics
consists of discrete steps of the motor and continuous Brownian motion of the
probe. Motivated by recent experiments on the F1-ATPase, we investigate three
types of efficiencies both in simulations and a Gaussian approximation.
Overall, we obtain good quantitative agreement with the experimental data. In
particular, we clarify the conditions under which one of these efficiencies
becomes larger than one.


  Liu et al recently proposed a minimum number of driver nodes needed to obtain
full structural controllability over a directed network. Driver nodes are
unmatched nodes, from which there are directed paths to all matched nodes.
Their most important assertion is that a system's controllability is to a great
extent encoded by the underlying network's degree distribution, $P(k_{in},
k_{out})$. Is the controllability of a network decided almost completely by the
immediate neighbourhood of a node, while, even slightly distant nodes play no
role at all? Motivated by the above question, in this communication, we argue
that an effective understanding of controllability in directed networks can be
reached using distance based measures of closeness centrality and betweenness
centrality and may not require the knowledge of local connectivity measures
like in-degree and out-degree.


  Epigenetic mechanisms of silencing via heritable chromatin modifications play
a major role in gene regulation and cell fate specification. We consider a
model of epigenetic chromatin silencing in budding yeast and study the
bifurcation diagram and characterize the bistable and the monostable regimes.
The main focus of this paper is to examine how the perturbations altering the
activity of histone modifying enzymes affect the epigenetic states. We analyze
the implications of having the total number of silencing proteins given by the
sum of proteins bound to the nucleosomes and the ones available in the ambient
to be constant. This constraint couples different regions of chromatin through
the shared reservoir of ambient silencing proteins. We show that the response
of the system to perturbations depends dramatically on the titration effect
caused by the above constraint. In particular, for a certain range of overall
abundance of silencing proteins, the hysteresis loop changes qualitatively with
certain jump replaced by continuous merger of different states. In addition, we
find a nonmonotonic dependence of gene expression on the rate of histone
deacetylation activity of Sir2. We discuss how these qualitative predictions of
our model could be compared with experimental studies of the yeast system under
anti-silencing drugs.


  The immune response to a pathogen has two basic features. The first is the
expansion of a few pathogen-specific cells to form a population large enough to
control the pathogen. The second is the process of differentiation of cells
from an initial naive phenotype to an effector phenotype which controls the
pathogen, and subsequently to a memory phenotype that is maintained and
responsible for long-term protection. The expansion and the differentiation
have been considered largely independently. Changes in cell populations are
typically described using ecologically based ordinary differential equation
models. In contrast, differentiation of single cells is studied within systems
biology and is frequently modeled by considering changes in gene and protein
expression in individual cells. Recent advances in experimental systems biology
make available for the first time data to allow the coupling of population and
high dimensional expression data of immune cells during infections. Here we
describe and develop population-expression models which integrate these two
processes into systems biology on the multicellular level. When translated into
mathematical equations, these models result in non-conservative, non-local
advection-diffusion equations. We describe situations where the
population-expression approach can make correct inference from data while
previous modeling approaches based on common simplifying assumptions would
fail. We also explore how model reduction techniques can be used to build
population-expression models, minimizing the complexity of the model while
keeping the essential features of the system. While we consider problems in
immunology in this paper, we expect population-expression models to be more
broadly applicable.


  Multimodal oncological strategies which combine chemotherapy or radiotherapy
with hyperthermia have a potential of improving the efficacy of the
non-surgical methods of cancer treatment. Hyperthermia engages the heat-shock
response mechanism (HSR), main component of which are heat-shock proteins
(HSP). Cancer cells have already partially activated HSR, thereby, hyperthermia
may be more toxic to them relative to normal cells. On the other hand, HSR
triggers thermotolerance, i.e. hyperthermia treated cells show an impairment in
their susceptibility to a subsequent heat-induced stress. This poses questions
about efficacy and optimal strategy of the anti-cancer therapy combined with
hyperthermia treatment.
  To address these questions, we adapt our previous HSR model and propose its
stochastic extension. We formalise the notion of a HSP-induced thermotolerance.
Next, we estimate the intensity and the duration of the thermotolerance.
Finally, we quantify the effect of a multimodal therapy based on hyperthermia
and a cytotoxic effect of bortezomib, a clinically approved proteasome
inhibitor. Consequently, we propose an optimal strategy for combining
hyperthermia and proteasome inhibition modalities.
  In summary, by a proof of concept mathematical analysis of HSR we are able to
support the common belief that the combination of cancer treatment strategies
increases therapy efficacy. thermotolerance.


  Proper lung functioning requires not only a correct structure of the
conducting airway tree, but also the simultaneous development of smooth muscles
and vasculature. Lung branching morphogenesis is strongly stereotyped and
involves the recursive use of only three modes of branching. We have previously
shown that the experimentally described interactions between Fibroblast growth
factor (FGF)10, Sonic hedgehog (SHH) and Patched (Ptc) can give rise to a
Turing mechanism that not only reproduces the experimentally observed wildtype
branching pattern but also, in part counterintuitive, patterns in mutant mice.
Here we show that, even though many proteins affect smooth muscle formation and
the expression of Vegfa, an inducer of blood vessel formation, it is sufficient
to add FGF9 to the FGF10/SHH/Ptc module to successfully predict simultaneously
the emergence of smooth muscles in the clefts between growing lung buds, and
Vegfa expression in the distal sub-epithelial mesenchyme. Our model reproduces
the phenotype of both wildtype and relevant mutant mice, as well as the results
of most culture conditions described in the literature.


  Using a network representation for real soil samples and mathematical models
for microbial spread, we show that the structural heterogeneity of the soil
habitat may have a very significant influence on the size of microbial
invasions of the soil pore space. In particular, neglecting the soil structural
heterogeneity may lead to a substantial underestimation of microbial invasion.
Such effects are explained in terms of a crucial interplay between
heterogeneity in microbial spread and heterogeneity in the topology of soil
networks. The main influence of network topology on invasion is linked to the
existence of long channels in soil networks that may act as bridges for
transmission of microorganisms between distant parts of soil.