  The epidemic spreading has been widely studied when each node may get
infected by an infected neighbor with the same rate. However, the infection
rate between a pair of nodes is usually heterogeneous and even correlated with
their nodal degrees in the contact network. We aim to understand how such
correlated heterogeneous infection rates influence the spreading on different
network topologies. Motivated by real-world datasets, we propose a correlated
heterogeneous Susceptible-Infected-Susceptible model which assumes that the
infection rate $\beta_{ij}(=\beta_{ji})$ between node $i$ and $j$ is correlated
with the degree of the two end nodes: $\beta_{ij}=c(d_id_j)^\alpha$, where
$\alpha$ indicates the strength of the correlation and $c$ is selected so that
the average infection rate is $1$. In order to understand the effect of such
correlation on epidemic spreading, we consider as well the corresponding
uncorrected but still heterogeneous infection rate scenario as a reference,
where the original correlated infection rates in our CSIS model are shuffled
and reallocated to the links of the same network topology. We compare these two
scenarios in the average fraction of infected nodes in the metastable state on
Erd{\"o}s-R{\'e}nyi (ER) and scale-free (SF) networks with a similar average
degree. Through the continuous-time simulations, we find that, when the
recovery rate is small, the negative correlation is more likely to help the
epidemic spread and the positive correlation prohibit the spreading; as the
recovery rate increases to be larger than a critical value, the positive but
not negative correlation tends to help the spreading. Our findings are further
analytically proved in a wheel network (one central node connects with each of
the nodes in a ring) and validated on real-world networks with correlated
heterogeneous interaction frequencies.


  The nucleotides guanosine tetraphosphate (ppGpp) and guanosine pentaphosphate
(pppGpp) bind to target proteins to promote bacterial survival (Corrigan et al.
2016). Thus, the binding of the nucleotides to RsgA, a GTPase, inhibits the
hydrolysis of GTP. The dose response, taken to be curvilinear with respect to
the logarithm of the inhibitor concentration, is instead much better (P<0.001
when the 6 experiments are combined) represented as multiphasic, with high to
exceedingly high absolute r values for the straight lines, and with transitions
in the form of non-contiguities (jumps). Profiles for the binding of
radiolabeled nucleotides to HprT and Gmk, GTP synthesis enzymes, were,
similarly, taken to be curvilinear with respect to the logarithm of the protein
concentration. However, the profiles are again much better represented as
multiphasic than as curvilinear (the P values range from 0.047 to <0.001 for
each of the 8 experiments for binding of ppGpp and pppGpp to HprT). The binding
of GTP to HprT and the binding of the three nucleotides to Gmk are also poorly
represented by curvilinear profiles, but well represented by multiphasic
profiles (straight and, in part, parallel lines).


  Catalytically active enzymes have recently been observed to exhibit enhanced
diffusion. In a recent work [C. Riedel et al., Nature 517, 227 (2015)], it has
been suggested that this phenomenon is correlated with the degree of
exothermicity of the reaction, and a mechanism was proposed to explain the
phenomenon based on channeling the released heat into the center of mass
kinetic energy of the enzyme. I addressed this question by comparing four
different mechanisms, and concluded that collective heating is the strongest
candidate out of those four to explain the phenomenon, and in particular,
several orders of magnitude stronger than the mechanism proposed by Riedel et
al. In a recent preprint (arXiv:1608.05433), K. Tsekouras, C. Riedel, R.
Gabizon, S. Marqusee, S. Presse, and C. Bustamante present a comment on my
paper [R. Golestanian, Phys. Rev. Lett. 115, 108102 (2015); arXiv:1508.03219],
which I address here in this reply.


  The living cell uses a variety of molecular receptors to read and process
chemical signals that vary in space and time. We model the dynamics of such
molecular level measurements as Markov processes in steady state, with a
coupling between the receptor and the signal. We prove exactly that, when the
the signal dynamics is not perturbed by the receptors, the free energy consumed
by the measurement process is lower bounded by a quantity proportional to the
mutual information. Our result is completely independent of the receptor
architecture and dependent on signal properties alone, and therefore holds as a
general principle for molecular information processing.


  The recent rapid growth of social media and online social networks (OSNs) has
raised interesting questions about the spread of ideas and fads within our
society. In the past year, several papers have drawn analogies between the rise
and fall in popularity of OSNs and mathematical models used to study infectious
disease. One such model, the irSIR model, made use of the idea of "infectious
recovery" to outperform the traditional SIR model in replicating the rise and
fall of MySpace and to predict a rapid drop in the popularity of Facebook. Here
we explore the irSIR model and two of its logical extensions and we
mathematically characterize the initial and long-run behavior of these
dynamical systems. In particular, while the original irSIR model always
predicts extinction of a social epidemic, we construct an extension of the
model that matches the exponential growth phase of the irSIR model while
allowing for the possibility of an arbitrary proportion of infections in the
long run.


  DNA replication is a process which is common to all domains of life yet
different replication mechanisms are seen among different organisms. The
mechanism of replication on such a structure is not yet understood. With this
bigger picture in mind, we have initiated the characterization involved in DNA
replication.


  Ultrasensitive response motifs, which are capable of converting graded
stimulus in binary responses, are very well-conserved in signal transduction
networks. Although it has been shown that a cascade arrangement of multiple
ultrasensitive modules can produce an enhancement of the system's
ultrasensitivity, how the combination of layers affects the cascade's
ultrasensitivity remains an open question for the general case. Here we
introduced a methodology that allowed us to determine the presence of
sequestration effects and to quantify the relative contribution of each module
to the overall cascade's ultrasensitivity. The proposed analysis framework
provides a natural link between global and local ultrasensitivity descriptors
and is particularly well-suited to characterize and better understand
mathematical models used to study real biological systems. As a case study we
considered three mathematical models introduced by O'Shaughnessy et al. to
study a tunable synthetic MAPK cascade, and showed how our methodology might
help modelers to better understand modeling alternatives.


  In this paper we consider dipole-mediated correlations between DNA and
enzymes in the context of their water environment. Such correlations emerge
from electric dipole-dipole interactions between aromatic ring structures in
DNA and in enzymes. We show that there are matching collective modes between
DNA and enzyme dipole fields, and that a dynamic time-averaged polarization
vanishes in the water dipole field only if either DNA, enzyme, or both are
absent from the sample. This persistent field may serve as the electromagnetic
image that, in popular colloquialisms about DNA biochemistry, allows enzymes to
"scan" or "read" the double helix. Topologically nontrivial configurations in
the coherent ground state requiring clamplike enzyme behavior on the DNA may
stem, ultimately, from spontaneously broken gauge symmetries.


  Standard neutral population genetics theory with a strictly fixed population
size has important limitations. An alternative model that allows independently
fluctuating population sizes and reproduces the standard neutral evolution is
reviewed. We then study a situation such that the competing species are neutral
at the equilibrium population size but population size fluctuations
nevertheless favor fixation of one species over the other. In this case, a
separation of timescales emerges naturally and allows adiabatic elimination of
a fast population size variable to deduce the fluctuations-induced selection
dynamics near the equilibrium population size. The results highlight the
incompleteness of the standard population genetics with a strictly fixed
population size.


  Understanding the evolution of dispersal is an important issue in
evolutionary ecology. For continuous time models in which individuals disperse
throughout their lifetime, it has been shown that a balanced dispersal
strategy, which results in an ideal free distribution, is evolutionary stable
in spatially varying but temporally constant environments. Many species,
however, primarily disperse prior to reproduction (natal dispersal) and less
commonly between reproductive events (breeding dispersal). As demographic and
dispersal terms combine in a multiplicative way for models of natal dispersal,
rather than the additive way for the previously studied models, we develop new
mathematical methods to study the evolution of natal dispersal for
continuous-time and discrete-time models. A fundamental ecological dichotomy is
identified for the non-trivial equilibrium of these models: (i) the per-capita
growth rates for individuals in all patches is equal to zero, or (ii)
individuals in some patches experience negative per-capita growth rates, while
individuals in other patches experience positive per-capita growth rates. The
first possibility corresponds to an ideal-free distribution, while the second
possibility corresponds to a "source-sink" spatial structure. We prove that
populations with a dispersal strategy leading to an ideal-free distribution
displace populations with dispersal strategy leading to a source-sink spatial
structure. When there are patches which can not sustain a population,
ideal-free strategies can be achieved by sedentary populations, and we show
that these populations can displace populations with any irreducible dispersal
strategy. Collectively, these results support that evolution selects for natal
or breeding dispersal strategies which lead to ideal-free distributions in
spatially heterogenous, but temporally homogenous, environments.


  A simple statistical mechanical model proposed by Wako and Sait$\hat{\rm o}$
has explained the aspects of protein folding surprisingly well. This model was
systematically applied to multiple proteins by Mu$\tilde{\rm n}$oz and Eaton
and has since been referred to as the Wako-Sait$\hat{\rm o}$-Mu$\tilde{\rm
n}$oz-Eaton (WSME) model. The success of the WSME model in explaining the
folding of many proteins has verified the hypothesis that the folding is
dominated by native interactions, which makes the energy landscape globally
biased toward native conformation. Using the WSME and other related models,
Sait$\hat{\rm o}$ emphasized the importance of the hierarchical pathway in
protein folding; folding starts with the creation of contiguous segments having
a native-like configuration and proceeds as growth and coalescence of these
segments. The $\phi$-values calculated for barnase with the WSME model
suggested that segments contributing to the folding nucleus are similar to the
structural modules defined by the pattern of native atomic contacts. The WSME
model was extended to explain folding of multi-domain proteins having a complex
topology, which opened the way to comprehensively understanding the folding
process of multi-domain proteins. The WSME model was also extended to describe
allosteric transitions, indicating that the allosteric structural movement does
not occur as a deterministic sequential change between two conformations but as
a stochastic diffusive motion over the dynamically changing energy landscape.
Statistical mechanical viewpoint on folding, as highlighted by the WSME model,
has been renovated in the context of modern methods and ideas, and will
continue to provide insights on equilibrium and dynamical features of proteins.


  Species or population that proliferate faster than others become dominant in
numbers. Catalysis allows catalytic sets within a molecular reaction network to
dominate the non catalytic parts of the network by processing most of the
available substrate. As a consequence one may consider a 'catalytic fitness' of
sets of molecular species. The fittest sets emerge as the expressed chemical
backbone or sub-network of larger chemical reaction networks employed by
organisms. However, catalytic fitness depends on the systemic context and the
stability of systemic dynamics. Unstable reaction networks would easily be
reshaped or destroyed by fluctuations of the chemical environment. In this
paper we therefore focus on recognizing systemic stability as an evolutionary
selection criterion. In fact, instabilities of regulatory systems dynamics
become predictive for associated evolutionary forces driving the emergence
large reaction networks that avoid or control inherent instabilities. Systemic
instabilities can be identified and analyzed using relatively simple
mathematical random networks models of complex regulatory systems. Using a
statistical ensemble approach one can identify fundamental causes of instable
dynamics, infer evolutionary preferred network properties, and predict
evolutionary emergent control mechanisms and their entanglement with cell
differentiation processes. Surprisingly, what systemic stability tells us here
is that cells (or other non-linear regulatory systems) never had to learn how
to differentiate, but rather how to avoid and control differentiation. For
example, in this framework we can predict that regulatory systems will
evolutionary favor networks where the number of catalytic enhancers is not
larger than the number of suppressors.


  We study the emergence of cell differentiation under the assumption of the
existence of a given number of tradeoffs between genes encoding different
functions. In the model the viability of colonies is determined by the
capability of their lower level units to perform different functions, which is
implicitly determined by external chemical stimuli. Due to the existence of
tradeoffs it can be evolutionarily advantageous to evolve the division of labor
whereby the cells can suppress their contributions to some of the activities
through the activation of regulatory genes, which in its turn inflicts a cost
in terms of fitness. Our simulation results show that cell differentiation is
more likely as the number of tradeoffs is increased but the outcome also
depends on their strength. We observe the existence of critical values for the
minimum number of tradeoffs and their strength beyond that maximum cell
differentiation can be attained. Remarkably, we observe the occurrence of a
maximum tradeoff strength beyond that the population is no longer viable
imposing an upper tolerable level of constraint at the system. This tolerance
is reduced as the number of tradeoffs grows.


  Spreading processes are often modelled as a stochastic dynamics occurring on
top of a given network with edge weights corresponding to the transmission
probabilities. Knowledge of veracious transmission probabilities is essential
for prediction, optimization, and control of diffusion dynamics. Unfortunately,
in most cases the transmission rates are unknown and need to be reconstructed
from the spreading data. Moreover, in realistic settings it is impossible to
monitor the state of each node at every time, and thus the data is highly
incomplete. We introduce an efficient dynamic message-passing algorithm, which
is able to reconstruct parameters of the spreading model given only partial
information on the activation times of nodes in the network. The method is
generalizable to a large class of dynamic models, as well to the case of
temporal graphs.


  We study the extinction of epidemics in a generalized
susceptible-infected-susceptible model, where a susceptible individual becomes
infected with the rate $\lambda$ when contacting $m$ infective individual(s)
simultaneously, and an infected individual spontaneously recovers with the rate
$\mu$. By employing the Wentzel-Kramers-Brillouin approximation for the master
equation, the problem is reduced to finding the zero-energy trajectories in an
effective Hamiltonian system, and the mean extinction time $\langle T\rangle$
depends exponentially on the associated action $\mathcal {S}$ and the size of
the population $N$, $\langle T\rangle \sim \exp(N\mathcal {S})$. Because of
qualitatively different bifurcation features for $m=1$ and $m\geq2$, we derive
independently the expressions of $\mathcal {S}$ as a function of the rescaled
infection rate $\lambda/\mu$. For the weak infection, $\mathcal {S}$ scales to
the distance to the bifurcation with an exponent $2$ for $m=1$ and $3/2$ for
$m\geq2$. Finally, a rare-event simulation method is used to validate the
theory.


  The amount of completely sequenced chloroplast genomes increases rapidly
every day, leading to the possibility to build large-scale phylogenetic trees
of plant species. Considering a subset of close plant species defined according
to their chloroplasts, the phylogenetic tree that can be inferred by their core
genes is not necessarily well supported, due to the possible occurrence of
problematic genes (i.e., homoplasy, incomplete lineage sorting, horizontal gene
transfers, etc.) which may blur the phylogenetic signal. However, a trustworthy
phylogenetic tree can still be obtained provided such a number of blurring
genes is reduced. The problem is thus to determine the largest subset of core
genes that produces the best-supported tree. To discard problematic genes and
due to the overwhelming number of possible combinations, this article focuses
on how to extract the largest subset of sequences in order to obtain the most
supported species tree. Due to computational complexity, a distributed Binary
Particle Swarm Optimization (BPSO) is proposed in sequential and distributed
fashions. Obtained results from both versions of the BPSO are compared with
those computed using an hybrid approach embedding both genetic algorithms and
statistical tests. The proposal has been applied to different cases of plant
families, leading to encouraging results for these families.


  We consider a two-dimensional model of three species in rock-paper-scissors
competition and study the self-organisation of the population into fascinating
spiraling patterns. Within our individual-based metapopulation formulation, the
population composition changes due to cyclic dominance (dominance-removal and
dominance-replacement), mutations, and pair-exchange of neighboring
individuals. Here, we study the influence of mobility on the emerging patterns
and investigate when the pair-exchange rate is responsible for spiral waves to
become elusive in stochastic lattice simulations. In particular, we show that
the spiral waves predicted by the system's deterministic partial equations are
found in lattice simulations only within a finite range of the mobility rate.
We also report that in the absence of mutations and dominance-replacement, the
resulting spiraling patterns are subject to convective instability and
far-field breakup at low mobility rate. Possible applications of these
resolution and far-field breakup phenomena are discussed.


  The eco-evolutionary dynamics of species are fundamentally linked to the
energetic constraints of its constituent individuals. Of particular importance
is the interplay between reproduction and the dynamics of starvation and
recovery. To elucidate this interplay, we introduce a nutritional
state-structured model that incorporates two classes of consumer: nutritionally
replete, reproducing consumers, and undernourished, non-reproducing consumers.
We obtain strong constraints on starvation and recovery rates by deriving
allometric scaling relationships and find that population dynamics are
typically driven to a steady state. Moreover, these rates fall within a
'refuge' in parameter space, where the probability of population extinction is
minimized. We also show that our model provides a natural framework to predict
maximum mammalian body size by determining the relative stability of an
otherwise homogeneous population to a competing population with altered percent
body fat. This framework provides a principled mechanism for a selective driver
of Cope's rule.


  Among all the proteins of Periplasmic C type cytochrome A (PPCA) family
obtained from cytochrome C7 found in Geobacter sulfurreducens, PPCA protein can
interact with Deoxycholate (DXCA), while its other homologs do not, as observed
from the crystal structures. Utilizing the concept of 'structure-function
relationship', an effort has been initiated towards understanding the driving
force for recognition of DXCA exclusively by PPCA among its homologs. Further,
a combinatorial analysis of the binding sequences (contiguous sequence of amino
acid residues of binding locations) is performed to build graph-theoretic
models, which show that PPCA differs from its homologues. Analysis of the
results suggests that the underlying impetus of recognition of DXCA by PPCA is
embedded in its primary sequence and 3D conformation.


  Historically, infectious diseases caused considerable damage to human
societies, and they continue to do so today. To help reduce their impact,
mathematical models of disease transmission have been studied to help
understand disease dynamics and inform prevention strategies. Vaccination - one
of the most important preventive measures of modern times - is of great
interest both theoretically and empirically. And in contrast to traditional
approaches, recent research increasingly explores the pivotal implications of
individual behavior and heterogeneous contact patterns in populations. Our
report reviews the developmental arc of theoretical epidemiology with emphasis
on vaccination, as it led from classical models assuming homogeneously mixing
(mean-field) populations and ignoring human behavior, to recent models that
account for behavioral feedback and/or population spatial/social structure.
Many of the methods used originated in statistical physics, such as lattice and
network models, and their associated analytical frameworks. Similarly, the
feedback loop between vaccinating behavior and disease propagation forms a
coupled nonlinear system with analogs in physics. We also review the new
paradigm of digital epidemiology, wherein sources of digital data such as
online social media are mined for high-resolution information on
epidemiologically relevant individual behavior. Armed with the tools and
concepts of statistical physics, and further assisted by new sources of digital
data, models that capture nonlinear interactions between behavior and disease
dynamics offer a novel way of modeling real-world phenomena, and can help
improve health outcomes. We conclude the review by discussing open problems in
the field and promising directions for future research.


  The description of the dynamics of a complex, high-dimensional system in
terms of a low-dimensional set of collective variables Y can be fruitful if the
low dimensional representation satisfies a Langevin equation with drift and
diffusion coefficients which depend only on Y. We present a computational
scheme to evaluate whether a given collective variable provides a faithful
low-dimensional representation of the dynamics of a high-dimensional system.
The scheme is based on the framework of finite-difference Langevin-equation,
similar to that used for molecular-dynamics simulations. This allows one to
calculate the drift and diffusion coefficients in any point of the
full-dimensional system. The width of the distribution of drift and diffusion
coefficients in an ensemble of microscopic points at the same value of Y
indicates to which extent the dynamics of Y is described by a simple Langevin
equation. Using a simple protein model we show that collective variables often
used to describe biopolymers display a non-negligible width both in the drift
and in the diffusion coefficients. We also show that the associated effective
force is compatible with the equilibrium free--energy calculated from a
microscopic sampling, but results in markedly different dynamical properties.


  In this article we build a mathematical model for forest growth and we
compare this model with a computer forest simulator named SORTIE. The main
ingredient taken into account in both models is the competition for light
between trees. The parameters of the mathematical model are estimated by using
SORTIE model, when the parameter values of SORTIE model correspond to the ones
previously evaluated for the Great Mountain Forest in USA. We construct a size
structured population dynamics model with one and two species and with spatial
structure.


  Accessing the network through which a propagation dynamics diffuse is
essential for understanding and controlling it. In a few cases, such
information is available through direct experiments or thanks to the very
nature of propagation data. In a majority of cases however, available
information about the network is indirect and comes from partial observations
of the dynamics, rendering the network reconstruction a fundamental inverse
problem. Here we show that it is possible to reconstruct the whole structure of
an interaction network and to simultaneously infer the complete time course of
activation spreading, relying just on single epoch (i.e. snapshot) or
time-scattered observations of a small number of activity cascades. The method
that we present is built on a Belief Propagation approximation, that has shown
impressive accuracy in a wide variety of relevant cases, and is able to infer
interactions in presence of incomplete time-series data by providing a detailed
modeling of the posterior distribution of trajectories conditioned to the
observations. Furthermore, we show by experiments that the information content
of full cascades is relatively smaller than that of sparse observations or
single snapshots.


  Network theorists have developed methods to characterize the complex
interactions in natural phenomena. The structure of the network of interactions
between proteins is important in the field of proteomics, and has been subject
to intensive research in recent years, as scientists have become increasingly
capable and interested in describing the underlying structure of interactions
in both normal and pathological biological processes. In this paper, we survey
the graph-theoretic characterization of protein-protein interaction networks
(PINs) in terms of structural features, and discuss its possible applications
in biomedical research. We also perform a brief revision of network theory's
classical literature and discuss modern statistical and computational
techniques to describe the structure of PINs


  Environmental fluctuations have important consequences in the organization of
ecological communities, and understanding how such a variability influences the
biodiversity of an ecosystem is a major question in ecology. In this paper, we
analyze the case of two species competing for the resources within the
framework of the neutral theory in the presence of environmental noise,
devoting special attention on how such a variability modulates species fitness.
The environment is dichotomous and stochastically alternates between periods
favoring one of the species while disfavoring the other one, preserving
neutrality on the long term. We study two different scenarios: in the first one
species fitness varies linearly with the environment, and in the second one the
effective fitness is re-scaled by the total fitness of the individuals
competing for the same resource. We find that, in the former case environmental
fluctuations always reduce the time of species coexistence, whereas such a time
can be enhanced or reduced in the latter case, depending on the correlation
time of the environment. This phenomenon can be understood as a direct
consequence of Chesson's storage effect.


  The cytoskeleton is an inhomogeneous network of semi-flexible filaments,
which are involved in a wide variety of active biological processes. Although
the cytoskeletal filaments can be very stiff and embedded in a dense and
cross-linked network, it has been shown that, in cells, they typically exhibit
significant bending on all length scales. In this work we propose a model of a
semi-flexible filament deformed by different types of cross-linkers for which
one can compute and investigate the bending spectrum. Our model allows to
couple the evolution of the deformation of the semi-flexible polymer with the
stochastic dynamics of linkers which exert transversal forces onto the
filament. We observe a $q^{-2}$ dependence of the bending spectrum for some
biologically relevant parameters and in a certain range of wavenumbers $q$.
However, generically, the spatially localized forcing and the non-thermal
dynamics both introduce deviations from the thermal-like $q^{-2}$ spectrum.


  Bone remodeling is a complex process involving cell-cell interactions,
biochemical signaling and mechanical stimuli. Early models of the biological
aspects of remodeling were non-spatial and focused on the local dynamics at a
fixed location in the bone. Several spatial extensions of these models have
been proposed, but they generally suffer from two limitations: first, they are
not amenable to analysis and are computationally expensive, and second, they
neglect the role played by bone-embedded osteocytes. To address these issues,
we developed a novel model of spatial remodeling based on the principles of
evolutionary game theory. The analytically tractable framework describes the
spatial interactions between zones of bone resorption, bone formation and
quiescent bone, and explicitly accounts for regulation of remodeling by
bone-embedded, mechanotransducing osteocytes. Using tools from the theory of
interacting particle systems we systematically classified the different dynamic
regimes of the spatial model and identified regions of parameter space that
allow for coexistence of resorption, formation and quiescence, as observed in
physiological remodeling. In coexistence scenarios, three-dimensional
simulations revealed the emergence of sponge-like bone clusters. Comparison
between spatial and non-spatial dynamics revealed substantial differences and
suggested a stabilizing role of space. Our findings emphasize the importance of
accounting for spatial structure and bone-embedded osteocytes when modeling the
process of bone remodeling. Thanks to the lattice-based framework, the proposed
model can easily be coupled to a mechanical model of bone loading.


  Recently exciting progress has been made on protein contact prediction, but
the predicted contacts for proteins without many sequence homologs is still of
low quality and not very useful for de novo structure prediction. This paper
presents a new deep learning method that predicts contacts by integrating both
evolutionary coupling (EC) and sequence conservation information through an
ultra-deep neural network formed by two deep residual networks. This deep
neural network allows us to model very complex sequence-contact relationship as
well as long-range inter-contact correlation. Our method greatly outperforms
existing contact prediction methods and leads to much more accurate
contact-assisted protein folding. Tested on three datasets of 579 proteins, the
average top L long-range prediction accuracy obtained our method, the
representative EC method CCMpred and the CASP11 winner MetaPSICOV is 0.47, 0.21
and 0.30, respectively; the average top L/10 long-range accuracy of our method,
CCMpred and MetaPSICOV is 0.77, 0.47 and 0.59, respectively. Ab initio folding
using our predicted contacts as restraints can yield correct folds (i.e.,
TMscore>0.6) for 203 test proteins, while that using MetaPSICOV- and
CCMpred-predicted contacts can do so for only 79 and 62 proteins, respectively.
Further, our contact-assisted models have much better quality than
template-based models. Using our predicted contacts as restraints, we can (ab
initio) fold 208 of the 398 membrane proteins with TMscore>0.5. By contrast,
when the training proteins of our method are used as templates, homology
modeling can only do so for 10 of them. One interesting finding is that even if
we do not train our prediction models with any membrane proteins, our method
works very well on membrane protein prediction. Finally, in recent blind CAMEO
benchmark our method successfully folded 5 test proteins with a novel fold.


  Motivation: Gene regulatory interactions are of fundamental importance to
various biological functions and processes. However, only a few previous
computational studies have claimed success in revealing genome-wide regulatory
landscapes from temporal gene expression data, especially for complex
eukaryotes like human. Moreover, recent work suggests that these methods still
suffer from the curse of dimensionality if network size increases to 100 or
higher.
  Result: We present a novel scalable algorithm for identifying genome-wide
regulatory network structures. The highlight of our method is that its superior
performance does not degenerate even for a network size on the order of $10^4$,
and is thus readily applicable to large-scale complex networks. Such a
breakthrough is achieved by considering both prior biological knowledge and
multiple topological properties (i.e., sparsity and hub gene structure) of
complex networks in the regularized formulation. We also illustrate the
application of our algorithm in practice using the time-course expression data
from an influenza infection study in respiratory epithelial cells.
  Availability and Implementation: The algorithm described in this article is
implemented in MATLAB$^\circledR$. The source code is freely available from
https://github.com/Hongyu-Miao/DMI.git.
  Contact: jliu@cs.rochester.edu; hongyu.miao@uth.tmc.edu
  Supplementary information: Supplementary data are available online.


  Neglected tropical diseases (NTD), particularly vector-borne diseases (VBD),
account for a large proportion of the global disease burden, and their control
faces several challenges including diminishing human and financial resources
for those distressed from such diseases. Visceral Leishmaniasis (VL), the
second-largest parasitic killer in the world (after malaria) affects poor
populations in endemic countries and causes considerable cost to the affected
individuals and their society. Mathematical models can serve as a critical tool
for understanding the driving mechanisms of a NTD such as VL. The WHO promotes
integrated control programs for VL but this policy is not well supported by
systematic quantitative and dynamic evidence and so potential benefits of the
policy are limited. Moreover, mathematical models can be readily developed and
used to understand the functioning of the VL system cheaply and systematically.
The focus of this research is three-fold: (i) to identify non-traditional but
critical mechanisms for ongoing VL transmission in resource limited regions,
(ii) to review mathematical models used for other infectious diseases that have
the potential to capture identified factors of VL, and (iii) to suggest novel
quantitative models for understanding VL dynamics and for evaluating control
programs in such frameworks for achieveing VL elimination goals.


  We study stochastic evolutionary game dynamics in a population of finite
size. Individuals in the population are divided into two dynamically evolving
groups. The structure of the population is formally described by a
Wright-Fisher type Markov chain with a frequency dependent fitness. In a strong
selection regime that favors one of the two groups, we obtain qualitatively
matching lower and upper bounds for the fixation probability of the
advantageous population. In the infinite population limit we obtain an exact
result showing that a single advantageous mutant can invade an infinite
population with a positive probability. We also give asymptotically sharp
bounds for the fixation time distribution.


  RNA-Seq is a widely-used method for studying the behavior of genes under
different biological conditions. An essential step in an RNA-Seq study is
normalization, in which raw data are adjusted to account for factors that
prevent direct comparison of expression measures. Errors in normalization can
have a significant impact on downstream analysis, such as inflated false
positives in differential expression analysis. An under-emphasized feature of
normalization is the assumptions upon which the methods rely and how the
validity of these assumptions can have a substantial impact on the performance
of the methods. In this paper, we explain how assumptions provide the link
between raw RNA-Seq read counts and meaningful measures of gene expression. We
examine normalization methods from the perspective of their assumptions, as an
understanding of methodological assumptions is necessary for choosing methods
appropriate for the data at hand. Furthermore, we discuss why normalization
methods perform poorly when their assumptions are violated and how this causes
problems in subsequent analysis. To analyze a biological experiment,
researchers must select a normalization method with assumptions that are met
and that produces a meaningful measure of expression for the given experiment.


  We study a theoretical model for the toxin-antitoxin (hok/sok) mechanism for
plasmid maintenance in bacteria. Toxin-antitoxin systems enforce the
maintenance of a plasmid through post-segregational killing of cells that have
lost the plasmid. Key to their function is the tight regulation of expression
of a protein toxin by an sRNA antitoxin. Here, we focus on the nonlinear nature
of the regulatory circuit dynamics of the toxin-antitoxin mechanism. The
mechanism relies on a transient increase in protein concentration rather than
on the steady state of the genetic circuit. Through a systematic analysis of
the parameter dependence of this transient increase, we confirm some known
design features of this system and identify new ones: for an efficient
toxin-antitoxin mechanism, the synthesis rate of the toxin's mRNA template
should be lower that of the sRNA antitoxin, the mRNA template should be more
stable than the sRNA antitoxin, and the mRNA-sRNA complex should be more stable
than the sRNA antitoxin. Moreover, a short half-life of the protein toxin is
also beneficial to the function of the toxin-antitoxin system. In addition, we
study a therapeutic scenario in which a competitor mRNA is introduced to
sequester the sRNA antitoxin, causing the toxic protein to be expressed.


  The field of cryo-electron microscopy has made astounding advancements in the
past few years, mainly due to advancements in electron detectors' technology.
Yet, one of the key open challenges of the field remains the processing of
heterogeneous data sets, produced from samples containing particles at several
different conformational states. For such data sets, the algorithms must
include some classification procedure to identify homogeneous groups within the
data, so that the images in each group correspond to the same underlying
structure. The fundamental importance of the heterogeneity problem in
cryo-electron microscopy has drawn many research efforts, and resulted in
significant progress in classification algorithms for heterogeneous data sets.
While these algorithms are extremely useful and effective in practice, they
lack rigorous mathematical analysis and performance guarantees.
  In this paper, we attempt to make the first steps towards rigorous
mathematical analysis of the heterogeneity problem in cryo-electron microscopy.
To that end, we present an algorithm for processing heterogeneous data sets,
and prove accuracy and stability bounds for it. We also suggest an extension of
this algorithm that combines the classification and reconstruction steps. We
demonstrate it on simulated data, and compare its performance to the
state-of-the-art algorithm in RELION.


  The epidemic spreading has been studied for years by applying the mean-field
approach in both homogeneous case, where each node may get infected by an
infected neighbor with the same rate, and heterogeneous case, where the
infection rates between different pairs of nodes are different. Researchers
have discussed whether the mean-field approaches could accurately describe the
epidemic spreading for the homogeneous cases but not for the heterogeneous
cases. In this paper, we explore under what conditions the mean-field approach
could perform well when the infection rates are heterogeneous. In particular,
we employ the Susceptible-Infected-Susceptible (SIS) model and compare the
average fraction of infected nodes in the metastable state obtained by the
continuous-time simulation and the mean-field approximation. We concentrate on
an individual-based mean-field approximation called the N-intertwined Mean
Field Approximation (NIMFA), which is an advanced approach considered the
underlying network topology. Moreover, we consider not only the independent and
identically distributed (i.i.d.) infection rate but also the infection rate
correlated with the degree of the two end nodes. We conclude that NIMFA is
generally more accurate when the prevalence of the epidemic is higher. Given
the same effective infection rate, NIMFA is less accurate when the variance of
the i.i.d.\ infection rate or the correlation between the infection rate and
the nodal degree leads to a lower prevalence. Moreover, given the same actual
prevalence, NIMFA performs better in the cases: 1) when the variance of the
i.i.d.\ infection rates is smaller (while the average is unchanged); 2) when
the correlation between the infection rate and the nodal degree is positive.


  Extracting associations that recur across multiple studies while controlling
the false discovery rate is a fundamental challenge. Here, we consider an
extension of Efron's single-study two-groups model to allow joint analysis of
multiple studies. We assume that given a set of p-values obtained from each
study, the researcher is interested in associations that recur in at least
$k>1$ studies. We propose new algorithms that differ in how the study
dependencies are modeled. We compared our new methods and others using various
simulated scenarios. The top performing algorithm, SCREEN (Scalable
Cluster-based REplicability ENhancement), is our new algorithm that is based on
three stages: (1) clustering an estimated correlation network of the studies,
(2) learning replicability (e.g., of genes) within clusters, and (3) merging
the results across the clusters using dynamic programming. We applied SCREEN to
two real datasets and demonstrated that it greatly outperforms the results
obtained via standard meta-analysis. First, on a collection of 29 case-control
large-scale gene expression cancer studies, we detected a large up-regulated
module of genes related to proliferation and cell cycle regulation. These genes
are both consistently up-regulated across many cancer studies, and are well
connected in known gene networks. Second, on a recent pan-cancer study that
examined the expression profiles of patients with or without mutations in the
HLA complex, we detected an active module of up-regulated genes that are
related to immune responses. Thanks to our ability to quantify the false
discovery rate, we detected thrice more genes as compared to the original
study. Our module contains most of the genes reported in the original study,
and many new ones. Interestingly, the newly discovered genes are needed to
establish the connectivity of the module.


  Organisms shape their own environment, which in turn affects their survival.
This feedback becomes especially important for communities containing a large
number of species; however, few existing approaches allow studying this regime,
except in simulations. Here, we use methods of statistical physics to
analytically solve a classic ecological model of resource competition
introduced by MacArthur in 1969. We show that the non-intuitive phenomenology
of highly diverse ecosystems includes a phase where the environment constructed
by the community becomes fully decoupled from the outside world.


  Chemical reactions often proceed through the formation and the consumption of
intermediate species. An example is the creation and subsequent degradation of
the substrate-enzyme complexes in an enzymatic reaction. In this paper we
provide a setting, based on ordinary differential equations, in which the
presence of intermediate species has little effect on the overall dynamics of a
biological system. The result provides a method to perform model reduction by
elimination of intermediate species. We study the problem in a multiscale
setting, where the species abundances as well a the reaction rates scale to
different orders of magnitudes. The different time and concentration scales are
parameterised by a single parameter $N$. We show that a solution to the
original reaction system is uniformly approximated on compact time intervals to
a solution of a reduced reaction system without intermediates and to a solution
of a certain limiting reaction systems, which does not depend on $N$. Known
approximation techniques such as the theorems by Tikhonov and Fenichel cannot
readily be used in this framework.


  Amyloid fibrillation is a protein self-assembly phenomenon that is intimately
related to well-known human neurodegenerative diseases. During the past few
decades, striking advances have been achieved in our understanding of the
physical origin of this phenomenon and they constitute the contents of this
review. Starting from a minimal model of amyloid fibrils, we explore
systematically the equilibrium and kinetic aspects of amyloid fibrillation in
both dilute and semi-dilute limits. We then incorporate further molecular
mechanisms into the analyses. We also discuss the mathematical foundation of
kinetic modeling based on chemical mass-action equations, the quantitative
linkage with experimental measurements, as well as the procedure to perform
global fitting.


  Control problem in a biological system is the problem of finding an
interventional policy for changing the state of the biological system from an
undesirable state, e.g. disease, into a desirable healthy state. Boolean
networks are utilized as mathematical model for gene regulatory networks. This
paper provides an algorithm to solve the control problem in Boolean networks.
The proposed algorithm is implemented and applied on two biological systems:
T-cell receptor network and Drosophila melanogaster network. Results show that
the proposed algorithm works faster in solving the control problem over these
networks, while having similar accuracy, in comparison to previous methods.


  Increased availability of data and accessibility of computational tools in
recent years have created unprecedented opportunities for scientific research
driven by statistical analysis. Inherent limitations of statistics impose
constrains on reliability of conclusions drawn from data but misuse of
statistical methods is a growing concern. Significance, hypothesis testing and
the accompanying P-values are being scrutinized as representing most widely
applied and abused practices. One line of critique is that P-values are
inherently unfit to fulfill their ostensible role as measures of scientific
hypothesis's credibility. It has also been suggested that while P-values may
have their role as summary measures of effect, researchers underappreciate the
degree of randomness in the P-value. High variability of P-values would suggest
that having obtained a small P-value in one study, one is, nevertheless, likely
to obtain a much larger P-value in a similarly powered replication study. Thus,
"replicability of P-value" is itself questionable. To characterize P-value
variability one can use prediction intervals whose endpoints reflect the likely
spread of P-values that could have been obtained by a replication study.
Unfortunately, the intervals currently in use, the P-intervals, are based on
unrealistic implicit assumptions. Namely, P-intervals are constructed with the
assumptions that imply substantial chances of encountering large values of
effect size in an observational study, which leads to bias. As an alternative
to P-intervals, we develop a method that gives researchers flexibility by
providing them with the means to control these assumptions. Unlike endpoints of
P-intervals, endpoints of our intervals are directly interpreted as
probabilistic bounds for replication P-values and are resistant to selection
bias contingent upon approximate prior knowledge of the effect size
distribution.


  Increasing accessibility of data to researchers makes it possible to conduct
massive amounts of statistical testing. Rather than follow a carefully crafted
set of scientific hypotheses with statistical analysis, researchers can now
test many possible relations and let P-values or other statistical summaries
generate hypotheses for them. Genetic epidemiology field is an illustrative
case in this paradigm shift. Driven by technological advances, testing a
handful of genetic variants in relation to a health outcome has been abandoned
in favor of agnostic screening of the entire genome, followed by selection of
top hits, e.g., by selection of genetic variants with the smallest association
P-values. At the same time, nearly total lack of replication of claimed
associations that has been shaming the field turned to a flow of reports whose
findings have been robustly replicating. Researchers may have adopted better
statistical practices by learning from past failures, but we suggest that a
steep increase in the amount of statistical testing itself is an important
factor. Regardless of whether statistical significance has been reached, an
increased number of tested hypotheses leads to enrichment of smallest P-values
with genuine associations. In this study, we quantify how the expected
proportion of genuine signals (EPGS) among top hits changes with an increasing
number of tests. When the rate of occurrence of genuine signals does not
decrease too sharply to zero as more tests are performed, the smallest P-values
are increasingly more likely to represent genuine associations in studies with
more tests.


  RNA junctions are important structural elements of RNA molecules. They are
formed when three or more helices come together in three-dimensional space.
Recent studies have focused on the annotation and prediction of coaxial helical
stacking (CHS) motifs within junctions. Here we exploit such predictions to
develop an efficient alignment tool to handle RNA secondary structures with CHS
motifs. Specifically, we build upon our Junction-Explorer software for
predicting coaxial stacking and RNAJAG for modelling junction topologies as
tree graphs to incorporate constrained tree matching and dynamic programming
algorithms into a new method, called CHSalign, for aligning the secondary
structures of RNA molecules containing CHS motifs. Thus, CHSalign is intended
to be an efficient alignment tool for RNAs containing similar junctions.
Experimental results based on thousands of alignments demonstrate that CHSalign
can align two RNA secondary structures containing CHS motifs more accurately
than other RNA secondary structure alignment tools. CHSalign yields a high
score when aligning two RNA secondary structures with similar CHS motifs or
helical arrangement patterns, and a low score otherwise. This new method has
been implemented in a web server, and the program is also made freely
available, at http://bioinformatics.njit.edu/CHSalign/.


  Allostery plays a fundament role in most biological processes. However,
little theory is available to describe it outside of two-state models. Here we
use a statistical mechanical approach to show that the allosteric coupling
between two collective variables is not a single number, but instead a
two-dimensional thermodynamic coupling function that is directly related to the
mutual information from information theory and the copula density function from
probability theory. On this basis, we demonstrate how to quantify the
contribution of specific energy terms to this thermodynamic coupling function,
enabling a decomposition that reveals the mechanism of allostery. We illustrate
the thermodynamic coupling function and its use by showing how allosteric
coupling in the alanine dipeptide molecule contributes to the overall shape of
the {\Phi}/{\Psi} free energy surface, and by identifying the interactions that
are necessary for this coupling.


  Purpose: Factor analysis (FA) has become part and parcel in metabolic
syndrome (MBS) research. Both exploration- and confirmation-driven factor
analyzes are rampant. However, factor analytic results on MBS differ widely. A
situation that is at least in part attributable to misapplication of FA. Here,
our purpose is (i) to review factor analytic efforts in the study of MBS with
emphasis on misusage of the FA model and (ii) to propose an alternative factor
analytic strategy.
  Methods: The proposed factor analytic strategy consists of four steps and
confronts weaknesses in application of the FA model. At its heart lies the
explicit separation of dimensionality and pattern selection as well as the
direct evaluation of competing inequality-constrained loading patterns. A
high-profile MBS data set with anthropometric measurements on overweight
children and adolescents is reanalyzed using this strategy.
  Results: The reanalysis implied a more parsimonious constellation of
pathophysiological domains underlying phenotypic expressions of MBS than the
original analysis (and many other analyzes). The results emphasize correlated
factors of impaired glucose metabolism and impaired lipid metabolism.
  Conclusions: Pathophysiological domains underlying phenotypic expressions of
MBS included in the analysis are driven by multiple interrelated metabolic
impairments. These findings indirectly point to the possible existence of a
multifactorial aetiology.


  Almost all processes -- highly correlated, weakly correlated, or correlated
not at all---exhibit statistical fluctuations. Often physical laws, such as the
Second Law of Thermodynamics, address only typical realizations -- as
highlighted by Shannon's asymptotic equipartition property and as entailed by
taking the thermodynamic limit of an infinite number of degrees of freedom.
Indeed, our interpretations of the functioning of macroscopic thermodynamic
cycles are so focused. Using a recently derived Second Law for information
processing, we show that different subsets of fluctuations lead to distinct
thermodynamic functioning in Maxwellian Demons. For example, while typical
realizations may operate as an engine -- converting thermal fluctuations to
useful work -- even "nearby" fluctuations (nontypical, but probable
realizations) behave differently, as Landauer erasers -- converting available
stored energy to dissipate stored information. One concludes that ascribing a
single, unique functional modality to a thermodynamic system, especially one on
the nanoscale, is at best misleading, likely masking an array of simultaneous,
parallel thermodynamic transformations. This alters how we conceive of cellular
processes, engineering design, and evolutionary adaptation.


  Recently, the structural controllability theory has been introduced to
analyze the Protein-Protein Interaction (PPI) network. The indispensable nodes,
which their removal increase the number of driver nodes to control the network,
are found essential in PPI network. However, the PPI network is far from
complete and there may exist many false-positive or false-negative
interactions, which promotes us to question: are these indispensable nodes
robust to structural change? Here we systematically investigate the robustness
of indispensable nodes of PPI network by removing and adding possible
interactions. We found that the indispensable nodes are sensitive to the
structural change and very few edges can change the type of many indispensable
nodes. The finding may promote our understanding to the control principle of
PPI network.


  The prospect of RNA nanotechnology is increasing because of its numerous
potential applications especially in medical science. The spherical
Nanodiamonds (NDs) are becoming popular because of their lesser toxicity,
desirable mechanical, optical properties, functionality and available surface
areas. On other hand RNAs are stable, flexible and easy to bind to the NDs. In
this work, we have studied the tRNA dynamics on ND surface by high-resolution
quasi-elastic neutron scattering spectroscopy and all atom molecular dynamics
simulation technique to understand how the tRNA motion is affected by the
presence of ND. The flexibly of the tRNA is analyzed by the Mean square
displacement analysis that shows tRNA have a sharp increase around 230K in its
hydrated form. The intermediate scattering function (ISF) representing the tRNA
dynamics follows the logarithmic decay as proposed by the Mode Coupling theory
(MCT). But most importantly the tRNA dynamics is found to be faster in presence
of ND within 220K to 310K compared to the freestanding ones. This is as we have
shown, is because of the swollen RNA molecule due the introduction of
hydrophilic ND surface.


  Over the last few decades, ecologists have come to appreciate that key
ecological patterns, which describe ecological communities at relatively large
spatial scales, are not only scale dependent, but also intimately intertwined.
The relative abundance of species, which informs us about the commonness and
rarity of species, changes its shape from small to large spatial scales. The
average number of species as a function of area has a steep initial increase,
followed by decreasing slopes at large scales. Finally, if we find a species in
a given location, it is more likely we find an individual of the same species
close-by, rather than farther apart. Such spatial turnover depends on the
geographical distribution of species, which often are spatially aggregated.
This reverberates on the abundances as well as the richness of species within a
region, but so far it has been difficult to quantify such relationships. Within
a neutral framework, which considers all individuals competitively equivalent,
we introduce a spatial stochastic model, which phenomenologically accounts for
birth, death, immigration and local dispersal of individuals. We calculate the
pair correlation function, which encapsulates spatial turnover, and the
conditional probability to find a species with a certain population within a
given circular area. Also, we calculate the macro-ecological patterns, which we
have referred to above, and compare the analytical formul{\ae} with the
numerical integration of the model. Finally, we contrast the model predictions
with the empirical data for two lowland tropical forest inventories, showing
always a good agreement.


  In this paper we consider the problem of learning the
genetic-interaction-map, i.e., the topology of a directed acyclic graph (DAG)
of genetic interactions from noisy double knockout (DK) data. Based on a set of
well established biological interaction models we detect and classify the
interactions between genes. We propose a novel linear integer optimization
program called the Genetic-Interactions-Detector (GENIE) to identify the
complex biological dependencies among genes and to compute the DAG topology
that matches the DK measurements best. Furthermore, we extend the GENIE-program
by incorporating genetic-interactions-profile (GI-profile) data to further
enhance the detection performance. In addition, we propose a sequential
scalability technique for large sets of genes under study, in order to provide
statistically stressable results for real measurement data. Finally, we show
via numeric simulations that the GENIE-program as well as the GI-profile data
extended GENIE (GI-GENIE)-program clearly outperform the conventional
techniques and present real data results for our proposed sequential
scalability technique.


  This review provides a pedagogic and self-contained introduction to master
equations and to their representation by path integrals. We discuss analytical
and numerical methods for the solution of master equations, keeping our focus
on methods that are applicable even when stochastic fluctuations are strong.
The reviewed methods include the generating function technique and the Poisson
representation, as well as novel ways of mapping the forward and backward
master equations onto linear partial differential equations (PDEs). Spectral
methods, WKB approximations, and a variational approach have been proposed for
the analysis of the PDE obeyed by the generating function. After outlining
these methods, we solve the derived PDEs in terms of two path integrals. The
path integrals provide distinct exact representations of the conditional
probability distribution solving the master equations. We exemplify both path
integrals in analysing elementary chemical reactions. Furthermore, we review a
method for the approximation of rare event probabilities and derive path
integral representations of Fokker-Planck equations. To make our review
accessible to a broad community, we have used the language of probability
theory rather than quantum (field) theory. The probabilistic structures
underpinning various technical concepts, such as coherent states, the
Doi-shift, and normal-ordered observables, are thereby made explicit.


  Many human males produce dysfunctional sperm. Various plants frequently abort
pollen. Hybrid matings often produce sterile males. Widespread male sterility
is puzzling. Natural selection prunes reproductive failure. Puzzling failure
implies something that we do not understand about how organisms are designed.
Solving the puzzle reveals the hidden processes of design.


  Human language emerged abruptly. Diverse body forms evolved suddenly.
Seed-bearing plants spread rapidly. How do complex evolutionary innovations
arise so quickly? Resolving alternative claims remains difficult. The great
events of the past happened a long time ago. Cancer provides a model to study
evolutionary innovation. A tumor must evolve many novel traits to become an
aggressive cancer. I use what we know or could study about cancer to describe
the key processes of innovation. In general, evolutionary systems form a
hierarchy of recursive processes. Those recursive processes determine the rates
at which innovations are generated, spread and transmitted. I relate the
recursive processes to abrupt evolutionary innovation.


  Genome assembly from the high-throughput sequencing (HTS) reads is a
fundamental yet challenging computational problem. An intrinsic challenge is
the uncertainty caused by the widespread repetitive elements. Here we get
around the uncertainty using the notion of uniquely mapped (UM) reads, which
motivated the design of a new assembler BAUM. It mainly consists of two types
of iterations. The first type of iterations constructs initial contigs from a
reference, say a genome of a species that could be quite distant, by adaptive
read mapping, filtration by the reference's unique regions, and reference
updating. A statistical test is proposed to split the layouts at possible
structural variation sites. The second type of iterations includes mapping,
scaffolding/contig-extension, and contig merging. We extend each contig by
locally assembling the reads whose mates are uniquely mapped to an end of the
contig. Instead of the de Bruijn graph method, we take the
overlap-layout-consensus (OLC) paradigm. The OLC is implemented by parallel
computation, and has linear complexity with respect to the number of contigs.
The adjacent extended contigs are merged if their alignment is confirmed by the
adjusted gap distance. Throughout the assembling, the mapping criterion is
selected by probabilistic calculations. These innovations can be used
complementary to the existing de novo assemblers. Applying this novel method to
the assembly of wild rice Oryza longistaminata genome, we achieved much
improved contig N50, 18.8k, compared with other assemblers. The assembly was
further validated by contigs constructed from an independent library of long
454 reads.


  This paper applies the Bayesian Model Averaging (BMA) statistical ensemble
technique to estimate small molecule solvation free energies. There is a wide
range of methods available for predicting solvation free energies, ranging from
empirical statistical models to ab initio quantum mechanical approaches. Each
of these methods is based on a set of conceptual assumptions that can affect
predictive accuracy and transferability. Using an iterative statistical
process, we have selected and combined solvation energy estimates using an
ensemble of 17 diverse methods from the fourth Statistical Assessment of
Modeling of Proteins and Ligands (SAMPL) blind prediction study to form a
single, aggregated solvation energy estimate. The ensemble design process
evaluates the statistical information in each individual method as well as the
performance of the aggregate estimate obtained from the ensemble as a whole.
Methods that possess minimal or redundant information are pruned from the
ensemble and the evaluation process repeats until aggregate predictive
performance can no longer be improved. We show that this process results in a
final aggregate estimate that outperforms all individual methods by reducing
estimate errors by as much as 91% to 1.2 kcal/mol accuracy. We also compare our
iterative refinement approach to other statistical ensemble approaches and
demonstrate that this iterative process reduces estimate errors by as much as
61%. This work provides a new approach for accurate solvation free energy
prediction and lays the foundation for future work on aggregate models that can
balance computational cost with prediction accuracy.


  While most studies emphasize on certain aspects of Pathogen-Host Interactions
(PHI), such as the preferential attachment of bacteria or virus to its human
receptor homolog, studies have attempted to methodically classify interactions
among pathogenic proteins and their host proteins. Here we have analyzed 182
pathogens from The Pathogen-Host Interaction Search Tool (PHISTO) and could
identify the proteins/protein coding genes that act on both virus and bacteria.
Importantly there were few proteins viz. P53 (Tumor protein p53), NFKB1
(Nuclear factor of kappa light polypeptide gene enhancer in B-cells 1), GBLP
(Guanine nucleotide-binding protein subunit beta-2-like-1), TOX4 (TOX high
mobility group box family member 4), PDIA1 (Protein disulfide-isomerase
precursor), MHY9 (Myosin 9), RAC1 (Ras-related C3 botulinum toxin substrate 1),
CCAR2 (Cell cycle and apoptosis regulator protein 2) and ILF3 (Interleukin
enhancer binding factor 3) that were more susceptible to both bacterial and
viral pathogens. Identification of such important interacting proteins (IIPs)
can elicit significant insights for better understanding the molecular
mechanisms of such pathogens that interact with the human host.


  A combined dynamics consisting of Brownian motion and L\'evy flights is
exhibited by a variety of biological systems performing search processes.
Assessing the search reliability of ever locating the target and the search
efficiency of doing so economically of such dynamics thus poses an important
problem. Here we model this dynamics by a one-dimensional fractional
Fokker-Planck equation combining unbiased Brownian motion and L\'evy flights.
By solving this equation both analytically and numerically we show that the
superposition of recurrent Brownian motion and L\'evy flights with stable
exponent $\alpha<1$, by itself implying zero probability of hitting a point on
a line, lead to transient motion with finite probability of hitting any point
on the line. We present results for the exact dependence of the values of both
the search reliability and the search efficiency on the distance between the
starting and target positions as well as the choice of the scaling exponent
$\alpha$ of the L\'evy flight component.


  Functional studies in Drosophila have been key for establishing a role for
the septin family of proteins in animal cell division and thus extending for
the first time observations from the budding yeast to animal cells. Visualizing
the distribution of specific septins in different Drosophila tissues and, in
particular, in the Drosophila embryo, together with biochemical and mutant
phenotype data, has contributed important advances to our understanding of
animal septin biology, suggesting roles in processes other than in cytokinesis.
Septin localization using immunofluorescence assays has been possible due to
the generation of antibodies against different Drosophila septins. The recent
availability of lines expressing fluorescent protein fusions of specific
septins further promises to facilitate studies on septin dynamics. Here, we
provide protocols for preparing early Drosophila embryos to visualize septins
using immunofluorescence assays and live fluorescence microscopy. The genetic
tractability of the Drosophila embryo together with its amenability to
high-resolution fluorescence microscopy promises to provide novel insights into
animal septin structure and function.


  Septins are guanine nucleotide-binding proteins that are conserved from fungi
to humans. Septins assemble into heterooligomeric complexes and higher-order
structures with key roles in various cellular functions including cell
migration and division. The mechanisms by which septins assemble and interact
with other cytoskeletal elements like actin remain elusive. A powerful approach
to address this question is by cell-free reconstitution of pu- rified
cytoskeletal proteins combined with fluorescence microscopy. Here, we describe
procedures for the purification of recombinant Drosophila and human septin
hexamers from Escherichia coli and reconstitution of actin-septin coassembly.
These procedures can be used to compare assembly of Drosophila and human
septins and their coassembly with the actin cytoskeleton by total internal
reflection fluorescence microscopy.


  Animals live in groups to defend against predation and to obtain food.
However, for some animals --- especially ones that spend long periods of time
feeding --- there are costs if a group chooses to move on before their
nutritional needs are satisfied. If the conflict between feeding and keeping up
with a group becomes too large, it may be advantageous to some animals to split
into subgroups of animals with similar nutritional needs. We model the costs
and benefits of splitting by a herd of cows using a cost function (CF) that
quantifies individual variation in hunger, desire to lie down, and predation
risk. We model the costs associated with hunger and lying desire as the
standard deviations of individuals within a group, and we model predation risk
as an inverse exponential function of group size. We minimize the cost function
over all plausible groups that can arise from a given herd and study the
dynamics of group splitting. We explore our model using two examples: (1) we
consider group switching and group fission in a herd of relatively homogeneous
cows; and (2) we examine a herd with an equal number of adult males (larger
animals) and adult females (smaller animals).


  In this paper, we consider the impact of time delay by media on the control
of the disease. We set up a class of SISM epidemic model with the time delay
and the cumulative density of awareness caused by media. The sufficient
condition of global asymptotic stability of disease-free equilibrium is
approved. We get the global stability of the epidemic equilibrium and the
existence conditions of Hopf bifurcation. Numerical simulations are presented
to illustrate the analytical results. Finally, we analyze the influence of
parameters on the control of infectious disease by combining the data of H1N1.
By shortening the time of media lag, increasing transmission rate of media and
the implementation rate of the media project, the spread of disease will be
controlled effectively.


  In this paper, we study a single-species population model with pulse toxicant
input to a small polluted environment. The intrinsic rate of population is
affected by environment and toxin in organisms. The toxin in organisms is
influenced by toxin in environment and the food chain. A new mathematical model
is derived. By the Pulse Compare Theorem, we find the surviving threshold of
the population and obtain the sufficient conditions of persistence and
extinction of the population.


  Many cell signaling pathways involve the diffusion of messengers that
bind/unbind to intracellular components. Quantifying their net transport rate
under different conditions, then requires having separate estimates of their
free diffusion coefficient and binding/unbinding rates. In this paper, we show
how performing sets of Fluorescence Correlation Spectroscopy (FCS) experiments
under different conditions, it is possible to quantify free diffusion
coefficients and on and off rates of reaction-diffusion systems. We develop the
theory and present a practical implementation for the case of the universal
second messenger, calcium (Ca$^{2+}$) and single-wavelength dyes that increase
their fluorescence upon Ca$^{2+}$ binding. We validate the approach with
experiments performed in aqueous solutions containing Ca$^{2+}$ and Fluo4
dextran (both in its High and Low Affinity versions). Performing FCS
experiments with tetramethylrhodamine-dextran in Xenopus laevis oocytes, we
infer the corresponding free diffusion coefficients in the cytosol of these
cells. Our approach can be extended to other physiologically relevant
reaction-diffusion systems to quantify biophysical parameters that determine
the dynamics of various variables of interest.


  The transmission dynamics of Tuberculosis (TB) involve complex
epidemiological and socio-economical interactions between individuals living in
highly distinct regional conditions. The level of exogenous reinfection and
first time infection rates within high-incidence settings may influence the
impact of control programs on TB prevalence. This study aims at enhancing the
understanding of TB dynamics via the study of scenarios, within {\it
simplified}, two patch, risk-defined environments, in the presence of short
term mobility and variations in reinfection and infection rates. The modeling
framework captures the role of individuals' `daily' dynamics within and between
places of residency, work or business via the proportion of time spent in
residence and as visitors to TB-risk environments (patches). As a result, the
{\it effective population size} of Patch $i$ (home of $i$-residents) at time
$t$ must account for visitors and residents of Patch $i$, at time $t$. The
impact that {\it effective population size} and the distribution of {\it
individuals' residence times} in different patches have on TB transmission and
control are studied using selected scenarios where risk is defined by the
estimated or perceive first time infection and/or exogenous re-infection rates.
Our results suggest that, under certain conditions, allowing infected
individuals to move from high to low TB prevalence areas (for example via the
sharing of treatment and isolation facilities) may lead to a reduction in the
total TB prevalence in the overall, here two-patch, population.


  Transposable elements, or transposons, are DNA sequences that can jump from
site to site in the genome during the life cycle of a cell, usually encoding
the very enzymes which perform their excision. However, some transposons are
parasitic, relying on the enzymes produced by the regular transposons. In this
case, we show that a stochastic model, which takes into account the small copy
numbers of the transposons in a cell, predicts noise-induced predator-prey
oscillations with a characteristic time scale that is much longer than the cell
replication time, indicating that the state of the predator-prey oscillator is
stored in the genome and transmitted to successive generations. Our work
demonstrates the important role of number fluctuations in the expression of
mobile genetic elements, and shows explicitly how ecological concepts can be
applied to the dynamics and fluctuations of living genomes.


  Spatially extended versions of the cyclic-dominance Rock-Paper-Scissors model
have traveling wave (in one dimension) and spiral (in two dimensions) behavior.
The far field of the spirals behave like traveling waves, which themselves have
profiles reminiscent of heteroclinic cycles. We compute numerically a nonlinear
dispersion relation between the wavelength and wavespeed of the traveling
waves, and, together with insight from heteroclinic bifurcation theory and
further numerical results from 2D simulations, we are able to make predictions
about the overall structure and stability of spiral waves in 2D cyclic
dominance models.


  The cytoskeleton protein actin assembles into large bundles when supporting
stresses in the cell, but grows into a fine branched network to induce cell
motion. Such self-organization processes are studied in artificial networks of
cytoskeleton proteins with thick actin bundles and large motor protein
aggregates to enable optical observation. The effect of motor aggregate size on
the cytoskeleton mechanical properties is studied here in networks comprised of
much smaller motor assemblies. Large motor protein clusters are known to
increase the stiffness of actin based networks by introducing tension and
additional cross-linking cites. We find that these effects are universal to
actin gels regardless of actin bundle thickness and motor aggregate size and
are relevant, therefore, to a wide range of cytoskeleton based cellular
processes. In contrast, motor induced active fluctuations depend significantly
on motor assembly size, featuring unique non-Gaussian statistics at high
concentrations of small assemblies.


  Flexibility-rigidity index (FRI) has been developed as a robust, accurate and
efficient method for macromolecular thermal fluctuation analysis and B-factor
prediction. The performance of FRI depends on its formulations of rigidity
index and flexibility index. In this work, we introduce alternative rigidity
and flexibility formulations. The structure of the classic Gaussian surface is
utilized to construct a new type of rigidity index, which leads to a new class
of rigidity densities with the classic Gaussian surface as a special case.
Additionally, we introduce a new type of flexibility index based on the domain
indicator property of normalized rigidity density. These generalized FRI (gFRI)
methods have been extensively validated by the B-factor predictions of 364
proteins. Significantly outperforming the classic Gaussian network model (GNM),
gFRI is a new generation of methodologies for accurate, robust and efficient
analysis of protein flexibility and fluctuation. Finally, gFRI based molecular
surface generation and flexibility visualization are demonstrated.


  In recent years the spread of the ectoparasitic mite Varroa destructor has
become the most serious threat to worldwide apiculture. In the model presented
here we extend the bee population dynamics with mite viral epidemiology
examined in an earlier paper by allowing a bee-dependent mite population size.
The results of the analysis match field observations well and give a clear
explanation of how Varroa affects the epidemiology of certain naturally
occurring bee viruses, causing considerable damages to colonies. The model
allows only four possible stable equilibria, using known field parameters. The
first one contains only the thriving healthy bees. Here the disease is
eradicated and also the mites are wiped out. Alternatively, we find the
equilibrium still with no mite population, but with endemic disease among the
thriving bee population. Thirdly, infected bees coexist with the mites in the
Varroa invasion scenario; in this situation the disease invades the hive,
driving the healthy bees to extinction and therefore affecting all the bees.
Coexistence is also possible, with both populations of bees and mites thriving
and with the disease endemically affecting both species. The analysis is in
line with field observations in natural honey bee colonies. Namely, these
diseases are endemic and if the mite population is present, necessarily the
viral infection occurs. Further, in agreement with the fact that the presence
of Varroa increments the viral transmission, the whole bee population may
become infected when the disease vector is present in the beehive. Also, a low
horizontal transmission rate of the virus among the honey bees will help in
protecting the bee colonies from Varroa infestation and viral epidemics.


  Food webs have been found to exhibit remarkable motif profiles, patterns in
the relative prevalences of all possible three-species sub-graphs, and this has
been related to ecosystem properties such as stability and robustness.
Analysing 46 food webs of various kinds, we find that most food webs fall into
one of two distinct motif families. The separation between the families is well
predicted by a global measure of hierarchical order in directed networks -
trophic coherence. We find that trophic coherence is also a good predictor for
the extent of omnivory, defined as the tendency of species to feed on multiple
trophic levels. We compare our results to a network assembly model that admits
tunable trophic coherence via a single free parameter. The model is able to
generate food webs in either of the two families by varying this parameter, and
correctly classifies almost all the food webs in our database. This establishes
a link between global order and local preying patterns in food webs.


  It is known that learning of players who interact in a repeated game can be
interpreted as an evolutionary process in a population of ideas. These
analogies have so far mostly been established in deterministic models, and
memory loss in learning has been seen to act similarly to mutation in
evolution. We here propose a representation of reinforcement learning as a
stochastic process in finite "populations of ideas". The resulting birth-death
dynamics has absorbing states and allows for the extinction or fixation of
ideas, marking a key difference to mutation-selection processes in finite
populations. We characterize the outcome of evolution in populations of ideas
for several classes of symmetric and asymmetric games.


  Advances in nanotechnology lead to an increasing interest in how
nanoparticles interact with biomembranes. Nanoparticles are wrapped
spontaneously by biomembranes if the adhesive interactions between the
particles and membranes compensate for the cost of membrane bending. In the
last years, the cooperative wrapping of spherical nanoparticles in membrane
tubules has been observed in experiments and simulations. For spherical
nanoparticles, the stability of the particle-filled membrane tubules strongly
depends on the range of the adhesive particle-membrane interactions. In this
article, we show via modeling and energy minimization that elongated and patchy
particles are wrapped cooperatively in membrane tubules that are highly stable
for all ranges of the particle-membrane interactions, compared to individual
wrapping of the particles. The cooperative wrapping of linear chains of
elongated or patchy particles in membrane tubules may thus provide an efficient
route to induce membrane tubulation, or to store such particles in membranes.


  Culture evolves, not just in the trivial sense that cultures change over
time, but also in the strong sense that such change is governed by Darwinian
principles. Both biological and cultural evolution are essentially cumulative
selection processes in which information (whether genetic or cultural) is
sieved, retained and then transmitted to the next generation. In both domains
such a process will result in recognizable lineages and tree-like phylogenies
so characteristic of Darwinian evolution. Because a principle of inheritance
(i.e., faithful replication of information) holds in both domains, we may trace
back particular transmission histories and identify the forces that influenced
them. The idea that culture evolves is quite old, but only in recent years
there has been a serious effort to turn this idea into science. This article
offers a concise analysis of how a rudimentary idea gradually evolved into a
thriving research program.


  There is an abundance of information about the genetic basis, physiological
and molecular mechanisms of bacterial pathogenesis. In contrast, relatively
little is known about population dynamic processes, by which bacteria colonize
hosts and invade tissues and cells and thereby cause disease. In an article
published in 1978, Moxon and Murphy presented evidence that, when inoculated
intranasally with a mixture streptomycin sensitive and resistant (Sm$^S$ and
Sm$^R$) and otherwise isogenic stains of Haemophilus influenzae type b (Hib),
neonatal rats develop a bacteremic infection that often is dominated by only
one strain, Sm$^S$ or Sm$^R$. After rulling out other possibilities through
years of related experiments, the field seems to have settled on a plausible
explanation for this phenomenon: the first bacterium to invade the host
activates the host immune response that `shuts the door' on the second invading
strain. To explore this hypothesis in a necessarily quantitative way, we
modeled this process with a set of mixed stochastic and deterministic
differential equations. Our analysis of the properties of this model with
realistic parameters suggests that this hypothesis cannot explain the
experimental results of Moxon and Murphy, and in particular the observed
relationship between the frequency of different types of blood infections
(bacteremias) and the inoculum size. We propose modifications to the model that
come closer to explaining these data. However, the modified and better fitting
model contradicts the common theory of independent action of individual
bacteria in establishing infections. We discuss the implications of these
results.


  Investigating the pleiotropic effects of genetic variants can increase
statistical power, provide important information to achieve deep understanding
of the complex genetic structures of disease, and offer powerful tools for
designing effective treatments with fewer side effects. However, the current
multiple phenotype association analysis paradigm lacks breadth (number of
phenotypes and genetic variants jointly analyzed at the same time) and depth
(hierarchical structure of phenotype and genotypes). A key issue for high
dimensional pleiotropic analysis is to effectively extract informative internal
representation and features from high dimensional genotype and phenotype data.
To explore multiple levels of representations of genetic variants, learn their
internal patterns involved in the disease development, and overcome critical
barriers in advancing the development of novel statistical methods and
computational algorithms for genetic pleiotropic analysis, we proposed a new
framework referred to as a quadratically regularized functional CCA (QRFCCA)
for association analysis which combines three approaches: (1) quadratically
regularized matrix factorization, (2) functional data analysis and (3)
canonical correlation analysis (CCA). Large-scale simulations show that the
QRFCCA has a much higher power than that of the nine competing statistics while
retaining the appropriate type 1 errors. To further evaluate performance, the
QRFCCA and nine other statistics are applied to the whole genome sequencing
dataset from the TwinsUK study. We identify a total of 79 genes with rare
variants and 67 genes with common variants significantly associated with the 46
traits using QRFCCA. The results show that the QRFCCA substantially outperforms
the nine other statistics.


  This paper compares a number of mathematical models for the Calvin cycle of
photosynthesis and presents theorems on the existence and stability of steady
states of these models. Results on five-variable models in the literature are
surveyed. Next a number of larger models related to one introduced by
Pettersson and Ryde-Pettersson are discussed. The mathematical nature of this
model is clarified, showing that it is naturally defined as a system of
differential-algebraic equations. It is proved that there are choices of
parameters for which this model admits more than one positive steady state.
This is done by analysing the limit where the storage of sugars from the cycle
as starch is shut down. There is also a discussion of the minimal models for
the cycle due to Hahn.


  HIV is a deadly virus transmitted either through having of unprotected sex,
mother to child transmission, sharing of unsterilized objects that is capable
of making cut or wounds on the body, through blood or bodily fluid
transmission. AIDS has no permanent cure but remedies that help in suppressing
the power effect of the virus are available. Previous studies has shown that
the epidemic claimed more lives via vertical transmission, mother-to-child
transmission, blood transfusion, sexual intercourse and injection drug users.
The associated models were derived using diagnosis and treatments records of
confirmed status of HIV/AIDS clients. A new model analysis were proposed that
could handle cases of HIV/AIDS routes of transmission and treatments via
vertical, mother-to-child, blood transfusion, sexual intercourse, injection
drug users that were not considered in previous studies. The new model was
solved using Fractional Differential Transformation by Caputo sense algorithm.
The new model is a major contribution for optimal assessment, monitoring,
evaluation, control and management of HIV/AIDS, with 75% accuracy.


  Predicting the 3D structure of a macromolecule, such as a protein or an RNA
molecule, is ranked top among the most difficult and attractive problems in
bioinformatics and computational biology. Its importance comes from the
relationship between the 3D structure and the function of a given protein or
RNA. 3D structures also help to find the ligands of the protein, which are
usually small molecules, a key step in drug design. Unfortunately, there is no
shortcut to accurately obtain the 3D structure of a macromolecule. Many
physical measurements of macromolecular 3D structures cannot scale up, due to
their large labor costs and the requirements for lab conditions.
  In recent years, computational methods have made huge progress due to advance
in computation speed and machine learning methods. These methods only need the
sequence information to predict 3D structures by employing various mathematical
models and machine learning methods. The success of computational methods is
highly dependent on a large database of the proteins and RNA with known
structures.
  However, the performance of computational methods are always expected to be
improved. There are several reasons for this. First, we are facing, and will
continue to face sparseness of data.Secondly, the 3D structure space is too
large for our computational capability.
  The two obstacles can be removed by knowledge-based methods, which combine
knowledge learned from the known structures and biologists' knowledge of the
folding process of protein or RNA. In the dissertation, I will present my
results in building a knowledge-based method by using machine learning methods
to tackle this problem. My methods include the knowledge constraints on
intermediate states, which can highly reduce the solution space of a protein or
RNA, in turn increasing the efficiency of the structure folding method and
improving its accuracy.


  How to best exploit patchy resources? This long-standing question belongs to
the extensively studied class of explore/exploit problems that arise in a wide
range of situations, from animal foraging, to robotic exploration, and to human
decision processes. Despite its broad relevance, the issue of optimal
exploitation has previously only been tackled through two paradigmatic limiting
models---patch-use and random search---that do not account for the interplay
between searcher motion within patches and resource depletion. Here, we bridge
this gap by introducing a minimal patch exploitation model that incorporates
this coupling: the searcher depletes the resources along its random-walk
trajectory within a patch and travels to a new patch after it takes
$\mathcal{S}$ consecutive steps without finding resources. We compute the
distribution of the amount of resources $F_t$ consumed by time $t$ for this
non-Markovian random walker and show that exploring multiple patches is
beneficial. In one dimension, we analytically derive the optimal strategy to
maximize $F_t$. We show that this strategy is robust with respect to the
distribution of resources within patches and the criterion for leaving a given
patch. We also show that $F_t$ can be optimized in the ecologically-relevant
case of two-dimensional patchy environments.


  Using force as a probe to map the folding landscapes of RNA molecules has
become a reality thanks to major advances in single molecule pulling
experiments. Although the unfolding pathways under tension are complicated to
predict studies in the context of proteins have shown that topology plays is
the major determinant of the unfolding landscapes. By building on this finding
we study the responses of RNA molecules to force by adapting Gaussian network
model (GNM) that represents RNAs using a bead-spring network with isotropic
interactions. Cross-correlation matrices of residue fluctuations, which are
analytically calculated using GNM even upon application of mechanical force,
show distinct allosteric communication as RNAs rupture. The model is used to
calculate the force-extension curves at full thermodynamic equilibrium, and the
corresponding unfolding pathways of four RNA molecules subject to a
quasi-statically increased force.Our study finds that the analysis using GNM
captures qualitatively the unfolding pathway of \emph{T}. ribozyme elucidated
by the optical tweezers measurement. However, the simple model is not
sufficient to capture subtle features, such as bifurcation in the unfolding
pathways or the ion effects, in the forced-unfolding of RNAs.


  Productive societies feature high levels of cooperation and strong
connections between individuals. Public Goods Games (PGGs) are frequently used
to study the development of social connections and cooperative behavior in
model societies. In such games, contributions to the public good are made only
by cooperators, while all players, including defectors, can reap public goods
benefits. Classic results of game theory show that mutual defection, as opposed
to cooperation, is the Nash Equilibrium of PGGs in well-mixed populations,
where each player interacts with all others. In this paper, we explore the
coevolutionary dynamics of a low information public goods game on a network
without spatial constraints in which players adapt to their environment in
order to increase individual payoffs. Players adapt by changing their
strategies, either to cooperate or to defect, and by altering their social
connections. We find that even if players do not know other players' strategies
and connectivity, cooperation can arise and persist despite large short-term
fluctuations.


  Living cells use readout molecules to record the state of receptor proteins,
similar to measurements or copies in typical computational devices. But is this
analogy rigorous? Can cells be optimally efficient, and if not, why? We show
that, as in computation, a canonical biochemical readout network generates
correlations; extracting no work from these correlations sets a lower bound on
dissipation. For general input, the biochemical network cannot reach this
bound, even with arbitrarily slow reactions or weak thermodynamic driving. It
faces an accuracy-dissipation trade-off that is qualitatively distinct from and
worse than implied by the bound, and more complex steady-state copy processes
cannot perform better. Nonetheless, the cost remains close to the thermodynamic
bound unless accuracy is extremely high. Additionally, we show that
biomolecular reactions could be used in thermodynamically optimal devices under
exogenous manipulation of chemical fuels, suggesting an experimental system for
testing computational thermodynamics.


  Cancer cells evolve through random somatic mutations. "Beneficial" mutations
which disrupt key pathways (e.g. cell cycle regulation) are subject to natural
selection. Multiple mutations may lead to the same "beneficial" effect, in
which case there is no selective advantage to having more than one of these
mutations. Hence we are interested in finding sets of genes whose mutations are
approximately mutually exclusive (anti-co-occurring) within the TCGA Pancancer
dataset. In principle, finding the best set is NP Hard. Nevertheless, we will
show how a new Mutation anti-co-OCcurrence Algorithm (MOCA) provides an
effective greedy search and testing algorithm with guaranteed control of the
familywise error rate or false discovery rate, by combining some
under-appreciated ideas from frequentist hypothesis testing. These ideas
include: (a) A novel exact conditional test for the tendency of multiple sets
to have a large/small union/intersection, which generalises Fisher's exact test
of 2x2 tables. (b) Randomised hypothesis tests for discrete distributions. (c)
Stouffer's method for combining p-values. (d) Weighted multiple hypothesis
testing. A new approach to setting a-priori weights which generates additional
implicit hypothesis tests is suggested, and allows us to preserve almost all
statistical power when testing pairs despite introducing a combinatorially
large number of additional hypotheses.