  We investigate dynamical coupling between water and amino acid side-chain
residues in solvation dynamics by selecting residues often used as natural
probes, namely tryptophan, tyrosine and histidine, located at different
positions on protein surface and having various degrees of solvent exposure.
Such differently placed residues are found to exhibit different timescales of
relaxation. The total solvation response, as measured by the probe is
decomposed in terms of its interactions with (i) protein core, (ii) side-chain
atoms and (iii) water molecules. Significant anti cross-correlations among
these contributions are observed as a result of side-chain assisted energy flow
between protein core and hydration layer, which is important for the proper
functionality of a protein. It is also observed that there are rotationally
faster as well as slower water molecules than that of bulk solvent, which are
considered to be responsible for the multitude of timescales that are observed
in solvation dynamics. We also establish that slow solvation derives a
significant contribution from protein side-chain fluctuations. When the motion
of the protein side-chains is forcefully quenched, solvation either becomes
faster or slower depending on the location of the probe.


  About half of human cancers show normal TP53 gene and aberrant overexpression
of Mdm2 and/or MdmX. This fact promotes a promising cancer therapeutic strategy
which targeting the interactions between p53 and Mdm2/MdmX. For developing the
inhibitors to disrupt the p53-Mdm2/MdmX interactions, we systematically
investigate structural and interaction characteristics of p53 and inhibitors
with Mdm2 and MdmX from atomistic level by exploiting stochastic molecular
dynamics simulations. We find that some specific $\alpha$ helices in Mdm2 and
MdmX structure play key role in their bindings with inhibitors and the hydrogen
bond formed by residue Trp23 of p53 with its counterpart in Mdm2/MdmX
determines dynamical competition processes of the disruption of Mdm2-p53
interaction and replacement of p53 from Mdm2-p53 complex {\it in vivo}. We hope
that the results reported in this paper provide basic information for designing
functional inhibitors and realizing cancer gene therapy.


  The principles behind the computation of protein-ligand binding free energies
by Monte Carlo integration are described in detail. The simulation provides
gas-phase binding free energies that can be converted to aqueous energies by
solvation corrections. The direct integration simulation has several
characteristics beneficial to free-energy calculations. One is that the number
of parameters that must be set for the simulation is small and can be
determined objectively, making the outcome more deterministic, with respect to
choice of input conditions, as compared to perturbation methods. Second, the
simulation is free from assumptions about the starting pose or nature of the
binding site. A final benefit is that binding free energies are a direct
outcome of the simulation, and little processing is required to determine them.
  The well-studied T4 lysozyme experimental free energy data and crystal
structures were used to evaluate the method.


  The rotary sequential hydrolysis of metabolic machine F1-ATPase is a
prominent feature to reveal high coordination among multiple chemical sites on
the stator F1 ring, which also contributes to tight coupling between the
chemical reaction and central {\gamma}-shaft rotation. High-speed AFM
experiments discovered that the sequential hydrolysis was maintained on the F1
ring even in the absence of the {\gamma} rotor. To explore how the intrinsic
sequential performance arises, we computationally investigated essential
inter-subunit couplings on the hexameric ring of mitochondrial and bacterial
F1. We first reproduced the sequential hydrolysis schemes as experimentally
detected, by simulating tri-site ATP hydrolysis cycles on the F1 ring upon
kinetically imposing inter-subunit couplings to substantially promote the
hydrolysis products release. We found that it is key for certain ATP binding
and hydrolysis events to facilitate the neighbor-site ADP and Pi release to
support the sequential hydrolysis. The kinetically feasible couplings were then
scrutinized through atomistic molecular dynamics simulations as well as
coarse-grained simulations, in which we enforced targeted conformational
changes for the ATP binding or hydrolysis. Notably, we detected the
asymmetrical neighbor-site opening that would facilitate the ADP release upon
the enforced ATP binding, and computationally captured the complete Pi release
through charge hopping upon the enforced neighbor-site ATP hydrolysis. The
ATP-hydrolysis triggered Pi release revealed in current TMD simulation confirms
a recent prediction made from statistical analyses of single molecule
experimental data in regard to the role ATP hydrolysis plays. Our studies,
therefore, elucidate both the concerted chemical kinetics and underlying
structural dynamics of the inter-subunit couplings that lead to the rotary
sequential hydrolysis of the F1 ring.


  We develop a variable population age-structured ODE model to investigate the
role of Intermittent Preventive Treatment (IPT) in averting malaria-induced
mortalities in children, as well as its related cost in promoting the spread of
anti-malarial drug resistance. IPT, a malaria control strategy in which a full
curative dose of an antimalarial medication is administered to vulnerable
asymptomatic individuals at specified intervals, has been shown to have a
positive impact on reducing malaria transmission and deaths in children and
pregnant women. However, it can also promote drug resistance spread. Our
mathematical model is used to explore IPT effects on drug resistance in
holoendemic malaria regions while quantifying the benefits in deaths averted.
Our model includes both drug-sensitive and drug-resistant strains of the
parasite as well as interactions between human hosts and mosquitoes. The basic
reproduction numbers for both strains as well as the invasion reproduction
numbers are derived and used to examine the role of IPT on drug resistance.
Numerical simulations show the individual and combined effects of IPT and
treatment of symptomatic infections on the prevalence levels of both parasite
strains and on the number of lives saved. The results suggest that while IPT
can indeed save lives, particularly in the high transmission region, certain
combinations of drugs used for IPT and drugs used to treat symptomatic
infection may result in more deaths when resistant parasite strains are
circulating. Moreover, the half-lives of the treatment and IPT drugs used play
an important role in the extent to which IPT may influence the rate of spread
of the resistant strain. A sensitivity analysis indicates the model outcomes
are most sensitive to the reduction factor of transmission for the resistant
strain, rate of immunity loss, and the clearance rate of sensitive infections.


  Population dynamics with spatial information is applied to understand the
spread of pests. We introduce a model describing how pests spread in discrete
space. The number of pest descendants at each site is controlled by local
information such as temperature, precipitation, and the density of pine trees.
Our simulation leads to a pest spreading pattern comparable to the real data
for pine needle gall midge in the past. We also simulate the model in two
different climate conditions based on two different representative
concentration pathways scenarios for the future. We observe that after an
initial stage of a slow spread of pests, a sudden change in the spreading speed
occurs, which is soon followed by a large-scale outbreak. We found that a
future climate change causes the outbreak point to occur earlier and that the
detailed spatio-temporal pattern of the spread depends on the source position
from which the initial pest infection starts.


  Falciparum malaria is a major parasitic disease causing widespread morbidity
and mortality globally. Artemisinin derivatives---the most effective and
widely-used antimalarials that have helped reduce the burden of malaria by 60%
in some areas over the past decade---have recently been found to induce growth
retardation of blood-stage Plasmodium falciparum when applied at clinically
relevant concentrations. To date, no model has been designed to quantify the
growth retardation effect and to predict the influence of this property on in
vivo parasite killing. Here we introduce a mechanistic model of parasite growth
from the ring to trophozoite stage of the parasite's life cycle, and by
modelling the level of staining with an RNA-binding dye, we demonstrate that
the model is able to reproduce fluorescence distribution data from in vitro
experiments using the laboratory 3D7 strain. We quantify the dependence of
growth retardation on drug concentration and demonstrate the model's utility as
a platform to propose experimentally-testable mechanisms of growth retardation.
Furthermore we illustrate that a drug-induced delay in growth may significantly
influence in vivo parasite dynamics, demonstrating the importance of
considering growth retardation in the design of optimal artemisinin-based
dosing regimens.


  Current treatment strategy for chronic myeloid leukemia (CML) mainly includes
inhibition of tyrosine kinase activity, which has dramatically improved the
prognosis of the disease but without cure. In addition some patients may become
drug resistant. Thus there is still the need for other therapies to avoid
resistance and if possible to cure the disease. Loss of p53 is known to play an
important role in the disease progression of CML and causes drug resistance.
Here I propose that in CML, inositol requiring enzyme 1 alpha (IRE1 alpha) may
cause abnormal degradation of p53 mRNA resulting in inhibition of apoptosis in
leukemic clonal cells, which has not been elucidated before. Hence, I propose
that inhibition of endoribonuclease activity of IRE1 alpha with small molecule
inhibitors may provide a novel strategy to enhance p53 function in CML leukemic
clones to overcome the limitations of current treatment regimens.


  Competition to bind microRNAs induces an effective positive crosstalk between
their targets, therefore known as `competing endogenous RNAs' or ceRNAs. While
such an effect is known to play a significant role in specific conditions,
estimating its strength from data and, experimentally, in physiological
conditions appears to be far from simple. Here we show that the susceptibility
of ceRNAs to different types of perturbations affecting their competitors (and
hence their tendency to crosstalk) can be encoded in quantities as intuitive
and as simple to measure as correlation functions. We confirm this scenario by
extensive numerical simulations and validate it by re-analyzing PTEN's
crosstalk pattern from TCGA breast cancer dataset. These results clarify the
links between different quantities used to estimate the intensity of ceRNA
crosstalk and provide new keys to analyze transcriptional datasets and
effectively probe ceRNA networks in silico.


  Cells contain elaborate and interconnected networks of protein polymers which
make up the cytoskeleton. The cytoskeleton governs the internal positioning and
movement of vesicles and organelles, and controls dynamic changes in cell
polarity, shape and movement. Many of these processes require tight control of
the size and shape of cytoskeletal structures, which is achieved despite rapid
turnover of their molecular components. Here we review mechanisms by which
cells control the size of filamentous cytoskeletal structures from the point of
view of simple quantitative models that take into account stochastic dynamics
of their assembly and disassembly. Significantly, these models make
experimentally testable predictions that distinguish different mechanisms of
length-control. While the primary focus of this review is on cytoskeletal
structures, we believe that the broader principles and mechanisms discussed
herein will apply to a range of other subcellular structures whose sizes are
tightly controlled and are linked to their functions.


  Motivated by applications in cancer genomics and following the work of
Hajirasouliha and Raphael (WABI 2014), Hujdurovi\'c et al. (IEEE TCBB, to
appear) introduced the minimum conflict-free row split (MCRS) problem: split
each row of a given binary matrix into a bitwise OR of a set of rows so that
the resulting matrix corresponds to a perfect phylogeny and has the minimum
possible number of rows among all matrices with this property. Hajirasouliha
and Raphael also proposed the study of a similar problem, in which the task is
to minimize the number of distinct rows of the resulting matrix. Hujdurovi\'c
et al. proved that both problems are NP-hard, gave a related characterization
of transitively orientable graphs, and proposed a polynomial-time heuristic
algorithm for the MCRS problem based on coloring cocomparability graphs.
  We give new, more transparent formulations of the two problems, showing that
the problems are equivalent to two optimization problems on branchings in a
derived directed acyclic graph. Building on these formulations, we obtain new
results on the two problems, including: (i) a strengthening of the heuristic by
Hujdurovi\'c et al. via a new min-max result in digraphs generalizing
Dilworth's theorem, which may be of independent interest, (ii) APX-hardness
results for both problems, (iii) approximation algorithms, and (iv)
exponential-time algorithms solving the two problems to optimality faster than
the na\"ive brute-force approach. Our work relates to several well studied
notions in combinatorial optimization: chain partitions in partially ordered
sets, laminar hypergraphs, and (classical and weighted) colorings of graphs.


  Signal transducers and activators of transcription (STATs) are key molecular
determinants of T cell fate and effector function. A number of inflammatory
diseases are characterized by an altered balance of T cell phenotypes and
cytokine secretion. STATs, therefore, represent viable therapeutic targets in
numerous pathologies. However, the underlying mechanisms of how the same STAT
proteins regulate both the development of different T cell phenotypes and their
plasticity during changes in extracellular conditions remain unclear. In this
study, we investigated the STAT mediated regulation of T cell phenotype
formation and plasticity using mathematical modeling and experimental data for
intracellular STAT signaling proteins. The close fit of our model predictions
to the experimental data for IFN-{\gamma} to IL-10 switching allows us to
propose a potential mechanism for T cell switching that regulates human Th1/Tr1
responses. According to this mechanism, T cell phenotype switching is due to
the relative redistribution of STAT dimer complexes caused by the extracellular
cytokine-dependent STAT competition effects. The proposed model is applicable
to a number of STAT signaling circuits.


  Dengue fever is increasing in geographical range, spread by invasion of its
vector mosquitoes. The trade in second-hand tires has been implicated as a
factor in this process as they act as mobile reservoirs of mosquito eggs and
larvae. Regional transportation of tires can create linkages between rural
areas with dengue to disease-free urban areas, potentially giving rise to
outbreaks even in areas with strong local control measures. In this work we
sought to model the dynamics of mosquito transportation via the tire trade, in
particular to predict its role in causing unexpected dengue outbreaks through
vertical transmission of the virus across generations of mosquitoes. We also
aimed to identify strategies for regulating the trade in second-hand tires,
improving disease control. We created a mathematical model which captures the
dynamics of dengue between rural and urban areas, taking into account the
movement, storage time of tires, and mosquito diapause. We simulate a series of
scenarios. First a mosquito population is introduced to a dengue-free area via
movement of tires, either as single or multiple events, increasing the
likelihood of a dengue outbreak. An endemic state can be induced regardless of
whether urban conditions for an outbreak are met, and an existing endemic state
can be enhanced by vector input. Finally we assess the potential for regulation
of tire processing as a means of reducing the transmission of dengue fever
using a specific case study from Puerto Rico. Our work demonstrates the
importance of the second-hand tire trade in modulating the spread of dengue
fever across regions, in particular its role in introducing dengue to
disease-free areas. We propose that regulation of tire storage and movement can
play a crucial role in containing outbreaks and dengue spread.


  Using a structure-based coarse-grained model of proteins, we study the
mechanism of unfolding of knotted proteins through heating. We find that the
dominant mechanisms of unfolding depend on the temperature applied and are
generally distinct from those identified for folding at its optimal
temperature. In particular, for shallowly knotted proteins, folding usually
involves formation of two loops whereas unfolding through high-temperature
heating is dominated by untying of single loops. Untying the knots is found to
generally precede unfolding unless the protein is deeply knotted and the
heating temperature exceeds a threshold value. We then use a phenomenological
model of the air-water interface to show that such an interface can untie
shallow knots, but it can also make knots in proteins that are natively
unknotted.


  We re-examined data from the classic Luria-Delbruck fluctuation experiment,
which is often credited with establishing a Darwinian basis for evolution. We
argue that, for the Lamarckian model of evolution to be ruled out by the
experiment, the experiment must favor pure Darwinian evolution over both the
Lamarckian model and a model that allows both Darwinian and Lamarckian
mechanisms. Analysis of the combined model was not performed in the original
1943 paper. The Luria-Delbruck paper also did not consider the possibility of
neither model fitting the experiment. Using Bayesian model selection, we find
that the Luria-Delbruck experiment, indeed, favors the Darwinian evolution over
purely Lamarckian. However, our analysis does not rule out the combined model,
and hence cannot rule out Lamarckian contributions to the evolutionary
dynamics.


  Since penicillin was discovered about 90 years ago, we have become used to
using drugs to eradicate unwanted pathogenic cells. However, using drugs to
kill bacteria, viruses or cancer cells has the serious side effect of selecting
for mutant types that survive the drug attack. A key question therefore is how
one could eradicate as many cells as possible for a given acceptable risk of
drug resistance evolution. We address this general question in a model of drug
resistance evolution in spatial drug gradients, which recent experiments and
theories have suggested as key drivers of drug resistance. Importantly, our
model takes into account the influence of convection, resulting for instance
from blood flow. Using stochastic simulations, we study the fates of individual
resistance mutations and quantify the trade-off between the killing of
wild-type cells and the rise of resistance mutations: shallow gradients and
convection into the antibiotic region promote wild-type death, at the cost of
increasing the establishment probability of resistance mutations. We can
explain these observed trends by modeling the adaptation process as a branching
random walk. Our analysis reveals that the trade-off between death and
adaptation depends on the relative length scales of the spatial drug gradient
and random dispersal, and the strength of convection. Our results show that
convection can have a momentous effect on the rate of establishment of new
mutations, and may heavily impact the efficiency of antibiotic treatment.


  Gene transcription is a critical step in gene expression. The currently
accepted physical model of transcription predicts the existence of a physical
limit on the maximal rate of transcription, which does not depend on the
transcribed gene. This limit appears as a result of polymerase "traffic jams"
forming in the bulk of the 1D DNA chain at high polymerase concentrations.
Recent experiments have, for the first time, allowed one to access live gene
expression dynamics in the Drosophila fly embryo in vivo under the conditions
of heavy polymerase load and test the predictions of the model. Our analysis of
the data shows that the maximal rate of transcription is indeed the same for
the Hunchback, Snail and Knirps gap genes, and modified gene constructs in
nuclear cycles 13, 14, but the experimentally observed value of the maximal
transcription rate corresponds to only 40 % of the one predicted by this model.
We argue that such a decrease must be due to a slower polymerase elongation
rate in the vicinity of the promoter region. This effectively shifts the
bottleneck of transcription from the bulk to the promoter region of the gene.
We suggest a quantitative explanation of the difference by taking into account
abortive transcription initiation. Our calculations based on the independently
measured abortive initiation constant in vitro confirm this hypothesis and find
quantitative agreement with MS2 fluorescence live imaging data in the early
fruit fly embryo. If our explanation is correct, then the transcription rate
cannot be increased by replacing "slow codons" in the bulk with synonymous
codons, and experimental efforts must be focused on the promoter region
instead. This study extends our understanding of transcriptional regulation,
re-examines physical constraints on the kinetics of transcription and
re-evaluates the validity of the standard TASEP model of transcription.


  We present a framework for computing the gating properties of ligand-gated
ion channel mutants using the Monod-Wyman-Changeux (MWC) model of allostery. We
derive simple analytic formulas for key functional properties such as the
leakiness, dynamic range, half-maximal effective concentration, and effective
Hill coefficient, and explore the full spectrum of phenotypes that are
accessible through mutations. Specifically, we consider mutations in the
channel pore of nicotinic acetylcholine receptor (nAChR) and the ligand binding
domain of a cyclic nucleotide-gated (CNG) ion channel, demonstrating how each
mutation can be characterized as only affecting a subset of the biophysical
parameters. In addition, we show how the unifying perspective offered by the
MWC model allows us, perhaps surprisingly, to collapse the plethora of
dose-response data from different classes of ion channels into a universal
family of curves.


  Despite recent advances in reputation technologies, it is not clear how
reputation systems can affect human cooperation in social networks. Although it
is known that two of the major mechanisms in the evolution of cooperation are
spatial selection and reputation-based reciprocity, theoretical study of the
interplay between both mechanisms remains almost uncharted. Here, we present a
new individual-based model for the evolution of reciprocal cooperation between
reputation and networks. We comparatively analyze four of the leading moral
assessment rules---shunning, image scoring, stern judging, and simple
standing---and base the model on the giving game in regular networks for
Cooperators, Defectors, and Discriminators. Discriminators rely on a proper
moral assessment rule. By using individual-based models, we show that the four
assessment rules are differently characterized in terms of how cooperation
evolves, depending on the benefit-to-cost ratio, the network-node degree, and
the observation and error conditions. Our findings show that the most tolerant
rule---simple standing---is the most robust among the four assessment rules in
promoting cooperation in regular networks.


  We study the coarsening of strongly microphase separated membrane domains in
the presence of recycling of material. We study the dynamics of the domain size
distribution under both scale-free and size-dependent recycling. Closed form
solutions to the steady state distributions and its associated central moments
are obtained in both cases. Moreover, for the size-independent case, the~time
evolution of the moments is analytically calculated, which provide us with
exact results for their corresponding relaxation times. Since these moments and
relaxation times are measurable quantities, the biophysically significant free
parameters in our model may be determined by comparison with experimental data.


  In Francis and Steel (2015), it was shown that there exists non-trivial
networks on $4$ leaves upon which the distance metric affords a metric on a
tree which is not the base tree of the network. In this paper we extend this
result in two directions. We show that for any tree $T$ there exists a family
of non-trivial HGT networks $N$ for which the distance metric $d_N$ affords a
metric on $T$. We additionally provide a class of networks on any number of
leaves upon which the distance metric affords a metric on a tree which is not
the base tree of the network.
  The family of networks are all "floating" networks, a subclass of a novel
family of networks introduced in this paper, and referred to as "versatile"
networks. Versatile networks are then characterised.
  Additionally, we find a lower bound for the number of `useful' HGT arcs in
such networks, in a sense explained in the paper. This lower bound is equal to
the number of HGT arcs required for each floating network in the main results,
and thus our networks are minimal in this sense.


  Cancer complexome comprises a heterogeneous and multifactorial milieu that
varies in cytology, physiology, signaling mechanisms and response to therapy.
The combined framework of network theory and spectral graph theory along with
the multilayer anal- ysis provides a comprehensive approach to analyze the
proteomic data of seven different cancers, namely, breast, oral, ovarian,
cervical, lung, colon and prostate. Our analysis demonstrates that the
protein-protein interaction networks of the normal and the cancerous tissues
associated with the seven cancers have overall similar structural and spectral
properties. However, few of these properties implicate unsystematic changes
from the normal to the disease networks depicting difference in the
interactions and highlighting changes in the complexity of different cancers.
Importantly, analysis of common proteins of all the cancer networks reveals few
proteins namely the sensors, which not only occupy significant position in all
the layers but also have direct involvement in causing cancer. The prediction
and analysis of miRNAs targeting these sensor proteins hint towards the
possible role of these proteins in tumorigenesis. This novel approach helps in
understanding cancer at the fundamental level and provides a clue to develop
promising and nascent concept of single drug therapy for multiple diseases as
well as personalized medicine.


  Evolutionary games on graphs describe how strategic interactions and
population structure determine evolutionary success, quantified by the
probability that a single mutant takes over a population. Graph structures,
compared to the well-mixed case, can act as amplifiers or suppressors of
selection by increasing or decreasing the fixation probability of a beneficial
mutant. Properties of the associated mean fixation times can be more intricate,
especially when selection is strong. The intuition is that fixation of a
beneficial mutant happens fast (in a dominance game), that fixation takes very
long (in a coexistence game), and that strong selection eliminates demographic
noise. Here we show that these intuitions can be misleading in structured
populations. We analyze mean fixation times on the cycle graph under strong
frequency-dependent selection for two different microscopic evolutionary update
rules (death-birth and birth-death). We establish exact analytical results for
fixation times under strong selection, and show that there are coexistence
games in which fixation occurs in time polynomial in population size. Depending
on the underlying game, we observe inherence of demographic noise even under
strong selection, if the process is driven by random death before selection for
birth of an offspring (death-birth update). In contrast, if selection for an
offspring occurs before random removal (birth-death update), strong selection
can remove demographic noise almost entirely.


  Chemical or enzymatic cross-linking of casein micelles (CMs) increases their
stability against dissociating agents. In this paper, a comparative study of
stability between native CMs and CMs cross-linked with genipin (CMs-GP) as a
function of pH is described. Stability to temperature and ethanol were
investigated in the pH range 2.0-7.0. The size and the charge
($\zeta$-potential) of the particles were determined by dynamic light
scattering. Native CMs precipitated below pH 5.5, CMs-GP precipitated from pH
3.5 to 4.5, whereas no precipitation was observed at pH 2.0-3.0 or pH 4.5-7.0.
The isoelectric point of CMs-GP was determined to be pH 3.7. Highest stability
against heat and ethanol was observed for CMs-GP at pH 2, where visible
coagulation was determined only after 800 s at 140 $^\circ$C or 87.5% (v/v) of
ethanol. These results confirmed the hypothesis that cross-linking by GP
increased the stability of CMs.


  While interdependent systems have usually been associated with increased
fragility, we show that strengthening the interdependence between dynamical
processes on different networks can make them more robust. By coupling the
dynamics of networks that in isolation exhibit catastrophic collapse with
extinction of nodal activity, we demonstrate system-wide persistence of
activity for an optimal range of interdependence between the networks. This is
related to the appearance of attractors of the global dynamics comprising
disjoint sets ("islands") of stable activity.


  Joint quantification of genetic and epigenetic effects on gene expression is
important for understanding the establishment of complex gene regulation
systems in living organisms. In particular, genomic imprinting and maternal
effects play important roles in the developmental process of mammals and
flowering plants. However, the influence of these effects on gene expression
are difficult to quantify because they act simultaneously with cis-regulatory
mutations. Here we propose a simple method to decompose cis-regulatory (i.e.,
allelic genotype, AG), genomic imprinting (i.e., parent-of-origin, PO), and
maternal (i.e., maternal genotype, MG) effects on allele-specific gene
expression using RNA-seq data obtained from reciprocal crosses. We evaluated
the efficiency of method using a simulated dataset and applied the method to
whole-body Drosophila and mouse trophoblast stem cell (TSC) and liver RNA-seq
data. Consistent with previous studies, we found little evidence of PO and MG
effects in adult Drosophila samples. In contrast, we identified dozens and
hundreds of mouse genes with significant PO and MG effects, respectively.
Interestingly, a similar number of genes with significant PO effect were detect
in mouse TSCs and livers, whereas more genes with significant MG effect were
observed in livers. Further application of this method will clarify how these
three effects influence gene expression levels in different tissues and
developmental stages, and provide novel insight into the evolution of gene
expression regulation.


  We introduce an in silico model for the initial spread of an aberrant
phenotype with Warburg-like overflow metabolism within a healthy homeostatic
tissue in contact with a nutrient reservoir (the blood), aimed at
characterizing the role of the microenvironment for aberrant growth. Accounting
for cellular metabolic activity, competition for nutrients, spatial diffusion
and their feedbacks on aberrant replication and death rates, we obtain a phase
portrait where distinct asymptotic whole-tissue states are found upon varying
the tissue-blood turnover rate and the level of blood-borne primary nutrient.
Over a broad range of parameters, the spreading dynamics is bistable as random
fluctuations can impact the final state of the tissue. Such a behaviour turns
out to be linked to the re-cycling of overflow products by non-aberrant cells.
Quantitative insight on the overall emerging picture is provided by a spatially
homogeneous version of the model.


  Coalescent models of evolution account for incomplete lineage sorting by
specifying a species tree parameter which determines a distribution on gene
trees. It has been shown that the unrooted topology of the species tree
parameter of the multispecies coalescent is generically identifiable. Moreover,
a statistically consistent reconstruction method called SVDQuartets has been
developed to recover this parameter. In this paper, we describe a modified
multispecies coalescent model that allows for varying effective population size
and violations of the molecular clock. We show that the unrooted topology of
the species tree for these models is generically identifiable and that
SVDQuartets is still a statistically consistent method for inferring this
parameter.


  An ecosystem made of nutrients, plants, detritus and dissolved oxygen is
presented. Its equilibria are established. Sufficient conditions for the
existence of the coexistence equilibrium are derived and its feasibility is
discussed in every detail.


  Bayesian networks can represent directed gene regulations and therefore are
favored over co-expression networks. However, hardly any Bayesian network study
concerns the false discovery control (FDC) of network edges, leading to low
accuracies due to systematic biases from inconsistent false discovery levels in
the same study. We design four empirical tests to examine the FDC of Bayesian
networks from three p-value based lasso regression variable selections --- two
existing and one we originate. Our method, lassopv, computes p-values for the
critical regularization strength at which a predictor starts to contribute to
lasso regression. Using null and Geuvadis datasets, we find that lassopv
obtains optimal FDC in Bayesian gene networks, whilst existing methods have
defective p-values. The FDC concept and tests extend to most network inference
scenarios and will guide the design and improvement of new and existing
methods. Our novel variable selection method with lasso regression also allows
FDC on other datasets and questions, even beyond network inference and
computational biology. Lassopv is implemented in R and freely available at
https://github.com/lingfeiwang/lassopv and
https://cran.r-project.org/package=lassopv


  Amino-acid substitutions are implicated in a wide range of human diseases,
many of which are lethal. Distinguishing such mutations from polymorphisms
without significant effect on human health is a necessary step in understanding
the etiology of such diseases. Computational methods can be used to select
interesting mutations within a larger set, to corroborate experimental findings
and to elucidate the cause of the deleterious effect. In this work, we show
that taking into account the sequence context in which the mutation appears
allows to improve the predictive and explanatory power of such methods. We
present an unsupervised approach based on the direct-coupling analysis of
homologous proteins. We show its capability to quantify mutations where methods
without context dependence fail. We highlight cases where the context
dependence is interpretable as functional or structural constraints and show
that our simple and unsupervised method has an accuracy similar to
state-of-the-art methods, including supervised ones.


  We consider a general branching population where the lifetimes of individuals
are i.i.d.\ with arbitrary distribution and where each individual gives birth
to new individuals at Poisson times independently from each other. In addition,
we suppose that individuals experience mutations at Poissonian rate $\theta$
under the infinitely many alleles assumption assuming that types are
transmitted from parents to offspring. This mechanism leads to a partition of
the population by type, called the allelic partition. The main object of this
work is the frequency spectrum $A(k,t)$ which counts the number of families of
size $k$ in the population at time $t$. The process $(A(k,t),\
t\in\mathbb{R}_+)$ is an example of non-Markovian branching process belonging
to the class of general branching processes counted by random characteristics.
In this work, we propose methods of approximation to replace the frequency
spectrum by simpler quantities. Our main goal is study the asymptotic error
made during these approximations through central limit theorems. In a last
section, we perform several numerical analysis using this model, in particular
to analyze the behavior of one of these approximations with respect to Sabeti's
Extended Haplotype Homozygosity [18].


  Ordered chains (such as chains of amino acids) are ubiquitous in biological
cells, and these chains perform specific functions contingent on the sequence
of their components. Using the existence and general properties of such
sequences as a theoretical motivation, we study the statistical physics of
systems whose state space is defined by the possible permutations of an ordered
list, i.e., the symmetric group, and whose energy is a function of how certain
permutations deviate from some chosen correct ordering. Such a non-factorizable
state space is quite different from the state spaces typically considered in
statistical physics systems and consequently has novel behavior in systems with
interacting and even non-interacting Hamiltonians. Various parameter choices of
a mean-field model reveal the system to contain five different physical regimes
defined by two transition temperatures, a triple point, and a quadruple point.
Finally, we conclude by discussing how the general analysis can be extended to
state spaces with more complex combinatorial properties and to other standard
questions of statistical mechanics models.


  Social conventions govern countless behaviors all of us engage in every day,
from how we greet each other to the languages we speak. But how can shared
conventions emerge spontaneously in the absence of a central coordinating
authority? The Naming Game model shows that networks of locally interacting
individuals can spontaneously self-organize to produce global coordination.
Here, we provide a gentle introduction to the main features of the model, from
the dynamics observed in homogeneously mixing populations to the role played by
more complex social networks, and to how slight modifications of the basic
interaction rules give origin to a richer phenomenology in which more
conventions can co-exist indefinitely.


  The history of gene families - which are equivalent to \emph{event-labeled}
gene trees - can be reconstructed from empirically estimated evolutionary
event-relations containing pairs of orthologous, paralogous or xenologous
genes. The question then arises as whether inferred event-labeled gene trees
are \emph{biologically feasible}, that is, if there is a possible true history
that would explain a given gene tree. In practice, this problem is boiled down
to finding a reconciliation map - also known as DTL-scenario - between the
event-labeled gene trees and a (possibly unknown) species tree.
  In this contribution, we first characterize whether there is a valid
reconciliation map for binary event-labeled gene trees $T$ that contain
speciation, duplication and horizontal gene transfer events and some unknown
species tree $S$ in terms of "informative" triples that are displayed in $T$
and provide information of the topology of $S$. These informative triples are
used to infer the unknown species tree $S$ for $T$. We obtain a similar result
for non-binary gene trees. To this end, however, the reconciliation map needs
to be further restricted. We provide a polynomial-time algorithm to decide
whether there is a species tree for a given event-labeled gene tree, and in the
positive case, to construct the species tree and the respective (restricted)
reconciliation map.
  However, informative triples as well as DTL-scenarios have its limitations
when they are used to explain the biological feasibility of gene trees. While
reconciliation maps imply biological feasibility, we show that the converse is
not true in general. Moreover, we show that informative triples do neither
provide enough information to characterize "relaxed" DTL-scenarios nor
non-restricted reconciliation maps for non-binary biologically feasible gene
trees.


  We deal with the study of the evolution of the allelic frequencies, at a
single locus, for a population distributed continuously over a bounded habitat.
We consider evolution which occurs under the joint action of selection and
arbitrary migration, that is independent of genotype, in absence of mutation
and random drift. The focus is on a conjecture, that was raised up in
literature of population genetics, about the possible uniqueness of polymorphic
equilibria, which are known as clines, under particular circumstances. We study
the number of these equilibria, making use of topological tools, and we give a
negative answer to that question by means of two examples. Indeed, we provide
numerical evidence of multiplicity of positive solutions for two different
Neumann problems satisfying the requests of the conjecture.


  We introduce a low dimensional function of the site frequency spectrum that
is tailor-made for distinguishing coalescent models with multiple mergers from
Kingman coalescent models with population growth, and use this function to
construct a hypothesis test between these model classes. The null and
alternative sampling distributions of the statistic are intractable, but its
low dimensionality renders them amenable to Monte Carlo estimation. We
construct kernel density estimates of the sampling distributions based on
simulated data, and show that the resulting hypothesis test dramatically
improves on the statistical power of a current state-of-the-art method. A key
reason for this improvement is the use of multi-locus data, in particular
averaging observed site frequency spectra across unlinked loci to reduce
sampling variance. We also demonstrate the robustness of our method to nuisance
and tuning parameters. Finally we show that the same kernel density estimates
can be used to conduct parameter estimation, and argue that our method is
readily generalisable for applications in model selection, parameter inference
and experimental design.


  News spread in internet media outlets can be seen as a contagious process
generating temporal networks representing the influence between published
articles. In this article we propose a methodology based on the application of
natural language analysis of the articles to reconstruct the spread network.
From the reconstructed network, we show that the dynamics of the news spread
can be approximated by a classical SIR epidemiological dynamics upon the
network. From the results obtained we argue that the methodology proposed can
be used to make predictions about media repercussion, and also to detect viral
memes in news streams.


  We study macroevolutionary dynamics by extending microevolutionary
competition models to long time scales. It has been shown that for a general
class of competition models, gradual evolutionary change in continuous
phenotypes (evolutionary dynamics) can be non-stationary and even chaotic when
the dimension of the phenotype space in which the evolutionary dynamics unfold
is high. It has also been shown that evolutionary diversification can occur
along non-equilibrium trajectories in phenotype space. We combine these lines
of thinking by studying long-term coevolutionary dynamics of emerging lineages
in multi-dimensional phenotype spaces. We use a statistical approach to
investigate the evolutionary dynamics of many different systems. We find: 1)
for a given dimension of phenotype space, the coevolutionary dynamics tends to
be fast and non-stationary for an intermediate number of coexisting lineages,
but tends to stabilize as the evolving communities reach a saturation level of
diversity; and 2) the amount of diversity at the saturation level increases
rapidly (exponentially) with the dimension of phenotype space. These results
have implications for theoretical perspectives on major macroevolutionary
patterns such as adaptive radiation, long-term temporal patterns of phenotypic
changes, and the evolution of diversity.


  Explicit accounting for copy number alterations can dramatically improve
mutation frequency estimates, leading to more accurate phylogeny
reconstructions and subclone characterizations.


  Summary: Counting all k-mers in a given dataset is a standard procedure in
many bioinformatics applications. We introduce KMC3, a significant improvement
of the former KMC2 algorithm together with KMC tools for manipulating k-mer
databases. Usefulness of the tools is shown on a few real problems.
Availability: Program is freely available at
http://sun.aei.polsl.pl/REFRESH/kmc. Contact: sebastian.deorowicz@polsl.pl


  Genome replication, a key process for a cell, relies on stochastic initiation
by replication origins, causing a variability of replication timing from cell
to cell. While stochastic models of eukaryotic replication are widely
available, the link between the key parameters and overall replication timing
has not been addressed systematically.We use a combined analytical and
computational approach to calculate how positions and strength of many origins
lead to a given cell-to-cell variability of total duration of the replication
of a large region, a chromosome or the entire genome.Specifically, the total
replication timing can be framed as an extreme-value problem, since it is due
to the last region that replicates in each cell. Our calculations identify two
regimes based on the spread between characteristic completion times of all
inter-origin regions of a genome. For widely different completion times, timing
is set by the single specific region that is typically the last to replicate in
all cells. Conversely, when the completion time of all regions are
comparable,an extreme-value estimate shows that the cell-to-cell variability of
genome replication timing has universal properties. Comparison with available
data shows that the replication program of three yeast species falls in this
extreme-value regime.


  Salmon farming has become a prosperous international industry over the last
decades. Along with growth in the production farmed salmon, however, an
increasing threat by pathogens has emerged. Of special concern is the
propagation and spread of the salmon louse, Lepeophtheirus salmonis. In order
to gain insight into this parasites population dynamics in large scale salmon
farming system, we present a fully mechanistic stage-structured population
model for the salmon louse, also allowing for complexities involved in the
hierarchical structure of full scale salmon farming. The model estimates
parameters controlling a wide range of processes, including temperature
dependent demographic rates, fish size and abundance effects on louse
transmission rates, effects sizes of various salmon louse control measures, and
distance based between farm transmission rates. Model parameters were estimated
from data including 32 salmon farms, except the last production months for five
farms which were used to evaluate model predictions. We used a Bayesian
estimation approach, combining the prior distributions and the data likelihood
into a joint posterior distribution for all model parameters. The model
generated expected values that fitted the observed infection levels of the
chalimus, adult female and other mobile stages of salmon lice, reasonably well.
Predictions for the time periods not used for fitting the model were also
consistent with the observational data. We argue that the present model for the
population dynamics of the salmon louse in aquaculture farm systems may
contribute to resolve the complexity of processes that drive that drive this
host-parasite relationship, and hence may improve strategies to control the
parasite in this production system.


  Identifying protein functional sites (PFSs) and protein-ligand interactions
(PLIs) are critically important in understanding the protein function and the
involved biochemical reactions. As large amount of unknown proteins are quickly
accumulated in this post-genome era, an urgent task arises to predict PFSs and
PLIs at residual level. Nowadays many knowledge-based methods have been well
developed for prediction of PFSs, however, accurate methods for PLI prediction
are still lacking. In this study, we have presented a new method for prediction
of PLIs and PFSs based on sequence of the inquiry protein. The key of the
method hinges on a function- and interaction-annotated protein domain profile
database, called fiDPD, which was built from the Structural Classification of
Proteins (SCOP) database, using a hidden Markov model program. The method was
applied to 13 target proteins from the recent Critical Assessment of Structure
Prediction (CASP10/11). Our calculations gave a Matthews correlation
coefficient (MCC) value of 0.66 for prediction of PFSs, and an 80% recall in
prediction of PLIs. Our method reveals that PLIs are conserved during the
evolution of proteins, and they can be reliably predicted from fiDPD. fiDPD can
be used as a complement to existent bioinformatics tools for protein function
annotation.


  Nonnative residual interactions have attracted increasing attention in recent
protein folding researches. Experimental and theoretical investigations had
been set out to catch nonnative contacts that might dominate key events in
protein folding. However, energetic frustrations caused by nonnative
inter-residue interactions are not systematically characterized, due to the
complicated folding mechanism. Recently, we studied the folding of a set of
homologous all-{\alpha} proteins and found that nonnative-contact-based
energetic frustrations are closely related to protein native-contact networks.
In this paper, we studied the folding of nine homologous immunoglobulin-like
(Ig-like) \b{eta}-sandwich proteins and examined energetic frustrations caused
by nonnative inter-residue interactions, based on analyses of residual
phi-values and contact maps of transition state ensembles. The proteins share
similar tertiary structures, thus minimize topology frustration differences in
protein folding and highlighting the effects of energetic frustrations as
caused by hydrophilic-hydrophobic mutations. Our calculations showed that
energetic frustrations have highly heterogeneous effects on the folding of most
secondary structures and on the folding correlations between different
folding-patches of \b{eta}-sandwich proteins. The simulations revealed a strong
interplay between energetic frustrations and native-contact networks in
\b{eta}-sandwich domains, which ensures that \b{eta}-sandwich domains follow a
common folding mechanism. Our results suggested that the folding processes of
\b{eta}-sandwich proteins might be redesigned by carefully manipulating
energetic frustrations at residue level.


  As high-throughput biological sequencing becomes faster and cheaper, the need
to extract useful information from sequencing becomes ever more paramount,
often limited by low-throughput experimental characterizations. For proteins,
accurate prediction of their functions directly from their primary amino-acid
sequences has been a long standing challenge. Here, machine learning using
artificial recurrent neural networks (RNN) was applied towards classification
of protein function directly from primary sequence without sequence alignment,
heuristic scoring or feature engineering. The RNN models containing
long-short-term-memory (LSTM) units trained on public, annotated datasets from
UniProt achieved high performance for in-class prediction of four important
protein functions tested, particularly compared to other machine learning
algorithms using sequence-derived protein features. RNN models were used also
for out-of-class predictions of phylogenetically distinct protein families with
similar functions, including proteins of the CRISPR-associated nuclease,
ferritin-like iron storage and cytochrome P450 families. Applying the trained
RNN models on the partially unannotated UniRef100 database predicted not only
candidates validated by existing annotations but also currently unannotated
sequences. Some RNN predictions for the ferritin-like iron sequestering
function were experimentally validated, even though their sequences differ
significantly from known, characterized proteins and from each other and cannot
be easily predicted using popular bioinformatics methods. As sequencing and
experimental characterization data increases rapidly, the machine-learning
approach based on RNN could be useful for discovery and prediction of
homologues for a wide range of protein functions.


  Positive density-dependence occurs when individuals experience increased
survivorship, growth, or reproduction with increased population densities.
Mechanisms leading to these positive relationships include mate limitation,
saturating predation risk, and cooperative breeding and foraging. Individuals
within these populations may differ in age, size, or geographic location and
thereby structure these populations. Here, I study structured population models
accounting for positive density-dependence and environmental stochasticity i.e.
random fluctuations in the demographic rates of the population. Under an
accessibility assumption (roughly, stochastic fluctuations can lead to
populations getting small and large), these models are shown to exhibit a
dynamical trichotomy: (i) for all initial conditions, the population goes
asymptotically extinct with probability one, (ii) for all positive initial
conditions, the population persists and asymptotically exhibits unbounded
growth, and (iii) for all positive initial conditions, there is a positive
probability of asymptotic extinction and a complementary positive probability
of unbounded growth. The main results are illustrated with applications to
spatially structured populations with an Allee effect and age-structured
populations experiencing mate limitation.


  Free energy perturbation (FEP) is frequently used to evaluate the free energy
change of a biological process, e.g. the drug binding free energy or the ligand
solvation free energy. Due to the sampling inefficiency, FEP is often employed
together with computationally expensive enhanced sampling methods. Here we show
that this sampling inefficiency, which stems from not accounting for the
environmental reorganization, is an inherent property of the single-ensemble
ansatz of FEP, and hence simply prolonging the MD simulation can hardly
alleviate the problem. Instead, we propose a new, multi-ensemble ansatz -- the
multi-layer free energy perturbation (MLFEP), which allows environmental
reorganization processes (relaxation) to occur automatically during the MD
sampling. Our study paves the way toward a fast but accurate free energy
calculation that can be employed in computer-aided drug design.


  Astrocytes affect neural transmission by a tight control via glutamate
transporters on glutamate concentrations in direct vicinity to the synaptic
cleft and by extracellular glutamate. Their relevance for information
representation has been supported by in-vivo studies in ferret and mouse
primary visual cortex. In ferret blocking glutamate transport pharmacologically
broadened tuning curves and enhanced the response at preferred orientation. In
knock-out mice with reduced expression of glutamate transporters sharpened
tuning was observed. It is however unclear how focal and ambient changes in
glutamate concentration affect stimulus representation. Here we develop a
computational framework, which allows the investigation of synaptic and
extrasynaptic effects of glutamate uptake on orientation tuning in recurrently
connected network models with pinwheel-domain (ferret) or salt-and-pepper
(mouse) organization. This model proposed that glutamate uptake shapes
information representation when it affects the contribution of excitatory and
inhibitory neurons to the network activity. Namely, strengthening the
contribution of excitatory neurons generally broadens tuning and elevates the
response. In contrast, strengthening the contribution of inhibitory neurons can
have a sharpening effect on tuning. In addition local representational topology
also plays a role: In the pinwheel-domain model effects were strongest within
domains - regions where neighboring neurons share preferred orientations.
Around pinwheels but also within salt-and-pepper networks the effects were less
strong. Our model proposes that the pharmacological intervention in ferret
increases the contribution of excitatory cells, while the reduced expression in
mouse increases the contribution of inhibitory cells to network activity.


  Networks can model real-world systems in a variety of domains. Network
alignment (NA) aims to find a node mapping that conserves similar regions
between compared networks. NA is applicable to many fields, including
computational biology, where NA can guide the transfer of biological knowledge
from well- to poorly-studied species across aligned network regions. Existing
NA methods can only align static networks. However, most complex real-world
systems evolve over time and should thus be modeled as dynamic networks. We
hypothesize that aligning dynamic network representations of evolving systems
will produce superior alignments compared to aligning the systems' static
network representations, as is currently done. For this purpose, we introduce
the first ever dynamic NA method, DynaMAGNA++. This proof-of-concept dynamic NA
method is an extension of a state-of-the-art static NA method, MAGNA++. Even
though both MAGNA++ and DynaMAGNA++ optimize edge as well as node conservation
across the aligned networks, MAGNA++ conserves static edges and similarity
between static node neighborhoods, while DynaMAGNA++ conserves dynamic edges
(events) and similarity between evolving node neighborhoods. For this purpose,
we introduce the first ever measure of dynamic edge conservation and rely on
our recent measure of dynamic node conservation. Importantly, the two dynamic
conservation measures can be optimized using any state-of-the-art NA method and
not just MAGNA++. We confirm our hypothesis that dynamic NA is superior to
static NA, under fair comparison conditions, on synthetic and real-world
networks, in computational biology and social network domains. DynaMAGNA++ is
parallelized and it includes a user-friendly graphical interface.


  An important problem in evolutionary biology is to reconstruct the
evolutionary history of a set $X$ of species. This history is often represented
as a phylogenetic network, that is, a connected graph with leaves labelled by
elements in $X$ (for example, an evolutionary tree), which is usually also
binary, i.e. all vertices have degree 1 or 3. A common approach used in
phylogenetics to build a phylogenetic network on $X$ involves constructing it
from networks on subsets of $X$. Here we consider the question of which
(unrooted) phylogenetic networks are leaf-reconstructible, i.e. which networks
can be uniquely reconstructed from the set of networks obtained from it by
deleting a single leaf (its $X$-deck). This problem is closely related to the
(in)famous reconstruction conjecture in graph theory but, as we shall show,
presents distinct challenges. We show that some large classes of phylogenetic
networks are reconstructible from their $X$-deck. This includes phylogenetic
trees, binary networks containing at least one non-trivial cut-edge, and binary
level-4 networks (the level of a network measures how far it is from being a
tree). We also show that for fixed $k$, almost all binary level-$k$
phylogenetic networks are leaf-reconstructible. As an application of our
results, we show that a level-3 network $N$ can be reconstructed from its
quarnets, that is, 4-leaved networks that are induced by $N$ in a certain
recursive fashion. Our results lead to several interesting open problems which
we discuss, including the conjecture that all phylogenetic networks with at
least five leaves are leaf-reconstructible.


  Phylogenetic networks are a generalization of evolutionary trees that are
used by biologists to represent the evolution of organisms which have undergone
reticulate evolution. Essentially, a phylogenetic network is a directed acyclic
graph having a unique root in which the leaves are labelled by a given set of
species. Recently, some approaches have been developed to construct
phylogenetic networks from collections of networks on 2- and 3-leaved networks,
which are known as binets and trinets, respectively. Here we study in more
depth properties of collections of binets, one of the simplest possible types
of networks into which a phylogenetic network can be decomposed. More
specifically, we show that if a collection of level-1 binets is compatible with
some binary network, then it is also compatible with a binary level-1 network.
Our proofs are based on useful structural results concerning lowest stable
ancestors in networks. In addition, we show that, although the binets do not
determine the topology of the network, they do determine the number of
reticulations in the network, which is one of its most important parameters. We
also consider algorithmic questions concerning binets. We show that deciding
whether an arbitrary set of binets is compatible with some network is at least
as hard as the well-known Graph Isomorphism problem. However, if we restrict to
level-1 binets, it is possible to decide in polynomial time whether there
exists a binary network that displays all the binets. We also show that to find
a network that displays a maximum number of the binets is NP-hard, but that
there exists a simple polynomial-time 1/3-approximation algorithm for this
problem. It is hoped that these results will eventually assist in the
development of new methods for constructing phylogenetic networks from
collections of smaller networks.


  The major biochemical networks of the living cell, the network of interacting
genes and the network of biochemical reactions, are highly interdependent,
however, they have been studied mostly as separate systems so far. In the last
years an appropriate theoretical framework for studying interdependent networks
has been developed in the context of statistical physics. Here we study the
interdependent network of gene regulation and metabolism of the model organism
Escherichia coli using the theoretical framework of interdependent networks. In
particular we aim at understanding how the biological system can consolidate
the conflicting tasks of reacting rapidly to (internal and external)
perturbations, while being robust to minor environmental fluctuations, at the
same time. For this purpose we study the network response to localized
perturbations and find that the interdependent network is sensitive to gene
regulatory and protein-level perturbations, yet robust against metabolic
changes. This first quantitative application of the theory of interdependent
networks to systems biology shows how studying network responses to localized
perturbations can serve as a useful strategy for analyzing a wide range of
other interdependent networks.


  Genetic oscillators are present in the cells of many organisms and control
several biological processes. The common feature of such oscillators is the
presence of a protein which represses the transcription of its own gene.
Recently, it has been shown that for many genes transcription is not a
continuous process, but that it proceeds in bursts. We study here the
relationship between bursty transcription and the robustness of protein
oscillations. We concentrate on the temporal profile of mRNA production by
studying regimes where this profile changes but the amount of mRNA produced is
kept fixed. For systems with different degrees of cooperativity we show that in
general bursts are associated with more robust oscillations, but when they are
too short and intense they can have the opposite effect. In other words, we
show that, in terms of the regularity of the oscillations generated, there is
an optimal value for the intensity of the bursts.


  Computational techniques are required for narrowing down the vast space of
possibilities to plausible prebiotic scenarios, since precise information on
the molecular composition, the dominant reaction chemistry, and the conditions
for that era are scarce. The exploration of large chemical reaction networks is
a central aspect in this endeavour. While quantum chemical methods can
accurately predict the structures and reactivities of small molecules, they are
not efficient enough to cope with large-scale reaction systems. The
formalization of chemical reactions as graph grammars provides a generative
system, well grounded in category theory, at the right level of abstraction for
the analysis of large and complex reaction networks. An extension of the basic
formalism into the realm of integer hyperflows allows for the identification of
complex reaction patterns, such as auto-catalysis, in large reaction networks
using optimization techniques.


  Background: MicroRNAs (miRNAs) play multiple roles in tumor biology [1].
Interestingly, reports from multiple groups suggest that miRNA targets may be
coupled through competitive stoichiometric sequestration [2]. Specifically,
computational models predicted [3, 4] and experimental assays confirmed [5, 6]
that miRNA activity is dependent on miRNA target abundance, and consequently,
changes to the abundance of some miRNA targets lead to changes to the
regulation and abundance of their other targets. The resulting indirect
regulatory influence between miRNA targets resembles competition and has been
dubbed competitive endogenous RNA (ceRNA) [5, 7, 8]. Recent studies have
questioned the physiological relevance of ceRNA interactions [9], researchers
ability to accurately predict these interactions [10], and the number of genes
that are impacted by ceRNA interactions in specific cellular contexts [11].
Results: To address these concerns, we reverse engineered ceRNA networks
(ceRNETs) in breast and prostate adenocarcinomas using context-specific TCGA
profiles [12-14], and tested whether ceRNA interactions can predict the effects
of RNAi-mediated gene silencing perturbations in PC3 and MCF7 cells. Our
results, based on tests of thousands of inferred ceRNA interactions that are
predicted to alter hundreds of cancer genes in each of the two tumor contexts,
confirmed statistically significant effects for half of the predicted targets.
Conclusions: Our results suggest that the expression of a significant fraction
of cancer genes may be regulated by ceRNA interactions in each of the two tumor
contexts.


  We introduce a new tuberculosis (TB) mathematical model, with $25$
state-space variables where $15$ are evolution disease states (EDSs), which
generalises previous models and takes into account the (seasonal) flux of
populations between a high incidence TB country (A) and a host country (B) with
low TB incidence, where (B) is divided into a community (G) with high
percentage of people from (A) plus the rest of the population (C). Contrary to
some beliefs, related to the fact that agglomerations of individuals increase
proportionally to the disease spread, analysis of the model shows that the
existence of semi-closed communities are beneficial for the TB control from a
global viewpoint. The model and techniques proposed are applied to a case-study
with concrete parameters, which model the situation of Angola (A) and Portugal
(B), in order to show its relevance and meaningfulness. Simulations show that
variations of the transmission coefficient on the origin country has a big
influence on the number of infected (and infectious) individuals on the
community and the host country. Moreover, there is an optimal ratio for the
distribution of individuals in (C) versus (G), which minimizes the reproduction
number $R_0$. Such value does not give the minimal total number of infected
individuals in all (B), since such is attained when the community (G) is
completely isolated (theoretical scenario). Sensitivity analysis and curve
fitting on $R_0$ and on EDSs are pursuit in order to understand the TB effects
in the global statistics, by measuring the variability of the relevant
parameters. We also show that the TB transmission rate $\beta$ does not act
linearly on $R_0$, as is common in compartment models where system feedback or
group interactions do not occur. Further, we find the most important parameters
for the increase of each EDS.


  In this paper, we show that many structured epidemic models may be described
using a straightforward product structure. Such products, derived from products
of directed graphs, may represent useful refinements including geographic and
demographic structure, age structure, gender, risk groups, or immunity status.
Extension to multi-strain dynamics, i.e. pathogen heterogeneity, is also shown
to be feasible in this framework. Systematic use of such products may aid in
model development and exploration, can yield insight, and could form the basis
of a systematic approach to numerical structural sensitivity analysis.


  Many biological populations, such as bacterial colonies, have developed
through evolution a protection mechanism, called bet-hedging, to increase their
probability of survival under stressful environmental fluctutation. In this
context, the concept of preadaptation refers to a common type of bet-hedging
protection strategy in which a relatively small number of individuals in a
population stochastically switch their phenotypes to a `dormant' metabolic
state in which they increase their probability of survival against potential
environmental shocks. Hence, if an environmental shock took place at some point
in time, preadapted organisms would be better adapted to survive and
proliferate once the shock is over. In many biological populations, the
mechanisms of preadaptation and proliferation present delays whose influence in
the fitness of the population are not well-understood. In this paper, we
propose a rigorous mathematical framework to analyze the role of delays in both
preadaptation and proliferation mechanisms in the survival of biological
populations, with an emphasis on bacterial colonies. Our theoretical framework
allows us to analytically quantify the average growth rate of a bet-hedging
bacterial colony with stochastically delayed reactions with arbitrary
precision. We verify the accuracy of the proposed method by numerical
simulations and conclude that the growth rate of a bet-hedging population shows
a non-trivial dependency on their preadaptation and proliferation delays.
Contrary to the current belief, our results show that faster reactions do not,
in general, increase the overall fitness of a biological population.


  The percentage and sequence of AT and GC base pairs and charges on the DNA
backbone contribute significantly to the stiffness of DNA. This elastic
property of DNA also changes with small interacting ligands. The
single-molecule force spectroscopy technique shows different interaction modes
by measuring the mechanical properties of DNA bound with small ligands. When a
ds-DNA molecule is overstretched in the presence of ligands, it undergoes a
co-operative structural transition based on the externally applied force, the
mode of binding of the ligands, the binding constant of the ligands to the DNA,
the concentration of the ligands and the ionic strength of the supporting
medium. This leads to the changes in the regions- upto 60 pN, cooperative
structural transition region and the overstretched region, compared to that of
the FEC in the absence of any binding ligand. The cooperative structural
transitions were studied by the extended and twistable worm-like chain model.
Here we have depicted these changes in persistence length and the elastic
modulus constant as a function of binding constant and the concentration of the
bound ligands, which vary with time. Therefore, besides ionic strength,
interacting proteins and content of AT and GC base pairs, the ligand binding or
intercalation with the ligands is an important parameter which changes the
stiffness of DNA.


  The sampling of scale-free networks in Molecular Biology is usually achieved
by growing networks from a seed using recursive algorithms with elementary
moves which include the addition and deletion of nodes and bonds. These
algorithms include the Barabasi-Albert algorithm. Later algorithms, such as the
Duplication-Divergence algorithm, the Sol\'e algorithm and the iSite algorithm,
were inspired by biological processes underlying the evolution of protein
networks, and the networks they produce differ essentially from networks grown
by the Barabasi-Albert algorithm. In this paper the mean field analysis of
these algorithms is reconsidered, and extended to variant and modified
implementations of the algorithms. The degree sequences of scale-free networks
decay according to a powerlaw distribution, namely $P(k) \sim k^{-\gamma}$,
where $\gamma$ is a scaling exponent. We derive mean field expressions for
$\gamma$, and test these by numerical simulations. Generally, good agreement is
obtained. We also found that some algorithms do not produce scale-free networks
(for example some variant Barabasi-Albert and Sol\'e networks).


  The past decade has seen a rapid growth in omics technologies. Genome-wide
association studies (GWAS) have uncovered susceptibility variants for a variety
of complex traits. However, the functional significance of most discovered
variants are still not fully understood. On the other hand, there is increasing
interest in exploring the role of epigenetic variations such as DNA methylation
in disease pathogenesis. In this work, we present a general framework for
epigenome-wide association study and integrative analysis with the
transcriptome based on GWAS summary statistics and data from methylation and
expression quantitative trait loci (QTL) studies. The framework is based on
Mendelian randomization, which is much less vulnerable to confounding and
reverse causation compared to conventional studies. The framework was applied
to five complex diseases. We first identified loci that are differentially
methylated due to genetic variations, and then developed several approaches for
joint testing with the GWAS-imputed transcriptome. We discovered a number of
novel candidate genes that are not implicated in the original GWAS studies. We
also observed strong evidence (lowest p = 2.01e-184) for differential
expression among the top genes mapped to methylation loci. The framework
proposed here opens a new way of analyzing GWAS summary data and will be useful
for gaining deeper insight into disease mechanisms.


  An essential quantity to ensure evolvability of populations is the
navigability of the genotype space. Navigability relies on the existence of
sufficiently large genotype networks, that is ensembles of sequences with the
same phenotype that guarantee an efficient random drift through sequence space.
The number of sequences compatible with a given structure (e.g. the number of
RNA sequences folding into a particular secondary structure, or the number of
DNA sequences coding for the same protein structure) is astronomically large in
all functional molecules investigated. However, an exhaustive experimental or
computational study of all RNA folds or all protein structures becomes
impossible even for moderately long sequences. Here, we analytically derive the
distribution of phenotype sizes for a hierarchy of models which successively
incorporate features of increasingly realistic sequence-to-structure
genotype-phenotype maps. The main feature of these models relies on the
characterization of each phenotype through a prototypical sequence whose sites
admit a variable fraction of letters of the alphabet. Our models interpolate
between two limit distributions: a powerlaw distribution, when the ordering of
sites in the prototypical sequence is strongly constrained, and a lognormal
distribution, as suggested for RNA, when different orderings of the same set of
sites yield different phenotypes. Our main result is the qualitative and
quantitative identification of those features of the sequence-to-structure map
that lead to different distributions of phenotype sizes.


  Improved understanding of molecular systems has only emphasised the
sophistication of networks within the cell. Simultaneously, the advance of
nucleic acid nanotechnology, a platform within which reactions can be
exquisitely controlled, has made the development of artificial architectures
and devices possible. Vital to this progress has been a solid foundation in the
thermodynamics of molecular systems. In this pedagogical review and
perspective, I will discuss how thermodynamics determines both the overall
potential of molecular networks, and the minute details of design. I will then
argue that, in turn, the need to understand molecular systems is helping to
drive the development of theories of thermodynamics at the microscopic scale.


  Optical biosensors are often used to measure kinetic rate constants
associated with chemical reactions. Such instruments operate in the
\textit{surface-volume} configuration, in which ligand molecules are convected
through a fluid-filled volume over a surface to which receptors are confined.
Currently, scientists are using optical biosenors to measure the kinetic rate
constants associated with DNA translesion synthesis--a process critical to DNA
damage repair. Biosensor experiments to study this process involve multiple
interacting components on the sensor surface. This multiple-component biosensor
experiment is modeled with a set of nonlinear integrodifferential equations
(IDEs). It is shown that in physically relevant asymptotic limits these
equations reduce to a much simpler set of Ordinary Differential Equations
(ODEs). To verify the validity of our ODE approximation, a numerical method for
the IDE system is developed and studied. Results from the ODE model agree with
simulations of the IDE model, rendering our ODE model useful for parameter
estimation.


  Cytochrome C oxidase and FoF1-ATP synthase constitute complex IV and V,
respectively, of the five membrane-bound enzymes in mitochondria comprising the
respiratory chain. These enzymes are located in the inner mitochondrial
membrane (IMM), which exhibits large invaginations called cristae. According to
recent cryo-tomography, FoF1-ATP synthases are located predominantly at the rim
of the cristae, while cytochrome C oxidases are likely distributed in planar
membrane areas of the cristae. Previous FLIM measurements (K. Busch and
coworkers) of complex II and III unravelled differences in the local
environment of the membrane enzymes in the cristae. Here, we tagged complex IV
and V with mNeonGreen and investigated their mitochondrial nano-environment by
FLIM and superresolution microscopy in living human cells. Different lifetimes
and anisotropy values were found and will be discussed.


  Protein synthesis rates are determined, at the translational level, by
properties of the transcript's sequence. The efficiency of an mRNA can be tuned
by varying the ribosome binding sites controlling the recruitment of the
ribosomes, or the codon usage establishing the speed of protein elongation. In
this work we propose transcript length as a further key determinant of
translation efficiency. Based on a physical model that considers the kinetics
of ribosomes advancing on the mRNA and diffusing in its surrounding, as well as
mRNA circularisation and ribosome drop-off, we explain how the transcript
length may play a central role in establishing ribosome recruitment and the
overall translation rate of an mRNA. According to our results, the proximity of
the 3' end to the ribosomal recruitment site of the mRNA could induce a
feedback in the translation process that would favour the recycling of
ribosomes. We also demonstrate how this process may be involved in shaping the
experimental ribosome density-gene length dependence. Finally, we argue that
cells could exploit this mechanism to adjust and balance the usage of its
ribosomal resources.


  Motivation: Epigenetic heterogeneity within a tumour can play an important
role in tumour evolution and the emergence of resistance to treatment. It is
increasingly recognised that the study of DNA methylation (DNAm) patterns along
the genome -- so-called `epialleles' -- offers greater insight into epigenetic
dynamics than conventional analyses which examine DNAm marks individually.
  Results: We have developed a Bayesian model to infer which epialleles are
present in multiple regions of the same tumour. We apply our method to reduced
representation bisulfite sequencing (RRBS) data from multiple regions of one
lung cancer tumour and a matched normal sample. The model borrows information
from all tumour regions to leverage greater statistical power. The total number
of epialleles, the epiallele DNAm patterns, and a noise hyperparameter are all
automatically inferred from the data. Uncertainty as to which epiallele an
observed sequencing read originated from is explicitly incorporated by
marginalising over the appropriate posterior densities. The degree to which
tumour samples are contaminated with normal tissue can be estimated and
corrected for. By tracing the distribution of epialleles throughout the tumour
we can infer the phylogenetic history of the tumour, identify epialleles that
differ between normal and cancer tissue, and define a measure of global
epigenetic disorder.


  Chemical reaction networks with generalized mass-action kinetics lead to
power-law dynamical systems. As a simple example, we consider the Lotka
reactions and the resulting planar ODE. We characterize the parameters
(positive coefficients and real exponents) for which the unique positive
equilibrium is a center.


  Recent experiments have shown that trivalent ion, spermidine$^{3+}$, can
provoke lateral microphase segregation in DNA brushes. Using molecular
simulations and simple theoretical arguments, we explore the effects of
trivalent counterions on polyelectrolyte brushes. At a proper range of grafting
density, polymer size, and ion concentration, the brush polymers collapse
heterogeneously into octopus-like surface micelles. Remarkably, the
heterogeneity in brush morphology is maximized and the relaxation dynamics of
chain and condensed ion are the slowest at the 1:3 stoichiometric concentration
of trivalent ions to polyelectrolyte charge. A further increase of trivalent
ion concentration conducive to a charge inversion elicits modest reswelling and
homogenizes the morphology of brush condensate. Our study provides a new
insight into the origin of the diversity in DNA organization in cell nuclei as
well as the ion-dependent morphological variation in polyelectrolyte brush
layer of biological membranes.


  We study a stochastic model of infection spreading on a network. At each time
step a node is chosen at random, along with one of its neighbors. If the node
is infected and the neighbor is susceptible, the neighbor becomes infected. How
many time steps $T$ does it take to completely infect a network of $N$ nodes,
starting from a single infected node? An analogy to the classic "coupon
collector" problem of probability theory reveals that the takeover time $T$ is
dominated by extremal behavior, either when there are only a few infected nodes
near the start of the process or a few susceptible nodes near the end. We show
that for $N \gg 1$, the takeover time $T$ is distributed as a Gumbel for the
star graph; as the sum of two Gumbels for a complete graph and an
Erd\H{o}s-R\'{e}nyi random graph; as a normal for a one-dimensional ring and a
two-dimensional lattice; and as a family of intermediate skewed distributions
for $d$-dimensional lattices with $d \ge 3$ (these distributions approach the
sum of two Gumbels as $d$ approaches infinity). Connections to evolutionary
dynamics, cancer, incubation periods of infectious diseases, first-passage
percolation, and other spreading phenomena in biology and physics are
discussed.


  Ribonucleic acid (RNA) is involved in many regulatory and catalytic processes
in the cell. The function of any RNA molecule is intimately related with its
structure. In-line probing experiments provide valuable structural datasets for
a variety of RNAs and are used to characterize conformational changes in
riboswitches. However, the structural determinants that lead to differential
reactivities in unpaired nucleotides have not been investigated yet. In this
work we used a combination of theoretical approaches, i.e., classical molecular
dynamics simulations, multiscale quantum mechanical/molecular mechanical
calculations, and enhanced sampling techniques in order to compute and
interpret the differential reactivity of individual residues in several RNA
motifs including members of the most important GNRA and UNCG tetraloop
families. Simulations on the multi ns timescale are required to converge the
related free-energy landscapes. The results for uGAAAg and cUUCGg tetraloops
and double helices are compared with available data from in-line probing
experiments and show that the introduced technique is able to distinguish
between nucleotides of the uGAAAg tetraloop based on their structural
predispositions towards phosphodiester backbone cleavage. For the cUUCGg
tetraloop, more advanced ab initio calculations would be required. This study
is the first attempt to computationally classify chemical probing experiments
and paves the way for an identification of tertiary structures based on the
measured reactivity of non-reactive nucleotides.


  Regional tuna fishery management organizations cannot provide specific advice
to local fishery managers in small island jurisdictions. The State of Hawaii
maintains time series of yellowfin tuna catches dating back to 1949, but these
data have never been formally applied to evaluating the effects of the
yellowfin fishery in the Main Hawaiian Islands on the local stock. I develop a
new approach utilizing these data that links the local stock dynamics to the
dynamics of the larger Pacific stock. This approach uses a state-space logistic
production model linked to the larger Pacific stock using an index of
abundance. The conclusion is that such a model is feasible, that the local
stock is not overfished and that local fisheries are fishing at acceptable
levels.


  Emergence of antibiotic resistance due to New Delhi Metallo $\beta$-lactamase
(NDM-1) bacterial enzymes is of great concern due to their ability to hydrolyze
wide range of antibiotics. Efforts are ongoing to obtain the atomistic details
of the hydrolysis mechanism in order to develop novel drugs and inhibitors
against NDM-1. Especially, it remains elusive how drug molecules of different
family of antibiotics are hydrolyzed by NDM-1 in an efficient manner. Here we
report the detailed molecular mechanism of NDM-1 catalyzed hydrolysis of
cephalexin, a cephalosporin family drug, and meropenem, a carbapenem family
drug. This study employs molecular dynamics (MD) simulations using hybrid
quantum mechanical/molecular mechanical (QM/MM) methods at the density
functional theory level, based on which reaction pathways and the associated
free energies are obtained. We find that the mechanism and the free energy
barrier for the ring-opening step are the same for both the drug molecules,
while the subsequent protonation step differs. In particular, we observe that
the mechanism of the protonation step depends on the R2 group of the drug
molecule. Our simulations show that allylic carbon protonation occurs in the
case of cephalexin drug molecule where Lys211 is the proton donor and the
proton transfer occurs via a water chain formed (only) at the ring-opened
intermediate structure. Based on the free energy profiles, the overall kinetics
of the drug hydrolysis is discussed. Finally, we show that the proposed
mechanisms and free energy profiles could explain various experimental
observations.


  During the asymmetric division of the Caenorhabditis elegans nematode zygote,
the polarity cues distribution and daughter cell fates depend on the correct
positioning of the mitotic spindle, which results from both centering and
cortical pulling forces. Revealed by anaphase spindle rocking, these pulling
forces are regulated by the force generator dynamics, which are in turn
consequent of mitotic progression. We found a novel, additional, regulation of
these forces by the spindle position. It controls astral microtubule
availability at the cortex, on which the active force generators can pull.
Importantly, this positional control relies on the polarity dependent LET-99
cortical band, which restricts or concentrates generators to a posterior
crescent. We ascribed this control to the microtubule dynamics at the cortex.
Indeed, in mapping the cortical contacts, we found a correlation between the
centrosome-cortex distance and the microtubule contact density. In turn, it
modulates pulling force generator activity. We modelled this control,
predicting and experimentally validating that the posterior crescent extent
controlled where the anaphase oscillations started, in addition to mitotic
progression. Finally, we propose that spatially restricting force generator to
a posterior crescent sets the spindle's final position, reflecting polarity
through the LET-99 dependent restriction of force generators to a posterior
crescent. This regulation superimposes that of force generator processivity.
This novel control confers a low dependence on microtubule and active force
generator exact numbers or dynamics, provided that they exceed the threshold
needed for posterior displacement. Interestingly, this robustness originates in
cell mechanics rather than biochemical networks.


  The study of Chemical Reaction Networks (CRN's) is a very active field.
Earlier well-known results \cite{Feinberg:def_01:87, Anderson:product_dist:10}
identify a topological quantity called deficiency, for any CRN, which, when
exactly equal to zero, leads to a unique factorized steady-state for these
networks. No results exist however for the steady states of non-zero-deficiency
networks. In this paper, we show how to write the full moment-hierarchy for any
non-zero-deficiency CRN obeying mass-action kinetics, in terms of equations for
the factorial moments (FM). Using these, we can recursively predict values for
lower moments from higher moments, reversing the procedure usually used to
solve moment hierarchies. We show, for non-trivial examples, that in this
manner we can predict to high accuracy, any moment of interest, for CRN's with
non-zero deficiency and non-factorizable steady states.


  One of the main aims in phylogenetics is the estimation of ancestral
sequences based on present-day data like, for instance, DNA alignments. One way
to estimate the data of the last common ancestor of a given set of species is
to first reconstruct a phylogenetic tree with some tree inference method and
then to use some method of ancestral state inference based on that tree. One of
the best-known methods both for tree inference as well as for ancestral
sequence inference is Maximum Parsimony (MP). In this manuscript, we focus on
this method and on ancestral state inference for fully bifurcating trees. In
particular, we investigate a conjecture published by Charleston and Steel in
1995 concerning the number of species which need to have a particular state,
say $a$, at a particular site in order for MP to unambiguously return $a$ as an
estimate for the state of the last common ancestor. We prove the conjecture for
all even numbers of character states, which is the most relevant case in
biology. We also show that the conjecture does not hold in general for odd
numbers of character states, but also present some positive results for this
case.


  Inverse problems in statistical physics are motivated by the challenges of
`big data' in different fields, in particular high-throughput experiments in
biology. In inverse problems, the usual procedure of statistical physics needs
to be reversed: Instead of calculating observables on the basis of model
parameters, we seek to infer parameters of a model based on observations. In
this review, we focus on the inverse Ising problem and closely related
problems, namely how to infer the coupling strengths between spins given
observed spin correlations, magnetisations, or other data. We review
applications of the inverse Ising problem, including the reconstruction of
neural connections, protein structure determination, and the inference of gene
regulatory networks. For the inverse Ising problem in equilibrium, a number of
controlled and uncontrolled approximate solutions have been developed in the
statistical mechanics community. A particularly strong method,
pseudolikelihood, stems from statistics. We also review the inverse Ising
problem in the non-equilibrium case, where the model parameters must be
reconstructed based on non-equilibrium statistics.


  We define and examine a model of epidemic propagation for a virus such as
Hepatitis C (with HIV co-infection) on a network of networks, namely the
network of French urban areas. One network level is that of the individual
interactions inside each urban area. The second level is that of the areas
themselves, linked by individuals travelling between these areas and
potentially helping the epidemic spread from one city to another. We choose to
encode the second level of the network as extra, special nodes in the first
level. We observe that such an encoding leads to sensible results in terms of
the extent and speed of propagation of an epidemic, depending on its source
point.


  The 70 kDa Heat Shock Proteins Hsp70 have several essential functions in
living systems, such as protecting cells against protein aggregation, assisting
protein folding, remodeling protein complexes and driving the translocation
into organelles. These functions require high affinity for non-specific
amino-acid sequences that are ubiquitous in proteins. It has been recently
shown that this high affinity, called ultra-affinity, depends on a process
driven out of equilibrium by ATP hydrolysis. Here we establish the
thermodynamic bounds for ultra-affinity, and further show that the same
reaction scheme can in principle be used both to strengthen and to weaken
affinities (leading in this case to infra-affinity). We show that cofactors are
essential to achieve affinity beyond the equilibrium range. Finally, biological
implications are discussed.


  DNA read mapping is a ubiquitous task in bioinformatics, and many tools have
been developed to solve the read mapping problem. However, there are two trends
that are changing the landscape of readmapping: First, new sequencing
technologies provide very long reads with high error rates (up to 15%). Second,
many genetic variants in the population are known, so the reference genome is
not considered as a single string over ACGT, but as a complex object containing
these variants. Most existing read mappers do not handle these new
circumstances appropriately.
  We introduce a new read mapper prototype called VATRAM that considers
variants. It is based on Min-Hashing of q-gram sets of reference genome
windows. Min-Hashing is one form of locality sensitive hashing. The variants
are directly inserted into VATRAMs index which leads to a fast mapping process.
Our results show that VATRAM achieves better precision and recall than
state-of-the-art read mappers like BWA under certain cirumstances. VATRAM is
open source and can be accessed at
https://bitbucket.org/Quedenfeld/vatram-src/.


  PHAST is a software package written in standard Fortran, with MPI and CUDA
extensions, able to efficiently perform parallel multicanonical Monte Carlo
simulations of single or multiple heteropolymeric chains, as coarse-grained
models for proteins. The outcome data can be straightforwardly analyzed within
its microcanonical Statistical Thermodynamics module, which allows for
computing the entropy, caloric curve, specific heat and free energies. As a
case study, we investigate the aggregation of heteropolymers bioinspired on
$A\beta_{25-33}$ fragments and their cross-seeding with $IAPP_{20-29}$
isoforms. Excellent parallel scaling is observed, even under numerically
difficult first-order like phase transitions, which are properly described by
the built-in fully reconfigurable force fields. Still, the package is free and
open source, this shall motivate users to readily adapt it to specific
purposes.


  As a cheap and safe antimalarial agent, chloroquine (CQ) has been used in the
battle against malaria for more than half century. However, the mechanism of CQ
action and resistance in Plasmodium falciparum remains elusive. Based on
further analysis of our published experimental results, we propose that the
mechanism of CQ action and resistance might be closely linked with
cell-cycle-associated amplified genomic-DNA fragments (CAGFs, singular form =
CAGF) as CQ induces CAGF production in P. falciparum, which could affect
multiple biological processes of the parasite, and thus might contribute to
parasite death and CQ resistance. Recently, we found that CQ induced one of
CAGFs, UB1- CAGF, might downregulate a probable P. falciparum cystine
transporter (Pfct) gene expression, which could be used to understand the
mechanism of CQ action and resistance in P. falciparum.


  Computer experiments that mirror the evolutionary dynamics of sexual and
asexual organisms as they occur in nature, tested features proposed to explain
the evolution of sexual recombination. Results show that this evolution is
better described as a network of interactions between possible sexual forms,
including diploidy, thelytoky, facultative sex, assortation, bisexuality, and
division of labor between the sexes, rather than a simple transition from
parthenogenesis to sexual recombination. Diploidy was shown to be fundamental
for the evolution of sex; bisexual reproduction emerged only among anisogamic
diploids with a synergistic division of reproductive labor; and facultative sex
was more likely to evolve among haploids practicing assortative mating. Looking
at the evolution of sex as a complex system through individual based
simulations, explains better the diversity of sexual strategies known to exist
in nature, compared to classical analytical models