  Intracellular bidirectional transport of cargo on microtubule filaments is
achieved by the collective action of oppositely directed dynein and kinesin
motors. Experiments have found that in certain cases, inhibiting the activity
of one type of motor results in an overall decline in the motility of the
cellular cargo in both directions. This counter-intuitive observation, referred
to as {\em paradox of codependence} is inconsistent with the existing paradigm
of a mechanistic tug-of-war between oppositely directed motors. Unlike kinesin,
dynein motors exhibit catchbonding, wherein the unbinding rates of these motors
decrease with increasing force on them. Incorporating this catchbonding
behavior of dynein in a theoretical model, we show that the functional
divergence of the two motors species manifests itself as an internal regulatory
mechanism, and leads to codependent transport behaviour in biologically
relevant regimes. Using analytical methods and stochastic simulations, we
analyse the processivity characteristics and probability distribution of run
times and pause times of transported cellular cargoes. We show that
catchbonding can drastically alter the transport characteristics and also
provide a plausible resolution of the paradox of codependence.


  The notions of taxis and kinesis are introduced and used to describe two
types of behavior of an organism in non-uniform conditions: (i) Taxis means the
guided movement to more favorable conditions; (ii) Kinesis is the
non-directional change in space motion in response to the change of conditions.
Migration and dispersal of animals has evolved under control of natural
selection. In a simple formalisation, the strategy of dispersal should increase
Darwinian fitness. We introduce new models of purposeful kinesis with diffusion
coefficient dependent on fitness. The local and instant evaluation of Darwinian
fitness is used, the reproduction coefficient. New models include one
additional parameter, intensity of kinesis, and may be considered as the {\em
minimal models of purposeful kinesis}. The properties of models are explored by
a series of numerical experiments. It is demonstrated how kinesis could be
beneficial for assimilation of patches of food or of periodic fluctuations.
Kinesis based on local and instant estimations of fitness is not always
beneficial: for species with the Allee effect it can delay invasion and
spreading. It is proven that kinesis cannot modify stability of positive
homogeneous steady states.


  We propose a Caputo type fractional-order mathematical model for the
transmission dynamics of tuberculosis (TB). Uniform asymptotic stability of the
unique endemic equilibrium of the fractional-order TB model is proved, for any
$\alpha \in (0, 1)$. Numerical simulations for the stability of the endemic
equilibrium are provided.


  Dispersal is ubiquitous throughout the tree of life: factors selecting for
dispersal include kin competition, inbreeding avoidance and spatiotemporal
variation in resources or habitat suitability. These factors differ in whether
they promote male and female dispersal equally strongly, and often selection on
dispersal of one sex depends on how much the other disperses. For example, for
inbreeding avoidance it can be sufficient that one sex disperses away from the
natal site. Attempts to understand sex-specific dispersal evolution have
created a rich body of theoretical literature, which we review here. We
highlight an interesting gap between empirical and theoretical literature. The
former associates different patterns of sex-biased dispersal with mating
systems, such as female-biased dispersal in monogamous birds and male-biased
dispersal in polygynous mammals. The predominant explanation is traceable back
to Greenwood's (1980) ideas of how successful philopatric or dispersing
individuals are at gaining mates or resources required to attract them. Theory,
however, has developed surprisingly independently of these ideas: predominant
ideas in theoretical work track how immigration and emigration change
relatedness patterns and alleviate competition for limiting resources,
typically considered sexually distinct, with breeding sites and fertilisable
females limiting reproductive success for females and males, respectively. We
show that the link between mating system and sex-biased dispersal is far from
resolved: there are studies showing that mating systems matter, but the
oft-stated association between polygyny and male-biased dispersal is not a
straightforward theoretical expectation... (full abstract in the PDF)


  We study the effect of contingent movement on the persistence of cooperation
on complex networks with empty nodes. Each agent plays Prisoner's Dilemma game
with its neighbors and then it either updates the strategy depending on the
payoff difference with neighbors or it moves to another empty node if not
satisfied with its own payoff. If no neighboring node is empty, each agent
stays at the same site. By extensive evolutionary simulations, we show that the
medium density of agents enhances cooperation where the network flow of mobile
agents is also medium. Moreover, if the movements of agents are more frequent
than the strategy updating, cooperation is further promoted. In scale-free
networks, the optimal density for cooperation is lower than other networks
because agents get stuck at hubs. Our study suggests that keeping a smooth
network flow is significant for the persistence of cooperation in ever-changing
societies.


  Computational protein design has a wide variety of applications. Despite its
remarkable success, designing a protein for a given structure and function is
still a challenging task. On the other hand, the number of solved protein
structures is rapidly increasing while the number of unique protein folds has
reached a steady number, suggesting more structural information is being
accumulated on each fold. Deep learning neural network is a powerful method to
learn such big data set and has shown superior performance in many machine
learning fields. In this study, we applied the deep learning neural network
approach to computational protein design for predicting the probability of 20
natural amino acids on each residue in a protein. A large set of protein
structures was collected and a multi-layer neural network was constructed. A
number of structural properties were extracted as input features and the best
network achieved an accuracy of 38.3%. Using the network output as residue type
restraints was able to improve the average sequence identity in designing three
natural proteins using Rosetta. Moreover, the predictions from our network show
~3% higher sequence identity than a previous method. Results from this study
may benefit further development of computational protein design methods.


  Previously, we developed a population model incorporating the Allee effect
and periodic environmental fluctuations, in which organisms alternate between
nomadic and colonial behaviours. This switching strategy is regulated by
biological clocks and the abundance of environmental resources, and can lead to
population persistence despite both behaviours being individually losing. In
the present study, we consider stochastic noise models in place of the original
periodic ones, thereby allowing a wider range of environmental fluctuations to
be modelled. The theoretical framework is generalized to account for resource
depletion by both nomadic and colonial sub-populations, and an ecologically
realistic population size-dependent switching scheme is proposed. We
demonstrate the robustness of the modified switching scheme to stochastic
noise, and we also present the intriguing possibility of consecutive
subsidence-recovery cycles within the resulting population dynamics. Our
results have relevance in biological and physical systems.


  In condensed matter physics, simplified descriptions are obtained by
coarse-graining the features of a system at a certain characteristic length,
defined as the typical length beyond which some properties are no longer
correlated. From a physics standpoint, in vitro DNA has thus a characteristic
length of 300 base pairs (bp), the Kuhn length of the molecule beyond which
correlations in its orientations are typically lost. From a biology standpoint,
in vivo DNA has a characteristic length of 1000 bp, the typical length of
genes. Since bacteria live in very different physico-chemical conditions and
since their genomes lack translational invariance, whether larger, universal
characteristic lengths exist is a non-trivial question. Here, we examine this
problem by leveraging the large number of fully sequenced genomes available in
public databases. By analyzing GC content correlations and the evolutionary
conservation of gene contexts (synteny) in hundreds of bacterial chromosomes,
we conclude that a fundamental characteristic length around 10-20 kb can be
defined. This characteristic length reflects elementary structures involved in
the coordination of gene expression, which are present all along the genome of
nearly all bacteria. Technically, reaching this conclusion required us to
implement methods that are insensitive to the presence of large idiosyncratic
genomic features, which may co-exist along these fundamental universal
structures.


  DNA and protein microarrays are a high-throughput technology that allow the
simultaneous quantification of tens of thousands of different biomolecular
species. The mediocre sensitivity and dynamic range of traditional fluorescence
microarrays compared to other techniques have been the technology's Achilles'
Heel, and prevented their adoption for many biomedical and clinical diagnostic
applications. Previous work to enhance the sensitivity of microarray readout to
the single-molecule ('digital') regime have either required signal amplifying
chemistry or sacrificed throughput, nixing the platform's primary advantages.
Here, we report the development of a digital microarray which extends both the
sensitivity and dynamic range of microarrays by about three orders of
magnitude. This technique uses functionalized gold nanorods as single-molecule
labels and an interferometric scanner which can rapidly enumerate individual
nanorods by imaging them with a 10x objective lens. This approach does not
require any chemical enhancement such as silver deposition, and scans arrays
with a throughput similar to commercial fluorescence devices. By combining
single-nanoparticle enumeration and ensemble measurements of spots when the
particles are very dense, this system achieves a dynamic range of about one
million directly from a single scan.


  Viral zoonoses have emerged as the key drivers of recent pandemics. Human
infection by zoonotic viruses are either spillover events -- isolated
infections that fail to cause a widespread contagion -- or species jumps, where
successful adaptation to the new host leads to a pandemic. Despite expensive
bio-surveillance efforts, historically emergence response has been reactive,
and post-hoc. Here we use machine inference to demonstrate a high accuracy
predictive bio-surveillance capability, designed to pro-actively localize an
impending species jump via automated interrogation of massive sequence
databases of viral proteins. Our results suggest that a jump might not purely
be the result of an isolated unfortunate cross-infection localized in space and
time; there are subtle yet detectable patterns of genotypic changes
accumulating in the global viral population leading up to emergence. Using tens
of thousands of protein sequences simultaneously, we train models that track
maximum achievable accuracy for disambiguating host tropism from the primary
structure of surface proteins, and show that the inverse classification
accuracy is a quantitative indicator of jump risk. We validate our claim in the
context of the 2009 swine flu outbreak, and the 2004 emergence of H5N1
subspecies of Influenza A from avian reservoirs; illustrating that
interrogation of the global viral population can unambiguously track a near
monotonic risk elevation over several preceding years leading to eventual
emergence.


  The Letter presents a novel way to connect random walks, stochastic
differential equations, and evolutionary game theory. We introduce a new
concept of potential function for discrete-space stochastic systems. It is
based on a correspondence between one-dimensional stochastic differential
equations and random walks, which may be exact not only in the continuous limit
but also in finite-state spaces. Our method is useful for computation of
fixation probabilities in discrete stochastic dynamical systems with two
absorbing states. We apply it to evolutionary games, formulating two simple and
intuitive criteria for evolutionary stability of pure Nash equilibria in finite
populations. In particular, we show that the $1/3$ law of evolutionary games,
introduced by Nowak et al [Nature, 2004], follows from a more general
mean-potential law.


  Diffusion processes are governed by external triggers and internal dynamics
in complex systems. Timely and cost-effective control of infectious disease
spread critically relies on uncovering the underlying diffusion mechanisms,
which is challenging due to invisible causality between events and their
time-evolving intensity. We infer causal relationships between infections and
quantify the reflexivity of a meta-population, the level of feedback on event
occurrences by its internal dynamics (likelihood of a regional outbreak
triggered by previous cases). These are enabled by our new proposed model, the
Latent Influence Point Process (LIPP) which models disease spread by
incorporating macro-level internal dynamics of meta-populations based on human
mobility. We analyse 15-year dengue cases in Queensland, Australia. From our
causal inference, outbreaks are more likely driven by statewide global
diffusion over time, leading to complex behavior of disease spread. In terms of
reflexivity, precursory growth and symmetric decline in populous regions is
attributed to slow but persistent feedback on preceding outbreaks via
inter-group dynamics, while abrupt growth but sharp decline in peripheral areas
is led by rapid but inconstant feedback via intra-group dynamics. Our proposed
model reveals probabilistic causal relationships between discrete events based
on intra- and inter-group dynamics and also covers direct and indirect
diffusion processes (contact-based and vector-borne disease transmissions).


  Biofilms are microbial collectives that occupy a diverse array of surfaces.
The function and evolution of biofilms are strongly influenced by the spatial
arrangement of different strains and species within them, but how
spatiotemporal distributions of different genotypes in biofilm populations
originate is still underexplored. Here, we study the origins of biofilm genetic
structure by combining model development, numerical simulations, and
microfluidic experiments using the human pathogen Vibrio cholerae. Using
spatial correlation functions to quantify the differences between emergent cell
lineage segregation patterns, we find that strong adhesion often, but not
always, maximizes the size of clonal cell clusters on flat surfaces.
Counterintuitively, our model predicts that, under some conditions, investing
in adhesion can reduce rather than increase clonal group size. Our results
emphasize that a complex interaction of fluid flow and cell adhesiveness can
underlie emergent patterns of biofilm genetic structure. This structure, in
turn, has an outsize influence on how biofilm-dwelling populations function and
evolve.


  We explore structural stability of weighted and unweighted networks of
positively interacting agents against a negative external field. We study how
the agents support the activity of each other to confront the negative field,
which suppresses the activity of agents and can lead to a collapse of the whole
network. The competition between the interactions and the field shape the
structure of stable states of the system. In unweighted networks (uniform
interactions) the stable states have the structure of $k$-cores of the
interaction network. The interplay between the topology and the distribution of
weights (heterogeneous interactions) impacts strongly the structural stability
against a negative field, especially in the case of fat-tailed distributions of
weights. We show that apart critical slowing down there is also a critical
change in the system structure that precedes the network collapse. The change
can serve as early warning of the critical transition. In order to characterize
changes of network structure we develop a method based on statistical analysis
of the $k$-core organization and so-called `corona' clusters belonging to the
$k$-cores.


  We present an importance sampling algorithm that can produce realisations of
Markovian epidemic models that exactly match observations, taken to be the
number of a single event type over a period of time. The importance sampling
can be used to construct an efficient particle filter that targets the states
of a system and hence estimate the likelihood to perform Bayesian parameter
inference. When used in a particle marginal Metropolis Hastings scheme, the
importance sampling provides a large speed-up in terms of the effective sample
size per unit of computational time, compared to simple bootstrap sampling. The
algorithm is general, with minimal restrictions, and we show how it can be
applied to any discrete-state continuous-time Markov chain where we wish to
exactly match the number of a single event type over a period of time.


  The technological and economic benefits of engineered nanomaterials may be
offset by their adverse effects on living organisms. One of the highly produced
nanomaterials under such scrutiny is amorphous silica nanoparticles, which are
known to have an appreciable, although reversible, inflammatory potential. This
is due to their selective toxicity toward macrophages, and it is thus important
to study the cellular responses of this cell type to silica nanoparticles to
better understand the direct or indirect adverse effects of nanosilica. We have
here studied the responses of the RAW264.7 murine macrophage cells and of the
control MPC11 plasma cells to subtoxic concentrations of nanosilica, using a
combination of pro-teomic and targeted approaches. This allowed us to document
alterations in the cellular cytoskeleton, in the phagocytic capacity of the
cells as well as their ability to respond to bacterial stimuli. More
surprisingly, silica nanoparticles also induce a greater sensitivity of
macrophages to DNA alkylating agents, such as styrene oxide, even at doses
which do not induce any appreciable cell death.


  The impact of random fluctuations on the dynamical behavior a complex
biological systems is a longstanding issue, whose understanding would shed
light on the evolutionary pressure that nature imposes on the intrinsic noise
levels and would allow rationally designing synthetic networks with controlled
noise. Using the It\=o stochastic differential equation formalism, we performed
both analytic and numerical analyses of several model systems containing
different molecular species in contact with the environment and interacting
with each other through mass-action kinetics. These systems represent for
example biomolecular oligomerization processes, complex-breakage reactions,
signaling cascades or metabolic networks. For chemical reaction networks with
zero deficiency values, which admit a detailed- or complex-balanced steady
state, all molecular species are uncorrelated. The number of molecules of each
species follow a Poisson distribution and their Fano factors, which measure the
intrinsic noise, are equal to one. Systems with deficiency one have an
unbalanced non-equilibrium steady state and a non-zero S-flux, defined as the
flux flowing between the complexes multiplied by an adequate stoichiometric
coefficient. In this case, the noise on each species is reduced if the flux
flows from the species of lowest to highest complexity, and is amplified is the
flux goes in the opposite direction. These results are generalized to systems
of deficiency two, which possess two independent non-vanishing S-fluxes, and we
conjecture that a similar relation holds for higher deficiency systems.


  Phylogenetic trees in genetics and biology in general are all binary. We make
an attempt to answer one fundamental question: Is such binary branching from
the coarsest to the finest scales sustained by data? We convert this question
into an equivalent one: where is the structural information of tree in a data
matrix? Results from this conceptual as well as computing issue afford us to
conclude a negative answer: Each branch being split into two at each inter-node
of tree from the top to bottom levels is a man-made structure. The data-driven
computing paradigm Data Mechanics is employed here to reveal that information
of tree is composed of a set of selected temperatures (or scales), each of
which has a clustering composition strictly regulated by a temperature-specific
cluster-sharing probability matrix. The resultant Data Cloud Geometry (DCG)
tree on the space of species is proposed as the authentic structure contained
in data. Particularly each core clusters on the finest scale, the bottom level,
of DCG tree should not be further partitioned because of uniformity. Beyond the
finest scale, the branching of DCG tree is primarily based on probability,
which induces an Ultrametric satisfying super triangular inequality property.
This Ultrametric property differentiates DCG tree from all popular trees based
on Hierarchical clustering (HC) algorithm, which typically employs an
empirical, often ad hoc distance measure. Since this measure is regulated by
the triangular inequality, it is not capable of producing a "flat" branch, in
which all its members (more than two) have equal distances to each others. We
demonstrate such information content on an illustrative zoo data first, and
then on two genomic data.


  Microorganisms, such as bacteria, are one of the first targets of
nanoparticles in the environment. In this study, we tested the effect of two
nanoparticles, ZnO and TiO2, with the salt ZnSO4 as the control, on the
Gram-positive bacterium Bacillus subtilis by 2D gel electrophoresis-based
proteomics. Despite a significant effect on viability (LD50), TiO2 NPs had no
detectable effect on the proteomic pattern, while ZnO NPs and ZnSO4
significantly modified B. subtilis metabolism. These results allowed us to
conclude that the effects of ZnO observed in this work were mainly attributable
to Zn dissolution in the culture media. Proteomic analysis highlighted twelve
modulated proteins related to central metabolism: MetE and MccB (cysteine
metabolism), OdhA, AspB, IolD, AnsB, PdhB and YtsJ (Krebs cycle) and XylA,
YqjI, Drm and Tal (pentose phosphate pathway). Biochemical assays, such as free
sulfhydryl, CoA-SH and malate dehydrogenase assays corroborated the observed
central metabolism reorientation and showed that Zn stress induced oxidative
stress, probably as a consequence of thiol chelation stress by Zn ions. The
other patterns affected by ZnO and ZnSO4 were the stringent response and the
general stress response. Nine proteins involved in or controlled by the
stringent response showed a modified expression profile in the presence of ZnO
NPs or ZnSO4: YwaC, SigH, YtxH, YtzB, TufA, RplJ, RpsB, PdhB and Mbl. An
increase in the ppGpp concentration confirmed the involvement of the stringent
response during a Zn stress. All these metabolic reorientations in response to
Zn stress were probably the result of complex regulatory mechanisms including
at least the stringent response via YwaC.


  What role do asymptomatically infected individuals play in the transmission
dynamics? There are many diseases, such as norovirus and influenza, where some
infected hosts show symptoms of the disease while others are asymptomatically
infected, i.e. do not show any symptoms. The current paper considers a class of
epidemic models following an SEIR (Susceptible $\to$ Exposed $\to$ Infectious
$\to$ Recovered) structure that allows for both symptomatic and asymptomatic
cases. The following question is addressed: what fraction $\rho$ of those
individuals getting infected are infected by symptomatic (asymptomatic) cases?
This is a more complicated question than the related question for the beginning
of the epidemic: what fraction of the expected number of secondary cases of a
typical newly infected individual, i.e. what fraction of the basic reproduction
number $R_0$, is caused by symptomatic individuals? The latter fraction only
depends on the type-specific reproduction numbers, while the former fraction
$\rho$ also depends on timing and hence on the probabilistic distributions of
latent and infectious periods of the two types (not only their means). Bounds
on $\rho$ are derived for the situation where these distributions (and even
their means) are unknown. Special attention is given to the class of Markov
models and the class of continuous-time Reed-Frost models as two classes of
distribution functions. We show how these two classes of models can exhibit
very different behaviour.


  A Y-linked two-sex branching process with mutations and blind choice of males
is a suitable model for analyzing the evolution of the number of carriers of an
allele and its mutations of a Y-linked gene. Considering a two-sex monogamous
population, in this model each female chooses her partner from among the male
population without caring about his type (i.e., the allele he carries). In this
work, we deal with the problem of estimating the main parameters of such model
developing the Bayesian inference in a parametric framework. Firstly, we
consider, as sample scheme, the observation of the total number of females and
males up to some generation as well as the number of males of each genotype at
last generation. Later, we introduce the information of the mutated males only
in the last generation obtaining in this way a second sample scheme. For both
samples, we apply the Approximate Bayesian Computation (ABC) methodology to
approximate the posterior distributions of the main parameters of this model.
The accuracy of the procedure based on these samples is illustrated and
discussed by way of simulated examples.


  A system made up of N interacting species is considered. Self-reaction terms
are assumed of the logistic type. Pairwise interactions take place among
species according to different modalities, thus yielding a complex asymmetric
disordered graph. A mathematical procedure is introduced and tested to
stabilise the ecosystem via an {\it ad hoc} rewiring of the underlying
couplings. The method implements minimal modifications to the spectrum of the
Jacobian matrix which sets the stability of the fixed point and traces these
changes back to species-species interactions. Resilience of the equilibrium
state appear to be favoured by predator-prey interactions.


  Development from egg to embryo to adult is a fascinating instance of
biological self-organization for which genetics has supplied us with a parts
list. It remains to find the principles organizing the assembly of those parts.
In the last decade embryonic stem cells (ESC) have provided the material from
which to build the mammalian embryo. This review, for a quantitative audience,
explains why colonies of ESC are an ideal system with which to peal back the
multiple layers of regulation that make embryonic development such a robust
process. It formed the basis of a presentation at the 27th Solvay Conference on
the Physics of Living Matter 2017 edited by Boris Shraiman.


  We analyze a class of chemical reaction networks under mass-action kinetics
and involving multiple time-scales, whose deterministic and stochastic models
display qualitative differences. The networks are inspired by gene-regulatory
networks, and consist of a slow-subnetwork, describing conversions among the
different gene states, and fast-subnetworks, describing biochemical
interactions involving the gene products. We show that the long-term dynamics
of such networks can consist of a unique attractor at the deterministic level
(unistability), while the long-term probability distribution at the stochastic
level may display multiple maxima (multimodality). The dynamical differences
stem from a novel phenomenon we call noise-induced mixing, whereby the
probability distribution of the gene products is a linear combination of the
probability distributions of the fast-subnetworks which are `mixed' by the
slow-subnetworks. The results are applied in the context of systems biology,
where noise-induced mixing is shown to play a biochemically important role,
producing phenomena such as stochastic multimodality and oscillations.


  Neurotensin receptor 1 (NTSR1) is a G protein-coupled receptor that is
important for signaling in the brain and the gut. Its agonist ligand
neurotensin (NTS), a 13-amino-acid peptide, binds with nanomolar affinity from
the extracellular side to NTSR1 and induces conformational changes that trigger
intracellular signaling processes. Our goal is to monitor the conformational
dynamics of single fluorescently labeled NTSR1. For this, we fused the
fluorescent protein mNeonGreen to the C terminus of NTSR1, purified the
receptor fusion protein from E. coli membranes, and reconstituted NTSR1 into
liposomes with E. coli polar lipids. Using single-molecule anisotropy
measurements, NTSR1 was found to be monomeric in liposomes, with a small
fraction being dimeric and oligomeric, showing homoFRET. Similar results were
obtained for NTSR1 in detergent solution. Furthermore, we demonstrated agonist
binding to NTSR1 by time-resolved single-molecule F\"orster resonance energy
transfer (smFRET), using neurotensin labeled with the fluorophore ATTO594.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.


  We introduce the exit time finite state projection (ETFSP) scheme, a
truncation-based method that yields approximations to the exit distribution and
occupation measure associated with the time of exit from a domain (i.e., the
time of first passage to the complement of the domain) of time-homogeneous
continuous-time Markov chains. We prove that: (i) the computed approximations
bound the measures from below; (ii) the total variation distances between the
approximations and the measures decrease monotonically as states are added to
the truncation; and (iii) the scheme converges, in the sense that, as the
truncation tends to the entire state space, the total variation distances tend
to zero. Furthermore, we give a computable bound on the total variation
distance between the exit distribution and its approximation, and we delineate
the cases in which the bound is sharp. We also revisit the related finite state
projection scheme and give a comprehensive account of its theoretical
properties. We demonstrate the use of the ETFSP scheme by applying it to two
biological examples: the computation of the first passage time associated with
the expression of a gene, and the fixation times of competing species subject
to demographic noise.


  Cell division site positioning is precisely regulated to generate correctly
sized and shaped daughters. We uncover a novel strategy to position the FtsZ
cytokinetic ring at midcell in the social bacterium Myxococcus xanthus. PomX,
PomY and the nucleoid-binding ParA/MinD ATPase PomZ self-assemble forming a
large nucleoid-associated complex that localizes at the division site before
FtsZ to directly guide and stimulate division. PomXYZ localization is generated
through self-organized biased random motion on the nucleoid towards midcell and
constrained motion at midcell. Experiments and theory show that PomXYZ motion
is produced by diffusive PomZ fluxes on the nucleoid into the complex. Flux
differences scale with the intracellular asymmetry of the complex and are
converted into a local PomZ concentration gradient across the complex with
translocation towards the higher PomZ concentration. At midcell, fluxes
equalize resulting in constrained motion. Flux-based mechanisms may represent a
general paradigm for positioning of macromolecular structures in bacteria.


  We use Markov state models (MSMs) to analyze the dynamics of a
$\beta$-hairpin-forming peptide in Monte Carlo (MC) simulations with
interacting protein crowders, for two different types of crowder proteins
[bovine pancreatic trypsin inhibitor (BPTI) and GB1]. In these systems, at the
temperature used, the peptide can be folded or unfolded and bound or unbound to
crowder molecules. Four or five major free-energy minima can be identified. To
estimate the dominant MC relaxation times of the peptide, we build MSMs using a
range of different time resolutions or lag times. We show that stable
relaxation-time estimates can be obtained from the MSM eigenfunctions through
fits to autocorrelation data. The eigenfunctions remain sufficiently accurate
to permit stable relaxation-time estimation down to small lag times, at which
point simple estimates based on the corresponding eigenvalues have large
systematic uncertainties. The presence of the crowders have a stabilizing
effect on the peptide, especially with BPTI crowders, which can be attributed
to a reduced unfolding rate $k_\text{u}$, while the folding rate $k_\text{f}$
is left largely unchanged.


  A state wide Human Respiratory Syncytial Virus (HRSV) surveillance system was
implemented in Florida in 1999 to support clinical decision-making for
prophylaxis of premature infants. The research presented in this paper
addresses the problem of fitting real data collected by the Florida HRSV
surveillance system by using a periodic SEIRS mathematical model. A sensitivity
and cost-effectiveness analysis of the model is done and an optimal control
problem is formulated and solved with treatment as the control variable.


  The performance of a molecular motor, characterized by its power output and
energy efficiency, is investigated in the motor design space spanned by the
stepping rate function and the motor-track interaction potential. Analytic
results and simulations show that a gating mechanism that restricts forward
stepping in a narrow window in configuration space is needed for generating
high power at physiologically relevant loads. By deriving general
thermodynamics laws for nonequilibrium motors, we find that the maximum torque
(force) at stall is less than its theoretical limit for any realistic
motor-track interactions due to speed fluctuations. Our study reveals a
tradeoff for the motor- track interaction: while a strong interaction generates
a high power output for forward steps, it also leads to a higher probability of
wasteful spontaneous back steps. Our analysis and simulations show that this
tradeoff sets a fundamental limit to the maximum motor efficiency in the
presence of spontaneous back steps, i.e., loose-coupling. Balancing this
tradeoff leads to an optimal design of the motor-track interaction for
achieving a maximum efficiency close to 1 for realistic motors that are not
perfectly coupled with the energy source.Comparison with existing data and
suggestions for future experiments are discussed.


  We apply matrix completion methods for haplotype assembly from NGS reads to
develop the new HapSVT, HapNuc, and HapOPT algorithms. This is performed by
applying a mathematical model to convert the reads to an incomplete matrix and
estimating unknown components. This process is followed by quantizing and
decoding the completed matrix in order to estimate haplotypes. These algorithms
are compared to the state-of-the-art algorithms using simulated data as well as
the real fosmid data. It is shown that the SNP missing rate and the haplotype
block length of the proposed HapOPT are better than those of HapCUT2 with
comparable accuracy in terms of reconstruction rate and switch error rate. A
program implementing the proposed algorithms in MATLAB is freely available at
https://github.com/smajidian/HapMC.


  Reaction-diffusion systems with a Lotka-Volterra-type reaction term, also
known as competition-diffusion systems, have been used to investigate the
dynamics of the competition among $m$ ecological species for a limited resource
necessary to their survival and growth. Notwithstanding their rather simple
mathematical structure, such systems may display quite interesting behaviours.
In particular, while for $m=2$ no coexistence of the two species is usually
possible, if $m \ge 3$ we may observe coexistence of all or a subset of the
species, sensitively depending on the parameter values. Such coexistence can
take the form of very complex spatio-temporal patterns and oscillations.
  Unfortunately, at the moment there are no known tools for a complete
analytical study of such systems for $m \ge 3$. This means that establishing
general criteria for the occurrence of coexistence appears to be very hard. In
this paper we will instead give some criteria for the non-coexistence of
species, motivated by the ecological problem of the invasion of an ecosystem by
an exotic species. We will show that when the environment is very favourable to
the invading species the invasion will always be successful and the native
species will be driven to extinction. On the other hand, if the environment is
not favourable enough, the invasion will always fail.


  Recent work indicates that twist-bend coupling plays an important role in DNA
micromechanics. Here we investigate its effect on bent DNA. We provide an
analytical solution of the minimum-energy shape of circular DNA, showing that
twist-bend coupling induces sinusoidal twist waves. This solution is in
excellent agreement with both coarse-grained simulations of minicircles and
nucleosomal DNA data, which is bent and wrapped around histone proteins in a
superhelical conformation. Our analysis shows that the observed twist
oscillation in nucleosomal DNA, so far attributed to the interaction with the
histone proteins, is an intrinsic feature of free bent DNA, and should be
observable in other protein-DNA complexes.


  Features of the CRISPR-Cas system, in which bacteria integrate small segments
of phage genome (spacers) into their DNA to neutralize future attacks, suggest
that its effect is not limited to individual bacteria but may control the fate
and structure of whole populations. Emphasizing the population-level impact of
the CRISPR-Cas system, recent experiments show that some bacteria regulate
CRISPR-associated genes via the quorum sensing (QS) pathway. Here we present a
model that shows that from the highly stochastic dynamics of individual spacers
under QS control emerges a rank-abundance distribution of spacers that is
time-invariant, a surprising prediction that we test with dynamic
spacer-tracking data from literature. This distribution depends on the state of
the competing phage-bacteria population, which due to QS-based regulation may
coexist in multiple stable states that vary significantly in their
phage-to-bacterium ratio, a widely used ecological measure to characterize
microbial systems.


  A variety of complex systems exhibit different types of relationships
simultaneously that can be modeled by multiplex networks. A typical problem is
to determine the community structure of such systems that, in general, depend
on one or more parameters to be tuned. In this study we propose one measure,
grounded on information theory, to find the optimal value of the relax rate
characterizing Multiplex Infomap, the generalization of the Infomap algorithm
to the realm of multilayer networks. We evaluate our methodology on synthetic
networks, to show that the most representative community structure can be
reliably identified when the most appropriate relax rate is used. Capitalizing
on these results, we use this measure to identify the most reliable meso-scale
functional organization in the human protein-protein interaction multiplex
network and compare the observed clusters against a collection of independently
annotated gene sets from the Molecular Signatures Database (MSigDB). Our
analysis reveals that modules obtained with the optimal value of the relax rate
are biologically significant and, remarkably, with higher functional content
than the ones obtained from the aggregate representation of the human proteome.
Our framework allows us to characterize the meso-scale structure of those
multilayer systems whose layers are not explicitly interconnected each other --
as in the case of edge-colored models -- the ones describing most biological
networks, from proteomes to connectomes.


  The effective representation of proteins is a crucial task that directly
affects the performance of many bioinformatics problems. Related proteins
usually bind to similar ligands. Chemical characteristics of ligands are known
to capture the functional and mechanistic properties of proteins suggesting
that a ligand based approach can be utilized in protein representation. In this
study, we propose SMILESVec, a SMILES-based method to represent ligands and a
novel method to compute similarity of proteins by describing them based on
their ligands. The proteins are defined utilizing the word-embeddings of the
SMILES strings of their ligands. The performance of the proposed protein
description method is evaluated in protein clustering task using TransClust and
MCL algorithms. Two other protein representation methods that utilize protein
sequence, BLAST and ProtVec, and two compound fingerprint based protein
representation methods are compared. We showed that ligand-based protein
representation, which uses only SMILES strings of the ligands that proteins
bind to, performs as well as protein-sequence based representation methods in
protein clustering. The results suggest that ligand-based protein description
can be an alternative to the traditional sequence or structure based
representation of proteins and this novel approach can be applied to different
bioinformatics problems such as prediction of new protein-ligand interactions
and protein function annotation.


  Gillespie's direct method for stochastic simulation of chemical kinetics is a
staple of computational systems biology research. However, the algorithm
requires explicit enumeration of all reactions and all chemical species that
may arise in the system. In many cases, this is not feasible due to the
combinatorial explosion of reactions and species in biological networks.
Rule-based modeling frameworks provide a way to exactly represent networks
containing such combinatorial complexity, and generalizations of Gillespie's
direct method have been developed as simulation engines for rule-based modeling
languages. Here, we provide both a high-level description of the algorithms
underlying the simulation engines, termed network-free simulation algorithms,
and how they have been applied in systems biology research. We also define a
generic rule-based modeling framework and describe a number of technical
details required for adapting Gillespie's direct method for network-free
simulation. Finally, we briefly discuss potential avenues for advancing
network-free simulation and the role they continue to play in modeling
dynamical systems in biology.


  Tree balance plays an important role in different research areas like
theoretical computer science and mathematical phylogenetics. For example, it
has long been known that under the Yule model, a pure birth process, imbalanced
trees are more likely than balanced ones. Also, concerning ordered search
trees, more balanced ones allow for more efficient data structuring than
imbalanced ones. Therefore, different methods to measure the balance of trees
were introduced. The Sackin index is one of the most frequently used measures
for this purpose. In many contexts, statements about the minimal and maximal
values of this index have been discussed, but formal proofs have only been
provided for some of them, and only in the context of ordered binary (search)
trees, not for general rooted trees. Moreover, while the number of trees with
maximal Sackin index as well as the number of trees with minimal Sackin index
when the number of leaves is a power of 2 are relatively easy to understand,
the number of trees with minimal Sackin index for all other numbers of leaves
has been completely unknown. In this manuscript, we extend the findings on
trees with minimal and maximal Sackin indices from the literature on ordered
trees and subsequently use our results to provide formulas to explicitly
calculate the numbers of such trees. We also extend previous studies by
analyzing the case when the underlying trees need not be binary. Finally, we
use our results to contribute both to the phylogenetic as well as the computer
scientific literature by using the new findings on Sackin minimal and maximal
trees in order to derive formulas to calculate the number of both minimal and
maximal phylogenetic trees as well as minimal and maximal ordered trees both in
the binary and non-binary settings. All our results have been implemented in
the Mathematica package SackinMinimizer, which has been made publicly
available.


  We propose an optimized parameter set for protein secondary structure
prediction using three layer feed forward back propagation neural network. The
methodology uses four parameters viz. encoding scheme, window size, number of
neurons in the hidden layer and type of learning algorithm. The input layer of
the network consists of neurons changing from 3 to 19, corresponding to
different window sizes. The hidden layer chooses a natural number from 1 to 20
as the number of neurons. The output layer consists of three neurons, each
corresponding to known secondary structural classes viz. alpha helix, beta
strands and coils respectively. It also uses eight different learning
algorithms and nine encoding schemes. Exhaustive experiments were performed
using non-homologues dataset. The experimental results were compared using
performance measures like Q3, sensitivity, specificity, Mathew correlation
coefficient and accuracy. The paper also discusses the process of obtaining a
stabilized cluster of 2530 records from a collection of 11340 records. The
graphs of these stabilized clusters of records with respect to accuracy are
concave, convergence is monotonic increasing and rate of convergence is
uniform. The paper gives BLOSUM62 as the encoding scheme, 19 as the window
size, 19 as the number of neurons in the hidden layer and One- Step Secant as
the learning algorithm with the highest accuracy of 78%. These parameter values
are proposed as the optimized parameter set for the three layer feed forward
back propagation neural network for the protein secondary structure predictionv


  Since social interactions have been shown to lead to symmetric clusters, we
propose here that symmetries play a key role in epidemic modeling. Mathematical
models on d-ary tree graphs were recently shown to be particularly effective
for modeling epidemics in simple networks [Seibold & Callender, 2016]. To
account for symmetric relations, we generalize this to a new type of networks
modeled on d-cliqued tree graphs, which are obtained by adding edges to regular
d-trees to form d-cliques. This setting gives a more realistic model for
epidemic outbreaks originating, for example, within a family or classroom and
which could reach a population by transmission via children in schools.
Specifically, we quantify how an infection starting in a clique (e.g. family)
can reach other cliques through the body of the graph (e.g. public places).
Moreover, we propose and study the notion of a safe zone, a subset that has a
negligible probability of infection.


  Feedback regulation in biochemical systems is fundamental to homeostasis,
with failure causing disease or death. Recent work has found that cooperation
between feedback and buffering---the use of reservoirs of molecules to maintain
molecular concentrations---is often critical for biochemical regulation, and
that buffering can act as a derivative or lead controller. However, buffering
differs from derivative feedback in important ways: it is not typically limited
by stability constraints on the parallel feedback loop, for some signals it
acts instead as a low-pass filter, and it can change the location of
disturbances in the closed-loop system. Here, we propose a frequency-domain
framework for studying the regulatory properties of buffer-feedback systems. We
determine standard single-output closed-loop transfer functions and discuss
loop-shaping properties. We also derive novel fundamental limits for
buffer-feedback regulation, which show that buffering and removal processes can
reduce the fundamental limits on feedback regulation. We apply the framework to
study the regulation for glycolysis (anaerobic metabolism) with creatine
phosphate buffering.


  Transient compartments have been recently shown to be able to maintain
functional replicators in the context of prebiotic studies. Motivated by this
experiment, we show that a broad class of selection dynamics is able to achieve
this goal. We identify two key parameters, the relative amplification of
non-active replicators (parasites) and the size of compartments. Since the
basic ingredients of our model are the competition between a host and its
parasite, and the diversity generated by small size compartments, our results
are relevant to various phage-bacteria or virus-host ecology problems.


  We propose the $S$-leaping algorithm for the acceleration of Gillespie's
stochastic simulation algorithm that combines the advantages of the two main
accelerated methods; the $\tau$-leaping and $R$-leaping algorithms. These
algorithms are known to be efficient under different conditions; the
$\tau$-leaping is efficient for non-stiff systems or systems with partial
equilibrium, while the $R$-leaping performs better in stiff system thanks to an
efficient sampling procedure. However, even a small change in a system's set up
can critically affect the nature of the simulated system and thus reduce the
efficiency of an accelerated algorithm. The proposed algorithm combines the
efficient time step selection from the $\tau$-leaping with the effective
sampling procedure from the $R$-leaping algorithm. The $S$-leaping is shown to
maintain its efficiency under different conditions and in the case of large and
stiff systems or systems with fast dynamics, the $S$-leaping outperforms both
methods. We demonstrate the performance and the accuracy of the $S$-leaping in
comparison with the $\tau$-leaping and $R$-leaping on a number of benchmark
systems involving biological reaction networks.


  A common problem in phylogenetics is to try to infer a species phylogeny from
gene trees. We consider different variants of this problem. The first variant,
called Unrestricted Minimal Episodes Inference, aims at inferring a species
tree based on a model with speciation and duplication where duplications are
clustered in duplication episodes. The goal is to minimize the number of such
episodes. The second variant, Parental Hybridization, aims at inferring a
species \emph{network} based on a model with speciation and reticulation. The
goal is to minimize the number of reticulation events. It is a variant of the
well-studied Hybridization Number problem with a more generous view on which
gene trees are consistent with a given species network. We show that these
seemingly different problems are in fact closely related and can, surprisingly,
both be solved in polynomial time, using a structure we call "beaded trees".
However, we also show that methods based on these problems have to be used with
care because the optimal species phylogenies always have a restricted form. To
mitigate this problem, we introduce a new variant of Unrestricted Minimal
Episodes Inference that minimizes the duplication episode depth. We prove that
this new variant of the problem can also be solved in polynomial time


  One of the ultimate goals in biology is to understand the design principles
of biological systems. Such principles, if they exist, can help us better
understand complex, natural biological systems and guide the engineering of de
novo ones. Towards deciphering design principles, in silico evolution of
biological systems with proper abstraction is a promising approach. Here, we
demonstrate the application of in silico evolution combined with rule-based
modelling for exploring design principles of cellular signaling networks. This
application is based on a computational platform, called BioJazz, which allows
in silico evolution of signaling networks with unbounded complexity. We provide
a detailed introduction to BioJazz architecture and implementation and describe
how it can be used to evolve and/or design signaling networks with defined
dynamics. For the latter, we evolve signaling networks with switch-like
response dynamics and demonstrate how BioJazz can result in new biological
insights on network structures that can endow bistable response dynamics. This
example also demonstrated both the power of BioJazz in evolving and designing
signaling networks and its limitations at the current stage of development.


  We establish a biophysical model for the dynamics of lipid vesicles exposed
to surfactants. The solubilization of the lipid membrane due to the insertion
of surfactant molecules induces a reduction of membrane surface area at almost
constant vesicle volume. This results in a rate-dependent increase of membrane
tension and leads to the opening of a micron-sized pore. We show that
solubilization kinetics due to surfactants can determine the regimes of pore
dynamics: either the pores open and reseal within a second (short-lived pore),
or the pore stays open up to a few minutes (long-lived pore). First, we
validate our model with previously published experimental measurements of pore
dynamics. Then, we investigate how the solubilization kinetics and membrane
properties affect the dynamics of the pore and construct a phase diagram for
short and long-lived pores. Finally, we examine the dynamics of sequential pore
openings and show that cyclic short-lived pores occur at a period inversely
proportional to the solubilization rate. By deriving a theoretical expression
for the cycle period, we provide an analytic tool to measure the solubilization
rate of lipid vesicles by surfactants. Our findings shed light on some
fundamental biophysical mechanisms that allow simple cell-like structures to
sustain their integrity against environmental stresses, and have the potential
to aid the design of vesicle-based drug delivery systems.


  Evolution occurs in populations of reproducing individuals. Reproduction
depends on the payoff a strategy receives. The payoff depends on the
environment that may change over time, on intrinsic uncertainties, and on other
sources of randomness. These temporal variations in the payoffs can affect
which traits evolve. Understanding evolutionary game dynamics that are affected
by varying payoffs remains difficult. Here we study the impact of arbitrary
amplitudes and covariances of temporally varying payoffs on the dynamics. The
evolutionary dynamics may be "unfair", meaning that, on average, two coexisting
strategies may persistently receive different payoffs. This mechanism can
induce an anomalous coexistence of cooperators and defectors in the Prisoner's
Dilemma, and an unexpected selection reversal in the Hawk-Dove game.


  The use of drug combinations, termed polypharmacy, is common to treat
patients with complex diseases and co-existing conditions. However, a major
consequence of polypharmacy is a much higher risk of adverse side effects for
the patient. Polypharmacy side effects emerge because of drug-drug
interactions, in which activity of one drug may change if taken with another
drug. The knowledge of drug interactions is limited because these complex
relationships are rare, and are usually not observed in relatively small
clinical testing. Discovering polypharmacy side effects thus remains an
important challenge with significant implications for patient mortality. Here,
we present Decagon, an approach for modeling polypharmacy side effects. The
approach constructs a multimodal graph of protein-protein interactions,
drug-protein target interactions, and the polypharmacy side effects, which are
represented as drug-drug interactions, where each side effect is an edge of a
different type. Decagon is developed specifically to handle such multimodal
graphs with a large number of edge types. Our approach develops a new graph
convolutional neural network for multirelational link prediction in multimodal
networks. Decagon predicts the exact side effect, if any, through which a given
drug combination manifests clinically. Decagon accurately predicts polypharmacy
side effects, outperforming baselines by up to 69%. We find that it
automatically learns representations of side effects indicative of
co-occurrence of polypharmacy in patients. Furthermore, Decagon models
particularly well side effects with a strong molecular basis, while on
predominantly non-molecular side effects, it achieves good performance because
of effective sharing of model parameters across edge types. Decagon creates
opportunities to use large pharmacogenomic and patient data to flag and
prioritize side effects for follow-up analysis.


  Advancements in genomic research such as high-throughput sequencing
techniques have driven modern genomic studies into "big data" disciplines. This
data explosion is constantly challenging conventional methods used in genomics.
In parallel with the urgent demand for robust algorithms, deep learning has
succeeded in a variety of fields such as vision, speech, and text processing.
Yet genomics entails unique challenges to deep learning since we are expecting
from deep learning a superhuman intelligence that explores beyond our knowledge
to interpret the genome. A powerful deep learning model should rely on
insightful utilization of task-specific knowledge. In this paper, we briefly
discuss the strengths of different deep learning models from a genomic
perspective so as to fit each particular task with a proper deep architecture,
and remark on practical considerations of developing modern deep learning
architectures for genomics. We also provide a concise review of deep learning
applications in various aspects of genomic research, as well as pointing out
potential opportunities and obstacles for future genomics applications.


  Biomedical sciences are increasingly recognising the relevance of gene
co-expression-networks for analysing complex-systems, phenotypes or diseases.
When the goal is investigating complex-phenotypes under varying conditions, it
comes naturally to employ comparative network methods. While approaches for
comparing two networks exist, this is not the case for multiple networks. Here
we present a method for the systematic comparison of an unlimited number of
networks: Co-expression Differential Network Analysis (CoDiNA) for detecting
links and nodes that are common, specific or different to the networks.
Applying CoDiNA to a neurogenesis study identified genes for neuron
differentiation. Experimentally overexpressing one candidate resulted in
significant disturbance in the underlying neurogenesis' gene regulatory
network. We compared data from adults and children with active tuberculosis to
test for signatures of HIV. We also identified common and distinct network
features for particular cancer types with CoDiNA. These studies show that
CoDiNA successfully detects genes associated with the diseases.


  Lipid phase heterogeneity in the plasma membrane is thought to be crucial for
many aspects of cell signaling, but the physical basis of participating
membrane domains such as "lipid rafts" remains controversial. Here we consider
a lattice model yielding a phase diagram that includes several states proposed
to be relevant for the cell membrane, including microemulsion - which can be
related to membrane curvature - and Ising critical behavior. Using a neural
network-based machine learning approach, we compute the full phase diagram of
this lattice model. We analyze selected regions of this phase diagram in the
context of a signaling initiation event in mast cells: recruitment of the
membrane-anchored tyrosine kinase Lyn to a cluster of transmembrane of
IgE-Fc{\epsilon}RI receptors. We find that model membrane systems in
microemulsion and Ising critical states can mediate roughly equal levels of
kinase recruitment (binding energy ~ -0.6 kBT), whereas a membrane near a
tricritical point can mediate much stronger kinase recruitment (-1.7 kBT). By
comparing several models for lipid heterogeneity within a single theoretical
framework, this work points to testable differences between existing models. We
also suggest the tricritical point as a new possibility for the basis of
membrane domains that facilitate preferential partitioning of signaling
components.


  In this paper, we study the interplay between individual behaviors and
epidemic spreading in a dynamical network. We distribute agents on a
square-shaped region with periodic boundary conditions. Every agent is regarded
as a node of the network and a wireless link is established between two agents
if their geographical distance is less than a certain radius. At each time,
every agent assesses the epidemic situation and make decisions on whether it
should stay in or leave its current place. An agent will leave its current
place with a speed if the number of infected neighbors reaches or exceeds a
critical value $E$. Owing to the movement of agents, the network's structure is
dynamical. Interestingly, we find that there exists an optimal value of $E$
leading to the maximum epidemic threshold. This means that epidemic spreading
can be effectively controlled by risk-averse migration. Besides, we find that
the epidemic threshold increases as the recovering rate increases, decreases as
the contact radius increases, and is maximized by an optimal moving speed. Our
findings offer a deeper understanding of epidemic spreading in dynamical
networks.


  We present a protocol for building, validating and simulating models of
signal transduction networks. These networks are challenging modelling targets
due to the combinatorial complexity and sparse data, which have made it a major
challenge even to formalise the current knowledge. To address this, the
community has developed methods to model biomolecular reaction networks based
on site dynamics. The strength of this approach is that reactions and states
can be defined at variable resolution, which makes it possible to adapt the
model resolution to the empirical data. This improves both scalability and
accuracy, making it possible to formalise large models of signal transduction
networks. Here, we present a method to build and validate large models of
signal transduction networks. The workflow is based on rxncon, the
reaction-contingency language. In a five-step process, we create a mechanistic
network model, convert it into an executable Boolean model, use the Boolean
model to evaluate and improve the network, and finally export the rxncon model
into a rule based format. We provide an introduction to the rxncon language and
an annotated, step-by-step protocol for the workflow. Finally, we create a
small model of the insulin signalling pathway to illustrate the protocol,
together with some of the challenges - and some of their solutions - in
modelling signal transduction.


  Organization and maintenance of the chromosomal DNA in living cells strongly
depends on the DNA interactions with a plethora of DNA-binding proteins.
Single-molecule studies show that formation of nucleoprotein complexes on DNA
by such proteins is frequently subject to force and torque constraints applied
to the DNA. Although the existing experimental techniques allow to exert these
type of mechanical constraints on individual DNA biopolymers, their exact
effects in regulation of DNA-protein interactions are still not completely
understood due to the lack of systematic theoretical methods able to
efficiently interpret complex experimental observations. To fill this gap, we
have developed a general theoretical framework based on the transfer-matrix
calculations that can be used to accurately describe behaviour of DNA-protein
interactions under force and torque constraints. Potential applications of the
constructed theoretical approach are demonstrated by predicting how these
constraints affect the DNA-binding properties of different types of
architectural proteins. Obtained results provide important insights into
potential physiological functions of mechanical forces in the chromosomal DNA
organization by architectural proteins as well as into single-DNA manipulation
studies of DNA-protein interactions.


  The purpose of this paper is to provide an algorithm for detecting and
tracking astral MTs in a fully automated way and supply a description of their
dynamic behaviour. For the algorithm testing, a dataset of stacks (i.e.
time-lapse image sequences), acquired with a confocal microscope, has been
employed. Cells were treated with two different drugs, nocodazole and taxol, in
order to explore their effect on microtubule dynamic instability.


  We develop a robust coarse-grained model for single and double stranded DNA
by representing each nucleotide by three interaction sites (TIS) located at the
centers of mass of sugar, phosphate, and base. The resulting TIS model includes
base-stacking, hydrogen bond, and electrostatic interactions as well as
bond-stretching and bond angle potentials that account for the polymeric nature
of DNA. The choices of force constants for stretching and the bending
potentials were guided by a Boltzmann inversion procedure using a large
representative set of DNA structures extracted from the Protein Data Bank. Some
of the parameters in the stacking interactions were calculated using a learning
procedure, which ensured that the experimentally measured melting temperatures
of dimers are faithfully reproduced. Without any further adjustments, the
calculations based on the TIS model reproduces the experimentally measured salt
and sequence dependence of the size of single stranded DNA (ssDNA), as well as
the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon
application of mechanical force the extension of poly(dA) exhibits a plateau,
which we trace to the formation of stacked helical domains. In contrast, the
force-extension curve (FEC) of poly(dT) is entropic in origin, and could be
described by a standard polymer model. We also show that the persistence length
of double stranded DNA is consistent with the prediction based on the worm-like
chain. The persistence length, which decreases with increasing salt
concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for
stiff polyelectrolyte chains near the rod limit. The range of applications,
which did not require adjusting any parameter after the initial construction
based solely on PDB structures and melting profiles of dimers, attests to the
transferability and robustness of the TIS model for ssDNA and dsDNA.


  The work reported here aims to address the effects of time-dependent
parameters and stochasticity on decision-making in biological systems. We
achieve this by extending previous studies that resorted to simple normal
forms. Yet, we focus primarily on the issue of the system's sensitivity to
initial conditions in the presence of different noise distributions. In
addition, we assess the impact of two-way sweeping through the critical region
of a canonical Pitchfork bifurcation with a constant external asymmetry. The
parallel with decision-making in bio-circuits is performed on this simple
system since it is equivalent in its available states and dynamics to more
complex genetic circuits. Overall, we verify that rate-dependent effects are
specific to particular initial conditions. Information processing for each
starting state is affected by the balance between sweeping speed through
critical regions, and the type of fluctuations added. For a heavy-tail noise,
forward-reverse dynamic bifurcations are more efficient in processing the
information contained in external signals, when compared to the system relying
on escape dynamics, if it starts at an attractor not favoured by the asymmetry
and, in conjunction, if the sweeping amplitude is large.


  We develop open source infection risk calculators for patients and healthcare
professionals as apps for hospital acquired infections (during child-delivery)
and sexually transmitted infections (like HIV). Advanced versions of ehealth in
non-communicable diseases do not apply to epidemiology much. There is, however,
no infection risk calculator in the Polish Internet so far, despite the
existence of data that may be applied to create such a tool.
  The algorithms involve data from Information Systems (like HIS in hospitals)
and surveys by applying mathematical modelling, Bayesian inference, logistic
regressions, covariance analysis and social network analysis. Finally, user may
fill or import data from Information System to obtain risk assessment and test
different settings to learn overall risk.
  The most promising risk calculator is developed for Healthcare-associated
infections in modes for patient hospital sanitary inspection. The most extended
version for hospital epidemiologists may include many layers of hospital
interactions by agent-based modeling. Simplified version of calculator is
dedicated to patients that require personalized hospitalization history of
pregnancy described by questions represented by quantitative and qualitative
variables. Patients receive risk assessment from interactive web application
with additional description about modifiable risk factors.
  We also provide solution for sexually transmitted infections like HIV. The
results of calculations with meaningful description and percentage chances are
presented in real-time to interested users. Finally, user fills the form to
obtain risk assessment for given settings.


  Many studies have shown that cells use temporal dynamics of signalling
molecules to encode information. One particular class of temporal dynamics is
persistent and transient signals, i.e. signals of long and short durations
respectively. It has been shown that the coherent type-1 feedforward loop with
an AND logic at the output (or C1-FFL for short) can be used to discriminate a
persistent input signal from a transient one. This has been done by modelling
the C1-FFL, and then use the model to show that persistent and transient input
signals give, respectively, a non-zero and zero output. Instead of assuming the
structure of C1-FFL, this paper shows that it is possible to deduce the C1-FFL
model from the requirement of discriminating a persistent signal. We do this by
first formulating a statistical detection problem of distinguishing persistent
signals from transient ones. The solution of the detection problem is to
compute the log-likelihood ratio of observing a persistent signal to a
transient signal. We show that, if this log-likelihood ratio is positive, which
happens when the signal is likely to be persistent, then it can be
approximately computed by a C1-FFL. Although the capability of C1-FFL to
discriminate persistent signals is known, this paper adds an information
processing interpretation on how a C1-FFL works as a detector of persistent
signals.


  In many animals parental gametes unite to form a zygote that develops into an
adult with gonads that, in turn, produce gametes. Interruption of this germinal
cycle by prezygotic or postzygotic reproductive barriers can result in two
independent cycles, each with the potential to evolve into a new species. When
the speciation process is complete, members of each species are fully
reproductively isolated from those of the other. During speciation a primary
barrier may be supported and eventually superceded by a later appearing
secondary barrier. For those holding certain cases of prezygotic isolation to
be primary (e.g. elephant cannot copulate with mouse), the onus is to show that
they had not been preceded over evolutionary time by periods of postzygotic
hybrid inviability (genically determined) or sterility (genically or
chromosomally determined). Likewise, the onus is upon those holding cases of
hybrid inviability to be primary (e.g. Dobzhansky-Muller epistatic
incompatibilities), to show that they had not been preceded by periods, however
brief, of hybrid sterility. The latter, when acting as a sympatric barrier
causing reproductive isolation, can only be primary. In many cases, hybrid
sterility may result from incompatibilities between parental chromosomes that
attempt to pair during meiosis in the gonad of their offspring
(Winge-Crowther-Bateson incompatibilities). While WCB incompatibilities have
long been observed on a microscopic scale, there is growing evidence for a role
of dispersed finer DNA sequence differences.


  During the initiation stage of protein synthesis, a ribosomal initiation
complex (IC) is assembled on a messenger RNA (mRNA) template. In bacteria, the
speed and accuracy of this assembly process are regulated by the complementary
activities of three essential initiation factors (IFs). Selection of an
authentic N-formylmethionyl-transfer RNA (fMet-tRNA\textsuperscript{fMet}) and
the canonical, triplet-nucleotide mRNA start codon are crucial events during
assembly of a canonical, ribosomal 70S IC. Mis-initiation due to the aberrant
selection of an elongator tRNA or a non-canonical start codon are rare events
that result in the assembly of a pseudo 70S IC or a non-canonical 70S IC,
respectively. Here, we have developed a theoretical model for the stochastic
kinetics of canonical-, pseudo-, and non-canonical 70S IC assembly that
includes all of the major steps of the IC assembly process that have been
observed and characterized in ensemble kinetic-, single-molecule kinetic-, and
structural studies of the fidelity of translation initiation. Specifically, we
use the rates of the individual steps in the IC assembly process and the
formalism of first-passage times to derive exact analytical expressions for the
probability distributions for the assembly of canonical-, pseudo- and
non-canonical 70S ICs. In order to illustrate the power of this analytical
approach, we compare the theoretically predicted first-passage time
distributions with the corresponding computer simulation data. We also compare
the mean times required for completion of these assemblies with experimental
estimates. In addition to generating new, testable hypotheses, our theoretical
model can also be easily extended as new experimental 70S IC assembly data
become available, thereby providing a versatile tool for interpreting these
data and developing advanced models of the mechanism and regulation of
translation initiation.


  PEBPs (PhosphatidylEthanolamine Binding Proteins) form a protein family
widely present in the living world since they are encountered in
microorganisms, plants and animals. In all organisms PEBPs appear to regulate
important mechanisms that govern cell cycle, proliferation, differentiation and
motility. In humans, three PEBPs have been identified, namely PEBP1, PEBP2 and
PEBP4. PEBP1 and PEBP4 are the most studied as they are implicated in the
development of various cancers. PEBP2 is specific of testes in mammals and was
essentially studied in rats and mice where it is very abundant. A lot of
information has been gained on PEBP1 also named RKIP (Raf Kinase Inhibitory
protein) due to its role as a metastasis suppressor in cancer. PEBP1 was also
demonstrated to be implicated in Alzheimers disease, diabetes and
nephropathies. Furthermore, PEBP1 was described to be involved in many cellular
processes, among them are signal transduction, inflammation, cell cycle,
proliferation, adhesion, differentiation, apoptosis, autophagy, circadian
rhythm and mitotic spindle checkpoint. On the molecular level, PEBP1 was shown
to regulate several signaling pathways such as Raf/MEK/ERK, NFkB,
PI3K/Akt/mTOR, p38, Notch and Wnt. PEBP1 acts by inhibiting most of the kinases
of these signaling cascades. Moreover, PEBP1 is able to bind to a variety of
small ligands such as ATP, phospholipids, nucleotides, flavonoids or drugs.
Considering PEBP1 is a small cytoplasmic protein (21kDa), its involvement in so
many diseases and cellular mechanisms is amazing. The aim of this review is to
highlight the molecular systems that are common to all these cellular
mechanisms in order to decipher the specific role of PEBP1. Recent discoveries
enable us to propose that PEBP1 is a modulator of molecular interactions that
control signal transduction during membrane and cytoskeleton reorganization.


  We study the spreading of cooperative infections in an empirical temporal
network of contacts between people, including health care workers and patients,
in a hospital. The system exhibits a phase transition leading to one or several
endemic branches, depending on the connectivity pattern and the temporal
correlations. There are two endemic branches in the original setting and the
non-cooperative case. However, the cooperative interaction between infections
reinforces the upper branch, leading to a smaller epidemic threshold and a
higher probability for having a big outbreak. We show the microscopic
mechanisms leading to these differences, characterize three different risks,
and use the influenza features as an example for this dynamics.


  Multiple scales in metapopulations can give rise to paradoxical behaviour: in
a conceptual model for a public goods game, the species associated with a
fitness cost due to the public good production can be stabilised in the
well-mixed limit due to the mere existence of these scales. The scales in this
model involve a length scale corresponding to separate patches, coupled by
mobility, and separate time scales for reproduction and interaction with a
local environment. Contrary to the well-mixed high mobility limit, we find that
for low mobilities, the interaction rate progressively stabilises this species
due to stochastic effects, and that the formation of spatial patterns is not
crucial for this stabilisation.


  The dissolved organic matter (DOM) composition in river ecosystems could
reflect the human impacts on the river ecosystem, and plays an important role
in the carbon cycling process. We collected water and phytoplankton samples at
107 sites in the Dongjiang River in two seasons to assess the impact of the
sub-catchments land use structure on the DOM composition. The results showed
that (1) the forested sub-catchments had higher humic-like C1 (16.45%) and C2
(25.04%) and lower protein-like C3 (22.57%) and C4 (35.95%) than urbanized and
mixed forest-agriculture sub-catchments, while the urbanized sub-catchments
showed an inverse trend (4.54%, 15.51%, 33.97% and 45.98%, respectively). (2)
The significant variation in the proportion of C1 and C4 between the dry and
rainy seasons was recorded in both the forested and the mixed
forest-agriculture sub-catchments (p<0.01), but only C4 showed an obvious
seasonal variation in the urbanized sub-catchments (p<0.01). While the DOM
composition was mainly related to the proportion of urbanized land and forested
land year-round (p<0.01), it had stronger correlation with forested land in the
dry season and urbanized land in the rainy season. (3) No significant
correlation between the DOM composition and the agricultural land proportion
was found in either season (p>0.05). Our findings indicated that the DOM
composition was strongly dependent on the land use structure of the
sub-catchments and varied seasonally, but the seasonal variation pattern could
be disturbed by human activities in the extensively urbanized catchments.
Notably, the higher C4 proportion compared with those in temperate rivers
implied a larger amount of CO2 released from the subtropical rivers into the
atmosphere when considering bioavailability.


  Inspired by DNA data of the human cytomegalovirus we propose a model of a
two-type parasite population distributed over its hosts. The parasite is
capable to persist in its host till the host dies, and to reinfect other hosts.
To maintain type diversity within a host, balancing selection is assumed.
  For a suitable parameter regime we show that in the limit of large host and
parasite populations the host state frequencies follow a dynamical system with
a globally stable equilibrium, guaranteeing that both types are maintained in
the parasite population for a long time on the host time scale.


  Small molecules are often found to exhibit extraordinarily diverse biological
activities. Metformin is one of them. It is widely used as anti-diabetic drug
for type-two diabetes. In addition to that, metformin hydrochloride shows
anti-tumour activities and increases the survival rate of patients suffering
from certain types of cancer namely colorectal, breast, pancreas and prostate
cancer. However, theoretical studies of structure and dynamics of metformin
have not yet been fully explored. In this work, we investigate the
characteristic structural and dynamical features of three mono-protonated forms
of metformin hydrochloride with the help of experiments, quantum chemical
calculations and atomistic molecular dynamics simulations. We validate our
force field by comparing simulation results to that of the experimental
findings. Nevertheless, we discover that the non-planar tautomeric form is the
most stable. Metformin forms strong hydrogen bonds with surrounding water
molecules and its solvation dynamics show unique features. Because of an
extended positive charge distribution, metformin possesses features of being a
permanent cationic partner toward several targets. We study its interaction and
binding ability with DNA using UV spectroscopy, circular dichroism, fluorimetry
and metadynamics simulation. We find a non-intercalating mode of interaction.
Metformin feasibly forms a minor/major groove-bound state within a few tens of
nanoseconds, preferably with AT rich domains. A significant decrease in the
free-energy of binding is observed when it binds to a minor groove of DNA.


  The so-called binary perfect phylogeny with persistent characters has
recently been thoroughly studied in computational biology as it is less
restrictive than the well known binary perfect phylogeny. Here, we focus on the
notion of (binary) persistent characters, i.e. characters that can be realized
on a phylogenetic tree by at most one $0 \rightarrow 1$ transition followed by
at most one $1 \rightarrow 0$ transition in the tree, and analyze these
characters under different aspects. First, we illustrate the connection between
persistent characters and Maximum Parsimony, where we characterize persistent
characters in terms of the first phase of the famous Fitch algorithm.
Afterwards we focus on the number of persistent characters for a given
phylogenetic tree. We show that this number solely depends on the balance of
the tree. To be precise, we develop a formula for counting the number of
persistent characters for a given phylogenetic tree based on an index of tree
balance, namely the Sackin index. Lastly, we consider the question of how many
(carefully chosen) binary characters together with their persistence status are
needed to uniquely determine a phylogenetic tree and provide an upper bound for
the number of characters needed.


  Histamine is a critical mediator of IgE/ cell-mediated anaphylaxis, a
neurotransmitter and a regulator of gastric acid secretion. Histamine is a
monoamine synthesized from the amino acid histidine through a reaction
catalyzed by the enzyme histidine decarboxylase (HDC), which removes carboxyl
group from histidine. Despite the importance of histamine, transcriptional
regulation of HDC gene expression in mammals is still poorly understood. In
this Review, we focus on discussing advances in the understanding of molecular
regulation of mammalian histamine synthesis.


  Albeit epidemic models have evolved into powerful predictive tools for the
spread of diseases and opinions, most assume memoryless agents and independent
transmission channels. We develop an infection mechanism that is endowed with
memory of past exposures and simultaneously incorporates the joint effect of
multiple infectious sources. Analytic equations and simulations of the
susceptible-infected-susceptible model in unstructured substrates reveal the
emergence of an additional phase that separates the usual healthy and endemic
ones. This intermediate phase shows fundamentally distinct characteristics, and
the system exhibits either excitability or an exotic variant of bistability.
Moreover, the transition to endemicity presents hybrid aspects. These features
are the product of an intricate balance between two memory modes and indicate
that non-Markovian effects significantly alter the properties of spreading
processes.


  Proteins are commonly used by biochemical industry for numerous processes.
Refining these proteins' properties via mutations causes stability effects as
well. Accurate computational method to predict how mutations affect protein
stability are necessary to facilitate efficient protein design. However,
accuracy of predictive models is ultimately constrained by the limited
availability of experimental data. We have developed mGPfusion, a novel
Gaussian process (GP) method for predicting protein's stability changes upon
single and multiple mutations. This method complements the limited experimental
data with large amounts of molecular simulation data. We introduce a Bayesian
data fusion model that re-calibrates the experimental and in silico data
sources and then learns a predictive GP model from the combined data. Our
protein-specific model requires experimental data only regarding the protein of
interest and performs well even with few experimental measurements. The
mGPfusion models proteins by contact maps and infers the stability effects
caused by mutations with a mixture of graph kernels. Our results show that
mGPfusion outperforms state-of-the-art methods in predicting protein stability
on a dataset of 15 different proteins and that incorporating molecular
simulation data improves the model learning and prediction accuracy.


  Background and Purpose: This study was designed to compare the effects of
high-intensity interval training (HIT) and moderate-intensity continuous
training (MOD) on functional recovery and cerebral plasticity during the first
2 weeks following cerebral ischemia. Methods: Rats were randomized as follows:
Control (n=15), SHAM (n=9), MCAO (n=13), MCAO-D1 (n=7), MOD (n=13) and HIT
(n=13). Incremental tests were performed at day 1 (D1) and 14 (D14) to identify
the running speed associated with the lactate threshold (SLT) and the maximal
speed (Smax). Functional tests were performed at D1, D7 and D14. Microglia
form, cytokines, p75NTR, KCC2 and NKCC1 expression were made at D15.
Results-HIT was more effective to improve the endurance performance than MOD
and induced a fast recovery of the impaired forelimb grip force. The Iba-1
positive cells with amoeboid form and the pro- and anti-inflammatory cytokine
expression were lower in HIT group, mainly in the ipsilesional hemisphere. A
p75NTR overexpression is observed on the ipsilesional side together with a
restored NKCC1/KCC2 ratio on the contralesional side. Conclusions-Low-volume
HIT based on lactate threshold appears to be more effective after cerebral
ischemia than work-matched MOD to improve aerobic fitness, grip strength and
might promote cerebral plasticity.


  Previously, a seven-cluster pattern claiming to be a universal one in
bacterial genomes has been reported. Keeping in mind the most popular theory of
chloroplast origin, we checked whether a similar pattern is observed in
chloroplast genomes. Surprisingly, eight cluster structure has been found, for
chloroplasts. The pattern observed for chloroplasts differs rather
significantly, from bacterial one, and from that latter observed for
cyanobacteria. The structure is provided by clustering of the fragments of
equal length isolated within a genome so that each fragment is converted in
triplet frequency dictionary with non-overlapping triplets with no gaps in
frame tiling. The points in 63-dimensional space were clustered due to elastic
map technique. The eight cluster found in chloroplasts comprises the fragments
of a genome bearing tRNA genes and exhibiting excessively high
$\mathsf{GC}$-content, in comparison to the entire genome.


  Experimentalists have observed phenotypic variability in isogenic bacteria
populations. We explore the hypothesis that in fluctuating environments this
variability is tuned to maximize a bacterium's expected log growth rate,
potentially aided by epigenetic markers that store information about past
environments. We show that, in a complex, memoryful environment, the maximal
expected log growth rate is linear in the instantaneous predictive
information---the mutual information between a bacterium's epigenetic markers
and future environmental states. Hence, under resource constraints, optimal
epigenetic markers are causal states---the minimal sufficient statistics for
prediction. This is the minimal amount of information about the past needed to
predict the future as well as possible. We suggest new theoretical
investigations into and new experiments on bacteria phenotypic bet-hedging in
fluctuating complex environments.


  We propose a mathematical model for collective sensing in a population
growing in a stochastically varying environment. In the population, individuals
use an information channel for sensing the environment, and two channels for
signal production and comprehension to communicate among themselves. We show
that existence of such system has a positive effect on population growth, hence
can have a positive evolutionary effect. We show that the gain in growth due to
the collective sensing is related to information theoretic entities, which can
be considered as the information content of this system from the environment.
We further show that heterogeneity in communication resulted from network or
spatial structure increases growth. We compute the growth rate of a population
residing on a lattice and show that growth rate near the maximum noise level in
observation or communication, increases exponentially as noise decreases. This
exponential effect makes the emergence of collective observation an easy
outcome in an evolutionary process. Furthermore, we are able to quantify
interesting effects such as accelerated growth, and simplification of decision
making due to information amplification by communication. Finally, we show that
an amount of noise in representation formation has more disadvantageous effect
compared to the same noise in signal production.


  Frequency-dependent selection reflects the interaction between different
species as they battle for limited resources in their environment. In a
stochastic evolutionary game the species relative fitnesses guides the
evolutionary dynamics which fluctuate due to random drift. Dependence of
species selection advantages on the environment introduces additional
possibilities for the evolutionary dynamics. We analyse a simple model in which
a random environment allows competing species to coexist for a long time before
a fixation of a single species happens. In our analysis we use stability in a
linear combination of competing species to approximate the stochastic dynamics
of the system by a diffusion on a one dimensional co-existence region. Our
method significantly simplifies calculating the probability of first extinction
and its expected time, and demonstrates a rigorous model reduction technique
for evaluating quasistationary properties of a stochastic evolutionary model.


  Chemotherapeutic response of cancer cells to a given compound is one of the
most fundamental information one requires to design anti-cancer drugs. Recent
advances in producing large drug screens against cancer cell lines provided an
opportunity to apply machine learning methods for this purpose. In addition to
cytotoxicity databases, considerable amount of drug-induced gene expression
data has also become publicly available. Following this, several methods that
exploit omics data were proposed to predict drug activity on cancer cells.
However, due to the complexity of cancer drug mechanisms, none of the existing
methods are perfect. One possible direction, therefore, is to combine the
strengths of both the methods and the databases for improved performance. We
demonstrate that integrating a large number of predictions by the proposed
method improves the performance for this task. The predictors in the ensemble
differ in several aspects such as the method itself, the number of tasks method
considers (multi-task vs. single-task) and the subset of data considered
(sub-sampling). We show that all these different aspects contribute to the
success of the final ensemble. In addition, we attempt to use the drug screen
data together with two novel signatures produced from the drug-induced gene
expression profiles of cancer cell lines. Finally, we evaluate the method
predictions by in vitro experiments in addition to the tests on data sets.The
predictions of the methods, the signatures and the software are available from
\url{http://mtan.etu.edu.tr/drug-response-prediction/}.


  Motivated by the problem of domain formation in chromosomes, we studied a
co--polymer model where only a subset of the monomers feel attractive
interactions. These monomers are displaced randomly from a regularly-spaced
pattern, thus introducing some quenched disorder in the system. Previous work
has shown that in the case of regularly-spaced interacting monomers this chain
can fold into structures characterized by multiple distinct domains of
consecutive segments. In each domain, attractive interactions are balanced by
the entropy cost of forming loops. We show by advanced replica-exchange
simulations that adding disorder in the position of the interacting monomers
further stabilizes these domains. The model suggests that the partitioning of
the chain into well-defined domains of consecutive monomers is a spontaneous
property of heteropolymers. In the case of chromosomes, evolution could have
acted on the spacing of interacting monomers to modulate in a simple way the
underlying domains for functional reasons.


  We propose a theoretical framework for the study of epidemics in structured
metapopulations, with heterogeneous agents, subjected to recurrent mobility
patterns. We propose to represent the heterogeneity in the composition of the
metapopulations as layers in a multiplex network, where nodes would correspond
to geographical areas and layers account for the mobility patterns of agents of
the same class. We analyze both the classical Susceptible-Infected-Susceptible
and the Susceptible-Infected-Removed epidemic models within this framework, and
compare macroscopic and microscopic indicators of the spreading process with
extensive Monte Carlo simulations. Our results are in excellent agreement with
the simulations. We also derive an exact expression of the epidemic threshold
on this general framework revealing a non-trivial dependence on the mobility
parameter. Finally, we use this new formalism to address the spread of diseases
in real cities, specifically in the city of Medellin, Colombia, whose
population is divided into six socio-economic classes, each one identified with
a layer in this multiplex formalism.


  Understanding the emergence of biological structures and their changes is a
complex problem. On a biochemical level, it is based on gene regulatory
networks (GRN) consisting on interactions between the genes responsible for
cell differentiation and coupled in a greater scale with external factors. In
this work we provide a systematic methodological framework to construct
Waddington's epigenetic landscape of the GRN involved in cellular determination
during the early stages of development of angiosperms. As a specific example we
consider the flower of the plant \textit{Arabidopsis thaliana}. Our model,
which is based on experimental data, recovers accurately the spatial
configuration of the flower during cell fate determination, not only for the
wild type, but for its homeotic mutants as well. The method developed in this
project is general enough to be used in the study of the relationship between
genotype-phenotype in other living organisms.


  Applying the method of moments to the chemical master equation (CME)
appearing in stochastic chemical kinetics often leads to the so-called closure
problem. Recently, several authors showed that this problem can be partially
overcome using moment-based semidefinite programs (SDPs). In particular, they
showed that moment-based SDPs can be used to calculate rigorous bounds on
various descriptions of the stochastic chemical kinetic system's stationary
distribution(s) -- for example, mean molecular counts, variances in these
counts, and so on. In this paper, we show that these ideas can be extended to
the corresponding dynamic problem, calculating time-varying bounds on the same
descriptions.