  Network medicine is an emerging area of research dealing with molecular and
genetic interactions, network biomarkers of disease, and therapeutic target
discovery. Large-scale biomedical data generation offers a unique opportunity
to assess the effect and impact of cellular heterogeneity and environmental
perturbations on the observed phenotype. Marrying the two, network medicine
with biomedical data provides a framework to build meaningful models and
extract impactful results at a network level. In this review, we survey
existing network types and biomedical data sources. More importantly, we delve
into ways in which the network medicine approach, aided by phenotype-specific
biomedical data, can be gainfully applied. We provide three paradigms, mainly
dealing with three major biological network archetypes: protein-protein
interaction, expression-based, and gene regulatory networks. For each of these
paradigms, we discuss a broad overview of philosophies under which various
network methods work. We also provide a few examples in each paradigm as a test
case of its successful application. Finally, we delineate several opportunities
and challenges in the field of network medicine. Taken together, the
understanding gained from combining biomedical data with networks can be useful
for characterizing disease etiologies and identifying therapeutic targets,
which, in turn, will lead to better preventive medicine with translational
impacts on personalized healthcare.


  The body is home to a diverse microbiota, mainly in the gut. Resistant
bacteria are selected for by antibiotic treatments, and once resistance becomes
widespread in a population of hosts, antibiotics become useless. Here, we
develop a multiscale model of the interaction between antibiotic use and
resistance spread in a host population, focusing on an important aspect of
within-host immunity. Antibodies secreted in the gut enchain bacteria upon
division, yielding clonal clusters of bacteria. We demonstrate that
immunity-driven bacteria clustering can hinder the spread of a novel resistant
bacterial strain in a host population. We quantify this effect both in the case
where resistance pre-exists and in the case where acquiring a new resistance
mutation is necessary for the bacteria to spread. We further show that the
reduction of spread by clustering can be countered when immune hosts are silent
carriers, and are less likely to get treated, and/or have more contacts. We
demonstrate the robustness of our findings to including stochastic within-host
bacterial growth, a fitness cost of resistance, and its compensation. Our
results highlight the importance of interactions between immunity and the
spread of antibiotic resistance, and argue in the favor of vaccine-based
strategies to combat antibiotic resistance.


  We identify hardware that is optimal to produce molecular dynamics
trajectories on Linux compute clusters with the GROMACS 2018 simulation
package. Therefore, we benchmark the GROMACS performance on a diverse set of
compute nodes and relate it to the costs of the nodes, which may include their
lifetime costs for energy and cooling. In agreement with our earlier
investigation using GROMACS 4.6 on hardware of 2014, the performance to price
ratio of consumer GPU nodes is considerably higher than that of CPU nodes.
However, with GROMACS 2018, the optimal CPU to GPU processing power balance has
shifted even more towards the GPU. Hence, nodes optimized for GROMACS 2018 and
later versions enable a significantly higher performance to price ratio than
nodes optimized for older GROMACS versions. Moreover, the shift towards GPU
processing allows to cheaply upgrade old nodes with recent GPUs, yielding
essentially the same performance as comparable brand-new hardware.


  This paper reports on the design of wildlife crossing structures (WCSs) along
a new expressway in China, which exemplifies the country's increasing efforts
on wildlife protection in infrastructure projects. The expert knowledge and
field surveys were used to determine the target species in the study area and
the quantity, locations, size, and type of the WCSs. The results on relative
abundance index and encounter rate showed that the ibex (\textit{Capra ibex}),
argali sheep (Ovis ammon), and goitered gazelle (Gazella subgutturosa) are the
main ungulates in the study area. Among them, the goitered gazelle is the most
widely distributed species. WCSs were proposed based on the estimated crossing
hotspots. The mean deviation distance between those hotspots and their nearest
proposed WCSs is around 341m. In addition, those 16 proposed underpass WCSs
have a width of no less than 12m and height of no lower than 3.5m, which is
believed to be sufficient for ungulates in the area. Given the limited
availability of high-resolution movement data and wildlife-vehicle collision
data during road's early design stage, the approach demonstrated in this paper
facilitates practical spatial planning and provides insights into designing
WCSs in a desert landscape.


  Many infectious diseases can lead to re-infection. We examined the
relationship between the prevalence of repeat infection and the basic
reproductive number (R0). First we solved a generic, deterministic
compartmental model of re-infection to derive an analytic solution for the
relationship. We then numerically solved a disease specific model of syphilis
transmission that explicitly tracked re-infection. We derived a generic
expression that reflects a non-linear and monotonically increasing relationship
between proportion re-infection and R0 and which is attenuated by entry/exit
rates and recovery (i.e. treatment). Numerical simulations from the syphilis
model aligned with the analytic relationship. Re-infection proportions could be
used to understand how far regions are from epidemic control, and should be
included as a routine indicator in infectious disease surveillance.


  Controlling recurrent infectious diseases is a vital yet complicated problem.
In this paper, we propose a novel active screening model (ACTS) and algorithms
to facilitate active screening for recurrent diseases (no permanent immunity)
under infection uncertainty. Our contributions are: (1) A new approach to
modeling multi-round network-based screening/contact tracing under uncertainty,
which is a common real-life practice in a variety of diseases; (2) Two novel
algorithms, Full- and Fast-REMEDY. Full-REMEDY considers the effect of future
actions and finds a policy that provides high solution quality, where
Fast-REMEDY scales linearly in the size of the network; (3) We evaluate Full-
and Fast-REMEDY on several real-world datasets which emulate human contact and
find that they control diseases better than the baselines. To the best of our
knowledge, this is the first work on multi-round active screening with
uncertainty for diseases with no permanent immunity.


  Ecological resilience refers to the ability of a system to retain its state
when subject to state variables perturbations or parameter changes. While
understanding and quantifying resilience is crucial to anticipate the possible
regime shifts, characterizing the influence of the system parameters on
resilience is the first step towards controlling the system to avoid
undesirable critical transitions. In this paper, we apply tools of qualitative
theory of differential equations to study the resilience of competing
populations as modeled by the classical Lotka-Volterra system. Within the high
interspecific competition regime, such model exhibits bistability, and the
boundary between the basins of attraction corresponding to exclusive survival
of each population is the stable manifold of a saddle-point. Studying such
manifold and its behavior in terms of the model parameters, we characterized
the populations resilience: while increasing competitiveness leads to higher
resilience, it is not always the case with respect to reproduction. Within a
pioneering context where both populations initiate with few individuals,
increasing reproduction leads to an increase in resilience; however, within an
environment previously dominated by one population and then invaded by the
other, an increase in resilience is obtained by decreasing the reproduction
rate. Besides providing interesting insights for the dynamics of competing
population, this work brings near to each other the theoretical concepts of
ecological resilience and the mathematical methods of differential equations
and stimulates the development and application of new mathematical tools for
ecological resilience.


  Bursting and non-renewal processes are common phenomena in birth-death
process, yet no theory can quantitatively describe a non-renewal birth process
with bursting. Here, we present a theoretical model that yields the product
number counting statistics of product creation occurring in bursts and of a
non-renewal creation process. When product creation is a stationary process,
our model confirms that product number fluctuation decreases with an increase
in the product lifetime fluctuation, originating from the non-Poisson
degradation dynamics, a result obtained in previous work. Our model
additionally demonstrates that the dependence of product number fluctuation on
product lifetime fluctuation varies with time, when product creation is a
non-stationary process. We find that bursting increases product number
fluctuation, compared to birth-processes without bursting. At time zero, in a
burst-less birth process, product number fluctuation is unsurprisingly found to
be zero, but we discover that, in a bulk creation process characterized by
bursting, product number fluctuation is a finite value at time zero. The
analytic expressions we obtain are applicable to many fields related to the
study system population, such as queueing models and gene expression.


  In many chemical reactions, reaction rate fluctuation is inevitable. Reaction
rates are different whenever chemical reaction occurs due to their dependence
on the number of reaction events or the product number. As such, understanding
the impact of rate fluctuation on product number counting statistics is of the
utmost importance when developing a quantitative explanation of chemical
reactions. In this work, we present a master equation that describes reaction
rates as a function of product number and time. Our equal reveals the
relationship between the reaction rate and product number fluctuation. Product
number counting statistics uncovers a stochastic property of the product
number; product number directly manipulates the reaction rate. Specifically, we
find that product number shows super-Poisson characteristics when the product
number increases, inducing an increase in the reaction rate. While, on the
other hand, when the product number shows sub-Poisson characteristics with an
increase in the product number, this is induced by a decrease in the reaction
rate. Furthermore, our analysis exploits reaction rate fluctuation, enabling
the quantification of the deviation of an elementary reaction process from a
renewal process.


  Existing strategies for determining the optimal treatment or monitoring
strategy typically assume unlimited access to resources. However, when a health
system has resource constraints, such as limited funds, access to medication,
or monitoring capabilities, medical decisions must balance impacts on both
individual and population health outcomes. That is, decisions should account
for competition between individuals in resource usage. One simple solution is
to estimate the (counterfactual) resource usage under the possible
interventions and choose the optimal strategy for which resource usage is
within acceptable limits. We propose a method to identify the optimal dynamic
intervention strategy that leads to the best expected health outcome accounting
for a health system's resource constraints. We then apply this method to
determine the optimal dynamic monitoring strategy for people living with HIV
when resource limits on monitoring exist using observational data from the
HIV-CAUSAL Collaboration.


  We present a new conjecture about a necessary condition that a (bio)chemical
network has to satisfy for it to exhibit multistationarity. According to a
Theorem of Feliu and Wiuf [27, 12], the conjecture is known for strictly
monotonic kinetics. We give several examples illustrating our conjecture.


  $\Lambda$-Wright--Fisher processes provide a robust framework to describe the
type-frequency evolution of an infinite neutral population. We add a polynomial
drift to the corresponding stochastic differential equation to incorporate
frequency-dependent selection. A decomposition of the drift allows us to
approximate the solution of the stochastic differential equation by a sequence
of Moran models. The genealogical structure underlying the Moran model leads in
the large population limit to a generalisation of the ancestral selection graph
of Krone and Neuhauser. Building on this object, we construct a continuous-time
Markov chain and relate it to the forward process via a new form of duality,
which we call Bernstein duality. We adapt classical methods based on the moment
duality to determine the time to absorption and criteria for the accessibility
of the boundaries; this extends a recent result by Gonz\'alez Casanova and
Span\`o. An intriguing feature of the construction is that the same forward
process is compatible with multiple backward models. In this context we
introduce suitable notions for minimality among the ancestral processes and
characterise the corresponding parameter sets. In this way we recover classic
ancestral structures as minimal ones.


  By constructing chemical reaction networks (CRNs), this paper proposes a
method of synthesizing polar decoder using belief propagation (BP) algorithm
and successive cancellation (SC) algorithm, respectively. Theoretical analysis
and simulation results have validated the feasibility of the method. Reactions
in the proposed design could be experimentally implemented with DNA strand
displacement reactions, making the proposed polar decoders promising for wide
application in nanoscale devices.


  Living organisms need to be sensitive to a changing environment while also
ignoring uninformative environmental fluctuations. Here, we show that the
circadian clock in \textit{Synechococcus elongatus} can naturally tune its
environmental sensitivity, through a clock-metabolism coupling quantified in
recent experiments. The metabolic coupling can detect mismatch between clock
predictions and the day-night light cycle, and temporarily raise the clock's
sensitivity to light changes and thus entrain faster. We also analyze analogous
behavior in recent experiments on switching between slow and fast osmotic
stress response pathways in yeast. In both cases, cells can raise their
sensitivity to new external information in epochs of frequent challenging
stress, much like a Kalman filter with adaptive gain in signal processing. Our
work suggests a new class of experiments that probe the history-dependence of
environmental sensitivity in biophysical sensing mechanisms.


  The original Hardenberg's model of biomass patterns in arid and semi-arid
regions is revisited to extend it to more general non flat regions. It is
proposed a technique to study these more generalized (non-flat) regions using
both a conservation criterion and a explicit spatial dependent function $\nu
(x)$. In this paper a study of dynamical stability around system's fixed points
made. Under the idea of predictability via air images a fitted relationship
among dynamical variables at stable fixed points is stablished. Also, is
presented a discrete version of the model, in the form of Cellular Automata
techniques, that allows to neglect the spatial scale and reproduces realistic
stable spatial patterns.


  The primary building block of the body is collagen, which is found in the
extracellular matrix and in many stress-bearing tissues such as tendon and
cartilage. It provides elasticity and support to cells and tissues while
influencing biological pathways including cell signaling, motility and
differentiation. Collagen's unique triple helical structure is thought to
impart mechanical stability. However, detailed experimental studies on its
molecular mechanics have been only recently emerging. Here, we review the
treatment of the triple helix as a homogeneous flexible rod, including bend
(standard worm-like chain model), twist, and stretch deformations, and the
assumption of backbone linearity. Additionally, we discuss protein-specific
properties of the triple helix including sequence dependence, and relate
single-molecule mechanics to collagen's physiological context.


  AGel amyloidosis, formerly known as familial amyloidosis of the Finnish-type,
is caused by pathological aggregation of proteolytic fragments of plasma
gelsolin. So far, four mutations in the gelsolin gene have been reported as
responsible for the disease. Although D187N is the first identified variant and
the best characterized, its structure has been hitherto elusive. Exploiting a
recently-developed nanobody targeting gelsolin, we were able to stabilize the
G2 domain of the D187N protein and obtained, for the first time, its
high-resolution crystal structure. In the nanobody-stabilized conformation, the
main effect of the D187N substitution is the impairment of the calcium binding
capability, leading to a destabilization of the C-terminal tail of G2. However,
molecular dynamics simulations show that in the absence of the nanobody,
D187N-mutated G2 further misfolds, ultimately exposing its hydrophobic core and
the furin cleavage site. The nanobody's protective effect is based on the
enhancement of the thermodynamic stability of different G2 mutants (D187N,
G167R and N184K). In particular, the nanobody reduces the flexibility of
dynamic stretches, and most notably decreases the conformational entropy of the
C-terminal tail, otherwise stabilized by the presence of the Ca2+ ion. A
Caenorhabditis elegans-based assay was also applied to quantify the proteotoxic
potential of the mutants and determine whether nanobody stabilization
translates into a biologically relevant effect. Successful protection from G2
toxicity in vivo points to the use of C. elegans as a tool for investigating
the mechanisms underlying AGel amyloidosis and rapidly screen new therapeutics.


  It has been suggested that the fundamental haploid-diploid cycle of
eukaryotic sex exploits a rudimentary form of the Baldwin effect. This paper
uses the well-known NKCS model to explore the effects of coevolution upon the
behaviour of eukaryotes. It is shown how varying fitness landscape size,
ruggedness and connectedness can vary the conditions under which eukaryotic sex
proves beneficial over asexual reproduction in haploids in a coevolutionary
context. Moreover, eukaryotic sex is shown to be more sensitive to the relative
rate of evolution exhibited by its partnering species than asexual haploids.


  Network-based interventions against epidemic spread are most powerful when
the full network structure is known. However, in practice, resource constraints
require decisions to be made based on partial network information. We
investigated how the accuracy of network data available at individual and
village levels affected network-based vaccination effectiveness. We simulated a
Susceptible-Infected-Recovered process on empirical social networks from 75
villages. First, we used regression to predict the percentage of individuals
ever infected based on village-level network. Second, we simulated vaccinating
10 percent of each of the 75 empirical village networks at baseline, selecting
vaccinees through one of five network-based approaches: random individuals;
random contacts of random individuals; random high-degree individuals; highest
degree individuals; or most central individuals. The first three approaches
require only sample data; the latter two require full network data. We also
simulated imposing a limit on how many contacts an individual can nominate
(Fixed Choice Design, FCD), which reduces the data collection burden but
generates only partially observed networks. We found mean and standard
deviation of the degree distribution to strongly predict cumulative incidence.
In simulations, the Nomination method reduced cumulative incidence by one-sixth
compared to Random vaccination; full network methods reduced infection by
two-thirds. The High Degree approach had intermediate effectiveness.
Surprisingly, FCD truncating individuals' degrees at three was as effective as
using complete networks. Using even partial network information to prioritize
vaccines at either the village or individual level substantially improved
epidemic outcomes. Such approaches may be feasible and effective in outbreak
settings, and full ascertainment of network structure may not be required.


  In this paper, we compared cancerous and normal cell according to their
protein-protein interaction network. Cancer is one of the complicated diseases
and experimental investigations have been showed that protein interactions have
an important role in the growth of cancer. We calculated some graph related
parameters such as Number of Vertices, Number of Edges, Closeness, Graph
Diameter, Graph Radius, Index of Aggregation, Connectivity, Number of Edges
divided by the Number of Vertices, Degree, Cluster Coefficient, Subgraph
Centrality, and Betweenness. Furthermore, the number of motifs and hubs in
these networks have been measured. In this paper bone, breast, colon, kidney
and liver benchmark datasets have been used for experiments. The experimental
results show that Graph Degree Mean, Subgraph Centrality, Betweenness, and Hubs
have higher values in the cancer cells and can be used as a measure to
distinguish between normal and cancerous networks. The cancerous tissues of the
five studied samples are denser in the interaction networks.


  Emerging research suggests that the extent to which activity spaces -- the
collection of an individual's routine activity locations -- overlap provides
important information about the functioning of a city and its neighborhoods. To
study patterns of overlapping activity spaces, we draw on the notion of an
ecological network, a type of two-mode network with the two modes being
individuals and the geographic locations where individuals perform routine
activities. We describe a method for detecting "ecological communities" within
these networks based on shared activity locations among individuals.
Specifically, we identify latent activity pattern profiles, which, for each
community, summarize its members' probability distribution of going to each
location, and community assignment vectors, which, for each individual,
summarize his/her probability distribution of belonging to each community.
Using data from the Adolescent Health and Development in Context (AHDC) Study,
we employ latent Dirichlet allocation (LDA) to identify activity pattern
profiles and communities. We then explore differences across neighborhoods in
the strength, and within-neighborhood consistency of community assignment. We
hypothesize that these aspects of the neighborhood structure of ecological
community membership capture meaningful dimensions of neighborhood functioning
likely to co-vary with economic and racial composition. We discuss the
implications of a focus on ecological communities for the conduct of
"neighborhood effects" research more broadly.


  To understand the effects of fluctuations on achieving homochirality, we
employ a Monte-Carlo method where autocatalysis and enantiomeric
cross-inhibition, as well as racemization and deracemization reactions are
included. The results of earlier work either without autocatalysis or without
cross-inhibition are reproduced. Bifurcation diagrams and the dependencies of
the number of reaction steps on parameters are studied. In systems with 30,000
molecules, for example, up to a billion reaction steps may be needed to achieve
homochirality without autocatalysis.


  Reciprocal best matches play an important role in numerous applications in
computational biology, in particular as the basis of many widely used tools for
orthology assessment. Nevertheless, very little is known about their
mathematical structure. Here, we investigate the structure of reciprocal best
match graphs (RBMGs). In order to abstract from the details of measuring
distances, we define reciprocal best matches here as pairwise most closely
related leaves in a gene tree, arguing that conceptually this is the notion
that is pragmatically approximated by distance- or similarity-based heuristics.
We start by showing that a graph $G$ is an RBMG if and only if its quotient
graph w.r.t.\ a certain thinness relation is an RBMG. Furthermore, it is
necessary and sufficient that all connected components of $G$ are RBMGs. The
main result of this contribution is a complete characterization of RBMGs with 3
colors/species that can be checked in polynomial time. For 3 colors, there are
three distinct classes of trees that are related to the structure of the
phylogenetic trees explaining them. We derive an approach to recognize RBMGs
with an arbitrary number of colors; it remains open however, whether a
polynomial-time for RBMG recognition exists. In addition, we show that RBMGs
that at the same time are cographs (co-RBMGs) can be recognized in polynomial
time. Co-RBMGs are characterized in terms of hierarchically colored cographs, a
particular class of vertex colored cographs that is introduced here. The (least
resolved) trees that explain co-RBMGs can be constructed in polynomial time.


  When analysing in vitro data, growth kinetics of influenza strains are often
compared by computing their growth rates, which are sometimes used as proxies
for fitness. However, analogous to mechanistic epidemic models, the growth rate
can be defined as a function of two parameters: the basic reproduction number
(the average number of cells each infected cell infects) and the mean
generation time (the average length of a replication cycle). Using a
mechanistic model, previously published data from experiments in human lung
cells, and newly generated data, we compared estimates of all three parameters
for six influenza A strains. Using previously published data, we found that the
two human-adapted strains (pre-2009 seasonal H1N1, and pandemic H1N1) had a
lower basic reproduction number, shorter mean generation time and slower growth
rate than the two avian-adapted strains (H5N1 and H7N9). These same differences
were then observed in data from new experiments where two strains were
engineered to have different internal proteins (pandemic H1N1 and H5N1), but
the same surface proteins (PR8), confirming our initial findings and implying
that differences between strains were driven by internal genes. Also, the model
predicted that the human-adapted strains underwent more replication cycles than
the avian-adapted strains by the time of peak viral load, potentially
accumulating mutations more quickly. These results suggest that the in vitro
reproduction number, generation time and growth rate differ between
human-adapted and avian-adapted influenza strains, and thus could be used to
assess host adaptation of internal proteins to inform pandemic risk assessment.


  Principal component analysis (PCA) represents a standard approach to identify
collective variables $\{x_i\}\!=\!\boldsymbol{x}$, which can be used to
construct the free energy landscape $\Delta G(\boldsymbol{x})$ of a molecular
system. While PCA is routinely applied to equilibrium molecular dynamics (MD)
simulations, it is less obvious how to extend the approach to nonequilibrium
simulation techniques. This includes, e.g., the definition of the statistical
averages employed in PCA, as well as the relation between the equilibrium free
energy landscape $\Delta G(\boldsymbol{x})$ and energy landscapes $\Delta{\cal
G} (\boldsymbol{x})$ obtained from nonequilibrium MD. As an example for a
nonequilibrium method, `targeted MD' is considered which employs a moving
distance constraint to enforce rare transitions along some biasing coordinate
$s$. The introduced bias can be described by a weighting function $P(s)$, which
provides a direct relation between equilibrium and nonequilibrium data, and
thus establishes a well-defined way to perform PCA on nonequilibrium data.
While the resulting distribution ${\cal P}(\boldsymbol{x})$ and energy
$\Delta{\cal G} \propto \ln {\cal P}$ will not reflect the equilibrium state of
the system, the nonequilibrium energy landscape $\Delta{\cal G}
(\boldsymbol{x})$ may directly reveal the molecular reaction mechanism. Applied
to targeted MD simulations of the unfolding of decaalanine, for example, a PCA
performed on backbone dihedral angles is shown to discriminate several
unfolding pathways. Although the formulation is in principle exact, its
practical use depends critically on the choice of the biasing coordinate $s$,
which should account for a naturally occurring motion between two well-defined
end-states of the system.


  In biology, there are several processes in which unfolded protein chains are
transported along narrow-tube channels. Normally, without such a severe
configurational constraint, unfolded polypeptides would not bind to each other.
However, when chain entropy is much reduced in the narrow channel, we find that
polypeptide chains have a propensity to bind, even if there is no great
potential energy gain in doing so. We find the average length of binding m*
(the number of monomers at the chain ends that form bonds) and the critical
tube diameter at which such constrained binding occurs. We carry out Brownian
dynamics simulations of tightly confined chains, demonstrating their binding
over the characteristic length m*, changing in tubes of different diameter.


  In this study, we describe different modeling approaches for ovarian follicle
population dynamics, based on either ordinary (ODE), partial (PDE) or
stochastic (SDE) differential equations, and accounting for interactions
between follicles. We put a special focus on representing the population-level
feedback exerted by growing ovarian follicles onto the activation of quiescent
follicles. We take advantage of the timescale difference existing between the
growth and activation processes to apply model reduction techniques in the
framework of singular perturbations. We first study the linear versions of the
models to derive theoretical results on the convergence to the limit models. In
the nonlinear cases, we provide detailed numerical evidence of convergence to
the limit behavior. We reproduce the main semi-quantitative features
characterizing the ovarian follicle pool, namely a bimodal distribution of the
whole population, and a slope break in the decay of the quiescent pool with
aging.


  Various kinetic Monte Carlo algorithms become inefficient when some of the
population sizes in a system are large, which gives rise to a large number of
reaction events per unit time. Here, we present a new acceleration algorithm
based on adaptive and heterogeneous scaling of reaction rates and
stoichiometric coefficients. The algorithm is conceptually related to the
commonly used idea of accelerating a stochastic simulation by considering a
sub-volume $\lambda \Omega$ ($0<\lambda<1$) within a system of interest, which
reduces the number of reaction events per unit time occurring in a simulation
by a factor $1/\lambda$ at the cost of greater error in unbiased estimates of
first moments and biased overestimates of second moments. Our new approach
offers two unique benefits. First, scaling is adaptive and heterogeneous, which
eliminates the pitfall of overaggressive scaling. Second, there is no need for
an \emph{a priori} classification of populations as discrete or continuous (as
in a hybrid method), which is problematic when discreteness of a chemical
species changes during a simulation. The method requires specification of only
a single algorithmic parameter, $N_c$, a global critical population size above
which populations are effectively scaled down to increase simulation
efficiency. The method, which we term partial scaling, is implemented in the
open-source BioNetGen software package. We demonstrate that partial scaling can
significantly accelerate simulations without significant loss of accuracy for
several published models of biological systems. These models characterize
activation of the mitogen-activated protein kinase ERK, prion protein
aggregation, and T-cell receptor signaling.


  The representation problem of finite-dimensional Markov matrices in Markov
semigroups is revisited, with emphasis on concrete criteria for matrix
subclasses of theoretical or practical relevance, such as equal-input,
circulant, symmetric or doubly stochastic matrices. Here, we pay special
attention to various algebraic properties of the embedding problem, and discuss
the connection with the centraliser of a Markov matrix.


  We study a variant of the voter model with multiple opinions; individuals can
imitate each other and also change their opinion randomly in mutation events.
We focus on the case of a population with all-to-all interaction. A
noise-driven transition between regimes with multi-modal and unimodal
stationary distributions is observed. In the former, the population is mostly
in consensus states; in the latter opinions are mixed. We derive an effective
death-birth process, describing the dynamics from the perspective of one of the
opinions, and use it to analytically compute marginals of the stationary
distribution. These calculations are exact for models with homogeneous
imitation and mutation rates, and an approximation if rates are heterogeneous.
Our approach can be used to characterize the noise-driven transition and to
obtain mean switching times between consensus states.


  Morphogenesis is the biological process that causes the emergence and changes
of patterns (tissues and organs) in living organisms. It is a robust,
self-organising mechanism, governed by Genetic Regulatory Networks (GRN), that
hasn't been thoroughly understood. In this work we propose Epigenetic Forests
as a tool to study morphogenesis and to extract valuable information from GRN.
Our method unfolds the richness and structure within the GRN.
  As a case study, we analyze the GRN during cell fate determination during the
early stages of development of the flower Arabidopsis thaliana and its spatial
dynamics. By using a genetic algorithm we optimize cell differentiation in our
model and correctly recover the architecture of the flower.


  Darwin's theory of evolution by natural selection does not predict long-term
progress or advancement, nor does it provide a useful way to define or
understand these concepts. Nevertheless, the history of life is marked by major
trends that appear progressive, and seemingly more advanced forms of life have
appeared. To reconcile theory and fact, evolutionists have proposed novel
theories that extend natural selection to levels and time frames not justified
by the original structure of Darwin's theory. To extend evolutionary theory
without violating the most basic tenets of Darwinism, I here identify a
separate struggle and an alternative evolutionary force. Owing to the abundant
free energy in our universe, there is a struggle for supremacy that naturally
rewards those that are first to invent novelties that allow exploitation of
untapped resources. This natural reward comes in form of a temporary monopoly,
which is granted to those who win a competitive race to innovate. By analogy to
human economies, natural selection plays the role of nature's inventor,
gradually fashioning inventions to the situation at hand, while natural reward
plays the role of nature's entrepreneur, choosing which inventions to first
disseminate to large markets. Natural reward leads to progress through a
process of invention-conquest macroevolution, in which the dual forces of
natural selection and natural reward create and disseminate major innovations.
Over vast time frames, natural reward drives the advancement of life by a
process of extinction-replacement megaevolution that releases constraints on
progress and increases the innovativeness of life.


  The emergence and maintenance of cooperation has attracted intensive
scholarly interest and has been analysed within the framework of evolutionary
game theory. The role of innovation, which introduces novel strategies into the
population, is a relatively understudied aspect of evolutionary game theory.
Here, we investigate the effects of two sources of innovation---mutation and
heterogeneous updating rules. These mechanisms allow agents to adopt strategies
that do not rely on the imitation of other individuals. The model
introduces---in addition to canonical imitation-based strategy
updating---aspiration-based updating, whereby agents switch their strategy by
referring solely to the performance of their own strategy; mutation also
introduces novel strategies into the population. Our simulation results show
that the introduction of aspiration-based rules into a population of imitators
leads to the deterioration of cooperation. In addition, mutation, in
combination with heterogeneous updating rules, also diminishes cooperators.
This phenomenon is prominent when a large proportion of the population consists
of imitators rather than adopters of aspiration-based updating. Nevertheless, a
high mutation rate, in combination with a low aspiration level, has positive
nonlinear effects, and a heterogeneous population achieves a higher level of
cooperation than the weighted average of homogeneous populations. Our results
demonstrate the profound role of innovation in the evolution of cooperation.


  We suggest a natural approach that leads to a modification of classical
quasispecies models and incorporates the possibility of population extinction
in addition to growth. The resulting modified models are called open. Their
essential properties, regarding in particular equilibrium behavior, are
investigated both analytically and numerically. The hallmarks of the
quasispecies dynamics, viz. the heterogeneous quasispecies distribution itself
and the error threshold phenomenon, can be observed in our models, along with
extinction. In order to demonstrate the flexibility of the introduced
framework, we study the inverse problem of fitness allocation under the
biologically motivated criterion of steady-state fitness maximization. Having
in mind the complexity of numerical investigation of high-dimensional
quasispecies problems and the fact that the actual number of genotypes or
alleles involved in a studied process can be extremely large, we also build
continuous-time distributed open quasispecies models. The obtained results may
serve as an initial step to developing mathematical models that involve
directed therapy against various pathogens.


  Although non-canonical residues, caps, crosslinks, and nicks play an
important role in the function of many DNA, RNA, proteins, and complexes, we do
not fully understand how networks of non-canonical macromolecules generate
behavior. One barrier is our limited formats, such as IUPAC, for abstractly
describing macromolecules. To overcome this barrier, we developed BpForms and
BcForms, a toolkit of ontologies, grammars, and software for abstracting the
primary structure of polymers and complexes as combinations of residues, caps,
crosslinks, and nicks. The toolkit can help quality control, exchange, and
integrate information about the primary structure of macromolecules into
fine-grained global networks of intracellular biochemistry.


  A model simulating transport of dense core vesicles (DCVs) in type II axonal
terminals of Drosophila motoneurons has been developed. The morphology of type
II terminals is characterized by the large number of en passant boutons. The
lack of both scaled up DCV transport and scaled down DCV capture in boutons
results in a less efficient supply of DCVs to distal boutons. Furthermore, the
large number of boutons that DCVs pass as they move anterogradely, until they
reach the most distal bouton, may lead to the capture of a majority of DCVs
before they turn around in the most distal bouton to move in the retrograde
direction. This may lead to a reduced retrograde flux of DCVs and a lack of DCV
circulation in type II terminals. The developed model simulates DCV
concentrations in boutons, DCV fluxes between the boutons, age density
distributions of DCVs, and the mean age of DCVs in various boutons. Unlike
published experimental observations, our model predicts DCV circulation in type
II terminals after these terminals are filled to saturation. This disagreement
is likely because experimentally observed terminals were not at steady-state,
but rather were accumulating DCVs for later release. Our estimates show that
the number of DCVs in the transiting state is much smaller than that in the
resident state. DCVs traveling in the axon, rather than DCVs transiting in the
terminal, may provide a reserve of DCVs for replenishing boutons after a
release. The techniques for modeling transport of DCVs developed in our paper
can be used to model the transport of other organelles in axons.


  The circadian clock is the molecular mechanism responsible for the adaptation
to daily rhythms in living organisms. Oscillations and fluctuations in
environmental conditions regulate the circadian clock through signaling
pathways. We study the response to continuous photic perturbations in a minimal
molecular network model of the circadian clock, composed of 5 nonlinear delay
differential equations with multiple feedbacks. We model the perturbation as a
stationary stochastic process, and we consider the resulting irreversibility of
trajectories as a key effect of the interaction. In particular we adopt a
measure of mutual mapping irreversibility in the time series thermodynamics
framework, and we find 12 hours harmonics.


  Ecological connectivity in coastal oceanic waters is mediated by dispersion
of the early life stages of marine organisms and conditions the structure of
biological communities and the provision of ecosystem services. Integrated
management strategies aimed at ensuring long-term service provision to society
do not currently consider the importance of dispersal and larval connectivity.
A spatial optimization model is introduced to maximise the potential provision
of ecosystem services in coastal areas by accounting for the role of dispersal
and larval connectivity. The approach combines a validated coastal circulation
model that reproduces realistic patterns of larval transport along the coast,
which ultimately conditions the biological connectivity and productivity of an
area, with additional spatial layers describing potential ecosystem services.
The spatial optimization exercise was tested along the coast of Central Chile,
a highly productive area dominated by the Humboldt Current. Results show it is
unnecessary to relocate existing management areas, as increasing no-take areas
by 10% could maximise ecosystem service provision, while improving the spatial
representativeness of protected areas and minimizing social conflicts. The
location of protected areas was underrepresented in some sections of the study
domain, principally due to the restriction of the model to rocky subtidal
habitats. Future model developments should encompass the diversity of coastal
ecosystems and human activities to inform integrative spatial management.
Nevertheless, the spatial optimization model is innovative not only for its
integrated ecosystem perspective, but also because it demonstrates that it is
possible to incorporate time-varying biophysical connectivity within the
optimization problem, thereby linking the dynamics of exploited populations
produced by the spatial management regime.


  In this paper, we provide a straightforward approach to defining and deriving
the key epidemiological quantity, the basic reproduction number, $R_0$, for
Markovian epidemics in structured populations. The methodology derived is
applicable to, and demonstrated on, both $SIR$ and $SIS$ epidemics and allows
for population as well as epidemic dynamics. The approach taken is to consider
the epidemic process as a multitype process by identifying and classifying the
different types of infectious units along with the infections from, and the
transitions between, infectious units. For the household model, we show that
our expression for $R_0$ agrees with earlier work despite the alternative
nature of the construction of the mean reproductive matrix, and hence, the
basic reproduction number.


  Cell competition is a quality control mechanism in tissues that results in
the elimination of less fit cells. Over the past decade, the phenomenon of cell
competition has been identified in many physiological and pathological
contexts, driven either by biochemical signaling or by mechanical forces within
the tissue. In both cases, competition has generally been characterized based
on the elimination of loser cells at the population level, but significantly
less attention has been focused on determining how single-cell dynamics and
interactions regulate population-wide changes. In this review, we describe
quantitative strategies and outline the outstanding challenges in understanding
the single cell rules governing tissue-scale competition dynamics. We propose
quantitative metrics to characterize single cell behaviors in competition and
use them to distinguish the types and outcomes of competition. We describe how
such metrics can be measured experimentally using a novel combination of
high-throughput imaging and machine learning algorithms. We outline the
experimental challenges to quantify cell fate dynamics with high-statistical
precision, and describe the utility of computational modeling in testing
hypotheses not easily accessible in experiments. In particular, cell-based
modeling approaches that combine mechanical interaction of cells with
decision-making rules for cell fate choices provide a powerful framework to
understand and reverse-engineer the diverse rules of cell competition.


  Population structure affects the outcome of natural selection. Static
population structures can be described by graphs, where individuals occupy the
nodes, and interactions occur along the edges. General conditions for
evolutionary success on any weighted graph were recently derived, for weak
selection, in terms of coalescence times of random walks. Here we show that for
a special class of graphs, the conditions for success take a particularly
simple form, in which all effects of graph structure are described by the
graph's "effective degree"---a measure of the effective number of neighbors per
individual. This result holds for all weighted graphs that are isothermal,
meaning that the sum of edge weights is the same at each node. Isothermal
graphs encompass a wide variety of underlying topologies, and arise naturally
from supposing that each individual devotes the same amount of time to
interaction. Cooperative behavior is favored on a large isothermal graph if the
benefit-to-cost ratio exceeds the effective degree. We relate the effective
degree of a graph to its spectral gap, thereby providing a link between
evolutionary dynamics and the theory of expander graphs. As a surprising
example, we report graphs of infinite average degree that are nonetheless
highly conducive for promoting cooperation.


  Influenza and respiratory syncytial virus (RSV) are the leading etiological
agents of seasonal acute respiratory infections (ARI) around the world. Medical
doctors typically base the diagnosis of ARI on patients' symptoms alone and do
not always conduct virological tests necessary to identify individual viruses,
which limits the ability to study the interaction between multiple pathogens
and make public health recommendations. We consider a stochastic kinetic model
(SKM) for two interacting ARI pathogens circulating in a large population and
an empirically motivated background process for infections with other pathogens
causing similar symptoms. An extended marginal sampling approach based on the
Linear Noise Approximation to the SKM integrates multiple data sources and
additional model components. We infer the parameters defining the pathogens'
dynamics and interaction within a Bayesian hierarchical model and explore the
posterior trajectories of infections for each illness based on aggregate
infection reports from six epidemic seasons collected by the state health
department, and a subset of virological tests from a sentinel program at a
general hospital in San Luis Potos\'i, M\'exico. We interpret the results based
on real and simulated data and make recommendations for future data collection
strategies. Supplementary materials and software are provided online.


  Angiotensin II (AngII) type 1 receptor (AT1R) is a G protein-coupled receptor
known for its role in numerous physiological processes and its implication in
many vascular diseases. Its functions are mediated through G protein dependent
and independent signaling pathways. AT1R has several endogenous peptidic
agonists, all derived from angiotensinogen, as well as several synthetic
ligands known to elicit biased signaling responses. Here, surface plasmon
resonance (SPR) was used as a cell- based and label-free technique to quantify,
in real time, the response of HEK293 cells stably expressing the human AT1R.
The goal was to take advantage of the integrative nature of this assay to
identify specific signaling pathways in the features of the response profiles
generated by numerous endogenous and synthetic ligands of AT1R. First, we
assessed the contributions of Gq, G12/13, Gi, Gbg, ERK1/2 and \b{eta}-
arrestins pathways in the cellular responses measured by SPR where Gq,
G12/Rho/ROCK together with \b{eta}-arrestins and ERK1/2 were found to play
significant roles. More specifically, we established a major role for G12 in
the early events of the AT1R-dependent response, which was followed by a robust
ERK1/2 component associated to the later phase of the signal. Interestingly,
endogenous AT1R ligands (AngII, AngIII and AngIV) exhibited distinct responses
signatures with a significant increase of the ERK1/2-like components for both
AngIII and AngIV, which points toward possibly distinct physiological roles for
the later. We also tested AT1R biased ligands, all of which affected both the
early and later events. Our results support SPR-based integrative cellular
assays as a powerful approach to delineate the contribution of specific
signaling pathways for a given cell response and reveal response differences
associated with ligands with distinct pharmacological properties.


  G protein coupled receptors (GPCRs) produce pleiotropic effects by their
capacity to engage numerous signaling pathways once activated. Functional
selectivity (also called biased signaling), where specific compounds can bring
GPCRs to adopt conformations that enable selective receptor coupling to
distinct signaling pathways, continues to be significantly investigated.
However, an important but often overlooked aspect of functional selectivity is
the capability of ligands such as angiotensin II (AngII) to adopt specific
conformations that may preferentially bind to selective GPCRs structures.
Understanding both receptor and ligand conformation is of the utmost importance
for the design of new drugs targeting GPCRs. In this study, we examined the
properties of AngII cyclic analogs to impart biased agonism on the angiotensin
type 1 receptor (AT1R). Positions 3 and 5 of AngII were substituted for
cysteine and homocysteine residues ([Sar1Hcy3,5]AngII, [Sar1Cys3Hcy5]AngII and
[Sar1Cys3,5]AngII) and the resulting analogs were evaluated for their capacity
to activate the Gq/11, G12, Gi2, Gi3, Gz, ERK and \b{eta}-arrestin (\b{eta}arr)
signaling pathways via AT1R. Interestingly, [Sar1Hcy3,5]AngII exhibited potency
and full efficacy on all pathways tested with the exception of the Gq pathway.
Molecular dynamic simulations showed that the energy barrier associated with
the insertion of residue Phe8 of AngII within the hydrophobic core of AT1R,
associated with Gq/11 activation, is increased with [Sar1Hcy3,5]AngII. These
results suggest that constraining the movements of molecular determinants
within a given ligand by introducing cyclic structures may lead to the
generation of novel ligands providing more efficient biased agonism.


  Recently, persistent homology has had tremendous success in biomolecular data
analysis. It works by examining the topological relationship or connectivity of
a group of atoms in a molecule at a variety of scales, then rendering a family
of topological representations of the molecule. However, persistent homology is
rarely {employed} for the analysis of atomic properties, such as biomolecular
flexibility analysis or B factor prediction. This work introduces atom-specific
persistent homology to provide a local atomic level representation of a
molecule via a global topological tool. This is achieved through the
construction of a pair of conjugated sets of atoms and corresponding conjugated
simplicial complexes, as well as conjugated topological spaces. The difference
between the topological invariants of the pair of conjugated sets is measured
by Bottleneck and Wasserstein metrics and leads to an atom-specific topological
representation of individual atomic properties in a molecule. Atom-specific
topological features are integrated with various machine learning algorithms,
including gradient boosting trees and convolutional neural network for protein
thermal fluctuation analysis and B factor prediction. Extensive numerical
results indicate the proposed method provides a powerful topological tool for
analyzing and predicting localized information.


  GTPases regulate a wide range of cellular processes, such as intracellular
vesicular transport, signal transduction, and protein translation. These
hydrolase enzymes operate as biochemical switches by toggling between an active
guanosine triphosphate (GTP)-bound state and an inactive guanosine diphosphate
(GDP)-bound state. We compare two network motifs, a single-species switch and
an interlinked cascade that consists of two species coupled through positive
and negative feedback loops. We find that interlinked cascades are closer to
the ideal all-or-none switch and are more robust against fluctuating signals.
While the single-species switch can only achieve bistability, interlinked
cascades can be converted into oscillators by tuning the cofactor
concentrations, which catalyse the activity of the cascade. These regimes can
only be achieved with sufficient chemical driving provided by GTP hydrolysis.
In this study, we present a thermodynamically consistent model that can achieve
bistability and oscillations with the same feedback motif.


  Social scientists have long appreciated that relationships between
individuals cannot be described from observing a single domain, and that the
structure across domains of interaction can have important effects on outcomes
of interest (e.g., cooperation).1 One debate explicitly about this surrounds
food sharing. Some argue that failing to find reciprocal food sharing means
that some process other than reciprocity must be occurring, whereas others
argue for models that allow reciprocity to span domains in the form of trade.2
Multilayer networks, high-dimensional networks that allow us to consider
multiple sets of relationships at the same time, are ubiquitous and have
consequences, so processes giving rise to them are important social phenomena.
The analysis of multi-dimensional social networks has recently garnered the
attention of the network science community.3 Recent models of these processes
show how ignoring layer interdependencies can lead one to miss why a layer
formed the way it did, and/or draw erroneous conclusions.6 Understanding the
structuring processes that underlie multiplex networks will help understand
increasingly rich datasets, giving more accurate and complete pictures of
social interactions.


  The challenge in controlling stochastic systems in which low-probability
events can set the system on catastrophic trajectories is to develop a robust
ability to respond to such events without significantly compromising the
optimality of the baseline control policy. This paper presents CelluDose, a
stochastic simulation-trained deep reinforcement learning adaptive feedback
control prototype for automated precision drug dosing targeting stochastic and
heterogeneous cell proliferation. Drug resistance can emerge from random and
variable mutations in targeted cell populations; in the absence of an
appropriate dosing policy, emergent resistant subpopulations can proliferate
and lead to treatment failure. Dynamic feedback dosage control holds promise in
combatting this phenomenon, but the application of traditional control
approaches to such systems is fraught with challenges due to the complexity of
cell dynamics, uncertainty in model parameters, and the need in medical
applications for a robust controller that can be trusted to properly handle
unexpected outcomes. Here, training on a sample biological scenario identified
single-drug and combination therapy policies that exhibit a 100% success rate
at suppressing cell proliferation and responding to diverse system
perturbations while establishing low-dose no-event baselines. These policies
were found to be highly robust to variations in a key model parameter subject
to significant uncertainty and unpredictable dynamical changes.


  Stem cell heterogeneity is essential for the homeostasis in tissue
development. This paper established a general formulation for understanding the
dynamics of stem cell regeneration with cell heterogeneity and random
transitions of epigenetic states. The model generalizes the classical G0 cell
cycle model, and incorporates the epigenetic states of stem cells that are
represented by a continuous multidimensional variable and the kinetic rates of
cell behaviors, including proliferation, differentiation, and apoptosis, that
are dependent on their epigenetic states. Moreover, the random transition of
epigenetic states is represented by an inheritance probability that can be
described as a conditional beta distribution. This model can be extended to
investigate gene mutation-induced tumor development. The proposed formula is a
generalized formula that helps us to understand various dynamic processes of
stem cell regeneration, including tissue development, degeneration, and
abnormal growth.


  The concentration of the population in coastal regions, in addition to the
direct human use, is leading to an accelerated process of change and
deterioration of the marine ecosystems. Human activities such as fishing
together with environmental drivers (e.g. climate change) are triggering major
threats to marine biodiversity, and impact directly the services they provide.
In the South and Southwest coasts of Portugal, the deep-water crustacean trawl
fishery is not exemption. This fishery is recognized to have large effects on a
number of species while generating high rates of unwanted catches. However,
taking into account an ecosystem-based perspective, the fishing impacts along
the food web accounting for biological interactions between and among species
caught remains poorly understood. These impacts are particularly troubling and
are a cause of concern given the cascading effects that might arise. Facing the
main policies and legislative instruments for the restoration and conservation
of the marine environment, times are calling for implementing ecosystem-based
approaches to fisheries management. To this end, we use a food web modelling
(Ecopath with Ecosim) approach to assess the fishing impacts of this particular
fishery on the marine ecosystem of southern and southwestern Portugal. In
particular, we describe the food web structure and functioning, identify the
main keystone species and/or groups, quantify the major trophic and energy
flows, and ultimately assess the impact of fishing on the target species but
also on the ecosystem by means of ecological and ecosystem-based indicators.
Finally, we examine limitations and weaknesses of the model for potential
improvements and future research directions.


  Following the ideas of Davydov's soliton theory, we study the bio-energy
transport in protein molecules. By using a quantum Brownian motion model for a
phonon dressed vibrational exciton, we calculate the time-dependence on the
mean square distance, diffusion coefficient, and energy of the vibrational
exciton. We find the time-dependence by solving the quantum Langevin equation
and find oscillatory behaviors due to the super-diffusive non-ohmic
dissipation. We find that the vibrational exciton gains an overall energy due
to the coupling to the phonon bath; it also dissipates its energy to the
environment as it propagates. The amount of energy gain and the oscillatory
features depend on both temperature and the phonon-vibron coupling.


  Accumulation of plastic waste is a major environmental problem. Enzymes,
particularly esterases, play an important role in the biodegradation of
polyesters. These enzymes are usually only active on aliphatic polyesters, but
a few have showed catalytic activity for semi-aromatic polyesters. Due to the
importance of these processes, an atomic level characterization of how common
polyesters are degraded by esterases is necessary. Hereby, we present a
Molecular dynamics (MD) and Quantum Mechanics/Molecular Mechanics (QM/MM) MD
study of the hydrolysis of a model of polycaprolactone (PCL), one of the most
widely used biomaterials, by the thermophilic esterase from the archaeon
Archaeoglobus fulgidus (AfEST). This enzyme is particularly interesting because
it can withstand temperatures well above the glass transition of many
polyesters. Our insights about the reaction mechanism are important for the
design of customized enzymes able to degrade different synthetic polyesters.


  Theory predicts rapid genetic drift during invasions, yet many expanding
populations maintain high genetic diversity. We find that genetic drift is
dramatically suppressed when dispersal rates increase with the population
density because many more migrants from the diverse, high-density regions
arrive at the expansion edge. When density-dependence is weak or negative, the
effective population size of the front scales only logarithmically with the
carrying capacity. The dependence, however, switches to a sublinear power law
and then to a linear increase as the density-dependence becomes strongly
positive. We develop a unified framework revealing that the transitions between
different regimes of diversity loss are controlled by a single, universal
parameter: the ratio of the expansion velocity to the geometric mean of
dispersal and growth rates at expansion edge. Our results suggest that positive
density-dependence could dramatically alter evolution in expanding populations
even when its contributions to the expansion velocity is small.


  The Colless index is one of the oldest and most widely used balance indices
for rooted bifurcating trees. Despite its popularity, its minimum value on the
space $\mathcal{T}_n$ of rooted bifurcating trees with $n$ leaves is only known
when $n$ is a power of 2. In this paper we fill this gap in the literature, by
providing a formula that computes, for each $n$, the minimum Colless index on
$\mathcal{T}_n$, and characterizing those trees where this minimum value is
reached.


  Coalescent theory combined with statistical modeling allows us to estimate
effective population size fluctuations from molecular sequences of individuals
sampled from a population of interest. When sequences are sampled serially
through time and the distribution of the sampling times depends on the
effective population size, explicit statistical modeling of sampling times
improves population size estimation. Previous work assumed that the genealogy
relating sampled sequences is known and modeled sampling times as an
inhomogeneous Poisson process with log-intensity equal to a linear function of
the log-transformed effective population size. We improve this approach in two
ways. First, we extend the method to allow for joint Bayesian estimation of the
genealogy, effective population size trajectory, and other model parameters.
Next, we improve the sampling time model by incorporating additional sources of
information in the form of time-varying covariates. We validate our new
modeling framework using a simulation study and apply our new methodology to
analyses of population dynamics of seasonal influenza and to the recent Ebola
virus outbreak in West Africa.


  The fitness landscape encodes the mapping of genotypes to fitness and
provides a succinct representation of possible trajectories followed by an
evolving population. Evolutionary accessibility is quantified by the existence
of fitness-monotonic paths connecting far away genotypes. Studies of
accessibility percolation use probabilistic fitness landscape models to explore
the emergence of such paths as a function of the initial fitness, the
parameters of the landscape or the structure of the genotype graph. This
chapter reviews these studies and discusses their implications for the
predictability of evolutionary processes.


  DNA and RNA are generally regarded as central molecules in molecular biology.
Recent advancements in the field of DNA/RNA nanotechnology successfully used
DNA/RNA as programmable molecules to construct molecular machines and
nanostructures with predefined shapes and functions. The key mechanism for
dynamic control of the conformations of these DNA/RNA nanodevices is a reaction
called strand displacement, in which one strand in a formed duplex is replaced
by a third invading strand. While DNA/RNA strand displacement has mainly been
used to de-novo design molecular devices, we argue in this review that this
reaction is also likely to play a key role in multiple cellular events such as
gene recombination, CRISPR-based genome editing, and RNA cotranscriptional
folding. We introduce the general mechanism of strand displacement reaction,
give examples of its use in the construction of molecular machines, and finally
review natural processes having characteristic which suggest that strand
displacement is occurring.


  The analysis of biological sequencing data has been one of the biggest
applications of string algorithms. The approaches used in many such
applications are based on the analysis of k-mers, which are short fixed-length
strings present in a dataset. While these approaches are rather diverse,
storing and querying a k-mer set has emerged as a shared underlying component.
A set of k-mers has unique features and applications that, over the last ten
years, have resulted in many specialized approaches for its representation. In
this survey, we give a unified presentation and comparison of the data
structures that have been proposed to store and query a k-mer set. We hope this
survey will serve as a resource for researchers in the field as well as make
the area more accessible to researchers outside the field.


  In this review, we focus on the kinesin-3 family molecular motor protein
UNC-104 and its regulatory protein ARL-8. UNC-104, originally identified in
Caenorhabditis elegans (C. elegans), has a primary role transporting synaptic
vesicle precursors (SVPs). Although in vitro single-molecule experiments have
been performed to primarily investigate the kinesin motor domain, these have
not addressed the in vivo reality of the existence of regulatory proteins, such
as ARL-8, that control kinesin attachment to/detachment from cargo vesicles,
which is essential to the overall transport efficiency of cargo vesicles. To
quantitatively understand the role of the regulatory protein, we review the in
vivo physical parameters of UNC-104-mediated SVP transport, including force,
velocity, run length and run time, derived from wild-type and arl-8-deletion
mutant C. elegans. Our future aim is to facilitate the construction of a
consensus physical model to connect SVP transport with pathologies related to
deficient synapse construction caused by the deficient UNC-104 regulation. We
hope that the physical parameters of SVP transport summarized in this review
become a useful guide for the development of such model.


  The World Health Organisation currently recommends pre-screening for past
infection prior to administration of the only licensed dengue vaccine, CYD-TDV.
Using a bounding analysis, we show that despite additional testing costs, this
approach can improve the economic viability of CYD-TDV: effective testing
reduces unnecessary vaccination costs while increasing the health benefit for
vaccine recipients. When testing is cheap enough, those trends outweigh
additional screening costs and make test-then-vaccinate strategies
net-beneficial in many settings.
  We derived these results using a general approach for determining price
thresholds for testing and vaccination, as well as indicating optimal start and
end ages of routine test-then-vaccinate programs. This approach only requires
age-specific seroprevalence and a cost estimate for second infections. We
demonstrate this approach across settings commonly used to evaluate CYD-TDV
economics, and highlight implications of our simple model for more detailed
studies. We found trends showing test-then-vaccinate strategies are generally
more beneficial starting at younger ages, and that in some settings multiple
years of testing can be more beneficial than only testing once, despite
increased investment in testing.


  Several methods now exist to solvate lipid bilayer discoids at the scale of
tens of nanometres. Due to their size, such nanodiscoids have a comparatively
large boundary-to-area ratio, making them unusually well-suited to probing the
effects of Gaussian curvature. Arguing that fluctuations in discoid size and
shape are quenched on formation, we quantify the stability, in terms of size
and shape, of near-solvation discoid-like flaps that are subject to thermal
fluctuations. Using cryo-Electron Microscopy images of Styrene Maleic Acid
stabilised discoids, we deduce that stable, saddle-like discoids (with high
Gaussian curvature) can likely be solvated from bulk lamellar ($L_\alpha$)
phase at moderate-to-high surface tensions ($>10^{-4}$ N/m). We then describe
how such tension-controlled solvation can be used for both measuring, and
fractionating membrane components according-to, the modulus of Gaussian
rigidity $\bar{\kappa}$. Opportunities for investigating the effects of
Gaussian curvature on membrane-embedded proteins, which can be co-solvated
during the formation process, are also discussed.


  The Price equation describes the change in populations. Change concerns some
value, such as biological fitness, information or physical work. The Price
equation reveals universal aspects for the nature of change, independently of
the meaning ascribed to values. By understanding those universal aspects, we
can see more clearly why fundamental mathematical results in different
disciplines often share a common form. We can also interpret more clearly the
meaning of key results within each discipline. For example, the mathematics of
natural selection in biology has a form closely related to information theory
and physical entropy. Does that mean that natural selection is about
information or entropy? Or do natural selection, information and entropy arise
as interpretations of a common underlying abstraction? The Price equation
suggests the latter. The Price equation achieves its abstract generality by
partitioning change into two terms. The first term naturally associates with
the direct forces that cause change. The second term naturally associates with
the changing frame of reference. In the Price equation's canonical form, total
change remains zero because the conservation of total probability requires that
all probabilities invariantly sum to one. Much of the shared common form for
the mathematics of different disciplines may arise from that seemingly trivial
invariance of total probability, which leads to the partitioning of total
change into equal and opposite components of the direct forces and the changing
frame of reference.


  Identifying stabilizing factors in foodwebs is a long standing challenge with
wide implications for community ecology and conservation. Here, we investigate
the stability of spatially resolved meta-foodwebs with far-ranging
super-predators for whom the whole meta-foodwebs appears to be a single
habitat. By using a combination of generalised modeling with a master stability
function approach, we are able to efficiently explore the asymptotic stability
of large classes of realistic many-patch meta-foodwebs. We show that
meta-foodwebs with far-ranging top predators are more stable than those with
localized top predators. Moreover, adding far-ranging generalist top predators
to a system can have a net stabilizing effect, despite increasing the food web
size. These results highlight the importance of top predator conservation.


  The theory of repeated games analyzes the long-term relationship of
interacting players and mathematically reveals the condition of how cooperation
is achieved, which is not achieved in a one-shot game. In the repeated
prisoner's dilemma (RPD) game with no errors, zero-determinant (ZD) strategies
allow a player to unilaterally set a linear relationship between the player's
own payoff and the opponent's payoff regardless of the strategy that the
opponent implements. In contrast, unconditional strategies such as ALLD and
ALLC also unilaterally set a linear payoff relationship. Errors often happen
between players in the real world. However, little is known about the existence
of such strategies in the RPD game with errors. Here, we analytically search
for all strategies that enforce a linear payoff relationship under observation
errors in the RPD game. As a result, we found that, even in the case with
observation errors, the only strategy sets that enforce a linear payoff
relationship are either ZD strategies or unconditional strategies and that no
other strategies can enforce it, which were numerically confirmed.


  RNA-RNA interaction (RRI) is ubiquitous and has complex roles in the cellular
functions. In human health studies, miRNA-target and lncRNAs are among an elite
class of RRIs that have been extensively studied. Bacterial ncRNA-target and
RNA interference are other classes of RRIs that have received significant
attention. In recent studies, mRNA-mRNA interaction instances have been
observed, where both partners appear in the same pathway without any direct
link between them, or any prior knowledge about their relationship. Those
recently discovered cases suggest that RRI scope is much wider than those
aforementioned elite classes.
  We revisit our RNA-RNA interaction partition function algorithm, piRNA, which
computes the partition function, base-pairing probabilities, and structure for
the comprehensive Turner energy model using 96 different dynamic programming
tables. In this study, we strategically retreat from sophisticated
thermodynamic models to the much simpler base pair counting model. That might
seem counter-intuitive at the first glance; our idea is to benefit from the
advantages of such simple models in terms of running time and memory footprint
and compensate for the associated information loss by adding machine learning
components in the future.
  Here, simple weighted base pair counting is considered to obtain BPPart for
Base-pair Partition function and BPMax for Base-pair Maximization, which use 9
and 2 tables respectively. They are empirically 225 and 1350 fold faster than
piRNA. A correlation of 0.855 and 0.836 was achieved between piRNA and BPPart
and between piRNA and BPMax, respectively, in 37 degrees, and 0.920 and 0.904
in -180 degrees. We also discover two partner RNAs, SNORD3D and TRAF3, and
hypothesize their potential roles in genetic diseases. We envision fusion of
machine learning methods with the proposed algorithms in the future.


  Translating the vast data generated by genomic platforms into reliable
predictions of clinical outcomes remains a critical challenge in realizing the
promise of genomic medicine largely due to small number of independent samples.
In this paper, we show that neural networks can be trained to predict clinical
outcomes using heterogeneous genomic data sources via multi-task learning and
adversarial representation learning, allowing one to combine multiple cohorts
and outcomes in training. We compare our proposed method to two baselines and
demonstrate that it can be used to help mitigate the data scarcity and clinical
outcome censorship in cancer genomics learning problems.


  Although the key promoter elements necessary to drive transcription in
Escherichia coli have long been understood, we still cannot predict the
behavior of arbitrary novel promoters, hampering our ability to characterize
the myriad of sequenced regulatory architectures as well as to design novel
synthetic circuits. This work builds upon a beautiful recent experiment by
Urtecho et al. who measured the gene expression of over 10,000 promoters
spanning all possible combinations of a small set of regulatory elements. Using
this data, we demonstrate that a central claim in energy matrix models of gene
expression - that each promoter element contributes independently and
additively to gene expression - contradicts experimental measurements. We
propose that a key missing ingredient from such models is the avidity between
the -35 and -10 RNA polymerase binding sites and develop what we call a refined
energy matrix model that incorporates this effect. We show that this the
refined energy matrix model can characterize the full suite of gene expression
data and explore several applications of this framework, namely, how
multivalent binding at the -35 and -10 sites can buffer RNAP kinetics against
mutations and how promoters that bind overly tightly to RNA polymerase can
inhibit gene expression. The success of our approach suggests that avidity
represents a key physical principle governing the interaction of RNA polymerase
to its promoter.


  Ligands entering a protein binding pocket essentially compete with water
molecules for binding to the protein. Hence, the location and thermodynamic
properties of water molecules in protein structures have gained increased
attention in the drug design community. Including corresponding data into 3D
pharmacophore modeling is essential for efficient high throughput virtual
screening. Here, we present PyRod, a free and open-source python software that
allows for visualization of pharmacophoric binding pocket characteristics,
identification of hot spots for ligand binding and subsequent generation of
pharmacophore features for virtual screening. The implemented routines analyze
the protein environment of water molecules in molecular dynamics (MD)
simulations and can differentiate between hydrogen bonded waters as well as
waters in a protein environment of hydrophobic, charged or aromatic atom
groups. The gathered information is further processed to generate dynamic
molecular interaction fields (dMIFs) for visualization and pharmacophoric
features for virtual screening. The described software was applied to 5
therapeutically relevant drug targets and generated pharmacophores were
evaluated using DUD-E benchmarking sets. The best performing pharmacophore was
found for the HIV1 protease with an early enrichment factor of 54.6. PyRod adds
a new perspective to structure-based screening campaigns by providing
easy-to-interpret dMIFs and purely protein-based 3D pharmacophores that are
solely based on tracing water molecules in MD simulations. Since structural
information about co-crystallized ligands is not needed, screening campaigns
can be followed, for which less or no ligand information is available. PyRod is
freely available at https://github.com/schallerdavid/pyrod.


  Living cells encode and transmit information in the temporal dynamics of
biochemical components. Gaining a detailed understanding of the input-output
relationship in biological systems therefore requires quantitative measures
that capture the interdependence between complete time trajectories of
biochemical components. Mutual information provides such a measure but its
calculation in the context of stochastic reaction networks is associated with
mathematical challenges. Here we show how to estimate the mutual information
between complete paths of two molecular species that interact with each other
through biochemical reactions. We demonstrate our approach using three simple
case studies.


  By constructing a Landau-like energy functional, we investigate the molecular
organization of a three-component mixture in cell membranes. In the strongly
interacting limit, we model the interaction between molecules using pseudospin
variables and convert them into non-interacting variables using a mean-field
theory. Next, we construct the two-order parameter Landau-type energy
functional through the Helmholtz free energy. By analyzing the Landau free
energy, we map out the phase diagram focusing on homogeneous and various phase
separated states on the cell membrane.


  We demonstrate that the conclusions drawn by Lefcheck et al. (2019) regarding
the positive effects of fish diversity on coral reef ecosystem functioning
across scales are flawed because of a series of conceptual and statistical
issues that include spurious correlations, the conflation of population size
and species diversity effects and a failure to recognize that observing a
biodiversity effect at multiple sites is not equivalent to observing it at
multiple scales.


  Most mechanistic predator-prey modelling has involved either parameterization
from process rate data or inverse modelling. Here, we take a median road: we
aim at identifying the potential benefits of combining datasets, when both
population growth and predation processes are viewed as stochastic. We fit a
discrete-time, stochastic predator-prey model of the Leslie type to simulated
time series of densities and kill rate data. Our model has both environmental
stochasticity in the growth rates and interaction stochasticity, i.e., a
stochastic functional response. We examine what the kill rate data brings to
the quality of the estimates, and whether estimation is possible (for various
time series lengths) solely with time series of population counts or biomass
data. Both Bayesian and frequentist estimation are performed, providing
multiple ways to check model identifiability. The Fisher Information Matrix
suggests that models with and without kill rate data are all identifiable,
although correlations remain between parameters that belong to the same
functional form. However, our results show that if the attractor is a fixed
point in the absence of stochasticity, identifying parameters in practice
requires kill rate data as a complement to the time series of population
densities, due to the relatively flat likelihood. Only noisy limit cycle
attractors can be identified directly from population count data (as in inverse
modelling), although even in this case, adding kill rate data - including in
small amounts - can make the estimates much more precise. Overall, we show that
under process stochasticity in interaction rates, interaction data might be
essential to obtain identifiable dynamical models for multiple species. These
results may extend to other biotic interactions than predation, for which
similar models combining interaction rates and population counts could be
developed.


  The double-helical structure of DNA results from canonical base pairing and
stacking interactions. However, variations from steady-state conformations
result from mechanical perturbations in cells. These different topologies have
physiological relevance but their dependence on sequence remains unclear. Here,
we use molecular dynamics simulations to show that sequence differences result
in markedly different structural motifs upon physiological twisting and
stretching. We simulated overextension on four different sequences of DNA
((AA)12, (AT)12, (GG)12 and (GC)12) with supercoiling densities within the
physiological range. We found that DNA denatures in the majority of stretching
simulations, surprisingly including those with overtwisted DNA. GC-rich
sequences were observed to be more stable than AT-rich, with the specific
response dependent on base pair ordering. Furthermore, we found that (AT)12
forms stable periodic structures with non-canonical hydrogen bonds in some
regions and non-canonical stacking in others, whereas (GC)12 forms a stacking
motif of four base pairs independent of supercoiling density. Our results
demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around
and significantly above cellular supercoiling, and that the DNA sequence is
crucial for understanding structural changes under mechanical stress. Our
findings have important implications for the activities of protein machinery
interacting with DNA in all cells.


  Species tree estimation is a complex problem, due to the fact that different
parts of the genome can have different evolutionary histories than the genome
itself. One of the causes for this discord is incomplete lineage sorting (also
called deep coalescence), which is a population-level process that produces
gene trees that differ from the species tree. The last decade has seen a large
number of new methods developed to estimate species trees from multi-locus
datasets, specifically addressing this cause of gene tree heterogeneity. In
this paper, we review these methods, focusing mainly on issues that relate to
analyses of datasets containing large numbers of species or loci (or both). We
also discuss divide-and-conquer strategies for enabling species tree estimation
methods to run on large datasets, including new approaches that are based on
algorithms (such as TreeMerge) for the Disjoint Tree Merger problem.


  In 1972, Robert May triggered a worldwide research program studying
ecological communities using random matrix theory. Yet, it remains unclear if
and when we can treat real communities as random ecosystems. Here, we draw on
recent progress in random matrix theory and statistical physics to extend May's
approach to generalized consumer-resource models. We show that in diverse
ecosystems adding even modest amounts of noise to consumer preferences results
in a transition to "typicality" where macroscopic ecological properties of
communities are indistinguishable from those of random ecosystems, even when
resource preferences have prominent designed structures. We test these ideas
using numerical simulations on a wide variety of ecological models. Our work
offers an explanation for the success of random consumer-resource models in
reproducing experimentally observed ecological patterns in microbial
communities and highlights the difficulty of scaling up bottom-up approaches in
synthetic ecology to diverse communities.


  The multiple futile cycle is a phosphorylation system in which a molecular
substrate might be phosphorylated sequentially n times by means of an enzymatic
mechanism. The system has been studied mathematically using reaction network
theory and ordinary differential equations. It is known that the system might
have at least as many as 2[n/2]+1 steady states (where [x] is the integer part
of x) for particular choices of parameters. Furthermore, for the simple and
dual futile cycles (n=1,2) the stability of the steady states has been
determined in the sense that the only steady state of the simple futile cycle
is globally stable, while there exist parameter values for which the dual
futile cycle admits two asymptotically stable and one unstable steady state.
For general n, evidence that the possible number of asymptotically stable
steady states increases with $n$ has been given, which has led to the
conjecture that parameter values can be chosen such that [n/2]+1 out of
2[n/2]+1 steady states are asymptotically stable and the remaining steady
states are unstable.
  We prove this conjecture here by first reducing the system to a smaller one,
for which we find a choice of parameter values that give rise to a unique
steady state with multiplicity 2[n/2]+1. Using arguments from geometric
singular perturbation theory, and a detailed analysis of the centre manifold of
this steady state, we achieve the desired result.


  Living in the age of the digital media explosion, the amount of data that is
being stored increases dramatically. However, even if existing storage systems
suggest efficiency in capacity, they are lacking in durability. Hard disks,
flash, tape or even optical storage have limited lifespan in the range of 5 to
20 years. Interestingly, recent studies have proven that it was possible to use
synthetic DNA for the storage of digital data, introducing a strong candidate
to achieve data longevity. The DNA's biological properties allows the storage
of a great amount of information into an extraordinary small volume while also
promising efficient storage for centuries or even longer with no loss of
information. However, encoding digital data onto DNA is not obvious, because
when decoding, we have to face the problem of sequencing noise robustness.
Furthermore, synthesizing DNA is an expensive process and thus, controlling the
compression ratio by optimizing the rate-distortion trade-off is an important
challenge we have to deal with. This work proposes a coding solution for the
storage of digital images onto synthetic DNA. We developed a new encoding
algorithm which generates a DNA code robust to biological errors coming from
the synthesis and the sequencing processes. Furthermore, thanks to an optimized
allocation process the solution is able to control the compression ratio and
thus the length of the synthesized DNA strand. Results show an improvement in
terms of coding potential compared to previous state-of-the-art works.


  We have developed a global analysis model for randomly oriented, fully
hydrated inverted hexagonal (H$_\text{II}$) phases formed by many amphiphiles
in aqueous solution, including membrane lipids. The model is based on a
structure factor for hexagonally packed rods and a compositional model for the
scattering length density (SLD) enabling also the analysis of positionally
weakly correlated H$_\text{II}$ phases. For optimization of the adjustable
parameters we used Bayesian probability theory, which allows to retrieve
parameter correlations in much more detail than standard analysis techniques,
and thereby enables a realistic error analysis. The model was applied to
different phosphatidylethanolamines including previously not reported
H$_\text{II}$ data for diC14:0 and diC16:1 phosphatidylethanolamine. The
extracted structural features include intrinsic lipid curvature, hydrocarbon
chain length and area per lipid at the position of the neutral plane.


  Can insects transmit information about a threat through a physical barrier
they cannot cross? While the emission and detection of alarm pheromones is the
universally accepted mechanism of panic communication, there are still open
questions about the specific pheromones and type of ant-to-ant contacts
involved in specific scenarios, and how relevant is the role of non-chemical
communication. Here we report experiments suggesting that panicked ants can
transmit danger information through a physical barrier to nestmates located in
a safe arena when danger is simulated by using insect repellent. As ants in the
safe arena get danger information from nestmates trapped in an area
contaminated with a chemical, they retreat from the barrier in average. Visual
clues and pure mutual antennation seem to play a minor role as mechanisms of
danger information transmission through the obstacle. A detailed quantification
of the ant behavior reveals that examination of the legs of ants in the danger
arena by nestmates in the safe one may be relevant for danger information
transmission.


  A feature of human creativity is the ability to take a subset of existing
items (e.g. objects, ideas, or techniques) and combine them in various ways to
give rise to new items, which, in turn, fuel further growth. Occasionally, some
of these items may also disappear (extinction). We model this process by a
simple stochastic birth--death model, with non-linear combinatorial terms in
the growth coefficients to capture the propensity of subsets of items to give
rise to new items. In its simplest form, this model involves just two
parameters $(P, \alpha)$. This process exhibits a characteristic 'hockey-stick'
behaviour: a long period of relatively little growth followed by a relatively
sudden 'explosive' increase. We provide exact expressions for the mean and
variance of this time to explosion and compare the results with simulations. We
then generalise our results to allow for more general parameter assignments,
and consider possible applications to data involving human productivity and
creativity.


  A multicompartment mathematical model is presented with the goal of studying
the role of dose-dense protocols in the context of combination cancer
chemotherapy. Dose-dense protocols aim at reducing the period between courses
of chemotherapy from three to two weeks or less, in order to avoid the regrowth
of the tumor during the meantime and achieve maximum cell kill at the end of
the treatment. Inspired by clinical trials, we carry out a randomized
computational study to systematically compare a variety of protocols using two
drugs of different specificity. Our results suggest that cycle specific drugs
can be administered at low doses between courses of treatment to arrest the
relapse of the tumor. This might be a better strategy than reducing the period
between cycles.


  Most known proteins in nature consist of multiple domains. Interactions
between domains may lead to unexpected folding and misfolding phenomena. This
study of human {\gamma}D-crystallin, a two-domain protein in the eye lens,
revealed one such surprise: conformational catalysis of misfolding via
intermolecular domain interface ''stealing''. An intermolecular interface
between the more stable domains outcompetes the native intramolecular domain
interface. Loss of the native interface in turn promotes misfolding and
subsequent aggregation, especially in cataract-related {\gamma}D-crystallin
variants. This phenomenon is likely a contributing factor in the development of
cataract disease, the leading worldwide cause of blindness. However, interface
stealing likely occurs in many proteins composed of two or more interacting
domains.


  We demonstrate that the conclusions drawn by Bernhard et al. (2018) regarding
the ability of nonlinear averaging to accurately predict organismal performance
under fluctuating temperatures are flawed because of a series of experimental
and statistical issues that include the presence of a hidden treatment effect,
the use of a single low frequency temperature fluctuation that could easily be
tracked by the fast growing organism, and the decision to quantify performance
via population growth rate, a metric that can mask significant variation in
population size.


  Summary: The advent of Web-based tools that assist in the analysis and
visualization of macromolecules require application programming interfaces
(APIs) designed for modern web frameworks. To this end, we have developed a
Node.js module pdbmine that allows any user to generate faster data-request
queries to the RCSB Protein Data Bank (PDB). This JavaScript API acts as a
layer over the XML-based RCSB PDB RESTful API. The relatively simple nature of
the function calls within this module allows the user to easily implement and
integrate pdbmine into larger Node.js web applications.
  Availability: This module can be installed via the Node Package Manager (NPM)
at https://www.npmjs.com/package/pdbmine/, and is hosted on GitHub under the
open-source MIT license at https://github.com/nnj1/pdbmine/. Relevant
documentation is detailed at https://nnj1.github.io/pdbmine/


  ASTRAL is a method for reconstructing species trees after inferring a set of
gene trees and is increasingly used in phylogenomic analyses. It is
statistically consistent under the multi-species coalescent model, is scalable,
and has shown high accuracy in simulated and empirical studies. This chapter
discusses practical considerations in using ASTRAL, starting with a review of
published results and pointing to the strengths and weaknesses of species tree
estimation using ASTRAL. It then continues to detail the best ways to prepare
input gene trees, interpret ASTRAL outputs, and perform follow-up analyses.


  Motivation: Network visualizations of complex biological datasets usually
result in 'hairball' images, which do not discriminate network modules.
Results: We present the EntOptLayout Cytoscape plug-in based on a recently
developed network representation theory. The plug-in provides an efficient
visualization of network modules, which represent major protein complexes in
protein-protein interaction and signalling networks. Importantly, the tool
gives a quality score of the network visualization by calculating the
information loss between the input data and the visual representation showing a
3- to 25-fold improvement over conventional methods. Availability and
implementation: The plug-in (running on Windows, Linux, or Mac OS) and its
tutorial (both in written and video forms) can be downloaded freely under the
terms of the MIT license from: http://apps.cytoscape.org/apps/entoptlayout.
Supplementary data are available at Bioinformatics online. Contact:
csermely.peter@med.semmelweis-univ.hu