  Bridging the exponentially growing gap between the numbers of unlabeled and
labeled protein sequences, several studies adopted semi-supervised learning for
protein sequence modeling. In these studies, models were pre-trained with a
substantial amount of unlabeled data, and the representations were transferred
to various downstream tasks. Most pre-training methods solely rely on language
modeling and often exhibit limited performance. In this paper, we introduce a
novel pre-training scheme called PLUS, which stands for Protein sequence
representations Learned Using Structural information. PLUS consists of masked
language modeling and a complementary protein-specific pre-training task,
namely same-family prediction. PLUS can be used to pre-train various model
architectures. In this work, we use PLUS to pre-train a bidirectional recurrent
neural network and refer to the resulting model as PLUS-RNN. Our experiment
results demonstrate that PLUS-RNN outperforms other models of similar size
solely pre-trained with the language modeling in six out of seven widely used
protein biology tasks. Furthermore, we present the results from our qualitative
interpretation analyses to illustrate the strengths of PLUS-RNN. PLUS provides
a novel way to exploit evolutionary relationships among unlabeled proteins and
is broadly applicable across a variety of protein biology tasks. We expect that
the gap between the numbers of unlabeled and labeled proteins will continue to
grow exponentially, and the proposed pre-training method will play a larger
role.


  In the last decade, advances in molecular dynamics (MD) and Markov State
Model (MSM) methodologies have made possible accurate and efficient estimation
of kinetic rates and reactive pathways for complex biomolecular dynamics
occurring on slow timescales. A promising approach to enhanced sampling of MSMs
is to use so-called "adaptive" methods, in which new MD trajectories are
"seeded" preferentially from previously identified states. Here, we investigate
the performance of various MSM estimators applied to reseeding trajectory data,
for both a simple 1D free energy landscape, and for mini-protein folding MSMs
of WW domain and NTL9(1-39). Our results reveal the practical challenges of
reseeding simulations, and suggest a simple way to reweight seeding trajectory
data to better estimate both thermodynamic and kinetic quantities.


  An environment far from equilibrium is thought to be a necessary condition
for the origin and persistence of life. In this context we report open-flow
simulations of a non-enzymic proto-metabolic system, in which hydrogen peroxide
acts both as oxidant and driver of thermochemical cycling. We find that a
Gaussian perturbed input produces a non-Boltzmann output fluctuation
distribution around the mean oscillation maximum. Our main result is that net
biosynthesis can occur under fluctuating cyclical but not steady drive.
Consequently we may revise the necessary condition to "dynamically far from
equilibrium".


  Growing efforts to measure fitness landscapes in molecular and microbial
systems are premised on a tight relationship between landscape topography and
evolutionary trajectories. This relationship, however, is far from being
straightforward: depending on their mutation rate, Darwinian populations can
climb the closest fitness peak (survival of the fittest), settle in lower
regions with higher mutational robustness (survival of the flattest), or fail
to adapt altogether (error catastrophes). These bifurcations highlight that
evolution does not necessarily drive populations "from lower peak to higher
peak", as Wright imagined. The problem therefore remains: how exactly does a
complex landscape topography constrain evolution, and can we predict where it
will go next? Here I introduce a generalization of quasispecies theory which
identifies metastable evolutionary states as minima of an effective potential.
From this representation I derive a coarse-grained, Markov state model of
evolution, which in turn forms a basis for evolutionary predictions across a
wide range of mutation rates. Because the effective potential is related to the
ground state of a quantum Hamiltonian, my approach could stimulate fruitful
interactions between evolutionary dynamics and quantum many-body theory


  We used machine learning methods to predict NaV1.7 inhibitors and found the
model RF-CDK that performed best on the imbalanced dataset. Using the RF-CDK
model for screening drugs, we got effective compounds K1. We use the cell patch
clamp method to verify K1. However, because the model evaluation method in this
article is not comprehensive enough, there is still a lot of research work to
be performed, such as comparison with other existing methods.
  The target protein has multiple active sites and requires our further
research. We need more detailed models to consider this biological process and
compare it with the current results, which is an error in this article.
  So we want to withdraw this article.


  The remarkable accuracy and versatility of single-molecule techniques make
possible new measurements that are not feasible in bulk assays. Among these,
the precise estimation of folding free energies using fluctuation theorems in
nonequilibrium pulling experiments has become a benchmark in modern biophysics.
In practice, the use of fluctuation relations to determine free energies
requires a thorough evaluation of the usually large energetic contributions
caused by the elastic deformation of the different elements of the experimental
setup (such as the optical trap, the molecular linkers and the
stretched-unfolded polymer). We review and describe how to optimally estimate
such elastic energy contributions to extract folding free energies, using DNA
and RNA hairpins as model systems pulled by laser optical tweezers. The
methodology is generally applicable to other force-spectroscopy techniques and
molecular systems.


  Assessing the risk of epidemic spread on networks and developing strategies
for its containment is of tremendous practical importance, both due to direct
effects in public health and its impact on economies. In this work we present
the numerical results of a stochastic, event-driven, hierarchical agent-based
model designed to reproduce the infectious dynamics of the cattle disease
called Bovine Viral Diarrhea (BVD), for which the corresponding network of
movements is the main route of spreading. For the farm-node dynamics, the model
takes into account a vast number of breeding, infectious and animal movement
mechanisms via a susceptible-infected-recovered (SIR) type of dynamics with an
additional permanently infectious class. The interaction between the farms is
described by a supply and demand farm manager mechanism governing the network
structure and dynamics. We discuss the disease and breeding dynamics, study
numerous mitigation strategies of present and past government regulations
taking Germany as a case-study country and perform a sensitivity analysis on
key parameters. We argue that the model, subject to calibration, has universal
predictive potential, can be extended to diseases beyond BVD and demonstrate
that appropriate measures can indeed lead to eradication regarding BVD. We
further present the results of the time-series of the model and conduct a
statistical analysis of and among the different mitigation strategies.


  Transport of ions and small molecules across the cell membrane against
electrochemical gradients is catalyzed by integral membrane proteins that use a
source of free energy to drive the energetically uphill flux of the transported
substrate. Secondary active transporters couple the spontaneous influx of a
"driving" ion such as Na$^+$ or H$^+$ to the flux of the substrate. The
thermodynamics of such cyclical non-equilibrium systems are well understood and
recent work has focused on the molecular mechanism of secondary active
transport. The fact that these transporters change their conformation between
an inward-facing and outward-facing conformation in a cyclical fashion, called
the alternating access model, is broadly recognized as the molecular framework
in which to describe transporter function. High resolution structures and
detailed computer simulations lead to the recognition of common molecular-level
principles between disparate transporter families. Inverted repeat symmetry in
secondary active transporters has shed light on how protein structures can
encode a bi-stable two-state system. Three broad classes of alternating access
transitions have been described as rocker-switch, rocking-bundle, and elevator
mechanisms. Transporters can be understood as gated pores with at least two
coupled gates that map to distinct parts of the transporter protein.
Enumerating all distinct gate states naturally includes occluded states in the
alternating access picture and suggests observable protein conformations. By
connecting the possible conformational states and ion/substrate bound states in
a kinetic model, a unified picture emerges in which symporter, antiporter, and
uniporter function are extremes in a continuum of functionality. We briefly
discuss how biological complexity may be integrated in quantitative kinetic
models to provide a bridge from structure to function.


  Predicting protein structure from the amino acid sequence has been a
challenge with theoretical and practical significance in biophysics. Despite
the recent progresses elicited by improved residue-residue contact prediction,
contact-based structure prediction has gradually reached the performance
ceiling. New methods have been proposed to predict the residue-residue
distance, but unanimously by simplifying the real-valued distance prediction
into a multiclass classification problem. Here we show a regression-based
distance prediction method, which adopts the generative adversarial network to
capture the delicate geometric relationship between residue pairs and thus
could predict the continuous, real-valued residue-residue distance
satisfactorily. The predicted residue distance map allows rapid structure
modeling by the CNS suite, and the constructed models approach at least the
same level of quality as the other state-of-the-art protein structure
prediction methods when tested on available CASP13 targets. Moreover, this
method can be used directly for the structure prediction of membrane proteins
without transfer learning.


  Temporal environmental variations are ubiquitous in nature, yet most of the
theoretical works in population genetics and evolution assume fixed
environment. Here we analyze the effect of variations in carrying capacity on
the fate of a mutant type. We consider a two-state Moran model, where selection
intensity at equilibrium may differ (in amplitude and in sign) from selection
during periods of sharp growth and sharp decline. Using Kimura's diffusion
approximation we present simple formulae for effective population size and
effective selection, and use it to calculate the chance of ultimate fixation,
the time to fixation and the time to absorption (either fixation or loss). Our
analysis shows perfect agreement with numerical solutions for neutral,
beneficial and deleterious mutant. The contributions of different processes to
the mean and the variance of abundance variations are additive and commutative.
As a result, when selection intensity $s$ is weak such that ${\cal O}(s^2)$
terms are negligible, periodic or stochastic environmental variations yield
identical results.


  Dynamic instability in the cytoskeleton underlies a great many vital cellular
processes. In neurons, it is thought to play a role in information processing
and the establishment of synaptic connections. It is from this interneuronal
connectivity that the large-scale organization and function of the brain
ultimately emerges. While brain networks above the cellular scale have been
extensively investigated with medical imaging, no research has yet examined the
contribution of the dynamic architecture subserving the components of these
networks to their functionality. Simultaneous functional imaging at molecular
and whole-brain scales may yield more comprehensive knowledge of this
entailment. Procedures of statistical analysis based on a recently developed
nonlinear model of microtubule dynamics are proposed to verify predictions
about the spatiotemporal relationships between various microtubule-associated
protein activities and functional connectivity. Known effects of brain cancer
on functional connectivity and peritumoral complexity serve as the basis for a
comparative study of these phenomena in regions of interest with particularly
high and low complexity.


  Cells are able to generate phenotypic diversity both during development and
in response to stressful and changing environments, aiding survival. The
biologically and medically vital process of a cell assuming a functionally
important fate from a range of phenotypic possibilities can be thought of as a
cell decision. To make these decisions, a cell relies on energy dependent
pathways of signalling and expression. However, energy availability is often
overlooked as a modulator of cellular decision-making. As cells can vary
dramatically in energy availability, this limits our knowledge of how this key
biological axis affects cell behaviour. Here, we consider the energy dependence
of a highly generalisable decision-making regulatory network, and show that
energy variability changes the sets of decisions a cell can make and the ease
with which they can be made. Increasing intracellular energy levels can
increase the number of stable phenotypes it can generate, corresponding to
increased decision-making capacity. For this decision-making architecture, a
cell with intracellular energy below a threshold is limited to a singular
phenotype, potentially forcing the adoption of a specific cell fate. We suggest
that common energetic differences between cells may explain some of the
observed variability in cellular decision-making, and demonstrate the
importance of considering energy levels in several diverse biological
decision-making phenomena.


  The complexity of human cancer often results in significant heterogeneity in
response to treatment. Precision medicine offers potential to improve patient
outcomes by leveraging this heterogeneity. Individualized treatment rules
(ITRs) formalize precision medicine as maps from the patient covariate space
into the space of allowable treatments. The optimal ITR is that which maximizes
the mean of a clinical outcome in a population of interest. Patient-derived
xenograft (PDX) studies permit the evaluation of multiple treatments within a
single tumor and thus are ideally suited for estimating optimal ITRs. PDX data
are characterized by correlated outcomes, a high-dimensional feature space, and
a large number of treatments. Existing methods for estimating optimal ITRs do
not take advantage of the unique structure of PDX data or handle the associated
challenges well. In this paper, we explore machine learning methods for
estimating optimal ITRs from PDX data. We analyze data from a large PDX study
to identify biomarkers that are informative for developing personalized
treatment recommendations in multiple cancers. We estimate optimal ITRs using
regression-based approaches such as Q-learning and direct search methods such
as outcome weighted learning. Finally, we implement a superlearner approach to
combine a set of estimated ITRs and show that the resulting ITR performs better
than any of the input ITRs, mitigating uncertainty regarding user choice of any
particular ITR estimation methodology. Our results indicate that PDX data are a
valuable resource for developing individualized treatment strategies in
oncology.


  Often in bacterial regulatory networks, small noncoding RNAs (sRNA) interact
with several mRNA species. The competition among mRNAs for binding to the same
pool of sRNA might lead to crosstalk between the mRNAs. This is similar to the
competing endogenous RNA (ceRNA) effect wherein the competition to bind to the
same pool of miRNA in Eukaryotes leads to miRNA mediated crosstalk resulting in
subtle and complex gene regulation with stabilised gene expression. We study an
sRNA-driven feedforward loop (sFFL) where the top-tier regulator, an sRNA
(RprA), translationally activates the target protein (RicI) directly and also,
indirectly, via up-regulation of its transcriptional activator (RpoS/sigma^s).
We show that the sRNA-mediated crosstalk between the two mRNA species leads to
maximum target protein synthesis for low synthesis rates of RpoS-mRNA. This
indicates the possibility of an optimal target protein synthesis with efficient
utilisation of RpoS-mRNA which is typically associated with various other
stress response activities inside the cell. Since gene expression is inherently
stochastic due to the probabilistic nature of various molecular interactions
associated with it, we next quantify the fluctuations in the target protein
level using generating function-based approach and stochastic simulations. The
coefficient of variation that provides a measure of fluctuations in the
concentration shows a minimum under conditions that also correspond to optimal
target protein synthesis. This prompts us to conclude that, in sFFL, the
crosstalk leads to optimal target protein synthesis with minimal noise and
efficient utilisation of RpoS-mRNA.


  We analyze the Conradi-Kahle Algorithm for detecting binomiality. We present
experiments using two implementations of the algorithm in Macaulay2 and Maple
on biological models and assess the performance of the algorithm on these
models. We compare the two implementations with each other and with Gr\"obner
bases computations up to their performance on these biological models.


  In arXiv:1911.06042v1 M. Casiulis et al. study a Hamiltonian model in which
rigid rotations of moving clusters give rise to scale-free correlations of
velocity and speed. M. Casiulis et al. compare correlations in their model to
those observed in real flocks of birds and claim that rigid-body rotations
provide an explanation that stands in contrast with, and it is simpler than,
previously proposed explanations of correlations in bird flocks, namely
Goldstone modes in the velocity orientations and marginal (or near-critical)
modes in the speed. Here, we show that the rigid rotation scenario is
completely inconsistent with a large body of well-established experimental
evidence on real flocks of birds, and it therefore does not provide an
appropriate explanation for the observed phenomenology.


  In ecological systems heterogeneous interactions between pathogens take place
simultaneously. This occurs, for instance, when two pathogens cooperate, while
at the same time multiple strains of these pathogens co-circulate and compete.
Notable examples include the cooperation of HIV with antibiotic-resistant and
susceptible strains of tuberculosis, or some respiratory infections with
Streptococcus pneumoniae strains. Models focusing on competition or cooperation
separately fail to describe how these concurrent interactions shape the
epidemiology of such diseases. We studied this problem considering two
cooperating pathogens, where one pathogen is further structured in two strains.
The spreading follows a susceptible-infected-susceptible process and the
strains differ in transmissibility and extent of cooperation with the other
pathogen. We combined a mean-field stability analysis with stochastic
simulations on networks considering both well-mixed and structured populations.
We observed the emergence of a complex phase diagram, where the conditions for
the less transmissible, but more cooperative strain to dominate are
non-trivial, e.g. non-monotonic boundaries and bistability. Coupled with
community structure, the presence of the cooperative pathogen enables the
co-existence between strains by breaking the spatial symmetry and dynamically
creating different ecological niches. These results shed light on ecological
mechanisms that may impact the epidemiology of diseases of public health
concern.


  Phosphate groups are naturally present in starch at C3- or C6-position of the
glucose residues and impact the structure of starch granules. Their precise
quantification is necessary for understanding starch physicochemical properties
and metabolism. Nevertheless, reliable quantification of Glc-3-P remains
laborious and time consuming. Here we describe a capillary electrophoresis
method for simultaneous measurement of both Glc-6-P and Glc-3-P after acid
hydrolysis of starch. The sensitivity threshold was estimated at the fg scale,
which is compatible with the analysis of less than a $\mu$g of sample. The
method was validated by analyzing antisense potato lines deficient in SBEs, GWD
or GBSS. We show that Glc-3-P content is altered in the latter and that these
variations do not correlate with modifications in Glc-6-P content. We
anticipate the method reported here to be an efficient tool for high throughput
study of starch phosphorylation at both C3- and C6-position.


  The glutamate metabotropic receptor 1 (GRM1) drives oncogenesis when
aberrantly activated in melanoma and several other cancers. Metabolomics
reveals that patient-derived xenografts with GRM1-positive melanoma tumors
exhibit elevated plasma glutamate levels associated with metastatic melanoma in
vivo. Stable isotope tracing and GCMS analysis determined that cells expressing
GRM1 fuel a substantial fraction of glutamate from glycolytic carbon.
Stimulation of GRM1 by glutamate leads to activation of mitogenic signaling
pathways, which in turn increases the production of glutamate, fueling
autocrine feedback. Implementing a rational drug-targeting strategy, we
critically evaluate metabolic bottlenecks in vitro and in vivo. Combined
inhibition of glutamate secretion and biosynthesis is an effective rational
drug targeting strategy suppressing tumor growth and restricting tumor
bioavailability of glutamate.


  A fitness landscape is a mapping from a space of discrete genotypes to the
real numbers. A path in a fitness landscape is a sequence of genotypes
connected by single mutational steps. Such a path is said to be accessible if
the fitness values of the genotypes encountered along the path increase
monotonically. We study accessible paths on random fitness landscapes of the
House-of-Cards type, on which fitness values are independent, identically and
continuously distributed random variables. The genotype space is taken to be a
Cartesian power graph $\mathcal{A}^L$, where $L$ is the number of genetic loci
and the allele graph $\mathcal{A}$ encodes the possible allelic states and
mutational transitions on one locus. The probability of existence of accessible
paths between two genotypes at a distance linear in $L$ displays a transition
from 0 to a positive value at a threshold $\beta_c$ for the fitness difference
between the initial and final genotype. We derive a lower bound on $\beta_c$
for general $\mathcal{A}$ and show that this bound is tight for a large class
of allele graphs. Our results generalize previous results for accessibility
percolation on the biallelic hypercube, and compare favorably to published
numerical results for multiallelic Hamming graphs.


  Background: Chronic myeloid leukemia (CML) results from hematopoietic stem
cell transformation by the bcr-abl chimeric oncogene, encoding a 210 kDa
protein with constitutive tyrosine kinase activity. In spite of the efficiency
of tyrosine kinase inhibitors (TKI; Imatinib), other strategies are explored to
eliminate CML leukemia stem cells, such as calcium pathways. Results: In this
work, we showed that Store-Operated Calcium Entry (SOCE) and thrombin induced
calcium influx were decreased in Bcr-Abl expressing 32d cells (32d-p210). The
32d-p210 cells showed modified Orai1/STIM1 ratio and reduced TRPC1 expression
that could explain SOCE reduction. Decrease in SOCE and thrombin induced
calcium entry was associated to reduced Nuclear Factor of Activated T cells
(NFAT) nucleus translocation in 32d-p210 cells. We demonstrated that SOCE
blockers enhanced cell mobility of 32d-p210 cells and reduced the proliferation
rate in both 32d cell lines. TKI treatment slightly reduced the
thrombin-induced response, but imatinib restored SOCE to the wild type level.
Bcr-Abl is also known to deregulate Protein Kinase C (PKC), which was described
to modulate calcium entries. We showed that PKC enhances SOCE and thrombin
induced calcium entries in control cells while this effect is lost in
Bcr-Abl-expressing cells. Conclusion: The tyrosine kinase activity seems to
regulate calcium entries probably not directly but through a global cellular
reorganization involving a PKC pathway. Altogether, calcium entries are
deregulated in Bcr-Abl-expressing cells and could represent an interesting
therapeutic target in combination with TKI.


  We review recent progress in the understanding of the role of multiple- and
simultaneous multiple merger coalescents as models for the genealogy in
idealised and real populations with exceptional reproductive behaviour. In
particular, we discuss models with `skewed offspring distribution' (or under
other non-classical evolutionary forces) which lead in the single locus haploid
case to multiple merger coalescents, and in the multi-locus diploid case to
simultaneous multiple merger coalescents. Further, we discuss inference methods
under the infinitely-many sites model which allow both model selection and
estimation of model parameters under these coalescents.


  Biopolymer self-assembly pathways are central to biological activity, but are
complicated by the ability of the monomeric subunits of biopolymers to adopt
different conformational states. As a result, biopolymer nucleation often
involves a two-step mechanism where the monomers first condense to form a
metastable intermediate, and this then converts to a stable polymer by
conformational rearrangement of its constituent monomers. While existing
mathematical models neglect the dynamics by which intermediates convert to
stable polymers, experiments and simulations show that these dynamics
frequently occur on comparable timescales to condensation of intermediates and
growth of mature polymers, and thus cannot be ignored. Moreover, nucleation
intermediates are responsible for cell toxicity in pathologies such as
Alzheimer's, Parkinson's, and prion diseases. Due to the relationship between
conformation and biological function, the slow conversion dynamics of these
species will strongly affect their toxicity. In this study, we present a
modified Oosawa model which explicitly accounts for simultaneous assembly and
conversion. To describe the conversion dynamics, we propose an experimentally
motivated initiation-propagation (IP) mechanism in which the stable phase
arises locally within the intermediate, and then spreads through additional
conversion events induced by nearest-neighbor interactions, analogous to
one-dimensional Glauber dynamics. Our mathematical analysis shows that the
competing timescales of assembly and conversion result in a nonequilibrium
critical point, separating a regime where intermediates are kinetically
unstable from one where conformationally mixed intermediates can accumulate.
Our work provides the first general model of two-step biopolymer nucleation,
which can be used to quantitatively predict the concentration and composition
of biologically crucial intermediates.


  A challenge in developing machine learning regression models is that it is
difficult to know whether maximal performance has been reached on a particular
dataset, or whether further model improvement is possible. In biology this
problem is particularly pronounced as sample labels (response variables) are
typically obtained through experiments and therefore have experiment noise
associated with them. Such label noise puts a fundamental limit to the
performance attainable by regression models. We address this challenge by
deriving a theoretical upper bound for the coefficient of determination (R2)
for regression models. This theoretical upper bound depends only on the noise
associated with the response variable in a dataset as well as its variance. The
upper bound estimate was validated via Monte Carlo simulations and then used as
a tool to bootstrap performance of regression models trained on biological
datasets, including protein sequence data, transcriptomic data, and genomic
data. Although we study biological datasets in this work, the new upper bound
estimates will hold true for regression models from any research field or
application area where response variables have associated noise.


  Infection by distinct Dengue virus serotypes and host immunity are
intricately linked. In particular, certain levels of cross-reactive antibodies
in the host may actually enhance infection severity leading to Dengue
hemorrhagic fever (DHF). The coupled immunological and epidemiological dynamics
of Dengue calls for a multi-scale modeling approach. In this work, we formulate
a within-host model which mechanistically recapitulates characteristics of
antibody dependent enhancement (ADE) in Dengue infection. The within-host scale
is then linked to epidemiological spread by a vector-host partial differential
equation model structured by host antibody level. The coupling allows for
dynamic population-wide antibody levels to be tracked through primary and
secondary infections by distinct Dengue strains, along with waning of
cross-protective immunity after primary infection. Analysis of both the
within-host and between-host systems are conducted. Stability results in the
epidemic model are formulated via basic and invasion reproduction numbers as a
function of immunological variables. Additionally, we develop numerical methods
in order to simulate the multi-scale model and assess the influence of
parameters on disease spread and DHF prevalence in the population.


  A current challenge for disease modeling and public health is understanding
pathogen dynamics across scales since their ecology and evolution ultimately
operate on several coupled scales. This is particularly true for vector-borne
diseases, where within-vector, within-host, and between vector-host populations
all play crucial roles in diversity and distribution of the pathogen. Despite
recent modeling efforts to determine the effect of within-host virus-immune
response dynamics on between-host transmission, the role of within-vector viral
dynamics on disease spread is overlooked. Here we formulate an
age-since-infection structured epidemic model coupled to nonlinear ordinary
differential equations describing within-host immune-virus dynamics and
within-vector viral kinetics, with feedbacks across these scales. We first
define the \emph{within-host viral-immune response and within-vector viral
kinetics dependent} basic reproduction number $\mathcal R_0.$ Then we prove
that whenever $\mathcal R_0<1,$ the disease free equilibrium is locally
asymptotically stable, and under certain biologically interpretable conditions,
globally asymptotically stable. Otherwise if $\mathcal R_0>1,$ it is unstable
and the system has a unique positive endemic equilibrium. In the special case
of constant vector to host inoculum size, we show the positive equilibrium is
locally asymptotically stable and the disease is weakly uniformly persistent.
Furthermore numerical results suggest that within-vector-viral kinetics and
dynamic inoculum size may play a substantial role in epidemics. Finally, we
address how the model can be utilized to better predict the success of control
strategies such as vaccination and drug treatment.


  Understanding the causes and effects of spatial vegetation patterns is a
fundamental problem in ecology, especially because these can be used as early
predictors of catastrophic shifts such as desertification processes. Empirical
studies of the vegetation cover in some areas such as drylands and semiarid
regions have revealed the existence of vegetation patches of broadly diverse
sizes. In particular, the probability distribution of patch sizes can be fitted
by a power law, i.e. vegetation patches are approximately scale free up to some
maximum size. Different explanatory mechanisms, such as plant-plant
interactions and plant-water feedback loops have been proposed to rationalize
the emergence of such scale-free patterns, yet a full understanding has not
been reached. Using a simple model for vegetation dynamics, we show that
environmental temporal variability -- a well-recognized feature of semiarid
environments -- promotes in a robust way (i.e. for a wide range of parameter
values) the emergence of vegetation patches with broadly distributed cluster
sizes. Furthermore, this result is related to a percolation phenomenon that
occurs in an intermittent or fluctuating way. The model also reveals that the
power-law exponents fitting the tails of the probability distributions depend
on the overall vegetation-cover density, in agreement with empirical
observations. This supports the idea that environmental variability plays a key
role in the formation of scale-free vegetation patterns. From a practical
viewpoint, this may be of importance to predict the effects that changes in
environmental conditions may have in real ecosystems. From a theoretical side,
our study sheds new light on a novel type of percolation phenomena occurring
under temporally-varying external conditions, that still needs further work to
be fully characterized.


  In this paper, experimental data is presented and a simple model is developed
for the time evolution of a F-amp \textit{E. Coli} culture population. In
general, the bacteria life cycle as revealed by monitoring a culture's
population consists of the lag phase, the growth (or exponential) phase, the
log (or stationary) phase, and finally the death phase. As the name suggests,
in the stationary phase, the population of the bacteria ceases to grow
exponentially and reaches a plateau before beginning the death phase. High
temporal resolution experimental observations using a unique light-scattering
technique in this work reveal all the expected phases in detail as well as an
oscillatory population behavior in the stationary phase. This unambiguous
oscillation behavior has been suggested previously using traditional surveys of
aliquots from a given population culture. An attempt is made to model these
experimental results by developing a differential equation that accounts for
the spatial distribution of the individual cells and the presence of the
self-organizing forces of competition and dispersion. The main phases are well
represented, and the oscillating behavior is attributed to intra-species
mixing. It is also observed, that the convective motion arising out of
intra-species mixing while plays a key role in limiting population growth,
scales as $t^{-\alpha}$, where $\alpha$ is bounded by model parameters.


  There is increasing evidence that protein binding to specific sites along DNA
can activate the reading out of genetic information without coming into direct
physical contact with the gene. There also is evidence that these distant but
interacting sites are embedded in a liquid droplet of proteins which condenses
out of the surrounding solution. We argue that droplet-mediated interactions
can account for crucial features of gene regulation only if the droplet is
poised at a non-generic point in its phase diagram. We explore a minimal model
that embodies this idea, show that this model has a natural mechanism for
self-tuning, and suggest direct experimental tests.


  We study a hybrid impulsive reaction-advection-diffusion model given by a
reaction-advection-diffusion equation composed with a discrete-time map in
space dimension $n\in\mathbb N$. The reaction-advection-diffusion equation
takes the form \begin{equation*}\label{} u^{(m)}_t = \text{div}(A\nabla
u^{(m)}-q u^{(m)}) + f(u^{(m)}) \quad \text{for} \ \ (x,t)\in\mathbb R^n \times
(0,1] , \end{equation*} for some function $f$, a drift $q$ and a diffusion
matrix $A$. When the discrete-time map is local in space we use $N_m(x)$ to
denote the density of population at a point $x$ at the beginning of
reproductive season in the $m$th year and when the map is nonlocal we use
$u_m(x)$. The local discrete-time map is \begin{eqnarray*}\label{}\left\{
\begin{array}{lcl} u^{(m)}(x,0) = g(N_m(x)) \quad \text{for} \ \ x\in \mathbb
R^n , \\ N_{m+1}(x):=u^{(m)}(x,1) \quad \text{for} \ \ x\in \mathbb R^n ,
\end{array}\right. \end{eqnarray*} for some function $g$. The nonlocal discrete
time map is \begin{eqnarray*}\label{}\left\{ \begin{array}{lcl} u^{(m)}(x,0) =
u_{m}(x) \quad \text{for} \ \ x\in \mathbb R^n , \\ \label{mainb2} u_{m+1}(x)
:= g\left(\int_{\mathbb R^n} K(x-y)u^{(m)}(y,1) dy\right) \quad \text{for} \ \
x\in \mathbb R^n, \end{array}\right. \end{eqnarray*} when $K$ is a nonnegative
normalized kernel. SEE THE ARTICLE FOR COMPLETE ABSTRACT.


  As genome sequencing tools and techniques improve, researchers are able to
incrementally assemble more accurate reference genomes, which enable
sensitivity in read mapping and downstream analysis such as variant calling. A
more sensitive downstream analysis is critical for a better understanding of
the genome donor (e.g., health characteristics). Therefore, read sets from
sequenced samples should ideally be mapped to the latest available reference
genome that represents the most relevant population. Unfortunately, the
increasingly large amount of available genomic data makes it prohibitively
expensive to fully re-map each read set to its respective reference genome
every time the reference is updated. There are several tools that attempt to
accelerate the process of updating a read data set from one reference to
another (i.e., remapping). However, if a read maps to a region in the old
reference that does not appear with a reasonable degree of similarity in the
new reference, the read cannot be remapped. We find that, as a result of this
drawback, a significant portion of annotations are lost when using
state-of-the-art remapping tools. To address this major limitation in existing
tools, we propose AirLift, a fast and comprehensive technique for remapping
alignments from one genome to another. Compared to the state-of-the-art method
for remapping reads (i.e., full mapping), AirLift reduces the overall execution
time to remap read sets between two reference genome versions by up to 27.4x.
We validate our remapping results with GATK and find that AirLift provides high
accuracy in identifying ground truth SNP/INDEL variants.


  Cancer diagnosis, prognosis, and therapeutic response predictions are based
on morphological information from histology slides and molecular profiles from
genomic data. However, most deep learning-based objective outcome prediction
and grading paradigms are based on histology or genomics alone and do not make
use of the complementary information in an intuitive manner. In this work, we
propose Pathomic Fusion, an interpretable strategy for end-to-end multimodal
fusion of histology image and genomic (mutations, CNV, RNA-Seq) features for
survival outcome prediction. Our approach models pairwise feature interactions
across modalities by taking the Kronecker product of unimodal feature
representations and controls the expressiveness of each representation via a
gating-based attention mechanism. Following supervised learning, we are able to
interpret and saliently localize features across each modality, and understand
how feature importance shifts when conditioning on multimodal input. We
validate our approach using glioma and clear cell renal cell carcinoma datasets
from the Cancer Genome Atlas (TCGA), which contains paired whole-slide image,
genotype, and transcriptome data with ground truth survival and histologic
grade labels. In a 15-fold cross-validation, our results demonstrate that the
proposed multimodal fusion paradigm improves prognostic determinations from
ground truth grading and molecular subtyping, as well as unimodal deep networks
trained on histology and genomic data alone. The proposed method establishes
insight and theory on how to train deep networks on multimodal biomedical data
in an intuitive manner, which will be useful for other problems in medicine
that seek to combine heterogeneous data streams for understanding diseases and
predicting response and resistance to treatment.


  The role of rotational molecular motors of the ATP synthase class is integral
to the metabolism of cells. Yet the function of FliI6-FliJ complex - homologous
to the F1 ATPase motor - within the flagellar export apparatus remains unclear.
We use a simple two-state model adapted from studies of linear molecular motors
to identify key features of this motor. The two states are the 'locked' ground
state where the FliJ coiled coil filament experiences fluctuations in an
asymmetric torsional potential, and a 'free' excited state in which FliJ
undergoes rotational diffusion. Michaelis-Menten kinetics was used to treat
transitions between these two states, and obtain the average angular velocity
of the FliJ filament within the FliI6 stator: ~9.0 rps. The motor was then
studied under external counter torque conditions in order to ascertain its
maximal power output: Pmax ~42 kBT/s, and the stall torque: ~3 kBT/rad. Two
modes of action within the flagellar export apparatus are proposed, in which
the motor performs useful work either by continuously 'grinding' through the
resistive environment, or by exerting equal and opposite stall force on it. In
both cases, the resistance is provided by flagellin subunits entering the
flagellar export channel prior to their unfolding. We therefore propose that
the function of the FliI6-FliJ complex is to lower the energy barrier and
therefore assist in unfolding of the flagellar proteins before feeding them
into the transport channel.


  Inferring concerted changes among biological traits along an evolutionary
history remains an important yet challenging problem. Besides adjusting for
spurious correlation induced from the shared history, the task also requires
sufficient flexibility and computational efficiency to incorporate multiple
continuous and discrete traits as data size increases. To accomplish this, we
jointly model mixed-type traits by assuming latent parameters for binary
outcome dimensions at the tips of an unknown tree informed by molecular
sequences. This gives rise to a phylogenetic multivariate probit model. With
large sample sizes, posterior computation under this model is problematic, as
it requires repeated sampling from a high-dimensional truncated normal
distribution. Current best practices employ multiple-try rejection sampling
that suffers from slow-mixing and a computational cost that scales
quadratically in sample size. We develop a new inference approach that exploits
1) the bouncy particle sampler (BPS) based on piecewise deterministic Markov
processes to simultaneously sample all truncated normal dimensions, and 2)
novel dynamic programming that reduces the cost of likelihood and gradient
evaluations for BPS to linear in sample size. In an application with 535 HIV
viruses and 24 traits that necessitates sampling from a 12,840-dimensional
truncated normal, our method makes it possible to estimate the across-trait
correlation and detect factors that affect the pathogen's capacity to cause
disease. This inference framework is also applicable to a broader class of
covariance structures beyond comparative biology.


  Fluctuating environments are situations where the spatio-temporal
stochasticity plays a significant role in the evolutionary dynamics. The study
of the evolution of cooperation in these environments typically assumes a
homogeneous, well mixed population, whose constituents are endowed with
identical capabilities. In this paper, we generalize these results by
developing a systematic study for the cooperation dynamics in fluctuating
environments under the consideration of structured, heterogeneous populations
with individual entities subjected to general behavioral rules. Considering
complex network topologies, and a behavioral rule based on generalized
reciprocity, we perform a detailed analysis of the effect of the underlying
interaction structure on the evolutionary stability of cooperation. We find
that, in the presence of environmental fluctuations, the cooperation dynamics
can lead to the creation of multiple network components, each with distinct
evolutionary properties. This is paralleled to the freezing state in the Random
Energy Model. We utilize this result to examine the applicability of our
generalized reciprocity behavioral rule in a variety of settings. We thereby
show that the introduced rule leads to steady state cooperative behavior that
is always greater than or equal to the one predicted by the evolutionary
stability analysis of unconditional cooperation. As a consequence, the
implementation of our results may go beyond explaining the evolution of
cooperation. In particular, they can be directly applied in domains that deal
with the development of artificial systems able to adequately mimic reality,
such as reinforcement learning.


  The mechanisms underlying eco-evolutionary dynamics (the feedback between
ecological and evolutionary processes) are often unknown. Here, we propose that
classical theory from behavioral ecology can provide a greater understanding of
the mechanisms underlying eco-evolutionary dynamics, and thus improve
predictions about the outcomes of these dynamics.


  We provide easy and readable GNU Octave/MATLAB code for the simulation of
mathematical models described by ordinary differential equations and for the
solution of optimal control problems through Pontryagin's maximum principle.
For that, we consider a normalized HIV/AIDS transmission dynamics model based
on the one proposed in our recent contribution (Silva, C.J.; Torres, D.F.M. A
SICA compartmental model in epidemiology with application to HIV/AIDS in Cape
Verde. Ecol. Complex. 2017, 30, 70--75), given by a system of four ordinary
differential equations. An HIV initial value problem is solved numerically
using the ode45 GNU Octave function and three standard methods implemented by
us in Octave/MATLAB: Euler method and second-order and fourth-order Runge-Kutta
methods. Afterwards, a control function is introduced into the normalized HIV
model and an optimal control problem is formulated, where the goal is to find
the optimal HIV prevention strategy that maximizes the fraction of uninfected
HIV individuals with the least HIV new infections and cost associated with the
control measures. The optimal control problem is characterized analytically
using the Pontryagin Maximum Principle, and the extremals are computed
numerically by implementing a forward-backward fourth-order Runge-Kutta method.
Complete algorithms, for both uncontrolled initial value and optimal control
problems, developed under the free GNU Octave software and compatible with
MATLAB are provided along the article.


  Securities markets are quintessential complex adaptive systems in which
heterogeneous agents compete in an attempt to maximize returns. Species of
trading agents are also subject to evolutionary pressure as entire classes of
strategies become obsolete and new classes emerge. Using an agent-based model
of interacting heterogeneous agents as a flexible environment that can
endogenously model many diverse market conditions, we subject deep neural
networks to evolutionary pressure to create dominant trading agents. After
analyzing the performance of these agents and noting the emergence of anomalous
superdiffusion through the evolutionary process, we construct a method to turn
high-fitness agents into trading algorithms. We backtest these trading
algorithms on real high-frequency foreign exchange data, demonstrating that
elite trading algorithms are consistently profitable in a variety of market
conditions---even though these algorithms had never before been exposed to real
financial data. These results provide evidence to suggest that developing
\textit{ab initio} trading strategies by repeated simulation and evolution in a
mechanistic market model may be a practical alternative to explicitly training
models with past observed market data.


  Species distribution models (SDMs) are important statistical tools for
ecologists to understand and predict species range. However, standard SDMs do
not explicitly incorporate dynamic processes like dispersal. This limitation
may lead to bias in inference about species distribution. Here, we adopt the
theory of ecological diffusion that has recently been introduced in statistical
ecology to incorporate spatio-temporal processes in ecological models. As a
case study, we considered the wolf (Canis lupus) that has been recolonizing
Eastern France naturally through dispersal from the Apennines since the early
90's. Using partial differential equations for modelling species diffusion and
growth in a fragmented landscape, we develop a mechanistic-statistical
spatio-temporal model accounting for ecological diffusion, logistic growth and
imperfect species detection. We conduct a simulation study and show the ability
of our model to i) estimate ecological parameters in various situations with
contrasted species detection probability and number of surveyed sites and ii)
forecast the distribution into the future. We found that the growth rate of the
wolf population in France was explained by the proportion of forest cover, that
diffusion was influenced by human density and that species detectability
increased with increasing survey effort. Using the parameters estimated from
the 2007-2015 period, we then forecasted wolf distribution in 2016 and found
good agreement with the actual detections made that year. Our approach may be
useful for managing species that interact with human activities to anticipate
potential conflicts.


  We study here the extreme statistics of Brownian particles escaping from a
cusp funnel: the fastest Brownian particles among $n$ follow an ensemble of
optimal trajectories located near the shortest path from the source to the
target. For the time of such first arrivers, we derive an asymptotic formula
that differs from the classical narrow escape and dire strait obtained for the
mean first passage time. Consequently, when particles are initially distributed
at a given distance from a cusp, the fastest do see some properties
characterizing the cusp geometry. Therefore, when many particles diffuse around
impermeable obstacles, the geometry plays a role in the time to reach a target.
In the biological context of cellular transduction with signalling molecules,
having to escape such cusp-like domains slows down fast signaling. To conclude,
generating multiple copies of the same molecule helps bypass a crowded
environment to transmit a molecular signal quickly.


  We show that weakly reversible mass-action systems can have a continuum of
positive steady states, coming from the zeroes of a multivariate polynomial.
Moreover, the same is true of systems whose underlying reaction network is
reversible and has a single connected component. In our construction, we relate
operations on the reaction network to the multivariate polynomial occurring as
a common factor in the system of differential equations.


  Bacteriophages are viruses infecting bacteria and archaea. Many phage species
cause infections which lead to the certain death of the infected prokaryotic
host cell and the release of a large batch of phage progeny, yet they have been
able to stably coexist with their bacterial hosts over the eons in nature, as
well as in the majority of laboratory experiments reported in the literature.
This possibility of a stable coexistence between populations of bacteria and
bacteriophages can be suspected to critically reduce the chances of a
successful therapeutic application of bacteriophages for combatting pathogenic
bacterial infections in vivo. Here, we extend an established differential
equation model describing the interaction of bacteria, bacteriophages, and the
host immune system, modelling different degrees of spatial heterogeneity of the
host organism by introducing a scaling parameter which alters the encounter
rates of the different cell populations. We demonstrate by rigorous
mathematical analysis that, depending on the degree of spatial heterogeneity,
the system will either converge to the desired state of bacterial and phage
extinction, to a state of persisting bacterial infection with a completely
eliminated bacteriophage population or to an equally undesirable state of a
stable long-term bacteria-bacteriophage coexistence. We additionally provide
numerical solutions of the model to illustrate the emerging dynamics and
discuss some implications for bacteriophage therapy.


  The integral control of positive systems using nonnegative control input is
an important problem arising, among others, in biochemistry, epidemiology and
ecology. An immediate solution is to use an ON-OFF nonlinearity between the
controller and the system. However, this solution is only available when
controllers are implemented in computer systems. When this is not the case,
like in biology, alternative approaches need to be explored. Based on recent
research in the control of biological systems, we propose to develop a theory
for the integral control of positive systems using nonnegative controls based
on the so-called \emph{antithetic integral controller} and two \emph{positively
regularized integral controllers}, the so-called \emph{exponential integral
controller} and \emph{logistic integral controller}. For all these controllers,
we establish several qualitative results, which we connect to standard results
on integral control. We also obtain additional results which are specific to
the type of controllers. For instance, we show an interesting result
stipulating that if the gain of the antithetic integral controller is suitably
chosen, then the local stability of the equilibrium point of the closed-loop
system does not depend on the choice for the coupling parameter, an additional
parameter specific to this controller. Conversely, we also show that if the
coupling parameter is suitably chosen, then the equilibrium point of the
closed-loop system is locally stable regardless the value of the gain. For the
exponential integral controller, we can show that the local stability of the
equilibrium point of the closed-loop system is independent of the gain of the
controller and the gain of the system. The stability only depends on the
exponential rate of the controller, again a parameter that is specific to this
type of controllers. Several examples are given for illustration.


  The huge amount of data acquired by high-throughput sequencing requires data
reduction for effective analysis. Here we give a clustering algorithm for
genome-wide open chromatin data using a new data reduction method. This method
regards the genome as a string of $1$s and $0$s based on a set of peaks and
calculates the Hamming distances between the strings. This algorithm with the
systematically optimized set of peaks enables us to quantitatively evaluate
differences between samples of hematopoietic cells and classify cell types,
potentially leading to a better understanding of leukemia pathogenesis.


  Resistance to insecticide is considered nowadays one of the major threats to
insect control, as its occurrence reduces drastically the efficiency of
chemical control campaigns, and may also perturb the application of other
control methods, like biological and genetic control. In order to account for
the emergence and spread of such phenomenon as an effect of exposition to
larvicide and/or adulticide, we develop in this paper a general time-continuous
population model with two life phases, subsequently simplified through slow
manifold theory. The derived models present density-dependent recruitment and
mortality rates in a non-conventional way. We show that in absence of
selection, they evolve in compliance with Hardy-Weinberg law; while in presence
of selection and in the dominant or codominant cases, convergence to the
fittest genotype occurs. The proposed mathematical models should allow for the
study of several issues of importance related to the use of insecticides and
other adaptive phenomena.


  Motivation: Laboratory gene regulatory data for a species are sporadic.
Despite the abundance of gene regulatory network algorithms that employ single
data sets, few algorithms can combine the vast but disperse sources of data and
extract the potential information. With a motivation to compensate for this
shortage, we developed an algorithm called GENEREF that can accumulate
information from multiple types of data sets in an iterative manner, with each
iteration boosting the performance of the prediction results.
  Results: The algorithm is examined extensively on data extracted from the
quintuple DREAM4 networks and DREAM5's Escherichia coli and Saccharomyces
cerevisiae networks and sub-networks. Many single-dataset and multi-dataset
algorithms were compared to test the performance of the algorithm. Results show
that GENEREF surpasses non-ensemble state-of-the-art multi-perturbation
algorithms on the selected networks and is competitive to present
multiple-dataset algorithms. Specifically, it outperforms dynGENIE3 and is on
par with iRafNet. Also, we argued that a scoring method solely based on the
AUPR criterion would be more trustworthy than the traditional score.
  Availability: The Python implementation along with the data sets and results
can be downloaded from github.com/msaremi/GENEREF


  Identifying protein-protein interactions is crucial for a systems-level
understanding of the cell. Recently, algorithms based on inverse statistical
physics, e.g. Direct Coupling Analysis (DCA), have allowed to use
evolutionarily related sequences to address two conceptually related inference
tasks: finding pairs of interacting proteins, and identifying pairs of residues
which form contacts between interacting proteins. Here we address two
underlying questions: How are the performances of both inference tasks related?
How does performance depend on dataset size and the quality? To this end, we
formalize both tasks using Ising models defined over stochastic block models,
with individual blocks representing single proteins, and inter-block couplings
protein-protein interactions; controlled synthetic sequence data are generated
by Monte-Carlo simulations. We show that DCA is able to address both inference
tasks accurately when sufficiently large training sets are available, and that
an iterative pairing algorithm (IPA) allows to make predictions even without a
training set. Noise in the training data deteriorates performance. In both
tasks we find a quadratic scaling relating dataset quality and size that is
consistent with noise adding in square-root fashion and signal adding linearly
when increasing the dataset. This implies that it is generally good to
incorporate more data even if its quality is imperfect, thereby shedding light
on the empirically observed performance of DCA applied to natural protein
sequences.


  We study the spread of an SIRS-type epidemic with vaccination on network.
Starting from an exact Markov description of the model, we investigate the mean
epidemic lifetime by providing a sufficient condition for fast extinction that
depends on the model parameters and the topology of the network. Then, we pass
to consider a firstorder mean-field approximation of the exact model and its
stability properties, by relying on the graph-theoretical notion of equitable
partition. In the case of graphs possessing this kind of partition, we prove
that the endemic equilibrium can be computed by using a lower-dimensional
dynamical system. Finally, in the special case of regular graphs, we
investigate the domain of attraction of the endemic equilibrium.


  Functional or non-coding RNAs are attracting more attention as they are now
potentially considered valuable resources in the development of new drugs
intended to cure several human diseases. The identification of drugs targeting
the regulatory circuits of functional RNAs depends on knowing its family, a
task which is known as RNA sequence classification. State-of-the-art small
noncoding RNA classification methodologies take secondary structural features
as input. However, in such classification, feature extraction approaches only
take global characteristics into account and completely oversight co-relative
effect of local structures. Furthermore secondary structure based approaches
incorporate high dimensional feature space which proves computationally
expensive. This paper proposes a novel Robust and Precise ConvNet (RPC-snRC)
methodology which classifies small non-coding RNAs sequences into their
relevant families by utilizing the primary sequence of RNAs. RPC-snRC
methodology learns hierarchical representation of features by utilizing
positioning and occurrences information of nucleotides. To avoid exploding and
vanishing gradient problems, we use an approach similar to DenseNet in which
gradient can flow straight from subsequent layers to previous layers. In order
to assess the effectiveness of deeper architectures for small non-coding RNA
classification, we also adapted two ResNet architectures having different
number of layers. Experimental results on a benchmark small non-coding RNA
dataset show that our proposed methodology does not only outperform existing
small non-coding RNA classification approaches with a significant performance
margin of 10% but it also outshines adapted ResNet architectures.


  In this paper we study self-similar and fractal networks from the
combinatorial perspective. We establish analogues of topological (Lebesgue) and
fractal (Hausdorff) dimensions for graphs and demonstrate that they are
naturally related to known graph-theoretical characteristics: rank dimension
and product (or Prague or Ne\v{s}et\v{r}il-R\"odl) dimension. Our approach
reveals how self-similarity and fractality of a network are defined by a
pattern of overlaps between densely connected network communities. It allows us
to identify fractal graphs, explore the relations between graph fractality,
graph colorings and graph Kolmogorov complexity, and analyze the fractality of
several classes of graphs and network models, as well as of a number of
real-life networks. We demonstrate the application of our framework to
evolutionary studies by revealing the growth of self-organization of
heterogeneous viral populations over the course of their intra-host evolution,
thus suggesting mechanisms of their gradual adaptation to the host's
environment. As far as the authors know, the proposed approach is the first
theoretical framework for study of network fractality within the combinatorial
paradigm. The obtained results lay a foundation for studying fractal properties
of complex networks using combinatorial methods and algorithms.


  Artificial neurons built on synthetic gene networks have potential
applications ranging from complex cellular decision-making to bioreactor
regulation. Furthermore, due to the high information throughput of natural
systems, it provides an interesting candidate for biologically-based
supercomputing and analog simulations of traditionally intractable problems. In
this paper, we propose an architecture for constructing multicellular neural
networks and programmable nonlinear systems. We design an artificial neuron
based on gene regulatory networks and optimize its dynamics for modularity.
Using gene expression models, we simulate its ability to perform arbitrary
linear classifications from multiple inputs. Finally, we construct a two-layer
neural network to demonstrate scalability and nonlinear decision boundaries and
discuss future directions for utilizing uncontrolled neurons in computational
tasks.


  Molecular fingerprints are the workhorse in ligand-based drug discovery. In
recent years, an increasing number of research papers reported fascinating
results on using deep neural networks to learn 2D molecular representations as
fingerprints. It is anticipated that the integration of deep learning would
also contribute to the prosperity of 3D fingerprints. Here, we unprecedentedly
introduce deep learning into 3D small molecule fingerprints, presenting a new
one we termed as the three-dimensional force fields fingerprint (TF3P). TF3P is
learned by a deep capsular network whose training is in no need of labeled
datasets for specific predictive tasks. TF3P can encode the 3D force fields
information of molecules and demonstrates the stronger ability to capture 3D
structural changes, to recognize molecules alike in 3D but not in 2D, and to
identify similar targets inaccessible by other 2D or 3D fingerprints based on
only ligands similarity. Furthermore, TF3P is compatible with both statistical
models (e.g. similarity ensemble approach) and machine learning models.
Altogether, we report TF3P as a new 3D small molecule fingerprint with a
promising future in ligand-based drug discovery. All codes are written in
Python and available at https://github.com/canisw/tf3p.


  Personalizing drug prescriptions in cancer care based on genomic information
requires associating genomic markers with treatment effects. This is an
unsolved challenge requiring genomic patient data in yet unavailable volumes as
well as appropriate quantitative methods. We attempt to solve this challenge
for an experimental proxy for which sufficient data is available: 42 drugs
tested on 1018 cancer cell lines. Our goal is to develop a method to identify
the drug that is most promising based on a cell line's genomic information. For
this, we need to identify for each drug the machine learning method, choice of
hyperparameters and genomic features for optimal predictive performance. We
extensively compare combinations of gene sets (both curated and random),
genetic features, and machine learning algorithms for all 42 drugs. For each
drug, the best performing combination (considering only the curated gene sets)
is selected. We use these top model parameters for each drug to build and
demonstrate a Drug Recommendation System (Dr.S). Insights resulting from this
analysis are formulated as best practices for developing drug recommendation
systems. The complete software system, called the Cell Line Analyzer, is
written in Python and available on github.


  In this perspective article, we present a multidisciplinary approach for
characterizing protein structure networks. We first place our approach in its
historical context and describe the manner in which it synthesizes concepts
from quantum chemistry, biology of polymer conformations, matrix mathematics,
and percolation theory. We then explicitly provide the method for constructing
the protein structure network in terms of non-covalently interacting amino acid
side chains and show how a mine of information can be obtained from the graph
spectra of these networks. Employing suitable mathematical approaches, such as
the use of a weighted, Laplacian matrix to generate the spectra, enables us to
develop rigorous methods for network comparison and to identify crucial nodes
responsible for the network integrity through a perturbation approach. Our
scoring methods have several applications in structural biology that are
elusive to conventional methods of analyses. Here, we discuss the instances of:
(a) Protein structure comparison that include the details of side chain
connectivity, (b) The contribution to node clustering as a function of bound
ligand, explaining the global effect of local changes in phenomena such as
allostery and (c) The identification of crucial amino acids for structural
integrity, derived purely from the spectra of the graph. We demonstrate how our
method enables us to obtain valuable information on key proteins involved in
cellular functions and diseases such as GPCR and HIV protease, and discuss the
biological implications. We then briefly describe how concepts from percolation
theory further augment our analyses. In our concluding perspective for future
developments, we suggest a further unifying approach to protein structure
analyses and a judicious choice of questions to employ our methods for larger,
more complex networks, such as metabolic and disease networks.


  The DNA double helix is a polyanionic macromolecule that in water solutions
is neutralized by metal ions (counterions). The property of the counterions to
stabilize the water network (positive hydration) or to make it friable
(negative hydration) is important in terms of the physical mechanisms of
stabilization of the DNA double helix. In the present research, the effects of
positive hydration of Na$^{+}$ counterions and negative hydration of K$^{+}$
and Cs$^{+}$ counterions, incorporated into the hydration shell of the DNA
double helix have been studied using molecular dynamics simulations. The
results have shown that the dynamics of the hydration shell of counterions
depends on region of the double helix: minor groove, major groove, and outside
the macromolecule. The longest average residence time has been observed for
water molecules contacting with the counterions, localized in the minor groove
of the double helix (about 50 ps for Na$^{+}$, and lower than 10 ps for K$^{+}$
and Cs$^{+}$). The estimated potentials of mean force for the hydration shells
of the counterions show that the water molecules are constrained too strong,
and, consequently, the effect of negative hydration for K$^{+}$ and Cs$^{+}$
counterions has not been observed in the simulations. The analysis has shown
that the effects of counterion hydration can be described more accurately with
water models having lower dipole moments.


  In multicellular organisms, cells differentiate into several distinct types
during early development. Determination of each cellular state, along with the
ratio of each cell type, as well as the developmental course during cell
differentiation are highly regulated processes that are robust to noise and
environmental perturbations throughout development. Waddington metaphorically
depicted this robustness as the epigenetic landscape in which the robustness of
each cellular state is represented by each valley in the landscape. This
robustness is now conceptualized as an approach toward an attractor in a
gene-expression dynamical system. However, there is still an incomplete
understanding of the origin of landscape change, which is accompanied by
branching of valleys that corresponds to the differentiation process. Recent
progress in developmental biology has unveiled the molecular processes involved
in epigenetic modification, which will be a key to understanding the nature of
slow landscape change. Nevertheless, the contribution of the interplay between
gene expression and epigenetic modification to robust landscape changes, known
as homeorhesis, remains elusive. Here, we introduce a theoretical model that
combines epigenetic modification with gene expression dynamics driven by a
regulatory network. In this model, epigenetic modification changes the
feasibility of expression, i.e., the threshold for expression dynamics, and a
slow positive-feedback process from expression to the threshold level is
introduced. Under such epigenetic feedback, several fixed-point attractors with
distinct expression patterns are generated hierarchically shaping the
epigenetic landscape with successive branching of valleys. This theory provides
a quantitative framework for explaining homeorhesis in development as
postulated by Waddington, based on dynamical-system theory with slow feedback
reinforcement.


  Agriculture operates within a global ecosystem for which planetary boundaries
have recently been defined. Efficiency in nitrogen use is essential for
maximizing the benefits of agriculture for humanity and minimizing adverse
socio-ecological impacts. The set of variables that support global system
efficiency also determine the food production boundaries of agriculture, which
govern the maximum supportable human population. Food production boundaries,
nitrogen loss and nitrogen self-sufficiency are combined here into the nitrogen
operating space of world food production. We position world regions and the
world trajectory (1961-2013) within the nitrogen operating space and show that
the maximum supportable human population ranges from 6 to almost 17 billion
people according to the share of grain used as feed and the nitrogen
fertilization regime. All UN population projections for the 21st century can
only be conditionally achieved. We discuss the growth rate requirements in
production and efficiency to meet food production boundaries and the nitrogen
planetary boundary by 2050.


  Mining gene expression profiles has proven valuable for identifying
signatures serving as surrogates of cancer phenotypes. However, the
similarities of such signatures across different cancer types have not been
strong enough to conclude that they represent a universal biological mechanism
shared among multiple cancer types. Here we describe a network-based approach
that explores gene-to-gene connections in multiple cancer datasets while
maximizing the overall association of the subnetwork with clinical outcomes.
With the dataset of The Cancer Genome Atlas (TCGA), we studied the
characteristics of common gene expression of three types of cancers: Rectum
adenocarcinoma (READ), Breast invasive carcinoma (BRCA) and Colon
adenocarcinoma (COAD). By analyzing several pairs of highly correlated genes
after filtering and clustering work, we found that the co-expressed genes
across multiple types of cancers point to particular biological mechanisms
related to cancer cell progression , suggesting that they represent important
attributes of cancer in need of being elucidated for potential applications in
diagnostic, prognostic and therapeutic products applicable to multiple cancer
types.


  For many stochastic models of interest in systems biology, such as those
describing biochemical reaction networks, exact quantification of parameter
uncertainty through statistical inference is intractable. Likelihood-free
computational inference techniques enable parameter inference when the
likelihood function for the model is intractable but the generation of many
sample paths is feasible through stochastic simulation of the forward problem.
The most common likelihood-free method in systems biology is approximate
Bayesian computation that accepts parameters that result in low discrepancy
between stochastic simulations and measured data. However, it can be difficult
to assess how the accuracy of the resulting inferences are affected by the
choice of acceptance threshold and discrepancy function. The pseudo-marginal
approach is an alternative likelihood-free inference method that utilises a
Monte Carlo estimate of the likelihood function. This approach has several
advantages, particularly in the context of noisy, partially observed,
time-course data typical in biochemical reaction network studies. Specifically,
the pseudo-marginal approach facilitates exact inference and uncertainty
quantification, and may be efficiently combined with particle filters for low
variance, high-accuracy likelihood estimation. In this review, we provide a
practical introduction to the pseudo-marginal approach using inference for
biochemical reaction networks as a series of case studies. Implementations of
key algorithms and examples are provided using the Julia programming language;
a high performance, open source programming language for scientific computing.


  Structural information about protein-protein interactions, often missing at
the interactome scale, is important for mechanistic understanding of cells and
rational discovery of therapeutics. Protein docking provides a computational
alternative to predict such information. However, ranking near-native docked
models high among a large number of candidates, often known as the scoring
problem, remains a critical challenge. Moreover, estimating model quality, also
known as the quality assessment problem, is rarely addressed in protein
docking.
  In this study the two challenging problems in protein docking are regarded as
relative and absolute scoring, respectively, and addressed in one
physics-inspired deep learning framework. We represent proteins' and encounter
complexes' 3D structures as intra- and inter-molecular residue contact graphs
with atom-resolution node and edge features. And we propose a novel graph
convolutional kernel that pool interacting nodes' features through edge
features so that generalized interaction energies can be learned directly from
graph data. The resulting energy-based graph convolutional networks (EGCN) with
multi-head attention are trained to predict intra- and inter-molecular
energies, binding affinities, and quality measures (interface RMSD) for
encounter complexes. Compared to a state-of-the-art scoring function for model
ranking, EGCN has significantly improved ranking for a CAPRI test set involving
homology docking; and is comparable for Score_set, a CAPRI benchmark set
generated by diverse community-wide docking protocols not known to training
data. For Score_set quality assessment, EGCN shows about 27% improvement to our
previous efforts. Directly learning from 3D structure data in graph
representation, EGCN represents the first successful development of graph
convolutional networks for protein docking.


  Empirical studies of graphs have contributed enormously to our understanding
of complex systems. Known today as network science, what was originally a
theoretical study of graphs has grown into a more scientific exploration of
communities spanning the physical, biological, and social. However, as the
quantity and types of networks have grown so has their heterogeneity in quality
and specificity. This has hampered efforts to develop general network theory
capable of inferring functioning and predicting dynamics across study systems.
We have successfully approached this challenge by aligning networks to each
other rather than comparing parameter estimates from individually fitted models
or properties of edge topologies. By comparing the predictability of
statistical dynamics originating from each network's constituent nodes we were
able to build a functional classifier that distinguished underlying processes
in both synthetic and real-world network data spanning the entire biological
scale from cellular machinery to ecosystems.


  Predicting compound-protein affinity is critical for accelerating drug
discovery. Recent progress made by machine learning focuses on accuracy but
leaves much to be desired for interpretability. Through molecular contacts
underlying affinities, our large-scale interpretability assessment finds
commonly-used attention mechanisms inadequate. We thus formulate a hierarchical
multi-objective learning problem whose predicted contacts form the basis for
predicted affinities. We further design a physics-inspired deep relational
network, DeepRelations, with intrinsically explainable architecture.
Specifically, various atomic-level contacts or "relations" lead to
molecular-level affinity prediction. And the embedded attentions are
regularized with predicted structural contexts and supervised with partially
available training contacts. DeepRelations shows superior interpretability to
the state-of-the-art: without compromising affinity prediction, it boosts the
AUPRC of contact prediction 9.5, 16.9, 19.3 and 5.7-fold for the test,
compound-unique, protein-unique, and both-unique sets, respectively. Our study
represents the first dedicated model development and systematic model
assessment for interpretable machine learning of compound-protein affinity.


  Characterization of proteins with intrinsic or unfolded state disorder
comprises a new frontier in structural biology, requiring the characterization
of diverse and dynamic structural ensembles. We introduce a comprehensive
Bayesian framework, the Extended Experimental Inferential Structure
Determination (X-EISD) method, that calculates the maximum log-likelihood of a
protein structural ensemble by accounting for the uncertainties of a wide range
of experimental data and back-calculation models from structures, including NMR
chemical shifts, J-couplings, Nuclear Overhauser Effects, paramagnetic
relaxation enhancements, residual dipolar couplings, and hydrodynamic radii,
single molecule fluorescence F\"orster resonance energy transfer efficiencies
and small angle X-ray scattering intensity curves. We apply X-EISD to the drkN
SH3 unfolded state domain and show that certain experimental data types are
more influential than others for both eliminating structural ensemble models,
while also finding equally probable disordered ensembles that have alternative
structural properties that will stimulate further experiments to discriminate
between them.


  Metal nanoclusters (NCs), typically consisting of a few to tens of metal
atoms, bridge the gap between organometallic compounds and crystalline metal
nanoparticles. As their size approaches the Fermi wavelength of electrons,
metal NCs exhibit discrete energy levels, which in turn results in the
emergence of intriguing physical and chemical (or physicochemical) properties,
especially strong fluorescence. In the past few decades, dramatic growth has
been witnessed in the development of different types of noble metal NCs (mainly
AuNCs and AgNCs). However, compared with noble metals, copper is a relatively
earth-abundant and cost-effective metal. Theoretical and experimental studies
have shown that copper NCs (CuNCs) possess unique catalytic and
photoluminescent properties. In this context, CuNCs are emerging as a new class
of nontoxic, economic, and effective phosphors and catalysts, drawing
significant interest across the life and medical sciences. To highlight these
achievements, this review begins by providing an overview of a multitude of
factors that play central roles in the fluorescence of CuNCs. Additionally, a
critical perspective of how the aggregation of CuNCs can efficiently improve
the florescent stability, tunability, and intensity is also discussed.
Following, we present representative applications of CuNCs in detection and
bioimaging. Finally, we outline current challenges and our perspective on the
development of CuNCs.


  Prolonged lag time can be induced by starvation contributing to the
antibiotic tolerance of bacteria. We analyze the optimal lag time to survive
and grow the iterative and stochastic application of antibiotics. A simple
model shows that the optimal lag time exhibits a discontinuous transition when
the severeness of the antibiotic is increased. This suggests the possibility of
reducing tolerant bacteria by controlled usage of antibiotics application. When
the bacterial populations are able to have two phenotypes with different lag
times, the fraction of the second phenotype that has different lag time shows a
continuous transition. We then present a generic framework to investigate the
optimal lag time distribution for total population fitness for a given
distribution of the antibiotic application duration. The obtained optimal
distributions have multiple peaks for a wide range of the antibiotic
application duration distributions, including the case where the latter is
monotonically decreasing. The analysis supports the advantage in evolving
multiple, possibly discrete phenotypes in lag time for bacterial long-term
fitness.


  In neurons, neuropeptides are synthesized in the soma and are then
transported along the axon in dense core vesicles (DCVs). DCVs are captured in
varicosities located along the axon terminal called en passant boutons, which
are active terminal sites that accumulate and release neurotransmitters.
Recently developed experimental techniques allow for the estimation of the age
of DCVs in various locations in the axon terminal. Accurate simulation of the
mean age of DCVs in boutons requires the development of a model that would
account for resident, transiting-anterograde, and transiting-retrograde DCV
populations. In this paper, such a model is developed. The model is applied to
simulating DCV transport in Drosophila type II motoneurons. The model simulates
DCV transport and capture in the axon terminals and makes it possible to
predict the age density distribution of DCVs in en passant boutons as well as
DCV's mean age in boutons. The predicted prevalence of older organelles in
distal boutons may explain the "dying back" pattern of axonal degeneration
observed in dopaminergic neurons in Parkinson's disease. The predicted
difference of two hours between the age of older DCVs residing in distal
boutons and the age of younger DCVs residing in proximal boutons is consistent
with an approximate estimate of age difference deduced from experimental
observations. The age density of resident DCVs is found to be bimodal, which is
because DCVs are captured from two transiting states: the anterograde
transiting state that contains younger DCVs and the retrograde transiting state
that contains older DCVs.


  RNA structure prediction is a challenging problem, especially with
pseudoknots. Recently, there has been a shift from the classical minimum free
energy-based methods (MFE) to partition function-based ones that assemble
structures using base-pairing probabilities. Two examples of the latter group
are the popular maximum expected accuracy (MEA) method and the ProbKnot method.
ProbKnot is a fast heuristic that pairs nucleotides that are reciprocally most
probable pairing partners, and unlike MEA, can also predict structures with
pseudoknots. However, ProbKnot's full potential has been largely overlooked. In
particular, when introduced, it did not have an MEA-like hyperparameter that
can balance between positive predictive value (PPV) and sensitivity. We show
that a simple thresholded version of ProbKnot, which we call ThreshKnot, leads
to more accurate overall predictions by filtering out unlikely pairs whose
probabilities fall under a given threshold. We also show that on three
widely-used folding engines (RNAstructure, Vienna RNAfold, and CONTRAfold),
ThreshKnot always outperforms the much more involved MEA algorithm in (1) its
higher structure prediction accuracy, (2) its capability to predict
pseudoknots, and (3) its faster runtime and easier implementation. This
suggests that ThreshKnot should replace MEA as the default partition
function-based structure prediction algorithm. ThreshKnot is already available
in the widely used RNAstructure software package version 6.2 (released November
27, 2019): https://rna.urmc.rochester.edu/RNAstructure.html


  Host-parasitoid population dynamics is often probed using a
semi-discrete/hybrid modeling framework. Here, the update functions in the
discrete-time model connecting year-to-year changes in the population densities
are obtained by solving ordinary differential equations that mechanistically
describe interactions when hosts become vulnerable to parasitoid attacks. We
use this semi-discrete formalism to study two key spatial effects: local
movement (migration) of parasitoids between patches during the vulnerable
period; and yearly redistribution of populations across patches outside the
vulnerable period. Our results show that in the absence of any redistribution,
constant density-independent migration and parasitoid attack rates are unable
to stabilize an otherwise unstable host-parasitoid population dynamics.
Interestingly, inclusion of host redistribution (but not parasitoid
redistribution) before the start of the vulnerable period can lead to stable
coexistence of both species. Next, we consider a Type-III functional response
(parasitoid attack rate increases with host density), where the absence of any
spatial effects leads to a neutrally stable host-parasitoid equilibrium. As
before, density-independent parasitoid migration by itself is again
insufficient to stabilize the population dynamics and host redistribution
provides a stabilizing influence. Finally, we show that a Type-III functional
response combined with density-dependent parasitoid migration leads to stable
coexistence, even in the absence of population redistributions. In summary, we
have systematically characterized parameter regimes leading to stable/unstable
population dynamics with different forms of spatial heterogeneity coupled to
the parasitoid's functional response using mechanistically formulated
semi-discrete models.


  Efficient protein synthesis depends on the availability of charged tRNA
molecules. With 61 different codons, shifting the balance among the tRNA
abundances can lead to large changes in the protein synthesis rate. Previous
theoretical work has asked about the optimization of these abundances, and
there is some evidence that regulatory mechanisms bring cells close to this
optimum, on average. We formulate the tradeoff between the precision of control
and the efficiency of synthesis, asking for the maximum entropy distribution of
tRNA abundances consistent with a desired mean rate of protein synthesis. Our
analysis, using data from E. coli, indicates that reasonable synthesis rates
are consistent only with rather low entropies, so that the cell's regulatory
mechanisms must encode a large amount of information about the "correct" tRNA
abundances.


  Recent studies have shown that even in the absence of extrinsic stress, the
morphologically symmetrically dividing model bacteria Escherichia coli do not
generate offspring with equal reproductive fitness. Instead, daughter cells
exhibit asymmetric division times that converge to two distinct growth states.
This represents a limited senescence / rejuvenation process derived from
asymmetric division that is stable for hundreds of generations. It remains
unclear why the bacteria do not continue the senescence beyond this asymptote.
Although there are inherent fitness benefits for heterogeneity in population
growth rates, the two growth equilibria are surprisingly similar, differing by
a few percent. In this work we derive an explicit model for the growth of a
bacterial population with two growth equilibria, based on a generalized
Fibonacci recurrence, in order to quantify the fitness benefit of a limited
senescence process and examine costs associated with asymmetry that could
generate the observed behavior. We find that with simple saturating effects of
asymmetric partitioning of subcellular components, two distinct but similar
growth states may be optimal while providing evolutionarily significant fitness
advantages.


  A novel method to classify human cells is presented in this work based on the
transform-domain method on DNA methylation data. DNA methylation profile
variations are observed in human cells with the progression of disease stages,
and the proposal is based on this DNA methylation variation to classify normal
and disease cells including cancer cells. The cancer cell types investigated in
this work cover hepatocellular (sample size n = 40), colorectal (n = 44), lung
(n = 70) and endometrial (n = 87) cancer cells. A new pipeline is proposed
integrating the DNA methylation intensity measurements on all the CpG islands
by the transformation of Walsh-Hadamard Transform (WHT). The study reveals the
three-step properties of the DNA methylation transform-domain data and the step
values of association with the cell status. Further assessments have been
carried out on the proposed machine learning pipeline to perform classification
of the normal and cancer tissue cells. A number of machine learning classifiers
are compared for whole sequence and WHT sequence classification based on public
Whole-Genome Bisulfite Sequencing (WGBS) DNA methylation datasets. The
WHT-based method can speed up the computation time by more than one order of
magnitude compared with whole original sequence classification, while
maintaining comparable classification accuracy by the selected machine learning
classifiers. The proposed method has broad applications in expedited disease
and normal human cell classifications by the epigenome and genome datasets.


  RNA secondary structure prediction is widely used to understand RNA function.
Recently, there has been a shift away from the classical minimum free energy
(MFE) methods to partition function-based methods that account for folding
ensembles and can therefore estimate structure and base pair probabilities.
However, the classical partition function algorithm scales cubically with
sequence length, and is therefore a slow calculation for long sequences. This
slowness is even more severe than cubic-time MFE-based methods due to a larger
constant factor in runtime. Inspired by the success of our recently proposed
LinearFold algorithm that predicts the approximate MFE structure in linear
time, we design a similar linear-time heuristic algorithm, LinearPartition, to
approximate the partition function and base pairing probabilities, which is
shown to be orders of magnitude faster than Vienna RNAfold and CONTRAfold
(e.g., 2.5 days vs. 1.3 minutes on a sequence with length 32,753 nt). More
interestingly, the resulting base pairing probabilities are even better
correlated with the ground truth structures. LinearPartition also leads to a
small accuracy improvement when used for downstream structure prediction on
families with the longest length sequences (16S and 23S rRNA), as well as a
substantial improvement on long-distance base pairs (500+ nt apart).


  Controlling an evolving population is an important task in modern molecular
genetics, including directed evolution for improving the activity of molecules
and enzymes, in breeding experiments in animals and in plants, and in devising
public health strategies to suppress evolving pathogens. An optimal
intervention to direct evolution should be designed by considering its impact
over an entire stochastic evolutionary trajectory that follows. As a result, a
seemingly suboptimal intervention at a given time can be globally optimal as it
can open opportunities for desirable actions in the future. Here, we propose a
feedback control formalism to devise globally optimal artificial selection
protocol to direct the evolution of molecular phenotypes. We show that
artificial selection should be designed to counter evolutionary tradeoffs among
multi-variate phenotypes to avoid undesirable outcomes in one phenotype by
imposing selection on another. Control by artificial selection is challenged by
our ability to predict molecular evolution. We develop an information
theoretical framework and show that molecular time-scales for evolution under
natural selection can inform how to monitor a population in order to acquire
sufficient predictive information for an effective intervention with artificial
selection. Our formalism opens a new avenue for devising artificial selection
methods for directed evolution of molecular functions.


  Abnormal hemoglobins can have major consequences for tissue delivery of
oxygen. Correct diagnosis of hemoglobinopathies with altered oxygen affinity
requires a determination of hemoglobin oxygen dissociation curve (ODC), which
relates the hemoglobin oxygen saturation to the partial pressure of oxygen in
the blood. Determination of the ODC of human hemoglobin is typically carried
out under conditions in which hemoglobin is in equilibrium with O2 at each
partial pressure. However, in the human body due to the fast transit of RBCs
through tissues hemoglobin oxygen exchanges occur under non-equilibrium
conditions. We describe the determination of non-equilibrium ODC, and show that
under these conditions Hb cooperativity has two apparent components in the
Adair, Perutz, and MWC models of Hb. The first component, which we call
sequential cooperativity, accounts for ~70% of Hb cooperativity, and emerges
from the constraint of sequential binding that is shared by the three models.
The second component, which we call conformational cooperativity, accounts for
~30% of Hb cooperativity, and is due either to a conformational equilibrium
between low affinity and high affinity tetramers (as in the MWC model), or to a
conformational change from low to high affinity once two of the tetramer sites
are occupied (Perutz model).


  A simplified model of clonal plant growth is formulated, motivated by
observations of spatial structures in Posidonia oceanica meadows in the
Mediterranean Sea. Two levels of approximation are considered for the
scale-dependent feedback terms. Both take into account mortality and clonal, or
vegetative, growth as well as competition and facilitation, but the first
version is nonlocal in space while the second is local. Study of the two
versions of the model in the one-dimensional case reveals that both cases
exhibit qualitatively similar behavior (but quantitative differences) and
describe the competition between three spatially extended states, the bare soil
state, the populated state, and a pattern state, and the associated spatially
localized structures. The latter are of two types, holes in the populated state
and vegetation patches on bare ground, and are organized within distinct
snaking bifurcation diagrams. Fronts between the three extended states are
studied and a transition between pushed and pulled fronts identified. Numerical
simulations in one spatial dimension are used to determine front speeds and
confront the predictions from the marginal stability condition for pulled
fronts.