  Interest in equilibrium-based sampling methods has grown with recent advances
in computational hardware and Markov state modeling (MSM) methods, yet
outstanding questions remain that hinder widespread adoption. Namely, how do
sampling strategies explore conformational space and how might this influence
predictions? Here, we seek to answer these questions for four commonly used
sampling methods: 1) a long simulation, 2) many short simulations, 3) adaptive
sampling, and 4) FAST. We first develop a theoretical framework for
analytically calculating the probability of discovering states and uncover the
drastic effects of varying the number and length of simulations. We then use
kinetic Monte Carlo simulations on a variety of physically inspired landscapes
to characterize state discovery and transition pathways. Consistently, we find
that FAST simulations discover target states with the highest probability and
traverse realistic pathways. Furthermore, we uncover the pathology that short
parallel simulations sometimes predict an incorrect transition pathway by
crossing large energy barriers that long simulations would typically
circumnavigate, which we refer to as pathway tunneling. To protect against
tunneling, we introduce FAST-string, which samples along the highest-flux
transition paths to refine an MSMs transition probabilities and discriminate
between competing pathways. Additionally, we compare MSM estimators in
describing thermodynamics and kinetics. For adaptive sampling, we recommend
normalizing the transition counts out of each state after adding pseudo-counts
to avoid creating sources or sinks. Lastly, we evaluate our insights from
simple landscapes with all-atom molecular dynamics simulations of the folding
of the {\lambda}-repressor protein. We find that FAST-contacts predicts the
same folding pathway as long simulations but with orders of magnitude less
simulation time.


  The mechanical properties of the cell depend crucially on the tension of its
cytoskeleton, a biopolymer network that is put under stress by active motor
proteins. While the fibrous nature of the network is known to strongly affect
the transmission of these forces to the cellular scale, our understanding of
this process remains incomplete. Here we investigate the transmission of forces
through the network at the individual filament level, and show that active
forces can be geometrically amplified as a transverse motor-generated force
force "plucks" the fiber and induces a nonlinear tension. In stiff and densely
connnected networks, this tension results in large network-wide tensile
stresses that far exceed the expectation drawn from a linear elastic theory.
This amplification mechanism competes with a recently characterized
network-level amplification due to fiber buckling, suggesting that that fiber
networks provide several distinct pathways for living systems to amplify their
molecular forces.


  Radix Bupleuri is a hepatoprotective traditional Chinese medicine (TCM) used
for thousands of years in clinical, which was reported to be linked with liver
damage. Previous studies have revealed that saikosaponins are the major types
of components that contribute to the hepatotoxicity of Radix Bupleuri. However
the underlying molecular mechanism is far from being understood. In order to
clarify whether these structural analogues exert toxicity effects through the
same molecular targets, a systematic comparison study was done in this paper.
The effects of SSa, b2, c, and d on isolated rat liver mitochondria and human
hepatocyte L02 cells were explored, respectively. The collective results
indicated that although saikosaponins share the similar structures but they
have quite different mechanisms. SSb2 and SSd showed most serious damage on the
function of mitochondria and survival rate of cell, respectively. SSb2 could
cause mitochondrial permeability transition pore (mPTP) opening and collapse of
mitochondrial membrane potential ({\Delta}{\Psi}m) by impairing the
mitochondrial respiratory chain complex III. While SSd destroyed plasma
membrane and led to the release of lactate dehydrogenase (LDH) mainly through
activating caspase-1. Furthermore, the combine index (CI) demonstrated that the
combined hepatotoxicity of SSb2 and SSd could be additive. This finding might
yield more in depth understanding of hepatotoxicity of Radix Bupleuri possess
many different saikosaponins.


  The outbreak of an infectious disease in a human population can lead to
individuals responding with preventive measures in an attempt to avoid getting
infected. This leads to changes in contact patterns. However, as we show in
this paper, rational behaviour at the individual level, such as social
distancing from infectious contacts, may not always be beneficial for the
population as a whole. We use epidemic network models to demonstrate the
potential negative consequences at the population level. We take into account
the social structure of the population through several network models. As the
epidemic evolves, susceptible individuals may distance themselves from their
infectious contacts. Some individuals replace their lost social connections by
seeking new ties. We show that social distancing can worsen the disease outcome
both in the initial phase of an outbreak and the final epidemic size. Moreover,
the same negative effect can arise in real-world networks. Our results suggest
that one needs to be careful when targeting behavioural changes as they could
potentially worsen the epidemic outcome. Furthermore, network structure
crucially influences the way that individual-level measures impact the epidemic
at the population level. These findings highlight the importance of careful
analysis of preventive measures in epidemic models.


  A new detailed mathematical model for dynamics of immune response to
hepatitis B is proposed, which takes into account contributions from innate and
adaptive immune responses, as well as cytokines. Stability analysis of
different steady states is performed to identify parameter regions where the
model exhibits clearance of infection, maintenance of a chronic infection, or
periodic oscillations. Effects of nucleoside analogues and interferon
treatments are analysed, and the critical drug efficiency is determined.


  Gene duplication is a major mechanism through which new genetic material is
generated. Although numerous methods have been developed to differentiate the
ortholog and paralogs, very few differentiate the "Parent-Daughter"
relationship among paralogous pairs. As coined by the Mira et al, we refer the
"Parent" copy as the paralogous copy that stays at the original genomic
position of the "original copy" before the duplication event, while the
"Daughter" copy occupies a new genomic locus.
  Here we present a novel method which combines the phylogenetic reconstruction
of duplications at different evolutionary periods and the synteny evidence
collected from the preserved homologous gene orders. We reconstructed for the
first time a deep evolutionary history of "Parent-Daughter" relationships among
genes that were descendants from 2 rounds of whole genome duplications (2R
WGDs) at early vertebrates and were further duplicated in later ceancestors
like early Mammalia and early Primates.
  Our analysis reveals that the "Parent" copy has significantly fewer
accumulated mutations compared with the "Daughter" copy since their divergence
after the duplication event. More strikingly, we found that the "Parent" copy
in a duplication event continues to be the "Parent" of the younger successive
duplication events which lead to "grand-daughters".


  The evolution of the genome has led to very sophisticated and complex
regulation. Because of the abundance of non-coding RNA (ncRNA) in the cell,
different species will promiscuously associate with each other, suggesting
collective dynamics similar to artificial neural networks. Here we present a
simple mechanism allowing ncRNA to perform computations equivalent to neural
network algorithms such as Boltzmann machines and the Hopfield model. The
quantities analogous to the neural couplings are the equilibrium constants
between different RNA species. The relatively rapid equilibration of RNA
binding and unbinding is regulated by a slower process that degrades and
creates new RNA. The model requires that the creation rate for each species be
an increasing function of the ratio of total to unbound RNA. Similar mechanisms
have already been found to exist experimentally for ncRNA regulation. With the
overall concentration of RNA regulated, equilibrium constants can be chosen to
store many different patterns, or many different input-output relations. The
network is also quite insensitive to random mutations in equilibrium constants.
Therefore one expects that this kind of mechanism will have a much higher
mutation rate than ones typically regarded as being under evolutionary
constraint.


  Explaining biodiversity in nature is a fundamental problem in ecology. An
outstanding challenge is embodied in the so-called Competitive Exclusion
Principle: two species competing for one limiting resource cannot coexist at
constant population densities, or more generally, the number of consumer
species in steady coexistence cannot exceed that of resources. The fact that
competitive exclusion is rarely observed in natural ecosystems has not been
fully understood. Here we show that by forming chasing triplets among the
consumers and resources in the consumption process, the Competitive Exclusion
Principle can be naturally violated. The modeling framework developed here is
broadly applicable and can be used to explain the biodiversity of many
consumer-resource ecosystems and hence deepens our understanding of
biodiversity in nature.


  Our understanding of evolution is shaped strongly by how we conceive of its
fundamental causes. In the original Modern Synthesis, evolution was defined as
a process of shifting the frequencies of available alleles at many loci
affecting a trait under selection. Events of mutation that introduce novelty
were not considered evolutionary causes, but proximate causes acting at the
wrong level. Today it is clear that long-term evolutionary dynamics depend on
the dynamics of mutational introduction. Yet, the implications of this
dependency remain unfamiliar, and have not yet penetrated into high-level
debates over evolutionary theory. Modeling the influence of biases in the
introduction process reveals behavior previously unimagined, as well as
behavior previously considered impossible. Quantitative biases in the
introduction of variation can impose biases on the outcome of evolution without
requiring high mutation rates or neutral evolution. Mutation-biased adaptation,
a possibility not previously imagined, has been observed among diverse taxa.
Directional trends are possible under a sustained bias. Biases that are
developmental in origin may have an effect analogous to mutational biases.
Structuralist arguments invoking the relative accessibility of forms in
state-space can be understood as references to the role of biases in the
introduction of variation. That is, the characteristic concerns of molecular
evolution, evo-devo and structuralism can be interpreted to implicate a kind of
causation absent from the original Modern Synthesis.


  The spread of infectious diseases, rumors, fashions, innovations are complex
contagion processes, embedded both in networked and spatial contexts. Here we
investigate the pattern dynamics of a complex contagion, where two agents, say
$A$ and $B$, interact with each other and diffuse simultaneously in the
geographic space. The contagion process for each follows the classical
susceptible-infected-susceptible kinetics, and their interaction introduces a
potential change in the secondary infection propensity compared to the baseline
reproduction ratio $R_0$. We show that nontrivial spatial infection patterns
arise, when the susceptible move faster than the infected and the interaction
between the two agents is neither too competitive nor too cooperative.
Interestingly, the system exhibits pattern hysteresis phenomena that quite
different parameter regions allowing for patterns exist in the direction of
increasing $R_0$ and in the direction of eradication by its reduction. The
latter shows a remarkable enhancement in the contagion prevalence, meaning that
the infection eradication now becomes extremely difficult compared to the
single-agent scenario and to the coinfection without space. Linearization
analysis supports our observations, and we identified the required elements and
dynamical mechanism behind the emergence of a pattern. These findings call for
further investigation for their close relevance, both in biological and social
contagions.


  We propose a novel method for numerical modeling of spatially inhomogeneous
moment dynamics of populations with nonlocal dispersal and competition in
continuous space. It is based on analytically solvable decompositions of the
time evolution operator for a coupled set of master equations. This has allowed
us -- for the first time in the literature -- to perform moment dynamics
simulations of spatially inhomogeneous systems beyond the mean-field approach
and to calculate the inhomogeneous pair correlation function using the Kirkwood
superposition ansatz. As a result, we revealed a number of new subtle effects,
possible in real populations. Namely, for systems with short-range dispersal
and mid-range competition, strong clustering of entities at small distances
followed by their deep disaggregation at larger separations are observed in the
wavefront of density propagation. For populations in which the competition
range is much shorter than that of dispersal, the pair correlation function
exhibits a long-tail behavior. Remarkably, the latter effect takes place only
due to the spatial inhomogeneity and thus was completely unknown before.
Moreover, both effects get stronger in the direction of propagation. All these
types of behavior are interpreted as a trade-off between the dispersal and
competition in the coexistence of reproductive pair correlations and the
inhomogeneity of the density of the system.


  Experiments measuring currents through single protein channels show unstable
currents, a phenomena called the gating of a single channel. Channels switch
between an 'open' state with a well defined single amplitude of current and
'closed' states with nearly zero current. The existing mean-field theory of ion
channels focuses almost solely on the open state. The physical modeling of the
dynamical features of ion channels is still in its infancy, and does not
describe the transitions between open and closed states, nor the distribution
of the duration times of open states. One hypothesis is that gating corresponds
to noise-induced fast transitions between multiple steady (equilibrium) states
of the underlying system. In this work, we aim to test this hypothesis.
Particularly, our study focuses on the (high order) steric
Poisson-Nernst-Planck-Cahn-Hilliard model since it has been successful in
predicting permeability and selectivity of ionic channels in their open state,
and since it gives rise to multiple steady states. We show that this system
gives rise to a gating-like behavior, but that important features of this
switching behavior are different from the defining features of gating in
biological systems. Furthermore, we show that noise prohibits switching in the
system of study. The above phenomena are expected to occur in other PNP-type
models, strongly suggesting that one has to go beyond over-damped (gradient
flow) Nernst-Planck type dynamics to explain the spontaneous gating of single
channels.


  Recent genetic studies and whole-genome sequencing projects have greatly
improved our understanding of human variation and clinically actionable genetic
information. Smaller ethnic populations, however, remain underrepresented in
both individual and large-scale sequencing efforts and hence present an
opportunity to discover new variants of biomedical and demographic
significance. This report describes the sequencing and analysis of a genome
obtained from an individual of Serbian origin, introducing tens of thousands of
previously unknown variants to the currently available pool. Ancestry analysis
places this individual in close proximity of the Central and Eastern European
populations; i.e., closest to Croatian, Bulgarian and Hungarian individuals
and, in terms of other Europeans, furthest from Ashkenazi Jewish, Spanish,
Sicilian, and Baltic individuals. Our analysis confirmed gene flow between
Neanderthal and ancestral pan-European populations, with similar contributions
to the Serbian genome as those observed in other European groups. Finally, to
assess the burden of potentially disease-causing/clinically relevant variation
in the sequenced genome, we utilized manually curated genotype-phenotype
association databases and variant-effect predictors. We identified several
variants that have previously been associated with severe early-onset disease
that is not evident in the proband, as well as variants that could yet prove to
be clinically relevant to the proband over the next decades. The presence of
numerous private and low-frequency variants along with the observed and
predicted disease-causing mutations in this genome exemplify some of the global
challenges of genome interpretation, especially in the context of understudied
ethnic groups.


  The predictive ability of stochastic chemical reactions is currently limited
by the lack of closed form solutions to the governing chemical master equation.
To overcome this limitation, this paper proposes a computational method capable
of predicting mathematically rigorous upper and lower bounds of transient
moments for reactions governed by the law of mass action. We first derive an
equation that transient moments must satisfy based on the moment equation.
Although this equation is underdetermined, we introduce a set of semidefinite
constraints known as moment condition to narrow the feasible set of the
variables in the equation. Using these conditions, we formulate a semidefinite
program that efficiently and rigorously computes the bounds of transient moment
dynamics. The proposed method is demonstrated with illustrative numerical
examples and is compared with related works to discuss advantages and
limitations.


  Markov jump process models have many applications across science. Often,
these models are defined on a state-space of product form and only one of the
components of the process is of direct interest. In this paper, we extend the
marginal process framework, which provides a marginal description of the
component of interest, to the case of fully coupled processes. We use entropic
matching to obtain a finite-dimensional approximation of the filtering
equation, which governs the transition rates of the marginal process. The
resulting equations can be seen as a combination of two projection operations
applied to the full master equation, so that we obtain a principled model
reduction framework. We demonstrate the resulting reduced description on the
totally asymmetric exclusion process. An important class of Markov jump
processes are stochastic reaction networks, which have applications in chemical
and biomolecular kinetics, ecological models and models of social networks. We
obtain a particularly simple instantiation of the marginal process framework
for mass-action systems by using product-Poisson distributions for the
approximate solution of the filtering equation. We investigate the resulting
approximate marginal process analytically and numerically.


  Data is scaling exponentially in fields ranging from genomics to neuroscience
to economics. A central question is: can modern machine learning methods be
applied to construct predictive models of natural systems like cells and brains
based on large data sets? In this paper, we examine how inference is impacted
when training data is generated by the statistical behavior of a physical
system, and hence outside direct control by the experimentalist. We develop an
information-theoretic analysis for the canonical problem of spin-network
inference. Our analysis reveals the essential role that the physical properties
of the spin network and its environment play in determining the difficulty of
the underlying machine learning problem. Specifically, stochastic fluctuations
drive a system to explore a range of configurations providing `raw' information
for a learning algorithm to construct an accurate model; yet they also blur
energetic differences between network states and thereby degrade information.
This competition leads spin networks to generically have an intrinsic optimal
temperature at which stochastic spin fluctuations provide maximal information
for discriminating among competing models, maximizing inference efficiency. We
demonstrate a simple active learning protocol that optimizes network
temperature to boost inference efficiency and dramatically increases the
efficiency of inference on a neural circuit reconstruction task. Our results
reveal a fundamental link between physics and information and show how the
physical environment can be tuned to optimize the efficiency of machine
learning.


  The study of high-throughput genomic profiles from a pharmacogenomics
viewpoint has provided unprecedented insights into the oncogenic features
modulating drug response. A recent screening of ~1,000 cancer cell lines to a
collection of anti-cancer drugs illuminated the link between genotypes and
vulnerability. However, due to essential differences between cell lines and
tumors, the translation into predicting drug response in tumors remains
challenging. Here we proposed a DNN model to predict drug response based on
mutation and expression profiles of a cancer cell or a tumor. The model
contains a mutation and an expression encoders pre-trained using a large
pan-cancer dataset to abstract core representations of high-dimension data,
followed by a drug response predictor network. Given a pair of mutation and
expression profiles, the model predicts IC50 values of 265 drugs. We trained
and tested the model on a dataset of 622 cancer cell lines and achieved an
overall prediction performance of mean squared error at 1.96 (log-scale IC50
values). The performance was superior in prediction error or stability than two
classical methods and four analog DNNs of our model. We then applied the model
to predict drug response of 9,059 tumors of 33 cancer types. The model
predicted both known, including EGFR inhibitors in non-small cell lung cancer
and tamoxifen in ER+ breast cancer, and novel drug targets. The comprehensive
analysis further revealed the molecular mechanisms underlying the resistance to
a chemotherapeutic drug docetaxel in a pan-cancer setting and the anti-cancer
potential of a novel agent, CX-5461, in treating gliomas and hematopoietic
malignancies. Overall, our model and findings improve the prediction of drug
response and the identification of novel therapeutic options.


  A fundamental question in biology is how cell populations evolve into
different subtypes based on homogeneous processes at the single cell level.
Here we show that population bimodality can emerge even when biological
processes are homogenous at the cell level and the environment is kept
constant. Our model is based on the stochastic partitioning of a cell component
with an optimal copy number. We show that the existence of unimodal or bimodal
distributions depends on the variance of partition errors and the growth rate
tolerance around the optimal copy number. In particular, our theory provides a
consistent explanation for the maintenance of aneuploid states in a population.
The proposed model can also be relevant for other cell components such as
mitochondria and plasmids, whose abundances affect the growth rate and are
subject to stochastic partition at cell division.


  Phylogenetic networks are rooted directed acyclic graphs that represent
evolutionary relationships between species whose past includes reticulation
events such as hybridisation and horizontal gene transfer. To search the space
of phylogenetic networks, the popular tree rearrangement operation rooted
subtree prune and regraft (rSPR) was recently generalised to phylogenetic
networks. This new operation - called subnet prune and regraft (SNPR) - induces
a metric on the space of all phylogenetic networks as well as on several
widely-used network classes. In this paper, we investigate several problems
that arise in the context of computing the SNPR-distance. For a phylogenetic
tree $T$ and a phylogenetic network $N$, we show how this distance can be
computed by considering the set of trees that are embedded in $N$ and then use
this result to characterise the SNPR-distance between $T$ and $N$ in terms of
agreement forests. Furthermore, we analyse properties of shortest
SNPR-sequences between two phylogenetic networks $N$ and $N'$, and answer the
question whether or not any of the classes of tree-child, reticulation-visible,
or tree-based networks isometrically embeds into the class of all phylogenetic
networks under SNPR.


  Bioinformatics tools have been developed to interpret gene expression data at
the gene set level, and these gene set based analyses improve the biologists'
capability to discover functional relevance of their experiment design. While
elucidating gene set individually, inter gene sets association is rarely taken
into consideration. Deep learning, an emerging machine learning technique in
computational biology, can be used to generate an unbiased combination of gene
set, and to determine the biological relevance and analysis consistency of
these combining gene sets by leveraging large genomic data sets. In this study,
we proposed a gene superset autoencoder (GSAE), a multi-layer autoencoder model
with the incorporation of a priori defined gene sets that retain the crucial
biological features in the latent layer. We introduced the concept of the gene
superset, an unbiased combination of gene sets with weights trained by the
autoencoder, where each node in the latent layer is a superset. Trained with
genomic data from TCGA and evaluated with their accompanying clinical
parameters, we showed gene supersets' ability of discriminating tumor subtypes
and their prognostic capability. We further demonstrated the biological
relevance of the top component gene sets in the significant supersets. Using
autoencoder model and gene superset at its latent layer, we demonstrated that
gene supersets retain sufficient biological information with respect to tumor
subtypes and clinical prognostic significance. Superset also provides high
reproducibility on survival analysis and accurate prediction for cancer
subtypes.


  Differences between men and women have intrigued generations of social
scientists, who have found that the two sexes behave differently in settings
requiring competition, risk taking, altruism, honesty, as well as many others.
Yet, little is known about whether there are gender differences in cooperative
behavior. Previous evidence is mixed and inconclusive. Here I shed light on
this topic by analyzing the totality of studies that my research group has
conducted since 2013. This is a dataset of 10,951 observations coming from
7,322 men and women living in the US, recruited through Amazon Mechanical Turk,
and who passed four comprehension questions to make sure they understand the
cooperation problem (a one-shot prisoner's dilemma). The analysis demonstrates
that women are more cooperative than men. The effect size is small (about 4
percentage points, and this might explain why previous studies failed to detect
it) but highly significant (p<.0001).


  This preprint has been reviewed and recommended by Peer Community In Ecology
(https://dx.doi.org/10.24072/pci.ecology.100001). In a recent article,
Bicca-Marques and Calegaro-Marques [2016] discussed the putative assumptions
related to an interpretation we provided regarding an observed positive
relationship between weekly averaged parasite richness of a group of mandrills
(Mandrillus sphinx) and their daily path lengths (DPL), published earlier in
the same journal (Brockmeyer et al., 2015). In our article, we proposed, inter
alia, that "the daily travels of mandrills could be seen as a way to escape
contaminated habitats on a local scale". In their article, Bicca-Marques and
Calegaro-Marques [2016] proposed an alternative mechanism that they considered
to be more parsimonious. In their view, increased DPL also increases exposure
to novel parasites from the environment. In other words, while we proposed that
elevated DPL may be a consequence of elevated parasite richness, they viewed it
as a cause. We are happy to see that our study attracted so much interest that
it evoked a public comment. We are also grateful to Bicca-Marques and
Calegaro-Marques [2016] for pointing out an obvious alternative scenario that
we failed to discuss and for laying out several key factors and assumptions
that should be addressed by future studies examining the links between parasite
risk and group ranging. We use this opportunity to advance this discourse by
responding to some of the criticisms raised in their discussion of our article.
In this reply, we briefly contextualize the main object of criticism. We then
discuss the putative parsimony of the two competing scenarios.


  Previous experiments have explored the effects of gender and cognitive
reflection on dishonesty separately. To the best of our knowledge, no studies
have investigated potential interactions between these two factors. Exploring
this interaction is important because previous work found that males tend to be
both more deliberative than females. Therefore, it is possible that the gender
effect on dishonesty is moderated by cognitive reflection. Here we report a
large online experiment (N = 766) where subjects first have a chance to lie for
their benefit and then take a Cognitive Reflection Test (CRT). We find a
significant interaction between gender and CRT score such that lack of
deliberation promotes honesty for men but not for women. Additional analyses
highlight that this effect is mainly driven by men whose answers in the CRT are
neither intuitive nor deliberative, who happen to be particularly honest in our
deception game.


  Learning language of protein sequences, which captures non-local interactions
between amino acids close in the spatial structure, is a long-standing
bioinformatics challenge, which requires at least context-free grammars.
However, complex character of protein interactions impedes unsupervised
learning of context-free grammars. Using structural information to constrain
the syntactic trees proved effective in learning probabilistic natural and RNA
languages. In this work, we establish a framework for learning probabilistic
context-free grammars for protein sequences from syntactic trees partially
constrained using amino acid contacts obtained from wet experiments or
computational predictions, whose reliability has substantially increased
recently. Within the framework, we implement the maximum-likelihood and
contrastive estimators of parameters for simple yet practical grammars. Tested
on samples of protein motifs, grammars developed within the framework showed
improved precision in recognition and higher fidelity to protein structures.
The framework is applicable to other biomolecular languages and beyond wherever
knowledge of non-local dependencies is available.


  Parameterizing the approximate posterior of a generative model with neural
networks has become a common theme in recent machine learning research. While
providing appealing flexibility, this approach makes it difficult to impose or
assess structural constraints such as conditional independence. We propose a
framework for learning representations that relies on Auto-Encoding Variational
Bayes and whose search space is constrained via kernel-based measures of
independence. In particular, our method employs the $d$-variable
Hilbert-Schmidt Independence Criterion (dHSIC) to enforce independence between
the latent representations and arbitrary nuisance factors. We show how to apply
this method to a range of problems, including the problems of learning
invariant representations and the learning of interpretable representations. We
also present a full-fledged application to single-cell RNA sequencing
(scRNA-seq). In this setting the biological signal is mixed in complex ways
with sequencing errors and sampling effects. We show that our method
out-performs the state-of-the-art in this domain.


  Genome-to-genome comparisons require designating anchor points, which are
given by Maximum Exact Matches (MEMs) between their sequences. For large
genomes this is a challenging problem and the performance of existing
solutions, even in parallel regimes, is not quite satisfactory. We present a
new algorithm, copMEM, that allows to sparsely sample both input genomes, with
sampling steps being coprime. Despite being a single-threaded implementation,
copMEM computes all MEMs of minimum length 100 between the human and mouse
genomes in less than 2 minutes, using less than 10 GB of RAM memory.


  With recent developments in remote sensing technologies, plot-level forest
resources can be predicted utilizing airborne laser scanning (ALS). The
prediction is often assisted by mostly vertical summaries of the ALS point
clouds. We present a spatial analysis of the point cloud by studying the
horizontal distribution of the pulse returns through canopy height models
thresholded at different height levels. The resulting patterns of patches of
vegetation and gabs on each layer are summarized to spatial ALS features. We
propose new features based on the Euler number, which is the number of patches
minus the number of gaps, and the empty-space function, which is a spatial
summary function of the gab space. The empty-space function is also used to
describe differences in the gab structure between two different layers. We
illustrate usefulness of the proposed spatial features for predicting different
forest variables that summarize the spatial structure of forests or their
breast height diameter distribution. We employ the proposed spatial features,
in addition to commonly used features from literature, in the well-known k-nn
estimation method to predict the forest variables. We present the methodology
on the example of a study site in Central Finland.


  In the present work, a non-linear stochastic model is presented to study the
effect of time variation of transmission rates on the co-evolution of epidemics
and its corresponding awareness over a two layered multiplex network. In this
model, the infection transmission rate of a given node in the epidemic layer
depends upon its awareness probability in the awareness layer. Similarly, the
infection information transmission rate of a node in the awareness layer
depends upon its infection probability in the epidemic layer. The spread of
disease resulting from physical contacts is described in terms of SIS
(Susceptible Infected Susceptible) process over the epidemic layer and the
spread of information about the disease outbreak is described in terms of UAU
(Unaware Aware Unaware) process over the virtual interaction mediated awareness
layer. The time variation of the transmission rates and the resulting
co-evolution of these mutually competing processes is studied in terms of a
network topology depend parameter({\alpha}). Using a second order linear theory
it has been shown that in the continuous time limit, the co-evolution of these
processes can be described in terms of damped and driven harmonic oscillator
equations. From the results of the Monte-Carlo simulation, it is shown that for
the suitable choice of parameter({\alpha}), the two process can either exhibit
sustained oscillatory or damped dynamics. The damped dynamics corresponds to
the endemic state. Further, for the case of endemic state it is shown that the
inclusion of awareness layer significantly lowers the disease transmission rate
and reduces the size of epidemic. The endemic state infection probability of a
given node corresponding to the damped dynamics is found to have dependence
upon both the transmission rates as well as on both absolute intra-layer and
relative inter-layer degree of the individual nodes.


  Proteins employ the information stored in the genetic code and translated
into their sequences to carry out well-defined functions in the cellular
environment. The possibility to encode for such functions is controlled by the
balance between the amount of information supplied by the sequence and that
left after that the protein has folded into its structure. We developed a
computational algorithm to evaluate the amount of information necessary to
specify the protein structure, keeping into account the thermodynamic
properties of protein folding. We thus show that the information remaining in
the protein sequence after encoding for its structure (the 'information gap')
is very close to what needed to encode for its function and interactions. Then,
by predicting the information gap directly from the protein sequence, we show
that it may be possible to use these insights from information theory to
discriminate between ordered and disordered proteins, to identify unknown
functions, and to optimize designed proteins sequences.


  Recent studies suggest that the emergence of cooperative behavior can be
explained by generalized reciprocity, a behavioral mechanism based on the
principle of "help anyone if helped by someone". In complex systems, the
cooperative dynamics is largely determined by the network structure which
dictates the interactions among neighboring individuals. These interactions
often exhibit multidimensional features, either as relationships of different
types or temporal dynamics, both of which may be modeled as a "multiplex"
network. Against this background, here we advance the research on cooperation
models inspired by generalized reciprocity by considering a multidimensional
networked society. Our results reveal that a multiplex network structure may
enhance the role of generalized reciprocity in promoting cooperation, whereby
some of the network dimensions act as a latent support for the others. As a
result, generalized reciprocity forces the cooperative contributions of the
individuals to concentrate in the dimension which is most favorable for the
existence of cooperation.


  In multi-agent systems, cooperative behavior is largely determined by the
network structure which dictates the interactions among neighboring agents.
These interactions often exhibit multidimensional features, either as
relationships of different types or temporal dynamics, both of which may be
modeled as a "multiplex" network. Against this background, here we advance the
research on cooperation models inspired by generalized reciprocity, a simple
pay-it-forward behavioral mechanism, by considering a multidimensional
networked society. Our results reveal that a multiplex network structure can
act as an enhancer of the role of generalized reciprocity in promoting
cooperation by acting as a latent support, even when the parameters in some of
the separate network dimensions suggest otherwise (i.e. favor defection). As a
result, generalized reciprocity forces the cooperative contributions of the
individual agents to concentrate in the dimension which is most favorable for
the existence of cooperation.


  Environmental contaminant exposure can pose significant risks to human
health. Therefore, evaluating the impact of this exposure is of great
importance; however, it is often difficult because both the molecular mechanism
of disease and the mode of action of the contaminants are complex. We used
network biology techniques to quantitatively assess the impact of environmental
contaminants on the human interactome and diseases with a particular focus on
seven major contaminant categories: persistent organic pollutants (POPs),
dioxins, polycyclic aromatic hydrocarbons (PAHs), pesticides,
perfluorochemicals (PFCs), metals, and pharmaceutical and personal care
products (PPCPs). We integrated publicly available data on toxicogenomics, the
diseasome, protein-protein interactions (PPIs), and gene essentiality and found
that a few contaminants were targeted to many genes, and a few genes were
targeted by many contaminants. The contaminant targets were hub proteins in the
human PPI network, whereas the target proteins in most categories did not
contain abundant essential proteins. Generally, contaminant targets and
disease-associated proteins were closely associated with the PPI network, and
the closeness of the associations depended on the disease type and chemical
category. Network biology techniques were used to identify environmental
contaminants with broad effects on the human interactome and
contaminant-sensitive biomarkers. Moreover, this method enabled us to quantify
the relationship between environmental contaminants and human diseases, which
was supported by epidemiological and experimental evidence. These methods and
findings have facilitated the elucidation of the complex relationship between
environmental exposure and adverse health outcomes.


  The Protein Data Bank (PDB) contains more than 135 000 entries today. From
these, relatively few amyloid structures can be identified, since amyloids are
insoluble in water. Therefore, mostly solid state NMR-recorded amyloid
structures are deposited in the PDB. Based on the geometric analysis of these
deposited structures we have prepared an automatically updated webserver, which
generates the list of the deposited amyloid structures, and, additionally,
those globular protein entries, which have amyloid-like substructures of a
given size and characteristics. We have found that applying only the properly
chosen geometric conditions, it is possible to identify the deposited amyloid
structures, and a number of globular proteins with amyloid-like substructures.
We have analyzed these globular proteins and have found that many of them are
known to form amyloids more easily than many other globular proteins. Our
results relate to the method of (Stankovic, I. et al. (2017): Construction of
Amyloid PDB Files Database. Transactions on Internet Research. 13 (1): 47-51),
who have applied a hybrid textual-search and geometric approach for finding
amyloids in the PDB.
  If one intends to identify a subset of the PDB for some applications, the
identification algorithm needs to be re-run periodically, since in 2017, on
average, every day 30 new entries were deposited in the data bank. Our
webserver is updated regularly and automatically, and the identified amyloid-
and partial amyloid structures can be viewed or their list can be downloaded
from the site https://pitgroup.org/amyloid.


  Previous experiments have found mixed results on whether honesty is intuitive
or requires deliberation. Here we add to this literature by building on prior
work of Capraro (2017). We report a large study (N=1,389) manipulating time
pressure vs time delay in a deception game. We find that, in this setting,
people are more honest under time pressure, and that this result is not driven
by confounds present in earlier work.


  In 1932, Paul Erdos asked whether a random walk constructed from a binary
sequence can achieve the lowest possible deviation (lowest discrepancy), for
the sequence itself and for all its subsequences formed by homogeneous
arithmetic progressions. Although avoiding low discrepancy is impossible for
infinite sequences, as recently proven by Terence Tao, attempts were made to
construct such sequences with finite lengths. We recognize that such
constructed sequences (we call these "Erdos sequences") exhibit certain
hallmarks of randomness at the local level: they show roughly equal frequencies
of subsequences, and at the same time exclude the trivial periodic patterns.
For the human DNA we examine the frequency of a set of Erdos motifs of
length-10 using three nucleotides-to-binary mappings. The particular length-10
Erdos sequence is derived by the length-11 Mathias sequence and is identical
with the first 10 digits of the Thue-Morse sequence, underscoring the fact that
both are deficient in periodicities. Our calculations indicate that: (1) the
purine (A and G)/pyridimine (C and T) based Erdos motifs are greatly
underrepresented in the human genome, (2) the strong(G and C)/weak(A and T)
based Erdos motifs are slightly overrepresented, (3) the densities of the two
are negatively correlated, (4) the Erdos motifs based on all three mappings
being combined are slightly underrepresented, and (5) the strong/weak based
Erdos motifs are greatly overrepresented in the human messenger RNA sequences.


  The use of mathematical methods for the analysis of chemical reaction systems
has a very long history, and involves many types of models: deterministic
versus stochastic, continuous versus discrete, and homogeneous versus spatially
distributed. Here we focus on mathematical models based on deterministic
mass-action kinetics. These models are systems of coupled nonlinear
differential equations on the positive orthant. We explain how mathematical
properties of the solutions of mass-action systems are strongly related to key
properties of the networks of chemical reactions that generate them, such as
specific versions of reversibility and feedback interactions.


  The gravity model (GM) analogous to Newton's law of universal gravitation has
successfully described the flow between different spatial regions, such as
human migration, traffic flows, international economic trades, etc. This simple
but powerful approach relies only on the 'mass' factor represented by the scale
of the regions and the 'geometrical' factor represented by the geographical
distance. However, when the population has a subpopulation structure
distinguished by different attributes, the estimation of the flow solely from
the coarse-grained geographical factors in the GM causes the loss of
differential geographical information for each attribute. To exploit the full
information contained in the geographical information of subpopulation
structure, we generalize the GM for population flow by explicitly harnessing
the subpopulation properties characterized by both attributes and geography. As
a concrete example, we examine the marriage patterns between the bride and the
groom clans of Korea in the past. By exploiting more refined geographical and
clan information, our generalized GM properly describes the real data, a part
of which could not be explained by the conventional GM. Therefore, we would
like to emphasize the necessity of using our generalized version of the GM,
when the information on such nongeographical subpopulation structures is
available.


  A simple, but ``classical``, stochastic model for epidemic spread in a
finite, but large, population is studied. The progress of the epidemic can be
divided into three different phases that requires different tools to analyse.
Initially the process is approximated by a branching process. It is discussed
for how long time this approximation is valid. When a non-negligible proportion
of the population is already infected the process can be studied using
differential equations. In a final phase the spread will fade out. The results
are used to investigate what happens if two strains of infectious agents, with
different potential for spread, are simultaneously introduced in a totally
susceptible population. It is assumed that an infection causes immunity, and
that a person can only be infected by one strain. The two epidemics will
initially develop approximately as independent branching processes. However, if
both strains causes large epidemics they will, due to immunity, eventually
interact. We will mainly be interested in the final outcome of the spread,
i.e., how large proportion of the population is infected by the different
strains.


  How to define and use the concept of inclusive fitness is a contentious topic
in evolutionary theory. Inclusive fitness can be used to calculate selection on
a focal gene, but it is also applied to whole organisms. Individuals are then
predicted to appear designed as if to maximise their inclusive fitness,
provided that certain conditions are met (formally when interactions between
individuals are 'additive'). Here we argue that applying the concept of
inclusive fitness to organisms is justified under far broader conditions than
previously shown, but only if it is appropriately defined. Specifically, we
propose that organisms should maximise the sum of their offspring (including
any accrued due to the behaviour/phenotype of relatives), plus any effects on
their relatives' offspring production, weighted by relatedness. In contrast,
most theoreticians have argued that a focal individual's inclusive fitness
should exclude any offspring accrued due to the behaviour of relatives. Our
approach is based on the notion that long-term evolution follows the genome's
'majority interest' of building coherent bodies that are efficient 'vehicles'
for gene propagation. A gene favoured by selection that reduces the propagation
of unlinked genes at other loci (e.g. meiotic segregation distorters that lower
sperm production) is eventually neutralised by counter-selection throughout the
rest of the genome. Most phenotypes will therefore appear as if designed to
maximise the propagation of any given gene in a focal individual and its
relatives.


  Biological systems process information under noisy environment. Sensory
adaptation model of E. coli is suitable for investigation because of its
simplicity. To understand the adaptation processing quantitatively, stochastic
thermodynamic approach has been attempted. Information processing can be
assumed as state transition of a system that consists of signal transduction
molecules using thermodynamic approach, and efficiency can be measured as
thermodynamic cost. Recently, using information geometry and stochastic
thermodynamics, a relationship between speed of the transition and the
thermodynamic cost has been investigated for a chemical reaction model. Here,
we introduce this approach to sensory adaptation model of E. coli, and examined
a relationship between adaptation speed and the thermodynamic cost, and
efficiency of the adaptation speed. For increasing external noise level in
stimulation, the efficiency decreased, but the efficiency was highly robust to
external stimulation strength. Moreover, we demonstrated that there is the best
noise to achieve the adaptation in the aspect of thermodynamic efficiency. Our
quantification method provides a framework to understand the adaptation speed
and the thermodynamic cost for various biological systems.


  It is crucial to measure position and conformational changes of a
membrane-interacting protein relative to the membrane surface. This is however
challenging because the thickness of a membrane is usually only about 4 nm. We
developed a fluorescence method which makes use of the principle of FRET
between a fluorophore and a cloud of quenchers encapsulated in a liposome,
hence the name LipoFRET. LipoFRET can readily locate a fluorophore in different
depths inside and at different heights above the membrane. We applied LipoFRET
to study {\alpha}-synuclein, a key player in the pathology of Parkinson's
disease. Our approach yielded quantita-tive information about the dynamics of
different regions of {\alpha}-syn in lipid membranes, which has never been
explored before.


  Phylogenetic tree inference using deep DNA sequencing is reshaping our
understanding of rapidly evolving systems, such as the within-host battle
between viruses and the immune system. Densely sampled phylogenetic trees can
contain special features, including "sampled ancestors" in which we sequence a
genotype along with its direct descendants, and "polytomies" in which multiple
descendants arise simultaneously. These features are apparent after identifying
zero-length branches in the tree. However, current maximum-likelihood based
approaches are not capable of revealing such zero-length branches. In this
paper, we find these zero-length branches by introducing adaptive-LASSO-type
regularization estimators to phylogenetics, deriving their properties, and
showing regularization to be a practically useful approach for phylogenetics.


  Biological oscillators are very diverse but can be classified based on
dynamical motifs such as the types of feedback loops present. The S. Elongatus
circadian clock is a remarkable phosphorylation-based oscillator that can be
reconstituted in vitro with only 3 different purified proteins: the clock
proteins KaiA, KaiB, and KaiC. Despite a growing body of knowledge about the
biochemistry of the Kai proteins, basic questions about how their interactions
lead to sustained oscillations remain unanswered. Here, we compare models of
this system that make opposing assumptions about whether KaiA sequestration
introduces a positive or a negative feedback loop. We find that the two
different feedback mechanisms can be distinguished experimentally by the
introduction of a protein that binds competitively with KaiA. Understanding the
dynamical mechanism responsible for oscillations in the Kai system may shed
light on the broader question of what clock architectures have been selected by
evolution and why.


  Protein structure reconstruction from Nuclear Magnetic Resonance (NMR)
experiments largely relies on computational algorithms. Recently, some
effective low-rank matrix completion (MC) methods, such as ASD and ScaledASD,
have been successfully applied to image processing, which inspires us to apply
the methods to reconstruct protein structures. In this paper, we present an
efficient method to determine protein structures based on experimental NMR
NOESY distances. ScaledASD algorithm is used in the method with several
post-procedures including chirality refinement, distance lower (upper) bound
refinement, force field-based energy minimization (EM) and water refinement. By
comparing several metrics in the conformation evaluation on our results with
Protein Data Bank (PDB) structures, we conclude that our method is consistent
with the popularly used methods. In particular, our results show higher
validities in Procheck dihedral angles G-factor. Furthermore, we compare our
calculation results with PDB structures by examining the structural similarity
to X-ray crystallographic structures in a special dataset. The software and its
MATLAB source codes are available in https://github.com/xubiaopeng/PRASD


  In structured populations the spatial arrangement of cooperators and
defectors on the interaction graph together with the structure of the graph
itself determines the game dynamics and particularly whether or not fixation of
cooperation (or defection) is favored. For a single cooperator (and a single
defector) and a network described by a regular graph the question of fixation
can be addressed by a single parameter, the structure coefficient. As this
quantity is generic for any regular graph, we may call it the generic structure
coefficient. For two and more cooperators (or several defectors) fixation
properties can also be assigned by structure coefficients. These structure
coefficients, however, depend on the arrangement of cooperators and defectors
which we may interpret as a configuration of the game. Moreover, the
coefficients are specific for a given interaction network modeled as regular
graph, which is why we may call them specific structure coefficients. In this
paper, we study how specific structure coefficients vary over interaction
graphs and link the distributions obtained over different graphs to spectral
properties of interaction networks. We also discuss implications for the
benefit-to-cost ratios of donation games.


  The sociality of cattle facilitates the maintenance of herd cohesion and
synchronisation, making these species the ideal choice for domestication as
livestock for humans. However, livestock populations are not self-regulated,
and farmers transfer individuals across different groups throughout their lives
for reasons such as genetic mixing, reproduction and pastureland management.
Individuals consequently have to adapt to different group compositions during
their lives rather than choose their own herd mates, as they would do in the
wild. These changes may lead to social instability and stress, entailing
potentially negative effects on animal welfare. In this study, we assess how
the transfer of Highland cattle (Bos taurus) impacts individual and group
social network measures. Four groups with nine different compositions and 18
individual transfers were studied to evaluate 1) the effect of group
composition on individual social centralities and 2) the effect of group
composition changes on these centralities. As shown in previous works, dyadic
associations are stronger between individuals with similar age and dominance
rank. This study reveals that the relative stability of dyadic spatial
relationships between changes in group composition or enclosure is due to the
identities of transferred individuals more than the quantity of individuals
that are transferred. Older cattle had higher network centralities than other
individuals. The centrality of individuals was also affected by their sex and
the number of familiar individuals in the group. This study reveals the
necessity of understanding the social structure of a group to predict social
instability following the transfer of individuals between groups. The
developing of guidelines for the modification of group composition could
improve livestock management and reduce stress for the animals concerned.


  In this investigation we study the effects of climate change on predator-prey
population dynamics. The interaction of predator and prey is described by the
Variable Territory (VT) model with Allee effects in an Ordinary Differential
Equation (ODE) framework. Climate influence is modeled in this ODE framework as
an exogenous driver which affects the prey growth rate. Mathematically, this
exogenous driver is a function of time, the climate function, that describes
how favorable is the climate, on a scale from $-1$ to $+1$. Climate change has
a number of effects; in this work, we focus on changes in climate variability,
in particular, a decrease in the switching rate between good years and bad.
With regard to the predator-prey populations, we focus on the conditions that
lead to extinction. We developed a high performance software tool to allow the
systematic exploration of a wide range of climate scenarios at a high level of
detail. Our results show that the predator-prey system is sensitive to the
rapid drop in growth rate that accompanies a switch from good growth years to
bad growth years. The amplitude of this negative variation as well as its
duration are key factors, but more importantly, it is the moment when this
negative change occurs during the predator-prey cycle that is critical for the
survival of the species.


  Followership is generally defined as a strategy that evolved to solve social
coordination problems, and particularly those involved in group movement.
Followership behaviour is particularly interesting in the context of
road-crossing behaviour because it involves other principles such as
risk-taking and evaluating the value of social information. This study sought
to identify the cognitive mechanisms underlying decision-making by pedestrians
who follow another person across the road at the green or at the red light in
two different countries (France and Japan). We used agent-based modelling to
simulate the road-crossing behaviours of pedestrians. This study showed that
modelling is a reliable means to test different hypotheses and find the exact
processes underlying decision-making when crossing the road. We found that two
processes suffice to simulate pedestrian behaviours. Importantly, the study
revealed differences between the two nationalities and between sexes in the
decision to follow and cross at the green and at the red light. Japanese
pedestrians are particularly attentive to the number of already departed
pedestrians and the number of waiting pedestrians at the red light, whilst
their French counterparts only consider the number of pedestrians that have
already stepped off the kerb, thus showing the strong conformism of Japanese
people. Finally, the simulations are revealed to be similar to observations,
not only for the departure latencies but also for the number of crossing
pedestrians and the rates of illegal crossings. The conclusion suggests new
solutions for safety in transportation research.


  We study a stochastic model of biosynthesis of proteins in generic bacterial
operons. The stochasticity arises from two different processes, namely from
`bursting' production of either mRNA and/or protein (in the
transcription/translation process) and from standard diffusive fluctuations.
The amount of protein follows the Feller diffusion, while the bursting
introduces random jumps between trajectories of the diffusion process. The
combined effect leads to a process commonly known as a diffusion process with
jumps. We study existence of invariant densities and the long time behavior of
distributions of the corresponding Markov process, proving asymptotic stability
in the evolution of the density.


  We formulate a nonequilibrium thermodynamic description for open chemical
reaction networks (CRN) described by a chemical master equation. The
topological properties of the CRN and its conservation laws are shown to play a
crucial role. They are used to decompose the entropy production into a
potential change and two work contributions, the first due to time dependent
changes in the externally controlled chemostats concentrations and the second
due to flows maintained across the system by nonconservative forces. These two
works jointly satisfy a Jarzynski and Crooks fluctuation theorem. In absence of
work, the potential is minimized by the dynamics as the system relaxes to
equilibrium and its equilibrium value coincides with the maximum entropy
principle. A generalized Landauer's principle also holds: the minimal work
needed to create a nonequilibrium state is the relative entropy of that state
to its equilibrium value reached in absence of any work.


  Evolutionary modeling applications are the best way to provide full
information to support in-depth understanding of evaluation of organisms. These
applications mainly depend on identifying the evolutionary history of existing
organisms and understanding the relations between them, which is possible
through the deep analysis of their biological sequences. Multiple Sequence
Alignment (MSA) is considered an important tool in such applications, where it
gives an accurate representation of the relations between different biological
sequences. In literature, many efforts have been put into presenting a new MSA
algorithm or even improving existing ones. However, little efforts on
optimizing parallel MSA algorithms have been done. Nowadays, large datasets
become a reality, and big data become a primary challenge in various fields,
which should be also a new milestone for new bioinformatics algorithms. This
survey presents four of the state-of-the-art parallel MSA algorithms, TCoffee,
MAFFT, MSAProbs, and M2Align. We provide a detailed discussion of each
algorithm including its strengths, weaknesses, and implementation details and
the effectiveness of its parallel implementation compared to the other
algorithms, taking into account the MSA accuracy on two different datasets,
BAliBASE and OXBench.


  Scientists are attempting to use models of ever increasing complexity,
especially in medicine, where gene-based diseases such as cancer require better
modeling of cell regulation. Complex models suffer from uncertainty and
experiments are needed to reduce this uncertainty. Because experiments can be
costly and time-consuming it is desirable to determine experiments providing
the most useful information. If a sequence of experiments is to be performed,
experimental design is needed to determine the order. A classical approach is
to maximally reduce the overall uncertainty in the model, meaning maximal
entropy reduction. A recently proposed method takes into account both model
uncertainty and the translational objective, for instance, optimal structural
intervention in gene regulatory networks, where the aim is to alter the
regulatory logic to maximally reduce the long-run likelihood of being in a
cancerous state. The mean objective cost of uncertainty (MOCU) quantifies
uncertainty based on the degree to which model uncertainty affects the
objective. Experimental design involves choosing the experiment that yields the
greatest reduction in MOCU. This paper introduces finite-horizon dynamic
programming for MOCU-based sequential experimental design and compares it to
the greedy approach, which selects one experiment at a time without
consideration of the full horizon of experiments. A salient aspect of the paper
is that it demonstrates the advantage of MOCU-based design over the widely used
entropy-based design for both greedy and dynamic-programming strategies and
investigates the effect of model conditions on the comparative performances.


  Lattice-based stochastic simulators are commonly used to study biological
reaction-diffusion processes. Some of these schemes that are based on the
reaction-diffusion master equation (RDME), can simulate for extended spatial
and temporal scales but cannot directly account for the microscopic effects in
the cell such as volume exclusion and diffusion-influenced reactions.
Nonetheless, schemes based on the high-resolution microscopic lattice method
(MLM) can directly simulate these effects by representing each finite-sized
molecule explicitly as a random walker on fine lattice voxels. The theory and
consistency of MLM in simulating diffusion-influenced reactions have not been
clarified in detail. Here, we examine MLM in solving diffusion-influenced
reactions in 3D space by employing the Spatiocyte simulation scheme. Applying
the random walk theory, we construct the general theoretical framework
underlying the method and obtain analytical expressions for the total rebinding
probability and the effective reaction rate. By matching Collins-Kimball and
lattice-based rate constants, we obtained the exact expressions to determine
the reaction acceptance probability and voxel size. We found that the size of
voxel should be about 2% larger than the molecule. MLM is validated by
numerical simulations, showing good agreement with the off-lattice
particle-based method, eGFRD. MLM run time is more than an order of magnitude
faster than eGFRD when diffusing macromolecules with typical concentrations in
the cell. MLM also showed good agreements with eGFRD and mean-field models in
case studies of two basic motifs of intracellular signaling, the protein
production-degradation process and the dual phosphorylation cycle. Moreover,
when a reaction compartment is populated with volume-excluding obstacles, MLM
captures the non-classical reaction kinetics caused by anomalous diffusion of
reacting molecules.


  Phylogenetic trees are the fundamental mathematical representation of
evolutionary processes in biology. They are also objects of interest in pure
mathematics, such as algebraic geometry and combinatorics, due to their
discrete geometry. Although they are important data structures, they face the
significant challenge that sets of trees form a non-Euclidean phylogenetic tree
space, which means that standard computational and statistical methods cannot
be directly applied. In this work, we explore the statistical feasibility of a
pure mathematical representation of the set of all phylogenetic trees based on
tropical geometry for both descriptive and inferential statistics, and
unsupervised and supervised machine learning. Our exploration is both
theoretical and practical. We show that the tropical geometric phylogenetic
tree space endowed with a generalized Hilbert projective metric exhibits
analytic, geometric, and topological properties that are desirable for
theoretical studies in probability and statistics and allow for well-defined
questions to be posed. We illustrate the statistical feasibility of the
tropical geometric perspective for phylogenetic trees with an example of both a
descriptive and inferential statistical task. Moreover, this approach exhibits
increased computational efficiency and statistical performance over the current
state-of-the-art, which we illustrate with a real data example on seasonal
influenza. Our results demonstrate the viability of the tropical geometric
setting for parametric statistical and probabilistic studies of sets of
phylogenetic trees.


  A fundamental goal of microbial ecology is to understand what determines the
diversity, stability, and structure of microbial ecosystems. The microbial
context poses special conceptual challenges because of the strong mutual
influences between the microbes and their chemical environment through the
consumption and production of metabolites. By analyzing a generalized consumer
resource model that explicitly includes cross-feeding, stochastic colonization,
and thermodynamics, we show that complex microbial communities generically
exhibit a transition as a function of available energy fluxes from a
"resource-limited" regime where community structure and stability is shaped by
energetic and metabolic considerations to a diverse regime where the dominant
force shaping microbial communities is the overlap between species' consumption
preferences. These two regimes have distinct species abundance patterns,
different functional profiles, and respond differently to environmental
perturbations. Our model reproduces large-scale ecological patterns observed
across multiple experimental settings such as nestedness and differential beta
diversity patterns along energy gradients. We discuss the experimental
implications of our results and possible connections with disorder-induced
phase transitions in statistical physics.


  Scanning fluorescence correlation spectroscopy (SFCS) with a scan path
perpendicular to the membrane plane was introduced to measure diffusion and
interactions of fluorescent components in free standing biomembranes. Using a
confocal laser scanning microscope (CLSM) the open detection volume is moved
laterally with kHz frequency through the membrane and the photon events are
continuously recorded and stored in a file. While the accessory hardware
requirements for a conventional CLSM are minimal, data evaluation can pose a
bottleneck. The photon events must be assigned to each scan, in which the
maximum signal intensities have to be detected, binned, and aligned between the
scans, in order to derive the membrane related intensity fluctuations of one
spot. Finally, this time-dependent signal must be correlated and evaluated by
well known FCS model functions. Here we provide two platform independent, open
source software tools (PyScanFCS and PyCorrFit) that allow to perform all of
these steps and to establish perpendicular SFCS in its one- or two-focus as
well as its single- or dual-colour modality.


  A double-stranded DNA unravels thermally through intermediate denatured
bubble segments. Intrinsically, fluctuations ensue at the bubble boundaries
from non-equilibrium (NE) energy exchanges with the environment. However, such
details gets obscured by large population kinetics at the macroscale,
associating equilibrium pathway to the unravelling landscape. In this work, we
capture evidence of fluctuation energetics with picoliter samples in a
microfluidic cantilever. We exploit nanomechanical resonance to measure the NE
energy exchanges through dissipation, revealing a crucial pre-melting
transition at T~42C . This signifies that unravelling possibly proceeds via
intermediate collapsed-bubble conformations releasing energy, sufficient to
unbind bubble ends, assisting further unbinding. Fluctuation theorem explains
the observations opening further avenues to investigate fluctuation kinetics in
other biological phenomena that also proceed through similar NE energetics.


  Stochastic models that incorporate birth, death and immigration (also called
birth-death and innovation models) are ubiquitous and applicable to many
research topics such as quantifying species sizes in ecological populations,
describing gene family sizes, modeling lymphocyte evolution in the body, and
modeling the evolution of firm sizes. Many of these applications involve the
immigration of new species into the system. We develop the full
high-dimensional stochastic process associated with multispecies
birth-death-immigration processes and present a number of exact and asymptotic
results for the steady-state solutions to these types of processes. We
considered models that include carrying capacity and random mutations and find
analytic and approximate results for the statistics of the total number of
individuals, the total number of species, the species size distribution, and
various diversity indices. Our results include a rigorous analysis of the
behavior of these systems in the fast immigration limit which shows that of the
different diversity indices, the species richness is best able to distinguish
different types of birth-death-immigration models. Surprisingly, we also find
that detailed balance is preserved in all of our neutral
birth-death-immigration models, except for those that incorporate carrying
capacity through the birth rate.


  We construct a pathwise formulation for a multi-type age-structured
population dynamics, which involves an age-dependent cell replication and
transition of gene- or phenotypes. By employing the formulation, we derive a
variational representation of the stationary population growth rate; the
representation comprises a trade-off relation between growth effects and a
single-cell intrinsic dynamics described by a semi-Markov process. This
variational representation leads to a response relation of the stationary
population growth rate, in which statistics on a retrospective history work as
the response coefficients. These results contribute to predicting and
controlling growing populations based on experimentally observed cell-lineage
information.


  Phenotypic variability in a population of cells can work as the bet-hedging
of the cells under an unpredictably changing environment, the typical example
of which is the bacterial persistence. To understand the strategy to control
such phenomena, it is indispensable to identify the phenotype of each cell and
its inheritance. Although recent advancements in microfluidic technology offer
us useful lineage data, they are insufficient to directly identify the
phenotypes of the cells. An alternative approach is to infer the phenotype from
the lineage data by latent-variable estimation. To this end, however, we must
resolve the bias problem in the inference from lineage called survivorship
bias. In this work, we clarify how the survivor bias distorts statistical
estimations. We then propose a latent-variable estimation algorithm without the
survivorship bias from lineage trees based on an expectation-maximization (EM)
algorithm, which we call Lineage EM algorithm (LEM). LEM provides a statistical
method to identify the traits of the cells applicable to various kinds of
lineage data.


  A chemostat enables long-term, continuous, exponential-phase growth in an
environment limited as prescribed by the researcher. It is thus a potent tool
for laboratory evolution - selecting for strains with desired phenotypes.
However, despite the apparently simple design governed by a limited set of
rules, analysis of chemostat dynamics shows that they display counter-intuitive
properties. For example, the concentration of limiting substrate in the
chemostat is independent of the concentration in the influx and only dependent
on the dilution rate and the strain parameters. Moreover, choosing optimal
operational parameters (dilution rate, volume size, etc.) can be challenging.
There are conflicting requirements in the experimental design, such as a need
for relatively fast growth conditions for mutation accumulation on the one hand
versus slow dilution for a large fitness advantage for mutants to take over the
population quickly on the other.In this study, we provide analytic and
computational tools to help understand and predict chemostat dynamics, and
choose suitable operational parameters. We refer to five stages of the process:
(A) parameter choice and setup, (B) basic steady state growth, (C) mutation,
(D) single takeover and (E) successive takeovers. We present a qualitative and
quantitative framework to answer the questions confronted in each of these
stages. We provide a set of simulations which support the quantitative results,
and a graphical user interface to give a hands-on opportunity to experience and
visualize the analytic results. We detail conditions that produce ineffectual
selection regimes, and find that when avoided, the selection time is relatively
robust, and usually varies by less than an order of magnitude. Finally, we
suggest rules of thumb to help ensure that the chosen parameters lead to
effective selection and minimize the duration of the selection process.


  We demonstrate the phenomenon of cumulative inertia in intracellular
transport involving multiple motor proteins in human epithelial cells by
measuring the empirical survival probability of cargoes on the microtubule and
their detachment rates. We found the longer a cargo moves along a microtubule,
the less likely it detaches from it. As a result, the movement of cargoes is
non-Markovian and involves a memory. We observe memory effects on the scale of
up to 2 seconds. We provide a theoretical link between the measured detachment
rate and the super-diffusive Levy walk-like cargo movement.


  Group living animals form aggregations and flocks that remain cohesive in
spite of internal movements of individuals. This is possible because individual
group members repeatedly adjust their position and motion in response to the
position and motion of other group members. Here we develop a theoretical
approach to address the question, what general features -- if any -- underlie
the interaction rules that mediate group stability in animals of all species?
We do so by considering how the spatial organisation of a group would change in
the complete absence of interactions. Without interactions, a group would
disperse in a way that can be easily characterised in terms of Fick's diffusion
equations. We can hence address the inverse theoretical problem of finding the
individual-level interaction responses that are required to counterbalance
diffusion and to preserve group stability. We show that an individual-level
response to neighbour densities in the form of Weber's law (a 'universal' law
describing the functioning of the sensory systems of animals of all species)
results in an 'anti-diffusion' term at the group level. On short time scales,
this anti-diffusion restores the initial group configuration in a way which is
reminiscent of methods for image deblurring in image processing. We also show
that any non-homogeneous, spatial density distribution can be preserved over
time if individual movement patterns have the form of a Weber's law response.
Weber's law describes the fundamental functioning of perceptual systems. Our
study indicates that it is also a necessary -- but not sufficient -- feature of
collective interactions in stable animal groups.


  The low-frequency collective vibrational modes in proteins as well as the
protein-water interface have been suggested as dominant factors controlling the
efficiency of biochemical reactions and biological energy transport. It is thus
crucial to uncover the mystery of hydration structure and dynamics as well as
their coupling to collective motions of proteins in aqueous solutions. Here we
report dielectric properties of aqueous bovine serum albumin protein solutions
as a model system using an extremely sensitive dielectric spectrometer with
frequencies spanning from megahertz to terahertz. The dielectric relaxation
spectra reveal several polarization mechanisms at the molecular level with
different time constants and dielectric strengths, reflecting the complexity of
protein-water interactions. Combining the effective-medium approximation and
molecular dynamics simulations, we have determined collective vibrational modes
at terahertz frequencies and the number of water molecules in the tightly-bound
and loosely-bound hydration layers. High-precision measurements of the number
of hydration water molecules indicate that the dynamical influence of proteins
extends beyond the first solvation layer, to around 7 {\AA} distance from the
protein surface, with the largest slowdown arising from water molecules
directly hydrogen-bonded to the protein. Our results reveal critical
information of protein dynamics and protein-water interfaces, which determine
biochemical functions and reactivity of proteins.


  We study the coevolutionary dynamics of the diversity of phenotype expression
and the evolution of cooperation in the Prisoner's Dilemma game. Rather than
pre-assigning zero-or-one interaction rate, we diversify the rate of
interaction by associating it with the phenotypes shared in common. Individuals
each carry a set of potentially expressible phenotypes and expresses a certain
number of phenotypes at a cost proportional to the number. The number of
expressed phenotypes and thus the rate of interaction is an evolvable trait.
Our results show that nonnegligible cost of expressing phenotypes restrains
phenotype expression, and the evolutionary race mainly proceeds on between
cooperative strains and defective strains who express a very few phenotypes. It
pays for cooperative strains to express a very few phenotypes. Though such a
low level of expression weakens reciprocity between cooperative strains, it
decelerates rate of interaction between cooperative strains and defective
strains to a larger degree, leading to the predominance of cooperative strains
over defective strains. We also find that evolved diversity of phenotype
expression can occasionally destabilize due to the invasion of defective
mutants, implying that cooperation and diversity of phenotype expression can
mutually reinforce each other. Therefore, our results provide new insights into
better understanding the coevolution of cooperation and the diversity of
phenotype expression.


  Representation learning is at the heart of what makes deep learning
effective. In this work, we introduce a new framework for representation
learning that we call "Holographic Neural Architectures" (HNAs). In the same
way that an observer can experience the 3D structure of a holographed object by
looking at its hologram from several angles, HNAs derive Holographic
Representations from the training set. These representations can then be
explored by moving along a continuous bounded single dimension. We show that
HNAs can be used to make generative networks, state-of-the-art regression
models and that they are inherently highly resistant to noise. Finally, we
argue that because of their denoising abilities and their capacity to
generalize well from very few examples, models based upon HNAs are particularly
well suited for biological applications where training examples are rare or
noisy.


  Gene expression dynamics satisfying given input-output relationships were
investigated by evolving the networks for an optimal response. We found three
types of networks and corresponding dynamics, depending on the sensitivity of
gene expression dynamics: direct response with straight paths, amplified
response by a feed-forward network, and cooperative response with a complex
network. When the sensitivity of each gene's response is low and expression
dynamics is sloppy, the last type is selected, in which many genes respond
collectively to inputs, with local-excitation and global-inhibition structures.
The result provides an insight into how a reliable response is achieved with
unreliable units, and on why complex networks with many genes are adopted in
cells.


  The stability of a complex system generally decreases with increasing system
size and interconnectivity, a counterintuitive result of widespread importance
across the physical, life, and social sciences. Despite recent interest in the
relationship between system properties and stability, the effect of variation
in response rate across system components remains unconsidered. Here I vary the
component response rates ($\boldsymbol{\gamma}$) of randomly generated complex
systems. I use numerical simulations to show that when component response rates
vary, the potential for system stability increases. These results are robust to
common network structures, including small-world and scale-free networks, and
cascade food webs. Variation in $\boldsymbol{\gamma}$ is especially important
for stability in highly complex systems, in which the probability of stability
would otherwise be negligible. At such extremes of simulated system complexity,
the largest stable complex systems would be unstable if not for variation in
$\boldsymbol{\gamma}$. My results therefore reveal a previously unconsidered
aspect of system stability that is likely to be pervasive across all realistic
complex systems.


  Probability modelling for DNA sequence evolution is well established and
provides a rich framework for understanding genetic variation between samples
of individuals from one or more populations. We show that both classical and
more recent models for coalescence (with or without recombination) can be
described in terms of the so-called phase-type theory, where complicated and
tedious calculations are circumvented by the use of matrices. The application
of phase-type theory consists of describing the stochastic model as a Markov
model by appropriately setting up a state space and calculating the
corresponding intensity and reward matrices. Formulae of interest are then
expressed in terms of these aforementioned matrices. We illustrate this by a
few examples calculating the mean, variance and even higher order moments of
the site frequency spectrum in the multiple merger coalescent models, and by
analysing the mean and variance for the number of segregating sites for
multiple samples in the two-locus ancestral recombination graph. We believe
that phase-type theory has great potential as a tool for analysing probability
models in population genetics. The compact matrix notation is useful for
clarification of current models, in particular their formal manipulation
(calculation), but also for further development or extensions.


  Eicosanoids and related species are critical, small bioactive mediators of
human physiology and inflammation. While ~1100 distinct eicosanoids have been
predicted to exist, to date, less than 150 of these molecules have been
measured in humans, limiting our understanding of eicosanoids and their role in
human biology. Using a directed non-targeted mass spectrometry approach in
conjunction with computational chemical networking of spectral fragmentation
patterns, we find over 500 discrete chemical signals highly consistent with
known and putative eicosanoids in human plasma, including 46 putative novel
molecules not previously described, thereby greatly expanding the breath of
prior analytical strategies. In plasma samples from 1500 individuals, we find
members of this expanded eicosanoid library hold close association with markers
of inflammation, as well as clinical characteristics linked with inflammation,
including advancing age and obesity. These experimental and computational
approaches enable discovery of new chemical entities and will shed important
insight into the role of bioactive molecules in human disease.


  It has been shown recently that changing the fluidic properties of a drug can
improve its efficacy in ablating solid tumors. We develop a modeling framework
for tumor ablation, and present the simplest possible model for drug diffusion
in a spherical tumor with leaky boundaries and assuming cell death eventually
leads to ablation of that cell effectively making the two quantities
numerically equivalent. The death of a cell after a given exposure time depends
on both the concentration of the drug and the amount of oxygen available to the
cell. Higher oxygen availability leads to cell death at lower drug
concentrations. It can be assumed that a minimum concentration is required for
a cell to die, effectively connecting diffusion with efficacy. The
concentration threshold decreases as exposure time increases, which allows us
to compute dose-response curves. Furthermore, these curves can be plotted at
much finer time intervals compared to that of experiments, which is used to
produce a dose-threshold-response surface giving an observer a complete picture
of the drug's efficacy for an individual. In addition, since the diffusion,
leak coefficients, and the availability of oxygen is different for different
individuals and tumors, we produce artificial replication data through
bootstrapping to simulate error. While the usual data-driven model with
Sigmoidal curves use 12 free parameters, our mechanistic model only has two
free parameters, allowing it to be open to scrutiny rather than forcing
agreement with data. Even so, the simplest model in our framework, derived
here, shows close agreement with the bootstrapped curves, and reproduces well
established relations, such as Haber's rule.


  Urbanisation introduces deep changes in habitats, eventually creating new
urban ecosystems where ecological functions are driven by human activities. The
higher frequency of some phenotypes in urban vs rural/wild areas has led to the
assumption that directional selection in urban habitats occurs, which may
thereby favour some behavioural and physiological traits in urban animal
populations compared to rural ones. However, empirical evidence of directional
selection on phenotypic traits in urban areas remains scarce. In this study we
tested whether eumelanin-based colouration could be linked to survival in two
urban populations of the feral pigeon Columba livia. A number of studies in
different cities pointed out a higher frequency of darker individuals in more
urbanised areas compared to rural ones. To investigate whether directional
selection through survival on this highly heritable trait could explain such
patterns, we conducted mark-recapture studies on two populations of feral
pigeons in highly urbanized areas. We predicted that darker coloured
individuals would exhibit higher survival and/or philopatry (integrated into
'local survival') than paler coloured ones. No difference in local survival was
found between adults of different colouration intensities. However, on one
site, we found that darker juveniles had a higher local survival probability
than light ones. Juvenile local survival on that site was also negatively
correlated with the number of chicks born. This suggests the existence of
colour- and/or density-dependent selection processes acting on juvenile feral
pigeons in urban environments, acting through differential mortality and/or
dispersal.


  The backbone of most proteins forms an open curve. To study their
entanglement, a common strategy consists in searching for the presence of knots
in their backbones using topological invariants. However, this approach
requires to close the curve into a loop, which alters the geometry of curve.
Knoto-ID allows evaluating the entanglement of open curves without the need to
close them, using the recent concept of knotoids which is a generalization of
the classical knot theory to open curves. Knoto-ID can analyse the global
topology of the full chain as well as the local topology by exhaustively
studying all subchains or only determining the knotted core. Knoto-ID permits
to localize topologically non-trivial protein folds that are not detected by
informatics tools detecting knotted protein folds.


  Adoptive Cell Transfer therapy of cancer is currently in full development and
mathematical modeling is playing a critical role in this area. We study a
stochastic model developed by Baar et al. in 2015 for modeling immunotherapy
against melanoma skin cancer. First, we estimate the parameters of the
deterministic limit of the model based on biological data of tumor growth in
mice. A Nonlinear Mixed Effects Model is estimated by the Stochastic
Approximation Expectation Maximization algorithm. With the estimated
parameters, we head back to the stochastic model and calculate the probability
that the T cells all get exhausted during the treatment. We show that for some
relevant parameter values, an early relapse is due to stochastic fluctuations
(complete T cells exhaustion) with a non negligible probability. Then, focusing
on the relapse related to the T cell exhaustion, we propose to optimize the
treatment plan (treatment doses and restimulation times) by minimizing the T
cell exhaustion probability in the parameter estimation ranges.


  Recent advances in molecular simulations allow the evaluation of previously
unattainable observables, such as rate constants for protein folding. However,
these calculations are usually computationally expensive and even significant
computing resources may result in a small number of independent estimates
spread over many orders of magnitude. Such small-sample, high "log-variance"
data are not readily amenable to analysis using the standard uncertainty (i.e.,
"standard error of the mean") because unphysical negative limits of confidence
intervals result. Bootstrapping, a natural alternative guaranteed to yield a
confidence interval within the minimum and maximum values, also exhibits a
striking systematic bias of the lower confidence limit in log space. As we
show, bootstrapping artifactually assigns high probability to improbably low
mean values. A second alternative, the Bayesian bootstrap strategy, does not
suffer from the same deficit and is more logically consistent with the type of
confidence interval desired. The Bayesian bootstrap provides uncertainty
intervals that are more reliable than those from the standard bootstrap method,
but must be used with caution nevertheless. Neither standard nor Bayesian
bootstrapping can overcome the intrinsic challenge of under-estimating the mean
from small-size, high log-variance samples. Our conclusions are based on
extensive analysis of model distributions and re-analysis of multiple
independent atomistic simulations. Although we only analyze rate constants,
similar considerations will apply to related calculations, potentially
including highly non-linear averages like the Jarzynski relation.


  In this letter, we analyze a genetic toggle switch recently studied in the
literature where the expression of two repressor proteins can be tuned by
controlling two different inputs, namely the concentration of two inducer
molecules in the growth medium of the cells. Specifically, we investigate the
dynamics of this system when subject to pulse-width modulated (PWM) input. We
provide an analytical model that captures qualitatively the experimental
observations reported in the literature and approximates its asymptotic
behavior. We also discuss the effect that the system parameters have on the
prediction accuracy of the model. Moreover, we propose a possible external
control strategy to regulate the mean value of the fluorescence of the reporter
proteins when the cells are subject to such periodic forcing.


  The origins of new genes are among the most fundamental questions in
evolutionary biology. Our understanding of the ways that new genetic material
appears and how that genetic material shapes population variation remains
incomplete. De novo genes and duplicate genes are a key source of new genetic
material on which selection acts. To better understand the origins of these new
gene sequences, we explored the ways that structural variation might alter
expression patterns and form novel transcripts. We provide evidence that
chromosomal rearrangements are a source of novel genetic variation that
facilitates the formation of de novo genes in Drosophila. We identify 52 cases
of de novo gene formation created by chromosomal rearrangements in 14 strains
of D. yakuba. These new genes inherit transcription start signals and open
reading frames when the 5' end of existing genes are combined with previously
untranscribed regions. Such new genes would appear with novel peptide
sequences, without the necessity for secondary transitions from non-coding RNA
to protein. This mechanism of new peptide formations contrasts with canonical
theory of de novo gene progression requiring non-coding intermediaries that
must acquire new mutations prior to loss via pseudogenization. Hence, these
mutations offer a means to de novo gene creation and protein sequence formation
in a single mutational step, answering a long standing open question concerning
new gene formation. We further identify gene expression changes to 134 existing
genes, indicating that these mutations can alter gene regulation. Population
variability for chromosomal rearrangements is considerable, with 2368
rearrangements observed across 14 inbred lines. More rearrangements were
identified on the X chromosome than any of the autosomes, suggesting the X is
more susceptible to chromosome alterations.


  Mitochondrial oxidative phosphorylation (mOxPhos) makes ATP, the energy
currency of life. Chemiosmosis, a proton centric mechanism, advocates that
Complex V harnesses a transmembrane potential (TMP) for ATP synthesis. This
perception of cellular respiration requires oxygen to stay tethered at Complex
IV (an association inhibited by cyanide) and diffusible reactive oxygen species
(DROS) are considered wasteful and toxic products. With new mechanistic
insights on heme and flavin enzymes, an oxygen or DROS centric explanation
(called murburn concept) was recently proposed for mOxPhos. In the new
mechanism, TMP is not directly harnessed, protons are a rate limiting reactant
and DROS within matrix serve as the chemical coupling agents that directly link
NADH oxidation with ATP synthesis. Herein, we report multiple ADP binding sites
and solvent accessible DROS channels in respiratory proteins, which validate
the oxygen or DROS centric power generation (ATP synthesis) system in mOxPhos.
Since cyanide's heme binding Kd is high (mM), low doses (uM) of cyanide is
lethal because cyanide disrupts DROS dynamics in mOxPhos. The critical study
also provides comprehensive arguments against Mitchell's and Boyer's
explanations and extensive support for murburn concept based holistic
perspectives for mOxPhos.