  Non-coding RNAs are ubiquitous, but the discovery of new RNA gene sequences
far outpaces research on their structure and functional interactions. We mine
the evolutionary sequence record to derive precise information about function
and structure of RNAs and RNA-protein complexes. As in protein structure
prediction, we use maximum entropy global probability models of sequence
co-variation to infer evolutionarily constrained nucleotide-nucleotide
interactions within RNA molecules, and nucleotide-amino acid interactions in
RNA-protein complexes. The predicted contacts allow all-atom blinded 3D
structure prediction at good accuracy for several known RNA structures and
RNA-protein complexes. For unknown structures, we predict contacts in 160
non-coding RNA families. Beyond 3D structure prediction, evolutionary couplings
help identify important functional interactions, e.g., at switch points in
riboswitches and at a complex nucleation site in HIV. Aided by accelerating
sequence accumulation, evolutionary coupling analysis can accelerate the
discovery of functional interactions and 3D structures involving RNA.


  This work exploits the zeptoliter sensing volume of electrolyte-filled
nanopores to determine, simultaneously and in real time, the approximate shape,
volume, charge, rotational diffusion coefficient, and dipole moment of
individual proteins. We have developed the theory for a quantitative
understanding and analysis of modulations in ionic current that arise from
rotational dynamics of single proteins as they move through the electric field
inside a nanopore. The resulting multi-parametric information raises the
possibility to characterize, identify, and quantify individual proteins and
protein complexes in a mixture. This approach interrogates single proteins in
solution and determines parameters such as the approximate shape and dipole
moment, which are excellent protein descriptors and cannot be obtained
otherwise from single protein molecules in solution. Taken together, this
five-dimensional characterization of biomolecules at the single particle level
has the potential for instantaneous protein identification, quantification, and
possibly sorting with implications for structural biology, proteomics,
biomarker detection, and routine protein analysis.


  Over the past decades, molecular dynamics (MD) simulations of biomolecules
have become a mainstream biophysics technique. As the length and time scales
amenable to the MD method increase, shortcomings of the empirical force
fields---which have been developed and validated using relatively short
simulations of small molecules---become apparent. One common artifact is
aggregation of water-soluble biomolecules driven by artificially strong
charge--charge interactions. Here, we report a systematic refinement of
Lennard-Jones parameters (NBFIX) describing amine--carboxylate and
amine--phosphate interactions, which brings MD simulations of basic
peptide-mediated nucleic acids assemblies and lipid bilayer membranes in better
agreement with experiment. As our refinement neither affects the existing
parameterization of bonded interaction nor does it alter the solvation free
energies, it improves realism of an MD simulation without introducing
additional artifacts.


  It has been shown that a random-effects framework can be used to test the
association between a gene's expression level and the number of DNA copies of a
set of genes. This gene-set modelling framework was later applied to find
associations between mRNA expression and microRNA expression, by defining the
gene sets using target prediction information.
  Here, we extend the model introduced by Menezes et al (2009) to consider the
effect of not just copy number, but also of other molecular profiles such as
methylation changes and loss-of-heterozigosity (LOH), on gene expression
levels. We will consider again sets of measurements, to improve robustness of
results and increase the power to find associations. Our approach can be used
genome-wide to find associations, yields a test to help separate true
associations from noise and can include confounders.
  We apply our method to colon and to breast cancer samples, for which
genome-wide copy number, methylation and gene expression profiles are
available. Our findings include interesting gene expression-regulating
mechanisms, which may involve only one of copy number or methylation, or both
for the same samples. We even are able to find effects due to different
molecular mechanisms in different samples.
  Our method can equally well be applied to cases where other types of
molecular (high-dimensional) data are collected, such as LOH, SNP genotype and
microRNA expression data. Computationally efficient, it represents a flexible
and powerful tool to study associations between high-dimensional datasets. The
method is freely available via the SIM BioConductor package.


  Phylogenetic inference can potentially result in a more accurate tree using
data from multiple loci. However, if the loci are incongruent--due to events
such as incomplete lineage sorting or horizontal gene transfer--it can be
misleading to infer a single tree. To address this, many previous contributions
have taken a mechanistic approach, by modelling specific processes.
Alternatively, one can cluster loci without assuming how these incongruencies
might arise. Such "process-agnostic" approaches typically infer a tree for each
locus and cluster these. There are, however, many possible combinations of tree
distance and clustering methods; their comparative performance in the context
of tree incongruence is largely unknown. Furthermore, because standard model
selection criteria such as AIC cannot be applied to problems with a variable
number of topologies, the issue of inferring the optimal number of clusters is
poorly understood. Here, we perform a large-scale simulation study of
phylogenetic distances and clustering methods to infer loci of common
evolutionary history. We observe that the best-performing combinations are
distances accounting for branch lengths followed by spectral clustering or
Ward's method. We also introduce two statistical tests to infer the optimal
number of clusters and show that they strongly outperform the silhouette
criterion, a general-purpose heuristic. We illustrate the usefulness of the
approach by (i) identifying errors in a previous phylogenetic analysis of yeast
species and (ii) identifying topological incongruence among newly sequenced
loci of the globeflower fly genus Chiastocheta. We release treeCl, a new
program to cluster genes of common evolutionary history (http://git.io/treeCl).


  Some natural proteins display recurrent structural patterns. Despite being
highly similar at the tertiary structure level, repetitions within a single
repeat protein can be extremely variable at the sequence level. We propose a
mathematical definition of a repeat and investigate the occurrences of these in
different protein families. We found that long stretches of perfect repetitions
are infrequent in individual natural proteins, even for those which are known
to fold into structures of recurrent structural motifs. We found that natural
repeat proteins are indeed repetitive in their families, exhibiting abundant
stretches of 6 amino acids or longer that are perfect repetitions in the
reference family. We provide a systematic quantification for this
repetitiveness, and show that this form of repetitiveness is not exclusive of
repeat proteins, but also occurs in globular domains. A by-product of this work
is a fast classifier of proteins into families, which yields likelihood value
about a given protein belonging to a given family.


  Genetic interactions confer robustness on cells in response to genetic
perturbations. This often occurs through molecular buffering mechanisms that
can be predicted using, among other features, the degree of coexpression
between genes, commonly estimated through marginal measures of association such
as Pearson or Spearman correlation coefficients. However, marginal correlations
are sensitive to indirect effects and often partial correlations are used
instead. Yet, partial correlations convey no information about the (linear)
influence of the coexpressed genes on the entire multivariate system, which may
be crucial to discriminate functional associations from genetic interactions.
To address these two shortcomings, here we propose to use the edge weight
derived from the covariance decomposition over the paths of the associated gene
network. We call this new quantity the networked partial correlation and use it
to analyze genetic interactions in yeast.


  Molecular recognition between two double stranded (ds) DNA with homologous
sequences may not seem compatible with the B-DNA structure because the sequence
information is hidden when it is used for joining the two strands.
Nevertheless, it has to be invoked to account for various biological data.
Using quantum chemistry, molecular mechanics, and hints from recent genetics
experiments I show here that direct recognition between homologous dsDNA is
possible through formation of short quadruplexes due to direct complementary
hydrogen bonding of major groove surfaces in parallel alignment. The
constraints imposed by the predicted structures of the recognition units
determine the mechanism of complexation between long dsDNA. This mechanism and
concomitant predictions agree with available experimental data and shed light
upon the sequence effects and the possible involvement of topoisomerase II in
the recognition.


  In higher eukaryotes, alternative splicing is usually regulated by protein
factors, which bind to the pre-mRNA and affect the recognition of splicing
signals. There is recent evidence that the secondary structure of the pre-mRNA
may also play an important role in this process, either by facilitating or by
hindering the interaction with factors and small nuclear ribonucleoproteins
(snRNPs) that regulate splicing. Moreover, the secondary structure could play a
fundamental role in the splicing of yeast species, which lack many of the
regulatory splicing factors present in metazoans. This review describes the
steps in the analysis of the secondary structure of the pre-mRNA and its
possible relation to splicing. As a working example, we use the case of yeast
and the problem of the recognition of the 3-prime splice site.


  SMAR1 is a sensitive signaling molecule in p53 regulatory network which can
drive p53 network dynamics to three distinct states, namely, stabilized (two),
damped and sustain oscillation states. In the interaction of p53 network with
SMAR1, p53 network sees SMAR1 as a sub-network with its new complexes formed by
SMAR1, where SMAR1 is the central node, and fluctuations in SMAR1 concentration
is propagated as a stress signal throughout the network. Excess stress induced
by SMAR1 can drive p53 network dynamics to amplitude death scenario which
corresponds to apoptotic state. The permutation entropy calculated for normal
state is minimum indicating self-organized behavior, whereas for apoptotic
state, the value is maximum showing breakdown of self-organization. We also
show that the regulation of SMAR1 togather with other signaling molecules p300
and HDAC1 in the p53 regulatory network can be engineered to extend the range
of stress such that the system can be save from apoptosis.


  Protein structure prediction is a challenging and unsolved problem in
computer science. Proteins are the sequence of amino acids connected together
by single peptide bond. The combinations of the twenty primary amino acids are
the constituents of all proteins. In-vitro laboratory methods used in this
problem are very time-consuming, cost-intensive, and failure-prone. Thus,
alternative computational methods come into play. The protein structure
prediction problem is to find the three-dimensional native structure of a
protein, from its amino acid sequence. The native structure of a protein has
the minimum free energy possible and arguably determines the function of the
protein. In this study, we present the preliminaries of proteins and their
structures, protein structure prediction problem, and protein models. We also
give a brief overview on experimental and computational methods used in protein
structure prediction. This study will provide a fundamental knowledge to the
computer scientists who are intending to pursue their future research on
protein structure prediction problem.


  Erwin Schrodinger famously and presciently ascribed the vehicle transmitting
the hereditary information underlying life to an `aperiodic crystal'. We
compare and contrast this, only later discovered to be stored in the linear
biomolecule DNA, with the information bearing, layered quasi-one-dimensional
materials investigated by the emerging field of chaotic crystallography.
Despite differences in functionality, the same information measures capture
structure and novelty in both, suggesting an intimate coherence between the
information character of biotic and abiotic matter---a broadly applicable
physics of information. We review layered solids and consider three examples of
how information- and computation-theoretic techniques are being applied to
understand their structure. In particular, (i) we review recent efforts to
apply new kinds of information measures to quantify disordered crystals; (ii)
we discuss the structure of ice I in information-theoretic terms; and (iii) we
recount recent experimental results on tris(bicyclo[2.1.1]hexeno)benzene TBHB),
showing how an information-theoretic analysis yields additional insight into
its structure. We then illustrate a new Second Law of Thermodynamics that
describes information processing in active low-dimensional materials, reviewing
Maxwell's Demon and a new class of molecular devices that act as information
catalysts. Lastly, we conclude by speculating on how these ideas from
informational materials science may impact biology.


  Deep convolutional neural networks comprise a subclass of deep neural
networks (DNN) with a constrained architecture that leverages the spatial and
temporal structure of the domain they model. Convolutional networks achieve the
best predictive performance in areas such as speech and image recognition by
hierarchically composing simple local features into complex models. Although
DNNs have been used in drug discovery for QSAR and ligand-based bioactivity
predictions, none of these models have benefited from this powerful
convolutional architecture. This paper introduces AtomNet, the first
structure-based, deep convolutional neural network designed to predict the
bioactivity of small molecules for drug discovery applications. We demonstrate
how to apply the convolutional concepts of feature locality and hierarchical
composition to the modeling of bioactivity and chemical interactions. In
further contrast to existing DNN techniques, we show that AtomNet's application
of local convolutional filters to structural target information successfully
predicts new active molecules for targets with no previously known modulators.
Finally, we show that AtomNet outperforms previous docking approaches on a
diverse set of benchmarks by a large margin, achieving an AUC greater than 0.9
on 57.8% of the targets in the DUDE benchmark.


  Studies of the genetic loci that contribute to variation in gene expression
frequently identify loci with broad effect on gene expression: expression
quantitative trait locus (eQTL) hotspots. We describe a set of exploratory
graphical methods as well as a formal likelihood-based test for assessing
whether a given hotspot is due to one or multiple polymorphisms. We first look
at the pattern of effects of the locus on the expression traits that map to the
locus: the direction of the effects, as well as the degree of dominance. A
second technique is to focus on the individuals that exhibit no recombination
event in the region, apply dimensionality reduction (such as with linear
discriminant analysis) and compare the phenotype distribution in the
non-recombinants to that in the recombinant individuals: If the recombinant
individuals display a different expression pattern than the non-recombinants,
this indicates the presence of multiple causal polymorphisms. In the formal
likelihood-based test, we compare a two-locus model, with each expression trait
affected by one or the other locus, to a single-locus model. We apply our
methods to a large mouse intercross with gene expression microarray data on six
tissues.


  Protein assembly plays an important role in the regulation of biological
systems. The cytoskeleton assembly activity is provided by the binding
cofactors GTP (guanidine triphosphate) or ATP(adenosine triphosphate) to
monomeric protein, and is initiated by assembling the monomeric proteins. The
binding GTP or ATP is hydrolyzed to GDP (guanidine diphosphate) or ADP
(adenosine diphosphate) by the monomeric enzymatic activity. This self-limited
assembly is characteristic of the cytoskeleton. To quantitatively evaluate the
assembly kinetics, we propose a nonlinear and non-equilibrium kinetic model,
with the nonlinearity provided by the fluctuation in monomer concentrations
during the diffusion. Numerical simulations suggest that the assembly and
disassembly oscillates in a chaos-like manner. We use a kinetic analysis of the
center manifold around the critical point to show that minimal increases in
ATP/GTP concentrations may lead to some attenuation in the amplification of
these fluctuations. The present model and our application of center manifold
theory illustrate a unique feature of protein assemblies, and our stability
analysis provides an analytical methodology for the biological reaction system.


  Cannibalism, which is the act of killing and at least partial consumption of
conspecifics, is ubiquitous in nature. Mathematical models have considered
cannibalism in the predator primarily, and show that predator cannibalism in
two species ODE models provides a strong stabilizing effect. There is strong
ecological evidence that cannibalism exists among prey as well, yet this
phenomenon has been much less investigated. In the current manuscript, we
investigate both the ODE and spatially explicit forms of a Holling-Tanner
model, with ratio dependent functional response. We show that cannibalism in
the predator provides a stabilizing influence as expected. However, when
cannibalism in the prey is considered, we show that it cannot stabilise the
unstable interior equilibrium in the ODE case, but can destabilise the stable
interior equilibrium. In the spatially explicit case, we show that in certain
parameter regime, prey cannibalism can lead to pattern forming Turing dynamics,
which is an impossibility without it. Lastly we consider a stochastic prey
cannibalism rate, and find that it can alter both spatial patterns, as well as
limit cycle dynamics.


  Renaturation of complementary single strands of DNA is one of the important
processes that requires better understanding in the view of molecular biology
and biological physics. Here we develop a stochastic dynamical model on the DNA
renaturation. According to our model there are at least three steps in the
renaturation process viz. incorrect-contact formation, correct-contact
formation and nucleation, and zipping. Most of the earlier two-state models
combined nucleation with incorrect-contact formation step. In our model we
suggest that it is considerably meaningful when we combine the nucleation with
the zipping since nucleation is the initial step of zipping and the nucleated
and zipping molecules are indistinguishable. Incorrect-contact formation step
is a pure three-dimensional diffusion controlled collision process. Whereas
nucleation involves several rounds of one-dimensional slithering dynamics of
one single strand of DNA on the other complementary strand in the process of
searching for the correct-contact and then initiate nucleation. Upon
nucleation, the stochastic zipping follows to generate a fully renatured double
stranded DNA. It seems that the square-root dependency of the overall
renaturation rate constant on the length of reacting single strands originates
mainly from the geometric constraints in the diffusion controlled
incorrect-contact formation step. Further the inverse scaling of the
renaturation rate on the viscosity of the reaction medium also originates from
the incorrect-contact formation step. On the other hand the inverse scaling of
the renaturation rate with the sequence complexity originates from the
stochastic zipping which involves several rounds of crossing over the
free-energy barrier at microscopic levels.


  Biochemical processes typically involve many chemical species, some in
abundance and some in low molecule numbers. Here we first identify the rate
constant limits under which the concentrations of a given set of species will
tend to infinity (the abundant species) while the concentrations of all other
species remains constant (the non-abundant species). Subsequently we prove that
in this limit, the fluctuations in the molecule numbers of non-abundant species
are accurately described by a hybrid stochastic description consisting of a
chemical master equation coupled to deterministic rate equations. This is a
reduced description when compared to the conventional chemical master equation
which describes the fluctuations in both abundant and non-abundant species. We
show that the reduced master equation can be solved exactly for a number of
biochemical networks involving gene expression and enzyme catalysis, whose
conventional chemical master equation description is analytically impenetrable.
We use the linear noise approximation to obtain approximate expressions for the
difference between the variance of fluctuations in the non-abundant species as
predicted by the hybrid approach and by the conventional chemical master
equation. Furthermore we show that surprisingly, irrespective of any separation
in the mean molecule numbers of various species, the conventional and hybrid
master equations exactly agree for a class of chemical systems.


  The quantitative characterization of mutational landscapes is a task of
outstanding importance in evolutionary and medical biology: It is, e.g., of
central importance for our understanding of the phenotypic effect of mutations
related to disease and antibiotic drug resistance. Here we develop a novel
inference scheme for mutational landscapes, which is based on the statistical
analysis of large alignments of homologs of the protein of interest. Our method
is able to capture epistatic couplings between residues, and therefore to
assess the dependence of mutational effects on the sequence context where they
appear. Compared to recent large-scale mutagenesis data of the beta-lactamase
TEM-1, a protein providing resistance against beta-lactam antibiotics, our
method leads to an increase of about 40% in explicative power as compared to
approaches neglecting epistasis. We find that the informative sequence context
extends to residues at native distances of about 20{\AA} from the mutated site,
reaching thus far beyond residues in direct physical contact.


  Cilia and flagella often exhibit synchronized behavior; this includes phase
locking, as seen in {\it Chlamydomonas}, and metachronal wave formation in the
respiratory cilia of higher organisms. Since the observations by Gray and
Rothschild of phase synchrony of nearby swimming spermatozoa, it has been a
working hypothesis that synchrony arises from hydrodynamic interactions between
beating filaments. Recent work on the dynamics of physically separated pairs of
flagella isolated from the multicellular alga {\it Volvox} has shown that
hydrodynamic coupling alone is sufficient to produce synchrony. However, the
situation is more complex in unicellular organisms bearing few flagella. We
show that flagella of {\it Chlamydomonas} mutants deficient in filamentary
connections between basal bodies display markedly different synchronization
from the wild type. We perform micromanipulation on configurations of flagella
and conclude that a mechanism, internal to the cell, must provide an additional
flagellar coupling. In naturally occurring species with 4, 8, or even 16
flagella, we find diverse symmetries of basal body positioning and of the
flagellar apparatus that are coincident with specific gaits of flagellar
actuation, suggesting that it is a competition between intracellular coupling
and hydrodynamic interactions that ultimately determines the precise form of
flagellar coordination in unicellular algae.


  Despite the biological importance of non-coding RNA, their structural
characterization remains challenging. Making use of the rapidly growing
sequence databases, we analyze nucleotide coevolution across homologous
sequences via Direct-Coupling Analysis to detect nucleotide-nucleotide
contacts. For a representative set of riboswitches, we show that the results of
Direct-Coupling Analysis in combination with a generalized Nussinov algorithm
systematically improve the results of RNA secondary structure prediction beyond
traditional covariance approaches based on mutual information. Even more
importantly, we show that the results of Direct-Coupling Analysis are enriched
in tertiary structure contacts. By integrating these predictions into molecular
modeling tools, systematically improved tertiary structure predictions can be
obtained, as compared to using secondary structure information alone.


  Astrophysical ionizing radiation events such as supernovae, gamma-ray bursts,
and solar proton events have been recognized as a potential threat to life on
Earth, primarily through depletion of stratospheric ozone and subsequent
increase in solar UV radiation at Earth's surface and in the upper levels of
the ocean. Other work has also considered the potential impact of nitric acid
rainout, concluding that no significant threat is likely. Not yet studied
to-date is the potential impact of ozone produced in the lower atmosphere
following an ionizing radiation event. Ozone is a known irritant to organisms
on land and in water and therefore may be a significant additional hazard.
Using previously completed atmospheric chemistry modeling we have examined the
amount of ozone produced in the lower atmosphere for the case of a gamma-ray
burst and find that the values are too small to pose a significant additional
threat to the biosphere. These results may be extended to other ionizing
radiation events, including supernovae and extreme solar proton events.


  It is well known that cyclins are a family of proteins that control
cell-cycle progression by activating cyclin-dependent kinase. Based on our
experimental results, we propose here a novel hypothesis that certain amplified
genomic-DNA fragments (AGFs) may also be required for the cell cycle
progression of eukaryotic cells and thus can be named as cell-cycle-associated
AGFs (CAGFs). Like fluctuation in cyclin levels during cell cycle progression,
these CAGFs are amplified and degraded at different points of the cell cycle.
The functions of CAGFs are unknown, but we speculate that CAGFs might be
involved in regulation of gene expression, genome protection, and formation of
certain macromolecular complexes required for the dynamic genome architecture
during cell cycle progression. Our experimental results also show that
chloroquine induces the production of CAGFs in Plasmodium falciparum,
suggesting that targeting cell cycle progression can be the primary mechanism
of chloroquine's antimalarial, anticancer, and immunomodulatory actions.


  Effective and controlled drug delivery systems with on-demand release and
targeting abilities have received enormous attention for biomedical
applications. Here, we describe a novel enzyme-based cap system for mesoporous
silica nanoparticles (MSNs) that is directly combined with a targeting ligand
via bio-orthogonal click chemistry. The capping system is based on the
pH-responsive binding of an aryl-sulfonamide-functionalized MSN and the enzyme
carbonic anhydrase (CA). An unnatural amino acid (UAA) containing a norbornene
moiety was genetically incorporated into CA. This UAA allowed for the
site-specific bio-orthogonal attachment of even very sensitive targeting
ligands such as folic acid and anandamide. This leads to specific
receptor-mediated cell and stem cell uptake. We demonstrate the successful
delivery and release of the chemotherapeutic agent Actinomycin D to KB cells.
This novel nanocarrier concept provides a promising platform for the
development of precisely controllable and highly modular theranostic systems.


  Gene assembly is an important step in functional analysis of shotgun
metagenomic data. Nonetheless, strain aware assembly remains a challenging
task, as current assembly tools often fail to distinguish among strain variants
or require closely related reference genomes of the studied species to be
available. We have developed Snowball, a novel strain aware and reference-free
gene assembler for shotgun metagenomic data. It uses profile hidden Markov
models (HMMs) of gene domains of interest to guide the assembly. Our assembler
performs gene assembly of individual gene domains based on read overlaps and
error correction using read quality scores at the same time, which result in
very low per-base error rates. The software runs on a user-defined number of
processor cores in parallel, runs on a standard laptop and is freely available
for installation under Linux or OS X on:
https://github.com/algbioi/snowball/wiki


  Pathway Tools is a bioinformatics software environment with a broad set of
capabilities. The software provides genome-informatics tools such as a genome
browser, sequence alignments, a genome-variant analyzer, and
comparative-genomics operations. It offers metabolic-informatics tools, such as
metabolic reconstruction, quantitative metabolic modeling, prediction of
reaction atom mappings, and metabolic route search. Pathway Tools also provides
regulatory-informatics tools, such as the ability to represent and visualize a
wide range of regulatory interactions. The software creates and manages a type
of organism-specific database called a Pathway/Genome Database (PGDB), which
the software enables database curators to interactively edit. It supports web
publishing of PGDBs and provides a large number of query, visualization, and
omics-data analysis tools. Scientists around the world have created more than
9,800 PGDBs by using Pathway Tools, many of which are curated databases for
important model organisms. Those PGDBs can be exchanged using a peer-to-peer
database-sharing system called the PGDB Registry.


  To study how plants respond to conditions under which the carbohydrates
needed for growth and development are exhausted, this project focuses on
cloning and functional characterization of a cytosolic protein kinase, STY46
that may potentially be involved in the sugar starvation response (SSR).
Generating A. thaliana STY46 knockout (KO) and overexpression (OE) lines is the
main method in this project to determine the function of the protein. A SALK
line (SALK116340) was confirmed to contain a TDNA insert in the STY46 gene and
the aim was to prove that this reduced the expression of this gene. Highly
intact RNA was isolated and reverse transcription PCR performed to determine
the level of transcript of STY46 in STY46KO. However the results were
inconclusive. This may be due to not enough amount of RNA use for the cDNA
reverse transcription. To develop an overexpression line, STY46 was cloned from
wild type Arabidopsis genotype. Primers were designed to flank the open read
frame of the STY46. They also contained restriction enzyme recognition sites
and a myctag engineered at their 5' ends to facilitate subsequent synthesis of
a plant transformation construct. A cDNA fragment of the expected size was
successfully amplified by PCR.


  We analyze the translocation of a charged particle across an Alpha-Hemolysin
(aHL) pore in the framework of a driven diffusion over an extended energy
barrier generated by the electrical charges of the aHL. A one-dimensional
electrostatic potential is extracted from the full 3D solution of the Poisson's
equation. We characterize the particle transport under the action of a constant
forcing by studying the statistics of the translocation time. We derive an
analytical expression of translocation time average that compares well with the
results from Brownian dynamic simulations of driven particles over the
electrostatic potential. Moreover, we show that the translocation time
distributions can be perfectly described by a simple theory which replaces the
true barrier by an equivalent structureless square barrier. Remarkably our
approach maintains its accuracy also for low-applied voltage regimes where the
usual inverse-Gaussian approximation fails. Finally we discuss how the
comparison between the simulated time distributions and their theoretical
prediction results to be greatly simplified when using the notion of the
empirical Laplace transform technique.


  We developed a strategy to control p53 levels when the p53 concentration
falls below a desired threshold . Based on an established p53-MDM2 mathematical
model that also incorporates microRNA regulators (p53-MDM2-miRNA), we applied a
control action in closed-loop with respect to this model. In particular, we
employ various architectures of synthetic circuit with integral feedback
control, illustrating varying degree in adaptive response to step-wise
perturbation of p53.


  We present a framework based on conditional entropy and the Dirichlet
distribution for classifying chromosomes based on the degree to which they use
synonymous codons uniformly or preferentially, that is, whether or not codons
that code for an amino acid appear with the same relative frequency. Applying
the approach to a large collection of annotated bacterial chromosomes reveals
three distinct groups of bacteria.


  Gene regulatory network (GRN) modeling is a well-established theoretical
framework for the study of cell-fate specification during developmental
processes. Recently, dynamical models of GRNs have been taken as a basis for
formalizing the metaphorical model of Waddington's epigenetic landscape,
providing a natural extension for the general protocol of GRN modeling. In this
contribution we present in a coherent framework a novel implementation of two
previously proposed general frameworks for modeling the Epigenetic Attractors
Landscape associated with boolean GRNs: the inter-attractor and inter-state
transition approaches. We implement novel algorithms for estimating
inter-attractor transition probabilities without necessarily depending on
intensive single-event simulations. We analyze the performance and sensibility
to parameter choices of the algorithms for estimating inter-attractor
transition probabilities using three real GRN models. Additionally, we present
a side-by-side analysis of downstream analysis tools such as the attractors'
temporal and global ordering in the EAL. Overall, we show how the methods
complement each other using a real case study: a cellular-level GRN model for
epithelial carcinogenesis. We expect the toolkit and comparative analyses put
forward here to be a valuable additional re- source for the systems biology
community interested in modeling cellular differentiation and reprogramming
both in normal and pathological developmental processes.


  Shortly after the determination of the first protein x-ray crystal
structures, researchers analyzed their cores and reported packing fractions
$\phi \approx 0.75$, a value that is similar to close packing equal-sized
spheres. A limitation of these analyses was the use of `extended atom' models,
rather than the more physically accurate `explicit hydrogen' model. The
validity of using the explicit hydrogen model is proved by its ability to
predict the side chain dihedral angle distributions observed in proteins. We
employ the explicit hydrogen model to calculate the packing fraction of the
cores of over $200$ high resolution protein structures. We find that these
protein cores have $\phi \approx 0.55$, which is comparable to random
close-packing of non-spherical particles. This result provides a deeper
understanding of the physical basis of protein structure that will enable
predictions of the effects of amino acid mutations and design of new functional
proteins.


  A great variety of biologically relevant monolayers present phase coexistence
characterized by domains formed by lipids in an ordered phase state dispersed
in a continuous, disordered phase. The difference in surface densities between
these phases originates inter-domain dipolar interactions, which are relevant
for the determination of the spacial distribution of domains, as well as their
dynamics. In this work, we propose a novel manner of estimating the dipolar
repulsion using a passive method that involves the analysis of images of the
monolayer with phase coexistence. The method is based on the comparison of the
pair correlation function obtained from experiments with that obtained from
Brownian dynamics simulations of a model system. As an example, we determined
the difference in dipolar density of a binary monolayer of DSPC/DMPC at the
air-water interface from the analysis of the radial distribution of domains,
and the results are compared with those obtained by surface potential
determinations. A systematic analysis for experimentally relevant parameter
range is given, which may be used as a working curve for obtaining the dipolar
repulsion in different systems.


  The movie Ben-Hur highlights the dynamics of contagion associated with
leprosy, a pattern of forced aggregation driven by the emergence of symptoms
and the fear of contagion. The 2014 Ebola outbreaks reaffirmed the dynamics of
redistribution among symptomatic and asymptomatic or non-infected individuals
as a way to avoid contagion. In this manuscript, we explore the establishment
of clusters of infection via density-dependence avoidance (diffusive
instability). We illustrate this possibility in two ways: using a
phenomenological driven model where disease incidence is assumed to be a
decreasing function of the size of the symptomatic population and with a model
that accounts for the deliberate movement of individuals in response to a
gradient of symptomatic infectious individuals. The results in this manuscript
are preliminary but indicative of the role that behavior, here modeled in crude
simplistic ways, may have on disease dynamics, particularly on the spatial
redistribution of epidemiological classes.


  In this paper, we study the selectivity of the potassium channel KcsA by a
recently developed image-charge solvation method(ICSM) combined with molecular
dynamics simulations. The hybrid solvation model in the ICSM is able to
demonstrate atomistically the function of the selectivity filter of the KcsA
channel when potassium and sodium ions are considered and their distributions
inside the filter are simulated. Our study also shows that the reaction field
effect, explicitly accounted for through image charge approximation in the ICSM
model, is necessary in reproducing the correct selectivity property of the
potassium channels.


  The use of illicit drugs has been on the rise in United States. It is very
detrimental on society, as fatal overdose is the fourth leading cause of death
in the United States, which is about the same as motor vehicle crashes. Of all
illicit drugs, one drug that has an severe adverse effect on a community as a
whole is heroin. This paper will discuss two mathematical models- the White and
Comiskey model and a newly introduced model proposed by the author, describing
heroin use within a fixed community. We will show the existence of stable
equilibrium from both models, suggesting both a situation where heroin use is
eradicated and one where it remains an endemic.


  Modern biomedicine is challenged to predict the effects of genetic variation.
Systematic functional assays of point mutants of proteins have provided
valuable empirical information, but vast regions of sequence space remain
unexplored. Fortunately, the mutation-selection process of natural evolution
has recorded rich information in the diversity of natural protein sequences.
Here, building on probabilistic models for correlated amino-acid substitutions
that have been successfully applied to determine the three-dimensional
structures of proteins, we present a statistical approach for quantifying the
contribution of residues and their interactions to protein function, using a
statistical energy, the evolutionary Hamiltonian. We find that these
probability models predict the experimental effects of mutations with
reasonable accuracy for a number of proteins, especially where the selective
pressure is similar to the evolutionary pressure on the protein, such as
antibiotics.


  Model reduction of biochemical networks relies on the knowledge of slow and
fast variables. We provide a geometric method, based on the Newton polytope, to
identify slow variables of a biochemical network with polynomial rate
functions. The gist of the method is the notion of tropical equilibration that
provides approximate descriptions of slow invariant manifolds. Compared to
extant numerical algorithms such as the intrinsic low dimensional manifold
method, our approach is symbolic and utilizes orders of magnitude instead of
precise values of the model parameters. Application of this method to a large
collection of biochemical network models supports the idea that the number of
dynamical variables in minimal models of cell physiology can be small, in spite
of the large number of molecular regulatory actors.


  Motivated by nucleation and molecular aggregation in physical, chemical and
biological settings, we present an extension to a thorough analysis of the
stochastic self-assembly of a fixed number of identical particles in a finite
volume. We study the statistic of times it requires for maximal clusters to be
completed, starting from a pure-monomeric particle configuration. For finite
volume, we extend previous analytical approaches to the case of arbitrary
size-dependent aggregation and fragmentation kinetic rates. For larger volume,
we develop a scaling framework to study the behavior of the first assembly time
as a function of the total quantity of particles.
  We find that the mean time to first completion of a maximum-sized cluster may
have surprisingly a very weak dependency on the total number of particles. We
highlight how the higher statistic (variance, distribution) of the first
passage time may still help to infer key parameters (such as the size of the
maximum cluster) from data. And last but not least, we present a framework to
quantify the formation of cluster of macroscopic size, whose formation is
(asymptotically) very unlikely and occurs as a large deviation phenomenon from
the mean-field limit. We argue that this framework is suitable to describe
phase transition phenomena, as inherent infrequent stochastic processes, in
contrast to classical nucleation theory.


  Biochemical networks related to similar functional pathways are often
correlated due to cross-talk among the homologous proteins in the different
networks. Using a stochastic framework, we address the functional significance
of the cross-talk between two pathways. Our theoretical analysis on generic
MAPK pathways reveals cross-talk is responsible for developing coordinated
fluctuations between the pathways. The extent of correlation evaluated in terms
of the information theoretic measure provides directionality to net information
propagation. Stochastic time series and scattered plot suggest that the
cross-talk generates synchronization within a cell as well as in a cellular
population. Depending on the number of input and output, we identify signal
integration and signal bifurcation motif that arise due to inter-pathway
connectivity in the composite network. Analysis using partial information
decomposition quantifies the net synergy in the information propagation through
these branched pathways.


  Spatial spread of infectious diseases among populations via the mobility of
humans is highly stochastic and heterogeneous. Accurate forecast/mining of the
spread process is often hard to be achieved by using statistical or mechanical
models. Here we propose a new reverse problem, which aims to identify the
stochastically spatial spread process itself from observable information
regarding the arrival history of infectious cases in each subpopulation. We
solved the problem by developing an efficient optimization algorithm based on
dynamical programming, which comprises three procedures: i, anatomizing the
whole spread process among all subpopulations into disjoint componential
patches; ii, inferring the most probable invasion pathways underlying each
patch via maximum likelihood estimation; iii, recovering the whole process by
assembling the invasion pathways in each patch iteratively, without burdens in
parameter calibrations and computer simulations. Based on the entropy theory,
we introduced an identifiability measure to assess the difficulty level that an
invasion pathway can be identified. Results on both artificial and empirical
metapopulation networks show the robust performance in identifying actual
invasion pathways driving pandemic spread.


  Flexible information routing fundamentally underlies the function of many
biological and artificial networks. Yet, how such systems may specifically
communicate and dynamically route information is not well understood. Here we
identify a generic mechanism to route information on top of collective
dynamical reference states in complex networks. Switching between collective
dynamics induces flexible reorganization of information sharing and routing
patterns, as quantified by delayed mutual information and transfer entropy
measures between activities of a network's units. We demonstrate the power of
this generic mechanism specifically for oscillatory dynamics and analyze how
individual unit properties, the network topology and external inputs coact to
systematically organize information routing. For multi-scale, modular
architectures, we resolve routing patterns at all levels. Interestingly, local
interventions within one sub-network may remotely determine non-local
network-wide communication. These results help understanding and designing
information routing patterns across systems where collective dynamics co-occurs
with a communication function.


  The ectoparasitic mite Varroa destructor has become one of the major
worldwide threats for apiculture. Varroa destructor attacks the honey bee Apis
mellifera weakening its host by sucking hemolymph. However, the damage to bee
colonies is not strictly related to the parasitic action of the mite but it
derives, above all, from its action as vector increasing the trasmission of
many viral diseases such as acute paralysis (ABPV) and deformed wing viruses
(DWV), that are considered among the main causes of CCD (Colony Collapse
Disorder). In this work we discuss an SI model that describes how the presence
of the mite affects the epidemiology of these viruses on adult bees. We
characterize the system behavior, establishing that ultimately either only
healthy bees survive, or the disease becomes endemic and mites are wiped out.
Another dangerous alternative is the Varroa invasion scenario with the
extinction of healthy bees. The final possible configuration is the coexistence
equilibrium in which honey bees share their infected hive with mites. The
analysis is in line with some observed facts in natural honey bee colonies.
Namely, these diseases are endemic. Further, if the mite population is present,
necessarily the viral infection occurs. The findings of this study indicate
that a low horizontal transmission rate of the virus among honey bees in
beehives will help in protecting bee colonies from Varroa infestation and viral
epidemics.


  Many of the observed species interactions embedded in ecological communities
are not permanent, but are characterized by temporal changes that are observed
along with abiotic and biotic variations. While work has been done describing
and quantifying these changes, little is known about their consequences for
species coexistence. Here, we investigate the extent to which changes of
species composition impacts the likelihood of persistence of the predator-prey
community in the highly seasonal Bialowieza Primeval Forest (NE Poland), and
the extent to which seasonal changes of species interactions (predator diet)
modulate the expected impact. This likelihood is estimated extending recent
developments on the study of structural stability in ecological communities. We
find that the observed species turnover strongly varies the likelihood of
community persistence between summer and winter. Importantly, we demonstrate
that the observed seasonal interaction changes minimize the variation in the
likelihood of persistence associated with species turnover across the year. We
find that these community dynamics can be explained as the coupling of
individual species to their environment by minimizing both the variation in
persistence conditions and the interaction changes between seasons. Our results
provide a homeostatic explanation for seasonal species interactions, and
suggest that monitoring the association of interactions changes with the level
of variation in community dynamics can provide a good indicator of the response
of species to environmental pressures.


  Polymer translocation through a nanometer-scale pore assisted by chaperones
binding to the polymer is a process encountered in vivo for proteins. Studying
the relevant models by computer simulations is computationally demanding.
Accordingly, previous studies are either for stiff polymers in three dimensions
or flexible polymers in two dimensions. Here, we study chaperone-assisted
translocation of flexible polymers in three dimensions using Langevin dynamics.
We show that differences in binding mechanisms, more specifically, whether a
chaperone can bind to a single or multiple sites on the polymer, lead to
substantial differences in translocation dynamics in three dimensions. We show
that the single-binding mode leads to dynamics that is very much like that in
the constant-force driven translocation and accordingly mainly determined by
tension propagation on the cis side. We obtain $\beta \approx 1.26$ for the
exponent for the scaling of the translocation time with polymer length. This
fairly low value can be explained by the additional friction due to binding
particles. The multiple-site binding leads to translocation whose dynamics is
mainly determined by the trans side. For this process we obtain $\beta \approx
1.36$. This value can be explained by our derivation of $\beta = 4/3$ for
constant-bias translocation, where translocated polymer segments form a globule
on the trans side. Our results pave the way for understanding and utilizing
chaperone-assisted translocation where variations in microscopic details lead
to rich variations in the emerging dynamics.


  General genealogical processes such as $\Lambda$- and $\Xi$-coalescents,
which respectively model multiple and simultaneous mergers, have important
applications in studying marine species, strong positive selection, recurrent
selective sweeps, strong bottlenecks, large sample sizes, and so on. Recently,
there has been significant progress in developing useful inference tools for
such general models. In particular, inference methods based on the site
frequency spectrum (SFS) have received noticeable attention. Here, we derive a
new formula for the expected SFS for general $\Lambda$- and $\Xi$-coalescents,
which leads to an efficient algorithm. For time-homogeneous coalescents, the
runtime of our algorithm for computing the expected SFS is $O(n^2)$, where $n$
is the sample size. This is a factor of $n^2$ faster than the state-of-the-art
method. Furthermore, in contrast to existing methods, our method generalizes to
time-inhomogeneous $\Lambda$- and $\Xi$-coalescents with measures that
factorize as $\Lambda(dx)/\zeta(t)$ and $\Xi(dx)/\zeta(t)$, respectively, where
$\zeta$ denotes a strictly positive function of time. The runtime of our
algorithm in this setting is $O(n^3)$. We also obtain general theoretical
results for the identifiability of the $\Lambda$ measure when $\zeta$ is a
constant function, as well as for the identifiability of the function $\zeta$
under a fixed $\Xi$ measure.


  The concept of fitness as a measure for a species's success in natural
selection is central to the theory of evolution. We here investigate how
reproduction rates which are not constant but vary in response to environmental
fluctuations, influence a species' prosperity and thereby its fitness.
Interestingly, we find that not only larger growth rates but also reduced
sensitivities to environmental changes substantially increase the fitness.
Thereby, depending on the noise level of the environment, it might be an
evolutionary successful strategy to minimize this sensitivity rather than to
optimize the reproduction speed. Also for neutral evolution, where species with
exactly the same properties compete, variability in the growth rates plays a
crucial role. The time for one species to fixate is strongly reduced in the
presence of environmental noise. Hence, environmental fluctuations constitute a
possible explanation for effective population sizes inferred from genetic data
that often are much smaller than the census population size.


  Proteins are the basic building blocks of life. They usually perform
functions by folding to a particular structure. Understanding the folding
process could help the researchers to understand the functions of proteins and
could also help to develop supplemental proteins for people with deficiencies
and gain more insight into diseases associated with troublesome folding
proteins. Experimental methods are both expensive and time consuming. In this
thesis I introduce a new machine learning based method to predict the protein
structure. The new method improves the performance from two directions:
creating accurate protein alignments and predicting accurate protein contacts.
First, I present an alignment framework MRFalign which goes beyond
state-of-the-art methods and uses Markov Random Fields to model a protein
family and align two proteins by aligning two MRFs together. Compared to other
methods, that can only model local-range residue correlation, MRFs can model
long-range residue interactions and thus, encodes global information in a
protein. Secondly, I present a Group Graphical Lasso method for contact
prediction that integrates joint multi-family Evolutionary Coupling analysis
and supervised learning to improve accuracy on proteins without many sequence
homologs. Different from single-family EC analysis that uses residue
co-evolution information in only the target protein family, our joint EC
analysis uses residue co-evolution in both the target family and its related
families, which may have divergent sequences but similar folds. Our method can
also integrate supervised learning methods to further improve accuracy. We
evaluate the performance of both methods including each of its components on
large public benchmarks. Experiments show that our methods can achieve better
accuracy than existing state-of-the-art methods under all the measurements on
most of the protein classes.


  In this paper we present a mode that incorporates two levels of dynamic
structuring of the coffee rust disease (caused by the rust fungus, Hemileia
vastatrix). First, two distinct spatial scales of transmission dynamics are
interrogated with respect to the resultant structure of catastrophic
transitions of epidemics. Second the effect of management style, especially the
well-known issue of shade management, is added to the base-line model. The
final structure includes a simple one-to-one functional structure of disease
incidence as a function as well as a classical critical transition structure,
and finally a hysteretic loop. These qualitative structures accord well with
the recent history of the coffee rust disease.


  Many organisms possess both a cell cycle to control DNA replication and a
circadian clock to anticipate changes between day and night. In some cases,
these two rhythmic systems are known to be coupled by specific,
cross-regulatory interactions. Here, we use mathematical modeling to show that,
additionally, the cell cycle generically influences circadian clocks in a
non-specific fashion: The regular, discrete jumps in gene-copy number arising
from DNA replication during the cell cycle cause a periodic driving of the
circadian clock, which can dramatically alter its behavior and impair its
function. A clock built on negative transcriptional feedback either phase locks
to the cell cycle, so that the clock period tracks the cell division time, or
exhibits erratic behavior. We argue that the cyanobacterium Synechococcus
elongatus has evolved two features that protect its clock from such
disturbances, both of which are needed to fully insulate it from the cell cycle
and give it its observed robustness: a phosphorylation-based protein
modification oscillator, together with its accompanying push-pull read-out
circuit that responds primarily to the ratios of the different phosphoforms,
makes the clock less susceptible to perturbations in protein synthesis; and the
presence of multiple, asynchronously replicating copies of the same chromosome
diminishes the effect of replicating any single copy of a gene.


  Recent data support the notion that a group of key transcriptional regulators
involved in tumorigenesis, including MYC, p53, E2F1, and BMI1, share an
intriguing capacity to simultaneously regulate metabolism and cell cycle. Here,
we show that another factor, the multifunctional protein E4F1, directly
controls genes involved in mitochondria functions and cell-cycle checkpoints,
including Chek1, a major component of the DNA damage response. Coordination of
these cellular functions by E4F1 appears essential for the survival of
p53-deficient transformed cells. Acute inactivation of E4F1 in these cells
results in CHK1-dependent checkpoint deficiency and multiple mitochondrial
dysfunctions that lead to increased ROS production, energy stress, and
inhibition of de novo pyrimidine synthesis. This deadly cocktail leads to the
accumulation of uncompensated oxidative damage to proteins and extensive DNA
damage, ending in cell death. This supports the rationale of therapeutic
strategies simultaneously targeting mitochondria and CHK1 for selective killing
of p53-deficient cancer cells.


  The question of what determines genetic diversity both between and within
species has long remained unsolved by the modern evolutionary theory (MET).
However, it has not deterred researchers from producing interpretations of
genetic diversity by using MET. We here examine the two key experimental
observations of genetic diversity made in the 1960s, one between species and
the other within a population of a species, that directly contributed to the
development of MET. The interpretations of these observations as well as the
assumptions by MET are widely known to be inadequate. We review the recent
progress of an alternative framework, the maximum genetic diversity (MGD)
hypothesis, that uses axioms and natural selection to explain the vast majority
of genetic diversity as being at optimum equilibrium that is largely determined
by organismal complexity. The MGD hypothesis fully absorbs the proven virtues
of MET and considers its assumptions relevant only to a much more limited
scope. This new synthesis has accounted for the much overlooked phenomenon of
progression towards higher complexity, and more importantly, been instrumental
in directing productive research into both evolutionary and biomedical
problems.


  The hypothesis that living systems can benefit from operating at the vicinity
of critical points has gained momentum in recent years. Criticality may confer
an optimal balance between exceedingly ordered and too noisy states. We here
present a model, based on information theory and statistical mechanics,
illustrating how and why a community of agents aimed at understanding and
communicating with each other converges to a globally coherent state in which
all individuals are close to an internal critical state, i.e. at the borderline
between order and disorder. We study --both analytically and computationally--
the circumstances under which criticality is the best possible outcome of the
dynamical process, confirming the convergence to critical points under very
generic conditions. Finally, we analyze the effect of cooperation (agents try
to enhance not only their fitness, but also that of other individuals) and
competition (agents try to improve their own fitness and to diminish those of
competitors) within our setting. The conclusion is that, while competition
fosters criticality, cooperation hinders it and can lead to more ordered or
more disordered consensual solutions.


  Bacteriocins are peptide-derived molecules produced by bacteria, whose
recently-discovered functions include virulence factors and signalling
molecules as well as their better known roles as antibiotics. To date, close to
five hundred bacteriocins have been identified and classified. Recent
discoveries have shown that bacteriocins are highly diverse and widely
distributed among bacterial species. Given the heterogeneity of bacteriocin
compounds, many tools struggle with identifying novel bacteriocins due to their
vast sequence and structural diversity. Many bacteriocins undergo
post-translational processing or modifications necessary for the biosynthesis
of the final mature form. Enzymatic modification of bacteriocins as well as
their export is achieved by proteins whose genes are often located in a
discrete gene cluster proximal to the bacteriocin precursor gene, referred to
as \textit{context genes} in this study. Although bacteriocins themselves are
structurally diverse, context genes have been shown to be largely conserved
across unrelated species. Using this knowledge, we set out to identify new
candidates for context genes which may clarify how bacteriocins are
synthesized, and identify new candidates for bacteriocins that bear no sequence
similarity to known toxins. To achieve these goals, we have developed a
software tool, Bacteriocin Operon and gene block Associator (BOA) that can
identify homologous bacteriocin associated gene clusters and predict novel
ones. We discover that several phyla have a strong preference for bactericon
genes, suggesting distinct functions for this group of molecules. Availability:
https://github.com/idoerg/BOA


  Two-locus sampling probabilities have played a central role in devising an
efficient composite likelihood method for estimating fine-scale recombination
rates. Due to mathematical and computational challenges, these sampling
probabilities are typically computed under the unrealistic assumption of a
constant population size, and simulation studies have shown that resulting
recombination rate estimates can be severely biased in certain cases of
historical population size changes. To alleviate this problem, we develop here
new methods to compute the sampling probability for variable population size
functions that are piecewise constant. Our main theoretical result, implemented
in a new software package called LDpop, is a novel formula for the sampling
probability that can be evaluated by numerically exponentiating a large but
sparse matrix. This formula can handle moderate sample sizes ($n \leq 50$) and
demographic size histories with a large number of epochs ($\mathcal{D} \geq
64$). In addition, LDpop implements an approximate formula for the sampling
probability that is reasonably accurate and scales to hundreds in sample size
($n \geq 256$). Finally, LDpop includes an importance sampler for the posterior
distribution of two-locus genealogies, based on a new result for the optimal
proposal distribution in the variable-size setting. Using our methods, we study
how a sharp population bottleneck followed by rapid growth affects the
correlation between partially linked sites. Then, through an extensive
simulation study, we show that accounting for population size changes under
such a demographic model leads to substantial improvements in fine-scale
recombination rate estimation. LDpop is freely available for download at
https://github.com/popgenmethods/ldpop


  Single Particle Tracking (SPT) can aid in understanding complex
spatio-temporal processes. However, quantifying diffusivity and forces from
individual live cell trajectories is complicated by inter- & intra-trajectory
kinetic heterogeneity, thermal fluctuations, and statistical temporal
dependence inherent to the underlying molecule's time correlated confined
dynamics experienced in the cell. Experimental artifacts such as localization
uncertainty and motion blur also obscure the data. We introduce a new maximum
likelihood estimation (MLE) technique that decouples the above noise sources
and systematically treats temporal correlation via a likelihood function
(permitting more reliable extraction of effective forces from position vs. time
data). Our estimator is demonstrated to be consistent over a wide range of
exposure times, diffusion coefficients, and confinement "radii". The algorithm
and corresponding software can reliably extract motion parameters independent
of exposure time in trajectories exhibiting confined and/or non-stationary
dynamics and will aid in directly comparing trajectories obtained from
different imaging modalities.


  We study coagulation-fragmentation equations inspired by a simple model
proposed in fisheries science to explain data for the size distribution of
schools of pelagic fish. Although the equations lack detailed balance and admit
no $H$-theorem, we are able to develop a rather complete description of
equilibrium profiles and large-time behavior, based on recent developments in
complex function theory for Bernstein and Pick functions. In the
large-population continuum limit, a scaling-invariant regime is reached in
which all equilibria are determined by a single scaling profile. This universal
profile exhibits power-law behavior crossing over from exponent $-\frac23$ for
small size to $-\frac32$ for large size, with an exponential cut-off.


  The CRISPR (clustered regularly interspaced short palindromic repeats)
mechanism allows bacteria to adaptively defend against phages by acquiring
short genomic sequences (spacers) that target specific sequences in the viral
genome. We propose a population dynamical model where immunity can be both
acquired and lost. The model predicts regimes where bacterial and phage
populations can co-exist, others where the populations exhibit damped
oscillations, and still others where one population is driven to extinction.
Our model considers two key parameters: (1) ease of acquisition and (2) spacer
effectiveness in conferring immunity. Analytical calculations and numerical
simulations show that if spacers differ mainly in ease of acquisition, or if
the probability of acquiring them is sufficiently high, bacteria develop a
diverse population of spacers. On the other hand, if spacers differ mainly in
their effectiveness, their final distribution will be highly peaked, akin to a
"winner-take-all" scenario, leading to a specialized spacer distribution.
Bacteria can interpolate between these limiting behaviors by actively tuning
their overall acquisition probability.


  Cooperative behaviors are common in humans and are fundamental to our
society. Theoretical and experimental studies have modeled environments in
which the behaviors of humans, or agents, have been restricted to analyze their
social behavior. However, it is important that such studies are generalized to
less restrictive environments to understand human society. Social network games
(SNGs) provide a particularly powerful tool for the quantitative study of human
behavior. In SNGs, numerous players can behave more freely than in the
environments used in previous studies; moreover, their relationships include
apparent conflicts of interest and every action can be recorded. We focused on
reciprocal altruism, one of the mechanisms that generate cooperative behavior.
This study aims to investigate cooperative behavior based on reciprocal
altruism in a less restrictive environment. For this purpose, we analyzed the
social behavior underlying such cooperative behavior in an SNG. We focused on a
game scenario in which the relationship between the players was similar to that
in the Leader game. We defined cooperative behaviors by constructing a payoff
matrix in the scenario. The results showed that players maintained cooperative
behavior based on reciprocal altruism, and cooperators received more advantages
than noncooperators. We found that players constructed reciprocal relationships
based on two types of interactions, cooperative behavior and unproductive
communication.


  Although significant insight has been gained into biochemical, genetic and
structural features of oxidative phosphorylation (OXPHOS) at the single-enzyme
level, relatively little was known of how the component complexes function
together in time and space until recently. Several pioneering single-molecule
studies have emerged over the last decade in particular, which have illuminated
our knowledge of OXPHOS, most especially on model bacterial systems. Here, we
discuss these recent findings of bacterial OXPHOS, many of which generate
time-resolved information of the OXPHOS machinery with the native physiological
context intact. These new investigations are transforming our knowledge of not
only the molecular arrangement of OXPHOS components in live bacteria, but also
of the way components dynamically interact with each other in a functional
state. These new discoveries have important implications towards putative
supercomplex formation in bacterial OXPHOS in particular.


  We present a quantum system incorporating qualitative aspects of enzyme
action in which the possibility of quantum superposition of several
conformations of the enzyme-substrate complex is investigated. We present
numerical results showing quantum effects that transcend the case of a
statistical mixture of conformations.


  Influenza viruses enter a cell via endocytosis after binding to the surface.
During the endosomal journey, acidification triggers a conformational change of
the virus spike protein hemagglutinin (HA) that results in escape of the viral
genome from the endosome to the cytoplasm. A quantitative understanding of the
processes involved in HA mediated fusion with the endosome is still missing. We
develop here a stochastic model to estimate the change of conformation of HAs
inside the endosome nanodomain. Using a Markov-jump process to model the
arrival of protons to HA binding sites, we compute the kinetics of their
accumulation and the mean first time for HAs to be activated. This analysis
reveals that HA proton binding sites possess a high chemical barrier, ensuring
a stability of the spike protein at sub-acidic pH. Finally, we predict that
activating more than 3 adjacent HAs is necessary to prevent a premature fusion.


  Cellular internalization of nanoparticles requires the full wrapping of the
nanoparticles by the cell membrane. This wrapping process can occur
spontaneously if the adhesive interactions between the nanoparticles and the
membranes are sufficiently strong to compensate for the cost of membrane
bending. In this article, we show that the membrane curvature prior to wrapping
plays a key role for the wrapping process, besides the size and shape of the
nanoparticles that have been investigated in recent years. For membrane
segments that initially bulge away from nanoparticles by having a mean
curvature of the same sign as the mean curvature of the particle surface, we
find strongly stable partially wrapped states that can prevent full wrapping.
For membrane segments that initially bulge towards the nanoparticles, in
contrast, partially wrapped states can constitute a significant energetic
barrier for the wrapping process.


  Owing to the complexity and variability of metagenomic studies, modern
machine learning approaches have seen increased usage to answer a variety of
question encompassing the full range of metagenomic NGS data analysis. We
review here the contribution of machine learning techniques for the field of
metagenomics, by presenting known successful approaches in a unified framework.
This review focuses on five important metagenomic problems: OTU-clustering,
binning, taxonomic profling and assignment, comparative metagenomics and gene
prediction. For each of these problems, we identify the most prominent methods,
summarize the machine learning approaches used and put them into perspective of
similar methods. We conclude our review looking further ahead at the challenge
posed by the analysis of interactions within microbial communities and
different environments, in a field one could call "integrative metagenomics".


  Leaving the joint enterprise when defection is unveiled is always a viable
option to avoid being exploited. Although loner strategy helps the population
not to be trapped into the tragedy of the commons state, it could offer only a
modest income for non-participants. In this paper we demonstrate that showing
some tolerance toward defectors could not only save cooperation in harsh
environments, but in fact results in a surprisingly high average payoff for
group members in public goods games. Phase diagrams and the underlying spatial
patterns reveal the high complexity of evolving states where cyclic dominant
strategies or two-strategy alliances can characterize the final state of
evolution. We identify microscopic mechanisms which are responsible for the
superiority of global solutions containing tolerant players. This phenomenon is
robust and can be observed both in well-mixed and in structured populations
highlighting the importance of tolerance in our everyday life.


  We investigate the dependence of the site frequency spectrum (SFS) on the
topological structure of genealogical trees. We show that basic population
genetic statistics - for instance estimators of $\theta$ or neutrality tests
such as Tajima's $D$ - can be decomposed into components of waiting times
between coalescent events and of tree topology. Our results clarify the
relative impact of the two components on these statistics. We provide a
rigorous interpretation of positive or negative values of an important class of
neutrality tests in terms of the underlying tree shape. In particular, we show
that values of Tajima's $D$ and Fay and Wu's $H$ depend in a direct way on a
peculiar measure of tree balance which is mostly determined by the root balance
of the tree. We present a new test for selection in the same class as Fay and
Wu's $H$ and discuss its interpretation and power. Finally, we determine the
trees corresponding to extreme expected values of these neutrality tests and
present formulae for these extreme values as a function of sample size and
number of segregating sites.


  The direct-coupling analysis is a powerful method for protein contact
prediction, and enables us to extract "direct" correlations between distant
sites that are latent in "indirect" correlations observed in a protein
multiple-sequence alignment. I show that the direct correlation can be obtained
by using a formulation analogous to the Ornstein-Zernike integral equation in
liquid theory. This formulation intuitively illustrates how the indirect or
apparent correlation arises from an infinite series of direct correlations, and
provides interesting insights into protein structure prediction.


  Community involvement and the preventive behavior of households are
considered to be at the heart of vector-control strategies. In this work, we
consider a simple theoretical model that enables us to take into account human
behaviors that may interfere with vector control. The model reflects the
trade-off between perceived costs and observed efficacy. Our theoretical
results emphasize that households may reduce their protective behavior in
response to mechanical elimination techniques piloted by a public agent,
leading to an increase of the total number of mosquitoes in the surrounding
environment and generating a barrier for vector-borne diseases control. Our
study is sufficiently generic to be applied to different arboviral diseases. It
also shows that vector-control models and strategies have to take into account
human behaviors.


  Bacteria pose unique challenges for genome-wide association studies (GWAS)
because of strong structuring into distinct strains and substantial linkage
disequilibrium across the genome. While methods developed for human studies can
correct for strain structure, this risks considerable loss- of-power because
genetic differences between strains often contribute substantial phenotypic
variability. Here we propose a new method that captures lineage-level
associations even when locus-specific associations cannot be fine-mapped. We
demonstrate its ability to detect genes and genetic variants underlying
resistance to 17 antimicrobials in 3144 isolates from four taxonomically
diverse clonal and recombining bacteria: Mycobacterium tuberculosis,
Staphylococcus aureus, Escherichia coli and Klebsiella pneumoniae. Strong
selection, recombination and penetrance confer high power to recover known
antimicrobial resistance mechanisms, and reveal a candidate association between
the outer membrane porin nmpC and cefazolin resistance in E. coli. Hence our
method pinpoints locus-specific effects where possible, and boosts power by
detecting lineage-level differences when fine-mapping is intractable.


  Systems composed of large numbers of interacting agents often admit an
effective coarse-grained description in terms of a multidimensional stochastic
dynamical system, driven by small-amplitude intrinsic noise. In applications to
biological, ecological, chemical and social dynamics it is common for these
models to posses quantities that are approximately conserved on short
timescales, in which case system trajectories are observed to remain close to
some lower-dimensional subspace. Here, we derive explicit and general formulae
for a reduced-dimension description of such processes that is exact in the
limit of small noise and well-separated slow and fast dynamics. The
Michaelis-Menten law of enzyme-catalyzed reactions, and the link between the
Lotka-Voltera and Wright-Fisher processes are explored as a simple worked
examples. Extensions of the method are presented for infinite dimensional
systems and processes coupled to non-Gaussian noise sources.


  In theoretical studies, the most commonly used measure of ecological
stability is resilience: ecosystems asymptotic rate of return to equilibrium
after a pulse-perturbation $-$or shock. A complementary notion of growing
popularity is reactivity: the strongest initial response to shocks. On the
other hand, empirical stability is often quantified as the inverse of temporal
variability, directly estimated on data, and reflecting ecosystems response to
persistent and erratic environmental disturbances. It is unclear whether and
how this empirical measure is related to resilience and reactivity. Here, we
establish a connection by introducing two variability-based stability measures
belonging to the theoretical realm of resilience and reactivity. We call them
intrinsic, stochastic and deterministic invariability; respectively defined as
the inverse of the strongest stationary response to white-noise and to
single-frequency perturbations. We prove that they predict ecosystems worst
response to broad classes of disturbances, including realistic models of
environmental fluctuations. We show that they are intermediate measures between
resilience and reactivity and that, although defined with respect to persistent
perturbations, they can be related to the whole transient regime following a
shock, making them more integrative notions than reactivity and resilience. We
argue that invariability measures constitute a stepping stone, and discuss the
challenges ahead to further unify theoretical and empirical approaches to
stability.


  In daily life, subjects often face a social dilemma in two stages. In Stage
1, they recognize the social dilemma structure of the decision problem at hand
(a tension between personal interest and collective interest); in Stage 2, they
have to choose between gathering additional information to learn the exact
payoffs corresponding to each of the two options or making a choice without
looking at the payoffs. While previous theoretical research suggests that the
mere act of considering one's strategic options in a social dilemma will be met
with distrust, no experimental study has tested this hypothesis. What does
"looking at payoffs" signal in observers? Do observers' beliefs actually match
decision makers' intentions? Experiment 1 shows that the actual action of
looking at payoffs signals selfish behavior, but it does not actually mean so.
Experiments 2 and 3 show that, when the action of looking at payoffs is
replaced by a self-report question asking the extent to which participants look
at payoffs in their everyday lives, subjects in high looking mode are indeed
more selfish than those in low looking mode, and this is correctly predicted by
observers. These results support Rand and colleagues' Social Heuristics
Hypothesis and the novel "cooperate without looking" model by Yoeli, Hoffman,
and Nowak. However, Experiment 1 shows that actual looking may lead to
different results, possibly caused by the emergence of a moral cleansing
effect.


  An inverse problem framework for constructing reaction systems with
prescribed properties is presented. Kinetic transformations are defined and
analysed as a part of the framework, allowing an arbitrary polynomial ordinary
differential equation to be mapped to the one that can be represented as a
reaction network. The framework is used for construction of specific two- and
three-dimensional bistable reaction systems undergoing a supercritical
homoclinic bifurcation, and the topology of their phase spaces is discussed.


  Recent studies has revealed that Microtubules (MTs) exhibit three transition
states of growth, shrinkage and pause. In this paper, we first introduce a
three states random evolution model as a framework for studying MTs dynamics in
three transition states of growth, pause and shrinkage. Then, we introduce a
non-traditional stack run encoding scheme with 5 symbols for detecting
transition states as well as to encode MT experimental data. The peak detection
is carried out in the wavelet domain to effectively detect these three
transition states. One of the added advantages of including peak information
while encoding being that it enables to detect the peaks efficiently and
encodes them simultaneously in the wavelet domain without having the need to do
further processing after the decoding stage. Experimental results show that
using this form of non-traditional stack run encoding has better compression
and reconstruction performance as opposed to traditional stack run encoding and
run length encoding schemes. Parameters for MTs modeled in the three states are
estimated and is shown to closely approximate original MT data for lower
compression rates. As the compression rate increases, we may end up throwing
away details that are required to detect transition states of MTs. Thus,
choosing the right compression rate is a trade-off between admissible level of
error in signal reconstruction, its parameter estimation and considerable rate
of compression of MT data.


  An interesting conundrum in biological control questions the efficiency of
generalist predators as biological control agents. Theory suggests, generalist
predators are poor agents for biological control, primarily due to mutual
interference. However field evidence shows they are actually quite effective in
regulating pest densities. In this work we provide a plausible answer to this
paradox. We analyze a three species model, where a generalist top predator is
introduced into an ecosystem as a biological control, to check the population
of a middle predator, that in turn is depredating on a prey species. We show
that the inclusion of predator interference alone, can cause the solution of
the top predator equation to blow-up in finite time, while there is global
existence in the no interference case. This result shows that interference
could actually cause a population explosion of the top predator, enabling it to
control the target species, thus corroborating recent field evidence. Our
results might also partially explain the population explosion of certain
species, introduced originally for biological control purposes, such as the
cane toad (\emph{Bufo marinus}) in Australia, which now functions as a
generalist top predator. We also show both Turing instability and
spatio-temporal chaos in the model. Lastly we investigate time delay effects.


  To verify the Global Burden of Disease Study and the provided healthy life
expectancy (HALE) estimates from the World Health Organization (WHO) we propose
a very simple model based on the mortality {\mu}x of a population provided in a
classical life table and a mortality diagram. We use the abridged life tables
provided by WHO. Our estimates are compared with the HALE estimates for the
World territories and the WHO countries. Even more we have developed the
related simple program in Excel which provides immediately the Life Expectancy,
the Loss of Healthy Life Years and the Healthy Life Expectancy estimate. We
also apply the health state function theory to have more estimates and
comparisons. The results suggest improved WHO estimates in recent years for the
majority of the cases. Keywords: Health state function, Healthy life
expectancy, Mortality Diagram, Loss of healthy years, LHLY, HALE, DALE, World
Health Organization, WHO, Global burden of Disease, Health status.


  We explore connections between metagenomic read assignment and the
quantification of transcripts from RNA-Seq data. In particular, we show that
the recent idea of pseudoalignment introduced in the RNA-Seq context is
suitable in the metagenomics setting. When coupled with the
Expectation-Maximization (EM) algorithm, reads can be assigned far more
accurately and quickly than is currently possible with state of the art
software.


  We formulate a two-patch mathematical model for Ebola Virus Disease dynamics
in order to evaluate the effectiveness of \textit{cordons sanitaires},
mandatory movement restrictions between communities while exploring their role
on disease dynamics and final epidemic size. Simulations show that severe
restrictions in movement between high and low risk areas of closely linked
communities may have a deleterious impact on the overall levels of infection in
the total population.


  Cells use temporal dynamical patterns to transmit information via signaling
pathways. As optimality with respect to the environment plays a fundamental
role in biological systems, organisms have evolved optimal ways to transmit
information. Here, we use optimal control theory to obtain the dynamical signal
patterns for the optimal transmission of information, in terms of efficiency
(low energy) and reliability (low uncertainty). Adopting an
activation-deactivation decoding network, we reproduce several dynamical
patterns found in actual signals, such as steep, gradual, and overshooting
dynamics. Notably, when minimizing the energy of the input signal, the optimal
signals exhibit overshooting, which is a biphasic pattern with transient and
steady phases; this pattern is prevalent in actual dynamical patterns. We also
identify conditions in which these three patterns (steep, gradual, and
overshooting) confer advantages. Our study shows that cellular signal
transduction is governed by the principle of minimizing free energy dissipation
and uncertainty; these constraints serve as selective pressures when designing
dynamical signaling patterns.


  We introduce exclusion models of two distinguishable species of hard rods
with their distinct sites of entry and exit under open boundary conditions. In
the first model both species of rods move in the same direction whereas in the
other two models they move in the opposite direction. These models are
motivated by the biological phenomenon known as Transcriptional Interference.
Therefore, the rules for the kinetics of the models, particularly the rules for
the outcome of the encounter of the rods, are also formulated to mimic those
observed in Transcriptional Interference. By a combination of mean-field theory
and computer simulation of these models we demonstrate how the flux of one
species of rods is completely switched off by the other. Exploring the
parameter space of the model we also establish the conditions under which
switch-like regulation of two fluxes is possible; from the extensive analysis
we discover more than one possible mechanism of this phenomenon.


  The degree to which population fluctuations arise from variable adult
survival relative to variable recruitment has been debated widely for marine
organisms. Disentangling these effects remains challenging because data
generally are not sufficient to evaluate if and how adult survival rates are
regulated by stochasticity and/or population density. Using unique time-series
for a largely unexploited reef fish, we found both population density and
stochastic food supply impacted adult survival. The estimated effect of
variable survival on adult abundance (both mean and variability) rivaled that
of variable recruitment. Moreover, we show density dependent adult survival can
dampen impacts of stochastic recruitment. Thus, food variability may alter
population fluctuations by simultaneously regulating recruitment and
compensatory adult survival. These results provide an additional mechanism for
why intensified density independent mortality (via harvest or other means)
amplifies population fluctuations and emphasizes need for research evaluating
the causes and consequences of variability in adult survival.


  Gaussian network model(GNM) and anisotropic network model(ANM) are some of
the most popular methods for the study of protein flexibility and related
functions. In this work, we propose generalized GNM(gGNM) and ANM methods and
show that the GNM Kirchhoff matrix can be built from the ideal low-pass filter,
which is a special case of a wide class of correlation functions underpinning
the linear scaling flexibility-rigidity index(FRI) method. Based on the
mathematical structure of correlation functions, we propose a unified framework
to construct generalized Kirchhoff matrices whose matrix inverse leads to
gGNMs, whereas, the direct inverse of its diagonal elements gives rise to FRI
method.With this connection,we further introduce two multiscale elastic network
models, namely, multiscale GNM(mGNM) and multiscale ANM(mANM), which are able
to incorporate different scales into the generalized Kirchkoff matrices or
generalized Hessian matrices.We validate our new multiscale methods with
extensive numerical experiments. We illustrate that gGNMs outperform the
original GNM method in the B-factor prediction of a set of 364 proteins.We
demonstrate that for a given correlation function, FRI and gGNM methods provide
essentially identical B-factor predictions when the scale value in the
correlation function is sufficiently large.More importantly,we reveal intrinsic
multiscale behavior in protein structures. The proposed mGNM and mANM are able
to capture this multiscale behavior and thus give rise to a significant
improvement of more than 11% in B-factor predictions over the original GNM and
ANM methods. We further demonstrate benefit of our mGNM in the B-factor
predictions on many proteins that fail the original GNM method. We show that
the present mGNM can also be used to analyze protein domain separations.
Finally, we showcase the ability of our mANM for the simulation of protein
collective motions.


  Protein-nucleic acid complexes are important for many cellular processes
including the most essential function such as transcription and translation.
For many protein-nucleic acid complexes, flexibility of both macromolecules has
been shown to be critical for specificity and/or function. Flexibility-rigidity
index (FRI) has been proposed as an accurate and efficient approach for protein
flexibility analysis. In this work, we introduce FRI for the flexibility
analysis of protein-nucleic acid complexes. We demonstrate that a multiscale
strategy, which incorporates multiple kernels to capture various length scales
in biomolecular collective motions, is able to significantly improve the state
of art in the flexibility analysis of protein-nucleic acid complexes. We take
the advantage of the high accuracy and ${\cal O}(N)$ computational complexity
of our multiscale FRI method to investigate the flexibility of large ribosomal
subunits, which is difficult to analyze by alternative approaches. An
anisotropic FRI approach, which involves localized Hessian matrices, is
utilized to study the translocation dynamics in an RNA polymerase.


  Multiple dynamic pathways always exist in biological networks, but their
robustness against internal fluctuations and relative stability have not been
well recognized and carefully analyzed yet. Here we try to address these issues
through an illustrative example, namely the Siah-1/beta-catenin/p14/19 ARF loop
of protein p53 dynamics. Its deterministic Boolean network model predicts that
two parallel pathways with comparable magnitudes of attractive basins should
exist after the protein p53 is activated when a cell becomes harmfully
disturbed. Once the low but non-neglectable intrinsic fluctuations are
incorporated into the model, we show that a phase transition phenomenon is
emerged: in one parameter region the probability weights of the normal pathway,
reported in experimental literature, are comparable with the other pathway
which is seemingly abnormal with the unknown functions, whereas, in some other
parameter regions, the probability weight of the abnormal pathway can even
dominate and become globally attractive. The theory of exponentially perturbed
Markov chains is applied and further generalized in order to quantitatively
explain such a phase transition phenomenon, in which the nonequilibrium
"activation energy barriers" along each transiting trajectory between the
parallel pathways and the number of "optimal transition paths" play a central
part. Our theory can also determine how the transition time and the number of
optimal transition paths between the parallel pathways depend on each
interaction's strength, and help to identify those possibly more crucial
interactions in the biological network.


  Under prolonged starvation, the nucleoid of E. coli bacteria forms
nanocrystalline complexes with the protein Dps. By one hypothesis, this effect
is considered to be a manifestation of a general strategy of self-preservation
of living organisms under adverse conditions by forming stable ordered
structures which do not require consumption of energy for maintaining the
order. The effect of nucleoid biocrystallization in starved bacteria has been
investigated in a number of studies using electron microscopy. However, this
method did not answer all questions regarding the structure of the
nanocrystals, and, primarily, the question of the true conformation of DNA in
these complexes. This issue is particularly important because it is related to
a more general, actively researched topic of organization of the genetic
material, which affects gene expression and is thus critical for the
functioning of organisms. The intracellular DNA-Dps nanocrystals have the size
in the range of hundreds of nanometers and are thus particularly suitable for
study at X-ray free electron lasers (XFELs). In this paper, we discuss the
possible configurations of XFEL experiments aimed at determining the
conformation of DNA in the DNA-Dps nanocrystals grown in vivo and in vitro, as
well as the specific features and challenges associated with this type of
sample. We conclude that the XFEL technique has a high potential for uncovering
the structural details of bacterial nucleoid restructuring under stress.


  The epidemic threshold of a spreading process indicates the condition for the
occurrence of the wide spreading regime, thus representing a predictor of the
network vulnerability to the epidemic. Such threshold depends on the natural
history of the disease and on the pattern of contacts of the network with its
time variation. Based on the theoretical framework introduced in (Valdano et
al. PRX 2015) for a susceptible-infectious-susceptible model, we formulate here
an infection propagator approach to compute the epidemic threshold accounting
for more realistic effects regarding a varying force of infection per contact,
the presence of immunity, and a limited time resolution of the temporal
network. We apply the approach to two temporal network models and an empirical
dataset of school contacts. We find that permanent or temporary immunity do not
affect the estimation of the epidemic threshold through the infection
propagator approach. Comparisons with numerical results show the good agreement
of the analytical predictions. Aggregating the temporal network rapidly
deteriorates the predictions, except for slow diseases once the heterogeneity
of the links is preserved. Weight-topology correlations are found to be the
critical factor to be preserved to improve accuracy in the prediction.


  We solve a class of optimization problems for (phylogenetic) $X$-trees or
their shapes. These problems have recently appeared in different contexts, e.g.
in the context of the impact of tree shapes on the size of TBR neighborhoods,
but so far these problems have not been characterized and solved in a
systematic way. In this work we generalize the concept and also present several
applications. Moreover, our results give rise to a nice notion of balance for
trees. Unsurprisingly, so-called caterpillars are the most unbalanced tree
shapes, but it turns out that balanced tree shapes cannot be described so
easily as they need not even be unique.