  Species sharing a prey or a predator species may go extinct due to
exploitative or apparent competition. We examine whether evolution of the
shared species acts as a coexistence mechanism and to what extent the answer
depends on the genetic architecture underlying trait evolution. In our models
of exploitative and apparent competition, the shared species evolves its
defense or prey use. Evolving species are either haploid or diploid. A single
locus pleiotropically determines prey nutritional quality and predator attack
rates. When pleiotropy is sufficiently antagonistic (e.g. nutritional prey are
harder to capture), eco-evolutionary assembly culminates in one of two stable
states supporting only two species. When pleiotropy is weakly antagonistic or
synergistic, assembly is intransitive: species-genotype pairs are cyclically
displaced by rare invasions of the missing genotypes or species. This
intransitivity allows for coexistence if, along its equilibria, the geometric
mean of recovery rates exceeds the geometric mean of loss rates of the rare
genotypes or species. By affecting these rates, synergistic pleiotropy can
mediate coexistence, while antagonistic pleiotropy does not. For diploid
populations experiencing weak antagonistic pleiotropy, superadditive allelic
contributions to fitness can mitigate coexistence via an eco-evolutionary
storage effect. Density-dependence and mutations also promote coexistence.
These results highlight how the efficacy of evolution as a coexistence
mechanism may depend on the underlying genetic architecture.


  In many important cellular processes, including mRNA translation, gene
transcription, phosphotransfer, and intracellular transport, biological
"particles" move along some kind of "tracks". The motion of these particles can
be modeled as a one-dimensional movement along an ordered sequence of sites.
The biological particles (e.g., ribosomes, RNAPs, phosphate groups, motor
proteins) have volume and cannot surpass one another. In some cases, there is a
preferred direction of movement along the track, but in general the movement
may be two-directional, and furthermore the particles may attach or detach from
various regions along the tracks (e.g. ribosomes may drop off the mRNA molecule
before reaching a stop codon). We derive a new deterministic mathematical model
for such transport phenomena that may be interpreted as the dynamic mean-field
approximation of an important model from mechanical statistics called the
asymmetric simple exclusion process (ASEP) with Langmuir kinetics. Using tools
from the theory of monotone dynamical systems and contraction theory we show
that the model admits a unique globally asymptotically stable equilibrium. This
means that the occupancy in all the sites along the lattice converges to a
steady-state value that depends on the parameters but not on the initial
conditions. We also show that the model entrains (or phase locks) to periodic
excitations in any of its forward, backward, attachment, or detachment rates.
We demonstrate an application of this phenomenological transport model for
analyzing the effect of ribosome drop off in mRNA translation. One may perhaps
expect that drop off from a jammed site may increase the total flow by reducing
congestion. Our results show that this is not true. Drop off has a substantial
effect on the flow, yet always leads to a reduction in the steady-state protein
production rate.


  The quasispecies model was introduced in 1971 by Manfred Eigen to discuss the
first stages of life on Earth. It provides an appealing mathematical framework
to study the evolution of populations in biology, for instance viruses. We
present briefly the model and we focus on its stationary solutions. These
formulae have a surprisingly rich combinatorial structure, involving for
instance the Eulerian and Stirling numbers, as well as the up--down
coefficients of permutations.


  Introduction: Zika epidemic in America was declared a public health emergency
of international concern after the rapid spread in the region. Stratification
of the potential transmission of the disease is needed to address the efforts
surveillance and disease control. The goal of this research is to compare the
basic reproductive number of Zika in different municipalities, from an SIR
model with implicit vector dynamics, based on the daily case reporting data of
Antioquia, Colombia, the second most affected country after Brazil. Methods: An
simple SIR model with implicit vector dynamics was derived and used. The
approximate solution of the model in terms of individuals recovered in each
time unit, allowed to obtain estimate of the model parameters including the
basic reproduction number (Ro) and its 95% confidence intervals. These
parameters were estimated via fitting the solution of the model to the daily
reported cumulative cases data from the regional surveillance system. Results:
Ro was estimated for 20 municipalities, all located at less than 2200 meters
above sea level. From January to April 2016, between 17 and 347 cases were
reported from these municipalities. Of these, 15 municipalities had a high
potential for transmission (Ro>1) and 5 had less potential (Ro<1), although in
3 of these, transmission was later possible because its upper confidence
interval of Ro was greater than one. Conclusion: Surveillance and control of
Zika should be directed primarily to municipalities with Ro>1. Furthermore,
strategies that will strengthen the detection and management of cases in the
remaining municipalities should be beneficial to the whole Antioquia region.


  Protein binding often involves conformational changes. Important questions
are whether a conformational change occurs prior to a binding event
('conformational selection') or after a binding event ('induced fit'), and how
conformational transition rates can be obtained from experiments. In this
article, we present general results for the chemical relaxation rates of
conformational-selection and induced-fit binding processes that hold for all
concentrations of proteins and ligands and, thus, go beyond the standard
pseudo-first-order approximation of large ligand concentration. These results
allow to distinguish conformational-selection from induced-fit processes - also
in cases in which such a distinction is not possible under pseudo-first-order
conditions - and to extract conformational transition rates of proteins from
chemical relaxation data.


  The aim of this study is to investigate the relation that can be found
between the phylogeny of a large set of complete chloroplast genomes, and the
evolution of gene content inside these sequences. Core and pan genomes have
been computed on \textit{de novo} annotation of these 845 genomes, the former
being used for producing well-supported phylogenetic tree while the latter
provides information regarding the evolution of gene contents over time. It
details too the specificity of some branches of the tree, when specificity is
obtained on accessory genes. After having detailed the material and methods, we
emphasize some remarkable relation between well-known events of the chloroplast
history, like endosymbiosis, and the evolution of gene contents over the
phylogenetic tree.


  The existence of phenotypic heterogeneity in single-species bacterial
biofilms is well-established in the published literature. However, the modeling
of population dynamics in biofilms from the viewpoint of social interactions,
i.e. interplay between heterotypic strains, and the analysis of this kind using
control theory are not addressed significantly. Therefore, in this paper, we
theoretically analyze the population dynamics model in microbial biofilms with
non-participating strains (coexisting with public goods producers and
non-producers) in the context of evolutionary game theory and nonlinear
dynamics. Our analysis of the replicator dynamics model is twofold: first
without the inclusion of spatial pattern, and second with the consideration of
degree of assortment. In the first case, Lyapunov stability analysis of the
stable equilibrium point of the proposed replicator system determines (1,0)
('full dominance of cooperators') as a global asymptotic stable equilibrium
whenever the return exceeds the metabolic cost of cooperation. Hence, the
global asymptotic stable nature of (1,0) in the context of non-consideration of
spatial pattern helps to justify mathematically the adversity in the
eradication of "cooperative enterprise" that is an infectious biofilm. In the
second case, we found non-existence of global asymptotic stability in the
system, and it unveils two additional phenomena - bistability and coexistence.
In this context, two inequality conditions are derived for the 'full dominance
of cooperators' and coexistence. Therefore, the inclusion of spatial pattern in
biofilms with non-competing strains intends conditional dominance of pathogenic
(with respect to the hosts) public goods producers which can be an effective
strategy towards the control of an infectious biofilm with the drug-dependent
regulation of degree of segregation.


  Gene expression is a readily-observed quantification of transcriptional
activity and cellular state that enables the recovery of the relationships
between regulators and their target genes. Reconstructing transcriptional
regulatory networks from gene expression data is a problem that has attracted
much attention, but previous work often makes the simplifying (but unrealistic)
assumption that regulator activity is represented by mRNA levels. We use a
latent tree graphical model to analyze gene expression without relying on
transcription factor expression as a proxy for regulator activity. The latent
tree model is a type of Markov random field that includes both observed gene
variables and latent (hidden) variables, which factorize on a Markov tree.
Through efficient unsupervised learning approaches, we determine which groups
of genes are co-regulated by hidden regulators and the activity levels of those
regulators. Post-processing annotates many of these discovered latent variables
as specific transcription factors or groups of transcription factors. Other
latent variables do not necessarily represent physical regulators but instead
reveal hidden structure in the gene expression such as shared biological
function. We apply the latent tree graphical model to a yeast stress response
dataset. In addition to novel predictions, such as condition-specific binding
of the transcription factor Msn4, our model recovers many known aspects of the
yeast regulatory network. These include groups of co-regulated genes,
condition-specific regulator activity, and combinatorial regulation among
transcription factors. The latent tree graphical model is a general approach
for analyzing gene expression data that requires no prior knowledge of which
possible regulators exist, regulator activity, or where transcription factors
physically bind.


  Life expectancies at birth are routinely computed from period life tables.
Such period life expectancies may be distorted by selection when comparing
countries where the living conditions improved earlier (like Norway and Sweden)
with countries where they improved later (like Italy and Japan). One way to get
a fair comparison between the countries, is to use cohort data and consider the
expected number of years lost before a given age a. Contrary to the results
based on period data, one then finds that Italian women may expect to lose more
years than women in Norway and Sweden, while there are no indications that
Japanese women will lose fewer years than Scandinavian women.


  Molecular profiling data (e.g., gene expression) has been used for clinical
risk prediction and biomarker discovery. However, it is necessary to integrate
other prior knowledge like biological pathways or gene interaction networks to
improve the predictive ability and biological interpretability of biomarkers.
Here, we first introduce a general regularized Logistic Regression (LR)
framework with regularized term $\lambda \|\bm{w}\|_1 +
\eta\bm{w}^T\bm{M}\bm{w}$, which can reduce to different penalties, including
Lasso, elastic net, and network-regularized terms with different $\bm{M}$. This
framework can be easily solved in a unified manner by a cyclic coordinate
descent algorithm which can avoid inverse matrix operation and accelerate the
computing speed. However, if those estimated $\bm{w}_i$ and $\bm{w}_j$ have
opposite signs, then the traditional network-regularized penalty may not
perform well. To address it, we introduce a novel network-regularized sparse LR
model with a new penalty $\lambda \|\bm{w}\|_1 + \eta|\bm{w}|^T\bm{M}|\bm{w}|$
to consider the difference between the absolute values of the coefficients. And
we develop two efficient algorithms to solve it. Finally, we test our methods
and compare them with the related ones using simulated and real data to show
their efficiency.


  Physical interactions between proteins are often difficult to decipher. The
aim of this paper is to present an algorithm designed to recognize binding
patches and supporting structural scaffolds of interacting heterodimer protein
chains using the Gaussian Network Model (GNM). The recognition is based on the
(self)adjustable identification of kinetically hot residues, i.e., residues
with highest contributions to the weighted sum of fastest modes per chain
extracted via GNM, and their connection to possible binding scaffolds. The
algorithm adjusts the number of modes used in the GNM's weighted sum
calculation using the ratio of predicted and expected numbers of target
residues (contact and first layer residues). This approach produces very good
results when applied to chains forming heterodimers, especially with dimers
with high chain length ratios. The protocol's ability to recognize near native
decoys was compared to the ability of the statistical potential of Lu and
Skolnick using the Sternberg and Vakser decoy dimers sets. The statistical
potential produced better overall results, but in a number of cases its
predicting ability was comparable, or even worse than the ability of the
adjustable GNM approach. The results presented in this paper suggest that in
heterodimers, at least one partnering chain has interacting scaffold determined
by the immovable kinetically hot residues. In many cases interacting chains
(especially if being of noticeably different sizes), either behave as rigid
lock and key, or exhibit opposite dynamic behaviors. While the binding surface
of one of the chains is rigid and stable, its partner's interacting scaffold is
more flexible and adaptable.
  Authors note: The approach described here was initially given as a rough
draft in 2013 [1]. The next manuscript will describe the behavior of protein
dimers incorrectly characterized with the present approach.


  The uniqueness issue of SDE decomposition theory proposed by Ao and his
co-workers has recently been discussed. A comprehensive study to investigate
connections among different landscape theories [J. Chem. Phys. 144, 094109
(2016)] has pointed out that the decomposition is generally not unique, while
Ao et al. (arXiv:1603.07927v1) argues that such conclusions are "incorrect"
because of the missing boundary conditions. In this comment, we will combine
literatures research and concrete examples to show that the concrete and
effective boundary conditions have not been proposed to guarantee the
uniqueness, hence the arguments in [arXiv:1603.07927v1] are not sufficient.
Moreover, we show that the "uniqueness" of the O-U process decomposition
referred by YTA paper is unable to serve as a counterexample to ZL's result
since additional assumptions have been made implicitly beyond the original SDE
decomposition framework, which cannot be applied to general nonlinear cases.
Some other issues such as the failure of gradient expansion method will also be
discussed. Our demonstration contributes to better understanding of the
relevant papers as well as the SDE decomposition theory.


  In cryo-electron microscopy (EM), molecular structures are determined from
large numbers of projection images of individual particles. To harness the full
power of this single-molecule information, we use the Bayesian inference of EM
(BioEM) formalism. By ranking structural models using posterior probabilities
calculated for individual images, BioEM in principle addresses the challenge of
working with highly dynamic or heterogeneous systems not easily handled in
traditional EM reconstruction. However, the calculation of these posteriors for
large numbers of particles and models is computationally demanding. Here we
present highly parallelized, GPU-accelerated computer software that performs
this task efficiently. Our flexible formulation employs CUDA, OpenMP, and MPI
parallelization combined with both CPU and GPU computing. The resulting BioEM
software scales nearly ideally both on pure CPU and on CPU+GPU architectures,
thus enabling Bayesian analysis of tens of thousands of images in a reasonable
time. The general mathematical framework and robust algorithms are not limited
to cryo-electron microscopy but can be generalized for electron tomography and
other imaging experiments.


  Understanding if and how mutants reach fixation in populations is an
important question in evolutionary biology. We study the impact of population
growth has on the success of mutants. To systematically understand the effects
of growth we decouple competition from reproduction; competition follows a
birth--death process and is governed by an evolutionary game, while growth is
determined by an externally controlled branching rate. In stochastic
simulations we find non-monotonic behaviour of the fixation probability of
mutants as the speed of growth is varied; the right amount of growth can lead
to a higher success rate. These results are observed in both coordination and
coexistence game scenarios, and we find that the 'one-third law' for
coordination games can break down in the presence of growth. We also propose a
simplified description in terms of stochastic differential equations to
approximate the individual-based model.


  Energy evaluation using fast Fourier transforms enables sampling billions of
putative complex structures and hence revolutionized rigid protein-protein
docking. However, in current methods efficient acceleration is achieved only in
either the translational or the rotational subspace. Developing an efficient
and accurate docking method that expands FFT based sampling to 5 rotational
coordinates is an extensively studied but still unsolved problem. The algorithm
presented here retains the accuracy of earlier methods but yields at least
tenfold speedup. The improvement is due to two innovations. First, the search
space is treated as the product manifold $\mathbf{SO(3)x(SO(3)\setminus S^1)}$,
where $\mathbf{SO(3)}$ is the rotation group representing the space of the
rotating ligand, and $\mathbf{(SO(3)\setminus S^1)}$ is the space spanned by
the two Euler angles that define the orientation of the vector from the center
of the fixed receptor toward the center of the ligand. This representation
enables the use of efficient FFT methods developed for $\mathbf{SO(3)}$.
Second, we select the centers of highly populated clusters of docked
structures, rather than the lowest energy conformations, as predictions of the
complex, and hence there is no need for very high accuracy in energy
evaluation. Therefore it is sufficient to use a limited number of spherical
basis functions in the Fourier space, which increases the efficiency of
sampling while retaining the accuracy of docking results. A major advantage of
the method is that, in contrast to classical approaches, increasing the number
of correlation function terms is computationally inexpensive, which enables
using complex energy functions for scoring.


  Phenotypic plasticity and its evolution may help evolutionary rescue in a
novel and stressful environment, especially if environmental novelty reveals
cryptic genetic variation that enables the evolution of increased plasticity.
However, the environmental stochasticity ubiquitous in natural systems may
alter these predictions because high plasticity may amplify
phenotype-environment mismatches. Although previous studies have highlighted
this potential detrimental effect of plasticity in stochastic environments,
they have not investigated how it affects extinction risk in the context of
evolutionary rescue and with evolving plasticity. We investigate this question
here by integrating stochastic demography with quantitative genetic theory in a
model with simultaneous change in the mean and predictability (temporal
autocorrelation) of the environment. We develop an approximate prediction of
long-term persistence under the new pattern of environmental fluctuations, and
compare it with numerical simulations for short- and long-term extinction risk.
We find that reduced predictability increases extinction risk and reduces
persistence because it increases stochastic load during rescue. This
understanding of how stochastic demography, phenotypic plasticity, and
evolution interact when evolution acts on cryptic genetic variation revealed in
a novel environment can inform expectations for invasions, extinctions, or the
emergence of chemical resistance in pests.


  We present an algorithm which explores permutation symmetries to describe the
time evolution of agent-based epidemic models. The main idea to improve
computation times relies on restricting the stochastic process to one sector of
the vector space, labeled by a single permutation eigenvalue. In this scheme,
the transition matrix reduces to block diagonal form, enhancing computational
performance.


  Cellulosomes are complex multi-enzyme machineries which efficiently degrade
plant cell-wall polysaccharides. The multiple domains of the cellulosome
proteins are often tethered together by intrinsically disordered regions. The
properties and functions of these disordered linkers are not well understood.
In this work, we study endoglucanase Cel8A, which is a relevant enzymatic
component of the cellulosomes of Clostridium thermocellum. We use both all-atom
and coarse-grained simulations to investigate how the equilibrium conformations
of the catalytic domain of Cel8A are affected by the disordered linker at its C
terminus. We find that when the endoglucanase is bound to its substrate, the
effective stiffness of the linker can influence the distances between groups of
amino-acid residues throughout the entire enzymatic domain. In particular,
variations in the linker stiffness can lead to small changes in the geometry of
the active-site cleft. We suggest that such geometrical changes may, in turn,
have an effect on the catalytic activity of the enzyme.


  In times of plenty expectations rise, just as in times of crisis they fall.
This can be mathematically described as a Win-Stay-Lose-Shift strategy with
dynamic aspiration levels, where individuals aspire to be as wealthy as their
average neighbor. Here we investigate this model in the realm of evolutionary
social dilemmas on the square lattice and scale-free networks. By using the
master equation and Monte Carlo simulations, we find that cooperators coexist
with defectors in the whole phase diagram, even at high temptations to defect.
We study the microscopic mechanism that is responsible for the striking
persistence of cooperative behavior and find that cooperation spreads through
second-order neighbors, rather than by means of network reciprocity that
dominates in imitation-based models. For the square lattice the master equation
can be solved analytically in the large temperature limit of the Fermi
function, while for other cases the resulting differential equations must be
solved numerically. Either way, we find good qualitative agreement with the
Monte Carlo simulation results. Our analysis also reveals that the evolutionary
outcomes are to a large degree independent of the network topology, including
the number of neighbors that are considered for payoff determination on
lattices, which further corroborates the local character of the microscopic
dynamics. Unlike large-scale spatial patterns that typically emerge due to
network reciprocity, here local checkerboard-like patterns remain virtually
unaffected by differences in the macroscopic properties of the interaction
network.


  I study two mechanisms based on punishment to promote cooperation in the
two-population snowdrift game. The first mechanism follows the traditional
approach in the literature and is based on the inclusion of a third additional
strategy in the payoff matrix of the stage-game. The second mechanism consists
of letting cooperators to punish defectors with a given exogenous frequency.
While both mechanisms share the same result regarding the minimum required
level of punishment in order to eliminate defectors in both populations,
stability in the mechanism following the second approach is more robust in the
sense that extinction of defectors is a globally asymptotically stable state
for any interior initial conditions in the phase space. Results were obtained
analytically through non-linear differential equations and also using an
agent-based simulation. There was a good level of agreement between both
approaches with respect to the evolutionary pattern over time and the possible
steady-states.


  We propose a new type of traveling wave pattern, one that can adapt to the
size of physical system in which it is embedded. Such a system arises when the
initial state has an instability that extends down to zero wavevector,
connecting at that point to two symmetry modes of the underlying dynamical
system. The Min system of proteins in E. coli is such as system with the
symmetry emerging from the global conservation of two proteins, MinD and MinE.
For this and related systems, traveling waves can adiabatically deform as the
system is increased in size without the increase in node number that would be
expected for an oscillatory version of a Turing instability containing an
allowed wavenumber band with a finite minimum.


  Motivation: New long read sequencers promise to transform sequencing and
genome assembly by producing reads tens of kilobases long. However their high
error rate significantly complicates assembly and requires expensive correction
steps to layout the reads using standard assembly engines.
  Results: We present an original and efficient spectral algorithm to layout
the uncorrected nanopore reads, and its seamless integration into a
straightforward overlap/layout/consensus (OLC) assembly scheme. The method is
shown to assemble Oxford Nanopore reads from several bacterial genomes into
good quality (~99% identity to the reference) genome-sized contigs, while
yielding more fragmented assemblies from a Sacharomyces cerevisiae reference
strain.
  Availability and implementation: http://github.com/antrec/spectrassembler
  Contact: antoine.recanati@inria.fr


  RNA can be used as a high-density medium for data storage and transmission;
however, an important RNA process -- replication -- is noisy. This paper
presents an error analysis for RNA as a data transmission medium, analyzing how
deletion errors increase in a collection of replicated DNA strands over time.


  During last years theoretical works shed new light and proposed new
hypothesis on the mechanisms which regulate the time behaviour of biological
populations in different natural systems. Despite of this, the role of
environmental variables in ecological systems is still an open question.
Filling this gap of knowledge is a crucial task for a deeper comprehension of
the dynamics of biological populations in real ecosystems. In this work we
study how the dynamics of food spoilage bacteria influences the sensory
characteristics of fresh fish specimens. This topic is crucial for a better
understanding of the role played by the bacterial growth on the organoleptic
properties, and for the quality evaluation and risk assessment of food
products. We therefore analyze the time behaviour, in fresh fish specimens, of
sensory characteristics starting from the growth curves of two spoilage
bacterial communities. The theoretical study, initially based on a
deterministic model, exploits experimental temperature profiles. As a first
step, a model of predictive microbiology is used to reproduce the experimental
behaviour of the two bacterial populations. Afterwards, the theoretical
bacterial growths are converted, through suitable differential equations, into
"sensory" scores, based on the Quality Index Method (QIM), a scoring system for
freshness and quality sensory estimation of fishery products. As a third step,
the theoretical curves of QIM scores are compared with the experimental data
obtained by sensory analysis. Finally, the differential equations for QIM
scores are modified by adding terms of multiplicative white noise, which mimics
the effects of uncertainty and variability in sensory analysis. A better
agreement between experimental and theoretical QIM scores is observed, in some
cases, in the presence of suitable values of noise intensity respect to the
deterministic analysis.


  Gene expression is inherently a noisy process which manifests as cell-to-cell
variability in time evolution of proteins. Consequently, events that trigger at
critical threshold levels of regulatory proteins exhibit stochasticity in their
timing. An important contributor to the noise in gene expression is translation
bursts which correspond to randomness in number of proteins produced in a
single mRNA lifetime. Modeling timing of an event as a first-passage time (FPT)
problem, we explore the effect of burst size distribution on event timing.
Towards this end, the probability density function of FPT is computed for a
gene expression model with burst size drawn from a generic non-negative
distribution. Analytical formulas for FPT moments are provided in terms of
known vectors and inverse of a matrix. The effect of burst size distribution is
investigated by looking at how the feedback regulation strategy that minimizes
noise in timing around a given time deviates from the case when burst is
deterministic. Interestingly, results show that the feedback strategy for
deterministic burst case is quite robust to change in burst size distribution,
and deviations from it are confined to about 20% of the optimal value. These
findings facilitate an improved understanding of noise regulation in event
timing.


  The structure of social networks is a key determinant in fostering
cooperation and other altruistic behavior among naturally selfish individuals.
However, most real social interactions are temporal, being both finite in
duration and spread out over time. This raises the question of whether stable
cooperation can form despite an intrinsically fragmented social fabric. Here we
develop a framework to study the evolution of cooperation on temporal networks
in the setting of the classic Prisoner's Dilemma. By analyzing both real and
synthetic datasets, we find that temporal networks generally facilitate the
evolution of cooperation compared to their static counterparts. More
interestingly, we find that the intrinsic human interactive pattern like bursty
behavior impedes the evolution of cooperation. Finally, we introduce a measure
to quantify the temporality present in networks and demonstrate that there is
an intermediate level of temporality that boosts cooperation most. Our results
open a new avenue for investigating the evolution of cooperation in more
realistic structured populations.


  Here we give a detailed description of MOLS 2.0.The source code and binary of
MOLS 2.0 are distributed free (under GNU Lesser General Public License) for the
scientific community. They are freely available for download at
https://sourceforge.net/projects/mols2-0/files/.


  Many chemical reactions in biological cells occur at very low concentrations
of constituent molecules. Thus, transcriptional gene-regulation is often
controlled by poorly expressed transcription-factors, such as E.coli lac
repressor with few tens of copies. Here we study the effects of inherent
concentration fluctuations of substrate-molecules on the seminal
Michaelis-Menten scheme of biochemical reactions. We present a universal
correction to the Michaelis-Menten equation for the reaction-rates. The
relevance and validity of this correction for enzymatic reactions and
intracellular gene-regulation is demonstrated. Our analytical theory and
simulation results confirm that the proposed variance-corrected
Michaelis-Menten equation predicts the rate of reactions with remarkable
accuracy even in the presence of large non-equilibrium concentration
fluctuations. The major advantage of our approach is that it involves only the
mean and variance of the substrate-molecule concentration. Our theory is
therefore accessible to experiments and not specific to the exact source of the
concentration fluctuations.


  We study the degree of success of a single predator hunting a herd of prey on
a two dimensional square lattice landscape. We explicitly consider the self
volume of the prey restraining their dynamics on the lattice. The movement of
both predator and prey is chosen to include an intelligent, decision making
step based on their respective sighting ranges, the radius in which they can
detect the other species (prey cannot recognise each other besides the self
volume interaction): after spotting each other the motion of prey and predator
turns from a nearest neighbour random walk into direct escape or chase,
respectively. We consider a large range of prey densities and sighting ranges
and compute the mean first passage time for a predator to catch a prey as well
as characterise the effective dynamics of the hunted prey. We find that the
prey's sighting range dominates their life expectancy and the predator profits
more from a bad eyesight of the prey than from his own good eye sight. We
characterise the dynamics in terms of the mean distance between the predator
and the nearest prey. It turns out that effectively the dynamics of this
distance coordinate can be captured in terms of a simple Ornstein-Uhlenbeck
picture. Reducing the many-body problem to a simple two-body problem by
imagining predator and nearest prey to be connected by a Hookean bond, all
features of the model such as prey density and sighting ranges merge into the
effective binding constant.


  The volume occupied by unconstrained bacterial DNA in physiological solutions
exceeds 1000 times the volume of the cell. Still, it is confined to a well
defined region of the cell called the nucleoid, which occupies only a fraction
of the cell volume. There is still no general agreement on the mechanism
leading to the compaction of the DNA and the formation of the nucleoid.
However, advances in in vivo sub-wavelength resolution microscopy techniques
have recently allowed the observation of the nucleoid at an unprecedented level
of detail. In particular, these observations show that the compaction of the
nucleoid is not static but is instead a highly dynamic feature, which depends
on several factors, like the richness of the nutrient, the cell cycle stage,
temperature, the action of an osmotic shock or antibiotics, etc. After a short
description of the electrolyte content of the cytosol and a brief overview of
the different mechanisms that may lead to the formation of the nucleoid, this
paper reviews some of the most fascinating recent results of in vivo
sub-wavelength resolution microscopy. It is furthermore argued that these
observations provide converging indications in favor of a model that describes
the cytosol as an aqueous electrolyte solution containing several
macromolecular species, where demixing and segregative phase separation occur
between DNA and RNA (essentially rRNA and mRNA involved in translation
complexes, but also the large amounts of rRNA synthesized at the rrn operons of
cells growing in rich media). It is also pointed out that crowding may play a
crucial role through its synergy with electrostatic forces. By constraining
macromolecules to remain at short distances from one another and feel
electrostatic interactions in spite of the strong screening exerted by
electrolyte species, crowding favors stronger DNA/RNA demixing and nucleoid
compaction.


  We report the folding thermodynamics of ccUUCGgg and ccGAGAgg RNA tetraloops
using atomistic molecular dynamics simulations. We obtain a previously
unreported estimation of the folding free energy using parallel tempering in
combination with well-tempered metadynamics. A key ingredient is the use of a
recently developed metric distance, eRMSD, as a biased collective variable. We
find that the native fold of both tetraloops is not the global free energy
minimum using the Amber\c{hi}OL3 force field. The estimated folding free
energies are 30.2kJ/mol for UUCG and 7.5 kJ/mol for GAGA, in striking
disagreement with experimental data. We evaluate the viability of all possible
one-dimensional backbone force field corrections. We find that disfavoring the
gauche+ region of {\alpha} and {\zeta} angles consistently improves the
existing force field. The level of accuracy achieved with these corrections,
however, cannot be considered sufficient by judging on the basis of available
thermodynamic data and solution experiments.


  The thermodynamic uncertainty relation provides an inequality relating any
mean current, the associated dispersion and the entropy production rate for
arbitrary non-equilibrium steady states. Applying it here to a general model of
a molecular motor running against an external force or torque, we show that the
thermodynamic efficiency of such motors is universally bounded by an expression
involving only experimentally accessible quantities. For motors pulling cargo
through a viscous fluid, a universal bound for the corresponding Stokes
efficiency follows as a variant. A similar result holds if mechanical force is
used to synthesize molecules of high chemical potential. Crucially, no
knowledge of the detailed underlying mechano-chemical mechanism is required for
applying these bounds.


  In order to model ion transport across protocell membranes in Hadean
hydrothermal vents, we consider both theoretically and experimentally the
planar growth of a precipitate membrane formed at the interface between two
parallel fluid streams in a two-dimensional microfluidic reactor. The growth
rate of the precipitate is found to be proportional to the square root of time,
which is characteristic of diffusive transport. However, the dependence of the
growth rate on the concentrations of hydroxide and metal ions is approximately
linear and quadratic, respectively. We show that such a difference in ionic
transport dynamics arises from the enhanced transport of metal ions across a
thin gel layer present at the surface of the precipitate. The fluctuations in
transverse velocity in this wavy porous gel layer allow an enhanced transport
of the cation, so that the effective diffusivity is about an order of magnitude
higher than that expected from molecular diffusion alone. Our theoretical
predictions are in excellent agreement with our laboratory measurements of the
growth of a manganese hydroxide membrane in a microfluidic channel, and this
enhanced transport is thought to have been needed to account for the
bioenergetics of the first single-celled organisms.


  The emergence of diseases such as Zika and Ebola has highlighted the need to
understand the role of sexual transmission in the spread of diseases with a
primarily non-sexual transmission route. In this paper we develop a number of
low-dimensional models which are appropriate for a range of assumptions for how
a disease will spread if it has sexual transmission through a sexual contact
network combined with some other transmission mechanism, such as direct contact
or vectors. The equations derived provide exact predictions for the dynamics of
the corresponding simulations in the large population limit.


  Intrinsically disordered proteins (IDPs) are typically low in
nonpolar/hydrophobic but relatively high in polar, charged, and aromatic amino
acid compositions. Some IDPs undergo liquid-liquid phase separation in the
aqueous milieu of the living cell. The resulting phase with enhanced IDP
concentration can function as a major component of membraneless organelles
that, by creating their own IDP-rich microenvironments, stimulate critical
biological functions. IDP phase behaviors are governed by their amino acid
sequences. To make progress in understanding this sequence-phase relationship,
we report further advances in a recently introduced application of
random-phase-approximation (RPA) heteropolymer theory to account for
sequence-specific electrostatics in IDP phase separation. Here we examine
computed variations in phase behavior with respect to block length and charge
density of model polyampholytes of alternating equal-length charge blocks to
gain insight into trends observed in IDP phase separation. As a real-life
example, the theory is applied to rationalize/predict binodal and spinodal
phase behaviors of the 236-residue N-terminal disordered region of RNA helicase
Ddx4 and its charge-scrambled mutant for which experimental data are available.
Fundamental differences are noted between the phase diagrams predicted by RPA
and those predicted by mean-field Flory-Huggins and
Overbeek-Voorn/Debye-H\"uckel theories. In the RPA context, a physically
plausible dependence of relative permittivity on protein concentration can
produce a cooperative effect in favor of IDP-IDP attraction and thus a
significant increased tendency to phase separate. Ramifications of these
findings for future development of IDP phase separation theory are discussed.


  The computer-aided folding of biomolecules, particularly RNAs, is one of the
most difficult challenges in computational structural biology. RNA tetraloops
are fundamental RNA motifs playing key roles in RNA folding and RNA-RNA and
RNA-protein interactions. Although state-of-the-art Molecular Dynamics (MD)
force fields correctly describe the native state of these tetraloops as a
stable free-energy basin on the microsecond time scale, enhanced sampling
techniques reveal that the native state is not the global free energy minimum,
suggesting yet unidentified significant imbalances in the force fields. Here,
we tested our ability to fold the RNA tetraloops in various force fields and
simulation settings. We employed three different enhanced sampling techniques,
namely, temperature replica exchange MD (T-REMD), replica exchange with solute
tempering (REST2), and well-tempered metadynamics (WT-MetaD). We aimed to
separate problems caused by limited sampling from those due to force-field
inaccuracies. We found that none of the contemporary force fields is able to
correctly describe folding of the 5'-GAGA-3' tetraloop over a range of
simulation conditions. We thus aimed to identify which terms of the force field
are responsible for this poor description of TL folding. We showed that at
least two different imbalances contribute to this behavior, namely,
overstabilization of base-phosphate and/or sugar-phosphate interactions and
underestimated stability of the hydrogen bonding interaction in base pairing.
The first artifact stabilizes the unfolded ensemble, while the second one
destabilizes the folded state. The former problem might be partially alleviated
by reparametrization of the van der Waals parameters of the phosphate oxygens
suggested by Case et al., while in order to overcome the latter effect we
suggest local potentials to better capture hydrogen bonding interactions.


  Numerous studies have attempted to model the effect of mass media on the
transmission of diseases such as influenza, however quantitative data on media
engagement has until recently been difficult to obtain. With the recent
explosion of "big data" coming from online social media and the like, large
volumes of data on a population's engagement with mass media during an epidemic
are becoming available to researchers. In this study we combine an online data
set comprising millions of shared messages relating to influenza with
traditional surveillance data on flu activity to suggest a functional form for
the relationship between the two. Using this data we present a simple
deterministic model for influenza dynamics incorporating media effects, and
show that such a model helps explain the dynamics of historical influenza
outbreaks. Furthermore, through model selection we show that the proposed media
function fits historical data better than other media functions proposed in
earlier studies.


  We consider a cyclically competing species model on a ring with global mixing
at finite rate, which corresponds to the well-known Lotka-Volterra equation in
the limit of infinite mixing rate. Within a perturbation analysis of the model
from the infinite mixing rate, we provide analytical evidence that extinction
occurs deterministically at sufficiently large but finite values of the mixing
rate for any species number $N\ge3$. Further, by focusing on the cases of
rather small species numbers, we discuss numerical results concerning the
trajectories toward such deterministic extinction, including global
bifurcations caused by changing the mixing rate.


  The discerning behavior of living systems relies on accurate interactions
selected from the lot of molecular collisions occurring in the cell. To ensure
the reliability of interactions, binding partners are classically envisioned as
finely preadapted molecules, evolutionarily selected on the basis of their
affinity in one-step associations. But the counterselection of inappropriate
interactions can in fact be much more efficiently obtained through difficult
multi-step adjustment, whose final high energy state is locked by a fluctuation
ratchet. The progressive addition of molecular bonds during stereo-adjustment
can be modeled as a predominantly backward random walk whose first arrival is
frozen by a micro-irreversible transition. A new criterion of ligand
specificity is presented, that is based on the ratio rejection/incorporation.
In addition to its role in the selectivity of interactions, this generic recipe
can underlie other important biological phenomena such as the regular synthesis
at low level of supramolecular complexes, monostable kinetic bimodality,
substrate concentration thresholds or the preparation of rapidly
depolymerizable structureswith stored energy, like microtubules.


  A long standing debate in biochemistry is to determine whether the
conformational changes observed during biomolecular interactions proceed
through conformational selection (of preexisting isoforms) or induced fit
(ligand-induced 3D reshaping). The latter mechanism had been invoked in certain
circumstances, for example to explain the non-Michaelian activity of monomeric
enzymes like glucokinase. But the relative importance of induced fit has been
recently depreciated in favor of conformational selection, assumed to be always
sufficient, predominant in general and in particular for glucokinase. The
relative contributions of conformational selection and induced fit are
reconsidered here in and out of equilibrium, in the light of earlier concepts
such as the cyclic equilibrium rule and the turning wheel of Wyman, using
single molecule state probability, one way fluxes and net fluxes. The
conditions for a switch from conformational selection to induced fit at a given
ligand concentration are explicitly determined. Out of equilibrium, the
inspection of the enzyme states circuit shows that conformational selection
alone would give a Michaelian reaction rate but not the established nonlinear
behaviour of glucokinase. Moreover, when induced fit and conformational
selection coexist and allow kinetic cooperativity, the net flux emerging in the
linkage cycle necessarily corresponds to the induced fit path.


  Epigenetics is most often reduced to chromatin marking in the current
literature, whereas this notion was initially defined in a more general
context. This restricted view ignores that epigenetic memories are in fact more
robustly ensured in living systems by steady state mechanisms with permanent
molecule renewal. This misconception is likely to result from misleading
intuitions and insufficient dialogues between traditional and quantitative
biologists. To demystify dynamic epigenetics, its most famous image, a
Waddington landscape and its attractors, are explicitly drawn. The simple
example provided is sufficient to highlight the main requirements and
characteristics of dynamic gene networks, underlying cellular differentiation,
de-differentiation and trans-differentiation.


  Networks of contacts capable of spreading infectious diseases are often
observed to be highly heterogeneous, with the majority of individuals having
fewer contacts than the mean, and a significant minority having relatively very
many contacts. We derive a two-dimensional diffusion model for the full
temporal behavior of the stochastic susceptible-infectious-recovered (SIR)
model on such a network, by making use of a time-scale separation in the
deterministic limit of the dynamics. This low-dimensional process is an
accurate approximation to the full model in the limit of large populations,
even for cases when the time-scale separation is not too pronounced, provided
the maximum degree is not of the order of the population size.


  The susceptible-infected-susceptible (SIS) agent-based model is usually
employed in the investigation of epidemics. The model describes a Markov
process for a single communicable disease among susceptible (S) and infected
(I) agents. However, the disease spreading forecasting is often restricted to
numerical simulations, while analytic formulations lack both general results
and perturbative approaches since they are subjected to asymmetric time
generators. Here, we discuss perturbation theory, approximations and
application of many-body techniques in epidemic models in the framework for
squared norm of probability vector $|P(t)| ^2$, in which asymmetric time
generators are replaced by their symmetric counterparts.


  Range expansion and range shifts are crucial population responses to climate
change. Genetic consequences are not well understood but are clearly coupled to
ecological dynamics that, in turn, are driven by shifting climate conditions.
We model a population with a deterministic reaction-- diffusion model coupled
to a heterogeneous environment that develops in time due to climate change. We
decompose the resulting travelling wave solution into neutral genetic
components to analyse the spatio-temporal dynamics of its genetic structure.
Our analysis shows that range expansions and range shifts under slow climate
change preserve genetic diversity. This is because slow climate change creates
range boundaries that promote spatial mixing of genetic components.
Mathematically , the mixing leads to so-called pushed travelling wave
solutions. This mixing phenomenon is not seen in spatially homogeneous
environments, where range expansion reduces genetic diversity through gene
surfing arising from pulled travelling wave solutions. However, the
preservation of diversity is diminished when climate change occurs too quickly.
Using diversity indices, we show that fast expansions and range shifts erode
genetic diversity more than slow range expansions and range shifts. Our study
provides analytical insight into the dynamics of travelling wave solutions in
heterogeneous environments.


  The data for the Ebola outbreak that occurred in 2014-2016 in three countries
of West Africa are analysed within a common framework. The analysis is made
using the results of an agent based Susceptible-Infected-Removed (SIR) model on
a Euclidean network, where nodes at a distance $l$ are connected with
probability $P(l) \propto l^{-\delta }$, $\delta$ determining the range of the
interaction, in addition to nearest neighbors. The cumulative (total) density
of infected population here has the form $R(t) =
\frac{a\exp(t/T)}{1+c\exp(t/T)}$, where the parameters depend on $\delta$ and
the infection probability $q$. This form is seen to fit well with the data.
Using the best fitting parameters, the time at which the peak is reached is
estimated and is shown to be consistent with the data. We also show that in the
Euclidean model, one can choose $\delta$ and $q$ values which reproduce the
data for the three countries qualitatively. These choices are correlated with
population density, control schemes and other factors. Comparing the real data
and the results from the model one can also estimate the size of the actual
population susceptible to the disease. Rescaling the real data a reasonably
good quantitative agreement with the simulation results is obtained.


  We consider evolution of a large population, where fitness of each organism
is defined by many phenotypical traits. These traits result from expression of
many genes. We propose a new model of gene regulation, where gene expression is
controlled by a gene network with a threshold mechanism and there is a feedback
between that threshold and gene expression. We show that this regulation is
very powerful: depending on parameters we can obtain any functional connection
between thresholds and genes. Under general assumptions on fitness we prove
that such model organisms are capable, to some extent, to recognize the fitness
landscape. That fitness landscape learning sharply reduces the number of
mutations necessary for adaptation and thus accelerates of evolution. Moreover,
this learning increases phenotype robustness with respect to mutations.
However, this acceleration leads to an additional risk since learning procedure
can produce errors. Finally evolution acceleration reminds races on a rugged
highway: when you speed up, you have more chances to crash. These results
explain recent experimental data on anticipation of environment changes by some
organisms.


  In recent years, a number of methods have been developed to infer complex
demographic histories, especially historical population size changes, from
genomic sequence data. Coalescent Hidden Markov Models have proven to be
particularly useful for this type of inference. Due to the Markovian structure
of these models, an essential building block is the joint distribution of local
genealogical trees, or statistics of these genealogies, at two neighboring loci
in populations of variable size. Here, we present a novel method to compute the
marginal and the joint distribution of the total length of the genealogical
trees at two loci separated by at most one recombination event for samples of
arbitrary size. To our knowledge, no method to compute these distributions has
been presented in the literature to date. We show that they can be obtained
from the solution of certain hyperbolic systems of partial differential
equations. We present a numerical algorithm, based on the method of
characteristics, that can be used to efficiently and accurately solve these
systems and compute the marginal and the joint distributions. We demonstrate
its utility to study the properties of the joint distribution. Our flexible
method can be straightforwardly extended to handle an arbitrary fixed number of
recombination events, to include the distributions of other statistics of the
genealogies as well, and can also be applied in structured populations.


  Bovine tuberculosis (TB) poses a serious threat for agricultural industry in
several countries, it involves potential interactions between wildlife and
cattle and creates societal problems in terms of human-wildlife conflict. This
study addresses connectedness network analysis, the spatial, and temporal
dynamics of TB between cattle in farms and the European badger (Meles meles)
using a large dataset generated by a calibrated agent based model. Results
showed that infected network connectedness was lower in badgers than in cattle.
The contribution of an infected individual to the mean distance of disease
spread over time was considerably lower for badger than cattle; badgers mainly
spread the disease locally while cattle infected both locally and across longer
distances. The majority of badger-induced infections occurred when individual
badgers leave their home sett, and this was positively correlated with badger
population growth rates. Point pattern analysis indicated aggregation in the
spatial pattern of TB prevalence in badger setts across all scales. The spatial
distribution of farms that were not TB free was aggregated at different scales
than the spatial distribution of infected badgers and became random at larger
scales. The spatial cross correlation between infected badger setts and
infected farms revealed that generally infected setts and farms do not coexist
except at few scales. Temporal autocorrelation detected a two year infection
cycle for badgers, while there was both within the year and longer cycles for
infected cattle. Temporal cross correlation indicated that infection cycles in
badgers and cattle are negatively correlated. The implications of these results
for understanding the dynamics of the disease are discussed.


  Many different proteins self-aggregate into insoluble fibrils growing
apically by reversible addition of elementary building blocks. But beyond this
common principle, the modalities of fibril formation are very disparate, with
various intermediate forms which can be reshuffled by minor modifications of
physico-chemical conditions or amino-acid sequences. To bypass this complexity,
the multifaceted phenomenon of fibril formation is reduced here to its most
elementary principles defined for a linear prototype of fibril. Selected
generic features, including nucleation, elongation and conformational
recruitment, are modeled using minimalist hypotheses and tools, by separating
equilibrium from kinetic aspects and in vitro from in vivo conditions. These
reductionist approaches allow to bring out known and new rudiments, including
the kinetic and equilibrium effects of nucleation, the dual influence of
elongation on nucleation, the kinetic limitations on nucleation and fibril
numbers and the accumulation of complexes in vivo by rescue from degradation.
Overlooked aspects of these processes are also pointed: the exponential
distribution of fibril lengths can be recovered using various models because it
is attributable to randomness only. It is also suggested that the same term
"critical concentration" is used for different things, involved in either
nucleation or elongation.


  Reconstructing who infected whom is a central challenge in analysing
epidemiological data. Recently, advances in sequencing technology have led to
increasing interest in Bayesian approaches to inferring who infected whom using
genetic data from pathogens. The logic behind such approaches is that isolates
that are nearly genetically identical are more likely to have been recently
transmitted than those that are very different. A number of methods have been
developed to perform this inference. However, testing their convergence,
examining posterior sets of transmission trees and comparing methods'
performance are challenged by the fact that the object of inference - the
transmission tree - is a complicated discrete structure. We introduce a metric
on transmission trees to quantify distances between them. The metric can
accommodate trees with unsampled individuals, and highlights differences in the
source case and in the number of infections per infector. We illustrate its
performance on simple simulated scenarios and on posterior transmission trees
from a TB outbreak. We find that the metric reveals where the posterior is
sensitive to the priors, and where collections of trees are composed of
distinct clusters. We use the metric to define median trees summarising these
clusters. Quantitative tools to compare transmission trees to each other will
be required for assessing MCMC convergence, exploring posterior trees and
benchmarking diverse methods as this field continues to mature.


  We consider a continuous-time Markov chain model of SIR disease dynamics with
two levels of mixing. For this so-called stochastic households model, we
provide two methods for inferring the model parameters---governing
within-household transmission, recovery, and between-household
transmission---from data of the day upon which each individual became
infectious and the household in which each infection occurred, as would be
available from first few hundred studies. Each method is a form of Bayesian
Markov Chain Monte Carlo that allows us to calculate a joint posterior
distribution for all parameters and hence the household reproduction number and
the early growth rate of the epidemic. The first method performs exact Bayesian
inference using a standard data-augmentation approach; the second performs
approximate Bayesian inference based on a likelihood approximation derived from
branching processes. These methods are compared for computational efficiency
and posteriors from each are compared. The branching process is shown to be an
excellent approximation and remains computationally efficient as the amount of
data is increased.


  To model discriminative, i.e. competition induced, self-thinning in even-aged
forest stands a concept has been explored that discriminative mortality alters
spatial arrangement of trees which in turn alters the mortality. Function of
density was suggested to be a balance of initial density minus mortality that
was dependent on initial spatial arrangement and mean horizontal size in the
forest stand. Assuming initial spatial arrangement to be aggregation and
performing normalizations gave the self-thinning function depending on only two
parameters, initial stem density and maximal horizontal size (represented by
stem diameter). Normalizations included integral of probability of trees to
collide to be unity and stem density at maximal size to be zero. The
self-thinning function obtained has been shown to successfully capture
non-linear self-thinning dynamics in Douglas-fir long-term experiments in which
competition-induced mortality prevailed.


  Recently, there has been a lot of effort in understanding sRNA mediated
regulation of gene expression and how this mode of regulation differs from
transcriptional regulation.In E.coli, in the presence of oxidative stress, the
synthesis of sigma^s is regulated through an interesting mechanism involving
both transcriptional and sRNA-mediated translational regulation. The key
regulatory factors involved in transcriptional and translational regulation are
ArcA and ArcB proteins and ArcZ sRNA, respectively. Phosphorylated ArcA, in a
feedforward mechanism, represses the transcriptions sigma^s and ArcZ sRNA with
the latter being a post-transcriptional activator for sigma^s. Through a
feedback mechanisms, ArcZ sRNA destabilises ArcB mRNA and regulates the level
of ArcB protein, a kinase that phosphorylates ArcA. The oxygen and energy
availability is expected to influence the ArcA phosphorylation rate and, as a
consequence, in equilibrium, the system is likely to be in either a high ArcB
(low ArcZ) or a low ArcB (high ArcZ) state. Kinetic modelling studies suggest
that the rate of destabilisation of ArcB mRNA by ArcZ sRNA must be
appropriately tuned for achieving the desired state. In particular, for a high
phosphorylation rate, the transition from a low to a high ArcZ synthesis regime
with the increase in sRNA-mRNA interaction is similar to the threshold-linear
response observed earlier. Further, we show that the mRNA destabilisation by
sRNA might be, in particular, beneficial in the low phosphorylation state for
having the right concentration levels of ArcZ and ArcB. Stochastic simulation
results suggest that as the ArcZ-ArcB binding affinity is increased, the
probability distribution for the number of ArcZ molecules becomes flatter
indicating frequently occurring transcriptional bursts of varying strengths.


  We present an algebraic solution for the Susceptible-Infective-Removed (SIR)
model originally presented by Kermack-McKendrick in 1927. Starting from the
differential equation for the removed subjects presented by them in the
original paper, we re-write it in a slightly different form in order to derive
formally the solution, unless one integration. Then, using algebraic techniques
and some well justified numerical assumptions we obtain an analytic solution
for the integral. Finally, we compare the numerical solution of the
differential equations of the SIR model with the analytically solution here
proposed, showing an excellent agreement.


  Understanding how species are distributed across landscapes over time is a
fundamental question in biodiversity research. Unfortunately, most species
distribution models only target a single species at a time, despite strong
ecological evidence that species are not independently distributed. We propose
Deep Multi-Species Embedding (DMSE), which jointly embeds vectors corresponding
to multiple species as well as vectors representing environmental covariates
into a common high-dimensional feature space via a deep neural network. Applied
to bird observational data from the citizen science project \textit{eBird}, we
demonstrate how the DMSE model discovers inter-species relationships to
outperform single-species distribution models (random forests and SVMs) as well
as competing multi-label models. Additionally, we demonstrate the benefit of
using a deep neural network to extract features within the embedding and show
how they improve the predictive performance of species distribution modelling.
An important domain contribution of the DMSE model is the ability to discover
and describe species interactions while simultaneously learning the shared
habitat preferences among species. As an additional contribution, we provide a
graphical embedding of hundreds of bird species in the Northeast US.


  In this communication we introduce an efficient implementation of adaptive
biasing that greatly improves the speed of free energy computation in molecular
dynamics simulations. We investigated the use of accelerated simulations to
inform on compound design using a recently reported and clinically relevant
inhibitor of the chromatin regulator BRD4. Benchmarking on our local compute
cluster, our implementation achieves up to 2.5 times more force calls per day
than plumed2. Results of five 1{\mu}second-long simulations are presented,
which reveal a conformational switch in the BRD4 inhibitor between a binding
competent and incompetent state. Stabilization of the switch led to a -3
kcal/mol improvement of absolute binding free energy. These studies suggest an
unexplored ligand design principle and offer new actionable hypotheses for
medicinal chemistry efforts against this druggable epigenetic target class.


  Accurate computational identification of promoters remains a challenge as
these key DNA regulatory regions have variable structures composed of
functional motifs that provide gene specific initiation of transcription. In
this paper we utilize Convolutional Neural Networks (CNN) to analyze sequence
characteristics of prokaryotic and eukaryotic promoters and build their
predictive models. We trained the same CNN architecture on promoters of four
very distant organisms: human, plant (Arabidopsis), and two bacteria
(Escherichia coli and Mycoplasma pneumonia). We found that CNN trained on
sigma70 subclass of Escherichia coli promoter gives an excellent classification
of promoters and non-promoter sequences (Sn=0.90, Sp=0.96, CC=0.84). The
Bacillus subtilis promoters identification CNN model achieves Sn=0.91, Sp=0.95,
and CC=0.86. For human and Arabidopsis promoters we employ CNNs for
identification of two well-known promoter classes (TATA and non-TATA
promoters). CNNs models nicely recognize these complex functional regions. For
human Sn/Sp/CC accuracy of prediction reached 0.95/0.98/0,90 on TATA and
0.90/0.98/0.89 for non-TATA promoter sequences, respectively. For Arabidopsis
we observed Sn/Sp/CC 0.95/0.97/0.91 (TATA) and 0.94/0.94/0.86 (non-TATA)
promoters. Thus, the developed CNN models (implemented in CNNProm program)
demonstrated the ability of deep learning with grasping complex promoter
sequence characteristics and achieve significantly higher accuracy compared to
the previously developed promoter prediction programs. As the suggested
approach does not require knowledge of any specific promoter features, it can
be easily extended to identify promoters and other complex functional regions
in sequences of many other and especially newly sequenced genomes. The CNNProm
program is available to run at web server http://www.softberry.com.


  The presence of membrane tubules in living cells is essential to many
biological processes. In cells, one mechanism to form nanosized lipid tubules
is via molecular motor induced bilayer extraction. In this paper, we describe a
simple experimental model to investigate the forces required for lipid tube
formation using kinesin motors anchored to DOPC vesicles. Previous related
studies have used molecular motors actively pulling on the membrane to extract
a nanotube. Here, we invert the system geometry; molecular motors are used as
static anchors linking DOPC vesicles to a two-dimensional microtubule network
and an external flow is introduced to generate nanotubes facilitated by the
drag force. We found that a drag force of ~7 pN was sufficient for tubule
extraction for vesicles ranging from 1 to 2 um in radius. By our method, we
found that the force generated by a single molecular motor was sufficient for
membrane tubule extraction from a spherical lipid vesicle.


  Aggregation of amphiphiles through the action of hydrophobic interactions is
a common feature in soft condensed matter systems and is of particular
importance in the context of biophysics as it underlies both the generation of
functional biological machinery as well as the formation of pathological
misassembled states of proteins. Here we explore the aggregation behaviour of
amphiphilic polymers using lattice Monte-Carlo calculations and show that the
distribution of hydrophobic residues within the polymer sequence determines the
facility with which dry/wet interfaces can be created and that such interfaces
drive the aggregation process.


  Within living cells, the transport of cargo is accomplished by groups of
molecular motors. Such collective transport could utilize mechanisms which
emerge from inter-motor interactions in ways that are yet to be fully
understood. Here we combined experimental measurements of two-kinesin transport
with a theoretical framework to investigate the functional ramifications of
inter-motor interactions on individual motor function and collective cargo
transport. In contrast to kinesin's low sidestepping frequency when present as
a single motor, with exactly two kinesins per cargo, we observed substantial
motion perpendicular to the microtubule. Our model captures a
surface-associated mode of kinesin, which is only accessible via inter-motor
interference in groups, in which kinesin diffuses along the microtubule surface
and rapidly "hops" between protofilaments without dissociating from the
microtubule. Critically, each kinesin transitions dynamically between the
active stepping mode and this weak surface-associated mode enhancing local
exploration of the microtubule surface, possibly enabling cellular cargos to
overcome macromolecular crowding and to navigate obstacles along microtubule
tracks without sacrificing overall travel distance.


  Native extracellular matrices (ECMs), such as those of the human brain and
other neural tissues, exhibit networks of molecular interactions between
specific matrix proteins and other tissue components. Guided by these naturally
self-assembling supramolecular systems, we have designed a matrix-derived
protein chimera that contains a laminin globular-like (LG) domain fused to an
elastin-like polypeptide (ELP). All-atom, classical molecular dynamics
simulations of our designed laminin-elastin fusion protein reveal
temperature-dependent conformational changes, in terms of secondary structure
composition, solvent accessible surface area, hydrogen bonding, and surface
hydration. These properties illuminate the phase behavior of this fusion
protein, via the emergence of $\beta$-sheet character in
physiologically-relevant temperature ranges.


  There are two contrasting explanations of sleep: as a proximate, essential
physiological function or as an adaptive state of inactivity and these
hypotheses remain widely debated. To investigate the adaptive significance of
sleep, we develop an evolutionary argument formulated as a tractable partial
differential equation model. We allow demographic parameters such as birth and
mortality rates to vary through time in both safe and vulnerable sleeping
environments. From this model we analytically calculate population growth rate
(fitness) for sleeping and non-sleeping strategies. We find that, in a
temporally heterogeneous environment, sleeping always achieves a higher fitness
than not sleeping. As organisms do not exist in constant environments, we
conclude that the evolution of sleep is inevitable.


  Recent genome and transcriptome sequencing projects have unveiled a plethora
of highly structured RNA molecules as central mediators of cellular function.
Single molecule Forster Resonance Energy Transfer (smFRET) is a powerful tool
for analyzing the temporal evolution of the global structure of individual RNA
molecules, in pursuit of understanding their essential
structure-dynamics-function relationships. In contrast to enzymatic and
chemical footprinting, NMR spectroscopy and X-ray crystallography, smFRET
yields temporally resolved, quantitative information about single molecules
rather than only time and ensemble averages of entire populations. This enables
unique observations of transient and rare conformations under both equilibrium
and non-equilibrium conditions.


  It is known that steric blocking by peptide sidechains of hydrogen bonding,
HB, between water and peptide groups, PGs, in beta sheets accords with an amino
acid intrinsic beta sheet preference. The present observations with Quantum
Molecular Dynamics, QMD, simulation with quantum mechanical treatment of every
water molecule solvating a beta sheet that would be transient in nature suggest
that this steric blocking is not applicable in a hydrophobic region unless a
cation is present, so that the amino acid beta sheet preference due to this
steric blocking is only effective in the presence of a cation. We observed
backbone hydration in a polyalanine and to a lesser extent polyvaline alpha
helix without a cation being present, but a cation could increase the strength
of these HBs.
  Parallel beta sheets have a greater tendency than antiparallel beta sheets of
equivalent small size to retain regular structure in solvated QMD, and a 4
strand 4 inter-PG HB chain parallel beta sheet was used. Stability was
reinforced by one surface being polyvaline, which buttressed the opposite
surface which was used for experimentation. A single Ca2+ ion was used for
investigation of individual binding events rather than bulk properties. No
direct binding between Ca2+ and the PG oxygen was observed in these
simulations, but perhaps it occurs at longer time scales as the transient beta
sheet unfolds.
  When linear scaling QMD methods that are accurate for peptide resonance,
Resonance-Assisted Hydrogen Bonding and the properties of water become
available, more extensive experiments having multiple ions of multiple types
could be performed at acceptable computational cost. It important that such
investigations be performed on protein secondary structures rather than model
amides so that sidechain limitation of backbone hydration and hence intrinsic
amino acid propensity is captured.


  Determining the three-dimensional structure of proteins and protein complexes
at atomic resolution is a fundamental task in structural biology. Over the last
decade, remarkable progress has been made using "single particle" cryo-electron
microscopy (cryo-EM) for this purpose. In cryo-EM, hundreds of thousands of
two-dimensional images are obtained of individual copies of the same particle,
each held in a thin sheet of ice at some unknown orientation. Each image
corresponds to the noisy projection of the particle's electron-scattering
density. The reconstruction of a high-resolution image from this data is
typically formulated as a nonlinear, non-convex optimization problem for
unknowns which encode the angular pose and lateral offset of each particle.
Since there are hundreds of thousands of such parameters, this leads to a very
CPU-intensive task---limiting both the number of particle images which can be
processed and the number of independent reconstructions which can be carried
out for the purpose of statistical validation.
  Here, we propose a deterministic method for high-resolution reconstruction
that operates in an ab initio manner---that is, without the need for an initial
guess. It requires a predictable and relatively modest amount of computational
effort, by marching out radially in the Fourier domain from low to high
frequency, increasing the resolution by a fixed increment at each step.


  Romania's forest cover is at this time below, but quite close to the European
average. Despite that, recent forestry activities tend to be more oriented
towards forest exploitation and not to increase the national forested area. The
general public perception is that afforestation activities are limited only to
NGO's projects and media actions and numerous reports and statements made by
some of these organizations are rather sensational and confusing. This paper
tries to cast more light on this controversial issue, presenting accurate data
and considering the whole context of afforestation-from the history of
afforestation activities in Romania to the recent statements and reports
submitted by the government agencies involved in these particular issues. The
necessity and opportunity of afforestation is also substantiated from the
perspective of reintroduction into production of marginal and degraded land. An
analysis of the situation of funding through the National Rural Development
Programme (2007-2013) measures is performed and possible future directions
regarding afforestation programmes are discussed.


  Network theoretic approach has been used to model and study the flow of
ecological information, growth and connectivity on landscape level of
anemochory plant species Abied pindrow, Betula utilis and Taxus wallichiana in
the Western Himalaya region. A network is formally defined and derived for seed
dispersion model of aforementioned species where vertices represent habitat
patches which are connected by an edge if the distance between the patches is
less than a threshold distance. We define centrality of a network and
computationally identify the habitat patches that are central to the process of
seed dispersion to occur across the network. We find that the network of
habitat patches is a scale free network and at the same time it also displays
small world property characterized by high clustering and low average shortest
path length. Due to high clustering, the spread of species is locally even as
seed disperse mutually among the member vertices of a cluster. Also since every
vertex is only a short number of steps away from every other vertex, the
species rapidly covers all the habitat patches in the component. Also due to
presence of hubs in the network the spread of species is greatly boosted
whenever the species establish and thrive in a hub patch and disperse to
adjacent patches. However, the network is not modular due to geographical
constraints, and is is negatively assortative as the high degree vertices are
connected to vertices of low degree.


  Most of the existing research in assembly pathway prediction/analysis of
virus cap- sids makes the simplifying assumption that the configuration of the
intermediate states can be extracted directly from the final configuration of
the entire capsid. This assump- tion does not take into account the
conformational changes of the constituent proteins as well as minor changes to
the binding interfaces that continues throughout the assembly process until
stabilization. This paper presents a statistical-ensemble based approach which
provides sufficient samples of the configurational space for each monomer and
the relative local orientation between monomers, to capture the uncertainties
in their binding and conformations. Furthermore, instead of using larger
capsomers (trimers, pentamers) as building blocks, we allow all possible
sub-assemblies to bind in all pos- sible combinations. We represent this
assembly graph in two different ways. First, we use the Wilcoxon signed rank
measure to compare the distributions of binding free energy computed on the
sampled conformations to predict likely pathways. Second, we represent chemical
equilibrium aspects of the transitions as a Bayesian Factor graph where both
associations and dissociations are modeled based on concentrations and the
binding free energies. Results from both of these experiments showed
significant departure from those one would obtain if only the static
configurations of the proteins were considered. Hence, we establish the
importance of an uncertainty-aware protocol for pathway analysis, and provide a
statistical framework as an important first step towards assembly pathway
prediction with high statistical confidence.


  Selective autophagy must not only select the correct type of organelle, but
also must discriminate between individual organelles of the same kind so that
some but not all of the organelles are removed. We propose that physical
clustering of autophagy receptor proteins on the organelle surface can provide
an appropriate all-or-none signal for organelle degradation. We explore this
proposal using a computational model restricted to peroxisomes and the
relatively well characterized pexophagy receptor proteins NBR1 and p62. We find
that larger peroxisomes nucleate NBR1 clusters first and lose them last through
competitive coarsening. This results in significant size-selectivity that
favors large peroxisomes, and can explain the increased catalase signal that
results from siRNA inhibition of p62. Excess ubiquitin, resulting from damaged
organelles, suppresses size-selectivity but not cluster formation. Our proposed
selectivity mechanism thus allows all damaged organelles to be degraded, while
otherwise selecting only a portion of organelles for degradation.


  Information transmission in biological signaling circuits has often been
described using the metaphor of a noise filter. Cellular systems need accurate,
real-time data about their environmental conditions, but the biochemical
reaction networks that propagate, amplify, and process signals work with noisy
representations of that data. Biology must implement strategies that not only
filter the noise, but also predict the current state of the environment based
on information delayed due to the finite speed of chemical signaling. The idea
of a biochemical noise filter is actually more than just a metaphor: we
describe recent work that has made an explicit mathematical connection between
signaling fidelity in cellular circuits and the classic theories of optimal
noise filtering and prediction that began with Wiener, Kolmogorov, Shannon, and
Bode. This theoretical framework provides a versatile tool, allowing us to
derive analytical bounds on the maximum mutual information between the
environmental signal and the real-time estimate constructed by the system. It
helps us understand how the structure of a biological network, and the response
times of its components, influences the accuracy of that estimate. The theory
also provides insights into how evolution may have tuned enzyme kinetic
parameters and populations to optimize information transfer.


  One of the long-expected goals of genome-scale metabolic modeling is to
evaluate the influence of the perturbed enzymes to the flux distribution. Both
ordinary differential equation (ODE) models and the constraint-based models,
like Flux balance analysis (FBA), lack of the room of performing metabolic
control analysis (MCA) for large-scale networks. In this study, we developed a
Hyper-Cube Shrink Algorithm (HCSA) to incorporate the enzymatic properties to
the FBA model by introducing a pseudo reaction constrained by enzymatic
parameters. Our algorithm was able to handle not only prediction of knockout
strains but also strains with quantitative adjustment of expression level or
activity. We first demonstrate the concept by applying HCSA to a simplest
three-node network. Then we validate its prediction by comparing with ODE and
with a synthetic network in Saccharomyces cerevisiae producing voilacein and
analogues. Finally we show its capability of predicting the flux distribution
in genome-scale networks by applying it to the sporulation in yeast.


  We present a model for the dynamics of a population of bacteria with a
continuum of traits, who compete for resources and exchange horizontally
(transfer) an otherwise vertically inherited trait with possible mutations.
Competition influences individual demographics, affecting population size,
which feeds back on the dynamics of transfer. We consider a stochastic
individual-based pure jump process taking values in the space of point
measures, and whose jump events describe the individual reproduction, transfer
and death mechanisms. In a large population scale, the stochastic process is
proved to converge to the solution of a nonlinear integro-differential
equation. When there are only two different traits and no mutation, this
equation reduces to a non-standard two-dimensional dynamical system. We show
how crucial the forms of the transfer rates are for the long-term behavior of
its solutions. We describe the dynamics of invasion and fixation when one of
the two traits is initially rare, and compute the invasion probabilities. Then,
we study the process under the assumption of rare mutations. We prove that the
stochastic process at the mutation time scale converges to a jump process which
describes the successive invasions of successful mutants. We show that the
horizontal transfer can have a major impact on the distribution of the
successive mutational fixations, leading to dramatically different behaviors,
from expected evolution scenarios to evolutionary suicide. Simulations are
given to illustrate these phenomena.


  With their longest dimension typically being less than 100 nm, molecular
motors are significantly below the optical-resolution limit. Despite
substantial advances in fluorescence-based imaging methodologies, labeling with
beads remains critical for optical-trapping-based investigations of molecular
motors. A key experimental challenge in bead-based assays is that the number of
motors on a bead is not well defined. Particularly for single-molecule
investigations, the probability of single versus multiple-motor events has not
been experimentally investigated. Here, we used bead travel distance as an
indicator of multiple-motor transport and determined the lower-bound
probability of bead transport by two or more motors. We limited the ATP
concentration to increase our detection sensitivity for multiple- versus
single-kinesin transport. Surprisingly, for all but the lowest motor number
examined, our measurements exceeded estimations of a previous model by R2-fold.
To bridge this apparent gap between theory and experiment, we derived a
closed-form expression for the probability of bead transport by multiple
motors, and constrained the only free parameter in this model using our
experimental measurements. Our data indicate that kinesin extends to ~57 nm
during bead transport, suggesting that kinesin exploits its conformational
flexibility to interact with microtubules at highly curved interfaces such as
those present for vesicle transport in cells. To our knowledge, our findings
provide the first experimentally constrained guide for estimating the
probability of multiple-motor transport in optical trapping studies. The
experimental approach utilized here (limiting ATP concentration) may be
generally applicable to studies in which molecular motors are labeled with
cargos that are artificial or are purified from cellular extracts.


  Microtubules are protein polymers that form "molecular highways" for
long-range transport within living cells. Molecular motors actively step along
microtubules to shuttle cellular materials between the nucleus and the cell
periphery; this transport is critical for the survival and health of all
eukaryotic cells. Structural defects in microtubules exist, but whether these
defects impact molecular motor-based transport remains unknown. Here, we report
a new, to our knowledge, approach that allowed us to directly investigate the
impact of such defects. Using a modified optical-trapping method, we examined
the group function of a major molecular motor, conventional kinesin, when
transporting cargos along individual microtubules. We found that microtubule
defects influence kinesin-based transport in vitro. The effects depend on motor
number: cargos driven by a few motors tended to unbind prematurely from the
microtubule, whereas cargos driven by more motors tended to pause. To our
knowledge, our study provides the first direct link between microtubule defects
and kinesin function. The effects uncovered in our study may have physiological
relevance in vivo.


  Although the disease-relevant microtubule-associated protein tau is known to
severely inhibit kinesin-based transport in vitro, the potential mechanisms for
reversing this detrimental effect to maintain healthy transport in cells remain
unknown. Here we report the unambiguous upregulation of multiple-kinesin travel
distance despite the presence of tau, via decreased single-kinesin velocity.
Interestingly, the presence of tau also modestly reduced cargo velocity in
multiple-kinesin transport, and our stochastic simulations indicate that the
tau-mediated reduction in single-kinesin travel underlies this observation.
Taken together, our observations highlight a nontrivial interplay between
velocity and travel distance for kinesin transport, and suggest that
single-kinesin velocity is a promising experimental handle for tuning the
effect of tau on multiple-kinesin travel distance.


  Microtubule-based molecular motors often work in small groups to transport
cargos in cells. A key question in understanding transport (and its regulation
in vivo) is to identify the sensitivity of multiple-motor-based motion to
various single molecule properties. Whereas both single-motor travel distance
and microtubule binding rate have been demonstrated to contribute to cargo
travel, the role of single-motor velocity is yet to be explored. Here, we
recast a previous theoretical study, and make explicit a potential contribution
of velocity to cargo travel. We test this possibility experimentally, and
demonstrate a strong negative correlation between single-motor velocity and
cargo travel for transport driven by two motors. Our study thus discovers a
previously unappreciated role of single-motor velocity in regulating
multiple-motor transport.


  We analyze the governing partial differential equations of a model of
pole-to-pole oscillations of the MinD protein in a bacterial cell. The
sensitivity to extrinsic noise in the parameters of the model is explored. Our
analysis shows that overall, the oscillations are robust to extrinsic
perturbations in the sense that small perturbations in reaction coefficients
result in small differences in the frequency and in the amplitude. However, a
combination of analysis and simulation also reveals that the oscillations are
more sensitive to some extrinsic time-scales than to others.


  The VHHs are antigen-binding region/domain of camelid heavy chain antibodies
(HCAb). They have many interesting biotechnological and biomedical properties
due to their small size, high solubility and stability, and high affinity and
specificity for their antigens. HCAb and classical IgGs are evolutionary
related and share a common fold. VHHs are composed of regions considered as
constant, called the frameworks (FRs) connected by Complementarity Determining
Regions (CDRs), a highly variable region that provide interaction with the
epitope. Actually, no systematic structural analyses had been performed on VHH
structures despite a significant number of structures. This work is the first
study to analyse the structural diversity of FRs of VHHs. Using a structural
alphabet that allows approximating the local conformation, we show that each of
the four FRs do not have a unique structure but exhibit many structural variant
patterns. Moreover, no direct simple link between the local conformational
change and amino acid composition can be detected. These results indicate that
long-range interactions affect the local conformation of FRs and impact the
building of structural models.


  Antibiotic-resistant bacteria has posed a grave threat to public health by
causing a number of nosocomial infections in hospitals. Mathematical models
have been used to study the transmission of antibiotic-resistant bacteria
within a hospital and the measures to control antibiotic resistance in
nosocomial pathogens. Studies presented in \cite{LBL,LB} have shown great value
in understanding the transmission of antibiotic-resistant bacteria in a
hospital. However, their results are limited to numerical simulations of a few
different scenarios without analytical analysis of the models in all
biologically feasible parameter regions. Bifurcation analysis and
identification of the global stability conditions are necessary to assess the
interventions which are proposed to limit nosocomial infection and stem the
spread of antibiotic-resistant bacteria. In this paper we study the global
dynamics of the mathematical model of antibiotic resistance in hospitals in
\cite{LBL,LB}. The invasion reproduction number $\mathcal R_{ar}$ of
antibiotic-resistant bacteria is introduced. We give the relationship of
$\mathcal R_{ar}$ and two control reproduction numbers of sensitive bacteria
and resistant bacteria ($\mathcal R_{sc}$ and $\mathcal R_{rc}$). More
importantly, we prove that a backward bifurcation may occur at $\mathcal
R_{ar}=1$ when the model includes superinfection which is not mentioned in
\cite{LB}. That is, there exists a new threshold $\mathcal R_{ar}^c$, and if
$\mathcal R_{ar}^c<\mathcal R_{ar}<1$, then the system can have two interior
equilibria and it supports an interesting bistable phenomenon. This provides
critical information on controlling the antibiotic-resistance in a hospital.


  In this research, methods and computational results based on statistical
analysis and mathematical modelling, data collection in situ in order to make a
hazard map of Hanta Virus infection in the region of Araucania, Chile are
presented. The development of this work involves several elements such as
Landsat satellite images, biological information regarding seropositivity of
Hanta Virus and information concerning positive cases of infection detected in
the region. All this information has been processed to find a function that
models the danger of contagion in the region, through logistic regression
analysis and Artificial Neural Networks


  We have recently introduced Forman's discretization of Ricci curvature to the
realm of complex networks. Forman curvature is an edge-based measure whose
mathematical definition elegantly encapsulates the weights of nodes and edges
in a complex network. In this contribution, we perform a comparative analysis
of Forman curvature with other edge-based measures such as edge betweenness,
embeddedness and dispersion in diverse model and real networks. We find that
Forman curvature in comparison to embeddedness or dispersion is a better
indicator of the importance of an edge for the large-scale connectivity of
complex networks. Based on the definition of the Forman curvature of edges,
there are two natural ways to define the Forman curvature of nodes in a
network. In this contribution, we also examine these two possible definitions
of Forman curvature of nodes in diverse model and real networks. Based on our
empirical analysis, we find that in practice the unnormalized definition of the
Forman curvature of nodes with the choice of combinatorial node weights is a
better indicator of the importance of nodes in complex networks.


  The phylogenetic network emerged in the 1990s as a new model to represent the
evolution of species in the case where coexisting species transfer genetic
information through hybridization, recombination, lateral gene transfer, etc.
As is true for many rapidly evolving fields, there is considerable
fragmentation and diversity in methodologies, standards and vocabulary in
phylogenetic network research, thus creating the need for an integrated
database of articles, authors, techniques, keywords and software. We describe
such a database, "Who is Who in Phylogenetic Networks", available at
http://phylnet.univ-mlv.fr. "Who is Who in Phylogenetic Networks" comprises
more than 600 publications and 500 authors interlinked with a rich set of more
than 200 keywords related to phylogenetic networks. The database is integrated
with web-based tools to visualize authorship and collaboration networks and
analyze these networks using common graph and social network metrics such as
centrality (betweenness, eigenvector, degree and closeness) and clustering. We
provide downloads of raw information about entries in the database, and a
facility to suggest modifications and contribute new information to the
database. We also present in this article common use cases of the database and
identify trends in the research on phylogenetic networks using the information
in the database and textual analysis.


  We examined the sensitivity of microtubule spools to transport velocity.
Perhaps surprisingly, we determined that the steady-state number and size of
spools remained constant over a seven-fold range of velocities. Our data on the
kinetics of spool assembly further suggest that the main mechanisms underlying
spool growth vary during assembly.