  Microbial ecosystems are remarkably diverse, stable, and often consist of a
balanced mixture of core and peripheral species. Here we propose a conceptual
model exhibiting all these emergent properties in quantitative agreement with
real ecosystem data, specifically species' abundance and prevalence
distributions. Resource competition and metabolic commensalism drive stochastic
ecosystem assembly in our model. We demonstrate that even when supplied with
just one resource, ecosystems can exhibit high diversity, increasing stability,
and partial reproducibility between samples.


  We introduce the Dynamic Influence Network (DIN), a novel visual analytics
technique for representing and analyzing rule-based models of protein-protein
interaction networks. Rule-based modeling has proved instrumental in developing
biological models that are concise, comprehensible, easily extensible, and that
mitigate the combinatorial complexity of multi-state and multi-component
biological molecules. Our technique visualizes the dynamics of these rules as
they evolve over time. Using the data produced by KaSim, an open source
stochastic simulator of rule-based models written in the Kappa language, DINs
provide a node-link diagram that represents the influence that each rule has on
the other rules. That is, rather than representing individual biological
components or types, we instead represent the rules about them (as nodes) and
the current influence of these rules (as links). Using our interactive DIN-Viz
software tool, researchers are able to query this dynamic network to find
meaningful patterns about biological processes, and to identify salient aspects
of complex rule-based models. To evaluate the effectiveness of our approach, we
investigate a simulation of a circadian clock model that illustrates the
oscillatory behavior of the KaiC protein phosphorylation cycle.


  The Na$^+$/K$^+$ ATPase is an essential component of cardiac
electrophysiology, maintaining physiological Na$^+$ and K$^+$ concentrations
over successive heart beats. Terkildsen et al. (2007) developed a model of the
ventricular myocyte Na$^+$/K$^+$ ATPase to study extracellular potassium
accumulation during ischaemia, demonstrating the ability to recapitulate a wide
range of experimental data, but unfortunately there was no archived code
associated with the original manuscript. Here we detail an updated version of
the model and provide CellML and MATLAB code to ensure reproducibility and
reusability. We note some errors within the original formulation which have
been corrected to ensure that the model is thermodynamically consistent, and
although this required some reparameterisation, the resulting model still
provides a good fit to experimental measurements that demonstrate the
dependence of Na$^+$/K$^+$ ATPase pumping rate upon membrane voltage and
metabolite concentrations. To demonstrate thermodynamic consistency we also
developed a bond graph version of the model. We hope that these models will be
useful for community efforts to assemble a whole-cell cardiomyocyte model which
facilitates the investigation of cellular energetics.


  Cooperation is ubiquitous at every level of living organisms. It is known
that spatial (network) structure is a viable mechanism for cooperation to
evolve. A recently proposed numerical metric, average gradient of selection
(AGoS), a useful tool for interpreting and visualizing evolutionary dynamics on
networks, allows simulation results to be visualized on a one-dimensional phase
space. However, stochastic mutation of strategies was not considered in the
analysis of AGoS. Here we extend AGoS so that it can analyze the evolution of
cooperation where mutation may alter strategies of individuals on networks. We
show that our extended AGoS correctly visualizes the final states of
cooperation with mutation in the individual-based simulations. Our analyses
revealed that mutation always has a negative effect on the evolution of
cooperation regardless of the payoff functions, fraction of cooperators, and
network structures. Moreover, we found that scale-free networks are the most
vulnerable to mutation and thus the dynamics of cooperation are altered from
bistability to coexistence on those networks, undergoing an imperfect pitchfork
bifurcation.


  Motivation: Seed location filtering is critical in DNA read mapping, a
process where billions of DNA fragments (reads) sampled from a donor are mapped
onto a reference genome to identify genomic variants of the donor.
State-of-the-art read mappers 1) quickly generate possible mapping locations
for seeds (i.e., smaller segments) within each read, 2) extract reference
sequences at each of the mapping locations, and 3) check similarity between
each read and its associated reference sequences with a
computationally-expensive algorithm (i.e., sequence alignment) to determine the
origin of the read. A seed location filter comes into play before alignment,
discarding seed locations that alignment would deem a poor match. The ideal
seed location filter would discard all poor match locations prior to alignment
such that there is no wasted computation on unnecessary alignments.
  Results: We propose a novel seed location filtering algorithm, GRIM-Filter,
optimized to exploit 3D-stacked memory systems that integrate computation
within a logic layer stacked under memory layers, to perform
processing-in-memory (PIM). GRIM-Filter quickly filters seed locations by 1)
introducing a new representation of coarse-grained segments of the reference
genome, and 2) using massively-parallel in-memory operations to identify read
presence within each coarse-grained segment. Our evaluations show that for a
sequence alignment error tolerance of 0.05, GRIM-Filter 1) reduces the false
negative rate of filtering by 5.59x--6.41x, and 2) provides an end-to-end read
mapper speedup of 1.81x--3.65x, compared to a state-of-the-art read mapper
employing the best previous seed location filtering algorithm.
  Availability: The code is available online at:
https://github.com/CMU-SAFARI/GRIM


  The first of $N$ identical independently distributed (i.i.d.) Brownian
trajectories that arrives to a small target, sets the time scale of activation,
which in general is much faster than the arrival to the target of only a single
trajectory. Analytical asymptotic expressions for the minimal time is
notoriously difficult to compute in general geometries. We derive here
asymptotic laws for the probability density function of the first and second
arrival times of a large number of i.i.d. Brownian trajectories to a small
target in 1,2, and 3 dimensions and study their range of validity by stochastic
simulations. The results are applied to activation of biochemical pathways in
cellular transduction.


  We report on the successful synthesis and hyperpolarization of N unprotected
{\alpha} amino acid ethyl acrylate esters and extensively, on an alanine
derivative hyperpolarized by PHIP (4.4$\pm$1% $^{13}$C-polarization), meeting
required levels for in vivo detection. Using water as solvent increases
biocompatibility and the absence of N-protection is expected to maintain
biological activity.


  MicroRNAs (miRNAs) are small non-coding RNAs that control gene expression at
the post-transcriptional level through complementary base pairing with the
target mRNA, leading to mRNA degradation and blocking translation process. Any
dysfunctions of these small regulatory molecules have been linked with the
development and progression of several diseases. Therefore, it is necessary to
reliably predict potential miRNA targets. A large number of computational
prediction tools have been developed which provide a faster way to find
putative miRNA targets, but at the same time their results are often
inconsistent. Hence, finding a reliable, functional miRNA target is still a
challenging task. Also, each tool is equipped with different algorithms, and it
is difficult for the biologists to know which tool is the best choice for their
study. This paper briefly describes fundamental of miRNA target prediction
algorithms, discuss frequently used prediction tools, and further, the
performance of frequently used prediction tools have been assessed using
experimentally validated high confident mature miRNAs and their targets for two
organisms Human and Drosophila Melanogaster. Both Drosophila Melanogaster and
Human supported miRNA target prediction tools have been evaluated separately to
find out best performing tool for each of these two organisms. In the human
dataset, TargetScan showed the best results amongst the other predictors
followed by the miRmap and microT, whereas in the D. Melanogaster dataset,
MicroT tool showed the best performance followed by the TargetScan in the
comparison of other tools.


  The elucidation of entropic contribution to enzyme catalysis has been debated
over decades. The recent experimentally measured activation enthalpy and
entropy, for chorismate rearrangement reaction in PchB brings up a hotly
debated issue whether the chorismate mutase catalyzed reaction is
entropy-driven reaction. Extensive configurational sampling combined with
quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) provides an
approach to calculate entropic contribution in condensed phase reactions.
Complete reaction pathway is exploited by QM/MM MD simulations at DFT and
SCC-DFTB levels. The overall entropy change calculated at SCC-DFTB level QM/MM
MD simulations, is close agreement with the experimental value. Conformation
analysis indicates that the self-ordering of chorismate in the active site of
PchB also contributes to total entropy change. This entropy penalty including
conformational transformation entropy and activation entropy cannot be
intuitively speculated from the crystal structure that only acts as a
stationary state along the reaction pathway of PchB catalyzed reaction. This is
the first time to use QM/MM MD simulations to calculate the activation entropy
from the temperature dependence of reliable free energy profiles with extensive
simulation time. The reasonable insight in enthalpy/entropy scheme clarifies
the detailed entropy change and provides a quantitative tool to the
contradicted experimental results.


  Methods for efficiently controlling dynamics propagated on networks are
usually based on identifying the most influential nodes. Knowledge of these
nodes can be used for the targeted control of dynamics such as epidemics, or
for modifying biochemical pathways relating to diseases. Similarly they are
valuable for identifying points of failure to increase network resilience in,
for example, social support networks and logistics networks. Many measures,
often termed `centrality', have been constructed to achieve these aims. Here we
consider Katz centrality and provide a new interpretation as a steady-state
solution to continuous-time dynamics. This enables us to implement a
sensitivity analysis which is similar to metabolic control analysis used in the
analysis of biochemical pathways. The results yield a centrality which
quantifies, for each node, the net impact of its absence from the network. It
also has the desirable property of requiring a node with a high centrality to
play a central role in propagating the dynamics of the system by having the
capacity to both receive flux from others and then to pass it on. This new
perspective on Katz centrality is important for a more comprehensive analysis
of directed networks.


  Topology affects physical and biological properties of DNA and impacts
fundamental cellular processes, such as gene expression, genome replication,
chromosome structure and segregation. In all organisms DNA topology is
carefully modulated and the supercoiling degree of defined genome regions may
change according to physiological and environmental conditions. Elucidation of
structural properties of DNA molecules with different topology may thus help to
better understand genome functions. Whereas a number of structural studies have
been published on highly negatively supercoiled DNA molecules, only preliminary
observations of highly positively supercoiled are available, and a description
of DNA structural properties over the full range of supercoiling degree is
lacking. Atomic Force Microscopy (AFM) is a powerful tool to study DNA
structure at single molecule level. We here report a comprehensive analysis by
AFM of DNA plasmid molecules with defined supercoiling degree, covering the
full spectrum of biologically relevant topologies, under different observation
conditions. Our data, supported by statistical and biochemical analyses,
revealed striking differences in the behavior of positive and negative plasmid
molecules.


  Cell adhesion complexes (CACs), which are activated by ligand binding, play
key roles in many cellular functions ranging from cell cycle regulation to
mediation of cell extracellular matrix adhesion. Inspired by single molecule
pulling experiments on leukocyte function-associated antigen-1 (LFA-1),
expressed in T-cells, bound to intercellular adhesion molecules (ICAM), we
performed constant loading rate ($r_f$) and constant force ($F$) simulations
using the Self-Organized Polymer (SOP) model to describe the mechanism of
ligand rupture from CACs. The simulations reproduce the major experimental
finding on the kinetics of the rupture process, namely, the dependence of the
most probable rupture forces ($f^*$s) on $\ln r_f$ ($r_f$ is the loading rate)
exhibits two distinct linear regimes. The first, at low $r_f$, has a shallow
slope whereas the slope at high $r_f$ is much larger, especially for
LFA-1/ICAM-1 complex with the transition between the two occurring over a
narrow $r_f$ range. Locations of the two transition states (TSs), extracted
from the simulations show an abrupt change from a high value at low $r_f$ or
$F$ to a low value at high $r_f$ or $F$. The unusual behavior in which the CACs
switch from one brittle (TS position is a constant over a range of forces)
state to another brittle state is not found in forced-rupture in other protein
complexes. We explain this novel behavior by constructing the free energy
profiles, $F(\Lambda)$s, as a function of a collective reaction coordinate
($\Lambda$), involving many key charged residues and a critical metal ion. The
TS positions in F($\Lambda) change abruptly at a critical force, demonstrating
that it, rather than the molecular extension is a good reaction coordinate. We
reveal a new mechanism for the two loading regimes observed in the rupture
kinetics in CACs.


  Through the lens of game theory, cooperation is frequently considered an
unsustainable strategy: if an entire population is cooperating, each indi-
vidual can increase its overall fitness by choosing not to cooperate, thereby
still receiving all the benefit of its cooperating neighbors while no longer
expending its own energy. Observable cooperation in naturally-occurring public
goods games is consequently of great interest, as such systems offer insight
into both the emergence and sustainability of cooperation. Here we consider a
population that obeys a public goods game on a network of discrete regions
(that we call nests), between any two of which individuals are free to migrate.
We construct a system of piecewise-smooth ordinary differential equations that
couple the within-nest population dynamics and the between-nest migratory
dynamics. Through a combination of analytical and numerical methods, we show
that if the workers within the population migrate sufficiently fast relative to
the cheaters, the network loses stability first through a Hopf bifurcation,
then a torus bifurcation, after which one or more nests collapse. Our results
indicate that fast moving cheaters can act to stabilize worker-cheater
coexistence within network that would otherwise collapse. We end with a
comparison of our results with the dynamics observed in colonies of the ant
species Pristomyrmex punctatus and in those of the Cape honeybee Apis mellifera
capensis, and argue that they qualitatively agree.


  We make a mathematical analysis of an age structured HIV infection model with
both virus-to-cell and cell-to-cell transmissions to understand the dynamical
behavior of HIV infection in vivo. In the model, we consider the proliferation
of uninfected CD4+ T cells by a logistic function and the infected CD4+ T cells
are assumed to have an infection-age structure. Our main results concern the
Hopf bifurcation of the model by using the theory of integrated semigroup and
the Hopf bifurcation theory for semilinear equations with nondense domain.
Bifurcation analysis indicates that there exist some parameter values such that
this HIV infection model has a non-trivial periodic solution which bifurcates
from the positive equilibrium. The numerical simulations are also carried out.


  This work reports on two related investigations of stochastic simulations
which are widely used to study biodiversity and other related issues. We first
deal with the behavior of the Hamming distance under the increase of the number
of species and the size of the lattice, and then investigate how the mobility
of the species contributes to jeopardize biodiversity. The investigations are
based on the standard rules of reproduction, mobility and predation or
competition, which are described by specific rules, guided by generalization of
the rock-paper-scissors game, valid in the case of three species. The results
on the Hamming distance indicate that it engenders universal behavior,
independently of the number of species and the size of the square lattice. The
results on the mobility confirm the prediction that it may destroy diversity,
if it is increased to higher and higher values.


  In this paper we derive several quasi steady-state approximations (QSSAs) to
the stochastic reaction network describing the Michaelis-Menten enzyme
kinetics. We show how the different assumptions about chemical species
abundance and reaction rates lead to the standard QSSA (sQSSA), the total QSSA
(tQSSA), and the reverse QSSA (rQSSA) approximations. These three QSSAs have
been widely studied in the literature in deterministic ordinary differential
equation (ODE) settings and several sets of conditions for their validity have
been proposed. By using multiscaling techniques introduced in Kang (and Kurtz
2013) and Ball et al. (2006) we show that these conditions for deterministic
QSSAs largely agree with the ones for QSSAs in the large volume limits of the
underlying stochastic enzyme kinetic network.


  Mutualistic interactions, which are beneficial for both interacting species,
are recurrently present in ecosystems. Observations of natural systems showed
that, if we draw mutualistic relationships as binary links between species, the
resulting bipartite network of interactions displays a widespread particular
ordering called nestedness. On the other hand, theoretical works have shown
that a nested structure has a positive impact on a number of relevant features
ranging from species coexistence to a higher structural stability of
communities and biodiversity. However, how nestedness emerges and what are its
determinants, are still open challenges that have led to multiple debates to
date. Here, we show, by applying a theoretical approach to the analysis of 167
real mutualistic networks, that nestedness is not an irreducible feature, but a
consequence of the degree sequences of both guilds of the mutualistic network.
Remarkably, we find that an outstanding majority of the analyzed networks does
not show statistical significant nestedness. These findings point to the need
of revising previous claims about the role of nestedness and might contribute
to expand our understanding of how evolution shapes mutualistic interactions
and communities by placing the focus on the local properties rather than on
global quantities.


  Global species richness is a key biodiversity metric. Despite recent efforts
to estimate global species richness, the resulting estimates have been highly
uncertain and often logically inconsistent. Estimates lower down either the
taxonomic or geographic hierarchies are often larger than those above. Further,
these estimates have been represented in a wide variety of forms, including
intervals (a, b), point estimates with no uncertainty, and point estimates with
either symmetrical or asymmetrical bounds, making it difficult to combine
information across different estimates. Here, we develop a Bayesian
hierarchical approach to estimate the global species richness from published
studies. It allows us to recover interval estimates at each level of the
hierarchy, even when data are partially or wholly unobserved, while respecting
logical constraints, and to determine the effects of estimation on the whole
hierarchy of obtaining future estimates anywhere within it


  The software program BAMM has been widely used to study the dynamics of
speciation, extinction, and phenotypic evolution on phylogenetic trees. The
program implements a model-based clustering algorithm to identify clades that
share common macroevolutionary rate dynamics and to estimate rate parameters. A
recent simulation study published in Evolution (2017) by Meyer and Wiens (M&W)
claimed that (i) simple ("MS") estimators of diversification rates perform much
better than BAMM, and (ii) evolutionary rates inferred with BAMM are weakly
correlated with the true rates in the generating model. I demonstrate that
their assessment suffers from two major conceptual errors that invalidate both
primary conclusions. These statistical considerations are not specific to BAMM
and apply to all methods for estimating parameters from empirical data where
the true grouping structure of the data is unknown. First, M&W\'s comparisons
between BAMM and MS estimators suffer from false equivalency because the MS
estimators are given perfect prior knowledge of the locations of rate shifts on
the simulated phylogenies. BAMM is given no such information and must
simultaneously estimate the number and location of rate shifts from the data,
thus resulting in a massive degrees-of-freedom advantage for the MS
estimators.When both methods are given equivalent information, BAMM
dramatically outperforms the MS estimators. Second, M&W\'s experimental design
is unable to assess parameter reliability because their analyses conflate small
effect sizes across treatment groups with error in parameter estimates. Nearly
all model-based frameworks for partitioning data are susceptible to the
statistical mistakes in M&W, including popular clustering algorithms in
population genetics, phylogenetics, and comparative methods.


  Protein adsorption/desorption upon nanoparticle surfaces is an important
process to understand for developing new nanotechnology involving biomaterials,
while atomistic picture of the process and its coupling with protein
conformational change is lacking. Here we report our study on the adsorption of
protein GB1 upon a polystyrene nanoparticle surface using atomistic molecular
dynamic simulations. Enabled by metadynamics, we explored the relevant phase
space and identified three protein states; each protein state involved both the
adsorbed and desorbed states. We also studied the change of secondary and
tertiary structures of GB1 during adsorption, and the dominant interactions
between protein and surface in different adsorbing stages. From the simulation
results we obtained a scenario that is more rational and complete than the
conventional one. We believe the new scenario is more appropriate as a
theoretical model in understanding and explaining experimental signals.
Introduction


  The introduction of Wolbachia-carrying mosquitoes into the population has
recently been proposed as an alternative strategy against dengue. Although
laboratory experiments have shown that the Wolbachia bacterium can reduce the
levels of dengue virus in mosquitoes, it is also important to assess the
performance of Wolbachia in reducing the incidence of dengue in human
populations.In this thesis, deterministic mathematical models of human and
mosquito populations in which either one or two dengue serotypes circulate are
developed. We adapt these models to enable the investigation of dengue disease
dynamics in the absence and presence of Wolbachia in order to assess the
performance of Wolbachia as a strategy to reduce human dengue incidence. When
studying the situation in which a single dengue serotype is present in the
population, we consider scenarios where dengue is introduced into the human
population once and multiple times.


  We introduce the \emph{frugal foraging} model in which a forager performs a
discrete-time random walk on a lattice, where each site initially contains
$\mathcal{S}$ food units. The forager metabolizes one unit of food at each step
and starves to death when it last ate $\mathcal{S}$ steps in the past. Whenever
the forager decides to eat, it consumes all food at its current site and this
site remains empty (no food replenishment). The crucial property of the forager
is that it is \emph{frugal} and eats only when encountering food within at most
$k$ steps of starvation. We compute the average lifetime analytically as a
function of frugality threshold and show that there exists an optimal strategy,
namely, a frugality threshold $k^*$ that maximizes the forager lifetime.


  Opioid addiction has become a global epidemic and a national health crisis in
recent years, with the number of opioid overdose fatalities steadily increasing
since the 1990s. In contrast to the dynamics of a typical illicit drug or
disease epidemic, opioid addiction has its roots in legal, prescription
medication - a fact which greatly increases the exposed population and provides
additional drug accessibility for addicts. In this paper, we present a
mathematical model for prescription drug addiction and treatment with
parameters and validation based on data from the opioid epidemic. Key dynamics
considered include addiction through prescription, addiction from illicit
sources, and treatment. Through mathematical analysis, we show that no
addiction-free equilibrium can exist without stringent control over how opioids
are administered and prescribed, effectively transforming the dynamics of the
opioid epidemic into those found in a purely illicit drug model. Numerical
sensitivity analysis suggests that relatively low states of endemic addiction
can be obtained by primarily focusing on medical prevention followed by
aggressive treatment of remaining cases - even when the probability of relapse
from treatment remains high. Further empirical study focused on understanding
the rate of illicit drug dependence versus overdose risk, along with the
current and changing rates of opioid prescription and treatment, would shed
significant light on optimal control efforts and feasible outcomes for this
epidemic and drug epidemics in general.


  In this paper, we study the distribution of the number of internal equilibria
of a multi-player two-strategy random evolutionary game. Using techniques from
the random polynomial theory, we obtain a closed formula for the probability
that the game has a certain number of internal equilibria. In addition, by
employing Descartes' rule of signs and combinatorial methods, we provide useful
estimates for this probability. Finally, we also compare our analytical results
with those obtained from samplings.


  Antibodies, an essential part of our immune system, develop through an
intricate process to bind a wide array of pathogens. This process involves
randomly mutating DNA sequences encoding these antibodies to find variants with
improved binding, though mutations are not distributed uniformly across
sequence sites. Immunologists observe this nonuniformity to be consistent with
"mutation motifs", which are short DNA subsequences that affect how likely a
given site is to experience a mutation. Quantifying the effect of motifs on
mutation rates is challenging: a large number of possible motifs makes this
statistical problem high dimensional, while the unobserved history of the
mutation process leads to a nontrivial missing data problem. We introduce an
$\ell_1$-penalized proportional hazards model to infer mutation motifs and
their effects. In order to estimate model parameters, our method uses a Monte
Carlo EM algorithm to marginalize over the unknown ordering of mutations. We
show that our method performs better on simulated data compared to current
methods and leads to more parsimonious models. The application of proportional
hazards to mutation processes is, to our knowledge, novel and formalizes the
current methods in a statistical framework that can be easily extended to
analyze the effect of other biological features on mutation rates.


  Background: Genome-scale metabolic network models and constraint-based
modeling techniques have become important tools for analyzing cellular
metabolism. Thermodynamically infeasible cy-cles (TICs) causing unbounded
metabolic flux ranges are often encountered. TICs satisfy the mass balance and
directionality constraints but violate the second law of thermodynamics.
Current prac-tices involve implementing additional constraints to ensure not
only optimal but also loopless flux distributions. However, the mixed integer
linear programming problems required to solve become computationally
intractable for genome-scale metabolic models. Results: We aimed to identify
the fewest needed constraints sufficient for optimality under the loop-less
requirement. We found that loopless constraints are required only for the
reactions that share elementary flux modes representing TICs with reactions
that are part of the objective function. We put forth the concept of localized
loopless constraints (LLCs) to enforce this minimal required set of loopless
constraints. By combining with a novel procedure for minimal null-space
calculation, the computational time for loopless flux variability analysis is
reduced by a factor of 10-150 compared to the original loopless constraints and
by 4-20 times compared to the currently fastest method Fast-SNP with the
percent improvement increasing with model size. Importantly, LLCs offer a
scalable strategy for loopless flux calculations for
multi-compartment/multi-organism models of very large sizes (e.g. >104
reactions) not feasible before. Matlab functions are available at
https://github.com/maranasgroup/lll-FVA.


  While linear mixed model (LMM) has shown a competitive performance in
correcting spurious associations raised by population stratification, family
structures, and cryptic relatedness, more challenges are still to be addressed
regarding the complex structure of genotypic and phenotypic data. For example,
geneticists have discovered that some clusters of phenotypes are more
co-expressed than others. Hence, a joint analysis that can utilize such
relatedness information in a heterogeneous data set is crucial for genetic
modeling.
  We proposed the sparse graph-structured linear mixed model (sGLMM) that can
incorporate the relatedness information from traits in a dataset with
confounding correction. Our method is capable of uncovering the genetic
associations of a large number of phenotypes together while considering the
relatedness of these phenotypes. Through extensive simulation experiments, we
show that the proposed model outperforms other existing approaches and can
model correlation from both population structure and shared signals. Further,
we validate the effectiveness of sGLMM in the real-world genomic dataset on two
different species from plants and humans. In Arabidopsis thaliana data, sGLMM
behaves better than all other baseline models for 63.4% traits. We also discuss
the potential causal genetic variation of Human Alzheimer's disease discovered
by our model and justify some of the most important genetic loci.


  Do biological cells sense time by the number of their divisions, a process
that ends at senescence? We consider the question "can the cell's perception of
time be expressed through the length of the shortest telomere?" The answer is
that the absolute time before senescence cannot be expressed by the telomere's
length and that a cell can survive many more divisions than intuitively
expected. This apparent paradox is due to shortening and elongation of the
telomere, which suggests a random walk model of the telomere's length. The
model indicates two phases, first, a determinist drift of the length toward a
quasi-equilibrium state, and second, persistence of the length near an
attracting state for the majority of divisions prior to senescence. The measure
of stability of the latter phase is the expected number of divisions at the
attractor ("lifetime") prior to crossing a threshold to senescence. The
telomerase regulates stability by creating an effective potential barrier that
separates statistically the shortest lifetime from the next shortest. The
random walk has to overcome the barrier in order to extend the range of the
first regime. The model explains how random telomere dynamics underlies the
extension of cell survival time.


  The worldwide surge of multiresistant microbial strains has propelled the
search for alternative treatment options. The study of Protein-Protein
Interactions (PPIs) has been a cornerstone in the clarification of complex
physiological and pathogenic processes, thus being a priority for the
identification of vital components and mechanisms in pathogens. Despite the
advances of laboratorial techniques, computational models allow the screening
of protein interactions between entire proteomes in a fast and inexpensive
manner. Here, we present a supervised machine learning model for the prediction
of PPIs based on the protein sequence. We cluster amino acids regarding their
physicochemical properties, and use the discrete cosine transform to represent
protein sequences. A mesh of classifiers was constructed to create
hyper-specialised classifiers dedicated to the most relevant pairs of molecular
function annotations from Gene Ontology. Based on an exhaustive evaluation that
includes datasets with different configurations, cross-validation and
out-of-sampling validation, the obtained results outscore the state-of-the-art
for sequence-based methods. For the final mesh model using SVM with RBF, a
consistent average AUC of 0.84 was attained.


  There are nowadays a huge load of publications about dengue epidemic models,
which mostly employ deterministic differential equations. The analytical
properties of deterministic models are always of particular interest by many
experts, but their validity - if they can indeed track some empirical data - is
an increasing demand by many practitioners. In this view, the data can tell to
which figure the solutions yielded from the models should be; they drift all
the involving parameters towards the most appropriate values. By prior
understanding of the population dynamics, some parameters with inherently
constant values can be estimated forthwith; some others can sensibly be
guessed. However, solutions from such models using sets of constant parameters
most likely exhibit, if not smoothness, at least noise-free behavior; whereas
the data appear very random in nature. Therefore, some parameters cannot be
constant as the solutions to seemingly appear in a high correlation with the
data. We were aware of impracticality to solve a deterministic model many times
that exhaust all trials of the parameters, or to run its stochastic version
with Monte Carlo strategy that also appeals for a high number of solving
processes. We were also aware that those aforementioned non-constant parameters
can potentially have particular relationships with several extrinsic factors,
such as meteorology and socioeconomics of the human population. We then study
an estimation of time-space-varying parameters within the framework of
variational calculus and investigate how some parameters are related to some
extrinsic factors. Numerical results for some examples with real data are shown
and discussed in details.


  The amplification cycle of many replicators (natural or artificial) involves
the usage of a host compartment, inside of which the replicator express
phenotypic compounds necessary to carry out its genetic replication. For
example, viruses infect cells, where they express their own proteins and
replicate. In this process, the host cell boundary limits the diffusion of the
viral protein products, thereby ensuring that phenotypic compounds, such as
proteins, promote the replication of the genes that encoded them. This role of
maintaining spatial co-localization, also called genotype-phenotype linkage, is
a critical function of compartments in natural selection. In most cases
however, individual replicating elements do not distribute systematically among
the hosts, but are randomly partitioned. Depending on the replicator-to-host
ratio, more than one variant may thus occupy some compartments, blurring the
genotype-phenotype linkage and affecting the effectiveness of natural
selection. We derive selection equations for a variety of such random multiple
occupancy situations, in particular considering the effect of replicator
population polymorphism and internal replication dynamics. We conclude that the
deleterious effect of random multiple occupancy on selection is relatively
benign, and may even completely vanish is some specific cases. In addition,
given that higher mean occupancy allows larger populations to be channeled
through the selection process, and thus provide a better exploration of
phenotypic diversity, we show that it may represent a valid strategy in both
natural and technological cases.


  The fixation probability of a mutant in the evolutionary dynamics of Moran
process is calculated by the Monte-Carlo method on a few families of
clique-based graphs. It is shown that the complete suppression of fixation can
be realized with the generalized clique-wheel graph in the limit of small
wheel-clique ratio and infinite size. The family of clique-star is an
amplifier, and clique-arms graph changes from amplifier to suppressor as the
fitness of the mutant increases. We demonstrate that the overall structure of a
graph can be more important to determine the fixation probability than the
degree or the heat heterogeneity. The dependence of the fixation probability on
the position of the first mutant is discussed.


  Recent development in unmanned UAV technology paved the way for numerous
applications in diverse cross discipline fields. One of the main feature of UAV
s are their portability in terms of size that allows them to navigate through
fairly hostile environments and collect data . This data collection leads to
the interpretation of the behavior and predictability according to the analysis
as presented by data science. The application of UAV to monitor the population
and climate geography is well documented. But the usage of UAV to study the
germs in the atmosphere is not well documented or absent. As air remains one of
main medium of transmission of germs so there must be some kind of signature
specific for a particular kind of germ. Using this as cue in this present
communication a hypothetical model to study the spread of disease is presented.
This model can help the epidemiologists to understand the mechanism of
microbial traffic like for example flu getting transferred within the same
species or cross species ,spatial diffusion like for example human traveling
pattern and newly recognized diseases for example various type of flu and
vector borne diseases like malaria , dengue etc. This model also covers some
relevant scenarios like global climate change, political ecologic emergences of
aerial transmitted diseases.


  Network analyses, such as of gene co-expression networks, metabolic networks
and ecological networks have become a central approach for the systems-level
study of biological data. Several software packages exist for generating and
analyzing such networks, either from correlation scores or the absolute value
of a transformed score called weighted topological overlap (wTO). However,
since gene regulatory processes can up- or down-regulate genes, it is of great
interest to explicitly consider both positive and negative correlations when
constructing a gene co-expression network. Here, we present an R package for
calculating the wTO, that, in contrast to existing packages, explicitly
addresses the sign of the wTO values, and is thus especially valuable for the
analysis of gene regulatory networks. The package includes the calculation of
p-values (raw and adjusted) for each pairwise gene score. Our package also
allows the calculation of networks from time series (without replicates). Since
networks from independent datasets (biological repeats or related studies) are
not the same due to technical and biological noise in the data, we
additionally, incorporated a novel method for calculating a consensus network
(CN) from two or more networks into our R package. We compare our new wTO
package to state of art packages and demonstrate the application of the wTO and
CN functions using 3 independently derived datasets from healthy human
pre-frontal cortex samples. To showcase an example for the time series
application we utilized a metagenomics data set. In this work, we developed a
software package that allows the computation of wTO networks, CNs and a
visualization tool in the R statistical environment. It is publicly available
on CRAN repositories under the GPL-2 Open Source License
(https://cran.r-project.org/web/packages/wTO/).


  Cancer is a heterogeneous disease with diverse molecular etiologies and
outcomes. The Cancer Genome Atlas (TCGA) has released a large compendium of
over 10,000 tumors with RNA-seq gene expression measurements. Gene expression
captures the diverse molecular profiles of tumors and can be interrogated to
reveal differential pathway activations. Deep unsupervised models, including
Variational Autoencoders (VAE) can be used to reveal these underlying patterns.
We compare a one-hidden layer VAE to two alternative VAE architectures with
increased depth. We determine the additional capacity marginally improves
performance. We train and compare the three VAE architectures to other
dimensionality reduction techniques including principal components analysis
(PCA), independent components analysis (ICA), non-negative matrix factorization
(NMF), and analysis of gene expression by denoising autoencoders (ADAGE). We
compare performance in a supervised learning task predicting gene inactivation
pan-cancer and in a latent space analysis of high grade serous ovarian cancer
(HGSC) subtypes. We do not observe substantial differences across algorithms in
the classification task. VAE latent spaces offer biological insights into HGSC
subtype biology.


  We propose and analyze a mathematical model of cargo transport by
non-processive molecular motors. In our model, the motors change states by
random discrete events (corresponding to stepping and binding/unbinding), while
the cargo position follows a stochastic differential equation (SDE) that
depends on the discrete states of the motors. The resulting system for the
cargo position is consequently an SDE that randomly switches according to a
Markov jump process governing motor dynamics. To study this system we (1) cast
the cargo position in a renewal theory framework and generalize the renewal
reward theorem and (2) decompose the continuous and discrete sources of
stochasticity and exploit a resulting pair of disparate timescales. With these
mathematical tools, we obtain explicit formulas for experimentally measurable
quantities, such as cargo velocity and run length. Analyzing these formulas
then yields some predictions regarding so-called non-processive clustering, the
phenomenon that a single motor cannot transport cargo, but that two or more
motors can. We find that having motor stepping, binding, and unbinding rates
depend on the number of bound motors, due to geometric effects, is necessary
and sufficient to explain recent experimental data on non-processive motors.


  This paper introduces the application of game theory to understand noisy
real-time signalling and the resulting behavioural dynamics in microscopic
populations such as bacteria and other cells. It presents a bridge between the
fields of molecular communication and microscopic game theory. Molecular
communication uses conventional communication engineering theory and techniques
to study and design systems that use chemical molecules as information
carriers. Microscopic game theory models interactions within and between
populations of cells and microorganisms. Integrating these two fields provides
unique opportunities to understand and control microscopic populations that
have imperfect signal propagation. Two examples, namely bacteria quorum sensing
and tumour cell signalling, are presented with potential games to demonstrate
the application of this approach. Finally, a case study of bacteria resource
sharing demonstrates how noisy signalling can alter the distribution of
behaviour.


  Phylogenetic networks are a generalisation of phylogenetic trees that allow
for more complex evolutionary histories that include hybridisation-like
processes. It is of considerable interest whether a network can be considered
`tree-like' or not, which lead to the introduction of \textit{tree-based}
networks in the rooted, binary context. Tree-based networks are those networks
which can be constructed by adding additional edges into a given phylogenetic
tree, called the \textit{base tree}. Previous extensions have considered
extending to the binary, unrooted case and the nonbinary, rooted case. We
extend tree-based networks to the context of unrooted, nonbinary networks in
three ways, depending on the types of additional edges that are permitted. A
phylogenetic network in which every embedded tree is a base tree is termed a
\textit{fully tree-based} network. We also extend this concept to unrooted,
nonbinary phylogenetic networks and classify the resulting networks. We also
derive some results on the colourability of tree-based networks, which can be
useful to determine whether a network is tree-based.


  We perform individual-based Monte Carlo simulations in a community consisting
of two predator species competing for a single prey species, with the purpose
of studying biodiversity stabilization in this simple model system. Predators
are characterized with predation efficiency and death rates, to which Darwinian
evolutionary adaptation is introduced. Competition for limited prey abundance
drives the populations' optimization with respect to predation efficiency and
death rates. We study the influence of various ecological elements on the final
state, finding that both indirect competition and evolutionary adaptation are
insufficient to yield a stable ecosystem. However, stable three-species
coexistence is observed when direct interaction between the two predator
species is implemented.


  CABLE is a land surface model (LSM) that can be applied stand-alone, as well
as providing for land surface-atmosphere exchange within the Australian
Community Climate and Earth System Simulator (ACCESS). We describe critical new
developments that extend the applicability of CABLE for regional and global
carbon-climate simulations, accounting for vegetation response to biophysical
and anthropogenic forcings. A land-use and land-cover change module, driven by
gross land-use transitions and wood harvest area was implemented, tailored to
the needs of the Coupled Model Intercomparison Project-6 (CMIP6). Novel aspects
include the treatment of secondary woody vegetation, which benefits from a
tight coupling between the land-use module and the Population Orders Physiology
(POP) module for woody demography and disturbance-mediated landscape
heterogeneity. Land-use transitions and harvest associated with secondary
forest tiles modify the annually-resolved patch age distribution within
secondary-vegetated tiles, in turn affecting biomass accumulation and turnover
rates and hence the magnitude of the secondary forest sink. Additionally, we
implemented a novel approach to constrain modelled GPP consistent with the
Co-ordination Hypothesis, predicted by evolutionary theory, which suggests that
electron transport and Rubisco-limited rates adjust seasonally and across
biomes to be co-limiting. We show that the default prior assumption - common to
CABLE and other LSMs - of a fixed ratio of electron transport to carboxylation
capacity is at odds with this hypothesis, and implement an alternative
algorithm for dynamic optimisation of this ratio, such that co-ordination is
achieved as an outcome of fitness maximisation. Results have significant
implications the magnitude of the simulated CO2 fertilisation effect on
photosynthesis in comparison to alternative estimates and observational
proxies.


  The conflict between pro-self and pro-social behaviour is at the core of many
key problems of our time, as, for example, the reduction of air pollution and
the redistribution of scarce resources. For the well-being of our societies, it
is thus crucial to find mechanisms to promote pro-social choices over egoistic
ones. Particularly important, because cheap and easy to implement, are those
mechanisms that can change people's behaviour without forbidding any options or
significantly changing their economic incentives, the so-called "nudges".
Previous research has found that moral nudges (e.g., making norms salient) can
promote pro-social behaviour. However, little is known about whether their
effect persists over time and spills across context. This question is key in
light of research showing that pro-social actions are often followed by selfish
actions, thus suggesting that some moral manipulations may backfire. Here we
present a class of simple moral nudges that have a great positive impact on
pro-sociality. In Studies 1-4 (total N = 1,400), we use economic games to
demonstrate that asking subjects to self-report "what they think is the morally
right thing to do" does not only increase pro-sociality in the choice
immediately after, but also in subsequent choices, and even when the social
context changes. In Study 5, we explore whether moral nudges promote charity
donations to humanitarian organisations in a large (N=1,800) crowdfunding
campaign. We find that, in this context, moral nudges increase donations by
about 44 percent.


  We propose a new mathematical model for the spread of Zika virus. Special
attention is paid to the transmission of microcephaly. Numerical simulations
show the accuracy of the model with respect to the Zika outbreak occurred in
Brazil.


  Opioid related deaths are increasing dramatically in recent years, and opioid
epidemic is worsening in the United States. Combating opioid epidemic becomes a
high priority for both the U.S. government and local governments such as New
York State. Analyzing patient level opioid related hospital visits provides a
data driven approach to discover both spatial and temporal patterns and
identity potential causes of opioid related deaths, which provides essential
knowledge for governments on decision making. In this paper, we analyzed opioid
poisoning related hospital visits using New York State SPARCS data, which
provides diagnoses of patients in hospital visits. We identified all patients
with primary diagnosis as opioid poisoning from 2010-2014 for our main studies,
and from 2003-2014 for temporal trend studies. We performed demographical based
studies, and summarized the historical trends of opioid poisoning. We used
frequent item mining to find co-occurrences of diagnoses for possible causes of
poisoning or effects from poisoning. We provided zip code level spatial
analysis to detect local spatial clusters, and studied potential correlations
between opioid poisoning and demographic and social-economic factors.


  The collection of immunoglobulin genes in an individual's germline, which
gives rise to B cell receptors via recombination, is known to vary
significantly across individuals. In humans, for example, each individual has
only a fraction of the several hundred known V alleles. Furthermore, the
currently-accepted set of known V alleles is both incomplete (particularly for
non-European samples), and contains a significant number of spurious alleles.
The resulting uncertainty as to which immunoglobulin alleles are present in any
given sample results in inaccurate B cell receptor sequence annotations, and in
particular inaccurate inferred naive ancestors. In this paper we first show
that the currently widespread practice of aligning each sequence to its closest
match in the full set of IMGT alleles results in a very large number of
spurious alleles that are not in the sample's true set of germline V alleles.
We then describe a new method for inferring each individual's germline gene set
from deep sequencing data, and show that it improves upon existing methods by
making a detailed comparison on a variety of simulated and real data samples.
This new method has been integrated into the partis annotation and clonal
family inference package, available at https://github.com/psathyrella/partis,
and is run by default without affecting overall run time.


  We examine a conservation problem in which the recovery of an endangered
species depends on a captive breeding and reintroduction program. The model is
applied to the case of the black-footed ferret (Mustela nigripes), an
endangered species in North America reliant on captive breeding for survival.
The timing of reintroduction is an important concern in these programs as there
is a tradeoff between the duration (and therefore the cost) of the captive
breeding program and the period the population spends in recovery and in the
wild. In this paper, we develop a stylized bioeconomic model to determine the
optimal reintroduction time, in which the objective is to minimize the cost of
reintroduction while providing a viably-sized population in the wild. Our
control variable is the timing of reintroduction, which departs from a large
body of work in bioeconomics that focuses on adjustable controls that directly
affect the target population. Generally, we find it is optimal to reintroduce
ferrets early in a reintroduction program, although this result is contingent
on species interactions and provisioning services.


  I extend the traditional SAR, which has achieved status of ecological law and
plays a critical role in global biodiversity assessment, to the general (alpha-
or beta-diversity in Hill numbers) diversity area relationship (DAR). The
extension was motivated to remedy the limitation of traditional SAR that only
address one aspect of biodiversity scaling, i.e., species richness scaling over
space. The extension was made possible by the fact that all Hill numbers are in
units of species (referred to as the effective number of species or as species
equivalents), and I postulated that Hill numbers should follow the same or
similar pattern of SAR. I selected three DAR models, the traditional power law
(PL), PLEC (PL with exponential cutoff) and PLIEC (PL with inverse exponential
cutoff). I defined three new concepts and derived their quantifications: (i)DAR
profile: z-q series where z is the PL scaling parameter at different diversity
order (q); (ii)PDO (pair-wise diversity overlap) profile: g-q series where g is
the PDO corresponding to q; (iii) MAD (maximal accrual diversity) profile:
Dmax-q series where Dmax is the MAD corresponding to q. Furthermore, the PDO-g
is quantified based on the self-similarity property of the PL model, and Dmax
can be estimated from the PLEC parameters. The three profiles constitute a
novel DAR approach to biodiversity scaling. I verified the postulation with the
American gut microbiome project (AGP) dataset of 1473 healthy North American
individuals (the largest human dataset from a single project to date). The PL
model was preferred due to its simplicity and established ecological properties
such as self-similarity (necessary for establishing PDO profile), and PLEC has
an advantage in establishing the MAD profile. All three profiles for the AGP
dataset were successfully quantified and compared with existing SAR parameters
in the literature whenever possible.


  Quantifying how many people are or will be sick, and where, is a critical
ingredient in reducing the burden of disease because it helps the public health
system plan and implement effective outbreak response. This process of disease
surveillance is currently based on data gathering using clinical and laboratory
methods; this distributed human contact and resulting bureaucratic data
aggregation yield expensive procedures that lag real time by weeks or months.
The promise of new surveillance approaches using internet data, such as web
event logs or social media messages, is to achieve the same goal but faster and
cheaper. However, prior work in this area lacks a rigorous model of information
flow, making it difficult to assess the reliability of both specific approaches
and the body of work as a whole.
  We model disease surveillance as a Shannon communication. This new framework
lets any two disease surveillance approaches be compared using a unified
vocabulary and conceptual model. Using it, we describe and compare the
deficiencies suffered by traditional and internet-based surveillance, introduce
a new risk metric called deceptiveness, and offer mitigations for some of these
deficiencies. This framework also makes the rich tools of information theory
applicable to disease surveillance. This better understanding will improve the
decision-making of public health practitioners by helping to leverage
internet-based surveillance in a way complementary to the strengths of
traditional surveillance.


  Pandemic influenza has the epidemic potential to kill millions of people.
While various preventive measures exist (i.a., vaccination and school
closures), deciding on strategies that lead to their most effective and
efficient use remains challenging. To this end, individual-based
epidemiological models are essential to assist decision makers in determining
the best strategy to curb epidemic spread. However, individual-based models are
computationally intensive and it is therefore pivotal to identify the optimal
strategy using a minimal amount of model evaluations. Additionally, as
epidemiological modeling experiments need to be planned, a computational budget
needs to be specified a priori. Consequently, we present a new sampling
technique to optimize the evaluation of preventive strategies using fixed
budget best-arm identification algorithms. We use epidemiological modeling
theory to derive knowledge about the reward distribution which we exploit using
Bayesian best-arm identification algorithms (i.e., Top-two Thompson sampling
and BayesGap). We evaluate these algorithms in a realistic experimental setting
and demonstrate that it is possible to identify the optimal strategy using only
a limited number of model evaluations, i.e., 2-to-3 times faster compared to
the uniform sampling method, the predominant technique used for epidemiological
decision making in the literature. Finally, we contribute and evaluate a
statistic for Top-two Thompson sampling to inform the decision makers about the
confidence of an arm recommendation.


  Explosive growth in the amount of genomic data is matched by increasing power
of consumer-grade computers. Even applications that require powerful servers
can be quickly tested on desktop or laptop machines if we can generate
representative samples from large data sets. I describe a fast and
memory-efficient implementation of an on-line sampling method developed for
tape drives 30 years ago. Focusing on genotype files, I test the performance of
this technique on modern solid-state and spinning hard drives, and show that it
performs well compared to a simple sampling scheme. I illustrate its utility by
developing a method to quickly estimate genome-wide patterns of linkage
disequilibrium (LD) decay with distance. I provide open-source software that
samples loci from several variant format files, a separate program that
performs LD decay estimates, and a C++ library that lets developers incorporate
these methods into their own projects.


  In this article, we quantitatively study, through stochastic models, the
efects of several intracellular phenomena, such as cell volume growth, cell
division, gene replication as well as fuctuations of available RNA polymerases
and ribosomes. These phenomena are indeed rarely considered in classic models
of protein production and no relative quantitative comparison among them has
been performed. The parameters for a large and representative class of proteins
are determined using experimental measures. The main important and surprising
conclusion of our study is to show that despite the signifcant fuctuations of
free RNA polymerases and free ribosomes, they bring little variability to
protein production contrary to what has been previously proposed in the
literature. After verifying the robustness of this quite counter-intuitive
result, we discuss its possible origin from a theoretical view, and interpret
it as the result of a mean-feld efect.


  We study the phase transition from the persistence phase to the extinction
phase for the SIRS (susceptible/ infected/ refractory/ susceptible) model of
diseases spreading on the networks. We derive an analytical expression of the
probability for the descendants nodes to re-infect their ancestors nodes. We
find that, in the case of the recovery time $\tau_R$ is larger than the
infection time $\tau_I$, the infection will flow directionally from the
ancestors to the descendants however, the descendants will not able to reinfect
their ancestors during their infection time. This behavior leads us to deduce
that, for this case and when the infection rate $\lambda$ is high enough in
such that, any infected node on the network infects all of its neighbors during
its infection time, SIRS model on the network evolves to extinction state,
where all the nodes on the network become susceptible. Moreover, we assert
that, in order to the infection occurs repeatedly inside the network, the loops
on the network are necessary, which means the clustering coefficient will play
an important role for this model. Hence, unlike the other models such as SIS
model and SIR model, SIRS model has a two critical threshold which separate the
persistence phase from the extinction phase when $\tau_I<\tau_R$. That means,
for fixed values of $\tau_I$ and $\tau_R$ there are a two critical points for
infection rate $\lambda_1$ and $\lambda_2$, where epidemic persists in between
of those two points. We confirm those results numerically by the simulation of
regular one dimensional SIRS system.


  An important problem in phylogenetics is the construction of phylogenetic
trees. One way to approach this problem, known as the supertree method,
involves inferring a phylogenetic tree with leaves consisting of a set $X$ of
species from a collection of trees, each having leaf-set some subset of $X$. In
the 1980's characterizations, certain inference rules were given for when a
collection of 4-leaved trees, one for each 4-element subset of $X$, can all be
simultaneously displayed by a single supertree with leaf-set $X$. Recently, it
has become of interest to extend such results to phylogenetic networks. These
are a generalization of phylogenetic trees which can be used to represent
reticulate evolution (where species can come together to form a new species).
It has been shown that a certain type of phylogenetic network, called a level-1
network, can essentially be constructed from 4-leaved trees. However, the
problem of providing appropriate inference rules for such networks remains
unresolved. Here we show that by considering 4-leaved networks, called
quarnets, as opposed to 4-leaved trees, it is possible to provide such rules.
In particular, we show that these rules can be used to characterize when a
collection of quarnets, one for each 4-element subset of $X$, can all be
simultaneously displayed by a level-1 network with leaf-set $X$. The rules are
an intriguing mixture of tree inference rules, and an inference rule for
building up a cyclic ordering of $X$ from orderings on subsets of $X$ of size
4. This opens up several new directions of research for inferring phylogenetic
networks from smaller ones, which could yield new algorithms for solving the
supernetwork problem in phylogenetics.


  Because human serum transferrin (hTF) exists freely in serum, it is a
potential target for cancer treatment drugs and in curing iron-overloaded
conditions in patients via long-term transfusion therapy. The understanding of
the interactions between hTF and metal ions is very important for biological,
pharmalogical, toxicological, and other protein engineering purposes. In this
paper, a simple linear free energy correlation is proposed to predict the
binding strength between hTF protein and metal cations. The stability constants
for a family of metal-hTF complexes can be correlated to the non-solvation
energies and the radii of cations. The binding strength is determined by both
the physical properties (charge and size or ionic radius) and chemical
properties (non-solvation energy) of a given cation. The binding strengths of
either divalent and trivalent metals can then be predicted systematically.


  We present a result on the number of decoupled molecules for systems binding
two different types of ligands. In the case of $n$ and $2$ binding sites
respectively, we show that, generically, there are $2(n!)^{2}$ decoupled
molecules with the same binding polynomial. For molecules with more binding
sites for the second ligand, we provide computational results.


  Many animal groups are heterogeneous and may even consist of individuals of
different species, called mixed-species flocks. Mathematical and computational
models of collective animal movement behaviour, however, typically assume that
groups and populations consist of identical individuals. In this paper, using
the mathematical framework of the coagulation-fragmentation process, we develop
and analyse a model of merge and split group dynamics, also called
fission-fusion dynamics, for heterogeneous populations that contain two types
(or species) of individuals. We assume that more heterogeneous groups
experience higher split rates than homogeneous groups, forming two daughter
groups whose compositions are drawn uniformly from all possible partitions. We
analytically derive a master equation for group size and compositions and find
mean-field steady-state solutions. We predict that there is a critical group
size below which groups are more likely to be homogeneous and contain the
abundant type/species. Despite the propensity of heterogeneous groups to split
at higher rates, we find that groups are more likely to be heterogeneous but
only above the critical group size. Monte-Carlo simulation of the model show
excellent agreement with these analytical model results. Thus, our model makes
a testable prediction that composition of flocks are group-size dependent and
do not merely reflect the population level heterogeneity. We discuss the
implications of our results to empirical studies on flocking systems.


  We give an new proof of the well-known competitive exclusion principle in the
chemostat model with $n$ species competing for a single resource, for any set
of increasing growth functions. The proof is constructed by induction on the
number of the species, after being ordered. It uses elementary analysis and
comparisons of solutions of ordinary differential equations.


  We introduce a unifying and generalizing framework for complex and detailed
balanced steady states in chemical reaction network theory. To this end, we
generalize the graph commonly used to represent a reaction network.
Specifically, we introduce a graph, called a reaction graph, that has one edge
for each reaction but potentially multiple nodes for each complex. A special
class of steady states, called node balanced steady states, is naturally
associated with such a reaction graph. We show that complex and detailed
balanced steady states are special cases of node balanced steady states by
choosing appropriate reaction graphs. Further, we show that node balanced
steady states have properties analogous to complex balanced steady states, such
as uniqueness and asymptotical stability in each stoichiometric compatibility
class. Moreover, we associate an integer, called the deficiency, to a reaction
graph that gives the number of independent relations in the reaction rate
constants that need to be satisfied for a positive node balanced steady state
to exist.
  The set of reaction graphs (modulo isomorphism) is equipped with a partial
order that has the complex balanced reaction graph as minimal element. We
relate this order to the deficiency and to the set of reaction rate constants
for which a positive node balanced steady state exists.


  Aedes aegypti is the main vector of multiple diseases, such as Dengue, Zika,
and Chikungunya. Due to modifications in weather patterns, its geographical
range is continuously evolving. Temperature is a key factor for its expansion
into regions with cool winters, but rainfall can also have a strong impact on
the colonization of these regions, since larvae emerging after a rainfall are
likely to die at temperatures below $10^{\circ}$C. As climate change is
expected to affect rainfall regimes, with a higher frequency of heavy storms
and an increase in drought-affected areas, it is important to understand how
different rainfall scenarios may shape Ae. aegypti's range. We develop a model
for the population dynamics of Ae. aegypti, coupled with a rainfall model to
study the effect of the temporal distribution of rainfall on mosquito
abundance. Using a fracturing process, we then investigate the effect of a
higher variability in the daily rainfall. As an example, we show that rainfall
distribution is necessary to explain the geographic range of Ae. aegypti in
Taiwan, an island characterized by rainy winters in the north and dry winters
in the south. We also predict that a higher variability in the rainfall time
distribution will decrease the maximum abundance of Ae. aegypti during the
summer. An increase in daily rainfall variability will likewise enhance its
extinction probability. Finally, we obtain a nonlinear relationship between dry
season duration and extinction probability. These findings can have a
significant impact on our ability to predict disease outbreaks.


  The DNA sequencing is the process of identifying the exact order of
nucleotides within a given DNA molecule. The new portable and relatively
inexpensive DNA sequencers, such as Oxford Nanopore MinION, have the potential
to move DNA sequencing outside of laboratory, leading to faster and more
accessible DNA-based diagnostics. However, portable DNA sequencing and analysis
are challenging for mobile systems, owing to high data throughputs and
computationally intensive processing performed in environments with unreliable
connectivity and power.
  In this paper, we provide an analysis of the challenges that mobile systems
and mobile computing must address to maximize the potential of portable DNA
sequencing, and in situ DNA analysis. We explain the DNA sequencing process and
highlight the main differences between traditional and portable DNA sequencing
in the context of the actual and envisioned applications. We look at the
identified challenges from the perspective of both algorithms and systems
design, showing the need for careful co-design.


  This is a short review of two common approximations in stochastic chemical
and biochemical kinetics. It will appear as Chapter 6 in the book "Quantitative
Biology: Theory, Computational Methods and Examples of Models" edited by Brian
Munsky, Lev Tsimring and Bill Hlavacek (to be published in late 2017/2018 by
MIT Press). All chapter references in this article refer to chapters in the
aforementioned book.


  Organisms result from adaptive processes interacting across different time
scales. One such interaction is that between development and evolution. Models
have shown that development sweeps over several traits in a single agent,
sometimes exposing promising static traits. Subsequent evolution can then
canalize these rare traits. Thus, development can, under the right conditions,
increase evolvability. Here, we report on a previously unknown phenomenon when
embodied agents are allowed to develop and evolve: Evolution discovers body
plans robust to control changes, these body plans become genetically
assimilated, yet controllers for these agents are not assimilated. This allows
evolution to continue climbing fitness gradients by tinkering with the
developmental programs for controllers within these permissive body plans. This
exposes a previously unknown detail about the Baldwin effect: instead of all
useful traits becoming genetically assimilated, only traits that render the
agent robust to changes in other traits become assimilated. We refer to this as
differential canalization. This finding also has implications for the
evolutionary design of artificial and embodied agents such as robots: robots
robust to internal changes in their controllers may also be robust to external
changes in their environment, such as transferal from simulation to reality or
deployment in novel environments.


  The present work reports about the dynamics of a collection of randomly
distributed, and randomly oriented, oscillators in 3D space, coupled by an
interaction potential falling as $1/r^3$, where r stands for the inter-particle
distance. This model schematically represents a collection of identical
biomolecules, coherently vibrating at some common frequency, coupled with a
$1/r^3$ potential stemming from the electrodynamic interaction between
oscillating dipoles. The oscillating dipole moment of each molecule being a
direct consequence of its coherent (collective) vibration. By changing the
average distance among the molecules, neat and substantial changes in the power
spectrum of the time variation of a collective observable are found. As the
average intermolecular distance can be varied by changing the concentration of
the solvated molecules, and as the collective variable investigated is
proportional to the projection of the total dipole moment of the model
biomolecules on a coordinate plane, we have found a prospective experimental
strategy of spectroscopic kind to check whether the mentioned intermolecular
electrodynamic interactions can be strong enough to be detectable, and thus to
be of possible relevance to biology.


  Important cellular processes such as migration, differentiation, and
development often rely on precise timing. Yet, the molecular machinery that
regulates timing is inherently noisy. How do cells achieve precise timing with
noisy components? We investigate this question using a first-passage-time
approach, for an event triggered by a molecule that crosses an abundance
threshold and that is regulated by either an accumulating activator or a
diminishing repressor. We find that the optimal strategy corresponds to a
nonlinear increase in the amount of the target molecule over time. Optimality
arises from a tradeoff between minimizing the extrinsic timing noise of the
regulator, and minimizing the intrinsic timing noise of the target molecule
itself. Although either activation or repression outperforms an unregulated
strategy, when we consider the effects of cell division, we find that
repression outperforms activation if division occurs late in the process. Our
results explain the nonlinear increase and low noise of mig-1 gene expression
in migrating neuroblast cells during Caenorhabditis elegans development, and
suggest that mig-1 regulation is dominated by repression for maximal temporal
precision. These findings suggest that dynamic regulation may be a simple and
powerful strategy for precise cellular timing.


  Predicting evolution of expanding populations is critical to control
biological threats such as invasive species and cancer metastasis. Expansion is
primarily driven by reproduction and dispersal, but nature abounds with
examples of evolution where organisms pay a reproductive cost to disperse
faster. When does selection favor this 'survival of the fastest?' We searched
for a simple rule, motivated by evolution experiments where swarming bacteria
evolved into an hyperswarmer mutant which disperses $ \sim 100\%$ faster but
pays a growth cost of $\sim 10 \%$ to make many copies of its flagellum. We
analyzed a two-species model based on the Fisher equation to explain this
observation: the population expansion rate ($v$) results from an interplay of
growth ($r$) and dispersal ($D$) and is independent of the carrying capacity:
$v=2\sqrt{rD}$. A mutant can take over the edge only if its expansion rate
($v_2$) exceeds the expansion rate of the established species ($v_1$); this
simple condition ($v_2 > v_1$) determines the maximum cost in slower growth
that a faster mutant can pay and still be able to take over. Numerical
simulations and time-course experiments where we tracked evolution by imaging
bacteria suggest that our findings are general: less favorable conditions delay
but do not entirely prevent the success of the fastest. Thus, the expansion
rate defines a traveling wave fitness, which could be combined with trade-offs
to predict evolution of expanding populations.


  The field of structural bioinformatics has seen significant advances with the
use of Molecular Dynamics (MD) simulations of biological systems. The MD
methodology has allowed to explain and discover molecular mechanisms in a wide
range of natural processes. There is an impending need to readily share the
ever-increasing amount of MD data, which has been hindered by the lack of
specialized tools in the past. To solve this problem, we present HTMoL, a
state-of-the-art plug-in-free hardware-accelerated web application specially
designed to efficiently transfer and visualize raw MD trajectory files on a web
browser. Now, individual research labs can publish MD data on the Internet, or
use HTMoL to profoundly improve scientific reports by including supplemental MD
data in a journal publication. HTMoL can also be used as a visualization
interface to access MD trajectories generated on a high-performance computer
center directly.
  Availability: HTMoL is available free of charge for academic use. All major
browsers are supported. A complete online documentation including instructions
for download, installation, configuration, and examples is available at the
HTMoL website http://htmol.tripplab.com. Supplementary data are available
online. Corresponding author: mauricio.carrillo@cinvestav.mx


  A recent paper (Manceau and Lambert, 2016) developed a novel approach for
describing two well-defined notions of 'species' based on a phylogenetic tree
and a phenotypic partition. In this paper, we explore some further
combinatorial properties of this approach and describe an extension that allows
an arbitrary number of phenotypic partitions to be combined with a phylogenetic
tree for these two species notions.


  Switching of the direction of flagella rotations is the key control mechanism
governing the chemotactic activity of E. coli and many other bacteria.
Power-law distributions of switching times are most peculiar because their
emergence cannot be deduced from simple thermodynamic arguments. Recently it
was suggested that by adding finite-time correlations into Gaussian
fluctuations regulating the energy height of barrier between the two rotation
states, one can generate a power-law switching statistics. By using a simple
model of a regulatory pathway, we demonstrate that the required amount of
correlated `noise' can be produced by finite number fluctuations of reacting
protein molecules, a condition common to the intracellular chemistry. The
corresponding power-law exponent appears as a tunable characteristic controlled
by parameters of the regulatory pathway network such as equilibrium number of
molecules, sensitivities, and the characteristic relaxation time.


  The antithetic integral feedback motif recently introduced in Briat, Gupta &
Khammash (Cell Systems, 2017) is known to ensure robust perfect adaptation for
the mean dynamics of a given molecular species involved in a complex stochastic
biomolecular reaction network. However, it was observed that it also leads to a
higher variance in the controlled network than that obtained when using a
constitutive (i.e. open-loop) control strategy. This was interpreted as the
cost of the adaptation property and may be viewed as a performance
deterioration for the overall controlled network. To decrease this variance and
improve the performance, we propose to combine the antithetic integral feedback
motif with a negative feedback strategy. Both theoretical and numerical results
are obtained. The theoretical ones are based on a tailored moment closure
method allowing one to obtain approximate expressions for the stationary
variance for the controlled network and predict that the variance can indeed be
decreased by increasing the strength of the negative feedback. Numerical
results verify the accuracy of this approximation and show that the controlled
species variance can indeed be decreased, sometimes below its constitutive
level. Three molecular networks are considered in order to verify the wide
applicability of two types of negative feedback strategies. The main conclusion
is that there is a trade-off between the speed of the settling-time of the mean
trajectories and the stationary variance of the controlled species; i.e.
smaller variance is associated with larger settling-time.


  Often the analysis of time-dependent chemical and biophysical systems
produces high-dimensional time-series data for which it can be difficult to
interpret which individual features are most salient. While recent work from
our group and others has demonstrated the utility of time-lagged co-variate
models to study such systems, linearity assumptions can limit the compression
of inherently nonlinear dynamics into just a few characteristic components.
Recent work in the field of deep learning has led to the development of
variational autoencoders (VAE), which are able to compress complex datasets
into simpler manifolds. We present the use of a time-lagged VAE, or variational
dynamics encoder (VDE), to reduce complex, nonlinear processes to a single
embedding with high fidelity to the underlying dynamics. We demonstrate how the
VDE is able to capture nontrivial dynamics in a variety of examples, including
Brownian dynamics and atomistic protein folding. Additionally, we demonstrate a
method for analyzing the VDE model, inspired by saliency mapping, to determine
what features are selected by the VDE model to describe dynamics. The VDE
presents an important step in applying techniques from deep learning to more
accurately model and interpret complex biophysics.


  In this paper we study group-based Markov models of evolution and their
mixtures. In the algebreo-geometric setting, group-based phylogenetic tree
models correspond to toric varieties, while their mixtures correspond to secant
and join varieties. Determining properties of these secant and join varieties
can aid both in model selection and establishing parameter identifiability.
Here we explore the first natural geometric property of these varieties: their
dimension. The expected projective dimension of the join variety of a set of
varieties is one more than the sum of their dimensions. A join variety that
realizes the expected dimension is nondefective. Nondefectiveness is not only
interesting from a geometric point-of-view, but has been used to establish
combinatorial identifiability for several classes of phylogenetic mixture
models. In this paper, we focus on group-based models where the equivalence
classes of identified parameters are orbits of a subgroup of the automorphism
group of the group defining the model. In particular, we show that, for these
group-based models, the variety corresponding to the mixture of $r$ trees with
$n$ leaves is nondefective when $n \geq 2r+5$. We also give improved bounds for
claw trees and give computational evidence that 2-tree and 3-tree mixtures are
nondefective for small~$n$.


  Nanopore sequencing technology has the potential to render other sequencing
technologies obsolete with its ability to generate long reads and provide
portability. However, high error rates of the technology pose a challenge while
generating accurate genome assemblies. The tools used for nanopore sequence
analysis are of critical importance as they should overcome the high error
rates of the technology. Our goal in this work is to comprehensively analyze
current publicly available tools for nanopore sequence analysis to understand
their advantages, disadvantages, and performance bottlenecks. It is important
to understand where the current tools do not perform well to develop better
tools. To this end, we 1) analyze the multiple steps and the associated tools
in the genome assembly pipeline using nanopore sequence data, and 2) provide
guidelines for determining the appropriate tools for each step. We analyze
various combinations of different tools and expose the tradeoffs between
accuracy, performance, memory usage and scalability. We conclude that our
observations can guide researchers and practitioners in making conscious and
effective choices for each step of the genome assembly pipeline using nanopore
sequence data. Also, with the help of bottlenecks we have found, developers can
improve the current tools or build new ones that are both accurate and fast, in
order to overcome the high error rates of the nanopore sequencing technology.


  Phylogenetic networks are a type of leaf-labelled, acyclic, directed graph
used by biologists to represent the evolutionary history of species whose past
includes reticulation events. A phylogenetic network is tree-child if each
non-leaf vertex is the parent of a tree vertex or a leaf. Up to a certain
equivalence, it has been recently shown that, under two different types of
weightings, edge-weighted tree-child networks are determined by their
collection of distances between each pair of taxa. However, the size of these
collections can be exponential in the size of the taxa set. In this paper, we
show that, if we ignore redundant edges, the same results are obtained with
only a quadratic number of inter-taxa distances by using the shortest distance
between each pair of taxa. The proofs are constructive and give cubic-time
algorithms in the size of the taxa sets for building such weighted networks.


  Module network inference is a statistical method to reconstruct gene
regulatory networks, which uses probabilistic graphical models to learn modules
of coregulated genes and their upstream regulatory programs from genome-wide
gene expression and other omics data. Here we review the basic theory of module
network inference, present protocols for common gene regulatory network
reconstruction scenarios based on the Lemon-Tree software, and show, using
human gene expression data, how the software can also be applied to learn
differential module networks across multiple experimental conditions.


  We study the influence of a randomly switching reproduction-predation rate on
the survival behavior of the non-spatial cyclic Lotka-Volterra model, also
known as the zero-sum rock-paper-scissors game, used to metaphorically describe
the cyclic competition between three species. In large and finite populations,
demographic fluctuations (internal noise) drive two species to extinction in a
finite time, while the species with the smallest reproduction-predation rate is
the most likely to be the surviving one ("law of the weakest"). Here, we model
environmental (external) noise by assuming that the reproduction-predation rate
of the "strongest species" (the fastest to reproduce/predate) in a given static
environment randomly switches between two values corresponding to more and less
favorable external conditions. We study the joint effect of environmental and
demographic noise on the species survival probabilities and on the mean
extinction time. In particular, we investigate whether the survival
probabilities follow the law of the weakest and analyze their dependence of the
external noise intensity and switching rate. Remarkably, when, on average,
there is a finite number of switches prior to extinction, the survival
probability of the predator of the species whose reaction rate switches
typically varies non-monotonically with the external noise intensity (with
optimal survival about a critical noise strength). We also outline the
relationship with the case where all reaction rates switch on markedly
different time scales.


  Clusters of circulating tumor cells (CTCs), although rare, may account for
more than 95% of metastases. Inflammatory breast cancer (IBC) is a highly
aggressive subtype that chiefly metastasizes via CTC clusters. Theory suggests
that physical systems with hierarchical organization tend to be more adaptable
due to their ability to efficiently span the set of available states. We used
the cophenetic correlation coefficient (CCC) to quantify the hierarchical
organization in the expression of collective dissemination associated and IBC
associated genes, and found that the CCC of both gene sets was higher in (a)
epithelial cell lines as compared to mesenchymal cell lines and (b) IBC tumor
samples as compared to non-IBC breast cancer samples. A higher CCC of both
networks was also correlated with a higher rate of metastatic relapse in breast
cancer patients. Gene set enrichment analysis could not provide similar
insights, indicating that the CCC provides additional information regarding the
organizational complexity of gene expression. These results suggest that
retention of epithelial traits in disseminating tumor cells as IBC progresses
promotes successful metastasis and the CCC may be a prognostic factor for IBC.


  CRISPR-Cas is a genetic adaptive immune system unique to prokaryotic cells
used to combat phage and plasmid threats. The host cell adapts by incorporating
DNA sequences from invading phages or plasmids into its CRISPR locus as
spacers. These spacers are expressed as mobile surveillance RNAs that direct
CRISPR-associated (Cas) proteins to protect against subsequent attack by the
same phages or plasmids. The threat from mobile genetic elements inevitably
shapes the CRISPR loci of archaea and bacteria, and simultaneously the
CRISPR-Cas immune system drives evolution of these invaders. Here we highlight
our recent work, as well as that of others, that seeks to understand phage
mechanisms of CRISPR-Cas evasion and conditions for population coexistence of
phages with CRISPR-protected prokaryotes.


  The standard genetic code is well known to be optimized for minimizing the
phenotypic effects of single nucleotide substitutions, a property that was
likely selected for during the emergence of a universal code. Given the fitness
advantage afforded by high standing genetic diversity in a population in a
dynamic environment, it is possible that selection to explore a large fraction
of the space of functional proteins also occurred. To determine whether
selection for such a property played a role during the emergence of the nearly
universal genetic code, we investigated the number of functional variants of
the Escherichia coli PhoQ protein explored at different time scales under
translation using different genetic codes. We found that the standard genetic
code is highly optimal for exploring a large fraction of the space of
functional PhoQ variants at intermediate time scales as compared to random
codes. Environmental changes, in response to which genetic diversity in a
population provides a fitness advantage, are likely to have occurred at these
intermediate time scales. Our results indicate that the ability of the standard
code to explore a large fraction of the space of functional sequence variants
arises from a balance between robustness and flexibility and is largely
independent of the property of the standard code to minimize the phenotypic
effects of mutations. We propose that selection to explore a large fraction of
the functional sequence space while minimizing the phenotypic effects of
mutations contributed towards the emergence of the standard code as the
universal genetic code.


  We have introduced a Pareto sorting algorithm into Synopsis, a de novo design
program that generates synthesizable molecules with desirable properties. We
give a detailed description of the algorithm and illustrate its working in 2
different de novo design settings: the design of putative dual and selective
FGFR and VEGFR inhibitors, and the successful design of organic structure
determining agents (OSDAs) for the synthesis of zeolites. We show that the
introduction of Pareto sorting not only enables the simultaneous optimization
of multiple properties but also greatly improves the performance of the
algorithm to generate molecules with hard-to-meet constraints. This in turn
allows us to suggest approaches to address the problem of false positive hits
in de novo structure based drug design by introducing structural and
physicochemical constraints in the designed molecules, and by forcing essential
interactions between these molecules and their target receptor.


  Abnormal metabolism is an emerging hallmark of cancer. Cancer cells utilize
both aerobic glycolysis and oxidative phosphorylation (OXPHOS) for energy
production and biomass synthesis. Understanding the metabolic reprogramming in
cancer can help design therapies to target metabolism and thereby to improve
prognosis. We have previously argued that more malignant tumors are usually
characterized by a more modular expression pattern of cancer-associated genes.
In this work, we analyzed the expression patterns of metabolism genes in terms
of modularity for 371 hepatocellular carcinoma (HCC) samples from the Cancer
Genome Atlas (TCGA). We found that higher modularity significantly correlated
with glycolytic phenotype, later tumor stages, higher metastatic potential, and
cancer recurrence, all of which contributed to poorer overall prognosis. Among
patients that recurred, we found the correlation of greater modularity with
worse prognosis during early to mid-progression. Furthermore, we developed
metrics to calculate individual modularity, which was shown to be predictive of
cancer recurrence and patients' survival and therefore may serve as a
prognostic biomarker. Our overall conclusion is that more aggressive HCC
tumors, as judged by decreased host survival probability, had more modular
expression patterns of metabolic genes. These results may be used to identify
cancer driver genes and for drug design.


  The twisting deformation of mechanically stretched DNA molecules is studied
by a coarse grained Hamiltonian model incorporating the fundamental
interactions that stabilize the double helix and accounting for the radial and
angular base pair fluctuations. The latter are all the more important at short
length scales in which DNA fragments maintain an intrinsic flexibility. The
presented computational method simulates a broad ensemble of possible molecule
conformations characterized by a specific average twist and determines the
energetically most convenient helical twist by free energy minimization. As
this is done for any external load, the method yields the characteristic
twist-stretch profile of the molecule and also computes the changes in the
macroscopic helix parameters i.e. average diameter and rise distance. It is
predicted that short molecules under stretching should first over-twist and
then untwist by increasing the external load. Moreover, applying a constant
load and simulating a torsional strain which over-twists the helix, it is found
that the average helix diameter shrinks while the molecule elongates, in
agreement with the experimental trend observed in kilo-base long sequences. The
quantitative relation between percent relative elongation and superhelical
density at fixed load is derived. The proposed theoretical model and
computational method offer a general approach to characterize specific DNA
fragments and predict their macroscopic elastic response as a function of the
effective potential parameters of the mesoscopic Hamiltonian.


  Modern algorithms for de novo prediction of protein structures typically
output multiple full-length models (decoys) rather than a single solution.
Subsequent clustering of such decoys is used both to gauge the success of the
modelling and to decide on the most native-like conformation. At the same time,
partial protein models are sufficient for some applications such as
crystallographic phasing by molecular replacement (MR) in particular, provided
these models represent a certain part of the target structure with reasonable
accuracy. Here we propose a novel clustering algorithm that natively operates
in the space of partial models through an approach known as granular clustering
(GC). The algorithm is based on growing local similarities found in a pool of
initial decoys. We demonstrate that the resulting clusters of partial models
provide a substantially more accurate structural detail on the target protein
than those obtained upon a global alignment of decoys. As the result, the
partial models output by our GC algorithm are also much more effective towards
the MR procedure, compared to the models produced by existing software. The
source code is freely available at https://github.com/biocryst/gc


  Humans are the ultimate ecosystem engineers who have profoundly transformed
the world's landscapes in order to enhance their survival. Somewhat
paradoxically, however, sometimes the unforeseen effect of this ecosystem
engineering is the very collapse of the population it intended to protect. Here
we use a spatial version of a standard population dynamics model of ecosystem
engineers to study the colonization of unexplored virgin territories by a small
settlement of engineers. We find that during the expansion phase the population
density reaches values much higher than those the environment can support in
the equilibrium situation. When the colonization front reaches the boundary of
the available space, the population density plunges sharply and attains its
equilibrium value. The collapse takes place without warning and happens just
after the population reaches its peak number. We conclude that overpopulation
and the consequent collapse of an expanding population of ecosystem engineers
is a natural consequence of the nonlinear feedback between the population and
environment variables.