  Managing biodiversity extinction crisis requires wise decision-making
processes able to account for the limited of resources available. In most
decision problems in conservation biology, several conflicting objectives have
to be taken into account. Most methods used in conservation either provide
suboptimal solutions or use strong assumptions about the decision-maker's
preferences. Our paper reviews some of the existing approaches to solve
multi-objective decision problems and presents new multi-objective linear
programming formulations of two multi-objective optimization problems in
conservation, allowing the use of a reference point approach. Reference point
approaches solve multi-objective optimization problems by interactively
representing the preferences of the decision-maker with a point in the criteria
space, called the reference point. We modelled and solved the following two
conservation problems: a dynamic multi-species management problem under
uncertainty and a spatial allocation resource management problem. Results show
that the reference point method outperforms classic methods while illustrating
the use of an interactive methodology for solving combinatorial problems with
multiple objectives. The method is general and can be adapted to a wide range
of ecological combinatorial problems.


  The notion that transcription factors bind DNA only through specific,
consensus binding sites has been recently questioned. In a pioneering study by
Pugh and Venters no specific consensus motif for the positioning of the human
pre-initiation complex (PIC) has been identified. Here, we reveal that
nonconsensus, statistical, DNA triplet code provides specificity for the
positioning of the human PIC. In particular, we reveal a highly non-random,
statistical pattern of repetitive nucleotide triplets that correlates with the
genome-wide binding preferences of PIC measured by Chip-exo. We analyze the
triplet enrichment and depletion near the transcription start site (TSS) and
identify triplets that have the strongest effect on PIC-DNA nonconsensus
binding. Our results constitute a proof-of-concept for a new design principle
for protein-DNA recognition in the human genome, which can lead to a better
mechanistic understanding of transcriptional regulation.


  Establishing a quantitative connection between the population growth rate and
the generation times of single cells is a prerequisite for understanding
evolutionary dynamics of microbes. However, existing theories fail to account
for the experimentally observed correlations between mother-daughter generation
times that are unavoidable when cell size is controlled for - which is
essentially always the case. Here, we study population-level growth in the
presence of cell size control and corroborate our theory using experimental
measurements of single-cell growth rates. We derive a closed formula for the
population growth rate and demonstrate that it only depends on the single-cell
growth rate variability, not other sources of stochasticity. Our work provides
an evolutionary rationale for the narrow growth rate distributions often
observed in nature: when single-cell growth rates are less variable but have a
fixed mean, the population will exhibit an enhanced population growth rate, as
long as the correlations between the mother and daughter cells' growth rates
are not too strong.


  Studies on the role of fluctuations in signal propagation and on gene
regulation in monoclonal bacterial population have been extensively pursued
based on the machinery of two-component system. The bacterial two-component
system shows noise utilisation through its inherent plasticity. The
fluctuations propagation takes place using the phosphotransfer module and the
feedback mechanism during gene regulation. To delicately observe the noisy
kinetics the generic cascade needs stochastic investigation at the mRNA and
protein levels. To this end, we propose a theoretical framework to investigate
the noisy signal transduction in a generic two-component system. The model
shows reliability in information transmission through quantification of several
statistical measures. We further extend our analysis to observe the protein
distribution in a population of cells. Through numerical simulation, we
identify the regime of the kinetic parameter set that generates a stability
switch in the steady state distribution of proteins. The results of our
theoretical analysis show key features of the network. The noise permeation and
information propagation in the autoregulation module is feeble. However, the
phosphotransfer module compensates such weakness and plays a significant role
in information transmission. The bimodality due to fluctuations pampers the
emergence of persistence in an isogenic bacterial pool.


  We propose to develop mean field theory in combination with Glauber
algorithm, to model and interpret protein dynamics and structure formation in
small to wide angle x-ray scattering (S/WAXS) experiments. We develop the
methodology by analysing the Engrailed homeodomain protein as an example. We
demonstrate how to interpret S/WAXS data with a good precision and over an
extended temperature range. We explain experimentally observed phenomena in
terms of protein phase structure, and we make predictions for future
experiments how the scattering data behaves at different ambient temperature
values. We conclude that a combination of mean field theory with Glauber
algorithm has the potential to develop into a highly accurate, computationally
effective and predictive tool for analysing S/WAXS data. Finally, we compare
our results with those obtained previously in an all-atom molecular dynamics
simulation.


  Penguins are flightless, so they are forced to walk while on land. In
particular, they show rather specific behaviors in their homecoming, which are
interesting to observe and to describe analytically. In this paper, we present
a simple mathematical formulation to describe the little penguins parade in
Phillip Island. We observed that penguins have the tendency to waddle back and
forth on the shore to create a sufficiently large group and then walk home
compactly together. The mathematical framework that we introduce describes this
phenomenon, by taking into account "natural parameters" such as the eye-sight
of the penguins, their cruising speed and the possible "fear" of animals. On
the one hand, this favors the formation of conglomerates of penguins that
gather together, but, on the other hand, this may lead to the "panic" of
isolated and exposed individuals. The model that we propose is based on a set
of ordinary differential equations. Due to the discontinuous behavior of the
speed of the penguins, the mathematical treatment (to get existence and
uniqueness of the solution) is based on a "stop-and-go" procedure. We use this
setting to provide rigorous examples in which at least some penguins manage to
safely return home (there are also cases in which some penguins freeze due to
panic). To facilitate the intuition of the model, we also present some simple
numerical simulations that can be compared with the actual movement of the
penguins parade.


  The Susceptible-Infected-Recovered (SIR) model is studied in multilayer
networks with arbitrary number of links across the layers. By following the
mapping to bond percolation we give the analytical expression for the epidemic
threshold and the fraction of the infected individuals in arbitrary number of
layers. These results provide an exact prediction of the epidemic threshold for
infinite locally tree-like multilayer networks, and an lower bound of the
epidemic threshold for more general multilayer networks. The case of a
multilayer network formed by two interconnected networks is specifically
studied as a function of the degree distribution within and across the layers.
We show that the epidemic threshold strongly depends on the degree correlations
of the multilayer structure. Finally we relate our results to the results
obtained in the annealed approximation for the Susceptible-Infected-Susceptible
(SIS) model.


  We explain that Stegosaurus exhibited exterior chirality and observe that the
largest plate in particular of USNM 4394, USNM 4714, DMNS 2818 and NHMUK R36730
appears to have tilted to the right rather than to the left in each case.
Several instances in which Stegosaurus specimens have been confused with their
distinct, hypothetical mirror-image forms are highlighted. We believe our
findings to be consistent with the hypothesis that Stegosaurus's plates acted
primarily as display structures. A collection of more than one Stegosaurus
might be referred to henceforth as a 'handful' of Stegosaurus.


  Gonadotropin-releasing hormone (GnRH) is reported to control mammalian
reproductive processes. GnRH a neurohormone which is pulsatile released into
the pituitary portal blood by hypothalamic GnRH neurons. In the present study,
the phase synchronization among a population of identical neurons subjected to
a pool of coupling molecules GnRH in extracellular medium via mean-field
coupling mechanism is investigated. In the model of populated neurons, GnRH is
considered to be autocrine signaling molecule and is taken to be common to all
neurons to act as synchronizing agent. The rate of synchrony is estimated
qualitatively and quantitatively by measuring phase locking values, time
evolution of the phase differences and recurrence plots. Our numerical results
show a phase transition like behavior separating the synchronized and
desynchronized regimes. We also investigated long range communication or relay
information transfer for one dimensional array of such neurons.


  Alloreactivity following stem cell transplantation (SCT) is difficult to
predict in patients undergoing transplantation from HLA matched donors. In this
study we performed whole exome sequencing of SCT donor-recipient pairs (DRP).
This allowed determination of entire library of alloreactive peptide sequences
which would bind HLA class I molecules in each DRP. Utilizing the HLA binding
affinity (IC50) and tissue expression levels of the parent proteins, an
aggregate donor T cell response to the recipient alloreactive peptides was
calculated using a vector-operator dynamical system model. Marked variability
in the simulated CD8+ T cell responses was observed in all the donor recipient
pairs.


  How animals allocate energy to different body functions is still not
completely understood and a challenging topic until recently. Here, we
investigate in more detail the allocation of energy intake to growth,
reproduction or heat production by developing energy budget models for
ectothermic and endothermic vertebrates using a mathematical approach. We
calculated energy intake functions of ectotherms and endotherms derived from
their body mass growth. We show that our energy budget model produces energy
intake patterns and distributions as observed in ectothermic and endothermic
species. Our results comply consistently with some empirical studies that in
endothermic species, like birds and mammals, energy is used for heat production
instead of growth. Our model additionally offers an explanation on known
differences in absolute energy intake between ectothermic fish and reptiles and
endothermic birds and mammals. From a mathematical point of view, the model
comes in two equivalent formulations, a differential and an integral one. It is
derived from a discrete level approach, and it is shown to be well-posed and to
attain a unique solution for (almost) every parameter set. Numerically, the
integral formulation of the model is considered as an inverse problem with
unknown parameters that are estimated using a series of experiments/realistic
data.


  Players of coevolutionary games may update not only their strategies but also
their networks of interaction. Based on interpreting the payoff of players as
fitness, dynamic landscape models are proposed. The modeling procedure is
carried out for Prisoner's Dilemma (PD) and Snowdrift (SD) games that both use
either birth--death (BD) or death--birth (DB) strategy updating. The main focus
is on using dynamic fitness landscapes as a mathematical model of
coevolutionary game dynamics. Hence, an alternative tool for analyzing
coevolutionary games becomes available, and landscape measures such as
modality, ruggedness and information content can be computed and analyzed. In
addition, fixation properties of the games and quantifiers characterizing the
interaction networks are calculated numerically. Relations are established
between landscape properties expressed by landscape measures and quantifiers of
coevolutionary game dynamics such as fixation probabilities, fixation times and
network properties.


  In this paper, a delay differential equations (DDEs) model of leukemia is
introduced and its dynamical properties are investigated in comparison with the
modified fractional-order system where the Caputo's derivative is used. The
model takes into account three types of division that a stem-like cell can
undergo and cell competition between healthy and leukemia cell populations. The
action of the immune system on the leukemic cell populations is also
considered. The stability properties of the equilibrium points are established
through numerical results and the differences between the two types of
approaches are discussed. Medical conclusions are drawn in view of the obtained
numerical simulations.


  The modeling of complex reaction-diffusion processes in, for instance,
cellular biochemical networks or self-assembling soft matter can be
tremendously sped up by employing a multiscale algorithm which combines the
mesoscopic Green's Function Reaction Dynamics (GFRD) method with explicit
stochastic Brownian, Langevin, or deterministic Molecular Dynamics to treat
reactants at the microscopic scale [A. Vijaykumar, P.G. Bolhuis and P.R. ten
Wolde, J. Chem. Phys. {\bf 43}, 21: 214102 (2015)]. Here we extend this
multiscale BD-GFRD approach to include the orientational dynamics that is
crucial to describe the anisotropic interactions often prevalent in
biomolecular systems. We illustrate the novel algorithm using a simple patchy
particle model. After validation of the algorithm we discuss its performance.
The rotational BD-GFRD multiscale method will open up the possibility for large
scale simulations of e.g. protein signalling networks.


  The changes on abiotic features of ecosystems have rarely been taken into
account by population dynamics models, which typically focus on trophic and
competitive interactions between species. However, understanding the population
dynamics of organisms that must modify their habitats in order to survive, the
so-called ecosystem engineers, requires the explicit incorporation of abiotic
interactions in the models. Here we study a model of ecosystem engineers that
is discrete both in space and time, and where the engineers and their habitats
are arranged in patches fixed to the sites of regular lattices. The growth of
the engineer population is modeled by Ricker equation with a density-dependent
carrying capacity that is given by the number of modified habitats. A diffusive
dispersal stage ensures that a fraction of the engineers move from their birth
patches to neighboring patches. We find that dispersal influences the
metapopulation dynamics only in the case that the local or single-patch
dynamics exhibits chaotic behavior. In that case, it can suppress the chaotic
behavior and avoid extinctions in the regime of large intrinsic growth rate of
the population.


  Self-assembly of protein monomers into distinct membrane protein oligomers
provides a general mechanism for diversity in the molecular architectures, and
resulting biological functions, of membrane proteins. We develop a general
physical framework describing the thermodynamic competition between different
oligomeric states of membrane proteins. Using the mechanosensitive channel of
large conductance as a model system, we show how the dominant oligomeric states
of membrane proteins emerge from the interplay of protein concentration in the
cell membrane, protein-induced lipid bilayer deformations, and direct
monomer-monomer interactions. Our results suggest general physical mechanisms
and principles underlying regulation of protein function via control of
membrane protein oligomeric state.


  In this work we show that there is a direct relationship between a graph's
topology and the free energy of a spin system on the graph. We develop a method
of separating topological and enthalpic contributions to the free energy, and
find that considering the topology is sufficient to qualitatively compare the
free energies of different graph systems at high temperature, even when the
energetics are not fully known. This method was applied to the metal lattice
system with defects, and we found that it partially explains why point defects
are more stable than high-dimensional defects. Given the energetics, we can
even quantitatively compare free energies of different graph structures via a
closed form of linear graph contributions. The closed form is applied to
predict the sequence space free energy of lattice proteins, which is a key
factor determining the designability of a protein structure.


  Single molecule detection, sequencing and conformational mapping of aptamers
are important for improving medical and biosensing technologies and for better
understanding of biological processes at the molecular level. We obtain
vibrational signals of single aptamers immobilized on gold substrates using
tip-enhanced Raman spectroscopy (TERS). We compare topographic and optical
signals and investigate the fluctuations of the position-dependent TERS
spectra. TERS mapping provides information about the chemical composition and
conformation of aptamers, and paves the way to future single-molecule
label-free sequencing.


  It is well known that many biochemical processes in the cell such as gene
regulation, growth signals and activation of ion channels, rely on mechanical
stimuli. However, the mechanism by which mechanical signals propagate through
cells is not as well understood. In this review we focus on stress propagation
in a minimal model for cell elasticity, actomyosin networks, which are
comprised of a sub-family of cytoskeleton proteins. After giving an overview of
th actomyosin network components, structure and evolution we review stress
propagation in these materials as measured through the correlated motion of
tracer beads. We also discuss the possibility to extract structural features of
these networks from the same experiments. We show that stress transmission
through these networks has two pathways, a quickly dissipative one through the
bulk, and a long ranged weakly dissipative one through the pre-stressed actin
network.


  The dynamics of two competing species in a finite size community is one of
the most studied problems in population genetics and community ecology.
Stochastic fluctuations lead, inevitably, to the extinction of one of the
species, but the relevant timescale depends on the underlying dynamics. The
persistence time of the community has been calculated for neutral models, where
the only drive of the system is drift (demographic stochasticity) and for
models with strong selection. Following recent analyses that stress the
importance of environmental stochasticity in empirical systems, we present here
a general theory of persistence time of two-species community where drift,
environmental variations and time independent selective advantage are all taken
into account.


  Predators may attack isolated or grouped prey in a cooperative, collective
way. Whether a gregarious behavior is advantageous to each species depends on
several conditions and game theory is a useful tool to deal with such a
problem. We here extend the Lett-Auger-Gaillard model [Theor. Pop. Biol. 65
(2004) 263] to spatially distributed populations and compare the resulting
behavior with their mean-field predictions for the coevolving densities of
predator and prey strategies. Besides its richer behavior in the presence of
spatial organization, we also show that the coexistence phase in which
collective and individual strategies for each group are present is stable
because of an effective, cyclic dominance mechanism similar to a well-studied
generalization of the Rock-Paper-Scissors game with four species, a further
example of how ubiquitous this coexistence mechanism is.


  Characterization of B-cell protein epitope and developing critical parameters
for its identification is one of the long standing interests. Using Layers
algorithm, we introduced the concept of anchor residues to identify epitope. We
have shown that majority of the epitope is composed of anchor residues and have
significant bias in epitope for these residues. We optimized the search space
reduction for epitope identification. We used Layers to non-randomly sample the
antigen surface reducing the molecular surface to an average of 75 residues
while preserving 50% of the epitope in the sample surface. To facilitate the
comparison of favorite methods of researchers we compared the popular
techniques used to identify epitope with their sampling performance and
evaluation. We proposed an optimum Sr of 16 {\AA} to sample the antigen
molecules to reduce the search space, in which epitope is identified using
buried surface area method. We used the combinations of molecular surface
sampling, anchor residue intensity in surface, secondary structure and sequence
information to predict epitope at an accuracy of 89%. A web application is made
available at
http://www.csb.iitkgp.ernet.in/applications/b_cell_epitope_pred/main.


  We apply a Gaussian variational approximation to model reduction in large
biochemical networks of unary and binary reactions. We focus on a small subset
of variables (subnetwork) of interest, e.g. because they are accessible
experimentally, embedded in a larger network (bulk). The key goal is to write
dynamical equations reduced to the subnetwork but still retaining the effects
of the bulk. As a result, the subnetwork-reduced dynamics contains a memory
term and an extrinsic noise term with non-trivial temporal correlations. We
first derive expressions for this memory and noise in the linearized (Gaussian)
dynamics and then use a perturbative power expansion to obtain first order
nonlinear corrections. For the case of vanishing intrinsic noise, our
description is explicitly shown to be equivalent to projection methods up to
quadratic terms, but it is applicable also in the presence of stochastic
fluctuations in the original dynamics. An example from the Epidermal Growth
Factor Receptor (EGFR) signalling pathway is provided to probe the increased
prediction accuracy and computational efficiency of our method.


  The reconstruction of phylogenetic trees from discrete character data
typically relies on models that assume the characters evolve under a
continuous-time Markov process operating at some overall rate $\lambda$. When
$\lambda$ is too high or too low, it becomes difficult to distinguish a short
interior edge from a polytomy (the tree that results from collapsing the edge).
In this note, we investigate the rate that maximizes the expected
log-likelihood ratio (i.e. the Kullback--Leibler separation) between the
four-leaf unresolved (star) tree and a four-leaf binary tree with interior edge
length $\epsilon$. For a simple two-state model, we show that as $\epsilon$
converges to $0$ the optimal rate also converges to zero when the four pendant
edges have equal length. However, when the four pendant branches have unequal
length, two local optima can arise, and it is possible for the globally optimal
rate to converge to a non-zero constant as $\epsilon \rightarrow 0$. Moreover,
in the setting where the four pendant branches have equal lengths and either
(i) we replace the two-state model by an infinite-state model or (ii) we retain
the two-state model and replace the Kullback--Leibler separation by Euclidean
distance as the maximization goal, then the optimal rate also converges to a
non-zero constant.


  The pathogenesis and progression of many tumors, including hematologic
malignancies is highly dependent on enhanced lipogenesis. De novo fatty-acid
synthesis permits accelerated proliferation of tumor cells by providing
structural components to build the membranes. It may also lead to alterations
of physicochemical properties of the formed membranes, which can have an impact
on signaling or even increase resistance to drugs in cancer cells. Cancer
type-specific lipid profiles would allow understanding the actual effects of
lipid changes and therefore could potentially serve as fingerprints for
individual tumors and be explored as diagnostic markers. We have used shotgun
MS approach to identify lipid patterns in different types of acute myeloid
leukemia (AML) patients that either show no karyotype changes or belong to
t(8;21) or inv16 types. The observed differences in lipidomes of t(8;21) and
inv(16) patients, as compared to AML patients without karyotype changes,
concentrate mostly on substantial modulation of ceramides/sphingolipids
synthesis. Also significant changes in the physicochemical properties of the
membranes, between the t(8;21) and the other patients, were noted that were
related to a marked alteration in the saturation levels of lipids. The revealed
differences in lipid profiles of various AML types increase our understanding
of the affected biochemical pathways and can potentially serve as diagnostic
tools.


  Context-free and context-sensitive formal grammars are often regarded as more
appropriate to model proteins than regular level models such as finite state
automata and Hidden Markov Models. In theory, the claim is well founded in the
fact that many biologically relevant interactions between residues of protein
sequences have a character of nested or crossed dependencies. In practice,
there is hardly any evidence that grammars of higher expressiveness have an
edge over old good HMMs in typical applications including recognition and
classification of protein sequences. This is in contrast to RNA modeling, where
CFG power some of the most successful tools. There have been proposed several
explanations of this phenomenon. On the biology side, one difficulty is that
interactions in proteins are often less specific and more "collective" in
comparison to RNA. On the modeling side, a difficulty is the larger alphabet
which combined with high complexity of CF and CS grammars imposes considerable
trade-offs consisting on information reduction or learning sub-optimal
solutions. Indeed, some studies hinted that CF level of expressiveness brought
an added value in protein modeling when CF and regular grammars where
implemented in the same framework. However, there have been no systematic study
of explanatory power provided by various grammatical models. The first step to
this goal is define objective criteria of such evaluation. Intuitively, a
decent explanatory grammar should generate topology, or the parse tree,
consistent with topology of the protein, or its secondary and/or tertiary
structure. In this piece of research we build on this intuition and propose a
set of measures to compare topology of the parse tree of a grammar with
topology of the protein structure.


  Defeat and success of the competitive invasion of a populated area is
described with a standard Lotka-Volterra competition model. The resident is
adapted to the heterogeneous living conditions, i.e., its motion is modelled as
space-dependent, so-called Fokker-Planck diffusion. The invader's diffusion is
taken as neutral Fickian. Furthermore, it is studied how multiplicative
environmental noise fosters or hinders the invasion.


  We propose a stochastic model in evolutionary game theory where individuals
(or subpopulations) can mutate changing their strategies randomly (but rarely)
and explore the external environment. This environment affects the selective
pressure by modifying the payoff arising from the interactions between
strategies. We derive a Fokker-Plank integro-differential equation and provide
Monte-Carlo simulations for the Hawks vs Doves game. In particular we show
that, in some cases, taking into account the external environment favors the
persistence of the low-fitness strategy.


  There is general agreement that the former (breeding) macrohabitat of the
Asian tiger mosquito, Aedes albopictus (synonym: Stegomyia albopicta) has been
described as phytotelmata in the forested areas of Southeast Asia, however,
this has changed in the last four decades as it has adapted to more urban areas
and antrotelmata. Capable of producing eggs with a certain dry resistance as
well as cold hardiness, populations of Ae. albopictus became distributed around
the globe. The invasion of the Ae. albopictus across all over the world is
thought to be the most rapid spread of any insect species in the last four
decades in tropical, subtropical and temperate climate zones. In addition, Ae.
albopictus is a potential vector for at least 27 viruses as well as for several
parasites and plays a major role in the global transmission of dengue virus and
chikungunya virus; and its contribution to the rapid spread of zika virus is
currently under investigation. Therefore, this species is considered a serious
threat to public health worldwide. This work broaches the issue of three
investigations about the potential to adapt and establish stabile populations
of the invasive Asian tiger mosquito Ae. albopictus. As a result of these
investigations, a high potential for this species to become a threat to public
health in many more countries, especially in temperate climatic zones, can be
seen. Due to the fact on the one hand that Ae. albopictus can gain fitness
advantages due to misapplied vector control measures and on the other hand of
the high epigenetically adaptation potential, it is in summary recommended to
focus further research on the development of vaccinations for viruses and other
pathogens. By doing so, citizens will be protected without putting ecosystems
and their services in danger and it is the even more economic solution.


  Forests are complex systems, and it is necessary to include this
characteristic in every forest definition, in order to consider the restriction
that this imposes in terms of prediction and control. This lost of
predictability and controllability should be incorporated in every
Environmental Impact Assessment or management program. We present two
case-studies located in Mexico and one in the US to illustrate three relevant
indicators of complexity. First, we introduce an informational framework to
measure the Zoquiapan forest systemic complexity. Then, we analyze complexity
changes among different types of forest and management systems, related with
spatial distributions, using data from a floristic study in the Montes Azules
National Park. Finally, we analyze time series of $CO_{2}$ fluctuations taken
from AMERIFLUX data bases. Our results show firstly that it is possible to
measure the systemic complexity of different forests, characteized by a
criticality state (1/f noise) which has been proposed as a finger print of
complexity. And secondly, that this characteristic can be used as a proxy of
their state of conservation, where the lowest complexity values are found in
perturbed areas showing the relevance of the concept and its measurement for
forest conservation and management.


  Abundance estimates from animal point-count surveys require accurate
estimates of detection probabilities. The standard model for estimating
detection from removal-sampled point-count surveys assumes that organisms at a
survey site are detected at a constant rate; however, this assumption is often
not justified. We consider a class of N-mixture models that allows for
detection heterogeneity over time through a flexibly defined time-to-detection
distribution (TTDD) and allows for fixed and random effects for both abundance
and detection. Our model is thus a combination of survival time-to-event
analysis with unknown-N, unknown-p abundance estimation. We specifically
explore two-parameter families of TTDDs, e.g. gamma, that can additionally
include a mixture component to model increased probability of detection in the
initial observation period. We find that modeling a TTDD by using a
two-parameter family is necessary when data have a chance of arising from a
distribution of this nature. In addition, models with a mixture component can
outperform non-mixture models even when the truth is non-mixture. Finally, we
analyze an Overbird data set from the Chippewa National Forest using mixed
effect models for both abundance and detection. We demonstrate that the effects
of explanatory variables on abundance and detection are consistent across
mixture TTDDs but that flexible TTDDs result in lower estimated probabilities
of detection and therefore higher estimates of abundance.


  Artemisinin resistance constitutes a major threat to the continued success of
control programs for malaria. With alternative antimalarial drugs not yet
available, improving our understanding of how artemisinin-based drugs act and
how resistance manifests is essential to enable optimisation of dosing regimens
in order to prolong the lifespan of current first-line treatment options. Here,
through introduction of a novel model of the dynamics of the parasites'
response to drug, we explore how artemisinin-based therapies may be adjusted to
maintain efficacy and how artemisinin resistance may manifest and be overcome.
We introduce a dynamic mathematical model, extending on the traditional
pharmacokinetic-pharmacodynamic framework, to capture the time-dependent
development of a stress response in parasites. We fit the model to in vitro
data and establish that the parasites' stress response explains the recently
identified complex interplay between drug concentration, exposure time and
parasite viability. Our model demonstrates that the previously reported
hypersensitivity of early ring stage parasites of the 3D7 strain to
dihydroartemisinin (DHA) is primarily due to the rapid development of stress,
rather than any change in the maximum achievable killing rate. Of direct
clinical relevance, we demonstrate that the complex temporal features of
artemisinin action observed in vitro have a significant impact on predictions
of in vivo parasite clearance using PK-PD models. Given the important role that
such models play in the design and evaluation of clinical trials for
alternative drug dosing regimens, our model contributes an enhanced predictive
platform for the continued efforts to minimise the burden of malaria.


  The Y-chromosome degeneration is still an intriguing mechanism and comprises
the very origin of sex. We present a coupled version of the well known logistic
map and the logistic equation describing the evolution of XY chromosomes.
Although chaos was found in X, Y chromosomes do not evolve chaotically. A
mathematical constraint is shown as the responsible for this behaviour. In
addition, analytical solutions are presented for the differential equations
herein.


  Accurately predicting drug responses to cancer is an important problem
hindering oncologists' efforts to find the most effective drugs to treat
cancer, which is a core goal in precision medicine. The scientific community
has focused on improving this prediction based on genomic, epigenomic, and
proteomic datasets measured in human cancer cell lines. Real-world cancer cell
lines contain noise, which degrades the performance of machine learning
algorithms. This problem is rarely addressed in the existing approaches. In
this paper, we present a noise-filtering approach that integrates techniques
from numerical linear algebra and information retrieval targeted at filtering
out noisy cancer cell lines. By filtering out noisy cancer cell lines, we can
train machine learning algorithms on better quality cancer cell lines. We
evaluate the performance of our approach and compare it with an existing
approach using the Area Under the ROC Curve (AUC) on clinical trial data. The
experimental results show that our proposed approach is stable and also yields
the highest AUC at a statistically significant level.


  We propose a new mathematical model for the transmission dynamics of the
human immunodeficiency virus (HIV). Global stability of the unique endemic
equilibrium is proved. Then, based on data provided by the "Progress Report on
the AIDS response in Cape Verde 2015", we calibrate our model to the cumulative
cases of infection by HIV and AIDS from 1987 to 2014 and we show that our model
predicts well such reality. Finally, a sensitivity analysis is done for the
case study in Cape Verde. We conclude that the goal of the United Nations to
end the AIDS epidemic by 2030 is a nontrivial task.


  Compartments are ubiquitous throughout biology, yet their importance
stretches back to the origin of cells. In the context of origin of life, we
assume that a protocell, a compartment enclosing functional components,
requires $N$ components to be evolvable. We take interest in the timescale in
which a minimal evolvable protocell is produced. We show that when protocells
fuse and share information, the time to produce an evolvable protocell scales
algebraically in $N$, in contrast to an exponential scaling in the absence of
fusion. We discuss the implications of this result for origins of life, as well
as other biological processes.


  Selection, the tendency of some traits to become more frequent than others in
a population under the influence of some (natural or artificial) agency, is a
key component of Darwinian evolution and countless other natural and social
phenomena. Yet a general theory of selection, analogous to the
Fisher-Tippett-Gnedenko theory of extreme events, is lacking. Here we introduce
a probabilistic definition of selection and show that selected values are
attracted to a universal family of limiting distributions. The universality
classes and scaling exponents are determined by the tail thickness of the
random variable under selection. Our results are supported by data from
molecular biology, agriculture and sport.


  A global optimization method called Greedy Neighborhood Search (GNS) and a
novel conformational sampling method using a spherical distribution is proposed
to find the minimum energy conformation of a protein-like heteropolymer model
called AB model. The AB model consists of hydrophobic (A) and hydrophilic (B)
monomers analogous to the real proteins. The AB model in three-dimensional
space is represented by simple bead-rod chain system which is identical to the
one-bead protein model. The minimum energy conformations of four different
sequences consisting of 13, 21, 34, and 55 monomers are obtained by the GNS
method. The minimum energies found are lower than those obtained by other
methods. Also the minimum energy conformations found have a similarity with the
real proteins forming a single hydrophobic core.


  Antimicrobial resistance is an important public health concern that has
implications in the practice of medicine worldwide. Accurately predicting
resistance phenotypes from genome sequences shows great promise in promoting
better use of antimicrobial agents, by determining which antibiotics are likely
to be effective in specific clinical cases. In healthcare, this would allow for
the design of treatment plans tailored for specific individuals, likely
resulting in better clinical outcomes for patients with bacterial infections.
In this work, we present the recent work of Drouin et al. (2016) on using Set
Covering Machines to learn highly interpretable models of antibiotic resistance
and complement it by providing a large scale application of their method to the
entire PATRIC database. We report prediction results for 36 new datasets and
present the Kover AMR platform, a new web-based tool allowing the visualization
and interpretation of the generated models.


  Gene expression is a noisy process that leads to regime shift between
alternative steady states among individual living cells, inducing phenotypic
variability. The effects of white noise on the regime shift in bistable systems
have been well characterized, however little is known about such effects of
colored noise (noise with non-zero correlation time). Here, we show that noise
correlation time, by considering a genetic circuit of autoactivation, can have
significant effect on the regime shift in gene expression. We demonstrate this
theoretically, using stochastic potential, stationary probability density
function and first-passage time based on the Fokker-Planck description, where
the Ornstein-Uhlenbeck process is used to model colored noise. We find that
increase in noise correlation time in degradation rate can induce a regime
shift from low to high protein concentration state and enhance the bistable
regime, while increase in noise correlation time in basal rate retain the
bimodal distribution. We then show how cross-correlated colored noises in basal
and degradation rates can induce regime shifts from low to high protein
concentration state, but reduce the bistable regime. In addition, we show that
early warning indicators can also be used to predict shifts between distinct
phenotypic states in gene expression. Predictions that a cell is about to shift
to a harmful phenotype could improve early therapeutic intervention in complex
human diseases.


  Fission yeast G2/M transition is regulated by a biochemical reaction networks
which contains four components: Cdc13, Cdc2, Wee1, and Cdc25. This circuit is
characterized by the ultrasensitive responses of Wee1 or Cdc25 to Cdc13/Cdc2
activity, and the bistability of Cdc2 activation. Previous work has shown that
this bistability is governed by phosphorylation energy. In this article, we
developed the kinetic model of this circuit and conducted further thermodynamic
analysis on the role of phosphorylation energy (&[Delta]G). We showed that
level &[Delta]G shapes the response curves of Wee1 or Cdc25 to Cdc2 and governs
the intrinsic noise level of Cdc2 activity. More importantly, the mutually
antagonistic chemical fluxes around the PdP cycles in G2/M circuit were shown
to act as a stabilizer of Cdc2 activity against &[Delta]G fluctuations. These
results suggests the fundamental role of free energy and chemical fluxes on the
sensitivity, bistability and robustness of G2/M transition.


  Epidemic models are increasingly used in real-world networks to understand
diffusion phenomena (such as the spread of diseases, emotions, innovations,
failures) or the transport of information (such as news, memes in social
on-line networks). A new analysis of the prevalence, the expected number of
infected nodes in a network, is presented and physically interpreted. The
analysis method is based on spectral decomposition and leads to a universal,
analytic curve, that can bound the time-varying prevalence in any finite time
interval. Moreover, that universal curve also applies to various types of
Susceptible-Infected-Susceptible (SIS) (and Susceptible-Infected-Removed (SIR))
infection processes, with both homogenous and heterogeneous infection
characteristics (curing and infection rates), in temporal and even disconnected
graphs and in SIS processes with and without self-infections. The accuracy of
the universal curve is comparable to that of well-established mean-field
approximations.


  We propose to combine a mean field approach with all atom molecular dynamics
into a multistage algorithm that can model protein folding and dynamics over
very long time periods yet with atomic level precision. As an example we
investigate an isolated monomeric Myc oncoprotein that has been implicated in
carcinomas including those in colon, breast and lungs. Under physiological
conditions a monomeric Myc is presumed to be an example of intrinsically
disordered proteins, that pose a serious challenge to existing modelling
techniques. We argue that a room temperature monomeric Myc is in a dynamical
state, it oscillates between different conformations that we identify. For this
we adopt the C-alpha backbone of Myc in a crystallographic heteromer as an
initial Ansatz for the monomeric structure. We construct a multisoliton of the
pertinent Landau free energy, to describe the C-alpha profile with ultra high
precision. We use Glauber dynamics to resolve how the multisoliton responds to
repeated increases and decreases in ambient temperature. We confirm that the
initial structure is unstable in isolation. We reveal a highly degenerate
ground state landscape, an attractive set towards which Glauber dynamics
converges in the limit of vanishing ambient temperature. We analyse the thermal
stability of this Glauber attractor using room temperature molecular dynamics.
We identify and scrutinise a particularly stable subset in which the two
helical segments of the original multisoliton align in parallel, next to each
other. During the MD time evolution of a representative structure from this
subset, we observe intermittent quasiparticle oscillations along the C-terminal
alpha-helix, some of which resemble a translating Davydov's Amide-I soliton. We
propose that the presence of oscillatory motion is in line with the expected
intrinsically disordered character of Myc.


  Cancer has become one of the most widespread diseases in the world.
Specifically, breast cancer is diagnosed more often than any other type of
cancer. However, breast cancer patients and their individual tumors are often
unique. Identifying the underlying genetic phenotype can lead to precision
(personalized) medicine. Tailoring medical treatment strategies to best fit the
needs of individual patients can dramatically improve their health. Such an
approach requires sufficient knowledge of the patients and the diseases, which
is currently unavailable to practitioners. This study focuses on breast cancer
and proposes a novel two-stage clustering method to partition patients into
hierarchical groups. The first stage is broad grouping, which is based on
phenotypes such as demographic information and clinical features. The second
stage is fine grouping based on genomic characteristics, such as copy number
variation and somatic mutation, of patients in a subgroup resulting from the
first stage. Generally, this framework offers a mechanism to mix multiple forms
of data, both phenotypic and genomic, to most effectively define individual
patients for personalized predictions. This method provides the ability to
detect correlation among all factors.


  In addition to their unusually long life cycle, periodical cicadas, {\it
Magicicada} spp., provide an exceptional example of spatially synchronized life
stage phenology in nature. Within regions ("broods") spanning 50,000 to 500,000
km$^2$, adults emerge synchronously every 13 or 17 years. While satiation of
avian predators is believed to be a key component of the ability of these
populations to reach high densities, it is not clear why populations at a
single location remain entirely synchronized. We develop nonlinear Leslie
matrix-type models of periodical cicadas that include predation-driven Allee
effects and competition in addition to reproduction and survival. Using both
analytical and numerical techniques, we demonstrate the observed presence of a
single brood critically depends on the relationship between fecundity,
competition, and predation. We analyze the single-brood, two-brood and
all-brood equilibria in the large life-span limit using a tractable hybrid
approximation to the Leslie matrix model with continuous time competition in
between discrete reproduction events. Within the hybrid model we prove that the
single-brood equilibrium is the only stable equilibrium. This hybrid model
allows us to quantitatively predict population sizes and the range of
parameters for which the stable single-brood and unstable two-brood and
all-brood equilibria exist. The hybrid model yields a good approximation to the
numerical results for the Leslie matrix model for the biologically relevant
case of a 17-year lifespan.


  Stochasticity can play an important role in the dynamics of biologically
relevant populations. These span a broad range of scales: from intra-cellular
populations of molecules to population of cells and then to groups of plants,
animals and people. Large deviations in stochastic population dynamics - such
as those determining population extinction, fixation or switching between
different states - are presently in a focus of attention of statistical
physicists. We review recent progress in applying different variants of
dissipative WKB approximation (after Wentzel, Kramers and Brillouin) to this
class of problems. The WKB approximation allows one to evaluate the mean time
and/or probability of population extinction, fixation and switches resulting
from either intrinsic (demographic) noise, or a combination of the demographic
noise and environmental variations, deterministic or random. We mostly cover
well-mixed populations, single and multiple, but also briefly consider
populations on heterogeneous networks and spatial populations. The spatial
setting also allows one to study large fluctuations of the speed of biological
invasions. Finally, we briefly discuss possible directions of future work.


  Cross-protection, which refers to a process whereby artificially inoculating
a plant with a mild strain provides protection against a more aggressive
isolate of the virus, is known to be an effective tool of disease control in
plants. In this paper we derive and analyse a new mathematical model of the
interactions between two competing viruses with particular account for RNA
interference. Our results show that co-infection of the host can either
increase or decrease the potency of individual infections depending on the
levels of cross-protection or cross-enhancement between different viruses.
Analytical and numerical bifurcation analyses are employed to investigate the
stability of all steady states of the model in order to identify parameter
regions where the system exhibits synergistic or antagonistic behaviour between
viral strains, as well as different types of host recovery. We show that not
only viral attributes but also the propagating component of RNA-interference in
plants can play an important role in determining the dynamics.


  RNA interference (RNAi) is a fundamental cellular process that inhibits gene
expression through cleavage and destruction of target mRNA. It is responsible
for a number of important intracellular functions, from being the first line of
immune defence against pathogens to regulating development and morphogenesis.
In this paper we consider a mathematical model of RNAi with particular emphasis
on time delays associated with two aspects of primed amplification: binding of
siRNA to aberrant RNA, and binding of siRNA to mRNA, both of which result in
the expanded production of dsRNA responsible for RNA silencing. Analytical and
numerical stability analyses are performed to identify regions of stability of
different steady states and to determine conditions on parameters that lead to
instability. Our results suggest that while the original model without time
delays exhibits a bi-stability due to the presence of a hysteresis loop, under
the influence of time delays, one of the two steady states with the high
(default) or small (silenced) concentration of mRNA can actually lose its
stability via a Hopf bifurcation. This leads to the co-existence of a stable
steady state and a stable periodic orbit, which has a profound effect on the
dynamics of the system.


  Recent theoretical work has demonstrated that Neighbor Joining applied to
concatenated DNA sequences is a statistically consistent method of species tree
reconstruction. This brief note compares the accuracy of this approach to other
popular statistically consistent species tree reconstruction algorithms
including ASTRAL-II Neighbor Joining using average gene-tree internode
distances (NJst) and SVD-Quartets+PAUP*, as well as concatenation using maximum
likelihood (RaxML). We find that the faster Neighbor Joining, applied to
concatenated sequences, is among the most effective of these methods for
accurate species tree reconstruction.


  A major problem for evolutionary theory is understanding the so called {\em
open-ended} nature of evolutionary change, from its definition to its origins.
Open-ended evolution (OEE) refers to the unbounded increase in complexity that
seems to characterise evolution on multiple scales. This property seems to be a
characteristic feature of biological and technological evolution and is
strongly tied to the generative potential associated with combinatorics, which
allows the system to grow and expand their available state spaces.
Interestingly, many complex systems presumably displaying OEE, from language to
proteins, share a common statistical property: the presence of Zipf's law.
Given an inventory of basic items (such as words or protein domains) required
to build more complex structures (sentences or proteins) Zipf's law tells us
that most of these elements are rare whereas a few of them are extremely
common. Using Algorithmic Information Theory, in this paper we provide a
fundamental definition for open-endedness, which can be understood as {\em
postulates}. Its statistical counterpart, based on standard Shannon Information
theory, has the structure of a variational problem which is shown to lead to
Zipf's law as the expected consequence of an evolutionary process displaying
OEE. We further explore the problem of information conservation through an OEE
process and we conclude that statistical information (standard Shannon
information) is not conserved, resulting into the paradoxical situation in
which the increase of information content has the effect of erasing itself. We
prove that this paradox is solved if we consider non-statistical forms of
information. This last result implies that standard information theory may not
be a suitable theoretical framework to explore the persistence and increase of
the information content in OEE systems.


  Mutualistic networks have been shown to involve complex patterns of
interactions among animal and plant species. The architecture of these webs
seems to pervade some of their robust and fragile behaviour. Recent work
indicates that there is a strong correlation between the patterning of
animal-plant interactions and their phylogenetic organisation. Here we show
that such pattern and other reported regularities from mutualistic webs can be
properly explained by means of a very simple model of speciation and
divergence. This model also predicts a co-extinction dynamics under species
loss consistent with the presence of an evolutionary signal. The agreement
between observed and model networks suggests that some patterns displayed by
real mutualistic webs might actually represent evolutionary spandrels.


  We propose a new approach for modelling the process of RNA folding as a graph
transformation guided by the global value of free energy. Since the folding
process evolves towards a configuration in which the free energy is minimal,
the global behaviour resembles the one of a self-adaptive system. Each RNA
configuration is a graph and the evolution of configurations is constrained by
precise rules that can be described by a graph grammar.


  Geometric, topological and graph theory modeling and analysis of biomolecules
are of essential importance in the conceptualization of molecular structure,
function, dynamics, and transport. On the one hand, geometric modeling provides
molecular surface and structural representation, and offers the basis for
molecular visualization, which is crucial for the understanding of molecular
structure and interactions. On the other hand, it bridges the gap between
molecular structural data and theoretical/mathematical models. Topological
analysis and modeling give rise to atomic critical points and connectivity, and
shed light on the intrinsic topological invariants such as independent
components (atoms), rings (pockets) and cavities. Graph theory analyzes
biomolecular interactions and reveals biomolecular structure-function
relationship. In this paper, we review certain geometric, topological and graph
theory apparatuses for biomolecular data modeling and analysis. These
apparatuses are categorized into discrete and continuous ones. For discrete
approaches, graph theory, Gaussian network model, anisotropic network model,
normal mode analysis, quasi-harmonic analysis, flexibility and rigidity index,
molecular nonlinear dynamics, spectral graph theory, and persistent homology
are discussed. For continuous mathematical tools, we present discrete to
continuum mapping, high dimensional persistent homology, biomolecular geometric
modeling, differential geometry theory of surfaces, curvature evaluation,
variational derivation of minimal molecular surfaces, atoms in molecule theory
and quantum chemical topology. Four new approaches, including analytical
minimal molecular surface, Hessian matrix eigenvalue map, curvature map and
virtual particle model, are introduced for the first time to bridge the gaps in
biomolecular modeling and analysis.


  Migration is a key mechanism for expansion of communities. In spatially
heterogeneous environments, rapidly gaining knowledge about the local
environment is key to the evolutionary success of a migrating population. For
historical human migration, environmental heterogeneity was naturally
asymmetric in the north-south (NS) and east-west (EW) directions. We here
consider the human migration process in the Americas, modeled as random,
asymmetric, modularly correlated environments. Knowledge about the environments
determines the fitness of each individual. We present a phase diagram for
asymmetry of migration as a function of carrying capacity and fitness
threshold. We find that the speed of migration is proportional to the inverse
complement of the spatial environmental gradient, and in particular we find
that north-south migration rates are lower than east-west migration rates when
the environmental gradient is higher in the north-south direction.
Communication of knowledge between individuals can help to spread beneficial
knowledge within the population. The speed of migration increases when
communication transmits pieces of knowledge that contribute in a modular way to
the fitness of individuals. The results for the dependence of migration rate on
asymmetry and modularity are consistent with existing archaeological
observations. The results for asymmetry of genetic divergence are consistent
with patterns of human gene flow.


  While confirming the long held view that viruses do not closely imitate the
use of their host's codon catalogue, Esposito and coworkers nevertheless
consider it surprising that, despite having the ability to infect the same
host, many mycobacteriophages share little or no genetic similarity (i.e.
similarity in their GC contents and codon utilization patterns). Arguing
correctly that efficient translation of a phage's proteins within a host is
likely to be optimized by the phage's ability to match the host's codon usage
pattern, it is concluded that the preferred host of many mycobacteriophages is
not Mycobacterium smegmatis, despite their having been isolated on that
organism. Thus, a virus and its elusive preferred hosts would have had similar
GC percentages and codon usages, but the same virus could still infect a
less-preferred host (Mycobacterium smegmatis), where the virus-host similarity
would be less evident. However, there is another evolutionary interpretation.


  Recently the physical characterization of a number of biological processes
has proven indispensable for a full understanding of natural phenomena. One
such example is the mechanical properties of transcription, which have been
shown to have significant effects in gene expression. In this letter we
introduce a simple description of the basic physical elements of transcription
where RNA elongation, RNA polymerase rotation and DNA supercoiling are coupled.
The resulting framework describes the relative amount of RNA polymerase
rotation and DNA supercoiling that occurs during RNA elongation. Asymptotic
behavior is derived and can be used to experimentally extract unknown
mechanical parameters of transcription. Incorporation of mechanical limits to
RNA polymerase is accomplished yielding an equation of motion for DNA
supercoiling and RNA elongation with transcriptional stalling. Important
implications for gene expression, chromatin structure and genome organization
are discussed.


  Signal processing in biological systems is delicately executed by specialised
networks, which are modular assemblies of network motifs. The motifs are
independently functional circuits found in enormous numbers in any living cell.
A very common network motif is the feed-forward loop (FFL), which regulates a
downstream node by an upstream one in a direct and an indirect way within the
network. If the direct and indirect regulations go antagonistic, the motif is
known as an incoherent FFL (ICFFL). The current study is aimed at exploring the
reason for the variation in the evolutionary selection of the four types of
ICFFLs. As comparative measures, I compute sensitivity amplification,
adaptation precision and efficiency from the temporal dynamics and mutual
information between the input-output nodes of the motifs at steady state. The
ICFFL II performs very efficiently in adaptation but poor in information
processing. On the other hand, ICFFL I and III are better in information
transmission compared to adaptation efficiency. Which is the fittest among them
under the pressure of natural selection? To sort out this puzzle, I take help
from the multi-objective Pareto efficiency. The results, found in the Pareto
task space, are in good agreement with the reported abundance level of all the
types in eukaryotes as well as prokaryotes.


  The principal pacemaker of the circadian clock of the cyanobacterium S.
elongatus is a protein phosphorylation cycle consisting of three proteins,
KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of
two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the
CII domain can be phosphorylated, at two residues, in a well-defined sequence.
While this system has been studied extensively, how the clock is driven
thermodynamically has remained elusive. Inspired by recent experimental
observations and building on ideas from previous mathematical models, we
present a new, thermodynamically consistent, statistical-mechanical model of
the clock. At its heart are two main ideas: i) ATP hydrolysis in the CI domain
provides the thermodynamic driving force for the clock, switching KaiC between
an active conformational state in which its phosphorylation level tends to rise
and an inactive one in which it tends to fall; ii) phosphorylation of the CII
domain provides the timer for the hydrolysis in the CI domain. The model also
naturally explains how KaiA, by acting as a nucleotide exchange factor, can
stimulate phosphorylation of KaiC, and how the differential affinity of KaiA
for the different KaiC phosphoforms generates the characteristic temporal order
of KaiC phosphorylation. As the phosphorylation level in the CII domain rises,
the release of ADP from CI slows down, making the inactive conformational state
of KaiC more stable. In the inactive state, KaiC binds KaiB, which not only
stabilizes this state further, but also leads to the sequestration of KaiA, and
hence to KaiC dephosphorylation. Using a dedicated kinetic Monte Carlo
algorithm, we show that the model can describe a wealth of experimental data.


  We report an experimental study of the influences of the fixed charge and
bulk ionic concentrations on the conduction of biological ion channels, and we
consider the results within the framework of the ionic Coulomb blockade model
of permeation and selectivity. Voltage clamp recordings were used to
investigate the Na$^+$/Ca$^{2+}$ anomalous mole fraction effect (AMFE)
exhibited by the bacterial sodium channel NaChBac and its mutants.
Site-directed mutagenesis was used to study the effect of either increasing or
decreasing the fixed charge in their selectivity filters for comparison with
the predictions of the Coulomb blockade model. The model was found to describe
well some aspects of the experimental (divalent blockade and AMFE) and
simulated (discrete multi-ion conduction and occupancy band) phenomena,
including a concentration-dependent shift of the Coulomb staircase. These
results substantially extend the understanding of ion channel selectivity and
may also be applicable to biomimetic nanopores with charged walls.


  Computational approaches to drug discovery can reduce the time and cost
associated with experimental assays and enable the screening of novel
chemotypes. Structure-based drug design methods rely on scoring functions to
rank and predict binding affinities and poses. The ever-expanding amount of
protein-ligand binding and structural data enables the use of deep machine
learning techniques for protein-ligand scoring.
  We describe convolutional neural network (CNN) scoring functions that take as
input a comprehensive 3D representation of a protein-ligand interaction. A CNN
scoring function automatically learns the key features of protein-ligand
interactions that correlate with binding. We train and optimize our CNN scoring
functions to discriminate between correct and incorrect binding poses and known
binders and non-binders. We find that our CNN scoring function outperforms the
AutoDock Vina scoring function when ranking poses both for pose prediction and
virtual screening.


  The dominant paradigm in origin of life research is that of an RNA world.
However, despite experimental progress towards the spontaneous formation of
RNA, the RNA world hypothesis still has its problems. Here, we introduce a
novel computational model of chemical reaction networks based on RNA secondary
structure and analyze the existence of autocatalytic sub-networks in random
instances of this model, by combining two well-established computational tools.
Our main results are that (i) autocatalytic sets are highly likely to exist,
even for very small reaction networks and short RNA sequences, and (ii)
sequence diversity seems to be a more important factor in the formation of
autocatalytic sets than sequence length. These findings could shed new light on
the probability of the spontaneous emergence of an RNA world as a network of
mutually collaborative ribozymes.


  Microbes produce metabolic resources that are important for cell growth yet
leak across membranes into the extracellular environment. Other microbes in the
same environment can use these resources and adjust their own metabolic
production accordingly---causing other resources to leak into the environment.
The combined effect of these processes is an economy in which organismal growth
and metabolic production are coupled to others in the community. We propose a
model for the co-evolving dynamics of metabolite concentrations, production
regulation, and population frequencies for the case of two cell types, each
requiring and capable of producing two metabolites. In this model, beneficial
trade relations emerge without any coordination, via individual-level
production decisions that maximize each cell's growth rate given its perceived
environment. As we vary production parameters of the model, we encounter three
paradoxical behaviors, where a change that should intuitively benefit some cell
type, actually harms it. (1) If a cell type is more efficient than its
counterpart at producing a metabolite and becomes even more efficient, its
frequency in the population can decrease. (2) If a cell type is less efficient
than its counterpart at producing a metabolite but becomes less inefficient,
the growth rate of the population can decrease. (3) Finally, if a cell type
controls its counterpart's production decisions so as to maximize its own
growth rate, the ultimate growth rate it achieves can be lower than if the two
cell types each maximized their own growth. These three paradoxes highlight the
complex and counter-intuitive dynamics that emerge in simple microbial
economies.


  In a previous paper, we examined a class of possible conformations for
helically patterned filaments in contact with a bonding surface. In particular,
we investigated geometries where contact between the pattern and the surface
was improved through a periodic twisting and lifting of the filament. A
consequence of this lifting is that the total length of the filament projected
onto the surface decreases after bonding. When the bonding character of the
surface is actuated, this phenomenon can lead to both lifelike "inchworm"
behavior of the filaments and ensemble movement. We illustrate, through
simulation, how pattern formation may be achieved through this mechanism.


  Transcriptional profiling on microarrays to obtain gene expressions has been
used to facilitate cancer diagnosis. We propose a deep generative machine
learning architecture (called DeepCancer) that learn features from unlabeled
microarray data. These models have been used in conjunction with conventional
classifiers that perform classification of the tissue samples as either being
cancerous or non-cancerous. The proposed model has been tested on two different
clinical datasets. The evaluation demonstrates that DeepCancer model achieves a
very high precision score, while significantly controlling the false positive
and false negative scores.


  Despite a low global prevalence, infections with Sarcoptes scabiei, or
scabies, are still common in remote communities such as in northern Australia
and the Solomon Islands. Mass drug administration (MDA) has been utilised in
these communities, and although prevalence drops substantially initially, these
reductions have not been sustained. We develop a compartmental model of scabies
infection dynamics and incorporate both ovicidal and non-ovicidal treatment
regimes. By including the dynamics of mass drug administration, we are able to
reproduce the phenomena of an initial reduction in prevalence, followed by the
recrudescence of infection levels in the population. We show that even under a
`perfect' two-round MDA, eradication of scabies under a non-ovicidal treatment
scheme is almost impossible. We then go on to consider how the probability of
elimination varies with the number of treatment rounds delivered in an MDA. We
find that even with infeasibly large numbers of treatment rounds, elimination
remains challenging.


  Zika virus (ZIKV), a disease spread primarily through the Aedes aegypti
mosquito, was identified in Brazil in 2015 and was declared a global health
emergency by the World Health Organization (WHO). Epidemiologists often use
common state-level attributes such as population density and temperature to
determine the spread of disease. By applying techniques from topological data
analysis, we believe that epidemiologists will be able to better predict how
ZIKV will spread. We use the Vietoris-Rips filtration on high-density mosquito
locations in Brazil to create simplicial complexes, from which we extract
homology group generators. Previously epidemiologists have not relied on
topological data analysis to model disease spread. Evaluating our model on ZIKV
case data in the states of Brazil demonstrates the value of these techniques
for the improved assessment of vector-borne diseases.


  Current models of animal evolution focus on selection of individuals,
ignoring the much faster selection of symbiotic bacteria. Here we take
host-symbiont interactions into account by introducing a Population
Genetics-like model of holobionts exposed to toxic stress. The stress can be
alleviated by selection of resistant individuals (host and bacteria) and by
secretion of a detoxification agent ("detox"). By defining a new measure,
termed the Lamarckian, we show that selection of resistant bacteria over one
generation of hosts leads to stress-dependent increase in the tolerance of the
host's offspring. This benefit is mediated by co-alleviation of toxic and
physiologic stress. Prolonged exposure leads to further adaptation by 'group
selection' of bacterial communities with higher detox per bacterium. These
findings show that Lamarckian adaptation can arise via interactions between two
levels of Darwinian selection within a holobiont system. The conclusions and
modelling framework are applicable to diverse types of holobiont systems.


  We model a social-encounter network where linked nodes match for reproduction
in a manner depending probabilistically on each node`s attractiveness. The
developed model reveals that increasing either the network`s mean degree or the
``choosiness`` exercised during pair-formation increases the strength of
positive assortative mating. That is, we note that attractiveness is correlated
among mated nodes. Their total number also increases with mean degree and
selectivity during pair-formation. By iterating over model mapping of parents
onto offspring across generations, we study the evolution of attractiveness.
Selection mediated by exclusion from reproduction increases mean
attractiveness, but is rapidly balanced by skew in the offspring distribution
of highly attractive mated pairs.


  We introduce a model, based on the Evolutionary Game Theory, for studying the
dynamics of group formation. The latter constitutes a relevant phenomenon
observed in different animal species, whose individuals tend to cluster
together forming groups of different size. Results of previous investigations
suggest that this phenomenon might have similar reasons across different
species, such as improving the individual safety (e.g. from predators), and
increasing the probability to get food resources. Remarkably, the group size
might strongly vary from species to species, and sometimes even within the same
species. In the proposed model, an agent population tries to form homogeneous
groups. The homogeneity of a group is computed according to a spin vector, that
characterizes each agent, and represents a set a features (e.g. physical
traits). We analyze the formation of groups of different size, on varying a
parameter named 'individual payoff'. The latter represents the gain one agent
would receive acting individually. In particular, the agents choose whether to
form a group (receiving a 'group payoff'), or if to play individually
(receiving an 'individual payoff'). The phase diagram representing the
equilibria of our population shows a sharp transition between the 'group phase'
and the 'individual phase', in correspondence of a critical 'individual
payoff'. In addition, we found that forming (homogeneous) small groups is
easier than forming big groups. To conclude, we deem that our model and the
related results supports the hypothesis that the phenomenon of group formation
has evolutionary roots.


  The study of populations subject to the phenomenon of loss and fragmentation
of habitat, transforming continuous areas into small ones, usually surrounded
by anthropogenic matrices, has been the focus of many researches within the
scope of conservation. The objective of this study was to develop a computer
model by introducing modifications to the renowned Penna Model for biological
aging, in order to evaluate the behavior of populations subjected to the
effects of fragmented environments. As an object of study, it was used
biological data of the robust tufted capuchin (Sapajus robustus), an endangered
primate species whose geographical distribution within the Atlantic Rain Forest
is part of the backdrop of intense habitat fragmentation. The simulations
showed the expected behavior based on the three main aspects that affects
populations under intense habitat fragmentation: the population density, area
and conformation of the fragments and deleterious effects due the low genetic
variability in small and isolated populations. The model showed itself suitable
to describe changes in viability and population dynamics of the species crested
capuchin considering critical levels of survival in a fragmented environment
and also, actions in order to preserve the species should be focused not only
on increasing available area but also in dispersion dynamics.


  We propose a numerical approach to study the invasion fitness of a mutant and
to determine evolutionary singular strategies in evolutionary structured models
in which the competitive exclusion principle holds. Our approach is based on a
dual representation, which consists of the modelling of the small size mutant
population by a stochastic model and the computation of its corresponding
deterministic model. The use of the deterministic model greatly facilitates the
numerical determination of the feasibility of invasion as well as the
convergence-stability of the evolutionary singular strategy. Our approach
combines standard adaptive dynamics with the link between the mutant survival
criterion in the stochastic model and the sign of the eigenvalue in the
corresponding deterministic model. We present our method in the context of a
mass-structured individual-based chemostat model. We exploit a previously
derived mathematical relationship between stochastic and deterministic
representations of the mutant population in the chemostat model to derive a
general numerical method for analyzing the invasion fitness in the stochastic
models. Our method can be applied to the broad class of evolutionary models for
which a link between the stochastic and deterministic invasion fitnesses can be
established.


  Understanding the spread of any disease is a highly complex and
interdisciplinary exercise as biological, social, geographic, economic, and
medical factors may shape the way a disease moves through a population and
options for its eventual control or eradication. Disease spread poses a serious
threat in animal and plant health and has implications for ecosystem
functioning and species extinctions as well as implications in society through
food security and potential disease spread in humans. Space-time epidemiology
is based on the concept that various characteristics of the pathogenic agents
and the environment interact in order to alter the probability of disease
occurrence and form temporal or spatial patterns. Epidemiology aims to identify
these patterns and factors, to assess the relevant uncertainty sources, and to
describe disease in the population. Thus disease spread at the population level
differs from the approach traditionally taken by veterinary practitioners that
are principally concerned with the health status of the individual. Patterns of
disease occurrence provide insights into which factors may be affecting the
health of the population, through investigating which individuals are affected,
where are these individuals located and when did they become infected. With the
rapid development of smart sensors, social networks, as well as digital maps
and remotely-sensed imagery spatio-temporal data are more ubiquitous and richer
than ever before. The availability of such large datasets (Big data) poses
great challenges in data analysis. In addition, increased availability of
computing power facilitates the use of computationally-intensive methods for
the analysis of such data. Thus new methods as well as case studies are needed
to understand veterinary and ecological epidemiology. A special issue aimed to
address this topic.


  This article presents an algorithm for the evaluation of organelles'
movements inside of an unmodified live cell. We used a time-lapse image series
obtained using wide-field bright-field photon transmission microscopy as an
algorithm input. The benefit of the algorithm is the application of the R\'enyi
information entropy, namely a variable called a point information gain, which
enables to highlight the borders of the intracellular organelles and to
localize the organelles' centers of mass with the precision of one pixel.


  Using a microfluidic trap, we study the behavior of individual phospholipid
vesicles in contact with fatty acids. We show that spontaneous fatty acids
insertion inside the bilayer is controlled by the vesicle size, osmotic
pressure difference across the membrane and fatty acids concentration in the
external bath. Depending on these parameters, vesicles can grow spherically or
become unstable and fragment into several daughter vesicles. We establish the
phase diagram for vesicle growth and we derive a simple thermodynamic model
that reproduces the time evolution of the vesicle volume. Finally, we show that
stable growth can be achieved on an artificial cell expressing a simple set of
bacterial cytoskeletal proteins, paving the way toward artificial cell
reproduction.


  Motivation: Biological data and knowledge bases increasingly rely on Semantic
Web technologies and the use of knowledge graphs for data integration,
retrieval and federated queries. In the past years, feature learning methods
that are applicable to graph-structured data are becoming available, but have
not yet widely been applied and evaluated on structured biological knowledge.
Results: We develop a novel method for feature learning on biological knowledge
graphs. Our method combines symbolic methods, in particular knowledge
representation using symbolic logic and automated reasoning, with neural
networks to generate embeddings of nodes that encode for related information
within knowledge graphs. Through the use of symbolic logic, these embeddings
contain both explicit and implicit information. We apply these embeddings to
the prediction of edges in the knowledge graph representing problems of
function prediction, finding candidate genes of diseases, protein-protein
interactions, or drug target relations, and demonstrate performance that
matches and sometimes outperforms traditional approaches based on manually
crafted features. Our method can be applied to any biological knowledge graph,
and will thereby open up the increasing amount of Semantic Web based knowledge
bases in biology to use in machine learning and data analytics. Availability
and Implementation:
https://github.com/bio-ontology-research-group/walking-rdf-and-owl Contact:
robert.hoehndorf@kaust.edu.sa


  The understanding of the macroscopic phenomenological models of the
population growth at a microscopic level is important to predict the population
behaviors emerged from the interactions between the individuals. In this work
we consider the influence of the cell-cell interaction on the population growth
rate $R$ in a tumor system, and show that, in most cases especially small
proliferative probabilities, the regulative role of the interaction will be
strengthened with the decline of the intrinsic proliferative probabilities. For
the high replication rates of an individual and the cooperative interactions,
the proliferative probability almost has no effect. We compute the dependences
of $R$ on the interactions between the cells under the approximation of the
nearest neighbor in the rim of an avascular tumor. Our results are helpful to
qualitatively understand the influence of the interactions between the
individuals on the growth rate in population systems.


  A two-types, discrete-time population model with finite, constant size is
constructed, allowing for a general form of frequency-dependent selection and
skewed offspring distribution. Selection is defined based on the idea that
individuals first choose a (random) number of $\textit{potential}$ parents from
the previous generation and then, from the selected pool, they inherit the type
of the fittest parent. The probability distribution function of the number of
potential parents per individual thus parametrises entirely the selection
mechanism. Using sampling- and moment-duality, weak convergence is then proved
both for the allele frequency process of the selectively weak type and for the
population's ancestral process. The scaling limits are, respectively, a
two-types $\Xi$-Fleming-Viot jump-diffusion process with frequency-dependent
selection, and a branching-coalescing process with general branching and
simultaneous multiple collisions. Duality also leads to a characterisation of
the probability of extinction of the selectively weak allele, in terms of the
ancestral process' ergodic properties.


  We generalize a model of growth over a disordered environment, to a large
class of It\=o processes. In particular, we study how the microscopic
properties of the noise influence the macroscopic growth rate. The present
model can account for growth processes in large dimensions, and provides a bed
to understand better the trade-off between exploration and exploitation. An
additional mapping to the Schr\"ordinger equation readily provides a set of
disorders for which this model can be solved exactly. This mean-field approach
exhibits interesting features, such as a freezing transition and an optimal
point of growth, that can be studied in details, and gives yet another
explanation for the occurrence of the $\textit{Zipf law}$ in complex,
well-connected systems.


  We study evolutionary multi-player games in finite populations, subject to
fluctuating environments. The population undergoes a birth-death process with
absorbing states, and the environment follows a Markovian process, resulting in
a fluctuating payoff matrix for the evolutionary game. Our focus is on the
fixation or extinction of a single mutant in a population of wildtypes. We show
that the nonlinear nature of fitnesses in multi-player games gives rise to an
intricate interplay of selection, genetic drift and environmental fluctuations.
This generates effects not seen in simpler two-player games. To analyse
trajectories towards fixation we analytically calculate sojourn times for
general birth-death processes in populations of two types of individuals and in
fluctuating environments.


  Theoretical analysis, which maps single molecule time trajectories of a
molecular motor onto unicyclic Markov processes, allows us to evaluate the heat
dissipated from the motor and to elucidate its dependence on the mean velocity
and diffusivity. Unlike passive Brownian particles in equilibrium, the velocity
and diffusion constant of molecular motors are closely inter-related to each
other. In particular, our study makes it clear that the increase of diffusivity
with the heat production is a natural outcome of active particles, which is
reminiscent of the recent experimental premise that the diffusion of an
exothermic enzyme is enhanced by the heat released from its own catalytic
turnover. Compared with freely diffusing exothermic enzymes, kinesin-1 whose
dynamics is confined on one-dimensional tracks is highly efficient in
transforming conformational fluctuations into a locally directed motion, thus
displaying a significantly higher enhancement in diffusivity with its turnover
rate. Putting molecular motors and freely diffusing enzymes on an equal
footing, our study offers thermodynamic basis to understand the heat enhanced
self-diffusion of exothermic enzymes.


  Revealing how a biological network is organized to realize its function is
one of the main topics in systems biology. The functional backbone network,
defined as the primary structure of the biological network, is of great
importance in maintaining the main function of the biological network. We
propose a new algorithm, the tinker algorithm, to determine this core structure
and apply it in the cell-cycle system. With this algorithm, the backbone
network of the cell-cycle network can be determined accurately and efficiently
in various models such as the Boolean model, stochastic model, and ordinary
differential equation model. Results show that our algorithm is more efficient
than that used in the previous research. We hope this method can be put into
practical use in relevant future studies.


  Background: Guanine quadruplexes (GQs) play vital roles in many cellular
processes and are of much interest as drug targets. In contrast to the
availability of many structural studies, there is still limited knowledge on GQ
folding.
  Scope of review: We review recent molecular dynamics (MD) simulation studies
of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We
explain the basic principles and limitations of all types of MD methods used to
study unfolding and folding in a way accessible to non-specialists. We discuss
the potential role of G-hairpin, G-triplex and alternative GQ intermediates in
the folding process. We argue that, in general, folding of GQs is fundamentally
different from funneled folding of small fast-folding proteins, and can be best
described by a kinetic partitioning (KP) mechanism. KP is a competition between
at least two (but often many) well-separated and structurally different
conformational ensembles.
  Major conclusions: The KP mechanism is the only plausible way to explain
experiments reporting long time-scales of GQ folding and the existence of
long-lived sub-states. A significant part of the natural partitioning of the
free energy landscape of GQs comes from the ability of the GQ-forming sequences
to populate a large number of syn-anti patterns in their G-tracts. The extreme
complexity of the KP of GQs typically prevents an appropriate description of
the folding landscape using just a few order parameters or collective
variables.
  General significance: We reconcile available computational and experimental
studies of GQ folding and formulate basic principles characterizing GQ folding
landscapes


  It has been hypothesized that one of the main reasons evolution has been able
to produce such impressive adaptations is because it has improved its own
ability to evolve -- "the evolution of evolvability". Rupert Riedl, for
example, an early pioneer of evolutionary developmental biology, suggested that
the evolution of complex adaptations is facilitated by a developmental
organization that is itself shaped by past selection to facilitate evolutionary
innovation. However, selection for characteristics that enable future
innovation seems paradoxical: natural selection cannot favor structures for
benefits they have not yet produced, and favoring characteristics for benefits
that have already been produced does not constitute future innovation. Here we
resolve this paradox by exploiting a formal equivalence between the evolution
of evolvability and learning systems. We use the conditions that enable simple
learning systems to generalize, i.e., to use past experience to improve
performance on previously unseen, future test cases, to demonstrate conditions
where natural selection can systematically favor developmental organizations
that benefit future evolvability. Using numerical simulations of evolution on
highly epistatic fitness landscapes, we illustrate how the structure of evolved
gene regulation networks can result in increased evolvability capable of
avoiding local fitness peaks and discovering higher fitness phenotypes. Our
findings support Riedl's intuition: Developmental organizations that "mimic"
the organization of constraints on phenotypes can be favored by short-term
selection and also facilitate future innovation. Importantly, the conditions
that enable the evolution of such surprising evolvability follow from the same
non-mysterious conditions that permit generalization in learning systems.


  In the context of bacteria and models of their evolution under genome
rearrangement, we explore a novel application of group representation theory to
the inference of evolutionary history. Our contribution is to show, in a very
general maximum likelihood setting, how to use elementary matrix algebra to
sidestep intractable combinatorial computations and convert the problem into
one of eigenvalue estimation amenable to standard numerical approximation
techniques.


  Plasmodium vivax populations are more resistant to malaria control strategies
than Plasmodium falciparum, maintaining high genetic diversity and gene flow
even at low transmission. To quantify the impact of declining transmission on
P. vivax populations, we investigated population genetic structure over time
during intensified control efforts and over a wide range of transmission
intensities and spatial scales in the Southwest Pacific. Analysis of 887 P.
vivax microsatellite haplotypes (Papua New Guinea, PNG = 443, Solomon Islands =
420, Vanuatu =24) revealed substantial population structure among countries and
modestly declining diversity as transmission decreases over space and time. In
the Solomon Islands, which has had sustained control efforts for 20 years,
significant population structure was observed on different spatial scales down
to the sub-village level. Up to 37% of alleles were partitioned between
populations and significant multilocus linkage disequilibrium was observed
indicating substantial inbreeding. High levels of haplotype relatedness around
households and within a range of 300m are consistent with a focal and clustered
infections suggesting that restricted local transmission occurs within the
range of vector movement and that subsequent focal inbreeding may be a key
factor contributing to the observed population structure. We conclude that
unique transmission strategies, including relapse allows P. vivax populations
to withstand pressure from control efforts for longer than P. falciparum.
However sustained control efforts do eventually impact parasite population
structure and with further control pressure, populations may eventually
fragment into clustered foci that could be targeted for elimination.


  We investigate the influence of intrinsic noise on stable states of a
one-dimensional dynamical system that shows in its deterministic version a
saddle-node bifurcation between monostable and bistable behaviour. The system
is a modified version of the Schl\"ogl model, which is a chemical reaction
system with only one type of molecule. The strength of the intrinsic noise is
varied without changing the deterministic description by introducing bursts in
the autocatalytic production step. We study the transitions between monostable
and bistable behavior in this system by evaluating the number of maxima of the
stationary probability distribution. We find that changing the size of bursts
can destroy and even induce saddle-node bifurcations. This means that a bursty
production of molecules can qualitatively change the dynamics of a chemical
reaction system even when the deterministic description remains unchanged.


  N-methyl-D-aspartate receptors (NMDARs) are glycoproteins in the brain
central to learning and memory. The effects of glycosylation on the structure
and dynamics of NMDARs are largely unknown. In this work, we use extensive
molecular dynamics simulations of GluN1 and GluN2B ligand binding domains
(LBDs) of NMDARs to investigate these effects. Our simulations predict that
intra-domain interactions involving the glycan attached to residue GluN1-N440
stabilize closed-clamshell conformations of the GluN1 LBD. The glycan on
GluN2B-N688 shows a similar, though weaker, effect. Based on these results, and
assuming the transferability of the results of LBD simulations to the full
receptor, we predict that glycans at GluN1-N440 might play a potentiator role
in NMDARs. To validate this prediction, we perform electrophysiological
analysis of full-length NMDARs with a glycosylation-preventing GluN1-N440Q
mutation, and demonstrate an increase in the glycine EC50 value. Overall, our
results suggest an intramolecular potentiating role of glycans on NMDA
receptors.


  Species experience both internal feedbacks with endogenous factors such as
trait evolution and external feedbacks with exogenous factors such as weather.
These feedbacks can play an important role in determining whether populations
persist or communities of species coexist. To provide a general mathematical
framework for studying these effects, we develop a theorem for coexistence for
ecological models accounting for internal and external feedbacks. Specifically,
we use average Lyapunov functions and Morse decompositions to develop
sufficient and necessary conditions for robust permanence, a form of
coexistence robust to large perturbations of the population densities and small
structural perturbations of the models. We illustrate how our results can be
applied to verify permanence in non-autonomous models, structured population
models, including those with frequency-dependent feedbacks, and models of
eco-evolutionary dynamics. In these applications, we discuss how our results
relate to previous results for models with particular types of feedbacks.


  Reconstructing the causal network in a complex dynamical system plays a
crucial role in many applications, from sub-cellular biology to economic
systems. Here we focus on inferring gene regulation networks (GRNs) from
perturbation or gene deletion experiments. Despite their scientific merit, such
perturbation experiments are not often used for such inference due to their
costly experimental procedure, requiring significant resources to complete the
measurement of every single experiment. To overcome this challenge, we develop
the Robust IDentification of Sparse networks (RIDS) method that reconstructs
the GRN from a small number of perturbation experiments. Our method uses the
gene expression data observed in each experiment and translates that into a
steady state condition of the system's nonlinear interaction dynamics. Applying
a sparse optimization criterion, we are able to extract the parameters of the
underlying weighted network, even from very few experiments. In fact, we
demonstrate analytically that, under certain conditions, the GRN can be
perfectly reconstructed using $K = \Omega (d_{max})$ perturbation experiments,
where $d_{max}$ is the maximum in-degree of the GRN, a small value for
realistic sparse networks, indicating that RIDS can achieve high performance
with a scalable number of experiments. We test our method on both synthetic and
experimental data extracted from the DREAM5 network inference challenge. We
show that the RIDS achieves superior performance compared to the
state-of-the-art methods, while requiring as few as ~60% less experimental
data. Moreover, as opposed to almost all competing methods, RIDS allows us to
infer the directionality of the GRN links, allowing us to infer empirical GRNs,
without relying on the commonly provided list of transcription factors.