  Lineage tracing, the determination and mapping of progeny arising from single
cells, is an important approach enabling the elucidation of mechanisms
underlying diverse biological processes ranging from development to disease. We
developed a dynamic sequence-based barcode for lineage tracing and have
demonstrated its performance in C. elegans, a model organism whose lineage tree
is well established. The strategy we use creates lineage trees based upon the
introduction of specific mutations into cells and the propagation of these
mutations to daughter cells at each cell division. We present an experimental
proof of concept along with a corresponding simulation and analytical model for
deeper understanding of the coding capacity of the system. By introducing
mutations in a predictable manner using CRISPR/Cas9, our technology will enable
more complete investigations of cellular processes.


  Hemoglobin (Hgb) forms tetramers (dimerized dimers), which enhance its
globular stability and may also facilitate small gas molecule transport, as
shown by recent all-atom Newtonian solvated simulations. Hydropathic
bioinformatic scaling reveals many wave-like features of strained Hgb
structures at the coarse-grained amino acid level, while distinguishing between
these features thermodynamically. Strain fields localized near hemes interfere
with extended strain fields associated with dimer interfacial misfit, resulting
in wave-length dependent dimer correlation function antiresonances.


  Genetic regulatory networks with adaptive responses are widely studied in
biology. Usually, models consisting only of a few nodes have been considered.
They present one input receptor for activation and one output node where the
adaptive response is computed. In this work, we design genetic regulatory
networks with many receptors and many output nodes able to produce delayed
adaptive responses. This design is performed by using an evolutionary algorithm
of mutations and selections that minimizes an error function defined by the
adaptive response in signal shapes. We present several examples of network
constructions with a predefined required set of adaptive delayed responses. We
show that an output node can have different kinds of responses as a function of
the activated receptor. Additionally, complex network structures are presented
since processing nodes can be involved in several input-output pathways.


  The dynamical organization of membrane-bound organelles along intracellular
transport pathways relies on vesicular exchange between organelles and on
biochemical maturation of the organelle content by specific enzymes. The
relative importance of each mechanism in controlling organelle dynamics remains
controversial, in particular for transport through the Golgi apparatus. Using a
stochastic model, we show that full maturation of membrane-bound compartments
can be seen as the stochastic escape from a steady-state in which export is
dominated by vesicular exchange. We show that full maturation can contribute a
significant fraction of the total out-flux for small organelles such as
endosomes and Golgi cisternae.


  We describe modeling approaches to a "network" of connected enzyme-catalyzed
reactions, with added (bio)chemical processes that introduce biochemical
filtering steps into the functioning of such a biocatalytic cascade.
Theoretical expressions are derived that allow simple, few-parameter modeling
of processes concatenated in such cascades, both with and without filtering.
The modeling approach captures and explains features identified in earlier
studies of enzymatic processes considered as potential "network components" for
multi-step information/signal processing systems.


  Processes involving multi-input multi-step reaction cascades are used in
developing novel biosensing, biocomputing, and decision making systems. In
various applications different changes in responses of the constituent
processing steps (reactions) in a cascade are desirable in order to allow
control of the system's response. Here we consider conversion of convex
response to sigmoid by "intensity filtering," as well as "threshold filtering,"
and we offer a general overview of this field of research. Specifically, we
survey rate equation modelling that has been used for enzymatic reactions. This
allows us to design modified biochemical processes as "network components" with
responses desirable in applications.


  The propensity of short DNA sequences to convert to the circular form is
studied by a mesoscopic Hamiltonian method which incorporates both the bending
of the molecule axis and the intrinsic twist of the DNA strands. The base pair
fluctuations with respect to the helix diameter are treated as path
trajectories in the imaginary time path integral formalism. The partition
function for the sub-ensemble of closed molecules is computed by imposing chain
ends boundary conditions both on the radial fluctuations and on the angular
degrees of freedom. The cyclization probability, the J-factor, proves to be
highly sensitive to the stacking potential, mostly to its nonlinear parameters.
We find that the J-factor generally decreases by reducing the sequence length (
N ) and, more significantly, below N = 100 base pairs. However, even for very
small molecules, the J-factors remain sizeable in line with recent experimental
indications. Large bending angles between adjacent base pairs and anharmonic
stacking appear as the causes of the helix flexibility at short length scales.


  Following the interest of L Messori and A Merlino 2016 Coordination Chemistry
Reviews in the platinum ions coordination geometries in our PDB entries 4dd4
and 4dd6 we have extended our original analyses.


  Iterated games provide a framework to describe social interactions among
groups of individuals. Recent work stimulated by the discovery of
"zero-determinant" strategies has rapidly expanded our ability to analyze such
interactions. This body of work has primarily focused on games in which players
face a simple binary choice, to "cooperate" or "defect". Real individuals,
however, often exhibit behavioral diversity, varying their input to a social
interaction both qualitatively and quantitatively. Here we explore how access
to a greater diversity of behavioral choices impacts the evolution of social
dynamics in finite populations. We show that, in public goods games, some
two-choice strategies can nonetheless resist invasion by all possible
multi-choice invaders, even while engaging in relatively little punishment. We
also show that access to greater behavioral choice results in more "rugged "
fitness landscapes, with populations able to stabilize cooperation at multiple
levels of investment, such that choice facilitates cooperation when returns on
investments are low, but hinders cooperation when returns on investments are
high. Finally, we analyze iterated rock-paper-scissors games, whose
non-transitive payoff structure means unilateral control is difficult and
zero-determinant strategies do not exist in general. Despite this, we find that
a large portion of multi-choice strategies can invade and resist invasion by
strategies that lack behavioral diversity -- so that even well-mixed
populations will tend to evolve behavioral diversity.


  Background: In 2015, the Zika arbovirus (ZIKV) began circulating in the
Americas, rapidly expanding its global geographic range in explosive outbreaks.
Unusual among mosquito-borne diseases, ZIKV has been shown to also be sexually
transmitted, although sustained autochthonous transmission due to sexual
transmission alone has not been observed, indicating the reproduction number
(R0) for sexual transmission alone is less than 1. Critical to the assessment
of outbreak risk, estimation of the potential attack rates, and assessment of
control measures, are estimates of the basic reproduction number, R0. Methods:
We estimated the R0 of the 2015 ZIKV outbreak in Barranquilla, Colombia,
through an analysis of the exponential rise in clinically identified ZIKV cases
(n = 359 to the end of November, 2015). Findings: The rate of exponential rise
in cases was rho=0.076 days-1, with 95 percent CI [0.066,0.087] days-1. We used
a vector-borne disease model with additional direct transmission to estimate
the R0; assuming the R0 of sexual transmission alone is less than 1, we
estimated the total R0 = 3.8 [2.4,5.6], and that the fraction of cases due to
sexual transmission was 0.23 [0.01,0.47] with 95 percent confidence.
Interpretation: This is among the first estimates of R0 for a ZIKV outbreak in
the Americas, and also among the first quantifications of the relative impact
of sexual transmission.


  A central issue in cell biology is the physico-chemical basis of organelle
biogenesis in intracellular trafficking pathways, its most impressive
manifestation being the biogenesis of Golgi cisternae. At a basic level, such
morphologically and chemically distinct compartments should arise from an
interplay between the molecular transport and chemical maturation. Here, we
formulate analytically tractable, minimalist models, that incorporate this
interplay between transport and chemical progression in physical space, and
explore the conditions for de novo biogenesis of distinct cisternae. We propose
new quantitative measures that can discriminate between the various models of
transport in a qualitative manner- this includes measures of the dynamics in
steady state and the dynamical response to perturbations of the kind amenable
to live-cell imaging.


  Understanding and even defining what constitutes animal interactions remains
a challenging problem. Correlational tools may be inappropriate for detecting
communication between a set of many agents exhibiting nonlinear behavior. A
different approach is to define coordinated motions in terms of an information
theoretic channel of direct causal information flow. In this work, we consider
time series data obtained by an experimental protocol of optical tracking of
the insect species Chironomus riparius. The data constitute reconstructed 3-D
spatial trajectories of the insects' flight trajectories and kinematics. We
present an application of the optimal causation entropy (oCSE) principle to
identify direct causal relationships or information channels among the insects.
The collection of channels inferred by oCSE describes a network of information
flow within the swarm. We find that information channels with a long spatial
range are more common than expected under the assumption that causal
information flows should be spatially localized. The tools developed herein are
general and applicable to the inference and study of intercommunication
networks in a wide variety of natural settings.


  Why is the spine of a neuron so small that only small numbers of molecules
can exist and reactions inevitably become stochastic? Despite such noisy
conditions, we previously showed that the spine exhibits robust, sensitive and
efficient features of information transfer using probability of Ca$^{2+}$
increase; however, their mechanisms remains unknown. Here we show that the
small-volume effect enables robust, sensitive and efficient information
transfer in the spine volume, but not in the cell volume. In the spine volume,
intrinsic noise in reactions becomes larger than extrinsic noise of input,
making robust information transfer against input fluctuation. Stochastic
facilitation of Ca$^{2+}$ increase occurs in the spine volume, making higher
sensitivity to lower intensity of input. Volume-dependency of information
transfer enables efficient information transfer per input in the spine volume.
Thus, we propose that the small-volume effect is the functional reasons why the
spine has to be so small.


  We derive structural and binding energy trends for twenty amino acids, their
dipeptides, and their interactions with the divalent cations Ca$^{2+}$,
Ba$^{2+}$, Sr$^{2+}$, Cd$^{2+}$, Pb$^{2+}$, and Hg$^{2+}$. The underlying data
set consists of 45,892 first-principles predicted conformers with relative
energies up to about 4 eV (about 400kJ/mol). We show that only very few
distinct backbone structures of isolated amino acids and their dipeptides
emerge as lowest-energy conformers. The isolated amino acids predominantly
adopt structures that involve an acidic proton shared between the carboxy and
amino function. Dipeptides adopt one of two intramolecular-hydrogen bonded
conformations C$_5$ or equatorial C$_7$. Upon complexation with a divalent
cation, the accessible conformational space shrinks and intramolecular hydrogen
bonding is prevented due to strong electrostatic interaction of backbone and
side chain functional groups with cations. Clear correlations emerge from the
binding energies of the six divalent ions with amino acids and dipeptides.
Cd$^{2+}$ and Hg$^{2+}$ show the largest binding energies - a potential
correlation with their known high acute toxicities. Ca$^{2+}$ and Pb$^{2+}$
reveal almost identical binding energies across the entire series of amino
acids and dipeptides. This observation validates past indications that
ion-mimicry of calcium and lead should play an important role in a
toxicological context.


  Cytoskeletal motor proteins are involved in major intracellular transport
processes which are vital for maintaining appropriate cellular function. The
motor exhibits distinct states of motility: active motion along filaments, and
effectively stationary phase in which it detaches from the filaments and
performs passive diffusion in the vicinity of the detachment point due to
cytoplasmic crowding. The transition rates between motion and pause phases are
asymmetric in general, and considerably affected by changes in environmental
conditions which influences the efficiency of cargo delivery to specific
targets. By considering the motion of molecular motor on a single filament as
well as a dynamic filamentous network, we present an analytical model for the
dynamics of self-propelled particles which undergo frequent pause phases. The
interplay between motor processivity, structural properties of filamentous
network, and transition rates between the two states of motility drastically
changes the dynamics: multiple transitions between different types of anomalous
diffusive dynamics occur and the crossover time to the asymptotic diffusive or
ballistic motion varies by several orders of magnitude. We map out the phase
diagrams in the space of transition rates, and address the role of initial
conditions of motion on the resulting dynamics.


  We discuss how to use evolutionary game theory (EGT) as a framework for
studying how cultural dynamics and structural properties can influence the
evolution of norms and behaviors within a society. We provide a brief tutorial
on how EGT works, and discuss what kinds of insights it can provide. We then
describe three published studies in which we have developed EGT models that
help explain how structural and external conditions in a society affect the
emergence of social norms.


  Multicellular organisms consist of cells that interact via elaborate adhesion
complexes. Desmosomes are membrane-associated adhesion complexes that
mechanically tether the cytoskeletal intermediate filaments (IFs) between two
adjacent cells, creating a network of tough connections in tissues such as skin
and heart. Desmoplakin (DP) is the key desmosomal protein that binds IFs, and
the DP-IF association poses a quandary: desmoplakin must stably and tightly
bind IFs to maintain the structural integrity of the desmosome. Yet, newly
synthesized DP must traffick along the cytoskeleton to the site of nascent
desmosome assembly without 'sticking' to the IF network, implying weak or
transient DP--IF contacts. Recent work reveals that these contacts are
modulated by post-translational modifications (PTMs) in DP's C-terminal tail.
Using molecular dynamics simulations, we have elucidated the structural basis
of these PTM-induced effects. Our simulations, nearing 2 microseconds in
aggregate, indicate that phosphorylation of S2849 induces an 'arginine claw' in
desmoplakin's C-terminal tail (DPCTT). If a key arginine, R2834, is methylated,
the DPCTT preferentially samples conformations that are geometrically
well-suited as substrates for processive phosphorylation by the cognate kinase
GSK3. We suggest that DPCTT is a molecular switch that modulates, via its
conformational dynamics, DP's efficacy as a substrate for GSK3. Finally, we
show that the fluctuating DPCTT can contact other parts of DP, suggesting a
competitive binding mechanism for the modulation of DP--IF interactions.


  We present a generalized Landau-Brazovskii theory for the solidification of
chiral molecules on a spherical surface. With increasing sphere radius one
encounters first intervals where robust achiral density modulations appear with
icosahedral symmetry via first-order transitions. Next, one en- counters
intervals where fragile but stable icosahedral structures still can be
constructed but only by superposition of multiple irreducible representations.
Chiral icoshedral structures appear via continuous or very weakly first-order
transitions. Outside these parameter intervals, icosahedral symmetry is broken
along a three-fold axis or a five-fold axis. The predictions of the theory are
compared with recent numerical simulations.


  Neutral models, in which individual agents with equal fitness undergo a
birth-death-mutation process, are very popular in population genetics and
community ecology. Usually these models are applied to populations and
communities with spatial structure, but the analytic results presented so far
are limited to well-mixed or mainland-island scenarios. Here we present a new
technique, based on interface dynamics analysis, and apply it to the neutral
dynamics in one, two and three spatial dimensions. New results are derived for
the correlation length and for the main characteristics of the community, like
total biodiversity and the species abundance distribution above the correlation
length. Our results are supported by extensive numerical simulations, and
provide qualitative and quantitative insights that allow for a rigorous
comparison between model predictions and empirical data.


  We consider continuous time Hopfield-like recurrent networks as dynamical
models for gene regulation and neural networks. We are interested in networks
that contain n high-degree nodes preferably connected to a large number of Ns
weakly connected satellites, a property that we call n/Ns-centrality. If the
hub dynamics is slow, we obtain that the large time network dynamics is
completely defined by the hub dynamics. Moreover, such networks are maximally
flexible and switchable, in the sense that they can switch from a globally
attractive rest state to any structurally stable dynamics when the response
time of a special controller hub is changed. In particular, we show that a
decrease of the controller hub response time can lead to a sharp variation in
the network attractor structure: we can obtain a set of new local attractors,
whose number can increase exponentially with N, the total number of nodes of
the nework. These new attractors can be periodic or even chaotic. We provide an
algorithm, which allows us to design networks with the desired switching
properties, or to learn them from time series, by adjusting the interactions
between hubs and satellites. Such switchable networks could be used as models
for context dependent adaptation in functional genetics or as models for
cognitive functions in neuroscience.


  Metabolomic data sets provide a direct read-out of cellular phenotypes and
are increasingly generated to study biological questions. Our previous work
revealed the potential of analyzing extracellular metabolomic data in the
context of the metabolic model using constraint-based modeling. Through this
work, which consists of a protocol, a toolbox, and tutorials of two use cases,
we make our methods available to the broader scientific community. The protocol
describes, in a step-wise manner, the workflow of data integration and
computational analysis. The MetaboTools comprise the Matlab code required to
complete the workflow described in the protocol. Tutorials explain the
computational steps for integration of two different data sets and demonstrate
a comprehensive set of methods for the computational analysis of metabolic
models and stratification thereof into different phenotypes. The presented
workflow supports integrative analysis of multiple omics data sets.
Importantly, all analysis tools can be applied to metabolic models without
performing the entire workflow. Taken together, this protocol constitutes a
comprehensive guide to the intra-model analysis of extracellular metabolomic
data and a resource offering a broad set of computational analysis tools for a
wide biomedical and non-biomedical research community.


  It is common in phylogenetics to have some, perhaps partial, information
about the overall evolutionary tree of a group of organisms and wish to find an
evolutionary tree of a specific gene for those organisms. There may not be
enough information in the gene sequences alone to accurately reconstruct the
correct "gene tree." Although the gene tree may deviate from the "species tree"
due to a variety of genetic processes, in the absence of evidence to the
contrary it is parsimonious to assume that they agree. A common statistical
approach in these situations is to develop a likelihood penalty to incorporate
such additional information. Recent studies using simulation and empirical data
suggest that a likelihood penalty quantifying concordance with a species tree
can significantly improve the accuracy of gene tree reconstruction compared to
using sequence data alone. However, the consistency of such an approach has not
yet been established, nor have convergence rates been bounded. Because
phylogenetics is a non-standard inference problem, the standard theory does not
apply. In this paper, we propose a penalized maximum likelihood estimator for
gene tree reconstruction, where the penalty is the square of the
Billera-Holmes-Vogtmann geodesic distance from the gene tree to the species
tree. We prove that this method is consistent, and derive its convergence rate
for estimating the discrete gene tree structure and continuous edge lengths
(representing the amount of evolution that has occurred on that branch)
simultaneously. We find that the regularized estimator is "adaptive fast
converging," meaning that it can reconstruct all edges of length greater than
any given threshold from gene sequences of polynomial length. Our method does
not require the species tree to be known exactly; in fact, our asymptotic
theory holds for any such guide tree.


  Gene regulatory networks can be successfully modeled as Boolean networks. A
much discussed hypothesis says that such model networks reproduce empirical
findings the best if they are tuned to operate at criticality, i.e. at the
borderline between their ordered and disordered phases. Critical networks have
been argued to lead to a number of functional advantages such as maximal
dynamical range, maximal sensitivity to environmental changes, as well as to an
excellent trade off between stability and flexibility. Here, we study the
effect of noise within the context of Boolean networks trained to learn complex
tasks under supervision. We verify that quasi-critical networks are the ones
learning in the fastest possible way --even for asynchronous updating rules--
and that the larger the task complexity the smaller the distance to
criticality. On the other hand, when additional sources of intrinsic noise in
the network states and/or in its wiring pattern are introduced, the optimally
performing networks become clearly subcritical. These results suggest that in
order to compensate for inherent stochasticity, regulatory and other type of
biological networks might become subcritical rather than being critical, all
the most if the task to be performed has limited complexity.


  Genetic data are frequently categorical and have complex dependence
structures that are not always well understood. For this reason, clustering and
classification based on genetic data, while highly relevant, are challenging
statistical problems. Here we consider a highly versatile U-statistics based
approach built on dissimilarities between pairs of data points for
nonparametric clustering. In this work we propose statistical tests to assess
group homogeneity taking into account the multiple testing issues, and a
clustering algorithm based on dissimilarities within and between groups that
highly speeds up the homogeneity test. We also propose a test to verify
classification significance of a sample in one of two groups. A Monte Carlo
simulation study is presented to evaluate power of the classification test,
considering different group sizes and degree of separation. Size and power of
the homogeneity test are also analyzed through simulations that compare it to
competing methods. Finally, the methodology is applied to three different
genetic datasets: global human genetic diversity, breast tumor gene expression
and Dengue virus serotypes. These applications showcase this statistical
framework's ability to answer diverse biological questions while adapting to
the specificities of the different datatypes.


  This paper explores the impact of surface area, volume, curvature and
Lennard-Jones potential on solvation free energy predictions. Rigidity surfaces
are utilized to generate robust analytical expressions for maximum, minimum,
mean and Gaussian curvatures of solvent-solute interfaces, and define a
generalized Poisson-Boltzmann (GPB) equation with a smooth dielectric profile.
Extensive correlation analysis is performed to examine the linear dependence of
surface area, surface enclosed volume, maximum curvature, minimum curvature,
mean curvature and Gaussian curvature for solvation modeling. It is found that
surface area and surfaces enclosed volumes are highly correlated to each
others, and poorly correlated to various curvatures for six test sets of
molecules. Different curvatures are weakly correlated to each other for six
test sets of molecules, but are strongly correlated to each other within each
test set of molecules. Based on correlation analysis, we construct twenty six
nontrivial nonpolar solvation models. Our numerical results reveal that the
Lennard-Jones (LJ) potential plays a vital role in nonpolar solvation modeling,
especially for molecules involving strong van der Waals interactions. It is
found that curvatures are at least as important as surface area or surface
enclosed volume in nonpolar solvation modeling. In conjugation with the GPB
model, various curvature based nonpolar solvation models are shown to offer
some of the best solvation free energy predictions for a wide range of test
sets. For example, root mean square errors from a model constituting surface
area, volume, mean curvature and LJ potential are less than 0.42 kcal/mol for
all test sets.


  The flexibility of short DNA fragments is studied by a Hamiltonian model
which treats the inter-strand and intra-strand forces at the level of the base
pair. The elastic response of a set of homogeneous helices to externally
applied forces is obtained by computing the average bending angles between
adjacent base pairs along the molecule axis. The ensemble averages are
performed over a room temperature equilibrium distribution of base pair
separations and bending fluctuations. The analysis of the end-to-end distances
and persistence lengths shows that even short sequences with less than $100$
base pairs maintain a significant bendability ascribed to thermal fluctuational
effects and kinks with large bending angles. The discrepancies between the
outcomes of the discrete model and those of the worm-like-chain model are
examined pointing out the inadequacy of the latter on short length scales.


  Sm proteins were discovered nearly 20 years ago as a group of small antigenic
proteins ($\approx$ 90-120 residues). Since then, an extensive amount of
biochemical and genetic data have illuminated the crucial roles of these
proteins in forming ribonucleoprotein (RNP) complexes that are used in RNA
processing, e.g., spliceosomal removal of introns from pre-mRNAs. Spliceosomes
are large macromolecular machines that are comparable to ribosomes in size and
complexity, and are composed of uridine-rich small nuclear RNPs (U snRNPs).
Various sets of seven different Sm proteins form the cores of most snRNPs.
Despite their importance, very little is known about the atomic-resolution
structure of snRNPs or their Sm cores. As a first step towards a
high-resolution image of snRNPs and their hierarchic assembly, we have
determined the crystal structures of archaeal homologs of Sm proteins, which we
term Sm-like archaeal proteins (SmAPs).


  This document attempts to clarify potential confusion regarding
electrostatics calculations, specifically in the context of biomolecular
structure and specifically as regards the units typically used to
contour/visualize isopotential surfaces, potentials mapped onto molecular
solvent-accessible surfaces, etc.


  In this paper I have given a mathematical model of Cell reprogramming from a
different contexts. Here I considered there is a delay in differential
regulator rate equations due to intermediate regulator's regulations. At first
I gave some basic mathematical models by Ferell Jr.[2] of reprogramming and
after that I gave mathematical model of cell reprogramming by Mithun Mitra[4].
In the last section I contributed a mathematical model of cell reprogramming
from intermediate steps regulations and tried to find the critical point of
pluripotent cell.


  The dynamics of populations is frequently subject to intrinsic noise. At the
same time unknown interaction networks or rate constants can present quenched
uncertainty. Existing approaches often involve repeated sampling of the
quenched disorder and then running the stochastic birth-death dynamics on these
samples. In this paper we take a different view, and formulate an effective
jump process, representative of the ensemble of quenched interactions as a
whole. Using evolutionary games with random payoff matrices as an example, we
develop an algorithm to simulate this process, and we discuss diffusion
approximations in the limit of weak intrinsic noise.


  This paper presents an intracellular feedback control strategy, to regulate
the gene expression process dynamics. For this purpose, two types of genetic
circuits are designed in order to compare concentrations of the input
transcription factor (the desired input) and the protein produced by the
expression of target gene (process output). Genetic Relay Switch and Shifted
Subtractor are proposed and employed to design Discrete Comparator. Also,
another type of genetic comparator is presented named Bistable Comparator.
These synthetic comparators are used to produce the required amount of the
protein that activates the expression of the target gene. Numerical simulation
results, demonstrate the effectiveness of the designed controllers.


  Identification and characterization of large protein complexes is a mainstay
of biochemical toolboxes. Utilization of cross-linking chemicals can facilitate
the capture and identification of transient or weak interactions of a transient
nature. Here we describe a detailed methodology for cell culture based
proteomic approach. We describe the generation of cells stably expressing green
fluorescent protein (GFP)- tagged proteins under the tetracycline-inducible
promoter and subsequent proteomic analysis of GFP-interacting proteins. We
include a list of proteins that were identified as interactors of GFP.


  Assembly and stability of mitotic spindle is governed by the interplay of
various intra-cellular forces, e.g. the forces generated by motor proteins by
sliding overlapping anti-parallel microtubules (MTs) polymerized from the
opposite centrosomes, the interaction of kinetochores with MTs, and the
interaction of MTs with the chromosomes arms. We study the mechanical behavior
and stability of spindle assembly within the framework of a minimal model which
includes all these effects. For this model, we derive a closed--form analytical
expression for the force acting between the centrosomes as a function of their
separation distance and we show that an effective potential can be associated
with the interactions at play. We obtain the stability diagram of spindle
formation in terms of parameters characterizing the strength of motor sliding,
repulsive forces generated by polymerizing MTs, and the forces arising out of
interaction of MTs with kinetochores. The stability diagram helps in
quantifying the relative effects of the different interactions and elucidates
the role of motor proteins in formation and inhibition of spindle structures
during mitotic cell division. We also predict a regime of bistability for
certain parameter range, wherein the spindle structure can be stable for two
different finite separation distances between centrosomes. This occurrence of
bistability also suggests mechanical versatility of such self-assembled spindle
structures.


  The adhesion of cell membranes is mediated by the binding of
membrane-anchored receptor and ligand proteins. In this article, we review
recent results from simulations and theory that lead to novel insights on how
the binding equilibrium and kinetics of these proteins is affected by the
membranes and by the membrane anchoring and molecular properties of the
proteins. Simulations and theory both indicate that the binding equilibrium
constant K2D and the on- and off-rate constants of anchored receptors and
ligands in their 'two-dimensional' (2D) membrane environment strongly depend on
the membrane roughness from thermally excited shape fluctuations on nanoscales.
Recent theory corroborated by simulations provides a general relation between
K2D} and the binding constant K3D of soluble variants of the receptors and
ligands that lack the membrane anchors and are free to diffuse in three
dimensions (3D).


  Escherichia coli swims using flagella activated by rotary motors. The
direction of rotation of the motors is indirectly regulated by the binding of a
single messenger protein. The conformational spread model has been shown to
accurately describe the equilibrium properties as well as the dynamics of the
flagellar motor. In this paper we study this model from an analytic point of
view. By exploiting the separation of timescales observed in experiments, we
show how to reduce the conformational spread model to a coarse-grained,
cooperative binding model. We show that this simplified model reproduces very
well the dynamics of the motor switch.


  A substantial period of life after reproduction ends, known as
postreproductive lifespan (PRLS), is at odds with classical life history
theory. Prolonged PRLS has been confirmed in only two non-human mammals, both
odontocete cetaceans. We investigate the evidence for PRLS in a third species,
the false killer whale, Pseudorca crassidens, using a quantitative measure of
PRLS and morphological evidence from reproductive tissue. We examined specimens
from false killer whales from combined strandings (South Africa, 1981) and
harvest (Japan 1979-80) and found morphological evidence of changes in the
activity of the ovaries in relation to age. Ovulation had ceased in 50% of
whales over 45 years, and all whales over 55 years old had ovaries classified
as postreproductive. We also calculated a measure of PRLS, known as
postreproductive representation (PrR) as an indication of the effect of
inter-population demographic variability. PrR for the combined sample was 0.14,
whereas the mean of the simulated distribution for PrR under the null
hypothesis of no PRLS was 0.02. The 99th percentile of the simulated
distribution was 0.08 and no simulated value exceeded 0.13. These results
suggest that PrR was convincingly different from the measures simulated under
the null hypothesis. We found morphological and statistical evidence for PRLS
in South African and Japanese pods of false killer whales, suggesting that this
species is the third non-human mammal in which this phenomenon has been
demonstrated in wild populations. Our estimates for PrR in false killer whales
(0.12-0.37) spanned the single values available for the short-finned pilot
whale (0.28) and the killer whale (0.22) and are comparable to estimates for
historical or hunter-gather human populations (0.3-0.47).


  Circular RNAs (circRNAs) are a class of RNA with an important role in micro
RNA (miRNA) regulation recently discovered in Human and various other
eukaryotes as well as in archaea. Here, we have analyzed RNA-seq data obtained
from {\it Enterococcus faecalis} and {\it Escherichia coli} in a way similar to
previous studies performed on eukaryotes. We report observations of circRNAs in
RNA-seq data that are reproducible across multiple experiments performed with
different protocols or growth conditions.


  Incidence of whooping cough, an infection caused by Bordetella pertussis and
Bordetella parapertussis, has been on the rise since the 1980s in many
countries. Immunological interactions, such as immune boosting and
cross-immunity between pathogens, have been hypothesised to be important
drivers of epidemiological dynamics. We present a two-pathogen model of
transmission which examines how immune boosting and cross-immunity can
influence the timing and severity of epidemics. We use a combination of
numerical simulations and bifurcation techniques to study the dynamical
properties of the system, particularly the conditions under which stable
periodic solutions are present. We derive analytic expressions for the steady
state of the single-pathogen model, and give a condition for the presence of
periodic solutions. A key result from our two-pathogen model is that, while
studies have shown that immune boosting at relatively strong levels can
independently generate periodic solutions, cross-immunity allows for the
presence of periodic solutions even when the level of immune boosting is weak.
Asymmetric cross-immunity can produce striking increases in the incidence and
period. Our study underscores the importance of developing a better
understanding of the immunological interactions between pathogens in order to
improve model-based interpretations of epidemiological data.


  One of the most important problems in development is how epigenetic domains
can be first established, and then maintained, within cells. To address this
question, we propose a framework which couples 3D chromatin folding dynamics,
to a "recolouring" process modelling the writing of epigenetic marks. Because
many intra-chromatin interactions are mediated by bridging proteins, we
consider a "two-state" model with self-attractive interactions between two
epigenetic marks which are alike (either active or inactive). This model
displays a first-order-like transition between a swollen, epigenetically
disordered, phase, and a compact, epigenetically coherent, chromatin globule.
If the self-attraction strength exceeds a threshold, the chromatin dynamics
becomes glassy, and the corresponding interaction network freezes. By modifying
the epigenetic read-write process according to more biologically-inspired
assumptions, our polymer model with recolouring recapitulates the
ultrasensitive response of epigenetic switches to perturbations, and accounts
for multi-domain conformations, strikingly similar to the
topologically-associating-domains observed in eukaryotic chromosomes.


  Previous research shows that women are more altruist than men in dictator
game experiments. Yet, little is known whether women are expected to be more
altruist than men. Here we elicit third-parties' beliefs about dictators'
donations conditional on knowing the gender of the dictator. Our data provide
evidence of three main findings: (i) women are expected to be more altruist
than men; (ii) both men and women have correct beliefs about the level of
altruism among men; and (iii) both men and women overestimate the level of
altruism among women. In doing so, our results uncover a perception gap
according to which, although women are more altruist than men, they are
expected to be even more altruist than they actually are.


  The amplitude of fluctuation-induced patterns might be expected to be
proportional to the strength of the driving noise, suggesting that such
patterns would be difficult to observe in nature. Here, we show that a large
class of spatially-extended dynamical systems driven by intrinsic noise can
exhibit giant amplification, yielding patterns whose amplitude is comparable to
that of deterministic Turing instabilities. The giant amplification results
from the interplay between noise and non-orthogonal eigenvectors of the linear
stability matrix, yielding transients that grow with time, and which, when
driven by the ever-present intrinsic noise, lead to persistent large amplitude
patterns. This mechanism provides a robust basis for fluctuation-induced
biological pattern formation based on the Turing mechanism, without requiring
fine tuning of diffusion constants.


  In the present work we study the relationship between population allocation
and the combined effects of urban size and energy consumption, for two given
areas and through a major part of the twentieth century. Along these lines a
general application model is laid down which relates city-growth rates to
initial inhabitants and to exosomatic energy increment, the deviations from it
showing order in space and time; as shown in a series of maps which hint at
unaccounted socioeconomic factors. The study of the maps by means of spectral
analysis allows finding patterns which reinforce over time, in such a manner
that spatial frequencies can be determined whose weight increases up so
granting surface evolution estimation.


  The nonlinear recombination equation from population genetics has a long
history and is notoriously difficult to solve, both in continuous and in
discrete time. This is particularly so if one aims at full generality, thus
also including degenerate parameter cases. Due to recent progress for the
continuous time case via the identification of an underlying stochastic
fragmentation process, it became clear that a direct general solution at the
level of the corresponding ODE itself should also be possible. This paper shows
how to do it, and how to extend the approach to the discrete-time case as well.


  Proteins are an example of heteropolymers able to self-assemble in specific
target structures. The self-assembly of designed artificial heteropolymers is
still, to the best of our knowledge, not possible with control over the single
chain self-assembling properties comparable to what natural proteins can
achieve. What artificial heteropolymers lacks compared to bio-heteropolymers
that grants the latter such a versatility? Is the geometry of the protein
skeleton the only a particular choice to be designable? Here we introduce a
general criteria to discriminate which polymer backbones can be designed to
adopt a predetermined structure. With our approach we can explore different
polymer backbones and different amino acids alphabets. By comparing the radial
distribution functions of designable and not-designable scenarios we identify
as designability criteria the presence of a particular peak in the radial
distribution function that dominates over the random packing of the
heteropolymer. We show that the peak is a universal feature of all designable
heteropolymers, as it is dominating also the radial distribution function of
natural proteins. Our finding can help in understanding the key features that
make proteins a highly designable system. The criteria that we present can be
applied to engineer new types of self-assembling modular polymers that will
open new applications for polymer based material science.


  Next Generation Sequencing (NGS), a recently evolved technology, have served
a lot in the research and development sector of our society. This novel
approach is a newbie and has critical advantages over the traditional Capillary
Electrophoresis (CE) based Sanger Sequencing. The advancement of NGS has led to
numerous important discoveries, which could have been costlier and time taking
in case of traditional CE based Sanger sequencing. NGS methods are highly
parallelized enabling to sequence thousands to millions of molecules
simultaneously. This technology results into huge amount of data, which need to
be analysed to conclude valuable information. Specific data analysis algorithms
are written for specific task to be performed. The algorithms in group, act as
a tool in analysing the NGS data. Analysis of NGS data unravels important clues
in quest for the treatment of various life-threatening diseases; improved crop
varieties and other related scientific problems related to human welfare. In
this review, an effort was made to address basic background of NGS
technologies, possible applications, computational approaches and tools
involved in NGS data analysis, future opportunities and challenges in the area.


  Though the problem of sequence-reversed protein folding is largely
unexplored, one might speculate that reversed native protein sequences should
be significantly more foldable than purely random heteropolymer sequences. In
this article, we investigate how the reverse-sequences of native proteins might
fold by examining a series of small proteins of increasing structural
complexity ({\alpha}-helix, \b{eta}-hairpin, {\alpha}-helix bundle, and
{\alpha}/\b{eta}-protein). Employing a tandem protein structure prediction
algorithmic and molecular dynamics simulation approach, we find that the
ability of reverse sequences to adopt native-like folds is strongly in
influenced by protein size and the flexibility of the native hydrophobic core.
For \b{eta}-hairpins with reverse-sequences that fail to fold, we employ a
simple mutational strategy for guiding stable hairpin formation that involves
the insertion of amino acids into the \b{eta}-turn region. This systematic look
at reverse sequence duality sheds new light on the problem of protein
sequence-structure mapping and may serve to inspire new protein design and
protein structure prediction protocols.


  In gene expression, various kinds of proteins need to bind to specific locus
of DNA. It is still not clear how these proteins find their target locus. In
this study, the mean first-passage time (FPT) of protein binding to its target
locus on DNA chain is discussed by a chain-space coupled model. Our results
show that the 1-dimensional diffusion constant has a critical value, with which
the mean time spent by a protein to find its target locus is almost independent
of the binding rate of protein to DNA chain and the detachment rate from DNA
chain. Which implies that, the frequency of protein binding to DNA and the
sliding time on DNA chain have little influence on the search efficiency, and
therefore whether or not the 1-dimensional sliding on DNA chain increases the
search efficiency depends on the 1-dimensional diffusion constant of the
protein on DNA chain. This study also finds that only protein bindings to DNA
loci which are close to the target locus help to increase the search
efficiency, while bindings to those loci which are far from the target locus
might delay the target binding process. As expected, the mean FPT increases
with the distance between the initial position of protein in cell space and its
target locus on DNA chain. The direct binding probability, which can be
regarded as one index to describe if the 1-dimensional sliding along DNA chain
is helpful to increase the search efficiency is calculated. Our results show
that the influence of 1-dimensional sliding along DNA chain on the search
process depends on both diffusion constants of protein in cell space and on the
1-dimensional DNA chain.


  The runtime of evolutionary algorithms (EAs) depends critically on their
parameter settings, which are often problem-specific. Automated schemes for
parameter tuning have been developed to alleviate the high costs of manual
parameter tuning. Experimental results indicate that self-adaptation, where
parameter settings are encoded in the genomes of individuals, can be effective
in continuous optimisation. However, results in discrete optimisation have been
less conclusive. Furthermore, a rigorous runtime analysis that explains how
self-adaptation can lead to asymptotic speedups has been missing. This paper
provides the first such analysis for discrete, population-based EAs. We apply
level-based analysis to show how a self-adaptive EA is capable of fine-tuning
its mutation rate, leading to exponential speedups over EAs using fixed
mutation rates.


  Emergence of new protein structures has proved difficult to trace in nature
and engineer in the laboratory. However, one aspect of structure evolution has
proved immensely helpful for determining the three-dimensional structure of
proteins from their sequences: in the vast majority of cases, proteins that
share more than 30% sequence identity have similar structures. Below this mark
is the "twilight zone" where proteins may have identical or very different
structures. These observations form the foundational intuition behind structure
homology modeling. Despite their importance, however, they have never received
a comprehensive biophysical justification. Here we show that the onset of the
twilight zone is more gradual for proteins with low contact density, a proxy
for low thermodynamic stability, than proteins with high contact density. Then
we present an analytical model that treats divergent fold evolution as an
activated process, in analogy to chemical kinetics, where sequence evolution
must overcome thermodynamically unstable evolutionary intermediates to discover
new folds. This model explains the existence of a twilight zone and explains
why its onset is more abrupt for some classes of proteins than for others. We
test the assumptions of the model and characterize the dynamics of fold
evolution using evolutionary simulations of model proteins and cell
populations. Overall these results show how fundamental biophysical constraints
directed evolutionary dynamics leading to the Universe of modern protein
structures and sequences.


  Understanding why strains with different metabolic pathways that compete for
a single limiting resource coexist is a challenging issue within a theoretical
perspective. Previous investigations rely on mechanisms such as group or
spatial structuring to achieve a stable coexistence between competing metabolic
strategies. Nevertheless, coexistence has been experimentally reported even in
situations where it cannot be attributed to spatial effects [Heredity {\bf
100}, 471 (2008)]. According to that study a toxin expelled by one of the
strains can be responsible for the stable maintenance of the two strain types.
We propose a resource-based model in which an efficient strain with a slow
metabolic rate competes with a second strain type which presents a fast but
inefficient metabolism. Moreover, the model assumes that the inefficient strain
produces a toxin as a byproduct. This toxin affects the growth rate of both
strains with different strength. Through an extensive exploration of the
parameter space we determine the situations at which the coexistence of the two
strains is possible. Interestingly, we observe that the resource influx rate
plays a key role in the maintenance of the two strain types. In a scenario of
resource scarcity the inefficient is favored, though as the resource influx
rate is augmented the coexistence becomes possible and its domain is enlarged.


  D. discoideum has been recently suggested as an example of bet-hedging. Upon
starvation a population of unicellular amoebae splits between aggregators,
which form a fruiting body made of a stalk and resistant spores, and
non-aggregators. Spores are favored by long starvation periods, but vegetative
cells can exploit resources in fast-recovering environments. This partition can
be understood as a bet-hedging strategy that evolves in response to stochastic
starvation times. A genotype is defined by a different balance between each
type of cells. In this framework, if the ecological conditions are defined in
terms of the mean starvation time (i.e. time between onset of starvation and
the arrival of a new food pulse), a single genotype dominates each environment,
which is inconsistent with the huge genetic diversity observed in nature. We
investigate whether seasonality, represented by a periodic alternation in the
mean starvation times, allows the coexistence of several strategies. We use a
non-spatial (well-mixed) setting where different strains compete for a pulse of
resources. We find that seasonality, which we model via two seasons, induces a
temporal storage effect that can promote the stable coexistence of multiple
genotypes. Two conditions need to be met. First, the distributions of
starvation times in each season cannot overlap in order to create two well
differentiated habitats within the year. Second, numerous growth-starvation
cycles have to occur during each season to allow well-adapted strains to grow
and survive the subsequent unfavorable period. Additional tradeoffs among
life-history traits can expand the range of coexistence and increase the number
of coexisting strategies, contributing towards explaining the genetic diversity
observed in D. discoideum


  Studies of cooperation in microbes often focus on one fitness component, with
little information about or attention to the ecological context, and this can
lead to paradoxical results. The life cycle of the social amoeba Dictyostelium
discoideum includes a multicellular stage in which not necessarily clonal
amoebae aggregate upon starvation to form a possibly chimeric (genetically
heterogeneous) fruiting body made of dead stalk and spores. The lab-measured
reproductive skew in the spores of chimeras indicates strong social antagonism;
this should result in low genotypic diversity, which is inconsistent with
observations from nature. Two studies have suggested that this inconsistency
stems from the one-dimensional assessment of fitness (spore production) and
that the solution lies in tradeoffs between multiple traits, e.g.: spore size
versus viability; and staying vegetative versus becoming dormant. We
theoretically explore different tradeoff-implementing mechanisms and provide a
unifying ecological framework in which the two tradeoffs above, as well as
novel ones, arise collectively in response to characteristics of the
environment. We find that spore production comes at the expense of vegetative
cell production, time to development, and, depending on the experimental setup,
spore size and viability. Furthermore, we find that all existing experimental
results regarding chimeric mixes can be qualitatively recapitulated without
needing to invoke social interactions, which allows for simple resolutions to
previously paradoxical results. We conclude that the complexities of life
histories, including social behavior and multicellularity, can only be
understood in the appropriate multidimensional ecological context.


  Time-series of allele frequencies are a useful and unique set of data to
determine the strength of natural selection on the background of genetic drift.
Technically, the selection coefficient is estimated by means of a likelihood
function built under the hypothesis that the available trajectory spans a
sufficiently large portion of the fitness landscape. Especially for ancient
DNA, however, often only one single such trajectories is available and the
coverage of the fitness landscape is very limited. In fact, one single
trajectory is more representative of a process conditioned both in the initial
and in the final condition than of a process free to visit the available
fitness landscape. Based on two models of population genetics, here we show how
to build a likelihood function for the selection coefficient that takes the
statistical peculiarity of single trajectories into account. We show that this
conditional likelihood delivers a precise estimate of the selection coefficient
also when allele frequencies are close to fixation whereas the unconditioned
likelihood fails. Finally, we discuss the fact that the traditional,
unconditioned likelihood always delivers an answer, which is often
unfalsifiable and appears reasonable also when it is not correct.


  Phylogenetic inference-the derivation of a hypothesis for the common
evolutionary history of a group of species- is an active area of research at
the intersection of biology, computer science, mathematics, and statistics. One
assumes the data contains a phylogenetic signal that will be recovered with
varying accuracy due to the quality of the method used, and the quality of the
data.
  The input for distance-based inference methods is an element of a Euclidean
space with coordinates indexed by the pairs of organisms. For several
algorithms there exists a subdivision of this space into polyhedral cones such
that inputs in the same cone return the same tree topology. The geometry of
these cones has been used to analyze the inference algorithms. In this chapter,
we model how input data points drawn from DNA sequences are distributed
throughout Euclidean space in relation to the space of tree metrics, which in
turn can also be described as a collection of polyhedral cones.


  With the resurgence of tick-borne diseases such as Lyme disease and the
emergence of new pathogens such as Powassan virus, understanding what
distinguishes vector from non-vector species, and predicting undiscovered tick
vectors is an important step towards mitigating human disease risk. We apply
generalized boosted regression to interrogate over 90 features for over 240
species of Ixodes ticks. Our model predicted vector status with ~97% accuracy
and implicated 14 tick species whose intrinsic trait profiles confer high
probabilities (~80%) that they are capable of transmitting infections from
animal hosts to humans. Distinguishing characteristics of zoonotic tick vectors
include several anatomical structures that facilitate efficient host seeking
and blood-feeding from a wide variety of host species. Boosted regression
analysis produced both actionable predictions to guide ongoing surveillance as
well as testable hypotheses about the biological underpinnings of vectorial
capacity across tick species.


  Woolly mammoths (Mammuthus primigenius) populated Siberia, Beringia, and
North America during the Pleistocene and early Holocene. Recent breakthroughs
in ancient DNA sequencing have allowed for complete genome sequencing for two
specimens of woolly mammoths (Palkopoulou et al. 2015). One mammoth specimen is
from a mainland population ~45,000 years ago when mammoths were plentiful. The
second, a 4300 yr old specimen, is derived from an isolated population on
Wrangel island where mammoths subsisted with small effective population size
more than 43-fold lower than previous populations. These extreme differences in
effective population size offer a rare opportunity to test nearly neutral
models of genome architecture evolution within a single species. Using these
previously published mammoth sequences, we identify deletions, retrogenes, and
non-functionalizing point mutations. In the Wrangel island mammoth, we identify
a greater number of deletions, a larger proportion of deletions affecting gene
sequences, a greater number of candidate retrogenes, and an increased number of
premature stop codons. This accumulation of detrimental mutations is consistent
with genomic meltdown in response to low effective population sizes in the
dwindling mammoth population on Wrangel island. In addition, we observe high
rates of loss of olfactory receptors and urinary proteins, either because these
loci are non-essential or because they were favored by divergent selective
pressures in island environments. Finally, at the locus of FOXQ1 we observe two
independent loss-of-function mutations, which would confer a satin coat
phenotype in this island woolly mammoth.


  We propose a stochastic model for gene transcription coupled to DNA
supercoiling, where we incorporate the experimental observation that
polymerases create supercoiling as they unwind the DNA helix, and that these
enzymes bind more favourably to regions where the genome is unwound. Within
this model, we show that when the transcriptionally induced flux of
supercoiling increases, there is a sharp crossover from a regime where
torsional stresses relax quickly and gene transcription is random, to one where
gene expression is highly correlated and tightly regulated by supercoiling. In
the latter regime, the model displays transcriptional bursts, waves of
supercoiling, and up-regulation of divergent or bidirectional genes. It also
predicts that topological enzymes which relax twist and writhe should provide a
pathway to down-regulate transcription. This article has been published in
Physical Review Letters, May 2016.


  The speciation model proposed by Derrida and Higgs demonstrated that a
sexually reproducing population can split into different species in the absence
of natural selection or any type of geographic isolation, provided that mating
is assortative and the number of genes involved in the process is infinite.
Here we revisit this model and simulate it for finite genomes, focusing on the
question of how many genes it actually takes to trigger neutral sympatric
speciation. We find that, for typical parameters used in the original model, it
takes of the order of $10^5$ genes. We compare the results with a similar
spatially explicit model where about 100 genes suffice for speciation. We show
that when the number of genes is small the species that emerge are strongly
segregated in space. For larger number of genes, on the other hand, the spatial
structure of the population is less important and the species distribution
overlap considerably.


  We report a new mechanism for allelic dominance in regulatory genetic
interactions that we call binding dominance. We investigated a biophysical
model of gene regulation, where the fractional occupancy of a transcription
factor (TF) on the cis-regulated promoter site it binds to is determined by
binding energy (-{\Delta}G) and TF dosage. Transcription and gene expression
proceed when the TF is bound to the promoter. In diploids, individuals may be
heterozygous at the cis-site, at the TF's coding region, or at the TF's own
promoter, which determines allele-specific dosage. We find that when the TF's
coding region is heterozygous, TF alleles compete for occupancy at the cis
sites and the tighter-binding TF is dominant in proportion to the difference in
binding strength. When the TF's own promoter is heterozygous, the TF produced
at the higher dosage is also dominant. Cis-site heterozygotes have additive
expression and therefore codominant phenotypes. Binding dominance propagates to
affect the expression of downstream loci and it is sensitive in both magnitude
and direction to genetic background, but its detectability often attenuates.
While binding dominance is inevitable at the molecular level, it is difficult
to detect in the phenotype under some biophysical conditions, more so when TF
dosage is high and allele-specific binding affinities are similar. A body of
empirical research on the biophysics of TF binding demonstrates the
plausibility of this mechanism of dominance, but studies of gene expression
under competitive binding in heterozygotes in a diversity of genetic
backgrounds are needed.


  Replicative senescence, induced by telomere shortening, exhibits considerable
asynchrony and heterogeneity, the origins of which remain unclear. Here, we
formally study how telomere shortening mechanisms impact on senescence kinetics
and define two regimes of senescence, depending on the initial telomere length
variance. We provide analytical solutions to the model, highlighting a
non-linear relationship between senescence onset and initial telomere length
distribution. This study reveals the complexity of the collective behavior of
telomeres as they shorten, leading to senescence heterogeneity.


  Many natural processes rely on optimizing the success ratio of a search
process. We use an experimental setup consisting of a simple online game in
which players have to find a target hidden on a board, to investigate the how
the rounds are influenced by the detection of cues. We focus on the search
duration and the statistics of the trajectories traced on the board. The
experimental data are explained by a family of random-walk-based models and
probabilistic analytical approximations. If no initial information is given to
the players, the search is optimized for cues that cover an intermediate
spatial scale. In addition, initial information about the extension of the cues
results, in general, in faster searches. Finally, strategies used by informed
players turn into non-stationary processes in which the length of each
displacement evolves to show a well-defined characteristic scale that is not
found in non-informed searches.


  In this paper we study properties of topological RNA structures, i.e.~RNA
contact structures with cross-serial interactions that are filtered by their
topological genus. RNA secondary structures within this framework are
topological structures having genus zero. We derive a new bivariate generating
function whose singular expansion allows us to analyze the distributions of
arcs, stacks, hairpin- , interior- and multi-loops. We then extend this
analysis to H-type pseudoknots, kissing hairpins as well as $3$-knots and
compute their respective expectation values. Finally we discuss our results and
put them into context with data obtained by uniform sampling structures of
fixed genus.


  We compare the phylogenetic tensors for various trees and networks for two,
three and four taxa. If the probability spaces between one tree or network and
another are not identical then there will be phylogenetic tensors that could
have arisen on one but not the other. We call these two trees or networks
distinguishable from each other. We show that for the binary symmetric model
there are no two-taxon trees and networks that are distinguishable from each
other, however there are three-taxon trees and networks that are
distinguishable from each other.
  We compare the time parameters for the phylogenetic tensors for various taxon
label permutations on a given tree or network. If the time parameters on one
taxon label permutation in terms of the other taxon label permutation are all
non-negative then we say that the two taxon label permutations are not network
identifiable from each other. We show that some taxon label permutations are
network identifiable from each other.
  We show that some four-taxon networks satisfy the four-point condition. Of
the two "shapes" of four-taxon rooted trees, one is defined by the cluster,
b,c,d, labelling taxa alphabetically from left to right. The network with this
shape and convergence between the two taxa with the root as their most recent
common ancestor satisfies the four-point condition.
  The phylogenetic tensors contain polynomial equations that cannot be easily
solved for four-taxon or higher trees or networks. We show how methods from
algebraic geometry, such as Gr\"obner bases, can be used to solve the
polynomial equations. We show that some four-taxon trees and networks can be
distinguished from each other.


  The question of how stiff polymers are able to pack into small containers is
particularly relevant to the study of DNA packaging in viruses. A reduced
version of the problem based on coarse-grained representations of the main
components of the system -- the DNA polymer and the spherical viral capsid --
has been studied by molecular dynamics simulation. The results, involving
longer polymers than in earlier work, show that as polymers become more rigid
there is an increasing tendency to self-organize as spools that wrap from the
inside out, rather than the inverse direction seen previously. In the final
state, a substantial part of the polymer is packed into one or more coaxial
spools, concentrically layered with different orientations, a form of packaging
achievable without twisting the polymer.


  Vascular endothelial growth factor receptors (VEGFRs) are receptor tyrosine
kinases (RTKs) that regulate proliferation, migration, angiogenesis and
vascular permeability of endothelial cells. VEGFR1 and VEGFR2 bind vascular
endothelial growth factors (VEGFs), inducing receptor dimerisation and
activation, characterised by phosphorylation of tyrosine residues in their
cytoplasmic domain. Although experimental evidence suggests that RTK signalling
occurs both on the plasma membrane and intra-cellularly, and reveals the role
of endocytosis in RTK signal transduction, we still lack knowledge of VEGFR
phosphorylation-site use and of the spatiotemporal regulation of VEGFR
signalling. In this paper, we introduce four stochastic mathematical models to
study the binding kinetics of vascular endothelial growth factor VEGF-A to
VEGFR1 and VEGFR2, and phosphorylation. The formation of phosphorylated dimers
on the cell surface is a two-step process: diffusive transport and binding. The
first two of our models only consider VEGFR2, which allows us to introduce new
stochastic descriptors making use of a matrix-analytic approach. The two
remaining models describe the competition of VEGFR1 and VEGFR2 for ligand
availability, and are analysed making use of Gillespie simulations and the van
Kampen approximation. Under the hypothesis that bound phosphorylated receptor
dimers are the signalling units, we study the time to reach a threshold number
of such complexes. Our results indicate that the presence of VEGFR1 does not
only affect the timescale to reach a given signalling threshold, but it also
affects the maximum attainable threshold. This result is consistent with the
conjectured role of VEGFR1 as a decoy receptor, that prevents VEGF-A binding to
VEGFR2, and thus, VEGFR2 attaining suitable phosphorylation levels. We identify
an optimum range of ligand concentration for sustained dimer phosphorylation.


  Phylogenetic networks model reticulate evolutionary histories. The last two
decades have seen an increased interest in establishing mathematical results
and developing computational methods for inferring and analyzing these
networks. A salient concept underlying a great majority of these developments
has been the notion that a network displays a set of trees and those trees can
be used to infer, analyze, and study the network.
  In this paper, we show that in the presence of coalescence effects, the set
of displayed trees is not sufficient to capture the network. We formally define
the set of parental trees of a network and make three contributions based on
this definition. First, we extend the notion of anomaly zone to phylogenetic
networks and report on anomaly results for different networks. Second, we
demonstrate how coalescence events could negatively affect the ability to infer
a species tree that could be augmented into the correct network. Third, we
demonstrate how a phylogenetic network can be viewed as a mixture model that
lends itself to a novel inference approach via gene tree clustering.
  Our results demonstrate the limitations of focusing on the set of trees
displayed by a network when analyzing and inferring the network. Our findings
can form the basis for achieving higher accuracy when inferring phylogenetic
networks and open up new venues for research in this area, including new
problem formulations based on the notion of a network's parental trees.


  We characterize the equation of state for a simple three-dimensional DNA
hairpin model using a Metropolis Monte Carlo algorithm. This algorithm was run
at constant temperature and fixed separation between the terminal ends of the
strand. From the equation of state, we compute the compressibility, thermal
expansion coefficient, and specific heat along with adiabatic path.


  The Ebola virus in West Africa has infected almost 30,000 and killed over
11,000 people. Recent models of Ebola Virus Disease (EVD) have often made
assumptions about how the disease spreads, such as uniform transmissibility and
homogeneous mixing within a population. In this paper, we test whether these
assumptions are necessarily correct, and offer simple solutions that may
improve disease model accuracy. First, we use data and models of West African
migration to show that EVD does not homogeneously mix, but spreads in a
predictable manner. Next, we estimate the initial growth rate of EVD within
country administrative divisions and find that it significantly decreases with
population density. Finally, we test whether EVD strains have uniform
transmissibility through a novel statistical test, and find that certain
strains appear more often than expected by chance.


  We study force generation by a set of parallel actin filaments growing
against a non-rigid obstacle, in presence of an external load. The filaments
polymerize by either moving the whole obstacle, with a large energy cost, or by
causing local distortion in its shape which costs much less energy. The
non-rigid obstacle also has local thermal fluctuations due to which its shape
can change with time and we describe this using fluctuations in the height
profile of a one dimensional interface with Kardar-Parisi-Zhang dynamics. We
find the shape fluctuations of the barrier strongly affects the force
generation mechanism. The qualitative nature of the force-velocity curve is
crucially determined by the relative time-scale of filament and barrier
dynamics. The height profile of the barrier also shows interesting variation
with the external load. Our analytical calculation within mean-field theory
shows reasonable agreement with our simulation results.


  Mathematical modeling and analysis can provide insight on the dynamics of
ecosystems which maintain biodiversity in the face of competitive and
prey-predator interactions. Of primary interests are the underlying structure
and features which stabilize diverse ecological networks. Recently Korytowski
and Smith [17] proved that a perfectly nested infection network, along with
appropriate life history trade-offs, leads to coexistence and persistence of
bacteria-phage communities in a chemostat model. In this article, we generalize
their model in order to apply it to the within-host dynamics virus and immune
response, in particular HIV and CTL (Cytotoxic T Lymphocyte) cells. Our model
can produce a diverse hierarchy of viral and immune populations, built through
sequential viral escape from dominant immune responses and rise in subdominant
immune responses, consistent with observed patterns of HIV/CTL evolution. We
find a Lyapunov function for the system which leads to rigorous
characterization of persistent viral and immune variants, and informs upon
equilibria stability and global dynamics. Results are interpreted in the
context of within-host HIV/CTL evolution and numerical simulations are
provided.


  Social interactions shape the patterns of spreading processes in a
population. Techniques such as diaries or proximity sensors allow to collect
data about encounters and to build networks of contacts between individuals.
The contact networks obtained from these different techniques are however
quantitatively different. Here, we first show how these discrepancies affect
the prediction of the epidemic risk when these data are fed to numerical models
of epidemic spread: low participation rate, under-reporting of contacts and
overestimation of contact durations in contact diaries with respect to sensor
data determine indeed important differences in the outcomes of the
corresponding simulations {with for instance an enhanced sensitivity to initial
conditions}. Most importantly, we investigate if and how information gathered
from contact diaries can be used in such simulations in order to yield an
accurate description of the epidemic risk, assuming that data from sensors
represent the ground truth. The contact networks built from contact sensors and
diaries present indeed several structural similarities: this suggests the
possibility to construct, using only the contact diary network information, a
surrogate contact network such that simulations using this surrogate network
give the same estimation of the epidemic risk as simulations using the contact
sensor network. We present and {compare} several methods to build such
surrogate data, and show that it is indeed possible to obtain a good agreement
between the outcomes of simulations using surrogate and sensor data, as long as
the contact diary information is complemented by publicly available data
describing the heterogeneity of the durations of human contacts.


  Methylation and hydroxylation of cytosines to form 5-methylcytosine (5mC) and
5-droxymethylcytosine (5hmC) belong to the most important epigenetic
modifications and their vital role in the regulation of gene expression has
been widely recognized. Recent experimental techniques allow to infer
methylation and hydroxylation levels at CpG dinucleotides but require a
sophisticated statistical analysis to achieve accurate estimates.


  Microtubule dynamic instability allows search and capture of kinetochores
during spindle formation, an important process for accurate chromosome
segregation during cell division. Recent work has found that microtubule
rotational diffusion about minus-end attachment points contributes to
kinetochore capture in fission yeast, but the relative contributions of dynamic
instability and rotational diffusion are not well understood. We have developed
a biophysical model of kinetochore capture in small fission-yeast nuclei using
hybrid Brownian dynamics/kinetic Monte Carlo simulation techniques. With this
model, we have studied the importance of dynamic instability and microtubule
rotational diffusion for kinetochore capture, both to the lateral surface of a
microtubule and at or near its end. Over a range of biologically relevant
parameters, microtubule rotational diffusion decreased capture time, but made a
relatively small contribution compared to dynamic instability. At most,
rotational diffusion reduced capture time by 25%. Our results suggest that
while microtubule rotational diffusion can speed up kinetochore capture, it is
unlikely to be the dominant physical mechanism. In addition, we found that when
microtubules undergo dynamic instability, lateral captures predominate even in
the absence of rotational diffusion. Counterintuitively, adding rotational
diffusion to a dynamic microtubule increases the probability of end-on capture.


  Stochastic epidemic models describe the dynamics of an epidemic as a disease
spreads through a population. Typically, only a fraction of cases are observed
at a set of discrete times. The absence of complete information about the time
evolution of an epidemic gives rise to a complicated latent variable problem in
which the state space size of the epidemic grows large as the population size
increases. This makes analytically integrating over the missing data infeasible
for populations of even moderate size. We present a data augmentation Markov
chain Monte Carlo (MCMC) framework for Bayesian estimation of stochastic
epidemic model parameters, in which measurements are augmented with
subject-level disease histories. In our MCMC algorithm, we propose each new
subject-level path, conditional on the data, using a time-inhomogeneous
continuous-time Markov process with rates determined by the infection histories
of other individuals. The method is general, and may be applied, with minimal
modifications, to a broad class of stochastic epidemic models. We present our
algorithm in the context of multiple stochastic epidemic models in which the
data are binomially sampled prevalence counts, and apply our method to data
from an outbreak of influenza in a British boarding school.


  Down regulation of mRNA translation is an important problem in various
bio-medical domains ranging from developing effective medicines for tumors and
for viral diseases to developing attenuated virus strains that can be used for
vaccination. Here, we study the problem of down regulation of mRNA translation
using a mathematical model called the ribosome flow model (RFM). In the RFM,
the mRNA molecule is modeled as a chain of $n$ sites. The flow of ribosomes
between consecutive sites is regulated by $n+1$ transition rates. Given a set
of feasible transition rates, that models the outcome of all possible
mutations, we consider the problem of maximally down regulating the translation
rate by altering the rates within this set of feasible rates. Under certain
conditions on the feasible set, we show that an optimal solution can be
determined efficiently. We also rigorously analyze two special cases of the
down regulation optimization problem. Our results suggest that one must focus
on the position along the mRNA molecule where the transition rate has the
strongest effect on the protein production rate. However, this rate is not
necessarily the slowest transition rate along the mRNA molecule. We discuss
some of the biological implications of these results.


  We study the relation between the age-specific prevalence, incidence and
mortality in an illness-death model consisting of the three states Healthy,
Ill, Dead. The dependency on three different time scales (age, calendar time,
disease duration) is considered. It is shown that Keiding's equation published
in 1991 is a generalisation of the solution of Brunet and Struchiner's partial
differential equation from 1999. In a special case, we propose a particularly
simple estimate of the incidence from prevalence data.


  We study DNA denaturation by integrating elasticity -- as described by the
Gaussian network model -- with bond binding energies, distinguishing between
different base-pair and stacking energies. We use exact calculation, within the
model, of the Helmholtz free-energy of any partial denaturation state, which
implies that the entropy of all formed bubbles ("loops") is accounted for.
Considering base-pair bond removal single events, the bond designated for
opening is chosen by minimizing the free-energy difference for the process,
over all remaining base-pair bonds. Despite of its great simplicity, for
several known DNA sequences our results are in accord with available
theoretical and experimental studies. Moreover, we report free-energy profiles
along the denaturation pathway, which allow to detect stable or meta-stable
partial denaturation states, composed of "bubbles", as local free-energy minima
separated by barriers. Our approach allows to study very long DNA strands with
commonly available computational power, as we demonstrate for a few random
sequences in the range 200-800 base-pairs. For the latter we also elucidate the
self-averaging property of the system. Implications for the well known
breathing dynamics of DNA are elucidated.


  Aedes aegypti is a known vector of Dengue, Chikungunya and Zika and the goal
of this study is to propose the first mathematical model to describe the
dynamic transmission of these three diseases. We present two preliminary models
that consist of the SEIR model for the human populations and an SEI model for
the vector to describe (a) the single transmission dynamics of dengue,
Chikungunya or Zika, and (b) any possible coinfection between two diseases in
the same population. In order to do that, we obtain an analytical solution of
the system of 17 and 30 coupled differential equations for each model
respectively, and later obtain the eigenvalues by analyzing the Jacobian matrix
in order to begin the development of a surveillance system to prevent the
spread of these three diseases.


  The model for cAMP-dependent synaptic plasticity relates the characterization
of a noradrenaline-stimulated adenylyl cyclase with DNA unzipping. Specific
proteins condition cAMP insertion in specific sites of the DNA structure to
direct cellular and synaptic differentiation in brain tissues, also providing
coding. Metabolic-dependent ATP binding of Mg2+, could control feedback by
inactivating AC dependent formation of cAMP. The level of cAMP and cGMP, which
could be assayed in red cells and cerebrospinal fluid, allows a clinical lab
diagnostic improvement. Also, provides a relationship of best fitting to cAMP
control by binding to DNA. The cAMP level allows the prediction of an inverse
relationship between neurodegeneration and cancer. The latter, could be
characterized by uncontrolled proliferation, whereas metabolic dominance by
stress over a long period of time, may deplete cerebral cAMP.


  Coral reefs are dynamic systems whose composition is highly influenced by
unpredictable biotic and abiotic factors. Understanding the spatial scale at
which long-term predictions of reef composition can be made will be crucial for
guiding conservation efforts. Using a 22-year time series of benthic
composition data from 20 reefs on the Kenyan and Tanzanian coast, we studied
the long-term behaviour of Bayesian vector autoregressive state-space models
for reef dynamics, incorporating among-site variability. We estimate that if
there were no among-site variability, the total long-term variability would be
approximately one third of its current value. Thus among-site variability
contributes more to long-term variability in reef composition than does
temporal variability. Individual sites are more predictable than previously
thought, and predictions based on current snapshots are informative about
long-term properties. Our approach allowed us to identify a subset of possible
climate refugia sites with high conservation value, where the long-term
probability of coral cover <= 0.1 was very low. Analytical results show that
this probability is most strongly influenced by among-site variability and by
interactions among benthic components within sites. These findings suggest that
conservation initiatives might be successful at the site scale as well as the
regional scale.


  For many types of learning, spaced training that involves repeated long
inter-trial intervals (ITIs) leads to more robust memory formation than does
massed training that involves short or no intervals. Several cognitive theories
have been proposed to explain this superiority, but only recently has data
begun to delineate the underlying cellular and molecular mechanisms of spaced
training. We review these theories and data here. Computational models of the
implicated signaling cascades have predicted that spaced training with
irregular ITIs can enhance learning. This strategy of using models to predict
optimal spaced training protocols, combined with pharmacotherapy, suggests
novel ways to rescue impaired synaptic plasticity and learning.


  Indoor skin-contact surfaces of public fitness centers may serve as
reservoirs of potential human transmission of methicillin-resistant
Staphylococcus aureus (MRSA). We found a high prevalence of multi-drug
resistant (MDR)-MRSA of CC59 lineage harboring a variety of extracellular toxin
genes from surface swab samples collected from inanimate surfaces of fitness
centers in Memphis metropolitan area, USA. Our findings underscore the role of
inanimate surfaces as potential sources of transmission of MDR-MRSA strains
with considerable genetic diversity.


  Although viral spreading processes taking place in networks are often
analyzed using Markovian models in which both the transmission and the recovery
times follow exponential distributions, empirical studies show that, in many
real scenarios, the distribution of these times are not necessarily
exponential. To overcome this limitation, we first introduce a generalized
susceptible-infected-susceptible (SIS) spreading model that allows transmission
and recovery times to follow phase-type distributions. In this context, we
derive a lower bound on the exponential decay rate towards the infection-free
equilibrium of the spreading model without relying on mean-field
approximations. Based on our results, we illustrate how the particular shape of
the transmission/recovery distribution influences the exponential rate of
convergence towards the equilibrium.


  A previous report claimed no evidence of transgenerational epigenetic
inheritance in a mouse model of in utero environmental exposure, based on the
observation that gene expression changes observed in the germ cells of G1 and
G2 male fetus were not in the same direction. A subsequent data reanalysis
however showed a statistically significant overlap between G1 and G2 genes
irrespective of direction, leading to the suggestion that, as phenotypic
variability in epigenetic transmission has been observed in several other
examples also, the above report provided evidence in favor of, not against,
transgenerational inheritance. This criticism has recently been questioned.
Here, it is shown that the questions raised are based not only on incorrect
statistical calculations but also on wrong premise that gene expression changes
do not constitute a phenotype.


  In living cells, chemical reactions are connected by sharing their products
and substrates, and form complex networks, e.g. metabolic pathways. Here we
developed a theory to predict the sensitivity, i.e. the responses of
concentrations and fluxes to perturbations of enzymes, from network structure
alone. Responses turn out to exhibit two characteristic patterns,
$localization$ and $hierarchy$. We present a general theorem connecting
sensitivity with network topology that explains these characteristic patterns.
Our results imply that network topology is an origin of biological robustness.
Finally, we suggest a strategy to determine real networks from experimental
measurements.


  Anomalous non-exponential relaxation in hydrated biomolecules is commonly
attributed to the complexity of the free-energy landscapes, similarly to
polymers and glasses. It was found recently that the hydrogen-bond breathing of
terminal DNA base pairs exhibits a slow power-law relaxation attributable to
weak Hamiltonian chaos, with parameters similar to experimental data. Here, the
relationship is studied between this motion and spectroscopic signals measured
in DNA with a small molecular photoprobe inserted into the base-pair stack. To
this end, the earlier computational approach in combination with an analytical
theory is applied to the experimental DNA fragment. It is found that the
intensity of breathing dynamics is strongly increased in the internal base
pairs that flank the photoprobe, with anomalous relaxation quantitatively close
to that in terminal base pairs. A physical mechanism is proposed to explain the
coupling between the relaxation of base-pair breathing and the experimental
response signal. It is concluded that the algebraic relaxation observed
experimentally is very likely a manifestation of weakly chaotic dynamics of
hydrogen-bond breathing in the base pairs stacked to the photoprobe, and that
the weak nanoscale chaos can represent an ubiquitous hidden source of
non-exponential relaxation in ultrafast spectroscopy.


  Oxford Nanopore MinION sequencer is currently the smallest sequencing device
available. While being able to produce very long reads (reads of up to 100~kbp
were reported), it is prone to high sequencing error rates of up to 30%. Since
most of these errors are insertions or deletions, it is very difficult to adapt
popular seed-based algorithms designed for aligning data sets with much lower
error rates.
  Base calling of MinION reads is typically done using hidden Markov models. In
this paper, we propose to represent each sequencing read by an ensemble of
sequences sampled from such a probabilistic model. This approach can improve
the sensitivity and false positive rate of seeding an alignment compared to
using a single representative base call sequence for each read.


  We introduce delays in a tuberculosis (TB) model, representing the time delay
on the diagnosis and commencement of treatment of individuals with active TB
infection. The stability of the disease free and endemic equilibriums is
investigated for any time delay. Corresponding optimal control problems, with
time delays in both state and control variables, are formulated and studied.
Although it is well-known that there is a delay between two to eight weeks
between TB infection and reaction of body's immune system to tuberculin, delays
for the active infected to be detected and treated, and delays on the treatment
of persistent latent individuals due to clinical and patient reasons, which
clearly justifies the introduction of time delays on state and control
measures, our work seems to be the first to consider such time-delays for TB
and apply time-delay optimal control to carry out the optimality analysis.


  The proper sorting of membrane components by regulated exchange between
cellular organelles is crucial to intra-cellular organisation. This process
relies on the budding and fusion of transport vesicles, and should be strongly
influenced by stochastic fluctuations considering the relatively small size of
many organelles. We identify the perfect sorting of two membrane components
initially mixed in a single compartment as a first passage process, and we show
that the mean sorting time exhibits two distinct regimes as a function of the
ratio of vesicle fusion to budding rates. Low ratio values lead to fast
sorting, but results in a broad size distribution of sorted compartments
dominated by small entities. High ratio values result in two well defined
sorted compartments but is exponentially slow. Our results suggest an optimal
balance between vesicle budding and fusion for the rapid and efficient sorting
of membrane components, and highlight the importance of stochastic effects for
the steady-state organisation of intra-cellular compartments.


  Many naturally occurring networks have a power-law degree distribution as
well as a non-zero degree correlation. Despite this, most studies analyzing the
robustness to random node-deletion and vulnerability to targeted node-deletion
have concentrated only on power-law degree distribution and ignored degree
correlation. This study looks specifically at the effect degree-correlation has
on robustness and vulnerability in scale-free networks. Our results confirm
Newman's finding that positive degree-correlation increases robustness and
decreases vulnerability. However, we found that networks with positive
degree-correlation are more vulnerable to random node-deletion than to targeted
deletion methods that utilize knowledge of initial node-degree only. Targeted
deletion sufficiently alters the topology of the network to render this method
less effective than uniform random methods unless changes in topology are
accounted for. This result indicates the importance of degree correlation in
certain network applications.


  We study ODE models of epidemic spreading with a preventive behavioral
response that is triggered by awareness of the infection. Previous studies of
such models have mostly focused on the impact of the response on the initial
growth of an outbreak and the existence and location of endemic equilibria.
Here we study the question whether this type of response is sufficient to
prevent future flare-ups from low endemic levels if awareness is assumed to
decay over time. In the ODE context, such flare-ups would translate into
sustained oscillations with significant amplitudes.
  Our results show that such oscillations are ruled out in
Susceptible-Aware-Infectious-Susceptible models with a single compartment of
aware hosts, but can occur if we consider two distinct compartments of aware
hosts who differ in their willingness to alert other susceptible hosts.


  We propose a general theory to describe the distribution of protein-folding
transition paths. We show that transition paths follow a predictable sequence
of high-free-energy transient states that are separated by free-energy
barriers. Each transient state corresponds to the assembly of one or more
discrete, cooperative units, which are determined directly from the native
structure. We show that the transition state on a folding pathway is reached
when a small number of critical contacts are formed between a specific set of
substructures, after which folding proceeds downhill in free energy. This
approach suggests a natural resolution for distinguishing parallel folding
pathways and provides a simple means to predict the rate-limiting step in a
folding reaction. Our theory identifies a common folding mechanism for proteins
with diverse native structures and establishes general principles for the
self-assembly of polymers with specific interactions.


  Motivation: Exome or targeted sequencing data exerts analytical challenge to
test single nucleotide polymorphisms (SNPs) with extremely small minor allele
frequency (MAF). Various rare variant tests were proposed to increase power by
aggregating SNPs within a fixed genomic region, such as a gene or pathway.
However, a gene could contain from several to thousands of markers, and not all
of them may be related to the phenotype. Combining functional and
non-functional SNPs in arbitrary genomic region could impair the testing power.
Results: We propose a Zoom-Focus algorithm (ZFA) to locate the optimal testing
region within a given genomic region, as a wrapper function to be applied in
conjunction with rare variant association tests. In the first Zooming step, a
given genomic region is partitioned by order of two, and the best partition is
located within all partition levels. In the next Focusing step, boundaries of
the zoomed region are refined. Simulation studies showed that ZFA substantially
enhanced the statistical power of rare variant tests by over 10 folds,
including the WSS, SKAT and W-test. The algorithm is applied on real exome
sequencing data of hypertensive disorder, and identified biologically relevant
genetic markers to metabolic disorder that are undiscoverable by testing using
full gene. The proposed algorithm is an efficient and powerful tool to increase
the effectiveness of rare variant association tests for exome sequencing
datasets of complex disorder.