  Icosahedral viruses have their infectious genome encapsulated by a shell
assembled by a multiscale process, starting from an integer multiple of 60
viral capsid or coat protein (VP) monomers.
  We predict and validate inter-atomic hotspot interactions between VP monomers
that are important for the assembly of 3 icosahedral viral capsids: Adeno
Associated Virus serotype 2 (AAV2) and Minute Virus of Mice (MVM), both T=1
single stranded DNA viruses, and Bromo Mosaic Virus (BMV), a T=3 single
stranded RNA virus. Experimental validation is by in-vitro, site-directed
mutagenesis data found in literature.
  We combine ab-initio predictions at two scales: at the interface-scale, we
predict the importance (cruciality) of an interaction for successful
subassembly across each interface between VP monomers; and at the capsid-scale,
we predict the cruciality of an interface for successful capsid assembly.
  At the interface-scale, we measure cruciality by changes in the capsid
free-energy landscape partition function when an interaction is removed. The
partition function computation uses atlases of interface subassembly
landscapes, rapidly generated by a novel geometric method and curated
opensource software EASAL (efficient atlasing and search of assembly
landscapes). At the capsid-scale, cruciality of an interface for successful
assembly of the capsid is based on combinatorial entropy.
  Our study goes from resource-light, multiscale computational predictions of
crucial hotspot inter-atomic interactions to validation using data on
site-directed mutagenesis' effect on capsid assembly. By reliably and rapidly
narrowing down target interactions, (no more than 1.5 hours per interface on a
laptop with Intel Core i5-2500K 3.2Ghz CPU and 8GB of RAM) our predictions can
inform and reduce time-consuming in-vitro and in-vivo experiments, or more
computationally intensive in-silico analyses.


  Background: Although recorded cancer mortality should include both deaths
from cancer and deaths from cancer treatment, there is evidence suggesting that
the measure may be incomplete. To investigate the completeness of recorded
prostate cancer mortality, we compared other-cause (non-prostate cancer)
mortality in men found and not found to have prostate cancer following a needle
biopsy.
  Methods: We linked Medicare claims data to SEER data to analyze survival in
the population of men aged 65+ enrolled in Medicare who resided in a SEER area
and received a needle biopsy in 1993-2001. We compared other-cause mortality in
men found to have prostate cancer (n=53,462) to that in men not found to have
prostate cancer (n=103,659).
  Results: The age-race adjusted other-cause mortality rate was 471 per 10,000
person-years in men found to have prostate cancer vs. 468 per 10,000 in men not
found to have prostate cancer (RR = 1.01;95% CI:0.98-1.03). The effect was
modified, however, by age. The RR declined in a stepwise fashion from 1.08 (95%
CI:1.03-1.14) in men age 65-69 to 0.89 (95% CI:0.83-0.95) in men age 85 and
older. If the excess (or deficit) in other-cause mortality were added to the
recorded prostate cancer mortality, prostate cancer mortality would rise 23% in
the youngest age group (from 90 to 111 per 10,000) and would fall 30% in the
oldest age group (from 551 to 388 per 10,000).
  Conclusion: Although recorded prostate cancer mortality appears to be an
accurate measure overall, it systematically underestimates the mortality
associated with prostate cancer diagnosis and treatment in younger men and
overestimates it in the very old. We surmise that in younger men
treatment-related deaths are incompletely captured in recorded prostate cancer
mortality, while in older men the diagnosis "sticks"-- once diagnosed, they are
more likely to be said to have died from the disease.


  A detailed examination of the effect of harvesting on a population has been
carried out by extending the standard cubic deterministic model by considering
a population under Allee effect with a quadratic function representing
harvesting. Weak and strong Allee effect transitions, carrying capacity, and
Allee threshold change according to harvesting is first discussed in the
deterministic model. A Fokker Planck equation has been obtained starting from a
Langevin equation subject to correlated Gaussian white noise with zero mean,
and an Approximate Fokker Planck Equation has been obtained from a Langevin
equation subject to correlated Gaussian colored noise with zero mean. This
allowed to calculate the stationary probability distributions of populations,
and thus to discuss the effects of linear and nonlinear (Holling type-II)
harvesting for populations under Allee effect and subject to white and colored
noises respectively.


  The abstract mathematical rules of artificial neural network (ANN) are
implemented through computation using electronic computers, photonics and
in-vitro DNA computation. Here we demonstrate the physical realization of ANN
in living bacterial cells. We created a single layer ANN using engineered
bacteria, where a single bacterium works as an artificial neuron and
demonstrated a 2-to-4 decoder and a 1-to-2 de-multiplexer for processing
chemical signals. The inputs were extracellular chemical signals, which
linearly combined and got processed through a non-linear log-sigmoid activation
function to produce fluorescent protein outputs. The activation function was
generated by synthetic genetic circuits, and for each artificial neuron, the
weight and bias values were adjusted manually by engineering the molecular
interactions within the bacterial neuron to represent a specific logical
function. The artificial bacterial neurons were connected as ANN architectures
to implement a 2-to-4 chemical decoder and a 1-to-2 chemical de-multiplexer. To
our knowledge, this is the first ANN created by artificial bacterial neurons.
Thus, it may open up a new direction in ANN research, where engineered
biological cells can be used as ANN enabled hardware.


  The cell cycle is tightly regulated by cyclins and their catalytic moieties,
the cyclin-dependent kinases (CDKs). Cyclin D1, in association with CDK4/6,
acts as a mitogenic sensor and integrates extracellular mitogenic signals and
cell cycle progression. When deregulated (overexpressed, accumulated,
inappropriately located), cyclin D1 becomes an oncogene and is recognized as a
driver of solid tumors and hemopathies. Recent studies on the oncogenic roles
of cyclin D1 reported non-canonical functions dependent on the partners of
cyclin D1 and its location within tumor cells or tissues. Support for these new
functions was provided by various mouse models of oncogenesis. Finally,
proteomic and transcriptomic data identified complex cyclin D1 networks. This
review focuses on these aspects of cyclin D1 pathophysiology, which may be
crucial for targeted therapy.


  We consider a nonautonomous eco-epidemiological model with general functions
for predation on infected and uninfected preys as well as general functions
associated to the vital dynamics of the susceptible prey and predator
populations. We obtain persistence and extinction results for the infected prey
based on assumptions on systems related to the dynamics in the absence of
infected preys. We apply our results to eco-epidemiological models constructed
from predator-prey models existent in the literature. Some illustrative
simulation is undertaken.


  Many of the commonly used methods for orthology detection start from mutually
most similar pairs of genes (reciprocal best hits) as an approximation for
evolutionary most closely related pairs of genes (reciprocal best matches).
This approximation of best matches by best hits becomes exact for ultrametric
dissimilarities, i.e., under the Molecular Clock Hypothesis. It fails, however,
whenever there are large lineage specific rate variations among paralogous
genes. In practice, this introduces a high level of noise into the input data
for best-hit-based orthology detection methods. If additive distances between
genes are known, then evolutionary most closely related pairs can be identified
by considering certain quartets of genes provided that in each quartet the
outgroup relative to the remaining three genes is known. \emph{A priori}
knowledge of underlying species phylogeny greatly facilitates the
identification of the required outgroup. Although the workflow remains a
heuristic since the correct outgroup cannot be determined reliably in all
cases, simulations with lineage specific biases and rate asymmetries show that
nearly perfect results can be achieved. In a realistic setting, where distances
data have to be estimated from sequence data and hence are noisy, it is still
possible to obtain highly accurate sets of best matches. Improvements of
tree-free orthology assessment methods can be expected from a combination of
the accurate inference of best matches reported here and recent mathematical
advances in the understanding of (reciprocal) best match graphs and orthology
relations.


  In recent years, the advances in single-cell RNA-seq techniques have enabled
us to perform large-scale transcriptomic profiling at single-cell resolution in
a high-throughput manner. Unsupervised learning such as data clustering has
become the central component to identify and characterize novel cell types and
gene expression patterns. In this study, we review the existing single-cell
RNA-seq data clustering methods with critical insights into the related
advantages and limitations. In addition, we also review the upstream
single-cell RNA-seq data processing techniques such as quality control,
normalization, and dimension reduction. We conduct performance comparison
experiments to evaluate several popular single-cell RNA-seq clustering
approaches on two single-cell transcriptomic datasets.


  The ability to consistently distinguish real protein structures from
computationally generated model decoys is not yet a solved problem. One route
to distinguish real protein structures from decoys is to delineate the
important physical features that specify a real protein. For example, it has
long been appreciated that the hydrophobic cores of proteins contribute
significantly to their stability. As a dataset of decoys to compare with real
protein structures, we studied submissions to the bi-annual CASP competition
(specifically CASP11, 12, and 13), in which researchers attempt to predict the
structure of a protein only knowing its amino acid sequence. Our analysis
reveals that many of the submissions possess cores that do not recapitulate the
features that define real proteins. In particular, the model structures appear
more densely packed (because of energetically unfavorable atomic overlaps),
contain too few residues in the core, and have improper distributions of
hydrophobic residues throughout the structure. Based on these observations, we
developed a deep learning method, which incorporates key physical features of
protein cores, to predict how well a computational model recapitulates the real
protein structure without knowledge of the structure of the target sequence. By
identifying the important features of protein structure, our method is able to
rank decoys from the CASP competitions equally well, if not better than,
state-of-the-art methods that incorporate many additional features.


  Motivated by models of cancer formation in which cells need to acquire $k$
mutations to become cancerous, we consider a spatial population model in which
the population is represented by the $d$-dimensional torus of side length $L$.
Initially, no sites have mutations, but sites with $i-1$ mutations acquire an
$i$th mutation at rate $\mu_i$ per unit area. Mutations spread to neighboring
sites at rate $\alpha$, so that $t$ time units after a mutation, the region of
individuals that have acquired the mutation will be a ball of radius $\alpha
t$. We calculate, for some ranges of the parameter values, the asymptotic
distribution of the time required for some individual to acquire $k$ mutations.
Our results, which build on previous work of Durrett, Foo, and Leder, are
essentially complete when $k = 2$ and when $\mu_i = \mu$ for all $i$.


  Comparing and aligning large datasets is a pervasive problem occurring across
many different knowledge domains. We introduce and study MREC, a recursive
decomposition algorithm for computing matchings between data sets. The basic
idea is to partition the data, match the partitions, and then recursively match
the points within each pair of identified partitions. The matching itself is
done using black box matching procedures that are too expensive to run on the
entire data set. Using an absolute measure of the quality of a matching, the
framework supports optimization over parameters including partitioning
procedures and matching algorithms. By design, MREC can be applied to extremely
large data sets. We analyze the procedure to describe when we can expect it to
work well and demonstrate its flexibility and power by applying it to a number
of alignment problems arising in the analysis of single cell molecular data.


  Deep learning is playing a vital role in every field which involves data. It
has emerged as a strong and efficient framework that can be applied to a broad
spectrum of complex learning problems which were difficult to solve using
traditional machine learning techniques in the past. In this study we focused
on classification of protein sequences with deep learning techniques. The study
of amino acid sequence is vital in life sciences. We used different word
embedding techniques from Natural Language processing to represent the amino
acid sequence as vectors. Our main goal was to classify sequences to four group
of classes, that are DNA, RNA, Protein and hybrid. After several tests we have
achieved almost 99% of train and test accuracy. We have experimented on CNN,
LSTM, Bidirectional LSTM, and GRU.


  The resurgence of measles is largely attributed to the decline in vaccine
adoption and the increase in mobility. Although the vaccine for measles is
readily available and highly successful, its current adoption is not adequate
to prevent epidemics. Vaccine adoption is directly affected by individual
vaccination decisions, and has a complex interplay with the spatial spread of
disease shaped by an underlying mobility (travelling) network. In this paper,
we model the travelling connectivity as a scale-free network, and investigate
dependencies between the network's assortativity and the resultant epidemic and
vaccination dynamics. In doing so we extend an SIR-network model with
game-theoretic components, capturing the imitation dynamics under a voluntary
vaccination scheme. Our results show a correlation between the epidemic
dynamics and the network's assortativity, highlighting that networks with high
assortativity tend to suppress epidemics under certain conditions. In highly
assortative networks, the suppression is sustained producing an early
convergence to equilibrium. In highly disassortative networks, however, the
suppression effect diminishes over time due to scattering of non-vaccinating
nodes, and frequent switching between the predominantly vaccinating and
non-vaccinating phases of the dynamics.


  The problem of matching a query string to a directed graph, whose vertices
are labeled by strings, has application in different fields, from data mining
to computational biology. Several variants of the problem have been considered,
depending on the fact that the match is exact or approximate and, in this
latter case, which edit operations are considered and where are allowed. In
this paper we present results on the complexity of the approximate matching
problem, where edit operations are symbol substitutions and are allowed only on
the graph labels or both on the graph labels and the query string. We introduce
a variant of the problem that asks whether there exists a path in a graph that
represents a query string with any number of edit operations and we show that
is is NP-complete, even when labels have length one and in the case the
alphabet is binary. Moreover, when it is parameterized by the length of the
input string and graph labels have length one, we show that the problem is
fixed-parameter tractable and it is unlikely to admit a polynomial kernel. The
NP-completeness of this problem leads to the inapproximability (within any
factor) of the approximate matching when edit operations are allowed only on
the graph labels. Moreover, we show that the variants of approximate string
matching to graph we consider are not fixed-parameter tractable, when the
parameter is the number of edit operations, even for graphs that have distance
one from a DAG. The reduction for this latter result allows us to prove the
inapproximability of the variant where edit operations can be applied both on
the query string and on graph labels.


  The ergodicity and the output-controllability of stochastic reaction networks
have been shown to be essential properties to fulfill to enable their control
using, for instance, antithetic integral control. We propose here to extend
those properties to the case of uncertain networks. To this aim, the notions of
interval, robust, sign, and structural ergodicity/output-controllability are
introduced. The obtained results lie in the same spirit as those obtained in
[Briat, Gupta & Khammash, Cell Systems, 2016] where those properties are
characterized in terms of control theoretic concepts, linear algebraic
conditions, linear programs, and graph-theoretic/algebraic conditions. An
important conclusion is that all those properties can be characterized by
linear programs. Two examples are given for illustration.


  Hydrogenase occupy a central place in the energy metabolism of anaerobic
bacteria. Although the structure of mitochondrial complex I is similar to that
of hydrogenase, whether it has hydrogen metabolic activity remain unclear.
Here, we show that a H2 evolving activity exists in higher plants mitochondria
and is closely related to complex I, especially around ubiquinone binding site.
The H2 production could be inhibited by rotenone and ubiquinone. Hypoxia could
simultaneously promote H2 evolution and succinate accumulation. Redox
properties of quinone pool, adjusted by NADH or succinate according to oxygen
concentration, acts as a valve to control the flow of protons and electrons and
the production of H2. The coupling of H2 evolving activity of mitochondrial
complex I with metabolic regulation reveals a more effective redox homeostasis
regulation mechanism. Considering the ubiquity of mitochondria in eukaryotes,
H2 metabolism might be the innate function of higher organisms. This may serve
to explain, at least in part, the broad physiological effects of H2.


  The genetic composition of a naturally developing population is considered as
due to mutation, selection, genetic drift and recombination. Selection is
modeled as single-locus terms (additive fitness) and two-loci terms (pairwise
epistatic fitness). The problem is posed to infer epistatic fitness from
population-wide whole-genome data from a time series of a developing
population. We generate such data in silico, and show that in the Quasi-Linkage
Equilibrium (QLE) phase of Kimura, Neher and Shraiman, that pertains at high
enough recombination rates and low enough mutation rates, epistatic fitness can
be quantitatively correctly inferred using inverse Ising/Potts methods.


  The calculation of thermodynamic properties of biochemical systems typically
requires the use of resource-intensive molecular simulation methods. One
example thereof is the thermodynamic profiling of hydration sites, i.e.
high-probability locations for water molecules on the protein surface, which
play an essential role in protein-ligand associations and must therefore be
incorporated in the prediction of binding poses and affinities. To replace
time-consuming simulations in hydration site predictions, we developed two
different types of deep neural-network models aiming to predict hydration site
data. In the first approach, meshed 3D images are generated representing the
interactions between certain molecular probes placed on regular 3D grids,
encompassing the binding pocket, with the static protein. These molecular
interaction fields are mapped to the corresponding 3D image of hydration
occupancy using a neural network based on an U-Net architecture. In a second
approach, hydration occupancy and thermodynamics were predicted point-wise
using a neural network based on fully-connected layers. In addition to direct
protein interaction fields, the environment of each grid point was represented
using moments of a spherical harmonics expansion of the interaction properties
of nearby grid points. Application to structure-activity relationship analysis
and protein-ligand pose scoring demonstrates the utility of the predicted
hydration information.


  We perform a detailed analysis of the behaviour of a non-autonomous
prey-predator model where age based growth with age discriminatory harvesting
in prey and predator's reliance upon alternative food in the absence of that
particular prey are considered. We begin by deriving certain sufficient
conditions for permanence and positive invariance and then proceed to construct
a Lyapunov function to derive some constraints for global attractivity. With
the help of continuation theorem we arrive at the best fit criterion to prove
the occurrence of a positive periodic solution. Moreover, using Arzela-Ascoli
theorem we formulate a proof for a unique positive solution to be almost
periodic and we carry out numerical simulation to verify the analytical
findings. With the aid of graphs and tables we show the nature of the
prey-predator system in response to the alternative food and delays.


  A stochastic model of autoregulated bursty gene expression by Kumar et al.
[Phys. Rev. Lett. 113, 268105 (2014)] has been exactly solved in steady-state
conditions under the implicit assumption that protein numbers are sufficiently
large such that fluctuations in protein numbers due to reversible
protein-promoter binding can be ignored. Here we derive an alternative model
that takes into account these fluctuations and hence can be used to study low
protein number effects. The exact steady-state protein number distributions is
derived as a sum of Gaussian hypergeometric functions. We use the theory to
study how promoter switching rates and the type of feedback influence the size
of protein noise and noise-induced bistability. Furthermore we show that our
model predictions for the protein number distribution are significantly
different from those of Kumar et al. when the protein mean is small, gene
switching is fast, and protein binding is faster than unbinding.


  It is becoming widely accepted that very early in the origin of life, even
before the emergence of genetic encoding, reaction networks of diverse small
chemicals might have manifested key properties of life, namely self-propagation
and adaptive evolution. To explore this possibility, we formalize the dynamics
of chemical reaction networks within the framework of chemical ecosystem
ecology. To capture the idea that life-like chemical systems are maintained out
of equilibrium by fluxes of energy-rich food chemicals, we model chemical
ecosystems in well-mixed containers that are subject to constant dilution by a
solution with a fixed concentration of food chemicals. Modelling all chemical
reactions as fully reversible, we show that seeding an autocatalytic cycle (AC)
with tiny amounts of one or more of its member chemicals results in logistic
growth of all member chemicals in the cycle. This finding justifies drawing an
instructive analogy between an AC and the population of a biological species.
We extend this finding to show that pairs of ACs can show competitive,
predator-prey, or mutualistic associations just like biological species.
Furthermore, when there is stochasticity in the environment, particularly in
the seeding of ACs, chemical ecosystems can show complex dynamics that can
resemble evolution. The evolutionary character is especially clear when the
network architecture results in ecological precedence (survival of the first),
which makes the path of succession historically contingent on the order in
which cycles are seeded. For all its simplicity, the framework developed here
is helpful for visualizing how autocatalysis in prebiotic chemical reaction
networks can yield life-like properties. Furthermore, chemical ecosystem
ecology could provide a useful foundation for exploring the emergence of
adaptive dynamics and the origins of polymer-based genetic systems.


  Microscopies have become pillars of our characterization tools to observe
biological systems and assemblies. Correlative and synchronous use of different
microscopies relies on the fundamental assumption of non-interference during
images acquisitions. In this work, by exploring the correlative use of Atomic
Force Microscopy and confocal-Fluorescence-Lifetime Imaging Microscopy
(AFM-FLIM), we quantify cross-talk effects occurring during synchronous
acquisition. We characterize and minimize optomechanical forces on different
AFM cantilevers interfering with normal AFM operation as well as spurious
luminescence from the tip and cantilever affecting time-resolved fluorescence
detection. By defining non-interfering experimental imaging parameters, we show
accurate real-time acquisition and two-dimensional mapping of interaction
force, fluorescence lifetime and intensity characterizing morphology (AFM) and
local viscosity (FLIM) of gel and fluid phases separation of supported lipid
model membranes. Finally, as proof of principle by means of synchronous force
and fluorescence spectroscopies, we precisely tune the lifetime of a
fluorescent nanodiamond positioned on the AFM tip by controlling its distance
from a metallic surface. This opens up a novel pathway of quench sensing to
image soft biological samples such as membranes since it does not require
tip-sample mechanical contact in contrast with conventional AFM in liquid.


  Proteins and nucleic acids can spontaneously self-assemble into membraneless
droplet-like compartments, both in vitro and in vivo. A key component of these
droplets are multi-valent proteins that possess several adhesive domains with
specific interaction partners (whose number determines total valency of the
protein) separated by disordered regions. Here, using multi-scale simulations
we show that such proteins self-organize into micro-phase separated droplets of
various sizes as opposed to the Flory-like macro-phase separated equilibrium
state of homopolymers or equilibrium physical gels. We show that the
micro-phase separated state is a dynamic outcome of the interplay between two
competing processes: a diffusion-limited encounter between proteins, and the
dynamics within small clusters that results in exhaustion of available
valencies whereby all specifically interacting domains find their interacting
partners within smaller clusters, leading to arrested phase separation. We
first model these multi-valent chains as bead-spring polymers with multiple
adhesive domains separated by semi-flexible linkers and use Langevin Dynamics
(LD) to assess how key timescales depend on the molecular properties of
associating polymers. Using the time-scales from LD simulations, we develop a
coarse-grained kinetic model to study this phenomenon at longer times.
Consistent with LD simulations, the macro-phase separated state was only
observed at high concentrations and large interaction valencies. Further, in
the regime where cluster sizes approach macro-phase separation, the condensed
phase becomes dynamically solid-like, suggesting that it might no longer be
biologically functional. Therefore, the micro-phase separated state could be a
hallmark of functional droplets formed by proteins with the sticker-spacer
architecture.


  BACKGROUND. Epidemic models of STIs are often used to characterize the
contribution of risk groups to overall transmission by projecting the
transmission population attributable fraction (tPAF) of unmet prevention and
treatment needs within risk groups. However, evidence suggests that STI risk is
dynamic over an individual's sexual life course, which manifests as turnover
between risk groups. We sought to examine the mechanisms by which turnover
influences modelled projections of the tPAF of high risk groups. METHODS. We
developed a unifying, data-guided framework to simulate risk group turnover in
deterministic, compartmental transmission models. We applied the framework to
an illustrative model of an STI and examined the mechanisms by which risk group
turnover influenced equilibrium prevalence across risk groups. We then fit a
model with and without turnover to the same risk-stratified STI prevalence
targets and compared the inferred level of risk heterogeneity and tPAF of the
highest risk group projected by the two models. RESULTS. The influence of
turnover on group-specific prevalence was mediated by three main phenomena:
movement of infectious individuals between risk groups; changes to herd
immunity; and changes in discordant partnerships. Faster turnover led to a
smaller ratio of STI prevalence between the highest and lowest risk groups.
Compared to the fitted model without turnover, the fitted model with turnover
inferred greater risk heterogeneity and consistently projected a larger tPAF of
the highest risk group over time. IMPLICATIONS. If turnover is not captured in
epidemic models, the projected contribution of high risk groups, and thus, the
potential impact of prioritizing interventions to address their needs, could be
underestimated. To aid the next generation of tPAF models, data collection
efforts to parameterize risk group turnover should be prioritized.


  The diversity of T-cell receptor (TCR) repertoires is achieved by a
combination of two intrinsically stochastic steps: random receptor generation
by VDJ recombination, and selection based on the recognition of random
self-peptides presented on the major histocompatibility complex. These
processes lead to a large receptor variability within and between individuals.
However, the characterization of the variability is hampered by the limited
size of the sampled repertoires. We introduce a new software tool SONIA to
facilitate inference of individual-specific computational models for the
generation and selection of the TCR beta chain (TRB) from sequenced repertoires
of 651 individuals, separating and quantifying the variability of the two
processes of generation and selection in the population. We find not only that
most of the variability is driven by the VDJ generation process, but there is a
large degree of consistency between individuals with the inter-individual
variance of repertoires being about 2% of the intra-individual variance. Known
viral-specific TCRs follow the same generation and selection statistics as all
TCRs.


  The tremendous boost in the next generation sequencing and in the omics
technologies makes it possible to characterize human gut microbiome (the
collective genomes of the microbial community that reside in our
gastrointestinal tract). While some of these microorganisms are considered as
essential regulators of our immune system, some others can cause several
diseases such as Inflammatory Bowel Diseases (IBD), diabetes, and cancer. IBD,
is a gut related disorder where the deviations from the healthy gut microbiome
are considered to be associated with IBD. Although existing studies attempt to
unveal the composition of the gut microbiome in relation to IBD diseases, a
comprehensive picture is far from being complete. Due to the complexity of
metagenomic studies, the applications of the state of the art machine learning
techniques became popular to address a wide range of questions in the field of
metagenomic data analysis. In this regard, using IBD associated metagenomics
dataset, this study utilizes both supervised and unsupervised machine learning
algorithms, i) to generate a classification model that aids IBD diagnosis, ii)
to discover IBD associated biomarkers, iii) to find subgroups of IBD patients
using k means and hierarchical clustering. To deal with the high dimensionality
of features, we applied robust feature selection algorithms such as Conditional
Mutual Information Maximization (CMIM), Fast Correlation Based Filter (FCBF),
min redundancy max relevance (mRMR) and Extreme Gradient Boosting (XGBoost). In
our experiments with 10 fold cross validation, XGBoost had a considerable
effect in terms of minimizing the microbiota used for the diagnosis of IBD and
thus reducing the cost and time. We observed that compared to the single
classifiers, ensemble methods such as kNN and logitboost resulted in better
performance measures for the classification of IBD.


  This paper is concerned with an analysis of the dynamics of a non-autonomous,
single population age based growth model with harvesting formulation. First, we
derive sufficient conditions for permanence and positive invariance. Then, by
constructing a scalar function, namely the Lyapunov function, we arrive at a
suitable criterion for global attractivity. With the help of Brouwer fixed
point and continuation theorems, we obtain constraints for the existence of a
positive periodic solution. Then we prove that there exists only one solution
which is almost periodic in nature that is distinct from every other solution.
Further, we carry out a numerical simulation to support the analytical
findings.


  Traditional approaches to elucidation of protein structures by NMR
spectroscopy rely on distance restraints also know as nuclear Overhauser
effects (NOEs). The use of NOEs as the primary source of structure
determination by NMR spectroscopy is time consuming and expensive. Residual
Dipolar Couplings (RDCs) have become an alternate approach for structure
calculation by NMR spectroscopy. In this work we report our results for
structure calculation of the novel protein PF2048.1 from RDC data and establish
the minimum data requirement for successful structure calculation using the
software package REDCRAFT. Our investigations start with utilizing four sets of
synthetic RDC data in two alignment media and proceed by reducing the RDC data
to the final limit of {CN, NH} and {NH} from two alignment media respectively.
Our results indicate that structure elucidation of this protein is possible
with as little as {CN, NH} and {NH} to within 0.533{\AA} of the target
structure.


  High-throughput cDNA sequencing (RNA-seq) is a very powerful technique to
quantify gene expression in an unbiased way. The Crustacean family is among the
groups of organisms sparsely represented in current genomic databases. Here we
present transcriptome data from Uca minax (red-jointed fiddler crab) as an
opportunity to extend our knowledge. Next generation sequencing was performed
on six tissue samples from Uca minax using the Illumina HiSeq system. Six
Transcriptome libraries were created using Trinity; a free, open-source
software tool for de novo transcriptome assembly of high-throughput mRNA
sequencing (RNA-seq) data with the absence of a reference genome. In addition,
several tools that aid in management of data were used, such as RSEM, Bowtie,
Blast, and IGV; a tool for visualizing RNA-seq analysis results. Fast quality
control (FastQC) analysis of the raw sequenced files revealed that both adapter
and PCR primer sequences were prevalently present, which may require a
preprocessing step.


  The arrival of Modern Humans (MHs) in Europe between 50 ka and 36 ka
coincides with significant changes in human behaviour, regarding the production
of tools, the exploitation of resources and the systematic use of ornaments and
colouring substances. The emergence of the so-called modern behaviours is
usually associated with MHs, although in these last decades findings relating
to symbolic thinking of pre-Sapiens groups have been claimed. In this paper we
present a synthesis of the Italian evidence concerning bone manufacturing and
the use of ornaments and pigments in the time span encompassing the demise of
Neandertals and their replacement by MHs. Current data show that Mousterian
bone tools are mostly obtained from bone fragments used as is. Conversely an
organized production of fine shaped bone tools is characteristic of the
Uluzzian and the Protoaurignacian, when the complexity inherent in the
manufacturing processes suggests that bone artefacts are not to be considered
as expedient resources. Some traces of symbolic activities are associated to
Neandertals in Northern Italy. Ornaments (mostly tusk shells) and pigments used
for decorative purposes are well recorded during the Uluzzian. Their features
and distribution witness to an intriguing cultural homogeneity within this
technocomplex. The Protoaurignacian is characterized by a wider archaeological
evidence, consisting of personal ornaments (mostly pierced gastropods),
pigments and artistic items.


  In this paper, we investigate potential biases in datasets used to make drug
binding predictions using machine learning. We investigate a recently published
metric called the Asymmetric Validation Embedding (AVE) bias which is used to
quantify this bias and detect overfitting. We compare it to a slightly revised
version and introduce a new weighted metric. We find that the new metrics allow
to quantify overfitting while not overly limiting training data and produce
models with greater predictive value.


  Fungal Biosynthetic Gene Clusters (BGCs) of secondary metabolites are
clusters of genes capable of producing natural products, compounds that play an
important role in the production of a wide variety of bioactive compounds,
including antibiotics and pharmaceuticals. Identifying BGCs can lead to the
discovery of novel natural products to benefit human health. Previous work has
been focused on developing automatic tools to support BGC discovery in plants,
fungi, and bacteria. Data-driven methods, as well as probabilistic and
supervised learning methods have been explored in identifying BGCs. Most
methods applied to identify fungal BGCs were data-driven and presented limited
scope. Supervised learning methods have been shown to perform well at
identifying BGCs in bacteria, and could be well suited to perform the same task
in fungi. But labeled data instances are needed to perform supervised learning.
Openly accessible BGC databases contain only a very small portion of previously
curated fungal BGCs. Making new fungal BGC datasets available could motivate
the development of supervised learning methods for fungal BGCs and potentially
improve prediction performance compared to data-driven methods. In this work we
propose new publicly available fungal BGC datasets to support the BGC discovery
task using supervised learning. These datasets are prepared to perform binary
classification and predict candidate BGC regions in fungal genomes. In addition
we analyse the performance of a well supported supervised learning tool
developed to predict BGCs.


  In modern age, the increasing complexity of computation and communication
technology is leading us towards the necessity of new paradigm. As a result,
unconventional approach like DNA coding theory is gaining considerable
attention. The storage capacity, information processing and transmission
properties of DNA molecules stimulate the notion of DNA coding theory as well
as DNA cryptography. In this paper we generate DNA codeword using DNA linear
block codes which ensures the secure transmission of information. In the
proposed code design strategy DNA-based XOR operation (DNAX) is applied for
effective construction of DNA codewords which are quadruples generated over the
set of alphabets consisting of four DNA bases adenine, thymine, guanine, and
cytosine. By worked out examples we explain the use of generator matrix and
parity check matrix in encryption and decryption of coded data in the form of
short single stranded DNA sequences. The newly developed technique can detect
as well as correcting error in transmission of DNA codewords through biological
channels from sender to the intended receiver. Through DNA coding theory we are
expanding the paths towards data compression and error correction in the form
of DNA strands. This leads us towards a broader domain of DNA cryptography.


  A number of discrete time, finite population size models in genetics
describing the dynamics of allele frequencies are known to converge (subject to
suitable scaling) to a diffusion process in the infinite population limit,
termed the Wright-Fisher diffusion. In this article we show that the diffusion
is ergodic uniformly in the selection and mutation parameters, and that the
measures induced by the solution to the stochastic differential equation are
uniformly locally asymptotically normal. Subsequently these two results are
used to analyse the statistical properties of the Maximum Likelihood and
Bayesian estimators for the selection parameter, when both selection and
mutation are acting on the population. In particular, it is shown that these
estimators are uniformly over compact sets consistent, display uniform in the
selection parameter asymptotic normality and convergence of moments over
compact sets, and are asymptotically efficient for a suitable class of loss
functions.


  The adoption of deep learning techniques in genomics has been hindered by the
difficulty of mechanistically interpreting the models that these techniques
produce. In recent years, a variety of post-hoc attribution methods have been
proposed for addressing this neural network interpretability problem in the
context of gene regulation. Here we describe a complementary way of approaching
this problem. Our strategy is based on the observation that two large classes
of biophysical models of cis-regulatory mechanisms can be expressed as deep
neural networks in which nodes and weights have explicit physiochemical
interpretations. We also demonstrate how such biophysical networks can be
rapidly inferred, using modern deep learning frameworks, from the data produced
by certain types of massively parallel reporter assays (MPRAs). These results
suggest a scalable strategy for using MPRAs to systematically characterize the
biophysical basis of gene regulation in a wide range of biological contexts.
They also highlight gene regulation as a promising venue for the development of
scientifically interpretable approaches to deep learning.


  Properties of non-canonical GC Base Pairs and their relation with
mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement
of the Transient GC Wobble Base Pairs in the mechanisms of the mechanochemical
cleavage of DNA and epigenetic mechanisms with participation of
5-methylcytosine is proposed. The hypothesis explains the increase in the
frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines,
asymmetric character of these breaks, and an increase in the frequency of
breaks in CpG after cytosine methylation.


  The mechanisms and properties of synchronization of oscillating ecological
populations attract attention because it is a fairly common phenomenon and
because spatial synchrony may elevate a risk of extinction and may lead to
other environmental impacts. Conditions for stable synchronization in a system
of linearly coupled predator-prey oscillators have been considered in the past.
However, the spatial dispersion coupling may be relatively weak and may not
necessarily lead to a stable, complete synchrony. If the coupling between
oscillators is too weak to induce a stable synchrony, oscillators may be
engaged into intermittent synchrony, when episodes of synchronized dynamics are
interspersed with the episodes of desynchronized dynamics. In the present study
we consider the temporal patterning of this kind of intermittent synchronized
dynamics in a system of two dispersal-coupled Rosenzweig-MacArthur
predator-prey oscillators. We consider the properties of the distributions of
durations of desynchronized intervals and their dependence on the model
parameters. We show that the temporal patterning of synchronous dynamics (an
ecological network phenomenon) may depend on the properties of individual
predator-prey patch (individual oscillator) and may vary independently of the
strength of dispersal. We also show that if the dynamics of predator is slow
relative to the dynamics of the prey (a situation that may promote brief but
large outbreaks), dispersal-coupled predator-prey oscillating populations
exhibit numerous short desynchronizations (as opposed to few long
desynchronizations) and may require weaker dispersal in order to reach strong
synchrony.


  Radio-frequency identification (RFID) is an increasingly popular wireless
technology that allows researchers to monitor wild bird populations from fixed
locations in the field. Our lab has developed an RFID-equipped birdfeeder based
on the Raspberry Pi Zero W, a low-cost single-board computer, that collects
continuous visitation data from birds tagged with passive integrated
transponder (PIT) tags. Each birdfeeder has a perch antenna connected to an
RFID reader board on a Raspberry Pi powered by a portable battery. When a
tagged bird lands on the perch to eat from the feeder, its unique code is
stored with the date and time on the Raspberry Pi. These birdfeeders require
only basic soldering and coding skills to assemble, and can be easily outfitted
with additional hardware like video cameras and microphones. We outline the
process of assembling the hardware and setting up the operating system for the
birdfeeders. Then, we describe an example implementation of the birdfeeders to
track house finches (Haemorhous mexicanus) on the campus of Queens College in
New York City.


  The inherent behavioral variability exhibited by stochastic biochemical
systems makes it a challenging task for human experts to manually analyze them.
Computational modeling of such systems helps in investigating and predicting
the behaviors of the underlying biochemical processes but at the same time
introduces the presence of several unknown parameters. A key challenge faced in
this scenario is to determine the values of these unknown parameters against
known behavioral specifications. The solutions that have been presented so far
estimate the parameters of a given model against a single specification whereas
a correct model is expected to satisfy all the behavioral specifications when
instantiated with a single set of parameter values. We present a new method,
BioMETA, to address this problem such that a single set of parameter values
causes a parameterized stochastic biochemical model to satisfy all the given
probabilistic temporal logic behavioral specifications simultaneously. Our
method is based on combining a multiple hypothesis testing based statistical
model checking technique with simulated annealing search to look for a single
set of parameter values so that the given parameterized model satisfies
multiple probabilistic behavioral specifications. We study two stochastic
rule-based models of biochemical receptors, namely, Fc$\epsilon$RI and T-cell
as our benchmarks to evaluate the usefulness of the presented method. Our
experimental results successfully estimate $26$ parameters of Fc$\epsilon$RI
and $29$ parameters of T-cell receptor model against three probabilistic
temporal logic behavioral specifications each.


  Spatio-temporal analysis typically performed in botany and in statistical
physics reveals persistent correlations of olive yields which range depend on
the size of the averaging region. Mapping spatially correlated regions reveals
areas which mimic historical spread of Olea europaea. These yield patterns are
remarkable given the intensive nature of modern agriculture, and cannot be
attributed to weather due to inability of weather variables to predict yields.
Long-range correlations between olive trees may indicate long-range
communications.


  Recently, Mikheyev et al. (2019) have produced a preprint study describing
the genomes of nine Khazars archeologically dated from the 7th to the 9th
centuries found in the Rostov county in modern-day Russia. Skull morphology
indicated a mix of "Caucasoid" and "Mongoloid" shapes. The authors compared the
samples to ancient and contemporary samples to study the genetic makeup of the
Khazars and their genetic legacy and addressed the question of the
relationships between the Khazar and Ashkenazic Jews. A careful examination
reveals grave concerns regarding all the aspects of the study from the
identification of the "Khazar" samples, the choice of environment for ancient
DNA sequencing, and the analyses. The authors did not disclose the data used in
their study, and their methodology is incoherent. We demonstrate that their
analyses yield nonsensical results and argue that none of the claims made in
this study are supported by the data unequivocally. Provided the destruction of
the bone samples and the irreproducibility of the analyses, even by the
forgivable standards of the field, this study is irreplicable, wasteful, and
misleading.


  A key question in the current diversity crisis is how diversity has been
maintained throughout evolution and how to preserve it. Modern coexistence
theories suggest that a high invasion rate of rare new types is directly
related to diversity. We show that adding almost any mechanism of catastrophes
to a stochastic birth, death, and mutation process with limited carrying
capacity induces a novel phase transition characterized by a positive invasion
rate but a low diversity. In this phase, new types emerge and grow rapidly, but
the resulting growth of very large types decreases diversity. This model also
resolves two major drawbacks of neutral evolution models: their failure to
explain balancing selection without resorting to fitness differences and the
unrealistic time required for the creation of the observed large types. We test
this model on a classical case of genetic polymorphism: the HLA locus.


  Motivation: Predicting the secondary structure of an RNA sequence is useful
in many applications. Existing algorithms (based on dynamic programming) suffer
from a major limitation: their runtimes scale cubically with the RNA length,
and this slowness limits their use in genome-wide applications.
  Results: We present a novel alternative $O(n^3)$-time dynamic programming
algorithm for RNA folding that is amenable to heuristics that make it run in
$O(n)$ time and $O(n)$ space, while producing a high-quality approximation to
the optimal solution. Inspired by incremental parsing for context-free grammars
in computational linguistics, our alternative dynamic programming algorithm
scans the sequence in a left-to-right (5'-to-3') direction rather than in a
bottom-up fashion, which allows us to employ the effective beam pruning
heuristic. Our work, though inexact, is the first RNA folding algorithm to
achieve linear runtime (and linear space) without imposing constraints on the
output structure. Surprisingly, our approximate search results in even higher
overall accuracy on a diverse database of sequences with known structures. More
interestingly, it leads to significantly more accurate predictions on the
longest sequence families in that database (16S and 23S Ribosomal RNAs), as
well as improved accuracies for long-range base pairs (500+ nucleotides apart),
both of which are well known to be challenging for the current models.
  Availability: Our source code is available at
https://github.com/LinearFold/LinearFold, and our webserver is at
http://linearfold.org (sequence limit: 100,000nt).


  Coarse-graining of fully atomistic molecular dynamics simulations is a
long-standing goal in order to allow the description of processes occurring on
biologically relevant timescales. For example, the prediction of pathways,
rates and rate-limiting steps in protein-ligand unbinding is crucial for modern
drug discovery. To achieve the enhanced sampling, we first perform
dissipation-corrected targeted molecular dynamics simulations, which yield free
energy and friction profiles of the molecular process under consideration. In a
second step, we use these fields to perform temperature-boosted Langevin
simulations which account for the desired molecular kinetics occurring on
multisecond timescales and beyond. Adopting the dissociation of solvated sodium
chloride as well as trypsin-benzamidine and Hsp90-inhibitor protein-ligand
complexes as test problems, we are able to reproduce rates from molecular
dynamics simulation and experiments within a factor of 2-20, and dissociation
constants within a factor of 1-4. Analysis of the friction profiles reveals
that binding and unbinding dynamics are mediated by changes of the surrounding
hydration shells in all investigated systems.


  We derive the asymptotic behavior of the total, active and inactive branch
lengths of the seed bank coalescent, when the size of the initial sample grows
to infinity. Those random variables have important applications for populations
evolving under some seed bank effects, such as plants and bacteria, and for
some cases of structured populations like metapopulations. The proof relies on
the study of the tree at a stopping time corresponding to the first time that a
deactivated lineage reactivates. We also give conditional sampling formulas for
the random partition and we study the system at the time of the first
deactivation of a lineage. All these results provide a good picture of the
different regimes and behaviors of the block-counting process of the seed bank
coalescent.


  In this work, a machine learning approach for identifying the multi-omics
metabolic regulatory control circuits inside the pathways is described.
Therefore, the identification of bacterial metabolic pathways that are more
regulated than others in term of their multi-omics follows from the analysis of
these circuits . This is a consequence of the alternation of the omic values of
codon usage and protein abundance along with the circuits. In this work, the
E.Coli's Glycolysis and its multi-omic circuit features are shown as an
example.


  Stochastic epidemic models (SEMs) fit to incidence data are critical to
elucidating outbreak dynamics, shaping response strategies, and preparing for
future epidemics. SEMs typically represent counts of individuals in discrete
infection states using Markov jump processes (MJPs), but are computationally
challenging as imperfect surveillance, lack of subject-level information, and
temporal coarseness of the data obscure the true epidemic. Analytic integration
over the latent epidemic process is impossible, and integration via Markov
chain Monte Carlo (MCMC) is cumbersome due to the dimensionality and
discreteness of the latent state space. Simulation-based computational
approaches can address the intractability of the MJP likelihood, but are
numerically fragile and prohibitively expensive for complex models. A linear
noise approximation (LNA) that approximates the MJP transition density with a
Gaussian density has been explored for analyzing prevalence data in
large-population settings, but requires modification for analyzing incidence
counts without assuming that the data are normally distributed. We demonstrate
how to reparameterize SEMs to appropriately analyze incidence data, and fold
the LNA into a data augmentation MCMC framework that outperforms deterministic
methods, statistically, and simulation-based methods, computationally. Our
framework is computationally robust when the model dynamics are complex and
applies to a broad class of SEMs. We evaluate our method in simulations that
reflect Ebola, influenza, and SARS-CoV-2 dynamics, and apply our method to
national surveillance counts from the 2013--2015 West Africa Ebola outbreak.


  Noise-modulating chemicals can synergize with transcriptional activators in
reactivating latent HIV to eliminate latent HIV reservoirs. To understand the
underlying biomolecular mechanism, we investigate a previous two-gene-state
model and identify two necessary conditions for the synergy: an assumption of
inhibition effect of transcription activators on noise enhancers; and frequent
transitions to the gene non-transcription-permissive state. We then develop a
loop-four-gene-state model with Tat transcription/translation and find that
drug synergy is mainly determined by the magnitude and direction of energy
input into the genetic regulatory kinetics of the HIV promoter. The inhibition
effect of transcription activators is actually a phenomenon of energy
dissipation in the nonequilibrium gene transition system. Overall, the
loop-four-state model demonstrates that energy dissipation plays a crucial role
in HIV latency reactivation, which might be useful for improving drug effects
and identifying other synergies on lentivirus latency reactivation.


  In eukaryotic cells, KDEL receptors (KDELRs) facilitate the retrieval of
endoplasmic reticulum (ER) luminal proteins from the Golgi compartment back to
the ER. Apart from the well-documented retention function, recent findings
reveal that the cellular KDELRs have more complex roles, e.g. in cell
signalling, protein secretion, cell adhesion and tumorigenesis. Furthermore,
several studies suggest that a sub-population of KDELRs is located at the cell
surface, where they could form and internalize KDELR/cargo clusters after
K/HDEL-ligand binding. However, so far it has been unclear whether there are
cell-type- or species-specific differences in KDELR clustering. By comparing
ligand-induced KDELR clustering in different mouse and human cell lines via
live cell imaging, we show that macrophage cell lines from both species do not
develop any clusters. Using RT-qPCR experiments and numerical analysis, we
address the role of KDELR expression as well as endocytosis and exocytosis
rates on the receptor clustering at the plasma membrane and discuss how the
efficiency of directed transport to preferred docking sites on the membrane
influences the exponent of the power-law distribution of the cluster size.


  A hybrid network--deterministic model for simulation of multiresistant
pathogen spread in a healthcare system is presented. The model accounts for two
paths of pathogen transmission between the healthcare facilities:
inter-hospital patient transfers (direct transfers) and readmission of
colonized patients (indirect transfers). In the latter case, the patients may
be readmitted to the same facility or to a different one. Intra-hospital
pathogen transmission is governed by a SIS model expressed by a system of
ordinary differential equations.
  Using a network model created for a Lower Saxony region (Germany), we showed
that the proposed model reproduces the basic properties of
healthcare-associated pathogen spread. Moreover, it shows the important
contribution of the readmission of colonized patients on the prevalence of
individual hospitals as well as of overall healthcare system: it can increase
the overall prevalence by the factor of 4 as compared to inter-hospital
transfers only. The final prevalence in individual healthcare facilities was
shown to depend on average length of stay by a non-linear concave function.
  Finally, we demonstrated that the network parameters of the model may be
derived from administrative admission/discharge records. In particular, they
are sufficient to obtain inter-hospital transfer probabilities, and to express
the patients' transfer as a Markov process.


  A tree-based network $N$ on $X$ is called universal if every phylogenetic
tree on $X$ is a base tree for $N$. Recently, binary universal tree-based
networks have attracted great attention in the literature and their existence
has been analyzed in various studies. In this note, we extend the analysis to
non-binary networks and show that there exist both a rooted and an unrooted
non-binary universal tree-based network with $n$ leaves for all positive
integers $n$.


  A major challenge in biomedical data science is to identify the causal genes
underlying complex genetic diseases. Despite the massive influx of genome
sequencing data, identifying disease-relevant genes remains difficult as
individuals with the same disease may share very few, if any, genetic variants.
Protein-protein interaction networks provide a means to tackle this
heterogeneity, as genes causing the same disease tend to be proximal within
networks. Previously, network propagation approaches have spread signal across
the network from either known disease genes or genes that are newly putatively
implicated in the disease (e.g., found to be mutated in exome studies or linked
via genome-wide association studies). Here we introduce a general framework
that considers both sources of data within a network context. Specifically, we
use prior knowledge of disease-associated genes to guide random walks initiated
from genes that are newly identified as perhaps disease-relevant. In
large-scale testing across 24 cancer types, we demonstrate that our approach
for integrating both prior and new information not only better identifies
cancer driver genes than using either source of information alone but also
readily outperforms other state-of-the-art network-based approaches. To
demonstrate the versatility of our approach, we also apply it to genome-wide
association data to identify genes functionally relevant for several complex
diseases. Overall, our work suggests that guided network propagation approaches
that utilize both prior and new data are a powerful means to identify disease
genes.


  Liquid-liquid phase separation occurs not only in bulk liquid, but also on
surfaces. In physiology, the nature and function of condensates on cellular
structures remain unexplored. Here, we study how the condensed protein TPX2
behaves on microtubules to initiate branching microtubule nucleation, which is
critical for spindle assembly in eukaryotic cells. Using fluorescence,
electron, and atomic force microscopies and hydrodynamic theory, we show that
TPX2 on a microtubule reorganizes according to the Rayleigh-Plateau
instability, like dew droplets patterning a spider web. After uniformly coating
microtubules, TPX2 forms regularly spaced droplets from which branches
nucleate. Droplet spacing increases with greater TPX2 concentration. A
stochastic model shows that droplets make branching nucleation more efficient
by confining the space along the microtubule where multiple necessary factors
colocalize to nucleate a branch.


  In the last decade, machine learning and artificial intelligence applications
have received a significant boost in performance and attention in both academic
research and industry. The success behind most of the recent state-of-the-art
methods can be attributed to the latest developments in deep learning. When
applied to various scientific domains that are concerned with the processing of
non-tabular data, for example, image or text, deep learning has been shown to
outperform not only conventional machine learning but also highly specialized
tools developed by domain experts. This review aims to summarize AI-based
research for GPCR bioactive ligand discovery with a particular focus on the
most recent achievements and research trends. To make this article accessible
to a broad audience of computational scientists, we provide instructive
explanations of the underlying methodology, including overviews of the most
commonly used deep learning architectures and feature representations of
molecular data. We highlight the latest AI-based research that has led to the
successful discovery of GPCR bioactive ligands. However, an equal focus of this
review is on the discussion of machine learning-based technology that has been
applied to ligand discovery in general and has the potential to pave the way
for successful GPCR bioactive ligand discovery in the future. This review
concludes with a brief outlook highlighting the recent research trends in deep
learning, such as active learning and semi-supervised learning, which have
great potential for advancing bioactive ligand discovery.


  The multispecies coalescent process models the genealogical relationships of
genes sampled from several species, enabling useful predictions about phenomena
such as the discordance between the gene tree and the species phylogeny due to
incomplete lineage sorting. Conversely, knowledge of large collections of gene
trees can inform us about several aspects of the species phylogeny, such as its
topology and ancestral population sizes. A fundamental open problem in this
context is how to efficiently compute the probability of a gene tree topology,
given the species phylogeny. Although a number of algorithms for this task have
been proposed, they either produce approximate results, or, when they are
exact, they do not scale to large data sets. In this paper, we present some
progress towards exact and efficient computation of the probability of a gene
tree topology. We provide a new algorithm that, given a species tree and the
number of genes sampled for each species, calculates the probability that the
gene tree topology will be concordant with the species tree. Moreover, we
provide an algorithm that computes the probability of any specific gene tree
topology concordant with the species tree. Both algorithms run in polynomial
time and have been implemented in Python. Experiments show that they are able
to analyse data sets where thousands of genes are sampled, in a matter of
minutes to hours.


  In this paper, we study the psychological effect in a SIS epidemic model. The
basic reproduction number is obtained. However, the disease free equilibrium is
always asymptotically stable, which doesn't depends on the basic reproduction
number. The system has a saddle-node bifurcation appear and displays bistable
behavior, which is a new phenomenon in epidemic dynamics and different from the
backward bifurcation behavior.


  Stochasticity (or noise) at cellular and molecular levels has been observed
extensively as a universal feature for living systems. However, how living
systems deal with noise while performing desirable biological functions remains
a major mystery. Regulatory network configurations, such as their topology and
timescale, are shown to be critical in attenuating noise, and noise is also
found to facilitate cell fate decision. Here we review major recent findings on
noise attenuation through regulatory control, the benefit of noise via
noise-induced cellular plasticity during developmental patterning, and
summarize key principles underlying noise control.


  Humans, domestic animals, orchard crops, and ornamental plants are commonly
treated with antibiotics in response to bacterial infection. By curing
infectious individuals, antibiotic therapy might limit the spread of contagious
disease among hosts. But antibiotic suppression of within-host pathogen density
might also reduce the probability that the host is otherwise removed from
infectious status before recovery. When rates of both recovery and removal
(isolation or mortality) depend directly on within-host density, antibiotic
therapy can relax the removal rate and so increase between-host disease
transmission. In this paper a deterministic within-host dynamics drives the
infectious host's probability of infection transmission, as well as the host's
time-dependent probability of surviving to recovery. The model varies (1)
inoculum size, (2) the time elapsing between infection and initiation of
therapy, (2) antibiotic efficacy, and (3) the size/susceptibility of groups
encountered by an infectious host. Results identify conditions where antibiotic
treatment simultaneously increases host survival and increases the expected
number of new infections. That is, antibiotics might convert a rare, serious
bacterial disease into a common, but treatable infection.


  Biosyntheses of proteins, nucleotides and fatty acids, are essential for the
malignant proliferation and survival of cancer cells. Cumulating research
findings show that amino acid restrictions are potential strategies for cancer
interventions. Meanwhile, dietary strategies are popular among cancer patients.
However, there is still lacking solid rationale to clarify what is the best
strategy, why and how it is. Here, integrated analyses and comprehensive
summaries for the abundances, signalling and functions of amino acids in
proteomes, metabolism, immunity and food compositions, suggest that,
intermittent fasting or intermittent dietary lysine restriction with normal
maize as an intermittent staple food for days or weeks, might have the value
and potential for cancer prevention or therapy. Moreover, dietary supplements
were also discussed for cancer cachexia including dietary immunomodulatory.


  Genomic data sets are growing dramatically as the cost of sequencing
continues to decline and small sequencing devices become available. Enormous
community databases store and share this data with the research community, but
some of these genomic data analysis problems require large scale computational
platforms to meet both the memory and computational requirements. These
applications differ from scientific simulations that dominate the workload on
high end parallel systems today and place different requirements on programming
support, software libraries, and parallel architectural design. For example,
they involve irregular communication patterns such as asynchronous updates to
shared data structures. We consider several problems in high performance
genomics analysis, including alignment, profiling, clustering, and assembly for
both single genomes and metagenomes. We identify some of the common
computational patterns or motifs that help inform parallelization strategies
and compare our motifs to some of the established lists, arguing that at least
two key patterns, sorting and hashing, are missing.


  We have developed a method of animal localisation that detects the angle from
a sensor towards the direction of an animal call. The method is as simple to
use as deploying a conventional static sound recorder, but provides tracking
information as well as sound recordings. The principal of operation is to
detect the phase difference between microphones positioned closely together.
The phase is detected by converting the signals to their analytic form with a
Hilbert transform. The angle is then calculated from the phase difference,
frequency and microphone separation. Angular measurements provide flight paths
above the sensor, and can give details of activity and behaviour that are not
possible with a single channel static recorder. We recorded flight paths for 5
bat species on a single night at a site in Denmark (Pipistrellus nathusii,
Pipistrellus pygmaeus, Eptesicus serotinus, Myotis daubentonii and Nyctalus
noctula). The median error in angular measurement for the species was between 3
and 7 degrees. Calls at high angles from normal, corresponding with a poor
signal-to-noise ratio, had larger errors compared to calls recorded in the
centre of the field of view. Locations in space could be estimated by combining
angular measurements from two or more sensors.


  Biological networks exhibit intricate architectures deemed to be crucial for
their functionality. In particular, gene regulatory networks, which play a key
role in information processing in the cell, display non-trivial architectural
features such as scale-free degree distributions, high modularity, and low
average distance between connected genes. Such networks result from complex
evolutionary and adaptive processes difficult to track down empirically. On the
other hand, there exists detailed information on the developmental (or
evolutionary) stages of open-software networks that result from self-organized
growth across versions. Here, we study the evolution of the Debian GNU/Linux
software network, focusing on the changes of key structural and statistical
features over time. Our results show that evolution has led to a network
structure in which the out-degree distribution is scale-free and the in-degree
distribution is a stretched exponential. In addition, while modularity,
directionality of information flow, and average distance between elements grew,
vulnerability decreased over time. These features resemble closely those
currently shown by gene regulatory networks, suggesting the existence of common
adaptive pathways for the architectural design of information-processing
networks. Differences in other hierarchical aspects point to system-specific
solutions to similar evolutionary challenges.


  Quantitative predictions about the processes that promote species coexistence
are a subject of active research in ecology. In particular, competitive
interactions are known to shape and maintain ecological communities, and
situations where some species out-compete or dominate over some others are key
to describe natural ecosystems. Here we develop ecological theory using a
stochastic, synthetic framework for plant community assembly that leads to
predictions amenable to empirical testing. We propose two stochastic
continuous-time Markov models that incorporate competitive dominance through a
hierarchy of species heights. The first model, which is spatially implicit,
predicts both the expected number of species that survive and the conditions
under which heights are clustered in realized model communities. The second one
allows spatially-explicit interactions of individuals and alternative
mechanisms that can help shorter plants overcome height-driven competition, and
it demonstrates that clustering patterns remain not only locally but also
across increasing spatial scales. Moreover, although plants are actually
height-clustered in the spatially-explicit model, it allows for plant species
abundances not necessarily skewed to taller plants.


  The phenomenon of molecular binding, where two molecules, referred to as a
receptor and a ligand, bind together to form a ligand-receptor complex, is
ubiquitous in biology and essential for the accurate functioning of all
life-sustaining processes. The probability of a single receptor forming a
complex with any one of L surrounding ligand molecules at thermal equilibrium
can be derived from a partition function obtained from the Gibbs-Boltzmann
distribution. We extend this approach to a system consisting of R receptors and
L ligands to derive the probability density function p(r;R,L) to find r bound
receptor-ligand complexes at thermal equilibrium. This extension allows us to
illustrate two aspects of this problem which are not apparent in the single
receptor problem, namely, a) a symmetry to be expected in the equilibrium
distribution of the number of bound complexes under exchange of R and L and b)
the number of bound complexes obtained from chemical kinetic equations has an
exact correspondence to the maximum probable value of r from the expression for
p(r;R,L). We derive the number fluctuations of r and present a practically
relevant molecular sensing application that benefits from the knowledge of
p(r;R,L).


  Advances in DNA sequencing have revolutionized our ability to read genomes.
However, even in the most well-studied of organisms, the bacterium ${\it
Escherichia coli}$, for $\approx$ 65$\%$ of the promoters we remain completely
ignorant of their regulation. Until we have cracked this regulatory Rosetta
Stone, efforts to read and write genomes will remain haphazard. We introduce a
new method (Reg-Seq) linking a massively-parallel reporter assay and mass
spectrometry to produce a base pair resolution dissection of more than 100
promoters in ${\it E. coli}$ in 12 different growth conditions. First, we show
that our method recapitulates regulatory information from known sequences.
Then, we examine the regulatory architectures for more than 80 promoters in the
${\it E. coli}$ genome which previously had no known regulation. In many cases,
we also identify which transcription factors mediate their regulation. The
method introduced here clears a path for fully characterizing the regulatory
genome of model organisms, with the potential of moving on to an array of other
microbes of ecological and medical relevance.


  The emergence of acquired drug resistance in cancer represents a major
barrier to treatment success. While research has traditionally focused on
genetic sources of resistance, recent findings suggest that cancer cells can
acquire transient resistant phenotypes via epigenetic modifications and other
non-genetic mechanisms. Although these resistant phenotypes are eventually
relinquished by individual cells, they can temporarily 'save' the tumor from
extinction and enable the emergence of more permanent resistance mechanisms.
These observations have generated interest in the potential of epigenetic
therapies for long-term tumor control or eradication. In this work, we develop
a mathematical model to study how phenotypic switching at the single-cell level
affects resistance evolution in cancer. We highlight unique features of
non-genetic resistance, probe the evolutionary consequences of epigenetic drugs
and explore potential therapeutic strategies. We find that even short-term
epigenetic modifications and stochastic fluctuations in gene expression can
drive long-term drug resistance in the absence of any bona fide resistance
mechanisms. We also find that an epigenetic drug that slightly perturbs the
average retention of the resistant phenotype can turn guaranteed treatment
failure into guaranteed success. Lastly, we find that combining an epigenetic
drug with an anti-cancer agent can significantly outperform monotherapy, and
that treatment outcome is heavily affected by drug sequencing.


  Magombedze and Mulder in 2013 studied the gene regulatory system of
Mycobacterium Tuberculosis (Mtb) by partitioning this into three subsystems
based on putative gene function and role in dormancy/latency development. Each
subsystem, in the form of S-system, is represented by an embedded chemical
reaction network (CRN), defined by a species subset and a reaction subset
induced by the set of digraph vertices of the subsystem. For the embedded
networks of S-system, we showed interesting structural properties and proved
that all S-system CRNs (with at least two species) are discordant. Analyzing
the subsystems as subnetworks, where arcs between vertices belonging to
different subsystems are retained, we formed a digraph homomorphism from the
corresponding subnetworks to the embedded networks. Lastly, we explored the
modularity concept of CRN in the context of digraph.


  The progesterone receptor (PR) mediates progesterone regulation of female
reproductive physiology, as well as gene transcription in non-reproductive
tissues, such as brain, bone, lung and vasculature, in both women and men. An
unusual property of progesterone is its high affinity for the mineralocorticoid
receptor (MR), which regulates electrolyte transport in the kidney in humans
and other terrestrial vertebrates. In humans, rats, alligators and frogs,
progesterone antagonizes activation of the MR by aldosterone, the physiological
mineralocorticoid in terrestrial vertebrates. In contrast, in elephant shark,
ray-finned fishes and chickens, progesterone activates the MR. Interestingly,
cartilaginous fishes and ray-finned fishes do not synthesize aldosterone,
raising the question of which steroid(s) activate the MR in cartilaginous
fishes and ray-finned fishes. The simpler synthesis of progesterone, compared
to cortisol and other corticosteroids, makes progesterone a candidate
physiological activator of the MR in elephant sharks and ray-finned fishes.
Elephant shark and ray-finned fish MRs are expressed in diverse tissues,
including heart, brain and lung, as well as, ovary and testis, two reproductive
tissues that are targets for progesterone, which together suggests a
multi-faceted physiological role for progesterone activation of the MR in
elephant shark and ray-finned fish. The functional consequences of progesterone
as an antagonist of some terrestrial vertebrate MRs and as an agonist of fish
and chicken MRs are not fully understood. Indeed, little is known of
physiological activities of progesterone via any vertebrate MR.


  Single Particle Tracking (SPT) is a powerful class of tools for analyzing the
dynamics of individual biological macromolecules moving inside living cells.
The acquired data is typically in the form of a sequence of camera images that
are then post-processed to reveal details about the motion. In this work, we
develop an algorithm for jointly estimating both particle trajectory and motion
model parameters from the data. Our approach uses Expectation Maximization (EM)
combined with an Unscented Kalman filter (UKF) and an Unscented
Rauch-Tung-Striebel smoother (URTSS), allowing us to use an accurate, nonlinear
model of the observations acquired by the camera. Due to the shot noise
characteristics of the photon generation process, this model uses a Poisson
distribution to capture the measurement noise inherent in imaging. In order to
apply a UKF, we first must transform the measurements into a model with
additive Gaussian noise. We consider two approaches, one based on variance
stabilizing transformations (where we compare the Anscombe and Freeman-Tukey
transforms) and one on a Gaussian approximation to the Poisson distribution.
Through simulations, we demonstrate efficacy of the approach and explore the
differences among these measurement transformations.


  Species tree estimation from multi-locus datasets is statistically
challenging for multiple reasons, including gene tree heterogeneity across the
genome due to incomplete lineage sorting (ILS). Species tree estimation methods
have been developed that operate by estimating gene trees and then using those
gene trees to estimate the species tree. Several of these methods (e.g.,
ASTRAL, ASTRID, and NJst) are provably statistically consistent under the
multi-species coalescent (MSC) model, provided that the gene trees are
estimated correctly, and there is no missing data. Recently, Nute et al. (BMC
Genomics 2018) addressed the question of whether these methods remain
statistically consistent under random models of taxon deletion, and asserted
that they do so. Here we provide a counterexample to one of these theorems, and
establish that ASTRID and NJst are not statistically consistent under an i.i.d.
model of taxon deletion.


  Active inference is a normative framework for explaining behaviour under the
free energy principle -- a theory of self-organisation originating in
neuroscience. It specifies neuronal dynamics for state-estimation in terms of a
descent on (variational) free energy -- a measure of the fit between an
internal (generative) model and sensory observations. The free energy gradient
is a prediction error -- plausibly encoded in the average membrane potentials
of neuronal populations. Conversely, the expected probability of a state can be
expressed in terms of neuronal firing rates. We show that this is consistent
with current models of neuronal dynamics and establish face validity by
synthesising plausible electrophysiological responses. We then show that these
neuronal dynamics approximate natural gradient descent, a well-known
optimisation algorithm from information geometry that follows the steepest
descent of the objective in information space. We compare the information
length of belief updating in both schemes, a measure of the distance traveled
in information space that has a direct interpretation in terms of metabolic
cost. We show that neural dynamics under active inference are metabolically
efficient and suggest that neural representations in biological agents may
evolve by approximating steepest descent in information space towards the point
of optimal inference.


  Multistationarity in molecular systems underlies switch-like responses in
cellular decision making. Determining whether and when a system displays
multistationarity is in general a difficult problem. In this work we completely
determine the set of kinetic parameters that enable multistationarity in a
ubiquitous motif involved in cell signaling, namely a dual phosphorylation
cycle. In addition we show that the regions of multistationarity and
monostationarity are both path connected.
  We model the dynamics of the concentrations of the proteins over time by
means of a parametrized polynomial ordinary differential equation (ODE) system
arising from the mass-action assumption. Since this system has three linear
first integrals defined by the total amounts of the substrate and the two
enzymes, we study for what parameter values the ODE system has at least two
positive steady states after suitably choosing the total amounts. We employ a
suite of techniques from (real) algebraic geometry, which in particular concern
the study of the signs of a multivariate polynomial over the positive orthant
and sums of nonnegative circuit polynomials.


  When analyzing communities of microorganisms from their sequenced DNA, an
important task is taxonomic profiling: enumerating the presence and relative
abundance of all organisms, or merely of all taxa, contained in the sample.
This task can be tackled via compressive-sensing-based approaches, which favor
communities featuring the fewest organisms among those consistent with the
observed DNA data. Despite their successes, these parsimonious approaches
sometimes conflict with biological realism by overlooking organism
similarities. Here, we leverage a recently developed notion of biological
diversity that simultaneously accounts for organism similarities and retains
the optimization strategy underlying compressive-sensing-based approaches. We
demonstrate that minimizing biological diversity still produces sparse
taxonomic profiles and we experimentally validate superiority to existing
compressive-sensing-based approaches. Despite showing that the objective
function is almost never convex and often concave, generally yielding NP-hard
problems, we exhibit ways of representing organism similarities for which
minimizing diversity can be performed via a sequence of linear programs
guaranteed to decrease diversity. Better yet, when biological similarity is
quantified by $k$-mer co-occurrence (a popular notion in bioinformatics),
minimizing diversity actually reduces to one linear program that can utilize
multiple $k$-mer sizes to enhance performance. In proof-of-concept experiments,
we verify that the latter procedure can lead to significant gains when
taxonomically profiling a metagenomic sample, both in terms of reconstruction
accuracy and computational performance. Reproducible code is available at
https://github.com/dkoslicki/MinimizeBiologicalDiversity.


  The TADPOLE Challenge compares the performance of algorithms at predicting
the future evolution of individuals at risk of Alzheimer's disease. TADPOLE
Challenge participants train their models and algorithms on historical data
from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. Participants
are then required to make forecasts of three key outcomes for ADNI-3 rollover
participants: clinical diagnosis, ADAS-Cog 13, and total volume of the
ventricles -- which are then compared with future measurements. Strong points
of the challenge are that the test data did not exist at the time of
forecasting (it was acquired afterwards), and that it focuses on the
challenging problem of cohort selection for clinical trials by identifying fast
progressors. The submission phase of TADPOLE was open until 15 November 2017;
since then data has been acquired until April 2019 from 219 subjects with 223
clinical visits and 150 Magnetic Resonance Imaging (MRI) scans, which was used
for the evaluation of the participants' predictions. Thirty-three teams
participated with a total of 92 submissions. No single submission was best at
predicting all three outcomes. For diagnosis prediction, the best forecast
(team Frog), which was based on gradient boosting, obtained a multiclass area
under the receiver-operating curve (MAUC) of 0.931, while for ventricle
prediction the best forecast (team EMC1), which was based on disease
progression modelling and spline regression, obtained mean absolute error of
0.41% of total intracranial volume (ICV). For ADAS-Cog 13, no forecast was
considerably better than the benchmark mixed effects model (BenchmarkME),
provided to participants before the submission deadline. Further analysis can
help understand which input features and algorithms are most suitable for
Alzheimer's disease prediction and for aiding patient stratification in
clinical trials.


  Dynamics of nucleosomes, the building blocks of the chromatin, has crucial
effects on expression, replication and repair of genomes in eukaryotes. Beside
constant movements of nucleosomes by thermal fluctuations, ATP-dependent
chromatin remodelling complexes cause their active displacements. Here we
propose a theoretical analysis of dinucleosome wrapping and unwrapping dynamics
in the presence of an external force. We explore the energy landscape and
configurations of dinucleosome in different unwrapped states. Moreover, using a
dynamical Monte-Carlo simulation algorithm, we demonstrate the dynamical
features of the system such as the unwrapping force for partial and full
wrapping processes. Furthermore, we show that in the short length of linker DNA
($\sim 10 - 90$ bp), the asymmetric unwrapping occurs. These findings could
shed some light on chromatin dynamics and gene accessibility.


  Cancer disease occurs because of a disorder in the cellular regulatory
mechanism, Which causes cellular malformation. The genes that start the
malformation are called Cancer driver genes (CDGs) . Numerous computational
methods have been introduced to identify cancer driver genes that use the
concept of mutation.Regarding abnormalities spread in human cell and tumor
development, CDGs are likely to be the potential types of gene with high
influence in the network. This increases the importance of influence diffusion
concept for the identification of CDGs.recently developed a method based on
influence maximization for identifying cancer driver genes. One of the
challenges in these types of networks is to find the power of regulatory
interaction between edges.The current study developed a technique to identify
cancer driver gene and predict the impact of regulatory interactions in a
transcriptional regulatory network. This technique utilizes the concept of
influence diffusion and optimizes the Hyperlink-Induced Topic Search algorithm
based on the influence diffusion. The results suggest the better performance of
our proposed technique than the other computational and network-based
approaches.


  Cold atmospheric plasma (CAP) is a promising new agent for (selective) cancer
treatment, but the underlying cause of the anti-cancer effect of CAP is not
well understood yet. Among different theories and observations, one theory in
particular has been postulated in great detail and consists of a very complex
network of reactions that are claimed to account for the anti-cancer effect of
CAP. Here, the key concept is a reactivation of two specific apoptotic cell
signaling pathways through catalase inactivation caused by CAP. Thus, it is
postulated that the anti-cancer effect of CAP is due to its ability to
inactivate catalase, either directly or indirectly. A theoretical investigation
of the proposed theory, especially the role of catalase inactivation, can
contribute to the understanding of the underlying cause of the anti-cancer
effect of CAP. In the present study, we develop a mathematical model to analyze
the proposed catalase-dependent anti-cancer effect of CAP. Our results show
that a catalase-dependent reactivation of the two apoptotic pathways of
interest is unlikely to contribute to the observed anti-cancer effect of CAP.
Thus, we believe that other theories of the underlying cause should be
considered and evaluated to gain knowledge about the principles of CAP-induced
cancer cell death.


  The study of motor protein dynamics within cytoskeletal networks is of high
interest to physicists and biologists to understand how the dynamics and
properties of individual motors lead to cooperative effects and control of
overall network behaviour. Here, we report a method to detect and track
muscular myosin II filaments within an actin network tethered to supported
lipid bilayers. Based on the characteristic shape of myosin II filaments, this
automated tracking routine allowed us to follow the position and orientation of
myosin II filaments over time, and to reliably classify their dynamics into
segments of diffusive and processive motion based on the analysis of
displacements and angular changes between time steps. This automated, high
throughput method will allow scientists to efficiently analyse motor dynamics
in different conditions, and will grant access to more detailed information
than provided by common tracking methods, without any need for time consuming
manual tracking or generation of kymographs.


  Outer hair cell (OHC) nonlinear capacitance (NLC) represents voltage sensor
charge movements of prestin (SLC26a5), the protein responsible for OHC
electromotility. Previous measures of NLC frequency response have employed
methods which did not assess the influence of dielectric loss (sensor charge
movements out of phase with voltage) that may occur, and such loss conceivably
may influence the frequency dependent activity of prestin. Here we evaluate
complex capacitance of prestin out to 30 kHz and find that its frequency
response determined using this approach coincides with all previous estimates.
We also show that membrane tension has no effect on the frequency response of
prestin, despite substantial shifts in its voltage operating range, indicating
that prestin transition rate alterations do not account for the shifts. The
magnitude roll-off of prestin activity across frequency surpasses the
reductions of NLC caused by salicylate treatments that are known to abolish
cochlear amplification. Such roll-off must therefore limit the effectiveness if
prestin in contributing to cochlear amplification at the very high acoustic
frequencies processed by some mammals.


  We live and cooperate in networks. However, links in networks only allow for
pairwise interactions, thus making the framework suitable for dyadic games, but
not for games that are played in groups of more than two players. Here, we
study the evolutionary dynamics of a public goods game in social systems with
higher-order interactions. First, we show that the game on uniform hypergraphs
corresponds to the replicator dynamics in the well-mixed limit, providing a
formal theoretical foundation to study cooperation in networked groups.
Secondly, we unveil how the presence of hubs and the coexistence of
interactions in groups of different sizes affects the evolution of cooperation.
Finally, we apply the proposed framework to extract the actual dependence of
the synergy factor on the size of a group from real-world collaboration data in
science and technology. Our work provides a way to implement informed actions
to boost cooperation in social groups.


  The data of the 2009 A(H1N1) epidemic in Istanbul, Turkey is unique in terms
of the collected data, which include not only the hospitalization but also the
fatality information recorded during the pandemic. The analysis of this data
displayed an unexpected time shift between the hospital referrals and
fatalities. This time shift, which does not conform to the SIR and SEIR models,
was explained by multi-stage SIR and SEIR models [21]. In this study we prove
that the delay for these models is half of the infectious period within a
quadratic approximation, and we determine the epidemic parameters $R_0$, $T$
and $I_0$ of the 2009 A(H1N1) Istanbul and Netherlands epidemics.These epidemic
parameters were estimated by comparing the normalized cumulative fatality data
with the solutions of the SIR model. Two different error criteria, the $L_2$
norms of the error over the whole observation period and over the initial
portion of the data, were used in order to obtain the best-fitting models. It
was observed that, with respect to both criteria, the parameters of "good"
models were agglomerated along a line in the $T$-$R_0$ plane, instead of being
scattered uniformly around a "best" model. As this fact indicates the existence
of a nearly invariant quantity, interval estimates for the parameters were
given. As the initial phase of the epidemics were less influenced by the
effects of medical interventions, the error norm based on the initial portion
of the data was preferred. However, the presented parameter ranges are well out
of the range for the usual influenza epidemic parameter values. To confirm our
observations on the Istanbul data, the same error criteria were also used for
the 2009 A(H1N1) epidemic for the Netherlands, which has a similar population
density as in Istanbul. As in the Istanbul case, the parameter ranges do not
match the usual influenza epidemic parameter values.


  Objectives.--To estimate the basic reproduction number of the Wuhan novel
coronavirus (2019-nCoV). Methods.--Based on the
susceptible-exposed-infected-removed (SEIR) compartment model and the
assumption that the infectious cases with symptoms occurred before January 25,
2020 are resulted from free propagation without intervention, we estimate the
basic reproduction number of 2019-nCoV according to the reported confirmed
cases and suspected cases, as well as the theoretical estimated number of
infected cases by other research teams, together with some epidemiological
determinants learned from the severe acute respiratory syndrome. Results The
basic reproduction number falls between 2.8 to 3.3 by using the real-time
reports on the number of 2019-nCoV infected cases from People's Daily in China,
and falls between 3.2 and 3.9 on the basis of the predicted number of infected
cases from colleagues. Conclusions.--The early transmission ability of
2019-nCoV is closed to or slightly higher than SARS. It is a controllable
disease with moderate-high transmissibility. Timely and effective control
measures are needed to suppress the further transmissions. Notes Added.--Using
a newly reported epidemiological determinants for early 2019-nCoV, the
estimated basic reproduction number is in the range [2.2,3.0].