  Sexual contacts are the main spreading route of HIV. This puts sex workers at
higher risk of infection even in populations where HIV prevalence is moderate
or low. Alongside condom use, Pre-Exposure Prophylaxis (PrEP) is an effective
tool for sex workers to reduce their risk of HIV acquisition. However, PrEP
provides no direct protection against sexually transmitted infections (STIs)
other than HIV, unlike condoms. We use an empirical network of sexual contacts
among female sex workers (FSWs) and clients to simulate the spread of HIV and
gonorrhea. We then investigate the effect of PrEP adoption and adherence, on
both HIV and gonorrhea prevalence. We also study the effect of a potential
increase in condomless acts due to lowered risk perception with respect of the
no-PrEP scenario (risk compensation). We find that when HIV is the only disease
circulating, PrEP is effective in reducing HIV prevalence, even with high risk
compensation. Instead, the complex interplay between the two diseases shows
that different levels of risk compensation require different intervention
strategies. Finally, we find that providing PrEP only to the most active FSWs
is less effective than uniform PrEP adoption. Our work shows that the effects
emerging from the complex interactions between these diseases and the available
prophylactic measures need to be accounted for, to devise effective
intervention strategies.


  Using a stochastic individual-based modelling approach, we examine the role
that Delta-Notch signalling plays in the regulation of a robust and reliable
somite segmentation clock. We find that not only can Delta-Notch signalling
synchronise noisy cycles of gene expression in adjacent cells in the presomitic
mesoderm (as is known), but it can also amplify and increase the coherence of
these cycles. We examine some of the shortcomings of deterministic approaches
to modelling these cycles and demonstrate how intrinsic noise can play an
active role in promoting sustained oscillations, giving rise to noise-induced
quasi-cycles. Finally, we explore how translational/transcriptional delays can
result in the cycles in neighbouring cells oscillating in anti-phase and we
study how this effect relates to the propagation of noise-induced stochastic
waves.


  With a view towards molecular communication systems and molecular multi-agent
systems, we propose the Chemical Baum-Welch Algorithm, a novel reaction network
scheme that learns parameters for Hidden Markov Models (HMMs). Each reaction in
our scheme changes only one molecule of one species to one molecule of another.
The reverse change is also accessible but via a different set of enzymes, in a
design reminiscent of futile cycles in biochemical pathways. We show that every
fixed point of the Baum-Welch algorithm for HMMs is a fixed point of our
reaction network scheme, and every positive fixed point of our scheme is a
fixed point of the Baum-Welch algorithm. We prove that the "Expectation" step
and the "Maximization" step of our reaction network separately converge
exponentially fast. We simulate mass-action kinetics for our network on an
example sequence, and show that it learns the same parameters for the HMM as
the Baum-Welch algorithm.


  Developing efficient computational methods to assess the impact of external
interventions on the dynamics of a network model is an important problem in
systems biology. This paper focuses on quantifying the global changes that
result from the application of an intervention to produce a desired effect,
which we define as the total effect of the intervention. The type of
mathematical models that we will consider are discrete dynamical systems which
include the widely used Boolean networks and their generalizations. The
potential interventions can be represented by a set of nodes and edges that can
be manipulated to produce a desired effect on the system. We use a class of
regulatory rules called nested canalizing functions that frequently appear in
published models and were inspired by the concept of canalization in
evolutionary biology. In this paper, we provide a polynomial normal form based
on the canalizing properties of regulatory functions. Using this polynomial
normal form, we give a set of formulas for counting the maximum number of
transitions that will change in the state space upon an edge deletion in the
wiring diagram. These formulas rely on the canalizing structure of the target
function since the number of changed transitions depends on the canalizing
layer that includes the input to be deleted. We also present computations on
random networks to compare the exact number of changes with the upper bounds
provided by our formulas. Finally, we provide statistics on the sharpness of
these upper bounds in random networks.


  Motivated by recent experimental work, we define and study a deterministic
model of the complex miRNA-based regulatory circuit that putatively controls
the early stage of myogenesis in human. We aim in particular at a quantitative
understanding of (i) the roles played by the separate and independent miRNA
biosynthesis channels (one involving a miRNA-decoy system regulated by an
exogenous controller, the other given by transcription from a distinct genomic
locus) that appear to be crucial for the differentiation program, and of (ii)
how competition to bind miRNAs can efficiently control molecular levels in such
an interconnected architecture. We show that optimal static control via the
miRNA-decoy system constrains kinetic parameters in narrow ranges where the
channels are tightly cross-linked. On the other hand, the alternative locus for
miRNA transcription can ensure that the fast concentration shifts required by
the differentiation program are achieved, specifically via non-linear response
of the target to even modest surges in the miRNA transcription rate. While
static, competition-mediated regulation can be achieved by the miRNA-decoy
system alone, both channels are essential for the circuit's overall
functionality, suggesting that that this type of joint control may represent a
minimal optimal architecture in different contexts.


  Interacting proteins coevolve at multiple but interconnected scales, from the
residue-residue over the protein-protein up to the family-family level. The
recent accumulation of enormous amounts of sequence data allows for the
development of novel, data-driven computational approaches. Notably, these
approaches can bridge scales within a single statistical framework. While being
currently applied mostly to isolated problems on single scales, their immense
potential for an evolutionary informed, structural systems biology is steadily
emerging.


  Using coupled biological-physical model based on NPZD-type biological model
and 3D coastal ocean model (ROMS) we studied dissolved oxygen (DO) dynamics and
hypoxia development on Oregon shelf during April-August of 2002, 2006, and
2008. We found that shelf hypoxia existed during summer months of all three
years. It was characterized by variable severity, horizontal and vertical
extent, duration and timing, and it was more pronounced in 2002 and 2006. By
the means of numerical sensitivity analysis we found out that: inadequate
initial DO and NO3 conditions in late-spring 2002 prevented or delayed hypoxia
development; offshore and especially northern DO and NO3 boundary conditions
are important to simulating hypoxia on the Oregon shelf, this was especially
critical for early bottom hypoxia on the shelf north of 45N in 2006; hypoxia
occurred earlier in the north in 2006 and in the south (Heceta Bank) in 2002,
perhaps, due to different northern boundary conditions for these years; the DO
and NO3 conditions at western open boundary located some 400 km offshore are
unimportant for DO dynamics in spring-summer. Although DO production due to
biological processes is large, physical processes, mostly horizontal advection
and diffusion, are responsible for net DO reduction in spring-summer and
hypoxia onset in summer on the Oregon shelf. The physical mechanism most
responsible for Oregon shelf hypoxia is the coastal upwelling. Diffusive fluxes
of NO3 and DO are negligible at northern and southern boundaries of the Oregon
shelf and appreciable at the western boundary. In 2006, about two thirds of
total April-August DO loss happened in April-May as a result of strong and
long-lasting upwelling event.


  Which social decisions are influenced by intuitive processes? Which by
deliberative processes? The dual-process approach to human sociality has
emerged in the last decades as a vibrant and exciting area of research. Yet, a
perspective that integrates empirical and theoretical work is lacking. This
review and meta-analysis synthesizes the existing literature on the cognitive
basis of cooperation, altruism, truth-telling, positive and negative
reciprocity, and deontology, and develops a framework that organizes the
experimental regularities. The meta-analytic results suggest that intuition
favours a set of heuristics that are related to the instinct for
self-preservation: people avoid being harmed, avoid harming others (especially
when there is a risk of harm to themselves), and are averse to disadvantageous
inequalities. Finally, this paper highlights some key research questions to
further advance our understanding of the cognitive foundations of human
sociality.


  We consider a measure of cooperativity based on the minimal absolute
interaction required to generate an observed titration behavior. We describe
the corresponding algebraic optimization problem and show how it can be solved
using the nonlinear algebra tool \texttt{SCIP}. Moreover, we compute the
minimal absolute interactions for various binding polynomials that describe the
oxygen binding of various hemoglobins under different conditions. While
calculated minimal absolute interactions are consistent with the expected
outcome of the chemical modifications, it ranks the cooperativity of the
molecules differently than the maximal Hill slope.


  The strength or weakness of an algorithm is ultimately governed by the
confidence of its result. When the domain of the problem is large (e.g.
traversal of a high-dimensional space), a perfect solution cannot be obtained,
so approximations must be made. These approximations often lead to a reported
quantity of interest (QOI) which varies between runs, decreasing the confidence
of any single run. When the algorithm further computes this final QOI based on
uncertain or noisy data, the variability (or lack of confidence) of the final
QOI increases. Unbounded, these two sources of uncertainty (algorithmic
approximations and uncertainty in input data) can result in a reported
statistic that has low correlation with ground truth.
  In biological applications, this is especially applicable, as the search
space is generally approximated at least to some degree (e.g. a high percentage
of protein structures are invalid or energetically unfavorable) and the
explicit conversion from continuous to discrete space for protein
representation implies some uncertainty in the input data. This research
applies uncertainty quantification techniques to the difficult protein-protein
docking problem, first showing the variability that exists in existing
software, and then providing a method for computing probabilistic certificates
in the form of Chernoff-like bounds. Finally, this paper leverages these
probabilistic certificates to accurately bound the uncertainty in docking from
two docking algorithms, providing a QOI that is both robust and statistically
meaningful.


  It is hard to identify nonlinear biological models strictly from data, with
results that are often sensitive to experimental conditions. Automated
experimental workflows and liquid handling enables unprecedented throughput, as
well as the capacity to generate extremely large datasets. We seek to develop
generalized identifiability conditions for informing the design of automated
experiments to discover predictive nonlinear biological models. For linear
systems, identifiability is characterized by persistence of excitation
conditions. For nonlinear systems, no such persistence of excitation conditions
exist. We use the input-Koopman operator method to model nonlinear systems and
derive identifiability conditions for open-loop systems initialized from a
single initial condition. We show that nonlinear identifiability is
intrinsically tied to the rank of a given dataset's power spectral density,
transformed through the lifted Koopman observable space. We illustrate these
identifiability conditions with a simulated synthetic gene circuit model, the
repressilator. We illustrate how rank degeneracy in datasets results in
overfitted nonlinear models of the repressilator, resulting in poor predictive
accuracy. Our findings provide novel experimental design criteria for discovery
of globally predictive nonlinear models of biological phenomena.


  The problem of determining which nucleotides of an RNA sequence are paired or
unpaired in the secondary structure of an RNA, which we call RNA state
inference, can be studied by different machine learning techniques. Successful
state inference of RNA sequences can be used to generate auxiliary information
for data-directed RNA secondary structure prediction. Bidirectional long
short-term memory (LSTM) neural networks have emerged as a powerful tool that
can model global nonlinear sequence dependencies and have achieved
state-of-the-art performances on many different classification problems. This
paper presents a practical approach to RNA secondary structure inference
centered around a deep learning method for state inference. State predictions
from a deep bidirectional LSTM are used to generate synthetic SHAPE data that
can be incorporated into RNA secondary structure prediction via the Nearest
Neighbor Thermodynamic Model (NNTM). This method produces predicted secondary
structures for a diverse test set of 16S ribosomal RNA that are, on average, 25
percentage points more accurate than undirected MFE structures. These
improvements range from several percentage points for some sequences to nearly
50 percentage points for others. Accuracy is highly dependent on the success of
our state inference method, and investigating the global features of our state
predictions reveals that accuracy of both our state inference and structure
inference methods are highly dependent on the similarity of the sequence to the
dataset. This paper presents a deep learning state inference tool, trained and
tested on 16S ribosomal RNA. Converting these state predictions into synthetic
SHAPE data with which to direct NNTM can result in large improvements in
secondary structure prediction accuracy, as shown on a test set of 16S rRNA.


  Using the continuous-time susceptible-infected-susceptible (SIS) model on
networks, we investigate the problem of inferring the class of the underlying
network when epidemic data is only available at population-level (i.e. the
number of infected individuals at a finite set of discrete times of a single
realisation of the epidemic), the only information likely to be available in
real world settings. To tackle this, epidemics on networks are approximated by
a Birth-and-Death process which keeps track of the number of infected nodes at
population level. The rates of this surrogate model encode both the structure
of the underlying network and disease dynamics. We use extensive simulations
over Regular, Erd\H{o}s-R\'enyi and Barab\'asi-Albert networks to build network
class-specific priors for these rates. % show that different well-known network
classes map onto distinct regions of the parameter space of this model. We then
use Bayesian model selection to recover the most likely underlying network
class, based only on a single realisation of the epidemic. We show that the
proposed methodology yields good results on both synthetic and real-world
networks.


  Sequence set is a widely-used type of data source in a large variety of
fields. A typical example is protein structure prediction, which takes an
multiple sequence alignment (MSA) as input and aims to infer structural
information from it. Almost all of the existing approaches exploit MSAs in an
indirect fashion, i.e., they transform MSAs into position-specific scoring
matrices (PSSM) that represent the distribution of amino acid types at each
column. PSSM could capture column-wise characteristics of MSA, however, the
column-wise characteristics embedded in each individual component sequence were
nearly totally neglected.
  The drawback of PSSM is rooted in the fact that an MSA is essentially an
unordered sequence set rather than a matrix. Specifically, the interchange of
any two sequences will not affect the whole MSA. In contrast, the pixels in an
image essentially form a matrix since any two rows of pixels cannot be
interchanged. Therefore, the traditional deep neural networks designed for
image processing cannot be directly applied on sequence sets. Here, we proposed
a novel deep neural network framework (called Seq-SetNet) for sequence set
processing. By employing a {\it symmetric function} module to integrate
features calculated from preceding layers, Seq-SetNet are immune to the order
of sequences in the input MSA. This advantage enables us to directly and fully
exploit MSAs by considering each component protein individually. We evaluated
Seq-SetNet by using it to extract structural information from MSA for protein
secondary structure prediction. Experimental results on popular benchmark sets
suggests that Seq-SetNet outperforms the state-of-the-art approaches by 3.6% in
precision. These results clearly suggest the advantages of Seq-SetNet in
sequence set processing and it can be readily used in a wide range of fields,
say natural language processing.


  In Nature Microbiology, Palleja and colleagues studied the changes in gut
microbiome composition in twelve healthy men over a period of six months
following an antibiotic intervention. The authors argued that the 'gut
microbiota of the subjects recovered to near-baseline composition within 1.5
months' and only exhibited a 'mild yet long-lasting imprint following
antibiotics exposure.' We here present a series of re-analyses of their
original data which demonstrate a significant loss of microbial taxa even after
the complete study period of 180 days. Additionally we show that the
composition of the microbiomes after the complete study period only moderately
correlates with the initial baseline states. Taken together with the lack of
significant compositional differences between day 42 and day 180, we think that
these findings suggest the convergence of the microbiomes to another stable
composition, which is different from the pre-treatment states, instead of a
recovery of the baseline state. Given the accumulating evidence of the role of
microbiome perturbations in a variety of infectious and non-infectious
diseases, as well as the crucial role antibiotics play in modern medicine, we
consider these differences in compositional states worthy of further
investigation.


  Dynamic properties of spreading infection through a heterogeneous population
are studied numerically and analytically using a dynamic variant of Watts and
Strogatz Small World Network-based stochastic
Susceptible-Exposed-Infectious-Removed (SEIR) epidemic model. This model
includes the main realistic parameters usually characterizing transmissible
diseases, such as the force of infection, latency and infection times. As far
as the latency time remains smaller than that of infection, the proportion of
infected individuals increases exponentially with time, otherwise an
oscillatory behavior appears. This may explain the periodic behaviors in time
observed by the health prevention services. It is also shown that periodic
epidemiological surveys overestimate or underestimate the dynamics of infection
if the survey periods do not exactly correspond to the characteristic times of
the infection. Further discussion is provided on the diffusion and relaxation
processes involved in this model, and their relation to the infection
characteristic time.


  Dynorphin is a neuropeptide involved in pain, addiction and mood regulation.
It exerts its activity by binding to the kappa opioid receptor (KOP) which
belongs to the large family of G-protein coupled receptors. The dynorphin
peptide was discovered in 1975, while its receptor was cloned in 1993. This
review will describe: a) the activities and physiological functions of
dynorphin and its receptor, b) early structure-activity relationship studies
performed before cloning of the receptor (mostly pharmacological and
biophysical studies of peptide analogues), c) structure-activity relationship
studies performed after cloning of the receptor via receptor mutagenesis and
the development of recombinant receptor expression systems, d) structural
biology of the opiate receptors culminating in X-ray structures of the four
opioid receptors in their inactive state and structures of MOP and KOP
receptors in their active state. X-ray and EM structures are combined with NMR
data, which gives complementary insight into receptor and peptide dynamics.
Molecular modelling greatly benefited from the availability of atomic
resolution 3D structures of receptor-ligand complexes and an example of the
strategy used to model a dynorphin-KOP receptor complex using NMR data will be
described. These achievements have led to a better understanding of the complex
dynamics of KOP receptor activation and to the development of new ligands and
drugs.


  Horizontal gene Transfer (HT) denotes the transmission of genetic material
between two living organisms, while the vertical transmission refers to a DNA
transfer from parents to their offspring. Consistent experimental evidence
report that this phenomenon plays an essential role in the evolution of certain
bacterias. In particular, HT is believed to be the main instrument of
developing the antibiotic resistance. In this work, we consider several models
which describe this phenomenon: a stochastic jump process (individual-based)
and the deterministic nonlinear integrod-ifferential equation obtained as a
limit for large populations. We also consider a Hamilton-Jacobi equation,
obtained as a limit of the deterministic model under the assumption of small
mutations. The goal of this paper is to compare these models with the help of
numerical simulations. More specifically, our goal is to understand to which
extent the Hamilton-Jacobi model reproduces the qualitative behavior of the
stochastic model and the phenomenon of evolutionary rescue in particular.


  For $\Lambda$-$n$-coalescents with mutation, we analyse the size $O_n$ of the
partition block of $i\in\{1,\ldots,n\}$ at the time where the first mutation
appears on the tree that affects $i$ and is shared with any other
$j\in\{1,\ldots,n\}$. We provide asymptotics of $O_n$ for $n\to\infty$ and a
recursion for all moments of $O_n$ for finite $n$. This variable gives an upper
bound for the minimal clade size [2], which is not observable in real data. In
applications to genetics, it has been shown to be useful to lower
classification errors in genealogical model selection [10].


  Natural environments can present diverse challenges, but some genotypes
remain fit across many environments. Such `generalists' can be hard to evolve,
out-competed by specialists fitter in any particular environment. Here,
inspired by the search for broadly-neutralising antibodies during B-cell
affinity maturation, we demonstrate that environmental changes on an
intermediate timescale can reliably evolve generalists, even when faster or
slower environmental changes are unable to do so. We find that changing
environments on timescales comparable to evolutionary transients in a
population enhances the rate of evolving generalists from specialists, without
enhancing the reverse process. The yield of generalists is further increased in
more complex dynamic environments, such as a `chirp' of increasing frequency.
Our work offers design principles for how non-equilibrium fitness `seascapes'
can dynamically funnel populations to genotypes unobtainable in static
environments.


  One of the fundamental questions of cultural evolutionary research is how
individual-level processes scale up to generate population-level patterns.
Previous studies in music have revealed that frequency-based bias (e.g.
conformity and novelty) drives large-scale cultural diversity in different ways
across domains and levels of analysis. Music sampling is an ideal research
model for this process because samples are known to be culturally transmitted
between collaborating artists, and sampling events are reliably documented in
online databases. The aim of the current study was to determine whether
frequency-based bias has played a role in the cultural transmission of music
sampling traditions, using a longitudinal dataset of sampling events across
three decades. Firstly, we assessed whether turn-over rates of popular samples
differ from those expected under neutral evolution. Next, we used agent-based
simulations in an approximate Bayesian computation framework to infer what
level of frequency-based bias likely generated the observed data. Despite
anecdotal evidence of novelty bias, we found that sampling patterns at the
population-level are most consistent with conformity bias.


  The human body is able to generate a diverse set of high affinity antibodies,
the soluble form of B cell receptors (BCRs), that bind to and neutralize
invading pathogens. The natural development of BCRs must be understood in order
to design vaccines for highly mutable pathogens such as influenza and HIV. BCR
diversity is induced by naturally occurring combinatorial "V(D)J"
rearrangement, mutation, and selection processes. Most current methods for BCR
sequence analysis focus on separately modeling the above processes. Statistical
phylogenetic methods are often used to model the mutational dynamics of BCR
sequence data, but these techniques do not consider all the complexities
associated with B cell diversification such as the V(D)J rearrangement process.
In particular, standard phylogenetic approaches assume the DNA bases of the
progenitor (or "naive") sequence arise independently and according to the same
distribution, ignoring the complexities of V(D)J rearrangement. In this paper,
we introduce a novel approach to Bayesian phylogenetic inference for BCR
sequences that is based on a phylogenetic hidden Markov model (phylo-HMM). This
technique not only integrates a naive rearrangement model with a phylogenetic
model for BCR sequence evolution but also naturally accounts for uncertainty in
all unobserved variables, including the phylogenetic tree, via posterior
distribution sampling.


  The relationship between division of labor and individuals' spatial behavior
in social insect colonies provides a useful context to study how social
interactions influence the spreading of agent (which could be information or
virus) across distributed agent systems. In social insect colonies, spatial
heterogeneity associated with variations of individual task roles, affects
social contacts, and thus the way in which agent moves through social contact
networks. We used an Agent Based Model (ABM) to mimic three realistic scenarios
of agent spreading in social insect colonies. Our model suggests that
individuals within a specific task interact more with consequences that agent
could potentially spread rapidly within that group, while agent spreads slower
between task groups. Our simulations show a strong linear relationship between
the degree of spatial heterogeneity and social contact rates, and that the
spreading dynamics of agents follow a modified nonlinear logistic growth model
with varied transmission rates for different scenarios. Our work provides an
important insights on the dual-functionality of physical contacts. This
dual-functionality is often driven via variations of individual spatial
behavior, and can have both inhibiting and facilitating effects on agent
transmission rates depending on environment. The results from our proposed
model not only provide important insights on mechanisms that generate spatial
heterogeneity, but also deepen our understanding of how social insect colonies
balance the benefit and cost of physical contacts on the agents' transmission
under varied environmental conditions.


  The total phenols and flavonoids content of different polar solvent extracts
from Aurea Helianthus flowers, and their antioxidant activity were determined.
The ethanol extract of Aurea Helianthus flowers were suspended in water and
fractionated using different polar solvents; hexane, chloroform, ethyl acetate,
butanol and water. The parameters of each extract mentioned above were
determined using Floin-ciocalteu reagent(FCR) method, AlCl3 colorimetry method,
ferric reducing ability of plasma(FRAP) assay, total antioxidant activity(TAA)
assay and DPPH radical scavenging assay. The highest total phenols
content(516.21 mg GAE/g) and flavonoids content(326.06 mg QCE/g) were obtained
in ethyl acetate extract, the correlation between TPC and TFC assay was founded
to be 0.967. All polar solvent extracts of Aurea Helianthus flowers showed
significant antioxidant effects, the hightest inhibition was obtained in ethyl
acetate and choroform extracts and the lowest inhibition in the water extract.
There is a good correlation of total phenols and flavonoids content with
antioxidant activity. This work indicated that the polar solvent extracts of
Aurea Helianthus flowers contain high phenols and flavonoids content and
exhibited antioxidant activities in vitro, therefore, could be candidates for
use as natural antioxidant.


  The effects of purification methods and its hypolipidemic function on the
total flavonoids of Aurea Helianthus flower were investigated. Liquid-liquid
extraction of ethanol extract from Aurea Helianthus flower was carried out by
using different polar solvents. The extract with the highest total flavonoid
content was selected, and the optimal conditions for purification of total
flavonoids were determined by purification with macroporous resin. The human
digestive environment was simulated in vitro, and the binding ability of
different flavonoid samples to three kinds of cholate was compared. The results
showed that the purity of total flavonoids in ethanol extract was 27.8%, the
purity of total flavonoids in ethyl acetate extract was 46.4%, and the purity
was increased by 18.6%. Subsequent purification with AB-8 macroporous resin;
loading of total flavonoids at a concentration of 5.5 mg/mL, flow rate of 1.5
mL/min, 110 mL; use of 75% ethanol, 80 mL as eluent at a flow rate of 1.5 mL
The elution at /min resulted in a total flavonoid purity of 83.5 % and an
increase of 37.1%, and a good purification effect was obtained. The binding
rate of total flavonoids purified by AB-8 macroporous resin to sodium
taurocholate, sodium glycocholate and sodium cholate was 88.2%, 73.2% and 75.8
%, respectively. The binding ability was the strongest, and the others were
ethyl acetate. Extract, ethanol extract. The purity of total flavonoids showed
a good correlation with the binding capacity of cholate, and the correlation
coefficient was between 0.963 and 0.988. The total flavonoids of Aurea
Helianthus flower have good bile acid binding ability and can be used as the
focus of natural hypolipidemic substances.


  We use a two-sex partial differential equation (PDE) model based on the
Grandmother Hypothesis. We build on an earlier model by Kim et al. [27] by
allowing for evolution in both longevity and age at last birth, and also
assuming that post-fertile females support only their daughters' fertility.
Similarly to Kim et al. [27], we find that only two locally stable equilibria
exist: one corresponding to great ape-like longevities and the other
corresponding to hunter-gatherer longevities. Our results show that
grandmothering enables the transition between these two equilibria, without
extending the end of fertility. Moreover, sensitivity analyses of the model
show that male competition, arising from a skew in the mating sex ratio towards
males, plays a significant role in determining whether the transition from
great ape-like longevities to higher longevities is possible and the
equilibrium value of the average adult lifespan. Whereas grandmothering effects
have a significant impact on the equilibrium value of the average age at last
birth and enable the transition to higher longevities, they have an
insignificant impact on the equilibrium value of the average adult lifespan.


  Great apes, the other living members of our hominid family, become decrepit
before the age of forty and rarely outlive their fertile years. In contrast,
women - even in high mortality hunter-gatherer populations - usually remain
healthy and productive well beyond menopause. The grandmother hypothesis aims
to account for the evolution of this distinctive feature of human life history.
Our previous mathematical simulations of that hypothesis fixed the end of
female fertility at the age of 45, based on the similarities among living
hominids, and then modeled the evolution of human-like longevity from an
ancestral state, like that of the great apes, due only to grandmother effects.
A major modification here allows the age female fertility ends to vary as well,
directly addressing a version of the question, influentially posed by GC
Williams six decades ago: Why isn't menopause later in humans? Our model is an
agent-based model (ABM) that accounts for the coevolution of both expected
adult lifespan and end of female fertility as selection maximizes reproductive
value. We find that grandmother effects not only drive the population from an
equilibrium representing a great ape-like longevity to a new human-like
longevity, they also maintain the observed termination of women's fertility
before the age of 50.


  Galled trees are studied as a recombination model in theoretic population
genetics. This class of phylogenetic networks has been generalized to
tree-child networks, normal networks and tree-based networks by relaxing a
structural condition. Although these networks are simple, their topological
structures have yet to be fully understood. It is well-known that all
phylogenetic trees on $n$ taxa can be generated by the insertion of the $n$-th
taxa to each edge of all the phylogenetic trees on $n-1$ taxa. We prove that
all tree-child networks with $k$ reticulate nodes on $n$ taxa can be uniquely
generated via three operations from all the tree-child networks with $k-1$ or
$k$ reticulate nodes on $n-1$ taxa . An application of this result is found in
counting tree-child networks and normal networks. In particular, a simple
formula is given for the number of rooted phylogenetic networks with one
reticulate node.


  Modern biological techniques such as Hi-C permit to measure probabilities
that different chromosomal regions are close in space. These probabilities can
be visualised as matrices called contact maps. In this paper, we introduce a
multifractal analysis of chromosomal contact maps. Our analysis reveals that
Hi-C maps are bifractal, i.e. complex geometrical objects characterized by two
distinct fractal dimensions. To rationalize this observation, we introduce a
model that describes chromosomes as a hierarchical set of nested domains and we
solve it exactly. The predicted multifractal spectrum is in excellent
quantitative agreement with experimental data. Moreover, we show that our
theory yields to a more robust estimation of the scaling exponent of the
contact probability than existing methods. By applying this method to
experimental data, we detect subtle conformational changes among chromosomes
during differentiation of human stem cells.


  In this work we addressed the problem of stability analysis for an uncertain
piecewise affine model of a genetic regulatory network. In particular we
considered polytopic parameter uncertainties on the proteins production rate
functions, giving conditions for the existence of a piecewise quadratic
Lyapunov function for any possible realization of the system. In the spirit of
other works in literature, the resulting conditions will be given on the
vertices of the parameter polytope, while still taking into consideration the
piecewise nature of the Lyapunov function and the presence, in general, of
sliding modes solutions. An example is shown to prove the validity and
applicability of the theoretical results.


  In the past years, many computational methods have been developed to infer
the structure of gene regulatory networks from time-series data. However, the
applicability and accuracy presumptions of such algorithms remain unclear due
to experimental heterogeneity. This paper assesses the performance of recent
and successful network inference strategies under a novel, multifactorial
evaluation framework in order to highlight pragmatic tradeoffs in experimental
design. The effects of data quantity and systems perturbations are addressed,
thereby formulating guidelines for efficient resource management.
  Realistic data were generated from six widely used benchmark models of
rhythmic and non-rhythmic gene regulatory systems with random perturbations
mimicking the effect of gene knock-out or chemical treatments. Then,
time-series data of increasing lengths were provided to five state-of-the-art
network inference algorithms representing distinctive mathematical paradigms.
The performances of such network reconstruction methodologies are uncovered
under various experimental conditions. We report that the algorithms do not
benefit equally from data increments. Furthermore, for rhythmic systems, it is
more profitable for network inference strategies to be run on long time-series
rather than short time-series with multiple perturbations. By contrast, for the
non-rhythmic systems, increasing the number of perturbation experiments yielded
better results than increasing the sampling frequency. We expect that future
benchmark and algorithm design would integrate such multifactorial
considerations to promote their widespread and conscientious usage.


  The use of model experimental systems and mathematical models is important to
further understanding of infectious disease dynamics and strategize disease
mitigation. Gyrodactylids are helminth ectoparasites of teleost fish which have
many dynamical characteristics of microparasites but offer the advantage that
they can be quantified and tracked over time, allowing further insight into
within-host and epidemic dynamics. In this paper, we design a model to describe
host-parasite dynamics of the well-studied guppy-Gyrodactylus turnbulli system,
using experimental data to estimate parameters and validate it. We estimate the
basic reproduction number (R_0), for this system. Sensitivity analysis reveals
that parasite growth rate, and the rate at which the guppy mounts an immune
response have the greatest impact on outbreak peak and timing both for initial
outbreaks and on longer time scales. These findings highlight guppy population
resistance and parasite virulence as key factors in disease control, and future
work should focus on incorporating heterogeneity in host resistance into
disease models and extrapolating to other host-parasite systems.


  The current state of cancer therapeutics has been moving away from
one-size-fits-all cytotoxic chemotherapy, and towards a more individualized and
specific approach involving the targeting of each tumor's genetic
vulnerabilities. Different tumors, even of the same type, may be more reliant
on certain cellular pathways more than others. With modern advancements in our
understanding of cancer genome sequencing, these pathways can be discovered.
Investigating each of the millions of possible small molecule inhibitors for
each kinase in vitro, however, would be extremely expensive and time consuming.
This project focuses on predicting the inhibition activity of small molecules
targeting 8 different kinases using multiple deep learning models. We trained
fingerprint-based MLPs and simplified molecular-input line-entry specification
(SMILES)-based recurrent neural networks (RNNs) and molecular graph
convolutional networks (GCNs) to accurately predict inhibitory activity
targeting these 8 kinases.


  Population boundary is a classic indicator of climatic response in ecology.
In addition to known challenges, the spatial and dynamical characteristics of
the boundary are not only affected by the spatial gradient in the environmental
factors, but also by local heterogeneities in the regional characteristics.
Here, we capture the effects of quenched heterogeneities on the ecological
boundary with the disordered contact process in one and two dimensions with a
linear spatial trend in the local control parameter. We apply the
strong-disorder renormalization group method to calculate the sites occupied
with an $O(1)$ probability in the stationary state, readily yielding the
population front's position as the outermost site locally as well as globally
for the entire boundary. We show that under a quasistatic change of the global
environment, mimicking climate change, the front advances intermittently: long
quiescent periods are interrupted by rare but long jumps. The characteristics
of this intermittent dynamics are found to obey universal scaling laws in terms
of the gradient, conjectured to be related to the correlation-length exponent
of the model. Our results suggest that current observations might misleadingly
show little to no climate response for an extended period of time, concealing
the long-term effects of climate change.


  Machine learning is often used in virtual screening to find compounds that
are pharmacologically active on a target protein. The weave module is a type of
graph convolutional deep neural network that uses not only features focusing on
atoms alone (atom features) but also features focusing on atom pairs (pair
features); thus, it can consider information of nonadjacent atoms. However, the
correlation between the distance on the graph and the three-dimensional
coordinate distance is uncertain. In this paper, we propose three improvements
for modifying the weave module. First, the distances between ring atoms on the
graph were modified to bring the distances on the graph closer to the
coordinate distance. Second, different weight matrices were used depending on
the distance on the graph in the convolution layers of the pair features.
Finally, a weighted sum, by distance, was used when converting pair features to
atom features. The experimental results show that the performance of the
proposed method is slightly better than that of the weave module, and the
improvement in the distance representation might be useful for compound
activity prediction.


  Local search is widely used to solve combinatorial optimisation problems and
to model biological evolution, but the performance of local search algorithms
on different kinds of fitness landscapes is poorly understood. Here we consider
how fitness landscapes can be represented using valued constraints, and
investigate what the structure of such representations reveals about the
complexity of local search.
  First, we show that for fitness landscapes representable by binary Boolean
valued constraints there is a minimal necessary constraint graph that can be
easily computed. Second, we consider landscapes as equivalent if they allow the
same (improving) local search moves; we show that a minimal constraint graph
still exists, but is NP-hard to compute.
  We then develop several techniques to bound the length of any sequence of
local search moves. We show that such a bound can be obtained from the
numerical values of the constraints in the representation, and show how this
bound may be tightened by considering equivalent representations. In the binary
Boolean case, we prove that a degree 2 or tree-structured constraint graph
gives a quadratic bound on the number of improving moves made by any local
search; hence, any landscape that can be represented by such a model will be
tractable for any form of local search.
  Finally, we build two families of examples to show that the conditions in our
tractability results are essential. With domain size three, even just a path of
binary constraints can model a landscape with an exponentially long sequence of
improving moves. With a treewidth-two constraint graph, even with a maximum
degree of three, binary Boolean constraints can model a landscape with an
exponentially long sequence of improving moves.


  Microbial communities routinely have several possible species compositions or
community states observed for the same environmental parameters. Changes in
these parameters can trigger abrupt and persistent transitions (regime shifts)
between such community states. Yet little is known about the main determinants
and mechanisms of multistability in microbial communities. Here we introduce
and study a resource-explicit model in which microbes compete for two types of
essential nutrients. We adapt game-theoretical methods of the stable matching
problem to identify all possible species compositions of a microbial community.
We then classify them by their resilience against three types of perturbations:
fluctuations in nutrient supply, invasions by new species, and small changes of
abundances of existing ones. We observe multistability and explore an intricate
network of regime shifts between stable states in our model. Our results
suggest that multistability requires microbial species to have different
stoichiometries of essential nutrients. We also find that balanced nutrient
supply promote multistability and species diversity yet make individual
community states less stable.


  The strong bending of polymers is poorly understood. We propose a general
quantitative framework of polymer bending that includes both the weak and
strong bending regimes on the same footing, based on a single general physical
principle. As the bending deformation increases beyond a certain
(polymer-specific) point, the change in the convexity properties of the
effective bending energy of the polymer makes the harmonic deformation
energetically unfavorable: in this strong bending regime the energy of the
polymer varies linearly with the average bending angle as the system follows
the convex hull of the deformation energy function. For double-stranded DNA,
the effective bending deformation energy becomes non-convex for bends greater
than about 2 degrees per base-pair, equivalent to the curvature of a closed
circular loop of about 160 base pairs. A simple equation is derived for the
polymer loop energy that covers both the weak and strong bending regimes. The
theory shows quantitative agreement with recent DNA cyclization experiments on
short DNA fragments, while maintaining the expected agreement with experiment
in the weak bending regime. Counter-intuitively, cyclization probability
(j-factor) of very short DNA loops is predicted to increase with decreasing
loop length; the j-factor reaches its minimum for loops of 45 base pairs.
Atomistic simulations reveal that the attractive component of the short-range
Lennard-Jones interaction between the backbone atoms can explain the underlying
non-convexity of the DNA effective bending energy, leading to the linear
bending regime. Applicability of the theory to protein-DNA complexes, including
the nucleosome, is discussed.


  We present a mathematical simplification for the evolutionary dynamics of a
heritable trait within a two-sex population. This trait is assumed to control
the timing of sex-specific life-history events, such as the age of sexual
maturity and end of female fertility, and each sex has a distinct fitness
tradeoff associated with the trait. We provide a formula for the fitness
landscape of the population and show a natural extension of the result to an
arbitrary number of heritable traits. Our method can be viewed as a dynamical
systems generalisation of the Price Equation to include two sexes{, age
structure} and multiple traits. We use this formula to examine the effect of
grandmothering, whereby post-fertile females subsidize their daughter's
fertility by provisioning grandchildren. Grandmothering can drive a shift
towards higher sex ratios due to lengthening female post-fertile longevity,
leading to changes in fitness for both sexes. For males, increased longevity is
accompanied by a substantially longer fertile lifespan resulting in higher sex
ratios in the fertile ages. Our fitness landscapes show a net increase in
fitness for both males and females at longer lifespans, and as a result, we
find that grandmothering alone provides an evolutionary trajectory to higher
longevities.


  Transposable Elements (TEs) or jumping genes are the DNA sequences that have
an intrinsic capability to move within a host genome from one genomic location
to another. Studies show that the presence of a TE within or adjacent to a
functional gene may alter its expression. TEs can also cause an increase in the
rate of mutation and can even mediate duplications and large insertions and
deletions in the genome, promoting gross genetic rearrangements. Thus, the
proper classification of the identified jumping genes is essential to
understand their genetic and evolutionary effects in the genome. While
computational methods have been developed that perform either binary
classification or multi-label classification of TEs, few studies have focused
on their hierarchical classification. The state-of-the-art machine learning
classification method utilizes a Multi-Layer Perceptron (MLP), a class of
neural network, for hierarchical classification of TEs. However, the existing
methods have limited accuracy in classifying TEs. A more effective classifier,
which can explain the role of TEs in germline and somatic evolution, is needed.
In this study, we examine the performance of a variety of machine learning (ML)
methods. And eventually, propose a robust approach for the hierarchical
classification of TEs, with higher accuracy, using Support Vector Machines
(SVM).


  We consider three classical models of biological evolution: (i) the Moran
process, an example of a reducible Markov Chain; (ii) the Kimura Equation, a
particular case of a degenerated Fokker-Planck Diffusion; (iii) the Replicator
Equation, a paradigm in Evolutionary Game Theory. While these approaches are
not completely equivalent, they are intimately connected, since (ii) is the
diffusion approximation of (i), and (iii) is obtained from (ii) in an
appropriate limit. It is well known that the Replicator Dynamics for two
strategies is a gradient flow with respect to the celebrated Shahshahani
distance. We reformulate the Moran process and the Kimura Equation as gradient
flows and in the sequel we discuss conditions such that the associated gradient
structures converge: (i) to (ii) and (ii) to (iii). This provides a geometric
characterisation of these evolutionary processes and provides a reformulation
of the above examples as time minimization of free energy functionals.


  Human social hierarchy has the unique characteristic of existing in two
forms. Firstly, as an informal hierarchy where leaders and followers are
implicitly defined by their personal characteristics, and secondly, as an
institutional hierarchy where leaders and followers are explicitly appointed by
group decision. Although both forms can reduce the time spent in organising
collective tasks, institutional hierarchy imposes additional costs. It is
therefore natural to question why it emerges at all. The key difference lies in
the fact that institutions can create hierarchy with only a single leader,
which is unlikely to occur in unregulated informal hierarchy. To investigate if
this difference can affect group decision-making and explain the evolution of
institutional hierarchy, we first build an opinion-formation model that
simulates group decision making. We show that in comparison to informal
hierarchy, a single-leader hierarchy reduces (i) the time a group spends to
reach consensus, (ii) the variation in consensus time, and (iii) the rate of
increase in consensus time as group size increases. We then use this model to
simulate the cost of organising a collective action which produces resources,
and integrate this into an evolutionary model where individuals can choose
between informal or institutional hierarchy. Our results demonstrate that
groups evolve preferences towards institutional hierarchy, despite the cost of
creating an institution, as it provides a greater organisational advantage
which is less affected by group size and inequality.


  Smooth cordgrass Spartina alterniflora is a grass species commonly found in
tidal marshes. It is an ecosystem engineer, capable of modifying the structure
of its surrounding environment through various feedbacks. The scale-dependent
feedback between marsh grass and sediment volume is particularly of interest.
Locally, the marsh vegetation attenuates hydrodynamic energy, enhancing
sediment accretion and promoting further vegetation growth. In turn, the
diverted water flow promotes the formation of erosion troughs over longer
distances. This scale-dependent feedback may explain the characteristic
spatially varying marsh shoreline, commonly observed in nature. We propose a
mathematical framework to model grass-sediment dynamics as a system of
reaction-diffusion equations with an additional nonlocal term quantifying the
short-range positive and long-range negative grass-sediment interactions. We
use a Mexican-hat kernel function to model this scale-dependent feedback. We
perform a steady state biharmonic approximation of our system and derive
conditions for the emergence of spatial patterns, corresponding to a spatially
varying marsh shoreline. We find that the emergence of such patterns depends on
the spatial scale and strength of the scale-dependent feedback, specified by
the width and amplitude of the Mexican-hat kernel function.


  Cooperators forgo their interest to benefit others. Thus cooperation should
not be favored by natural selection. It challenges the evolutionists, since
cooperation is widespread. As one of the resolutions, information spreading has
been revealed to play a key role in the emergence of cooperation. Individuals,
however, are typically assumed to be passive in the information spreading. Here
we assume that individuals are active to spread the information via
self-recommendation. Individuals with higher intensities of self-recommendation
are likely to have more neighbors. We find that i) eloquent cooperators are
necessary to promote cooperation; ii) individuals need to be open to the
self-recommendation to enhance cooperation level; iii) the cost-to-benefit
ratio should be smaller than one minus the ratio between self-recommendation
intensities of defector and cooperator, which qualitatively measures the
viscosity of the population. Our results highlight the importance of active
information spreading on cooperation.


  Aim: In the present work, we aimed to evaluate a comprehensive radiomics
framework that enabled prediction of EGFR and KRAS mutation status in NSCLC
cancer patients based on PET and CT multi-modalities radiomic features and
machine learning (ML) algorithms. Methods: Our study involved 211 NSCLC cancer
patient with PET and CTD images. More than twenty thousand radiomic features
from different image-feature sets were extracted Feature value was normalized
to obtain Z-scores, followed by student t-test students for comparison, high
correlated features were eliminated and the False discovery rate (FDR)
correction were performed Six feature selection methods and twelve classifiers
were used to predict gene status in patient and model evaluation was reported
on independent validation sets (68 patients). Results: The best predictive
power of conventional PET parameters was achieved by SUVpeak (AUC: 0.69,
P-value = 0.0002) and MTV (AUC: 0.55, P-value = 0.0011) for EGFR and KRAS,
respectively. Univariate analysis of radiomics features improved prediction
power up to AUC: 75 (q-value: 0.003, Short Run Emphasis feature of GLRLM from
LOG preprocessed image of PET with sigma value 1.5) and AUC: 0.71 (q-value
0.00005, The Large Dependence Low Gray Level Emphasis from GLDM in LOG
preprocessed image of CTD sigma value 5) for EGFR and KRAS, respectively.
Furthermore, the machine learning algorithm improved the perdition power up to
AUC: 0.82 for EGFR (LOG preprocessed of PET image set with sigma 3 with VT
feature selector and SGD classifier) and AUC: 0.83 for KRAS (CT image set with
sigma 3.5 with SM feature selector and SGD classifier). Conclusion: We
demonstrated that radiomic features extracted from different image-feature sets
could be used for EGFR and KRAS mutation status prediction in NSCLC patients,
and showed that they have more predictive power than conventional imaging
parameters.


  Animals use a wide variety of strategies to reduce or avoid aggression in
conflicts over resources. These strategies range from sharing resources without
outward signs of conflict to the development of dominance hierarchies, in which
initial fighting is followed by the submission of subordinates. Although models
have been developed to analyze specific strategies for resolving conflicts over
resources, little work has focused on trying to understand why particular
strategies are more likely to arise in certain situations. In this paper, we
use a model based on an iterated Hawk--Dove game to analyze how resource
holding potentials (RHPs) and other factors affect whether sharing, dominance
relationships, or other behaviours are evolutionarily stable. We find through
extensive numerical simulations that sharing is stable only when the cost of
fighting is low and the animals in a contest have similar RHPs, whereas
dominance relationships are stable in most other situations. We also explore
what happens when animals are unable to assess each other's RHPs without
fighting, and we compare a range of strategies for this problem using
simulations. We find (1) that the most successful strategies involve a limited
period of assessment followed by a stable relationship in which fights are
avoided and (2) that the duration of assessment depends both on the costliness
of fighting and on the difference between the animals' RHPs. Along with our
direct work on modeling and simulations, we develop extensive software to
facilitate further testing; it is available at
\url{https://bitbucket.org/CameronLHall/dominancesharingassessmentmatlab/}.


  Current methods for investigation of receptor - ligand interactions in drug
discovery are based on three-dimensional complementarity of receptor and ligand
surfaces, and they include pharmacophore modelling, QSAR, molecular docking
etc. Those methods only consider short-range molecular interactions (distances
<5A), and not include long-range interactions (distances >5A) which are
essential for kinetic of biochemical reactions because they influence the
number of productive collisions between interacting molecules. Previously was
shown that the electron-ion interaction potential (EIIP) represents the
physical property which determines the long-range properties of biological
molecules. This molecular descriptor served as a base for development of the
informational spectrum method (ISM), a virtual spectroscopy method for
investigation of protein-protein interactions. In this paper, we proposed a new
approach to treat small molecules as linear entities, allowing study of the
small molecule - protein interaction by ISM. We analyzed here 21 sets of KEGG
drug-protein interactions and showed that this new approach allows an efficient
discrimination between biologically active and inactive ligands, and
consistence with AA regions of their binding site on the target protein.


  Cellular functions are established through biological evolution, but are
constrained by the laws of physics. For instance, the physics of protein
folding limits the lengths of cellular polypeptide chains. Consequently, many
cellular functions are carried out not by long, isolated proteins, but rather
by multi-protein complexes. Protein complexes themselves do not escape physical
constraints, one of the most important being the difficulty to assemble
reliably in the presence of cellular noise. In order to lay the foundation for
a theory of reliable protein complex assembly, we study here an equilibrium
thermodynamic model of self-assembly that exhibits four distinct assembly
behaviors: diluted protein solution, liquid mixture, "chimeric assembly" and
"multifarious assembly". In the latter regime, different protein complexes can
coexist without forming erroneous chimeric structures. We show that two
conditions have to be fulfilled to attain this regime: (i) the composition of
the complexes needs to be sufficiently heterogeneous, and (ii) the use of the
set of components by the complexes has to be sparse. Our analysis of publicly
available databases of protein complexes indicates that cellular protein
systems might have indeed evolved so to satisfy both of these conditions.


  We train a neural network to predict human gene expression levels based on
experimental data for rat cells. The network is trained with paired human/rat
samples from the Open TG-GATES database, where paired samples were treated with
the same compound at the same dose. When evaluated on a test set of held out
compounds, the network successfully predicts human expression levels. On the
majority of the test compounds, the list of differentially expressed genes
determined from predicted expression levels agrees well with the list of
differentially expressed genes determined from actual human experimental data.


  To effectively forage in natural environments, organisms must adapt to
changes in the quality and yield of food sources across multiple timescales.
Individuals foraging in groups act based on both their private observations and
the opinions of their neighbors. How do these information sources interact in
changing environments? We address this problem in the context of honeybee
swarms, showing inhibitory social interactions help maintain adaptivity and
consensus needed for effective foraging. Individual and social interactions of
a mathematical swarm model shape the nutrition yield of a group foraging from
feeders with temporally switching food quality. Social interactions improve
foraging from a single feeder if temporal switching is fast or feeder quality
is low. When the swarm chooses from multiple feeders, the most effective form
of social interaction is direct switching, whereby bees flip the opinion of
nestmates foraging at lower yielding feeders. Model linearization shows that
effective social interactions increase the fraction of the swarm at the correct
feeder (consensus) and the rate at which bees reach that feeder (adaptivity).
Our mathematical framework allows us to compare a suite of social inhibition
mechanisms, suggesting experimental protocols for revealing effective swarm
foraging strategies in dynamic environments.


  Protein function and activity are a consequence of its three-dimensional
structure. Single particle analysis of cryogenic electron micrographs has
radically changed structural biology allowing atomic reconstruction of almost
any type of proteins. While such an approach provides snapshots of
three-dimensional structural information that can be correlated with function,
the new frontier of protein structural biology is in the fourth dimension,
time. Here we propose the use of liquid phase electron microscopy to expand
structural biology into dynamic studies. We apply here single particle analysis
algorithm to images of proteins in Brownian motion through time; thus, Brownian
single particle analysis (BSPA). BSPA enables to reduce the acquisition time
from hours, in cryo-EM, to seconds and achieve information on conformational
changes, hydration dynamics, and effects of thermal fluctuations. Yielding all
these previously neglected aspects, BSPA may lead to the verge of a new field:
dynamic structural biology.


  Most amino acids are encoded by multiple synonymous codons. For an amino
acid, some of its synonymous codons are used much more rarely than others.
Analyses of positions of such rare codons in protein sequences revealed that
rare codons can impact co-translational protein folding and that positions of
some rare codons are evolutionary conserved. Analyses of positions of rare
codons in proteins' 3-dimensional structures, which are richer in biochemical
information than sequences alone, might further explain the role of rare codons
in protein folding. We analyze a protein set recently annotated with codon
usage information, considering non-redundant proteins with sufficient
structural information. We model the proteins' structures as networks and study
potential differences between network positions of amino acids encoded by
evolutionary conserved rare, evolutionary non-conserved rare, and commonly used
codons. In 84% of the proteins, at least one of the three codon categories
occupies significantly more or less network-central positions than the other
codon categories. Different protein groups showing different codon centrality
trends (i.e., different types of relationships between network positions of the
three codon categories) are enriched in different biological functions,
implying the existence of a link between codon usage, protein folding, and
protein function.


  Understanding age-group dynamics of infectious diseases is a fundamental
issue for both scientific study and policymaking. Age-structure epidemic models
were developed in order to study and improve our understanding of these
dynamics. By fitting the models to incidence data of real outbreaks one can
infer estimates of key epidemiological parameters. However, estimation of the
transmission in an age-structured populations requires first to define the
age-groups of interest. Misspecification in representing the heterogeneity in
the age-dependent transmission rates can potentially lead to biased estimation
of parameters. We develop the first statistical, data-driven methodology for
deciding on the best partition of incidence data into age-groups. The method
employs a top-down hierarchical partitioning algorithm, with a metric distance
built for maximizing mathematical identifiability of the transmission matrix,
and a stopping criteria based on significance testing. The methodology is
tested using simulations showing good statistical properties. The methodology
is then applied to influenza incidence data of 14 seasons in order to extract
the significant age-group clusters in each season.


  Background and aim: This study aimed to predict methylation status of the O-6
methyl guanine-DNA methyl transferase (MGMT) gene promoter status by using MRI
radiomics features, as well as univariate and multivariate analysis.
  Material and Methods: Eighty-two patients who had a MGMT methylation status
were include in this study. Tumor were manually segmented in the four regions
of MR images, a) whole tumor, b) active/enhanced region, c) necrotic regions
and d) edema regions (E). About seven thousand radiomics features were
extracted for each patient. Feature selection and classifier were used to
predict MGMT status through different machine learning algorithms. The area
under the curve (AUC) of receiver operating characteristic (ROC) curve was used
for model evaluations.
  Results: Regarding univariate analysis, the Inverse Variance feature from
gray level co-occurrence matrix (GLCM) in Whole Tumor segment with 4.5 mm Sigma
of Laplacian of Gaussian filter with AUC: 0.71 (p-value: 0.002) was found to be
the best predictor. For multivariate analysis, the decision tree classifier
with Select from Model feature selector and LOG filter in Edema region had the
highest performance (AUC: 0.78), followed by Ada Boost classifier with Select
from Model feature selector and LOG filter in Edema region (AUC: 0.74).
  Conclusion: This study showed that radiomics using machine learning
algorithms is a feasible, noninvasive approach to predict MGMT methylation
status in GBM cancer patients
  Keywords: Radiomics, Radiogenomics, GBM, MRI, MGMT


  Various theoretical studies have been carried out to infer relevant
protein-protein interactions among pathogens and their hosts. Such studies are
generally based on preferential attachment of bacteria / virus to their human
receptor homologs. We have analyzed 17 pathogenic species mainly belonging to
bacterial and viral pathogenic classes, with the aim to identify the
interacting human proteins which are targeted by both bacteria and virus
specifically. Our study reveals that the TLRs play a crucial role between the
pathogen-associated molecular patterns (PAMPs) and the damage associated
molecular patterns (DAMPS). PAMPs include bacterial lipopolysaccharides (LPs),
endotoxins, bacterial flagellin, lipoteichoic acid, peptidoglycan in bacterial
organisms and nucleic acid variants associated with viral organisms, whereas
DAMPs are associated with host biomolecules that perpetuate non-infectious
inflammatory responses. We found out the presence of SOD1 and SOD2 (superoxide
dismutase) is crucial, as it acts as an anti-oxidant and helps in eliminating
oxidative stress by preventing damage from reactive oxygen species. Hence, such
strategies can be used as new therapies for anti-inflammatory diseases with
significant clinical outcomes.


  The classification of amino acids and their sequence analysis plays a vital
role in life sciences and is a challenging task. This article uses and compares
state-of-the-art deep learning models like convolution neural networks (CNN),
long short-term memory (LSTM), and gated recurrent units (GRU) to solve
macromolecule classification problems using amino acids. These models have
efficient frameworks for solving a broad spectrum of complex learning problems
compared to traditional machine learning techniques. We use word embedding to
represent the amino acid sequences as vectors. The CNN extracts features from
amino acid sequences, which are treated as vectors, then fed to the models
mentioned above to train a robust classifier. Our results show that word2vec as
embedding combined with VGG-16 performs better than LSTM and GRU. The proposed
approach gets an error rate of 1.5%.


  The Trojan Y Chromosome (TYC) Strategy is a promising eradication method that
attempts to manipulate the female to male ratio to promote the reduction of the
population of an invasive species. The manipulation stems from an introduction
of sex-reversed males, called supermales, into an ecosystem. The offspring of
the supermales is guaranteed to be male. Mathematical models have shown that
the population can be driven to extinction with a continuous supply of
supermales. In this paper, a new model of the TYC strategy is introduced and
analyzed that includes two important modeling characteristics, that are
neglected in all previous models. First, the new model includes intraspecies
competition for mates. Second, a strong Allee effect is included. Several
conclusions about the strategy via optimal control are established. These
results have large scale implications for the biological control of invasive
species.


  The adaptive immune system is a dynamical, self-organized multiscale system
that protects vertebrates from both pathogens and internal irregularities, such
as tumours. For these reason it fascinates physicists, yet the multitude of
different cells, molecules and sub-systems is often also petrifying. Despite
this complexity, as experiments on different scales of the adaptive immune
system become more quantitative, many physicists have made both theoretical and
experimental contributions that help predict the behaviour of ensembles of
cells and molecules that participate in an immune response. Here we review some
recent contributions with an emphasis on quantitative questions and
methodologies. We also provide a more general methods section that presents
some of the wide array of theoretical tools used in the field.


  Surface segregation of the low-molecular weight component in a polymeric
mixture leads to degradation of industrial formulations. We report a
simultaneous phase separation and surface migration phenomena in
oligomer-polymer and oligomer-gel systems following a temperature quench. We
compute equilibrium and time varying migrant density profiles and wetting layer
thickness using coarse grained molecular dynamics and mesoscale hydrodynamics
simulations to demonstrate that surface migration in oligomer-gel systems is
significantly reduced due to network elasticity. Further, phase separation
processes are significantly slowed in gels, modifying the
Lifshitz-Slyozov-Wagner (LSW) law $\ell(\tau) \sim \tau^{1/3}$. Our work allows
for rational design of polymer/gel-oligomer mixtures with predictable surface
segregation characteristics.


  The output of molecular dynamics simulations is high-dimensional, and the
degrees of freedom among the atoms are related in intricate ways. Therefore, a
variety of analysis frameworks have been introduced in order to distill complex
motions into lower-dimensional representations that model the system dynamics.
These dynamical models have been developed to optimally approximate the
system's global kinetics. However, the separate aims of optimizing global
kinetics and modeling a process of interest diverge when the process of
interest is not the slowest process in the system. Here, we introduce deflation
into state-of-the-art methods in molecular kinetics in order to preserve the
use of variational optimization tools when the slowest dynamical mode is not
the same as the one we seek to model and understand. First, we showcase
deflation for a simple toy system and introduce the deflated variational
approach to Markov processes (dVAMP). Using dVAMP, we show that nondominant
reaction coordinates produced using deflation are more informative than their
counterparts generated without deflation. Then, we examine a protein folding
system in which the slowest dynamical mode is not folding. Following a dVAMP
analysis, we show that deflation can be used to obscure this undesired slow
process from a kinetic model, in this case a VAMPnet. The incorporation of
deflation into current methods opens the door for enhanced sampling strategies
and more flexible, targeted model building.


  When studying multi-body protein complexes, biochemists use computational
tools that can suggest hundreds or thousands of their possible spatial
configurations. However, it is not feasible to experimentally verify more than
only a very small subset of them. In this paper, we propose a novel multiscale
visual drilldown approach that was designed in tight collaboration with
proteomic experts, enabling a systematic exploration of the configuration
space. Our approach takes advantage of the hierarchical structure of the data
-- from the whole ensemble of protein complex configurations to the individual
configurations, their contact interfaces, and the interacting amino acids. Our
new solution is based on interactively linked 2D and 3D views for individual
hierarchy levels and at each level, we offer a set of selection and filtering
operations enabling the user to narrow down the number of configurations that
need to be manually scrutinized. Furthermore, we offer a dedicated filter
interface, which provides the users with an overview of the applied filtering
operations and enables them to examine their impact on the explored ensemble.
This way, we maintain the history of the exploration process and thus enable
the user to return to an earlier point of the exploration. We demonstrate the
effectiveness of our approach on two case studies conducted by collaborating
proteomic experts.


  Chemical reaction networks (CRNs) provide a convenient language for modelling
a broad variety of biological systems. These models are commonly studied with
respect to the time series they generate in deterministic or stochastic
simulations. Their dynamic behaviours are then analysed, often by using
deterministic methods based on differential equations with a focus on the
steady states. Here, we propose a method for comparing CRNs with respect to
their behaviour in stochastic simulations. Our method is based on using the
flux graphs that are delivered by stochastic simulations as abstract
representations of their dynamic behaviour. This allows us to compare the
behaviour of any two CRNs for any time interval, and define a notion of
equivalence on them that overlaps with graph isomorphism at the lowest level of
representation. The similarity between the compared CRNs can be quantified in
terms of their distance. The results can then be used to refine the models or
to replace a larger model with a smaller one that produces the same behaviour
or vice versa.


  Treatment recommendations within Clinical Practice Guidelines (CPGs) are
largely based on findings from clinical trials and case studies, referred to
here as research studies, that are often based on highly selective clinical
populations, referred to here as study cohorts. When medical practitioners
apply CPG recommendations, they need to understand how well their patient
population matches the characteristics of those in the study cohort, and thus
are confronted with the challenges of locating the study cohort information and
making an analytic comparison. To address these challenges, we develop an
ontology-enabled prototype system, which exposes the population descriptions in
research studies in a declarative manner, with the ultimate goal of allowing
medical practitioners to better understand the applicability and
generalizability of treatment recommendations. We build a Study Cohort Ontology
(SCO) to encode the vocabulary of study population descriptions, that are often
reported in the first table in the published work, thus they are often referred
to as Table 1. We leverage the well-used Semanticscience Integrated Ontology
(SIO) for defining property associations between classes. Further, we model the
key components of Table 1s, i.e., collections of study subjects, subject
characteristics, and statistical measures in RDF knowledge graphs. We design
scenarios for medical practitioners to perform population analysis, and
generate cohort similarity visualizations to determine the applicability of a
study population to the clinical population of interest. Our semantic approach
to make study populations visible, by standardized representations of Table 1s,
allows users to quickly derive clinically relevant inferences about study
populations.


  Personalized models of the gut microbiome are valuable for disease prevention
and treatment. For this, one requires a mathematical model that predicts
microbial community composition and the emergent behavior of microbial
communities. We seek a modeling strategy that can capture emergent behavior
when built from sets of universal individual interactions. Our investigation
reveals that species-metabolite interaction modeling is better able to capture
emergent behavior in community composition dynamics than direct species-species
modeling.
  Using publicly available data, we examine the ability of species-species
models and species-metabolite models to predict trio growth experiments from
the outcomes of pair growth experiments. We compare quadratic species-species
interaction models and quadratic species-metabolite interaction models, and
conclude that only species-metabolite models have the necessary complexity to
to explain a wide variety of interdependent growth outcomes. We also show that
general species-species interaction models cannot match patterns observed in
community growth dynamics, whereas species-metabolite models can. We conclude
that species-metabolite modeling will be important in the development of
accurate, clinically useful models of microbial communities.


  Bird collisions with high-speed transport modes is a vital topic on vehicle
safety and wildlife protection, especially when high-speed trains, with an
average speed of 250km/h, have to run across the habitat of an endangered bird
species. This paper evaluates the bird-train collision risk associated with a
recent high-speed railway project in Qinling Mountains, China, for the crested
ibis (Nipponia nippon) and other local bird species. Using line transect
surveys and walking monitoring techniques, we surveyed the population
abundance, spatial-temporal distributions, and bridge-crossing behaviors of the
birds in the study area. The results show that: (1) The crested ibis and the
egret were the two most abundant waterfowl species in the study area. The RAI
of these two species were about 43.69% and 42.91%, respectively; (2) Crested
ibises overall habitat closer to the railway bridge. 91.63% of them were
firstly detected within the range of 0m to 25m of the vicinity of the bridge;
(3) the ratio between crossing over and under the railway bridge was about 7:3.
Crested ibises were found to prefer to fly over the railway bridge (89.29% of
the total crossing activities observed for this species). Egrets were more
likely to cross the railway below the bridge, and they accounted for 60.27% of
the total observations of crossing under the bridge. We recommend that, while
the collision risk of crested ibises could be low, barrier-like structures,
such as fences, should still be considered to promote the conservation of
multiple bird species in the area. This paper provides a practical case for
railway ecology studies in China. To our best knowledge, this is the first
high-speed railway project that takes protecting crested ibises as one of the
top priorities, and exemplifies the recent nationwide initiative towards the
construction of "eco-civilization" in the country.


  In host-pathogen interactions, often the host (attacked organism) defends
itself by some toxic compound and the parasite, in turn, responds by producing
an enzyme that inactivates that compound. In some cases, the host can respond
by producing an inhibitor of that enzyme, which can be considered as a
counter-counter defence. An example is provided by cephalosporins,
beta-lactamases and clavulanic acid (an inhibitor of beta-lactamases). Here, we
tackle the question under which conditions it pays, during evolution, to
establish a counter-counter defence rather than to intensify or widen the
defence mechanisms. We establish a mathematical model describing this
phenomenon, based on enzyme kinetics for competitive inhibition. We use an
objective function based on Haber's rule, which says that the toxic effect is
proportional to the time integral of toxin concentration. The optimal
allocation of defence and counter-counter defence can be calculated in an
analytical way despite the nonlinearity in the underlying differential
equation. The calculation provides a threshold value for the dissociation
constant of the inhibitor. Only if the inhibition constant is below that
threshold, that is, in the case of strong binding of the inhibitor, it pays to
have a counter-counter defence. This theoretical prediction accounts for the
observation that not for all defence mechanisms, a counter-counter defence
exists. Our results should be of interest for computing optimal mixtures of
beta-lactam antibiotics and beta-lactamase inhibitors such as sulbactam, as
well as for plant-herbivore and other molecular-ecological interactions and to
fight antibiotic resistance in general.


  Genomics is the critical key to enabling precision medicine, ensuring global
food security and enforcing wildlife conservation. The massive genomic data
produced by various genome sequencing technologies presents a significant
challenge for genome analysis. Because of errors from sequencing machines and
genetic variations, approximate pattern matching (APM) is a must for practical
genome analysis. Recent work proposes FPGA, ASIC and even
process-in-memory-based accelerators to boost the APM throughput by
accelerating dynamic-programming-based algorithms (e.g., Smith-Waterman).
However, existing accelerators lack the efficient hardware acceleration for the
exact pattern matching (EPM) that is an even more critical and essential
function widely used in almost every step of genome analysis including
assembly, alignment, annotation and compression.
  State-of-the-art genome analysis adopts the FM-Index that augments the
space-efficient BWT with additional data structures permitting fast EPM
operations. But the FM-Index is notorious for poor spatial locality and massive
random memory accesses. In this paper, we propose a ReRAM-based
process-in-memory architecture, FindeR, to enhance the FM-Index EPM search
throughput in genomic sequences. We build a reliable and energy-efficient
Hamming distance unit to accelerate the computing kernel of FM-Index search
using commodity ReRAM chips without introducing extra CMOS logic. We further
architect a full-fledged FM-Index search pipeline and improve its search
throughput by lightweight scheduling on the NVDIMM. We also create a system
library for programmers to invoke FindeR to perform EPMs in genome analysis.
Compared to state-of-the-art accelerators, FindeR improves the FM-Index search
throughput by $83\%\sim 30K\times$ and throughput per Watt by $3.5\times\sim
42.5K\times$.


  Proteins in photosynthetic membranes can organize into patterned arrays that
span the membrane's lateral size. Attractions between proteins in different
layers of a membrane stack can play a key role in this ordering, as was
suggested by microscopy and fluorescence spectroscopy and demonstrated by
computer simulations of a coarse-grained model. The architecture of thylakoid
membranes, however, also provides opportunities for inter-layer interactions
that instead disfavor the high protein densities of ordered arrangements. Here
we explore the interplay between these opposing driving forces, and in
particular the phase transitions that emerge in the periodic geometry of
stacked thylakoid membrane discs. We propose a lattice model that roughly
accounts for proteins' attraction within a layer and across the stromal gap,
steric repulsion across the lumenal gap, and regulation of protein density by
exchange with the stroma lamellae. Mean field analysis and computer simulation
reveal rich phase behavior for this simple model, featuring a broken-symmetry
striped phase that is disrupted at both high and low extremes of chemical
potential. The resulting sensitivity of microscopic protein arrangement to the
thylakoid's mesoscale vertical structure raises intriguing possibilities for
regulation of photosynthetic function.


  The rise of antimicrobial drug resistance is an imminent threat to global
health that has warranted, and duly received, considerable attention within the
medical, microbiological and modelling communities. Outbreaks of drug-resistant
pathogens are ignited by the emergence and transmission of mutant variants
descended from wild-type strains circulating in the community. In this work we
investigate the stochastic dynamics of the emergence of a novel disease strain,
introduced into a population in which it must compete with an existing endemic
strain. In analogy with past work on single-strain epidemic outbreaks, we apply
a branching process approximation to calculate the probability that the new
strain becomes established. As expected, a critical determinant of the survival
prospects of any invading strain is the magnitude of its reproduction number
relative to that of the background endemic strain. Whilst in most circumstances
this ratio must exceed unity in order for invasion to be viable, we show that
differential control scenarios can lead to less-fit novel strains invading
populations hosting a fitter endemic one. This analysis and the accompanying
findings will inform our understanding of the mechanisms that have led to past
instances of successful strain invasion, and provide valuable lessons for
thwarting future drug-resistant strain incursions.


  Host-microbe interactions play crucial roles in marine ecosystems, but we
still have very little understanding of the mechanisms that govern these
relationships, the evolutionary processes that shape them, and their ecological
consequences. The holobiont concept is a renewed paradigm in biology that can
help describe and understand these complex systems. It posits that a host and
its associated microbiota, living together in a long-lasting relationship, form
the holobiont, and have to be studied together, as a coherent biological and
functional unit, in order to understand the biology, ecology and evolution of
the organisms. Here we discuss critical concepts and opportunities in marine
holobiont research and identify key challenges in the field. We highlight the
potential economic, sociological, and environmental impacts of the holobiont
concept in marine biological, evolutionary, and environmental sciences with
comparisons to terrestrial science whenever appropriate. A deeper understanding
of such complex systems, however, will require further technological and
conceptual advances. The most significant challenge will be to bridge
functional research on simple and tractable model systems and global
approaches. This will require scientists to work together as an (inter)active
community in order to address, for instance, ecological and evolutionary
questions and the roles of holobionts in biogeochemical cycles.


  Measures of tree balance play an important role in the analysis of
phylogenetic trees. One of the oldest and most popular indices in this regard
is the Colless index for rooted bifurcating trees, introduced by Colless
(1982). While many of its statistical properties under different probabilistic
models for phylogenetic trees have already been established, little is known
about its minimum value and the trees that achieve it. In this manuscript, we
fill this gap in the literature. To begin with, we derive both recursive and
closed expressions for the minimum Colless index of a tree with $n$ leaves.
Surprisingly, these expressions show a connection between the minimum Colless
index and the so-called Blancmange curve, a fractal curve. We then fully
characterize the tree shapes that achieve this minimum value and we introduce
both an algorithm to generate them and a recurrence to count them. After
focusing on two extremal classes of trees with minimum Colless index (the
maximally balanced trees and the greedy from the bottom trees), we conclude by
showing that all trees with minimum Colless index also have minimum Sackin
index, another popular balance index.


  We consider the spread of a Susceptible-Infected-Recovered (SIR) disease
through finite populations and derive an expression for the final size
distribution. Our derivation allows arbitrary distributions of the number of
transmissions caused by an infected individual. We show how this calculation
can be used to infer parameters of the infectious disease through observations
in multiple small populations. The inference suffers from some identifiability
difficulties, and it requires many observations to distinguish between
parameter combinations that correspond to the same reproductive number.


  Rock-paper-scissors games metaphorically model cyclic dominance in ecology
and microbiology. In a static environment, these models are characterized by
fixation probabilities obeying two different "laws" in large and small
well-mixed populations. Here, we investigate the evolution of these
three-species models subject to a randomly switching carrying capacity modeling
the endless change between states of resources scarcity and abundance. Focusing
mainly on the zero-sum rock-paper-scissors game, equivalent to the cyclic
Lotka-Volterra model, we study how the ${\it coupling}$ of demographic and
environmental noise influences the fixation properties. More specifically, we
investigate which species is the most likely to prevail in a population of
fluctuating size and how the outcome depends on the environmental variability.
We show that demographic noise coupled with environmental randomness "levels
the field" of cyclic competition by balancing the effect of selection. In
particular, we show that fast switching effectively reduces the selection
intensity proportionally to the variance of the carrying capacity. We determine
the conditions under which new fixation scenarios arise, where the most likely
species to prevail changes with the rate of switching and the variance of the
carrying capacity. Random switching has a limited effect on the mean fixation
time that scales linearly with the average population size. Hence,
environmental randomness makes the cyclic competition more egalitarian, but
does not prolong the species coexistence. We also show how the fixation
probabilities of close-to-zero-sum rock-paper-scissors games can be obtained
from those of the zero-sum model by rescaling the selection intensity.


  Clinical biosensors with low detection limit hold significant promise in the
early diagnosis of debilitating diseases. Recent progress in sensor development
has led to the demonstration of detection capable of detecting target molecules
even down to single-molecule level. One crucial performance parameter which is
not adequately discussed is the issue of measurement fidelity in such sensors.
We define measurement fidelity based on the false positive rate of the system
as we expect systems with higher sensitivity to concomitantly respond more to
interfering molecules thus increasing the false-positive rates. We present a
model which allows us to estimate the limit of detection of a biosensor system
constrained by a specified false-positive rate. Two major results emerging from
our model is that a) there is a lower bound to the detection limit for a target
molecule determined by the variation in the concentration of background
molecules interfering with the molecular recognition process and b) systems
which use a secondary label, such as a fluorophore, can achieve lower detection
limits for a given false positive rate. We also present data collected from
literature to support our model. The insights from our model will be useful in
the systematic design of future clinical biosensors to achieve relevant
detection limits with assured fidelity.


  Artificially engineered biosensors are highly inefficient in accurately
measuring the concentration of biomarkers, particularly, during early diagnosis
of diseases. On the other hand, single cellular systems such as chemotactic
bacteria can sense their environment with extraordinary precision. Therefore,
one would expect that implementing the optimal cellular sensing strategies in
state-of-the-art artificial sensors can produce optimally precise biosensors.
However because of the presence of measurement noise, strategies that are
optimal in biological systems may not be optimal in artificial systems.
Therefore, mimicking biological strategies may not be the optimal path in case
of artificial sensing systems because of the presence of inherent measurement
noise.


  Single-particle reconstruction in cryo-electron microscopy (cryo-EM) is an
increasingly popular technique for determining the 3-D structure of a molecule
from several noisy 2-D projections images taken at unknown viewing angles. Most
reconstruction algorithms require a low-resolution initialization for the 3-D
structure, which is the goal of ab initio modeling. Suggested by Zvi Kam in
1980, the method of moments (MoM) offers one approach, wherein low-order
statistics of the 2-D images are computed and a 3-D structure is estimated by
solving a system of polynomial equations. Unfortunately, Kam's method suffers
from restrictive assumptions, most notably that viewing angles should be
distributed uniformly. Often unrealistic, uniformity entails the computation of
higher-order correlations, as in this case first and second moments fail to
determine the 3-D structure. In the present paper, we remove this hypothesis,
by permitting an unknown, non-uniform distribution of viewing angles in MoM.
Perhaps surprisingly, we show that this case is statistically easier than the
uniform case, as now first and second moments generically suffice to determine
low-resolution expansions of the molecule. In the idealized setting of a known,
non-uniform distribution, we find an efficient provable algorithm inverting
first and second moments. For unknown, non-uniform distributions, we use
non-convex optimization methods to solve for both the molecule and
distribution.


  A common trait of complex systems is that they can be represented by means of
a network of interacting parts. It is, in fact, the network organisation (more
than the parts) what largely conditions most higher-level properties, which are
not reducible to the properties of the individual parts. Can the topological
organisation of these webs provide some insight into their evolutionary
origins? Both biological and artificial networks share some common
architectural traits. They are often heterogeneous and sparse, and most exhibit
the small-world property or the presence of modular or hierarchical patterns.
These properties have often been attributed to the selection of functionally
meaningful traits. However, the proper formulation of generative network models
suggests a rather different picture. Against the standard
selection-optimisation argument, some networks reveal the inevitable generation
of complex patterns resulting from reuse and can be modelled using
duplication-rewiring rules. These give rise to the heterogeneous, scale-free
and modular architectures observed in the real case studies. Tinkering is a
universal mechanism that drives not only biological evolution but also the
large-scale dynamics of some technological designs. Here we examine the
evidence for tinkering in cellular, technological and ecological webs and its
impact in shaping their architecture and deeply affecting their functional
properties. Our analysis suggests to seriously reconsider the role played by
selection forces or design principles as main drivers of network evolution.


  One of the most progressive methods of cancer treatment is ion beam therapy.
The simulations of water radiolysis show that in the cell medium the most
long-living species are hydrogen peroxide (H$_2$O$_2$) molecules. But up to the
present time the role of H$_2$O$_2$ molecules in the deactivation of cancer
cells has not been determined yet. To understand the possible role of
H$_2$O$_2$ in ion beam therapy, the competitive interaction of H$_2$O$_2$ and
H$_2$O molecules with nucleic bases in a pair on the different stages of
genetic information transfer is studied in the present work. Atom-atomic
potential functions method is used for the calculations. It is shown, that some
configurations of AT and GC complementary pairs are stabilized much better by
H$_2$O$_2$ molecule compared to water molecule. Formation of such interaction
complexes can terminate the processes of DNA unzipping by enzymes and
consequently block genetic information transfer processes in cancer cells
during ion beam therapy treatment. Experimental verification method of hydrogen
peroxide interaction with nucleic base pairs is proposed.