  Understanding the mathematical properties of graphs underling biological
systems could give hints on the evolutionary mechanisms behind these
structures. In this article we perform a complete statistical analysis over
thousands of graphs representing metabolic and protein-protein interaction
(PPI) networks. First, we investigate the quality of fits obtained for the
nodes degree distributions to power-law functions. This analysis suggests that
a power-law distribution poorly describes the data except for the far right
tail in the case of PPI networks. Next we obtain descriptive statistics for the
main graph parameters and try to identify the properties that deviate from the
expected values had the networks been built by randomly linking nodes with the
same degree distribution. This survey identifies the properties of biological
networks which are not solely the result of their degree distribution, but
emerge from yet unidentified mechanisms other than those that drive these
distributions. The findings suggest that, while PPI networks have properties
that differ from their expected values in their randomized versions with great
statistical significance, the differences for metabolic networks have a smaller
statistical significance, though it is possible to identify some drift.


  In this report, we present an unsupervised machine learning method for
determining groups of molecular systems according to similarity in their
dynamics or structures using Ward's minimum variance objective function. We
first apply the minimum variance clustering to a set of simulated tripeptides
using the information theoretic Jensen-Shannon divergence between Markovian
transition matrices in order to gain insight into how point mutations affect
protein dynamics. Then, we extend the method to partition two chemoinformatic
datasets according to structural similarity to motivate a train/validation/test
split for supervised learning that avoids overfitting.


  Even the simplest organisms are too complex to have spontaneously arisen
fully-formed, yet precursors to first life must have emerged ab initio from
their environment. A watershed event was the appearance of the first entity
capable of evolution: the Initial Darwinian Ancestor. Here we suggest that
nucleopeptide reciprocal replicators could have carried out this important role
and contend that this is the simplest way to explain extant replication systems
in a mathematically consistent way. We propose short nucleic-acid templates on
which amino-acylated adapters assembled. Spatial localization drives peptide
ligation from activated precursors to generate phosphodiester-bond-catalytic
peptides. Comprising autocatalytic protein and nucleic acid sequences, this
dynamical system links and unifies several previous hypotheses and provides a
plausible model for the emergence of DNA and the operational code.


  A metapopulation structure in landscape ecology comprises a group of
interacting spatially separated subpopulations or patches of the same species
that may experience several local extinctions. This makes the investigation of
survivability (in the form of global oscillation) of a metapopulation on top of
diverse dispersal topologies extremely crucial. However, among various
dispersal topologies in ecological networks, which one can provide higher
metapopulation survivability under local extinction is still not well explored.
In this article, we scrutinize the robustness of an ecological network
consisting of prey-predator patches having Holling type I functional response,
against progressively extinct population patches. We present a comprehensive
study on this while considering global, small-world and scale-free dispersal of
the subpopulations. Furthermore, we extend our work in enhancing survivability
in the form of sustained global oscillation by introducing asymmetries in the
dispersal rates of the considered species. Our findings affirm that the
asynchrony among the patches plays an important role in the survivability of a
metapopulation. In order to demonstrate the model independence of the observed
phenomenon, we perform a similar analysis for patches exhibiting Holling type
II functional response. On the grounds of the obtained results, our work is
expected to provide a better perception of the influence of dispersal
arrangements on the global survivability of ecological networks.


  In this note, we investigate the relationship between probabilistic updating
mechanisms and discrete-time replicator-mutator dynamics. We consider the
recently shown connection between Bayesian updating and replicator dynamics and
extend it to the replicator-mutator dynamics by considering prediction and
filtering recursions in hidden Markov models (HMM). We show that it is possible
to understand the evolution of the frequency vector of a population under the
replicator-mutator equation as a posterior predictive inference procedure in an
HMM. This view enables us to derive a natural dual version of the
replicator-mutator equation, which corresponds to updating the filtering
distribution. Finally, we conclude with the implications of the interpretation
and with some comments related to the recent discussions about evolution and
learning.


  We consider three regions with different public health conditions. In the
absence of migration among these regions, the first two have good health
conditions and the disease free state is stable; for the third region, on the
other hand, the only stable state is the endemic one. When migration is
included in the model, we assume that the second region has a disease risk that
makes its inhabitants prone to accept to be vaccinated, while the population in
the first region tends to reject the vaccination, considered riskier that the
disease. Therefore, the second region is a "buffer zone" between the two
extremal regions. We study the basic reproductive ratio as a function of the
vaccination in all regions and migration among them. This problem is studied
numerically, showing explicit situations in which migration will have an
overall positive effect in the disease dynamics, with and without vaccinations.
We also find explicit formula in the limit of small ("closed borders") and high
migration ("open borders").


  Environmental variability greatly influences the eco-evolutionary dynamics of
a population, i.e. it affects how its size and composition evolve. Here, we
study a well-mixed population of finite and fluctuating size whose growth is
limited by a randomly switching carrying capacity. This models the
environmental fluctuations between states of resources abundance and scarcity.
The population consists of two strains, one growing slightly faster than the
other, competing under two scenarios: one in which competition is solely for
resources, and one in which the slow ("cooperating") strain produces a public
good that benefits also the fast ("freeriding") strain. We investigate how the
coupling of demographic and environmental (external) noise affects the
population's eco-evolutionary dynamics. By analytical and computational means,
we study the correlations between the population size and its composition, and
discuss the social-dilemma-like "eco-evolutionary game" characterizing the
public good production. We determine in what conditions it is best to produce a
public good; when cooperating is beneficial but outcompeted by freeriding, and
when the public good production is detrimental for cooperators. Within a linear
noise approximation to populations of varying size, we also accurately analyze
the coupled effects of demographic and environmental noise on the size
distribution.


  Molecular dynamics (MD) simulations can model the interactions between
macromolecules with high spatiotemporal resolution but at a high computational
cost. By combining high-throughput MD with Markov state models (MSMs), it is
now possible to obtain long-timescale behavior of small to intermediate
biomolecules and complexes. To model the interactions of many molecules at
large lengthscales, particle-based reaction-diffusion (RD) simulations are more
suitable but lack molecular detail. Thus, coupling MSMs and RD simulations
(MSM/RD) would be highly desirable, as they could efficiently produce
simulations at large time- and lengthscales, while still conserving the
characteristic features of the interactions observed at atomic detail. While
such a coupling seems straightforward, fundamental questions are still open:
Which definition of MSM states is suitable? Which protocol to merge and split
RD particles in an association/dissociation reaction will conserve the correct
bimolecular kinetics and thermodynamics? In this paper, we make the first step
towards MSM/RD by laying out a general theory of coupling and proposing a first
implementation for association/dissociation of a protein with a small ligand (A
+ B <--> C). Applications on a toy model and CO diffusion into the heme cavity
of myoglobin are reported.


  In early development, regulation of transcription results in precisely
positioned and highly reproducible expression patterns that specify cellular
identities. How transcription, a fundamentally noisy molecular process, is
regulated to achieve reliable embryonic patterning remains unclear. In
particular, it is unknown how gene-specific regulation mechanisms affect
kinetic rates of transcription, and whether there are common, global features
that govern these rates across a genetic network. Here, we measure nascent
transcriptional activity in the gap gene network of early Drosophila embryos
and characterize the variability in absolute activity levels across expression
boundaries. We demonstrate that boundary formation follows a common
transcriptional principle: a single control parameter determines the
distribution of transcriptional activity, regardless of gene identity, boundary
position, or enhancer-promoter architecture. By employing a minimalist model of
transcription, we infer kinetic rates of transcriptional bursting for these
patterning genes; we find that the key regulatory parameter is the fraction of
time a gene is in an actively transcribing state, while the rate of Pol II
loading appears globally conserved. These results point to a universal
simplicity underlying the apparently complex transcriptional processes
responsible for early embryonic patterning and indicate a path to general rules
in transcriptional regulation.


  A computational model of aquaporin regulation in cancer cells has been
constructed as a Qualitative Network in the software BioModelAnalyzer (BMA).
The model connects some important aquaporins expressed in human cancer to
common phenotypes via a number of fundamental, dysregulated signalling
pathways. Based on over 60 publications, this model can not only reproduce the
results reported in a discrete, qualitative manner, but also reconcile the
seemingly incompatible phenotype with research consensus by suggesting
molecular mechanisms accountable for it. Novel predictions have also been made
by mimicking real-life experiments in the model.


  By challenging E. coli with sublethal norfloxacin for 10 days, Henry Lee and
James Collins suggests the bacterial altruism leads to the population-wide
resistance. By detailedly analyzing experiment data, we suggest that bacterial
cooperation leads to population-wide resistance under norfloxacin pressure and
simultaneously propose the bacteria shield is the possible feedback mechanism
of less resistant bacteria. The bacteria shield is that the less resistant
bacteria sacrifice the large number of themselves to consume norfloxacin and
then to relieve the norfloxacin burden from highly resistant bacteria. Thus,
due to highly resistant bacteria and less resistant bacteria extracted from the
same bacteria population, bacterial cooperation leads to heteroresistance.


  This paper considers a mathematical model based on the transmission dynamics
of hepatitis C virus (HCV) infection. In addition to the usual compartments for
susceptible, exposed, and infected individuals, this model includes
compartments for individuals who are under treatment and those who have had
vaccination against HCV infection. It is assumed that the immunity provided by
the vaccine fades with time. The basic reproduction number, $R_0$, and the
equilibrium solutions of the model are determined. The model exhibits the
phenomenon of backward bifurcation where a stable disease-free equilibrium
co-exists with a stable endemic equilibrium whenever $R_0$ is less than unity.
It is shown that the use of only a perfect vaccine can eliminate backward
bifurcation completely. Furthermore, a unique endemic equilibrium of the model
is proved to be globally asymptotically stable under certain restrictions on
the parameter values. Numerical simulation results are given to support the
theoretical predictions. [epidemiological model; equilibrium solutions;
backward bifurcation; global asymptotic stability; Lyapunov function.]


  Mechanoelectrical transduction in the inner ear is a biophysical process
underlying the senses of hearing and balance. The key players involved in this
process are mechanosensitive ion channels. They are located in the stereocilia
of hair cells and opened by the tension in specialized molecular springs, the
tip links, connecting adjacent stereocilia. When channels open, the tip links
relax, reducing the hair-bundle stiffness. This gating compliance makes hair
cells especially sensitive to small stimuli. The classical explanation for the
gating compliance is that the conformational rearrangement of a single channel
directly shortens the tip link. However, to reconcile theoretical models based
on this mechanism with experimental data, an unrealistically large structural
change of the channel is required. Experimental evidence indicates that each
tip link is a dimeric molecule, associated on average with two channels at its
lower end. It also indicates that the lipid bilayer modulates channel gating,
although it is not clear how. Here, we design and analyze a model of
mechanotransduction where each tip link attaches to two channels, mobile within
the membrane. Their states and positions are coupled by membrane-mediated
elastic forces arising from the interaction between the channels' hydrophobic
cores and that of the lipid bilayer. This coupling induces cooperative opening
and closing of the channels. The model reproduces the main properties of
hair-cell mechanotransduction using only realistic parameters constrained by
experimental evidence. This work provides an insight into the fundamental role
that membrane-mediated ion-channel cooperativity can play in sensory
physiology.


  Stochastic simulation has been a powerful tool for studying the dynamics of
gene regulatory networks, particularly in terms of understanding how
cell-phenotype stability and fate-transitions are impacted by noisy gene
expression. However, gene networks often have dynamics characterized by
multiple attractors. Stochastic simulation is often inefficient for such
systems, because most of the simulation time is spent waiting for rare,
barrier-crossing events to occur. We present a rare-event simulation-based
method for computing epigenetic landscapes and phenotype-transitions in
metastable gene networks. Our computational pipeline was inspired by studies of
metastability and barrier-crossing in protein folding, and provides an
automated means of computing and visualizing essential stationary and dynamic
information that is generally inaccessible to conventional simulation. Applied
to a network model of pluripotency in Embryonic Stem Cells, our simulations
revealed rare phenotypes and approximately Markovian transitions among
phenotype-states, occurring with a broad range of timescales. The relative
probabilities of phenotypes and the transition paths linking pluripotency and
differentiation are sensitive to global kinetic parameters governing
transcription factor-DNA binding kinetics. Our approach significantly expands
the capability of stochastic simulation to investigate gene regulatory network
dynamics, which may help guide rational cell reprogramming strategies. Our
approach is also generalizable to other types of molecular networks and
stochastic dynamics frameworks.


  In the real world, most objects and data have multiple types of attributes
and inter-connections. Such data structures are named "Heterogeneous
Information Networks" (HIN) and have been widely researched. Biological systems
are also considered to be highly complicated HIN. In this work, we review
various applications of HIN methods to biological and chemical data, discuss
some advanced topics, and describe some future research directions.


  Coffee rust is one of the main diseases that affect coffee plantations
worldwide. This causes an important economic impact in the coffee production
industry in countries where coffee is an important part of the economy. A
common method for combating this disease is using copper hydroxide as a
fungicide, which can have damaging effects both on the coffee tree and on human
health. A novel method for biological control of coffee rust using bacteria has
been proven to be an effective alternative to copper hydroxide fungicides as
anti-fungal compounds. In this paper, we develop and explore a spatial
stochastic model for this interaction in a coffee plantation. We analyze
equilibria for specific control strategies, as well as compute the basic
reproductive number, R0, of individual coffee trees, conditions for local and
global stability under specific conditions, parameter estimation of key
parameters, as well as sensitivity analysis, and numerical experiments under
local and global control strategies for key scenarios.


  \noindent We formulate an age-structured three-staged nonlinear partial
differential equation model that features {\it nonlinear} recidivism to the
infected ({\it infectious}) class from the {\it temporarily} recovered class.
Equilibria are computed, as well as local and global stability of the {\it
infection-free} equilibrium. As a result, a backward-bifurcation exists under
necessary and sufficient conditions. A generalized numerical framework is
established and numerical experiments are explored for two positive solutions
to exist in the {\it infectious} class.


  Data of alien species presences are generally more readily available in
protected than non-protected areas due to higher sampling efforts inside
protected areas. Are the results and conclusions based on analyses of data
collected in protected areas representative of wider non-protected regions? We
address this question by analysing some recently published data of alien plants
in Greece. Mixed effects models were used with alien species presences in 8.25
x 8.25 km cells as dependent variable and the percentage of protected area, as
well as the agricultural and artificial land cover types richness (as
indicators of human presence) as independent variables. In addition, the
spatial cross-correlation between the percentage of protected area and alien
species richness was examined across scales. Results indicated that the
percentage of protected area per cell is a poor predictor of alien species
richness. Spatial analysis indicated that cells with higher percentage of
protected areas have slightly less alien species than cells with lower
percentage of protected areas. This result is likely to be driven by the
overall negative correlation between habitat protection and anthropogenic
activities. Thus, the conclusions deduced by data deriving from protected areas
are likely to hold true for patterns of alien species in non-protected areas
when the human pressures are accounted for.


  The CryoEM single particle imaging method has recently received broad
attention in the field of structural biology for determining the structures of
biological molecules. The structures can be resolved to near-atomic resolutions
after rending a large number of CryoEM images measuring molecules in different
orientations. However, the factors for model resolution need to be further
explored. Here, we provide a theoretical framework in conjunction with
numerical simulations to gauge the influence of several key factors that are
determinant in model resolution. We found that the number of measured
projection images and the quality of each measurement (quantified using average
signal-noise-ratio) can be combined to a single factor, which is dominant to
the constructed model resolution. Furthermore, the intrinsic thermal motion of
the molecules and the defocus levels of the electron microscope both have
significant effects on the model resolution. These effects can be
quantitatively summarized using an analytical formula that provides a
theoretical guideline on structure resolutions for given experimental
measurements.


  Invasion by the red imported fire ant, Solenopsis invicta Buren, has
destructive effects on native biodiversity, agriculture, and public health.
This ant's aggressive foraging behaviour and high reproductive capability have
enabled its establishment of wild populations in most regions into which it has
been imported. An important aspect of eradication is thorough nest monitoring
and destruction during early invasion to prevent range expansion. The question
is: How intense must monitoring be on temporal and spatial scales to eradicate
the fire ant? Assuming that the ant was introduced into a region and that
monitoring was conducted immediately after nest detection in an effort to
detect all other potentially established nests, we developed a mathematical
model to investigate detection rates. Setting the monitoring limit to three
years, the detection rate was maximized when monitoring was conducted shifting
bait trap locations and setting them at intervals of 30 m for each monitoring.
Monitoring should be conducted in a radius of at least 4 km around the source
nest, or wider --depending on how late a nest is found. For ease of
application, we also derived equations for finding the minimum bait interval
required in an arbitrary ant species for thorough monitoring.


  Considerable data and analysis support the detection of one or more
supernovae at a distance of about 50 pc, ~2.6 million years ago. This is
possibly related to the extinction event around that time and is a member of a
series of explosions which formed the Local Bubble in the interstellar medium.
We build on previous work, and propagate the muon flux from supernova-initiated
cosmic rays from the surface to the depths of the ocean. We find that the
radiation dose from the muons will exceed the total present surface dose from
all sources at depths up to a kilometer and will persist for at least the
lifetime of marine megafauna. It is reasonable to hypothesize that this
increase in radiation load may have contributed to a newly documented marine
megafaunal extinction at that time.


  The {\em double-cut-and-join} (DCJ) operation, introduced by Yancopoulos
\emph{et al.}, allows minimum edit distance to be computed by modeling all
possible classical rearrangement operations, such as inversions, fusions,
fissions, translocations, and transpositions, in linear-time between two
genomes. However, there is lack of visualization tool that can effectively
present DCJ operations that will help biologists to use DCJ operation. In this
paper, a new visualization program is introduced, DCJVis, to create a diagram
of each DCJ operation necessary to transform between the genomes of two
distinct organisms by describing a possible sequence of genome graphs based on
the selected gene adjacency on the source genome for the DCJ operation. Our
program is the first visualization tool for DCJ operations using circular
layout. Specifically, the genomes of \textit{Saccharomyces cerevisiae} and
\textit{Candida albicans} are used to demonstrate the functionality of this
program and provide an example of the type of problem this program can solve
for biologists.


  Mutual information and its causal variant, directed information, have been
widely used to quantitatively characterize the performance of biological
sensing and information transduction. However, once coupled with selection in
response to decision-making, the sensing signal could have more or less
evolutionary value than its mutual or directed information. In this work, we
show that an individually sensed signal always has a better fitness value, on
average, than its mutual or directed information. The fitness gain, which
satisfies fluctuation relations (FRs), is attributed to the selection of
organisms in a population that obtain a better sensing signal by chance. A new
quantity, similar to the coarse-grained entropy production in information
thermodynamics, is introduced to quantify the total fitness gain from
individual sensing, which also satisfies FRs. Using this quantity, the
optimizing fitness gain from individual sensing is shown to be related to
fidelity allocations for individual environmental histories. Our results are
supplemented by numerical verifications of FRs, and a discussion on how this
problem is linked to information encoding and decoding.


  The isoelectric point (IEP or pI) is the pH where the net charge on the
molecular ensemble of peptides and proteins is zero. This physical-chemical
property is dependent on protonable/deprotonable sidechains and their pKa
values. Here an pI prediction model is trained from a database of peptide
sequences and pIs using a recurrent neural network (RNN) with long short-term
memory (LSTM) cells. The trained model obtains an RMSE and R$^2$ of 0.28 and
0.95 for the external test set. The model is not based on pKa values, but
prediction of constructed test sequences show similar rankings as already known
pKa values. The prediction depends mostly on the existence of known acidic and
basic amino acids with fine-adjusted based on the neighboring sequence and
position of the charged amino acids in the peptide chain.


  Recently, we have shown that the age-specific prevalence of a disease can be
related to the transition rates in the illness-death model via a partial
differential equation (PDE). In case of a chronic disease, we show that the PDE
can be used to estimate excess mortality from prevalence and incidence.
Applicability of the new method is demonstrated in a simulation and claims data
about diabetes in German men.


  Modern problems of concept annotation associate an object of interest (gene,
individual, text document) with a set of interrelated textual descriptors
(functions, diseases, topics), often organized in concept hierarchies or
ontologies. Most ontologies can be seen as directed acyclic graphs, where nodes
represent concepts and edges represent relational ties between these concepts.
Given an ontology graph, each object can only be annotated by a consistent
subgraph; that is, a subgraph such that if an object is annotated by a
particular concept, it must also be annotated by all other concepts that
generalize it. Ontologies therefore provide a compact representation of a large
space of possible consistent subgraphs; however, until now we have not been
aware of a practical algorithm that can enumerate such annotation spaces for a
given ontology. In this work we propose an algorithm for enumerating consistent
subgraphs of directed acyclic graphs. The algorithm recursively partitions the
graph into strictly smaller graphs until the resulting graph becomes a rooted
tree (forest), for which a linear-time solution is computed. It then combines
the tallies from graphs created in the recursion to obtain the final count. We
prove the correctness of this algorithm and then apply it to characterize four
major biomedical ontologies. We believe this work provides valuable insights
into concept annotation spaces and predictability of ontological annotation.


  Clustered regularly interspaced short palindromic repeats (CRISPR) and
CRISPR-associated proteins (Cas) constitute a multi-functional, constantly
evolving immune system in bacteria and archaea cells. A heritable, molecular
memory is generated of phage, plasmids, or other mobile genetic elements that
attempt to attack the cell. This memory is used to recognize and interfere with
subsequent invasions from the same genetic elements. This versatile prokaryotic
tool has also been used to advance applications in biotechnology. Here we
review a large body of CRISPR-Cas research to explore themes of evolution and
selection, population dynamics, horizontal gene transfer, specific and
cross-reactive interactions, cost and regulation, non-immunological CRISPR
functions that boost host cell robustness, as well as applicable mechanisms for
efficient and specific genetic engineering. We offer future directions that can
be addressed by the physics community. Physical understanding of the CRISPR-Cas
system will advance uses in biotechnology, such as developing cell lines and
animal models, cell labeling and information storage, combatting antibiotic
resistance, and human therapeutics.


  Sustainable yields that are at least 80% of the maximum sustainable yield are
sometimes referred to as pretty good yield (PGY). The range of PGY harvesting
strategies is generally broad and thus leaves room to account for additional
objectives besides high yield. Here, we analyze stage-dependent harvesting
strategies that realize PGY with conservation as a second objective. We show
that (1) PGY harvesting strategies can give large conservation benefits and (2)
equal harvesting rates of juveniles and adults is often a good strategy. These
conclusions are based on trade-off curves between yield and four measures of
conservation that form in two established population models, one age-structured
and one stage-structured model, when considering different harvesting rates of
juveniles and adults. These conclusions hold for a broad range of parameter
settings, though our investigation of robustness also reveals that (3)
predictions of the age-structured model are more sensitive to variations in
parameter values than those of the stage-structured model. Finally, we find
that (4) measures of stability that are often quite difficult to assess in the
field (e.g.~basic reproduction ratio and resilience) are systematically
negatively correlated with impacts on biomass and impact on size structure, so
that these later quantities can provide integrative signals to detect possible
collapses.


  This paper studies lying in a novel context. Previous work has focused on
situations in which people are either fully aware of the economic consequences
of all available actions (e.g., die-under-cup paradigm), or they are uncertain,
but this uncertainty cannot be cleared in any way (e.g., sender-receiver game).
On the contrary, in reality, people oftentimes know that they will have a
chance to lie, they are initially uncertain about the economic consequences of
the available actions, but they can invest resources (e.g., time) to find them
out. Here we capture the essence of this type of situations by means of a novel
decision problem. Two experiments provide evidence of four empirical
regularities regarding the distribution of choices, and suggest that
participants vary along two dimensions: the moral cost of lying, and the cost
of investing time to find out the payoffs associated to the available actions.
Taking inspiration from these observations, we introduce a model that is
consistent with the main empirical results.


  This paper presents a novel optical processing approach for exploring genome
sequences built upon optical correlator for global alignment and extended
DV-curve method for local alignment. To overcome the problem of traditional
DV-curve method for presenting an accurate and simplified output, we propose
HAWPOD, built upon DV- curve method, to analyze genome sequences in five steps:
DNA coding, alignment, noise cancellation, simplification, and modification.
Moreover, all-optical implementation of the HAWPOD method is developed, while
its accuracy is validated through numerical simulations in LUMERICAL FDTD. The
results express the proposed method is much faster than its electrical
counterparts, such as Basic Local Alignment Search Tools.


  Variational auto-encoder frameworks have demonstrated success in reducing
complex nonlinear dynamics in molecular simulation to a single non-linear
embedding. In this work, we illustrate how this non-linear latent embedding can
be used as a collective variable for enhanced sampling, and present a simple
modification that allows us to rapidly perform sampling in multiple related
systems. We first demonstrate our method is able to describe the effects of
force field changes in capped alanine dipeptide after learning a model using
AMBER99. We further provide a simple extension to variational dynamics encoders
that allows the model to be trained in a more efficient manner on larger
systems by encoding the outputs of a linear transformation using time-structure
based independent component analysis (tICA). Using this technique, we show how
such a model trained for one protein, the WW domain, can efficiently be
transferred to perform enhanced sampling on a related mutant protein, the GTT
mutation. This method shows promise for its ability to rapidly sample related
systems using a single transferable collective variable and is generally
applicable to sets of related simulations, enabling us to probe the effects of
variation in increasingly large systems of biophysical interest.


  This paper gives a new flavor of what Peter Jagers and his co-authors call
`the path to extinction'. In a neutral population with constant size $N$, we
assume that each individual at time $0$ carries a distinct type, or allele. We
consider the joint dynamics of these $N$ alleles, for example the dynamics of
their respective frequencies and more plainly the nonincreasing process
counting the number of alleles remaining by time $t$. We call this process the
extinction process. We show that in the Moran model, the extinction process is
distributed as the process counting (in backward time) the number of common
ancestors to the whole population, also known as the block counting process of
the $N$-Kingman coalescent. Stimulated by this result, we investigate: (1)
whether it extends to an identity between the frequencies of blocks in the
Kingman coalescent and the frequencies of alleles in the extinction process,
both evaluated at jump times; (2) whether it extends to the general case of
$\Lambda$-Fleming-Viot processes.


  Recent progress in microdissection and in DNA sequencing has enabled
subsampling of multi-focal cancers in organs such as the liver in several
hundred spots, helping to determine the pattern of mutations in each of these
spots. This has led to the construction of genealogies of the primary,
secondary, tertiary and so forth, foci of the tumor. These studies have led to
diverse conclusions concerning the Darwinian (selective) or neutral evolution
in cancer. Mathematical models of development of multifocal tumors have been
developed to support these claims. We report a model of development of a
multifocal tumor, which is a mathematically rigorous refinement of a model of
Ling et al. (2015). Guided by numerical studies and simulations, we show that
the rigorous model, in the form of an infinite-type branching process, displays
distributions of tumors size which have heavy tails and moments that become
infinite in finite time. To demonstrate these points, we obtain bounds on the
tails of the distributions of the process and infinite-series expression for
the first moments. In addition to its inherent mathematical interest, the model
is corroborated by recent reports of apparent super-exponential growth in
cancer metastases.


  Life can be viewed as a localized chemical system that sits on, or in the
basin of attraction of, a metastable dynamical attractor state that remains out
of equilibrium with the environment. Such a view of life allows that new living
states can arise through chance changes in local chemical concentration
(=mutations) that move points in space into the basin of attraction of a life
state - the attractor being an autocatalytic sets whose essential (=keystone)
species are produced at a higher rate than they are lost to the environment by
diffusion, such that growth in expected. This conception of life yields several
new insights and conjectures. (1) This framework suggests that the first new
life states to arise are likely at interfaces where the rate of diffusion of
keystone species is tied to a low-diffusion regime, while precursors and waste
products diffuse at a higher rate. (2) There are reasons to expect that once
the first life state arises, most likely on a mineral surface, additional
mutations will generate derived life states with which the original state will
compete. (3) I propose that in the resulting adaptive process there is a
general tendency for higher complexity life states (i.e., ones that are further
from being at equilibrium with the environment) to dominate a given mineral
surface. (4) The framework suggests a simple and predictable path by which
cells evolve and provides pointers on why such cells are likely to acquire
particulate inheritance. Overall, the dynamical systems theoretical framework
developed provides an integrated view of the origin and early evolution of life
and supports novel empirical approaches.


  Second generation sequencing technologies are being increasingly used for
genetic association studies, where the main research interest is to identify
sets of genetic variants that contribute to various phenotype. The phenotype
can be univariate disease status, multivariate responses and even
high-dimensional outcomes. Considering the genotype and phenotype as two
complex objects, this also poses a general statistical problem of testing
association between complex objects. We here proposed a similarity-based test,
generalized similarity U (GSU), that can test the association between complex
objects. We first studied the theoretical properties of the test in a general
setting and then focused on the application of the test to sequencing
association studies. Based on theoretical analysis, we proposed to use
Laplacian kernel based similarity for GSU to boost power and enhance
robustness. Through simulation, we found that GSU did have advantages over
existing methods in terms of power and robustness. We further performed a whole
genome sequencing (WGS) scan for Alzherimer Disease Neuroimaging Initiative
(ADNI) data, identifying three genes, APOE, APOC1 and TOMM40, associated with
imaging phenotype. We developed a C++ package for analysis of whole genome
sequencing data using GSU. The source codes can be downloaded at
https://github.com/changshuaiwei/gsu.


  The duration, type and structure of connections between individuals in
real-world populations play a crucial role in how diseases invade and spread.
Here, we incorporate the aforementioned heterogeneities into a model by
considering a dual-layer static-dynamic multiplex network. The static network
layer affords tunable clustering and describes an individual's permanent
community structure. The dynamic network layer describes the transient
connections an individual makes with members of the wider population by
imposing constant edge rewiring. We follow the edge-based compartmental
modelling approach to derive equations describing the evolution of a
susceptible - infected - recovered (SIR) epidemic spreading through this
multiplex network of individuals. We derive the basic reproduction number,
measuring the expected number of new infectious cases caused by a single
infectious individual in an otherwise susceptible population. We validate model
equations by showing convergence to pre-existing edge-based compartmental model
equations in limiting cases and by comparison with stochastically simulated
epidemics. We explore the effects of altering model parameters and multiplex
network attributes on resultant epidemic dynamics. We validate the basic
reproduction number by plotting its value against associated final epidemic
sizes measured from simulation and predicted by model equations for a number of
setups. Further, we explore the effect of varying individual model parameters
on the basic reproduction number. We conclude with a discussion of the
significance and interpretation of the model and its relation to existing
research literature. We highlight intrinsic limitations and potential
extensions of the present model and outline future research considerations,
both experimental and theoretical.


  Protein pattern formation is essential for the spatial organization of many
intracellular processes like cell division, flagellum positioning, and
chemotaxis. A prominent example of intracellular patterns are the oscillatory
pole-to-pole oscillations of Min proteins in \textit{E. coli} whose biological
function is to ensure precise cell division. Cell polarization, a prerequisite
for processes such as stem cell differentiation and cell polarity in yeast, is
also mediated by a diffusion-reaction process. More generally, these functional
modules of cells serve as model systems for self-organization, one of the core
principles of life. Under which conditions spatio-temporal patterns emerge, and
how these patterns are regulated by biochemical and geometrical factors are
major aspects of current research. Here we review recent theoretical and
experimental advances in the field of intracellular pattern formation, focusing
on general design principles and fundamental physical mechanisms.


  Vegetation patterns are abundant in arid and semiarid ecosystems, but how
they form remains unclear. One of the most extended theories lies in the
existence of scale-dependent feedbacks (SDF) in plant-to-plant and plant-water
interactions. Short distances are dominated by facilitative interactions,
whereas competitive interactions dominate at larger scales. These feedbacks
shape spatially inhomogeneous distributions of water that ultimately drive the
emergence of patterns of vegetation. Even though the presence of facilitative
and competitive interactions is clear, they are often hard to disentangle in
the field, and therefore their relevance in vegetation pattern formation is
still disputable. Here, we review the biological processes that have been
proposed to explain pattern formation in arid ecosystems and how they have been
implemented in mathematical models. We conclude by discussing the existence of
similar structures in different biological and physical systems.


  Alzheimer's Disease (AD) is a neurodegenerative disorder that lacks effective
treatment options. Anti-amyloid beta (ABeta) antibodies are the leading drug
candidates to treat AD, but the results of clinical trials have been
disappointing. Introducing rational mutations into anti-ABeta antibodies to
increase their effectiveness is a way forward, but the path to take is unclear.
In this study, we demonstrate the use of computational fragment-based docking
and MMPBSA binding free energy calculations in the analysis of anti-ABeta
antibodies for rational drug design efforts. Our fragment-based docking method
successfully predicted the emergence of the common EFRH epitope, MD simulations
coupled with MMPBSA binding free energy calculations were used to analyze
scenarios described in prior studies, and we introduced rational mutations into
PFA1 to improve its calculated binding affinity towards the pE3-ABeta3-8 form
of ABeta. Two out of four proposed mutations stabilized binding. Our study
demonstrates that a computational approach may lead to an improved drug
candidate for AD in the future.


  In population genetics, diffusions on the unit interval are often used to
model the frequency path of an allele. In this setting we derive approximations
for fixation probabilities, expected hitting times and the expected
site-frequency-spectrum under small frequency-dependent selection.
Specifically, we rederive and extend the one-third rule of evolutionary game
theory (Nowak et al., 2004) and effects of stochastic slowdown (Altrock and
Traulsen, 2009). Since similar effects are of interest in other application
areas, we formulate our results for general one-dimensional diffusions.


  There is a growing need for unbiased clustering methods, ideally automated.
We have developed a topology-based analysis tool called Two-Tier Mapper (TTMap)
to detect subgroups in global gene expression datasets and identify their
distinguishing features. First, TTMap discerns and adjusts for highly variable
features in the control group and identifies outliers. Second, the deviation of
each test sample from the control group in a high-dimensional space is computed
and the test samples are clustered in a global and local network using a new
topological algorithm based on Mapper. Validation of TTMap on both synthetic
and biological datasets shows that it outperforms current clustering methods in
sensitivity and stability; clustering is not affected by removal of samples
from the control group, choice of normalization nor subselection of data. There
is no user induced bias because all parameters are data-driven. Datasets can
readily be combined into one analysis. TTMap reveals hitherto undetected gene
expression changes in mouse mammary glands related to hormonal changes during
the estrous cycle. This illustrates the ability to extract information from
highly variable biological samples and its potential for personalized medicine.


  Deterministic models are developed for the spatial spread of epidemic
diseases in geographical settings. The models are focused on outbreaks that
arise from a small number of infected hosts imported into sub-regions of the
geographical settings. The goal is to understand how spatial heterogeneity
influences the transmission dynamics of the susceptible and infected
populations. The models consist of systems of partial differential equations
with diffusion terms describing the spatial spread of the underlying microbial
infectious agents. The model is compared with real data from seasonal influenza
epidemics in Puerto Rico.


  Many living and non-living complex systems can be modeled and understood as
collective systems made of heterogeneous components that self-organize and
generate nontrivial morphological structures and behaviors. This chapter
presents a brief overview of our recent effort that investigated various
aspects of such morphogenetic collective systems. We first propose a
theoretical classification scheme that distinguishes four complexity levels of
morphogenetic collective systems based on the nature of their components and
interactions. We conducted a series of computational experiments using a
self-propelled particle swarm model to investigate the effects of (1)
heterogeneity of components, (2) differentiation/re-differentiation of
components, and (3) local information sharing among components, on the
self-organization of a collective system. Results showed that (a) heterogeneity
of components had a strong impact on the system's structure and behavior, (b)
dynamic differentiation/re-differentiation of components and local information
sharing helped the system maintain spatially adjacent, coherent organization,
(c) dynamic differentiation/re-differentiation contributed to the development
of more diverse structures and behaviors, and (d) stochastic re-differentiation
of components naturally realized a self-repair capability of self-organizing
morphologies. We also explored evolutionary methods to design novel
self-organizing patterns, using interactive evolutionary computation and
spontaneous evolution within an artificial ecosystem. These self-organizing
patterns were found to be remarkably robust against dimensional changes from 2D
to 3D, although evolution worked efficiently only in 2D settings.


  We study the behavior of five proteins at the air-water and oil-water
interfaces by all-atom molecular dynamics. The proteins are found to get
distorted when pinned to the interface. This behavior is consistent with the
phenomenological way of introducing the interfaces in a coarse-grained model
through a force that depends on the hydropathy indices of the residues.
Proteins couple to the oil-water interface stronger than to the air- water one.
They diffuse slower at the oil-water interface but do not depin from it,
whereas depinning events are observed at the other interface. The reduction of
the disulfide bonds slows the diffusion down.


  The probability that the frequency of a particular trait will eventually
become unity, the so-called fixation probability, is a central issue in the
study of population evolution. Its computation, once we are given a stochastic
finite population model without mutations and a (possibly frequency dependent)
fitness function, is straightforward and it can be done in several ways.
Nevertheless, despite the fact that the fixation probability is an important
macroscopic property of the population, its precise knowledge does not give any
clear information about the interaction patterns among individuals in the
population. Here we address the inverse problem: From a given fixation pattern
and population size, we want to infer what is the game being played by the
population. This is done by first exploiting the framework developed in FACC
Chalub and MO Souza, J. Math. Biol. 75: 1735, 2017., which yields a fitness
function that realises this fixation pattern in the Wright-Fisher model. This
fitness function always exists, but it is not necessarily unique. Subsequently,
we show that any such fitness function can be approximated, with arbitrary
precision, using $d$-player game theory, provided $d$ is large enough. The
pay-off matrix that emerges naturally from the approximating game will provide
useful information about the individual interaction structure that is not
itself apparent in the fixation pattern. We present extensive numerical support
for our conclusions.


  Under many in vitro conditions, some small viruses spontaneously encapsidate
a single stranded (ss) RNA into a protein shell called the capsid. While viral
RNAs are found to be compact and highly branched because of long distance
base-pairing between nucleotides, recent experiments reveal that in a
head-to-head competition between a ssRNA with no secondary or higher order
structure and a viral RNA, the capsid proteins preferentially encapsulate the
linear polymer! In this paper, we study the impact of genome stiffness on the
encapsidation free energy of the complex of RNA and capsid proteins. We show
that an increase in effective chain stiffness because of base-pairing could be
the reason why under certain conditions linear chains have an advantage over
branched chains when it comes to encapsidation efficiency. While branching
makes the genome more compact, RNA base-pairing increases the effective Kuhn
length of the RNA molecule, which could result in an increase of the free
energy of RNA confinement, that is, the work required to encapsidate RNA, and
thus less efficient packaging.


  Given multi-platform genome data with prior knowledge of functional gene
sets, how can we extract interpretable latent relationships between patients
and genes? More specifically, how can we devise a tensor factorization method
which produces an interpretable gene factor matrix based on gene set
information while maintaining the decomposition quality and speed? We propose
GIFT, a Guided and Interpretable Factorization for Tensors. GIFT provides
interpretable factor matrices by encoding prior knowledge as a regularization
term in its objective function. Experiment results demonstrate that GIFT
produces interpretable factorizations with high scalability and accuracy, while
other methods lack interpretability. We apply GIFT to the PanCan12 dataset, and
GIFT reveals significant relations between cancers, gene sets, and genes, such
as influential gene sets for specific cancer (e.g., interferon-gamma response
gene set for ovarian cancer) or relations between cancers and genes (e.g., BRCA
cancer - APOA1 gene and OV, UCEC cancers - BST2 gene).


  In recent years, the number of biomedical publications has steadfastly grown,
resulting in a rich source of untapped new knowledge. Most biomedical facts are
however not readily available, but buried in the form of unstructured text, and
hence their exploitation requires the time-consuming manual curation of
published articles. Here we present INtERAcT, a novel approach to extract
protein-protein interactions from a corpus of biomedical articles related to a
broad range of scientific domains in a completely unsupervised way. INtERAcT
exploits vector representation of words, computed on a corpus of domain
specific knowledge, and implements a new metric that estimates an interaction
score between two molecules in the space where the corresponding words are
embedded. We demonstrate the power of INtERAcT by reconstructing the molecular
pathways associated to 10 different cancer types using a corpus of
disease-specific articles for each cancer type. We evaluate INtERAcT using
STRING database as a benchmark, and show that our metric outperforms currently
adopted approaches for similarity computation at the task of identifying known
molecular interactions in all studied cancer types. Furthermore, our approach
does not require text annotation, manual curation or the definition of semantic
rules based on expert knowledge, and hence it can be easily and efficiently
applied to different scientific domains. Our findings suggest that INtERAcT may
increase our capability to summarize the understanding of a specific disease
using the published literature in an automated and completely unsupervised
fashion.


  In this paper a novel biclustering algorithm based on artificial intelligence
(AI) is introduced. The method called EBIC aims to detect biologically
meaningful, order-preserving patterns in complex data. The proposed algorithm
is probably the first one capable of discovering with accuracy exceeding 50%
multiple complex patterns in real gene expression datasets. It is also one of
the very few biclustering methods designed for parallel environments with
multiple graphics processing units (GPUs). We demonstrate that EBIC outperforms
state-of-the-art biclustering methods, in terms of recovery and relevance, on
both synthetic and genetic datasets. EBIC also yields results over 12 times
faster than the most accurate reference algorithms. The proposed algorithm is
anticipated to be added to the repertoire of unsupervised machine learning
algorithms for the analysis of datasets, including those from large-scale
genomic studies.


  A growing number of RNA sequences are now known to have distributions of
multiple stable sequences. Recent algorithms use the list of nucleotides in a
sequence and auxiliary experimental data to predict such distributions.
Although the algorithms are largely successful in identifying a distribution's
constituent structures, it remains challenging to recover their relative
weightings. In this paper, we quantify this issue using a total variation
distance. Then, we prove under a Nussinov-Jacobson model that a large
proportion of RNA structure pairs cannot be jointly reconstructed with low
total variation distance. Finally, we characterize the uncertainty in
predicting conformational ratios by analyzing the amount of information in the
auxiliary data.


  This paper presents a comprehensive list of the scientific articles of Giulio
Fermi (1936-1997), son of the Italian-American physicist Enrico Fermi,
published between 1962 and 1997. The initial research activity of Giulio was
concerned with virology and biological cybernetics while, from 1975 onward, his
work was completely devoted to protein crystallography. The crystallographic
research was carried out in collaboration with Nobel laureate Max Perutz at the
Medical Research Council (MRC) Laboratory of Molecular Biology in Cambridge
(United Kingdom). A short biography of Giulio (Judd) Fermi appears inside John
Finch's book A Nobel Fellow on Every Floor: A History of the Medical Research
Council Laboratory of Molecular Biology published by the MRC in 2008.


  In this work, we analyze and model the cholera epidemics that affected South
Sudan, the newest country in the world, during 2014 and 2015. South Sudan
possibly represents one of the most difficult context in which adapt the
deterministic mathematical cholera model, due to the unstable social and
political situation that clearly affects the fluxes of people and the sanitary
conditions, increasing the risk of large outbreaks. Despite the limitation of a
static gravity model in describing the chaotic human mobility of South Sudan,
the SIRB model, calibrated with a data assimilation technique (Ensemble Kalman
Filter), retrieves the epidemic dynamics in the counties with the largest
number of infected cases, showing the potentiality of the methodology in
forecasting future outbreaks.


  Genes are frequently regulated in complex manners, necessitating modelling
approaches which go beyond simple (linear) gene-to-gene interactions and
address the modularity of cis-regulatory regions and alternate transcription
initiation sites. In particular, sharp expression patterns (peaks or stripes)
indicate that gene regulation involves nonlinear transcription factor kinetics
(beyond first-order). We propose a methodology for approaching this problem,
using the example of the multiple cis-regulatory modules (CRMs) and two
transcripts (P1 and P2) found in the Drosophila hunchback (hb) and (CD1 and
CD2)in Kruppel (Kr) genes, the genes expressed along the anteroposterior axis
of the embryo. The positional characteristics of the mRNA expression patterns
are studied for their sensitivity to the variation of the input factors and
initial conditions.


  A central claim in modern network science is that real-world networks are
typically "scale free," meaning that the fraction of nodes with degree $k$
follows a power law, decaying like $k^{-\alpha}$, often with $2 < \alpha < 3$.
However, empirical evidence for this belief derives from a relatively small
number of real-world networks. We test the universality of scale-free structure
by applying state-of-the-art statistical tools to a large corpus of nearly 1000
network data sets drawn from social, biological, technological, and
informational sources. We fit the power-law model to each degree distribution,
test its statistical plausibility, and compare it via a likelihood ratio test
to alternative, non-scale-free models, e.g., the log-normal. Across domains, we
find that scale-free networks are rare, with only 4% exhibiting the
strongest-possible evidence of scale-free structure and 52% exhibiting the
weakest-possible evidence. Furthermore, evidence of scale-free structure is not
uniformly distributed across sources: social networks are at best weakly scale
free, while a handful of technological and biological networks can be called
strongly scale free. These results undermine the universality of scale-free
networks and reveal that real-world networks exhibit a rich structural
diversity that will likely require new ideas and mechanisms to explain.


  Resource diffusion is an ubiquitous phenomenon, but how it impacts epidemic
spreading has received little study. We propose a model that couples epidemic
spreading and resource diffusion in multiplex networks. The spread of disease
in a physical contact layer and the recovery of the infected nodes are both
strongly dependent upon resources supplied by their counterparts in the social
layer. The generation and diffusion of resources in the social layer are in
turn strongly dependent upon the state of the nodes in the physical contact
layer. Resources diffuse preferentially or randomly in this model. To quantify
the degree of preferential diffusion, a bias parameter that controls the
resource diffusion is proposed. We conduct extensive simulations and find that
the preferential resource diffusion can change phase transition type of the
fraction of infected nodes. When the degree of interlayer correlation is below
a critical value, increasing the bias parameter changes the phase transition
from double continuous to single continuous. When the degree of interlayer
correlation is above a critical value, the phase transition changes from
multiple continuous to first discontinuous and then to hybrid. We find
hysteresis loops in the phase transition. We also find that there is an optimal
resource strategy at each fixed degree of interlayer correlation where the
threshold reaches a maximum and under which the disease can be maximally
suppressed. In addition, the optimal controlling parameter increases as the
degree of inter-layer correlation increases.


  A stable population network is hard to interrupt without any ecological
consequences. A communication blockage between patches may destabilize the
populations in the ecological network. This work deals with the construction of
a safe cut passing through metapopulations habitat such that populations remain
stable. We combine the dynamical system stability analysis with graph
partitioning algorithms in our approach to the problem. It finds such a safe
construction, when one exists, provided the algebraic connectivity of the graph
components is stronger than all the spatially local instabilities in the
respective components. The dynamics of the populations on the spatially
discrete patches (graph nodes) and their spatial communication with other
patches is modeled as a reaction-diffusion system. By reversing the
Turing-instability idea the stability conditions of the partitioned system are
found to depend on local dynamics of the metapopulations and the Fiedler value
of the Laplacian matrix of the graph. This leads to the necessary and
sufficient conditions for removal of the graph edges subject to the stability
of the partitioned graph networks. An heuristic bisection graph partitioning
algorithm has been proposed and examples illustrate the theoretical result.


  The function of proteins arises from cooperative interactions and
rearrangements of their amino acids, which exhibit large-scale dynamical modes.
Long-range correlations have also been revealed in protein sequences, and this
has motivated the search for physical links between the observed genetic and
dynamic cooperativity. We outline here a simplified theory of protein, which
relates sequence correlations to physical interactions and to the emergence of
mechanical function. Our protein is modeled as a strongly-coupled amino acid
network whose interactions and motions are captured by the mechanical
propagator, the Green function. The propagator describes how the gene
determines the connectivity of the amino acids, and thereby the transmission of
forces. Mutations introduce localized perturbations to the propagator which
scatter the force field. The emergence of function is manifested by a
topological transition when a band of such perturbations divides the protein
into subdomains. We find that epistasis -- the interaction among mutations in
the gene -- is related to the nonlinearity of the Green function, which can be
interpreted as a sum over multiple scattering paths. We apply this mechanical
framework to simulations of protein evolution, and observe long-range epistasis
which facilitates collective functional modes.


  The transition from single-cell to multicellular behavior is important in
early development but rarely studied. The starvation-induced aggregation of the
social amoeba Dictyostelium discoideum into a multicellular slug is known to
result from single-cell chemotaxis towards emitted pulses of cyclic adenosine
monophosphate (cAMP). However, how exactly do transient short-range chemical
gradients lead to coherent collective movement at a macroscopic scale? Here, we
use a multiscale model verified by quantitative microscopy to describe
wide-ranging behaviors from chemotaxis and excitability of individual cells to
aggregation of thousands of cells. To better understand the mechanism of
long-range cell-cell communication and hence aggregation, we analyze cell-cell
correlations, showing evidence for self-organization at the onset of
aggregation (as opposed to following a leader cell). Surprisingly, cell
collectives, despite their finite size, show features of criticality known from
phase transitions in physical systems. Application of external cAMP
perturbations in our simulations near the sensitive critical point allows
steering cells into early aggregation and towards certain locations but not
once an aggregation center has been chosen.


  We give correct explicit formulas for the probabilities of rooted binary
trees and cladograms under Ford's $\alpha$-model.


  T cell receptors (TCRs) bind foreign or self-peptides attached to major
histocompatibility complex (MHC) molecules, and the strength of this
interaction determines T cell activation. Optimizing the ability of T cells to
recognize a diversity of foreign peptides yet be tolerant of self-peptides is
crucial for the adaptive immune system to properly function. This is achieved
by selection of T cells in the thymus, where immature T cells expressing
unique, stochastically generated TCRs interact with a large number of
self-peptide-MHC; if a TCR does not bind strongly enough to any
self-peptide-MHC, or too strongly with at least one self-peptide-MHC, the T
cell dies. Past theoretical work cast thymic selection as an extreme value
problem, and characterized the statistical enrichment or depletion of amino
acids in the post-selection TCR repertoire, showing how T cells are selected to
be able to specifically recognize peptides derived from diverse pathogens, yet
have limited self-reactivity. Here, we investigate how the degree of enrichment
is modified by nonuniform contacts that a TCR makes with peptide-MHC.
Specifically, we were motivated by recent experiments showing that amino acids
at certain positions of a TCR sequence have large effects on thymic selection
outcomes, and crystal structure data that reveal a nonuniform contact profile
between a TCR and its peptide-MHC ligand. Using a representative TCR contact
profile as an illustration, we show via simulations that the degree of
enrichment now varies by position according to the contact profile, and,
importantly, it depends on the implementation of nonuniform contacts during
thymic selection. We explain these nontrivial results analytically. Our study
has implications for understanding the selection forces that shape the
functionality of the post-selection TCR repertoire.


  Self-sustaining autocatalytic networks play a central role in living systems,
from metabolism at the origin of life, simple RNA networks, and the modern
cell, to ecology and cognition. A collectively autocatalytic network that can
be sustained from an ambient food set is also referred to more formally as a
`Reflexively Autocatalytic F-generated' (RAF) set. In this paper, we first
investigate a simplified setting for studying RAFs, which are nevertheless
relevant to real biochemistry and allows for a more exact mathematical analysis
based on graph-theoretic concepts. This, in turn, allows for the development of
efficient (polynomial-time) algorithms for questions that are computationally
NP-hard in the general RAF setting. We then show how this simplified setting
for RAF systems leads naturally to a more general notion of RAFs that are
`generative' (they can be built up from simpler RAFs) and for which efficient
algorithms carry over to this more general setting. Finally, we show how
classical RAF theory can be extended to deal with ensembles of catalysts as
well as the assignment of rates to reactions according to which catalysts (or
combinations of catalysts) are available.


  The mitogen-activated protein kinase (MAPK) signaling cascade, an
evolutionarily conserved motif present in all eukaryotic cells, is involved in
coordinating critical cell-fate decisions, regulating protein synthesis, and
mediating learning and memory. While the steady-state behavior of the pathway
stimulated by a time-invariant signal is relatively well-understood, we show
using a computational model that it exhibits a rich repertoire of transient
adaptive responses to changes in stimuli. When the signal is switched on, the
response is characterized by long-lived modulations in frequency as well as
amplitude. On withdrawing the stimulus, the activity decays over timescales
much longer than that of phosphorylation-dephosphorylation processes,
exhibiting reverberations characterized by repeated spiking in the activated
MAPK concentration. The long-term persistence of such post-stimulus activity
suggests that the cascade retains memory of the signal for a significant
duration following its removal, even in the absence of any explicit feedback or
cross-talk with other pathways. We find that the molecular mechanism underlying
this behavior is related to the existence of distinct relaxation rates for the
different cascade components. This results in the imbalance of fluxes between
different layers of the cascade, with the repeated reuse of activated kinases
as enzymes when they are released from sequestration in complexes leading to
one or more spike events following the removal of the stimulus. The persistent
adaptive response reported here, indicative of a cellular "short-term" memory,
suggests that this ubiquitous signaling pathway plays an even more central role
in information processing by eukaryotic cells.


  Gene regulatory networks are powerful abstractions of biological systems.
Since the advent of high-throughput measurement technologies in biology in the
late 90s, reconstructing the structure of such networks has been a central
computational problem in systems biology. While the problem is certainly not
solved in its entirety, considerable progress has been made in the last two
decades, with mature tools now available. This chapter aims to provide an
introduction to the basic concepts underpinning network inference tools,
attempting a categorisation which highlights commonalities and relative
strengths. While the chapter is meant to be self-contained, the material
presented should provide a useful background to the later, more specialised
chapters of this book.


  Global coevolutionary models of homologous protein families, as constructed
by direct coupling analysis (DCA), have recently gained popularity in
particular due to their capacity to accurately predict residue-residue contacts
from sequence information alone, and thereby to facilitate tertiary and
quaternary protein structure prediction. More recently, they have also been
used to predict fitness effects of amino-acid substitutions in proteins, and to
predict evolutionary conserved protein-protein interactions. These models are
based on two currently unjustified hypotheses: (a) correlations in the
amino-acid usage of different positions are resulting collectively from
networks of direct couplings; and (b) pairwise couplings are sufficient to
capture the amino-acid variability. Here we propose a highly precise inference
scheme based on Boltzmann-machine learning, which allows us to systematically
address these hypotheses. We show how correlations are built up in a highly
collective way by a large number of coupling paths, which are based on the
protein's three-dimensional structure. We further find that pairwise
coevolutionary models capture the collective residue variability across
homologous proteins even for quantities which are not imposed by the inference
procedure, like three-residue correlations, the clustered structure of protein
families in sequence space or the sequence distances between homologs. These
findings strongly suggest that pairwise coevolutionary models are actually
sufficient to accurately capture the residue variability in homologous protein
families.


  Certain retroviruses, including HIV, insert their DNA in a non-random
fraction of the host genome via poorly understood selection mechanisms. Here,
we develop a biophysical model for retroviral integrations as stochastic and
quasi-equilibrium topological reconnections between polymers. We discover that
physical effects, such as DNA accessibility and elasticity, play important and
universal roles in this process. Our simulations predict that integration is
favoured within nucleosomal and flexible DNA, in line with experiments, and
that these biases arise due to competing energy barriers associated with DNA
deformations. By considering a long chromosomal region in human T-cells during
interphase, we discover that at these larger scales integration sites are
predominantly determined by chromatin accessibility. Finally, we propose and
solve a reaction-diffusion problem that recapitulates the distribution of HIV
hot-spots within T-cells. With few generic assumptions, our model can
rationalise experimental observations and identifies previously unappreciated
physical contributions to retroviral integration site selection.


  The classical idea of evolutionarily stable strategy (ESS) modeling animal
behavior does not involve any spatial dependence. We considered a spatial
Hawk-Dove game played by animals in a patchy environment with wrap around
boundaries. We posit that each site contains the same number of individuals. An
evolution equation for analyzing the stability of the ESS is found as the mean
dynamics of the classical frequency dependent Moran process coupled via
migration and nonlocal payoff calculation in 1D and 2D habitats. The linear
stability analysis of the model is performed and conditions to observe spatial
patterns are investigated. For the nearest neighbor interactions (including von
Neumann and Moore neighborhoods in 2D) we concluded that it is possible to
destabilize the ESS of the game and observe pattern formation when the
dispersal rate is small enough. We numerically investigate the spatial patterns
arising from the replicator equations coupled via nearest neighbor payoff
calculation and dispersal.


  We study the solvation free energy of two different conformations (helix and
extended) of two different peptides (deca-alanine and deca-glycine) in two
different solvents (water and aqueous guanidinium chloride, GdmCl). The free
energies are obtained using the quasichemical organization of the potential
distribution theorem, an approach that naturally provides the repulsive
(solvophobic or cavity) and attractive (solvophilic) contributions to
solvation. The solvophilic contribution is further parsed into a chemistry
contribution arising from solute interaction with the solvent in the first
solvation shell and a long-range contribution arising from non-specific
interactions between the solute and the solvent beyond the first solvation
shell. The cavity contribution is obtained for two different envelopes,
$\Sigma_{SE}$ which theory identifies as the solvent excluded volume and a
larger envelope ($\Sigma_G$) beyond which solute-solvent interactions are
Gaussian. For both envelopes, the cavity contribution in water is proportional
to the surface area of the envelope. The same does not hold for GdmCl(aq),
revealing limitations of using molecular area to assess solvation energetics,
especially in mixed solvents. The $\Sigma_G$-cavity contribution predicts that
GdmCl(aq) should favor the more compact state, contrary to the role of GdmCl in
unfolding proteins. The chemistry contribution attenuates this effect, but
still the net local (chemistry plus $\Sigma_G$-packing) contribution is
inadequate in capturing the role of GdmCl. With the inclusion of the long-range
contribution, which is dominated by van~der~Waals interaction, aqueous GdmCl
favors the extended conformation over the compact conformation. Our finding
emphasizes the importance of weak, but attractive, long-range dispersion
interactions in protein solution thermodynamics.


  ATP synthases utilize a proton motive force to synthesize ATP. In reverse,
these membrane-embedded enzymes can also hydrolyze ATP to pump protons over the
membrane. To prevent wasteful ATP hydrolysis, distinct control mechanisms exist
for ATP synthases in bacteria, archaea, chloroplasts and mitochondria.
Single-molecule F\"orster resonance energy transfer (smFRET) demonstrated that
the C-terminus of the rotary subunit epsilon in the Escherichia coli enzyme
changes its conformation to block ATP hydrolysis. Previously we investigated
the related conformational changes of subunit F of the A1AO-ATP synthase from
the archaeon Methanosarcina mazei G\"o1. Here, we analyze the lifetimes of
fluorescence donor and acceptor dyes to distinguish between smFRET signals for
conformational changes and potential artefacts.


  This paper is a follow-up work about the artificial ecosystem model: number
soup (Liu and Sumpter, J. Royal Soc. Interface, 2017). It elaborates more
details about this model and points out future directions.


  Whereas the relationship between criticality of gene regulatory networks
(GRNs) and dynamics of GRNs at a single cell level has been vigorously studied,
the relationship between the criticality of GRNs and system properties at a
higher level has remained unexplored. Here we aim at revealing a potential role
of criticality of GRNs at a multicellular level which are hard to uncover
through the single-cell-level studies, especially from an evolutionary
viewpoint. Our model simulated the growth of a cell population from a single
seed cell. All the cells were assumed to have identical GRNs. We induced
genetic perturbations to the GRN of the seed cell by adding, deleting, or
switching a regulatory link between a pair of genes. From numerical
simulations, we found that the criticality of GRNs facilitated the formation of
nontrivial morphologies when the GRNs were critical in the presence of the
evolutionary perturbations. Moreover, the criticality of GRNs produced
topologically homogenous cell clusters by adjusting the spatial arrangements of
cells, which led to the formation of nontrivial morphogenetic patterns. Our
findings corresponded to an epigenetic viewpoint that heterogeneous and complex
features emerge from homogeneous and less complex components through the
interactions among them. Thus, our results imply that highly structured tissues
or organs in morphogenesis of multicellular organisms might stem from the
criticality of GRNs.


  Genetic robustness, the preservation of evolved phenotypes against genotypic
mutations, is one of the central concepts in evolution. In recent years a large
body of work has focused on the origins, mechanisms, and consequences of
robustness in a wide range of biological systems. In particular, research on
ncRNAs studied the ability of sequences to maintain folded structures against
single-point mutations. In these studies, the structure is merely a reference.
However, recent work revealed evidence that structure itself contributes to the
genetic robustness of ncRNAs. We follow this line of thought and consider
sequence-structure pairs as the unit of evolution and introduce the spectrum of
inverse folding rates (IFR-spectrum) as a measurement of genetic robustness.
Our analysis of the miRNA let-7 family captures key features of
structure-modulated evolution and facilitates the study of robustness against
multiple-point mutations.


  In this review we address to what extent computational techniques can augment
our ability to predict toxicity. The first section provides a brief history of
empirical observations on toxicity dating back to the dawn of Sumerian
civilization. Interestingly, the concept of dose emerged very early on, leading
up to the modern emphasis on kinetic properties, which in turn encodes the
insight that toxicity is not solely a property of a compound but instead
depends on the interaction with the host organism. The next logical step is the
current conception of evaluating drugs from a personalized medicine
point-of-view. We review recent work on integrating what could be referred to
as classical pharmacokinetic analysis with emerging systems biology approaches
incorporating multiple omics data. These systems approaches employ advanced
statistical analytical data processing complemented with machine learning
techniques and use both pharmacokinetic and omics data. We find that such
integrated approaches not only provide improved predictions of toxicity but
also enable mechanistic interpretations of the molecular mechanisms
underpinning toxicity and drug resistance. We conclude the chapter by
discussing some of the main challenges, such as how to balance the inherent
tension between the predictive capacity of models, which in practice amounts to
constraining the number of features in the models versus allowing for rich
mechanistic interpretability, i.e. equipping models with numerous molecular
features. This challenge also requires patient-specific predictions on
toxicity, which in turn requires proper stratification of patients as regards
how they respond, with or without adverse toxic effects. In summary, the
transformation of the ancient concept of dose is currently successfully
operationalized using rich integrative data encoded in patient-specific models.


  We study the phase diagram of a minority game where three classes of agents
are present. Two types of agents play a risk-loving game that we model by the
standard Snowdrift Game. The behaviour of the third type of agents is coded by
{\em indifference} w.r.t. the game at all: their dynamics is designed to
account for risk-aversion as an innovative behavioral gambit. From this point
of view, the choice of this solitary strategy is enhanced when innovation
starts, while is depressed when it becomes the majority option. This implies
that the payoff matrix of the game becomes dependent on the global awareness of
the agents measured by the relevance of the population of the indifferent
players. The resulting dynamics is non-trivial with different kinds of phase
transition depending on a few model parameters. The phase diagram is studied on
regular as well as complex networks.


  Biochemical reactions are fundamentally noisy at a molecular scale. This
limits the precision of reaction networks, but also allows fluctuation
measurements which may reveal the structure and dynamics of the underlying
biochemical network. Here, we study non-equilibrium reaction cycles, such as
the mechanochemical cycle of molecular motors, the phosphorylation cycle of
circadian clock proteins, or the transition state cycle of enzymes.
Fluctuations in such cycles may be measured using either of two classical
definitions of the randomness parameter, which we show to be equivalent in
general microscopically reversible cycles. We define a stochastic period for
reversible cycles and present analytical solutions for its moments.
Furthermore, we associate the two forms of the randomness parameter with the
thermodynamic uncertainty relation, which sets limits on the timing precision
of the cycle in terms of thermodynamic quantities. Our results should prove
useful also for the study of temporal fluctuations in more general networks.


  First-passage times in random walks have a vast number of diverse
applications in physics, chemistry, biology, and finance. In general,
environmental conditions for a stochastic process are not constant on the time
scale of the average first-passage time, or control might be applied to reduce
noise. We investigate moments of the first-passage time distribution under a
transient describing relaxation of environmental conditions. We solve the
Laplace-transformed (generalized) master equation analytically using a novel
method that is applicable to general state schemes. The first-passage time from
one end to the other of a linear chain of states is our application for the
solutions. The dependence of its average on the relaxation rate obeys a power
law for slow transients. The exponent $\nu$ depends on the chain length $N$
like $\nu=-N/(N+1)$ to leading order. Slow transients substantially reduce the
noise of first-passage times expressed as the coefficient of variation (CV),
even if the average first-passage time is much longer than the transient. The
CV has a pronounced minimum for some lengths, which we call resonant lengths.
These results also suggest a simple and efficient noise control strategy, and
are closely related to the timing of repetitive excitations, coherence
resonance and information transmission by noisy excitable systems. A resonant
number of steps from the inhibited state to the excitation threshold and slow
recovery from negative feedback provide optimal timing noise reduction and
information transmission.


  The metabolic fluxes in cells follow physical, biochemical, and economic
principles. Some flux balance analysis (FBA) methods trade flux benefit against
flux cost. However, if flux cost functions are linear and meant to describe
underlying enzyme costs, this entails that enzyme efficiencies are constant and
ignores the interplay between fluxes, metabolite concentrations and enzyme
levels in cells. Here I introduce realistic flux cost functions that describe
an "overhead cost", namely the minimum enzyme and metabolite cost associated
with the fluxes in a kinetic model. These flux cost functions have general
mathematical properties. Enzymatic flux cost functions, which represent enzyme
costs, scale proportionally with the flux profile and are concave on the flux
polytope. Kinetic flux cost functions represent the sum of enzyme and
metabolite costs. If two flux profiles are superimposed, their different
demands for metabolite concentrations cause an extra compromise cost, which
makes flux cost functions strictly concave in almost all cases. When fluxes
change their direction, the enzymatic cost jumps abruptly. Here I propose two
methods for flux modelling: Flux Cost Minimisation, a nonlinear variant of FBA
with flux minimisation, and Flux Benefit Optimisation, a nonlinear variant of
FBA with molecular crowding. The optimal flux profiles, at a given flux
benefit, are vertices of the flux polytope. Linear approximations of enzymatic
flux cost can be used in FBA. In contrast to flux costs chosen ad hoc, these
functions reflect the enzyme kinetics and extracellular concentrations in
realistic kinetic models. Based on enzymatic flux costs, we can describe the
cell growth rate as a convex function on the flux polytope and derive
growth-optimal metabolic states and statistical distributions for the fluxes in
cell populations.


  Despite their dynamic nature, certain chromatin marks must be maintained over
the long term. This is particulary true for histone 3 lysine 9 (H3K9)
trimethylation, that is involved in the maintenance of healthy differentiated
cellular states by preventing inappropriate gene expression, and has been
recently identified as the most efficient barrier to cellular reprogramming in
nuclear transfer experiments. We propose that the capacity of the enzymes
SUV39H1/2 to rebind to a minor fraction of their products, either directly or
via HP1, contributes to the solidity of this mark through (i) a positive
feedback involved in its establishment by the mutual enforcement of H3K9me3 and
SUV39H1/2 and then (ii) a negative feedback sufficient to strongly stabilize
H3K9me3 heterochromatin in post-mitotic cells by generating local enzyme
concentrations capable of counteracting transient bursts of demethylation. This
model does not require direct molecular interactions with adjacent nucleosomes
and is favoured by a series of additional mechanisms including (i) the
protection of chromatin-bound SUV39H1/2 from the turnovers of soluble proteins,
which can explain the uncoupling between the cellular contents in SUV39H1 mRNA
and protein; (ii) the cooperative dependence on the local density of the
H3K9me3 of HP1-dependent heterochomatin condensation and, dispensably (iii)
restricted enzyme exchanges with chromocenters confining the reactive bursts of
SUV39H1/2 in heterochromatin. This mechanism illustrates how seemingly static
epigenetic states can be firmly maintained by dynamic and reversible
modifications.


  Explaining the origin of life requires us to explain how self-replication
arises. To be specific, how can a self-replicating entity develop spontaneously
from a chemical reaction system in which no reaction is self-replicating?
Previously proposed mathematical models either supply an explicit framework for
a minimal living system or only consider catalyzed reactions, and thus fail to
provide a comprehensive theory. We set up a general model for chemical reaction
systems that properly accounts for energetics, kinetics and the conservation
law. We find that (1) some systems are collectively-catalytic where reactants
are transformed into end products with the assistance of intermediates (as in
the citric acid cycle), while some others are self-replicating where different
parts replicate each other and the system self-replicates as a whole (as in the
formose reaction); (2) many alternative chemical universes often contain one or
more such systems; (3) it is possible to construct a self-replicating system
where the entropy of some parts spontaneously decreases, in a manner similar to
that discussed by Schr\"odinger; (4) complex self-replicating molecules can
emerge spontaneously and relatively easily from simple chemical reaction
systems through a sequence of transitions. Together these results start to
explain the origins of prebiotic evolution.


  Commonly in biomedical research, studies collect data in which an outcome
measure contains informative excess zeros; for example when observing the
burden of neuritic plaques in brain pathology studies, those who show none
contribute to our understanding of neurodegenerative disease. The outcome may
be characterized by a mixture distribution with one component being the
`structural zero' and the other component being a Poisson distribution. We
propose a novel variance components score test of genetic association between a
set of genetic markers and a zero-inflated count outcome from a mixture
distribution. This test shares advantageous properties with SNP-set tests which
have been previously devised for standard continuous or binary outcomes, such
as the Sequence Kernel Association Test (SKAT). In particular, our method has
superior statistical power compared to competing methods, especially when there
is correlation within the group of markers, and when the SNPs are associated
with both the mixing proportion and the rate of the Poisson distribution. We
apply the method to Alzheimer's data from the Rush University Religious Orders
Study and Memory and Aging Project (ROSMAP), where as proof of principle we
find highly significant associations with the APOE gene, in both the
`structural zero' and `count' parameters, when applied to a zero-inflated
neuritic plaques count outcome.


  DNA rearrangement processes recombine gene segments that are organized on the
chromosome in a variety of ways. The segments can overlap, interleave or one
may be a subsegment of another. We use directed graphs to represent segment
organizations on a given locus where contigs containing rearranged segments
represent vertices and the edges correspond to the segment relationships. Using
graph properties we associate a point in a higher dimensional Euclidean space
to each graph such that cluster formations and analysis can be performed with
methods from topological data analysis. The method is applied to a recently
sequenced model organism \textit{Oxytricha trifallax}, a species of ciliate
with highly scrambled genome that undergoes massive rearrangement process after
conjugation. The analysis shows some emerging star-like graph structures
indicating that segments of a single gene can interleave, or even contain all
of the segments from fifteen or more other genes in between its segments. We
also observe that as many as six genes can have their segments mutually
interleaving or overlapping.


  The dynamics of a two-species community of $N$ competing individuals is
considered, with an emphasis on the role of environmental variations that
affect coherently the fitness of entire populations. The chance of fixation of
a mutant population (and, in particular, of a single mutant) is calculated as a
function of the mean relative fitness of the mutant, the amplitude of fitness
variations and their typical duration. We emphasize the distinction between the
case of local competition and the case of global competition; in the latter a
noise-induced stabilization mechanism yields a higher chance of fixation for a
single mutant. The distinction becomes dramatic in the weak selection phase,
where the chance of fixation for a single deleterious mutant is an
$N$-independent constant for global competition and decays like $(\ln N)^{-1}$
in the local competition case. The possibility of long-term persistence of
${\cal O} (N)$ suboptimal (and extinction-prone) populations is discussed, as
well as its relevance to stochastic tunneling between fitness peaks.


  Precise positioning of the cell division site is essential for the correct
segregation of the genetic material into the two daughter cells. In the
bacterium Myxococcus xanthus, the proteins PomX and PomY form a cluster on the
chromosome that performs a biased random walk to midcell and positively
regulates cell division there. PomZ, an ATPase, is necessary for tethering of
the cluster to the nucleoid and regulates its movement towards midcell. It has
remained unclear how the cluster dynamics change when the biochemical
parameters, such as the attachment rates of PomZ to the nucleoid and the
cluster, the ATP hydrolysis rate of PomZ or the mobility of PomZ dimers
interacting with the nucleoid and cluster, are varied. To answer these
questions, we investigate a one-dimensional model that includes the nucleoid,
the Pom cluster and the PomZ protein. We find that a mechanism based on the
diffusive PomZ fluxes on the nucleoid into the cluster can explain the latter's
midnucleoid localization for a broad parameter range. Furthermore, there is an
ATP hydrolysis rate that minimizes the time the cluster needs to reach
midnucleoid. If the dynamics of PomZ dimers on the nucleoid is slow relative to
the cluster's velocity, we observe oscillatory cluster movements around
midnucleoid. To understand midnucleoid localization, we developed a
semi-analytical approach that dissects the net movement of the cluster into its
components: the difference in PomZ fluxes into the cluster from either side,
the force exerted by a single PomZ dimer on the cluster and the effective
friction coefficient of the cluster. Importantly, we predict that the Pom
cluster oscillates around midnucleoid if the diffusivity of PomZ on the
nucleoid is reduced. A similar approach to that applied here may also prove
useful for cargo localization in ParABS systems.


  The computational prediction of a protein structure from its sequence
generally relies on a method to assess the quality of protein models. Most
assessment methods rank candidate models using heavily engineered structural
features, defined as complex functions of the atomic coordinates. However, very
few methods have attempted to learn these features directly from the data. We
show that deep convolutional networks can be used to predict the ranking of
model structures solely on the basis of their raw three-dimensional atomic
densities, without any feature tuning. We develop a deep neural network that
performs on par with state-of-the-art algorithms from the literature. The
network is trained on decoys from the CASP7 to CASP10 datasets and its
performance is tested on the CASP11 dataset. On the CASP11 stage 2 dataset, it
achieves a loss of 0.064, whereas the best performing method achieves a loss of
0.063. Additional testing on decoys from the CASP12, CAMEO, and 3DRobot
datasets confirms that the network performs consistently well across a variety
of protein structures. While the network learns to assess structural decoys
globally and does not rely on any predefined features, it can be analyzed to
show that it implicitly identifies regions that deviate from the native
structure.


  Zooplankton have evolved several mechanisms to deal with environmental
threats, such as ultraviolet radiation (UVR), and in order to identify
strategies inherent to organisms exposed to different UVR environments, we here
examine life-history traits of two lineages of Daphnia pulex. The lineages
differed in the UVR dose they had received at their place of origin from
extremely high UVR stress at high-altitude Bolivian lakes to low UVR stress
near the sea level in temperate Sweden. Nine life-history variables of each
lineage were analyzed in laboratory experiments in the presence and the absence
of sub-lethal doses of UVR (UV-A band), and we identified trade-offs among
variables through structural equation modeling (SEM). The UVR treatment was
detrimental to almost all life-history variables of both lineages; however, the
Daphnia historically exposed to higher doses of UVR (HighUV) showed a higher
overall fecundity than those historically exposed to lower doses of UVR
(LowUV). The total offspring and ephippia production, as well as the number of
clutches and number of offspring at first reproduction, was directly affected
by UVR in both lineages. Main differences between lineages involved indirect
effects that affected offspring production as the age at first reproduction. We
here show that organisms within the same species have developed different
strategies as responses to UVR, although no increased physiological tolerance
or plasticity was shown by the HighUV lineage. In addition to known tolerance
strategies to UVR, including avoidance, prevention, or repairing of damages, we
here propose a population strategy that includes early reproduction and high
fertility, which we show compensated for the fitness loss imposed by UVR
stress.


  We consider a stochastic model of population dynamics where each individual
is characterised by a trait in {0,1,...,L} and has a natural reproduction rate,
a logistic death rate due to age or competition and a probability of mutation
towards neighbouring traits at each reproduction event. We choose parameters
such that the induced fitness landscape exhibits a valley: mutant individuals
with negative fitness have to be created in order for the population to reach a
trait with positive fitness. We focus on the limit of large population and rare
mutations at several speeds. In particular, when the mutation rate is low
enough, metastability occurs: the exit time of the valley is random,
exponentially distributed.


  Astrocytes express a large variety of G~protein-coupled receptors (GPCRs)
which mediate the transduction of extracellular signals into intracellular
calcium responses. This transduction is provided by a complex network of
biochemical reactions which mobilizes a wealth of possible calcium-mobilizing
second messenger molecules. Inositol 1,4,5-trisphosphate is probably the best
known of these molecules whose enzymes for its production and degradation are
nonetheless calcium-dependent. We present a biophysical modeling approach based
on the assumption of Michaelis-Menten enzyme kinetics, to effectively describe
GPCR-mediated astrocytic calcium signals. Our model is then used to study
different mechanisms at play in stimulus encoding by shape and frequency of
calcium oscillations in astrocytes.