  Ultrasound is increasingly being used to modulate the properties of
biological membranes for applications in drug delivery and neuromodulation.
While various studies have investigated the mechanical aspect of the
interaction such as acoustic absorption and membrane deformation, it is not
clear how these effects transduce into biological functions, for example,
changes in the permeability or the enzymatic activity of the membrane. A
critical aspect of the activity of an enzyme is the thermal fluctuations of its
solvation or hydration shell. Thermal fluctuations are also known to be
directly related to membrane permeability. Here solvation shell changes of
lipid membranes subject to an acoustic impulse were investigated using a
fluorescence probe, Laurdan. Laurdan was embedded in multi-lamellar lipid
vesicles in water, which were exposed to broadband pressure impulses of the
order of 1MPa peak amplitude and 10{\mu}s pulse duration. An instrument was
developed to monitor changes in the emission spectrum of the dye at two
wavelengths with sub-microsecond temporal resolution. The experiments show that
changes in the emission spectrum, and hence the fluctuations of the solvation
shell, are related to the changes in the thermodynamic state of the membrane
and correlated with the compression and rarefaction of the incident sound wave.
The results suggest that acoustic fields affect the state of a lipid membrane
and therefore can potentially modulate the kinetics of channels and proteins
embedded in the membrane.


  The consensus that complexity begets stability in ecosystems was challenged
in the seventies, a result recently extended to ecologically-inspired networks.
The approaches assume the existence of a feasible equilibrium, i.e. with
positive abundances. However, this key assumption has not been tested. We
provide analytical results complemented by simulations which show that
equilibrium feasibility vanishes in species rich systems. This result leaves us
in the uncomfortable situation in which the existence of a feasible equilibrium
assumed in local stability criteria is far from granted. We extend our analyses
by changing interaction structure and intensity, and find that feasibility and
stability is warranted irrespective of species richness with weak interactions.
Interestingly, we find that the dynamical behaviour of ecologically inspired
architectures is very different and richer than that of unstructured systems.
Our results suggest that a general understanding of ecosystem dynamics requires
focusing on the interplay between interaction strength and network
architecture.


  The paper deals with using chaos to direct trajectories to targets and
analyzes ruggedness and fractality of the resulting fitness landscapes. The
targeting problem is formulated as a dynamic fitness landscape and four
different chaotic maps generating such a landscape are studied. By using a
computational approach, we analyze properties of the landscapes and quantify
their fractal and rugged characteristics. In particular, it is shown that
ruggedness measures such as correlation length and information content are
scale-invariant and self-similar.


  The properties and fundamental limits of chemical computers have recently
attracted significant interest as a model of computation, an unifying principle
of cellular organisation and in the context of bio-engineering. As of yet,
research in this topic is based on case-studies. There exists no generally
accepted criterion to distinguish between chemical processes that compute and
those that do not. Here, the concept of entropy driven computer (EDC) is
proposed as a general model of chemical computation. It is found that entropy
driven computation is subject to a trade-off between accuracy and entropy
production, but unlike many biological systems, there are no trade-offs
involving time. The latter only arise when it is taken into account that the
observation of the state of the EDC is not energy neutral, but comes at a cost.
The significance of this conclusion in relation to biological systems is
discussed. Three examples of biological computers, including an implementation
of a neural network as an EDC are given.


  We propose a model for the formation of chromatin loops based on the
diffusive sliding of a DNA-bound factor which can dimerise to form a molecular
slip-link. Our slip-links mimic the behaviour of cohesin-like molecules, which,
along with the CTCF protein, stabilize loops which organize the genome. By
combining 3D Brownian dynamics simulations and 1D exactly solvable
non-equilibrium models, we show that diffusive sliding is sufficient to account
for the strong bias in favour of convergent CTCF-mediated chromosome loops
observed experimentally. Importantly, our model does not require any
underlying, and energetically costly, motor activity of cohesin. We also find
that the diffusive motion of multiple slip-links along chromatin may be
rectified by an intriguing ratchet effect that arises if slip-links bind to the
chromatin at a preferred "loading site". This emergent collective behaviour is
driven by a 1D osmotic pressure which is set up near the loading point, and
favours the extrusion of loops which are much larger than the ones formed by
single slip-links.


  We investigate methods to vaccinate contact networks -- i.e. removing nodes
in such a way that disease spreading is hindered as much as possible -- with
respect to their cost-efficiency. Any real implementation of such protocols
would come with costs related both to the vaccination itself, and gathering of
information about the network. Disregarding this, we argue, would lead to
erroneous evaluation of vaccination protocols. We use the
susceptible-infected-recovered model -- the generic model for diseases making
patients immune upon recovery -- as our disease-spreading scenario, and analyze
outbreaks on both empirical and model networks. For different relative costs,
different protocols dominate. For high vaccination costs and low costs of
gathering information, the so-called acquaintance vaccination is the most cost
efficient. For other parameter values, protocols designed for query-efficient
identification of the network's largest degrees are most efficient.


  In this paper, we develop and analyze a minimal model for a 2D network of
cross-linked actin filaments and myosin motors, representing the cortical
cytoskeleton of eukaryotic cells. We implement coarse-grained representations
of force production by myosin motors and stress dissipation through an
effective cross-link friction and filament turnover. We use this model to
characterize how the sustained production of active stress, and the steady
dissipation of elastic stress, depend individually on motor activity, effective
cross-link friction and filament turnover. Then we combine these results to
gain insights into how microscopic network parameters control steady state flow
produced by asymmetric distributions of motor activity. Our results provide a
framework for understanding how local modulation of microscopic interactions
within contractile networks control macroscopic quantities like active stress
and effective viscosity to control cortical deformation and flow at cellular
scales.


  Combinatorial analysis of a certain abstract of RNA structures has been
studied to investigate their statistics. Our approach regards the backbone of
secondary structures as an alternate sequence of paired and unpaired sets of
nucleotides, which can be described by random matrix model. We obtain the
generating function of the structures using Hermitian matrix model with
Chebyshev polynomial of the second kind and analyze the statistics with respect
to the number of stems. To match the experimental findings of the statistical
behavior, we consider the structures in a grand canonical ensemble and find a
fugacity value corresponding to an appropriate number of stems.


  Nowadays different experimental techniques, such as single molecule or
relaxation experiments, can provide dynamic properties of biomolecular systems,
but the amount of detail obtainable with these methods is often limited in
terms of time or spatial resolution. Here we use state-of-the-art computational
techniques, namely atomistic molecular dynamics and Markov state models, to
provide insight into the rapid dynamics of short RNA oligonucleotides, in order
to elucidate the kinetics of stacking interactions. Analysis of multiple
microsecond-long simulations indicates that the main relaxation modes of such
molecules can consist of transitions between alternative folded states, rather
than between random coils and native structures. After properly removing
structures that are artificially stabilized by known inaccuracies of the
current RNA AMBER force field, the kinetic properties predicted are consistent
with the timescales of previously reported relaxation experiments.


  The view that the initiation of branching into two sympatric species may not
require natural selection emerged in Victorian times (Fleeming Jenkin, George
Romanes, William Bateson). In the 1980s paleontologist Steven Jay Gould gave a
theoretical underpinning of this non-genic "chromosomal" view, thus reinstating
Richard Goldschmidt's "heresy" of the 1930s. From modelling studies with
computer-generated "biomorphs," zoologist Richard Dawkins also affirmed
Goldschmidt, proclaiming the "evolution of evolvability." However, in the
1990s, while Gould and Dawkins were recanting, bioinformatic, biochemical and
cytological studies were providing a deeper underpinning. In 2001 this came
under attack from leaders in the field who favored Dawkins' genic emphasis.
Now, with growing evidence from multiple sources, we can reinstate again
Goldschmidt's view and clarify its nineteenth century roots.


  Mutations in a microbial population can increase the frequency of a genotype
not only by increasing its exponential growth rate, but also by decreasing its
lag time or adjusting the yield (resource efficiency). The contribution of
multiple life-history traits to selection is a critical question for
evolutionary biology as we seek to predict the evolutionary fates of mutations.
Here we use a model of microbial growth to show there are two distinct
components of selection corresponding to the growth and lag phases, while the
yield modulates their relative importance. The model predicts rich population
dynamics when there are tradeoffs between phases: multiple strains can coexist
or exhibit bistability due to frequency-dependent selection, and strains can
engage in rock-paper-scissors interactions due to non-transitive selection. We
characterize the environmental conditions and patterns of traits necessary to
realize these phenomena, which we show to be readily accessible to experiments.
Our results provide a theoretical framework for analyzing high-throughput
measurements of microbial growth traits, especially interpreting the pleiotropy
and correlations between traits across mutants. This work also highlights the
need for more comprehensive measurements of selection in simple microbial
systems, where the concept of an ordinary fitness landscape breaks down.


  Synthetic biology sets out to implement new functions in cells, and to
develop a deeper understanding of biological design principles. In 2000,
Elowitz and Leibler showed that by rational design of the reaction network, and
using existing biological components, they could create a network that exhibits
periodic gene expression, dubbed the repressilator (Elowitz and Leibler,
Nature, 2000). More recently, Stricker et al. presented another synthetic
oscillator, called the dual-feedback oscillator (Stricker et al., 2008), which
is more stable. How the stability of these oscillators is affected by the
intrinsic noise of the interactions between the components and the stochastic
expression of their genes, has been studied in considerable detail. However, as
all biological oscillators reside in growing and dividing cells, an important
question is how these oscillators are perturbed by the cell cycle. In previous
work we showed that the periodic doubling of the gene copy numbers due to DNA
replication can couple not only natural, circadian oscillators to the cell
cycle (Paijmans et al., PNAS, \textbf{113}, 4063, (2016)), but also these
synthetic oscillators. Here we expand this study. We find that the strength of
the locking between oscillators depends not only on the positions of the genes
on the chromosome, but also on the noise in the timing of gene replication:
noise tends to weaken the coupling. Yet, even in the limit of high levels of
noise in the replication times of the genes, both synthetic oscillators show
clear signatures of locking to the cell cycle. This work enhances our
understanding of the design of robust biological oscillators inside growing and
diving cells.


  We construct a model of cell reprogramming (the conversion of fully
differentiated cells to a state of pluripotency, known as induced pluripotent
stem cells, or iPSCs) which builds on key elements of cell biology viz. cell
cycles and cell lineages. Although reprogramming has been demonstrated
experimentally, much of the underlying processes governing cell fate decisions
remain unknown. This work aims to bridge this gap by modelling cell types as a
set of hierarchically related dynamical attractors representing cell cycles.
Stages of the cell cycle are characterised by the configuration of gene
expression levels, and reprogramming corresponds to triggering transitions
between such configurations. Two mechanisms were found for reprogramming in a
two level hierarchy: cycle specific perturbations and a noise-induced
switching. The former corresponds to a \emph{directed} perturbation that
induces a transition into a cycle state of a different cell type in the potency
hierarchy (mainly a stem cell) whilst the latter is a priori undirected and
could be induced, e.g., by a (stochastic) change in the cellular environment.
These reprogramming protocols were found to be effective in large regimes of
the parameter space and make specific predictions concerning reprogramming
dynamics which are broadly in line with experimental findings.


  Gene expression levels carry information about signals that have functional
significance for the organism. Using the gap gene network in the fruit fly
embryo as an example, we show how this information can be decoded, building a
dictionary that translates expression levels into a map of implied positions.
The optimal decoder makes use of graded variations in absolute expression
level, resulting in positional estimates that are precise to ~1% of the
embryo's length. We test this optimal decoder by analyzing gap gene expression
in embryos lacking some of the primary maternal inputs to the network. The
resulting maps are distorted, and these distortions predict, with no free
parameters, the positions of expression stripes for the pair-rule genes in the
mutant embryos.


  We introduce a tensor-based clustering method to extract sparse,
low-dimensional structure from high-dimensional, multi-indexed datasets. This
framework is designed to enable detection of clusters of data in the presence
of structural requirements which we encode as algebraic constraints in a linear
program. Our clustering method is general and can be tailored to a variety of
applications in science and industry. We illustrate our method on a collection
of experiments measuring the response of genetically diverse breast cancer cell
lines to an array of ligands. Each experiment consists of a cell line-ligand
combination, and contains time-course measurements of the early-signalling
kinases MAPK and AKT at two different ligand dose levels. By imposing
appropriate structural constraints and respecting the multi-indexed structure
of the data, the analysis of clusters can be optimized for biological
interpretation and therapeutic understanding. We then perform a systematic,
large-scale exploration of mechanistic models of MAPK-AKT crosstalk for each
cluster. This analysis allows us to quantify the heterogeneity of breast cancer
cell subtypes, and leads to hypotheses about the signalling mechanisms that
mediate the response of the cell lines to ligands.


  Circadian clocks must be able to entrain to time-varying signals to keep
their oscillations in phase with the day-night rhythm. On the other hand, they
must also exhibit input compensation: their period must remain about one day in
different constant environments. The post-translational oscillator of the Kai
system can be entrained by transient or oscillatory changes in the ATP
fraction, yet is insensitive to constant changes in this fraction. We study in
three different models of this system how these two seemingly conflicting
criteria are met: the Van Zon model (Van Zon et al., PNAS, 2007), the Rust
model (Phong et al., PNAS, 2013), and our new model (Paijmans et al.,
arXiv:1612.02715). We find that the new model exhibits the best trade-off
between input compensation and entrainability: on the footing of equal
phase-response curves, it exhibits the strongest input compensation. Performing
stochastic simulations at the level of individual hexamers allows us to
identify a new mechanism, which is employed by the new model to achieve input
compensation: At lower ATP fraction, the individual hexamers make a shorter
cycle in the phosphorylation state space, which compensates for the slower pace
at which they traverse the cycle.


  In this paper, we introduce multiscale persistent functions for biomolecular
structure characterization. The essential idea is to combine our multiscale
rigidity functions with persistent homology analysis, so as to construct a
series of multiscale persistent functions, particularly multiscale persistent
entropies, for structure characterization. To clarify the fundamental idea of
our method, the multiscale persistent entropy model is discussed in great
detail. Mathematically, unlike the previous persistent entropy or topological
entropy, a special resolution parameter is incorporated into our model. Various
scales can be achieved by tuning its value. Physically, our multiscale
persistent entropy can be used in conformation entropy evaluation. More
specifically, it is found that our method incorporates in it a natural
classification scheme. This is achieved through a density filtration of a
multiscale rigidity function built from bond and/or dihedral angle
distributions. To further validate our model, a systematical comparison with
the traditional entropy evaluation model is done. It is found that our model is
able to preserve the intrinsic topological features of biomolecular data much
better than traditional approaches, particularly for resolutions in the mediate
range. Moreover, our method can be successfully used in protein classification.
For a test database with around nine hundred proteins, a clear separation
between all-alpha and all-beta proteins can be achieved, using only the
dihedral and pseudo-bond angle information. Finally, a special protein
structure index (PSI) is proposed, for the first time, to describe the
"regularity" of protein structures. Essentially, PSI can be used to describe
the "regularity" information in any systems.


  Measures of wealth and production have been found to scale superlinearly with
the population of a city. Therefore, it makes economic sense for humans to
congregate together in dense settlements. A recent model of population dynamics
showed that population growth can become superexponential due to the
superlinear scaling of production with population in a city. Here, we
generalize this population dynamics model and demonstrate the existence of
multiple stable equilibrium points, showing how population growth can be
stymied by a poor economic environment. This occurs when the goods and services
produced by the city become less profitable due to a lack of diversification in
the city's economy. Then, relying on critical slowing down signals related to
the stability of an equilibrium point, we present an algorithm for engineering
regime shifts such that a city at a stable equilibrium point may continue to
grow again. The generality of the model and the algorithm used here implies
that the model and algorithm need not be restricted to urban systems; they are
easily applicable to other types of systems where the assumptions used are
valid.


  We present a computational model to reconstruct trees of ancestors for
animals with sexual reproduction. Through a recursive algorithm combined with a
random number generator, it is possible to reproduce the number of ancestors
for each generation and use it to constraint the maximum number of the
following generation. This new model allows to consider the reproductive
preferences of particular species and combine several trees to simulate the
behavior of a population. It is also possible to obtain a description
analytically, considering the simulation as a theoretical stochastic process.
Such process can be generalized in order to use an algorithm associated with it
to simulate other similar processes of stochastic nature. The simulation is
based in the theoretical model previously presented before.


  HIV-1 virions assemble as immature particles containing Gag polyproteins that
are processed by the viral protease into individual components, resulting in
the formation of mature infectious particles. There are two competing models
for the process of forming the mature HIV-1 core: the disassembly and de novo
reassembly model and the non-diffusional displacive model. To study the
maturation pathway, we simulate HIV-1 maturation in vitro by digesting immature
particles and assembled virus-like particles with recombinant HIV-1 protease
and monitor the process with biochemical assays and cryoEM structural analysis
in parallel. Processing of Gag in vitro is accurate and efficient and results
in both soluble capsid protein and conical or tubular capsid assemblies,
seemingly converted from immature Gag particles. Computer simulations further
reveal probable assembly pathways of HIV-1 capsid formation. Combining the
experimental data and computer simulations, our results suggest a sequential
combination of both displacive and disassembly/reassembly processes for HIV-1
maturation.


  Efficient replication and assembly of virus particles are integral to the
establishment of infection. In addition to the primary role of the capsid
protein (CP) in encapsidating the RNA progeny, experimental evidence on
positive sense single-stranded RNA viruses suggests that the CP also regulates
RNA synthesis. Here, we demonstrate that replication of Satellite tobacco
mosaic virus (STMV) is controlled by the cooperative interaction between STMV
CP and the helper virus (HV) Tobacco mosaic virus (TMV) replicase. We
identified that the STMV CP-HV replicase interaction requires a positively
charged residue at the third position (3R) in the N-terminal 13 amino acid (aa)
motif. Far-Northwestern blotting showed that STMV CP promotes binding between
HV-replicase and STMV RNA. An STMV CP variant having an arginine to alanine
substitution at position 3 in the N-terminal 13aa motif abolished replicase-CP
binding. The N-terminal 13aa motif of the CP bearing alanine substitutions for
positively charged residues located at positions 5, 7, 10 and 11 are defective
in packaging full-length STMV, but can package a truncated STMV RNA lacking the
3' terminal 150 nt region. These findings provide insights into the mechanism
underlying the regulation of STMV replication and packaging.


  Interaction with divalent cations is of paramount importance for RNA
structural stability and function. We here report a detailed molecular dynamics
study of all the possible binding sites for Mg$^{2+}$ on a RNA duplex,
including both direct (inner sphere) and indirect (outer sphere) binding. In
order to tackle sampling issues, we develop a modified version of bias-exchange
metadynamics which allows us to simultaneously compute affinities with
previously unreported statistical accuracy. Results correctly reproduce trends
observed in crystallographic databases. Based on this, we simulate a carefully
chosen set of models that allows us to quantify the effects of competition with
monovalent cations, RNA flexibility, and RNA hybridization. Our simulations
reproduce the decrease and increase of Mg$^{2+}$ affinity due to ion
competition and hybridization respectively, and predict that RNA flexibility
has a site dependent effect. This suggests a non trivial interplay between RNA
conformational entropy and divalent cation binding.


  We introduce a general framework for biological systems, called MESSI
systems, that describe Modifications of type Enzyme-Substrate or Swap with
Intermediates, and we prove general results based on the network structure.
Many post-translational modification networks are MESSI systems. For example:
the motifs in [Feliu and Wiuf (2012a)], sequential distributive and processive
multisite phosphorylation networks, most of the examples in [Angeli et al.
(2007)], phosphorylation cascades, two component systems as in [Kothamachu et
al. (2015)], the bacterial EnvZ/OmpR network in [Shinar and Feinberg (2010)],
and all linear networks. We show that, under mass-action kinetics, MESSI
systems are conservative. We simplify the study of steady states of these
systems by explicit elimination of intermediate complexes and we give
conditions to ensure an explicit rational parametrization of the variety of
steady states (inspired by [Feliu and Wiuf (2013a, 2013b), Thomson and
Gunawardena (2009)]). We define an important subclass of MESSI systems with
toric steady states [P\'erez Mill\'an et al. (2012)] and we give for MESSI
systems with toric steady states an easy algorithm to determine the capacity
for multistationarity. In this case, the algorithm provides rate constants for
which multistationarity takes place, based on the theory of oriented matroids.


  Fisher's geometric model was originally introduced to argue that complex
adaptations must occur in small steps because of pleiotropic constraints. When
supplemented with the assumption of additivity of mutational effects on
phenotypic traits, it provides a simple mechanism for the emergence of
genotypic epistasis from the nonlinear mapping of phenotypes to fitness. Of
particular interest is the occurrence of reciprocal sign epistasis, which is a
necessary condition for multipeaked genotypic fitness landscapes. Here we
compute the probability that a pair of randomly chosen mutations interacts sign
epistatically, which is found to decrease with increasing phenotypic dimension
$n$, and varies nonmonotonically with the distance from the phenotypic optimum.
We then derive expressions for the mean number of fitness maxima in genotypic
landscapes comprised of all combinations of $L$ random mutations. This number
increases exponentially with $L$, and the corresponding growth rate is used as
a measure of the complexity of the landscape. The dependence of the complexity
on the model parameters is found to be surprisingly rich, and three distinct
phases characterized by different landscape structures are identified. Our
analysis shows that the phenotypic dimension, which is often referred to as
phenotypic complexity, does not generally correlate with the complexity of
fitness landscapes and that even organisms with a single phenotypic trait can
have complex landscapes. Our results further inform the interpretation of
experiments where the parameters of Fisher's model have been inferred from
data, and help to elucidate which features of empirical fitness landscapes can
be described by this model.


  Extensive molecular dynamics simulations reveal that the interactions between
proteins and poly(ethylene glycol)(PEG) can be described in terms of the
surface composition of the proteins. PEG molecules accumulate around non-polar
residues while avoiding polar ones. A solvent-accessible-surface-area model of
protein adsorption on PEGylated nanoparticles accurately fits a large set of
data on the composition of the protein corona recently obtained by label-free
proteomic mass spectrometry.


  Game theory research on the snowdrift game has showed that gradual evolution
of the continuously varying level of cooperation in joint enterprises can
demonstrate evolutionary merging as well as evolutionary branching. However,
little is known about the consequences of changes in diversity at the
cooperation level. In the present study I consider effects of costly rewards on
the continuous snowdrift game. I show that not evolutionary merging but
evolutionary branching can promote the emergence of pool reward, which can then
enhance social welfare.


  Complex diseases can be modeled as damage to intracellular networks that
results in abnormal cell behaviors. Network-based dynamic models such as
Boolean models have been employed to model a variety of biological systems
including those corresponding to disease. Previous work designed compensatory
interactions to stabilize an attractor of a Boolean network after single node
damage. We generalize this method to a multinode damage scenario and to the
simultaneous stabilization of multiple steady state attractors. We classify the
emergent situations, with a special focus on combinatorial effects, and
characterize each class through simulation. We explore how the structural and
functional properties of the network affect its resilience and its possible
repair scenarios. We demonstrate the method's applicability to two
intracellular network models relevant to cancer. This work has implications in
designing prevention strategies for complex disease.


  The search for high-affinity aptamers for targets such as proteins, small
molecules, or cancer cells remains a formidable endeavor. Systematic Evolution
of Ligands by EXponential Enrichment (SELEX) offers an iterative process to
discover these aptamers through evolutionary selection of high-affinity
candidates from a highly diverse random pool. This randomness dictates an
unknown population distribution of fitness parameters, encoded by the binding
affinities, toward SELEX targets. Adding to this uncertainty, repeating SELEX
under identical conditions may lead to variable outcomes. These uncertainties
pose a challenge when tuning selection pressures to isolate high-affinity
ligands. Here, we present a novel stochastic hybrid model that describes the
evolutionary selection of aptamers in order to explore the impact of these
unknowns. To our surprise, we find that even single copies of high-affinity
ligands in a pool of billions can strongly influence population dynamics, yet
their survival is highly dependent on chance. We perform Monte Carlo
simulations to explore the impact of environmental parameters, such as the
target concentration, on selection efficiency in SELEX and identify new
strategies to control these uncertainties to ultimately improve the outcome and
speed of this time- and resource-intensive process.


  We show, within the context of the standard class of deterministic ODE
predator-prey mathematical models, that predator culling does not produce a
long term decrease in the predator population.


  In nature microbial populations are subject to fluctuating nutrient levels.
Nutrient fluctuations are important for evolutionary and ecological dynamics in
microbial communities since they impact growth rates, population sizes and
biofilm formation. Here we use automated continuous-culture devices and
high-throughput imaging to show that when populations of Escherichia coli are
subjected to cycles of nutrient excess (feasts) and scarcity (famine) their
abundance dynamics during famines depend on the frequency and amplitude of
feasts. We show that frequency and amplitude dependent dynamics in planktonic
populations arise from nutrient and history dependent rates of aggregation and
dispersal. A phenomenological model recapitulates our experimental
observations. Our results show that the statistical properties of environmental
fluctuations have substantial impacts on spatial structure in bacterial
populations driving large changes in abundance dynamics.


  Ecosystems are complex systems, currently experiencing several threats
associated with global warming, intensive exploitation, and human-driven
habitat degradation. Such threats are pushing ecosystems to the brink of
collapse. Because of a general presence of multiple stable states, including
states involving population extinction, and due to intrinsic nonlinearities
associated with feedback loops, collapse can occur in a catastrophic manner.
Such catastrophic shifts have been suggested to pervade many of the future
transitions affecting ecosystems at many different scales. Many studies have
tried to delineate potential warning signals predicting such ongoing shifts but
little is known about how such transitions might be effectively prevented. It
has been recently suggested that a potential path to prevent or modify the
outcome of these transitions would involve designing synthetic organisms and
synthetic ecological interactions that could push these endangered systems out
of the critical boundaries. Four classes of such ecological engineering designs
or {\em Terraformation motifs} have been defined in a qualitative way. Here we
develop the simplest mathematical models associated with these motifs, defining
the expected stability conditions and domains where the motifs shall properly
work.


  The probability distribution of sequences with maximum entropy that satisfies
a given amino acid composition at each site and a given pairwise amino acid
frequency at each site pair is a Boltzmann distribution with $\exp(-\psi_N)$,
where the total interaction $\psi_N$ is represented as the sum of one body and
pairwise interactions. A protein folding theory based on the random energy
model (REM) indicates that the equilibrium ensemble of natural protein
sequences is a canonical ensemble characterized by $\exp(-\Delta G_{ND}/k_B
T_s)$ or by $\exp(- G_{N}/k_B T_s)$ if an amino acid composition is kept
constant, meaning $\psi_N = \Delta G_{ND}/k_B T_s +$ constant, where $\Delta
G_{ND} \equiv G_N - G_D$, $G_N$ and $G_D$ are the native and denatured free
energies, and $T_s$ is the effective temperature of natural selection. Here, we
examine interaction changes ($\Delta \psi_N$) due to single nucleotide
nonsynonymous mutations, and have found that the variance of their $\Delta
\psi_N$ over all sites hardly depends on the $\psi_N$ of each homologous
sequence, indicating that the variance of $\Delta G_N (= k_B T_s \Delta
\psi_N)$ is nearly constant irrespective of protein families. As a result,
$T_s$ is estimated from the ratio of the variance of $\Delta \psi_N$ to that of
a reference protein, which is determined by a direct comparison between
$\Delta\Delta \psi_{ND} (\simeq \Delta \psi_N)$ and experimental $\Delta\Delta
G_{ND}$. Based on the REM, glass transition temperature $T_g$ and $\Delta
G_{ND}$ are estimated from $T_s$ and experimental melting temperatures ($T_m$)
for 14 protein domains. The estimates of $\Delta G_{ND}$ agree well with their
experimental values for 5 proteins, and those of $T_s$ and $T_g$ are all within
a reasonable range. This method is coarse-grained but much simpler in
estimating $T_s$, $T_g$ and $\Delta\Delta G_{ND}$ than previous methods.


  Assuming that mutation and fixation processes are reversible Markov
processes, we prove that the equilibrium ensemble of sequences obeys a
Boltzmann distribution with $\exp(4N_e m(1 - 1/(2N)))$, where $m$ is Malthusian
fitness and $N_e$ and $N$ are effective and actual population sizes. On the
other hand, the probability distribution of sequences with maximum entropy that
satisfies a given amino acid composition at each site and a given pairwise
amino acid frequency at each site pair is a Boltzmann distribution with
$\exp(-\psi_N)$, where $\psi_N$ is represented as the sum of one body and
pairwise potentials. A protein folding theory indicates that homologous
sequences obey a canonical ensemble characterized by $\exp(-\Delta G_{ND}/k_B
T_s)$ or by $\exp(- G_{N}/k_B T_s)$ if an amino acid composition is kept
constant, where $\Delta G_{ND} \equiv G_N - G_D$, $G_N$ and $G_D$ are the
native and denatured free energies, and $T_s$ is selective temperature. Thus,
$4N_e m (1 - 1/(2N))$, $-\Delta \psi_{ND}$, and $-\Delta G_{ND}/k_B T_s$ must
be equivalent to each other. Based on the analysis of the changes ($\Delta
\psi_N$) of $\psi_N$ due to single nucleotide nonsynonymous substitutions,
$T_s$, and then glass transition temperature $T_g$, and $\Delta G_{ND}$ are
estimated with reasonable values for 14 protein domains. In addition,
approximating the probability density function (PDF) of $\Delta \psi_N$ by a
log-normal distribution, PDFs of $\Delta \psi_N$ and $K_a/K_s$, which is the
ratio of nonsynonymous to synonymous substitution rate per site, in all and in
fixed mutants are estimated. It is confirmed that $T_s$ negatively correlates
with the mean of $K_a/K_s$. Stabilizing mutations are significantly fixed by
positive selection, and balance with destabilizing mutations fixed by random
drift. Supporting the nearly neutral theory, neutral selection is not
significant.


  Recently described stochastic models of protein evolution have demonstrated
that the inclusion of structural information in addition to amino acid
sequences leads to a more reliable estimation of evolutionary parameters. We
present a generative, evolutionary model of protein structure and sequence that
is valid on a local length scale. The model concerns the local dependencies
between sequence and structure evolution in a pair of homologous proteins. The
evolutionary trajectory between the two structures in the protein pair is
treated as a random walk in dihedral angle space, which is modelled using a
novel angular diffusion process on the two-dimensional torus. Coupling sequence
and structure evolution in our model allows for modelling both "smooth"
conformational changes and "catastrophic" conformational jumps, conditioned on
the amino acid changes. The model has interpretable parameters and is
comparatively more realistic than previous stochastic models, providing new
insights into the relationship between sequence and structure evolution. For
example, using the trained model we were able to identify an apparent
sequence-structure evolutionary motif present in a large number of homologous
protein pairs. The generative nature of our model enables us to evaluate its
validity and its ability to simulate aspects of protein evolution conditioned
on an amino acid sequence, a related amino acid sequence, a related structure
or any combination thereof.


  Cancer is known as a disease mainly caused by gene alterations. Discovery of
mutated driver pathways or gene sets is becoming an important step to
understand molecular mechanisms of carcinogenesis. However, systematically
investigating commonalities and specificities of driver gene sets among
multiple cancer types is still a great challenge, but this investigation will
undoubtedly benefit deciphering cancers and will be helpful for personalized
therapy and precision medicine in cancer treatment. In this study, we propose
two optimization models to \emph{de novo} discover common driver gene sets
among multiple cancer types (ComMDP) and specific driver gene sets of one
certain or multiple cancer types to other cancers (SpeMDP), respectively. We
first apply ComMDP and SpeMDP to simulated data to validate their efficiency.
Then, we further apply these methods to 12 cancer types from The Cancer Genome
Atlas (TCGA) and obtain several biologically meaningful driver pathways. As
examples, we construct a common cancer pathway model for BRCA and OV, infer a
complex driver pathway model for BRCA carcinogenesis based on common driver
gene sets of BRCA with eight cancer types, and investigate specific driver
pathways of the liquid cancer lymphoblastic acute myeloid leukemia (LAML)
versus other solid cancer types. In these processes more candidate cancer genes
are also found.


  Estimating the phylogeny of the genus Homo is entering a new phase of vastly
improved data and methodology. There is increasing evidence of 6 to 10
competing species/lineages at any point in the last half million years, making
the elucidation of the relationships of individual specimens particularly
important. Recent estimates of the phylogeny of key specimens include Waddell
(2013, 2014, 2015, 2016), and Mounier et al. (2016). These are made with quite
different data (3D skull shapes and discrete morphological characters,
respectively) and methods of analysis (unweighted least squares fitting of
distances, OLS+, and reweighted maximum parsimony, respectively). Initial
inspection of the trees in these articles might leave the impression of a great
deal of disagreement and confused results. Here it is shown this need not be
the case, and that these two types of data and analysis may be indicating a
very similar tree, one that is in good agreement also with subjective current
wisdom/expert opinions on particular parts of the phylogeny. The precise
location of the African LH18 specimen arises as key to a better understanding
of the likely form of the last common ancestors of H. sapiens and Neanderthals.
A diverse approach seems to bring forth much more agreement of trees than
otherwise perceived, and argues against being dogmatic about methods of
phylogenetic analysis particularly when working with difficult problems.


  We propose a bio-inspired, agent-based approach to describe the natural
phenomenon of group chasing in both two and three dimensions. Using a set of
local interaction rules we created a continuous-space and discrete-time model
with time delay, external noise and limited acceleration. We implemented a
unique collective chasing strategy, optimized its parameters and studied its
properties when chasing a much faster, erratic escaper. We show that collective
chasing strategies can significantly enhance the chasers' success rate. Our
realistic approach handles group chasing within closed, soft boundaries -
contrasting most of those published in the literature with periodic ones -- and
resembles several properties of pursuits observed in nature, such as the
emergent encircling or the escaper's zigzag motion.


  Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely
used in enzyme simulation. Over ten convergence studies of QM/MM methods have
revealed over the past several years that key energetic and structural
properties approach asymptotic limits with only very large (ca. 500-1000 atom)
QM regions. This slow convergence has been observed to be due in part to
significant charge transfer between the core active site and surrounding
protein environment, which cannot be addressed by improvement of MM force
fields or the embedding method employed within QM/MM. Given this slow
convergence, it becomes essential to identify strategies for the most
atom-economical determination of optimal QM regions and to gain insight into
the crucial interactions captured only in large QM regions. Here, we extend and
develop two methods for quantitative determination of QM regions. First, in the
charge shift analysis (CSA) method, we probe the reorganization of electron
density when core active site residues are removed completely, as determined by
large-QM region QM/MM calculations. Second, we introduce the
highly-parallelizable Fukui shift analysis (FSA), which identifies how
core/substrate frontier states are altered by the presence of an additional QM
residue on smaller initial QM regions. We demonstrate that the FSA and CSA
approaches are complementary and consistent on three test case enzymes:
catechol O-methyltransferase, cytochrome P450cam, and hen eggwhite lysozyme. We
also introduce validation strategies and test sensitivities of the two methods
to geometric structure, basis set size, and electronic structure methodology.
Both methods represent promising approaches for the systematic, unbiased
determination of quantum mechanical effects in enzymes and large systems that
necessitate multi-scale modeling.


  Motif finding in DNA, RNA and proteins plays an important role in life
science research. Recent patents concerning motif finding in the biomolecular
data are recorded in the DNA Patent Database which serves as a resource for
policy makers and members of the general public interested in fields like
genomics, genetics and biotechnology. In this paper we present a computational
approach to mining for RNA tertiary motifs in genomic sequences. Specifically
we describe a method, named CSminer, for finding RNA coaxial helical stackings
in genomes. A coaxial helical stacking occurs in an RNA tertiary structure
where two separate helical elements form a pseudocontiguous helix and provides
thermodynamic stability to the molecule as a whole. Experimental results
demonstrate the effectiveness of our approach.


  Protein motifs are conserved fragments occurred frequently in protein
sequences. They have significant functions, such as active site of an enzyme.
Search and clustering protein sequence motifs are computational intensive. Most
existing methods are not fast enough to analyze large data sets for motif
finding or achieve low accuracy for motif clustering. We present a new protein
sequence motif finding and clustering algorithm, called HSEARCH. It converts
fixed length protein sequences to data points in high dimensional space, and
applies locality-sensitive hashing to fast search homologous protein sequences
for a motif. HSEARCH is significantly faster than the brute force algorithm for
protein motif finding and achieves high accuracy for protein motif clustering.


  Biological organisms have to cope with stochastic variations in both the
external environment and the internal population dynamics. Theoretical studies
and laboratory experiments suggest that population diversification could be an
effective bet-hedging strategy for adaptation to varying environments. Here we
show that bet-hedging can also be effective against demographic fluctuations
that pose a trade-off between growth and survival for populations even in a
constant environment. A species can maximize its overall abundance in the long
term by diversifying into coexisting subpopulations of both "fast-growing" and
"better-surviving" individuals. Our model generalizes statistical physics
models of birth-death processes to incorporate dispersal, during which new
populations are founded, and can further incorporate variations of local
environments. In this way we unify different bet-hedging strategies against
demographic and environmental variations as a general means of adaptation to
both types of uncertainties in population growth.


  In order to function reliably, synthetic molecular circuits require
mechanisms that allow them to adapt to environmental disturbances. Least mean
squares (LMS) schemes, such as commonly encountered in signal processing and
control, provide a powerful means to accomplish that goal. In this paper we
show how the traditional LMS algorithm can be implemented at the molecular
level using only a few elementary biomolecular reactions. We demonstrate our
approach using several simulation studies and discuss its relevance to
synthetic biology.


  Computational prediction of origin of replication (ORI) has been of great
interest in bioinformatics and several methods including GC Skew, Z curve,
auto-correlation etc. have been explored in the past. In this paper, we have
extended the auto-correlation method to predict ORI location with much higher
resolution for prokaryotes. The proposed complex correlation method (iCorr)
converts the genome sequence into a sequence of complex numbers by mapping the
nucleotides to {+1,-1,+i,-i} instead of {+1,-1} used in the auto-correlation
method (here, 'i' is square root of -1). Thus, the iCorr method uses
information about the positions of all the four nucleotides unlike the earlier
auto-correlation method which uses the positional information of only one
nucleotide. Also, this earlier method required visual inspection of the
obtained graphs to identify the location of origin of replication. The proposed
iCorr method does away with this need and is able to identify the origin
location simply by picking the peak in the iCorr graph. The iCorr method also
works for a much smaller segment size compared to the earlier auto-correlation
method, which can be very helpful in experimental validation of the
computational predictions. We have also developed a variant of the iCorr method
to predict ORI location in eukaryotes and have tested it with the
experimentally known origin locations of S. cerevisiae with an average accuracy
of 71.76%.


  Biological functions are typically performed by groups of cells that express
predominantly the same genes, yet display a continuum of phenotypes. While it
is known how one genotype can generate such non-genetic diversity, it remains
unclear how different phenotypes contribute to the performance of biological
function at the population level. We developed a microfluidic device to
simultaneously measure the phenotype and chemotactic performance of tens of
thousands of individual, freely-swimming Escherichia coli as they climbed a
gradient of attractant. We discovered that spatial structure spontaneously
emerged from initially well-mixed wild type populations due to non-genetic
diversity. By manipulating the expression of key chemotaxis proteins, we
established a causal relationship between protein expression, non- genetic
diversity, and performance that was theoretically predicted. This approach
generated a complete phenotype-to-performance map, in which we found a
nonlinear regime. We used this map to demonstrate how changing the shape of a
phenotypic distribution can have as large of an effect on collective
performance as changing the mean phenotype, suggesting that selection could act
on both during the process of adaptation.


  How can tissues generate large numbers of cells, yet keep the divisional load
(the number of divisions along cell lineages) low in order to curtail the
accumulation of somatic mutations and reduce the risk of cancer? To answer the
question we consider a general model of hierarchically organized self-renewing
tissues and show that the lifetime divisional load of such a tissue is
independent of the details of the cell differentiation processes, and depends
only on two structural and two dynamical parameters. Our results demonstrate
that a strict analytical relationship exists between two seemingly disparate
characteristics of self-renewing tissues: divisional load and tissue
organization. Most remarkably, we find that a sufficient number of
progressively slower dividing cell types can be almost as efficient in
minimizing the divisional load, as non-renewing tissues. We argue that one of
the main functions of tissue-specific stem cells and differentiation
hierarchies is the prevention of cancer.


  To introduce selection into a model of coalescence, I explore the use of
modified integer partitions that allow the identification of a preferred
lineage. I show that a partition-partition transition matrix, along with Monte
Carlo discrete time kinetics, treats both the neutral case and a wide range of
positive and negative selection pressures for small population sizes. Selection
pressure causes multiple collisions per generation, short coalescence times,
increased lengths of terminal branches, increased tree asymmetry, and
dependence of coalescence times on the logarithm of population size. These
features are consistent with higher order coalescences that permit multiple
collisions per generation. While the treatment is exact in terms of the
simplified Wright-Fisher model used, it is not easily extended to large
population size. Keywords: Selection, Coalescence, Integer Partitions, Multiple
Collisions, Tree Asymmetry.


  Circular polarization spectroscopy has proven to be an indispensable tool in
photosynthesis research and (bio)-molecular research in general. Oxygenic
photosystems typically display an asymmetric Cotton effect around the
chlorophyll absorbance maximum with a signal $\leq 1 \%$. In vegetation, these
signals are the direct result of the chirality of the supramolecular
aggregates. The circular polarization is thus directly influenced by the
composition and architecture of the photosynthetic macrodomains, and is thereby
linked to photosynthetic functioning. Although ordinarily measured only on a
molecular level, we have developed a new spectropolarimetric instrument,
TreePol, that allows for both laboratory and in-the-field measurements. Through
spectral multiplexing, TreePol is capable of fast measurements with a
sensitivity of $\sim 1*10^{-4}$ and is therefore suitable of non-destructively
probing the molecular architecture of whole plant leaves. We have measured the
chiroptical evolution of \textit{Hedera helix} leaves for a period of 22 days.
Spectrally resolved circular polarization measurements (450-900 nm) on whole
leaves in transmission exhibit a strong decrease in the polarization signal
over time after plucking, which we accredit to the deterioration of chiral
macro-aggregates. Chlorophyll \textit{a} levels measured over the same period
by means of UV-Vis absorption and fluorescence spectroscopy showed a much
smaller decrease. With these results we are able to distinguish healthy from
deteriorating leaves. Hereby we indicate the potency of circular polarization
spectroscopy on whole and intact leaves as a nondestructive tool for structural
and plant stress assessment. Additionally, we underline the establishment of
circular polarization signals as remotely accessible means of detecting the
presence of extraterrestrial life.


  Recent technological advances and long-term data studies provide interaction
data that can be modelled through dynamic networks, i.e a sequence of different
snapshots of an evolving ecological network. Most often time is the parameter
along which these networks evolve but any other one-dimensional gradient
(temperature, altitude, depth, humidity, . . . ) could be considered.Here we
propose a statistical tool to analyse the underlying structure of these
networks and follow its evolution dynamics (either in time or any other
one-dimensional factor). It consists in extracting the main features of these
networks and summarise them into a high-level view.We analyse a dynamic animal
contact network and a seasonal food web and in both cases we show that our
approach allows for the identification of a backbone organisation as well as
interesting temporal variations at the individual level.Our method, implemented
into the R package dynsbm, can handle the largest ecological datasets and is a
versatile and promising tool for ecologists that study dynamic interactions.


  A cell has the ability to convert an environmental change into the expression
of genetic information through a chain of intracellular signal transduction
reactions. Here, we aimed to develop a method for quantifying this signal
transduction. We showed that the channel capacities of individual steps in a
given general model cascade were equivalent in an independent manner, and were
given by the entropy production rate. Signal transduction was transmitted by
fluctuation of the entropy production rate and quantified transduction was
estimated by the work done in individual steps. If the individual step
representing the modification to demodification of the signal molecules is
considered to be a Szilard engine, the maximal work done is equivalent to the
chemical potential change of the messenger that is consumed during the
modification reaction. Our method was applicable to calculate the channel
capacity of the MAPK cascade. In conclusion, our method is suitable for
quantitative analyses of signal transduction.


  Genetic and environmental perturbation experiments have been used to study
microbes in a bid to gain insight into transcriptional regulation, adaptive
evolution, and other cellular dynamics. These studies have potential in
enabling rational strain design. Unfortunately, experimentally determined
intracellular flux distribution are often inconsistent or incomparable due to
different experimental conditions and methodologies. Computational strain
design relies on constraint-based reconstruction and analysis (COBRA)
techniques to predict the effect of gene knockouts such as flux balance
analysis (FBA), regulatory on/off minimization(ROOM), minimization of metabolic
adjustment (MOMA), relative optimality in metabolic networks (RELATCH). Most of
these knock-out prediction methods are based on conserving inherent flux
patterns (between wild type and mutant) that are thought to be representative
of the cellular regulatory structure. However, it has been recently
demonstrated that these methods show poor agreement with experimental data. To
improve the fidelity of knockout predictions and subsequent computational
strain design, we developed REMEP, a metabolite-centric method. We demonstrate
the improved performance of REMEP by comparing the different methods on
experimental knockout data of E. coli, and S. cerevisiae grown in batch
cultures. REMEP retains most of the features of earlier algorithms but is much
more accurate in capturing cellular responses to genetic perturbations. A
primary reason for this is that REMEP relies on the assumption that cellular
regulatory structure leaves a signature on metabolite patterns and not just
flux patterns. REMEP will also prove useful in uncovering novel insights into
cellular regulation and control.


  Background: The Epithelial-Mesenchymal Transition (EMT) endows
epithelial-looking cells with enhanced migratory ability during embryonic
development and tissue repair. EMT can also be co-opted by cancer cells to
acquire metastatic potential and drug-resistance. Recent research has argued
that epithelial (E) cells can undergo either a partial EMT to attain a hybrid
epithelial/mesenchymal (E/M) phenotype that typically displays collective
migration, or a complete EMT to adopt a mesenchymal (M) phenotype that shows
individual migration. The core EMT regulatory network -
miR-34/SNAIL/miR-200/ZEB1 - has been identified by various studies, but how
this network regulates the transitions among the E, E/M, and M phenotypes
remains controversial. Two major mathematical models - ternary chimera switch
(TCS) and cascading bistable switches (CBS) - that both focus on the
miR-34/SNAIL/miR-200/ZEB1 network, have been proposed to elucidate the EMT
dynamics, but a detailed analysis of how well either or both of these two
models can capture recent experimental observations about EMT dynamics remains
to be done. Results: Here, via an integrated experimental and theoretical
approach, we first show that both these two models can be used to understand
the two-step transition of EMT - E-E/M-M, the different responses of SNAIL and
ZEB1 to exogenous TGF-b and the irreversibility of complete EMT. Next, we
present new experimental results that tend to discriminate between these two
models. We show that ZEB1 is present at intermediate levels in the hybrid E/M
H1975 cells, and that in HMLE cells, overexpression of SNAIL is not sufficient
to initiate EMT in the absence of ZEB1 and FOXC2. Conclusions: These
experimental results argue in favor of the TCS model proposing that
miR-200/ZEB1 behaves as a three-way decision-making switch enabling transitions
among the E, hybrid E/M and M phenotypes.


  Understanding the operation of biological molecular motors, nanoscale
machines that transduce electrochemical energy into mechanical work, is
enhanced by bottom-up strategies to synthesize novel motors.


  We prove a sample path Large Deviation Principle (LDP) for a class of jump
processes whose rates are not uniformly Lipschitz continuous in phase space.
Building on it we further establish the corresponding Wentzell-Freidlin (W-F)
(infinite time horizon) asymptotic theory. These results apply to jump Markov
processes that model the dynamics of chemical reaction networks under mass
action kinetics, on a microscopic scale. We provide natural sufficient
topological conditions for the applicability of our LDP and W-F results. This
then justifies the computation of non-equilibrium potential and exponential
transition time estimates between different attractors in the large volume
limit, for systems that are beyond the reach of standard chemical reaction
network theory.


  To be able to understand how infectious diseases spread on networks, it is
important to understand the network structure itself in the absence of
infection. In this text we consider dynamic network models that are inspired by
the (static) configuration network. The networks are described by
population-level averages such as the fraction of the population with $k$
partners, $k=0,1,2,\ldots$ This means that the bookkeeping contains information
about individuals and their partners, but no information about partners of
partners. Can we average over the population to obtain information about
partners of partners? The answer is `it depends', and this is where the mean
field at distance one assumption comes into play. In this text we explain that,
yes, we may average over the population (in the right way) in the static
network. Moreover, we provide evidence in support of a positive answer for the
network model that is dynamic due to partnership changes. If, however, we
additionally allow for demographic changes, dependencies between partners
arise. In earlier work we used the slogan `mean field at distance one' as a
justification of simply ignoring the dependencies. Here we discuss the
subtleties that come with the mean field at distance one assumption, especially
when demography is involved. Particular attention is given to the accuracy of
the approximation in the setting with demography. Next, the mean field at
distance one assumption is discussed in the context of an infection
superimposed on the network. We end with the conjecture that an extension of
the bookkeeping leads to an exact description of the network structure.


  Sedimentation velocity analytical ultracentrifugation with fluorescence
detection has emerged as a powerful method for the study of interacting systems
of macromolecules. It combines picomolar sensitivity with high hydrodynamic
resolution, and can be carried out with photoswitchable fluorophores for
multi-component discrimination, to determine the stoichiometry, affinity, and
shape of macromolecular complexes with dissociation equilibrium constants from
picomolar to micromolar. A popular approach for data interpretation is the
determination of the binding affinity by isotherms of weight-average
sedimentation coefficients, sw. A prevailing dogma in sedimentation analysis is
that the weight-average sedimentation coefficient from the transport method
corresponds to the signal- and population-weighted average of all species. We
show that this does not always hold true for systems that exhibit significant
signal changes with complex formation - properties that may be readily
encountered in practice, e.g., from a change in fluorescence quantum yield.
Coupled transport in the reaction boundary of rapidly reversible systems can
make significant contributions to the observed migration in a way that cannot
be accounted for in the standard population-based average. Effective particle
theory provides a simple physical picture for the reaction-coupled migration
process. On this basis we develop a more general binding model that converges
to the well-known form of sw with constant signals, but can account
simultaneously for hydrodynamic co-transport in the presence of changes in
fluorescence quantum yield. We believe this will be useful when studying
interacting systems exhibiting fluorescence quenching, enhancement or Forster
resonance energy transfer with transport methods.


  Agent-based models (ABMs) simulate interactions between autonomous agents in
constrained environments over time. ABMs are often used for modeling the spread
of infectious diseases. In order to simulate disease outbreaks or other
phenomena, ABMs rely on "synthetic ecosystems," or information about agents and
their environments that is representative of the real world. Previous
approaches for generating synthetic ecosystems have some limitations: they are
not open-source, cannot be adapted to new or updated input data sources, and do
not allow for alternative methods for sampling agent characteristics and
locations. We introduce a general framework for generating Synthetic
Populations and Ecosystems of the World (SPEW), implemented as an open-source R
package. SPEW allows researchers to choose from a variety of sampling methods
for agent characteristics and locations when generating synthetic ecosystems
for any geographic region. SPEW can produce synthetic ecosystems for any agent
(e.g. humans, mosquitoes, etc), provided that appropriate data is available. We
analyze the accuracy and computational efficiency of SPEW given different
sampling methods for agent characteristics and locations and provide a suite of
diagnostics to screen our synthetic ecosystems. SPEW has generated over five
billion human agents across approximately 100,000 geographic regions in about
70 countries, available online.


  Despite the obvious advantage of simple life forms capable of fast
replication, different levels of cognitive complexity have been achieved by
living systems in terms of their potential to cope with environmental
uncertainty. Against the inevitable cost associated to detecting environmental
cues and responding to them in adaptive ways, we conjecture that the potential
for predicting the environment can overcome the expenses associated to
maintaining costly, complex structures. We present a minimal formal model
grounded in information theory and selection, in which successive generations
of agents are mapped into transmitters and receivers of a coded message. Our
agents are guessing machines and their capacity to deal with environments of
different complexity defines the conditions to sustain more complex agents.


  Viruses are incapable of autonomous energy production. Although many
experimental studies make it clear that viruses are parasitic entities that
hijack the host's molecular resources, a detailed estimate for the energetic
cost of viral synthesis is largely lacking. To quantify the energetic cost of
viruses to their hosts, we enumerated the costs associated with two very
distinct but representative DNA and RNA viruses, namely T4 and influenza. We
found that for these viruses, translation of viral proteins is the most
energetically expensive process. Interestingly, the cost of building a T4 phage
and a single influenza virus are nearly the same. Due to influenza's higher
burst size, however, the overall cost of a T4 phage infection is only 2-3% of
the cost of an influenza infection. The costs of these infections relative to
their host's estimated energy budget during the infection reveal that a T4
infection consumes about a third of its host's energy budget, where as an
influenza infection consumes only 1%. Building on our estimates for T4, we show
how the energetic costs of double-stranded DNA viruses scale with virus size,
revealing that the dominant cost of building a virus can switch from
translation to genome replication above a critical virus size. Lastly, using
our predictions for the energetic cost of viruses, we provide estimates for the
strengths of selection and genetic drift acting on newly incorporated genetic
elements in viral genomes, under conditions of energy limitation.


  We explore the collective behaviours of 7 group sizes: 1, 2, 3, 5, 7, 10 and
20 AB zebrafish (Danio rerio) in a constraint environment composed of two
identical squared rooms connected by a corridor. This simple set-up is similar
to a natural patchy environment. We track the positions and the identities of
the fish and compute the metrics at the group and at the individual levels.
First, we show that the size of the population affects the behaviour of each
individual in a group, the cohesion of the groups, the preferential
interactions and the transition dynamics between the two rooms. Second, during
collective departures, we show that the rankings of exit correspond to the
topological organisations of the fish prior to their collective departure with
no leadership. This spatial organisation emerge in the group a few seconds
before a collective departure. These results provide new evidences on the
spatial organisation of the groups and the effect of the population size on
individual and collective behaviours in a patchy environment.


  We study the influence of the complex topology of scale-free graphs on the
dynamics of anti-coordination games (e.g. snowdrift games). These reference
models are characterized by the coexistence (evolutionary stable mixed
strategy) of two competing species, say "cooperators" and "defectors", and, in
finite systems, by metastability and large-fluctuation-driven fixation. In this
work, we use extensive computer simulations and an effective diffusion
approximation (in the weak selection limit) to determine under which
circumstances, depending on the individual-based update rules, the topology
drastically affects the long-time behavior of anti-coordination games. In
particular, we compute the variance of the number of cooperators in the
metastable state and the mean fixation time when the dynamics is implemented
according to the voter model (death-first/birth-second process) and the link
dynamics (birth/death or death/birth at random). For the voter update rule, we
show that the scale-free topology effectively renormalizes the population size
and as a result the statistics of observables depend on the network's degree
distribution. In contrast, such a renormalization does not occur with the link
dynamics update rule and we recover the same behavior as on complete graphs.


  This thesis focuses on the applications of mathematical tools and concepts
brought from nonequilibrium statistical physics to the modeling of ecological
problems.
  The first part provides a short introduction where the theoretical concepts
and mathematical tools that are going to be used in subsequent chapters are
presented. Firstly, the different levels of description usually employed in the
models are explained. Secondly, the mathematical relationships among them are
presented. Finally, the notation and terminology that will be used later on are
explained.
  The second part is devoted to studying vegetation pattern formation in
regions where precipitations are not frequent and resources for plant growth
are scarce. This part comprises two chapters.
  The third part of the thesis develops a series of mathematical models
describing the collective movement and behavior of some animal species. Its
primary objective is to investigate the effect that communication among
foragers has on searching times and the formation of groups. It consists of two
chapters.
  The fourth part covers the effect of stochastic temporal disorder, mimicking
climate and environmental variability, on systems formed by many interacting
particles. These models may serve as an example of ecosystems. The thesis ends
with a summary and devising future research lines.


  We show how active transport of ions can be interpreted as an entropy
facilitated process. In this interpretation, the pore geometry through which
substrates are transported can give rise to a driving force. This gives a
direct link between the geometry and the changes in Gibbs energy required.
Quantifying the size of this effect for several proteins we find that the
entropic contribution from the pore geometry is significant and we discuss how
the effect can be used to interpret variations in the affinity at the binding
site.


  Genetic drift is stochastic fluctuations of alleles frequencies in a
population due to sampling effects. We consider a model of drift in an
equilibrium population, with high mutation rates: few functional mutations per
generation. Such mutation rates are common in multicellular organisms including
humans, however they are not explicitly considered in most population genetics
models. Under these assumptions the drift shows properties distinct from the
classical drift models, which ignore realistic mutation rates: i) All
(non-lethal) variants of a site have a characteristic average frequencies,
which are independent of population size, however the magnitude of fluctuations
around these frequencies depends on population size. ii) There is no
"mutational meltdown" due to "low efficiency of selection" for small population
size. Population average fitness does not depend on population size. iii) Drift
(and molecular clock) can be represented as wandering by compensatory
mutations, postulate of neutral mutations is not necessary for explaining the
high rate of mutation accumulation. Our results, which adjust the meaning of
the neutral theory from the individual neutrality of the majority of mutations,
to the collective neutrality of compensatory mutations, are applicable to
investigations in phylogeny and coalescent and for GWAS design and analysis.


  Next generation sequencing allows the identification of genes consisting of
differentially expressed transcripts, a term which usually refers to changes in
the overall expression level. A specific type of differential expression is
differential transcript usage (DTU) and targets changes in the relative within
gene expression of a transcript. The contribution of this paper is to: (a)
extend the use of cjBitSeq to the DTU context, a previously introduced Bayesian
model which is originally designed for identifying changes in overall
expression levels and (b) propose a Bayesian version of DRIMSeq, a frequentist
model for inferring DTU. cjBitSeq is a read based model and performs fully
Bayesian inference by MCMC sampling on the space of latent state of each
transcript per gene. BayesDRIMSeq is a count based model and estimates the
Bayes Factor of a DTU model against a null model using Laplace's approximation.
The proposed models are benchmarked against the existing ones using a recent
independent simulation study as well as a real RNA-seq dataset. Our results
suggest that the Bayesian methods exhibit similar performance with DRIMSeq in
terms of precision/recall but offer better calibration of False Discovery Rate.


  We introduce diffusively coupled networks where the dynamical system at each
vertex is planar Hamiltonian. The problems we address are synchronisation and
an analogue of diffusion-driven Turing instability for time-dependent
homogeneous states. As a consequence of the underlying Hamiltonian structure
there exist unusual behaviours compared with networks of coupled limit cycle
oscillators or activator-inhibitor systems.


  Plants rarely occur in isolated systems. Bacteria can inhabit either the
endosphere, the region inside the plant root, or the rhizosphere, the soil
region just outside the plant root. Our goal is to understand if using genomic
data and media dependent metabolic model information is better for training
machine learning of predicting bacterial ecological niche than media
independent models or pure genome based species trees. We considered three
machine learning techniques: support vector machine, non-negative matrix
factorization, and artificial neural networks. In all three machine-learning
approaches, the media-based metabolic models and flux balance analyses were
more effective at predicting bacterial niche than the genome or PRMT models.
Support Vector Machine trained on a minimal media base with Mannose, Proline
and Valine was most predictive of all models and media types with an f-score of
0.8 for rhizosphere and 0.97 for endosphere. Thus we can conclude that
media-based metabolic modeling provides a holistic view of the metabolome,
allowing machine learning algorithms to highlight the differences between and
categorize endosphere and rhizosphere bacteria. There was no single media type
that best highlighted differences between endosphere and rhizosphere bacteria
metabolism and therefore no single enzyme, reaction, or compound that defined
whether a bacteria's origin was of the endosphere or rhizosphere.


  Monitoring the kinetics and conformational dynamics of single enzymes is
crucial in order to better understand their biological functions as these
motions and structural dynamics are usually unsynchronized among the molecules.
Detecting the enzyme-reactant interactions and associated conformational
changes of the enzyme on a single molecule basis, however, remain as a
challenge with established optical techniques due to the commonly required
labeling of the reactants or the enzyme itself. The labeling process is usually
non-trivial and the labels themselves might skew the physical properties of the
enzyme. Here we demonstrate an optical, label-free method capable of observing
enzymatic interactions and the associated conformational changes on the single
molecule level. We monitor polymerase/DNA interactions via the strong
near-field enhancement provided by plasmonic nanorods resonantly coupled to
whispering gallery modes in microcavities. Specifically, we employ two
different recognition schemes: one in which the kinetics of polymerase/DNA
interactions are probed in the vicinity of DNA-functionalized nanorods, and the
other in which these interactions are probed via the magnitude of
conformational changes in the polymerase molecules immobilized on nanorods. In
both approaches we find that low and high polymerase activities can be clearly
discerned via their characteristic signal amplitude and signal length
distributions. Furthermore, the thermodynamic study of the monitored
interactions suggests the occurrence of DNA polymerization. This work
constitutes a proof-of-concept study of enzymatic activities via plasmonically
enhanced microcavities and establishes an alternative and label-free method
capable of investigating structural changes in single molecules.


  Trapping nanoscopic objects to observe their dynamic behaviour for extended
periods of time is an ongoing quest. Particularly, sub-100nm transparent
objects are hard to catch and most techniques rely on immobilisation or
transient diffusion through a confocal laser focus. We present an Anti-Brownian
ELectrokinetic trap (pioneered by A. E. Cohen and W. E. Moerner) to hold
nanoparticles and individual FoF1-ATP synthase proteins in solution. We are
interested in the conformational dynamics of this membrane-bound rotary motor
protein that we monitor using single-molecule FRET. The ABELtrap is an active
feedback system cancelling the nano-object's Brownian motion by applying an
electric field. We show how the induced electrokinetic forces confine the
motion of nanoparticles and proteoliposomes to the centre of the trap.


  The segregation of plasmids in a bacterial population is investigated.
Hereby, a dynamical model is formulated in terms of a size-structured
population using a hyperbolic partial differential equation incorporating
non-local terms (the fragmentation equation). For a large class of parameter
functions this PDE can be re-written as an infinite system of ordinary
differential equations for the moments of its solution. We investigate the
influence of different plasmid production modes, kinetic parameters, and
plasmid segregation modes on the equilibrium plasmid distribution. In
particular, at small plasmid numbers the distribution is strongly influenced by
the production mode, while the kinetic parameters (cell growth rate resp. basic
plasmid reproduction rate) influence the distribution mainly at large plasmid
numbers. The plasmid transmission characteristics only gradually influence the
distribution, but may become of importance for biologically relevant cases. We
compare the theoretical findings with experimental results.


  A tri-block nanoparticle (TBN) comprising of an enzymatically cleavable
porous gelatin nanocore encapsulated with gefitinib (tyrosine kinase inhibitor
(TKI)) and surface functionalized with cetuximab-siRNA conjugate has been
synthesized. Targeted delivery of siRNA to undruggable KRAS mutated non-small
cell lung cancer cells would sensitize the cells to TKI drugs and offers an
efficient therapy for treating cancer; however, efficient delivery of siRNA and
releasing it in cytoplasm remains a major challenge. We have shown TBN can
efficiently deliver siRNA to cytoplasm of KRAS mutant H23 Non-Small Cell Lung
Cancer (NSCLC) cells for oncogene knockdown; subsequently, sensitizing it to
TKI. In the absence of TKI, the nanoparticle showed minimal toxicity suggesting
that the cells adapt a parallel GAB1 mediated survival pathway. In H23 cells,
activated ERK results in phosphorylation of GAB1 on serine and threonine
residues to form GAB1-p85 PI3K complex. In the absence of TKI, knocking down
the oncogene dephosphorylated ERK, and negated the complex formation. This
event led to tyrosine phosphorylation at Tyr627 domain of GAB1 that regulated
EGFR signaling by recruiting SHP2. In the presence of TKI, GAB1-SHP2
dissociation occurs, leading to cell death. The outcome of this study provides
a promising platform for treating NSCLC patients harboring KRAS mutation.


  For animals living in groups, one of the important questions is to understand
what are the decision-making mechanisms that lead to choosing a motion
direction or leaving an area while preserving group cohesion. Here, we analyse
the initiation of collective departure in zebrafish \textit{Danio rerio}. In
particular, we observed groups of 2, 3, 5, 7 and 10 zebrafish swimming in a two
resting sites arena and quantify the number of collective departure initiated
by each fish. While all fish initiated at least one departure, the probability
to be the first one to exit a resting site is not homogeneously distributed
with some individuals leading more departures than others. We show that this
number of initiation is linearly proportional to the number of attempts
performed and that all fish have the same success rate to lead the group out of
a resting sites after an attempt. In addition, by measuring the average
swimming speed of all fish, we highlight that the intra-group ranking of a fish
for its proportion of initiation is correlated to its intra-group ranking in
average speed. These results highlight that the initiation of collective
departure in zebrafish is a heterogeneously distributed process, even if all
individual have the same success rate after attempting a departure.


  In this paper, we study how the pro-social impact due to the vigilance by
other individuals is conditioned by both environmental and evolutionary
effects. To this aim, we consider a known model where agents play a Prisoner's
Dilemma Game (PDG) among themselves and the pay-off matrix of an individual
changes according to the number of neighbors that are "vigilant", i.e., how
many neighbors watch out for her behavior. In particular, the temptation to
defect decreases linearly with the number of vigilant neighbors. This model
proved to support cooperation in specific conditions, and here we check its
robustness with different topologies, microscopical update rules and initial
conditions. By means of many numerical simulations and few theoretical
considerations, we find in which situations the vigilance by the others is more
effective in favoring cooperative behaviors and when its influence is weaker.


  The emergence and survival of cooperation is one of the hardest problems
still open in science. Several factors such as the existence of punishment,
fluctuations in finite systems, repeated interactions and the formation of
prestige may all contribute to explain the counter-intuitive prevalence of
cooperation in natural and social systems. The characteristics of the
interaction networks have been also signaled as an element favoring the
persistence of cooperators. Here we consider the invasion dynamics of
cooperative behaviors in complex topologies. The invasion of a heterogeneous
network fully occupied by defectors is performed starting from nodes with a
given number of connections (degree) $k_0$. The system is then evolved within a
Prisoner's Dilemma game and the outcome is analyzed as a function of $k_0$ and
the degree $k$ of the nodes adopting cooperation. Carried out using both
numerical and analytical approaches, our results show that the invasion
proceeds following preferentially a hierarchical order in the nodes from those
with higher degree to those with lower degree. However, the invasion of
cooperation will succeed only when the initial cooperators are numerous enough
to form a cluster from which cooperation can spread. This implies that the
initial condition must be a suitable equilibrium between high degree and high
numerosity, which usually takes place, when possible, at intermediate values of
$k_0$. These findings have potential applications, as they suggest that, in
order to promote cooperative behavior on complex networks, one should infect
with cooperators \emph{high but not too high} degree nodes.


  The provision of intergenerational care, via the Grandmother Hypothesis, has
been implicated in the evolution of post-fertile longevity, particularly in
humans. However, if grandmothering does provide fitness benefits, a key
question is why has it evolved so infrequently? We investigate this question
with a combination of life-history and evolutionary game theory. We derive
simple eligibility and stability thresholds, both of which must be satisfied if
intergenerational care is first to evolve and then to persist in a population.
As one threshold becomes easier to fulfill, the other becomes more difficult,
revealing a conflict between the two. As such, we suggest that, in fact, we
should expect the evolution of grandmothering to be rare.


  Epidemic spreading on complex networks depends on the topological structure
as well as on the dynamical properties of the infection itself. Generally
speaking, highly connected individuals play the role of hubs and are crucial to
channel information across the network. On the other hand, static topological
quantities measuring the connectivity structure are independent on the
dynamical mechanisms of the infection. A natural question is therefore how to
improve the topological analysis by some kind of dynamical information that may
be extracted from the ongoing infection itself. In this spirit, we propose a
novel vaccination scheme that exploits information from the details of the
infection pattern at the moment when the vaccination strategy is applied.
Numerical simulations of the infection process show that the proposed
immunization strategy is effective and robust on a wide class of complex
networks.


  The concept of allostery in which macromolecules switch between two different
conformations is a central theme in biological processes ranging from gene
regulation to cell signaling to enzymology. Allosteric enzymes pervade
metabolic processes, yet a simple and unified treatment of the effects of
allostery in enzymes has been lacking. In this work, we take the first step
towards this goal by modeling allosteric enzymes and their interaction with two
key molecular players - allosteric regulators and competitive inhibitors. We
then apply this model to characterize existing data on enzyme activity, comment
on how enzyme parameters (such as substrate binding affinity) can be
experimentally tuned, and make novel predictions on how to control phenomena
such as substrate inhibition.


  While all organisms on Earth descend from a common ancestor, there is no
consensus on whether the origin of this ancestral self-replicator was a one-off
event or whether it was only the final survivor of multiple origins. Here we
use the digital evolution system Avida to study the origin of self-replicating
computer programs. By using a computational system, we avoid many of the
uncertainties inherent in any biochemical system of self-replicators (while
running the risk of ignoring a fundamental aspect of biochemistry). We
generated the exhaustive set of minimal-genome self-replicators and analyzed
the network structure of this fitness landscape. We further examined the
evolvability of these self-replicators and found that the evolvability of a
self-replicator is dependent on its genomic architecture. We studied the
differential ability of replicators to take over the population when competed
against each other (akin to a primordial-soup model of biogenesis) and found
that the probability of a self-replicator out-competing the others is not
uniform. Instead, progenitor (most-recent common ancestor) genotypes are
clustered in a small region of the replicator space. Our results demonstrate
how computational systems can be used as test systems for hypotheses concerning
the origin of life.


  The rapid expansion in the spectrum of two-dimensional (2D) materials has
driven the efforts of research on the fabrication of 2D composites and
heterostructures. Highly ordered structure of 2D materials provides an
excellent platform for controlling the ultimate structure and properties of the
composite material with precision. However, limited control over the structure
of the adherent material and its interactions with highly ordered 2D materials
results in defective composites with inferior performance. Here, we demonstrate
the successful synthesis, integration, and characterization of hybrid 2D
material systems consisting of tandem repeat (TR) proteins inspired by squid
ring teeth and graphene oxide (GO). The TR protein layer acts as a unique
programmable molecular spacer between GO layers. As an application, we further
demonstrate thermal actuation using bimorph molecular composite films. Bimorph
actuators made of molecular composite films (GO/TR) can lead to energy
efficiencies 18 times higher than regular bimorph actuators consisting of a GO
layer and a TR protein layer (i.e., conventional bulk composite of GO and TR).
Additionally, molecular composite bimorph actuators can reach curvature values
as high as 1.2 cm-1 by using TR proteins with higher molecular weight, which is
3 times higher than conventional GO and TR composites.


  The spindle checkpoint assembly (SAC) ensures genome fidelity by temporarily
delaying anaphase onset, until all chromosomes are properly attached to the
mitotic spindle. The SAC delays mitotic progression by preventing activation of
the ubiquitin ligase anaphase-promoting complex (APC/C) or cyclosome; whose
activation by Cdc20 is required for sister-chromatid separation marking the
transition into anaphase. The mitotic checkpoint complex (MCC), which contains
Cdc20 as a subunit, binds stably to the APC/C. Compelling evidence by Izawa and
Pines (Nature 2014; 10.1038/nature13911) indicates that the MCC can inhibit a
second Cdc20 that has already bound and activated the APC/C. Whether or not MCC
per se is sufficient to fully sequester Cdc20 and inhibit APC/C remains
unclear. Here, a dynamic model for SAC regulation in which the MCC binds a
second Cdc20 was constructed. This model is compared to the MCC, and the
MCC-and-BubR1 (dual inhibition of APC) core model variants and subsequently
validated with experimental data from the literature. By using ordinary
nonlinear differential equations and spatial simulations, it is shown that the
SAC works sufficiently to fully sequester Cdc20 and completely inhibit APC/C
activity. This study highlights the principle that a systems biology approach
is vital for molecular biology and could also be used for creating hypotheses
to design future experiments.


  We are living in an uncertain and dynamically changing world, where optimal
decision-making under uncertainty is directly linked to the survival of
species. However, evolutionary selection pressures that shape value-based
decision-making under uncertainty have thus far received limited attention.
Here, we demonstrate that fitness associated with different value-based
decision-making preferences is influenced by the value properties of the
environment, as well as the characteristics and the density of competitors in
the population. We show that risk-seeking tendencies will eventually dominate
the population, when there are a relatively large number of discrete strategies
competing in volatile value environments. These results may have important
implications for behavioural ecology: (i) to inform the prediction that species
which naturally exhibit risk-averse characteristics and live alongside
risk-seeking competitors may be selected against; (ii) to potentially improve
our understanding of day-traders value-based decision-making preferences in
volatile financial markets in terms of an environmental adaptation.


  Proteins are biological polymers that underlie all cellular functions. The
first high-resolution protein structures were determined by x-ray
crystallography in the 1960s. Since then, there has been continued interest in
understanding and predicting protein structure and stability. It is
well-established that a large contribution to protein stability originates from
the sequestration from solvent of hydrophobic residues in the protein core. How
are such hydrophobic residues arranged in the core? And how can one best model
the packing of these residues? Here we show that to properly model the packing
of residues in protein cores it is essential that amino acids are represented
by appropriately calibrated atom sizes, and that hydrogen atoms are explicitly
included. We show that protein cores possess a packing fraction of $\phi
\approx 0.56$, which is significantly less than the typically quoted value of
0.74 obtained using the extended atom representation. We also compare the
results for the packing of amino acids in protein cores to results obtained for
jammed packings from disrete element simulations composed of spheres, elongated
particles, and particles with bumpy surfaces. We show that amino acids in
protein cores pack as densely as disordered jammed packings of particles with
similar values for the aspect ratio and bumpiness as found for amino acids.
Knowing the structural properties of protein cores is of both fundamental and
practical importance. Practically, it enables the assessment of changes in the
structure and stability of proteins arising from amino acid mutations (such as
those identified as a result of the massive human genome sequencing efforts)
and the design of new folded, stable proteins and protein-protein interactions
with tunable specificity and affinity.


  We call change-point problem (CPP) the identification of changes in the
probabilistic behavior of a sequence of observations. Solving the CPP involves
detecting the number and position of such changes. In genetics the study of how
and what characteristics of a individual's genetic content might contribute to
the occurrence and evolution of cancer has fundamental importance in the
diagnosis and treatment of such diseases and can be formulated in the framework
of chage-point analysis. In this article we propose a modification to a
existing method of segmentation with the objective of producing a algorithm
that is robust to a variety of sampling distributions and that is adequate for
more recent method of accessing DNA copy-number which might require a
restriction on the minimum length of a altered segment.


  We present a generalized Landau-Brazovskii free energy for the solidification
of chiral molecules on a spherical surface in the context of the assembly of
viral shells. We encounter two types of icosahedral solidification transitions.
The first type is a conventional first-order phase transition from the uniform
to the icosahedral state. It can be described by a single icosahedral spherical
harmonic of even $l$. The chiral pseudo-scalar term in the free energy creates
secondary terms with chiral character but it does not affect the thermodynamics
of the transition. The second type, associated with icosahedral spherical
harmonics with odd $l$, is anomalous. Pure odd $l$ icosahedral states are
unstable but stability is recovered if admixture with the neighboring $l+1$
icosahedral spherical harmonic is included, generated by the non-linear terms.
This is in conflict with the principle of Landau theory that symmetry-breaking
transitions are characterized by only a \textit{single} irreducible
representation of the symmetry group of the uniform phase and we argue that
this principle should be removed from Landau theory. The chiral term now
directly affects the transition because it lifts the degeneracy between two
isomeric mixed-$l$ icosahedral states. A direct transition is possible only
over a limited range of parameters. Outside this range, non-icosahedral states
intervene. For the important case of capsid assembly dominated by $l=15$, the
intervening states are found to be based on octahedral symmetry.


  For decades, dengue virus has been a cause of major public health concern in
Costa Rica, due to its landscape and climatic conditions that favor the
circumstances in which the vector, Aedes aegypti, thrives. The emergence and
introduction throughout tropical and subtropical countries of the chikungunya
virus, as of 2014, challenged Costa Rican health authorities to provide a
correct diagnosis since it is also transmitted by the same vector and infected
hosts may share similar symptoms. We study the 2015-2016 dengue and chikungunya
outbreaks in Costa Rica while establishing how point estimates of epidemic
parameters for both diseases compare to one another. Longitudinal weekly
incidence reports of these outbreaks signal likely misdiagnosis of infected
individuals: underreporting of chikungunya cases, while overreporting cases of
dengue. Our comparative analysis is formulated with a single-outbreak
deterministic model that features an undiagnosed class. Additionally, we also
used a genetic algorithm in the context of weighted least squares to calculate
point estimates of key model parameters and initial conditions, while formally
quantifying misdiagnosis.


  In this paper, a brief review of delay population models and their
applications in ecology is provided. The inclusion of diffusion and nonlocality
terms in delay models has given more capabilities to these models enabling them
to capture several ecological phenomena such as the Allee effect, waves of
invasive species and spatio-temporal competitions of interacting species.
Moreover, recent advances in the studies of traveling and stationary wave
solutions of delay models are outlined. In particular, the existence of
stationary and traveling wave solutions of delay models, stability of wave
solutions, formation of wavefronts in the special domain, and possible outcomes
of delay models are discussed.