  Engineered biomedical nanoparticles (NP) administered via intravenous routes
are prone to associate to serum proteins. The protein corona can mask the NP
surface functionalization and hamper the delivery of the NP to its biological
target. The design of corona-free NPs relies on our understanding of the
chemical-physical features of the NP surface driving the interaction with serum
proteins. Here we address, by computational means, the interaction between
human serum albumin (HSA) and a prototypical monolayer-protected Au
nanoparticle. We show that both the chemical composition (charge,
hydrophobicity) and the conformational preferences of the ligands decorating
the NP surface affect the NP propensity to bind HSA.


  Two types of osmolytes, i.e., trimethylamin N-oxide (TMAO) and urea,
demonstrate dramatically different properties in a protein folding process.
Even with the great progresses in revealing the potential underlying mechanism
of these two osmolyte systems, many problems still remain unsolved. In this
paper, we propose to use the persistent homology, a newly-invented topological
method, to systematically study the osmolytes molecular aggregation and their
hydrogen-bonding network from a global topological perspective. It has been
found that, for the first time, TMAO and urea show two extremely different
topological behaviors, i.e., extensive network and local cluster. In general,
TMAO forms highly consistent large loop or circle structures in high
concentrations. In contrast, urea is more tightly aggregated locally. Moreover,
the resulting hydrogen-bonding networks also demonstrate distinguishable
features. With the concentration increase, TMAO hydrogen-bonding networks vary
greatly in their total number of loop structures and large-sized loop
structures consistently increase. In contrast, urea hydrogen-bonding networks
remain relatively stable with slight reduce of the total loop number. Moreover,
the persistent entropy (PE) is, for the first time, used in characterization of
the topological information of the aggregation and hydrogen-bonding networks.
The average PE systematically increases with the concentration for both TMAO
and urea, and decreases in their hydrogen-bonding networks. But their PE
variances have totally different behaviors. Finally, topological features of
the hydrogen-bonding networks are found to be highly consistent with those from
the ion aggregation systems, indicating that our topological invariants can
characterize intrinsic features of the "structure making" and "structure
breaking" systems.


  Statistical models for families of evolutionary related proteins have
recently gained interest: in particular pairwise Potts models, as those
inferred by the Direct-Coupling Analysis, have been able to extract information
about the three-dimensional structure of folded proteins, and about the effect
of amino-acid substitutions in proteins. These models are typically requested
to reproduce the one- and two-point statistics of the amino-acid usage in a
protein family, {\em i.e.}~to capture the so-called residue conservation and
covariation statistics of proteins of common evolutionary origin. Pairwise
Potts models are the maximum-entropy models achieving this. While being
successful, these models depend on huge numbers of {\em ad hoc} introduced
parameters, which have to be estimated from finite amount of data and whose
biophysical interpretation remains unclear. Here we propose an approach to
parameter reduction, which is based on selecting collective sequence motifs. It
naturally leads to the formulation of statistical sequence models in terms of
Hopfield-Potts models. These models can be accurately inferred using a mapping
to restricted Boltzmann machines and persistent contrastive divergence. We show
that, when applied to protein data, even 20-40 patterns are sufficient to
obtain statistically close-to-generative models. The Hopfield patterns form
interpretable sequence motifs and may be used to clusterize amino-acid
sequences into functional sub-families. However, the distributed collective
nature of these motifs intrinsically limits the ability of Hopfield-Potts
models in predicting contact maps, showing the necessity of developing models
going beyond the Hopfield-Potts models discussed here.


  Calculation of the log-likelihood stands as the computational bottleneck for
many statistical phylogenetic algorithms. Even worse is its gradient
evaluation, often used to target regions of high probability. Order ${\cal
O}\hspace{-0.2em}\left( N \right)$-dimensional gradient calculations based on
the standard pruning algorithm require ${\cal O}\hspace{-0.2em}\left( N^2
\right)$ operations where N is the number of sampled molecular sequences. With
the advent of high-throughput sequencing, recent phylogenetic studies have
analyzed hundreds to thousands of sequences, with an apparent trend towards
even larger data sets as a result of advancing technology. Such large-scale
analyses challenge phylogenetic reconstruction by requiring inference on larger
sets of process parameters to model the increasing data heterogeneity. To make
this tractable, we present a linear-time algorithm for ${\cal
O}\hspace{-0.2em}\left( N \right)$-dimensional gradient evaluation and apply it
to general continuous-time Markov processes of sequence substitution on a
phylogenetic tree without a need to assume either stationarity or
reversibility. We apply this approach to learn the branch-specific evolutionary
rates of three pathogenic viruses: West Nile virus, Dengue virus and Lassa
virus. Our proposed algorithm significantly improves inference efficiency with
a 126- to 234-fold increase in maximum-likelihood optimization and a 16- to
33-fold computational performance increase in a Bayesian framework.


  Naive human T cells are produced in the thymus, which atrophies abruptly and
severely in response to physical or psychological stress. To understand how an
instance of stress affects the size and "diversity" of the peripheral naive T
cell pool, we derive a mean-field autonomous ODE model of T cell replenishment
that allows us to track the clone abundance distribution (the mean number of
different TCRs each represented by a specific number of cells). We identify
equilibrium solutions that arise at different rates of T cell production, and
derive analytic approximations to the dominant eigenvalues and eigenvectors of
the problem linearized about these equilibria. From the forms of the
eigenvalues and eigenvectors, we estimate rates at which counts of clones of
different sizes converge to and depart from equilibrium values--that is, how
the number of clones of different sizes "adjust" to the changing rate of T cell
production. Under most physiologically realistic realizations of our model, the
dominant eigenvalue (representing the slowest dynamics of the clone abundance
distribution) scales as a power law in the thymic output for low output levels,
but saturates at higher T cell production rates. Our analysis provides a
framework for quantitatively understanding how the clone abundance
distributions evolve under small changes in the overall T cell production rate
by the thymus.


  We study a simple model of DNA evolution in a growing population of cells.
Each cell contains a nucleotide sequence which randomly mutates at cell
division. Cells divide according to a branching process. Following typical
parameter values in bacteria and cancer cell populations, we take the mutation
rate to zero and the final number of cells to infinity. We prove that almost
every site (entry of the nucleotide sequence) is mutated in only a finite
number of cells, and these numbers are independent across sites. However
independence breaks down for the rare sites which are mutated in a positive
fraction of the population. The model is free from the popular but disputed
infinite sites assumption. Violations of the infinite sites assumption are
widespread while their impact on mutation frequencies is negligible at the
scale of population fractions. Some results are generalised to allow for cell
death, selection, and site-specific mutation rates. For illustration we
estimate mutation rates in a lung adenocarcinoma.


  Upon contact with antigen presenting cells (APCs) cytotoxic T lymphocytes (T
cells) establish a highly organized contact zone denoted as immunological
synapse (IS). The formation of the IS implies relocation of the microtubule
organizing center (MTOC) towards the contact zone, which necessitates a proper
connection between MTOC and IS via dynamic microtubules (MTs). The efficiency
of the MTs finding the IS within relevant time scale is, however, still
illusive. We investigate how MTs search the three-dimensional constrained
cellular volume for the IS and bind upon encounter to dynein anchored at the IS
cortex. The search efficiency is estimated by calculating the time required for
the MTs to reach the dynein-enriched region of the IS. In this study, we
develop simple mathematical and numerical models incorporating relevant
components of a cell and propose an optimal search strategy. Using the
mathematical model, we have quantified the average search time for a wide range
of model parameters and proposed an optimized set of values leading to the
minimum capture time. Our results show that search times are minimal when the
IS formed at the nearest or at the farthest sites on the cell surface, with
respect to the perinuclear MTOC. The search time increases monotonically away
from these two specific sites and are maximal at an intermediate position near
the equator of the cell. We observed that search time strongly depends on the
number of searching MTs and distance of the MTOC from the nuclear surface.


  Synthesis of biopolymers such as DNA, RNA, and proteins are biophysical
processes aided by enzymes. Performance of these enzymes is usually
characterized in terms of their average error rate and speed. However, because
of thermal fluctuations in these single-molecule processes, both error and
speed are inherently stochastic quantities. In this paper, we study
fluctuations of error and speed in biopolymer synthesis and show that they are
in general correlated. This means that, under equal conditions, polymers that
are synthesized faster due to a fluctuation tend to have either better or worse
errors than the average. The error-correction mechanism implemented by the
enzyme determines which of the two cases holds. For example, discrimination in
the forward reaction rates tends to grant smaller errors to polymers with
faster synthesis. The opposite occurs for discrimination in monomer rejection
rates. Our results provide an experimentally feasible way to identify
error-correction mechanisms by measuring the error-speed correlations.


  Systems biology and whole-cell modelling are demanding increasingly
comprehensive mathematical models of cellular biochemistry. These models
require the development of simplified models of specific processes which
capture essential biophysical features but without unnecessarily complexity.
Recently there has been renewed interest in thermodynamically-based modelling
of cellular processes. Here we present an approach to developing of simplified
yet thermodynamically consistent (hence physically plausible) models which can
readily be incorporated into large scale biochemical descriptions but which do
not require full mechanistic detail of the underlying processes. We illustrate
the approach through development of a simplified, physically plausible model of
the mitochondrial electron transport chain and show that the simplified model
behaves like the full system.


  Designing new molecules with a set of predefined properties is a core problem
in modern drug discovery and development. There is a growing need for de-novo
design methods that would address this problem. We present MolecularRNN, the
graph recurrent generative model for molecular structures. Our model generates
diverse realistic molecular graphs after likelihood pretraining on a big
database of molecules. We perform an analysis of our pretrained models on
large-scale generated datasets of 1 million samples. Further, the model is
tuned with policy gradient algorithm, provided a critic that estimates the
reward for the property of interest. We show a significant distribution shift
to the desired range for lipophilicity, drug-likeness, and melting point
outperforming state-of-the-art works. With the use of rejection sampling based
on valency constraints, our model yields 100% validity. Moreover, we show that
invalid molecules provide a rich signal to the model through the use of
structure penalty in our reinforcement learning pipeline.


  The translocation of double-stranded DNA through a solid-state nanopore may
either decrease or increase the ionic current depending on the ionic
concentration of the surrounding solution. Below a certain crossover ionic
concentration, the current change inverts from a current blockade to current
enhancement. In this paper, we show that the crossover concentration for
bundled DNA nanostructures composed of multiple connected DNA double-helices is
lower than that of double-stranded DNA. Our measurements suggest that
counterion mobility in the vicinity of DNA is reduced depending on the
three-dimensional structure of the molecule. We further demonstrate that
introducing neutral polymers such as polyethylene glycol into the measurement
solution reduces electroosmotic outflow from the nanopore, allowing
translocation of large DNA structures at low salt concentrations. Our
experiments contribute to an improved understanding of ion transport in
confined DNA environments, which is critical for the development of nanopore
sensing techniques as well as synthetic membrane channels. Our salt-dependent
measurements of model DNA nanostructures will guide the development of
computational models of DNA translocation through nanopores.


  Boolean networks (BNs) are widely used to model the qualitative dynamics of
biological systems. Besides the logical rules determining the evolution of each
component with respect to the state of its regulators, the scheduling of
component updates can have a dramatic impact on the predicted behaviours. In
this paper, we explore the use of Read (contextual) Petri Nets (RPNs) to study
dynamics of BNs from a concurrency theory perspective. After showing
bi-directional translations between RPNs and BNs and analogies between results
on synchronism sensitivity, we illustrate that usual updating modes for BNs can
miss plausible behaviours, i.e., incorrectly conclude on the
absence/impossibility of reaching specific configurations. We propose an
encoding of BNs capitalizing on the RPN semantics enabling more behaviour than
the generalized asynchronous updating mode. The proposed encoding ensures a
correct abstraction of any multivalued refinement, as one may expect to achieve
when modelling biological systems with no assumption on its time features.


  Searching new molecules in areas like drug discovery often starts from the
core structures of candidate molecules to optimize the properties of interest.
The way as such has called for a strategy of designing molecules retaining a
particular scaffold as a substructure. On this account, our present work
proposes a scaffold-based molecular generative model. The model generates
molecular graphs by extending the graph of a scaffold through sequential
additions of vertices and edges. In contrast to previous related models, our
model guarantees the generated molecules to retain the given scaffold with
certainty. Our evaluation of the model using unseen scaffolds showed the
validity, uniqueness, and novelty of generated molecules as high as the case
using seen scaffolds. This confirms that the model can generalize the learned
chemical rules of adding atoms and bonds rather than simply memorizing the
mapping from scaffolds to molecules during learning. Furthermore, despite the
restraint of fixing core structures, our model could simultaneously control
multiple molecular properties when generating new molecules.


  A generalisation of reaction diffusion systems and their travelling solutions
to cases when the productive part of the reaction happens only on a surface in
space or on a line on plane but the degradation and the diffusion happen in
bulk are important for modelling various biological processes. These include
problems of invasive species propagation along boundaries of ecozones, problems
of gene spread in such situations, morphogenesis in cavities, intracellular
reaction etc. Piecewise linear approximations of reaction terms in
reaction-diffusion systems often result in exact solutions of propagation front
problems. This article presents an exact travelling solution for a
reaction-diffusion system with a piecewise constant production restricted to a
codimension-1 subset. The solution is monotone, propagates with the unique
constant velocity, and connects the trivial solution to a nontrivial
nonhomogeneous stationary solution of the problem. The properties of the
solution closely parallel the properties of monotone travelling solutions in
classical bistable reaction-diffusion systems.


  While deep learning has achieved great success in many fields, one common
criticism about deep learning is its lack of interpretability. In most cases,
the hidden units in a deep neural network do not have a clear semantic meaning
or correspond to any physical entities. However, model interpretability and
explainability are crucial in many biomedical applications. To address this
challenge, we developed the Factor Graph Neural Network model that is
interpretable and predictable by combining probabilistic graphical models with
deep learning. We directly encode biological knowledge such as Gene Ontology as
a factor graph into the model architecture, making the model transparent and
interpretable. Furthermore, we devised an attention mechanism that can capture
multi-scale hierarchical interactions among biological entities such as genes
and Gene Ontology terms. With parameter sharing mechanism, the unrolled Factor
Graph Neural Network model can be trained with stochastic depth and generalize
well. We applied our model to two cancer genomic datasets to predict target
clinical variables and achieved better results than other traditional machine
learning and deep learning models. Our model can also be used for gene set
enrichment analysis and selecting Gene Ontology terms that are important to
target clinical variables.


  Disease psoriasis occurs as chronic inflammation of skin and appears as scaly
red lesions on skin surface. Advent of several immunosuppressive drugs
established that the disease stems from immuno-pathogenic disorder in human
blood. Cell biological as well as clinical research on the disease reveals that
the helper T-cells and other Leucocytes, responsible for human immunity, may
lead to psoriasis pathogenesis if produced in plenty at locations close to the
dermal region. Research findings also showed that a complex, self-sustaining
(cytokine and related) proteins network play important role in disease
maturation by actually leading to a huge proliferation of epidermal
keratinocytes. Disease pathogenesis is identified with such hyperproliferation
leading to flaking of skin surface (psoriatic plaques). An excessive generation
of nitric oxide by proliferated keratinocytes, through a complex chain of
bio-chemical events, is causal to the scaliness of psoriatic plaques.
Considering these immunopathogenic mechanisms, we propose and analyse a
mathematical time differential model for the disease psoriasis. Outcomes of
analysis are consistent with existing cell biological and clinical findings
with some new predictions which could be tested further.


  Artificial analogues of the natural nucleic acids have attracted recent
interest as a diverse class of information storage molecules capable of
self-replication. In the present study, we use the computational potential
energy landscape framework to investigate the structural and dynamical
properties of xylo- and deoxyxylo-nucleic acids (XyNA and dXyNA), which are
derived from their respective RNA and DNA analogues by an inversion of
configuration at a single chiral center in the sugar moiety of the nucleotide
unit. The free energy landscapes of an octameric XyNA sequence and its dXyNA
analogue demonstrate the existence of a facile conformational transition
between a left-handed helix that is the global free energy minimum, and a
closely competing ladder-type structure with approximately zero helicity. The
separation of the competing conformational ensembles is better-defined for the
dXyNA system, whereas the XyNA analogue is inherently more flexible. The former
therefore appear more suitable candidates for a molecular switch. The
landscapes differ qualitatively from those reported in previous studies for
evolved biomolecules: they are significantly more frustrated, so that XyNAs
provide an example of an unnatural system for which the conditions constituting
the principle of minimal frustration are, as may be expected, violated.


  Optimal sensor placement is an important yet unsolved problem in control
theory. In biological organisms, genetic activity is often highly nonlinear,
making it difficult to design libraries of promoters to act as reporters of the
cell state. We make use of the Koopman observability gramian to develop an
algorithm for optimal sensor (or reporter) placement for discrete time
nonlinear dynamical systems to ease the difficulty of design of the promoter
library. This ease is enabled due to the fact that the Koopman operator
represents the evolution of a nonlinear system linearly by lifting the states
to an infinite-dimensional space of observables. The Koopman framework ideally
demands high temporal resolution, but data in biology are often sampled
sparsely in time. Therefore we compute what we call the temporally fine-grained
Koopman operator from the temporally coarse-grained Koopman operator, the
latter of which is identified from the sparse data. The optimal placement of
sensors then corresponds to maximizing the observability of the fine-grained
system. We demonstrate the algorithm on a simulation example of a circadian
oscillator.


  The spatial structure of an evolving population affects which mutations
become fixed. Some structures amplify selection, increasing the likelihood that
beneficial mutations become fixed while deleterious mutations do not. Other
structures suppress selection, reducing the effect of fitness differences and
increasing the role of random chance. This phenomenon can be modeled by
representing spatial structure as a graph, with individuals occupying vertices.
Births and deaths occur stochastically, according to a specified update rule.
We study death-Birth updating: An individual is chosen to die and then its
neighbors compete to reproduce into the vacant spot. Previous numerical
experiments suggested that amplifiers of selection for this process are either
rare or nonexistent. We introduce a perturbative method for this problem for
weak selection regime, meaning that mutations have small fitness effects. We
show that fixation probability under weak selection can be calculated in terms
of the coalescence times of random walks. This result leads naturally to a new
definition of effective population size. Using this and other methods, we
uncover the first known examples of transient amplifiers of selection (graphs
that amplify selection for a particular range of fitness values) for the
death-Birth process. We also exhibit new families of "reducers of fixation",
which decrease the fixation probability of all mutations, whether beneficial or
deleterious.


  From fake news to innovative technologies, many contagions spread via a
process of social reinforcement, where multiple exposures are distinct from
prolonged exposure to a single source. Contrarily, biological agents such as
Ebola or measles are typically thought to spread as simple contagions. Here, we
demonstrate that interacting simple contagions are indistinguishable from
complex contagions. In the social context, our results highlight the challenge
of identifying and quantifying mechanisms, such as social reinforcement, in a
world where an innumerable amount of ideas, memes and behaviors interact. In
the biological context, this parallel allows the use of complex contagions to
effectively quantify the non-trivial interactions of infectious diseases.


  Most graph kernels are an instance of the class of $\mathcal{R}$-Convolution
kernels, which measure the similarity of objects by comparing their
substructures. Despite their empirical success, most graph kernels use a naive
aggregation of the final set of substructures, usually a sum or average,
thereby potentially discarding valuable information about the distribution of
individual components. Furthermore, only a limited instance of these approaches
can be extended to continuously attributed graphs. We propose a novel method
that relies on the Wasserstein distance between the node feature vector
distributions of two graphs, which allows to find subtler differences in data
sets by considering graphs as high-dimensional objects, rather than simple
means. We further propose a Weisfeiler-Lehman inspired embedding scheme for
graphs with continuous node attributes and weighted edges, enhance it with the
computed Wasserstein distance, and thus improve the state-of-the-art prediction
performance on several graph classification tasks.


  The processes of interplant competition within a field are still poorly
understood. However, they explain a large part of the heterogeneity in a field
and may have longer-term consequences, especially in mixed stands. Modeling can
help to better understand these phenomena but requires simulating the
interactions between different individuals. In the case of large populations,
assessing the parameters of a heterogeneous population model from experimental
data is intractable computationally. This paper investigates the mean-field
approximation of large dynamical systems with random initial conditions and
individual parameters, and with interaction being represented by pairwise
potentials between individuals. Under this approximation, each individual is in
interaction with an infinitely-crowded population, summarized by a probability
measure, the mean-field limit distribution, being itself the weak solution of a
non-linear hyperbolic partial differential equation. In particular, the
phenomenon of chaos propagation implies that the individuals are independent
asymptotically when the size of the population tends towards infinity. This
result provides perspectives for a possible simplification of the inference
problem. The simulation of the mean-field distribution, consisting in a
semi-Lagrangian scheme with an interpolation step using Gaussian process
regression, is illustrated for a heterogeneous population model representing
plants in competition for light.


  Spatial many-species predator-prey systems have been shown to yield very rich
space-time patterns. This observation begs the question whether there exist
universal mechanisms for generating this type of emerging complex patterns in
non-equilibrium systems. In this work we investigate the possibility of
dynamically generated hierarchies in predator-prey systems. We analyze a
nine-species model with competing interactions and show that the studied
situation results in the spontaneous formation of spirals within spirals. The
parameter dependence of these intriguing nested spirals is elucidated. This is
achieved through the numerical investigation of various quantities (correlation
lengths, densities of empty sites, Fourier analysis of species densities,
interface fluctuations) that allows to gain a rather complete understanding of
the spatial arrangements and the temporal evolution of the system. A possible
generalization of the interaction scheme yielding dynamically generated
hierarchies is discussed. As cyclic interactions occur spontaneously in systems
with competing strategies, the mechanism discussed in this work should
contribute to our understanding of various social and biological systems.


  Phylogenetics is the study of the evolutionary relationships between
organisms. One of the main challenges in the field is to take biological data
for a group of organisms and to infer an evolutionary tree, a graph that
represents these relationships. Developing practical and efficient methods for
inferring phylogenetic trees has lead to a number of interesting mathematical
questions across a variety of fields. However, due to hybridization and gene
flow, a phylogenetic network may be a better representation of the evolutionary
history of some groups of organisms. In this chapter, we introduce some of the
basic concepts in phylogenetics, and present related undergraduate research
projects on phylogenetic networks that touch on areas of graph theory and
abstract algebra. In the first section, we describe several open research
questions related to the combinatorics of phylogenetic networks. In the second,
we describe problems related to understanding phylogenetic statistical models
as algebraic varieties.


  DNA surface-hybridization biosensors utilize the selective hybridization of
target sequences in solution to surface-immobilized probes. In this process,
the target is usually assumed to be in excess, so that its concentration does
not significantly vary while hybridizing to the surface-bound probes. If the
target is initially at low concentrations and/or if the number of probes is
very large and have high affinity for the target, the DNA in solution may get
depleted. In this paper we analyze the equilibrium and kinetics of
hybridization of DNA biosensors in the case of strong target depletion, by
extending the Langmuir adsorption model. We focus, in particular, on the
detection of a small amount of a single-nucleotide "mutant" sequence
(concentration $c_2$) in a solution, which differs by one or more nucleotides
from an abundant "wild-type" sequence (concentration $c_1 \gg c_2$). We show
that depletion can give rise to a strongly-enhanced sensitivity of the
biosensors. Using representative values of rate constants and hybridization
free energies, we find that in the depletion regime one could detect relative
concentrations $c_2/c_1$ that are up to three orders of magnitude smaller than
in the conventional approach. The kinetics is surprisingly rich, and exhibits a
non-monotonic adsorption with no counterpart in the no-depletion case. Finally,
we show that, alongside enhanced detection sensitivity, this approach offers
the possibility of sample enrichment, by substantially increasing the relative
amount of the mutant over the wild-type sequence.


  Amphiphilic molecules such as surfactants, lipids and block copolymers can be
assembled into bilayers and form vesicles. Fluorescent membrane labeling
methods require the use of dye molecules that can be inserted in the bilayers
at different stages of the synthesis. To our knowledge there is no generalized
method for labeling preformed vesicles. Herein we develop a versatile protocol
that is suitable to both surfactant and lipid preformed vesicles and requires
no separation or purification steps. Based on the lipophilic carbocyanine green
dye PKH67, the methodology is assessed on zwitterionic phosphatidylcholine
vesicles. To demonstrate its versatility, it is applied to dispersions of
anionic or cationic vesicles, such as a drug administrated to premature infants
with respiratory distress syndrome, or a vesicle formulation used as a fabric
softener for homecare applications. By means of fluorescence microscopy we then
visualize the interaction mechanisms of nanoparticles crossing live cell
membranes and of surfactant adsorbing on cotton fabric. These results highlight
the advantages of a membrane labeling technique that is simple and applicable
to a large number of soft matter systems.


  Microbial electrolysis cells (MECs) employ electroactive bacteria to perform
extracellular electron transfer, enabling hydrogen generation from
biodegradable substrates. In previous work, we developed and analyzed a
differential-algebraic equation (DAE) model for MECs. The model resembles a
chemostat with ordinary differential equations (ODEs) for concentrations of
substrate, microorganisms, and an extracellular mediator involved in electron
transfer. There is also an algebraic constraint for electric current and
hydrogen production. Our goal is to determine the outcome of competition
between methanogenic archaea and electroactive bacteria, because only the
latter contribute to electric current and resulting hydrogen production. We
investigate asymptotic stability in two industrially relevant versions of the
model. An important aspect of chemostats models is the principle of competitive
exclusion -- only microbes which grow at the lowest substrate concentration
will survive as $t\to\infty$. We show that if methanogens grow at the lowest
substrate concentration, then the equilibrium corresponding to competitive
exclusion by methanogens is globally asymptotically stable. The analogous
result for electroactive bacteria is not necessarily true. We show that local
asymptotic stability of exclusion by electroactive bacteria is not guaranteed,
even in a simplified version of the model. In this case, even if electroactive
bacteria can grow at the lowest substrate concentration, a few additional
conditions are required to guarantee local asymptotic stability. We also
provide numerical simulations supporting these arguments. Our results suggest
operating conditions that are most conducive to success of electroactive
bacteria and the resulting current and hydrogen production in MECs. This will
help identify when methane production or electricity and hydrogen production
are favored.


  Historically, the majority of statistical association methods have been
designed assuming availability of SNP-level information. However, modern
genetic and sequencing data present new challenges to access and sharing of
genotype-phenotype datasets, including cost management, difficulties in
consolidation of records across research groups, etc. These issues make methods
based on SNP-level summary statistics for a joint analysis of variants in a
group particularly appealing. The most common form of combining statistics is a
sum of SNP-level squared scores, possibly weighted, as in burden tests for rare
variants. The overall significance of the resulting statistic is evaluated
using its distribution under the null hypothesis. Here, we demonstrate that
this basic approach can be substantially improved by decorrelating scores prior
to their addition, resulting in remarkable power gains in situations that are
most commonly encountered in practice; namely, under heterogeneity of effect
sizes and diversity between pairwise LD. In these situations, the power of the
traditional test, based on the added squared scores, quickly reaches a ceiling,
as the number of variants increases. Thus, the traditional approach does not
benefit from information potentially contained in any additional SNPs, while
our decorrelation by orthogonal transformation (DOT) method yields steady gain
in power. We present theoretical and computational analyses of both approaches,
and reveal causes behind sometimes dramatic difference in their respective
powers. We showcase DOT by analyzing breast cancer data, in which our method
strengthened levels of previously reported associations and implied the
possibility of multiple new alleles that jointly confer breast cancer risk.


  Metrics on rooted phylogenetic trees are integral to a number of areas of
phylogenetic analysis. Cluster-similarity metrics have recently been introduced
in order to limit skew in the distribution of distances, and to ensure that
trees in the neighbourhood of each other have similar hierarchies. In the
present paper we introduce a new cluster-similarity metric on rooted
phylogenetic tree space that has an associated local operation, allowing for
easy calculation of neighbourhoods, a trait that is desirable for MCMC
calculations. The metric is defined by the distance on the Hasse diagram
induced by a partial order on the set of rooted phylogenetic trees, itself
based on the notion of a hierarchy-preserving map between trees. The partial
order we introduce is a refinement of the well-known refinement order on
hierarchies. Both the partial order and the hierarchy-preserving maps may also
be of independent interest.


  Computational drug discovery provides an efficient tool helping large scale
lead molecules screening. One of the major tasks of lead discovery is
identifying molecules with promising binding affinities towards a target, a
protein in general. The accuracies of current scoring functions which are used
to predict the binding affinity are not satisfactory enough. Thus, machine
learning (ML) or deep learning (DL) based methods have been developed recently
to improve the scoring functions. In this study, a deep convolutional neural
network (CNN) model (called OnionNet) is introduced and the features are based
on rotation-free element-pair specific contacts between ligands and protein
atoms, and the contacts were further grouped in different distance ranges to
cover both the local and non-local interaction information between the ligand
and the protein. The prediction power of the model is evaluated and compared
with other scoring functions using the comparative assessment of scoring
functions (CASF-2013) benchmark and the v2016 core set of PDBbind database.
When compared to a previous CNN-based scoring function, our model shows
improvements of 0.08 and 0.16 in the correlations (R) and standard deviations
(SD) of regression, respectively, between the predicted binding affinities and
the experimental measured binding affinities. The robustness of the model is
further explored by predicting the binding affinities of the complexes
generated from docking simulations instead of experimentally determined PDB
structures.


  The mitotic spindle lies at the heart of the spatio-temporal control over
cellular components during cell division. The spindle consists of microtubules,
which are not only crosslinked by motor proteins but also by passive binding
proteins. These passive crosslinkers stabilize the highly dynamic mitotic
spindle by generating friction forces between sliding filaments. However, it
remains unclear how the friction coefficient depends on the number of
crosslinkers and the size of the overlap between the microtubules. Here, we use
theory and computer simulations to study the friction between two filaments
that are crosslinked by passive proteins, which can hop between neighboring
binding sites while physically excluding each other. The simulations reveal
that the movement of one microtubule relative to the other is limited by
free-energy barrier crossings, causing rare and discrete jumps of the
microtubule that span the distance between adjacent crosslinker binding sites.
We derive an exact analytical expression for the free-energy landscape and
identify the reaction coordinate that governs the relative movement, which
allows us to determine the effective barrier height for the microtubule jumps.
Both through simulations and reaction rate theory, we make the experimentally
testable prediction that the friction between the microtubules increases
superexponentially with the density of crosslinkers.


  We once proposed that cell-type-associated chromatin configurations determine
cell types and that cancer cell type is determined by cancer-associated
chromatin configuration (CACC). In this paper, we hypothesize that flexible
cell-type-associated chromatin configuration is associated with cell potency
and has an advantage over inflexible one in regulating genome related
activities, such as DNA replication, DNA transcription, DNA repair, and DNA
mutagenesis. The reason why cancer is so difficult to treat is because CACC is
flexible, which enables cancer cells not only to produce heterogeneous
subclones through limited cell differentiation, but also to maximally and
efficiently use genome related resources to survive environmental changes.
Therefore, to beat cancer, more efforts should be made to restrict the
flexibility of CACC or to change CACC so that cancer cells can be turned back
to normal or become less malignant.


  The fixation probability of a single mutant invading a population of
residents is among the most widely-studied quantities in evolutionary dynamics.
Amplifiers of natural selection are population structures that increase the
fixation probability of advantageous mutants, compared to well-mixed
populations. Extensive studies have shown that many amplifiers exist for the
Birth-death Moran process, some of them substantially increasing the fixation
probability or even guaranteeing fixation in the limit of large population
size. On the other hand, no amplifiers are known for the death-Birth Moran
process, and computer-assisted exhaustive searches have failed to discover
amplification. In this work we resolve this disparity, by showing that any
amplification under death-Birth updating is necessarily \emph{bounded} and
\emph{transient}. Our boundedness result states that even if a population
structure does amplify selection, the resulting fixation probability is close
to that of the well-mixed population. Our transience result states that for any
population structure there exists a threshold $r^*$ such that the population
structure ceases to amplify selection if the mutant fitness advantage $r$ is
larger than $r^\star$. Finally, we also extend the above results to
$\delta$-death-Birth updating, which is a combination of Birth-death and
death-Birth updating. On the positive side, we identify population structures
that maintain amplification for a wide range of values $r$ and $\delta$. These
results demonstrate that amplification of natural selection depends on the
specific mechanisms of the evolutionary process.


  In the gonad of C. elegans worms, germ cells form the outer lining of the
central tube, the rachis, which carries a cytoplasmic flow.


  Recently several deep learning models have been used for DNA sequence based
classification tasks. Often such tasks require long and variable length DNA
sequences in the input. In this work, we use a sequence-to-sequence autoencoder
model to learn a latent representation of a fixed dimension for long and
variable length DNA sequences in an unsupervised manner. We evaluate both
quantitatively and qualitatively the learned latent representation for a
supervised task of splice site classification. The quantitative evaluation is
done under two different settings. Our experiments show that these
representations can be used as features or priors in closely related tasks such
as splice site classification. Further, in our qualitative analysis, we use a
model attribution technique Integrated Gradients to infer significant sequence
signatures influencing the classification accuracy. We show the identified
splice signatures resemble well with the existing knowledge.


  In the cell, proteins fold and perform complex functions through global
structural rearrangements. Function requires a protein to be at the brink of
stability to be susceptible to small environmental fluctuations, yet stable
enough to maintain structural integrity. These apparently conflicting behaviors
are exhibited by systems near a critical point, where distinct phases merge $-$
a concept beyond previous studies indicating proteins have a well-defined
folded/unfolded phase boundary in the pressure-temperature plane. Here, by
modeling the protein phosphoglycerate kinase (PGK) on the temperature (T),
pressure (P), and crowding volume-fraction ($\phi$) phase diagram, we
demonstrate a critical transition where phases merge, and PGK exhibits large
structural fluctuations. Above the critical temperature (Tc), the difference
between the intermediate and unfolded phases disappears. When $\phi$ increases,
the Tc moves to a lower T. We verify the calculations with experiments mapping
the T-P-$\phi$ space, which likewise reveal a critical point at 305 K and 170
MPa that moves to a lower T as $\phi$ increases. Crowding places PGK near a
critical line in its natural parameter space, where large conformational
changes can occur without costly free energy barriers. Specific structures are
proposed for each phase based on simulation.


  Evolution is the principal theoretical framework in Pacific Island
archaeology. Progressive cultural evolutionary stages, particularly the
chiefdom, provided the earliest framework. These ideas were increasingly linked
to islands as relatively isolated laboratories of cultural evolution and
diversification, albeit their populations phylogenetically linked in a series
of ancestral and descendant societies. With more sophisticated environmental
data, much research has investigated cultural adaptations and their homologous
or analogous origins. Most recently, Darwinian mechanisms of transmission,
selection, and drift comprise an important approach. The history of evolution
in Pacific archaeology suggests it will continue as a useful theoretical
framework, perhaps more so than in other world regions.


  Organisms have evolved a variety of mechanisms to cope with the
unpredictability of environmental conditions, and yet mainstream models of
metabolic regulation are typically based on strict optimality principles that
do not account for uncertainty. This paper introduces a dynamic metabolic
modelling framework that is a synthesis of recent ideas on resource allocation
and the powerful optimal control formulation of Ramkrishna and colleagues. In
particular, their work is extended based on the hypothesis that cellular
resources are allocated among elementary flux modes according to the principle
of maximum entropy. These concepts both generalise and unify prior approaches
to dynamic metabolic modelling by establishing a smooth interpolation between
dynamic flux balance analysis and dynamic metabolic models without regulation.
The resulting theory is successful in describing `bet-hedging' strategies
employed by cell populations dealing with uncertainty in a fluctuating
environment, including heterogenous resource investment, accumulation of
reserves in growth-limiting conditions, and the observed behaviour of yeast
growing in batch and continuous cultures. The maximum entropy principle is also
shown to yield an optimal control law consistent with partitioning resources
between elementary flux mode families, which has important practical
implications for model reduction, selection, and simulation.


  We study families of chemical reaction networks whose positive steady states
are toric, and therefore can be parameterized by monomials. Families are
constructed algorithmically from a core network; we show that if a family
member is multistationary, then so are all subsequent networks in the family.
Further, we address the questions of model selection and experimental design
for families by investigating the algebraic dependencies of the chemical
concentrations using matroids. Given a family with toric steady states and a
constant number of conservation relations, we construct a matroid that encodes
important information regarding the steady state behaviour of the entire
family. Among other things, this gives necessary conditions for the
distinguishability of families of reaction networks with respect to a data set
of measured chemical concentrations. We illustrate our results using multi-site
phosphorylation networks.


  We study the evolution of cooperation as a birth-death process in spatially
extended populations. The benefit from the altruistic behavior of a cooperator
is implemented by decreasing the death rate of its direct neighbors. The cost
of cooperation is the increase of a cooperator's death rate proportional to the
number of its neighbors. When cooperation has higher cost than benefit,
cooperators disappear. Then the dynamics reduces to the contact process with
defectors as the single particle type. Increasing the benefit-cost ratio above
1, the extinction transition of the contact process splits into a set of
nonequilibrium transitions between four regimes when increasing the baseline
death rate $p$ as a control parameter: (i) defection only, (ii) coexistence,
(iii) cooperation only, (iv) extinction. We investigate the transitions between
these regimes. As the main result, we find that full cooperation is established
at the extinction transition as long as benefit is strictly larger than cost.
Qualitatively identical phase diagrams are obtained for populations on square
lattices and in pair approximation. Spatial correlations with nearest neighbors
only are thus sufficient for sustained cooperation.


  Determining which proteins interact together is crucial to a systems-level
understanding of the cell. Recently, algorithms based on Direct Coupling
Analysis (DCA) pairwise maximum-entropy models have allowed to identify
interaction partners among paralogous proteins from sequence data. This success
of DCA at predicting protein-protein interactions could be mainly based on its
known ability to identify pairs of residues that are in contact in the
three-dimensional structure of protein complexes and that coevolve to remain
physicochemically complementary. However, interacting proteins possess similar
evolutionary histories. What is the role of purely phylogenetic correlations in
the performance of DCA-based methods to infer interaction partners? To address
this question, we employ controlled synthetic data that only involve phylogeny
and no interactions or contacts. We find that DCA accurately identifies the
pairs of synthetic sequences that share evolutionary history. While
phylogenetic correlations confound the identification of contacting residues by
DCA, they are thus useful to predict interacting partners among paralogs. We
find that DCA performs as well as phylogenetic methods to this end, and
slightly better than them with large and accurate training sets. Employing DCA
or phylogenetic methods within an Iterative Pairing Algorithm (IPA) allows to
predict pairs of evolutionary partners without a training set. We demonstrate
the ability of these various methods to correctly predict pairings among real
paralogous proteins with genome proximity but no known physical interaction,
illustrating the importance of phylogenetic correlations in natural data.
However, for physically interacting and strongly coevolving proteins, DCA and
mutual information outperform phylogenetic methods. We discuss how to
distinguish physically interacting proteins from those only sharing
evolutionary history.


  Rearrangement operations transform a phylogenetic tree into another one and
hence induce a metric on the space of phylogenetic trees. Popular operations
for unrooted phylogenetic trees are NNI (nearest neighbour interchange), SPR
(subtree prune and regraft), and TBR (tree bisection and reconnection).
Recently, these operations have been extended to unrooted phylogenetic
networks, which are generalisations of phylogenetic trees that can model
reticulated evolutionary relationships. Here, we study global and local
properties of spaces of phylogenetic networks under these three operations. In
particular, we prove connectedness and asymptotic bounds on the diameters of
spaces of different classes of phylogenetic networks, including tree-based and
level-k networks. We also examine the behaviour of shortest TBR-sequence
between two phylogenetic networks in a class, and whether the TBR-distance
changes if intermediate networks from other classes are allowed: for example,
the space of phylogenetic trees is an isometric subgraph of the space of
phylogenetic networks under TBR. Lastly, we show that computing the
TBR-distance and the PR-distance of two phylogenetic networks is NP-hard.


  Given native 2D contact map, protein 3D structure could be reconstructed with
accuracy of 2A or better, and such reconstruction is a feasible computational
approach for protein folding problem. The prediction accuracy from traditional
methods is generally too poor to useful, but the recent deep learning model has
significantly improved the accuracy. In this study, we proposed a neural
network model comprising a bi-directional recurrent neural network and
artificial neural network. Over the non-redundant database of all available
protein 3D structures in Protein Data Bank, this deep learning model achieved
an accuracy of 0.80, much higher than those of previous models. This study
represents a major breakthrough in protein 2D contact map prediction and likely
a major step forward for the protein folding problem.


  In the context of the ongoing Ebola epidemic in DRC, active conflict and
community distrust are undermining control efforts, including vaccination
strategies. In this paper, we employed an individual-level stochastic
structured transmission model to assess the impact of vaccination strategies on
epidemic control in the context of variable levels of household
inaccessibility. We found that a ring vaccination strategy of close contacts
would not be effective for containing the epidemic in the context of
significant delays to vaccinating contacts even for low levels of household
inaccessibility and evaluate the impact of a supplemental community vaccination
strategy. For lower levels of inaccessibility, the probability of epidemic
containment increases over time. For higher levels of inaccessibility, even the
combined ring and community vaccination strategies are not expected to contain
the epidemic even though they help lower incidence levels, which saves lives,
makes the epidemic easier to contain and reduces spread to other communities.
We found that ring vaccination is effective for containing an outbreak until
the levels of inaccessibility exceeds approximately 10%, a combined ring and
community vaccination strategy is effective until the levels of inaccessibility
exceeds approximately 50%. Our findings underscore the need to enhance
community engagement to public health interventions.


  Relaxed random walk (RRW) models of trait evolution introduce branch-specific
rate multipliers to modulate the variance of a standard Brownian diffusion
process along a phylogeny and more accurately model overdispersed biological
data. Increased taxonomic sampling challenges inference under RRWs as the
number of unknown parameters grows with the number of taxa. To solve this
problem, we present a scalable method to efficiently fit RRWs and infer this
branch-specific variation in a Bayesian framework. We develop a Hamiltonian
Monte Carlo (HMC) sampler to approximate the high-dimensional, correlated
posterior that exploits a closed-form evaluation of the gradient of the trait
data log-likelihood with respect to all branch-rate multipliers simultaneously.
Our gradient calculation achieves computational complexity that scales only
linearly with the number of taxa under study. We compare the efficiency of our
HMC sampler to the previously standard univariable Metropolis-Hastings approach
while studying the spatial emergence of the West Nile virus in North America in
the early 2000s. Our method achieves an over 300-fold speed-increase over the
univariable approach. Additionally, we demonstrate the scalability of our
method by applying the RRW to study the correlation between five mammalian life
history traits in a phylogenetic tree with 3650 tips.


  State-free reversible VAMPnets (SRVs) are a neural network-based framework
capable of learning the leading eigenfunctions of the transfer operator of a
dynamical system from trajectory data. In molecular dynamics simulations, these
data-driven collective variables (CVs) capture the slowest modes of the
dynamics and are useful for enhanced sampling and free energy estimation. In
this work, we employ SRV coordinates as a feature set for Markov state model
(MSM) construction. Compared to the current state of the art, MSMs constructed
from SRV coordinates are more robust to the choice of input features, exhibit
faster implied timescale convergence, and permit the use of shorter lagtimes to
construct higher kinetic resolution models. We apply this methodology to study
the folding kinetics and conformational landscape of the Trp-cage miniprotein.
Folding and unfolding mean first passage times are in good agreement with prior
literature, and a nine macrostate model is presented. The unfolded ensemble
comprises a central kinetic hub with interconversions to several metastable
unfolded conformations and which serves as the gateway to the folded ensemble.
The folded ensemble comprises the native state, a partially unfolded
intermediate "loop" state, and a previously unreported short-lived intermediate
that we were able to resolve due to the high time-resolution of the SRV-MSM. We
propose SRVs as an excellent candidate for integration into modern MSM
construction pipelines.


  We present a computational study on the folding and aggregation of proteins
in aqueous environment, as function of its concentration. We show how the
increase of the concentration of individual protein species can induce a
partial unfolding of the native conformation without the occurrence of
aggregates. A further increment of the protein concentration results in the
complete loss of the folded structures and induces the formation of protein
aggregates. We discuss the effect of the protein interface on the water
fluctuations in the protein hydration shell and their relevance in the
protein-protein interaction.


  The validity of a complex reaction pathway proposed to treat Inflammatory
Bowel Disease (IBD) was verified by a comprehensive time and frequency domain
analysis. The model was taken to the frequency domain to study the effect and
the significance of the negative feedback loop introduced by the reaction
pathways. It could be shown that such proposed probiotics have very interesting
potentials that could be used extensively in the near future.


  Markov models of character substitution on phylogenies form the foundation of
phylogenetic inference frameworks. Early models made the simplifying assumption
that the substitution process is homogeneous over time and across sites in the
molecular sequence alignment. While standard practice adopts extensions that
accommodate heterogeneity of substitution rates across sites, heterogeneity in
the process over time in a site-specific manner remains frequently overlooked.
This is problematic, as evolutionary processes that act at the molecular level
are highly variable, subjecting different sites to different selective
constraints over time, impacting their substitution behaviour. We propose
incorporating time variability through Markov-modulated models (MMMs) that
allow the substitution process (including relative character exchange rates as
well as the overall substitution rate) that models the evolution at an
individual site to vary across lineages. We implement a general MMM framework
in BEAST, a popular Bayesian phylogenetic inference software package, allowing
researchers to compose a wide range of MMMs through flexible XML specification.
Using examples from bacterial, viral and plastid genome evolution, we show that
MMMs impact phylogenetic tree estimation and can substantially improve model
fit compared to standard substitution models. Through simulations, we show that
marginal likelihood estimation accurately identifies the generative model and
does not systematically prefer the more parameter-rich MMMs. In order to
mitigate the increased computational demands associated with MMMs, our
implementation exploits recently developed updates to BEAGLE, a
high-performance computational library for phylogenetic inference.


  The discovery of genomic polymorphisms influencing gene expression (also
known as expression quantitative trait loci or eQTLs) can be formulated as a
sparse Bayesian multivariate/multiple regression problem. An important aspect
in the development of such models is the implementation of bespoke inference
methodologies, a process which can become quite laborious, when multiple
candidate models are being considered. We describe automatic, black-box
inference in such models using Stan, a popular probabilistic programming
language. The utilisation of systems like Stan can facilitate model prototyping
and testing, thus accelerating the data modelling process. The code described
in this chapter can be found at https://github.com/dvav/eQTLBookChapter.


  Reaction networks are often used to model interacting species in fields such
as biochemistry and ecology. When the counts of the species are sufficiently
large, the dynamics of their concentrations are typically modeled via a system
of differential equations. However, when the counts of some species are small,
the dynamics of the counts are typically modeled stochastically via a discrete
state, continuous time Markov chain.
  A key quantity of interest for such models is the probability mass function
of the process at some fixed time. Since paths of such models are relatively
straightforward to simulate, we can estimate the probabilities by constructing
an empirical distribution. However, the support of the distribution is often
diffuse across a high-dimensional state space, where the dimension is equal to
the number of species. Therefore generating an accurate empirical distribution
can come with a large computational cost.
  We present a new Monte Carlo estimator that fundamentally improves on the
"classical" Monte Carlo estimator described above. It also preserves much of
classical Monte Carlo's simplicity. The idea is basically one of conditional
Monte Carlo. Our conditional Monte Carlo estimator has two parameters, and
their choice critically affects the performance of the algorithm. Hence, a key
contribution of the present work is that we demonstrate how to approximate
optimal values for these parameters in an efficient manner. Moreover, we
provide a central limit theorem for our estimator, which leads to approximate
confidence intervals for its error.


  The outcome of an epidemic is closely related to the network of interactions
between the individuals. Likewise, protein functions depend on the 3D
arrangement of their residues and on the underlying energetic interaction
network. Borrowing ideas from the theoretical framework that has been developed
to address the spreading of real diseases, we study the diffusion of a
fictitious epidemic inside the protein non-bonded interaction network. Our
approach allowed to probe the overall stability and the capability to propagate
information in the complex 3D-structures and proved to be very efficient in
addressing different problems, from the assessment of thermal stability to the
identification of allosteric sites.


  We investigate a population-genetic model with a temporally-fluctuating
sawtooth fitness landscape. We numerically show that a counter-intuitive
behavior occurs where the rate of evolution of the system decreases as
selection pressure increases from zero. This phenomenon is understood by
analogy with absolute negative mobility in particle flow. A phenomenological
explanation about the direction of evolution is also provided.


  In this work, we have developed a Coxian distributed SEIR model in
incorporating an empirical incubation period. We show that the global dynamics
are completely determined by a basic reproduction number. An application of the
Coxian distributed SEIR model using data of an empirical incubation period is
explored. The model may be useful for resolving causing the realistic intrinsic
parts in classical epidemic models since Coxian distribution approximately
converges to any distribution.


  Many proteins have the potential to aggregate into amyloid fibrils, which are
associated with a wide range of human disorders including Alzheimer's and
Parkinson's disease. In contrast to that of folded proteins, the thermodynamic
stability of amyloid fibrils is not well understood: specifically the balance
between entropic and enthalpic terms, including the chain entropy and the
hydrophobic effect, are poorly characterised. Using simulations of a
coarse-grained protein model we delineate the enthalpic and entropic
contributions dominating amyloid fibril elongation, predicting a characteristic
temperature-dependent enthalpic signature. We confirm this thermodynamic
signature by performing calorimetric experiments and a meta-analysis over
published data. From these results, we can also elucidate the necessary
conditions to observe cold denaturation of amyloid fibrils. Overall, we show
that amyloid fibril elongation is associated with a negative heat capacity, the
magnitude of which correlates closely with the hydrophobic surface area that is
buried upon fibril formation, highlighting the importance of hydrophobicity for
fibril stability.


  We study a class of multi-species birth-and-death processes going almost
surely to extinction and admitting a unique quasi-stationary distribution (qsd
for short). When rescaled by $K$ and in the limit $K\to+\infty$, the
realizations of such processes get close, in any fixed finite-time window, to
the trajectories of a dynamical system whose vector field is defined by the
birth and death rates. Assuming that this dynamical has a unique attracting
fixed point, we analyzed in a previous work what happens for large but finite
$K$, especially the different time scales showing up. In the present work, we
are mainly interested in the following question: Observing a realization of the
process, can we determine the so-called engineering resilience? To answer this
question, we establish two relations which intermingle the resilience, which is
a macroscopic quantity defined for the dynamical system, and the fluctuations
of the process, which are microscopic quantities. Analogous relations are well
known in nonequilibrium statistical mechanics. To exploit these relations, we
need to introduce several estimators which we control for times between $\log
K$ (time scale to converge to the qsd) and $\exp(K)$ (time scale of mean time
to extinction).


  I survey some recent approaches to studying change in the lexicon,
particularly change in meaning across phylogenies. I briefly sketch an
evolutionary approach to language change and point out some issues in recent
approaches to studying semantic change that rely on temporally stratified word
embeddings. I draw illustrations from lexical cognate models in Pama-Nyungan to
identify meaning classes most appropriate for lexical phylogenetic inference,
particularly highlighting the importance of variation in studying change over
time.


  We consider a modified Rosenzweig-MacArthur predator-prey model, based on the
premise that the search rate of predators is dependent on the prey density,
rather than constant. A complete analysis of the global behavior of the model
is presented, and shows that the model exhibits a dichotomy similar to the
classical Rosenzweig-MacArthur model: either the coexistence steady state is
globally stable; or it is unstable, and then a unique, globally stable limit
cycle exists. We discuss the similarities, but also important differences
between our model and the Rosenweig-MacArthur model. The main differences are
that: 1. The paradox of enrichment which always occurs in the
Rosenzweig-MacArthur model, does not always occur here, and 2. Even when the
paradox of enrichment occurs, predators can adapt by lowering their search
rate, and effectively stabilize the system.


  Tree-based networks are a class of phylogenetic networks that attempt to
formally capture what is meant by "tree-like" evolution. A given non-tree-based
phylogenetic network, however, might appear to be very close to being
tree-based, or very far. In this paper, we formalise the notion of proximity to
tree-based for unrooted phylogenetic networks, with a range of proximity
measures. These measures also provide characterisations of tree-based networks.
One measure in particular, related to the nearest neighbour interchange
operation, allows us to define the notion of "tree-based rank". This provides a
subclassification within the tree-based networks themselves, identifying those
networks that are "very" tree-based. Finally, we prove results relating
tree-based networks in the settings of rooted and unrooted phylogenetic
networks, showing effectively that an unrooted network is tree-based if and
only if it can be made a rooted tree-based network by rooting it and orienting
the edges appropriately. This leads to a clarification of the contrasting
decision problems for tree-based networks, which are polynomial in the rooted
case but NP complete in the unrooted.


  Membrane neck formation is essential for scission, which, as recent
experiments on tubules have demonstrated, can be location dependent. The
diversity of biological machinery that can constrict a neck such as dynamin,
actin, ESCRTs and BAR proteins, and the range of forces and deflection over
which they operate, suggest that the constriction process is functionally
mechanical and robust to changes in biological environment. In this study, we
used a mechanical model of the lipid bilayer to systematically investigate the
influence of location, symmetry constraints, and helical forces on membrane
neck constriction. Simulations from our model demonstrated that the energy
barriers associated with constriction of a membrane neck are
location-dependent. Importantly, if symmetry restrictions are relaxed, then the
energy barrier for constriction is dramatically lowered and the membrane
buckles at lower values of forcing parameters. Our simulations also show that
constriction due to helical proteins further reduces the energy barrier for
neck formation compared to cylindrical proteins. These studies establish that
despite different molecular mechanisms of neck formation in cells, the
mechanics of constriction naturally leads to a loss of symmetry that can lower
the energy barrier to constriction.


  We develop a greedy algorithm that is fast and scalable in the detection of a
nested partition extracted from a dendrogram obtained from hierarchical
clustering of a multivariate series. Our algorithm provides a $p$-value for
each clade observed in the hierarchical tree. The $p$-value is obtained by
computing a number of bootstrap replicas of the dissimilarity matrix and by
performing a statistical test on each difference between the dissimilarity
associated with a given clade and the dissimilarity of the clade of its parent
node. We prove the efficacy of our algorithm with a set of benchmarks generated
by using a hierarchical factor model. We compare the results obtained by our
algorithm with those of Pvclust. Pvclust is a widely used algorithm developed
with a global approach originally motivated by phylogenetic studies. In our
numerical experiments we focus on the role of multiple hypothesis test
correction and on the robustness of the algorithms to inaccuracy and errors of
datasets. We also apply our algorithm to a reference empirical dataset. We
verify that our algorithm is much faster than Pvclust algorithm and has a
better scalability both in the number of elements and in the number of records
of the investigated multivariate set. Our algorithm provides a hierarchically
nested partition in much shorter time than currently widely used algorithms
allowing to perform a statistically validated cluster analysis detection in
very large systems.


  The abundance of different species in a community often follows the log
series distribution. Other ecological patterns also have simple forms. Why does
the complexity and variability of ecological systems reduce to such simplicity?
Common answers include maximum entropy, neutrality, and convergent outcome from
different underlying biological processes. This article proposes a more general
answer based on the concept of invariance, the property by which a pattern
remains the same after transformation. Invariance has a long tradition in
physics. For example, general relativity emphasizes the need for the equations
describing the laws of physics to have the same form in all frames of
reference. By bringing this unifying invariance approach into ecology, we show
that the log series pattern dominates when the consequences of processes acting
on abundance are invariant to the addition or multiplication of abundance by a
constant. The lognormal pattern dominates when the processes acting on net
species growth rate obey rotational invariance (symmetry) with respect to the
summing up of the individual component processes. Recognizing how these
invariances connect pattern to process leads to a synthesis of previous
approaches. First, invariance provides a simpler and more fundamental maximum
entropy derivation of the log series distribution. Second, invariance provides
a simple derivation of the key result from neutral theory: the log series at
the metacommunity scale and a clearer form of the skewed lognormal at the local
community scale. The invariance expressions are easy to understand because they
uniquely describe the basic underlying components that shape pattern.


  Sexual reproduction is not always synonymous with the existence of two
morphologically different sexes; isogamous species produce sex cells of equal
size, typically falling into multiple distinct self-incompatible classes,
termed mating types. A long-standing open question in evolutionary biology is:
what governs the number of these mating types across species? Simple
theoretical arguments imply an advantage to rare types, suggesting the number
of types should grow consistently, however, empirical observations are very
different. While some isogamous species exhibit thousands of mating types, such
species are exceedingly rare, and most have fewer than ten. In this paper, we
present a mathematical analysis to quantify the role of fitness variation -
characterised by different mortality rates - in determining the number mating
types emerging in simple evolutionary models. We predict that the number of
mating types decreases as the variance of mortality increases.


  Long non-coding RNA, microRNA, and messenger RNA enable key regulations of
various biological processes through a variety of diverse interaction
mechanisms. Identifying the interactions and cross-talk between these
heterogeneous RNA classes is essential in order to uncover the functional role
of individual RNA transcripts, especially for unannotated and newly-discovered
RNA sequences with no known interactions. Recently, sequence-based deep
learning and network embedding methods are becoming promising approaches that
can either predict RNA-RNA interactions from a sequence or infer missing
interactions from patterns that may exist in the network topology. However, the
majority of these methods have several limitations, eg, the inability to
perform inductive predictions, to distinguish the directionality of
interactions, or to integrate various sequence, interaction, and annotation
biological datasets. We proposed a novel deep learning-based framework,
rna2rna, which learns from RNA sequences to produce a low-dimensional embedding
that preserves the proximities in both the interactions topology and the
functional affinity topology. In this proposed embedding space, we have
designated a two-part" source and target contexts" to capture the targeting and
receptive fields of each RNA transcript, while encapsulating the heterogenous
cross-talk interactions between lncRNAs and miRNAs. From experimental results,
our method exhibits superior performance in AUPR rates compared to state-of-art
approaches at predicting missing interactions in different RNA-RNA interaction
databases and was shown to accurately perform link predictions to novel RNA
sequences not seen at training time, even without any prior information.
Additional results suggest that our proposed framework can capture a manifold
for heterogeneous RNA sequences to discover novel functional annotations.


  This paper studies the relationship between undirected (unrooted) and
directed (rooted) phylogenetic networks. We describe a polynomial-time
algorithm for deciding whether an undirected nonbinary phylogenetic network,
given the locations of the root and reticulation vertices, can be oriented as a
directed nonbinary phylogenetic network. Moreover, we characterize when this is
possible and show that, in such instances, the resulting directed nonbinary
phylogenetic network is unique. In addition, without being given the location
of the root and the reticulation vertices, we describe an algorithm for
deciding whether an undirected binary phylogenetic network $N$ can be oriented
as a directed binary phylogenetic network of a certain class. The algorithm is
fixed-parameter tractable (FPT) when the parameter is the level of $N$ and is
applicable to classes of directed phylogenetic networks that satisfy certain
conditions. As an example, we show that the well-studied class of binary
tree-child networks satisfies these conditions.


  Viruses evolve in the background of host immune systems that exert selective
pressure and drive viral evolutionary trajectories. This interaction leads to
different evolutionary patterns in antigenic space. Examples observed in nature
include the effectively one-dimensional escape characteristic of influenza A
and the prolonged coexistence of lineages in influenza B. Here we use an
evolutionary model for viruses in the presence of immune host systems with
finite memory to delineate parameter regimes of these patterns in a in
two-dimensional antigenic space. We find that for small effective mutation
rates and mutation jump ranges, a single lineage is the only stable solution.
Large effective mutation rates combined with large mutational jumps in
antigenic space lead to multiple stably co-existing lineages over prolonged
evolutionary periods. These results combined with observations from data
constrain the parameter regimes for the adaptation of viruses, including
influenza.


  The set of T cells that express the same T cell receptor (TCR) sequence
represents a T cell clone. The number of different naive T cell clones in an
organism reflects the number of different T cell receptors (TCRs) arising from
recombination of the V(D)J gene segments during T cell development in the
thymus. TCR diversity and more specifically, the clone abundance distribution
is an important factor in immune function. Specific recombination patterns
occur more frequently than others while subsequent interactions between TCRs
and self-antigens are known to trigger proliferation and sustain naive T cell
survival. These processes are TCR-dependent, leading to clone-dependent thymic
export and naive T cell proliferation rates. Using a mean-field approximation
to the solution of a regulated birth-death-immigration model and a modification
arising from sampling, we systematically quantify how TCR-dependent
heterogeneities in immigration and proliferation rates affect the shape of
clone abundance distributions (the number of different clones that are
represented by a specific number of cells, or "clone counts"). By comparing
predicted clone counts derived from our heterogeneous birth-death-immigration
model with experimentally sampled clone abundances, we show that although
heterogeneity in immigration rates causes very little change to predicted
clone-counts, significant heterogeneity in proliferation rates is necessary to
generate the observed abundances with reasonable physiological parameter
values. Our analysis provides constraints among physiological parameters that
are necessary to yield predictions that qualitatively match the data.
Assumptions of the model and potentially other important mechanistic factors
are discussed.


  Social discrimination seems to be a persistent phenomenon in many cultures.
It is important to understand the mechanisms that lead people to judge others
by the group to which they belong, rather than individual qualities. It was
recently shown that evolutionary (imitation) dynamics can lead to a
hierarchical discrimination between agents marked with observable, but
otherwise meaningless, labels. These findings suggest that it can give useful
insight, to describe the phenomenon of social discrimination in terms of
spontaneous symmetry breaking. The investigations so far have, however, only
considered binary labels. In this contribution we extend the investigations to
models with up to seven different labels. We find the features known from the
binary label model remain remarkably robust when the number of labels is
increased. We also discover a new feature, namely that it is more likely for
neighbours to have strategies which are similar, in the sense that they agree
on how to act towards a subset of the labels.


  Background Precise prediction of cancer types is vital for cancer diagnosis
and therapy. Important cancer marker genes can be inferred through predictive
model. Several studies have attempted to build machine learning models for this
task however none has taken into consideration the effects of tissue of origin
that can potentially bias the identification of cancer markers. Results In this
paper, we introduced several Convolutional Neural Network (CNN) models that
take unstructured gene expression inputs to classify tumor and non-tumor
samples into their designated cancer types or as normal. Based on different
designs of gene embeddings and convolution schemes, we implemented three CNN
models: 1D-CNN, 2D-Vanilla-CNN, and 2D-Hybrid-CNN. The models were trained and
tested on combined 10,340 samples of 33 cancer types and 731 matched normal
tissues of The Cancer Genome Atlas (TCGA). Our models achieved excellent
prediction accuracies (93.9-95.0%) among 34 classes (33 cancers and normal).
Furthermore, we interpreted one of the models, known as 1D-CNN model, with a
guided saliency technique and identified a total of 2,090 cancer markers (108
per class). The concordance of differential expression of these markers between
the cancer type they represent and others is confirmed. In breast cancer, for
instance, our model identified well-known markers, such as GATA3 and ESR1.
Finally, we extended the 1D-CNN model for prediction of breast cancer subtypes
and achieved an average accuracy of 88.42% among 5 subtypes. The codes can be
found at https://github.com/chenlabgccri/CancerTypePrediction.


  In single-molecule force spectroscopy experiments, a biomolecule is attached
to a force probe via polymer linkers, and the total extension -- of molecule
plus apparatus -- is monitored as a function of time. In a typical unfolding
experiment at constant force, the total extension jumps between two values that
correspond to the folded and unfolded states of the molecule. For several
biomolecular systems the committor, which is the probability to fold starting
from a given extension, has been used to extract the molecular activation
barrier (a technique known as "committor inversion"). In this work, we study
the influence of the force probe, which is much larger than the molecule being
measured, on the activation barrier obtained by committor inversion. We use a
two-dimensional framework in which the diffusion coefficient of the molecule
and of the pulling device can differ. We systematically study the free energy
profile along the total extension obtained from the committor, by numerically
solving the Onsager equation and using Brownian dynamics simulations. We
analyze the dependence of the extracted barrier on the linker stiffness,
molecular barrier height, and diffusion anisotropy, and thus, establish the
range of validity of committor inversion. Along the way, we showcase the
committor of 2-dimensional diffusive models and illustrate how it is affected
by barrier asymmetry and diffusion anisotropy.


  The ribosome is one of the largest and most complex macromolecular machines
in living cells. It polymerizes a protein in a step-by-step manner as directed
by the corresponding nucleotide sequence on the template messenger RNA (mRNA)
and this process is referred to as `translation' of the genetic message encoded
in the sequence of mRNA transcript. In each successful chemo-mechanical cycle
during the (protein) elongation stage, the ribosome elongates the protein by a
single subunit, called amino acid, and steps forward on the template mRNA by
three nucleotides called a codon. Therefore, a ribosome is also regarded as a
molecular motor for which the mRNA serves as the track, its step size is that
of a codon and two molecules of GTP and one molecule of ATP hydrolyzed in that
cycle serve as its fuel. What adds further complexity is the existence of
competing pathways leading to distinct cycles, branched pathways in each cycle
and futile consumption of fuel that leads neither to elongation of the nascent
protein nor forward stepping of the ribosome on its track. We investigate a
model formulated in terms of the network of discrete chemo-mechanical states of
a ribosome during the elongation stage of translation. The model is analyzed
using a combination of stochastic thermodynamic and kinetic analysis based on a
graph-theoretic approach. We derive the exact solution of the corresponding
master equations. We represent the steady state in terms of the cycles of the
underlying network and discuss the energy transduction processes. We identify
the various possible modes of operation of a ribosome in terms of its average
velocity and mean rate of GTP hydrolysis. We also compute entropy production as
functions of the rates of the interstate transitions and the thermodynamic cost
for accuracy of the translation process.


  Protein modeling is an increasingly popular area of machine learning
research. Semi-supervised learning has emerged as an important paradigm in
protein modeling due to the high cost of acquiring supervised protein labels,
but the current literature is fragmented when it comes to datasets and
standardized evaluation techniques. To facilitate progress in this field, we
introduce the Tasks Assessing Protein Embeddings (TAPE), a set of five
biologically relevant semi-supervised learning tasks spread across different
domains of protein biology. We curate tasks into specific training, validation,
and test splits to ensure that each task tests biologically relevant
generalization that transfers to real-life scenarios. We benchmark a range of
approaches to semi-supervised protein representation learning, which span
recent work as well as canonical sequence learning techniques. We find that
self-supervised pretraining is helpful for almost all models on all tasks, more
than doubling performance in some cases. Despite this increase, in several
cases features learned by self-supervised pretraining still lag behind features
extracted by state-of-the-art non-neural techniques. This gap in performance
suggests a huge opportunity for innovative architecture design and improved
modeling paradigms that better capture the signal in biological sequences. TAPE
will help the machine learning community focus effort on scientifically
relevant problems. Toward this end, all data and code used to run these
experiments are available at https://github.com/songlab-cal/tape.


  Monolayer-protected gold nanoparticles (Au NPs) are promising biomedical
tools with applications to diagnosis and therapy, thanks to their
biocompatibility and versatility. Here we show how the NP surface
functionalization can drive the mechanism of interaction with lipid membranes.
In particular, we show that the spontaneous protonation of anionic carboxylic
groups on the NP surface can make the NP-membrane interaction faster and less
disruptive.


  Non-Gaussian diffusion has been intensively studied in recent years, which
reflects the dynamic heterogeneity in the disordered media. The recent study on
the non-Gaussian diffusion in a static disordered landscape suggests novel
phenomena due to the quenched disorder. In this work, we further investigate
the random walk in this landscape under various effective temperature $\mu$,
which continuously modulates the dynamic heterogeneity. We show in the long
time limit, the trap dynamics on the landscape is equivalent to the quenched
trap model, in which sub-diffusion appears for $\mu<1$. The non-Gaussian
distribution of displacement has been analytically estimated for short $t$, of
which the stretched exponential tail is expected for $\mu\neq1$. Due to the
localization in the ensemble of trajectory segments, an additional peak arises
in $P(x,t)$ around $x=0$ even for $\mu>1$. Evolving in different time scales,
the peak and the tail of $P(x,t)$ are well split for a wide range of $t$. This
theoretical study reveals the connections among the sub-diffusion, non-Gaussian
diffusion, and the dynamic heterogeneity in the static disordered medium. It
also offers an insight on how the cell would benefit from the quasi-static
disordered structures.


  The first cell fate decisions in the developing fly embryo are made very
rapidly : hunchback genes decide in a few minutes whether a given nucleus
follows the anterior or the posterior developmental blueprint by reading out
the positional information encoded in the Bicoid morphogen. This developmental
system constitutes a prototypical instance of the broad spectrum of regulatory
decision processes that combine speed and accuracy. Traditional arguments based
on fixed-time sampling of Bicoid concentration indicate that an accurate
readout is not possible within the short times observed experimentally. This
raises the general issue of how speed-accuracy tradeoffs are achieved. Here, we
compare fixed-time sampling strategies to decisions made on-the-fly, which are
based on updating and comparing the likelihoods of being at an anterior or a
posterior location. We found that these more efficient schemes can complete
reliable cell fate decisions even within the very short embryological
timescales. We discuss the influence of promoter architectures on the mean
decision time and decision error rate and present concrete promoter
architectures that allow for the fast readout of the morphogen. Lastly, we
formulate explicit predictions for new experiments involving Bicoid mutants.


  An important dynamical property of biological interaction networks is
persistence, which intuitively means that "no species goes extinct". It has
been conjectured that dynamical system models of weakly reversible networks
(i.e., networks for which each reaction is part of a cycle) are persistent. The
property of persistence is also related to the well known global attractor
conjecture. An approach for the proof of the global attractor conjecture uses
an embedding of weakly reversible dynamical systems into toric differential
inclusions. We show that the larger class of endotactic dynamical systems can
also be embedded into toric differential inclusions. Moreover, we show that,
essentially, endotactic networks form the largest class of networks with this
property.


  The biological cell exhibits a fantastic range of behaviors, but ultimately
these are governed by a handful of physical and chemical principles. Here we
explore simple theory, known for decades and based on the simple thermodynamics
of mixtures of ideal gases, which illuminates several key functions performed
within the cell. Our focus is the free-energy-driven import and export of
molecules, such as nutrients and other vital compounds, via transporter
proteins. Complementary to a thermodynamic picture is a description of
transporters via "mass-action" chemical kinetics, which lends further insights
into biological machinery and free energy use. Both thermodynamic and kinetic
descriptions can shed light on the fundamental non-equilibrium aspects of
transport. On the whole, our biochemical-physics discussion will remain
agnostic to chemical details, but we will see how such details ultimately enter
a physical description through the example of the cellular fuel ATP.


  A small fraction of all protein structures characterized so far are
entangled. The challenge of understanding the properties of these knotted
proteins, and the why and the how of their natural folding process, has been
taken up in the past decade with different approaches, such as structural
characterization, in vitro experiments, and simulations of protein models with
varying levels of complexity. The simplest among these are the lattice G\=o
models, which belong to the class of structure-based models, i.e., models that
are biased to the native structure by explicitly including structural data. In
this review we highlight the contributions to the field made in the scope of
lattice G\=o models, putting them into perspective in the context of the main
experimental and theoretical results and of other, more realistic,
computational approaches.


  In recent decades, an increasing persistence of atmospheric circulation
patterns has been observed. In the course of the associated long-lasting
anticyclonic summer circulations, heat waves and drought spells often coincide,
leading to so-called hotter droughts. Previous hotter droughts caused a
decrease in agricultural yields and increase in tree mortality, and thus, had a
remarkable effect on carbon budgets and negative economic impacts.
Consequently, a quantification of ecosystem responses to hotter droughts and a
better understanding of the underlying mechanisms is crucial. In this context,
the European hotter drought of the year 2018 may be considered as a key event.
As a first step towards the quantification of its causes and consequences, we
here assess anomalies of atmospheric circulation patterns, temperature loads,
and climatic water balance as potential drivers of ecosystem responses as
quantified by remote sensing using the MODIS vegetation indices NDVI and EVI.
To place the drought of 2018 within a climatological context, we compare its
climatic features and ecosystem response with the extreme hot drought of 2003.
Our results indicated 2018 to be characterized by a climatic dipole, featuring
extremely hot and dry weather conditions north of the Alps but comparably cool
and moist conditions across large parts of the Mediterranean. Analyzing
ecosystem response of five dominant land-cover classes, we found significant
positive effects of April-July climatic water balance on ecosystem
productivity. Negative drought impacts appeared to affect a larger area in 2018
compared to 2003. We found a significantly higher sensitivity of pastures and
arable land to climatic water balance compared to forests in both years. This
study quantifies the drought of 2018 as a yet unprecedented event and provides
valuable insights into the heterogeneous drought responses of European
ecosystems.


  We develop a multiple compartment Susceptible-Infected-Recovered (SIR) model
to analyze the spread of several infectious diseases through different
geographic areas. Additionally, we propose a data-quality sensitive
optimization framework for fitting this model to observed data.
  We fit the model to the temporal profile of the number of people infected by
one of six influenza strains in Europe over $7$ influenza seasons. In addition
to describing the temporal and spatial spread of influenza, the model provides
an estimate of the inter-country and intra-country infection and recovery rates
of each strain and in each season. We find that disease parameters remain
relatively stable, with a correlation greater than $0.5$ over seasons and
stains. Clustering of influenza strains by the inferred disease parameters is
consistent with genome sub-types. Surprisingly, our analysis suggests that
inter-country human mobility plays a negligible role in the spread of influenza
in Europe. Finally, we show that the model allows the estimation of disease
load in countries with poor or none existent data from the disease load in
adjacent countries.
  Our findings reveal information on the spreading mechanism of influenza and
on disease parameters. These can be used to assist in disease surveillance and
in control of influenza as well as of other infectious pathogens in a
heterogenic environment.


  Cells estimate concentrations of chemical ligands in their environment using
a limited set of receptors. Recent work has shown that the temporal sequence of
binding and unbinding events on just a single receptor can be used to estimate
the concentrations of multiple ligands. Here, for a network of many ligands and
many receptors, we show that such temporal sequences can be used to estimate
the concentration of a few times as many ligand species as there are receptors.
Crucially, we show that the spectrum of the inverse covariance matrix of these
estimates has several universal properties, which we trace to properties of
Vandermonde matrices. We argue that this can be used by cells in realistic
biochemical decoding networks.


  Gene expression is a random or noisy process. The process consists of several
random events among which the reinitiation of transcription by RNAP is an
important one. The RNAP molecules can bind the gene only after the promoter
gets activated by transcription factors. Several transcription factors bind the
promoter to put the gene in the active state. The gene turns into inactive
state as the bound transcription factors leave the promoter. During the active
period of the gene, many RNAP molecules transcribe the gene to synthesize the
mRNAs. The binding event of RNAP to the active state of the gene is a
probabilistic process and therefore, introduces noise or fluctuations in the
mRNA and protein levels. In this paper, we analytically calculate the Fano
factor in mRNA and protein levels and also the probability distribution of mRNA
numbers exactly with the binding event of RNAPs in gene transcription process.
The analytically calculated expression of Fano factor of proteins shows
excellent agreement with an experimental result. Then we show that the Fano
factor in mRNA levels can be sub-Poissonian due to the reinitiation of
transcription by RNAP and the mean mRNA level can be increased without
increasing the Fano factor. Our study show that the Fano factor can also be
reduced keeping mRNA levels fixed. We find that the reinitiation of
transcription can behave as a fine-tuned control process to regulate the
mRNA/protein level in the cell.