  Mathematical models of cardiac action potentials have become increasingly
important in the study of heart disease and pharmacology, but concerns linger
over their robustness during long periods of simulation, in particular due to
issues such as model drift and non-unique steady states. Previous studies have
linked these to violation of conservation laws, but only explored those issues
with respect to charge conservation in specific models. Here, we propose a
general and systematic method of identifying conservation laws hidden in models
of cardiac electrophysiology by using bond graphs, and develop a bond graph
model of the cardiac action potential to study long-term behaviour. Bond graphs
provide an explicit energy-based framework for modelling physical systems,
which makes them well-suited for examining conservation within
electrophysiological models. We find that the charge conservation laws derived
in previous studies are examples of the more general concept of a "conserved
moiety". Conserved moieties explain model drift and non-unique steady states,
generalising the results from previous studies. The bond graph approach
provides a rigorous method to check for drift and non-unique steady states in a
wide range of cardiac action potential models, and can be extended to examine
behaviours of other excitable systems.


  Nano-communication has gained significant attention in the last few years, as
a means to establish information transfer between future nano-machines.
Comparing with other communication techniques for nano-scale (calcium ions
signaling, molecular or catalytic nanomotors, pheromones propagation,
bacteria-based communication), the phenomenon called Forster Resonance Energy
Transfer (FRET) offers significantly smaller propagation delays and high
channel throughput. In this paper, we report our recent experiments on
FRET-based nano-networks performed in the Laboratory of Cell Biophysics of the
Jagiellonian University, Krakow. We propose to use Alexa Fluor dyes as nano
transmitters and receivers, as they enable to create multiple-input
multi-output (MIMO) FRET communication channels and thus enhance FRET
efficiency. We measure FRET efficiency values, calculate bit error rates for
the measured scenarios and extend the calculations to consider a general case
of MIMO (n,m) FRET channels.


  The phenomenon of Forster Resonance Energy Transfer, commonly used to measure
the distances between fluorophore molecules and to study interactions between
fluorescent-tagged proteins in life sciences, can also be applied in
nanocommunication networks to transfer information bits. The mechanism offers a
relatively large throughput and very small delays, but at the same time the
channel bit error rate is too high and the transmission ranges are too limited
for communication purposes. In this paper, multiple donors at the transmitter
side and multiple acceptors at the receiver side are considered to decrease the
bit error rate. As nanoantennas, the DyLight fluorescent dyes, which are very
well suited to long range nanocommunication due to their large Forster
distances and high degrees of labeling, are proposed. The reported results of
the recent laboratory experiments confirm efficient communication on distances
over 10 nm.


  Nanocommunications, understood as communications between nanoscale devices,
is commonly regarded as a technology essential for cooperation of large groups
of nanomachines and thus crucial for development of the whole area of
nanotechnology. While solutions for point-to-point nanocommunications have been
already proposed, larger networks cannot function properly without routing. In
this article we focus on the nanocommunications via Forster Resonance Energy
Transfer (FRET), which was found to be a technique with a very high signal
propagation speed, and discuss how to route signals through nanonetworks. We
introduce five new routing mechanisms, based on biological properties of
specific molecules. We experimentally validate one of these mechanisms.
Finally, we analyze open issues showing the technical challenges for signal
transmission and routing in FRET-based nanocommunications.


  The paper is concerned with Forster Resonance Energy Transfer (FRET)
considered as a mechanism for communication between nanodevices. Two solved
issues are reported in the paper, namely: signal generation and signal storage
in FRET-based nanonetworks. First, luciferase molecules as FRET transmitters
which are able to generate FRET signals themselves, taking energy from chemical
reactions without any external light exposure, are proposed. Second,
channelrhodopsins as FRET receivers, as they can convert FRET signals into
voltage, are suggested. Further, medical in-body systems where both molecule
types might be successfully applied, are discussed. Luciferase-channelrhodopsin
communication is modeled and its performance is numerically validated,
reporting on its throughput, bit error rate, propagation delay and energy
consumption.


  In tropical regions, fires propagate readily in grasslands but typically
consume only edges of forest patches. Thus forest patches grow due to tree
propagation and shrink by fires in surrounding grasslands. The interplay
between these competing edge effects is unknown, but critical in determining
the shape and stability of individual forest patches, as well the
landscape-level spatial distribution and stability of forests. We analyze
high-resolution remote-sensing data from protected areas of the Brazilian
Cerrado and find that forest shapes obey a robust perimeter-area scaling
relation across climatic zones. We explain this scaling by introducing a
heterogeneous fire propagation model of tropical forest-grassland ecotones.
Deviations from this perimeter-area relation determine the stability of
individual forest patches. At a larger scale, our model predicts that the
relative rates of tree growth due to propagative expansion and long-distance
seed dispersal determine whether collapse of regional-scale tree cover is
continuous or discontinuous as fire frequency changes.


  In this paper, it is presented a methodology for implementing arbitrarily
constructed time-homogenous Markov chains with biochemical systems. Not only
discrete but also continuous-time Markov chains are allowed to be computed. By
employing chemical reaction networks (CRNs) as a programmable language,
molecular concentrations serve to denote both input and output values. One
reaction network is elaborately designed for each chain. The evolution of
species' concentrations over time well matches the transient solutions of the
target continuous-time Markov chain, while equilibrium concentrations can
indicate the steady state probabilities. Additionally, second-order Markov
chains are considered for implementation, with bimolecular reactions rather
that unary ones. An original scheme is put forward to compile unimolecular
systems to DNA strand displacement reactions for the sake of future physical
implementations. Deterministic, stochastic and DNA simulations are provided to
enhance correctness, validity and feasibility.


  Approaches to gene interactions based on sign epistasis have been highly
influential in recent time. Sign epistasis is useful for relating local and
global properties of fitness landscapes, as well as for analyzing evolutionary
trajectories and constraints. The geometric theory of gene interactions, on the
other hand, provides complete information on interactions in terms of minimal
dependence relations. We propose a new framework that combines aspects of both
approaches. In particular, we provide efficient tools for identifying sign
epistasis and related order perturbations in large genetic systems, with
applications to the malaria-causing parasite Plasmodium vivax. We found that
order perturbations beyond sign epistasis are prevalent in the drug-free
environment, which agrees well with the observation that reversed evolution
back to the ancestral type is difficult. As a theoretical application, we
investigate how rank orders of genotypes with respect to fitness relates to
additivity.


  A practical, biologically motivated case of protein complexes (immunoglobulin
G and FcRII receptors) moving on the surface of mastcells, that are common
parts of an immunological system, is investigated. Proteins are considered as
nanomachines creating a nanonetwork. Accurate molecular models of the proteins
and the fluorophores which act as their nanoantennas are used to simulate the
communication between the nanomachines when they are close to each other. The
theory of diffusion-based Brownian motion is applied to model movements of the
proteins. It is assumed that fluorophore molecules send and receive signals
using the Forster Resonance Energy Transfer. The probability of the efficient
signal transfer and the respective bit error rate are calculated and discussed.


  In this paper, we study through mathematical modelling the combined effect of
transcriptional and translational regulation by proteins and small noncoding
RNAs (sRNA) in a genetic feedback motif that has an important role in the
survival of E.coli under stress associated with oxygen and energy availability.
We show that subtle changes in this motif can bring in drastically different
effects on the gene expression. In particular, we show that a threshold
response in the gene expression changes to a bistable response as the
regulation on sRNA synthesis or degradation is altered. These results are
obtained under deterministic conditions. Next, we study how the gene expression
is altered by additive and multiplicative noise which might arise due to
probabilistic occurrences of different biochemical events. Using the
Fokker-Planck formulation, we obtain steady state probability distributions for
sRNA concentration for the network motifs displaying bistability. The
probability distributions are found to be bimodal with two peaks at low and
high concentrations of sRNAs. We further study the variations in the
probability distributions under different values of noise strength and
correlations. The results presented here might be of interest for designing
synthetic network for artificial control.


  This paper provides a global picture about the deployment of networked
processing services for genomic data sets. Many current research make an
extensive use genomic data, which are massive and rapidly increasing over time.
They are typically stored in remote databases, accessible by using Internet.
For this reason, a significant issue for effectively handling genomic data
through data networks consists of the available network services. A first
contribution of this paper consists of identifying the still unexploited
features of genomic data that could allow optimizing their networked
management. The second and main contribution of this survey consists of a
methodological classification of computing and networking alternatives which
can be used to offer what we call the Genomic-as-a-Service (GaaS) paradigm. In
more detail, we analyze the main genomic processing applications, and classify
not only the main computing alternatives to run genomics workflows in either a
local machine or a distributed cloud environment, but also the main software
technologies available to develop genomic processing services. Since an
analysis encompassing only the computing aspects would provide only a partial
view of the issues for deploying GaaS system, we present also the main
networking technologies that are available to efficiently support a GaaS
solution. We first focus on existing service platforms, and analyze them in
terms of service features, such as scalability, flexibility, and efficiency.
Then, we present a taxonomy for both wide area and datacenter network
technologies that may fit the GaaS requirements. It emerges that
virtualization, both in computing and networking, is the key for a successful
large-scale exploitation of genomic data, by pushing ahead the adoption of the
GaaS paradigm. Finally, the paper illustrates a short and long-term vision on
future research challenges in the field.


  Classical approaches to analyzing dynamical systems, including bifurcation
analysis, can provide invaluable insights into underlying structure of a
mathematical model, and the spectrum of all possible dynamical behaviors.
However, these models frequently fail to take into account population
heterogeneity, which, while critically important to understanding and
predicting the behavior of any evolving system, is a common simplification that
is made in analysis of many mathematical models of ecological systems. Attempts
to include population heterogeneity frequently result in expanding system
dimensionality, effectively preventing qualitative analysis. Reduction Theorem,
or Hidden keystone variable (HKV) method, allows incorporating population
heterogeneity while still permitting the use of previously existing classical
bifurcation analysis. A combination of these methods allows visualization of
evolutionary trajectories and making meaningful predictions about dynamics over
time of evolving populations. Here, we discuss three examples of combination of
these methods to augment understanding of evolving ecological systems. We
demonstrate what new meaningful questions can be asked through this approach,
and propose that the large existing literature of fully analyzed models can
reveal new and meaningful dynamical behaviors with the application of the
HKV-method, if the right questions are asked.


  Configurational arrangement of network architecture and interaction character
of individuals are two most influential factors on the mechanisms underlying
the evolutionary outcome of cooperation, which is explained by the
well-established framework of evolutionary game theory. In the current study,
not only qualitatively but also quantitatively, we measure Moran-evolution of
cooperation to support an analytical agreement based on the consequences of the
replicator equation in a finite population. The validity of the measurement has
been double-checked in the well-mixed network by the Langevin stochastic
differential equation and the Gillespie-algorithmic version of Moran-evolution,
while in a structured network, the measurement of accuracy is verified by the
standard numerical simulation. Considering the Birth-Death and Death-Birth
updating rules through diffusion of individuals, the investigation is carried
out in the wide range of game environments those relate to the various social
dilemmas where we are able to draw a new rigorous mathematical track to tackle
the heterogeneity of complex networks. The set of modified criteria reveals the
exact fact about the emergence and maintenance of cooperation in the structured
population. We find that in general, nature promotes the environment of
coexistent traits.


  Human pathogens transmitted through environmental pathways are subject to
stress and pressures outside of the host. These pressures may cause pathogen
pathovars to diverge in their environmental persistence and their infectivity
on an evolutionary time-scale. On a shorter time-scale, a single-genotype
pathogen population may display wide variation in persistence times and exhibit
biphasic decay. Using an infectious disease transmission modeling framework, we
demonstrate in both cases that fitness-preserving trade-offs have implications
for the dynamics of associated epidemics: less infectious, more persistent
pathogens cause epidemics to progress more slowly than more infectious, less
persistent (labile) pathogens, even when the overall risk is the same. Using
identifiability analysis, we show that the usual disease surveillance data does
not sufficiently inform these underlying pathogen population dynamics, even
with basic environmental monitoring. These results suggest directions for
future microbial research and environmental monitoring. In particular,
determining the relative infectivity of persistent pathogen subpopulations and
the rates of phenotypic conversion will help ascertain how much disease risk is
associated with the long tails of biphasic decay. Alternatively, risk can be
indirectly ascertained by developing methods to separately monitor labile and
persistent subpopulations. A better understanding of persistence--infectivity
trade-offs and associated dynamics can improve risk assessment and disease
control strategies.


  On 2018-01-17 two electron crystallography structures (with PDB entries 6AXZ,
6BTK) on a prion protofibril of bank vole PrP(168-176) (a segment in the PrP
$\beta$2-$\alpha$2 loop) were released into the PDB Bank. The paper published
by [Nat Struct Mol Biol 25(2):131-134 (2018)] reports some polar clasps for
these two crystal structures, and "an intersheet hydrogen bond between Tyr169
and the backbone carbonyl of Asn171 on an opposing strand." - this hydrogen
bond is not between the neighbouring Chain B and Chain A directly. In addition,
by revisiting the polar clasps, we found another two hydrogen bonds
(B.Asn171@H-A.Gln172@OE1, B.Tyr169@OH-A.Gln172@N) between the strand A of one
sheet and the opposing strand B of the mating sheet. For the neighbouring two
single $\beta$-sheets AB, the two new hydrogen bonds are completely different
from the experimental one (an intersheet hydrogen bond between Tyr169 and the
backbone carbonyl of Asn171 on an opposing strand) in [Nat Struct Mol Biol
25(2):131-134 (2018)].


  The relationship between the dynamics of a community and its constituent
pairwise interactions is a fundamental problem in ecology. Higher-order
ecological effects beyond pairwise interactions may be key to complex
ecosystems, but mechanisms to produce these effects remain poorly understood.
Here we show that higher-order effects can arise from variation in multiple
microbial growth traits, such as lag times and growth rates, on a single
limiting resource with no other interactions. These effects produce a range of
ecological phenomena: an unlimited number of strains can exhibit multistability
and neutral coexistence, potentially with a single keystone strain; strains
that coexist in pairs do not coexist all together; and the champion of all
pairwise competitions may not dominate in a mixed community. Since variation in
multiple growth traits is ubiquitous in microbial populations due to pleiotropy
and non-genetic variation, our results indicate these higher-order effects may
also be widespread, especially in laboratory ecology and evolution experiments.


  Han Chinese experienced substantial population migrations and admixture in
history, yet little is known about the evolutionary process of Chinese
dialects. Here, we used phylogenetic approaches and admixture inference to
explicitly decompose the underlying structure of the diversity of Chinese
dialects, based on the total phoneme inventories of 140 dialect samples from
seven traditional dialect groups: Mandarin, Wu, Xiang, Gan, Hakka, Min and Yue.
We found a north-south gradient of phonemic differences in Chinese dialects
induced from historical population migrations. We also quantified extensive
horizontal language transfers among these dialects, corresponding to the
complicated socio-genetic history in China. We finally identified that the
middle latitude dialects of Xiang, Gan and Hakka were formed by admixture with
other four dialects. Accordingly, the middle-latitude areas in China were a
linguistic melting pot of northern and southern Han populations. Our study
provides a detailed phylogenetic and historical context against family-tree
model in China.


  Here we address the challenge of profiling causal properties and tracking the
transformation of chemical compounds from an algorithmic perspective. We
explore the potential of applying a computational interventional calculus based
on the principles of algorithmic probability to chemical structure networks. We
profile the sensitivity of the elements and covalent bonds in a chemical
structure network algorithmically, asking whether reprogrammability affords
information about thermodynamic and chemical processes involved in the
transformation of different compound classes. We arrive at numerical results
suggesting a correspondence between some physical, structural and functional
properties. Our methods are capable of separating chemical classes that reflect
functional and natural differences without considering any information about
atomic and molecular properties. We conclude that these methods, with their
links to chemoinformatics via algorithmic, probability hold promise for future
research.


  Ring oscillators are biochemical circuits consisting of a ring of
interactions capable of sustained oscillations. The non-linear interactions
between genes hinder the analytical insight into their function, usually
requiring computational exploration. Here we show that, despite the apparent
complexity, the stability of the unique steady state in an incoherent feedback
ring depends only on the degradation rates and a single parameter summarizing
the feedback of the circuit. Concretely, we show that the range of regulatory
parameters that yield oscillatory behaviour, is maximized when the degradation
rates are equal. Strikingly, this results holds independently of the regulatory
functions used or number of genes. We also derive properties of the
oscillations as a function of the degradation rates and number of nodes forming
the ring. Finally, we explore the role of mRNA dynamics by applying the generic
results to the specific case with two naturally different degradation time
scales


  The ecologist H. T. Odum introduced a principle of physics, called Maximum
Empower, in order to explain self-organization in a system (e.g. physical,
biological, social, economical, mathematical, ...). The concept of empower
relies on emergy, which is a second notion introduced by Odum for comparing
energy systems on the same basis. The roots of these notions trace back to the
50's (with the work of H. T. Odum and R. C. Pinkerton) and is becoming now an
important sustainability indicator in the ecologist community. In 2012, Le
Corre and Truffet developed a recursive method, based on max-plus algebra, to
compute emergy of a system. Recently, using this max-plus algebra approach, it
has been shown that the Maximum Empower Principle can be formalized as a new
combinatorial optimization problem (called the Maximum Empower Problem).
  In this paper we show that the Maximum Empower Problem can be solved by
finding a maximum weighted clique in a cograph, which leads to an
exponential-time algorithm in the worst-case. We also provide a polynomial-time
algorithm when there is no cycle in the graph modeling the system. Finally, we
prove that the Maximum Empower Problem is #P-hard in the general case, i.e. it
is as hard as computing the permanent of a matrix.


  In this work, we improve a previous minimalistic tree-grass savanna model by
taking into account water availability, in addition to fire, since both factors
are known to be important for shaping savanna physiognomies along a climatic
gradient. As in our previous models, we consider two nonlinear functions of
grass and tree biomasses to respectively take into account grass-fire
feedbacks, and the response of trees to fire of a given intensity. The novelty
is that rainfall is taken into account in the tree and grass growth functions
and in the biomass carrying capacities. Then, we derive a qualitative analysis
of the ODE model, showing existence of equilibria, and studying their stability
conditions. We also construct a two dimension bifurcation diagram based on
rainfall and fire frequency. This led to summarize different scenarios for the
model including multi-stabilities that are proven possible. Next, to bring more
realism in the model, pulsed fire events are modelled as part of an IDE
(Impulsive differential Equations) system analogous to the ODE system.
Numerical simulations are provided and we discuss some important ecological
outcomes that our ODE and IDE models are able to predict. Notably, the
expansion of forest into tree-poor physiognomies (grassland and savanna) is
systematically predicted when fire return period increases, especially in mesic
and humid climatic areas.


  The Gene Ontology aims to define the universe of functions known for gene
products, at the molecular, cellular and organism levels. While the ontology is
designed to cover all aspects of biology in a "species independent manner", the
fact remains that many if not most biological functions are restricted in their
taxonomic range. This is simply because functions evolve, i.e. like other
biological characteristics they are gained and lost over evolutionary time.
Here we introduce a general method of representing the evolutionary gain and
loss of biological functions within the Gene Ontology. We then apply a variety
of techniques, including manual curation, logical reasoning over the ontology
structure, and previously published "taxon constraints" to assign evolutionary
gain and loss events to the majority of terms in the GO. These gain and loss
events now almost triple the number of terms with taxon constraints, and
currently cover a total of 76% of GO terms, including 40% of molecular function
terms, 78% of cellular component terms, and 89% of biological process terms.
  Database URL: GOTaxon is freely available at
https://github.com/haimingt/GOTaxonConstraint


  Gene gains and losses have shaped the gene repertoire of species since the
universal last common ancestor to species today. Genes in extant species were
gained at different historical times via de novo creation of new genes,
duplication of existing genes or transfer from genes of another species (HGT),
and get lost gradually. With the increasing number of sequenced genomes, some
comparative analyses have been done to quantify the evolutionary history of
gene gains and losses in restricted lineages like vertebrates, insects, fungi,
plants and so on. Here, we have constructed and analyzed over 10,000 gene
family trees to reconstruct the gene content of ancestral genomes at an
unprecedented scale, covering hundreds of genomes across all domains of life.
This is the most comprehensive genome-wide analysis of all events in gene
evolutionary histories. We find that our results are largely consistent with
earlier, less complete comparative studies on specific lineages such as the
vertebrates, but find significant differences especially in recent evolutionary
histories. We find that the rate of gene gain varies widely among branches of
the species tree, and find that some periods of rapid gene duplication are
associated with great extinctions in geological history.


  Genealogical networks, also known as family trees or population pedigrees,
are commonly studied by genealogists wanting to know about their ancestry, but
they also provide a valuable resource for disciplines such as digital
demography, genetics, and computational social science. These networks are
typically constructed by hand through a very time-consuming process, which
requires comparing large numbers of historical records manually. We develop
computational methods for automatically inferring large-scale genealogical
networks. A comparison with human-constructed networks attests to the accuracy
of the proposed methods. To demonstrate the applicability of the inferred
large-scale genealogical networks, we present a longitudinal analysis on the
mating patterns observed in a network. This analysis shows a consistent
tendency of people choosing a spouse with a similar socioeconomic status, a
phenomenon known as assortative mating. Interestingly, we do not observe this
tendency to consistently decrease (nor increase) over our study period of 150
years.


  Electrical forces are the background of all the interactions occurring in
biochemical systems. From here and by using a combination of ab-initio and
ad-hoc models, we introduce the first description of electric field profiles
with intrabond resolution to support a characterization of single bond forces
attending to its electrical origin. This fundamental issue has eluded a
physical description so far. Our method is applied to describe hydrogen bonds
(HB) in DNA base pairs. Numerical results reveal that base pairs in DNA could
be equivalent considering HB strength contributions, which challenges previous
interpretations of thermodynamic properties of DNA based on the assumption that
Adenine/Thymine pairs are weaker than Guanine/Cytosine pairs due to the sole
difference in the number of HB. Thus, our methodology provides solid
foundations to support the development of extended models intended to go deeper
into the molecular mechanisms of DNA functioning.


  An explosion of high-throughput DNA sequencing in the past decade has led to
a surge of interest in population-scale inference with whole-genome data.
Recent work in population genetics has centered on designing inference methods
for relatively simple model classes, and few scalable general-purpose inference
techniques exist for more realistic, complex models. To achieve this, two
inferential challenges need to be addressed: (1) population data are
exchangeable, calling for methods that efficiently exploit the symmetries of
the data, and (2) computing likelihoods is intractable as it requires
integrating over a set of correlated, extremely high-dimensional latent
variables. These challenges are traditionally tackled by likelihood-free
methods that use scientific simulators to generate datasets and reduce them to
hand-designed, permutation-invariant summary statistics, often leading to
inaccurate inference. In this work, we develop an exchangeable neural network
that performs summary statistic-free, likelihood-free inference. Our framework
can be applied in a black-box fashion across a variety of simulation-based
tasks, both within and outside biology. We demonstrate the power of our
approach on the recombination hotspot testing problem, outperforming the
state-of-the-art.


  Understanding how transcription factors (TFs) regulate mammalian gene
expression in space and time is a central topic in biology. To activate a gene,
a TF has first to diffuse in the available space of the nucleus until it
reaches a target DNA sequence or protein (target site). This eventually results
in the recruitment of the whole transcriptional machinery.
  All these processes take place in the mammalian nucleoplasm, a highly
organized and dynamic environment, in which some complexes transiently assemble
and break apart, whereas others appear more stable. This diversity of dynamic
behaviors arises from the number of biomolecules that make up the nucleoplasm
and their pairwise interactions. Indeed, interactions energies that span
several orders of magnitude, from covalent bounds to transient and dynamic
interactions can shape nuclear landscapes. Thus, the nuclear environment
determines how frequently and how fast a TF contacts its target site, and
indirectly gene expression. How exactly transient interactions are involved in
the regulation of TF diffusion is unclear, but are reflected by live cell
imaging techniques such as fluorescence correlation spectroscopy, fluorescence
recovery after photobleaching or single-particle tracking. Overall, the
macroscopic result of these microscopic interactions is almost always anomalous
diffusion, a phenomenon widely studied and modeled.
  Here, we review the connections between the anomalous diffusion of a TF and
the microscopic organization of the nucleus, including recently described
topologically associated domains and dynamic phase-separated compartments. We
propose that anomalous diffusion found in single particle tracking (SPT) data
result from weak and transient interactions with dynamic nuclear substructures,
and that SPT data analysis would benefit form a better description of such
structures.


  The Kinfold and KFOLD programs for RNA folding kinetics implement the
Gillespie algorithm to generate stochastic folding trajectories from an initial
structure s to a target structure t, in which each intermediate secondary
structure is obtained from its predecessor by the addition, removal or shift of
a single base pair. Define MS2 distance between secondary structures s and t to
be the minimum path length to refold s to t, where a move from MS2 is applied
in each step. We describe algorithms to compute the shortest MS2 folding
trajectory between any two given RNA secondary structures. These algorithms
include an optimal integer programming (IP) algorithm, an accurate and
efficient near-optimal algorithm, a greedy algorithm, a branch-and-bound
algorithm, and an optimal algorithm if one allows intermediate structures to
contain pseudoknots. Our optimal IP [resp. near-optimal IP] algorithm maximizes
[resp. approximately maximizes] the number of shifts and minimizes [resp.
approximately minimizes] the number of base pair additions and removals by
applying integer programming to (essentially) solve the minimum feedback vertex
set (FVS) problem for the RNA conflict digraph, then applies topological sort
to tether subtrajectories into the final optimal folding trajectory. We prove
NP-hardness of the problem to determine the minimum barrier energy over all
possible MS2 folding pathways, and conjecture that computing the MS2 distance
between arbitrary secondary structures is NP-hard. Since our optimal IP
algorithm relies on the FVS, known to be NP-complete for arbitrary digraphs, we
compare the family of RNA conflict digraphs with the following classes of
digraphs (planar, reducible flow graph, Eulerian, and tournament) for which FVS
is known to be either polynomial time computable or NP-hard. Source code
available at http://bioinformatics.bc.edu/clotelab/MS2distance/.


  B cells develop high affinity receptors during the course of affinity
maturation, a cyclic process of mutation and selection. At the end of affinity
maturation, a number of cells sharing the same ancestor (i.e. in the same
"clonal family") are released from the germinal center, their amino acid
frequency profile reflects the allowed and disallowed substitutions at each
position. These clonal-family-specific frequency profiles, called "substitution
profiles", are useful for studying the course of affinity maturation as well as
for antibody engineering purposes. However, most often only a single sequence
is recovered from each clonal family in a sequencing experiment, making it
impossible to construct a clonal-family-specific substitution profile. Given
the public release of many high-quality large B cell receptor datasets, one may
ask whether it is possible to use such data in a prediction model for
clonal-family-specific substitution profiles. In this paper, we present the
method "Substitution Profiles Using Related Families" (SPURF), a penalized
tensor regression framework that integrates information from a rich assemblage
of datasets to predict the clonal-family-specific substitution profile for any
single input sequence. Using this framework, we show that substitution profiles
from similar clonal families can be leveraged together with simulated
substitution profiles and germline gene sequence information to improve
prediction. We fit this model on a large public dataset and validate the
robustness of our approach on an external dataset. Furthermore, we provide a
command-line tool in an open-source software package
(https://github.com/krdav/SPURF) implementing these ideas and providing easy
prediction using our pre-fit models.


  Each living species carries a complex DNA sequence that determines their
unique features and functionalities. It is generally assumed that life started
from a random pool of oligonucleotides sequences, generated by a prebiotic
polymerization of nucleotides. The mechanism that initially facilitated the
emergence of sequences that code for the function of the first species from
such a random pool of sequences remains unknown. It is a central problem of the
origin of life. An interesting option would be a self-selection mechanism by
spontaneous symmetry breaking. Initial concentration fluctuations of specific
sequence motifs would have been amplified and outcompeted less abundant
sequences, enhancing the signal to noise to replicate and select functional
sequences. Here, we demonstrate with experimental and theoretical findings that
templated ligation would provide such a self-selection. In templated ligation,
two adjacent single sequences strands are chemically joined when a third
complementary strand sequence brought them in close proximity. This simple
mechanism was a likely side-product of a prebiotic polymerization chemistry
once the strands reach the length to form double stranded species. As shown
here, the ligation gave rise to a nonlinear replication process by the
cooperative ligation of matching sequences which self-promoted their own
elongation. This led to a cascade of enhanced template binding and faster
ligation reactions. A requirement was the reshuffling of the strands by thermal
cycling, enabled for example by microscale convection. Assuming that templated
ligation was driven by the same chemical mechanism that generated prebiotic
polymerization of oligonucleotides, the mechanism could function as a missing
link between polymerization and the self-stabilized replication, offering a
pathway to the autonomous emergence of Darwinian evolution for the origin of
life.


  Inferential methods can be used to integrate experimental informations and
molecular simulations. The maximum entropy principle provides a framework for
using equilibrium experimental data and it has been shown that replica-averaged
simulations, restrained using a static potential, are a practical and powerful
implementation of such principle. Here we show that replica-averaged
simulations restrained using a time-dependent potential are equivalent to the
principle of maximum caliber, the dynamic version of the principle of maximum
entropy, and thus may allow to integrate time-resolved data in molecular
dynamics simulations. We provide an analytical proof of the equivalence as well
as a computational validation making use of simple models and synthetic data.
Some limitations and possible solutions are also discussed.


  Exploiting others is beneficial individually but it could also be detrimental
globally. The reverse is also true: a higher cooperation level may change the
environment in a way that is beneficial for all competitors. To explore the
possible consequence of this feedback we consider a coevolutionary model where
the local cooperation level determines the payoff values of the applied
prisoner's dilemma game. We observe that the coevolutionary rule provides a
significantly higher cooperation level comparing to the traditional setup
independently of the topology of the applied interaction graph. Interestingly,
this cooperation supporting mechanism offers lonely defectors a high surviving
chance for a long period hence the relaxation to the final cooperating state
happens logarithmically slow. As a consequence, the extension of the
traditional evolutionary game by considering interactions with the environment
provides a good opportunity for cooperators, but their reward may arrive with
some delay.


  During the last decade coarse-grained nucleotide models have emerged that
allow us to DNA and RNA on unprecedented time and length scales. Among them is
oxDNA, a coarse-grained, sequence-specific model that captures the
hybridisation transition of DNA and many structural properties of single- and
double-stranded DNA. oxDNA was previously only available as standalone
software, but has now been implemented into the popular LAMMPS molecular
dynamics code. This article describes the new implementation and analyses its
parallel performance. Practical applications are presented that focus on
single-stranded DNA, an area of research which has been so far
under-investigated. The LAMMPS implementation of oxDNA lowers the entry barrier
for using the oxDNA model significantly, facilitates future code development
and interfacing with existing LAMMPS functionality as well as other
coarse-grained and atomistic DNA models.


  Dynamic patterning of specific proteins is essential for the spatiotemporal
regulation of many important intracellular processes in procaryotes,
eucaryotes, and multicellular organisms. The emergence of patterns generated by
interactions of diffusing proteins is a paradigmatic example for
self-organization. In this article we review quantitative models for
intracellular Min protein patterns in E. coli, Cdc42 polarization in S.
cerevisiae, and the bipolar PAR protein patterns found in C. elegans. By
analyzing the molecular processes driving these systems we derive a theoretical
perspective on general principles underlying self-organized pattern formation.
We argue that intracellular pattern formation is not captured by concepts such
as "activators"', "inhibitors", or "substrate-depletion". Instead,
intracellular pattern formation is based on the redistribution of proteins by
cytosolic diffusion, and the cycling of proteins between distinct
conformational states. Therefore, mass-conserving reaction-diffusion equations
provide the most appropriate framework to study intracellular pattern
formation. We conclude that directed transport, e.g. cytosolic diffusion along
an actively maintained cytosolic gradient, is the key process underlying
pattern formation. Thus the basic principle of self-organization is the
establishment and maintenance of directed transport by intracellular protein
dynamics.


  In this paper, I model and study the process of natural selection between all
possible mixed strategies in classical two-player two-strategy games. I derive
and solve an equation that is a natural generalization of the Taylor-Jonker
replicator equation that describes dynamics of pure strategy frequencies. I
then investigate the evolution of not only frequencies of pure strategies but
also of total distribution of mixed strategies. I show that the process of
natural selection of strategies for all games obeys the dynamical Principle of
minimum of information gain.


  The famous 2R hypothesis was first proposed by Susumu Ohno in 1970. It states
that the two whole genome duplications had shaped the genome of early
vertebrates. The most convincing evidence for 2R hypothesis comes from the 4:1
ratio chromosomal regions that have preserved both gene content and order in
vertebrates compared with closely related. However, due to the shortage of such
strict evidence, the 2R hypothesis is still under debates.
  Here, we present a combined perspective of phylogenetic and genomic homology
to revisit the hypothesis of 2R whole genome duplications. Ancestral vertebrate
genomes as well as ancient duplication events were created from 17 extant
vertebrate species. Extant descendants from the duplication events at early
vertebrates were extracted and reorganized to partial genomes. We then examined
the gene order based synteny, and projected back to phylogenetic gene trees for
examination of synteny evidence of the reconstructed early vertebrate genes. We
identified 7877 ancestral genes that were created from 3026 duplication events
at early vertebrates, and more than 50% of the duplication events show synteny
evidence. Thus, our reconstructions provide very strong evidence for the 2R
hypothesis.
  We also reconstructed the genome of early vertebrates, and built a model of
the gene gains and losses in early vertebrates. We estimated that there were
about 12,000 genes in early vertebrates before 2R, and the probability of a
random gene get lost after the first round of whole genome duplication is
around 0.45, and the probability of a random gene get lost after the second
round of whole genome duplication is around 0.55.
  This research provides convincing evidence for the 2R hypothesis, and may
provide further insights in vertebral evolution.
  Data availability: https://github.com/haimingt/Ohnologs-and-2R-WGD


  The Gene Ontology Consortium launched the GO-PAINT project (Phylogenetic
Annotation and INference Tool) 9 years ago and is currently being used in the
GO Reference Genome Annotation Project to support inference of GO function
terms (molecular function, cellular component and biological process) by
homology. PAINT uses a phylogenetic model to infer gene function by homology, a
process that requires manual curation of experienced biocurators. Tremendous
amount of time and efforts have been spent on the GO-PAINT project yielding
more than 4000 fully annotated phylogenetic families with more than 170,000
annotations. These preliminary data have thus enabled potential algorithmic
representation and automatic solvation of the additional 9000 unannoated
phylogenetic families. Here we present an automated pipeline for phylogenetic
annotation and inference, which simulates the standard annotation procedures of
curators and models the curators' decisions during the manual curation process.
The pipeline has been built into the newest version of PAINT software available
at http://www.pantherdb.org/downloads/index.jsp. The standalone automation
pipeline and datasets are available at
https://github.com/haimingt/GO-PAINT-automation


  Summary: TreeGrafter is a new software tool for annotating protein sequences
using annotated phylogenetic trees. Cur-rently, the tool provides annotations
to Gene Ontology terms, and PANTHER protein class, family and subfamily. The
ap-proach is generalizable to any annotations that have been made to internal
nodes of a reference phylogenetic tree. Tree-Grafter takes each input query
protein sequence, finds the best matching homologous family in a library of
pre-calculated, pre-annotated gene trees, and then grafts it to the best
location in the tree. It then annotates the sequence by propagating
anno-tations from its ancestral nodes in the reference tree. We show that
TreeGrafter outperforms subfamily HMM scoring for cor-rectly assigning
subfamily membership, and that it produces highly specific annotations of GO
terms based on annotated reference phylogenetic trees. This method will be
further inte-grated into InterProScan, enabling an even broader user
com-munity.
  Availability: TreeGrafter is freely available on the web at
https://github.com/haimingt/TreeGrafting.


  The stationary distribution of the diffusion limit of the 2-island, 2-allele
Wright-Fisher with small but otherwise arbitrary mutation and migration rates
is investigated. Following a method developed by Burden and Tang (2016, 2017)
for approximating the forward Kolmogorov equation, the stationary distribution
is obtained to leading order as a set of line densities on the edges of the
sample space, corresponding to states for which one island is bi-allelic and
the other island is non-segregating, and a set of point masses at the corners
of the sample space, corresponding to states for which both islands are
simultaneously non-segregating. Analytic results for the corner probabilities
and line densities are verified independently using the backward generator and
for the corner probabilities using the coalescent.


  The collective behavior of active semiflexible filaments is studied with a
model of tangentially driven self-propelled worm-like chains. The combination
of excluded-volume interactions and self-propulsion leads to several distinct
dynamic phases as a function of bending rigidity, activity, and aspect ratio of
individual filaments. We consider first the case of intermediate filament
density. For high-aspect-ratio filaments, we identify a transition with
increasing propulsion from a state of free-swimming filaments to a state of
spiraled filaments with nearly frozen translational motion. For lower aspect
ratios, this gas-of-spirals phase is suppressed with growing density due to
filament collisions; instead, filaments form clusters similar to self-propelled
rods, as activity increases. Finite bending rigidity strongly effects the
dynamics and phase behavior. Flexible filaments form small and transient
clusters, while stiffer filaments organize into giant clusters, similarly as
self-propelled rods, but with a reentrant phase behavior from giant to smaller
clusters as activity becomes large enough to bend the filaments. For high
filament densities, we identify a nearly frozen jamming state at low
activities, a nematic laning state at intermediate activities, and an
active-turbulence state at high activities. The latter state is characterized
by a power-law decay of the energy spectrum as a function of wave number. The
resulting phase diagrams encapsulate tunable non-equilibrium steady states that
can be used in the organization of living matter.


  The study of genome rearrangement has many flavours, but they all are somehow
tied to edit distances on variations of a multi-graph called the breakpoint
graph. We study a weighted 2-break distance on Eulerian 2-edge-colored
multi-graphs, which generalizes weighted versions of several Double Cut and
Join problems, including those on genomes with unequal gene content. We affirm
the connection between cycle decompositions and edit scenarios first discovered
with the Sorting By Reversals problem. Using this we show that the problem of
finding a parsimonious scenario of minimum cost on an Eulerian 2-edge-colored
multi-graph - with a general cost function for 2-breaks - can be solved by
decomposing the problem into independent instances on simple alternating
cycles. For breakpoint graphs, and a more constrained cost function, based on
coloring the vertices, we give a polynomial-time algorithm for finding a
parsimonious 2-break scenario of minimum cost, while showing that finding a
non-parsimonious 2-break scenario of minimum cost is NP-Hard.


  Summary
  1. Infectious disease outbreaks in plants threaten ecosystems, agricultural
crops and food trade. Currently, several fungal diseases are affecting forests
worldwide, posing a major risk to tree species, habitats and consequently
ecosystem decay. Prediction and control of disease spread are difficult, mainly
due to the complexity of the interaction between individual components
involved.
  2. In this work, we introduce a lattice-based epidemic model coupled with a
stochastic process that mimics, in a very simplified way, the interaction
between the hosts and pathogen. We studied the disease spread by measuring the
propagation velocity of the pathogen on the susceptible hosts. Quantitative
results indicate the occurrence of a critical transition between two stable
phases: local confinement and an extended epiphytotic outbreak that depends on
the density of the susceptible individuals.
  3. Quantitative predictions of epiphytotics are performed using the framework
early-warning indicators for impending regime shifts, widely applied on
dynamical systems. These signals forecast successfully the outcome of the
critical shift between the two stable phases before the system enters the
epiphytotic regime.
  4. Synthesis: Our study demonstrates that early-warning indicators could be
useful for the prediction of forest disease epidemics through mathematical and
computational models suited to more specific pathogen-host-environmental
interactions.


  Classical simulations of protein flexibility remain computationally
expensive, especially for large proteins. A few years ago, we developed a fast
method for predicting protein structure fluctuations that uses a single protein
model as the input. The method has been made available as the CABS-flex web
server and applied in numerous studies of protein structure-function
relationships. Here, we present a major update of the CABS-flex web server to
version 2.0. The new features include: extension of the method to significantly
larger and multimeric proteins, customizable distance restraints and simulation
parameters, contact maps and a new, enhanced web server interface. CABS-flex
2.0 is freely available at http://biocomp.chem.uw.edu.pl/CABSflex2


  A standard belief on emerging collective behavior is that it emerges from
simple individual rules. Most of the mathematical research on such collective
behavior starts from imperative individual rules, like always go to the center.
But how could an (optimal) individual rule emerge during a short period within
the group lifetime, especially if communication is not available. We argue that
such rules can actually emerge in a group in a short span of time via
collective (multi-agent) reinforcement learning, i.e learning via rewards and
punishments. We consider the gathering problem: several agents (social animals,
swarming robots...) must gather around a same position, which is not determined
in advance. They must do so without communication on their planned decision,
just by looking at the position of other agents. We present the first
experimental evidence that a gathering behavior can be learned without
communication in a partially observable environment. The learned behavior has
the same properties as a self-stabilizing distributed algorithm, as processes
can gather from any initial state (and thus tolerate any transient failure).
Besides, we show that it is possible to tolerate the brutal loss of up to 90\%
of agents without significant impact on the behavior.


  The stationary distribution of a sample taken from a Wright-Fisher diffusion
with general small mutation rates is found using a coalescent approach. The
approximation is equivalent to having at most one mutation in the coalescent
tree to the first order in the rates. The sample probabilities characterize an
approximation for the stationary distribution from the Wright-Fisher diffusion.
The approach is different from Burden and Tang (2016,2017) who use a
probability flux argument to obtain the same results from a forward diffusion
generator equation. The solution has interest because the solution is not known
when rates are not small. An analogous solution is found for the configuration
of alleles in a general exchangeable binary coalescent tree. In particular an
explicit solution is found for a pure birth process tree when individuals
reproduce at rate lambda.


  Timing is of fundamental importance in biology and our life. Borrowing ideas
from mechanism, we map our clock signals onto a gear system, in pursuit of
better depiction of a clock signal implemented with chemical reaction networks
(CRNs). On a chassis of gear theory, more quantitative descriptions are offered
for our method. Inspired by gears, our work to synthesize a tunable clock
signal could be divided into two parts. Part I, this paper, mainly focuses on
the implementation of clock signals with three types of duty cycles, namely
$1/2$, $1/N$ ($N > 2$), and $M/N$. Part II devotes itself in addressing
frequency alteration issues of clock signals. \textcolor{black}{Guaranteed by
existing literature, the experimental chassis can be taken care of by DNA
strand displacement reactions, which lay a solid foundation for the physical
implementation of nearly arbitrary CRNs.


  On a chassis of gear model, we have offered a quantitative description for
our method to synthesize a chemical clock signal with various duty cycles in
Part I. As Part II of the study, this paper devotes itself in proposing a
design methodology to handle frequency alteration issues for the chemical
clock, including both frequency division and frequency multiplication. Several
interesting examples are provided for a better explanation of our contribution.
All the simulation results verify and validate the correctness and efficiency
of our proposal.


  Many examples of cooperation exist in biology. In chemical systems however,
which can sometimes be quite complex, we do not appear to observe intricate
cooperative interactions. A key question for the origin of life, is then how
can molecular cooperation first arise in an abiotic system prior to the
emergence of biological replication. We postulate that selection at the
molecular level is a driving force behind the complexification of chemical
systems, particularly during the origins of life. In the theory of multilevel
selection the two selective forces are: within-group and between-group, where
the former tends to favor "selfish" replication of individuals and the latter
favor cooperation between individuals enhancing the replication of the group as
a whole. These forces can be quantified using the Price equation, which is a
standard tool used in evolutionary biology to quantify evolutionary change. Our
central claim is that replication and heredity in chemical systems are subject
to selection, and quantifiable using the multilevel Price equation. We
demonstrate this using the Graded Autocatalysis Replication Domain computer
model, describing simple protocell composed out of molecules and its
replication, which respectively analogue to the group and the individuals. In
contrast to previous treatments of this model, we treat the lipid molecules
themselves as replicating individuals and the protocells they form as groups of
individuals. Our goal is to demonstrate how evolutionary biology tools and
concepts can be applied in chemistry and we suggest that molecular cooperation
may arise as a result of group selection. Further, the biological relation of
parent-progeny is proposed to be analogue to the reactant-product relation in
chemistry, thus allowing for tools from evolutionary biology to be applied to
chemistry and would deepen the connection between chemistry and biology.


  Bakground: With the proliferation of available microarray and high throughput
sequencing experiments in the public domain, the use of meta-analysis methods
increases. In these experiments, where the sample size is often limited,
meta-analysis offers the possibility to considerably enhance the statistical
power and give more accurate results. For those purposes, it combines either
effect sizes or results of single studies in a appropriate manner. R packages
metaMA and metaRNASeq perform meta-analysis on microarray and NGS data,
respectively. They are not interchangeable as they rely on statistical modeling
specific to each technology.
  Results: SMAGEXP (Statistical Meta-Analysis for Gene EXPression) integrates
metaMA and metaRNAseq packages into Galaxy. We aim to propose a unified way to
carry out meta-analysis of gene expression data, while taking care of their
specificities. We have developed this tool suite to analyse microarray data
from Gene Expression Omnibus (GEO) database or custom data from affymetrix
microarrays. These data are then combined to carry out meta-analysis using
metaMA package. SMAGEXP also offers to combine raw read counts from Next
Generation Sequencing (NGS) experiments using DESeq2 and metaRNASeq package. In
both cases, key values, independent from the technology type, are reported to
judge the quality of the meta-analysis. These tools are available on the Galaxy
main tool shed. Source code, help and installation instructions are available
on github.
  Conclusion: The use of Galaxy offers an easy-to-use gene expression
meta-analysis tool suite based on the metaMA and metaRNASeq packages.


  Elastic network models (ENM) and constraint-based, topological rigidity
analysis are two distinct, coarse-grained approaches to study conformational
flexibility of macromolecules. In the two decades since their introduction,
both have contributed significantly to insights into protein molecular
mechanisms and function. However, despite a shared purpose of these approaches,
the topological nature of rigidity analysis, and thereby the absence of motion
modes, has impeded a direct comparison. Here, we present an alternative,
kinematic approach to rigidity analysis, which circumvents these drawbacks. We
introduce a novel protein hydrogen bond network spectral decomposition, which
provides an orthonormal basis for collective motions modulated by non-covalent
interactions, analogous to the eigenspectrum of normal modes, and decomposes
proteins into rigid clusters identical to those from topological rigidity. Our
kinematic flexibility analysis bridges topological rigidity theory and ENM, and
enables a detailed analysis of motion modes obtained from both approaches. Our
analysis reveals that collectivity of protein motions, reported by the Shannon
entropy, is significantly lower for rigidity theory versus normal mode
approaches. Strikingly, kinematic flexibility analysis suggests that the
hydrogen bonding network encodes a protein-fold specific, spatial hierarchy of
motions, which goes nearly undetected in ENM. This hierarchy reveals distinct
motion regimes that rationalize protein stiffness changes observed from
experiment and molecular dynamics simulations. A formal expression for changes
in free energy derived from the spectral decomposition indicates that motions
across nearly 40% of modes obey enthalpy-entropy compensation. Taken together,
our analysis suggests that hydrogen bond networks have evolved to modulate
protein structure and dynamics.


  We confront a network model of human aging and mortality in which nodes
represent health attributes that interact within a scale-free network topology,
with observational data that uses both clinical and laboratory (pre-clinical)
health deficits as network nodes. We find that individual health attributes
exhibit a wide range of mutual information with mortality and that, with a re-
construction of their relative connectivity, higher-ranked nodes are more
informative. Surprisingly, we find a broad and overlapping range of mutual
information of laboratory measures as compared with clinical measures. We
confirm similar behavior between most-connected and least-connected model
nodes, controlled by the nearest-neighbor connectivity. Furthermore, in both
model and observational data, we find that the least-connected (laboratory)
nodes damage earlier than the most-connected (clinical) deficits. A mean-field
theory of our network model captures and explains this phenomenon, which
results from the connectivity of nodes and of their connected neighbors. We
find that other network topologies, including random, small-world, and
assortative scale-free net- works, exhibit qualitatively different behavior.
Our disassortative scale-free network model behaves consistently with our
expanded phenomenology observed in human aging, and so is a useful tool to
explore mechanisms of and to develop new predictive measures for human aging
and mortality.


  Understanding the relationship between protein sequence, function, and
stability is a fundamental problem in biology. While high-throughput methods
have produced large numbers of sequence-function pairs, functional assays do
not distinguish whether mutations directly affect function or are destabilizing
the protein. Here, we introduce a statistical method to infer the underlying
biophysics from a high-throughput binding assay by combining information from
many mutated variants. We fit a thermodynamic model describing the bound,
unbound, and unfolded states to high quality data of protein G domain B1
binding to IgG-Fc. We infer an energy landscape with distinct folding and
binding energies for each substitution providing a detailed view of how
mutations affect binding and stability across the protein. We accurately infer
folding energy of each variant in physical units, validated by independent
data, whereas previous high-throughput methods could only measure indirect
changes in stability. While we assume an additive sequence-energy relationship,
the binding fraction is epistatic due its non-linear relation to energy.
Despite having no epistasis in energy, our model explains much of the observed
epistasis in binding fraction, with the remaining epistasis identifying
conformationally dynamic regions.


  During chronic infection, HIV-1 engages in a rapid coevolutionary arms race
with the host's adaptive immune system. While it is clear that HIV exerts
strong selection on the adaptive immune system, the characteristics of the
somatic evolution that shape the immune response are still unknown. Traditional
population genetics methods fail to distinguish chronic immune response from
healthy repertoire evolution. Here, we infer the evolutionary modes of B-cell
repertoires and identify complex dynamics with a constant production of better
B-cell receptor mutants that compete, maintaining large clonal diversity and
potentially slowing down adaptation. A substantial fraction of mutations that
rise to high frequencies in pathogen engaging CDRs of B-cell receptors (BCRs)
are beneficial, in contrast to many such changes in structurally relevant
frameworks that are deleterious and circulate by hitchhiking. We identify a
pattern where BCRs in patients who experience larger viral expansions undergo
stronger selection with a rapid turnover of beneficial mutations due to clonal
interference in their CDR3 regions. Using population genetics modeling, we show
that the extinction of these beneficial mutations can be attributed to the rise
of competing beneficial alleles and clonal interference. The picture is of a
dynamic repertoire, where better clones may be outcompeted by new mutants
before they fix.


  We investigate the expected time to extinction in the
susceptible-infectious-susceptible (SIS) model of disease spreading. Rather
than using stochastic simulations, or asymptotic calculations in network
models, we solve the extinction time exactly for all connected graphs with
three to eight vertices. This approach enables us to discover descriptive
relations that would be impossible with stochastic simulations. It also helps
us discovering graphs and configurations of S and I with anomalous behaviors
with respect to disease spreading. We find that for large transmission rates
the extinction time is independent of the configurations, just dependent on the
graph. In this limit, the number of vertices and edges determine the extinction
time very accurately (deviations primarily coming from the fluctuations in
degrees). We find that the rankings of configurations with respect to
extinction times at low and high transmission rates are correlated at low
prevalences and negatively correlated for high prevalences. The most important
structural factor determining this ranking is the degrees of the infectious
vertices.


  Steady state is an essential concept in reaction networks. Its stability
reflects fundamental characteristics of several biological phenomena such as
cellular signal transduction and gene expression. Because biochemical reactions
occur at the cellular level, they are affected by unavoidable fluctuations.
Although several methods have been proposed to detect and analyze the stability
of steady states for deterministic models, these methods cannot be applied to
stochastic reaction networks. In this paper, we propose an algorithm based on
algebraic computations to calculate parameter regions for constrained
steady-state distribution of stochastic reaction networks, in which the means
and variances satisfy some given inequality constraints. To evaluate our
proposed method, we perform computer simulations for three typical chemical
reactions and demonstrate that the results obtained with our method are
consistent with the simulation results.


  The availability of protein is an important factor for the determination of
the size of the mitotic spindle. Involved in spindle-size regulation is
kinesin-8, a molecular motor and microtubule (MT) depolymerase, which is known
to tightly control MT length. Here, we propose and analyze a theoretical model
in which kinesin-induced MT depolymerization competes with spontaneous
polymerization while supplies of both tubulin and kinesin are limited. In
contrast to previous studies where resources were unconstrained, we find that,
for a wide range of concentrations, MT length regulation is bistable. We test
our predictions by conducting in vitro experiments, and find that the bistable
behavior manifests in a bimodal MT length distribution.


  Despite the great advancement of experimental tools and theoretical models, a
quantitative characterization of the microscopic structures of ion aggregates
and its associated water hydrogen-bonding networks still remains a challenging
problem. In this paper, a newly-invented mathematical method called persistent
homology is introduced, for the first time, to quantitatively analyze the
intrinsic topological properties of ion aggregation systems and
hydrogen-bonding networks. Two most distinguishable properties of persistent
homology analysis of assembly systems are as follows. First, it does not
require a predefined bond length to construct the ion or hydrogen network.
Persistent homology results are determined by the morphological structure of
the data only. Second, it can directly measure the size of circles or holes in
ion aggregates and hydrogen-bonding networks. To validate our model, we
consider two well-studied systems, i.e., NaCl and KSCN solutions, generated
from molecular dynamics simulations. They are believed to represent two
morphological types of aggregation, i.e., local clusters and extended ion
network. It has been found that the two aggregation types have distinguishable
topological features and can be characterized by our topological model very
well. For hydrogen-bonding networks, KSCN systems demonstrate much more
dramatic variations in their local circle structures with the concentration
increase. A consistent increase of large-sized local circle structures is
observed and the sizes of these circles become more and more diverse. In
contrast, NaCl systems show no obvious increase of large-sized circles. Instead
a consistent decline of the average size of circle structures is observed and
the sizes of these circles become more and more uniformed with the
concentration increase.


  Single-cell RNA sequencing (scRNA-seq) is a fast growing approach to measure
the genome-wide transcriptome of many individual cells in parallel, but results
in noisy data with many dropout events. Existing methods to learn molecular
signatures from bulk transcriptomic data may therefore not be adapted to
scRNA-seq data, in order to automatically classify individual cells into
predefined classes. We propose a new method called DropLasso to learn a
molecular signature from scRNA-seq data. DropLasso extends the dropout
regularisation technique, popular in neural network training, to esti- mate
sparse linear models. It is well adapted to data corrupted by dropout noise,
such as scRNA-seq data, and we clarify how it relates to elastic net
regularisation. We provide promising results on simulated and real scRNA-seq
data, suggesting that DropLasso may be better adapted than standard regularisa-
tions to infer molecular signatures from scRNA-seq data.


  We use an evolutionary game model to study the interplay between corporate
environmental compliance and enforcement promoted by the policy maker in a
country facing a pollution trap, i.e., a scenario in which the vast majority of
firms do not internalize their pollution negative externality and auditors do
not inspect firms. The game conflict is due to the trade-off in which firms are
better-off when they pollute and are not inspected, while social welfare is
maximized when auditors do not need to inspect socially responsible
corporations that account for pollution in their production decisions regarding
technology used and emission level. Starting with a well-mixed two-population
game model, there is no long-run equilibrium and the shares of polluters and
shirking auditors keep oscillating over time. In contrast, when firms and
auditors are allocated in a spatial network, the game displays a rich dynamics
depending on the inspecting cost. While the oscillatory behaviour is still
possible, there is a set of parameters for which a long run robust equilibrium
is achieved with the country leaving the pollution trap. On the other hand, an
excessively high inspection cost leads to an ineffective auditing process where
the few compliant firms are driven out of the country.


  Species accumulation curves (SAC), i.e. the relationship between species
richness and the number of sampling units in a given community, can be used to
describe diversity patterns while accounting for the well-known
scale-dependence of species richness. Despite their value, the functional form
and the parameters of SAC, as well as their determinants, have barely been
investigated in plant communities, particularly in drylands. We characterized
the SAC of perennial plant communities from 233 dryland ecosystems from six
continents by comparing the fit of major functions (power-law, logarithmic and
Michaelis-Menten). We tested the theoretical prediction that the effects of
aridity and soil pH on SAC are mediated by vegetation attributes such as
evenness, cover, and spatial aggregation. We found that the logarithmic
relationship was the most common functional form, followed by Michaelis-Menten
and power-law. Functional form was mainly determined by evenness while the SAC
parameters (intercept and slope) were largely determined by spatial
aggregation. In addition, aridity decreased small scale richness (intercept of
SAC) but did not affect accumulation rate (slope of the SAC). Our results
highlight the role that attributes such as spatial aggregation and evenness
play as main mediators of the SAC of vegetation in drylands, the Earth's
largest biome.


  The COI mitochondrial gene is present in all animal phyla and in a few
others, and is the leading candidate for species identification through DNA
barcoding. Calculating a generalized form of total correlation on publicly
available data on the gene yields distinctive information-theoretic descriptors
of the phyla represented in the data. Moreover, performing principal component
analysis on standardized versions of these descriptors reveals a strong
correlation between the first principal component and the natural logarithm of
the number of known living species. The descriptors thus constitute clear
information-theoretic signatures of the processes whereby evolution has given
rise to current biodiversity.


  Sensing and reciprocating cellular systems (SARs) are important for the
operation of many biological systems. Production in interferon (IFN) SARs is
achieved through activation of the Jak-Stat pathway, and downstream
upregulation of IFN regulatory factor (IRF)-3 and IFN transcription, but the
role that high and low affinity IFNs play in this process remains unclear. We
present a comparative between a minimal spatio-temporal partial differential
equation (PDE) model and a novel spatio-structural-temporal (SST) model for the
consideration of receptor, binding, and metabolic aspects of SAR behaviour.
Using the SST framework, we simulate single- and multi-cluster paradigms of IFN
communication. Simulations reveal a cyclic process between the binding of IFN
to the receptor, and the consequent increase in metabolism, decreasing the
propensity for binding due to the internal feed-back mechanism. One observes
the effect of heterogeneity between cellular clusters, allowing them to
individualise and increase local production, and within clusters, where we
observe `sub popular quiescence'; a process whereby intra-cluster
subpopulations reduce their binding and metabolism such that other such
subpopulations may augment their production. Finally, we observe the ability
for low affinity IFN to communicate a long range signal, where high affinity
cannot, and the breakdown of this relationship through the introduction of cell
motility. Biological systems may utilise cell motility where environments are
unrestrictive and may use fixed system, with low affinity communication, where
a localised response is desirable.


  Considering all the PDB annotated allosteric proteins (from ASD - AlloSteric
Database) belonging to four different classes (kinases, nuclear receptors,
peptidases and transcription factors), this work has attempted to decipher
certain consistent patterns present in the residues constituting the allosteric
communication sub-system (ACSS). The thermal fluctuations of hydrophobic
residues in ACSSs were found to be significantly higher than those present in
the non-ACSS part of the same proteins, while polar residues showed the
opposite trend.
  The basic residues and hydroxyl residues were found to be slightly more
predominant than the acidic residues and amide residues in ACSSs, hydrophobic
residues were found extremely frequently in kinase ACSSs. Despite having
different sequences and different lengths of ACSS, they were found to be
structurally quite similar to each other - suggesting a preferred structural
template for communication. ACSS structures recorded low RMSD and high Akaike
Information Criterion(AIC) scores among themselves. While the ACSS networks for
all the groups of allosteric proteins showed low degree centrality and
closeness centrality, the betweenness centrality magnitudes revealed nonuniform
behavior. Though cliques and communities could be identified within the ACSS,
maximal-common-subgraph considering all the ACSS could not be generated,
primarily due to the diversity in the dataset. Barring one particular case, the
entire ACSS for any class of allosteric proteins did not demonstrate "small
world" behavior, though the sub-graphs of the ACSSs, in certain cases, were
found to form small-world networks.


  Trait variation and similarity among coexisting species can provide a window
into the mechanisms that maintain their coexistence. Recent theoretical
explorations suggest that competitive interactions will lead to groups, or
clusters, of species with similar traits. However, theoretical predictions
typically assume complete knowledge of the map between competition and measured
traits. These assumptions limit the plausible application of these patterns for
inferring competitive interactions in nature. Here we relax these restrictions
and find that the clustering pattern is robust to contributions of unknown or
unobserved niche axes. However, it may not be visible unless measured traits
are close proxies for niche strategies. We conclude that patterns along single
niche axes may reveal properties of interspecific competition in nature, but
detecting these patterns requires natural history expertise firmly tying traits
to niches.


  Autoreactive B cells have a central role in the pathogenesis of rheumatoid
arthritis (RA), and recent findings have proposed that anti-citrullinated
protein autoantibodies (ACPA) may be directly pathogenic. Herein, we
demonstrate the frequency of variable-region glycosylation in single-cell
cloned mAbs. A total of 14 ACPA mAbs were evaluated for predicted N-linked
glycosylation motifs in silico and compared to 452 highly-mutated mAbs from RA
patients and controls. Variable region N-linked motifs (N-X-S/T) were
strikingly prevalent within ACPA (100%) compared to somatically hypermutated
(SHM) RA bone marrow plasma cells (21%), and synovial plasma cells from
seropositive (39%) and seronegative RA (7%). When normalized for SHM, ACPA
still had significantly higher frequency of N-linked motifs compared to all
studied mAbs including highly-mutated HIV broadly-neutralizing and
malaria-associated mAbs. The Fab glycans of ACPA-mAbs were highly sialylated,
contributed to altered charge, but did not influence antigen binding. The
analysis revealed evidence of unusual B-cell selection pressure and
SHM-mediated decreased in surface charge and isoelectric point in ACPA. It is
still unknown how these distinct features of anti-citrulline immunity may have
an impact on pathogenesis. However, it is evident that they offer selective
advantages for ACPA+ B cells, possibly also through non-antigen driven
mechanisms.


  Selection of appropriate collective variables for enhancing sampling of
molecular simulations remains an unsolved problem in computational biophysics.
In particular, picking initial collective variables (CVs) is particularly
challenging in higher dimensions. Which atomic coordinates or transforms there
of from a list of thousands should one pick for enhanced sampling runs? How
does a modeler even begin to pick starting coordinates for investigation? This
remains true even in the case of simple two state systems and only increases in
difficulty for multi-state systems. In this work, we solve the initial CV
problem using a data-driven approach inspired by the filed of supervised
machine learning. In particular, we show how the decision functions in
supervised machine learning (SML) algorithms can be used as initial CVs
(SML_cv) for accelerated sampling. Using solvated alanine dipeptide and
Chignolin mini-protein as our test cases, we illustrate how the distance to the
Support Vector Machines' decision hyperplane, the output probability estimates
from Logistic Regression, the outputs from deep neural network classifiers, and
other classifiers may be used to reversibly sample slow structural transitions.
We discuss the utility of other SML algorithms that might be useful for
identifying CVs for accelerating molecular simulations.


  A common data mining task on networks is community detection, which seeks an
unsupervised decomposition of a network into structural groups based on
statistical regularities in the network's connectivity. Although many methods
exist, the No Free Lunch theorem for community detection implies that each
makes some kind of tradeoff, and no algorithm can be optimal on all inputs.
Thus, different algorithms will over or underfit on different inputs, finding
more, fewer, or just different communities than is optimal, and evaluation
methods that use a metadata partition as a ground truth will produce misleading
conclusions about general accuracy. Here, we present a broad evaluation of over
and underfitting in community detection, comparing the behavior of 16
state-of-the-art community detection algorithms on a novel and structurally
diverse corpus of 406 real-world networks. We find that (i) algorithms vary
widely both in the number of communities they find and in their corresponding
composition, given the same input, (ii) algorithms can be clustered into
distinct high-level groups based on similarities of their outputs on real-world
networks, and (iii) these differences induce wide variation in accuracy on link
prediction and link description tasks. We introduce a new diagnostic for
evaluating overfitting and underfitting in practice, and use it to roughly
divide community detection methods into general and specialized learning
algorithms. Across methods and inputs, Bayesian techniques based on the
stochastic block model and a minimum description length approach to
regularization represent the best general learning approach, but can be
outperformed under specific circumstances. These results introduce both a
theoretically principled approach to evaluate over and underfitting in models
of network community structure and a realistic benchmark by which new methods
may be evaluated and compared.


  Integration of transcriptomic and metabolomic data improves functional
interpretation of disease-related metabolomic phenotypes, and facilitates
discovery of putative metabolite biomarkers and gene targets. For this reason,
these data are increasingly collected in large cohorts, driving a need for the
development of novel methods for their integration. Of note,
clinical/translational studies typically provide snapshot gene and metabolite
profiles and, oftentimes, most metabolites are not identified. Thus, in these
types of studies, pathway/network approaches that take into account the
complexity of gene-metabolite relationships may neither be applicable nor
readily uncover novel relationships. With this in mind, we propose a simple
linear modeling approach to capture phenotype-specific gene-metabolite
associations, with the assumption that co-regulation patterns reflect
functionally related genes and metabolites. The proposed linear model,
metabolite ~ gene + phenotype + gene:phenotype, specifically evaluates whether
gene-metabolite relationships differ by phenotype, by testing whether the
relationship in one phenotype is significantly different from the relationship
in another phenotype (via an interaction gene:phenotype p-value). Interaction
p-values for all possible gene-metabolite pairs are computed and significant
pairs are clustered by the directionality of associations. We implemented our
approach as an R package, IntLIM, which includes a user-friendly Shiny app. We
applied IntLIM to two published datasets, collected in NCI-60 cell lines and in
human breast tumor and non-tumor tissue. We demonstrate that IntLIM captures
relevant tumor-specific gene-metabolite associations involved in cancer-related
pathways. and also uncover novel relationships that could be tested
experimentally. The IntLIM R package is publicly available in GitHub
(https://github.com/mathelab/IntLIM).


  We introduce a framework for studying social dilemmas in networked societies
where individuals follow a simple state-based behavioral mechanism based on
generalized reciprocity, which is rooted in the principle "help anyone if
helped by someone". Within this general framework, which applies to a wide
range of social dilemmas including, among others, public goods, donation and
snowdrift games, we study the cooperation dynamics on a variety of complex
network examples. By interpreting the studied model through the lenses of
nonlinear dynamical systems, we show that cooperation through generalized
reciprocity always emerges as the unique attractor in which the overall level
of cooperation is maximized, while simultaneously exploitation of the
participating individuals is prevented. The analysis elucidates the role of the
network structure, here captured by a local centrality measure which uniquely
quantifies the propensity of the network structure to cooperation, by dictating
the degree of cooperation displayed both at microscopic and macroscopic level.
We demonstrate the applicability of the analysis on a practical example by
considering an interaction structure that couples a donation process with a
public goods game.


  Transcription factors regulate gene expression, but how these proteins
recognize and specifically bind to their DNA targets is still debated. Machine
learning models are effective means to reveal interaction mechanisms. Here we
studied the ability of a quantum machine learning approach to predict binding
specificity. Using simplified datasets of a small number of DNA sequences
derived from actual binding affinity experiments, we trained a commercially
available quantum annealer to classify and rank transcription factor binding.
The results were compared to state-of-the-art classical approaches for the same
simplified datasets, including simulated annealing, simulated quantum
annealing, multiple linear regression, LASSO, and extreme gradient boosting.
Despite technological limitations, we find a slight advantage in classification
performance and nearly equal ranking performance using the quantum annealer for
these fairly small training data sets. Thus, we propose that quantum annealing
might be an effective method to implement machine learning for certain
computational biology problems.


  Ecosystems dynamics is often considered as driven by a coupling of species'
resource consumption and its population size dynamics. Such resource-population
dynamics is captured by MacArthur-type models. One biologically relevant
feature that would also need to be captured by such models is the introduction
of new and different species. Speciation introduces a stochastic component in
the otherwise deterministic MacArthur theory. We describe here how speciation
can be implemented to yield a model that is consistent with current theory on
equilibrium resource-consumer models, but also displays readily observable rank
diversity metric changes. The model also reproduces a priority effect. Adding
speciation to a MacArthur-style model so provides an attractively simple
extension to explore the rich dynamics in evolving ecosystems.


  Evolutionarily stable strategy (ESS) is an important solution concept in game
theory which has been applied frequently to biological models. Informally an
ESS is a strategy that if followed by the population cannot be taken over by a
mutation strategy that is initially rare. Finding such a strategy has been
shown to be difficult from a theoretical complexity perspective. We present an
algorithm for the case where mutations are restricted to pure strategies, and
present experiments on several game classes including random and a
recently-proposed cancer model. Our algorithm is based on a mixed-integer
non-convex feasibility program formulation, which constitutes the first general
optimization formulation for this problem. It turns out that the vast majority
of the games included in the experiments contain ESS with small support, and
our algorithm is outperformed by a support-enumeration based approach. However
we suspect our algorithm may be useful in the future as games are studied that
have ESS with potentially larger and unknown support size.


  Despite the extreme diversity of T cell repertoires, many identical T-cell
receptor (TCR) sequences are found in a large number of individual mice and
humans. These widely-shared sequences, often referred to as `public', have been
suggested to be over-represented due to their potential immune functionality or
their ease of generation by V(D)J recombination. Here we show that even for
large cohorts the observed degree of sharing of TCR sequences between
individuals is well predicted by a model accounting for by the known
quantitative statistical biases in the generation process, together with a
simple model of thymic selection. Whether a sequence is shared by many
individuals is predicted to depend on the number of queried individuals and the
sampling depth, as well as on the sequence itself, in agreement with the data.
We introduce the degree of publicness conditional on the queried cohort size
and the size of the sampled repertoires. Based on these observations we propose
a public/private sequence classifier, `PUBLIC' (Public Universal Binary
Likelihood Inference Classifier), based on the generation probability, which
performs very well even for small cohort sizes.


  In view of recent intense experimental and theoretical interests in the
biophysics of liquid-liquid phase separation (LLPS) of intrinsically disordered
proteins (IDPs), heteropolymer models with chain molecules configured as
self-avoiding walks on the simple cubic lattice are constructed to study how
phase behaviors depend on the sequence of monomers along the chains. To address
pertinent general principles, we focus primarily on two fully charged
50-monomer sequences with significantly different charge patterns. Each monomer
in our models occupies a single lattice site and all monomers interact via a
screened pairwise Coulomb potential. Phase diagrams are obtained by extensive
Monte Carlo sampling performed at multiple temperatures on ensembles of 300
chains in boxes of sizes ranging from $52\times 52\times 52$ to $246\times
246\times 246$ to simulate a large number of different systems with the overall
polymer volume fraction $\phi$ in each system varying from $0.001$ to $0.1$.
Phase separation in the model systems is characterized by the emergence of a
large cluster connected by inter-monomer nearest-neighbor lattice contacts and
by large fluctuations in local polymer density. The simulated critical
temperatures, $T_{\rm cr}$, of phase separation for the two sequences differ
significantly, whereby the sequence with a more "blocky" charge pattern
exhibits a substantially higher propensity to phase separate. The trend is
consistent with our sequence-specific random-phase-approximation (RPA) polymer
theory, but the variation of the simulated $T_{\rm cr}$ with a previously
proposed "sequence charge decoration" pattern parameter is milder than that
predicted by RPA. Ramifications of our findings for the development of
analytical theory and simulation protocols of IDP LLPS are discussed.


  In recent years non-demographic variability has been shown to greatly affect
dynamics of stochastic populations. For example, non-demographic noise in the
form of a bursty reproduction process with an a-priori unknown burst size, or
environmental variability in the form of time-varying reaction rates, have been
separately found to dramatically impact the extinction risk of isolated
populations. In this work we investigate the extinction risk of an isolated
population under the combined influence of these two types of non-demographic
variation. Using the so-called momentum-space WKB approach we arrive at a set
of time-dependent Hamilton equations. In order to account for the explicit time
dependence, we find the instanton of the time-perturbed Hamiltonian
numerically, where analytical expressions are presented in particular limits
using various perturbation techniques. We focus on two classes of time-varying
environments: periodically-varying rates corresponding to seasonal effects, and
a sudden decrease in the birth rate corresponding to a catastrophe. All our
theoretical results are tested against numerical Monte Carlo simulations with
time-dependent rates and also against a numerical solution of the corresponding
time-dependent Hamilton equations.


  The transient response to a stimulus and subsequent recovery to a steady
state are the fundamental characteristics of a living organism. Here we study
the relaxation kinetics of autoregulatory gene networks based on the chemical
master equation model of single-cell stochastic gene expression with nonlinear
feedback regulation. We report a novel relation between the rate of relaxation,
characterized by the spectral gap of the Markov model, and the feedback sign of
the underlying gene circuit. When a network has no feedback, the relaxation
rate is exactly the decaying rate of the protein. We further show that positive
feedback always slows down the relaxation kinetics while negative feedback
always speeds it up. Numerical simulations demonstrate that this relation
provides a possible method to infer the feedback topology of autoregulatory
gene networks by using time-series data of gene expression.


  Solanum melongena (eggplant) is one of the diversity of the Solanum family
which is grown and widely spread in Indonesia and widely used by the community.
This research explored the genetic diversity of four local Indonesian eggplant
species namely leuca, tekokak, gelatik and kopek by using RAPD (Random
Amplified Polymorphic DNA). The samples were obtained from Agricultural
Technology Assessment Institute (BPTP) Bogor, Indonesia. The result of data
observation was in the form of Solanum melongena plants DNA profile analyzed
descriptively and quantitatively. 30 DNA bands (28 polymorphic and 2
monomorphic) were successfully scored by using four primers (OPF01, OPF02,
OPF03 and OPF04). The Primers were used able to amplify all of the four
eggplant samples. The result of PCR-RAPD visualization produces bands of
300-1500 bp. The result of cluster analysis showed the existence of three
clusters (A, B, and C). Cluster A (coefficient of equal to 49%) consisted of a
gelatik, cluster B (coefficient of 65% equilibrium) consisted of TPU (Kopek)
and TK (Tekokak), and cluster C (55% equilibrium coefficient) consisted of LC
(Leunca). These results indicated that the closest proximity is found in
samples of TK (Tekokak) and TPU (Kopek)


  Threshold theorem is probably the most important development of mathematical
epidemic modelling. Unfortunately, some models may not behave according to the
threshold. In this paper, we will focus on the final outcome of SIR model with
demography. The behaviour of the model approached by deteministic and
stochastic models will be introduced, mainly using simulations. Furthermore, we
will also investigate the dynamic of susceptibles in population in absence of
infective. We have successfully showed that both deterministic and stochastic
models performed similar results when $R_0 \leq 1$. That is, the disease-free
stage in the epidemic. But when $R_0 > 1$, the deterministic and stochastic
approaches had different interpretations.


  Studies about epidemic modelling have been conducted since before 19th
century. Both deterministic and stochastiic model were used to capture the
dynamic of infection in the population. The purpose of this project is to
investigate the behaviour of the models when we set the basic reproduction
number, $R_0$. This quantity is defined as the expected number of contacts made
by a typical infective to susceptibles in the population. According to the
epidemic threshold theory, when $R_0 \leq 1$, minor epidemic occurs with
probability one in both approaches, but when $R_0 > 1$, the deterministic and
stochastic models have different interpretation. In the deterministic approach,
major epidemic occurs with probability one when $R_0 > 1$ and predicts that the
disease will settle down to an endemic equilibrium. Stochastic models, on the
other hand, identify that the minor epidemic can possibly occur. If it does,
then the epidemic will die out quickly. Moreover, if we let the population size
be large and the major epidemic occurs, then it will take off and then reach
the endemic level and move randomly around the deterministic's equilibrium.


  Biosystems contain an almost infinite amount of vital important details,
which together ensure their life. There are, however, some common structures
and reactions in the systems: the homochirality of carbohydrates and proteins,
the metabolism and the genetics. The Abiogenesis, or the origin of life, is
probably not a result of a series of single events, but rather the result of a
gradual process with increasing complexity of molecules and chemical reactions,
and the prebiotic synthesis of molecules might not have left a trace of the
establishment of structures and reactions at the beginning of the evolution.
But alternatively, one might be able to determine some order in the formation
of the chemical denominators in the Abiogenesis. Here we review experimental
results and present a model of the start of the Abionenesis, where the
spontaneous formation of homochirality in proteins is the precondition for the
establishment of homochirality of carbohydrates and for the metabolism at the
start of the Abiogenesis.