  In this work we study the stability of the equilibria reached by ecosystems
formed by a large number of species. The model we focus on are Lotka-Volterra
equations with symmetric random interactions. Our theoretical analysis,
confirmed by our numerical studies, shows that for strong and heterogeneous
interactions the system displays multiple equilibria which are all marginally
stable. This property allows us to obtain general identities between diversity
and single species responses, which generalize and saturate May's bound. By
connecting the model to systems studied in condensed matter physics, we show
that the multiple equilibria regime is analogous to a critical spin-glass
phase. This relation provides a new perspective as to why many systems in
several different fields appear to be poised at the edge of stability and also
suggests new experimental ways to probe marginal stability.


  The ability to precisely visualize the atomic geometry of the interactions
between a drug and its protein target in structural models is critical in
predicting the correct modifications in previously identified inhibitors to
create more effective next generation drugs. It is currently common practice
among medicinal chemists while attempting the above to access the information
contained in three-dimensional structures by using two-dimensional projections,
which can preclude disclosure of useful features. A more precise visualization
of the three-dimensional configuration of the atomic geometry in the models can
be achieved through the implementation of immersive virtual reality (VR). In
this work, we present a freely available software pipeline for visualising
protein structures through VR. New customer hardware, such as the HTC Vive and
the Oculus Rift utilized in this study, are available at reasonable prices.
Moreover, we have combined VR visualization with fast algorithms for simulating
intramolecular motions of protein flexibility, in an effort to further improve
structure-lead drug design by exposing molecular interactions that might be
hidden in the less informative static models.


  Actin is a key component of the cytoskeleton, which plays central roles in
cell motility, division, growth, and tensile strength. To enable this wide
range of transient mechanical processes and properties, networks of actin
filaments continuously disassemble and reassemble via active
de/re-polymerization. However, the question remains as to how
de/re-polymerization kinetics of individual actin filaments translate to
time-varying mechanics of dis/re-assembling networks. To address this question,
to ultimately elucidate the molecular mechanisms that enable cells to exhibit a
myriad of transitory mechanical properties, we couple time-resolved active
microrheology with microfluidics to measure the time-varying viscoelastic
moduli of entangled and crosslinked actin networks during chemically-triggered
network disassembly and reassembly. We also develop a corresponding
mathematical model that relates the time-evolution of filament length to the
resulting time-dependent storage moduli of evolving networks. During
disassembly, we find that the moduli exhibit two distinct exponential decays,
with experimental time constants of $\sim$169 and $\sim$47 min, which we show
arises from a phase transition from a rigid percolated network to a non-rigid
regime occurring after $\sim$90 min. During reassembly, measured moduli exhibit
linear increase with time for $\sim$90 min, after which steady-state values are
achieved. Our theoretical model shows that reassembly mechanics are dominated
by elongation kinetics, and that elongation rates are much slower than has been
recently assumed. Our measurements shed much needed light onto how
polymerization kinetics map to time-varying mechanics of actively evolving
cytoskeleton networks, and provide a powerful platform for studying other
active and dynamic systems currently under intense investigation.


  Most current methods for inferring species-level phylogenies under the
coalescent model assume that no gene flow occurs following speciation. While
some studies have examined the impact of gene flow on estimation accuracy for
certain methods, limited analytical work has been undertaken to directly assess
the potential effect of gene flow across a species phylogeny. In this paper, we
consider a three-taxon isolation-with-migration model that allows gene flow
between sister taxa for a brief period following speciation, as well as
variation in the effective population sizes across the tree. We derive the
probabilities of each of the three gene tree topologies under this model, and
show that for certain choices of the gene flow and effective population size
parameters, anomalous gene trees (i.e., gene trees that are discordant with the
species tree but that have higher probability than the gene tree concordant
with the species tree) exist. We characterize the region of parameter space
producing anomalous trees, and show that the probability of the gene tree that
is concordant with the species tree can be arbitrarily small. We then show that
the SVDQuartets method is theoretically valid under the model of gene flow
between sister taxa. We study its performance on simulated data and compare it
to two other commonly-used methods for species tree inference, ASTRAL and
MP-EST. The simulations show that ASTRAL and MP-EST can be statistically
inconsistent when gene flow is present, while SVDQuartets performs well, though
large sample sizes may be required for certain parameter choices.


  There are economic and physical limitations when applying prevention and
control strategies for urban vector borne diseases. Consequently, there are
increasing concerns and interest in designing efficient strategies and
regulations that health agencies can follow in order to reduce the imminent
impact of viruses like Dengue, Zika and Chikungunya. That includes fumigation,
abatization, reducing the hatcheries, picking up trash, information campaigns.
A basic question that arise when designing control strategies is about which
and where these ones should focus. In other words, one would like to know
whether preventing the contagion or decrease vector population, and in which
area of the city, is more efficient. In this work, we propose risk indexes
based on the idea of secondary cases from patch to patch. Thus, they take into
account human mobility and indicate which patch has more chance to be a
corridor for the spread of the disease and which is more vulnerable. They can
also indicate the neighborhood where hatchery control will reduce more the
number of potential cases. In order to illustrate the usefulness of these
indexes, we run a set of numerical simulations in a mathematical model that
takes into account the urban mobility and the differences in population density
among the areas of a city. If i is a particular neighborhood, the transmission
risk index TR_i measures the potential secondary cases caused by a host in that
neighborhood. The vector transmission risk index VTR_i measures the potential
secondary cases caused by a vector. Finally, the vulnerability risk index VR_i
measures the potential secondary cases in the neighborhood. Transmission
indexes can be used to give geographical priority to some neighborhoods when
applying prevention and control measures. On the other hand, the vulnerability
index can be useful to implement monitoring campaigns or public health
investment.


  The emergence of a predominant phenotype within a cell population is often
triggered by a rare accumulation of DNA mutations in a single cell. For
example, tumors may be initiated by a single cell in which multiple mutations
cooperate to bypass a cell's defense mechanisms. The risk of such an event is
thus determined by the extremal accumulation of mutations across tissue cells.
To address this risk, we study the statistics of the maximum mutation numbers
in a generic, but tested, model of a renewing cell population. By drawing an
analogy between the genealogy of a cell population and the theory of branching
random walks, we obtain analytical estimates for the probability of exceeding a
threshold number of mutations and determine how the statistical distribution of
maximum mutation numbers scales with age and cell population size.


  We reconsider the deterministic haploid mutation-selection equation with two
types. This is an ordinary differential equation that describes the type
distribution (forward in time) in a population of infinite size. This paper
establishes ancestral (random) structures inherent in this deterministic model.
In a first step, we obtain a representation of the deterministic equation's
solution (and, in particular, of its equilibrium) in terms of an ancestral
process called the killed ancestral selection graph. This representation allows
one to understand the bifurcations related to the error threshold phenomenon
from a genealogical point of view. Next, we characterise the ancestral type
distribution by means of the pruned lookdown ancestral selection graph and
study its properties at equilibrium. We also provide an alternative
characterisation in terms of a piecewise-deterministic Markov process.
Throughout, emphasis is on the underlying dualities as well as on explicit
results.


  Economic experiments have shown that punishment can increase public goods
game contributions over time. However, the effectiveness of punishment is
challenged by second-order freeriding and antisocial punishment. The latter
implies that non-cooperators punish cooperators, while the former implies
unwillingness to shoulder the cost of punishment. Here we extend the theory of
cooperation in the spatial public goods game by considering four competing
strategies, which are traditional cooperators and defectors, as well as
cooperators who punish defectors and defectors who punish cooperators. We show
that if the synergistic effects are high enough to sustain cooperation based on
network reciprocity alone, antisocial punishment does not deter public
cooperation. Conversely, if synergistic effects are low and punishment is
actively needed to sustain cooperation, antisocial punishment does act
detrimental, but only if the cost-to-fine ratio is low. If the costs are
relatively high, cooperation again dominates as a result of spatial pattern
formation. Counterintuitively, defectors who do not punish cooperators, and are
thus effectively second-order freeriding on antisocial punishment, form an
active layer around punishing cooperators, which protects them against
defectors that punish cooperators. A stable three-strategy phase that is
sustained by the spontaneous emergence of cyclic dominance is also possible via
the same route. The microscopic mechanism behind the reported evolutionary
outcomes can be explained by the comparison of invasion rates that determine
the stability of subsystem solutions. Our results reveal an unlikely
evolutionary escape from adverse effects of antisocial punishment, and they
provide a rationale for why second-order freeriding is not always an impediment
to the evolutionary stability of punishment.


  This paper investigates the effects of vaccination on the dynamics of
infectious disease, which is spreading in a population concurrently with
awareness. The model considers contributions to the overall awareness from a
global information campaign, direct contacts between unaware and aware
individuals, and reported cases of infection. It is assumed that there is some
time delay between individuals becoming aware and modifying their behaviour.
Vaccination is administered to newborns, as well as to aware individuals, and
it is further assumed that vaccine-induced immunity may wane with time.
Feasibility and stability of the disease-free and endemic equilibria are
studied analytically, and conditions for the Hopf bifurcation of the endemic
steady state are found in terms of system parameters and the time delay.
Analytical results are supported by numerical continuation of the Hopf
bifurcation and numerical simulations of the model to illustrate different
types of dynamical behaviour.


  To survive starvation, Bacillus subtilis forms durable spores. After
asymmetric cell division, the septum grows around the forespore in a process
called engulfment, but the mechanism of force generation is unknown. Here, we
derived a novel biophysical model for the dynamics of cell-wall remodeling
during engulfment based on a balancing of dissipative, active, and mechanical
forces. By plotting phase diagrams, we predict that sporulation is promoted by
a line tension from the attachment of the septum to the outer cell wall, as
well as by an imbalance in turgor pressures in the mother-cell and forespore
compartments. We also predict that significant mother-cell growth hinders
engulfment. Hence, relatively simple physical principles may guide this complex
biological process.


  Stochastic simulations of cyclic three-species spatial predator-prey models
are usually performed in square lattices with nearest neighbor interactions
starting from random initial conditions. In this Letter we describe the results
of off-lattice Lotka-Volterra stochastic simulations, showing that the
emergence of spiral patterns does occur for sufficiently high values of the
(conserved) total density of individuals. We also investigate the dynamics in
our simulations, finding an empirical relation characterizing the dependence of
the characteristic peak frequency and amplitude on the total density. Finally,
we study the impact of the total density on the extinction probability, showing
how a low population density may jeopardize biodiversity.


  We propose a probabilistic model for interpreting gene expression levels that
are observed through single-cell RNA sequencing. In the model, each cell has a
low-dimensional latent representation. Additional latent variables account for
technical effects that may erroneously set some observations of gene expression
levels to zero. Conditional distributions are specified by neural networks,
giving the proposed model enough flexibility to fit the data well. We use
variational inference and stochastic optimization to approximate the posterior
distribution. The inference procedure scales to over one million cells, whereas
competing algorithms do not. Even for smaller datasets, for several tasks, the
proposed procedure outperforms state-of-the-art methods like ZIFA and
ZINB-WaVE. We also extend our framework to take into account batch effects and
other confounding factors and propose a natural Bayesian hypothesis framework
for differential expression that outperforms tradition DESeq2.


  Despite being similar in structure, functioning, and size viral pathogens
enjoy very different mostly well-defined ways of life. They occupy their hosts
for a few days (influenza), for a few weeks (measles), or even lifelong (HCV),
which manifests in acute or chronic infections. The various transmission routes
(airborne, via direct contact, etc.), degrees of infectiousness (referring to
the load required for transmission), antigenic variation/immune escape and
virulence define further pathogenic lifestyles. To survive pathogens must
infect new hosts; the success determines their fitness. Infection happens with
a certain likelihood during contact of hosts, where contact can also be
mediated by vectors. Besides structural aspects of the host-contact network,
three parameters/concepts appear to be key: the contact rate and the
infectiousness during contact, which encode the mode of transmission, and third
the immunity of susceptible hosts. From here, what can be concluded about the
evolutionary strategies of viral pathogens? This is the biological question
addressed in this paper. The answer extends earlier results (Lange & Ferguson
2009, PLoS Comput Biol 5 (10): e1000536) and makes explicit connection to
another basic work on the evolution of pathogens (Grenfell et al. 2004, Science
303: 327-332). A mathematical framework is presented that models intra- and
inter-host dynamics in a minimalistic but unified fashion covering a broad
spectrum of viral pathogens, including those that cause flu-like infections,
childhood diseases, and sexually transmitted infections. These pathogens turn
out as local maxima of the fitness landscape. The models involve differential-
and integral equations, agent-based simulation, networks, and probability.


  The graph-related symmetries of a reaction network give rise to certain
special equilibria (such as complex balanced equilibria) in deterministic
models of dynamics of the reaction network. Correspondingly, in the stochastic
setting, when modeled as a continuous-time Markov chain, these symmetries give
rise to certain special stationary measures. Previous work by Anderson, Craciun
and Kurtz identified stationary distributions of a complex balanced network;
later Cappelletti and Wiuf developed the notion of complex balancing for
stochastic systems. We define and establish the relations between reaction
balanced measure, complex balanced measure, reaction vector balanced measure,
and cycle balanced measure and prove that with mild additional hypotheses, the
former two are stationary distributions. Furthermore, in spirit of earlier work
by Joshi, we give sufficient conditions under which detailed balance of the
stationary distribution of Markov chain models implies the existence of
positive detailed balance equilibria for the related deterministic reaction
network model. Finally, we provide a complete map of the implications between
balancing properties of deterministic and corresponding stochastic reaction
systems, such as complex balance, reaction balance, reaction vector balance and
cycle balance.


  For light harvesters with a reaction center complex (LH1-RC complex) of three
types, we propose an experiment to verify our analysis based upon antenna
theories that automatically include the required structural information. Our
analysis conforms to the current understanding of light-harvesting antennae in
that we can explain known properties of these complexes. We provide an
explanation for the functional roles of the notch at the light harvester, a
functional role of the polypeptide called PufX or W at the opening, a
functional role of the special pair, a reason that the cross section of the
light harvester must not be circular, a reason that the light harvester must
not be spherical, reasons for the use of dielectric bacteriochlorophylls
instead of conductors to make the light harvester, a mechanism to prevent
damage from excess sunlight, an advantage of the dimeric form, and reasons for
the modular design of nature. Based upon our analysis we provide a mechanism
for dimerization. We predict the dimeric form of light-harvesting complexes is
favoured under intense sunlight. We further comment upon the classification of
the dimeric or S-shape complexes. The S-shape complexes should not be
considered as the third type of light harvester but simply as a composite form.


  In large but finite populations, weak demographic stochasticity due to random
birth and death events can lead to population extinction. The process is
analogous to the escaping problem of trapped particles under random forces.
Methods widely used in studying such physical systems, for instance,
Wentzel-Kramers-Brillouin (WKB) and Fokker-Planck methods, can be applied to
solve similar biological problems. In this article, we comparatively analyse
applications of WKB and Fokker-Planck methods to some typical stochastic
population dynamical models, including the logistic growth, endemic SIR,
predator-prey, and competitive Lotka-Volterra models. The mean extinction time
strongly depends on the nature of the corresponding deterministic fixed
point(s). For different types of fixed points, the extinction can be driven
either by rare events or typical Gaussian fluctuations. In the former case, the
large deviation function that governs the distribution of rare events can be
well-approximated by the WKB method in the weak noise limit. In the later case,
the simpler Fokker-Planck approximation approach is also appropriate.


  The duration of the infectious period is a crucial determinant of the ability
of an infectious disease to spread. We consider an epidemic model that is
network based and non-Markovian, containing classic Kermack-McKendrick,
pairwise, message passing, and spatial models as special cases. For this model,
we prove a monotonic relationship between the variability of the infectious
period (with fixed mean) and the probability that the infection will reach any
given subset of the population by any given time. For certain families of
distributions, this result implies that epidemic severity is decreasing with
respect to the variance of the infectious period. The striking importance of
this relationship is demonstrated numerically. We then prove, with a fixed
basic reproductive ratio (R0), a monotonic relationship between the variability
of the posterior transmission probability (which is a function of the
infectious period) and the probability that the infection will reach any given
subset of the population by any given time. Thus again, even when R0 is fixed,
variability of the infectious period tends to dampen the epidemic. Numerical
results illustrate this but indicate the relationship is weaker. We then show
how our results apply to message passing, pairwise, and Kermack-McKendrick
epidemic models, even when they are not exactly consistent with the stochastic
dynamics. For Poisson contact processes, and arbitrarily distributed infectious
periods, we demonstrate how systems of delay differential equations and
ordinary differential equations can provide upper and lower bounds,
respectively, for the probability that any given individual has been infected
by any given time


  Molecular dynamics simulation is used to model the self-assembly of
polyhedral shells containing 180 trapezoidal particles that correspond to the
T=3 virus capsid. Three kinds of particle, differing only slightly in shape,
are used to account for the effect of quasi-equivalence. Bond formation between
particles is reversible and an explicit atomistic solvent is included. Under
suitable conditions the simulations are able to produce complete shells, with
the majority of unused particles remaining as monomers, and practically no
other clusters. There are also no incorrectly assembled clusters. The
simulations reveal details of intermediate structures along the growth pathway,
information that is relevant for interpreting experiment.


  We propose a targeted intervention protocol where recovery is restricted to
individuals that have the least number of infected neighbours. Our recovery
strategy is highly efficient on any kind of network, since epidemic outbreaks
are minimal when compared to the baseline scenario of spontaneous recovery. In
the case of spatially embedded networks, we find that an epidemic stays
strongly spatially confined with a characteristic length scale undergoing a
random walk. We demonstrate numerically and analytically that this dynamics
leads to an epidemic spot with a flat surface structure and a radius that grows
linearly with the spreading rate.


  Engineering the synthetic nanopores through lipid bilayer membrane to access
the interior of a cell is a long persisting challenge in biotechnology. Here,
we demonstrate the stability and dynamics of a tile-based 6-helix DNA nanotube
(DNT) embedded in POPC lipid bilayer using the analysis of 0.2 microsecond long
equilibrium MD simulation trajectories. We observe that the head groups of the
lipid molecules close to the lumen cooperatively tilt towards the hydrophilic
sugar-phosphate backbone of DNA and form a toroidal structure around the patch
of DNT protruding in the membrane. Further, we explore the effect of ionic
concentrations to the in-solution structure and stability of the lipid-DNT
complex. Transmembrane ionic current measurements for the constant electric
field MD simulation provide the I-V characteristics of the water filled DNT
lumen in lipid membrane. With increasing salt concentrations, the measured
values of transmembrane ionic conductance of the porous DNT lumen vary from 4.3
nS to 20.6 nS. Simulations of the DNTs with ssDNA and dsDNA overhangs at the
mouth of the pore show gating effect with remarkable difference in the
transmembrane ionic conductivities for open and close state nanopores.


  Discriminating between correct and incorrect substrates is a core process in
biology but how is energy apportioned between the conflicting demands of
accuracy ($\mu$), speed ($\sigma$) and total entropy production rate ($P$)?
Previous studies have focussed on biochemical networks with simple structure or
relied on simplifying kinetic assumptions. Here, we use the linear framework
for timescale separation to analytically examine steady-state probabilities
away from thermodynamic equilibrium for networks of arbitrary complexity. We
also introduce a method of scaling parameters that is inspired by Hopfield's
treatment of kinetic proofreading. Scaling allows asymptotic exploration of
high-dimensional parameter spaces. We identify in this way a broad class of
complex networks and scalings for which the quantity $\sigma\ln(\mu)/P$ remains
asymptotically finite whenever accuracy improves from equilibrium, so that
$\mu_{eq}/\mu \to 0$. Scalings exist, however, even for Hopfield's original
network, for which $\sigma\ln(\mu)/P$ is asymptotically infinite, illustrating
the parametric complexity. Outside the asymptotic regime, numerical
calculations suggest that, under more restrictive parametric assumptions,
networks satisfy the bound, $\sigma\ln(\mu/\mu_{eq})/P < 1$, and we discuss the
biological implications for discrimination by ribosomes and DNA polymerase. The
methods introduced here may be more broadly useful for analysing complex
networks that implement other forms of cellular information processing.


  Understanding the mechanisms governing population extinctions is of key
importance to many problems in ecology and evolution. Stochastic factors are
known to play a central role in extinction, but the interactions between a
population's demographic stochasticity and environmental noise remain poorly
understood. Here, we model environmental forcing as a stochastic fluctuation
between two states, one with a higher death rate than the other. We find that
in general there exists a rate of fluctuations that minimizes the mean time to
extinction, a phenomenon previously dubbed "resonant activation." We develop a
heuristic description of the phenomenon, together with a criterion for the
existence of resonant activation. Specifically the minimum extinction time
arises as a result of the system approaching a scenario wherein the severity of
rare events is balanced by the time interval between them. We discuss our
findings within the context of more general forms of environmental noise, and
suggest potential applications to evolutionary models.


  Reduction of information entropy along with ever-increasing complexity are
among the key signatures of living matter. Understanding the onset of such
behavior in early prebiotic world is essential for solving the problem of
origins of life. To elucidate this transition, we study a theoretical model of
information-storing heteropolymers capable of template-assisted ligation and
subjected to cyclic non-equilibrium driving forces. We discover that this
simple physical system undergoes a spontaneous reduction of the information
entropy due to the competition of chains for constituent monomers. This
natural-selection-like process ultimately results in the survival of a limited
subset of polymer sequences. Importantly, the number of surviving sequences
remains exponentially large, thus opening up the possibility of further
increase in complexity due to Darwinian evolution. We also propose potential
experimental implementations of our model using either biopolymers or
artificial nano-structures.


  Organisms across all domains of life regulate the size of their cells.
However, the means by which this is done is poorly understood. We study two
abstracted "molecular" models for size regulation: inhibitor dilution and
initiator accumulation. We apply the models to two settings: bacteria like
Escherichia coli, that grow fully before they set a division plane and divide
into two equally sized cells, and cells that form a bud early in the cell
division cycle, confine new growth to that bud, and divide at the connection
between that bud and the mother cell, like the budding yeast Saccharomyces
cerevisiae. In budding cells, delaying cell division until buds reach the same
size as their mother leads to very weak size control, with average cell size
and standard deviation of cell size increasing over time and saturating up to
100-fold higher than those values for cells that divide when the bud is still
substantially smaller than its mother. In budding yeast, both inhibitor
dilution or initiator accumulation models are consistent with the observation
that the daughters of diploid cells add a constant volume before they divide.
This adder behavior has also been observed in bacteria. We find that in
bacteria an inhibitor dilution model produces adder correlations that are not
robust to noise in the timing of DNA replication initiation or in the timing
from initiation of DNA replication to cell division (the C + D period). In
contrast, in bacteria an initiator accumulation model yields robust adder
correlations in the regime where noise in the timing of DNA replication
initiation is much greater than noise in the C + D period, as reported
previously [1]. In bacteria, division into two equally sized cells does not
broaden the size distribution.


  We theoretically investigate an exciton-coupled electron transfer (XCET)
process that is conversion of an exciton into a charge transfer state. This
conversion happens in an exciton transfer (XT) process, and the electron moves
away in an electron transfer(ET) process in multiple environments (baths). This
XCET process plays an essential role in the harvesting of solar energy in
biological and photovoltaic materials. We develop a practical theoretical model
to study the efficiency of XCET process that occurs either in consecutive or
concerted processes under the influence of non-Markovian baths. The role of
quantum coherence in the XT-ET system and the baths is investigated using
reduced hierarchal equations of motion (HEOM). This model includes independent
baths for each XT and ET state, in addition to a XCET bath for the conversion
process. We found that, while quantum system-bath coherence is important in the
XT and ET processes, coherence between the XT and ET processes must be
suppressed in order to realize efficient irreversible XCET process through the
weak off-diagonal interaction between the XT and ET bridge sites arises from a
XCET bath.


  Genetic regulatory circuits universally cope with different sources of noise
that limit their ability to coordinate input and output signals. In many cases,
optimal regulatory performance can be thought to correspond to configurations
of variables and parameters that maximize the mutual information between inputs
and outputs. Such optima have been well characterized in several biologically
relevant cases over the past decade. Here we use methods of statistical field
theory to calculate the statistics of the maximal mutual information (the
`capacity') achievable by tuning the input variable only in an ensemble of
regulatory motifs, such that a single controller regulates N targets. Assuming
(i) sufficiently large N, (ii) quenched random kinetic parameters, and (iii)
small noise affecting the input-output channels, we can accurately reproduce
numerical simulations both for the mean capacity and for the whole
distribution. Our results provide insight into the inherent variability in
effectiveness occurring in regulatory systems with heterogeneous kinetic
parameters.


  Should we input known genome sequence features or input sequence itself in
deep learning framework? As deep learning more popular in various applications,
researchers often come to question whether to generate features or use raw
sequences for deep learning. To answer this question, we study the prediction
accuracy of precursor miRNA prediction of feature-based deep belief network and
sequence-based convolution neural network. Tested on a variant of six-layer
convolution neural net and three-layer deep belief network, we find the raw
sequence input based convolution neural network model performs similar or
slightly better than feature based deep belief networks with best accuracy
values of 0.995 and 0.990, respectively. Both the models outperform existing
benchmarks models. The results shows us that if provided large enough data,
well devised raw sequence based deep learning models can replace feature based
deep learning models. However, construction of well behaved deep learning model
can be very challenging. In cased features can be easily extracted,
feature-based deep learning models may be a better alternative.


  The functions of RNA pseudoknots (PKs), which are minimal tertiary structural
motifs and an integral part of several ribozymes and ribonucleoprotein
complexes, are determined by their structure, stability and dynamics.
Therefore, it is important to elucidate the general principles governing their
thermodynamics/folding mechanisms. Here, we combine experiments and simulations
to examine the folding/unfolding pathways of the VPK pseudoknot, a variant of
the Mouse Mammary Tumor Virus (MMTV) PK involved in ribosomal frameshifting.
Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at
different stem/loop positions in the PK, and laser temperature-jump approaches
serve as local probes allowing us to monitor the order of assembly of VPK with
two helices with different intrinsic stabilities. The experiments and molecular
simulations show that at 50 mM KCl the dominant folding pathway populates only
the more stable partially folded hairpin. As the salt concentration is
increased a parallel folding pathway emerges, involving the less stable hairpin
structure as an alternate intermediate. Notably, the flux between the pathways
is modulated by the ionic strength. The findings support the principle that the
order of PK structure formation is determined by the relative stabilities of
the hairpins, which can be altered by sequence variations or salt
concentrations. Our study not only unambiguously demonstrates that PK folds by
parallel pathways, but also establishes that quantitative description of RNA
self-assembly requires a synergistic combination of experiments and
simulations.


  We extend a recent investigation by Meehan et al. (2019) regarding the global
stability properties of the general Kermack-McKendrick (renewal) model to the
multi-strain case. We demonstrate that the basic reproduction number of each
strain $R_{0j}$ represents a sharp threshold parameter such that when $R_{0j}
\leq 1$ for all $j$ each strain dies out and the infection-free equilibrium is
globally asymptotically stable; whereas for $R_{01} \equiv \mathrm{max}_j\,
R_{0j} > 1$ the endemic equilibrium point $\bar{P}^1$, at which only the
fittest strain (i.e. strain 1) remains in circulation, becomes globally
asymptotically stable.


  This paper presents a novel Differential Evolution algorithm for protein
folding optimization that is applied to a three-dimensional AB off-lattice
model. The proposed algorithm includes two new mechanisms. A local search is
used to improve convergence speed and to reduce the runtime complexity of the
energy calculation. For this purpose, a local movement is introduced within the
local search. The designed evolutionary algorithm has fast convergence speed
and, therefore, when it is trapped into the local optimum or a relatively good
solution is located, it is hard to locate a better similar solution. The
similar solution is different from the good solution in only a few components.
A component reinitialization method is designed to mitigate this problem. Both
the new mechanisms and the proposed algorithm were analyzed on well-known amino
acid sequences that are used frequently in the literature. Experimental results
show that the employed new mechanisms improve the efficiency of our algorithm
and that the proposed algorithm is superior to other state-of-the-art
algorithms. It obtained a hit ratio of 100% for sequences up to 18 monomers,
within a budget of $10^{11}$ solution evaluations. New best-known solutions
were obtained for most of the sequences. The existence of the symmetric
best-known solutions is also demonstrated in the paper.


  Kinesin is a molecular motor that transports cargo along microtubules. The
results of many {\it in vitro} experiments on kinesin-1 are described by
kinetic models \cite{Clancy11} in which one transition corresponds to the
forward motion and subsequent binding of the tethered motor head. We argue that
in a viscoelastic medium like the cytosol of a cell this step is not Markov and
has to be described by a non-exponential waiting time distribution. We
introduce a semi-Markov kinetic model for kinesin that takes this effect into
account. We calculate, for arbitrary waiting time distributions, the moment
generating function of the number of steps made, and determine from this the
average velocity and the diffusion constant of the motor. We illustrate our
results for the case of a waiting time distribution that is Weibull. We find
that for realistic parameter values, viscoelasticity decreases the velocity and
the diffusion constant, but increases the randomness (or Fano-factor).


  Thousands of risk variants underlying complex phenotypes (quantitative traits
and diseases) have been identified in genome-wide association studies (GWAS).
However, there are still two major challenges towards deepening our
understanding of the genetic architectures of complex phenotypes. First, the
majority of GWAS hits are in the non-coding region and their biological
interpretation is still unclear. Second, accumulating evidence from GWAS
suggests the polygenicity of complex traits, i.e., a complex trait is often
affected by many variants with small or moderate effects, whereas a large
proportion of risk variants with small effects remains unknown. The
availability of functional annotation data enables us to address the above
challenges. In this study, we propose a latent sparse mixed model (LSMM) to
integrate functional annotations with GWAS data. Not only does it increase
statistical power of the identification of risk variants, but also offers more
biological insights by detecting relevant functional annotations. To allow LSMM
scalable to millions of variants and hundreds of functional annotations, we
developed an efficient variational expectation-maximization (EM) algorithm for
model parameter estimation and statistical inference. We first conducted
comprehensive simulation studies to evaluate the performance of LSMM. Then we
applied it to analyze 30 GWAS of complex phenotypes integrated with 9 genic
category annotations and 127 tissue-specific functional annotations from the
Roadmap project. The results demonstrate that our method possesses more
statistical power over conventional methods, and can help researchers achieve
deeper understanding of genetic architecture of these complex phenotypes.


  In an earlier paper, we proved the validity of large deviations theory for
the particle approximation of quite general chemical reaction networks (CRNs).
In this paper, we extend its scope and present a more geometric insight into
the mechanism of that proof, exploiting the notion of spherical image of the
reaction polytope. This allows to view the asymptotic behavior of the vector
field describing the mass-action dynamics of chemical reactions as the result
of an interaction between the faces of this polytope in different dimensions.
We also illustrate some local aspects of the problem in a discussion of
Wentzell-Freidlin (WF) theory, together with some examples.


  The Moran process, as studied by [Lieberman, E., Hauert, C. and Nowak, M.
Evolutionary dynamics on graphs. Nature 433, pp. 312-316 (2005)], is a
stochastic process modeling the spread of genetic mutations in populations. In
this process, agents of a two-type population (i.e. mutants and residents) are
associated with the vertices of a graph. Initially, only one vertex chosen
uniformly at random is a mutant, with fitness $r > 0$, while all other
individuals are residents, with fitness $1$. In every step, an individual is
chosen with probability proportional to its fitness, and its state (mutant or
resident) is passed on to a neighbor which is chosen uniformly at random. In
this paper, we introduce and study a generalization of the model of Lieberman
et al. by assuming that different types of individuals perceive the population
through different graphs, namely $G_R(V,E_R)$ for residents and $G_M(V,E_M)$
for mutants. In this model, we study the fixation probability, i.e. the
probability that eventually only mutants remain in the population, for various
pairs of graphs. First, we transfer known results from the original
single-graph model of Lieberman et al. to our 2-graph model. Among them, we
provide a generalization of the Isothermal Theorem of Lieberman et al., that
gives sufficient conditions for a pair of graphs to have the same fixation
probability as a pair of cliques. Next, we give a 2-player strategic game view
of the process where player payoffs correspond to fixation and/or extinction
probabilities. In this setting, we attempt to identify best responses for each
player and give evidence that the clique is the most beneficial graph for both
players. Finally, we examine the possibility of efficient approximation of the
fixation probability and provide a FPRAS for the special case where the mutant
graph is complete.


  Docking is an important tool in computational drug discovery that aims to
predict the binding pose of a ligand to a target protein through a combination
of pose scoring and optimization. A scoring function that is differentiable
with respect to atom positions can be used for both scoring and gradient-based
optimization of poses for docking. Using a differentiable grid-based atomic
representation as input, we demonstrate that a scoring function learned by
training a convolutional neural network (CNN) to identify binding poses can
also be applied to pose optimization. We also show that an iteratively-trained
CNN that includes poses optimized by the first CNN in its training set performs
even better at optimizing randomly initialized poses than either the first CNN
scoring function or AutoDock Vina.


  The gene silencing effect of short interfering RNA (siRNA) is known to vary
strongly with the targeted position of the mRNA. A number of hypotheses have
been suggested to explain this phenomenon. We would like to test if this
positional effect is mainly due to the secondary structure of the mRNA at the
target site. We proposed that this structural factor can be characterized by a
single parameter called "the hydrogen bond (H-b) index", which represents the
average number of hydrogen bonds formed between nucleotides in the target
region and the rest of the mRNA. This index can be determined using a
computational approach. We tested the correlation between the H-b index and the
gene-silencing effects on three genes (Bcl-2, hTF and cyclin B1) using a
variety of siRNAs. We found that the gene-silencing effect is inversely
dependent on the H-b index, indicating that the local mRNA structure at the
targeted site is the main cause of the positional effect. Based on this
finding, we suggest that the H-b index can be a useful guideline for future
siRNA design.


  The distributions of the times to the first common ancestor t_mrca is
numerically studied for an ecological population model, the extended Moran
model. This model has a fixed population size N. The number of descendants is
drawn from a beta distribution Beta(alpha, 2-alpha) for various choices of
alpha. This includes also the classical Moran model (alpha->0) as well as the
uniform distribution (alpha=1). Using a statistical mechanics-based
large-deviation approach, the distributions can be studied over an extended
range of the support, down to probabilities like 10^{-70}, which allowed us to
study the change of the tails of the distribution when varying the value of
alpha in [0,2]. We find exponential distributions p(t_mrca)~ delta^{t_mrca} in
all cases, with systematically varying values for the base delta. Only for the
cases alpha=0 and alpha=1, analytical results are known, i.e.,
delta=\exp(-2/N^2) and delta=2/3, respectively. We recover these values,
confirming the validity of our approach. Finally, we also study the
correlations between t_mrca and the number of descendants.


  Classical rich-get-richer models have found much success in being able to
broadly reproduce the statistics and dynamics of diverse real complex systems.
These rich-get-richer models are based on classical urn models and unfold
step-by-step in discrete time. Here, we consider a natural variation acting on
a temporal continuum in the form of a partial differential equation (PDE). We
first show that the continuum version of Herbert Simon's canonical preferential
attachment model exhibits an identical size distribution. In relaxing Simon's
assumption of a linear growth mechanism, we consider the case of an arbitrary
growth kernel and find the general solution to the resultant PDE. We then
extend the PDE to multiple spatial dimensions, again determining the general
solution. Finally, we apply the model to size and wealth distributions of
firms. We obtain power law scaling for both to be concordant with simulations
as well as observational data.


  Phylogenetic networks generalise phylogenetic trees and allow for the
accurate representation of the evolutionary history of a set of present-day
species whose past includes reticulate events such as hybridisation and lateral
gene transfer. One way to obtain such a network is by starting with a (rooted)
phylogenetic tree $T$, called a base tree, and adding arcs between arcs of $T$.
The class of phylogenetic networks that can be obtained in this way is called
tree-based networks and includes the prominent classes of tree-child and
reticulation-visible networks. Initially defined for binary phylogenetic
networks, tree-based networks naturally extend to arbitrary phylogenetic
networks. In this paper, we generalise recent tree-based characterisations and
associated proximity measures for binary phylogenetic networks to arbitrary
phylogenetic networks. These characterisations are in terms of matchings in
bipartite graphs, path partitions, and antichains. Some of the generalisations
are straightforward to establish using the original approach, while others
require a very different approach. Furthermore, for an arbitrary tree-based
network $N$, we characterise the support trees of $N$, that is, the tree-based
embeddings of $N$. We use this characterisation to give an explicit formula for
the number of support trees of $N$ when $N$ is binary. This formula is written
in terms of the components of a bipartite graph.


  Non-normality can underlie pulse dynamics in many engineering contexts.
However, its role in pulses generated in biomolecular contexts is generally
unclear. Here, we address this issue using the mathematical tools of linear
algebra and systems theory on simple computational models of biomolecular
circuits. We find that non-normality is present in standard models of
feedforward loops. We used a generalized framework and pseudospectrum analysis
to identify non-normality in larger biomolecular circuit models, finding that
it correlates well with pulsing dynamics. Finally, we illustrate how these
methods can be used to provide analytical support to numerical screens for
pulsing dynamics as well as provide guidelines for design.


  Effective intervention strategies for epidemics rely on the identification of
their origin and on the robustness of the predictions made by network disease
models. We introduce a Bayesian uncertainty quantification framework to infer
model parameters for a disease spreading on a network of communities from
limited, noisy observations; the state-of-the-art computational framework
compensates for the model complexity by exploiting massively parallel computing
architectures. Using noisy, synthetic data, we show the potential of the
approach to perform robust model fitting and additionally demonstrate that we
can effectively identify the disease origin via Bayesian model selection. As
disease-related data are increasingly available, the proposed framework has
broad practical relevance for the prediction and management of epidemics.


  Signal transmission at the molecular level in many biological complexes
occurs through allosteric transitions. They describe the response a complex to
binding of ligands at sites that are spatially well separated from the binding
region. We describe the Structural Perturbation Method (SPM), based on phonon
propagation in solids, that can be used to determine the signal transmitting
allostery wiring diagram (AWD) in large but finite-sized biological complexes.
Applications to the bacterial chaperonin GroEL-GroES complex shows that the AWD
determined from structures also drive the allosteric transitions dynamically.
Both from a structural and dynamical perspective these transitions are largely
determined by formation and rupture of salt-bridges. The molecular description
of allostery in GroEL provides insights into its function, which is
quantitatively described by the Iterative Annealing Mechanism. Remarkably, in
this complex molecular machine, a deep connection is established between the
structures, reaction cycle during which GroEL undergoes a sequence of
allosteric transitions, and function in a self-consistent manner.


  Single-Particle Reconstruction (SPR) in Cryo-Electron Microscopy (cryo-EM) is
the task of estimating the 3D structure of a molecule from a set of noisy 2D
projections, taken from unknown viewing directions. Many algorithms for SPR
start from an initial reference molecule, and alternate between refining the
estimated viewing angles given the molecule, and refining the molecule given
the viewing angles. This scheme is called iterative refinement. Reliance on an
initial, user-chosen reference introduces model bias, and poor initialization
can lead to slow convergence. Furthermore, since no ground truth is available
for an unsolved molecule, it is difficult to validate the obtained results.
This creates the need for high quality ab initio models that can be quickly
obtained from experimental data with minimal priors, and which can also be used
for validation. We propose a procedure to obtain such an ab initio model
directly from raw data using Kam's autocorrelation method. Kam's method has
been known since 1980, but it leads to an underdetermined system, with missing
orthogonal matrices. Until now, this system has been solved only for special
cases, such as highly symmetric molecules or molecules for which a homologous
structure was already available. In this paper, we show that knowledge of just
two clean projections is sufficient to guarantee a unique solution to the
system. This system is solved by an optimization-based heuristic. For the first
time, we are then able to obtain a low-resolution ab initio model of an
asymmetric molecule directly from raw data, without 2D class averaging and
without tilting. Numerical results are presented on both synthetic and
experimental data.


  We review recent progress on ancestral processes related to
mutation-selection models, both in the deterministic and the stochastic
setting. We mainly rely on two concepts, namely, the killed ancestral selection
graph and the pruned lookdown ancestral selection graph. The killed ancestral
selection graph gives a representation of the type of a random individual from
a stationary population, based upon the individual's potential ancestry back
until the mutations that define the individual's type. The pruned lookdown
ancestral selection graph allows one to trace the ancestry of individuals from
a stationary distribution back into the distant past, thus leading to the
stationary distribution of ancestral types. We illustrate the results by
applying them to a prototype model for the error threshold phenomenon.


  Spindles are self-organized microtubule-based structures that segregate
chromosomes during cell division. The mass of the spindle is controlled by the
balance between microtubule turnover and nucleation. The mechanisms that
control the spatial regulation of microtubule nucleation remain poorly
understood. Previous work has found that microtubule nucleators bind to
microtubules in the spindle, but it is unclear if this binding regulates the
activity of those nucleators. Here we use a combination of experiments and
mathematical modeling to investigate this issue. We measure the concentration
of tubulin and microtubules in and around the spindle. We found a very sharp
decay in microtubules at the spindle interface, which is inconsistent with the
activity of microtubule nucleators being independent of their association with
microtubules and consistent with a model in which microtubule nucleators are
only active when bound to a microtubule. This strongly argues that the activity
of microtubule nucleators is greatly enhanced when bound to microtubules. Thus,
microtubule nucleators are both localized and activated by the microtubules
they generate.


  While a monomer of the ubiquitous hormone insulin is the biologically active
form in the human body, its hexameric assembly acts as an efficient storage
unit. However, the role of water molecules in the structure, stability and
dynamics of the insulin hexamer is poorly understood. Here we combine
experimental data with molecular dynamics simulations to investigate the shape,
structure and stability of an insulin hexamer focusing on the role of water
molecules. Both X-Ray analysis and computer simulations show that the core of
the hexamer cavity is barrel-shaped, holding, on an average, sixteen water
molecules. These encapsulated and constrained molecules impart structural
stability to the hexamer. Apart from the electrostatic interactions with Zn2+
ions, an intricate hydrogen bond network amongst cavity water and neighboring
protein residues stabilizes the hexameric association. These water molecules
solvate six glutamate residues inside the cavity decreasing electrostatic
repulsions amongst the negatively charged carboxylate groups. They also prevent
association between glutamate residues and Zn2+ ions and maintain the integrity
of the cavity. Simulations reveal that removal of these waters results in a
collapse of the cavity. Subsequent analyses also show that the hydrogen bond
network among these water molecules and protein residues that face the inner
side of the cavity is more rigid with a slower relaxation as compared to that
of the bulk solvent. Dynamics of cavity water reveal certain slow water
molecules which form the back bone of the stable hydrogen bond network. The
analysis presented here suggests a dominant role of structurally conserved
water molecules in maintaining the integrity of the hexameric assembly and
potentially modulating the dissociation of this assembly into the functional
monomeric form.


  Pluripotent embryonic stem cells are of paramount importance for biomedical
research thanks to their innate ability for self-renewal and differentiation
into all major cell lines. The fateful decision to exit or remain in the
pluripotent state is regulated by complex genetic regulatory network. Latest
advances in transcriptomics have made it possible to infer basic topologies of
pluripotency governing networks. The inferred network topologies, however, only
encode boolean information while remaining silent about the roles of dynamics
and molecular noise in gene expression. These features are widely considered
essential for functional decision making. Herein we developed a framework for
extending the boolean level networks into models accounting for individual
genetic switches and promoter architecture which allows mechanistic
interrogation of the roles of molecular noise, external signaling, and network
topology. We demonstrate the pluripotent state of the network to be a broad
attractor which is robust to variations of gene expression. Dynamics of exiting
the pluripotent state, on the other hand, is significantly influenced by the
molecular noise originating from genetic switching events which makes cells
more responsive to extracellular signals. Lastly we show that steady state
probability landscape can be significantly remodeled by global gene switching
rates alone which can be taken as a proxy for how global epigenetic
modifications exert control over stability of pluripotent states.


  The first proposed hypothesis about the mechanism of chloroquine (CQ) action
on malaria parasites is DNA intercalation hypothesis which indicates that the
site of CQ action is within the nucleus. Later on the interest of research was
shifted from nucleus to lysosome due to the report of CQ accumulation within
lysosome. The current opinions about CQ action and resistance are mainly based
on the results of more than 30-year studies on lysosome, which can be used to
explain some facts but still remains incomplete and controversial. Based on
recently published papers and our related data it is possible that the key CQ
target protein may exist in the nucleus. Development of CQ resistance is
probably mainly due to the alteration in the CQ target protein or certain
mechanism which prevents CQ from reaching its target protein in the nucleus. In
conclusion, the key site of CQ action may be in the nucleus though it has not
been well explored while CQ action in lysosome which has been well studied may
be secondarily important in CQ action and resistance.


  A fundamental problem in the fields of population genetics, evolution, and
community ecology, is the fate of a single mutant, or invader, introduced in a
finite population of wild types. For a fixed-size community of $N$ individuals,
with Markovian, zero-sum dynamics driven by stochastic birth-death events, the
mutant population eventually reaches either fixation or extinction. The
classical analysis, provided by Kimura and his coworkers, is focused on the
neutral case, [where the dynamics is only due to demographic stochasticity
(drift)], and on \emph{time-independent} selective forces
(deleterious/beneficial mutation). However, both theoretical arguments and
empirical analyses suggest that in many cases the selective forces fluctuate in
time (temporal environmental stochasticity). Here we consider a generic model
for a system with demographic noise and fluctuating selection. Our system is
characterized by the time-averaged (log)-fitness $s_0$ and zero-mean fitness
fluctuations. These fluctuations, in turn, are parameterized by their amplitude
$\gamma$ and their correlation time $\delta$. We provide asymptotic (large $N$)
formulas for the chance of fixation, the mean time to fixation and the mean
time to absorption. Our expressions interpolate correctly between the constant
selection limit $\gamma \to 0$ and the time-averaged neutral case $s_0=0$.


  With the help of a mathematical model, the metabolic process of the Krebs
cycle is studied. The autocatalytic processes resulting in both the formation
of the self-organization in the Krebs cycle and the appearance of a cyclicity
of its dynamics are determined. Some structural-functional connections creating
the synchronism of an autoperiodic functioning at the transport in the
respiratory chain and the oxidative phosphorylation are investigated. The
conditions for breaking the synchronization of processes, increasing the
multiplicity of a cyclicity, and for the appearance of chaotic modes are
analyzed. The phase-parametric diagram of a cascade of bifurcations showing the
transition to a chaotic mode by the Feigenbaum scenario is obtained. The
fractal nature of the revealed cascade of bifurcations is demonstrated. The
strange attractors formed as a result of the folding are obtained. The results
obtained give the idea of structural-functional connections, due to which the
self-organization appears in the metabolic process running in a cell. The
constructed mathematical model can be applied to the study of the toxic and
allergic effects of drugs and various substances on the metabolism of a cell.


  We investigate with computer simulations the critical radius of pores in a
lipid bilayer membrane. Ilton et al. (2016) recently showed that nucleated
pores in a homopolymer film can increase or decrease in size, depending on
whether they are larger or smaller than a critical size which scales linearly
with film thickness. Using dissipative particle dynamics, a particle-based
simulation method, we investigate the same scenario for a lipid bilayer
membrane whose structure is determined by lipid-water interactions. We simulate
a perforated membrane in which holes larger than a critical radius grow, while
holes smaller than the critical radius close, as in the experiment of Ilton et
al. (2016). By altering system parameters such as the number of particles per
lipid and the periodicity, we also describe scenarios in which pores of any
initial size can seal or even remain stable, showing a fundamental difference
in the behavior of lipid membranes from polymer films.


  Single-cell experiments show that gene expression is stochastic and bursty, a
feature that can emerge from slow switching between promoter states with
different activities. One source of long-lived promoter states is the slow
binding and unbinding kinetics of transcription factors to promoters, i.e. the
non-adiabatic binding regime. Here, we introduce a simple analytical framework,
known as a piecewise deterministic Markov process (PDMP), that accurately
describes the stochastic dynamics of gene expression in the non-adiabatic
regime. We illustrate the utility of the PDMP on a non-trivial dynamical system
by analyzing the properties of a titration-based oscillator in the
non-adiabatic limit. We first show how to transform the underlying Chemical
Master Equation into a PDMP where the slow transitions between promoter states
are stochastic, but whose rates depend upon the faster deterministic dynamics
of the transcription factors regulated by these promoters. We show that the
PDMP accurately describes the observed periods of stochastic cycles in
activator and repressor-based titration oscillators. We then generalize our
PDMP analysis to more complicated versions of titration-based oscillators to
explain how multiple binding sites lengthen the period and improve coherence.
Last, we show how noise-induced oscillation previously observed in a
titration-based oscillator arises from non-adiabatic and discrete binding
events at the promoter site.


  Chirality in shape and motility can evolve rapidly in microbes and cancer
cells. To determine how chirality affects cell fitness, we developed a model of
chiral growth in compact aggregates such as microbial colonies and solid
tumors. Our model recapitulates previous experimental findings and shows that
mutant cells can invade by increasing their chirality or switching their
handedness. The invasion results either in a takeover or stable coexistence
between the mutant and the ancestor depending on their relative chirality. For
large chiralities, the coexistence is accompanied by strong intermixing between
the cells, while spatial segregation occurs otherwise. We show that the
competition within the aggregate is mediated by bulges in regions where the
cells with different chiralities meet. The two-way coupling between aggregate
shape and natural selection is described by the chiral Kardar-Parisi-Zhang
equation coupled to the Burgers' equation with multiplicative noise. We solve
for the key features of this theory to explain the origin of selection on
chirality. Overall, our work suggests that chirality could be an important
ecological trait that mediates competition, invasion, and spatial structure in
cellular populations.


  The adaptation process of a species to a new environment is a significant
area of study in biology. As part of natural selection, adaptation is a
mutation process which improves survival skills and reproductive functions of
species. Here, we investigate this process by combining the idea of
incompetence with evolutionary game theory. In the sense of evolution,
incompetence and training can be interpreted as a special learning process.
With focus on the social side of the problem, we analyze the influence of
incompetence on behavior of species. We introduce an incompetence parameter
into a learning function in a single-population game and analyze its effect on
the outcome of the replicator dynamics. Incompetence can change the outcome of
the game and its dynamics, indicating its significance within what are
inherently imperfect natural systems.


  Objectives Influenza outbreaks have been widely studied. However, the
patterns between influenza and religious festivals remained unexplored. This
study examined the patterns of influenza and Hanukkah in Israel, and that of
influenza and Hajj in Bahrain, Egypt, Iraq, Jordan, Oman and Qatar. Method
Influenza surveillance data of these seven countries from 2009 to 2017 were
downloaded from the FluNet of the World Health Organization. Secondary data
were collected for the countries' population, and the dates of Hajj and
Hanukkah. We aggregated the weekly influenza A and B laboratory confirmations
for each country over the study period. Weekly influenza A patterns and
religious festival dates were further explored across the study period. Results
We found that influenza A peaks closely followed Hanukkah in Israel in six out
of seven years from 2010 to 2017. Aggregated influenza A peaks of the other six
Middle East countries also occurred right after Hajj every year during the
study period. Conclusions We predict that unless there is an emergence of new
influenza strain, such influenza patterns are likely to persist in future
years. Our results suggested that the optimal timing of mass influenza
vaccination should take into considerations of the dates of these religious
festivals.


  We present a dynamic coarse-graining technique that allows to simulate the
mechanical unfolding of biomolecules or molecular complexes on experimentally
relevant time scales. It is based on Markov state models (MSM), which we
construct from molecular dynamics simulations using the pulling coordinate as
an order parameter. We obtain a sequence of MSMs as a function of the
discretized pulling coordinate, and the pulling process is modeled by switching
among the MSMs according to the protocol applied to unfold the complex. This
way we cover seven orders of magnitude in pulling speed. In the region of rapid
pulling we additionally perform steered molecular dynamics simulations and find
excellent agreement between the results of the fully atomistic and the
dynamically coarse-grained simulations. Our technique allows the determination
of the rates of mechanical unfolding in a dynamical range from approximately
$10^{-8}$/ns to $1$/ns thus reaching experimentally accessible time regimes
without abandoning atomistic resolution.


  First, we revisit the stochastic Luria-Delbr\"uck model: a classic two-type
branching process which describes cell proliferation and mutation. We prove
limit theorems and exact results for the mutation times, clone sizes, and
number of mutants. Second, we extend the framework to consider mutations at
multiple sites along the genome. The number of mutants in the two-type model
characterises the mean site frequency spectrum in the multiple-site model. Our
predictions are consistent with previously published cancer genomic data.


  Bird species' migratory patterns have typically been studied through
individual observations and historical records. In recent years however, the
eBird citizen science project, which solicits observations from thousands of
bird watchers around the world, has opened the door for a data-driven approach
to understanding the large-scale geographical movements. Here, we focus on the
North American Tree Swallow (\textit{Tachycineta bicolor}) occurrence patterns
throughout the eastern United States. Migratory departure dates for this
species are widely believed by both ornithologists and casual observers to vary
substantially across years, but the reasons for this are largely unknown. In
this work, we present evidence that maximum daily temperature is a major factor
influencing Tree Swallow occurrence. Because it is generally understood that
species occurrence is a function of many complex, high-order interactions
between ecological covariates, we utilize the flexible modeling approach
offered by random forests. Making use of recent asymptotic results, we provide
formal hypothesis tests for predictive significance various covariates and also
develop and implement a permutation-based approach for formally assessing
interannual variations by treating the prediction surfaces generated by random
forests as functional data. Each of these tests suggest that maximum daily
temperature has a significant effect on migration patterns.


  The stochastic simulation of large-scale biochemical reaction networks is of
great importance for systems biology since it enables the study of inherently
stochastic biological mechanisms at the whole cell scale. Stochastic Simulation
Algorithms (SSA) allow us to simulate the dynamic behavior of complex kinetic
models, but their high computational cost makes them very slow for many
realistic size problems. We present a pilot service, named WebStoch, developed
in the context of our StochSoCs research project, allowing life scientists with
no high-performance computing expertise to perform over the internet stochastic
simulations of large-scale biological network models described in the SBML
standard format. Biomodels submitted to the service are parsed automatically
and then placed for parallel execution on distributed worker nodes. The workers
are implemented using multi-core and many-core processors, or FPGA accelerators
that can handle the simulation of thousands of stochastic repetitions of
complex biomodels, with possibly thousands of reactions and interacting
species. Using benchmark LCSE biomodels, whose workload can be scaled on
demand, we demonstrate linear speedup and more than two orders of magnitude
higher throughput than existing serial simulators.


  Bipartite networks provide an insightful representation of many systems,
ranging from mutualistic networks of species interactions to investment
networks in finance. The analysis of their topological structures has revealed
the ubiquitous presence of properties which seem to characterize many -
apparently different - systems. Nestedness, for example, has been observed in
plants-pollinator as well as in country-product trade networks. This has raised
questions about the significance of these patterns, which are often believed to
constitute a genuine signature of self-organization. Here, we review several
methods that have been developed for the analysis of such evidence. Due to the
interdisciplinary character of complex networks, tools developed in one field,
for example ecology, can greatly enrich other areas of research, such as
economy and finance, and vice versa. With this in mind, we briefly review
several entropy-based bipartite null models that have been recently proposed
and discuss their application to several real-world systems. The focus on these
models is motivated by the fact that they show three very desirable features:
analytical character, general applicability and versatility. In this respect,
entropy-based methods have been proven to perform satisfactorily both in
providing benchmarks for testing evidence-based null hypotheses and in
reconstructing unknown network configurations from partial information. On top
of that, entropy-based models have been successfully employed to analyze
ecological as well as economic systems, thus representing an ideal,
interdisciplinary tool to approach the study of bipartite complex systems.
[...]


  Recent studies reveal even the smallest genomes such as viruses evolve
through complex and stochastic processes, and the assumption of independent
alleles is not valid in most applications. Advances in sequencing technologies
produce multiple time-point whole-genome data, which enable potential
interactions between these alleles to be investigated empirically. To
investigate these interactions, we represent alleles as distributed vectors
that encode for relationships with other alleles in the course of evolution,
and apply artificial neural networks to time-sampled whole-genome datasets for
feature learning. We build this platform using methods and algorithms derived
from Natural Language Processing (NLP), and we denote it as the nucleotide
skip-gram neural network. We learn distributed vectors of alleles using the
changes in allele frequency of echovirus 11 in the presence or absence of the
disinfectant from the experimental evolution data. Results from the training
using a new open-source software TensorFlow show that the learned distributed
vectors can be clustered using Principal Component Analysis and Hierarchical
Clustering to reveal a list of non-synonymous mutations that arise on the
structural protein VP1 in connection to the candidate mutation for disinfectant
adaptation. Furthermore, this method can account for recombination rates by
setting the extent of interactions as a biological hyper-parameter, and the
results show the most realistic scenario of mid-range interactions across the
genome is most consistent with the previous studies.


  Recently, much attention has been given to understanding recombination events
along a chromosome in a variety of field. For instance, many population
genetics problems are limited by the inaccuracy of inferred evolutionary
histories of chromosomes sampled randomly from a population. This evolutionary
history differs among genomic locations as an artifact of recombination events
along a chromosome. Thus, much recent attention has been focused on identifying
these recombination points. However, many proposed methods either make
simplifying, but unrealistic, assumptions about recombination along a
chromosome, or are unable to scale to large genome-wide data like what has
become commonplace in statistical genetics. Here, we introduce a
\emph{phylogenetic derivative} to describe the relatedness of neighboring trees
along a chromosome. This phylogenetic derivative is a computationally
efficient, flexible metric that can be also be used assess the prevalence of
recombination across a chromosome. These proposed methods are tested and
perform well in analyzing both simulated data and a real mouse data set.


  This paper suggests that aging is influenced synthetically by pro-aging
factors such as ROS and anti-aging factors such as protective responses. The
anti-aging effect may be side effects of retrograde responses motivated against
adverse circumstances. ROS may be more closely correlated with metabolism
rather than aging.


  Molecular simulations produce very high-dimensional data-sets with millions
of data points. As analysis methods are often unable to cope with so many
dimensions, it is common to use dimensionality reduction and clustering methods
to reach a reduced representation of the data. Yet these methods often fail to
capture the most important features necessary for the construction of a Markov
model. Here we demonstrate the results of various dimensionality reduction
methods on two simulation data-sets, one of protein folding and another of
protein-ligand binding. The methods tested include a k-means clustering
variant, a non-linear auto encoder, principal component analysis and tICA. The
dimension-reduced data is then used to estimate the implied timescales of the
slowest process by a Markov state model analysis to assess the quality of the
projection. The projected dimensions learned from the data are visualized to
demonstrate which conformations the various methods choose to represent the
molecular process.


  Network translation has recently been used to establish steady state
properties of mass action systems by corresponding the given system to a
generalized one which is either dynamically or steady state equivalent. In this
work we further use network translation to identify network structures which
give rise to the well-studied property of absolute concentration robustness in
the corresponding mass action systems. In addition to establishing the capacity
for absolute concentration robustness, we show that network translation can
often provide a method for deriving the steady state value of the robust
species. We furthermore present a MILP algorithm for the identification of
translated chemical reaction networks that improves on previous approaches,
allowing for easier application of the theory.


  Progress in probabilistic generative models has accelerated, developing
richer models with neural architectures, implicit densities, and with scalable
algorithms for their Bayesian inference. However, there has been limited
progress in models that capture causal relationships, for example, how
individual genetic factors cause major human diseases. In this work, we focus
on two challenges in particular: How do we build richer causal models, which
can capture highly nonlinear relationships and interactions between multiple
causes? How do we adjust for latent confounders, which are variables
influencing both cause and effect and which prevent learning of causal
relationships? To address these challenges, we synthesize ideas from causality
and modern probabilistic modeling. For the first, we describe implicit causal
models, a class of causal models that leverages neural architectures with an
implicit density. For the second, we describe an implicit causal model that
adjusts for confounders by sharing strength across examples. In experiments, we
scale Bayesian inference on up to a billion genetic measurements. We achieve
state of the art accuracy for identifying causal factors: we significantly
outperform existing genetics methods by an absolute difference of 15-45.3%.


  Oxidation-associated malondialdehyde (MDA) modification of proteins can
generate immunogenic neo-epitopes that are recognized by autoantibodies. In
health, IgM antibodies to MDA-adducts are part of the natural antibody pool,
while elevated levels of IgG anti-MDA are associated with inflammatory
conditions. Yet, in human autoimmune disease IgG anti-MDA responses have not
been well characterized and their potential contribution to disease
pathogenesis is not known. Here, we investigate MDA-modifications and
anti-MDA-modified protein autoreactivity in rheumatoid arthritis (RA). While RA
is primarily associated with autoreactivity to citrullinated antigens, we also
observed increases in serum IgG anti-MDA in RA patients compared to controls.
IgG anti-MDA levels significantly correlated with disease activity by DAS28-ESR
and serum TNF-alpha, IL-6, and CRP. Mass spectrometry analysis of RA synovial
tissue identified MDA-modified proteins and revealed shared peptides between
MDA-modified and citrullinated actin and vimentin. Furthermore, anti-MDA
autoreactivity among synovial B cells was discovered when investigating
recombinant monoclonal antibodies (mAbs) cloned from single B cells. Several
clones were highly specific for MDA-modification with no cross-reactivity to
other antigen modifications. The mAbs recognized MDA-adducts in a variety of
proteins. Interestingly, the most reactive clone, originated from an
IgG1-bearing memory B cell, was encoded by germline variable genes, and showed
similarity to previously reported natural IgM. Other anti-MDA clones display
somatic hypermutations and lower reactivity. These anti-MDA antibodies had
significant in vitro functional properties and induced enhanced
osteoclastogenesis, while the natural antibody related high-reactivity clone
did not. We postulate that these may represent distinctly different facets of
anti-MDA autoreactive responses.


  Natural IgM autoantibodies have been proposed to convey protection from
autoimmune pathogenesis. Herein, we investigated the IgM responses in 396
systemic lupus erythematosus (SLE) patients, divided into subgroups based on
distinct autoantibody profiles. Depressed IgM levels were more common in SLE
than in matched population controls. Strikingly, an autoreactivity profile
defined by IgG anti-Ro/La was associated with reduced levels of specific
natural IgM anti-phosphorylcholine (PC) antigens and anti-malondialdehyde (MDA)
modified-protein, as well total IgM, while no differences were detected in SLE
patients with an autoreactivity profile defined by
anti-cardiolipin/Beta2glycoprotein-I. We also observed an association of
reduced IgM levels with the HLA-DRB1*03 allelic variant amongst SLE patients
and controls. Associations of low IgM anti-PC with cardiovascular disease were
primarily found in patients without antiphospholipid antibodies. These studies
further highlight the clinical relevance of depressed IgM. Our results suggest
that low IgM levels in SLE patients reflect immunological and genetic
differences between SLE subgroups.


  Individual cells exhibit specific proliferative responses to changes in
microenvironmental conditions. Whether such potential is constrained by the
cell density throughout the growth process is however unclear. Here, we
identify a theoretical framework that captures how the information encoded in
the initial density of cancer cell populations impacts their growth profile. By
following the growth of hundreds of populations of cancer cells, we found that
the time they need to adapt to the environment decreases as the initial cell
density increases. Moreover, the population growth rate shows a maximum at
intermediate initial densities. With the support of a mathematical model, we
show that the observed interdependence of adaptation time and growth rate is
significantly at odds both with standard logistic growth models and with the
Monod-like function that governs the dependence of the growth rate on nutrient
levels. Our results (i) uncover and quantify a previously unnoticed
heterogeneity in the growth dynamics of cancer cell populations; (ii) unveil
how population growth may be affected by single-cell adaptation times; (iii)
contribute to our understanding of the clinically-observed dependence of the
primary and metastatic tumor take rates on the initial density of implanted
cancer cells.


  The fast and affordable sequencing of large clinical and environmental
metagenomic datasets opens up new horizons in medical and biotechnological
applications. It is believed that today we have described only about 1\% of the
microorganisms on the Earth, therefore, metagenomic analysis mostly deals with
unknown species in the samples. Microbial communities in extreme environments
may contain genes with high biotechnological potential, and clinical
metagenomes, related to diseases, may uncover still unknown pathogens and
pathological mechanisms in known diseases. While the species-level
identification and description of the taxa in the samples does not seem to be
possible today, we can search for novel genes with known functions in these
samples, using numerous techniques, including artificial intelligence tools,
like the hidden Markov models (HMMs). Here we describe a simple-to-use
webserver, the MetaHMM, which is capable of homology-based automatic
model-building for the genes to be searched for, and it also finds the closest
matches in the metagenome. The webserver uses already highly successful
building blocks: it performs multiple alignment by applying Clustal Omega,
builds a hidden Markov model with HMMER components of hmmbuild and uses
hmmsearch for finding similar sequences to the specified model in the
metagenomes. The webserver is publicly available at
\url{https://metahmm.pitgroup.org}.


  Much research has focused on the deleterious effects of free-riding in public
goods games, and a variety of mechanisms that suppresses cheating behaviour.
Here we argue that under certain conditions cheating behaviour can be
beneficial to the population. In a public goods game, cheaters do not pay for
the cost of the public goods, yet they receive the benefit. Although this
free-riding harms the entire population in the long run, the success of
cheaters may aid the population when there is a common enemy that antagonizes
both cooperators and cheaters. Here we study models in which an immune system
antagonizes a cooperating pathogen. We investigate three population dynamics
models, and determine under what conditions the presence of cheaters help
defeat the immune system. The mechanism of action is that a polymorphism of
cheaters and altruists optimizes the average growth rate. Our results give
support for a possible synergy between cooperators and cheaters in ecological
public goods games.


  We model a situation in which a collection of species derive their fitnesses
via a rock-paper-scissors-type game; however, the precise payoffs are a
function of the environment. The new aspect of our model lies in adding a
feedback loop: the environment changes according to the relative fitnesses of
the species; in particular, it gives a boost to the species having small
populations. We cast our model in the setting of a differential game and we
show that for a certain setting of parameters, this dynamics cycles. Our model
is a natural one, since depletion of resources used by more frequent species
will shift the payoff matrix towards favoring less frequent ones. Since the
dynamics cycles, no species goes extinct and diversity is maintained.


  In this paper, we carry out a computational study using the spectral
decomposition of the fluctuations of a two-pathogen epidemic model around its
deterministic attractor, i.e., steady state or limit cycle, to examine the role
of partial vaccination and between-host pathogen interaction on early pathogen
replacement during seasonal epidemics of influenza and respiratory syncytial
virus.


  The structure of complex networks has been of interest in many scientific and
engineering disciplines over the decades. A number of studies in the field have
been focused on finding the common properties among different kinds of networks
such as heavy-tail degree distribution, small-worldness and modular structure
and they have tried to establish a theory of structural universality in complex
networks. However, there is no comprehensive study of network structure across
a diverse set of domains in order to explain the structural diversity we
observe in the real-world networks. In this paper, we study 986 real-world
networks of diverse domains ranging from ecological food webs to online social
networks along with 575 networks generated from four popular network models.
Our study utilizes a number of machine learning techniques such as random
forest and confusion matrix in order to show the relationships among network
domains in terms of network structure. Our results indicate that there are some
partitions of network categories in which networks are hard to distinguish
based purely on network structure. We have found that these partitions of
network categories tend to have similar underlying functions, constraints
and/or generative mechanisms of networks even though networks in the same
partition have different origins, e.g., biological processes, results of
engineering by human being, etc. This suggests that the origin of a network,
whether it's biological, technological or social, may not necessarily be a
decisive factor of the formation of similar network structure. Our findings
shed light on the possible direction along which we could uncover the hidden
principles for the structural diversity of complex networks.


  We propose and investigate a one-parameter probabilistic mixture of
one-dimensional elementary cellular automata under the guise of a model for the
dynamics of a single-species unstructured population with nonoverlapping
generations in which individuals have smaller probability of reproducing and
surviving in a crowded neighbourhood but also suffer from isolation and
dispersal. Remarkably, the first-order mean field approximation to the dynamics
of the model yields a cubic map containing terms representing both logistic and
weak Allee effects. The model has a single absorbing state devoid of
individuals, but depending on the reproduction and survival probabilities can
achieve a stable population. We determine the critical probability separating
these two phases and find that the phase transition between them is in the
directed percolation universality class of critical behaviour.


  An RNA sequence is a word over an alphabet on four elements $\{A,C,G,U\}$
called bases. RNA sequences fold into secondary structures where some bases
match one another while others remain unpaired. Pseudoknot-free secondary
structures can be represented as well-parenthesized expressions with additional
dots, where pairs of matching parentheses symbolize paired bases and dots,
unpaired bases. The two fundamental problems in RNA algorithmic are to predict
how sequences fold within some model of energy and to design sequences of bases
which will fold into targeted secondary structures. Predicting how a given RNA
sequence folds into a pseudoknot-free secondary structure is known to be
solvable in cubic time since the eighties and in truly subcubic time by a
recent result of Bringmann et al. (FOCS 2016). As a stark contrast, it is
unknown whether or not designing a given RNA secondary structure is a tractable
task; this has been raised as a challenging open question by Anne Condon (ICALP
2003). Because of its crucial importance in a number of fields such as
pharmaceutical research and biochemistry, there are dozens of heuristics and
software libraries dedicated to RNA secondary structure design. It is therefore
rather surprising that the computational complexity of this central problem in
bioinformatics has been unsettled for decades.
  In this paper we show that, in the simplest model of energy which is the
Watson-Crick model the design of secondary structures is NP-complete if one
adds natural constraints of the form: index $i$ of the sequence has to be
labeled by base $b$. This negative result suggests that the same lower bound
holds for more realistic models of energy. It is noteworthy that the additional
constraints are by no means artificial: they are provided by all the RNA design
pieces of software and they do correspond to the actual practice.


  Hepatitis B is a potentially life-threatening liver infection caused by the
hepatitis B virus (HBV). In this paper, the transmission dynamics of hepatitis
B is formulated with a mathematical model with considerations of different
classes of individuals, namely immunized, susceptible, latent,infected and
recovered class. The role of vaccination of new born babies against hepatitis B
and the treatment of both latently and actively infected individuals in
controlling the spread are factored into the model. The model in this study is
based on the standard SEIR model. The disease-free equilibrium state of the
model was established and its stability analyzed using the Routh-Hurwitz
theorem. The result of the analysis of the stability of the disease-free
equilibrium state shows that hepatitis B can totally be eradicated if effort is
made to ensure that the sum of the rate of recovery of the latent class, the
rate at which latently infected individuals become actively infected and the
rate of natural death must have a lower bound.


  We introduce a two-species exclusion model to describe the key features of
the conflict between the RNA polymerase (RNAP) motor traffic, engaged in the
transcription of a segment of DNA, concomitant with the progress of two DNA
replication forks on the same DNA segment. One of the species of particles
($P$) represents RNAP motors while the other ($R$) represents replication
forks. Motivated by the biological phenomena that this model is intended to
capture, a maximum of only two $R$ particles are allowed to enter the lattice
from two opposite ends whereas the unrestricted number of $P$ particles
constitute a totally asymmetric simple exclusion process (TASEP) in a segment
in the middle of the lattice. Consequently, the lattice consists of three
segments; the encounters of the $P$ particles with the $R$ particles are
confined within the middle segment (segment $2$) whereas only the $R$ particles
can occupy the sites in the segments $1$ and $3$. The model captures three
distinct pathways for resolving the co-directional as well as head-collision
between the $P$ and $R$ particles. Using Monte Carlo simulations and heuristic
analytical arguments that combine exact results for the TASEP with mean-field
approximations, we predict the possible outcomes of the conflict between the
traffic of RNAP motors ($P$ particles engaged in transcription) and the
replication forks ($R$ particles). The outcomes, of course, depend on the
dynamical phase of the TASEP of $P$ particles. In principle, the model can be
adapted to the experimental conditions to account for the data quantitatively.


  In living cells, chemical reactions form a complex network. Complicated
dynamics arising from such networks are the origins of biological functions. We
propose a novel mathematical method to analyze bifurcation behaviors of a
reaction system from the network structure alone. The whole network is
decomposed into subnetworks based on "buffering structures". For each
subnetwork, the bifurcation condition is studied independently, and the
parameters that can induce bifurcations and the chemicals that can exhibit
bifurcations are determined. We demonstrate our theory using hypothetical and
real networks.


  A sufficient statistic is a significant concept in statistics, which means a
probability variable that has sufficient information required for an inference
task. We investigate the roles of sufficient statistics and related quantities
in stochastic thermodynamics. Specifically, we prove that for general
continuous-time bipartite networks, the existence of a sufficient statistic
implies that an informational quantity called the sensory capacity takes the
maximum. Since the maximal sensory capacity imposes a constraint that the
energetic efficiency cannot exceed one-half, our result implies that the
existence of a sufficient statistic is inevitably accompanied by energetic
dissipation. We also show that, in a particular parameter region of linear
Langevin systems, there exists the optimal noise intensity, at which the
sensory capacity, the information-thermodynamic efficiency, and the total
entropy production are optimized at the same time. We apply our general result
to a model of sensory adaptation of E. Coli, and find that the sensory capacity
is nearly maximal with experimentally realistic parameters, suggesting that E.
Coli approximately realizes a sufficient statistic in signal transduction at
the cost of energetic dissipation.


  We investigated the potential biological impacts at Earth's surface of
stratospheric O3 depletion caused by nearby supernovae known to have occurred
about 2.5 and 8 million years ago at about 50 pc distance. New and previously
published atmospheric chemistry modeling results were combined with radiative
transfer modeling to determine changes in surface-level Solar irradiance and
biological responses. We find that UVB irradiance is increased by a factor of
1.1 to 2.8, with large variation in latitude, and seasonally at high latitude
regions. Changes in UVA and PAR (visible light) are much smaller. DNA damage
(in vitro) is increased by factors similar to UVB, while other biological
impacts (erythema, skin cancer, cataracts, marine phytoplankton photosynthesis
inhibition, and plant damage) are increased by smaller amounts. We conclude
that biological impacts due to increased UV irradiance in this SN case are not
mass-extinction level, but might be expected to contribute to changes in
species abundances; this result fits well with species turnover observed around
the Pliocene-Pleistocene boundary.


  We study unbinding of multivalent cationic ligands from oppositely charged
polymeric binding sites sparsely grafted on a flat neutral substrate. Our
molecular dynamics (MD) simulations are suggested by single-molecule studies of
protein-DNA interactions. We consider univalent salt concentrations spanning
roughly a thousandfold range, together with various concentrations of excess
ligands in solution. To reveal the ionic effects on unbinding kinetics of
spontaneous and facilitated dissociation mechanisms, we treat electrostatic
interactions both at a Debye-H\"{u}ckel (DH, or `implicit' ions, i.e., use of
an electrostatic potential with a prescribed decay length) level, as well as by
the more precise approach of considering all ionic species explicitly in the
simulations. We find that the DH approach systematically overestimates
unbinding rates, relative to the calculations where all ion pairs are present
explicitly in solution, although many aspects of the two types of calculation
are qualitatively similar. For facilitated dissociation (FD, acceleration of
unbinding by free ligands in solution) explicit ion simulations lead to
unbinding at lower free ligand concentrations. Our simulations predict a
variety of FD regimes as a function of free ligand and ion concentrations; a
particularly interesting regime is at intermediate concentrations of ligands
where non-electrostatic binding strength controls FD. We conclude that
explicit-ion electrostatic modeling is an essential component to quantitatively
tackle problems in molecular ligand dissociation, including
nucleic-acid-binding proteins.


  We introduce the SPlit-and-conQueR (SPQR) model, a coarse-grained
representation of RNA designed for structure prediction and refinement. In our
approach, the representation of a nucleotide consists of a point particle for
the phosphate group and an anisotropic particle for the nucleoside. The
interactions are, in principle, knowledge-based potentials inspired by the
ESCORE function, a base-centered scoring function. However, a special treatment
is given to base-pairing interactions and certain geometrical conformations
which are lost in a raw knowledge-base model. This results in a representation
able to describe planar canonical and non-canonical base pairs and
base-phosphate interactions and to distinguish sugar puckers and glycosidic
torsion conformations. The model is applied to the folding of several
structures, including duplexes with internal loops of non-canonical base pairs,
tetraloops, junctions and a pseudoknot. For the majority of these systems,
experimental structures are correctly predicted at the level of individual
contacts. We also propose a method for efficiently reintroducing atomistic
detail from the coarse-grained representation.