  Energy saving mechanisms are ubiquitous in nature. Aerodynamic and
hydrodynamic drafting, vortice uplift, Bernoulli suction, thermoregulatory
coupling, path following, physical hooks, synchronization, and cooperation are
only some of the better-known examples. While drafting mechanisms also appear
in non-biological systems such as sedimentation and particle vortices, the
broad spectrum of these mechanisms appears more diversely in biological systems
including bacteria, spermatozoa, various aquatic species, birds, land animals,
semi-fluid dwellers like turtle hatchlings, as well as human systems. We
present the thermodynamic framework for energy saving mechanisms, and we review
evidence in favor of the variation range hypothesis. This hypothesis posits
that, as an evolutionary process, the variation range between strongest and
weakest group members converges on the equivalent energy saving quantity that
is generated by the energy saving mechanism. We also review self-organized
structures that emerge due to energy saving mechanisms, including convective
processes that can be observed in many systems over both short and long time
scales, as well as high collective output processes in which a form of
collective position locking occurs.


  We consider a population dynamics model coupling cell growth to a diffusion
in the space of metabolic phenotypes as it can be obtained from realistic
constraints-based modelling. In the asymptotic regime of slow diffusion, that
coincides with the relevant experimental range, the resulting non-linear
Fokker-Planck equation is solved for the steady state in the WKB approximation
that maps it into the ground state of a quantum particle in an Airy potential
plus a centrifugal term. We retrieve scaling laws for growth rate fluctuations
and time response with respect to the distance from the maximum growth rate
suggesting that suboptimal populations can have a faster response to
perturbations.


  The Rosids is one of the largest groups of flowering plants, with 140
families and ~70,000 species. Previous phylogenetic studies of the rosids have
primarily utilized organelle genes that likely differ in evolutionary histories
from nuclear genes. To better understand the evolutionary history of rosids, it
is necessary to investigate their phylogenetic relationships using nuclear
genes. Here, we employed large-scale phylogenomic datasets composed of nuclear
genes, including 891 clusters of putative orthologous genes. Combined with
comprehensive taxon sampling covering 63 species representing 14 out of the 17
orders, we reconstructed the rosids phylogeny with coalescence and
concatenation methods, yielding similar tree topologies from all datasets.
However, these topologies did not agree on the placement of Zygophyllales.
Through comprehensive analyses, we found that missing data and gene tree
heterogeneity were potential factors that may mislead concatenation methods, in
particular, large amounts of missing data under high gene tree heterogeneity.
Our results provided new insights into the deep phylogenetic relationships of
the rosids, and demonstrated that coalescence methods may effectively resolve
the phylogenetic relationships of the rosids with missing data under high gene
tree heterogeneity.


  Known graphical conditions for the generic or global convergence to
equilibria of the dynamical system arising from a reaction network are shown to
be invariant under the so-called successive removal of intermediates, a
systematic procedure to simplify the network, making the graphical conditions
easier to check.


  We address the problem of determining the stationary distribution of the
multi-allelic, neutral-evolution Wright-Fisher model in the diffusion limit. A
full solution to this problem for an arbitrary K x K mutation rate matrix
involves solving for the stationary solution of a forward Kolmogorov equation
over a (K - 1)-dimensional simplex, and remains intractable. In most practical
situations mutations rates are slow on the scale of the diffusion limit and the
solution is heavily concentrated on the corners and edges of the simplex. In
this paper we present a practical approximate solution for slow mutation rates
in the form of a set of line densities along the edges of the simplex. The
method of solution relies on parameterising the general non-reversible rate
matrix as the sum of a reversible part and a set of (K - 1)(K - 2)/2
independent terms corresponding to fluxes of probability along closed paths
around faces of the simplex. The solution is potentially a first step in
estimat- ing non-reversible evolutionary rate matrices from observed allele
frequency spectra.


  New experimental results on bacterial growth inspire a novel top-down
approach to study cell metabolism, combining mass balance and proteomic
constraints to extend and complement Flux Balance Analysis. We introduce here
Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic
costs associated to growth are accounted for in an effective way through a
single additional genome-wide constraint. Its roots lie in the experimentally
observed pattern of proteome allocation for metabolic functions, allowing to
bridge regulation and metabolism in a transparent way under the principle of
growth-rate maximization. We provide a simple method to solve CAFBA efficiently
and propose an "ensemble averaging" procedure to account for unknown protein
costs. Applying this approach to modeling E. coli metabolism, we find that, as
the growth rate increases, CAFBA solutions cross over from respiratory,
growth-yield maximizing states (preferred at slow growth) to fermentative
states with carbon overflow (preferred at fast growth). In addition, CAFBA
allows for quantitatively accurate predictions on the rate of acetate excretion
and growth yield based on only 3 parameters determined by empirical growth
laws.


  Since the arrival of next-generation sequencing technologies the amount of
genetic sequencing data has increased dramatically. This has has fueled an
increase in human genetics research. At the same time, with the recent advent
of technologies in processing large data sets, lot of these technologies are
proving valuable and efficient in analyzing these huge datasets. In this paper
we use some of these technologies to analyze genetic sequencing data of 1000
Genomes Project,produce and evaluate a framework to process the sequencing data
thereof and look into structural variations with respect to population groups.


  This paper suggests that the fundamental haploid-diploid cycle of eukaryotic
sex exploits a rudimentary form of the Baldwin effect. With this explanation
for the basic cycle, the other associated phenomena can be explained as
evolution tuning the amount and frequency of learning experienced by an
organism. Using the well-known NK model of fitness landscapes it is shown that
varying landscape ruggedness varies the benefit of the haploid-diploid cycle,
whether based upon endomitosis or syngamy. The utility of pre-meiotic doubling
and recombination during the cycle are also shown to vary with landscape
ruggedness. This view is suggested as underpinning, rather than contradicting,
many existing explanations for sex.


  We show that, under certain assumptions, the fitness of almost all
quasi-species becomes independent of mutational probabilities and the initial
frequency distributions of the sequences in high dimensional sequence spaces.
This result is the consequence of the concentration of measure on a high
dimensional hypersphere and its extension to Lipschitz functions knows as the
Levy's Lemma. Therefore, evolutionary dynamics almost always yields the same
value for fitness of the quasi-species, independent of the mutational process
and initial conditions, and is quite robust to mutational changes and
fluctuations in initial conditions. Our results naturally extend to any
Lipschitz function whose input parameters are the frequencies of individual
constituents of the quasi-species. This suggests that the functional
capabilities of high dimensional quasi-species are robust to fluctuations in
the mutational probabilities and initial conditions.


  We introduce a model called Host-Pathogen game for studying biological
competitions. Notably, we focus on the invasive dynamics of external agents,
like bacteria, within a host organism. The former are mapped to a population of
defectors that aim to spread in the extracellular medium of the host. In turn,
the latter is composed of cells, mapped to a population of cooperators, that
aim to kill pathogens. The cooperative behavior of cells is fundamental for the
emergence of the living functions of the whole organism, since each one
contributes to a specific set of tasks. So, broadly speaking, their
contribution can be viewed as a form of energy. When bacteria are spatially
close to a cell, the latter can use a fraction of its energy to remove them. On
the other hand, when bacteria survive an attack, they absorb the received
energy, becoming stronger and more resistant to further attacks. In addition,
since bacteria play as defectors, their unique target is to increase their
wealth, without supporting their own kind. As in many living organisms, the
host temperature plays a relevant role in the host-pathogen equilibrium. For
instance, in animals like human beings, a neural mechanism triggers the
increasing of the body temperature in order to activate the immune system.
Here, cooperators succeed once bacteria are completely removed while, in the
opposite scenario, the host undergoes a deep invasive process, like a blood
poisoning. Results of numerical simulations show that the dynamics of the
proposed model allow to reach a variety of states. At a very high level of
abstraction, some of these states seem to be similar to those that can be
observed in some living systems. Therefore, to conclude, we deem that our model
might be exploited for studying further biological phenomena.


  A substantial body of work supports a teleconnection between the El Ni\~no
Southern Oscillation (ENSO) and cholera incidence in Bangladesh. In particular,
high positive anomalies during the winter (Dec-Feb) in Sea Surface Temperatures
(SST) in the Tropical Pacific have been shown to exacerbate the seasonal
outbreak of cholera following the monsoons from Aug to Nov, and climate studies
have indicated a role of regional precipitation over Bangladesh in mediating
this long-distance effect. Thus, the current strong El Ni\~no has the potential
to significantly increase cholera risk this year in Dhaka, Bangladesh, where
the last five years have experienced low seasons of the disease. To examine
this possibility and produce a forecast for the city, we considered two models
for the transmission dynamics of cholera: a statistical model previously
developed for the disease in this region, and a process-based model presented
here that includes the effect of SST anomalies in the force of infection and is
fitted to extensive cholera surveillance record between 1995 and 2010.
Prediction accuracy was evaluated with 'out-of-fit' data from the same
surveillance efforts, by comparing the total number of cholera cases observed
for the season to those predicted by model simulations 8 to 12 months ahead,
starting in January each year. Encouraged by accurate forecasts for the low
risk of cholera for this period, we then generated a prediction for this coming
season. An increase above the third quantile in cholera cases is expected for
the period of Aug - Dec 2016 with 92% and 87% probability respectively for the
two models. This alert warrants the preparedness of the public health system.
We discuss the possible limitations of our approach, including variations in
the impact of El Ni\~no events, and the importance of this large, warm event
for further informing an early-warning system for cholera in Dhaka


  We consider the problem of finding the Perron-Frobenius eigenvector of a
primitive matrix. Dividing each of the rows of the matrix by the sum of the
elements in the row, the resulting new matrix is stochastic. We give a formula
for the Perron-Frobenius eigenvector of the original matrix, in terms of a
realization of the Markov chain defined by the associated stochastic matrix.
This formula is a generalization of the classical formula for the invariant
probability measure of a Markov chain.


  SMISS is a novel web server for protein function prediction. Three different
predictors can be selected for different usage. It integrates different sources
to improve the protein function prediction accuracy, including the query
protein sequence, protein-protein interaction network, gene-gene interaction
network, and the rules mined from protein function associations. SMISS
automatically switch to ab initio protein function prediction based on the
query sequence when there is no homologs in the database. It takes fasta format
sequences as input, and several sequences can submit together without
influencing the computation speed too much. PHP and Perl are two primary
programming language used in the server. The CodeIgniter MVC PHP web framework
and Bootstrap front-end framework are used for building the server. It can be
used in different platforms in standard web browser, such as Windows, Mac OS X,
Linux, and iOS. No plugins are needed for our website. Availability:
http://tulip.rnet.missouri.edu/profunc/.


  The presence of knots has been observed in a small fraction of single-domain
proteins and related to their thermodynamic and kinetic properties. The
exchanging of identical structural elements, typical of domain-swapped
proteins, make such dimers suitable candidates to validate the possibility that
mutual entanglement between chains may play a similar role for protein
complexes. We suggest that such entanglement is captured by the linking number.
This represents, for two closed curves, the number of times that each curve
winds around the other. We show that closing the curves is not necessary, as a
novel parameter $G'$, termed Gaussian entanglement, is strongly correlated with
the linking number. Based on $110$ non redundant domain-swapped dimers, our
analysis evidences a high fraction of chains with a significant intertwining,
that is with $|G'| > 1$. We report that Nature promotes configurations with
negative mutual entanglement and surprisingly, it seems to suppress
intertwining in long protein dimers. Supported by numerical simulations of
dimer dissociation, our results provide a novel topology-based classification
of protein-swapped dimers together with some preliminary evidence of its impact
on their physical and biological properties.


  Geometric constraints impact the formation of a broad range of spatial
networks, from amino acid chains folding to proteins structures to rearranging
particle aggregates. How the network of interactions dynamically self-organizes
in such systems is far from fully understood. Here, we analyze a class of
spatial network formation processes by introducing a mapping from geometric to
graph-theoretic constraints. Combining stochastic and mean field analyses
yields an algebraic scaling law for the extent (graph diameter) of the
resulting networks with system size, in contrast to logarithmic scaling known
for networks without constraints. Intriguingly, the exponent falls between that
of self-avoiding random walks and that of space filling arrangements,
consistent with experimentally observed scaling (of the spatial radius of
gyration) for protein tertiary structures.


  In many virus families, tens to thousands of proteins assemble spontaneously
into a capsid (protein shell) while packaging the genomic nucleic acid. This
review summarizes recent advances in computational modeling of these dynamical
processes. We present an overview of recent technological and algorithmic
developments, which are enabling simulations to describe the large ranges of
length-and time-scales relevant to assembly, under conditions more closely
matched to experiments than in earlier work. We then describe two examples in
which computational modeling has recently provided an important complement to
experiments.


  The interplay between parasites and their hosts is found in all kinds of
species and plays an important role in understanding the principles of
evolution and coevolution. Usually, the different genotypes of hosts and
parasites oscillate in their abundances. The well-established theory of
oscillatory Red Queen dynamics proposes an ongoing change in frequencies of the
different types within each species. So far, it is unclear in which way Red
Queen dynamics persists with more than two types of hosts and parasites. In our
analysis, an arbitrary number of types within two species are examined in a
deterministic framework with constant or changing population size. This general
framework allows for analytical solutions for internal fixed points and their
stability. For more than two species, apparently chaotic dynamics has been
reported. Here we show that even for two species, once more than two types are
considered per species, irregular dynamics in their frequencies can be observed
in the long run. The nature of the dynamics depends strongly on the initial
configuration of the system; the usual regular Red Queen oscillations are only
observed in some parts of the parameter region.


  We propose a dynamical model depicting the interactions between DNA and a
specific binding protein involving long range transmissions. The dynamics rely
on the coupling between Hydrogen bonds formed between DNA and protein and
between the base pairs because they account for site specificity of the
binding. We adopt the Morse potential with coupling terms to construct the
Hamiltonian. This model gives rise to a breather excitation, corresponding to
the DNA bubble formation, which propagates as the carrier of genetic
information. We examine the various kind of possible coupling dynamics and
suggest the model feasibility in depicting the renaturation or hybridization
processes.


  We propose and study an optimal control model for malaria infection described
by system of fractional differential equations. The model is formulated in
terms of the left and right Caputo fractional derivatives. We determine the
necessary conditions for the optimality of the controlled dynamical system.The
forward-backward sweep method with generalized euler scheme is applied to
numerically compute the solutions of the optimality system.


  We report on the synthesis, characterization, stability and pharmacokinetics
of novel iron based contrast agents for magnetic resonance imaging (MRI).
Statistical copolymers combining multiple phosphonic acid groups and
poly(ethylene glycol) (PEG) were synthesized and used as coating agents for 10
nm iron oxide nanocrystals. In vitro, protein corona and stability assays show
that phosphonic acid PEG copolymers outperform all other coating types
examined, including low molecular weight anionic ligands and polymers. In vivo,
the particle pharmacokinetics is investigated by monitoring the MRI signal
intensity from mouse liver, spleen and arteries as a function of the time,
between one minute and seven days after injection. Iron oxide particles coated
with multi-phosphonic acid PEG polymers are shown to have a blood circulation
lifetime of 250 minutes, i.e. 10 to 50 times greater than that of recently
published PEGylated probes and benchmarks. The clearance from the liver takes
in average 2 to 3 days and is independent of the core size, coating and
particle stability. By comparing identical core particles with different
coatings, we are able to determine the optimum conditions for stealth MRI
probes.


  Mays celebrated theoretical work of the 70s contradicted the established
paradigm by demonstrating that complexity leads to instability in biological
systems. Here Mays random-matrix modelling approach is generalized to realistic
large-scale webs of species interactions, be they structured by networks of
competition, mutualism or both. Simple relationships are found to govern these
otherwise intractable models, and control the parameter ranges for which
biological systems are stable and feasible. Our analysis of model and real
empirical networks is only achievable upon introducing a simplifying
Google-matrix reduction scheme, which in the process, yields a practical
ecological eigenvalue stability index. These results provide an understanding
on how network topology, especially connectance, influences species stable
coexistence. Constraints controlling feasibility (positive equilibrium
populations) in these systems, are found more restrictive than those
controlling stability, helping explain the enigma of why many classes of
feasible ecological models are nearly always stable.


  After an introduction to the general topic of models for a given locus of a
diploid population whose quadratic dynamics is determined by a fitness
landscape, we consider more specifically the models that can be treated using
genetic (or train) algebras. In this setup, any quadratic offspring interaction
can produce any type of offspring and after the use of specific changes of
basis, we study the evolution and possible stability of some examples. We also
consider some examples that cannot be treated using the framework of genetic
algebras. Among these are bistochastic matrices.


  We derive some additional results on the Bienyam\'e-Galton-Watson branching
process with $\theta -$linear fractional branching mechanism, as studied in
\cite{Sag}. This includes: the explicit expression of the limit laws in both
the sub-critical cases and the super-critical cases with finite mean, the
long-run behavior of the population size in the critical case, limit laws in
the super-critical cases with infinite mean when the $\theta $-process is
either regular or explosive, results regarding the time to absorption, an
expression of the probability law of the $\theta $-branching mechanism
involving Bell polynomials, the explicit computation of the stochastic
transition matrix of the $\theta -$process, together with its powers.


  Motivation: Histone modifications are among the most important factors that
control gene regulation. Computational methods that predict gene expression
from histone modification signals are highly desirable for understanding their
combinatorial effects in gene regulation. This knowledge can help in developing
'epigenetic drugs' for diseases like cancer. Previous studies for quantifying
the relationship between histone modifications and gene expression levels
either failed to capture combinatorial effects or relied on multiple methods
that separate predictions and combinatorial analysis. This paper develops a
unified discriminative framework using a deep convolutional neural network to
classify gene expression using histone modification data as input. Our system,
called DeepChrome, allows automatic extraction of complex interactions among
important features. To simultaneously visualize the combinatorial interactions
among histone modifications, we propose a novel optimization-based technique
that generates feature pattern maps from the learnt deep model. This provides
an intuitive description of underlying epigenetic mechanisms that regulate
genes. Results: We show that DeepChrome outperforms state-of-the-art models
like Support Vector Machines and Random Forests for gene expression
classification task on 56 different cell-types from REMC database. The output
of our visualization technique not only validates the previous observations but
also allows novel insights about combinatorial interactions among histone
modification marks, some of which have recently been observed by experimental
studies.


  Considerable mechanistic insight has been gained into amyloid aggregation;
however, a large class of non-amyloid protein aggregates are considered
'amorphous,' and in most cases little is known about their mechanisms.
Amorphous aggregation of {\gamma}-crystallins in the eye lens causes a
widespread disease of aging, cataract. We combined simulations and experiments
to study the mechanism of aggregation of two {\gamma}D-crystallin mutants, W42R
and W42Q - the former a congenital cataract mutation, and the latter a mimic of
age-related oxidative damage. We found that formation of an internal disulfide
was necessary and sufficient for aggregation under physiological conditions.
Two-chain all-atom simulations predicted that one non-native disulfide in
particular, between Cys32 and Cys41, was likely to stabilize an unfolding
intermediate prone to intermolecular interactions. Mass spectrometry and
mutagenesis experiments confirmed the presence of this bond in the aggregates
and its necessity for oxidative aggregation under physiological conditions in
vitro. Mining the simulation data linked formation of this disulfide to
extrusion of the N-terminal \b{eta}-hairpin and rearrangement of the native
\b{eta}-sheet topology. Specific binding between the extruded hairpin and a
distal \b{eta}-sheet, in an intermolecular chain reaction similar to domain
swapping, is the most probable mechanism of aggregate propagation.


  This thesis is aimed at studying mutations, understood as trajectories in the
DNA configuration space. An evolutive model of mutations in terms of Levy
flights is proposed. The parameters of the model are estimated by means of data
from the Long-Term Evolution Experiment (LTEE) with {\it E. Coli} bacteria. The
results of simulations on competition of clones, mean fitness, etc are compared
with experimental data. We discuss the qualitative analogy found between the
bacterial mutator phenotype and the cancerous cells. The role of radiation as
source of mutations is analyzed. We focus on the case of Radon's decay in the
lungs in breathing.


  The somatic mutations in the pathways that drive cancer development tend to
be mutually exclusive across tumors, providing a signal for distinguishing
driver mutations from a larger number of random passenger mutations. This
mutual exclusivity signal can be confounded by high and highly variable
mutation rates across a cohort of samples. Current statistical tests for
exclusivity that incorporate both per-gene and per-sample mutational
frequencies are computationally expensive and have limited precision.
  We formulate a weighted exact test for assessing the significance of
mutational exclusivity in an arbitrary number of mutational events. Our test
conditions on the number of samples with a mutation as well as per-event,
per-sample mutation probabilities. We provide a recursive formula to compute
$p$-values for the weighted test exactly as well as a highly accurate and
efficient saddlepoint approximation of the test. We use our test to approximate
a commonly used permutation test for exclusivity that conditions on per-event,
per-sample mutation frequencies. However, our test is more efficient and it
recovers more significant results than the permutation test. We use our
Weighted Exclusivity Test (WExT) software to analyze hundreds of colorectal and
endometrial samples from The Cancer Genome Atlas, which are two cancer types
that often have extremely high mutation rates. On both cancer types, the
weighted test identifies sets of mutually exclusive mutations in cancer genes
with fewer false positives than earlier approaches.


  A certain class of viruses replicates inside a cell if they can enter the
nucleus through one of many small target pores, before being permanently
trapped or degraded. We adopt for viral motion a switching stochastic process
model and we estimate here the probability and the conditional mean first
passage time for a viral particle to attain alive the nucleus. The cell nucleus
is covered with thousands of small absorbing nuclear pores and the minimum
distance between them defines the smallest spatial scale that limits the
efficiency of stochastic simulations. Using the Neuman-Green's function method
to solve the steady-state Fokker-Planck equation, we derive asymptotic formula
for the probability and mean arrival time to a small window for various pores'
distributions, that agree with stochastic simulations. These formulas reveal
how key geometrical parameters defines the cytoplasmic stage of viral
infection.


  There has been a large number of reported cases of the occurrence of Zika in
different countries in 2016 and it is necessary to develop an early warning
system to initiate preventive campaigns against the disease. A potential early
warning system based on the rise in ocean temperature of the Pacific Ni\~no
Index is proposed. The efficacy is verified using data for the outbreak in
Colombia as obtained from Google Trends.


  The alignment of biological networks has the potential to teach us as much
about biology and disease as has sequence alignment. Sequence alignment can be
optimally solved in polynomial time. In contrast, network alignment is
$NP$-hard, meaning optimal solutions are impossible to find, and the quality of
found alignments depend strongly upon the algorithm used to create them. Every
network alignment algorithm consists of two orthogonal components: first, an
objective function or measure $M$ that is used to evaluate the quality of any
proposed alignment, and second, a search algorithm used to explore the
exponentially large set of possible alignments in an effort to find "good" ones
according to $M$. Objective functions fall into many categories, including
biological measures such as sequence similarity, as well as topological
measures like graphlet similarity and edge coverage (possibly weighted).
Algorithms to search the space of all possible alignments can be deterministic
or stochastic, and many possibilities have been offered over the past decade.
In this paper we introduce a new stochastic search algorithm called SANA:
Simulated Annealing Network Aligner. We test it on several popular objective
functions and demonstrate that it almost universally optimizes each one
significantly better than existing search algorithms. Finally, we compare
several topological objective functions using SANA. Software available at
http://sana.ics.uci.edu.


  The equilibrium structure of a Dinucleosome is studied using an elastic model
that takes into account the force and torque balance conditions. Using the
proper boundary conditions, it is found that the conformational energy of the
problem does not depend on the length of the linker DNA. In addition it is
shown that the two histone octamers are almost perpendicular to each other and
the linker DNA in short lengths is almost straight. These findings could shed
some light on the role of DNA elasticity in the chromatin structure.


  Phytoplankton photosynthesis in most natural waters is often inhibited by
ultraviolet (UV) and intense photosynthetically available radiation (PAR) but
the effects on ocean productivity have received little consideration aside from
polar areas subject to periodic enhanced UV-B due to depletion of stratospheric
ozone. Here we consider responses in the temperate and tropical mid-ocean
regions typically dominated by picophytoplankton including the prokaryotic
lineages, $\textit{Prochlorococcus}$ and $\textit{Synechococcus}$. Spectral
models of photosynthetic response for each lineage were constructed using model
strains cultured at different growth irradiances and temperatures. In the
model, inhibition becomes more severe once exposure exceeds a threshold (Emax)
related to repair capacity. Model parameters are presented for
$\textit{Prochlorococcus}$ adding to those previously presented for
$\textit{Synechococcus}$. The models were applied to the estimation of mid-day,
water-column photosynthesis based on an atmospheric model of spectral incident
radiation, satellite-derived spectral water transparency and sea surface
temperature. Comparing predictions integrated over the water column including
versus excluding inhibition, production was 7-28% lower due to inhibition
depending on strain and site conditions. Inhibition was consistently greater
for $\textit{Prochlorococcus}$ compared to two strains of
$\textit{Synechococcus}$. Considering only the surface mixed layer, production
was inhibited 7-73%. Weighted by lineage abundance and daily PAR exposure,
inhibition of full water column production averages around 20% for the modeled
region of the Pacific with two-thirds of the inhibition due to UV. The results
suggest more consideration is needed on inhibition effects, especially due to
UV, which are largely excluded from present models of global phytoplankton
production.


  Quantitative characterization of bacterial growth has gathered substantial
attention since Monod's pioneering study. Theoretical and experimental work has
uncovered several laws for describing the log growth phase, in which the number
of cells grows exponentially. However, microorganism growth also exhibits lag,
stationary, and death phases under starvation conditions, in which cell growth
is highly suppressed, while quantitative laws or theories for such phases are
underdeveloped. In fact, models commonly adopted for the log phase that consist
of autocatalytic chemical components, including ribosomes, can only show
exponential growth or decay in a population, and phases that halt growth are
not realized. Here, we propose a simple, coarse-grained cell model that
includes inhibitor molecule species in addition to the autocatalytic active
protein. The inhibitor forms a complex with active proteins to suppress the
catalytic process. Depending on the nutrient condition, the model exhibits the
typical transition among the lag, log, stationary, and death phases.
Furthermore, the lag time needed for growth recovery after starvation follows
the square root of the starvation time and is inverse to the maximal growth
rate, in agreement with experimental observations. Moreover, the distribution
of lag time among cells shows an exponential tail, also consistent with
experiments. Our theory further predicts strong dependence of lag time upon the
speed of substrate depletion, which should be examined experimentally. The
present model and theoretical analysis provide universal growth laws beyond the
log phase, offering insight into how cells halt growth without entering the
death phase.


  The HP model of protein folding, where the chain exists in a free medium, is
investigated using a parallel Monte Carlo scheme based upon Wang-Landau
sampling. Expanding on the work of Wust and Landau by introducing a lesser
known replica -exchange scheme between individual Wang- Landau samplers, the
problem of dynamical trapping (spiking in the density of states) was avoided
and an enhancement in the efficiency of traversing configuration space was
obtained. Highlighting dynamical trapping as an issue for lattice polymer
simulations for increasing lengths is explicitly done here for the first time.
The 1/t scheme is also integrated within this sophisticated Monte Carlo
methodology. A trial move set was developed which includes pull, bond
re-bridging, pivot, kink-flip and a newly invented and implemented fragment
random walk move which allowed rapid exploration of high and low temperature
configurations. A native state search was conducted leading to the attainment
of the native states of the benchmark sequences of 2D50 (-21), 2D60 (-36)and
2D64 (-42), whilst attaining minimum energies close to the native state for
2D85(-52 NATIVE= -53), 2D100a (-47 NATIVE= -48) and 2D100b (-49 NATIVE=-50).
Thermodynamic observables such as the specific heat, internal energy, entropy
and free energy were computed for 2D benchmark sequences at varying
temperatures and folding and unfolding behaviour was investigated. Lattice
polymers with monomeric hydrophobic structure were also studied in the same
manner with the recording of minimum energy values and thermodynamic behaviour.
The native results for the benchmark sequences and lattice polymers were
compared with varying computational methods.


  Real-world optimisation problems are often dynamic. Previously good solutions
must be updated or replaced due to changes in objectives and constraints. It is
often claimed that evolutionary algorithms are particularly suitable for
dynamic optimisation because a large population can contain different solutions
that may be useful in the future. However, rigorous theoretical demonstrations
for how populations in dynamic optimisation can be essential are sparse and
restricted to special cases.
  This paper provides theoretical explanations of how populations can be
essential in evolutionary dynamic optimisation in a general and natural
setting. We describe a natural class of dynamic optimisation problems where a
sufficiently large population is necessary to keep track of moving optima
reliably. We establish a relationship between the population-size and the
probability that the algorithm loses track of the optimum.


  The ability to directly record human face-to-face interactions increasingly
enables the development of detailed data-driven models for the spread of
directly transmitted infectious diseases at the scale of individuals. Complete
coverage of the contacts occurring in a population is however generally
unattainable, due for instance to limited participation rates or experimental
constraints in spatial coverage. Here, we study the impact of spatially
constrained sampling on our ability to estimate the epidemic risk in a
population using such detailed data-driven models. The epidemic risk is
quantified by the epidemic threshold of the
susceptible-infectious-recovered-susceptible model for the propagation of
communicable diseases, i.e. the critical value of disease transmissibility
above which the disease turns endemic. We verify for both synthetic and
empirical data of human interactions that the use of incomplete data sets due
to spatial sampling leads to the underestimation of the epidemic risk. The bias
is however smaller than the one obtained by uniformly sampling the same
fraction of contacts: it depends nonlinearly on the fraction of contacts that
are recorded and becomes negligible if this fraction is large enough. Moreover,
it depends on the interplay between the timescales of population and spreading
dynamics.


  The main chain dihedral angles play an important role to decide the protein
conformation. The native states of a protein can be regard as an ensemble of a
lot of similar conformations, in different conformations the main chain
dihedral angles vary in a certain range. Each dihedral angle value can be
described as a distribution, but only using the distribution can't describe the
real conformation space. The reason is that the dihedral angle has correlation
with others, especially the neighbor dihedral angles in primary sequence. In
our study we analysis extensive molecular dynamics (MD) simulation trajectories
of eleven proteins with different sizes and folds, we found that in stable
second structure the correlations only exist between the dihedrals near to each
other in primary sequence, long range correlations are rare. But in unstable
structures (loop) long range correlations exist. Further we observed some
characteristics of the short range correlations in different second structures
({\alpha}-helix, {\beta}-sheet) and we found that we can approximate good high
order dihedral angle distribution good only use three order distribution in
stable second structure which illustrates that high order correlations (over
three order) is small in stable second structure.


  Alchemical theory is emerging as a promising tool in the context of molecular
dynamics simulations for drug discovery projects. In this theoretical
contribution, I revisit the statistical mechanics foundation of non covalent
interactions in drug-receptor systems, providing a unifying treatment that
encompasses the most important variants in the alchemical approaches, from the
seminal Double Annihilation Method by Jorgensen and Ravimohan [W.L. Jorgensen
and C. Ravimohan, J. Chem. Phys. 83,3050, 1985], to the Gilson's Double
Decoupling Method [M. K. Gilson and J. A. Given and B. L. Bush and J. A.
McCammon, Biophys. J. 72, 1047 1997] and the Deng and Roux alchemical theory
[Y. Deng and B. Roux, J. Chem. Theory Comput., 2, 1255 2006]. Connections and
differences between the various alchemical approaches are highlighted and
discussed, and finally placed into the broader context of nonequilibrium
thermodynamics.


  One successful model of interacting biological systems is the Boolean
network. The dynamics of a Boolean network, controlled with Boolean functions,
is usually considered to be a Markovian (memory-less) process. However, both
self organizing features of biological phenomena and their intelligent nature
should raise some doubt about ignoring the history of their time evolution.
Here, we extend the Boolean network Markovian approach: we involve the effect
of memory on the dynamics. This can be explored by modifying Boolean functions
into non-Markovian functions, for example, by investigating the usual
non-Markovian threshold function, - one of the most applied Boolean functions.
By applying the non-Markovian threshold function on the dynamical process of a
cell cycle network, we discover a power law memory with a more robust dynamics
than the Markovian dynamics.


  Tuna are large pelagic fish whose populations are close to panmixia. In
addition, they are threatened species, so it is important for the maintenance
and monitoring of genetic diversity that genetic information at a genome level
be obtained. Here we report the draft assembly of the southern bluefin tuna
genome and the collection of genome-wide sequence data for five other tuna
species. We sampled five tuna species of the genus Thunnus, the northern and
southern bluefin, yellowfin, albacore, and bigeye, as well as the skipjack
(Katsuwonis pelamis), a tuna-like species. Genome assembly was facilitated at
k-mer=25 while k-mer=51 generated assembly artefacts. The estimated size of the
southern bluefin tuna genome was 795 Mb. We assembled two southern bluefin tuna
individuals independently using both paired end and mate pair sequence. This
resulted in scaffolds with N50>174,000 bp and maximum scaffold lengths>1.4 Mb.
Our estimate of the size of the assembled genome was the scaffolded sequences
in common to both assemblies, which amounted to 721 Mb of the 795 Mb of the
southern bluefin tuna genome sequence. Using BLAST, there were matches between
13,039 of 14,341 (91%) refseq mRNA of the zebrafish Danio rerio to the tuna
assembly indicating that most of a generic fish transcriptome was covered by
the assembly.


  The simultaneous expression of the hunchback gene in the numerous nuclei of
the developing fly embryo gives us a unique opportunity to study how
transcription is regulated in living organisms. A recently developed MS2-MCP
technique for imaging nascent messenger RNA in living Drosophila embryos allows
us to quantify the dynamics of the developmental transcription process. The
initial measurement of the morphogens by the hunchback promoter takes place
during very short cell cycles, not only giving each nucleus little time for a
precise readout, but also resulting in short time traces of transcription.
Additionally, the relationship between the measured signal and the promoter
state depends on the molecular design of the reporting probe. We develop an
analysis approach based on tailor made autocorrelation functions that overcomes
the short trace problems and quantifies the dynamics of transcription
initiation. Based on live imaging data, we identify signatures of bursty
transcription initiation from the hunchback promoter. We show that the
precision of the expression of the hunchback gene to measure its position along
the anterior-posterior axis is low both at the boundary and in the anterior
even at cycle 13, suggesting additional post-transcriptional averaging
mechanisms to provide the precision observed in fixed embryos.


  During mRNA translation, several ribosomes attach to the same mRNA molecule
simultaneously translating it into a protein. This pipelining increases the
protein production rate. A natural and important question is what ribosomal
density maximizes the protein production rate. Using mathematical models of
ribosome flow along both a linear and a circular mRNA molecule we prove that
typically the steady-state production rate is maximized when the ribosomal
density is one half of the maximal possible density. We discuss the
implications of our results to endogenous genes under natural cellular
conditions and also to synthetic biology.


  Untargeted metabolomic studies are revealing large numbers of naturally
occurring metabolites that cannot be characterized because their chemical
structures and MS/MS spectra are not available in databases. Here we present
iMet, a computational tool based on experimental tandem mass spectrometry that
could potentially allow the annotation of metabolites not discovered
previously. iMet uses MS/MS spectra to identify metabolites structurally
similar to an unknown metabolite, and gives a net atomic addition or removal
that converts the known metabolite into the unknown one. We validate the
algorithm with 148 metabolites, and show that for 89% of them at least one of
the top four matches identified by iMet enables the proper annotation of the
unknown metabolite. iMet is freely available at http://imet.seeslab.net.


  Chitin and protein are two main building blocks for many natural
biomaterials. The interaction between chitin and protein critically determines
the properties of the composite biological materials. As living organisms
usually encounter complex ambient conditions like water, pH and ions are
critical factors towards the structural integrity of biomaterials. It is
therefore essential to study the chitin-protein interface under different
environmental conditions. Here, an atomistic model consisting of a chitin
substrate and a protein filament is constructed, which is regarded as a
representative of the chitin-protein interface existing in many chitin-based
biomaterials. Based on this model, the mechanical properties of chitin-protein
interface under different moisture and pH values are investigated through
molecular dynamics simulations. The results reveal a weakening effect of water
towards the chitin-protein interface, as well as acidity, i.e. the protonated
protein forms a stronger adhesion to chitin than that in the alkaline
environment. In addition, the effect from side-chain of protein is studied and
it is found that certain kinds of amino acid can form hydrophobic connections
to chitin surface, which means that these peptides partly dodge the weakening
effect of water. Our observation indicates that terminuses and side-chains in
protein are of importance in forming interfacial hydrogen bonds. From our full
atomistic models, we can observe some molecular mechanisms about how protein
interacts with chitin in different conditions, which may spotlight the
engineering on biomaterials with similar interfaces.


  The Sheldon spectrum describes a remarkable regularity in aquatic ecosystems:
the biomass density as a function of logarithmic body mass is approximately
constant over many orders of magnitude. While size-spectrum models have
explained this phenomenon for assemblages of multicellular organisms, this
paper introduces a species-resolved size-spectrum model to explain the
phenomenon in unicellular plankton. A Sheldon spectrum spanning the cell-size
range of unicellular plankton necessarily consists of a large number of
coexisting species covering a wide range of characteristic sizes. The
coexistence of many phytoplankton species feeding on a small number of
resources is known as the Paradox of the Plankton. Our model resolves the
paradox by showing that coexistence is facilitated by the allometric scaling of
four physiological rates. Two of the allometries have empirical support, the
remaining two emerge from predator-prey interactions exactly when the
abundances follow a Sheldon spectrum. Our plankton model is a scale-invariant
trait-based size-spectrum model: it describes the abundance of phyto- and
zooplankton cells as a function of both size and species trait (the maximal
size before cell division). It incorporates growth due to resource consumption
and predation on smaller cells, death due to predation, and a flexible cell
division process. We give analytic solutions at steady state for both the
within-species size distributions and the relative abundances across species.


  In this paper we extend the familiar continuum electrostatic model with a
perturbation to the usual macroscopic boundary condition. The perturbation is
based on the mean spherical approximation (MSA), to derive a multiscale
hydration-shell boundary condition (HSBC). We show that the HSBC/MSA model
reproduces MSA predictions for Born ions in a variety of polar solvents,
including both protic and aprotic solvents. Importantly, the HSBC/MSA model
predicts not only solvation free energies accurately but also solvation
entropies, which standard continuum electrostatic models fail to predict. The
HSBC/MSA model depends only on the normal electric field at the dielectric
boundary, similar to our recent development of an HSBC model for charge-sign
hydration asymmetry, and the reformulation of the MSA as a boundary condition
enables its straightforward application to complex molecules such as proteins.


  The rock-paper-scissors game is a paradigmatic model for biodiversity, with
applications ranging from microbial populations to human societies. Research
has shown, however, that mobility jeopardizes biodiversity by promoting the
formation of spiral waves, especially if there is no conservation law in place
for the total number of competing players. Firstly, we show that even if such a
conservation law applies, mobility still jeopardizes biodiversity in the
spatial rock-paper-scissors game if only a small fraction of links of the
square lattice is randomly rewired. Secondly, we show that zealots are very
effective in taming the amplitude of oscillations that emerge due to mobility
and/or interaction randomness, and this regardless of whether the later is
quenched or annealed. While even a tiny fraction of zealots brings significant
benefits, at 5\% occupancy zealots practically destroy all oscillations
regardless of the intensity of mobility, and regardless of the type and
strength of randomness in the interaction structure. Interestingly, by annealed
randomness the impact of zealots is qualitatively the same as by mobility,
which highlights that fast diffusion does not necessarily destroy the
coexistence of species, and that zealotry thus helps to recover the stable
mean-field solution. Our results strengthen the important role of zealots in
models of cyclic dominance, and they reveal fascinating evolutionary outcomes
in structured populations that are a unique consequence of such uncompromising
behavior.


  A procedure is described for estimating evolutionary rate matrices from
observed site frequency data. The procedure assumes (1) that the data are
obtained from a constant size population evolving according to a stationary
Wright-Fisher model; (2) that the data consist of a multiple alignment of a
moderate number of sequenced genomes drawn randomly from the population; and
(3) that within the genome a large number of independent, neutral sites
evolving with with a common mutation rate matrix can be identified. No
restrictions are imposed on the scaled rate matrix other than that the
off-diagonal elements are positive and <<1, and that the rows sum to zero. In
particular the rate matrix is not assumed to be reversible. The key to the
method is an approximate stationary solution to the forward Kolmogorov equation
for the multi-allele neutral Wright-Fisher model in the limit of low mutation
rates.


  A quantitative characterization of the relationship between molecular
sequence and structure is essential to improve our understanding of how
function emerges. This particular genotype-phenotype map has been often studied
in the context of RNA sequences, with the folded configurations standing as a
proxy for the phenotype. Here, we count the secondary structures of circular
RNAs of length $n$ and calculate the asymptotic distributions of different
structural moieties, such as stems or hairpin loops, by means of symbolic
combinatorics. Circular RNAs differ in essential ways from their linear
counterparts. From the mathematical viewpoint, the enumeration of the
corresponding secondary structures demands the use of combinatorial techniques
additional to those used for linear RNAs. The asymptotic number of secondary
structures for circular RNAs grows as $a^nn^{-5/2}$, with a depending on
particular constraints applied to the secondary structure. The abundance of any
structural moiety is normally distributed in the limit $n\to\infty$, with a
mean and a variance that increase linearly with $n$.


  Inspired by recent works on evolutionary graph theory, an area of growing
interest in mathematical and computational biology, we present the first known
examples of undirected structures acting as suppressors of selection for any
fitness value $r > 1$. This means that the average fixation probability of an
advantageous mutant or invader individual placed at some node is strictly less
than that of this individual placed in a well-mixed population. This leads the
way to study more robust structures less prone to invasion, contrary to what
happens with the amplifiers of selection where the fixation probability is
increased on average for advantageous invader individuals. A few families of
amplifiers are known, although some effort was required to prove it. Here, we
use computer aided techniques to find an exact analytical expression of the
fixation probability for some graphs of small order (equal to $6$, $8$ and
$10$) proving that selection is effectively reduced for $r > 1$. Some numerical
experiments using Monte Carlo methods are also performed for larger graphs.


  Time- and state-discrete dynamical systems are frequently used to model
molecular networks. This paper provides a collection of mathematical and
computational tools for the study of robustness in Boolean network models. The
focus is on networks governed by $k$-canalizing functions, a recently
introduced class of Boolean functions that contains the well-studied class of
nested canalizing functions. The activities and sensitivity of a function
quantify the impact of input changes on the function output. This paper
generalizes the latter concept to $c$-sensitivity and provides formulas for the
activities and $c$-sensitivity of general $k$-canalizing functions as well as
canalizing functions with more precisely defined structure. A popular measure
for the robustness of a network, the Derrida value, can be expressed as a
weighted sum of the $c$-sensitivities of the governing canalizing functions,
and can also be calculated for a stochastic extension of Boolean networks.
These findings provide a computationally efficient way to obtain Derrida values
of Boolean networks, deterministic or stochastic, that does not involve
simulation.


  Developing accurate solvers for the Poisson Boltzmann (PB) model is the first
step to make the PB model suitable for implicit solvent simulation. Reducing
the grid size influence on the performance of the solver benefits to increasing
the speed of solver and providing accurate electrostatics analysis for solvated
molecules. In this work, we explore the accurate coarse grid PB solver based on
the Green's function treatment of the singular charges, matched interface and
boundary (MIB) method for treating the geometric singularities, and posterior
electrostatic potential field extension for calculating the reaction field
energy. We made our previous PB software, MIBPB, robust and provides almost
grid size independent reaction field energy calculation. Large amount of the
numerical tests verify the grid size independence merit of the MIBPB software.
The advantage of MIBPB software directly make the acceleration of the PB solver
from the numerical algorithm instead of utilization of advanced computer
architectures. Furthermore, the presented MIBPB software is provided as a free
online sever.


  The threading of a polymer chain through a small pore is a classic problem in
polymer dynamics and underlies nanopore sensing technology. However important
experimental aspects of the polymer motion in a solid-state nanopore, such as
an accurate measurement of the velocity variation during translocation, have
remained elusive. In this work we analysed the translocation through conical
quartz nanopores of a 7 kbp DNA double-strand labelled with six markers equally
spaced along its contour. These markers, constructed from DNA hairpins, give
direct experimental access to the translocation dynamics. On average we measure
a 5% reduction in velocity during the translocation. We also find a striking
correlation in velocity fluctuations with a decay constant of 100s of {\mu}s.
These results shed light on hitherto unresolved problems in the dynamics of DNA
translocation and provide guidance for experiments seeking to determine
positional information along a DNA strand.


  In large clonal populations, several clones generally compete which results
in complex evolutionary and ecological dynamics: experiments show successive
selective sweeps of favorable mutations as well as long-term coexistence of
multiple clonal strains. The mechanisms underlying either coexistence or
fixation of several competing strains have rarely been studied altogether.
Conditions for coexistence have mostly been studied by population and community
ecology, while rates of invasion and fixation have mostly been studied by
population genetics. In order to provide a global understanding of the
complexity of the dynamics observed in large clonal populations, we develop a
stochastic model where three clones compete. Competitive interactions can be
intransitive and we suppose that strains enter the population via mutations or
rare immigrations. We first describe all possible final states of the
population, including stable coexistence of two or three strains, or the
fixation of a single strain. Second, we give estimate of the invasion and
fixation times of a favorable mutant (or immigrant) entering the population in
a single copy. We show that invasion and fixation can be slower or faster when
considering complex competitive interactions. Third, we explore the parameter
space assuming prior distributions of reproduction, death and competitive rates
and we estimate the likelihood of the possible dynamics. We show that when
mutations can affect competitive interactions, even slightly, stable
coexistence is likely. We discuss our results in the context of the
evolutionary dynamics of large clonal populations.


  The assembly of ecological communities from a pool of species is central to
ecology, but the effect of this process on properties of community interaction
networks is still largely unknown. Here, we use a systematic analytical
framework to describe how assembly from a species pool gives rise to community
network properties that differ from those of the pool: Compared to the pool,
the community shows a bias towards higher carrying capacities, weaker
competitive interactions and stronger beneficial interactions. Moreover, even
if interactions between all pool species are completely random, community
networks are more structured, with correlations between interspecies
interactions, and between interactions and carrying capacities. Nonetheless, we
show that these properties are not sufficient to explain the coexistence of all
community species, and that it is a simple relation between interactions and
species abundances that is responsible for the diversity within a community.


  In the present work the spread of epidemic is studied over complex networks
which are characterized by power law degree distribution of links and
heterogeneous rate of disease transmission. The random allocation of epidemic
transmission rates to the nodes results in the heterogeneity, which in turn
causes the segregation of nodes in terms of various sub populations. The aim of
the study is to gain microscopic insight into the effect of interactions among
various sub populations in the spreading processes of disease over such
networks. The discrete time Markov chain method based upon the susceptible
infected susceptible (SIS) model of diseases transmission has been used to
describe the spreading of epidemic over the networks. The study is
parameterized in terms of variable $\lambda$, defined as the number of contacts
a node makes with the fraction of its neighboring nodes. From the simulation
results, it is found that the spread of epidemic on such networks is critical
in terms of number of minimum contacts made by a node below which there is no
outbreak of disease. The degree of infection in these networks is assessed from
the size of epidemic defined in terms of fraction of infected nodes of the
total number and their corresponding level of infection. The results of the
parametric study demonstrates the dependence of the epidemic size upon number
of concurrent contacts made by a node ($\lambda$ ) and the average number of
links per node. In both these cases, the size of the epidemic is found to
increase with the corresponding increase in respective parameters.


  When evaluating the ecological value of land use within a landscape,
investigators typically rely on measures of habitat selection and habitat
quality. Traditional measures of habitat selection and habitat quality require
data from resource intensive study designs (e.g., telemetry or mark-recapture).
Often, managers must evaluate ecological value despite only having data from
less resource intensive study designs. In this paper, we use occupancy data to
measure habitat quality and habitat selection response for the Puerto Rican
Vireo, an endemic songbird whose population growth is depressed by brood
parasitism from the Shiny Cowbird. We were interested in how vireo habitat
quality and vireo habitat selection varied among three land uses (forest,
shaded coffee plantations, and sun coffee) in Puerto Rico. We estimated vireo
occupancy probability as a measure of habitat selection, and the probability of
cowbird occurrence given vireo presence as a measure of habitat quality. To
estimate the latter, we explored different ways of modeling the occurrence of
the two species jointly, and compared these models to independent models using
measures of predictive performance. Vireos preferentially selected forested
sites and shaded coffee sites over sun coffee sites. By our measure of habitat
quality, either type of coffee plantation was poor quality, and the forested
sites were high quality. This suggests that shade coffee may be an ecological
trap for the vireo in our study area. One joint model performed best by our
measures of predictive ability, thus showing that the cowbird may not occur
independently of the vireo. Vireo population dynamics in our study area may
benefit from having large amounts of forest relative to coffee plantations.
Incorporating species interactions into occupancy models has the potential to
improve monitoring for conservation.


  Diffusive motion of regulatory enzymes on biopolymers with eventual capture
at a reaction site is a common feature in cell biology. Using a lattice gas
model we study the impact of diffusion and capture for a microtubule polymerase
and a depolymerase. Our results show that the capture mechanism localizes the
proteins and creates large-scale spatial correlations. We develop an analytic
approximation that globally accounts for relevant correlations and yields
results that are in excellent agreement with experimental data. Our results
show that diffusion and capture operates most efficiently at cellular enzyme
concentrations which points to in vivo relevance.


  General laws in ecological parasitology are scarce. Here we evaluate data
published by over 100 authors to determine whether the number of hosts in a
life cycle is associated with the degree of aggregation of fish parasites at
different stages. Parasite species were grouped taxonomically to produce 20 or
more data points per group as far as possible. Most parasites that remained at
one trophic level were less aggregated than those that had passed up a food
chain. We use a stochastic model to show that high parasite overdispersion in
predators can be solely the result of the accumulation of parasites in their
prey. The model is further developed to show that a change in the predators
feeding behaviour with age may further increase parasite aggregation.


  Nucleosomes form the basic unit of compaction within eukaryotic genomes and
their locations represent an important, yet poorly understood, mechanism of
genetic regulation. Quantifying the strength of interactions within the
nucleosome is a central problem in biophysics and is critical to understanding
how nucleosome positions influence gene expression. By comparing to
single-molecule experiments, we demonstrate that a coarse-grained molecular
model of the nucleosome can reproduce key aspects of nucleosome unwrapping.
Using detailed simulations of DNA and histone proteins, we calculate the
tension-dependent free energy surface corresponding to the unwrapping process.
The model reproduces quantitatively the forces required to unwrap the
nucleosome, and reveals the role played by electrostatic interactions during
this process. We then demonstrate that histone modifications and DNA sequence
can have significant effects on the energies of nucleosome formation. Most
notably, we show that histone tails are crucial for stabilizing the outer turn
of nucleosomal DNA.


  Almost 60 years ago Andrew Huxley with his seminal paper \cite{Huxley1957}
laid the foundation of modern muscle modeling, linking chemical events to
mechanical performance. He described mechanics and energetics of muscle
contraction through the cyclical attachment and detachment of myosin motors to
the actin filament with ad hoc assumptions on the dependence of the rate
constants on the strain of the myosin motors. That relatively simple hypothesis
is still present in recent models, even though with several modifications to
adapt the model to the different experimental constraints which became
subsequently available. However, already in that paper, one controversial
aspect of the model became clear. Relatively high attachment and detachment
rates of myosin to the actin filament were needed to simulate the high power
output at intermediate velocity of contraction. However, these rates were
incompatible with the relatively slow rise in tension after activation, despite
the rise should be generated by the same rate functions. This discrepancy has
not been fully solved till today, despite several hypotheses have been
forwarded to reconcile the two aspects. Here, using a conventional muscle
model, we show that the recently revealed mechanosensing mechanism of
recruitment of myosin motors \cite{Linarietal2015} can solve this long standing
problem without any further ad-hoc hypotheses.


  In this work, a systematic protocol is proposed to automatically parametrize
implicit solvent models with polar and nonpolar components. The proposed
protocol utilizes the classical Poisson model or the Kohn-Sham density
functional theory (KSDFT) based polarizable Poisson model for modeling polar
solvation free energies. For the nonpolar component, either the standard model
of surface area, molecular volume, and van der Waals interactions, or a model
with atomic surface areas and molecular volume is employed. Based on the
assumption that similar molecules have similar parametrizations, we develop
scoring and ranking algorithms to classify solute molecules. Four sets of
radius parameters are combined with four sets of charge force fields to arrive
at a total of 16 different parametrizations for the Poisson model. A large
database with 668 experimental data is utilized to validate the proposed
protocol. The lowest leave-one-out root mean square (RMS) error for the
database is 1.33k cal/mol. Additionally, five subsets of the database, i.e.,
SAMPL0-SAMPL4, are employed to further demonstrate that the proposed protocol
offers some of the best solvation predictions. The optimal RMS errors are 0.93,
2.82, 1.90, 0.78, and 1.03 kcal/mol, respectively for SAMPL0, SAMPL1, SAMPL2,
SAMPL3, and SAMPL4 test sets. These results are some of the best, to our best
knowledge.


  In this work we study the likelihood of survival of single-species in the
context of hostile and disordered environments. Population dynamics in this
environment, as modeled by the Fisher equation, is characterized by negative
average growth rate, except in some random spatially distributed patches that
may support life. In particular, we are interested in the phase diagram of the
survival probability and in the critical size problem, i.e., the minimum patch
size required for surviving in the long time dynamics. We propose a measure for
the critical patch size as being proportional to the participation ratio (PR)
of the eigenvector corresponding to the largest eigenvalue of the linearized
Fisher dynamics. We obtain the (extinction-survival) phase diagram and the
probability distribution function (PDF) of the critical patch sizes for two
topologies, namely, the one-dimensional system and the fractal Peano basin. We
show that both topologies share the same qualitative features, but the fractal
topology requires higher spatial fluctuations to guarantee species survival. We
perform a finite-size scaling and we obtain the associated scaling exponents.
In addition, we show that the PDF of the critical patch sizes has an universal
shape for the 1D case in terms of the model parameters (diffusion, growth rate,
etc.). In contrast, the diffusion coefficient has a drastic effect on the PDF
of the critical patch sizes of the fractal Peano basin, and it does not obey
the same scaling law of the 1D case.


  Phylogenetic networks are increasingly used in evolutionary biology to
represent the history of species that have undergone reticulate events such as
horizontal gene transfer, hybrid speciation and recombination. One of the most
fundamental questions that arise in this context is whether the evolution of a
gene with one copy in all species can be explained by a given network. In
mathematical terms, this is often translated in the following way: is a given
phylogenetic tree contained in a given phylogenetic network? Recently this tree
containment problem has been widely investigated from a computational
perspective, but most studies have only focused on the topology of the phylo-
genies, ignoring a piece of information that, in the case of phylogenetic
trees, is routinely inferred by evolutionary analyses: branch lengths. These
measure the amount of change (e.g., nucleotide substitutions) that has occurred
along each branch of the phylogeny. Here, we study a number of versions of the
tree containment problem that explicitly account for branch lengths. We show
that, although length information has the potential to locate more precisely a
tree within a network, the problem is computationally hard in its most general
form. On a positive note, for a number of special cases of biological
relevance, we provide algorithms that solve this problem efficiently. This
includes the case of networks of limited complexity, for which it is possible
to recover, among the trees contained by the network with the same topology as
the input tree, the closest one in terms of branch lengths.


  Background: Many mathematical models have now been employed across every area
of systems biology. These models increasingly involve large numbers of unknown
parameters, have complex structure which can result in substantial evaluation
time relative to the needs of the analysis, and need to be compared to observed
data. The correct analysis of such models usually requires a global parameter
search, over a high dimensional parameter space, that incorporates and respects
the most important sources of uncertainty. This can be an extremely difficult
task, but it is essential for any meaningful inference or prediction to be made
about any biological system. It hence represents a fundamental challenge for
the whole of systems biology.
  Results: Bayesian statistical methodology for the uncertainty analysis of
complex models is introduced, which is designed to address the high dimensional
global parameter search problem. Bayesian emulators that mimic the systems
biology model but which are extremely fast to evaluate are embedded within an
iterative history match: an efficient method to search high dimensional spaces
within a more formal statistical setting, while incorporating major sources of
uncertainty. The approach is demonstrated via application to two models of
hormonal crosstalk in Arabidopsis root development, which have 32 rate
parameters, for which we identify the sets of rate parameter values that lead
to acceptable matches to observed trend data. The biological consequences of
the resulting comparison, including the evaluation of gene functions, are
described.


  The bioprotective nature of disaccharides is hypothesized to derive from the
modification of the hydrogen bonding network of water which protects
biomolecules through lowered water activity at the protein interface. Using
ultrafast fluorescence spectroscopy we measured the relaxation of bulk water
dynamics around the induced dipole moment of two fluorescent probes (Lucifer
Yellow Ethylenediamine and Tryptophan). Our results indicate a reduction in
bulk water reorganization rate of approximately of 30%. We observe this
retardation in the low concentration regime measured at 0.1M and 0.25 M, far
below the onset of glassy dynamics. This reduction in water activity could be
significant in crowded biological systems, contributing to global change in
protein energy landscape, resulting in a significant enhancement of protein
stability under environmental stress. We observed similar dynamic reduction for
two disaccharide osmolytes, sucrose and trehalose, with trehalose being the
more effective dynamic reducer.


  Fluorescence microscopy reveals that the contents of many (membrane-free)
nuclear "bodies" exchange rapidly with the soluble pool whilst the underlying
structure persists; such observations await a satisfactory biophysical
explanation. To shed light on this, we perform large-scale Brownian dynamics
simulations of a chromatin fiber interacting with an ensemble of (multivalent)
DNA-binding proteins; these proteins switch between two states -- active
(binding) and inactive (non-binding). This system provides a model for any
DNA-binding protein that can be modified post-translationally to change its
affinity for DNA (e.g., like the phosphorylation of a transcription factor).
Due to this out-of-equilibrium process, proteins spontaneously assemble into
clusters of self-limiting size, as individual proteins in a cluster exchange
with the soluble pool with kinetics like those seen in photo-bleaching
experiments. This behavior contrasts sharply with that exhibited by
"equilibrium", or non-switching, proteins that exist only in the binding state;
when these bind to DNA non-specifically, they form clusters that grow
indefinitely in size. Our results point to post-translational modification of
chromatin-bridging proteins as a generic mechanism driving the self-assembly of
highly dynamic, non-equilibrium, protein clusters with the properties of
nuclear bodies. Such active modification also reshapes intra-chromatin contacts
to give networks resembling those seen in topologically-associating domains, as
switching markedly favors local (short-range) contacts over distant ones.


  Many of the processes that underly neural computation are carried out by ion
channels embedded in the plasma membrane, a two-dimensional liquid that
surrounds all cells. Recent experiments have demonstrated that this membrane is
poised close to a liquid-liquid critical point in the Ising universality class.
Here we use both exact and stochastic techniques on the lattice Ising model to
explore the ramifications of proximity to criticality for proteins that are
allosterically coupled to Ising composition modes. Owing to diverging
generalized susceptibilities, such a protein's activity becomes strongly
influenced by perturbations that influence the two relevant parameters of the
critical point, especially the critical temperature. In addition, the protein's
kinetics acquire a range of time scales from its surrounding membrane,
naturally leading to non-Markovian dynamics.


  Guided movement in response to slowly diffusing polymeric trails provides a
unique mechanism for self-organization of some microorganisms. To elucidate how
this signaling route leads to microcolony formation, we experimentally probe
the trajectory and orientation of Pseudomonas aeruginosa that propel themselves
on a surface using type IV pili motility appendages, which preferentially
attach to deposited exopolysaccharides. We construct a stochastic model by
analyzing single-bacterium trajectories, and show that the resulting
theoretical prediction for the many-body behavior of the bacteria is in
quantitative agreement with our experimental characterization of how cells
explore the surface via a power law strategy.


  A classic problem in computational biology is constructing a phylogenetic
tree given a set of distances between n species. In most cases, a tree
structure is too constraining. We consider a circular split network, a
generalization of a tree in which multiple parallel edges signify divergence. A
geometric space of such networks is introduced, forming a natural extension of
the work by Billera, Holmes, and Vogtmann on tree space. We explore properties
of this space, and show a natural embedding of the compactification of the real
moduli space of curves within it.


  Controlling stochastic reactions networks is a challenging problem with
important implications in various fields such as systems and synthetic biology.
Various regulation motifs have been discovered or posited over the recent
years, the most recent one being the so-called Antithetic Integral Control
(AIC) motif in Briat et al. (Cell Systems, 2016). Several favorable properties
for the AIC motif have been demonstrated for classes of reaction networks that
satisfy certain irreducibility, ergodicity and output controllability
conditions. Here we address the problem of verifying these conditions for large
sets of reaction networks with fixed topology using two different approaches.
The first one is quantitative and relies on the notion of interval matrices
while the second one is qualitative and is based on sign properties of
matrices. The obtained results lie in the same spirit as those obtained in
Briat et al. (Cell Systems, 2016) where properties of reaction networks are
independently characterized in terms of control theoretic concepts, linear
programming conditions and graph theoretic conditions.


  Protein molecules in cells are synthesized by macromolecular machines called
ribosomes. According to recent experimental data, we reduce the complexity of
the ribosome and propose a model to express its activity in six main states.
Using our model, we study the translation rate in different biological relevant
situations in the presence of external force, and translation through the RNA
double stranded region in the absence or presence of the external force. In the
present study, we give a quantitative theory for translation rate and show that
the ribosome behaves more like a Brownian Ratchet motor. Our findings could
shed some light on understanding behaviors of the ribosome in biological
conditions.


  Cells of the immune system are confronted with opposing pro- and
anti-inflammatory signals. Dendritic cells (DC) integrate these cues to make
informed decisions whether to initiate an immune response. Confronted with
exogenous microbial stimuli, DC endogenously produce both anti- (IL-10) and
pro-inflammatory (TNF$\alpha$) cues whose joint integration controls the cell's
final decision. We combine experimental measurements with theoretical modeling
to quantitatively describe the integration mode of these opposing signals. We
propose a two step integration model that modulates the effect of the two types
of signals: an initial bottleneck integrates both signals (IL-10 and
TNF$\alpha$), the output of which is later modulated by the anti-inflammatory
signal. We show that the anti-inflammatory IL-10 signaling is long ranged, as
opposed to the short-ranged pro-inflammatory TNF$\alpha$ signaling. The model
suggests that the population averaging and modulation of the pro-inflammatory
response by the anti-inflammatory signal is a safety guard against excessive
immune responses.


  In considering evolution of transcribed regions, regulatory modules, and
other genomic loci of interest, we are often faced with a situation in which
the number of allelic states greatly exceeds the population size. In this
limit, the population eventually adopts a steady state characterized by
mutation-selection-drift balance. Although new alleles continue to be explored
through mutation, the statistics of the population, and in particular the
probabilities of seeing specific allelic configurations in samples taken from a
population, do not change with time. In the absence of selection, probabilities
of allelic configurations are given by the Ewens sampling formula, widely used
in population genetics to detect deviations from neutrality. Here we develop an
extension of this formula to arbitrary, possibly epistatic, fitness landscapes.
Although our approach is general, we focus on the class of landscapes in which
alleles are grouped into two, three, or several fitness states. This class of
landscapes yields sampling probabilities that are computationally more
tractable, and can form a basis for the inference of selection signatures from
sequence data. We demonstrate that, for a sizeable range of mutation rates and
selection coefficients, the steady-state allelic diversity is not neutral.
Therefore, it may be used to infer selection coefficients, as well as other key
evolutionary parameters, using high-throughput sequencing of evolving
populations to collect data on locus polymorphisms. We also find that our
theory remains sufficiently accurate even if the assumptions such as the
infinite allele limit and the "full connectivity" assumption in which each
allele can mutate into any other allele are relaxed. Thus, our framework
establishes a theoretical foundation for inferring selection signatures from
samples of sequences produced by evolution on epistatic fitness landscapes.


  Recently, an elegant approach in phylogenetics was introduced by
Billera-Holmes-Vogtmann that allows a systematic comparison of different
evolutionary histories using the metric geometry of tree spaces. In many
problem settings one encounters heavily populated phylogenetic trees, where the
large number of leaves encumbers visualization and analysis in the relevant
evolutionary moduli spaces. To address this issue, we introduce tree
dimensionality reduction, a structured approach to reducing large phylogenetic
trees to a distribution of smaller trees. We prove a stability theorem ensuring
that small perturbations of the large trees are taken to small perturbations of
the resulting distributions.
  We then present a series of four biologically motivated applications to the
analysis of genomic data, spanning cancer and infectious disease. The first
quantifies how chemotherapy can disrupt the evolution of common leukemias. The
second examines a link between geometric information and the histologic grade
in relapsed gliomas, where longer relapse branches were specific to high grade
glioma. The third concerns genetic stability of xenograft models of cancer,
where heterogeneity at the single cell level increased with later mouse
passages. The last studies genetic diversity in seasonal influenza A virus. We
apply tree dimensionality reduction to 24 years of longitudinally collected
H3N2 hemagglutinin sequences, generating distributions of smaller trees
spanning between three and five seasons. A negative correlation is observed
between the influenza vaccine effectiveness during a season and the variance of
the distributions produced using preceding seasons' sequence data. We also show
how tree distributions relate to antigenic clusters and choice of influenza
vaccine. Our formalism exposes links between viral genomic data and clinical
observables such as vaccine selection and efficacy.


  Scientists use mathematical modelling to understand and predict the
properties of complex physical systems. In highly parameterised models there
often exist relationships between parameters over which model predictions are
identical, or nearly so. These are known as structural or practical
unidentifiabilities, respectively. They are hard to diagnose and make reliable
parameter estimation from data impossible. They furthermore imply the existence
of an underlying model simplification. We describe a scalable method for
detecting unidentifiabilities, and the functional relations defining them, for
generic models. This allows for model simplification, and appreciation of which
parameters (or functions thereof) cannot be estimated from data. Our algorithm
can identify features such as redundant mechanisms and fast timescale
subsystems, as well as the regimes in which such approximations are valid. We
base our algorithm on a novel quantification of regional parametric
sensitivity: multiscale sloppiness. Traditionally, the link between parametric
sensitivity and the conditioning of the parameter estimation problem is made
locally, through the Fisher Information Matrix. This is valid in the regime of
infinitesimal measurement uncertainty. We demonstrate the duality between
multiscale sloppiness and the geometry of confidence regions surrounding
parameter estimates made where measurement uncertainty is non-negligible.
Further theoretical relationships are provided linking multiscale sloppiness to
the Likelihood-ratio test. From this, we show that a local sensitivity analysis
(as typically done) is insufficient for determining the reliability of
parameter estimation, even with simple (non)linear systems. Our algorithm
provides a tractable alternative. We finally apply our methods to a
large-scale, benchmark Systems Biology model of NF-$\kappa$B, uncovering
previously unknown unidentifiabilities.


  Theoretical results regarding two-dimensional ordinary-differential equations
(ODEs) with second-degree polynomial right-hand sides are summarized, with an
emphasis on limit cycles, limit cycle bifurcations and multistability. The
results are then used for construction of two reaction systems, which are at
the deterministic level described by two-dimensional third-degree kinetic ODEs.
The first system displays a homoclinic bifurcation, and a coexistence of a
stable critical point and a stable limit cycle in the phase plane. The second
system displays a multiple limit cycle bifurcation, and a coexistence of two
stable limit cycles. The deterministic solutions (obtained by solving the
kinetic ODEs) and stochastic solutions (noisy time-series generating by the
Gillespie algorithm, and the underlying probability distributions obtained by
solving the chemical master equation (CME)) of the constructed systems are
compared, and the observed differences highlighted. The constructed systems are
proposed as test problems for statistical methods, which are designed to detect
and classify properties of given noisy time-series arising from biological
applications.


  We assume that the protein folding process follows two autonomous steps: the
conformational search for the native, mainly ruled by the hydrophobic effect;
and, the final adjustment stage, which eventually gives stability to the
native. Our main tool of investigation is a 3D lattice model provided with a
ten-letter alphabet, the stereochemical model. This model was conceived for
Monte Carlo (MC) simulations when one keeps in mind the kinetic behavior of
protein-like chains in solution. In order to characterize the folding
characteristic time ({\tau}) by two distinct sampling methods, first we present
two sets of 10^{3} MC simulations for a fast protein-like sequence. For these
sets of folding times, {\tau} and {\tau}_{q} were obtained with the application
of the standard Metropolis algorithm (MA), and a modified algorithm (M_{q}A).
The results for {\tau}_{q}reveal two things: i) the hydrophobic chain-solvent
interactions plus a set of inter-residues steric constraints are enough to
emulate the first stage of the process: for each one of the 10^{3} MC performed
simulations, the native is always found without exception, ii) the ratio
{\tau}_{q}/{\tau}~1/3 suggests that the effect of local thermal fluctuations,
encompassed by the Tsallis weight, provides an innate efficiency to the chain
escapes from energetic and steric traps. ...


  The Ornstein-Uhlenbeck (OU) process plays a major role in the analysis of the
evolution of phenotypic traits along phylogenies. The standard OU process
includes drift and stabilizing selection and assumes that species evolve
independently. However, especially in plants, there is ample evidence of
hybridization and introgression during evolution. In this work we present a
statistical approach with analytical solutions that allows for the inclusion of
adaptation and migration in a common phylogenetic framework. We furthermore
present a detailed simulation study that clearly indicates the adverse effects
of ignoring migration. Similarity between species due to migration could be
misinterpreted as very strong convergent evolution without proper correction
for these additional dependencies. Our model can also be useful for studying
local adaptation among populations within the same species. Finally, we show
that our model can be interpreted in terms of ecological interactions between
species, providing a general framework for the evolution of traits between
"interacting" species or populations.


  We present a mathematical model, based on ordinary differential equations,
for the evolution of solid tumors and their response to treatment. Specifically
the effects of a cytotoxic agent and a monoclonal antibody are included as
control term in the equations. The variables considered here are: the number of
cancerous cells sensitive to chemotherapy, the number of cancerous cells
resistant to chemotherapy, the degree of angiogenesis and the average intensity
of VEGF. The rules that govern the quantities mentioned above are based on a
geometrical argument: we approximate the tumor mass as a sphere and thus derive
basic formulae for the normoxic cells and for VEGF production. The monoclonal
antibody acts on VEGF and thus has in uence to the global degree of
angiogenesis. Numerical estimates on some of the parameters are performed in
order to match the main landmark in tumor progression and reaction to treatment
in the specific case of colorectal cancer.


  A multi-resolution bead-spring model for polymer dynamics is developed as a
generalization of the Rouse model. A polymer chain is described using beads of
variable sizes connected by springs with variable spring constants. A numerical
scheme which can use different timesteps to advance the positions of different
beads is presented and analyzed. The position of a particular bead is only
updated at integer multiples of the timesteps associated with its connecting
springs. This approach extends the Rouse model to a multiscale model on both
spatial and temporal scales, allowing simulations of localized regions of a
polymer chain with high spatial and temporal resolution, while using a coarser
modelling approach to describe the rest of the polymer chain. A method for
changing the model resolution on-the-fly is developed using the
Metropolis-Hastings algorithm. It is shown that this approach maintains key
statistics of the end-to-end distance and diffusion of the polymer filament and
makes computational savings when applied to a model for the binding of a
protein to the DNA filament.


  The theory of computer science is based around Universal Turing Machines
(UTMs): abstract machines able to execute all possible algorithms. Modern
digital computers are physical embodiments of UTMs. The nondeterministic
polynomial (NP) time complexity class of problems is the most significant in
computer science, and an efficient (i.e. polynomial P) way to solve such
problems would be of profound economic and social importance. By definition
nondeterministic UTMs (NUTMs) solve NP complete problems in P time. However,
NUTMs have previously been believed to be physically impossible to construct.
Thue string rewriting systems are computationally equivalent to UTMs, and are
naturally nondeterministic. Here we describe the physical design for a NUTM
that implements a universal Thue system. The design exploits the ability of DNA
to replicate to execute an exponential number of computational paths in P time.
Each Thue rewriting step is embodied in a DNA edit implemented using a novel
combination of polymerase chain reactions and site-directed mutagenesis. We
demonstrate that this design works using both computational modelling and in
vitro molecular biology experimentation. The current design has limitations,
such as restricted error-correction. However, it opens up the prospect of
engineering NUTM based computers able to outperform all standard computers on
important practical problems.


  Inferring evolutionary histories (phylogenetic trees) has important
applications in biology, criminology and public health. However, phylogenetic
trees are complex mathematical objects that reside in a non-Euclidean space,
which complicates their analysis. While our mathematical, algorithmic, and
probabilistic understanding of phylogenies in their metric space is mature,
rigorous inferential infrastructure is as yet undeveloped. In this manuscript
we unify recent computational and probabilistic advances to construct
tree--valued confidence sets. The procedure accounts for both centre and
multiple directions of tree--valued variability. We draw on block replicates to
improve testing, identifying the best supported most recent ancestor of the
Zika virus, and formally testing the hypothesis that a Floridian dentist with
AIDS infected two of his patients with HIV. The method illustrates connections
between variability in Euclidean and tree space, opening phylogenetic tree
analysis to techniques available in the multivariate Euclidean setting.


  The association between human cytomegalovirus (hCMV) reactivation and the
development of graft-versus-host-disease (GVHD) has been observed in stem cell
transplantation (SCT). Seventy seven SCT donor-recipient pairs (DRP) (HLA
matched unrelated donor (MUD), n=50; matched related donor (MRD), n=27)
underwent whole exome sequencing to identify single nucleotide polymorphisms
(SNPs) generating alloreactive peptide libraries for each DRP (9-mer
peptide-HLA complexes); Human CMV CROSS (Cross-Reactive Open Source Sequence)
Database was compiled from NCBI; HLA class I binding affinity for each DRPs HLA
was calculated by NetMHCpan 2.8 and hCMV- derived 9-mers algorithmically
compared to the alloreactive peptide-HLA complex libraries. Short consecutive
(6 or greater) amino acid (AA) sequence homology matching hCMV to recipient
peptides was considered for HLA-bound-peptide (IC50<500 nM) cross reactivity.
Of the 70,686 hCMV 9-mers contained within the hCMV CROSS database, 29,658.8
+/- 9038.5 were found to match MRD DRP alloreactive peptides and 52,910.2 +/-
16121.8 matched MUD DRP peptides (Student's T-test, p<0.001). In silico
analysis revealed multiple high affinity, immunogenic CMV-Human peptide matches
(IC50<500 nM) expressed in GVHD-affected tissue-specific manner (proteins
expressed at 10 RPKM or greater). hCMV+GVHD was found in 18 patients, 13
developing hCMV viremia before GVHD onset with a subset analysis of 7 instances
of hCMV viremia prior to acute GVHD onset (n=3), chronic GVHD (n=2) and acute +
chronic GVHD (n=2) indicating cross reactive peptide expression within affected
organs. We propose that based on our analysis and preliminary clinical
correlations that hCMV immune cross-reactivity may cause antigenic mimicry of
human alloreactive peptides triggering GVHD.


  Muscle uses Ca2+ as a messenger to control contraction and relies on ATP to
maintain the intracellular Ca2+ homeostasis. Mitochondria are the major
sub-cellular organelle of ATP production. With a negative inner membrane
potential, mitochondria take up Ca2+ from their surroundings, a process called
mitochondrial Ca2+ uptake. Under physiological conditions, Ca2+ uptake into
mitochondria promotes ATP production. Excessive uptake causes mitochondrial
Ca2+ overload, which activates downstream adverse responses leading to cell
dysfunction. Moreover, mitochondrial Ca2+ uptake could shape spatio-temporal
patterns of intracellular Ca2+ signaling. Malfunction of mitochondrial Ca2+
uptake is implicated in muscle degeneration. Unlike non-excitable cells,
mitochondria in muscle cells experience dramatic changes of intracellular Ca2+
levels. Besides the sudden elevation of Ca2+ level induced by action
potentials, Ca2+ transients in muscle cells can be as short as a few
milliseconds during a single twitch or as long as minutes during tetanic
contraction, which raises the question whether mitochondrial Ca2+ uptake is
fast and big enough to shape intracellular Ca2+ signaling during
excitation-contraction coupling and creates technical challenges for
quantification of the dynamic changes of Ca2+ inside mitochondria. This review
focuses on characterization of mitochondrial Ca2+ uptake in skeletal muscle and
its role in muscle physiology and diseases.


  Tolerance implies enduring trying circumstances with a fair and objective
attitude. To determine whether evolutionary advantages might be stemming from
diverse levels of tolerance in a population, we study a spatial public goods
game, where in addition to cooperators, defectors, and loners, tolerant players
are also present. Depending on the number of defectors within a group, a
tolerant player can either cooperate in or abstain from a particular instance
of the game. We show that the diversity of tolerance can give rise to
synergistic effects, wherein players with a different threshold in terms of the
tolerated number of defectors in a group compete most effectively against
defection and default abstinence. Such synergistic associations can stabilise
states of full cooperation where otherwise defection would dominate. We observe
complex pattern formation that gives rise to an intricate phase diagram, where
invisible yet stable strategy alliances require outmost care lest they are
overlooked. Our results highlight the delicate importance of diversity and
tolerance for the provisioning of public goods, and they reveal fascinating
subtleties of the spatiotemporal dynamics that is due to the competition of
subsystem solutions in structured populations.