  Encounter rates link movement strategies to intra- and inter-specific
interactions, and therefore translate individual movement behavior into
higher-level ecological processes. Indeed, a large body of interacting
population theory rests on the law of mass action, which can be derived from
assumptions of Brownian motion in an enclosed container with exclusively local
perception. These assumptions imply completely uniform space use, individual
home ranges equivalent to the population range, and encounter dependent on
movement paths actually crossing. Mounting empirical evidence, however,
suggests that animals use space non-uniformly, occupy home ranges substantially
smaller than the population range, and are often capable of nonlocal
perception. Here, we explore how these empirically supported behaviors change
pairwise encounter rates. Specifically, we derive novel analytical expressions
for encounter rates under Ornstein-Uhlenbeck motion, which features non-uniform
space use and allows individual home ranges to differ from the population
range. We compare OU-based encounter predictions to those of Reflected Brownian
Motion, from which the law of mass action can be derived. For both models, we
further explore how the interplay between the scale of perception and home
range size affects encounter rates. We find that neglecting realistic movement
and perceptual behaviors can systematically bias encounter rate predictions.


  The Billera-Holmes-Vogtmann (BHV) space of weighted trees can be embedded in
Euclidean space, but the extrinsic Euclidean mean often lies outside of
treespace. Sturm showed that the intrinsic Frechet mean exists and is unique in
treespace. This Frechet mean can be approximated with an iterative algorithm,
but bounds on the convergence of the algorithm are not known, and there is no
other known polynomial algorithm for computing the Frechet mean nor even the
edges present in the mean. We give the first necessary and sufficient
conditions for an edge to be in the Frechet mean. The conditions are in the
form of inequalities on the weights of the edges. These conditions provide a
pre-processing step for finding the treespace orthant containing the Frechet
mean. This work generalizes to orthant spaces.


  The Trojan Y Chromosome Strategy (TYC) is an extremely well investigated
biological control method for controlling invasive populations with an XX-XY
sex determinism. In \cite{GP12, WP14} various dynamical properties of the
system are analyzed, including well posedness, boundedness of solutions, and
conditions for extinction or recovery. These results are derived under the
assumption of positive solutions. In the current manuscript, we show that if
the introduction rate of trojan fish is zero, under certain large data
assumptions, negative solutions are possible for the male population, which in
turn can lead to finite time blow-up in the female and male populations. A
comparable result is established for \emph{any} positive initial condition if
the introduction rate of trojan fish is large enough. Similar finite time
blow-up results are obtained in a spatial temporal TYC model that includes
diffusion. Lastly, we investigate improvements to the TYC modeling construct
that may dampen the mechanisms to the blow-up phenomenon or remove the
negativity of solutions. The results draw into suspect the reliability of
current TYC models under certain situations.


  It has long been observed that trimethylamin N-oxide (TMAO) and urea
demonstrate dramatically different properties in a protein folding process.
Even with the enormous theoretical and experimental research work of the two
osmolytes, various aspects of their underlying mechanisms still remain largely
elusive. In this paper, we propose to use the weighted persistent homology to
systematically study the osmolytes molecular aggregation and their
hydrogen-bonding network from a local topological perspective. We consider two
weighted models, i.e., localized persistent homology (LPH) and interactive
persistent homology (IPH). From the localized persistent homology models, we
have found that TMAO and urea have very different local topology. TMAO shows
local network structures. With the concentration increase, the circle elements
in these networks show a clear increase in their total numbers and a decrease
in their relative sizes. In contrast, urea shows two types of local topological
patterns, i.e., local clusters around 6 \AA~ and a few global circle elements
at around 12 \AA. From the interactive persistent homology models, it has been
found that our persistent radial distribution function (PRDF) from the
global-scale IPH has same physical properties as the traditional radial
distribution function (RDF). Moreover, PRDFs from the local-scale IPH can also
be generated and used to characterize the local interaction information. Other
than the clear difference of the first peak value of PRDFs at filtration size
4\AA, TMAO and urea also shows very different behaviors at the second peak
region from filtration size 5\AA~ to 10 \AA.


  Foraging in uncertain environments requires balancing the risks associated
with finding alternative resources against potential gains. In aridland
environments characterized by extreme variation in the amount and seasonal
timing of primary production, consumer communities must weigh the risks
associated with foraging for preferred seeds that can be cached against
fallback foods of low nutritional quality (e.g., leaves) that must be consumed
immediately. Here we explore the influence of resource-scarcity, body size, and
seasonal uncertainty on the expected foraging behaviors of caching rodents in
the northern Chihuahaun Desert by integrating these elements with a Stochastic
Dynamic Program (SDP) to determine fitness-maximizing foraging strategies. We
demonstrate that resource-limited environments promote dependence on fallback
foods, reducing the likelihood of starvation while increasing future risk
exposure. Our results point to a qualitative difference in the use of fallback
foods and the fitness benefits of caching at the threshold body size of 50 g.
Above this threshold the fitness benefits are greater for consumers with
smaller caches, affirming empirical observations of cache use among rodents in
such dynamic environments. This suggests that larger-bodied consumers with
larger caches may be less sensitive to the future uncertainties in monsoonal
onset predicted by global climate scenarios.


  The spatial localization or sequestering of motile cargo and their dispersal
within cells is an important process in a number of physiological contexts. The
morphology of the cytoskeletal network, along which active, motor-driven
intracellular transport takes place, plays a critical role in regulating such
transport phases. Here, we use a computational model to address the existence
and sensitivity of dynamic sequestering and how it depends on the parameters
governing the cytoskeletal network geometry, with a focus on filament lengths
and polarization away or toward the periphery. Our model of intracellular
transport solves for the time evolution of a probability distribution of cargo
that is transported by passive diffusion in the bulk cytoplasm and driven by
motors on explicitly rendered, polar cytoskeletal filaments with random
orientations. We show that depending on the lengths and polarizations of
filaments in the network, dynamic sequestering regions can form in different
regions of the cell. Furthermore, we find that, for certain parameters, the
residence time of cargo is non-monotonic with increasing filament length,
indicating an optimal regime for dynamic sequestration that is potentially
tunable via filament length. Our results are consistent with {\it in vivo}
observations and suggest that the ability to tunably control cargo
sequestration via cytoskeletal network regulation could provide a general
mechanism to regulate intracellular transport phases.


  The simultaneous emergence of several abrupt disease outbreaks or the
extinction of some serotypes of multi-strain diseases are fingerprints of the
interaction between pathogens spreading within the same population. Here, we
propose a general and versatile benchmark to address the unfolding of both
cooperative and competitive interacting diseases. We characterize the explosive
transitions between the disease-free and the epidemic regimes arising from the
cooperation between pathogens, and show the critical degree of cooperation
needed for the onset of such abrupt transitions. For the competing diseases, we
characterize the mutually exclusive case and derive analytically the transition
point between the full-dominance phase, in which only one pathogen propagates,
and the coexistence regime. Finally, we use this framework to analyze the
behavior of the former transition point as the competition between pathogens is
relaxed.


  We investigate the role of order/disorder transitions in alchemical
simulations of protein-ligand absolute binding free energies. We show, in the
context of a potential of mean force description, that for a benchmarking
system (the complex between the L99A mutant of T4 lysozyme and 3-iodotoluene)
and for a more challenging system relevant for medicinal applications (the
complex of the farnesoid X receptor and inhibitor 26 from a recent D3R
challenge) that order/disorder transitions can significantly hamper Hamiltonian
replica exchange sampling efficiency and slow down the rate of equilibration of
binding free energy estimates. We further show that our analytical model of
alchemical binding combined with the formalism developed by Straub et al. for
the treatment of order/disorder transitions of molecular systems can be
successfully employed to analyze the transitions and help design alchemical
schedules and soft-core functions that avoid or reduce the adverse effects of
rare binding/unbinding transitions. The results of this work pave the way for
the application of these techniques to the alchemical estimation with explicit
solvation of hydration free energies and absolute binding free energies of
systems undergoing order/disorder transitions.


  In this paper we introduce a formal method for the derivation of a predator's
functional response from a system of fast state transitions of the prey or
predator on a time scale during which the total prey and predator densities
remain constant. Such derivation permits an explicit interpretation of the
structure and parameters of the functional response in terms of individual
behaviour. The same method is also used here to derive the corresponding
numerical response of the predator as well as of the prey.


  First passage time (FPT) theory is often used to estimate timescales in
cellular and molecular biology. While the overwhelming majority of studies have
focused on the time it takes a given single Brownian searcher to reach a
target, cellular processes are instead often triggered by the arrival of the
first molecule out of many molecules. In these scenarios, the more relevant
timescale is the FPT of the first Brownian searcher to reach a target from a
large group of independent and identical Brownian searchers. Though the
searchers are identically distributed, one searcher will reach the target
before the others and will thus have the fastest FPT. This fastest FPT depends
on extremely rare events and its mean can be orders of magnitude faster than
the mean FPT of a given single searcher. In this paper, we use rigorous
probabilistic methods to study this fastest FPT. We determine the asymptotic
behavior of all the moments of this fastest FPT in the limit of many searchers
in a general class of two- and three-dimensional domains. We establish these
results by proving that the fastest searcher takes an almost direct path to the
target.


  Here we propose that the upper bound marginal stability of proteins (7.4
kcal/mol) is a universal property that includes macro-molecular complexes and
is not affected by molecular changes such as mutations and Post-Translational
Modifications. Its existence is, essentially, a consequence of the Anfinsen
thermodynamic hypothesis rather than a result of an evolutive process. This
result enables us to conjecture that neutral evolution should also be, with
respect to protein stability, a universal phenomenon.


  Despite significant progress in dissecting the genetic architecture of
complex diseases by genome-wide association studies (GWAS), the signals
identified by association analysis may not have specific pathological relevance
to diseases so that a large fraction of disease causing genetic variants is
still hidden. Association is used to measure dependence between two variables
or two sets of variables. Genome-wide association studies test association
between a disease and SNPs (or other genetic variants) across the genome.
Association analysis may detect superficial patterns between disease and
genetic variants. Association signals provide limited information on the causal
mechanism of diseases. The use of association analysis as a major analytical
platform for genetic studies of complex diseases is a key issue that hampers
discovery of the mechanism of diseases, calling into question the ability of
GWAS to identify loci underlying diseases. It is time to move beyond
association analysis toward techniques enabling the discovery of the underlying
causal genetic strctures of complex diseases. To achieve this, we propose a
concept of a genome-wide causation studies (GWCS) as an alternative to GWAS and
develop additive noise models (ANMs) for genetic causation analysis. Type I
error rates and power of the ANMs to test for causation are presented. We
conduct GWCS of schizophrenia. Both simulation and real data analysis show that
the proportion of the overlapped association and causation signals is small.
Thus, we hope that our analysis will stimulate discussion of GWAS and GWCS.


  We consider the dynamics of spatially-distributed, diffusing populations of
organisms with antagonistic interactions. These interactions are found on many
length scales, ranging from kilometer-scale animal range dynamics with
selection against hybrids to micron-scale interactions between poison-secreting
microbial populations. We find that the dynamical line tension at the interface
between antagonistic organisms suppresses survival probabilities of small
clonal clusters: the line tension introduces a critical cluster size that an
organism with a selective advantage must achieve before deterministically
spreading through the population. We calculate the survival probability as a
function of selective advantage $\delta$ and antagonistic interaction strength
$\sigma$. Unlike a simple Darwinian selective advantage, the survival
probability depends strongly on the spatial diffusion constant $D_s$ of the
strains when $\sigma>0$, with suppressed survival when both species are more
motile. Finally, we study the survival probability of a single mutant cell at
the frontier of a growing spherical cluster of cells, such as the surface of an
avascular spherical tumor. Both the inflation and curvature of the frontier
significantly enhance the survival probability by changing the critical size of
the nucleating cell cluster.


  Aggregated metapopulation lifetime statistics has been used to access
stylized facts that might help identify the underlying patch-level dynamics.
For instance, the emergence of scaling laws in the aggregated probability
distribution of patch lifetimes can be associated to critical phenomena, in
which the correlation length among system units tends to diverge. Nevertheless,
an aggregated approach is biased by patch-level variability, a fact that can
blur the interpretation of the data. Here, I propose a weakly-coupled
metapopulation model to show how patch variability can solely trigger
qualitatively different lifetime probability distribution at the aggregated
level. In a generalized approach, I obtain a two-way connection between the
variability of a certain patch property (e.g. carrying capacity, environment
condition or connectivity) and the aggregated lifetime probability
distribution. Furthermore, for a particular case, assuming that scaling laws
are observed at the aggregated-level, I speculate the heterogeneity that could
be behind it, relating the qualitative features the variability (mean, variance
and concentration) to the scaling exponents. In this perspective, the
application points to the possibility of equivalence between heterogeneous
weakly-coupled metapopulations and homogeneous ones that exhibit critical
behavior.


  We present a method for estimating parameters in stochastic models of
biochemical reaction networks by fitting steady-state distributions using
Wasserstein distances. We simulate a reaction network at different parameter
settings and train a Gaussian process to learn the Wasserstein distance between
observations and the simulator output for all parameters. We then use Bayesian
optimization to find parameters minimizing this distance based on the trained
Gaussian process. The effectiveness of our method is demonstrated on the
three-stage model of gene expression and a genetic feedback loop for which
moment-based methods are known to perform poorly. Our method is applicable to
any simulator model of stochastic reaction networks, including Brownian
Dynamics.


  Immune repertoires rely on diversity of T-cell and B-cell receptors to
protect us against foreign threats. The ability to recognize a wide variety of
pathogens is linked to the number of different clonotypes expressed by an
individual. Out of the estimated $\sim 10^{12}$ different B and T cells in
humans, how many of them express distinct receptors? We review current and past
estimates for these numbers. We point out a fundamental limitation of current
methods, which ignore the tail of small clones in the distribution of clone
sizes. We show that this tail strongly affects the total number of clones, but
it is impractical to access experimentally. We propose that combining
statistical models with mechanistic models of lymphocyte clonal dynamics offers
possible new strategies for estimating the number of clones.


  There have been several studies suggesting that protein structures solved by
NMR spectroscopy and x-ray crystallography show significant differences. To
understand the origin of these differences, we assembled a database of
high-quality protein structures solved by both methods. We also find
significant differences between NMR and crystal structures---in the
root-mean-square deviations of the C$_{\alpha}$ atomic positions, identities of
core amino acids, backbone and sidechain dihedral angles, and packing fraction
of core residues. In contrast to prior studies, we identify the physical basis
for these differences by modelling protein cores as jammed packings of
amino-acid-shaped particles. We find that we can tune the jammed packing
fraction by varying the degree of thermalization used to generate the packings.
For an athermal protocol, we find that the average jammed packing fraction is
identical to that observed in the cores of protein structures solved by x-ray
crystallography. In contrast, highly thermalized packing-generation protocols
yield jammed packing fractions that are even higher than those observed in NMR
structures. These results indicate that thermalized systems can pack more
densely than athermal systems, which suggests a physical basis for the
structural differences between protein structures solved by NMR and x-ray
crystallography.


  Motivated by the problem of classifying individuals with a disease versus
controls using functional genomic attributes as input, we encode the input as a
string of 1s (presence) or 0s (absence) of the genomic attribute across the
genome. Blocks of physical regions in the subdivided genome serve as the
feature dimensions, which takes full advantage of remaining in the
computational basis of a quantum computer. Given that a natural distance
between two binary strings is the Hamming distance and that this distance
shares properties with an inner product between qubits, we developed two
Hamming-distance-like classifiers which apply two different kinds of inner
products ("active" and "symmetric") to directly and efficiently classify a test
input into either of two training classes. To account for multiple disease and
control samples, each train class can be composed of an arbitrary number of bit
strings (i.e., number of samples) that can be compressed into one input string
per class. Thus, our circuits require the same number of qubits regardless of
the number of training samples. These algorithms, which implement a training
bisection decision plane, were simulated and implemented on IBMQX4 and IBMQX16.
The latter allowed for encoding of $64$ training features across the genome for
$2$ (disease and control) training classes.


  We present the first fixed-parameter algorithm for constructing a tree-child
phylogenetic network that displays an arbitrary number of binary input trees
and has the minimum number of reticulations among all such networks. The
algorithm uses the recently introduced framework of cherry picking sequences
and runs in $O((8k)^k \mathrm{poly}(n, t))$ time, where $n$ is the number of
leaves of every tree, $t$ is the number of trees, and $k$ is the reticulation
number of the constructed network. Moreover, we provide an efficient parallel
implementation of the algorithm and show that it can deal with up to $100$
input trees on a standard desktop computer, thereby providing a major
improvement over previous phylogenetic network construction methods.


  Apis mellifera plays a crucial role as pollinator of the majority of crops
linked to food production and thus its presence is currently fundamental to our
health and survival. The composition and configuration of the landscape in
which Apis mellifera lives will likely determine the well-being of the hives
and the pollination service that their members can provide to the crops. Here
we present a spatially explicit model that predicts the spatial distribution of
visits by Apis mellifera to crops, by simulating daily trips of honey bees, the
demographical dynamic of each hive and their honey production. This model goes
beyond existing approaches by including 1) a flower resource affected by the
feedback interaction between nectar extraction, pollination, blossoming and
repeated visits, 2) a pollinators dynamic that allows competition through short
term resource depletion, 3) a probabilistic approach of the foraging behavior,
modeling the fact that the pollinators have only partial knowledge of the
resource on their surroundings, and 4) the specific and systematic foraging
behavior and strategies of Apis mellifera at the moment of choosing foraging
sites, as opposed to those adopted by solitary and wild pollinators. With a
balance between simplicity and realism we show the importance of keeping a
minimal fraction of natural habitat in an agricultural landscape. We also
evaluate the effects of the landscape's structure on pollination, and
demonstrate that there exists an optimal size of natural habitat patches that
maximizes the pollination service for a fixed fraction of natural habitat.


  Unlike any other field, the science of morality has drawn attention from an
extraordinarily diverse set of disciplines. An interdisciplinary research
program has formed in which economists, biologists, neuroscientists,
psychologists, and even philosophers have been eager to provide answers to
puzzling questions raised by the existence of human morality. Models and
simulations, for a variety of reasons, have played various important roles in
this endeavor. Their use, however, has sometimes been deemed as useless,
trivial and inadequate. The role of models in the science of morality has been
vastly underappreciated. This omission shall be remedied here, offering a much
more positive picture on the contributions modelers made to our understanding
of morality.


  There is mounting empirical evidence that many communities of living
organisms display key features which closely resemble those of physical systems
at criticality. We here introduce a minimal model framework for the dynamics of
a community of individuals which undergoes local birth-death, immigration and
local jumps on a regular lattice. We study its properties when the system is
close to its critical point. Even if this model violates detailed balance,
within a physically relevant regime dominated by fluctuations, it is possible
to calculate analytically the probability density function of the number of
individuals living in a given volume, which captures the close-to-critical
behavior of the community across spatial scales. We find that the resulting
distribution satisfies an equation where spatial effects are encoded in
appropriate functions of space, which we calculate explicitly. The validity of
the analytical formul{\ae} is confirmed by simulations in the expected regimes.
We finally discuss how this model in the critical-like regime is in agreement
with several biodiversity patterns observed in tropical rain forests.


  The propagation of infectious diseases and its impact on individuals play a
major role in disease dynamics, and it is important to incorporate population
heterogeneity into efforts to study diseases. As a simplistic but illustrative
example, we examine interactions between urban and rural populations in the
dynamics of disease spreading. Using a compartmental framework of
susceptible--infected--susceptible ($\mathrm{SI\widetilde{S}}$) dynamics with
some level of immunity, we formulate a model that allows nonlinear reinfection.
We investigate the effects of population movement in the simplest scenario: a
case with two patches, which allows us to model movement between urban and
rural areas. To study the dynamics of the system, we compute a basic
reproduction number for each population (urban and rural). We also compute
steady states, determine the local stability of the disease-free steady state,
and identify conditions for the existence of endemic steady states. From our
analysis and computational experiments, we illustrate that population movement
plays an important role in disease dynamics. In some cases, it can be rather
beneficial, as it can enlarge the region of stability of a disease-free steady
state.


  Many enzymes appear to diffuse faster in the presence of substrate and to
drift either up or down a concentration gradient of their substrate.
Observations of these phenomena, termed enhanced enzyme diffusion (EED) and
enzyme chemotaxis, respectively, lead to a novel view of enzymes as active
matter. Enzyme chemotaxis and EED may be important in biology, and they could
have practical applications in biotechnology and nanotechnology. They also are
of considerable biophysical interest; indeed, their physical mechanisms are
still quite uncertain. This review provides an analytic summary of experimental
studies of these phenomena and of the mechanisms that have been proposed to
explain them, and offers a perspective of future directions for the field.


  The number of ribosomes in a cell is considered as limiting, and gene
expression is thus largely determined by their cellular concentration. In this
work we develop a toy model to study the trade-off between the ribosomal supply
and the demand of the translation machinery, dictated by the composition of the
transcript pool. Our equilibrium framework is useful to highlight qualitative
behaviours and new means of gene expression regulation determined by the fine
balance of this trade-off. We also speculate on the possible impact of these
mechanisms on cellular physiology.


  While much has been learnt about the impacts of specific stressors on
individual marine organisms, considerable debate exists over the nature and
impact of multiple simultaneous stressors on both individual species and marine
ecosystems. We describe a modelling tool (OSIRIS) for integrating the effects
of multiple simultaneous stressors. The model is relatively computationally
light, and demonstrated using a coarse-grained, non-spatial and simplified
representation of a temperate marine ecosystem. This version is capable of
reproducing a wide range of dynamic responses.Results indicate the degree to
which interactions are synergistic is crucial in determining sensitivity to
forcing, particularly for the higher trophic levels, which can respond
non-linearly to stronger forcing. Stronger synergistic interactions sensitize
the system to variability in forcing, and combinations of stronger forcing,
noise and synergies between effects are particularly potent. This work also
underlines the significant potential risk incurred in treating stressors on
ecosystems as individual and additive.


  Marine species reproduce and compete while being advected by turbulent flows.
It is largely unknown, both theoretically and experimentally, how population
dynamics and genetics are changed by the presence of fluid flows. Discrete
agent-based simulations in continuous space allow for accurate treatment of
advection and number fluctuations, but can be computationally expensive for
even modest organism densities. In this report, we propose an algorithm to
overcome some of these challenges. We first provide a thorough validation of
the algorithm in one and two dimensions without flow. Next, we focus on the
case of weakly compressible flows in two dimensions. This models organisms such
as phytoplankton living at a specific depth in the three-dimensional,
incompressible ocean experiencing upwelling and/or downwelling events. We show
that organisms born at sources in a two-dimensional time-independent flow
experience an increase in fixation probability.


  The primary step in search of the gene prediction is an identification of the
coding region from genomic DNA sequence. Gene structure in the case of a
eukaryotic organism is composed of promoter, intron, start codon, exons, stop
codon, etc. Splice site prediction, which separates the junction between exon
and intron, though the sequence beside. The splice sites have huge
preservation, however, the precision of the tool exhibits less than 90%. The
main objective of this work to exhibits a hybrid technique that efficiently
improves the existing gene recognition technique. Therefore to enhance the
identification of splice sites, the respective algorithm needs to be improved.
Over the last decade, the researcher paid more attention to improve the
accuracy of a predicted model in this domain. Our proposed method, SpliceCombo
involves three stages. At initial stage, which considers the principal
Component Analysis, based on the feature extracted. In the intermediate stage,
i.e.,, the second stage Case- Based Reasoning is done, i.e., feature selection.
The third stage uses support vector machine based along with polynomial kernel
function for final classification. In comparison with other methods, the
proposed SpliceCombo model outperforms other prediction models with respect to
prediction accuracies. Particularly for donor splice site the methodology
exhibits sensitivity is 97.25% accurate and specificity is 97.46% accurate. For
acceptor Splice Site the sensitivity is 96.51% and Specificity is 94.48%
correct.


  I revisit two theories of cell differentiation in multicellular organisms
published a half-century ago, Stuart Kauffman's global gene regulatory dynamics
(GGRD) model and Roy Britten's and Eric Davidson's modular gene regulatory
network (MGRN) model, in light of newer knowledge of mechanisms of gene
regulation in the metazoans (animals). The two models continue to inform
hypotheses and computational studies of differentiation of lineage-adjacent
cell types. However, their shared notion (based on bacterial regulatory
systems) of gene switches and networks built from them, have constrained
progress in understanding the dynamics and evolution of differentiation. Recent
work has described unique write-read-rewrite chromatin-based expression
encoding in eukaryotes, as well metazoan-specific processes of gene activation
and silencing in condensed-phase, enhancer-recruiting regulatory hubs,
employing disordered proteins, including transcription factors, with
context-dependent identities. These findings suggest an evolutionary scenario
in which the origination of differentiation in animals, rather than depending
exclusively on adaptive natural selection, emerged as a consequence of a type
of multicellularity in which the novel metazoan gene regulatory apparatus was
readily mobilized to amplify and exaggerate inherent cell functions of
unicellular ancestors. The plausibility of this hypothesis is illustrated by
the evolution of the developmental role of Grainyhead-like in the formation of
epithelium.


  Given that biochemical circuits can process information by using analog
computation, a question is: What can biochemical circuits compute? This paper
considers the problem of using biochemical circuits to distinguish persistent
signals from transient ones. We define a statistical detection problem over a
reaction pathway consisting of three species: an inducer, a transcription
factor (TF) and a gene promoter, where the inducer can activate the TF and an
active TF can bind to the gene promoter. We model the pathway using the
chemical master equation so the counts of bound promoters over time is a
stochastic signal. We consider the problem of using the continuous-time
stochastic signal of the counts of bound promoters to infer whether the inducer
signal is persistent or not. We use statistical detection theory to derive the
solution to this detection problem, which is to compute the log-likelihood
ratio of observing a persistent signal to a transient one. We then show, using
time-scale separation and other assumptions, that this log-likelihood ratio can
be approximately computed by using the continuous-time signals of the number of
active TF molecules and the number of bound promoters when the input is
persistent. Finally, we show that the coherent feedforward gene circuits can be
used to approximately compute this log-likelihood ratio when the inducer signal
is persistent.


  We consider a hybrid method to simulate the return time to the initial state
in a critical-case birth--death process. The expected value of this return time
is infinite, but its distribution asymptotically follows a power-law. Hence,
the simulation approach is to directly simulate the process, unless the
simulated time exceeds some threshold and if it does, draw the return time from
the tail of the power law.


  Numerical solution of the chemical master equation for stochastic reaction
networks typically suffers from the state space explosion problem due to the
curse of dimensionality and from stiffness due to multiple time scales. The
dimension of the state space equals the number of molecular species involved in
the reaction network and the size of the system of differential equations
equals the number of states in the corresponding continuous-time Markov chain,
which is usually enormously huge and often even infinite. Thus, efficient
numerical solution approaches must be able to handle huge, possibly infinite
and stiff systems of differential equations efficiently. We present an
approximate numerical integration approach that combines a dynamical state
space truncation procedure with efficient numerical integration schemes for
systems of ordinary differential equations including adaptive step size
selection based on local error estimates. The efficiency and accuracy is
demonstrated by numerical examples.


  We present a class of SEIR Markov chain models for infectious diseases
observed over discrete time in a random human population living in a closed
environment. The population changes over time through random births, deaths,
and transitions between states of the population. The SEIR models consist of
random dynamical equations for each state (S, E, I and R) involving driving
events for the process. We characterize some special types of SEIR Markov chain
models in the class including: (1) when birth and death are zero or non-zero,
and (2) when the incubation and infectious periods are constant or random. A
detailed parameter estimation applying the maximum likelihood estimation
technique and expectation maximization algorithm are presented for this study.
Numerical simulation results are given to validate the epidemic models.


  How cells regulate the number of organelles is a fundamental question in cell
biology. While decades of experimental work have uncovered four fundamental
processes that regulate organelle biogenesis, namely, de novo synthesis,
fission, fusion and decay, a comprehensive understanding of how these processes
together control organelle abundance remains elusive. Recent fluorescence
microscopy experiments allow for the counting of organelles at the single-cell
level. These measurements provide information about the cell-to-cell
variability in organelle abundance in addition to the mean level. Motivated by
such measurements, we build upon a recent study and analyze a general
stochastic model of organelle biogenesis. We compute the exact analytical
expressions for the probability distribution of organelle numbers, their mean,
and variance across a population of single cells. It is shown that different
mechanisms of organelle biogenesis lead to distinct signatures in the
distribution of organelle numbers which allows us to discriminate between these
various mechanisms. By comparing our theory against published data for
peroxisome abundance measurements in yeast, we show that a widely believed
model of peroxisome biogenesis that involves de novo synthesis, fission, and
decay is inadequate in explaining the data. Also, our theory predicts
bimodality in certain limits of the model. Overall, the framework developed
here can be harnessed to gain mechanistic insights into the process of
organelle biogenesis.


  Estrogen receptors (ERs) are a group of proteins activated by
17$\beta$-estradiol. The endocrine-disrupting chemicals (EDCs) mimic estrogen
action by bind directly to the ligand binding domain of ER. From this
perspective, ER represent a good model for identifying and assessing the health
risk of potential EDCs. This ability is best reflected by the ligand-ER binding
energy. Multilayer fragment molecular orbital (MFMO) calculations were
performed which allowed us to obtain the binding energy using a calculation
scheme that considers the molecular interactions that occur on the following
model systems: the bound and free receptor, 17$\beta$-estradiol and a water
cluster. The bound and free receptor and 17$\beta$-estradiol were surrounded by
a water shell containing the same number of molecules as the water cluster. The
structures required for MFMO calculations were obtained from molecular dynamics
simulations and cluster analysis. Attractive dispersion interactions were
observed between 17$\beta$-estradiol and the binding site hydrophobic residues.
In addition, strong electrostatic interactions were found between
17$\beta$-estradiol and the following charged/polarized residues: Glu 353, His
524 and Arg 394. The FMO2-RHF/STO-3G:MP2/6-31G(d) weighted binding energy was
of -67.2 kcal/mol. We hope that the model developed in this study can be useful
for identifying and assessing the health risk of potential EDCs.


  Measuring the oligomerization of plasma membrane proteins is rife with
biophysical and biomedical implications. This is particularly true for GPCRs, a
large family of proteins representing the targets of over one third of all FDA
approved medications. Over the last thirty years, fluorescence microscopy has
been the leading approach to address this problem. However, in spite of a large
number of studies and approaches, for most GPCRs the results have remained
highly contentious, possibly due to the large spectrum of specific methods
employed.


  We here report on non-equilibrium targeted Molecular Dynamics simulations as
tool for the estimation of protein-ligand unbinding kinetics. Correlating
simulations with experimental data from SPR kinetics measurements and X-ray
crystallography on two small molecule compound libraries bound to the
N-terminal domain of the chaperone Hsp90, we show that the mean non-equilibrium
work computed in an ensemble of trajectories of enforced ligand unbinding is a
promising predictor for ligand unbinding rates. We furthermore investigate the
molecular basis determining unbinding rates within the compound libraries. We
propose ligand conformational changes and protein-ligand nonbonded interactions
to impact on unbinding rates. Ligands may remain longer at the protein if they
exhibit strong electrostatic and/or van der Waals interactions with the target.
In the case of ligands with rigid chemical scaffold that exhibit longer
residence times however, transient electrostatic interactions with the protein
appear to facilitate unbinding. Our results imply that understanding the
unbinding pathway and the protein-ligand interactions along this path is
crucial for the prediction of small molecule ligands with defined unbinding


  Epithelial-Mesenchymal Transition (EMT), and the corresponding reverse
process, Mesenchymal-Epithelial Transition (MET), are dynamic and reversible
cellular programs orchestrated by many changes at biochemical and morphological
levels. A recent surge in identifying the molecular mechanisms underlying
EMT/MET has led to the development of various mathematical models that have
contributed to our improved understanding of dynamics at single-cell and
population levels: a) multi-stability (how many phenotypes can cells attain en
route EMT/MET?), b) reversibility/irreversibility (what time and/or
concentration of an EMT inducer marks the 'tipping point' when cells induced to
undergo EMT cannot revert?), c) symmetry in EMT/MET (do cells take the same
path while reverting as they took during the induction of EMT?), and d)
non-cell autonomous mechanisms (how does a cell undergoing EMT alter the
tendency of its neighbors to undergo EMT?). These dynamical traits may
facilitate a heterogeneous response within a cell population undergoing
EMT/MET. Here, we present a few examples of designing different mathematical
models that can contribute to decoding EMT/MET dynamics.


  We revisit the problem of the predominance of the 'weakest' species in the
context of Lotka-Volterra and May-Leonard implementations of a spatial
stochastic rock-paper-scissors model in which one of the species has its
predation probability reduced by $0 < \mathcal{P}_w < 1$. We show that,despite
the different population dynamics and spatial patterns, these two
implementations lead to qualitatively similar results for the late time values
of the relative abundances of the three species (as a function of
$\mathcal{P}_w$), as long as the simulation lattices are sufficiently large for
coexistence to prevail --- the 'weakest' species generally having an advantage
over the others (specially over its predator). However, for smaller simulation
lattices, we find that the relatively large oscillations at the initial stages
of simulations with random initial conditions may result in a significant
dependence of the probability of species survival on the lattice size and total
simulation time.


  Gene regulatory network (GRN) plays a central role in system biology and
genomics. It provides a promising way to model and study complex biological
processes. Several computational methods have been developed for the
construction and analysis of GRN. In particular, Petri net and its variants
were introduced for GRN years ago. On the other hand, Petri net theory itself
has been rapidly advanced recently. Especially noteworthy is the combination or
treatment of Petri net with the mathematical framework of category theory
(categorization), which endows Petri net with the power of abstraction and
composability. Open Petri net is a state-of-art implementation of such
"categorized" Petri nets. Applying open Petri net to GRN may potentially
facilitate the modeling of large scale GRNs. In this manuscript, we took a
shallow step towards that direction.


  Gene regulatory networks (GRNs) are complex systems in which many genes
regulate mutually to adapt the cell state to environmental conditions. In
addition to function, the GRNs possess several kinds of robustness. This
robustness means that systems do not lose their functionality when exposed to
disturbances such as mutations or noise, and is widely observed at many levels
in living systems. Both function and robustness have been acquired through
evolution. In this respect, GRNs utilized in living systems are rare among all
possible GRNs. In this study, we explored the fitness landscape of GRNs and
investigated how robustness emerged in highly-fit GRNs. We considered a toy
model of GRNs with one input gene and one output gene. The difference in the
expression level of the output gene between two input states, "on" and "off",
was considered as fitness. Thus, the determination of the fitness of a GRN was
based on how sensitively it responded to the input. We employed the
multicanonical Monte Carlo method, which can sample GRNs randomly in a wide
range of fitness levels, and classified the GRNs according to their fitness. As
a result, the following properties were found: (1) Highly-fit GRNs exhibited
bistability for intermediate input between "on" and "off". This bistability
emerges necessarily as fitness increases. (2) These highly-fit GRNs were robust
against noise because of their bistability. (3) GRNs that were robust against
mutations were not extremely rare among the highly-fit GRNs. This implies that
mutational robustness is readily acquired through the evolutionary process.


  Germinal centres are anatomically defined lymphoid organ structures that
mediate B cell affinity maturation and affect the quality of humoral immune
responses. Mathematical models based on differential equations or agent-based
simulations have been widely used to deepen our understanding of the cellular
and molecular processes characterizing these complex dynamic systems. Along
with experimental studies, these tools have provided insights into the
spatio-temporal behavior of B cells in germinal center reactions and the
mechanisms underlying their clonal selection and cellular differentiation. More
recently, mathematical models have been used to define key parameters that
influence the quality of humoral vaccine responses such as vaccine composition
and vaccination schedule. Here, we review current mathematical and
computational models of the GC response, highlight interesting model features
and discuss how they helped to improve our understanding of B cell responses
specifically in immunization settings.


  The entanglement of population dynamics, evolution, and adaptive radiation
for species competing for resources is studied. For resource harvesting, we
modify the model used in Ref. Phys. Rev. Lett. 118 048103 and introduce new
resource contest principles. We realistically implement the effects of
beneficial and deleterious mutations on the coefficients in the equations
governing the population dynamics and consider the emergence of reproductive
isolation. The proposed model is shown to be in agreement with the competitive
exclusion principle and no vacant niche principle. We establish a mechanism
that contributes to preventing the accumulation of irreversible deleterious
mutations: competition between recently diverged species/subpopulations. The
proposed model is applicable for descriptions of more complex systems. In case
of many constants in time resources, one observes very rapid specialization, a
feature not reproducible by the common model.


  We show that the basic reproduction number of an SIS patch model with
standard incidence is either strictly decreasing and strictly convex with
respect to the diffusion coefficient of infected subpopulation if the patch
reproduction numbers of at least two patches in isolation are distinct or
constant otherwise. Biologically, it means that fast diffusion of the infected
people reduces the risk of infection. This completely solves and generalizes a
conjecture by Allen et al. ({\it SIAM J Appl Math}, 67: 1283-1309, 2007).
Furthermore, a substantially improved lower bound on the multipatch
reproduction number, a generalized monotone result on the spectral bound the
Jacobian matrix of the model system at the disease-free equilibrium, and the
limiting endemic equilibrium are obtained. The approach and results can be
applied to a class of epidemic patch models where only one class of infected
compartments migrate between patches and one transmission route is involved.


  Coevolutionary arms races form between interacting populations that
constitute each other's environment and respond to mutual changes. This
inherently far-from-equilibrium process finds striking manifestations in the
adaptive immune system, where highly variable antigens and a finite repertoire
of immune receptors coevolve on comparable timescales. This unique challenge to
the immune system motivates general questions: How do ecological and
evolutionary processes interplay to shape diversity? What determine the
endurance and fate of coevolution? Here, we take the perspective of responsive
environments and develop a phenotypic model of coevolution between receptors
and antigens that both exhibit cross-reactivity (one-to-many responses). The
theory predicts that the extent of asymmetry in cross-reactivity is a key
determinant of repertoire composition: small asymmetry supports persistent
large diversity, whereas strong asymmetry yields long-lived transients of
quasispecies in both populations. The latter represents a new type of Turing
mechanism. More surprisingly, patterning in the trait space feeds back on
population dynamics: spatial resonance between the Turing modes breaks the
dynamic balance, leading to antigen extinction or unrestrained growth. Model
predictions can be tested via combined genomic and phenotypic measurements. Our
work identifies cross-reactivity as an important regulator of diversity and
coevolutionary outcome, and reveals the remarkable effect of ecological
feedback in pattern-forming systems, which drives evolution toward non-steady
states different than the Red Queen persistent cycles.


  Selectivity and sensitivity are important figures of merit in the design and
optimization of electrochemical biosensors. The efficiency of the fabricated
immunosensing surface can easily be influenced by several factors, such as
detection limit, non-specific binding, and type of sensing platform. Here, we
demonstrate the effect of macro- and micro-sized planner working electrodes (4
mm and 400 um diameter, respectively) on the electrochemical behavior and the
performance of the developed biosensor to detect human chorionic gonadotropin
(hCG). The fabricated screen-printed sensor was constructed by modifying the
carbon macro- and micro-electrodes with a linker, 1-pyrenebutyric
acid-N-hydroxysuccinimide ester (PANHS) and immobilization of anti-hCG
antibodies to detect specifically the hCG protein. The characterization of the
developed electrodes was performed by cyclic voltammetry (CV) and square wave
voltammetry (SWV). Each immunesensing system has its unique electrochemical
behavior which might be attributed to arrangement of particles on the surface.
However, the smaller surface area of the micro-electrode is found to show
higher sensitivity (1 pg/mL) compared to the macro-electrode sensor with a
lower detection limit of 100 pg/mL. The proposed assay represents a promising
approach that is highly effective for specific detection of an analyte and can
be exploited to target biomarkers for a variety of point-of-care diagnostic
applications.


  The Gaia hypothesis considers the life-environment coupled system as a single
entity that acts to regulate and maintain habitable conditions on Earth. In
this paper we discuss three mechanisms which could potentially lead to Gaia:
Selection by Survival, Sequential Selection and Entropic Hierarchy. We use the
Tangled Nature Model of co-evolution as a common framework for investigating
all three, using an extended version of the standard model to elaborate on Gaia
as an example of an entropic hierarchy. This idea, which combines sequential
selection together with a reservoir of diversity that acts as a 'memory',
implies a tendency towards growth and increasing resilience of the Gaian system
over time. We then discuss how Gaian memory could be realised in practice via
the microbial seed bank, climate refugia and lateral gene transfer and conclude
by discussing testable implications of an entropic hierarchy for the study of
Earth history and the search for life in the universe. This paper adds to the
existing taxonomy of Gaia hypotheses to suggest an "Entropic Gaia" where we
argue that increasing biomass, complexity and enhanced habitability over time
is a statistically likely feature of a co-evolving system.


  Motivation: Next generation methods of DNA sequencing produce relatively high
rate of reading errors, which interfere with de novo genome assembly of newly
sequenced organisms and particularly affect the quality of SNP detection
important for diagnostics of many hereditary diseases. There exists a number of
programs developed for correcting errors in NGS reads. Such programs utilize
various approaches and are optimized for different specific tasks, but all of
them are far from being able to correct all errors, especially in sequencing
reads that crossing by repeats and DNA from di/polyploid eukaryotic genomes.
Results: This paper describes a novel method of error correction based on
clustering of alignments of similar reads. This method is implemented in
ReadsClean program, which is designed for cleaning Illumina HiSeq sequencing
reads. We compared ReadsClean to other reads cleaning programs recognized to be
the best by several publications. Our sequence assembly tests using actual and
simulated sequencing reads show superior results achieved by ReadsClean.
Availability and implementation: ReadsClean is implemented as a standalone C
code. It is incorporated in an error correction pipeline and is freely
available to academic users at Softberry web server www.softberry.com.


  Grate and co-workers at Pacific Northwest National Laboratory recently
developed high information content triazine-based sequence-defined polymers
that are robust by not having hydrolyzable bonds and can encode structure and
functionality by having various side chains. Through molecular dynamics (MD)
simulations, the triazine polymers have been shown to form particular
sequential stacks, have stable backbone-backbone interactions through hydrogen
bonding and $\pi$-$\pi$ interactions, and conserve their \emph{cis/trans}
conformations throughout the simulation. However, we do not know the effects of
having different side chains and backbone structures on the entire conformation
and whether the \emph{cis} or \emph{trans} conformation is more stable for the
triazine polymers. For this reason, we investigate the role of non-covalent
interactions for different side chains and backbone structures on the
conformation and assembly of triazine polymers in MD simulations. Since there
is a high energy barrier associated to the \emph{cis}-\emph{trans}
isomerization, we use replica exchange molecular dynamics (REMD) to sample
various conformations of triazine hexamers. To obtain rates and intermediate
conformations, we use the recently developed concurrent adaptive sampling (CAS)
algorithm for dimer of triazine trimers. We found that the hydrogen bonding
ability of the backbone structure is critical for the triazine polymers to
self-assemble into nanorod-like structures, rather than that of the side
chains, which can help researchers design more robust materials.


  Artificial life simulations are an important tool in the study of ecological
phenomena that can be difficult to examine directly in natural environments.
Recent work has established the soundscape as an ecologically important
resource and it has been proposed that the differentiation of animal
vocalizations within a soundscape is driven by the imperative of intraspecies
communication. The experiments in this paper test that hypothesis in a
simulated soundscape in order to verify the feasibility of intraspecies
communication as a driver of acoustic niche differentiation. The impact of
intraspecies communication is found to be a significant factor in the division
of a soundscape's frequency spectrum when compared to simulations where the
need to identify signals from conspecifics does not drive the evolution of
signalling. The method of simulating the effects of interspecies interactions
on the soundscape is positioned as a tool for developing artificial life agents
that can inhabit and interact with physical ecosystems and soundscapes.


  Innovations in Next-Generation Sequencing are enabling generation of DNA
sequence data at ever faster rates and at very low cost. Large sequencing
centers typically employ hundreds of such systems. Such high-throughput and
low-cost generation of data underscores the need for commensurate acceleration
in downstream computational analysis of the sequencing data. A fundamental step
in downstream analysis is mapping of the reads to a long reference DNA
sequence, such as a reference human genome. Sequence mapping is a
compute-intensive step that accounts for more than 30% of the overall time of
the GATK workflow. BWA-MEM is one of the most widely used tools for sequence
mapping and has tens of thousands of users.
  In this work, we focus on accelerating BWA-MEM through an efficient
architecture aware implementation, while maintaining identical output. The
volume of data requires distributed computing environment, usually deploying
multicore processors. Since the application can be easily parallelized for
distributed memory systems, we focus on performance improvements on a single
socket multicore processor. BWA-MEM run time is dominated by three kernels,
collectively responsible for more than 85% of the overall compute time. We
improved the performance of these kernels by 1) improving cache reuse, 2)
simplifying the algorithms, 3) replacing small fragmented memory allocations
with a few large contiguous ones, 4) software prefetching, and 5) SIMD
utilization wherever applicable - and massive reorganization of the source code
enabling these improvements.
  As a result, we achieved nearly 2x, 183x, and 8x speedups on the three
kernels, respectively, resulting in up to 3.5x and 2.4x speedups on end-to-end
compute time over the original BWA-MEM on single thread and single socket of
Intel Xeon Skylake processor. To the best of our knowledge, this is the highest
reported speedup over BWA-MEM.


  Climate has been an important factor in shaping the distribution and
incidence of dengue cases in tropical and subtropical countries. In Costa Rica,
a tropical country with distinctive micro-climates, dengue has been endemic
since its introduction in 1993, inflicting substantial economic, social, and
public health repercussions. Using the number of dengue reported cases and
climate data from 2007-2017, we fitted a prediction model applying a
Generalized Additive Model (GAM) and Random Forest (RF) approach, which allowed
us to retrospectively predict the relative risk of dengue in five
climatological diverse municipalities around the country.


  Naively, one expects the information communicated by an enzyme downstream
within a signaling network, in which the enzyme is embedded, to grow
monotonically with the enzyme's rate of product formation. However, here we
observe that this does not necessarily hold true for allosterically regulated
enzymes, often observed in signaling networks. In particular, we find that the
mutual information between the catalytic sites of an allosterically regulated
enzyme and a receiver protein downstream in the signaling pathway depends on
the transition kinetics between the different allosterically regulated states
of the enzyme and their respective rates of product formation. Surprisingly,
our work implies that allosteric down-regulation of an enzyme's rate of product
formation may not only be used as a way to silence itself, as one would
normally expect. Rather, down-regulation may also be used to increase the
information communicated by this enzyme to a receiver protein downstream in a
signaling pathway.


  Proteins are a matter of dual nature. As a physical object, a protein
molecule is a folded chain of amino acids with multifarious biochemistry. But
it is also an instantiation along an evolutionary trajectory determined by the
function performed by the protein within a hierarchy of interwoven interaction
networks of the cell, the organism and the population. A physical theory of
proteins therefore needs to unify both aspects, the biophysical and the
evolutionary. Specifically, it should provide a model of how the DNA gene is
mapped into the functional phenotype of the protein.
  We review several physical approaches to the protein problem, focusing on a
mechanical framework which treats proteins as evolvable condensed matter:
Mutations introduce localized perturbations in the gene, which are translated
to localized perturbations in the protein matter. A natural tool to examine how
mutations shape the phenotype are Green's functions. They map the evolutionary
linkage among mutations in the gene (termed epistasis) to cooperative physical
interactions among the amino acids in the protein. We discuss how the
mechanistic view can be applied to examine basic questions of protein evolution
and design.


  We experimentally study how redistribution choices are affected by positive
and negative information regarding the behaviour of a previous participant in a
dictator game with a taking option. We use the strategy method to identify
behavioural "types", and thus distinguish "conformists" from
"counter-conformists", and unconditional choosers. Unconditional choosers make
up the greatest proportion of types (about 80%) while only about 20% of
subjects condition their responses to social information. We find that both
conformity and counter-conformity are driven by a desire to be seen as moral
(the "symbolization" dimension of moral identity). The main difference is that,
conformity is also driven by a sensitivity to what others think ("attention to
social comparison"). Unconditional giving (about 30% of players) on the other
hand is mainly driven by the centrality of moral identity to the self (the
"internalization" dimension of moral identity). Social information thus seems
to mainly affect those who care about being seen to be moral. The direction of
effect however depends on how sensitive one is to what others think.


  In this article, we derive statistical tools to analyze and distinguish the
patterns of genetic variability produced by classical and recent population
genetic models related to seed banks. In particular, we are concerned with
models described by the Kingman coalescent (K), models exhibiting so-called
weak seed banks described by a time-changed Kingman coalescent (W), models with
so-called strong seed bank described by the seed bank coalescent (S) and the
classical two-island model by Wright, described by the structured coalescent
(TI). As the presence of a (strong) seed bank should stratify a population, we
expect it to produce a signal roughly comparable to the presence of population
structure.
  We begin with a brief analysis of Wright's $F_{ST}$, which is a classical but
crude measure for population structure, followed by a derivation of the
expected site frequency spectrum (SFS) in the infinite sites model based on
'phase-type distribution calculus' as recently discussed by Hobolth et al.
(2019). Both the $F_{ST}$ and the SFS can be readily computed under various
population models, they discard statistical signal. Hence we also derive exact
likelihoods for the full sampling probabilities, which can be achieved via
recursions and a Monte Carlo scheme both in the infinite alleles and the
infinite sites model. We employ a pseudo-marginal Metropolis-Hastings algorithm
of Andrieu and Roberts (2009) to provide a method for simultaneous model
selection and parameter inference under the so-called infinitely-many sites
model, which is the most relevant in real applications.
  It turns out that this full likelihood method can reliably distinguish among
the model classes (K, W), (S) and (TI) on the basis of simulated data even from
moderate sample sizes. It is also possible to infer mutation rates, and in
particular determine whether mutation is taking place in the (strong) seed
bank.


  Significant progress in computer hardware and software have enabled molecular
dynamics (MD) simulations to model complex biological phenomena such as protein
folding. However, enabling MD simulations to access biologically relevant
timescales (e.g., beyond milliseconds) still remains challenging. These
limitations include (1) quantifying which set of states have already been
(sufficiently) sampled in an ensemble of MD runs, and (2) identifying novel
states from which simulations can be initiated to sample rare events (e.g.,
sampling folding events). With the recent success of deep learning and
artificial intelligence techniques in analyzing large datasets, we posit that
these techniques can also be used to adaptively guide MD simulations to model
such complex biological phenomena. Leveraging our recently developed
unsupervised deep learning technique to cluster protein folding trajectories
into partially folded intermediates, we build an iterative workflow that
enables our generative model to be coupled with all-atom MD simulations to fold
small protein systems on emerging high performance computing platforms. We
demonstrate our approach in folding Fs-peptide and the $\beta\beta\alpha$ (BBA)
fold, FSD-EY. Our adaptive workflow enables us to achieve an overall root-mean
squared deviation (RMSD) to the native state of 1.6$~\AA$ and 4.4~$\AA$
respectively for Fs-peptide and FSD-EY. We also highlight some emerging
challenges in the context of designing scalable workflows when data intensive
deep learning techniques are coupled to compute intensive MD simulations.


  Recent years have witnessed a trend that advanced mathematical tools, such as
algebraic topology, differential geometry, graph theory, and partial
differential equations, have been developed for describing biological
macromolecules. These tools have considerably strengthened our ability to
understand the molecular mechanism of macromolecular function, dynamics and
transport from their structures. However, currently, there is no unified
mathematical theory to analyze, describe and characterize biological
macromolecular geometry, topology, flexibility and natural mode at a variety of
scales. We introduce the de Rham-Hodge theory, a landmark of 20th Century's
mathematics, as a unified paradigm for analyzing biological macromolecular
geometry and algebraic topology, for predicting macromolecular flexibility, and
for modeling macromolecular natural modes at a variety of scales. In this
paradigm, macromolecular geometric characteristic and topological invariants
are revealed by de Rham-Hodge spectral analysis. By using the Helmholtz-Hodge
decomposition, every macromolecular vector field is split into orthogonal
divergence-free, curl-free, and harmonic components with a distinct physical
interpretation. We organize the eigenvalues and eigenvectors of the 0-form
Laplace-de Rham operator into one of the most accurate protein B-factor
predictors. By combining the 1-form Laplace-de Rham operator and the
Helfrich-type curvature energy, we predict the natural modes of both X-ray
protein structures and cryo-EM maps. We construct accurate and efficient
three-dimensional discrete exterior calculus algorithms for the aforementioned
modeling and analysis of biological macromolecules. Using extensive
experiments, we validate that the proposed paradigm is one of the most
versatile and powerful tools for biological macromolecular studies.


  In this paper, we study two-dimensional, three-dimensional monotonic and
nonmonotonic immune responses in viral infection systems. Our results show that
the viral infection systems with monotonic immune response has no bistability
appear. However, the systems with nonmonotonic immune response has bistability
appear under some conditions. For immune intensity, we got two important
thresholds, post-treatment control threshold and elite control threshold. When
immune intensity is less than post-treatment control threshold, the virus will
be rebound. The virus will be under control when immune intensity is larger
than elite control threshold. While between the two thresholds is a bistable
interval. When immune intensity is in the bistable interval, the system can
have bistability appear. Select the rate of immune cells stimulated by the
viruses as a bifurcation parameter for nonmonotonic immune responses, we prove
the system exhibits saddle-node bifurcation and transcritical bifurcation.


  The recurrent geometric network (RGN), the first end-to-end differentiable
neural architecture for protein structure prediction, is a competitive
alternative to existing models. However, the RGN's use of recurrent neural
networks (RNNs) as internal representations results in long training time and
unstable gradients. And because of its sequential nature, it is less effective
at learning global dependencies among amino acids than existing transformer
architectures. We propose the Universal Transforming Geometric Network (UTGN),
an end-to-end differentiable model that uses the encoder portion of the
Universal Transformer architecture as an alternative for internal
representations. Our experiments show that compared to RGN, UTGN achieve a
$1.7$ \si{\angstrom} improvement on the free modeling portion and a $0.7$
\si{\angstrom} improvement on the template based modeling of the CASP12
competition.


  Common respiratory viruses cause seasonal epidemics in sequential patterns.
The underlying mechanisms for this pattern have been debated for some time. For
influenza, contenders include temperature, humidity, and vitamin D levels.
While such seasonal drivers may be sufficient to explain midwinter peaking it
is unclear if other respiratory viruses peaking in spring, summer, or fall have
individual environmental drivers. Here we present a dynamic model of
interacting diseases, in an effort to explain observed seasonal patterns
without requiring multiple seasonal drivers. Our model extends the classical
SIRS-model to include multiple diseases that interact via a short-lasting,
non-specific immunity component. This temporal protection is triggered by
exposure to any epidemic disease and lasts only a couple of weeks. We show that
the inhibiting disease-interaction allows recurrent epidemic behaviour without
any seasonal driving. In the presence of a single seasonal driver our model
parsimoniously predicts epidemic patterns such as sequentially occurring
epidemic diseases. We also present a two-disease simulation reproducing
multiple features observed in time series of parainfluenza strains PIV-3 and
biennial PIV-1 epidemic patterns. Complex epidemic patterns of seasonal
diseases may be explained by non-specific innate or cloub cell (T-cell)
mediated immune responses that result in interaction between unrelated
respiratory epidemic diseases.


  Cryo-electron microscopy is a revolutionary technique that can provide 3D
density maps at near-atomic resolution. However, map validation is still an
open issue in the field. Despite several efforts from the community, it is
possible to overfit the reconstructions to noisy data. Here, inspired by modern
statistics, we develop a novel methodology that uses a small independent
particle set to validate the 3D maps. The main idea is to monitor how the map
probability evolves over the control set during the refinement. The method is
complementary to the gold-standard procedure, which generates two
reconstructions at each iteration. We low-pass filter the two reconstructions
for different frequency cutoffs, and we calculate the probability of each
filtered map given the control set. For high-quality maps, the probability
should increase as a function of the frequency cutoff and of the refinement
iteration. We also compute the similarity between the probability distributions
of the two reconstructions. As higher frequencies are added to the maps, more
dissimilar are the distributions. We optimized the BioEM software package to
perform these calculations, and tested the method on several systems, some
which were overfitted. Our results show that our method is able to discriminate
the overfitted sets from the non-overfitted ones. We conclude that having a
control particle set, not used for the refinement, is essential for
cross-validating cryo-EM maps.


  The $r$-index is a tool for compressed indexing of genomic databases for
exact pattern matching, which can be used to completely align reads that
perfectly match some part of a genome in the database or to find seeds for
reads that do not. This paper shows how to download and install the programs
ri-buildfasta and ri-align; how to call ri-buildfasta on a FASTA file to build
an $r$-index for that file; and how to query that index with ri-align.
  Availability: The source code for these programs is released under GPLv3 and
available at https://github.com/alshai/r-index .


  As genomic scale datasets motivate research on species tree inference,
simulators of the multispecies coalescent (MSC) process are essential for the
testing and evaluation of new inference methods. However, the simulators
themselves must be tested to ensure they give valid samples from the coalescent
process. In this work we develop several statistical tools using summary
statistics to evaluate the fit of a simulated gene tree sample to the MSC
model. Using these tests on samples from four published simulators, we uncover
flaws in several. The tests are implemented as an R package, so that both
developers and users will be able to easily check proper performance of future
simulators.


  There are currently plenty of programs available for mapping short sequences
(reads) to a genome. Most of them, however, including such popular and actively
developed programs as Bowtie, BWA, TopHat and many others, are based on
Burrows-Wheeler Transform (BWT) algorithm. This approach is very effective for
mapping high-homology reads, but runs into problems when mapping reads with
high level of errors or SNP. Also it has problems with mapping RNASeq spliced
reads (such as reads that aligning with gaps corresponding intron sequences),
the kind that is essential for finding introns and alternative splicing gene
isoforms. Meanwhile, finding intron positions is the most important task for
determining the gene structure, and especially alternatively spliced variants
of genes. In this paper, we propose a new algorithm that involves hashing
reference genome. ReadsMap program, implementing such algorithm, demonstrate
very high-accuracy mapping of large number of short reads to one or more
genomic contigs. It is achieved mostly by better alignment of very short parts
of reads separated by long introns with accounting information from mapping
other reads containing the same intron inserted between bigger blocks.
Availability and implementation: ReadsMap is implemented in C. It is
incorporated in Fgenesh++ gene identification pipeline and is freely available
to academic users at Softberry web server www.softberry.com.


  Oxidative stress occurs when reactive oxygen species (ROS) exceed antioxidant
defences, which can have deleterious effects on cell function, health and
survival. Therefore, organisms are expected to finely regulate pro-oxidant and
antioxidant processes. ROS are mainly produced through aerobic metabolism and
vary in response to changes in energetic requirements, whereas antioxidants may
be enhanced, depleted or show no changes in response to changes in ROS levels.
We investigated the repeatability, within-individual variation and correlation
across different conditions of two plasmatic markers of the oxidative balance
in 1108 samples from 635 free-living adult collared flycatchers (Ficedula
albicollis). We sought to manipulate energy constraints by increasing wing load
in 2012 and 2013 and by providing additional food in 2014. We then tested the
relative importance of within- and between-individual variation on reactive
oxygen metabolites (ROMs), a marker of lipid and protein peroxidation, and on
non-enzymatic antioxidant defences (OXY test). We also investigated whether the
experimental treatments modified the correlation between markers. Antioxidant
defences were repeatable (range of repeatability estimates = 0.128--0.581),
whereas ROMs were not (0--0.061). Antioxidants varied neither between
incubation and nestling feeding nor between sexes. ROMs increased from
incubation to nestling feeding in females and were higher in females than
males. Antioxidant defences and ROM concentration were globally positively
correlated, but the correlation varied between experimental conditions and
between years. Hence, the management of oxidative balance in wild animals
appears flexible under variable environmental conditions, an observation which
should be confirmed over a wider range of markers.


  Galled trees are studied as a recombination model in population genetics.
This class of phylogenetic networks is generalized into tree-child, galled and
reticulation-visible network classes by relaxing a structural condition imposed
on galled trees. We count tree-child networks through enumerating their
component graphs. Explicit counting formulas are also given for galled trees
through their relationship to ordered trees, phylogenetic networks with few
reticulations and phylogenetic networks in which the child of each reticulation
is a leaf.


  The processes that determine the establishment of the complex morphology of
neurons during development are still poorly understood. We present experiments
that use live imaging to examine the role of vesicle transport and propose a
lattice-based model that shows symmetry breaking features similar to a neuron
during its polarization. In a otherwise symmetric situation our model predicts
that a difference in neurite length increases the growth potential of the
longer neurite indicating that vesicle transport can be regarded as a major
factor in neurite growth.


  Predicting the three-dimensional (3D) structure of a protein from its primary
sequence of amino acids is known as the protein folding (PF) problem. Due to
the central role of proteins' 3D structures in chemistry, biology and medicine
applications (e.g., in drug discovery) this subject has been intensively
studied for over half a century. Although classical algorithms provide
practical solutions, sampling the conformation space of small proteins, they
cannot tackle the intrinsic NP-hard complexity of the problem, even reduced to
its simplest Hydrophobic-Polar model. While fault-tolerant quantum computers
are still beyond reach for state-of-the-art quantum technologies, there is
evidence that quantum algorithms can be successfully used on Noisy
Intermediate-Scale Quantum (NISQ) computers to accelerate energy optimization
in frustrated systems. In this work, we present a model Hamiltonian with
$\mathcal{O}(N^4)$ scaling and a corresponding quantum variational algorithm
for the folding of a polymer chain with $N$ monomers on a tetrahedral lattice.
The model reflects many physico-chemical properties of the protein, reducing
the gap between coarse-grained representations and mere lattice models. We use
a robust and versatile optimisation scheme, bringing together variational
quantum algorithms specifically adapted to classical cost functions and
evolutionary strategies (genetic algorithms), to simulate the folding of the 10
amino acid Angiotensin peptide on 22 qubits. The same method is also
successfully applied to the study of the folding of a 7 amino acid neuropeptide
using 9 qubits on an IBM Q 20-qubit quantum computer. Bringing together recent
advances in building gate-based quantum computers with noise-tolerant hybrid
quantum-classical algorithms, this work paves the way towards accessible and
relevant scientific experiments on real quantum processors.


  Dengue is a serious health problem in the Philippines. In 2016, the
Department of Health of the country launched a dengue vaccination campaign
using Dengvaxia. However, the campaign was mired with controversy and the use
of Dengvaxia was banned in the country. This study proposes a mathematical
model that represents the dynamics of the transmission of dengue with its four
strains. Considering that the Dengvaxia vaccine was intended to be given only
to people aging from 9 to 45 years old, the human population is divided into
two age groups: from 9-45 years old and the rest of the population. Using this
model and optimal control theory we simulate what could have been the effect of
Dengvaxia in the number of dengue cases and then compare this with the other
usual dengue intervention strategies. Results show that the best implementation
of the usual strategies is better than that of Dengvaxia.


  Lately, experimental research on microbial symbioses based on nutrient
exchange and interdependence has yielded a number of interesting findings,
however an in-depth mathematical description of the exact underlying dynamics
of such symbiotic associations is still missing. Here, we derive and analyse a
mechanistic mathematical model of such a relationship in a continuous chemostat
culture based on five coupled differential equations. The influence of the
biological traits of the involved organisms on the position and stability of
the equilibrium states of the system is examined. We also demonstrate how
manipulating the external metabolite concentrations of the system can shift the
species interaction on a continuous spectrum ranging from mutualism over
commensalism to parasitism.


  The complexity of an ecological community can be distilled into a network,
where diverse interactions connect species in a web of dependencies. Species
interact not only with each other but indirectly through environmental effects,
however the role of these ecosystem engineers has not yet been considered in
models of ecological networks. Here we explore the dynamics of ecosystem
assembly, where the colonization and extinction of species within a community
depends on the constraints imposed by trophic, service, and engineering
dependencies. We show that our assembly model reproduces many key features of
ecological systems, such as the role of generalists during assembly, realistic
maximum trophic levels, and increased nestedness with higher frequencies of
mutualisms. We find that ecosystem engineering has large and nonlinear effects
on extinction rates, facilitating robustness by creating niche space, but at
the same time increasing the magnitude of extinction cascades. We emphasize the
importance of redundancies in engineered effects and show that such redundancy
lowers the barriers to colonization, promoting community diversity. Together,
our results suggest that ecological engineers may enhance community diversity
while increasing persistence by facilitating colonization and limiting
competitive exclusion.


  We construct an epidemic model for the transmission of dengue fever with an
early-life stage in the vector dynamics and age-structure within hosts. The
early-life stage of the vector is modeled via a general function that supports
multiple vector densities. The {\it basic reproductive number} and {\it vector
demographic threshold} are computed to study the local and global stability of
the infection-free state. A numerical framework is implemented and simulations
are performed.


  Coexistence of plants depends on their competition for common resources and
indirect interactions mediated by shared exploiters or mutualists. These
interactions are driven either by changes in animal abundance (density-mediated
interactions, e.g., apparent competition), or by changes in animal preferences
for plants (behaviorally-mediated interactions). This article studies effects
of behaviorally-mediated interactions on two plant population dynamics and
animal preference dynamics when animal densities are fixed. Animals can be
either adaptive exploiters or adaptive mutualists (e.g., herbivores or
pollinators) that maximize their fitness. Analysis of the model shows that
adaptive animal preferences for plants can lead to multiple outcomes of plant
coexistence with different levels of specialization or generalism for the
mediator animal species. In particular, exploiter generalism promotes plant
coexistence even when inter-specific competition is too strong to make plant
coexistence possible without exploiters, and mutualist specialization promotes
plant coexistence at alternative stable states when plant inter-specific
competition is weak. Introducing a new concept of generalized isoclines allows
us to fully analyze the model with respect to the strength of competitive
interactions between plants (weak or strong), and the type of interaction
between plants and animals (exploitation or mutualism).
  Keywords: behaviorally-mediated interactions, competition for preference,
differential inclusion, generalized isocline, switching, sliding and repelling
regimes.


  Many physical and natural systems, including the population of species,
evolve in habitats with spatial stochastic variations of the individuals'
motility. We study here the effect of those fluctuations on invasion and
genetic loss. A Langevin equation for the \textit{position} and \textit{border}
of the invasion front is obtained. A striking result is that small/large
fluctuations of diffusivity suppress/intensify genetic loss. Our findings
reveal the potential role of environmental fluctuations as a regulating factor
for genetic loss and provide a simple explanation for the regional differences
in the intensity of genetic drift observed during the final stages of human
evolution and in tumor mutational landscapes.


  In microbial communities, each species often has multiple, distinct
phenotypes, but studies of ecological stability have largely ignored this
subpopulation structure. Here, we show that such implicit averaging over
phenotypes leads to incorrect linear stability results. We then analyze the
effect of phenotypic switching in detail in an asymptotic limit and partly
overturn classical stability paradigms: abundant phenotypic variation is
linearly destabilizing but, surprisingly, a rare phenotype such as bacterial
persisters has a stabilizing effect. Finally, we extend these results by
showing how phenotypic variation modifies the stability of the system to large
perturbations such as antibiotic treatments.


  Seasonal variation in environmental variables, and in rates of contact among
individuals, are fundamental drivers of infectious disease dynamics. Unlike
most periodically-forced physical systems, for which the precise pattern of
forcing is typically known, underlying patterns of seasonal variation in
transmission rates can be estimated approximately at best, and only the period
of forcing is accurately known. Yet solutions of epidemic models depend
strongly on the forcing function, so dynamical predictions---such as changes in
epidemic patterns that can be induced by demographic transitions or mass
vaccination---are always subject to the objection that the underlying patterns
of seasonality are poorly specified. Here, we demonstrate that the key
bifurcations of the standard epidemic model are invariant to the shape of
seasonal forcing if the amplitude of forcing is appropriately adjusted.
Consequently, analyses applicable to real disease dynamics can be conducted
with a smooth, idealized sinusoidal forcing function, and qualitative changes
in epidemic patterns can be predicted without precise knowledge of the
underlying forcing pattern. We find similar invariance in a seasonally forced
predator-prey model, and conjecture that this phenomenon---and the associated
robustness of predictions---might be a feature of many other periodically
forced dynamical systems.


  The HANDY model of Motesharrei, Rivas, and Kalnay examines interactions with
the environment by human populations, both between poor and rich people, i.e.,
"Commoners" and "Elites". The Elites control the society's wealth and consume
it at a higher rate than Commoners, whose work produces the wealth. We say a
model is "Elite-dominated" when the Elites' per capita population change rate
is always at least as large as the Commoners'. We can show the HANDY model
always exhibits population crashes for all choices of parameter values for
which it is Elite-dominated. But any such model with explicit equations raises
questions of how the resulting behaviors depend on the details of the models.
How important are the particular design features codified in the differential
equations? In this paper, we first replace the explicit equations of HANDY with
differential equations that are only described conceptually or qualitatively -
using only conditions that can be verified for explicit systems. Next, we
discard the equations entirely, replacing them with qualitative conditions, and
we prove these conditions imply population collapse must occur. In particular,
one condition is that the model is Elite-dominated. We show that the HANDY
model with Elite-dominated parameters satisfies our hypotheses and thus must
undergo population collapse. Our approach of introducing qualitative
mathematical hypotheses can better show the underlying features of the model
that lead to collapse. We also ask how societies can avoid collapse.


  Dengue virus has caused major problems for public health officials for
decades in tropical and subtropical countries. We construct a compartmental
model that includes the risk of hospitalization and its impact on public health
policies. The basic reproductive number, $\mathcal{R}_0$, is computed, as well
as a sensitivity analysis on $\mathcal{R}_0$ parameters and discuss the
relevance in public health policies. The local and global stability of the
disease-free equilibrium is established. Numerical simulations are performed to
better determine future prevention/control strategies.


  Recent work highlights that cooperation in the one-shot Prisoner's dilemma
(PD) is primarily driven by moral preferences for doing the right thing, rather
than social preferences for equity or efficiency. By contrast, little is known
on what motivates cooperation in the Stag-Hunt Game (SHG). Cooperation in the
SHG fundamentally differs from cooperation in the PD in that it is not costly,
but risky: players have no temptation to deviate from the cooperative outcome,
but cooperation only pays off if the other player cooperates. Here, we provide
data from a large (N=436), pre-registered, experiment. Contrary to what has
been observed for the PD, we find that SHG cooperation is primarily driven by
preferences for efficiency, rather than preferences for doing the right thing.


  Cryo-Electron Microscopy (cryo-EM) has become an extremely powerful method
for resolving structural details of large biomolecular complexes. However,
challenging problems in single-particle methods remain open because of (1) the
low signal-to-noise ratio in EM; and (2) the potential anisotropy and lack of
coverage of projection directions relative to the body-fixed coordinate system
for some complexes. Whereas (1) is usually addressed by class averaging (and
increasingly due to rapid advances in microscope and sensor technology), (2) is
an artifact of the mechanics of interaction of biomolecular complexes and the
vitrification process. In the absence of tilt series, (2) remains a problem,
which is addressed here by supplementing EM data with Small-Angle X-Ray
Scattering (SAXS). Whereas SAXS is of relatively low resolution and contains
much lower information content than EM, we show that it is nevertheless
possible to use SAXS to fill in blind spots in EM in difficult cases where the
range of projection directions is limited.


  When can complex ecological interactions drive an entire ecosystem into a
persistent non-equilibrium state, where species abundances keep fluctuating
without going to extinction? We show that high-diversity spatially-extended
systems, in which conditions vary somewhat between spatial locations, can
exhibit chaotic dynamics which persist for extremely long times. We develop a
theoretical framework, based on dynamical mean-field theory, to quantify the
conditions under which these fluctuating states exist, and predict their
properties. We uncover parallels with the persistence of externally-perturbed
ecosystems, such as the role of perturbation strength, synchrony and
correlation time. But uniquely to endogenous fluctuations, these properties
arise from the species dynamics themselves, creating feedback loops between
perturbation and response. A key result is that the fluctuation amplitude and
species diversity are tightly linked, in particular fluctuations enable
dramatically more species to coexist than at equilibrium in the very same
system. Our findings highlight crucial differences between well-mixed and
spatially-extended systems, with implications for experiments and their ability
to reproduce natural dynamics. They shed light on the maintenance of
biodiversity, and the strength and synchrony of fluctuations observed in
natural systems.


  A hallmark of the adaptive immune response is the proliferation of
pathogen-specific lymphocytes that leave in their wake a long lived population
of cells that provide lasting immunity. A subject of ongoing investigation is
when during an adaptive immune response those memory cells are produced. In two
ground-breaking studies, Buchholz et al. (Science, 2013) and Gerlach et al.
(Science, 2013) employed experimental methods that allowed identification of
offspring from individual lymphocytes in vivo, which we call clonal data, at a
single time point. Through the development, application and fitting of a
mathematical model, Buchholz et al. (Science, 2013) concluded that, if memory
is produced during the expansion phase, memory cell precursors are made before
the effector cells that clear the original pathogen. We sought to determine the
general validity and power of the modeling approach introduced in Buchholz et
al. (Science, 2013) for quickly evaluating differentiation networks by adapting
it to make it suitable for drawing inferences from more readily available
non-clonal phenotypic proportion time-courses. We first established the method
drew consistent deductions when fit to the non-clonal data in Buchholz et al.
(Science, 2013) itself. We fit a variant of the model to data reported in
Badovinac et al. (J. Immun., 2007), Schlub et al. (Immun. & Cell Bio., 2010),
and Kinjo et al. (Nature Commun., 2015) with necessary simplifications to match
different reported data in these papers. The deduction from the model was
consistent with that in Buchholz et al. (Science, 2013), albeit with
questionable parameterizations. An alternative possibility, supported by the
data in Kinjo et al. (Nature Commun., 2015), is that memory precursors are
created after the expansion phase, which is a deduction not possible from the
mathematical methods provided in Buchholz et al. (Science, 2013).


  Gene regulation is one of the most important fundamental biological processes
in living cells. It involves multiple protein molecules that locate specific
sites on DNA and assemble gene initiation or gene repression multi-molecular
complexes. While the protein search dynamics for DNA targets has been
intensively investigated, the role of inter-molecular interactions during the
genetic activation or repression remains not well quantified. Here we present a
simple one-dimensional model of target search for two interacting molecules
that can reversibly form a dimer molecular complex, which also participates in
the search process. In addition, the proteins have finite residence times on
specific target sites, and the gene is activated or repressed when both
proteins are simultaneously present at the target. The model is analyzed using
first-passage analytical calculations and Monte Carlo computer simulations. It
is shown that the search dynamics exhibits a complex behavior depending on the
strength of inter-molecular interactions and on the target residence times. We
also found that the search time shows a non-monotonic behavior as a function of
the dissociation rate for the molecular complex. Physical-chemical arguments to
explain these observations are presented. Our theoretical approach highlights
the importance of molecular interactions in the complex process of gene
activation/repression by multiple transcription factor proteins.


  In evolutionary dynamics, a key measure of a mutant trait's success is the
probability that it takes over the population given some initial
mutant-appearance distribution. This "fixation probability" is difficult to
compute in general, as it depends on the mutation's effect on the organism as
well as the population's spatial structure, mating patterns, and other factors.
In this study, we consider weak selection, which means that the mutation's
effect on the organism is small. We obtain a weak-selection perturbation
expansion of a mutant's fixation probability, from an arbitrary initial
configuration of mutant and resident types. Our results apply to a broad class
of stochastic evolutionary models, in which the size and spatial structure are
arbitrary (but fixed). The problem of whether selection favors a given trait is
thereby reduced from exponential to polynomial complexity in the population
size, when selection is weak. We conclude by applying these methods to obtain
new results for evolutionary dynamics on graphs.


  Transcriptomics response of SK-N-AS cells to methamidophos (an acetylcholine
esterase inhibitor) exposure was measured at 10 time points between 0.5 and 48
h. The data was analyzed using a combination of traditional statistical methods
and novel machine learning algorithms for detecting anomalous behavior and
infer causal relations between time profiles. We identified several processes
that appeared to be upregulated in cells treated with methamidophos including:
unfolded protein response, response to cAMP, calcium ion response, and
cell-cell signaling. The data confirmed the expected consequence of
acetylcholine buildup. In addition, transcripts with potentially key roles were
identified and causal networks relating these transcripts were inferred using
two different computational methods: Siamese convolutional networks and time
warp causal inference. Two types of anomaly detection algorithms, one based on
Autoencoders and the other one based on Generative Adversarial Networks (GANs),
were applied to narrow down the set of relevant transcripts.