  The adaptive immune system is a natural diagnostic and therapeutic. It
recognizes threats earlier than clinical symptoms manifest and neutralizes
antigen with exquisite specificity. Recognition specificity and broad
reactivity is enabled via adaptive B- and T-cell receptors: the immune receptor
repertoire. The human immune system, however, is not omnipotent. Our natural
defense system sometimes loses the battle to parasites and microbes and even
turns against us in the case of cancer, autoimmune and inflammatory disease. A
long-standing dream of immunoengineers has been, therefore, to mechanistically
understand how the immune system sees, reacts and remembers antigens. Only very
recently, experimental and computational methods have achieved sufficient
quantitative resolution to start querying and engineering adaptive immunity
with great precision. In specific, these innovations have been applied with the
greatest fervency and success in immunotherapy, autoimmunity and vaccine
design. The work here highlights advances, challenges and future directions of
quantitative approaches which seek to advance the fundamental understanding of
immunological phenomena, and reverse engineer the immune system to produce
auspicious biopharmaceutical drugs and immunodiagnostics. Our review indicates
that the merger of fundamental immunology, computational immunology and
digital-biotechnology minimizes black box engineering, thereby advancing both
immunological knowledge and as well immunoengineering methodologies.


  Internal feedbacks are commonly present in biological populations and can
play a crucial role in the emergence of collective behavior. We consider a
generalization of Fisher-KPP equation to describe the temporal evolution of the
distribution of a single-species population. This equation includes the
elementary processes of random motion, reproduction and, importantly, nonlocal
interspecific competition, which introduces a spatial scale of interaction.
Furthermore, we take into account feedback mechanisms in diffusion and growth
processes, mimicked through density-dependencies controlled by exponents $\nu$
and $\mu$, respectively. These feedbacks include, for instance, anomalous
diffusion, reaction to overcrowding or to rarefaction of the population, as
well as Allee-like effects. We report that, depending on the dynamics in place,
the population can self-organize splitting into disconnected sub-populations,
in the absence of environment constraints. Through extensive numerical
simulations, we investigate the temporal evolution and stationary features of
the population distribution in the one-dimensional case. We discuss the crucial
role that density-dependency has on pattern formation, particularly on
fragmentation, which can bring important consequences to processes such as
epidemic spread and speciation.


  In growing plant cells, parallel ordering of microtubules (MTs) along the
inner surface of the cell membrane influences the direction of cell expansion
and thereby plant morphology. For correct expansion of organs that primarily
grow by elongating, such as roots and stems, MTs must bend in the
high-curvature direction along the cylindrically shaped cell membrane in order
to form the required circumferential arrays. Computational studies, which have
recapitulated the self-organization of these arrays, ignored MT mechanics and
assumed MTs follow geodesics of the cell surface. Here, we show, through
analysis of a derived Euler-Lagrange equation, that an elastic MT constrained
to a cylindrical surface will deflect away from geodesics and toward low
curvature directions to minimize bending energy. This occurs when the curvature
of the cell surface is relatively high for a given anchor density. In the limit
of infinite anchor density, MTs always follow geodesics. We compare our
analytical predictions to measured curvatures and anchor densities and find
that the regime in which cells are forming these cortical arrays straddles the
region of parameter space in which arrays must form under the antagonistic
influence of this mechanically induced deflection. Although this introduces a
potential obstacle to forming circumferentially orientated arrays that needs to
be accounted for in the models, it also raises the question of whether plants
use this mechanical phenomenon to regulate array orientation. The model also
constitutes an elegant generalization of the classical Euler-bucking
instability along with an intrinsic unfolding of the associated pitchfork
bifurcation.


  Using a coarse-grained model, self-organized assembly of proteins (e.g. CorA
and its inner segment iCorA) is studied by examining quantities such as contact
profile, radius of gyration, and structure factor as a function of protein
concentration at a range of low (native phase) to high (denature phase)
temperatures. Visual inspections show distinct structures, i.e. isolated
globular bundles to entangled network on multiple length scales in dilute to
crowded protein concentrations. In native phase, the radius of gyration of the
protein does not vary much with the protein concentration while that of its
inner segment increases systematically. In contrast, the radius of gyration of
the protein shows enormous growth with the concentration due to entanglement
while that of the inner segment remains almost constant in denatured phase. The
multi-scale morphology of the collective assembly is quantified by estimating
the effective dimension D of protein from scaling of the structure factor:
collective assembly from inner segments remains globular (D aroud 3) at almost
all length scales in its native phase while that from protein chains shows
sparsely distributed morphology with D around 2 in entire temperature range due
to entanglement except in crowded environment at low temperature where D around
2.6. Higher morphological response of chains with only the inner-segments due
to selective interactions in its native phase may be more conducive to
self-organizing mechanism than that of the remaining segments of the protein
chains.


  By combining analytical and numerical calculations, we investigate the
minimal-energy shape of short DNA loops of approximately $100$ base pairs (bp).
We show that in these loops the excess twist density oscillates as a response
to an imposed bending stress, as recently found in DNA minicircles and observed
in nucleosomal DNA. These twist oscillations, here referred to as twist waves,
are due to the coupling between twist and bending deformations, which in turn
originates from the asymmetry between DNA major and minor grooves. We introduce
a simple analytical variational shape, that reproduces the exact loop energy up
to the fourth significant digit, and is in very good agreement with shapes
obtained from coarse-grained simulations. We, finally, analyze the loop
dynamics at room temperature, and show that the twist waves are robust against
thermal fluctuations. They perform a normal diffusive motion, whose origin is
briefly discussed.


  Here, we show by solution nuclear magnetic resonance measurements that the
urea-unfolded protein apoazurin becomes elongated when the synthetic crowding
agent dextran 20 is present, in contrast to the prediction from the
macromolecular crowding effect based on the argument of volume exclusion. To
explore the complex interactions beyond volume exclusion, we employed
coarse-grained molecular dynamics simulations to explore the conformational
ensemble of apoazurin in a box of monodisperse crowders under strong chemically
denaturing conditions. The elongated conformation of unfolded apoazurin appears
to result from the interplay of the effective attraction between the protein
and crowders and the shape of the crowders. With a volume-conserving crowder
model, we show that the crowder shape provides an anisotropic direction of the
depletion force, in which a bundle of surrounding rod-like crowders stabilize
an elongated conformation of unfolded apoazurin in the presence of effective
attraction between the protein and crowders.


  Influenza epidemics result in a public health and economic burden around the
globe. Traditional surveillance techniques, which rely on doctor visits,
provide data with a delay of 1-2 weeks. A means of obtaining real-time data and
forecasting future outbreaks is desirable to provide more timely responses to
influenza epidemics. In this work, we present the first implementation of a
novel data set by demonstrating its ability to supplement traditional disease
surveillance at multiple spatial resolutions. We use Internet traffic data from
the Centers for Disease Control and Prevention (CDC) website to determine the
potential usability of this data source. We test the traffic generated by ten
influenza-related pages in eight states and nine census divisions within the
United States and compare it against clinical surveillance data. Our results
yield $r^2$ = 0.955 in the most successful case, promising results for some
cases, and unsuccessful results for other cases. These results demonstrate that
Internet data may be able to complement traditional influenza surveillance in
some cases but not in others. Specifically, our results show that the CDC
website traffic may inform national and division-level models but not models
for each individual state. In addition, our results show better agreement when
the data were broken up by seasons instead of aggregated over several years. In
the interest of scientific transparency to further the understanding of when
Internet data streams are an appropriate supplemental data source, we also
include negative results (i.e., unsuccessful models). We anticipate that this
work will lead to more complex nowcasting and forecasting models using this
data stream.


  We modified the Lotka-Volterra Equations with the assumption that two of the
original four constant parameters in the traditional equations are
time-dependent. In the first place, we assumed that the human population
(borrowed from the T-Function) plays the role as the prey while all lethal
factors that jeopardize the existence of the human race as the predator.
Although we could still calculate the time-dependent lethal function, the idea
of treating the lethal factors as the prey was too general to recognize the
meaning of them. Hence, in the second part of the modified Lotka-Volterra
Equations, we exchanged the roles between the prey and the predator. This time,
we treated the prey as the natural resources while the predator as the human
population (still borrowed from the T-Function). After carefully choosing
appropriate parameters to match the maximum carrying capacity with the
saturated number of the human population predicted by the T-Function, we
successfully calculated the natural resources as a function of time. Contrary
to our intuition, the carrying capacity is constant over time rather than a
time-varying function, with the constant value of 10.2 billion people.


  Cryptography is the science that secures data and communication over the
network by applying mathematics and logic to design strong encryption methods.
In the modern era of e-business and e-commerce the protection of
confidentiality, integrity and availability (CIA triad) of stored information
as well as of transmitted data is very crucial. Deoxyribonucleic acid (DNA) is
a genetic molecule consisting of two linked strands that wind around each other
to form a double helical structure. The backbone of each strand is made of
alternating deoxyribose sugar and phosphate groups. To each sugar one of four
bases are attached i.e., adenine (A), cytosine (C), guanine (G) and thymine
(T). DNA molecules, having the capacity to store, process and transmit
information, inspires the idea of DNA cryptography. It is the rapid emerging
unconventional techniques which combines the chemical characteristics of
biological DNA sequences with classical cryptography to ensure non-vulnerable
transmission of data. This innovative method is based on the notion of DNA
computing. The methodologies of DNA cryptography are not coded mathematically;
thus, it could be too secure to be cracked easily.


  The adaptive immune system of vertebrates can detect, respond to, and
memorize diverse pathogens from past experience. While the clonal selection of
T helper (Th) cells is the simple and established mechanism to better recognize
new pathogens, the question that still remains unexplored is how the Th cells
can acquire better ways to bias the responses of immune cells for eliminating
pathogens more efficiently by translating the recognized antigen information
into regulatory signals. In this work, we address this problem by associating
the adaptive immune network organized by the Th cells with reinforcement
learning (RL). By employing recent advancements of network-based RL, we show
that the Th immune network can acquire the association between antigen patterns
of and the effective responses to pathogens. Moreover, the clonal selection as
well as other inter-cellular interactions are derived as a learning rule of the
network. We also demonstrate that the stationary clone-size distribution after
learning shares characteristic features with those observed experimentally. Our
theoretical framework may contribute to revising and renewing our understanding
of adaptive immunity as a learning system.


  Techniques for leaf-cutting ant control have been investigated in literature
due to the importance of the damage they cause to agriculture. The
effectiveness of different forms of control is explored in researches aimed at
identifying the balance between pest control efficiency and environmental
damage caused by the treatments applied. Plantations with large territorial
extensions, which can be contiguous or not, are usually subdivided into local
administration that collects data to determine the frequencies of areas size
occupied by ant nests. The purpose of this work is to build a relationship of
similarities among different geographical regions using the frequency data of
nests size occurrence by applying Information Bottleneck (IB) method and
Principal Component Analysis (PCA). IB allows simultaneous clustering of each
region with the ant nest size distribution, while PCA is used to reduce the
variable dimensionalities into a three-dimensional representation of the
results. The approach was applied to data of leaf-cutting ants Atta spp.
(Hymenoptera: Formicidae) in cultivated Eucalyptus spp forests in Sao Paulo,
Brazil. The results suggest information acquired by the method can help
coordinate pest management, such as the allocation of baits, material and
personnel.


  The bacterial flagellar motor is an ion-powered transmembrane protein complex
which drives swimming in many bacterial species. The motor consists of a
cytoplasmic 'rotor' ring and a number of 'stator' units, which are bound to the
cell wall of the bacterium. Recently, it has been shown that the number of
functional torque-generating stator units in the motor depends on the external
load, and suggested that mechanosensing in the flagellar motor is driven via a
'catch bond' mechanism in the motor's stator units. We present a method that
allows us to measure -- on a single motor -- stator unit dynamics across a
large range of external loads, including near the zero-torque limit. By
attaching superparamagnetic beads to the flagellar hook, we can control the
motor's speed via a rotating magnetic field. We manipulate the motor to four
different speed levels in two different ion-motive force (IMF) conditions. This
framework allows for a deeper exploration into the mechanism behind
load-dependent remodelling by separating out motor properties, such as rotation
speed and energy availability in the form of IMF, that affect the motor torque.


  Mass spectrometry (MS) is used widely in biomolecular structural analysis and
is particularly dominant in the study of proteins. Despite its considerable
power, state-of-the-art protein MS frequently suffers from limited reliability
of spectrum-to-structure assignments. This could not be solved fully by the
dramatic increase in mass accuracy and resolution of modern MS instrumentation
or by the introduction of new fragmentation methods. Here we present a new kind
of two-dimensional mass spectrometry for high fidelity determination of a
biomolecular primary structure based on partial covariance mapping. Partial
covariance two-dimensional mass spectrometry (pC-2DMS) detects intrinsic
statistical correlations between biomolecular fragments originating from the
same or consecutive decomposition events. This enables identification of pairs
of ions produced along the same fragmentation pathway of a biomolecule across
its entire fragment mass spectrum. We demonstrate that the fragment-fragment
correlations revealed by pC-2DMS provide much more specific information on the
amino acid sequence and its covalent modifications than the individual fragment
mass-to-charge ratios on which standard one-dimensional MS is based. We
illustrate the power of pC-2DMS by using it to resolve structural isomers of
combinatorially modified histone peptides inaccessible to standard MS.


  This paper introduces a new way to define a genome rearrangement distance,
using the concept of mean first passage time from probability theory.
Crucially, this distance estimate provides a genuine metric on genome space. We
develop the theory and introduce a link to a graph-based zeta function. The
approach is very general and can be applied to a wide variety of
group-theoretic models of genome evolution.


  Genome editing allows scientists to change an organism's DNA. One promising
genome editing protocol, already validated in living organisms, is based on
clustered regularly interspaced short palindromic repeats (CRISPR)/Cas
protein-nucleic acid complexes. When the CRISPR/Cas approach was first
demonstrated in 2012, its advantages with respect to previously available
techniques, such as zinc-finger nucleases (ZFNs) and transcription
activator-like effector nucleases (TALENs), immediately got attention and the
method has seen a surge of experimental and computational investigations since
then. However, the molecular mechanisms involved in target DNA recognition and
cleavage are still not completely resolved and need further attention. The
large size and complex nature of CRISPR/Cas9 complexes has been a challenge for
computational studies, but some seed results exist and are illuminating on the
cleavage activity. In this short review, we present recent progress in studying
CRISPR/Cas9 systems by molecular dynamics simulations with coarse-grained and
atomistic descriptions, including enhanced sampling.


  The classical fold bifurcation is a paradigmatic example of a critical
transition. It has been used in a variety of contexts, including in particular
ecology and climate science, to motivate the role of slow recovery rates and
increased autocorrelations as early-warning signals of such transitions.
  We study the influence of external forcing on fold bifurcations and the
respective early-warning signals. Thereby, our prime examples are
single-species population dynamical models with Allee effect under the
influence of either quasiperiodic forcing or bounded random noise. We show that
the presence of these external factors may lead to so-called non-smooth fold
bifurcations, and thereby has a significant impact on the behaviour of the
Lyapunov exponents (and hence the recovery rates). In particular, it may lead
to the absence of critical slowing down prior to population collapse. More
precisely, unlike in the unforced case, the question whether slow recovery
rates can be observed or detected prior to the transition crucially depends on
the chosen time-scales and the size of the considered data set.


  Structured RNA plays many functionally relevant roles in molecular life.
Structural information, while required to understand the functional cycles in
detail, is challenging to gather. Computational methods promise to complement
experimental efforts by predicting three-dimensional RNA models. Here, we
provide a concise view of the state of the art methodologies with a focus on
the strengths and the weaknesses of the different approaches. Furthermore, we
analyzed the recent developments regarding the use of coevolutionary
information and how it can boost the prediction performances. We finally
discuss some open perspectives and challenges for the near future in the RNA
structural stability field.


  In recent years machine learning (ML) took bio- and cheminformatics fields by
storm, providing new solutions for a vast repertoire of problems related to
protein sequence, structure, and interactions analysis. ML techniques, deep
neural networks especially, were proven more effective than classical models
for tasks like predicting binding affinity for molecular complex. In this work
we investigated the earlier stage of drug discovery process - finding druggable
pockets on protein surface, that can be later used to design active molecules.
For this purpose we developed a 3D fully convolutional neural network capable
of binding site segmentation. Our solution has high prediction accuracy and
provides intuitive representations of the results, which makes it easy to
incorporate into drug discovery projects. The model's source code, together
with scripts for most common use-cases is freely available at
http://gitlab.com/cheminfIBB/kalasanty


  Fluctuations in intracellular reactions (intrinsic noise) reduce the
information transmitted from an extracellular input to a cellular response.
However, recent studies have demonstrated that the decrease in the transmitted
information with respect to extracellular input fluctuations (extrinsic noise)
is smaller when the intrinsic noise is larger. Therefore, it has been suggested
that robustness against extrinsic noise increases with the level of the
intrinsic noise. We call this phenomenon intrinsic noise-induced robustness
(INIR). As previous studies on this phenomenon have focused on complex
biochemical reactions, the relation between INIR and the input--output of a
system is unclear. Moreover, the mechanism of INIR remains elusive. In this
paper, we address these questions by analyzing simple models. We first analyze
a model in which the input--output relation is linear. We show that the
robustness against extrinsic noise increases with the intrinsic noise,
confirming the INIR phenomenon. Moreover, the robustness against the extrinsic
noise is more strongly dependent on the intrinsic noise when the variance of
the intrinsic noise is larger than that of the input distribution. Next, we
analyze a threshold model in which the output depends on whether the input
exceeds the threshold. When the threshold is equal to the mean of the input,
INIR is realized, but when the threshold is much larger than the mean, the
threshold model exhibits stochastic resonance, and INIR is not always apparent.
The robustness against extrinsic noise and the transmitted information can be
traded off against one another in the linear model and the threshold model
without stochastic resonance, whereas they can be simultaneously increased in
the threshold model with stochastic resonance.


  Despite the wealth of empirical and theoretical studies, the origin and
maintenance of cooperation is still an evolutionary riddle. In this context,
ecological life-history traits which affect the efficiency of selection may
play a role, though these are often ignored. We consider here species such as
bacteria, fungi, invertebrates and plants which exhibit resting stages in the
form of a quiescent state or a seedbank. When quiescent, individuals are
inactive and reproduce upon activation, while under seed bank parents produce
offspring remaining dormant for different amount of time. We assume weak
frequency-dependent selection modeled using game-theory and the prisoners
dilemma (cooperation/defect) as payoff matrix. The cooperators and defectors
are allowed to evolve different quiescence or dormancy times. By means of
singular perturbation theory we reduce the model to a one-dimensional equation
resembling the well known replicator equation, where the gain functions are
scaled with lumped parameters reflecting the time scale of the resting state of
the cooperators and defectors. If both time scales are identical cooperation
cannot persist in a homogeneous population. If, however, the time scale of the
cooperator is distinctively different from that of the defector, cooperation
may become a locally asymptotically stable strategy. Interestingly enough, in
the seedbank case the cooperator needs to be faster than the defector, while in
the quiescent case the cooperator has to be slower. We use adaptive dynamics to
identify situations where cooperation may evolve and form a convergent stable
ESS. We conclude by highlighting the relevance fo these results for many
non-model species and the maintenance of cooperation in microbial, invertebrate
or plant populations.


  We study a class of the stochastic May-Leonard models, with three species
dominating each other in a cyclic nonhierarchical way, according to the
rock-paper-scissors game. We introduce an unevenness in the system, by
considering that one of the species is weaker because of a lower selection
probability. The simulation results show that the pattern formation is
drastically affected by the presence of the weaker species, with no spiral
waves arising immediately from random initial conditions. Instead,
single-species spatial domains cyclically dominate the entire territory until a
region occupied by the weaker species is sufficiently narrow to be crossed by
individuals without being selected. This leads to the appearance of spatial
patterns responsible for the species coexistence. We verify that the asymmetry
in the selection probabilities leads to different spatial autocorrelation
function and average relative species abundances. Finally, we investigate the
coexistence probability and show that the surviving species depends on the
level of unevenness of the model and the mobility of individuals.


  Sequencing large number of candidate disease genes which cause diseases in
order to identify the relationship between them is an expensive and
time-consuming task. To handle these challenges, different computational
approaches have been developed. Based on the observation that genes associated
with similar diseases have a higher likelihood of interaction, a large class of
these approaches relay on analyzing the topological properties of biological
networks. However, the incomplete and noisy nature of biological networks is
known as an important challenge in these approaches. In this paper, we propose
a two-step framework for disease gene prioritization: (1) construction of a
reliable human FLN using sequence information and machine learning techniques,
(2) prioritizing the disease gene relations based on the constructed FLN. On
our framework, unlike other FLN based frameworks that using FLNs based on
integration of various low quality biological data, the sequence of proteins is
used as the comprehensive data to construct a reliable initial network. In
addition, the physicochemical properties of amino-acids are employed to
describe the functionality of proteins. All in all, the proposed approach is
evaluated and the results indicate the high efficiency and validity of the FLN
in disease gene prioritization.


  Systems of random linear equations may or may not have solutions with all
components being non-negative. The question is, e.g., of relevance when the
unknowns are concentrations or population sizes. In the present paper we show
that if such systems are large the transition between these two possibilities
occurs at a sharp value of the ratio between the number of unknowns and the
number of equations. We analytically determine this threshold as a function of
the statistical properties of the random parameters and show its agreement with
numerical simulations. We also make contact with two special cases that have
been studied before: the storage problem of a perceptron and the resource
competition model of MacArthur.


  Species inhabiting urban environments experience enormous anthropogenic
stress. Behavioural plasticity and flexibility of temperament are crucial to
successful urban-adaptation. Urban free-ranging dogs experience variable human
impact, from positive to negative and represent an ideal system to evaluate the
effects of human-induced stress on behavioural plasticity. We tested 600 adult
dogs from 60 sites across India, categorised as high - HF, low - LF, and
intermediate - IF human flux zones, to understand their sociability towards an
unfamiliar human. Dogs in the HF and IF zones were bolder and as compared to
their shy counterparts in LF zones. The IF zone dogs were the most sociable.
This is the first-ever study aimed to understand how the experiences of
interactions with humans in its immediate environment might shape the responses
of an animal to humans. This is very relevant in the context of human-animal
conflict induced by rapid urbanization and habitat loss across the world.


  Cryo-electron microscopy (EM) single particle reconstruction is an entirely
general technique for 3D structure determination of macromolecular complexes.
However, because the images are taken at low electron dose, it is extremely
hard to visualize the individual particle with low contrast and high noise
level. In this paper, we propose a novel approach called multi-frequency vector
diffusion maps (MFVDM) to improve the efficiency and accuracy of cryo-EM 2D
image classification and denoising. This framework incorporates different
irreducible representations of the estimated alignment between similar images.
In addition, we propose a graph filtering scheme to denoise the images using
the eigenvalues and eigenvectors of the MFVDM matrices. Through both simulated
and publicly available real data, we demonstrate that our proposed method is
efficient and robust to noise compared with the state-of-the-art cryo-EM 2D
class averaging and image restoration algorithms.


  Cancer cells can acquire a spectrum of stable hybrid epithelial/mesenchymal
(E/M) states during epithelial-mesenchymal transition (EMT). Cells in these
hybrid E/M phenotypes often combine epithelial and mesenchymal features and
tend to migrate collectively commonly as small clusters. Such collectively
migrating cancer cells play a pivotal role in seeding metastases and their
presence in cancer patients indicates an adverse prognostic factor. Moreover,
cancer cells in hybrid E/M phenotypes tend to be more associated with stemness
which endows them with tumor-initiation ability and therapy resistance. Most
recently, cells undergoing EMT have been shown to promote immune suppression
for better survival. A systematic understanding of the emergence of hybrid E/M
phenotypes and the connection of EMT with stemness and immune suppression would
contribute to more effective therapeutic strategies. In this review, we first
discuss recent efforts combining theoretical and experimental approaches to
elucidate mechanisms underlying EMT multi-stability (i.e. the existence of
multiple stable phenotypes during EMT) and the properties of hybrid E/M
phenotypes. Following we discuss non-cell-autonomous regulation of EMT by cell
cooperation and extracellular matrix. Afterwards, we discuss various metrics
that can be used to quantify EMT spectrum. We further describe possible
mechanisms underlying the formation of clusters of circulating tumor cells.
Last but not least, we summarize recent systems biology analysis of the role of
EMT in the acquisition of stemness and immune suppression.


  A great variety of molecular components is encapsulated in cells. Each of
these components is replicated for cell reproduction. To address an essential
role of the huge diversity of cellular components, we study a model of
protocells that convert resources into catalysts with the aid of a catalytic
reaction network. As the resources are limited, it is shown that diversity in
intracellular components is increased to allow the use of diverse resources for
cellular growth. Scaling relation is demonstrated between resource abundances
and molecular diversity. We then study how the molecule species diversify and
complex catalytic reaction networks develop through the evolutionary course. It
is shown that molecule species first appear, at some generations, as parasitic
ones that do not contribute to replication of other molecules. Later, the
species turn to be host species that support the replication of other species.
With this successive increase of host species, a complex joint network evolves.
The present study sheds new light on the origin of molecular diversity and
complex reaction network at the primitive stage of a cell.


  Many of the existing stochastic models of gene expression contain the
first-order decay reaction term that may describe active protein degradation or
dilution. If the model variable is interpreted as the molecule number, and not
concentration, the decay term may also approximate the loss of protein
molecules due to cell division as a continuous degradation process. The seminal
model of that kind leads to gamma distributions of protein levels, whose
parameters are defined by the mean frequency of protein bursts and mean burst
size. However, such models (whether interpreted in terms of molecule numbers or
concentrations) do not correctly account for the noise due to protein
partitioning between daughter cells.
  We present an exactly solvable stochastic model of gene expression in cells
dividing at random times, where we assume description in terms of molecule
numbers with a constant mean protein burst size. The model is based on a
population balance equation supplemented with protein production in random
bursts.
  If protein molecules are partitioned equally between daughter cells, we
obtain at steady state the analytical expressions for probability distributions
similar in shape to gamma distribution, yet with quite different values of mean
burst size and mean burst frequency than would result from fitting of the
classical continuous-decay model to these distributions. For random
partitioning of protein molecules between daughter cells, we obtain the moment
equations for the protein number distribution and thus the analytical formulas
for the squared coefficient of variation.


  Personalized cancer vaccines are envisioned as the next generation rational
cancer immunotherapy. The key step in developing personalized therapeutic
cancer vaccines is to identify tumor-specific neoantigens that are on the
surface of tumor cells. A promising method for this is through de novo peptide
sequencing from mass spectrometry data. In this paper we introduce DeepNovoV2,
the state-of-the-art model for peptide sequencing. In DeepNovoV2, a spectrum is
directly represented as a set of (m/z, intensity) pairs, therefore it does not
suffer from the accuracy-speed/memory trade-off problem. The model combines an
order invariant network structure (T-Net) and recurrent neural networks and
provides a complete end-to-end training and prediction framework to sequence
patterns of peptides. Our experiments on a wide variety of data from different
species show that DeepNovoV2 outperforms previous state-of-the-art methods,
achieving 13.01-23.95\% higher accuracy at the peptide level.


  Solutions of manually purified gastric mucins have been shown to be promising
lubricants for biomedical purposes, where they can efficiently reduce friction
and wear. However, so far, such mucin solutions have been mostly tested in
specific settings, and variations in the composition of the lubricating fluid
have not been systematically explored. We here fill this gap and determine the
viscosity, adsorption behavior, and lubricity of porcine gastric mucin
solutions on hydrophobic surfaces at different pH levels, mucin and salt
concentrations and in the presence of other proteins. We demonstrate that mucin
solutions provide excellent lubricity even at very low concentrations of 0.01 %
(w/v), over a broad range of pH levels and even at elevated ionic strength.
Furthermore, we provide mechanistic insights into mucin lubricity, which help
explain how certain variations in physiologically relevant parameters can limit
the lubricating potential of mucin solutions. Our results motivate that
solutions of manually purified mucin solutions can be powerful biomedical
lubricants, e.g. serving as eye drops, mouth sprays or as a personal lubricant
for intercourse.


  One of the simplest examples of a robust heteroclinic cycle involves three
saddle equilibria: each one is unstable to the next in turn, and connections
from one to the next occur within invariant subspaces. Such a situation can be
described by a third-order ordinary differential equation (ODE), and typical
trajectories approach each equilibrium point in turn, spending progressively
longer to cycle around the three points but never stopping. This cycle has been
invoked as a model of cyclic competition between populations adopting three
strategies, characterised as Rock, Paper and Scissors. When spatial
distribution and mobility of the populations is taken into account, waves of
Rock can invade regions of Scissors, only to be invaded by Paper in turn. The
dynamics is described by a set of partial differential equations (PDEs) that
has travelling wave (in one dimension) and spiral (in two dimensions)
solutions. In this paper, we explore how the robust heteroclinic cycle in the
ODE manifests itself in the PDEs. Taking the wavespeed as a parameter, and
moving into a travelling frame, the PDEs reduce to a sixth-order set of ODEs,
in which travelling waves are created in a Hopf bifurcation and are destroyed
in three different heteroclinic bifurcations, depending on parameters, as the
travelling wave approaches the heteroclinic cycle. We explore the three
different heteroclinic bifurcations, none of which have been observed in the
context of robust heteroclinic cycles previously. These results are an
important step towards a full understanding of the spiral patterns found in two
dimensions, with possible application to travelling waves and spirals in other
population dynamics models.


  Identifying which taxa in our microbiota are associated with traits of
interest is important for advancing science and health. However, the
identification is challenging because the measured vector of taxa counts (by
amplicon sequencing) is compositional, so a change in the abundance of one
taxon in the microbiota induces a change in the number of sequenced counts
across all taxa. The data is typically sparse, with zero counts present either
due to biological variance or limited sequencing depth (technical zeros). For
low abundance taxa, the chance for technical zeros is non-negligible. We show
that existing methods designed to identify differential abundance for
compositional data may have an inflated number of false positives due to
improper handling of the zero counts. We introduce a novel non-parametric
approach which provides valid inference even when the fraction of zero counts
is substantial. Our approach uses a set of reference taxa that are
non-differentially abundant, which can be estimated from the data or from
outside information. We show the usefulness of our approach via simulations, as
well as on three different data sets: a Crohn's disease study, the Human
Microbiome Project, and an experiment with 'spiked-in' bacteria.


  It has been known for nearly a decade that deterministically modeled reaction
networks that are weakly reversible and consist of a single linkage class have
trajectories that are bounded from both above and below by positive constants
(so long as the initial condition has strictly positive components). It is
conjectured that the stochastically modeled analogs of these systems are
positive recurrent. We prove this conjecture in the affirmative under the
following additional assumptions: (i) the system is binary and (ii) for each
species, there is a complex (vertex in the associated reaction diagram) that is
a multiple of that species. To show this result, a new proof technique is
developed in which we study the recurrence properties of the $n-$step embedded
discrete time Markov chain.


  In the drug repurposing approach, the chemically diverse and potentially safe
molecules can be explored as therapeutic potential out of those originally
targeted indications. However, the intellectual property rights, and
competitive over-heating issue from sharing the same drug pipeline information
are drawback of the approach. In current paradigm of drug discovery, despite
approval of new chemical entity, the finding of a new ring system for drug
discovery have stagnated. In addition, there are still rarely used ring systems
from FDA approved drug that can be optimized and used for repurposing purpose.
In the present work, these rarely used drug ring system were re-purposed with
our bidirectional target and ring screening (TR- dual screening) approach
within our defined objectives. The TR-dual screening suggested RARbeta and
cyproheptadine as a pair of target and ring system. The selected ring system
was virtually grown with predefined criterion, ranking of molecular docking and
top score poses made us select test compounds for synthesis and biological
proof. Among the tested compounds, 6c proved RARbeta agonist activity. The
concept of TR dual optimization can contribute to picking either a compound
having acceptable ADMET profile or a target for a selected drug scaffold.


  Natural microbial communities contain hundreds to thousands of interacting
species. For this reason, computational simulations are playing an increasingly
important role in microbial ecology. In this manuscript, we present a new
open-source, freely available Python package called Community Simulator for
simulating microbial population dynamics in a reproducible, transparent and
scalable way. The Community Simulator includes five major elements: tools for
preparing the initial states and environmental conditions for a set of samples,
automatic generation of dynamical equations based on a dictionary of modeling
assumptions, random parameter sampling with tunable levels of metabolic and
taxonomic structure, parallel integration of the dynamical equations, and
support for metacommunity dynamics with migration between samples. To
significantly speed up simulations using Community Simulator, our Python
package implements a new Expectation-Maximization (EM) algorithm for finding
equilibrium states of community dynamics that exploits a recently discovered
duality between ecological dynamics and convex optimization. We present data
showing that this EM algorithm improves performance by between one and two
orders compared to direct numerical integration of the corresponding ordinary
differential equations. We conclude by listing several recent applications of
the Community Simulator to problems in microbial ecology, and discussing
possible extensions of the package for directly analyzing microbiome
compositional data.


  Hi-C experiments are used to infer the contact probabilities between loci
separated by varying genome lengths. Contact probability should decrease as the
spatial distance between two loci increases. However, studies comparing Hi-C
and FISH data show that in some cases the distance between one pair of loci,
with larger Hi-C readout, is paradoxically larger compared to another pair with
a smaller value of the contact probability. Here, we show that the FISH-Hi-C
paradox can be resolved using a theory based on a Generalized Rouse Model for
Chromosomes (GRMC). The FISH-Hi-C paradox arises because the cell population is
highly heterogeneous, which means that a given contact is present in only a
fraction of cells. Insights from the GRMC is used to construct a theory,
without any adjustable parameters, to extract the distribution of
subpopulations from the FISH data, which quantitatively reproduces the Hi-C
data. Our results show that heterogeneity is pervasive in genome organization
at all length scales, reflecting large cell-to-cell variations.


  Machine learning provides a broad framework for addressing high-dimensional
prediction problems in classification and regression. While machine learning is
often applied for imaging problems in medical physics, there are many efforts
to apply these principles to biological data towards questions of radiation
biology. Here, we provide a review of radiogenomics modeling frameworks and
efforts towards genomically-guided radiotherapy. We first discuss medical
oncology efforts to develop precision biomarkers. We next discuss similar
efforts to create clinical assays for normal tissue or tumor radiosensitivity.
We then discuss modeling frameworks for radiosensitivity and the evolution of
machine learning to create predictive models for radiogenomics.


  Measures of tree balance play an important role in various research areas,
for example in phylogenetics. There they are for instance used to test whether
an observed phylogenetic tree differs significantly from a tree generated by
the Yule model of speciation. One of the most popular indices in this regard is
the Colless index, which measures the degree of balance for rooted binary
trees. While many statistical properties of the Colless index (e.g. asymptotic
results for its mean and variance under different models of speciation) have
already been discussed in different contexts, we focus on its extremal
properties. While it is relatively straightforward to characterize trees with
maximal Colless index, the analysis of the minimal value of the Colless index
and the characterization of trees that achieve it, are much more involved. In
this note, we therefore focus on the minimal value of the Colless index for any
given number of leaves. We derive both a recursive formula for this minimal
value, as well as an explicit expression, which shows a surprising connection
between the Colless index and the so-called Blancmange curve, a fractal curve
that is also known as the Takagi curve. Moreover, we characterize two classes
of trees that have minimal Colless index, consisting of the set of so-called
\emph{maximally balanced trees} and a class of trees that we call \emph{greedy
from the bottom trees}. Furthermore, we derive an upper bound for the number of
trees with minimal Colless index by relating these trees with trees with
minimal Sackin index (another well-studied index of tree balance).


  A metapopulations network is a multi-patch habitat system, where populations
live and interact in the habitat patches, and individuals disperse from one
patch to the other via dispersal connections. The loss of dispersal connections
among the habitat patches can impact the stability of the system. In this work,
we determine if there exist(s) set(s) of dispersal connections removal of which
causes partitioning(s) of the metapopulations network into dynamically stable
sub-networks. Our study finds that there exists a lower bound threshold Fiedler
value which guarantees the dynamical stability of the network dynamics.
Necessary and sufficient mathematical conditions for finding partitions that
result in sub-networks with the desired threshold Fiedler values have been
derived and illustrated with examples. Although posed and discussed in the
ecological context, it may be pointed out that such partitioning problems exist
across any spatially discrete but connected dynamical systems with
reaction-diffusion. Non-ecological examples are power distribution grids,
intra-cellular reaction pathway networks and high density nano-fluidic
lab-on-chip applications.


  In this paper, we study the number of equilibria of the replicator-mutator
dynamics for both deterministic and random multi-player two-strategy
evolutionary games. For deterministic games, using Decartes' rule of signs, we
provide a formula to compute the number of equilibria in multi-player games via
the number of change of signs in the coefficients of a polynomial. For
two-player social dilemmas (namely, the Prisoner's Dilemma, Snowdrift, Stag
Hunt, and Harmony), we characterize (stable) equilibrium points and
analytically calculate the probability of having a certain number of equilibria
when the payoff entries are uniformly distributed. For multi-player random
games whose payoffs are independently distributed according to a normal
distribution, by employing techniques from random polynomial theory, we compute
the expected or average number of internal equilibria. In addition, we perform
extensive simulations by sampling and averaging over a large number of possible
payoff matrices to compare with and illustrate analytical results. Numerical
simulations also suggest several interesting behaviour of the average number of
equilibria when the number of players is sufficiently large or when the
mutation is sufficiently small. In general, we observe that introducing
mutation results in a larger average number of internal equilibria than when
mutation is absent, implying that mutation leads to larger behavioural
diversity in dynamical systems. Interestingly, this number is largest when
mutation is rare rather than when it is frequent.


  A stochastic reaction-diffusion model is studied on a networked support. In
each patch of the network two species are assumed to interact following a
non-normal reaction scheme. When the interaction unit is replicated on a
directed linear lattice, noise gets amplified via a self-consistent process
which we trace back to the degenerate spectrum of the embedding support. The
same phenomenon holds when the system is bound to explore a quasi degenerate
network. In this case, the eigenvalues of the Laplacian operator, which governs
species diffusion, accumulate over a limited portion of the complex plane. The
larger the network, the more pronounced the amplification. Beyond a critical
network size, a system deemed deterministically stable, hence resilient, may
turn unstable, yielding seemingly regular patterns in the concentration amount.
Non-normality and quasi-degenerate networks may therefore amplify the inherent
stochasticity, and so contribute to altering the perception of resilience, as
quantified via conventional deterministic methods.


  Many questions that we have about the history and dynamics of organisms have
a geographical component: How many are there, and where do they live? How do
they move and interbreed across the landscape? How were they moving a thousand
years ago, and where were the ancestors of a particular individual alive today?
Answers to these questions can have profound consequences for our understanding
of history, ecology, and the evolutionary process. In this review, we discuss
how geographic aspects of the distribution, movement, and reproduction of
organisms are reflected in their pedigree across space and time. Because the
structure of the pedigree is what determines patterns of relatedness in modern
genetic variation, our aim is to thus provide intuition for how these processes
leave an imprint in genetic data. We also highlight some current methods and
gaps in the statistical toolbox of spatial population genetics.


  It has been a puzzling question why some organisms reproduce sexually. Fisher
and Muller hypothesized that reproducing by sex can speed up the evolution.
They explained that in the sexual reproduction, recombination can combine
beneficial alleles that lie on different chromosomes, which speeds up the time
that those beneficial alleles spread to the entire population. We consider a
population model of fixed size $N$, in which we will focus on two loci on a
chromosome. Each allele at each locus can mutate into a beneficial allele at
rate $\mu_N$. The individuals with 0, 1, and 2 beneficial alleles die at rates
$1, 1-s_N$ and $1-2s_N$ respectively. When an individual dies, with probability
$1-r_N$, the new individual inherits both alleles from one parent, chosen at
random from the population, while with probability $r_N$, recombination occurs,
and the new individual receives its two alleles from different parents. Under
certain assumptions on the parameters $N, \mu_N, s_N$ and $r_N$, we obtain an
asymptotic approximation for the time that both beneficial alleles spread to
the entire population. When the recombination probability is small, we show
that recombination does not speed up the time that the two beneficial alleles
spread to the entire population, while when the recombination probability is
large, we show that recombination decreases the time, which agrees with
Fisher-Muller hypothesis, and confirms the advantage of reproducing by sex.


  Matsushita et al (1992) have done an interesting finding. They observed that
the shape parameter of the Weibull model presented systematic changes over time
and age when applied to mortality data for males and females in Japan. They
have also estimated that this parameter was smaller in the 1891-1898 data in
Japan compared to the 1980 mortality data and they presented an illustrative
figure for females where the values of the shape parameter are illustrated on
the diagram close to the corresponding survival curves. However, they have not
provided an analytical explanation of this behavior of the shape parameter of
the Weibull model. The cumulative hazard of this model can express the additive
process of applying a force in a material for enough time before cracking. To
pass to the human data, the Weibull model and the cumulative hazard can express
the additive process which disabilities and diseases cause the human organism
during the life span leading to healthy life years lost. In this paper we
further analytically derive a more general model of survival-mortality in which
we estimate a parameter related to the Healthy Life Years Lost (HLYL) and
leading to the Weibull model and the corresponding shape parameter as a
specific case. We have also demonstrated that the results found for the general
HLYL parameter we have proposed provides results similar to those provided by
the World Health Organization for the Healthy Life Expectancy (HALE) and the
corresponding HLYL estimates. An analytic derivation of the mathematical
formulas is presented along with an easy to apply Excel program. This program
is an extension of the classical life table including four more columns to
estimate the cumulative mortality, the average mortality, the person life years
lost and finally the HLYL parameter bx. The latest versions of this program
appear in the Demographics2019 website


  The Oxford Nanopore Technologies's MinION is the first portable DNA
sequencing device. It is capable of producing long reads, over 100 kBp were
reported. However, it has significantly higher error rate than other methods.
In this study, we present MinCall, an end2end basecaller model for the MinION.
The model is based on deep learning and uses convolutional neural networks
(CNN) in its implementation. For extra performance, it uses cutting edge deep
learning techniques and architectures, batch normalization and Connectionist
Temporal Classification (CTC) loss. The best performing deep learning model
achieves 91.4% median match rate on E. Coli dataset using R9 pore chemistry and
1D reads.


  In this paper, we propose a semi-supervised deep learning method for
detecting the specific types of reads that impede the de novo genome assembly
process. Instead of dealing directly with sequenced reads, we analyze their
coverage graphs converted to 1D-signals. We noticed that specific signal
patterns occur in each relevant class of reads. Semi-supervised approach is
chosen because manually labelling the data is a very slow and tedious process,
so our goal was to facilitate the assembly process with as little labeled data
as possible. We tested two models to learn patterns in the coverage graphs:
M1+M2 and semi-GAN. We evaluated the performance of each model based on a
manually labeled dataset that comprises various reads from multiple reference
genomes with respect to the number of labeled examples that were used during
the training process. In addition, we embedded our detection in the assembly
process which improved the quality of assemblies.


  Chronic exposure to inorganic arsenic is a potential cause of carcinogenesis.
It elicits its potential by generation of ROS, leading to DNA, protein and
lipid damage. Therefore, the deleterious effect of arsenic can be mitigated by
quenching ROS using antioxidants. There is a homology between the protein
coding regions of mice and human. Effect of these alterations in human can be
mimicked in mice. Therefore to understand the underlying mechanism of arsenic
toxicity and its amelioration by black tea, studies have been conducted in mice
model. Long term exposure to iAs leads to tumour growth, which has been found
to be alleviated by black tea. Observations reveal that black tea has two
salutary effects on the growth of tumour: the rate of growth of damaged cells
was appreciably reduced and an early saturation of the level of damage is
achieved. To take the experimental findings further, the experimental data have
been modelled with simple dynamical equations. The curves obtained from
\textit{in vivo} studies have been fitted with the data obtained from the
model. The corresponding steady states and their stabilities are analyzed.


  We consider the problem of estimating the elapsed time since the most recent
common ancestor of a finite random sample drawn from a population which has
evolved through a Bienayme-Galton-Watson branching process. More specifically,
we are interested in the diffusion limit appropriate to a supercritical process
in the near-critical limit evolving over a large number of time steps. Our
approach differs from earlier analyses in that we assume the only known
information is the mean and variance of the number of offspring per parent, the
observed total population size at the time of sampling, and the size of the
sample. We obtain a formula for the probability that a finite random sample of
the population is descended from a single ancestor in the initial population,
and derive a confidence interval for the initial population size in terms of
the final population size and the time since initiating the process. We also
determine a joint likelihood surface from which confidence regions can be
determined for simultaneously estimating two parameters, (1) the population
size at the time of the most recent common ancestor, and (2) the time elapsed
since the existence of the most recent common ancestor.


  14-3-3 proteins are a group of seven dimeric adapter proteins that exert
their biological function by interacting with hundreds of phosphorylated
proteins, thus influencing their sub-cellular localization, activity or
stability in the cell. Due to this remarkable interaction network, 14-3-3
proteins have been associated with several pathologies and the protein-protein
interactions established with a number of partners are now considered promising
drug targets. The activity of 14-3-3 proteins is often isoform specific and to
our knowledge only one out of seven isoforms, 14-3-3$\zeta$, has been assigned.
Despite the availability of the crystal structures of all seven isoforms of
14-3-3, the additional NMR assignments of 14-3-3 proteins are important for
both biological mechanism studies and chemical biology approaches. Herein, we
present a robust backbone assignment of 14-3-3$\sigma$, which will allow
advances in the discovery of potential therapeutic compounds. This assignment
is now being applied to the discovery of both inhibitors and stabilizers of
14-3-3 protein-protein interactions.


  Explaining coexistence in species-rich communities of primary producers
remains a challenge for ecologists because of their likely competition for
shared resources. Following Hutchinson's seminal suggestion, many theoreticians
have tried to create diversity through a fluctuating environment, which impairs
or slows down competitive exclusion. However, fluctuating-environment models
often only produce a dozen of coexisting species at best. Here, we investigate
how to create richer communities in fluctuating environments, using an
empirically parameterized model. Building on the forced Lotka-Volterra model of
Scranton and Vasseur (Theor Ecol 9(3):353-363, 2016), inspired by phytoplankton
communities, we have investigated the effect of two coexistence mechanisms,
namely the storage effect and higher intra- than interspecific competition
strengths (i.e., strong self-regulation). We tuned the intra/inter competition
ratio based on empirical analyses, in which self-regulation dominates
interspecific interactions. Although a strong self-regulation maintained more
species (50%) than the storage effect (25%), we show that none of the two
coexistence mechanisms considered could ensure the coexistence of all species
alone. Realistic seasonal environments only aggravated that picture, as they
decreased persistence relative to a random environment. However, strong
self-regulation and the storage effect combined superadditively so that all
species could persist with both mechanisms at work. Our results suggest that
combining different coexistence mechanisms into community models might be more
fruitful than trying to find which mechanism best explains diversity. We
additionally highlight that while biomass-trait distributions provide some
clues regarding coexistence mechanisms, they cannot indicate unequivocally
which mechanisms are at play.


  In recent years, there has been growing interest in studying evolutionary
games with environmental feedback. Previous studies exclusively focus on
two-player games. However, extension to multi-player game is needed to study
problems such as microbial cooperation and crowdsourcing collaborations. Here,
we study coevolutionary public goods games where strategies coevolve with the
multiplication factors of group cooperation. Asymmetry can arise in such
environmental feedback, where games organized by focal cooperators may have a
different efficiency than the ones by defectors. Our analysis shows that
co-evolutionary dynamics with asymmetrical environmental feedback can yield
oscillatory convergence to persistent cooperation, if the relative changing
speed of cooperators' multiplication factor is above a certain threshold. Our
work provides useful insights into sustaining group cooperation in a changing
world.


  While ubiquitous, energy redistribution remains a poorly understood facet of
the nonequilibrium thermodynamics of biomolecules. At the molecular level,
finite-size effects, pronounced nonlinearities, and ballistic processes produce
behavior that diverges from the macroscale. Here, we show that transient
thermal transport reflects macromolecular energy landscape architecture through
the topological characteristics of molecular contacts and the nonlinear
processes that mediate dynamics. While the former determines transport pathways
via pairwise interactions, the latter reflects frustration within the landscape
for local conformational rearrangements. Unlike transport through
small-molecule systems, such as alkanes, nonlinearity dominates over coherent
processes at even quite short time- and length-scales. Our exhaustive all-atom
simulations and novel local-in-time and space analysis, applicable to both
theory and experiment, permit dissection of energy migration in biomolecules.
The approach demonstrates that vibrational energy transport can probe otherwise
inaccessible aspects of macromolecular dynamics and interactions that underly
biological function.


  Using a Multi-agent systems paradigm, the present project develops, validates
and exploits a computational $\textit{testbed}$ that simulates micro-biological
complex systems, namely the aggregation patterns of the social amoeba
$\textit{Dyctiostelium discoideum}$. We propose a new design and implementation
for managing discrete simulations with autonomous agents on a microscopic
scale, thus focusing on their social behavior and mutual interactions. Then,
the dependence on the main physical variables is tested, namely density and
number of amoebas; in addition, we analyze the robustness of the dynamics
against various noise sources. Along with these results, we suggest a
methodology for further studies that make use of our validated model.


  DNA supercoiling, the under or overwinding of DNA, is a key physical
mechanism both participating to compaction of bacterial genomes and making
genomic sequences adopt various structural forms. DNA supercoiling may lead to
the formation of braided superstructures (plectonemes), or it may locally
destabilize canonical B-DNA to generate denaturation bubbles, left-handed Z-DNA
and other functional alternative forms. Prediction of the relative fraction of
these structures has been limited because of a lack of predictive polymer
models that can capture the multiscale properties of long DNA molecules. In
this work, we address this issue by extending the self-avoiding rod-like chain
model of DNA so that every site of the chain is allocated with an additional
structural degree of freedom reflecting variations of DNA forms. Efficient
simulations of the model reveal its relevancy to capture multiscale properties
of long chains (here up to 21 kb) as reported in magnetic tweezers experiments.
Well-controlled approximations further lead to accurate analytical estimations
of thermodynamic properties in the high force regime, providing, in combination
with experiments, a simple, yet powerful framework to infer physical parameters
describing alternative forms. In this regard, using simulated data, we find
that extension curves at forces above 2 pN may lead, alone, to erroneous
parameter estimations as a consequence of an underdetermination problem. We
thus revisit published data in light of these findings and discuss the
relevancy of previously proposed sets of parameters for both denatured and
left-handed DNA forms. Altogether, our work paves the way for a scalable
quantitative model of bacterial DNA.


  It is well known that stochastically modeled reaction networks that are
complex balanced admit a stationary distribution that is a product of Poisson
distributions. In this paper, we consider the following related question:
supposing that the initial distribution of a stochastically modeled reaction
network is a product of Poissons, under what conditions will the distribution
remain a product of Poissons for all time? By drawing inspiration from Crispin
Gardiner's "Poisson representation" for the solution to the chemical master
equation, we provide a necessary and sufficient condition for such a
product-form distribution to hold for all time. Interestingly, the condition is
a dynamical "complex-balancing" for only those complexes that have multiplicity
greater than or equal to two (i.e. the higher order complexes that yield
non-linear terms to the dynamics). We term this new condition the "dynamical
and restricted complex balance" condition (DR for short).


  The distributive sequential n-site phosphorylation/dephosphorylation system
is an important building block in networks of chemical reactions arising in
molecular biology, which has been intensively studied. In the nice paper of
Wang and Sontag (2008) it is shown that for certain choices of the reaction
rate constants and total conservation constants, the system can have 2[n/2]+1
positive steady states (that is, n+1 positive steady states for n even and n
positive steady states for n odd). In this paper we give open parameter regions
in the space of reaction rate constants and total conservation constants that
ensure these number of positive steady states, while assuming in the modeling
that roughly only 1/4 of the intermediates occur in the reaction mechanism.
This result is based on the general framework developed by Bihan, Dickenstein,
and Giaroli (2018), which can be applied to other networks. We also describe
how to implement these tools to search for multistationarity regions in a
computer algebra system and present some computer aided results.


  Self-organization of kinesin-driven, microtubule-based 3D active fluids
relies on the collective dynamics of single microtubules. However, the
connection between macroscopic fluid flows and microscopic motion of
microtubules remains unclear. In this work, the motion of single microtubules
was characterized by means of 2D gliding assays and compared with the flows of
3D active fluids. While the scales of the two systems differ by 1,000x, both
were driven by processive, non-processive or an equal mixture of both molecular
motor proteins. To search for the dynamic correlation between both systems, the
motor activities were tuned by varying temperature and ATP concentration, and
the changes in both systems were compared. Motor processivity played an
important role in active fluid flows but only when the fluids were nearly
motionless; otherwise, flows were dominated by hydrodynamic resistance
controlled by sample size. Furthermore, while the motors' thermal reaction led
active fluids to flow faster with increasing temperature, such temperature
dependence could be reversed by introducing temperature-varying depletants,
emphasizing the potential role of the depletant in designing an active fluid's
temperature response. The temperature response of active fluids was nearly
immediate (<10 sec). Such a characteristic enables active fluids to be
controlled with a temperature switch. Overall, this work not only clarifies the
role of temperature in active fluid activity, but also sheds light on the
underlying principles of the relationship between the collective dynamics of
active fluids and the dynamics of their constituent single microtubules.


  A wide variety of problems in computational biology, most notably the
assessment of orthology, are solved with the help of reciprocal best matches.
Using an evolutionary definition of best matches that captures the intuition
behind the concept we clarify rigorously the relationships between reciprocal
best matches, orthology, and evolutionary events under the assumption of
duplication/loss scenarios. We show that the orthology graph is a subgraph of
the reciprocal best match graph (RBMG). We furthermore give conditions under
which an RBMG that is a cograph identifies the correct orthlogy relation. Using
computer simulations we find that most false positive orthology assignments can
be identified as so-called good quartets -- and thus corrected -- in the
absence of horizontal transfer. Horizontal transfer, however, may introduce
also false-negative orthology assignments.


  We first present a predator-prey model for two species and then extend the
model to three species where the two predator species engage in mutualistic
predation. Constant effort harvesting and the impact of by-catch issue are also
incorporated. Necessary sufficient conditions for the existence and stability
of positive equilibrium points are examined. It is shown that harvesting is
sustainable, and the memory concept of the fractional derivative damps out
oscillations in the population numbers so that the system as a whole settles on
an equilibrium quicker than it would with integer time derivatives. Finally,
some possible physical explanations are given for the obtained results. It is
shown that the stability requires the memory concept in the model.


  The advent of plant phenomics, coupled with the wealth of genotypic data
generated by next-generation sequencing technologies, provides exciting new
resources for investigations into and improvement of complex traits. However,
these new technologies also bring new challenges in quantitative genetics,
namely, a need for the development of robust frameworks that can accommodate
these high-dimensional data. In this chapter, we describe methods for the
statistical analysis of high-throughput phenotyping (HTP) data with the goal of
enhancing the prediction accuracy of genomic selection (GS). Following the
Introduction in Section 1, Section 2 discusses field-based HTP, including the
use of unmanned aerial vehicles and light detection and ranging, as well as how
we can achieve increased genetic gain by utilizing image data derived from HTP.
Section 3 considers extending commonly used GS models to integrate HTP data as
covariates associated with the principal trait response, such as yield.
Particular focus is placed on single-trait, multi-trait, and genotype by
environment interaction models. One unique aspect of HTP data is that phenomics
platforms often produce large-scale data with high spatial and temporal
resolution for capturing dynamic growth, development, and stress responses.
Section 4 discusses the utility of a random regression model for performing
longitudinal GS. The chapter concludes with a discussion of some standing
issues.


  During muscle contraction, myosin motors anchored to thick filaments bind to
and slide actin thin filaments. These motors rely on energy derived from ATP,
supplied, in part, by diffusion from the sarcoplasm to the interior of the
lattice of actin and myosin filaments. The radial spacing of filaments in this
lattice may change or remain constant during contraction. If the lattice is
isovolumetric, it must expand when the muscle shortens. If, however, the
spacing is constant or has a different pattern of axial and radial motion, then
the lattice changes volume during contraction, driving fluid motion and
assisting in the transport of molecules between the contractile lattice and the
surrounding intracellular space. We first create an
advective-diffusive-reaction flow model and show that the flow into and out of
the sarcomere lattice would be significant in the absence of lattice expansion.
Advective transport coupled to diffusion has the potential to substantially
enhance metabolite exchange within the crowded sarcomere. Using time-resolved
x-ray diffraction of contracting muscle, we next show that the contractile
lattice is neither isovolumetric nor constant in spacing. Instead, lattice
spacing is time-varying, depends on activation, and can manifest as an
effective time-varying Poisson ratio. The resulting fluid flow in the sarcomere
lattice of synchronous insect flight muscles is greater than expected for
constant lattice spacing conditions. Lattice spacing depends on a variety of
factors that produce radial force, including crossbridges, titin-like
molecules, and other structural proteins. Volume change and advective transport
varies with the phase of muscle stimulation but remains significant at all
conditions. Akin to "breathing," advective-diffusive transport in sarcomeres is
sufficient to promote metabolite exchange and may play a role in the regulation
of contraction itself.


  'Tree-based' phylogenetic networks proposed by Francis and Steel have
attracted much attention of theoretical biologists in the last few years. At
the heart of the definitions of tree-based phylogenetic networks is the notion
of 'support trees', about which there are numerous algorithmic problems that
are important for evolutionary data analysis. Recently, Hayamizu
(arXiv:1811.05849 [math.CO]) proved a structure theorem for tree-based
phylogenetic networks and obtained linear-time and linear-delay algorithms for
many basic problems on support trees, such as counting, optimisation, and
enumeration. In the present paper, we consider the following fundamental
problem in statistical data analysis: given a tree-based phylogenetic network
$N$ whose arcs are associated with probability, create the top-$k$ support tree
ranking for $N$ by their likelihood values. We provide a linear-delay (and
hence optimal) algorithm for the problem and thus reveal the interesting
property of tree-based phylogenetic networks that ranking top-$k$ support trees
is as computationally easy as picking $k$ arbitrary support trees.


  Background: Many genome-wide association studies have detected genomic
regions associated with traits, yet understanding the functional causes of
association often remains elusive. Utilizing systems approaches and focusing on
intermediate molecular phenotypes might facilitate biologic understanding.
Results: The availability of exome sequencing of two populations of
African-Americans and European-Americans from the Atherosclerosis Risk in
Communities study allowed us to investigate the effects of annotated
loss-of-function (LoF) mutations on 122 serum metabolites. To assess the
findings, we built metabolomic causal networks for each population separately
and utilized structural equation modeling. We then validated our findings with
a set of independent samples. By use of methods based on concepts of Mendelian
randomization of genetic variants, we showed that some of the affected
metabolites are risk predictors in the causal pathway of disease. For example,
LoF mutations in the gene KIAA1755 were identified to elevate the levels of
eicosapentaenoate (p-value=5E-14), an essential fatty acid clinically
identified to increase essential hypertension. We showed that this gene is in
the pathway to triglycerides, where both triglycerides and essential
hypertension are risk factors of metabolomic disorder and heart attack. We also
identified that the gene CLDN17, harboring loss-of-function mutations, had
pleiotropic actions on metabolites from amino acid and lipid pathways.
Conclusion: Using systems biology approaches for the analysis of metabolomics
and genetic data, we integrated several biological processes, which lead to
findings that may functionally connect genetic variants with complex diseases.


  Surveys of microbial biodiversity such as the Earth Microbiome Project (EMP)
and the Human Microbiome Project (HMP) have revealed robust ecological patterns
across different environments. A major goal in ecology is to leverage these
patterns to identify the ecological processes shaping microbial ecosystems. One
promising approach is to use minimal models that can relate mechanistic
assumptions at the microbe scale to community-level patterns. Here, we
demonstrate the utility of this approach by showing that the Microbial Consumer
Resource Model (MiCRM) -- a minimal model for microbial communities with
resource competition, metabolic crossfeeding and stochastic colonization -- can
qualitatively reproduce patterns found in survey data including compositional
gradients, dissimilarity/overlap correlations, richness/harshness correlations,
and nestedness of community composition. By using the MiCRM to generate
synthetic data with different environmental and taxonomical structure, we show
that large scale patterns in the EMP can be reproduced by considering the
energetic cost of surviving in harsh environments and HMP patterns may reflect
the importance of environmental filtering in shaping competition. We also show
that recently discovered dissimilarity-overlap correlations in the HMP likely
arise from communities that share similar environments rather than reflecting
universal dynamics. We identify ecologically meaningful changes in parameters
that alter or destroy each one of these patterns, suggesting new mechanistic
hypotheses for further investigation. These findings highlight the promise of
minimal models for microbial ecology.


  The humanity has been facing a plethora of challenges associated with
infectious diseases, which kill more than 6 million people a year. Although
continuous efforts have been applied to relieve the potential damages from such
misfortunate events, it is unquestionable that there are many persisting
challenges yet to overcome. One related issue we particularly address here is
the assessment and prediction of such epidemics. In this field of study,
traditional and ad-hoc models frequently fail to provide proper predictive
situation awareness (PSAW), characterized by understanding the current
situations and predicting the future situations. Comprehensive PSAW for
infectious disease can support decision making and help to hinder disease
spread. In this paper, we develop a computing system platform focusing on
collective intelligence causal modeling, in order to support PSAW in the domain
of infectious disease. Analyses of global epidemics require integration of
multiple different data and models, which can be originated from multiple
independent researchers. These models should be integrated to accurately assess
and predict the infectious disease in terms of holistic view. The system shall
provide three main functions: (1) collaborative causal modeling, (2) causal
model integration, and (3) causal model reasoning. These functions are
supported by subject-matter expert and artificial intelligence (AI), with
uncertainty treatment. Subject-matter experts, as collective intelligence,
develop causal models and integrate them as one joint causal model. The
integrated causal model shall be used to reason about: (1) the past, regarding
how the causal factors have occurred; (2) the present, regarding how the spread
is going now; and (3) the future, regarding how it will proceed. Finally, we
introduce one use case of predictive situation awareness for the Ebola virus
disease.


  Fluorescence microscopy has led to impressive quantitative models and new
insights gained from richer sets of biomedical imagery. However, there is a
dearth of rigorous and established bioimaging strategies for modeling
spatiotemporal behavior of diffuse, subcellular components such as mitochondria
or actin. In many cases, these structures are assessed by hand or with other
semi-quantitative measures. We propose to build descriptive and dynamic models
of diffuse subcellular morphologies, using the mitochondrial protein patterns
of cervical epithelial (HeLa) cells. We develop a parametric representation of
the patterns as a mixture of probability masses. This mixture is iteratively
perturbed over time to fit the evolving spatiotemporal behavior of the
subcellular structures. We convert the resulting trajectory into a series of
graph Laplacians to formally define a dynamic network. Finally, we demonstrate
how graph theoretic analyses of the trajectories yield biologically-meaningful
quantifications of the structures.


  Every now and then the cultural paradigm of a society changes. Human history
can be regarded as a sequence of long periods of cultural stasis punctuated by
paradigm shifts that transform culture upside-down over the turn of a few
generations. We propose here a population dynamics model devised to analyse
paradigm shifts. In this model individuals are defined by a vector of cultural
traits that can change mainly through imitation of other individuals' traits.
The novelty of the model is that cultural traits may interact reinforcing or
hindering each other. Imitation is then biased by the 'cultural fitness'
landscape thus defined. Our main result is that abrupt paradigm shifts occur,
as a response to weak changes in the landscape, only when cultural traits do
interact---whereas adaptation is smooth if there is no interaction. Borrowing
the genetic term, this interaction is called 'cultural epistasis'. The result
is robust to the way that epistasis is implemented, to whether imitation is
biased by homophily, or to changes in other model parameters. Finally, a
relevant consequence of this dynamics is the irreversible nature of paradigm
shifts: the old paradigm cannot be restored even if the external changes are
undone. Our model puts the phenomenon of paradigm shifts in cultural evolution
in the same category as catastrophic shifts in ecology or phase transitions in
physics.


  We introduce a method to generate synthetic protein sequences which are
predicted to be resistant to certain antibiotics. We did this using 6,023 genes
that were predicted to be resistant to antibiotics in the intestinal region of
the human gut and were fed as input to a Wasserstein generative adversarial
network (W-GAN) model a variant to the original generative adversarial model
which has been known to perform efficiently when it comes to mimicking the
distribution of the real data in order to generate new data which is similar in
style to the original data which was fed as the training data


  Understanding individual decisions in a world where communications and
information move instantly via cell phones and the internet, contributes to the
development and implementation of policies aimed at stopping or ameliorating
the spread of diseases. In this manuscript, the role of official social network
perturbations generated by public health officials to slow down or stop a
disease outbreak are studied over distinct classes of static social networks.
The dynamics are stochastic in nature with individuals (nodes) being assigned
fixed levels of education or wealth. Nodes may change their epidemiological
status from susceptible, to infected and to recovered. Most importantly, it is
assumed that when the prevalence reaches a pre-determined threshold level, P*,
information, called awareness in our framework, starts to spread, a process
triggered by public health authorities. Information is assumed to spread over
the same static network and whether or not one becomes a temporary informer, is
a function of his/her level of education or wealth and epidemiological status.
Stochastic simulations show that threshold selection P* and the value of the
average basic reproduction number impact the final epidemic size
differentially. For the Erdos-Renyi and Small-world networks, an optimal choice
for P* that minimize the final epidemic size can be identified under some
conditions while for Scale-free networks this is not case.


  In recent years the publication of genome sequences for the Chinese hamster
and Chinese hamster ovary (CHO) cell lines have facilitated study of these
biopharmaceutical cell factories with unprecedented resolution. Our
understanding of the CHO cell transcriptome, in particular, has rapidly
advanced through the application of next-generation sequencing (NGS) technology
to characterise RNA expression (RNA-Seq). In this chapter we present a
computational pipeline for the analysis of CHO cell RNA-Seq data from the
Illumina platform to identify differentially expressed genes. The example data
and bioinformatics workflow required to run this analysis are freely available
at www.cgcdb.org/rnaseq_analysis_protocol.html.


  At the scale of the individual cell, protein production is a stochastic
process with multiple time scales, combining quick and slow random steps with
discontinuous and smooth variation. Hybrid stochastic processes, in particular
piecewise-deterministic Markov processes (PDMP), are well adapted for
describing such situations. PDMPs approximate the jump Markov processes
traditionally used as models for stochastic chemical reaction networks.
Although hybrid modelling is now well established in biology, these models
remain computationally challenging. We propose several improved methods for
computing time dependent multivariate probability distributions (MPD) of PDMP
models of gene networks. In these models, the promoter dynamics is described by
a finite state, continuous time Markov process, whereas the mRNA and protein
levels follow ordinary differential equations (ODEs). The Monte-Carlo method
combines direct simulation of the PDMP with analytic solutions of the ODEs. The
push-forward method numerically computes the probability measure advected by
the deterministic ODE flow, through the use of analytic expressions of the
corresponding semigroup. Compared to earlier versions of this method, the
probability of the promoter states sequence is computed beyond the naive mean
field theory and adapted for non-linear regulation functions.


  To explain the stability of cooperation is a central task of evolutionary
theory. We investigate this question in the case of quorum sensing (QS)
bacteria, which regulate cooperative traits in response to population density.
Cooperation is modeled by the prisoner's dilemma, where individuals produce a
costly public good (PG) that equally benefits all members of a community
divided into multiple, distinct patches (multilevel selection). Cost and
benefit are non-linear functions of the PG production. The analysis of
evolutionary stability yields an optimization problem for the expression of PG
in dependency on the number of QS individuals within a colony. We find that the
optimal total PG production of the QS population mainly depends on the shape of
the benefit. A graded and a switch-like response is possible, in accordance
with earlier results. Interestingly, at the level of the individual cell, the
QS response is determined by the shape of the costs. All QS individuals respond
either homogeneously if cost are a convex function of the PG production rate,
or they respond heterogeneously with distinct ON/OFF responses if the costs are
concave. The latter finding is consistent with recent experimental findings,
and contradicts the usual interpretation of QS as a mechanism to establish a
uniform, synchronized response of a bacterial population.


  The biological function of a protein stems from its 3-dimensional structure,
which is thermodynamically determined by the energetics of interatomic forces
between its amino acid building blocks (the order of amino acids, known as the
sequence, defines a protein). Given the costs (time, money, human resources) of
determining protein structures via experimental means such as X-ray
crystallography, can we better describe and compare protein 3D structures in a
robust and efficient manner, so as to gain meaningful biological insights? We
begin by considering a relatively simple problem, limiting ourselves to just
protein secondary structural elements. Historically, many computational methods
have been devised to classify amino acid residues in a protein chain into one
of several discrete secondary structures, of which the most well-characterized
are the geometrically regular $\alpha$-helix and $\beta$-sheet; irregular
structural patterns, such as 'turns' and 'loops', are less understood. Here, we
present a study of Deep Learning techniques to classify the loop-like end cap
structures which delimit $\alpha$-helices. Previous work used highly empirical
and heuristic methods to manually classify helix capping motifs. Instead, we
use structural data directly--including (i) backbone torsion angles computed
from 3D structures, (ii) macromolecular feature sets (e.g., physicochemical
properties), and (iii) helix cap classification data (from CAPS-DB)--as the
ground truth to train a bidirectional long short-term memory (BiLSTM) model to
classify helix cap residues. We tried different network architectures and
scanned hyperparameters in order to train and assess several models; we also
trained a Support Vector Classifier (SVC) to use as a baseline. Ultimately, we
achieved 85% class-balanced accuracy with a deep BiLSTM model.


  We analyse the mechanism of enzyme-substrate catalysis from the perspective
of minimizing the load on the enzymes through sequestration, whilst maintaining
at least a minimum reaction flux. In particular, we ask: which binding free
energies of the enzyme-substrate and enzyme-product reaction intermediates
minimize the fraction of enzymes sequestered in complexes, while sustaining at
a certain minimal flux? Under reasonable biophysical assumptions, we find that
the optimal design will saturate the bound on the minimal flux, and reflects a
basic trade-off in catalytic operation. If both binding free energies are too
high, there is low sequestration, but the effective progress of the reaction is
hampered. If both binding free energies are too low, there is high
sequestration, and the reaction flux may also be suppressed in extreme cases.
The optimal binding free energies are therefore neither too high nor too low,
but in fact moderate. Moreover, the optimal difference in substrate and product
binding free energies, which contributes to the thermodynamic driving force of
the reaction, is in general strongly constrained by the intrinsic free-energy
difference between products and reactants. Both the strategies of using a
negative binding free-energy difference to drive the catalyst-bound reaction
forward, and of using a positive binding free-energy difference to enhance
detachment of the product, are limited in their efficacy.


  African Swine Fever (ASF) is viral infection which causes acute disease in
domestic pigs and wild boar. Although the virus does not cause disease in
humans, the impact it has on the economy, especially through trade and farming,
is substantial. Recent rapid propagation of the (ASF) from East to West of
Europe encouraged us to prepare risk assessment for Poland. The early growth
estimation can be easily done by matching incidence trajectory to the
exponential function, resulting in the approximation of the force of infection.
With these calculations the basic reproduction rate of the epidemic, the
effective outbreaks detection and elimination times could be estimated. In
regression mode, 380 Polish counties (poviats) have been analysed, where 18
(located in Northeast Poland) have been affected (until August 2017) for
spatial propagation (risk assessment for future). Mathematical model has been
applied by taking into account: swine amount significance, disease vectors
(wild boards) significance. We use pseudogravitational models of short and
longrange interactions referring to the socio-migratory behavior of wild boars
and the pork production chain significance. Spatial modeling in a certain range
of parameters proves the existence of a natural protective barrier within
boarders of the Congress Poland. The spread of the disease to the Greater
Poland should result in the accelerated outbreak of ASF production chain. In
the preliminary setup, we perform regression analysis, network outbreak
investigation, early epidemic growth estimation and simulate landscape-based
propagation.


  Most living systems rely on double-stranded DNA (dsDNA) to store their
genetic information and perpetuate themselves. This biological information has
been considered the main target of evolution. However, here we show that
symmetries and patterns in the dsDNA sequence can emerge from the physical
peculiarities of the dsDNA molecule itself and the maximum entropy principle
alone, rather than from biological or environmental evolutionary pressure. The
randomness justifies the human codon biases and context-dependent mutation
patterns in human populations. Thus, the DNA "exceptional symmetries", emerged
from the randomness, have to be taken into account when looking for the DNA
encoded information. Our results suggest that the double helix energy
constraints and, more generally, the physical properties of the dsDNA are the
hard drivers of the overall DNA sequence architecture, whereas the biological
selective processes act as soft drivers, which only under extraordinary
circumstances overtake the overall entropy content of the genome.


  We developed an integrated hierarchical Bayesian life cycle model that
simultaneously estimates the abundance of post-smolts at sea, post-smolt
survival rates, and proportions maturing as 1SW, for all SU in Northern Europe,
Southern Europe and North America. The model is an age- and stage-based life
cycle model that considers 1SW and 2SW life history strategies and harmonizes
the life history dynamics among SU in North America and Europe. The new
framework brought a major contribution to improve the scientific basis for
Atlantic salmon stock assessment. It is a benchmark for the assessment and
forecast models currently used by ICES for Atlantic salmon stock assessment in
the North Atlantic. ...


  Background. Non-specific Lipid Transfer Proteins (nsLTPs) are widely
distributed in the plant kingdom and constitute a superfamily of related
proteins. More than 800 different sequences have been characterized so far, but
their biological functions remain unclear. It has been clear for years that
they present a certain interest for agronomic and nutritional issues.
Deciphering their functions means collecting and analyzing a variety of data
from gene sequence to protein structure, from cellular localization to the
physiological role. As a huge and growing number of new protein sequences are
available nowadays, extracting meaningful knowledge from
sequence-structure-function relationships calls for the development of new
tools and approaches. As nsLTPs show high evolutionary divergence, but a
conserved common right-handed superhelix structural fold, and as they are
involved in a large number of key roles in plant development and defense, they
are a stimulating case study for validating such an approach.


  We present a series of SIR-network models, extended with a game-theoretic
treatment of imitation dynamics which result from regular population mobility
across residential and work areas and the ensuing interactions. Each considered
SIR-network model captures a class of vaccination behaviours influenced by
epidemic characteristics, interaction topology, and imitation dynamics. Our
focus is the eventual vaccination coverage, produced under voluntary
vaccination schemes, in response to these varying factors. Using the next
generation matrix method, we analytically derive and compare expressions for
the basic reproduction number $R_0$ for the proposed SIR-network models.
Furthermore, we simulate the epidemic dynamics over time for the considered
models, and show that if individuals are sufficiently responsive towards the
changes in the disease prevalence, then the more expansive travelling patterns
encourage convergence to the endemic, mixed equilibria. On the contrary, if
individuals are insensitive to changes in the disease prevalence, we find that
they tend to remain unvaccinated in all the studied models. Our results concur
with earlier studies in showing that residents from highly connected
residential areas are more likely to get vaccinated. We also show that the
existence of the individuals committed to receiving vaccination reduces $R_0$
and delays the disease prevalence, and thus is essential to containing
epidemics.


  Stochastic simulation is a widely used method for estimating quantities in
models of chemical reaction networks where uncertainty plays a crucial role.
However, reducing the statistical uncertainty of the corresponding estimators
requires the generation of a large number of simulation runs, which is
computationally expensive. To reduce the number of necessary runs, we propose a
variance reduction technique based on control variates. We exploit constraints
on the statistical moments of the stochastic process to reduce the estimators'
variances. We develop an algorithm that selects appropriate control variates in
an on-line fashion and demonstrate the efficiency of our approach on several
case studies.


  The dynamics of a SIVR model with power relationship incidence rates $(\beta
I^p S^q)$ is investigated. It is assumed an individual can be susceptible after
receiving the first dose of the vaccine, hence a second dose is required to
attain permanent immunity. The steady states conditions of the disease-free
equilibrium and the endemic equilibrium are critically presented. Numerical
simulations are carried out to determine the impact of the exponential
parameters $(p;q)$ on infection.