  The chromosome is a key player of cell physiology, and its dynamics provides
valuable information about its physical organization. In both prokaryotes and
eukaryotes, the short-time motion of chromosomal loci has been described as a
Rouse model in a simple or viscoelastic medium. However, little emphasis has
been put on the role played by the folded organization of chromosomes on the
local dynamics. Clearly, stress-propagation, and thus dynamics, must be
affected by such organization, but a theory allowing to extract such
information from data, e.g.\ of two-point correlations, is lacking. Here, we
describe a theoretical framework able to answer this general polymer dynamics
question, and we provide a general scaling analysis of the stress-propagation
time between two loci at a given arclength distance along the chromosomal
coordinate. The results suggest a precise way to detect folding information
from the dynamical coupling of chromosome segments. Additionally, we realize
this framework in a specific theoretical model of a polymer with variable-range
interactions in a viscoelastic medium characterized by a tunable scaling
exponent, where we derive analytical estimates of the correlation functions.


  Recent experiments and simulations have demonstrated that proteins can fold
on the ribosome. However, the extent and generality of fitness effects
resulting from co-translational folding remain open questions. Here we report a
genome-wide analysis that uncovers evidence of evolutionary selection for
co-translational folding. We describe a robust statistical approach to identify
loci within genes that are both significantly enriched in slowly translated
codons and evolutionarily conserved. Surprisingly, we find that domain
boundaries can explain only a small fraction of these conserved loci. Instead,
we propose that regions enriched in slowly translated codons are associated
with co-translational folding intermediates, which may be smaller than a single
domain. We show that the intermediates predicted by a native-centric model of
co-translational folding account for the majority of these loci across more
than 500 E. coli proteins. By making a direct connection to protein folding,
this analysis provides strong evidence that many synonymous substitutions have
been selected to optimize translation rates at specific locations within genes.
More generally, our results indicate that kinetics, and not just
thermodynamics, can significantly alter the efficiency of self-assembly in a
biological context.


  A stochastic version of the Barkai-Leibler model of chemotaxis receptors in
{\it E. coli} is studied here to elucidate the effects of intrinsic network
noise in their conformational dynamics. It was originally proposed to explain
the robust and near-perfect adaptation of {\it E. coli} observed across a wide
range of spatially uniform attractant/repellent (ligand) concentrations. A
receptor is either active or inactive and can stochastically switch between the
two states. Enzyme CheR methylates inactive receptors while CheB demethylates
active ones and the probability for it to be active depends on its level of
methylation and ligandation. A simple version of the model with two methylation
sites per receptor (M=2) shows zero-order ultrasensitivity (ZOU) akin to the
classical 2-state model of covalent modification studied by Goldbeter and
Koshland (GK). For extremely small and large ligand concentrations, the system
reduces to two 2-state GK modules. A quantitative measure of the spontaneous
fluctuations in activity (variance) estimated mathematically under linear noise
approximation (LNA) is found to peak near the ZOU transition. The variance is a
weak, non-monotonic and decreasing functions of ligand and receptor
concentrations. Gillespie simulations for M=2 show excellent agreement with
analytical results obtained under LNA. Numerical results for M=2, 3 and 4 show
ZOU in mean activity; the variance is found to be smaller for larger M. The
magnitude of receptor noise deduced from available experimental data is
consistent with our predictions. A simple analysis of the downstream signaling
pathway shows that this noise is large enough to have a beneficial effect on
the motility of the organism. The response of mean receptor activity to small
time-dependent changes in the external ligand concentration, computed within
linear response theory, is found to have a bilobe form.


  Protein-ligand binding is a fundamental biological process that is paramount
to many other biological processes, such as signal transduction, metabolic
pathways, enzyme construction, cell secretion, gene expression, etc. Accurate
prediction of protein-ligand binding affinities is vital to rational drug
design and the understanding of protein-ligand binding and binding induced
function. Existing binding affinity prediction methods are inundated with
geometric detail and involve excessively high dimensions, which undermines
their predictive power for massive binding data. Topology provides an ultimate
level of abstraction and thus incurs too much reduction in geometric
information. Persistent homology embeds geometric information into topological
invariants and bridges the gap between complex geometry and abstract topology.
However, it over simplifies biological information. This work introduces
element specific persistent homology (ESPH) to retain crucial biological
information during topological simplification. The combination of ESPH and
machine learning gives rise to one of the most efficient and powerful tools for
revealing protein-ligand binding mechanism and for predicting binding
affinities.


  Understanding the evolution of binary traits, which affects the birth and
survival of species and also the rate of molecular evolution, remains
challenging. A typical example is the evolution of mating systems in plant
species. In this work, we present a probabilistic modeling framework for binary
trait, random species trees, in which the number of species and their traits
are represented by a two-type, continuous time Markov branching process. We
develop our model by considering the impact of mating systems on dN/dS, the
ratio of nonsynonymous to synonymous substitutions. A methodology is introduced
which enables us to match model parameters with parameter estimates from
phylogenetic tree data. The properties obtained from the model are applied to
outcrossing and selfing species trees in the Geraniaceae and Solanaceae family.
This allows us to investigate not only the branching tree rates, but also the
mutation rates and the intensity of selection.


  The process of desertification in the semi-arid climatic zone is considered
by many as a catastrophic regime shift, since the positive feedback of
vegetation density on growth rates yields a system that admits alternative
steady states. Some support to this idea comes from the analysis of static
patterns, where peaks of the vegetation density histogram were associated with
these alternative states. Here we present a large-scale empirical study of
vegetation dynamics, aimed at identifying and quantifying directly the effects
of positive feedback. To do that, we have analyzed vegetation density across
$~2.5 \times 10^6 \ \rm{km}^2$ of the African Sahel region, with spatial
resolution of $30 \times 30$ meters, using three consecutive snapshots. The
results are mixed. The local vegetation density (measured at a single pixel)
moves towards the average of the corresponding rainfall line, indicating a
purely negative feedback. On the other hand, the chance of spatial clusters (of
many "green" pixels) to expand in the next census is growing with their size,
suggesting some positive feedback. We show that these apparently contradicting
results emerge naturally in a model with positive feedback and strong
demographic stochasticity, a model that allows for a catastrophic shift only in
a certain range of parameters. Static patterns, like the double peak in the
histogram of vegetation density, are shown to vary between censuses, with no
apparent correlation with the actual dynamical features.


  A free energy landscape estimation-method based on Bayesian inference is
presented and used for comparing the efficiency of thermally enhanced sampling
methods with respect to regular molecular dynamics, where the simulations are
carried out on two binding states of calmodulin. The proposed free energy
estimation method (the GM method) is compared to other estimators using a toy
model showing that the GM method provides a robust estimate not subject to
overfitting. The continuous nature of the GM method, as well as predictive
inference on the number of basis functions, provide better estimates on sparse
data. We find that the free energy diffusion proper- ties determine sampling
method effectiveness, such that the diffusion dominated apo-calmodulin is most
efficiently sampled by regular molecular dynamics, while the holo with its
rugged free energy landscape is better sampled by enhanced methods.


  UPGMA (Unweighted Pair Group Method with Arithmetic Mean) is a widely used
clustering method. Here we show that UPGMA is a greedy heuristic for the
normalized equidistant minimum evolution (NEME) problem, that is, finding a
rooted tree that minimizes the minimum evolution score relative to the
dissimilarity matrix among all rooted trees with the same leaf-set in which all
leaves have the same distance to the root. We prove that the NEME problem is
NP-hard. In addition, we present some heuristic and approximation algorithms
for solving the NEME problem, including a polynomial time algorithm that yields
a binary, rooted tree whose NEME score is within O(log^2 n) of the optimum. We
expect that these results to eventually provide further insights into the
behavior of the UPGMA algorithm.


  Aims: Ischaemic cardiomyopathy (ICM) leads to impaired contraction and
ventricular dysfunction causing high rates of morbidity and mortality.
Epigenomics allows the identification of epigenetic signatures in human
diseases. We analyse the differential epigenetic patterns of ASB gene family in
ICM patients and relate these alterations to their haemodynamic and functional
status. Methods and Results: Epigenomic analysis was carried out using 16 left
ventricular (LV) tissue samples, 8 from ICM patients undergoing heart
transplantation and 8 from control (CNT) subjects without cardiac disease. We
increased the sample size up to 13 ICM and 10 CNT for RNA-sequencing and to 14
ICM for pyrosequencing analyses. We found a hypermethylated profile
(cg11189868) in the ASB1 gene that showed a differential methylation of 0.26
beta difference, P < 0.05. This result was validated by pyrosequencing
technique (0.23 beta difference, P < 0.05). Notably, the methylation pattern
was strongly related to LV ejection fraction (r = -0.849, P = 0.008) stroke
volume (r = -0.929, P = 0.001) and end-systolic and diastolic LV diameters (r =
-0.743, P = 0.035 for both). ASB1 showed a down regulation in mRNA levels (-1.2
fold, P < 0.05). Conclusion: Our findings link a specific ASB1 methylation
pattern to LV structure and performance in end-stage ICM, opening new
therapeutic opportunities and providing new insights regarding which is the
functionally relevant genome in the ischemic failing myocardium. Keywords:
ischaemic cardiomyopathy; epigenomics; heart failure; left ventricular
dysfunction; stroke volume; ASB1.


  Epidemic spreading has been intensively studied in SIS epidemic model.
Although the mean-field theory of SIS model has been widely used in the
research, there is a lack of comparative results between different theoretical
calculations, and the differences between them should be systematically
explained. In this paper, we have compared different theoretical solutions for
mean-field theory and explained the underlying reason. We first describe the
differences between different equations for mean-field theory in different
networks. The results show that the difference between mean-field reaction
equations is due to the different probability consideration for the infection
process. This finding will help us to design better theoretical solutions for
epidemic models.


  Gene Ontology (GO) terms are frequently used to score alignments between
protein-protein interaction (PPI) networks. Methods exist to measure the GO
similarity between two proteins in isolation, but pairs of proteins in a
network alignment are not isolated: each pairing is implicitly dependent upon
every other pairing via the alignment itself. Current methods fail to take into
account the frequency of GO terms across the networks, and attempt to account
for common GO terms in an ad hoc fashion by imposing arbitrary rules on when to
"allow" GO terms based on their location in the GO hierarchy, rather than using
readily available frequency information in the PPI networks themselves. Here we
develop a new measure, NetGO, that naturally weighs infrequent, informative GO
terms more heavily than frequent, less informative GO terms, without requiring
arbitrary cutoffs. In particular, NetGO down-weights the score of frequent GO
terms according to their frequency in the networks being aligned. This is a
global measure applicable only to alignments, independent of pairwise GO
measures, in the same sense that the edge-based EC or S3 scores are global
measures of topological similarity independent of pairwise topological
similarities. We demonstrate the superiority of NetGO by creating alignments of
predetermined quality based on homologous pairs of nodes and show that NetGO
correlates with alignment quality much better than any existing GO-based
alignment measures. We also demonstrate that NetGO provides a measure of
taxonomic similarity between species, consistent with existing taxonomic
measures--a feature not shared with existing GO-based network alignment
measures. Finally, we re-score alignments produced by almost a dozen aligners
from a previous study and show that NetGO does a better job than existing
measures at separating good alignments from bad ones.


  Network alignment (NA) compares networks with the goal of finding a node
mapping that uncovers highly similar (conserved) network regions. Existing NA
methods are homogeneous, i.e., they can deal only with networks containing
nodes and edges of one type. Due to increasing amounts of heterogeneous network
data with nodes or edges of different types, we extend three recent
state-of-the-art homogeneous NA methods, WAVE, MAGNA++, and SANA, to allow for
heterogeneous NA for the first time. We introduce several algorithmic
novelties. Namely, these existing methods compute homogeneous graphlet-based
node similarities and then find high-scoring alignments with respect to these
similarities, while simultaneously maximizing the amount of conserved edges.
Instead, we extend homogeneous graphlets to their heterogeneous counterparts,
which we then use to develop a new measure of heterogeneous node similarity.
Also, we extend $S^3$, a state-of-the-art measure of edge conservation for
homogeneous NA, to its heterogeneous counterpart. Then, we find high-scoring
alignments with respect to our heterogeneous node similarity and edge
conservation measures. In evaluations on synthetic and real-world biological
networks, our proposed heterogeneous NA methods lead to higher-quality
alignments and better robustness to noise in the data than their homogeneous
counterparts. The software and data from this work is available upon request.


  Single-molecule FRET is widely used to study helicases by detecting distance
changes between a fluorescent donor and an acceptor anchored to overhangs of a
forked DNA duplex. However, it has lacked single-base pair (1-bp) resolution
required for revealing stepping dynamics in unwinding because FRET signals are
usually blurred by thermal fluctuations of the overhangs. We designed a
nanotensioner in which a short DNA is bent to exert a force on the overhangs,
just as in optical/magnetic tweezers. The strategy improved the resolution of
FRET to 0.5 bp, high enough to uncover the differences in DNA unwinding by
yeast Pif1 and E. coli RecQ whose unwinding behaviors cannot be differentiated
by currently practiced methods. We found that Pif1 exhibits 1-bp-stepping
kinetics, while RecQ breaks 1 bp at a time but questers the nascent nucleotides
and releases them randomly. The high-resolution data allowed us to propose a
three-parameter model to quantitatively interpret the apparently different
unwinding behaviors of the two helicases which belong to two superfamilies.


  The ability to measure the transcriptomes of single cells has only been
feasible for a few years, and is becoming an extremely popular assay. While
many types of analysis and questions can be answered using single cell
RNA-sequencing, a central focus is the ability to survey the diversity of cell
types within a sample. Unbiased and reproducible cataloging of distinct cell
types requires large numbers of cells. Technological developments and protocol
improvements have fuelled a consistent exponential increase in the numbers of
cells studied in single cell RNA-seq analyses. In this perspective, we will
highlight the key technological developments which have enabled this growth in
data.


  We prove that the probability substitution matrices obtained from a
continuous-time Markov chain form a multiplicatively closed set if and only if
the rate matrices associated to the chain form a linear space spanning a Lie
algebra. The key original contribution we make is to overcome an obstruction,
due to the presence of inequalities that are unavoidable in the probabilistic
application, that prevents free manipulation of terms in the
Baker-Campbell-Haursdorff formula.


  The notion of entropy is shared between statistics and thermodynamics, and is
fundamental to both disciplines. This makes statistical problems particularly
suitable for reaction network implementations. In this paper we show how to
perform a statistical operation known as Information Projection or E projection
with stochastic mass-action kinetics. Our scheme encodes desired conditional
distributions as the equilibrium distributions of reaction systems. To our
knowledge this is a first scheme to exploit the inherent stochasticity of
reaction networks for information processing. We apply this to the problem of
an artificial cell trying to infer its environment from partial observations.


  Constraints on changes in expression levels across all cell components
imposed by the steady growth of cells have recently been discussed both
experimentally and theoretically. By assuming a small environmental
perturbation and considering a linear response to it, a common proportionality
in such expression changes was derived and partially verified by experimental
data. Here, we examined global protein expression in ${\it Escherichia coli}$
under various environmental perturbations. Remarkably they are proportional
across components, even though these environmental changes are not small and
cover different types of stresses, while the proportion coefficient corresponds
to the change in growth rate. However, since such global proportionality is not
generic to all systems under a condition of steady growth, a new conceptual
framework is needed. We hypothesized that such proportionality is a result of
evolution. To validate this hypothesis, we analyzed a cell model with a huge
number of components that reproduces itself via a catalytic reaction network,
and confirmed that the common proportionality in the concentrations of all
components is shaped through evolutionary processes to maximize cell growth
(and therefore fitness) under a given environmental condition. Further, we
found that the concentration changes across all components in response to
environmental and evolutionary changes are constrained along a one-dimensional
major axis within a huge-dimensional state space. Based on these observations,
we propose a theory in which high-dimensional phenotypic changes after
evolution are constrained along a one-dimensional major axis that correlates
with the growth rate, which can explain broad experimental and numerical
results.


  High throughput sequencing is a technology that allows for the generation of
millions of reads of genomic data regarding a study of interest, and data from
high throughput sequencing platforms are usually count compositions. Subsequent
analysis of such data can yield information on tran- scription profiles,
microbial diversity, or even relative cellular abundance in culture. Because of
the high cost of acquisition, the data are usually sparse, and always contain
far fewer observations than variables. However, an under-appreciated pathology
of these data are their often unbalanced nature: i.e, there is often systematic
variation between groups simply due to presence or absence of features, and
this variation is important to the biological interpretation of the data. A
simple example would be comparing transcriptomes of yeast cells with and
without a gene knockout. This causes samples in the comparison groups to
exhibit widely varying centres. This work extends a previously described
log-ratio transformation method that allows for variable comparisons between
samples in a Bayesian compositional context. We demonstrate the pathology in
modelled and real unbalanced experimental designs to show how this dramatically
causes both false negative and false positive inference. We then introduce
several approaches to demonstrate how the pathologies can be addressed. An
extreme example is presented where only the use of a predefined basis is
appropriate. The transformations are implemented as an extension to a general
compositional data analysis tool known as ALDEx2 which is available on
Bioconductor.


  Topological phylogenetic trees can be assigned edge weights in several
natural ways, highlighting different aspects of the tree. Here the rooted
triple and quartet metrizations are introduced, and applied to formulate novel
fast methods of inferring large trees from rooted triple and quartet data.
These methods can be applied in new statistically consistent procedures for
inference of a species tree from gene trees under the multispecies coalescent
model.


  Phylogenetic networks are a generalization of phylogenetic trees that are
used to represent non-tree-like evolutionary histories that arise in organisms
such as plants and bacteria, or uncertainty in evolutionary histories. An
\emph{unrooted} phylogenetic network on a nonempty, finite set $X$ of taxa, or
\emph{network}, is a connected graph in which every vertex has degree 1 or 3
and whose leaf-set is $X$. It is called a \emph{phylogenetic tree} if the
underlying graph is a tree. In this paper we consider properties of
\emph{tree-based networks}, that is, networks that can be constructed by adding
edges into a phylogenetic tree. We show that although they have some properties
in common with their rooted analogues which have recently drawn much attention
in the literature, they have some striking differences in terms of both their
structural and computational properties. We expect that our results could
eventually have applications to, for example, detecting horizontal gene
transfer or hyrbridization which are important factors in the evolution of many
organisms.


  We consider excursions for a class of stochastic processes describing a
population of discrete individuals experiencing density-limited growth, such
that the population has a finite carrying capacity and behaves qualitatively
like the classical logistic model when the carrying capacity is large. Being
discrete and stochastic, however, our population nonetheless goes extinct in
finite time. We present results concerning the maximum of the population prior
to extinction in the large population limit, from which we obtain establishment
probabilities and upper bounds for the process, as well as estimates for the
waiting time to establishment and extinction. As a consequence, we show that
conditional upon establishment, the stochastic logistic process will with high
probability greatly exceed carrying capacity an arbitrary number of times prior
to extinction.


  We study an optimal control problem arising from a resource allocation
problem in cellular metabolism. A minimalistic model that describes the
production of enzymatic vs. non-enzymatic biomass components from a single
nutrient source is introduced. The basic growth modes with this model are
linear growth, where only the non-enzymatic component is produced, and
exponential growth with only enzymatic components being produced. Using
Pontryagin's maximum principle, we derive the optimal growth trajectory
depending on the model's parameters. It turns out that depending on the
parameters, either a single growth mode is optimal, or otherwise the optimal
solution is a concatenation of exponential growth with linear growth.
Importantly, on the short time scale, the choice of growth mode depends only on
catalytic rate constants and biomass weights of the two component types,
whereas on longer time scales, where the nutrient amount becomes limiting, also
the yield coefficients play a role.


  We investigate the impact of misinformation about the contact structure on
the ability to predict disease outbreaks. We base our study on 31 empirical
temporal networks and tune the frequencies in errors in the node identities or
timestamps of contacts. We find that for both these spreading scenarios, the
maximal misprediction of both the outbreak size and time to extinction follows
an stretched exponential convergence as a function of the error frequency. We
furthermore determine the temporal-network structural factors influencing the
parameters of this convergence.


  Transcription is regulated through binding factors to gene promoters to
activate or repress expression, however, the mechanisms by which factors find
targets remain unclear. Using single-molecule fluorescence microscopy, we
determined in vivo stoichiometry and spatiotemporal dynamics of a GFP tagged
repressor, Mig1, from a paradigm signaling pathway of Saccharomyces cerevisiae.
We find the repressor operates in clusters, which upon extracellular signal
detection, translocate from the cytoplasm, bind to nuclear targets and
turnover. Simulations of Mig1 configuration within a 3D yeast genome model
combined with a promoter-specific, fluorescent translation reporter confirmed
clusters are the functional unit of gene regulation. In vitro and structural
analysis on reconstituted Mig1 suggests that clusters are stabilized by
depletion forces between intrinsically disordered sequences. We observed
similar clusters of a co-regulatory activator from a different pathway,
supporting a generalized cluster model for transcription factors that reduces
promoter search times through intersegment transfer while stabilizing gene
expression.


  Cells receive signaling molecules by receptors and relay information via
sensory networks so that they can respond properly depending on the type of
signal. Recent studies have shown that cells can extract multi-dimensional
information from dynamical concentration patterns of signaling molecules. We
herein study how biochemical systems can process multi-dimensional information
embedded in dynamical patterns. We model the decoding networks by linear
response functions, and optimize the functions with the calculus of variations
to maximize the mutual information between patterns and output. We find that,
when the noise intensity is lower, decoders with different linear response
functions, i.e., distinct decoders, can extract much information. However, when
the noise intensity is higher, distinct decoders do not provide the maximum
amount of information. This indicates that, when transmitting information by
dynamical patterns, embedding information in multiple patterns is not optimal
when the noise intensity is very large. Furthermore, we explore the biochemical
implementations of these decoders using control theory and demonstrate that
these decoders can be implemented biochemically through the modification of
cascade-type networks, which are prevalent in actual signaling pathways.


  Desertification in dryland ecosystems is considered to be a major
environmental threat that may lead to devastating consequences. The concern
increases when the system admits two alternative steady states and the
transition is abrupt and irreversible (catastrophic shift). However, recent
studies show that the inherent stochasticity of the birth-death process, when
superimposed on the presence of an absorbing state, may lead to a continuous
(second order) transition even if the deterministic dynamics supports a
catastrophic transition. Following these works we present here a numerical
study of a one-dimensional stochastic desertification model, where the
deterministic predictions are confronted with the observed dynamics. Our
results suggest that a stochastic spatial system allows for a propagating front
only when its active phase invades the inactive (desert) one. In the extinction
phase one observes transient front propagation followed by a global collapse.
In the presence of a seed bank the vegetation state is shown to be more robust
against demographic stochasticity, but the transition in that case still
belongs to the directed percolation equivalence class.


  The medical research facilitates to acquire a diverse type of data from the
same individual for particular cancer. Recent studies show that utilizing such
diverse data results in more accurate predictions. The major challenge faced is
how to utilize such diverse data sets in an effective way. In this paper, we
introduce a multiple kernel based pipeline for integrative analysis of
high-throughput molecular data (somatic mutation, copy number alteration, DNA
methylation and mRNA) and clinical data. We apply the pipeline on Ovarian
cancer data from TCGA. After multiple kernels have been generated from the
weighted sum of individual kernels, it is used to stratify patients and predict
clinical outcomes. We examine the survival time, vital status, and neoplasm
cancer status of each subtype to verify how well they cluster. We have also
examined the power of molecular and clinical data in predicting dichotomized
overall survival data and to classify the tumor grade for the cancer samples.
It was observed that the integration of various data types yields higher
log-rank statistics value. We were also able to predict clinical status with
higher accuracy as compared to using individual data types.


  Following \cite{ipel1}, we consider a nonlinear SIS-type nonlocal system
describing the spread of epidemics on networks, assuming nonlimited
transmission, We prove local existence of a unique solution for any diffusion
coefficients and global existence in the case of equal diffusion coefficients.
Next we study the asymptotic behaviour of the solution and show that the
disease-free equilibrium (DFE) is linearly and globally asymptotically stable
when the total mean population is small. Finally, we prove that the solution of
the system converge to the $DFE$.


  Multiplex and multi-directional control of metabolic pathways is crucial for
metabolic engineering to improve product yield of fuels, chemicals, and
pharmaceuticals. To achieve this goal, artificial transcriptional regulators
such as CRISPR-based transcription regulators have been developed to
specifically activate or repress genes of interest. Here, we found that by
deploying guide RNAs to target on DNA sites at different locations of genetic
cassettes, we could use just one synthetic CRISPR-based transcriptional
regulator to simultaneously activate and repress gene expressions. By using the
pairwise datasets of guide RNAs and gene expressions, we developed a
data-driven predictive model to rationally design this system for fine-tuning
expression of target genes. We demonstrated that this system could achieve
programmable control of metabolic fluxes when using yeast to produce versatile
chemicals. We anticipate that this master CRISPR-based transcription regulator
will be a valuable addition to the synthetic biology toolkit for metabolic
engineering, speeding up the design-build-test cycle in industrial
biomanufacturing as well as generating new biological insights on the fates of
eukaryotic cells.


  Many enhanced sampling methods, such as Umbrella Sampling, Metadynamics or
Variationally Enhanced Sampling, rely on the identification of appropriate
collective variables. For proteins, even small ones, finding appropriate
collective variables has proven challenging. Here we suggest that the NMR $S^2$
order parameter can be used to this effect. We trace the validity of this
statement to the suggested relation between $S^2$ and entropy. Using the $S^2$
order parameter and a surrogate for the protein enthalpy in conjunction with
Metadynamics or Variationally Enhanced Sampling we are able to reversibly fold
and unfold a small protein and draw its free energy at a fraction of the time
that is needed in unbiased simulations. From a more conceptual point of view
this implies describing folding as a resulting from a trade off between entropy
and enthalpy. We also use $S^2$ in combination with the free energy flooding
method to compute the unfolding rate of this peptide. We repeat this
calculation at different temperatures to obtain the unfolding activation
energy.


  Hyperamylinemia, a condition characterized by above-normal blood levels of
the pancreas-derived peptide amylin, accompanies obesity and precedes type II
diabetes. Human amylin oligomerizes easily and can deposit in the pancreas,
brain, and heart, where they have been associated with calcium dysregulation.
In the heart, accumulating evidence suggests that human amylin oligomers form
modestly cation-selective, voltage-dependent ion channels that embed in the
cell sarcolemma (SL). The oligomers increase membrane conductance in a
dose-dependent manner, which is correlated with elevated cytosolic calcium
(Ca2+). Our core hypothesis therefore is that non-selective inward Ca2+
conduction afforded by human amylin oligomers increase cytosolic and
sarcoplasmic reticulum (SR) Ca2+ load, which thereby magnifies intracellular
Ca2+ transients. Questions remain however regarding the mechanism of
amylin-induced Ca2+ dysregulation, including whether enhanced SL Ca2+ influx is
sufficient to elevate cytosolic Ca2+ load, and if so, how might amplified Ca2+
transients perturb Ca2+-dependent cardiac pathways. To investigate these
questions, we modified a computational model of cardiomyocytes Ca2+ signaling
to reflect experimentally-measured changes in SL membrane permeation and
decreased Sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) function
stemming from acute and transgenic human amylin peptide exposure. With this
model, we confirmed the hypothesis that increasing SL permeation alone was
sufficient to enhance Ca2+ transient amplitudes. Our model indicated that
amplified cytosolic transients are driven by increased Ca2+ loading of the
sarcoplasmic reticulum and may contribute to the Ca2+-dependent activation of
calmodulin. These findings suggest that increased membrane permeation induced
by deposition of amylin oligomers contributes to Ca2+ dysregulation in
pre-diabetes.


  In a standard bifurcation of a dynamical system, the stationary points (or
more generally attractors) change qualitatively when varying a control
parameter. Here we describe a novel unusual effect, when the change of a
parameter, e.g. a growth rate, does not influence the stationary states, but
nevertheless leads to a qualitative change of dynamics. For instance, such a
dynamic transition can be between the convergence to a stationary state and a
strong increase without stationary states, or between the convergence to one
stationary state and that to a different state. This effect is illustrated for
a dynamical system describing two symbiotic populations, one of which exhibits
a growth rate larger than the other one. We show that, although the stationary
states of the dynamical system do not depend on the growth rates, the latter
influence the boundary of the basins of attraction. This change of the basins
of attraction explains this unusual effect of the quantitative change of
dynamics by growth rate variation.


  Gene drives have the potential to rapidly replace a harmful wild-type allele
with a gene drive allele engineered to have desired functionalities. However,
an accidental or premature release of a gene drive construct to the natural
environment could damage an ecosystem irreversibly. Thus, it is important to
understand the spatiotemporal consequences of the super-Mendelian population
genetics prior to potential applications. Here, we employ a reaction-diffusion
model for sexually reproducing diploid organisms to study how a locally
introduced gene drive allele spreads to replace the wild-type allele, even
though it possesses a selective disadvantage $s>0$. Using methods developed by
N. Barton and collaborators, we show that socially responsible gene drives
require $0.5<s<0.697$, a rather narrow range. In this "pushed wave" regime, the
spatial spreading of gene drives will be initiated only when the initial
frequency distribution is above a threshold profile called "critical
propagule", which acts as a safeguard against accidental release. We also study
how the spatial spread of the pushed wave can be stopped by making gene drives
uniquely vulnerable ("sensitizing drive") in a way that is harmless for a
wild-type allele. Finally, we show that appropriately sensitized drives in two
dimensions can be stopped even by imperfect barriers perforated by a series of
gaps.


  Polysaccharides (carbohydrates) are key regulators of a large number of cell
biological processes. However, precise biochemical or genetic manipulation of
these often complex structures is laborious and hampers experimental
structure-function studies. Molecular Dynamics (MD) simulations provide a
valuable alternative tool to generate and test hypotheses on saccharide
function. Yet, currently used MD force fields often overestimate the
aggregation propensity of polysaccharides, affecting the usability of those
simulations. Here we tested MARTINI, a popular coarse-grained (CG) force field
for biological macromolecules, for its ability to accurately represent
molecular forces between saccharides. To this end, we calculated a
thermodynamic solution property, the second virial coefficient of the osmotic
pressure ($B_{22}$). Comparison with light scattering experiments revealed a
non-physical aggregation of a prototypical polysaccharide in MARTINI, pointing
at an imbalance of the non-bonded solute-solute, solute-water, and water-water
interactions. This finding also applies to smaller oligosaccharides which were
all found to aggregate in simulations even at moderate concentrations, well
below their solubility limit. Finally, we explored the influence of the
Lennard-Jones (LJ) interaction between saccharide molecules and propose a
simple scaling of the LJ interaction strength that makes MARTINI more reliable
for the simulation of saccharides.


  The intrinsic stochasticity of gene expression can give rise to large
fluctuations and rare events that drive phenotypic variation in a population of
genetically identical cells. Characterizing the fluctuations that give rise to
such rare events motivates the analysis of large deviations in stochastic
models of gene expression. Recent developments in non-equilibrium statistical
mechanics have led to a framework for analyzing Markovian processes conditioned
on rare events and for representing such processes by conditioning-free driven
Markovian processes. We use this framework, in combination with approaches
based on queueing theory, to analyze a general class of stochastic models of
gene expression. Modeling gene expression as a Batch Markovian Arrival Process
(BMAP), we derive exact analytical results quantifying large deviations of
time-integrated random variables such as promoter activity fluctuations. We
find that the conditioning-free driven process can also be represented by a
BMAP that has the same form as the original process, but with renormalized
parameters. The results obtained can be used to quantify the likelihood of
large deviations, to characterize system fluctuations conditional on rare
events and to identify combinations of model parameters that can give rise to
dynamical phase transitions in system dynamics.


  Predicting the biological function of molecules, be it proteins or drug-like
compounds, from their atomic structure is an important and long-standing
problem. Function is dictated by structure, since it is by spatial interactions
that molecules interact with each other, both in terms of steric
complementarity, as well as intermolecular forces. Thus, the electron density
field and electrostatic potential field of a molecule contain the "raw
fingerprint" of how this molecule can fit to binding partners. In this paper,
we show that deep learning can predict biological function of molecules
directly from their raw 3D approximated electron density and electrostatic
potential fields. Protein function based on EC numbers is predicted from the
approximated electron density field. In another experiment, the activity of
small molecules is predicted with quality comparable to state-of-the-art
descriptor-based methods. We propose several alternative computational models
for the GPU with different memory and runtime requirements for different sizes
of molecules and of databases. We also propose application-specific
multi-channel data representations. With future improvements of training
datasets and neural network settings in combination with complementary
information sources (sequence, genomic context, expression level), deep
learning can be expected to show its generalization power and revolutionize the
field of molecular function prediction.


  We use optimal control theory with the purpose of finding the best spraying
policy with the aim of at least to minimize and possibly to eradicate the
number of parasites, i.e., the prey for the spiders living in an
agroecosystems. Two different optimal control problems are posed and solved,
and their implications discussed.


  Using topological summaries of gene trees as a basis for species tree
inference is a promising approach to obtain acceptable speed on genomic-scale
datasets, and to avoid some undesirable modeling assumptions. Here we study the
probabilities of splits on gene trees under the multispecies coalescent model,
and how their features might inform species tree inference. After investigating
the behavior of split consensus methods, we investigate split invariants ---
that is, polynomial relationships between split probabilities. These invariants
are then used to show that, even though a split is an unrooted notion, split
probabilities retain enough information to identify the rooted species tree
topology for trees of more than 5 taxa, with one possible 6-taxon exception.


  We investigate the rates of drug resistance acquisition in a natural
population using molecular epidemiological data from Bolivia. First, we study
the rate of direct acquisition of double resistance from the double sensitive
state within patients and compare it to the rates of evolution to single
resistance. In particular, we address whether or not double resistance can
evolve directly from a double sensitive state within a given host. Second, we
aim to understand whether the differences in mutation rates to rifampicin and
isoniazid resistance translate to the epidemiological scale. Third, we estimate
the proportion of MDR TB cases that are due to the transmission of MDR strains
compared to acquisition of resistance through evolution. To address these
problems we develop a model of TB transmission in which we track the evolution
of resistance to two drugs and the evolution of VNTR loci. However, the
available data is incomplete, in that it is recorded only {for a fraction of
the population and} at a single point in time. The likelihood function induced
by the proposed model is computationally prohibitive to evaluate and
accordingly impractical to work with directly. We therefore approach
statistical inference using approximate Bayesian computation techniques.


  Synergies between evolutionary game theory and statistical physics have
significantly improved our understanding of public cooperation in structured
populations. Multiplex networks, in particular, provide the theoretical
framework within network science that allows us to mathematically describe the
rich structure of interactions characterizing human societies. While research
has shown that multiplex networks may enhance the resilience of cooperation,
the interplay between the overlap in the structure of the layers and the
control parameters of the corresponding games has not yet been investigated.
With this aim, we consider here the public goods game on a multiplex network,
and we unveil the role of the number of layers and the overlap of links, as
well as the impact of different synergy factors in different layers, on the
onset of cooperation. We show that enhanced public cooperation emerges only
when a significant edge overlap is combined with at least one layer being able
to sustain some cooperation by means of a sufficiently high synergy factor. In
the absence of either of these conditions, the evolution of cooperation in
multiplex networks is determined by the bounds of traditional network
reciprocity with no enhanced resilience. These results caution against overly
optimistic predictions that the presence of multiple social domains may in
itself promote cooperation, and they help us better understand the complexity
behind prosocial behavior in layered social systems.


  Energy saving mechanisms in nature allow following organisms to expend less
energy than leaders. Queues, or ordered rows of individuals, may form when
organisms exploit the available energy saving mechanism while travelling at
near-maximal sustainable metabolic capacities; compact clusters form when group
members travel well below maximal sustainable metabolic capacities. The group
size range, given here as the ratio of the difference between the size of the
largest and smallest group members, and the size of the largest member (as a
per cent), has been hypothesized to correspond proportionately to the energy
saving quantity because weaker, smaller, individuals sustain speeds of
stronger, larger, individuals by exploiting the energy saving mechanism (as a
per cent). During migration, small individuals outside this range may perish,
or form sub-groups, or simply not participate in migratory behavior. We
approximate drag forces for leading and following individuals in queues of the
late Devonian (approx. 370 Ma) trilobite Trimerocephalus chopini. Applying data
from literature of R. herculea, a living crustacean, we approximate the
hypothetical walking speed and maximal speeds for T. chopini. Findings
reasonably support the hypothesis: among the population of fossilized queues of
T. chopini reported by Kin and Blazejowski (2013), trilobite size range was 75%
while the size range within queues, was 63%; this corresponds reasonably with
drag reductions in following positions that permit approx. 61.5% energy saving
for trilobites following others in optimal low-drag positions. We model
collective trilobite behavior associated with hydrodynamic drafting.


  In this conceptual paper we propose to explore the analogy between
ontic/epistemic description of quantum phenomena and interrelation between
dynamics of conformational and functional states of proteins. Another new idea
is to apply theory of automata to model the latter dynamics. In our model
protein's behavior is modeled with the aid of two dynamical systems, ontic and
epistemic, which describe evolution of conformational and functional states of
proteins, respectively. The epistemic automaton is constructed from the ontic
automaton on the basis of functional (observational) equivalence relation on
the space of ontic states. This reminds a few approaches to emergent quantum
mechanics in which a quantum (epistemic) state is treated as representing a
class of prequantum (ontic) states. This approach does not match to the
standard {\it protein structure-function paradigm.} However, it is perfect for
modeling of behavior of intrinsically disordered proteins. Mathematically space
of protein's ontic states (conformational states) is modeled with the aid of
$p$-adic numbers or more general ultrametric spaces encoding the internal
hierarchical structure of proteins. Connection with theory of $p$-adic
dynamical systems is briefly discussed.


  Even this saying itself is a variant of a similar statement attributed to
Bernard of Chartres in the 12th Century, and inspired the title for a book by
Steven Hawking and an album by Oasis. Creative ideas beget other creative ideas
and, as a result, modifications accumulate, and we see an overall increase in
the complexity of cultural novelty over time, a phenomenon sometimes referred
to as the ratchet effect (Tomasello, Kruger, & Ratner, 1993). Although we may
never meet the people or objects that creatively influence us, by assimilating
what we encounter around us and bringing to bear our own insights and
perspectives, we all contribute in our own way, however small, to a second
evolutionary process -- the evolution of culture. This chapter explores how we
can better understand culture by understanding the creative processes that fuel
it, and better understand creativity by examining it from its cultural context.
First, we look at some theoretical frameworks for how culture evolves and what
these frameworks imply for the role of creativity. Then we will see how
questions about the relationship between creativity and cultural evolution have
been addressed using an agent-based model. We will also discuss studies of how
creative outputs are influenced, in perhaps unexpected ways, by other ideas and
individuals, and how individual creative styles "peek through" cultural outputs
in different domains.


  Nonequilibrium energetics of single molecule translational motor kinesin was
investigated by measuring heat dissipation from the violation of the
fluctuation-response relation of a probe attached to the motor using optical
tweezers. The sum of the dissipation and work did not amount to the input free
energy change, indicating large hidden dissipation exists. Possible sources of
the hidden dissipation were explored by analyzing the Langevin dynamics of the
probe, which incorporates the two-state Markov stepper as a kinesin model. We
conclude that internal dissipation is dominant.


  Understanding gene regulation is a fundamental step towards understanding of
how cells function and respond to environmental cues and perturbations. An
important step in this direction is to infer the transcription factor-gene
regulatory network (GRN). However gene regulatory networks are typically
constructed disregarding the fact that regulatory programs are conditioned on
tissue type, developmental stage, sex, and other factors. Collecting multitude
of features required for a reliable construction of GRNs such as physical
features and functional features for every context of interest is costly.
Therefore we need methods that is able to use the knowledge of a
context-agnostic network for construction of a context specific regulatory
network.
  To address this challenge we developed a computational approach that uses
context specific expression data and a GRN constructed in a different but
related context to construct a context specific GRN. Our method, NetREX, is
inspired by network component analysis that estimates TF activities and their
influences on target genes given predetermined topology of a TF-gene network.
To predict a network under a different condition, NetREX removes the
restriction that the topology of the TF-gene network is fixed and allows for
adding and removing edges to that network. To solve the corresponding
optimization problem, we provide a general mathematical framework allowing use
of the recently proposed PALM technique and develop a convergent algorithm.
  We tested our NetREX on simulated data and subsequently applied it to gene
expression data in adult females from Drosophila deletion (DrosDel) panel. The
networks predicted by NetREX showed higher biological consistency than
alternative approaches. In addition, we used the list of recently identified
targets of the Doublesex (DSX) transcription factor to demonstrate the
predictive power of our method.


  Most infectious diseases including more than half of known human pathogens
are not restricted to just one host, yet much of the mathematical modeling of
infections has been limited to a single species. We investigate consequences of
a single epidemic propagating in multiple species and compare and contrast it
with the endemic steady state of the disease. We use the two-species
Susceptible-Infected-Recovered (SIR) model to calculate the severity of
post-epidemic collapses in populations of two host species as a function of
their initial population sizes, the times individuals remain infectious, and
the matrix of infection rates. We derive the criteria for a very large,
extinction-level, population collapse in one or both of the species. The main
conclusion of our study is that a single epidemic could drive a species with
high mortality rate to local or even global extinction provided that it is
co-infected with an abundant species. Such collapse-driven extinctions depend
on factors different than those in the endemic steady state of the disease.


  Precision and reliability of barcode-based biodiversity assessment can be
affected at several steps during acquisition and analysis of the data.
Identification of barcodes is one of the crucial steps in the process and can
be accomplished using several different approaches, namely, alignment-based,
probabilistic, tree-based and phylogeny-based. Number of identified sequences
in the reference databases affects the precision of identification. This paper
compares the identification of marine nematode barcodes using alignment-based,
tree-based and phylogeny-based approaches. Because the nematode reference
dataset is limited in its taxonomic scope, barcodes can only be assigned to
higher taxonomic categories, families. Phylogeny-based approach using
Evolutionary Placement Algorithm provided the largest number of positively
assigned metabarcodes and was least affected by erroneous sequences and
limitations of reference data, comparing to alignment-based and tree-based
approaches.


  The ionic environment of biomolecules strongly influences their structure,
conformational stability, and inter-molecular interactions.This paper
introduces GIBS, a grand-canonical Monte Carlo (GCMC) simulation program for
computing the thermodynamic properties of ion solutions and their distributions
around biomolecules. This software implements algorithms that automate the
excess chemical potential calculations for a given target salt concentration.
GIBS uses a cavity-bias algorithm to achieve high sampling acceptance rates for
inserting ions and solvent hard spheres in simulating dense ionic systems. In
the current version, ion-ion interactions are described using Coulomb,
hard-sphere, or Lennard-Jones (L-J) potentials; solvent-ion interactions are
described using hard-sphere, L-J and attractive square-well potentials; and,
solvent-solvent interactions are described using hard-sphere repulsions. This
paper and the software package includes examples of using GIBS to compute the
ion excess chemical potentials and mean activity coefficients of sodium
chloride as well as to compute the cylindrical radial distribution functions of
monovalent (Na$^+$, Rb$^+$), divalent (Sr$^{2+}$), and trivalent (CoHex$^{3+}$)
around fixed all-atom models of 25 base-pair nucleic acid duplexes. GIBS is
written in C++ and is freely available community use; it can be downloaded at
https://github.com/Electrostatics/GIBS.


  Biochemical reaction networks frequently consist of species evolving on
multiple timescales. Stochastic simulations of such networks are often
computationally challenging and therefore various methods have been developed
to obtain sensible stochastic approximations on the timescale of interest. One
of the rigorous and popular approaches is the multiscale approximation method
for continuous time Markov processes. In this approach, by scaling species
abundances and reaction rates, a family of processes parameterized by a scaling
parameter is defined. The limiting process of this family is then used to
approximate the original process. However, we find that such approximations
become inaccurate when combinations of species with disparate abundances either
constitute conservation laws or form virtual slow auxiliary species. To obtain
more accurate approximation in such cases, we propose here an appropriate
modification of the original method.


  We investigate global stability properties of a HIV/AIDS population model
with constant recruitment rate, mass action incidence, and variable population
size. Existence and uniqueness results for disease-free and endemic equilibrium
points are proved. Global stability of the equilibria is obtained through
Lyapunov's direct method and LaSalle's invariance principle.


  We investigate the dynamics of a greedy forager that moves by random walking
in an environment where each site initially contains one unit of food. Upon
encountering a food-containing site, the forager eats all the food there and
can subsequently hop an additional $\mathcal{S}$ steps without food before
starving to death. Upon encountering an empty site, the forager goes hungry and
comes one time unit closer to starvation. We investigate the new feature of
forager greed; if the forager has a choice between hopping to an empty site or
to a food-containing site in its nearest neighborhood, it hops preferentially
towards food. If the neighboring sites all contain food or are all empty, the
forager hops equiprobably to one of these neighbors. Paradoxically, the
lifetime of the forager can depend non-monotonically on greed, and the sense of
the non-monotonicity is opposite in one and two dimensions. Even more
unexpectedly, the forager lifetime in one dimension is substantially enhanced
when the greed is negative; here the forager tends to avoid food in its local
neighborhood. We also determine the average amount of food consumed at the
instant when the forager starves. We present analytic, heuristic, and numerical
results to elucidate these intriguing phenomena.


  Protein-ligand binding is essential to almost all life processes. The
understanding of protein-ligand interactions is fundamentally important to
rational drug design and protein design. Based on large scale data sets, we
show that protein rigidity strengthening or flexibility reduction is a pivoting
mechanism in protein-ligand binding. Our approach based solely on rigidity is
able to unveil a surprisingly long range contribution of four residue layers to
protein-ligand binding, which has a ramification for drug and protein design.
Additionally, the present work reveals that among various pairwise
interactions, the short range ones within the distance of the van der Waals
diameter are most important. It is found that the present approach outperforms
all the other state-of-the-art scoring functions for protein-ligand binding
affinity predictions of two benchmark data sets


  This paper analyses the dynamics of infectious disease with a concurrent
spread of disease awareness. The model includes local awareness due to contacts
with aware individuals, as well as global awareness due to reported cases of
infection and awareness campaigns. We investigate the effects of time delay in
response of unaware individuals to available information on the epidemic
dynamics by establishing conditions for the Hopf bifurcation of the endemic
steady state of the model. Analytical results are supported by numerical
bifurcation analysis and simulations.


  Ammonium assimilation in E. coli is regulated by two paralogous proteins
(GlnB and GlnK), which orchestrate interactions with regulators of gene
expression, transport proteins and metabolic pathways. Yet how they conjointly
modulate the activity of glutamine synthetase (GS), the key enzyme for nitrogen
assimilation, is poorly understood. We combine experiments and theory to study
the dynamic roles of GlnB and GlnK during nitrogen starvation and upshift. We
measure time-resolved in vivo concentrations of metabolites, total and
post-translationally modified proteins, and develop a concise biochemical model
of GlnB and GlnK that incorporates competition for active and allosteric sites,
as well as functional sequestration of GlnK. The model predicts the responses
of GS, GlnB and GlnK under time-varying external ammonium level in the wild
type and two genetic knock-outs. Our results show that GlnK is tightly
regulated under nitrogen-rich conditions, yet it is expressed during ammonium
run-out and starvation. This suggests a role for GlnK as a buffer of nitrogen
shock after starvation, and provides a further functional link between nitrogen
and carbon metabolisms.


  We used a microfluidic platform to address the problems of obtaining
diffraction quality crystals and crystal handling during transfer to the X-ray
diffractometer. We optimize crystallization conditions of a pharmaceutical
protein and collect X-ray data both in situ and ex situ.


  Boolean networks is a well-established formalism for modelling biological
systems. A vital challenge for analysing a Boolean network is to identify all
the attractors. This becomes more challenging for large asynchronous Boolean
networks, due to the asynchronous updating scheme. Existing methods are
prohibited due to the well-known state-space explosion problem in large Boolean
networks. In this paper, we tackle this challenge by proposing a SCC-based
decomposition method. We prove the correctness of our proposed method and
demonstrate its efficiency with two real-life biological networks.


  Gene regulation is a series of processes that control gene expression and its
extent. The connections among genes and their regulatory molecules, usually
transcription factors, and a descriptive model of such connections, are known
as gene regulatory networks (GRNs). Elucidating GRNs is crucial to understand
the inner workings of the cell and the complexity of gene interactions. To
date, numerous algorithms have been developed to infer gene regulatory
networks. However, as the number of identified genes increases and the
complexity of their interactions is uncovered, networks and their regulatory
mechanisms become cumbersome to test. Furthermore, prodding through
experimental results requires an enormous amount of computation, resulting in
slow data processing. Therefore, new approaches are needed to expeditiously
analyze copious amounts of experimental data resulting from cellular GRNs. To
meet this need, cloud computing is promising as reported in the literature.
Here we propose new MapReduce algorithms for inferring gene regulatory networks
on a Hadoop cluster in a cloud environment. These algorithms employ an
information-theoretic approach to infer GRNs using time-series microarray data.
Experimental results show that our MapReduce program is much faster than an
existing tool while achieving slightly better prediction accuracy than the
existing tool.


  GWAS in humans are revealing the genetic architecture of biomedical and
anthropomorphic traits, i.e., the frequencies and effect sizes of variants that
contribute to heritable variation in a trait. To interpret these findings, we
need to understand how genetic architecture is shaped by basic population
genetics processes-notably, by mutation, natural selection and genetic drift.
Because many quantitative traits are subject to stabilizing selection and
genetic variation that affects one trait often affects many others, we model
the genetic architecture of a focal trait that arises under stabilizing
selection in a multi-dimensional trait space. We solve the model for the
phenotypic distribution and allelic dynamics at steady state and derive robust,
closed form solutions for summary statistics of the genetic architecture. Our
results provide a simple interpretation for missing heritability and why it
varies among traits. They also predict that the distribution of variances
contributed by loci identified in GWAS is well approximated by a simple
functional form that depends on a single parameter: the expected contribution
to genetic variance of a strongly selected site affecting the trait. We test
this prediction against the results of GWAS for height and BMI and find that it
fits the data well, allowing us to make inferences about the degree of
pleiotropy and mutational target size for these traits. Our findings help to
explain why the GWAS for height explains more of the heritable variance than
similarly-sized GWAS for BMI, and to predict the increase in explained
heritability with study sample size. Considering the demographic history of
European populations, in which these GWAS were performed, we further find that
most of the associations they identified likely involve mutations that arose
during the out of Africa bottleneck at sites with selection coefficients around
$s=10^{-3}$.


  This paper considers a broadly biologically relevant question of a chain
(such as a protein) binding to a sequence of receptors with matching multiple
ligands distributed along the chain. This binding is critical in cell adhesion
events, and in protein self-assembly. Using a mean field approximation of
polymer dynamics, we first calculate the characteristic binding time for a
tethered ligand reaching for a specific binding site on the surface. This time
is determined by two separate entropic effects: an entropic barrier for the
chain to be stretched sufficiently to reach the distant target, and a
restriction on chain conformations near the surface. We then derive the
characteristic time for a sequence of single binding events, and find that it
is determined by the `zipper effect', optimizing the sequence of single and
multiple binding steps.


  Species tree reconstruction from genomic data is increasingly performed using
methods that account for sources of gene tree discordance such as incomplete
lineage sorting. One popular method for reconstructing species trees from
unrooted gene tree topologies is ASTRAL. In this paper, we derive theoretical
sample complexity results for the number of genes required by ASTRAL to
guarantee reconstruction of the correct species tree with high probability. We
also validate those theoretical bounds in a simulation study. Our results
indicate that ASTRAL requires $\mathcal{O}(f^{-2} \log n)$ gene trees to
reconstruct the species tree correctly with high probability where n is the
number of species and f is the length of the shortest branch in the species
tree. Our simulations, which are the first to test ASTRAL explicitly under the
anomaly zone, show trends consistent with the theoretical bounds and also
provide some practical insights on the conditions where ASTRAL works well.


  In recent years there has been a growing interest in the study of the
dynamics of stochastic populations. A key question in population biology is to
understand the conditions under which populations coexist or go extinct.
Theoretical and empirical studies have shown that coexistence can be
facilitated or negated by both biotic interactions and environmental
fluctuations. We study the dynamics of $n$ populations that live in a
stochastic environment and which can interact nonlinearly (through competition
for resources, predator-prey behavior, etc.). Our models are described by
$n$-dimensional Kolmogorov systems with white noise (stochastic differential
equations - SDE). We give sharp conditions under which the populations converge
exponentially fast to their unique stationary distribution as well as
conditions under which some populations go extinct exponentially fast.
  The analysis is done by a careful study of the properties of the invariant
measures of the process that are supported on the boundary of the domain. To
our knowledge this is one of the first general results describing the
asymptotic behavior of stochastic Kolmogorov systems in non-compact domains.
  We are able to fully describe the properties of many of the SDE that appear
in the literature. In particular, we extend results on two dimensional
Lotka-Volterra models, two dimensional predator-prey models, $n$ dimensional
simple food chains, and two predator and one prey models. We also show how one
can use our methods to classify the dynamics of any two-dimensional stochastic
Kolmogorov system satisfying some mild assumptions.


  Signaling in enzymatic networks is typically triggered by environmental
fluctuations, resulting in a series of stochastic chemical reactions, leading
to corruption of the signal by noise. For example, information flow is
initiated by binding of extracellular ligands to receptors, which is
transmitted through a {cascade involving} kinase-phosphatase stochastic
chemical reactions. For a class of such networks, we develop a general
field-theoretic approach in order to calculate the error in signal transmission
as a function of an appropriate control variable. Application of the theory to
a simple push-pull network, a module in the kinase-phosphatase cascade,
recovers the exact results for error in signal transmission previously obtained
using umbral calculus (Phys. Rev. X., {\bf 4}, 041017 (2014)). We illustrate
the generality of the theory by studying the minimal errors in noise reduction
in a reaction cascade with two connected push-pull modules. Such a cascade
behaves as an effective {three-species} network with a pseudo intermediate. In
this case, optimal information transfer, resulting in the smallest square of
the error between the input and output, occurs with a time delay, which is
given by the inverse of the decay rate of the pseudo intermediate. There are
substantial deviations between simulations and predictions of the linear theory
in error in signal propagation in an enzymatic push-pull network for a certain
range of parameters. Inclusion of second order perturbative corrections shows
that differences between simulations and theoretical predictions are minimized.
Our study establishes that {a} field theoretic formulation {of} stochastic
biological {signaling offers} a systematic way to understand error propagation
in networks of arbitrary complexity.


  Transforming error-prone immunosequencing datasets into antibody repertoires
is a fundamental problem in immunogenomics, and a prerequisite for studies of
immune responses. Although various repertoire reconstruction algorithms were
released in the last three years, it remains unclear how to benchmark them and
how to assess the accuracy of the reconstructed repertoires. We describe a
novel IgReC algorithm for constructing antibody repertoires from
high-throughput immunosequencing datasets and a new framework for assessing the
quality of reconstructed repertoires. Benchmarking IgReC against the existing
antibody repertoire reconstruction tools has demonstrated that it results in
highly accurate repertoire reconstructions. Surprisingly, antibody repertoires
constructed by IgReC from barcoded immunosequencing datasets in blind mode
(without using unique molecular identifiers information) improved upon the
repertoires constructed by the state-of-the-art tools that use barcoding. This
finding suggests that IgReC may alleviate the need to generate repertoires
using the barcoding technology (the workhorse of current immunogenomics
efforts) because our computational approach to error correction of
immunosequencing data ends up being nearly as powerful as the experimental
approach based on barcoding.


  Computational prediction of membrane protein (MP) structures is very
challenging partially due to lack of sufficient solved structures for homology
modeling. Recently direct evolutionary coupling analysis (DCA) sheds some light
on protein contact prediction and accordingly, contact-assisted folding, but
DCA is effective only on some very large-sized families since it uses
information only in a single protein family. This paper presents a deep
transfer learning method that can significantly improve MP contact prediction
by learning contact patterns and complex sequence-contact relationship from
thousands of non-membrane proteins (non-MPs). Tested on 510 non-redundant MPs,
our deep model (learned from only non-MPs) has top L/10 long-range contact
prediction accuracy 0.69, better than our deep model trained by only MPs (0.63)
and much better than a representative DCA method CCMpred (0.47) and the CASP11
winner MetaPSICOV (0.55). The accuracy of our deep model can be further
improved to 0.72 when trained by a mix of non-MPs and MPs. When only contacts
in transmembrane regions are evaluated, our method has top L/10 long-range
accuracy 0.62, 0.57, and 0.53 when trained by a mix of non-MPs and MPs, by
non-MPs only, and by MPs only, respectively, still much better than MetaPSICOV
(0.45) and CCMpred (0.40). All these results suggest that sequence-structure
relationship learned by our deep model from non-MPs generalizes well to MP
contact prediction. Improved contact prediction also leads to better
contact-assisted folding. Using only top predicted contacts as restraints, our
deep learning method can fold 160 and 200 of 510 MPs with TMscore>0.6 when
trained by non-MPs only and by a mix of non-MPs and MPs, respectively, while
CCMpred and MetaPSICOV can do so for only 56 and 77 MPs, respectively. Our
contact-assisted folding also greatly outperforms homology modeling.


  We study Eigen's quasispecies model in the asymptotic regime where the length
of the genotypes goes to infinity and the mutation probability goes to 0. A
limiting infinite system of differential equations is obtained. We prove the
convergence of the trajectories, as well as the convergence of the equilibrium
solutions. We give the analogous results for a discrete-time version of Eigen's
model, which coincides with a model proposed by Moran.


  Given a conjunctive Boolean network (CBN) with $n$ state-variables, we
consider the problem of finding a minimal set of state-variables to directly
affect with an input so that the resulting conjunctive Boolean control network
(CBCN) is controllable. We give a necessary and sufficient condition for
controllability of a CBCN; an $O(n^2)$-time algorithm for testing
controllability; and prove that nonetheless the minimal controllability problem
for CBNs is NP-hard.


  This paper proposes an approach to environmental accounting useful for
studying the feasibility of socio-economic systems in relation to the external
constraints posed by ecological compatibility. The approach is based on a
multi-scale analysis of the metabolic pattern of ecosystems and societies and
it provides an integrated characterization of the resulting interaction. The
text starts with a theoretical part explaining (i) the implicit epistemological
revolution implied by the notion of ecosystem metabolism and the fund-flow
model developed by Georgescu-Roegen applied to environmental accounting, and
(ii) the potentials of this approach to create indicators to assess ecological
integrity and environmental impacts. This revolution also makes it possible to
carry out a multi-scale integrated assessment of ecosystem and societal
metabolisms at the territorial level. In the second part, two applications of
this approach using an indicator of the negentropic cost show the possibility
to characterize in quantitative and qualitative terms degrees of alteration
(crop cultivation, tree plantations)for different biomes (tropical and boreal
forests). Also, a case study for land use scenarios has been included. The
proposed approach represents an integrated multi-scale tool for the analysis of
nature conservation scenarios and strategies.


  This paper is devoted to the analysis of a simple Lotka-Volterra food chain
evolving in a stochastic environment. It can be seen as the companion paper of
Hening and Nguyen (J. of Math. Biol. `18) where we have characterized the
persistence and extinction of such a food chain under the assumption that there
is no intraspecific competition among predators. In the current paper we focus
on the case when all the species experience intracompetition. The food chain we
analyze consists of one prey and $n-1$ predators. The $j$th predator eats the
$j-1$st species and is eaten by the $j+1$st predator; this way each species
only interacts with at most two other species - the ones that are immediately
above or below it in the trophic chain. We show that one can classify, based on
the invasion rates of the predators (which we can determine from the
interaction coefficients of the system via an algorithm), which species go
extinct and which converge to their unique invariant probability measure. We
obtain stronger results than in the case with no intraspecific competition
because in this setting we can make use of the general results of Hening and
Nguyen (Ann. of Appl. Probab.). Unlike most of the results available in the
literature, we provide an in depth analysis for both non-degenerate and
degenerate noise.
  We exhibit our general results by analysing trophic cascades in a
plant--herbivore--predator system and providing persistence/extinction criteria
for food chains of length $n\leq 3$.


  The relationship between the M-species stochastic Lotka-Volterra competition
(SLVC) model and the M-allele Moran model of population genetics is explored
via timescale separation arguments. When selection for species is weak and the
population size is large but finite, precise conditions are determined for the
stochastic dynamics of the SLVC model to be mappable to the neutral Moran
model, the Moran model with frequency-independent selection and the Moran model
with frequency-dependent selection (equivalently, a game-theoretic formulation
of the Moran model). We demonstrate how these mappings can be used to calculate
extinction probabilities and the times until a species' extinction in the SLVC
model.


  Model-based prediction of stochastic noise in biomolecular reactions often
resorts to approximation with unknown precision. As a result, unexpected
stochastic fluctuation causes a headache for the designers of biomolecular
circuits. This paper proposes a convex optimization approach to quantifying the
steady state moments of molecular copy counts with theoretical rigor. We show
that the stochastic moments lie in a convex semi-algebraic set specified by
linear matrix inequalities. Thus, the upper and the lower bounds of some
moments can be computed by a semidefinite program. Using a protein dimerization
process as an example, we demonstrate that the proposed method can precisely
predict the mean and the variance of the copy number of the monomer protein.


  The origin of population-scale coordination has puzzled philosophers and
scientists for centuries. Recently, game theory, evolutionary approaches and
complex systems science have provided quantitative insights on the mechanisms
of social consensus. However, the literature is vast and widely scattered
across fields, making it hard for the single researcher to navigate it. This
short review aims to provide a compact overview of the main dimensions over
which the debate has unfolded and to discuss some representative examples. It
focuses on those situations in which consensus emerges 'spontaneously' in
absence of centralised institutions and covers topic that include the
macroscopic consequences of the different microscopic rules of behavioural
contagion, the role of social networks, and the mechanisms that prevent the
formation of a consensus or alter it after it has emerged. Special attention is
devoted to the recent wave of experiments on the emergence of consensus in
social systems.


  We present a method of detecting sequence defects by supercoiling DNA with
magnetic tweezers. The method is sensitive to a single mismatched base pair in
a DNA sequence of several thousand base pairs. We systematically compare DNA
molecules with 0 to 16 adjacent mismatches at 1 M monovalent salt and 3.5 pN
force and show that, under these conditions, a single plectoneme forms and is
stably pinned at the defect. We use these measurements to estimate the energy
and degree of end-loop kinking at defects. From this, we calculate the relative
probability of plectoneme pinning at the mismatch under physiologically
relevant conditions. Based on this estimate, we propose that DNA supercoiling
could contribute to mismatch and damage sensing in vivo.


  Containing the recent West African outbreak of Ebola virus (EBOV) required
the deployment of substantial global resources. Operationally, health workers
and surveillance teams treated cases, collected genetic samples, and tracked
case contacts. Despite the substantial progress in analyzing and modeling EBOV
epidemiological data, a complete characterization of the spatiotemporal spread
of Ebola cases remains a challenge. In this work, we offer a novel perspective
on the EBOV epidemic that utilizes virus genome sequences to inform
population-level, spatial models. Calibrated to phylogenetic linkages, these
dynamic spatial models provide unique insight into the disease mobility of EBOV
in Sierra Leone. Further, we developed a model selection framework that
identifies important epidemiological variables influencing the spatiotemporal
propagation of EBOV. Consistent with other investigations, our results show
that the spread of EBOV during the beginning and middle portions of the
epidemic strongly depended on the size of and distance between populations. Our
analysis also revealed a substantial decline in the dependence on population
size at the end of the epidemic, coinciding with the large-scale intervention
campaign: Operation Western Area Surge. More generally, we believe this
framework, pairing molecular diagnostics with dynamic models, has the potential
to be a powerful forecasting tool along with offering operationally-relevant
guidance for surveillance and sampling strategies during an epidemic.


  The missing phase problem in X-ray crystallography is commonly solved using
the technique of molecular replacement, which borrows phases from a previously
solved homologous structure, and appends them to the measured Fourier
magnitudes of the diffraction patterns of the unknown structure. More recently,
molecular replacement has been proposed for solving the missing orthogonal
matrices problem arising in Kam's autocorrelation analysis for single particle
reconstruction using X-ray free electron lasers and cryo-EM. In classical
molecular replacement, it is common to estimate the magnitudes of the unknown
structure as twice the measured magnitudes minus the magnitudes of the
homologous structure, a procedure known as `twicing'. Mathematically, this is
equivalent to finding an unbiased estimator for a complex-valued scalar. We
generalize this scheme for the case of estimating real or complex valued
matrices arising in single particle autocorrelation analysis. We name this
approach "Anisotropic Twicing" because unlike the scalar case, the unbiased
estimator is not obtained by a simple magnitude isotropic correction. We
compare the performance of the least squares, twicing and anisotropic twicing
estimators on synthetic and experimental datasets. We demonstrate 3D homology
modeling in cryo-EM directly from experimental data without iterative
refinement or class averaging, for the first time.


  We consider the problem of inferring the probability distribution of flux
configurations in metabolic network models from empirical flux data. For the
simple case in which experimental averages are to be retrieved, data are
described by a Boltzmann-like distribution ($\propto e^{F/T}$) where $F$ is a
linear combination of fluxes and the `temperature' parameter $T\geq 0$ allows
for fluctuations. The zero-temperature limit corresponds to a Flux Balance
Analysis scenario, where an objective function ($F$) is maximized. As a test,
we have inverse modeled, by means of Boltzmann learning, the catabolic core of
Escherichia coli in glucose-limited aerobic stationary growth conditions.
Empirical means are best reproduced when $F$ is a simple combination of biomass
production and glucose uptake and the temperature is finite, implying the
presence of fluctuations. The scheme presented here has the potential to
deliver new quantitative insight on cellular metabolism. Our implementation is
however computationally intensive, and highlights the major role that effective
algorithms to sample the high-dimensional solution space of metabolic networks
can play in this field.


  The use of computers in statistical physics is common because the sheer
number of equations that describe the behavior of an entire system particle by
particle often makes it impossible to solve them exactly. Monte Carlo methods
form a particularly important class of numerical methods for solving problems
in statistical physics. Although these methods are simple in principle, their
proper use requires a good command of statistical mechanics, as well as
considerable computational resources. The aim of this paper is to demonstrate
how the usage of widely accessible graphics cards on personal computers can
elevate the computing power in Monte Carlo simulations by orders of magnitude,
thus allowing live classroom demonstration of phenomena that would otherwise be
out of reach. As an example, we use the public goods game on a square lattice
where two strategies compete for common resources in a social dilemma
situation. We show that the second-order phase transition to an absorbing phase
in the system belongs to the directed percolation universality class, and we
compare the time needed to arrive at this result by means of the main processor
and by means of a suitable graphics card. Parallel computing on graphics
processing units has been developed actively during the last decade, to the
point where today the learning curve for entry is anything but steep for those
familiar with programming. The subject is thus ripe for inclusion in graduate
and advanced undergraduate curricula, and we hope that this paper will
facilitate this process in the realm of physics education. To that end, we
provide a documented source code for an easy reproduction of presented results
and for further development of Monte Carlo simulations of similar systems.


  Metal ions, including copper and zinc, have been implicated in the
pathogenesis of Alzheimers disease through a variety of mechanisms including
increased amyloid \b{eta} affinity and redox effects. Recent reports have
demonstrated that the amyloid \b{eta} monomer does not necessarily travel
through a definitive intermediary en-route to a stable amyloid fibril
structure. Rather, amyloid \b{eta} misfolding may follow a variety of pathways
resulting in a fibrillar end-product or a variety of oligomeric end-products
with a diversity of structures and sizes. The presence of metal ions has been
demonstrated to alter the kinetic pathway of the amyloid \b{eta} peptide which
may lead to more toxic oligomeric end-products. In this work, we review the
contemporary literature supporting the hypothesis that metal ions alter the
reaction pathway of amyloid \b{eta} misfolding leading to more neurotoxic
species.


  Using atomic force microscopy (AFM) we investigated the interaction of
amyloid beta (Ab) (1 42) peptide with chemically modified surfaces in order to
better understand the mechanism of amyloid toxicity, which involves interaction
of amyloid with cell membrane surfaces. We compared the structure and density
of Ab fibrils on positively and negatively charged as well as hydrophobic
chemically modified surfaces at physiologically relevant conditions.


  Proposal to develop an Improved immunological assay employing primary IgG
antibodies and secondary IgM antibodies labeled with quantum dots to amplify
antigen detection.


  A central goal in cancer genomics is to identify the somatic alterations that
underpin tumor initiation and progression. This task is challenging as the
mutational profiles of cancer genomes exhibit vast heterogeneity, with many
alterations observed within each individual, few shared somatically mutated
genes across individuals, and important roles in cancer for both frequently and
infrequently mutated genes. While commonly mutated cancer genes are readily
identifiable, those that are rarely mutated across samples are difficult to
distinguish from the large numbers of other infrequently mutated genes. Here,
we introduce a method that considers per-individual mutational profiles within
the context of protein-protein interaction networks in order to identify small
connected subnetworks of genes that, while not individually frequently mutated,
comprise pathways that are perturbed across (i.e., "cover") a large fraction of
the individuals. We devise a simple yet intuitive objective function that
balances identifying a small subset of genes with covering a large fraction of
individuals. We show how to solve this problem optimally using integer linear
programming and also give a fast heuristic algorithm that works well in
practice. We perform a large-scale evaluation of our resulting method, nCOP, on
6,038 TCGA tumor samples across 24 different cancer types. We demonstrate that
our approach nCOP is more effective in identifying cancer genes than both
methods that do not utilize any network information as well as state-of-the-art
network-based methods that aggregate mutational information across individuals.
Overall, our work demonstrates the power of combining per-individual mutational
information with interaction networks in order to uncover genes functionally
relevant in cancers, and in particular those genes that are less frequently
mutated.


  This work introduces a new statistical physics lattice model of bacteria
interacting with anti-microbial drugs that can reproduce qualitative features
of resistance emergence and whose model parameters and outputs can be measured
with controlled \textit{in vitro} experiments. The lattice is inhabited by
agents modeled by Ising perceptrons. The results show the advantage of mixing
drugs among the population compared to other treatment protocols.


  We propose new activity-dependent adaptive Boolean networks inspired by the
cis-regulatory mechanism in gene regulatory networks. We analytically show that
our model can be solved for stationary in-degree distribution for a wide class
of update rules by employing the annealed approximation of Boolean network
dynamics and that evolved Boolean networks have a preassigned average
sensitivity that can be set independently of update rules if certain conditions
are satisfied. In particular, when it is set to 1, our theory predicts that the
proposed network rewiring algorithm drives Boolean networks towards
criticality. We verify that these analytic results agree well with numerical
simulations for four representative update rules. We also discuss the
relationship between sensitivity of update rules and stationary in-degree
distributions and compare it with that in real-world gene regulatory networks.


  Why do some females menstruate at all? Answering this question has
implications for understanding the tight links between reproductive function
and organismal immunity. Here we build on the growing evidence that
menstruation is the by-product of a "choosy" uterus to (i) make the theoretical
case for the idea that female immunity is cyclical in menstruating species,
(ii) evaluate the evidence for the menstrual modulation of immunity and health
in humans and (iii) speculate on the implications of cyclical female health for
female behaviour, male immunity and host-pathogen interactions. We argue that
an understanding of females' evolved reproductive system is foundational for
both tackling the future challenges of the global women's health agenda and
predicting eco-evolutionary dynamics in cyclically reproducing species.


  A computation-oriented representation of uncertain kinetic systems is
introduced and analysed in this paper. It is assumed that the monomial
coefficients of the ODEs belong to a polytopic set, which defines a set of
dynamical systems for an uncertain model. An optimization-based computation
model is proposed for the structural analysis of uncertain models. It is shown
that the so-called dense realization containing the maximum number of reactions
(directed edges) is computable in polynomial time, and it forms a
super-structure among all the possible reaction graphs corresponding to an
uncertain kinetic model, assuming a fixed set of complexes. The set of core
reactions present in all reaction graphs of an uncertain model is also studied.
Most importantly, an algorithm is proposed to compute all possible reaction
graph structures for an uncertain kinetic model.


  Since there is now a growing wish by referees to judge the underpinning data
for a submitted article it is timely to provide a summary of the data
evaluation checks required to be done by a referee. As these checks will vary
from field to field this article focuses on the needs of biological X-ray
crystallography articles, which is the predominantly used method leading to
depositions in the PDB. The expected referee checks of data underpinning an
article are described with examples. These checks necessarily include that a
referee checks the PDB validation report for each crystal structure
accompanying the article submission; this check whilst necessary is not
sufficient for a complete evaluation. A referee would be expected to undertake
one cycle of model refinement of the authors biological macromolecule
coordinates against the authors processed diffraction data and look at the
various validation checks of the model and Fo-Fc electron density maps in e.g.
Phenix_refine and in COOT. If the referee deems necessary the diffraction data
images should be reprocessed (e.g. to a different diffraction resolution than
the authors submission). This can be requested to be done by the authors or if
the referee prefers can be undertaken directly by the referee themselves. A
referee wishing to do these data checks may wish to receive a certificate that
they have command of these data science skills. The organisation of such
voluntary certification training can e.g. be via those crystallography
associations duly recognised by the IUCr to issue such certificates.


  The analysis of the mutational landscape of cancer, including mutual
exclusivity and co-occurrence of mutations, has been instrumental in studying
the disease. We hypothesized that exploring the interplay between
co-occurrence, mutual exclusivity, and functional interactions between genes
will further improve our understanding of the disease and help to uncover new
relations between cancer driving genes and pathways. To this end, we designed a
general framework, BeWith, for identifying modules with different combinations
of mutation and interaction patterns. We focused on three different settings of
the BeWith schema: (i) BeME-WithFun in which the relations between modules are
enriched with mutual exclusivity while genes within each module are
functionally related; (ii) BeME-WithCo which combines mutual exclusivity
between modules with co-occurrence within modules; and (iii) BeCo-WithMEFun
which ensures co-occurrence between modules while the within module relations
combine mutual exclusivity and functional interactions. We formulated the
BeWith framework using Integer Linear Programming (ILP), enabling us to find
optimally scoring sets of modules. Our results demonstrate the utility of
BeWith in providing novel information about mutational patterns, driver genes,
and pathways. In particular, BeME-WithFun helped identify functionally coherent
modules that might be relevant for cancer progression. In addition to finding
previously well-known drivers, the identified modules pointed to the importance
of the interaction between NCOR and NCOA3 in breast cancer. Additionally, an
application of the BeME-WithCo setting revealed that gene groups differ with
respect to their vulnerability to different mutagenic processes, and helped us
to uncover pairs of genes with potentially synergetic effects, including a
potential synergy between mutations in TP53 and metastasis related DCC gene.


  One of the causes of high fidelity of copying in biological systems is
kinetic discrimination. In this mechanism larger dissipation and copying
velocity result in improved copying accuracy. We consider a model of a
polymerase which simultaneously copies a single stranded RNA and opens a
single- to double-stranded junction serving as an obstacle. The presence of the
obstacle slows down the motor, resulting in a change of its fidelity, which can
be used to gain information about the motor and junction dynamics. We find that
the motor's fidelity does not depend on details of the motor-junction
interaction, such as whether the interaction is passive or active. Analysis of
the copying fidelity can still be used as a tool for investigating the junction
kinetics.


  The correct prediction of protein secondary structures is one of the key
issues in predicting the correct protein folded shape, which is used for
determining gene function. Existing methods make use of amino acids properties
as indices to classify protein secondary structures, but are faced with a
significant number of misclassifications. The paper presents a technique for
the classification of protein secondary structures based on protein
"signal-plotting" and the use of the Fourier technique for digital signal
processing. New indices are proposed to classify protein secondary structures
by analyzing hydrophobicity profiles. The approach is simple and
straightforward. Results show that the more types of protein secondary
structures can be classified by means of these newly-proposed indices.


  We present a conditional generative model to learn variation in cell and
nuclear morphology and the location of subcellular structures from microscopy
images. Our model generalizes to a wide range of subcellular localization and
allows for a probabilistic interpretation of cell and nuclear morphology and
structure localization from fluorescence images. We demonstrate the
effectiveness of our approach by producing photo-realistic cell images using
our generative model. The conditional nature of the model provides the ability
to predict the localization of unobserved structures given cell and nuclear
morphology.