  Binding of acridine derivative quinacrine (QA) to chicken erythrocyte DNA was
studied by methods of absorption and polarized fluorescent spectroscopy.
Measurements were carried out in aqueous buffered solutions (pH 6.9) of
different dye concentrations (QA concentration range from $10^{-6}$ till
$10^{-4}$ M) and ionic strengths ($Na^{+}$ concentration rang from $10^{-3}$
till 0.15 M) in a wide range of phosphate-to-dye molar ratios ($P/D$). It is
established that the minimum of fluorescent titration curve plotted as relative
fluorescence intensity $vs$ $P/D$ is conditioned by the competition between the
two types of QA binding to DNA which posses by different emission parameters:
(i) intercalative one dominating under high $P/D$ values, and (ii) outside
electrostatic binding dominating under low $P/D$ values, which is accompanied
by the formation of non-fluorescent dye associates on the DNA backbone.
Absorption and fluorescent characteristics of complexes formed were determined.
The method of calculation of different binding modes contribution to the
complex formation depending on $P/D$ value is presented. It was shown that the
size of binding site measured as the number of DNA base pairs per one QA
molecule bound in the case of the electrostatic interaction is 8 times less
than that for the intercalative one that determines the competitive ability of
the outside binding against the stronger intercalative binding mode.


  Mathematical models of epidemic dynamics offer significant insight into
predicting and controlling infectious diseases. The dynamics of a disease model
generally follow a susceptible, infected, and recovered (SIR) model, with some
standard modifications. In this paper, we extend the work of Munz et.al (2009)
on the application of disease dynamics to the so-called "zombie apocalypse",
and then apply the identical methods to influenza dynamics. Unlike Munz et.al
(2009), we include data taken from specific depictions of zombies in popular
culture films and apply Markov Chain Monte Carlo (MCMC) methods on improved
dynamical representations of the system. To demonstrate the usefulness of this
approach, beyond the entertaining example, we apply the identical methodology
to Google Trend data on influenza to establish infection and recovery rates.
Finally, we discuss the use of the methods to explore hypothetical intervention
policies regarding disease outbreaks.


  Because of the huge number of graphs possible even with a small number of
nodes, inference on network structure is known to be a challenging problem.
Generating large random directed graphs with prescribed probabilities of
occurrences of some meaningful patterns (motifs) is also difficult. We show how
to generate such random graphs according to a formal probabilistic
representation, using fast Markov chain Monte Carlo methods to sample them. As
an illustration, we generate realistic graphs with several hundred nodes
mimicking a gene transcription interaction network in Escherichia coli.


  Hepatitis C virus helicase is a molecular motor that splits nucleic acid
duplex structures during viral replication, therefore representing a promising
target for antiviral treatment. Hence, a detailed understanding of the
mechanism by which it operates would facilitate the development of efficient
drug-assisted therapies aiming to inhibit helicase activity. Despite extensive
investigations performed in the past, a thorough understanding of the activity
of this important protein was lacking since the underlying internal
conformational motions could not be resolved. Here we review investigations
that have been previously performed by us for HCV helicase. Using methods of
structure-based computational modelling it became possible to follow entire
operation cycles of this motor protein in structurally resolved simulations and
uncover the mechanism by which it moves along the nucleic acid and accomplishes
strand separation. We also discuss observations from that study in the light of
recent experimental studies that confirm our findings.


  The ratio of non-synonymous to synonymous substitutions
$\omega(=d_{N}/d_{S})$ has been widely used as a measure of adaptive evolution
in protein coding genes. Omega can be defined in terms of population genetics
parameters as the fixation ratio of selected vs. neutral mutants. Here it is
argued that approaches based on the infinite sites model are not appropriate to
define $\omega$ for single codon locations. Simple models of amino acid
substitution with reversible mutation and selection are analysed, and used to
define $\omega$ under several evolutionary scenarios. In most practical cases
$\omega<1$ when selection is constant throughout time. However, it is shown
that when the pattern of selection on amino acids changes, for example after an
environment shift, a temporary burst of adaptive evolution ($\omega\gg1$) can
be observed. The fixation probability of a novel mutant under frequency
dependent selection is calculated, and it is used to show why $\omega>1$ can be
sometimes expected for single locations at equilibrium. An example with
influenza data is discussed.


  The emergence of agriculture is suggested to have driven extensive human
population growths. However, genetic evidence from maternal mitochondrial
genomes suggests major population expansions began before the emergence of
agriculture. Therefore, role of agriculture that played in initial population
expansions still remains controversial. Here, we analyzed a set of globally
distributed whole Y chromosome and mitochondrial genomes of 526 male samples
from 1000 Genome Project. We found that most major paternal lineage expansions
coalesced in Neolithic Time. The estimated effective population sizes through
time revealed strong evidence for 10- to 100-fold increase in population growth
of males with the advent of agriculture. This sex-biased Neolithic expansion
might result from the reduction in hunting-related mortality of males.


  Boolean networks have been used successfully in modeling biological networks
and provide a good framework for theoretical analysis. However, the analysis of
large networks is not trivial. In order to simplify the analysis of such
networks, several model reduction algorithms have been proposed; however, it is
not clear if such algorithms scale well with respect to the number of nodes.
The goal of this paper is to propose and implement an algorithm for the
reduction of AND-NOT network models for the purpose of steady state
computation. Our method of network reduction is the use of "steady state
approximations" that do not change the number of steady states. Our algorithm
is designed to work at the wiring diagram level without the need to evaluate or
simplify Boolean functions. Also, our implementation of the algorithm takes
advantage of the sparsity typical of discrete models of biological systems. The
main features of our algorithm are that it works at the wiring diagram level,
it runs in polynomial time, and it preserves the number of steady states. We
used our results to study AND-NOT network models of gene networks and showed
that our algorithm greatly simplifies steady state analysis. Furthermore, our
algorithm can handle sparse AND-NOT networks with up to 1000000 nodes.


  A biomimetic minimalist model membrane was used to study the mechanism and
kinetics of cell-free in vitro HIV-1 Gag budding from a giant unilamellar
vesicle (GUV). Real time interaction of Gag, RNA and lipid leading to the
formation of mini-vesicles was measured using confocal microscopy. Gag forms
resolution limited punctae on the GUV lipid membrane. Introduction of the Gag
and urea to a GUV solution containing RNA led to the budding of mini-vesicles
on the inside surface of the GUV. The GUV diameter showed a linear decrease in
time due to bud formation. Both bud formation and decrease in GUV size were
proportional to Gag concentration. In the absence of RNA, addition of urea to
GUVs incubated with Gag also resulted in subvesicle formation but exterior to
the surface. These observations suggest the possibility that clustering of GAG
proteins leads to membrane invagination even in the absence of host cell
proteins. The method presented here is promising, and allows for systematic
study of the dynamics of assembly of immature HIV and help classify the
hierarchy of factors that impact the Gag protein initiated assembly of
retroviruses such as HIV.


  The mechanism allowing a protein to search of a target sequence on DNA is
currently described as an intermittent process composed of 3D diffusion in bulk
and 1D diffusion along the DNA molecule. Due to the relevant charge of protein
and DNA, electrostatic interaction should play a crucial role during this
search. In this paper, we explicitly derive the mean field theory allowing for
a description of the protein-DNA electrostatics in solution. This approach
leads to a unified model of the search process, where 1D and 3D diffusion
appear as a natural consequence of the diffusion on an extended interaction
energy profile.


  Despite the greater functional importance of protein levels, our knowledge of
gene expression evolution is based almost entirely on studies of mRNA levels.
In contrast, our understanding of how translational regulation evolves has
lagged far behind. Here we have applied ribosome profiling - which measures
both global mRNA levels and their translation rates - to two species of
Saccharomyces yeast and their interspecific hybrid in order to assess the
relative contributions of changes in mRNA abundance and translation to
regulatory evolution. We report that both cis and trans-acting regulatory
divergence in translation are abundant, affecting at least 35% of genes. The
majority of translational divergence acts to buffer changes in mRNA abundance,
suggesting a widespread role for stabilizing selection acting across regulatory
levels. Nevertheless, we observe evidence of lineage-specific selection acting
on a number of yeast functional modules, including instances of reinforcing
selection acting at both levels of regulation. Finally, we also uncover
multiple instances of stop-codon readthrough that are conserved between
species. Our analysis reveals the under-appreciated complexity of
post-transcriptional regulatory divergence and indicates that partitioning the
search for the locus of selection into the binary categories of 'coding' vs.
'regulatory' may overlook a significant source of selection, acting at multiple
regulatory levels along the path from genotype to phenotype.


  Enhancers play a central role in cell-type-specific gene expression and are
marked by H3K4me1/2. Active enhancers are further marked by H3K27ac. However,
the methyltransferases responsible for H3K4me1/2 on enhancers remain elusive.
Furthermore, how these enzymes function on enhancers to regulate
cell-type-specific gene expression is unclear. Here we identify MLL4 (KMT2D) as
a major mammalian H3K4 mono- and di-methyltransferase with partial functional
redundancy with MLL3 (KMT2C). Using adipogenesis and myogenesis as model
systems, we show that MLL4 exhibits cell-type- and
differentiation-stage-specific genomic binding and is predominantly localized
on enhancers. MLL4 co-localizes with lineage-determining transcription factors
(TFs) on active enhancers during differentiation. Deletion of MLL4 markedly
decreases H3K4me1/2, H3K27ac, Polymerase II and Mediator levels on enhancers
and leads to severe defects in cell-type-specific gene expression and cell
differentiation. Together, these findings identify MLL4 as a major mammalian
H3K4 mono- and di-methyltransferase essential for enhancer activation during
cell differentiation.


  In the past one hundred years, deterministic rate equations have been
successfully used to infer enzyme-catalysed reaction mechanisms and to estimate
rate constants from reaction kinetics experiments conducted in vitro. In recent
years, sophisticated experimental techniques have been developed that allow the
measurement of enzyme- catalysed and other biopolymer-mediated reactions inside
single cells at the single molecule level. Time course data obtained by these
methods are considerably noisy because molecule numbers within cells are
typically quite small. As a consequence, the interpretation and analysis of
single cell data requires stochastic methods, rather than deterministic rate
equations. Here we concisely review both experimental and theoretical
techniques which enable single molecule analysis with particular emphasis on
the major developments in the field of theoretical stochastic enzyme kinetics,
from its inception in the mid-twentieth century to its modern day status. We
discuss the differences between stochastic and deterministic rate equation
models, how these depend on enzyme molecule numbers and substrate inflow into
the reaction compartment and how estimation of rate constants from single cell
data is possible using recently developed stochastic approaches.


  Gene expression and regulation rely on an apparently finely tuned set of
reactions between some proteins and DNA. Such DNA-binding proteins have to find
specific sequences on very long DNA molecules and they mostly do so in absence
of any active process. It has been rapidly recognized that to achieve this task
these proteins should be efficient at both searching (i.e. sampling fast
relevant parts of DNA) and finding (i.e. recognizing the specific site). A
two-mode search and variants of it have been suggested since the 70s to explain
either a fast search or an efficient recognition. Combining these two
properties at a phenomenological level is however more difficult as they appear
to have antagonist roles. To overcome this difficulty, one may simply need to
drop the dihcotomic view inherent to the two-mode search and look more
thoroughly at the set of interactions between DNA-binding proteins and a given
DNA segment either specific or non-specific. This chapter demonstrates that, in
doing so in a very generic way, one may indeed find a potential reconciliation
between a fast search and an efficient recognition. Although a lot remains to
be done, this could be the time for a change of paradigm


  The statistical properties of intra-chromosome maps obtained by a genome-wide
chromosome conformation capture method (Hi-C) are described in the framework of
the hierarchical crumpling model of heteropolymer chain with quenched disorder
in the primary sequence. We conjecture that the observed Hi-C maps are
statistical averages over many different ways of hierarchical genome folding,
and show that the existence of quenched primary structure coupled with
hierarchical folding can induce the observed fine structure of intra-chromosome
contact maps.


  Prediction of the functions of riboswitches requires a quantitative
description of the folding landscape so that the barriers and time scales for
the conformational change in the switching region in the aptamer can be
estimated. Using a combination of all atom molecular dynamics and
coarse-grained model simulations we studied the response of adenine (A) binding
add and pbuE A-riboswitches to mechanical force. The two riboswitches contain a
structurally similar three-way junction formed by three paired helices, P1, P2,
and P3, but carry out different functions. Using pulling simulations, with
structures generated in MD simulations, we show that after P1 rips the dominant
unfolding pathway in add A-riboswitch is the rupture of P2 followed by
unraveling of P3. In the pbuE A-riboswitch, after P1 unfolds P3 ruptures ahead
of P2. The order of unfolding of the helices, which is in accord with single
molecule pulling experiments, is determined by the relative stabilities of the
individual helices. Our results show that the stability of isolated helices
determines the order of assembly and response to force in these non-coding
regions. We use the simulated free energy profile for pbuE A-riboswitch to
estimate the time scale for allosteric switching, which shows that this
riboswitch is under kinetic control lending additional support to the
conclusion based on single molecule pulling experiments. A consequence of the
stability hypothesis is that a single point mutation (U28C) in the P2 helix of
the add A-riboswitch, which increases the stability of P2, would make the
folding landscapes of the two riboswitches similar. This prediction can be
tested in single molecule pulling experiments.


  The inheritance of characteristics induced by the environment has often been
opposed to the theory of evolution by natural selection. Yet, while evolution
by natural selection requires new heritable traits to be produced and
transmitted, it does not prescribe, per se, the mechanisms by which this is
operated. The mechanisms of inheritance are not, however, unconstrained, since
they are themselves subject to natural selection. We introduce a general,
analytically solvable mathematical model to compare the adaptive value of
different schemes of inheritance. Our model allows for variations to be
inherited, randomly produced, or environmentally induced, and, irrespectively,
to be either transmitted or not during reproduction. The adaptation of the
different schemes for processing variations is quantified for a range of
fluctuating environments, following an approach that links quantitative
genetics with stochastic control theory.


  Reconstructing transcriptional regulatory networks is an important task in
functional genomics. Data obtained from experiments that perturb genes by
knockouts or RNA interference contain useful information for addressing this
reconstruction problem. However, such data can be limited in size and/or are
expensive to acquire. On the other hand, observational data of the organism in
steady state (e.g. wild-type) are more readily available, but their
informational content is inadequate for the task at hand. We develop a
computational approach to appropriately utilize both data sources for
estimating a regulatory network. The proposed approach is based on a three-step
algorithm to estimate the underlying directed but cyclic network, that uses as
input both perturbation screens and steady state gene expression data. In the
first step, the algorithm determines causal orderings of the genes that are
consistent with the perturbation data, by combining an exhaustive search method
with a fast heuristic that in turn couples a Monte Carlo technique with a fast
search algorithm. In the second step, for each obtained causal ordering, a
regulatory network is estimated using a penalized likelihood based method,
while in the third step a consensus network is constructed from the highest
scored ones. Extensive computational experiments show that the algorithm
performs well in reconstructing the underlying network and clearly outperforms
competing approaches that rely only on a single data source. Further, it is
established that the algorithm produces a consistent estimate of the regulatory
network.


  Over the last two decades, several methods have been proposed for stabilizing
the dynamics of biological populations. However, these methods have typically
been evaluated using different population dynamics models and in the context of
very different concepts of stability, which makes it difficult to compare their
relative efficiencies. Moreover, since the dynamics of populations are
dependent on the life-history of the species and its environment, it is
conceivable that the stabilizing effects of control methods would also be
affected by such factors, a complication that has typically not been
investigated. In this study we compare six different control methods with
respect to their efficiency at inducing a common level of enhancement (50%
increase) for two kinds of stability (constancy and persistence) under four
different life history/ environment combinations. Since these methods have been
analytically studied elsewhere, we focus on an intuitive understanding of
realistic simulations incorporating noise, extinction probability and lattice
effect. We show that for these six methods, even when the magnitude of
stabilization attained is the same, other aspects of the dynamics like
population size distribution can be very different. Consequently, correlated
aspects of stability, like the amount of persistence for a given degree of
constancy stability (and vice versa) or the corresponding effective population
size (a measure of resistance to genetic drift) vary widely among the methods.
Moreover, the number of organisms needed to be added or removed to attain
similar levels of stabilization also varies for these methods, a fact that has
economic implications. Finally, we compare the relative efficiency of these
methods through a composite index of various stability related measures. We
find that restocking to a constant lower threshold seems to be the optimal
method under most conditions.


  Predicting adaptive evolutionary trajectories is a primary goal of
evolutionary biology. One can differentiate between forward and backward
predictability, where forward predictability measures the likelihood of the
same adaptive trajectory occurring in independent evolutions and backward
predictability measures the likelihood of a particular adaptive path given the
knowledge of starting and final states. Recent studies have attempted to
measure both forward and backward predictability using experimental evolution
in asexual haploid microorganisms. Similar experiments in diploid organisms
have not been conducted. Here we simulate adaptive walks using Fisher's
Geometric Model in haploids and diploids and find that adaptive walks in
diploids are less forward- and more backward-predictable than adaptive walks in
haploids. We argue that the difference is due to the ability of diploids in our
simulations to generate transiently stable polymorphisms and to allow adaptive
mutations of larger phenotypic effect. As stable polymorphisms can be generated
in both haploid and diploid natural populations through a number of mechanisms,
we argue that inferences based on experiments in which adaptive walks proceed
through succession of monomorphic states might miss many of the key features of
adaptation.


  The standard approach to analyzing 16S tag sequence data, which relies on
clustering reads by sequence similarity into Operational Taxonomic Units
(OTUs), underexploits the accuracy of modern sequencing technology. We present
a clustering-free approach to multi-sample Illumina datasets that can identify
independent bacterial subpopulations regardless of the similarity of their 16S
tag sequences. Using published data from a longitudinal time-series study of
human tongue microbiota, we are able to resolve within standard 97% similarity
OTUs up to 20 distinct subpopulations, all ecologically distinct but with 16S
tags differing by as little as 1 nucleotide (99.2% similarity). A comparative
analysis of oral communities of two cohabiting individuals reveals that most
such subpopulations are shared between the two communities at 100% sequence
identity, and that dynamical similarity between subpopulations in one host is
strongly predictive of dynamical similarity between the same subpopulations in
the other host. Our method can also be applied to samples collected in
cross-sectional studies and can be used with the 454 sequencing platform. We
discuss how the sub-OTU resolution of our approach can provide new insight into
factors shaping community assembly.


  Social reinforcement and modular structure are two salient features observed
in the spreading of behavior through social contacts. In order to investigate
the interplay between these two features, we study the generalized epidemic
process on modular networks with equal-sized finite communities and adjustable
modularity. Using the analytical approach originally applied to clique-based
random networks, we show that the system exhibits a bond-percolation type
continuous phase transition for weak social reinforcement, whereas a
discontinuous phase transition occurs for sufficiently strong social
reinforcement. Our findings are numerically verified using the finite-size
scaling analysis and the crossings of the bimodality coefficient.


  Under the hypothesis that both influenza and respiratory syncytial virus
(RSV) are the two leading causes of acute respiratory infections (ARI), in this
paper we have used a standard two-pathogen epidemic model as a regressor to
explain, on a yearly basis, high season ARI data in terms of the contact rates
and initial conditions of the mathematical model. The rationale is that ARI
high season is a transient regime of a noisy system, e.g., the system is driven
away from equilibrium every year by fluctuations in variables such as humidity,
temperature, viral mutations and human behavior. Using the value of the
replacement number as a phenotypic trait associated to fitness, we provide
evidence that influenza and RSV coexists throughout the ARI high season through
superinfection.


  Evolutionary trajectories and phenotypic states available to cell populations
are ultimately dictated by intermolecular interactions between DNA, RNA,
proteins, and other molecular species. Here we study how evolution of gene
regulation in a single-cell eukaryote S. cerevisiae is affected by the
interactions between transcription factors (TFs) and their cognate genomic
sites. Our study is informed by high-throughput in vitro measurements of TF-DNA
binding interactions and by a comprehensive collection of genomic binding
sites. Using an evolutionary model for monomorphic populations evolving on a
fitness landscape, we infer fitness as a function of TF-DNA binding energy for
a collection of 12 yeast TFs, and show that the shape of the predicted fitness
functions is in broad agreement with a simple thermodynamic model of two-state
TF-DNA binding. However, the effective temperature of the model is not always
equal to the physical temperature, indicating selection pressures in addition
to biophysical constraints caused by TF-DNA interactions. We find little
statistical support for the fitness landscape in which each position in the
binding site evolves independently, showing that epistasis is common in
evolution of gene regulation. Finally, by correlating TF-DNA binding energies
with biological properties of the sites or the genes they regulate, we are able
to rule out several scenarios of site-specific selection, under which binding
sites of the same TF would experience a spectrum of selection pressures
depending on their position in the genome. These findings argue for the
existence of universal fitness landscapes which shape evolution of all sites
for a given TF, and whose properties are determined in part by the physics of
protein-DNA interactions.


  We study a susceptible-vaccinated--infected--recovered (SVIR)
epidemic-spreading model with diversity of infection rate of the individuals.
By means of analytical arguments as well as extensive computer simulations, we
demonstrate that the heterogeneity in infection rate can either impede or
accelerate the epidemic spreading, which depends on the amount of vaccinated
individuals introduced in the population as well as the contact pattern among
the individuals. Remarkably, as long as the individuals with different
capability of acquiring the disease interact with unequal frequency, there
always exist a cross point for the fraction of vaccinated, below which the
diversity of infection rate hinders the epidemic spreading and above which
expedites it. The overall results are robust to the SVIR dynamics defined on
different population models; the possible applications of the results are
discussed.


  We analyze protein-protein interaction networks for six different species
under the framework of random matrix theory. Nearest neighbor spacing
distribution of the eigenvalues of adjacency matrices of the largest connected
part of these networks emulate universal Gaussian orthogonal statistics of
random matrix theory. We demonstrate that spectral rigidity, which quantifies
long range correlations in eigenvalues, for all protein-protein interaction
networks follow random matrix prediction up to certain ranges indicating
randomness in interactions. After this range, deviation from the universality
evinces underlying structural features in network.


  This chapter aims to inform a practitioner about current methods for
predicting potential distributions of invasive species. It mostly addresses
single species models, covering the conceptual bases, touching on mechanistic
models, and then focusing on methods using species distribution records and
environmental data to predict distributions. The commentary in this last
section is oriented towards key issues that arise in fitting, and predicting
with, these models (which include CLIMEX, MaxEnt and other regression methods).
In other words, it is more about the process of thinking about the data and the
modelling problem (which is a challenging one) than it is about one technique
versus another. The discussion helps clarify the necessary steps and expertise
for predicting distributions. Some researchers are optimistic that correlative
models will predict with high precision; while that may be true for some
species at some scales of evaluation, I believe that the issues discussed in
this chapter show that substantial errors are reasonably likely. I am hopeful
that ongoing developments will produce models better suited to the task and
tools to help practitioners to better understand predictions and their
uncertainties.


  We theoretically investigate exciton relaxation dynamics in molecular
aggregates based on model photosynthetic complexes under various conditions of
incoherent excitation. We show that noise-induced quantum coherence is
generated between spatially-separated exciton states which belong to the same
or different energy transfer pathways, coupled via real and virtual transfer
processes. Such quantum coherence effects may be used to improve
light-harvesting efficiency and to reveal quantum phenomena in biology.


  Biomolecules are the prime information processing elements of living matter.
Most of these inanimate systems are polymers that compute their structures and
dynamics using as input seemingly random character strings of their sequence,
following which they coalesce and perform integrated cellular functions. In
large computational systems with a finite interaction-codes, the appearance of
conflicting goals is inevitable. Simple conflicting forces can lead to quite
complex structures and behaviors, leading to the concept of "frustration" in
condensed matter. We present here some basic ideas about frustration in
biomolecules and how the frustration concept leads to a better appreciation of
many aspects of the architecture of biomolecules, and how structure connects to
function. These ideas are simultaneously both seductively simple and perilously
subtle to grasp completely. The energy landscape theory of protein folding
provides a framework for quantifying frustration in large systems and has been
implemented at many levels of description. We first review the notion of
frustration from the areas of abstract logic and its uses in simple condensed
matter systems. We discuss then how the frustration concept applies
specifically to heteropolymers, testing folding landscape theory in computer
simulations of protein models and in experimentally accessible systems.
Studying the aspects of frustration averaged over many proteins provides ways
to infer energy functions useful for reliable structure prediction. We discuss
how frustration affects folding, how a large part of the biological functions
of proteins are related to subtle local frustration effects and how frustration
influences the appearance of metastable states, the nature of binding
processes, catalysis and allosteric transitions. We hope to illustrate how
Frustration is a fundamental concept in relating function to structural
biology.


  X-ray data indicate that complexes of ribosomes with cognate and near cognate
tRNAs are very similar structurally, and this was the ground for a suggestion
that the ribosome discriminates correct codon-anticodon pair because of its
higher stability. Here an alternative explanation of kinetic proofreading is
suggested, and intramolecular oscillations in tRNAs play a keystone role in it.
Resonance of the oscillations allows the cognate codon-anticodon pair to be
conserved due to fast energy transfer to other part of the tRNA molecule. This
mechanism can potentially discriminate correct pair from an incorrect one even
if they have similar stabilities.


  Two-dimensional electrophoresis has nurtured the birth of proteomics. It is
however no longer the exclusive setup used in proteomics, with the development
of shotgun proteomics techniques that appear more fancy and fashionable
nowadays.Nevertheless, 2D gel-based proteomics still has valuable features, and
sometimes unique ones, which make it often an attractive choice when a
proteomics strategy must be selected. These features are detailed in this
chapter, as is the rationale for selecting or not 2D gel-based proteomics as a
proteomic strategy.


  The macroscopic properties, the properties of individual components and how
those components interact with each other are three important aspects of a
composited structure. An understanding of the interplay between them is
essential in the study of complex systems. Using axonal cytoskeleton as an
example system, here we perform a theoretical study of slender structures that
can be coarse-grained as a simple smooth 3-dimensional curve. We first present
a generic model for such systems based on the fundamental theorem of curves. We
use this generic model to demonstrate the applicability of the well-known
worm-like chain (WLC) model to the network level and investigate the situation
when the system is stretched by strong forces (weakly bending limit). We
specifically studied recent experimental observations that revealed the
hitherto unknown periodic cytoskeleton structure of axons and measured the
longitudinal fluctuations. Instead of focusing on single molecules, we apply
analytical results from the WLC model to both single molecule and network
levels and focus on the relations between extensions and fluctuations. We show
how this approach introduces constraints to possible local dynamics of the
spectrin tetramers in the axonal cytoskeleton and finally suggests simple but
self-consistent dynamics of spectrins in which the spectrins in one spatial
period of axons fluctuate in-sync.


  In this short note we prove that, given two (not necessarily binary) rooted
phylogenetic trees T_1, T_2 on the same set of taxa X, where |X|=n, the
hybridization number of T_1 and T_2 can be computed in time O^{*}(2^n) i.e.
O(2^{n} poly(n)). The result also means that a Maximum Acyclic Agreement Forest
(MAAF) can be computed within the same time bound.


  Semiflexible polymers characterized by the contour length $L$ and persistent
length $\ell_p$ confined in a spatial region $D$ have been described as a
series of ``{\em spherical blobs}'' and ``{\em deflecting lines}'' by de Gennes
and Odjik for $\ell_p < D$ and $\ell_p \gg D$ respectively. Recently new
intermediate regimes ({\em extended de Gennes} and {\em Gauss-de Gennes}) have
been investigated by Tree {\em et al.} [Phys. Rev. Lett. {\bf 110}, 208103
(2013)]. In this letter we derive scaling relations to characterize these
transitions in terms of universal scaled fluctuations in $d$-dimension as a
function of $L,\ell_p$, and $D$, and show that the Gauss-de Gennes regime is
absent and extended de Gennes regime is vanishingly small for polymers confined
in a 2D strip. We validate our claim by extensive Brownian dynamics (BD)
simulation which also reveals that the prefactor $A$ used to describe the chain
extension in the Odjik limit is independent of physical dimension $d$ and is
the same as previously found by Yang {\em et al.}[Y. Yang, T. W. Burkhardt, G.
Gompper, Phys. Rev. E {\bf 76}, 011804 (2007)]. Our studies are relevant for
optical maps of DNA stretched inside a nano-strip.


  Bacterial flagellar motor (BFM) is one of the ion-driven molecular machines,
which drives the rotation of flagellar filaments and enable bacteria to swim in
viscous solutions. Understanding its mechanism is one challenge in biophysics.
Based on previous models and inspired by the idea used in description of motor
proteins, in this study one two-state model is provided. Meanwhile, according
to corresponding experimental data, mathematical relationship between BFM
membrane voltage and pH value of the environment, and relationship between
internal and external sodium concentrations are given. Therefore, with model
parameter values obtained by fitting theoretical results of torque-speed
relation to recent experimental data, many biophysical properties of bacterial
flagellar motor can be obtained for any pH values and any external sodium
concentrations. Including the rotation speed, stall torque (i.e. the torque
generated by BFM), rotation dispersion, and rotation randomness. In this study,
the single-stator BFM will be firstly analyzed, and then properties of
multiple-stator BFM are addressed briefly.


  In a molecular communication network, transmitters and receivers communicate
by using signalling molecules. At the receivers, the signalling molecules
react, via a chain of chemical reactions, to produce output molecules. The
counts of output molecules over time is considered to be the output signal of
the receiver. This output signal is used to detect the presence of signalling
molecules at the receiver. The output signal is noisy due to the stochastic
nature of diffusion and chemical reactions. The aim of this paper is to
characterise the properties of the output signals for two types of receivers,
which are based on two different types of reaction mechanisms. We derive
analytical expressions for the mean, variance and frequency properties of these
two types of receivers. These expressions allow us to study the properties of
these two types of receivers. In addition, our model allows us to study the
effect of the diffusibility of the receiver membrane on the performance of the
receivers.


  About 2% of human genetic polymorphisms have been hypothesized to arise via
multinucleotide mutations (MNMs), complex events that generate SNPs at multiple
sites in a single generation. MNMs have the potential to accelerate the pace at
which single genes evolve and to confound studies of demography and selection
that assume all SNPs arise independently. In this paper, we examine clustered
mutations that are segregating in a set of 1,092 human genomes, demonstrating
that MNMs become enriched as large numbers of individuals are sampled. We
leverage the size of the dataset to deduce new information about the allelic
spectrum of MNMs, estimating the percentage of linked SNP pairs that were
generated by simultaneous mutation as a function of the distance between the
affected sites and showing that MNMs exhibit a high percentage of transversions
relative to transitions. These findings are reproducible in data from multiple
sequencing platforms. Among tandem mutations that occur simultaneously at
adjacent sites, we find an especially skewed distribution of ancestral and
derived dinucleotides, with $\textrm{GC}\to \textrm{AA}$, $\textrm{GA}\to
\textrm{TT}$ and their reverse complements making up 36% of the total. These
same mutations dominate the spectrum of tandem mutations produced by the
upregulation of low-fidelity Polymerase $\zeta$ in mutator strains of S.
cerevisiae that have impaired DNA excision repair machinery. This suggests that
low-fidelity DNA replication by Pol $\zeta$ is at least partly responsible for
the MNMs that are segregating in the human population, and that useful
information about the biochemistry of MNM can be extracted from ordinary
population genomic data. We incorporate our findings into a mathematical model
of the multinucleotide mutation process that can be used to correct
phylogenetic and population genetic methods for the presence of MNMs.


  We present a computational model to study the spatiotemporal dynamics of the
epidermis homeostasis under normal and pathological conditions. The model
consists of a population kinetics model of the central transition pathway of
keratinocyte proliferation, differentiation and loss and an agent-based model
that propagates cell movements and generates the stratified epidermis. The
model recapitulates observed homeostatic cell density distribution, the
epidermal turnover time and the multilayered tissue structure. We extend the
model to study the onset, recurrence and phototherapy-induced remission of
psoriasis. The model considers the psoriasis as a parallel homeostasis of
normal and psoriatic keratinocytes originated from a shared stem-cell niche
environment and predicts two homeostatic modes of the psoriasis: a disease mode
and a quiescent mode. Interconversion between the two modes can be controlled
by interactions between psoriatic stem cells and the immune system and by the
normal and psoriatic stem cells competing for growth niches. The prediction of
a quiescent state potentially explains the efficacy of the multi-episode UVB
irradiation therapy and recurrence of psoriasis plaques, which can further
guide designs of therapeutics that specifically target the immune system and/or
the keratinocytes.


  The formula for the probability of fixation of a new mutation is widely used
in theoretical population genetics and molecular evolution. Here we derive a
series of identities, inequalities and approximations for the exact probability
of fixation of a new mutation under the Moran process (equivalent results hold
for the approximate probability of fixation for the Wright-Fisher process after
an appropriate change of variables). We show that the behavior of the logarithm
of the probability of fixation is particularly simple when the selection
coefficient is measured as a difference of Malthusian fitnesses, and we exploit
this simplicity to derive several inequalities and approximations. We also
present a comprehensive comparison of both existing and new approximations for
the probability of fixation, highlighting in particular approximations that
result in a reversible Markov chain when used to model the dynamics of
evolution under weak mutation.


  Understanding the patterns and causes of phenotypic divergence is a central
goal in evolutionary biology. Much work has shown that mRNA abundance is highly
variable between closely related species. However, the extent and mechanisms of
post-transcriptional gene regulatory evolution are largely unknown. Here we
used ribosome profiling to compare transcript abundance and translation
efficiency in two closely related yeast species (S. cerevisiae and S.
paradoxus). By comparing translation regulatory divergence to interspecies
differences in mRNA sequence features, we show that differences in transcript
leaders and codon bias substantially contribute to divergent translation.
Globally, we find that translation regulatory divergence often buffers
species-differences in mRNA abundance, such that ribosome occupancy is more
conserved than transcript abundance. We used allele-specific ribosome profiling
in interspecies hybrids to compare the relative contributions of cis- and
trans-regulatory divergence to species differences in mRNA abundance and
translation efficiency. The mode of gene regulatory divergence differs for
these processes, as trans-regulatory changes play a greater role in divergent
mRNA abundance than in divergent translation efficiency. Strikingly, most genes
with aberrant transcript abundance in F1 hybrids (either over- or
under-expressed compared to both parent species) did not exhibit aberrant
ribosome occupancy. Our results show that interspecies differences in
translation contribute substantially to the evolution of gene expression.
Compensatory differences in transcript abundance and translation efficiency may
increase the robustness of gene regulation.


  In this work we investigate the development of stable dynamical structures
along interfaces separating domains belonging to enemy partnerships, in the
context of cyclic predator-prey models with an even number of species $N \ge
8$. We use both stochastic and field theory simulations in one and two spatial
dimensions, as well as analytical arguments, to describe the association at the
interfaces of mutually neutral individuals belonging to enemy partnerships and
to probe their role in the development of the dynamical structures at the
interfaces. We identify an interesting behaviour associated to the symmetric or
asymmetric evolution of the interface profiles depending on whether $N/2$ is
odd or even, respectively. We also show that the macroscopic evolution of the
interface network is not very sensitive internal structure of the interfaces.
Although this work focus on cyclic predator prey-models with an even number of
species, we argue that the results are expected to be quite generic in the
context of spatial stochastic May-Leonard models.


  Motor enzymes are remarkable molecular machines that use the energy derived
from the hydrolysis of a nucleoside triphosphate to generate mechanical
movement, achieved through different steps that constitute their kinetic cycle.
These macromolecules, nowadays investigated with advanced experimental
techniques to unveil their molecular mechanisms and the properties of their
kinetic cycles, are implicated in many biological processes, ranging from
biopolymerisation (e.g. RNA polymerases and ribosomes) to intracellular
transport (motor proteins such as kinesins or dyneins). Although the kinetics
of individual motors is well studied on both theoretical and experimental
grounds, the repercussions of their stepping cycle on the collective dynamics
still remains unclear. Advances in this direction will improve our
comprehension of transport process in the natural intracellular medium, where
processive motor enzymes might operate in crowded conditions. In this work, we
therefore extend the current statistical kinetic analysis to study collective
transport phenomena of motors in terms of lattice gas models belonging to the
exclusion process class. Via numerical simulations, we show how to interpret
and use the randomness calculated from single particle trajectories in crowded
conditions. Importantly, we also show that time fluctuations and non-Poissonian
behavior are intrinsically related to spatial correlations and the emergence of
large, but finite, clusters of co-moving motors. The properties unveiled by our
analysis have important biological implications on the collective transport
characteristics of processive motor enzymes in crowded conditions.


  We show that, if truth assignments on $n$ variables reproduce through
recombination so that satisfaction of a particular Boolean function confers a
small evolutionary advantage, then a polynomially large population over
polynomially many generations (polynomial in $n$ and the inverse of the initial
satisfaction probability) will end up almost certainly consisting exclusively
of satisfying truth assignments. We argue that this theorem sheds light on the
problem of novelty in Evolution.


  Modern population genetics studies typically involve genome-wide genotyping
of individuals from a diverse network of ancestries. An important, unsolved
problem is how to formulate and estimate probabilistic models of observed
genotypes that allow for complex population structure. We formulate two general
probabilistic models, and we propose computationally efficient algorithms to
estimate them. First, we show how principal component analysis (PCA) can be
utilized to estimate a general model that includes the well-known
Pritchard-Stephens-Donnelly mixed-membership model as a special case. Noting
some drawbacks of this approach, we introduce a new "logistic factor analysis"
(LFA) framework that seeks to directly model the logit transformation of
probabilities underlying observed genotypes in terms of latent variables that
capture population structure. We demonstrate these advances on data from the
human genome diversity panel and 1000 genomes project, where we are able to
identify SNPs that are highly differentiated with respect to structure while
making minimal modeling assumptions.


  We consider an age-structured epidemic model with two basic public health
interventions: (i) identifying and isolating symptomatic cases, and (ii)
tracing and quarantine of the contacts of identified infectives. The dynamics
of the infected population are modeled by a nonlinear infection-age-dependent
partial differential equation, which is coupled with an ordinary differential
equation that describes the dynamics of the susceptible population. Theoretical
results about global existence and uniqueness of positive solutions are proved.
We also present two practical applications of our model: (1) we assess public
health guidelines about emergency preparedness and response in the event of a
smallpox bioterrorist attack; (2) we simulate the 2003 SARS outbreak in Taiwan
and estimate the number of cases avoided by contact tracing. Our model can be
applied as a rational basis for decision makers to guide interventions and
deploy public health resources in future epidemics.


  Microtubules (MTs) are cytoplasmic protein polymers that are essential for
fundamental cellular processes including the maintenance of cell shape,
organelle transport and formation of the mitotic spindle. Microtubule dynamic
instability is critical for these processes, but it remains poorly understood,
in part because the relationship between the structure of the MT tip and the
growth/depolymerization transitions is enigmatic. What are the functionally
significant aspects of a tip structure that is capable of promoting MT growth,
and how do changes in these characteristics cause the transition to
depolymerization (catastrophe)? Here we use computational models to investigate
the connection between cracks (laterally unbonded regions) between
protofilaments and dynamic instability. Our work indicates that it is not the
depth of the cracks per se that governs MT dynamic instability. Instead it is
whether the cracks terminate in GTP-rich or GDP-rich areas of the MT that
governs whether a particular MT tip structure is likely to grow, shrink, or
transition: the identity of the crack-terminating subunit pairs has a profound
influence on the moment-by-moment behavior of the MT.


  We use ideas from kinetic proofreading, an error-correcting mechanism in
biology, to identify new kinetic regimes in non-equilibrium systems. These
regimes are defined by the sensitivity of the occupancy of a state of the
system to a change in its energy. In biological contexts, higher sensitivity
corresponds to stronger discrimination between molecular substrates with
different energetics competing in the same reaction. We study this
discriminatory ability in systems with discrete states that are connected by a
general network of transitions. We find multiple regimes of different
discriminatory ability when the energy of a given state of the network is
varied. Interestingly, the occupancy of the state can even increase with its
energy, corresponding to an "anti-proofreading" regime. The number and
properties of such discriminatory regimes are limited by the topology of the
network. Finally, we find that discriminatory regimes can be changed without
modifying any "hard-wired" structural aspects of the system but rather by
simply changing external chemical potentials.


  Recent years have seen a rapid expansion of the model space explored in
statistical phylogenetics, emphasizing the need for new approaches to
statistical model representation and software development. Clear communication
and representation of the chosen model is crucial for: (1) reproducibility of
an analysis, (2) model development and (3) software design. Moreover, a
unified, clear and understandable framework for model representation lowers the
barrier for beginners and non-specialists to grasp complex phylogenetic models,
including their assumptions and parameter/variable dependencies.
  Graphical modeling is a unifying framework that has gained in popularity in
the statistical literature in recent years. The core idea is to break complex
models into conditionally independent distributions. The strength lies in the
comprehensibility, flexibility, and adaptability of this formalism, and the
large body of computational work based on it. Graphical models are well-suited
to teach statistical models, to facilitate communication among phylogeneticists
and in the development of generic software for simulation and statistical
inference.
  Here, we provide an introduction to graphical models for phylogeneticists and
extend the standard graphical model representation to the realm of
phylogenetics. We introduce a new graphical model component, tree plates, to
capture the changing structure of the subgraph corresponding to a phylogenetic
tree. We describe a range of phylogenetic models using the graphical model
framework and introduce modules to simplify the representation of standard
components in large and complex models. Phylogenetic model graphs can be
readily used in simulation, maximum likelihood inference, and Bayesian
inference using, for example, Metropolis-Hastings or Gibbs sampling of the
posterior distribution.


  The homogeneous reconstructed evolutionary process is a birth-death process
without observed extinct lineages. Each species evolves independently with the
same diversification rates (speciation rate $\lambda(t)$ and extinction rate
$\mu(t)$) that may change over time. The process is commonly applied to model
species diversification where the data are reconstructed phylogenies, e.g.,
trees reconstructed from present-day molecular data, and used to infer
diversification rates.
  In the present paper I develop the general probability density of a
reconstructed tree under any time-dependent birth-death process. I elaborate on
how to adapt this probability density if conditioned on survival of one or two
initial lineages, or having sampled $n$ species and show how to transform
between the probability density of a reconstructed and the probability density
of the speciation times.
  I demonstrate the use of the general time-dependent probability density
functions by deriving the probability density of a reconstructed tree under a
birth-death-shift model with explicit mass-extinction events. I enrich this
compendium by providing and discussing several special cases, including: the
pure birth process, the pure death process, the birth-death process and the
critical branching process. Thus, I provide here most of the commonly used
birth-death models in a unified framework (e.g., same condition and same data)
with common notation.


  We give a complete closed form description of the evolution of human blood
genotypes frequencies (in the ABO and Rh classification) after any (finite or
infinite) number of generations and for any initial distribution.


  It is the main purpose of this paper to introduce a graph-valued stochastic
process in order to model the spread of a communicable infectious disease. The
major novelty of the SIR model we promote lies in the fact that the social
network on which the epidemics is taking place is not specified in advance but
evolves through time, accounting for the temporal evolution of the interactions
involving infective individuals. Without assuming the existence of a fixed
underlying network model, the stochastic process introduced describes, in a
flexible and realistic manner, epidemic spread in non-uniformly mixing and
possibly heterogeneous populations. It is shown how to fit such a
(parametrised) model by means of Approximate Bayesian Computation methods based
on graph-valued statistics. The concepts and statistical methods described in
this paper are finally applied to a real epidemic dataset, related to the
spread of HIV in Cuba in presence of a contact tracing system, which permits
one to reconstruct partly the evolution of the graph of sexual partners
diagnosed HIV positive between 1986 and 2006.


  We scanned through the genomes of 29,141 African Americans, searching for
loci where the average proportion of African ancestry deviates significantly
from the genome-wide average. We failed to find any genome-wide significant
deviations, and conclude that any selection in African Americans since
admixture is sufficiently weak that it falls below the threshold of our power
to detect it using a large sample size. These results stand in contrast to the
findings of a recent study of selection in African Americans. That study, which
had 15 times fewer samples, reported six loci with significant deviations. We
show that the discrepancy is likely due to insufficient correction for multiple
hypothesis testing in the previous study. The same study reported 14 loci that
showed greater population differentiation between African Americans and
Nigerian Yoruba than would be expected in the absence of natural selection.
Four such loci were previously shown to be genome-wide significant and likely
to be affected by selection, but we show that most of the 10 additional loci
are likely to be false positives. Additionally, the most parsimonious
explanation for the loci that have significant evidence of unusual
differentiation in frequency between Nigerians and Africans Americans is
selection in Africa prior to their forced migration to the Americas.


  We consider the Moran model on the sharp peak landscape, in the asymptotic
regime studied in [3], where a quasispecies is formed. We find explicitly the
distribution of this quasispecies.


  Objectives: Re-infection after treatment for Chlamydia trachomatis or
Neisseria gonorrhoeae reduces the effect of control interventions. We explored
the impact of delays in partner treatment on the expected probability of
re-infection of index cases using a mathematical model. Methods: We used
previously reported parameter distributions to calculate the probability that
index cases would be re-infected by their untreated partners. We then assumed
different delays between index case and partner treatment to calculate the
probabilities of re-infection. Results: In the absence of partner treatment,
the medians of the expected re-infection probabilities are 19.4% (interquartile
range (IQR) 9.2-31.6%) for chlamydia and 12.5% (IQR 5.6-22.2%) for gonorrhoea.
If all current partners receive treatment three days after the index case, the
expected re-infection probabilities are 4.2% (IQR 2.1-6.9%) for chlamydia and
5.5% (IQR 2.6-9.5%) for gonorrhoea. Conclusions: Quicker partner referral and
treatment can substantially reduce re-infection rates for chlamydia and
gonorrhoea. The formula we used to calculate re-infection rates can be used to
inform the design of randomised controlled trials of novel partner notification
technologies like accelerated partner therapy.


  Biological systems produce outputs in response to variable inputs.
Input-output relations tend to follow a few regular patterns. For example, many
chemical processes follow the S-shaped Hill equation relation between input
concentrations and output concentrations. That Hill equation pattern
contradicts the fundamental Michaelis-Menten theory of enzyme kinetics. I use
the discrepancy between the expected Michaelis-Menten process of enzyme
kinetics and the widely observed Hill equation pattern of biological systems to
explore the general properties of biological input-output relations. I start
with the various processes that could explain the discrepancy between basic
chemistry and biological pattern. I then expand the analysis to consider
broader aspects that shape biological input-output relations. Key aspects
include the input-output processing by component subsystems and how those
components combine to determine the system's overall input-output relations.
That aggregate structure often imposes strong regularity on underlying
disorder. Aggregation imposes order by dissipating information as it flows
through the components of a system. The dissipation of information may be
evaluated by the analysis of measurement and precision, explaining why certain
common scaling patterns arise so frequently in input-output relations. I
discuss how aggregation, measurement and scale provide a framework for
understanding the relations between pattern and process. The regularity imposed
by those broader structural aspects sets the contours of variation in biology.
Thus, biological design will also tend to follow those contours. Natural
selection may act primarily to modulate system properties within those broad
constraints.


  A reaction network consists of a finite number of species, which interact
through predefined reaction channels. Traditionally such networks were modeled
deterministically, but it is now well-established that when reactant copy
numbers are small, the random timing of the reactions create internal noise
that can significantly affect the macroscopic properties of the system. To
understand the role of noise and quantify its effects, stochastic models are
necessary. In the stochastic setting, the population is described by a
probability distribution, which evolves according to a set of ordinary
differential equations known as the Chemical Master Equation (CME). This set is
infinite in most cases making the CME practically unsolvable. In many
applications, it is important to determine if the solution of a CME has a
globally attracting fixed point. This property is called ergodicity and its
presence leads to several important insights about the underlying dynamics. The
goal of this paper is to present a simple procedure to verify ergodicity in
stochastic reaction networks. We provide a set of simple linear-algebraic
conditions which are sufficient for the network to be ergodic. In particular,
our main condition can be cast as a Linear Feasibility Problem (LFP) which is
essentially the problem of determining the existence of a vector satisfying
certain linear constraints. The inherent scalability of LFPs make our approach
efficient, even for very large networks. We illustrate our procedure through an
example from systems biology.


  Evolutionary graph theory is a well established framework for modelling the
evolution of social behaviours in structured populations. An emerging consensus
in this field is that graphs that exhibit heterogeneity in the number of
connections between individuals are more conducive to the spread of cooperative
behaviours. In this article we show that such a conclusion largely depends on
the individual-level interactions that take place. In particular, averaging
payoffs garnered through game interactions rather than accumulating the payoffs
can altogether remove the cooperative advantage of heterogeneous graphs while
such a difference does not affect the outcome on homogeneous structures. In
addition, the rate at which game interactions occur can alter the evolutionary
outcome. Less interactions allow heterogeneous graphs to support more
cooperation than homogeneous graphs, while higher rates of interactions make
homogeneous and heterogeneous graphs virtually indistinguishable in their
ability to support cooperation. Most importantly, we show that common measures
of evolutionary advantage used in homogeneous populations, such as a comparison
of the fixation probability of a rare mutant to that of the resident type, are
no longer valid in heterogeneous populations. Heterogeneity causes a bias in
where mutations occur in the population which affects the mutant's fixation
probability. We derive the appropriate measures for heterogeneous populations
that account for this bias.


  Protein contacts contain important information for protein structure and
functional study, but contact prediction from sequence remains very
challenging. Both evolutionary coupling (EC) analysis and supervised machine
learning methods are developed to predict contacts, making use of different
types of information, respectively. This paper presents a group graphical lasso
(GGL) method for contact prediction that integrates joint multi-family EC
analysis and supervised learning. Different from existing 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. To implement joint EC analysis, we model a set of related protein
families using Gaussian graphical models (GGM) and then co-estimate their
precision matrices by maximum-likelihood, subject to the constraint that the
precision matrices shall share similar residue co-evolution patterns. To
further improve the accuracy of the estimated precision matrices, we employ a
supervised learning method to predict contact probability from a variety of
evolutionary and non-evolutionary information and then incorporate the
predicted probability as prior into our GGL framework. Experiments show that
our method can predict contacts much more accurately than existing methods, and
that our method performs better on both conserved and family-specific contacts.


  Here we present the first genome wide statistical test for recessive
selection. This test uses explicitly non-equilibrium demographic differences
between populations to infer the mode of selection. By analyzing the transient
response to a population bottleneck and subsequent re-expansion, we
qualitatively distinguish between alleles under additive and recessive
selection. We analyze the response of the average number of deleterious
mutations per haploid individual and describe time dependence of this quantity.
We introduce a statistic, $B_R$, to compare the number of mutations in
different populations and detail its functional dependence on the strength of
selection and the intensity of the population bottleneck. This test can be used
to detect the predominant mode of selection on the genome wide or regional
level, as well as among a sufficiently large set of medically or functionally
relevant alleles.


  We present a new proof, using the log-sum inequality, that the
pseudo-Helmholtz free energy function is a Lyapunov function for
complex-balanced mass-action systems. This proof is shorter and simpler than
previous proofs.


  As a complementary tool to nanofluidics, biomolecular based transport is
envisioned for nanotechnological devices. We report a new method for guiding
microtubule shuttles on multi-walled carbon nanotube tracks, aligned by
dielectrophoresis on a functionalized surface. In the absence of electric field
and in fluid flow, alignment is maintained. The directed translocation of
kinesin propelled microtubules has been investigated using fluorescence
microscopy. To our knowledge, this is the first demonstration of microtubules
gliding along carbon nanotubes.


  Most animals have separate sexes. The differential expression of gene
products, in particular that of gene regulators, is underlying sexual
dimorphism. Analyses of sex-biased expression have focused mostly in protein
coding genes. Several lines of evidence indicate that microRNAs, a class of
major gene regulators, are likely to have a significant role in sexual
dimorphism. This role has not been systematically explored so far. Here I study
the sex-biased expression pattern of microRNAs in the model species Drosophila
melanogaster. As with protein coding genes, sex biased microRNAs are associated
with the reproductive function. Strikingly, contrary to protein-coding genes,
male biased microRNAs are enriched in the X chromosome whilst female microRNAs
are mostly autosomal. I propose that the chromosomal distribution is a
consequence of high rates of de novo emergence, and a preference of new
microRNAs to be expressed in the testis. I also suggest that demasculinization
of the X chromosome may not affect microRNAs. Interestingly, female biased
microRNAs are often encoded within protein coding genes that are also expressed
in females. MicroRNAs with sex-biased expression do not preferentially target
sex-biased gene transcripts. These results strongly suggest that the sex-biased
expression of microRNAs is mainly a consequence of high rates of microRNA
emergence in the X (male bias) or hitch-hiked expression by host genes (female
bias).


  DNA methylation (DNAme) is a critical component of the epigenetic regulatory
machinery and aberrations in DNAme patterns occur in many diseases, such as
cancer. Mapping and understanding DNAme profiles offers considerable promise
for reversing the aberrant states. There are several approaches to analyze
DNAme, which vary widely in cost, resolution and coverage. Affinity capture and
high-throughput sequencing of methylated DNA strike a good balance between the
high cost of whole genome bisulphite sequencing (WGBS) and the low coverage of
methylation arrays. However, existing methods cannot adequately differentiate
between hypomethylation patterns and low capture efficiency, and do not offer
flexibility to integrate copy number variation (CNV). Furthermore, no
uncertainty estimates are provided, which may prove useful for combining data
from multiple protocols or propagating into downstream analysis. We propose an
empirical Bayes framework that uses a fully methylated (i.e. SssI treated)
control sample to transform observed read densities into regional methylation
estimates. In our model, inefficient capture can be distinguished from low
methylation levels by means of larger posterior variances. Furthermore, we can
integrate CNV by introducing a multiplicative offset into our Poisson model
framework. Notably, our model offers analytic expressions for the mean and
variance of the methylation level and thus is fast to compute. Our algorithm
outperforms existing approaches in terms of bias, mean-squared error and
coverage probabilities as illustrated on multiple reference datasets. Although
our method provides advantages even without the SssI-control, considerable
improvement is achieved by its incorporation. Our method can be applied to
methylated DNA affinity enrichment assays (e.g MBD-seq, MeDIP-seq) and a
software implementation is available in the Bioconductor Repitools package.


  PLUMED-GUI is an interactive environment to develop and test complex PLUMED
scripts within the Visual Molecular Dynamics (VMD) environment. Computational
biophysicists can take advantage of both PLUMED's rich syntax to define
collective variables (CVs) and VMD's chemically-aware atom selection language,
while working within a natural point-and-click interface. Pre-defined templates
and syntax mnemonics facilitate the definition of well-known reaction
coordinates. Complex CVs, e.g. involving reference snapshots used for RMSD or
native contacts calculations, can be built through dialogs that provide a
synoptic view of the available options. Scripts can be either exported for use
in simulation programs, or evaluated on the currently loaded molecular
trajectories. Development of scripts takes place without leaving VMD, thus
enabling an incremental try-see-modify development model for molecular metrics.


  Several studies question the adaptive value of female preferences for older
males. Theory and evidence show that older males carry more deleterious
mutations in their sperm than younger males carry. These mutations are not
visible to females choosing mates. Germ-line mutations could oppose preferences
for "good genes." Choosy females run the risk that offspring of older males
will be no more attractive or healthy than offspring of younger males.
Germ-line mutations could pose a particular problem when females can only judge
male trait size, rather than assessing age directly. I ask whether or not
females will prefer extreme traits, despite reduced offspring survival due to
age-dependent mutation. I use a quantitative genetic model to examine the
evolution of female preferences, an age-dependent male trait, and overall
health ("condition"). My dynamical equation includes mutation bias that depends
on the generation time of the population. I focus on the case where females
form preferences for older males because male trait size depends on male age.
My findings agree with good genes theory. Females at equilibrium always select
above-average males. The trait size preferred by females directly correlates
with the direct costs of the preference. Direct costs can accentuate the
equilibrium preference at a higher rate than mutational parameters. Females can
always offset direct costs by mating with older, more ornamented males.
Age-dependent mutation in condition maintains genetic variation in condition
and thereby maintains the selective value of female preferences. Rather than
eliminating female preferences, germ-line mutations provide an essential
ingredient in sexual selection.


  We study the features on the unzipping process of a modified version of the
Peyrard-Bishop- Dauxois model. We show that the inclusion of a barrier in the
on-site potential allows for the existence of stable domain wall solutions
between open and closed regions of the DNA chain. We analyze the linear
stability of such solutions and study their relevance in the dynamical behavior
of DNA under ac forces.


  Single nucleotide polymorphisms (SNPs) represent an important type of dynamic
sites within the human genome. These common variants often locally correlate
into more complex multi-SNP haploblocks that are maintained throughout
generations in a stable population. The information encoded in the structure of
common SNPs and SNP haploblock variation can be characterized through a
normalized information content (NIC) metric. Such an intrinsic measure allows
disparate regions of individual genomes and the genomes of various populations
to be quantitatively compared in a meaningful way.
  Using our defined measures of genomic information, the interplay of
maintained statistical variations due to the environmental baths within which
stable populations exist can be interrogated. We develop the analogous
"thermodynamics" characterizing the state variables for genomic populations
that are stable under stochastic environmental stresses. Since living systems
have not been found to develop in the absence of environmental influences, we
focus on describing the analogous genomic free energy measures in this
development.
  The intensive parameter describing how an environment drives genomic
diversity is found to depend inversely upon the NIC of the genome of a stable
population within that environment. Once this environmental potential has been
determined from the whole genome of a population, additive state variables can
be directly related to the probabilities of the occurrence of given viable SNP
based units (alleles) within that genome. This formulation allows the
determination of both population averaged state variables as well as the
genomic energies of individual alleles and their combinations. The
determination of individual allelic potentials then should allow the
parameterization of specific environmental influences upon shared alleles
across populations in varying environments.


  Recent H5N1 influenza research has revived the debate on the storage and
manipulation of potentially harmful pathogens. In the last two decades, new
high biosafety (BSL-4) laboratories entered into operation, raising strong
concerns from the public. The probability of an accidental release of a
pathogen from a BSL-4 laboratory is extremely low, but the corresponding risk
-- defined as the probability of occurrence multiplied by its impact -- could
be significant depending on the pathogen specificities and the population
potentially affected. A list of BSL-4 laboratories throughout the world, with
their location and date of first activity, was established from publicly
available sources. This database was used to estimate the total population
living within a daily commuting distance of BSL-4 laboratories, and to quantify
how this figure changed over time. We show that from 1990 to present, the
population living within the commuting belt of BSL-4 laboratories increased by
a factor of 4 to reach up to 1.8% of the world population, owing to an increase
in the number of facilities and their installation in cities. Europe is
currently hosting the largest population living in the direct vicinity of BSL-4
laboratories, while the recent building of new facilities in Asia suggests that
an important increase of the population living close to BSL-4 laboratories will
be observed in the next decades. We discuss the potential implications in term
of global risk, and call for better pathogen-specific quantitative assessment
of the risk of outbreaks resulting from the accidental release of potentially
pandemic pathogens


  We metaphorically call "Grammar of Biology" a small field of genomic
research, whose main objective is to search for intrinsic DNA sequence
properties. Erwin Chargaff inaugurated it back in 50s, but since then little
progress has been made. It remained almost neglected until early 90s, when
Vinayakumar V. Prabhu made a major contribution determining the Symmetry
Principle. Remarkably, different sciences have contributed for its development,
for instance, Chargaff used his Chemistry background to discover the so-called
Chargaff's rules; taking advantage of several publicly available genomic
sequences, and through Computational and Statistical analyses, Prabhu
identified the Symmetry Principle and, recently, using a Mathematical approach,
we have discovered four new Generalized Chargaff's rules. Our work has expanded
the "Grammar of Biology", and created the conceptual and theoretical framework
necessary to further developments.


  The effect of disturbance on a model ecosystem of sessile and mutually
competitive species [Mathiesen et al. Phys. Rev. Lett. 107, 188101 (2011);
Mitarai et al. Phys. Rev. E 86, 011929 (2012) ] is studied. The disturbance
stochastically removes individuals from the system, and the created empty sites
are re-colonized by neighbouring species. We show that the stable
high-diversity state, maintained by occasional cyclic species interactions that
create isolated patches of meta-populations, is robust against small
disturbance. We further demonstrate that finite disturbance can accelerate the
transition from the low- to high-diversity state by helping creation of small
patches through diffusion of boundaries between species with stand-off
relation.


  The study of the evolution of cooperative behaviours --which provide benefits
to others-- and altruism --which provides benefits to others at a cost to
oneself-- has been on the core of the evolutionary game theoretical framework
since its foundation. The fast development of the theory during the last years
has improved our knowledge of the issue, but carried attached a diversification
of concepts which affected communication between scientists. Furthermore, the
main root of conflict in the struggle for life identified by Darwin, the
limited amount of resources present in any ecosystem, which is assumed to keep
a constant population size in most game theoretical studies, has only recently
been taken into account as explicitly influencing the evolutionary process.
This review concerns about both issues, the conceptual diversification during
the last years and the new results of the resource dependent models. In
extenso: After a historical introduction, a review of the most important
concepts is carried out. Then it is shown that pairwise interactions and
additive fitness determine prisoner's dilemmas (PDs) or harmony games, that two
altruists interacting together may determine a PD, and that the interaction
environment of the most cooperative and less selfish individual in any
population is always a PD. After that, it is shown that in addition to
altruists versus free-riders, the combination of free-riders and parasites
determines a fundamentally different PD. Computer simulations are then carried
out to show that random exploration of parasitism, free-riding and altruism
enables coexistence of the three strategies without the need of reciprocating,
punishing or rewarding strategies. To finish, the problem of the limitation of
resources is reviewed, showing that...


  In the study of the evolution of cooperation, resource limitations are
usually assumed just to provide a finite population size. Recently, however, it
has been pointed out that resource limitation may also generate dynamical
payoffs able to modify the original structure of the games. Here we study
analytically a phase transition from a homogeneous population of defectors when
resources are abundant to the survival of unconditional cooperators when
resources reduce below a threshold. To this end, we introduce a model of simple
agents, with no memory or ability of recognition, interacting in well-mixed
populations. The result might shed light on the role played by resource
constraints on the origin of multicellularity.


  Complex systems theory pays much attention to simple mechanisms producing
nontrivial patterns, especially power laws. However, power laws with exponent
close to one also result from complex mixtures of mechanisms that, in
isolation, would not necessarily give this type of distribution. Probably, both
paths to the power law are relevant in nature. The second gives a plausible
explanation for some instances of power laws emerging in extremely complex
systems, such as ecosystems.


  Among the different computational approaches modelling the dynamics of
isogenic cell populations, discrete stochastic models can describe with
sufficient accuracy the evolution of small size populations. However, for a
systematic and efficient study of their long-time behaviour over a wide range
of parameter values, the performance of solely direct temporal simulations
requires significantly high computational time. In addition, when the dynamics
of the cell populations exhibit non-trivial bistable behaviour, such an
analysis becomes a prohibitive task, since a large ensemble of initial states
need to be tested for the quest of possibly co-existing steady state solutions.
In this work, we study cell populations which carry the {\it lac} operon
network exhibiting solution multiplicity over a wide range of extracellular
conditions (inducer concentration). By adopting ideas from the so-called
``equation-free'' methodology, we perform systems-level analysis, which
includes numerical tasks such as the computation of {\it coarse} steady state
solutions, {\it coarse} bifurcation analysis, as well as {\it coarse} stability
analysis. Dynamically stable and unstable macroscopic (population level) steady
state solutions are computed by means of bifurcation analysis utilising short
bursts of fine-scale simulations, and the range of bistability is determined
for different sizes of cell populations. The results are compared with the
deterministic cell population balance (CPB) model, which is valid for large
populations, and we demonstrate the increased effect of stochasticity in small
size populations with asymmetric partitioning mechanisms.


  Degradation of the myelin sheath is a common pathology underlying
demyelinating neurological diseases from Multiple Sclerosis to
Leukodistrophies. Although large malformations of myelin ultrastructure in the
advanced stages of Wallerian degradation is known, its subtle structural
variations at early stages of demyelination remains poorly characterized. This
is partly due to the lack of suitable and non-invasive experimental probes
possessing sufficient resolution to detect the degradation. Here we report the
feasibility of the application of an innovative non-invasive local structure
experimental approach for imaging the changes of statistical structural
fluctuations in the first stage of myelin degeneration. Scanning micro X-ray
diffraction, using advances in synchrotron x-ray beam focusing, fast data
collection, paired with spatial statistical analysis, has been used to unveil
temporal changes in the myelin structure of dissected nerves following
extraction of the Xenopus laevis sciatic nerve. The early myelin degeneration
is a specific ordered compacted phase preceding the swollen myelin phase of
Wallerian degradation. Our demonstration of the feasibility of the statistical
analysis of SmXRD measurements using biological tissue paves the way for
further structural investigations of degradation and death of neurons and other
cells and tissues in diverse pathological states where nanoscale structural
changes may be uncovered.


  One way in which the human immunodeficiency virus (HIV-1) replicates within a
host is by infecting activated CD4+ T-cells, which then produce additional
copies of the virus. Even with the introduction of antiretroviral drug therapy,
which has been very successful over the past decade, a large obstacle to the
complete eradication of the virus is the presence of viral reservoirs in the
form of latently infected CD4+ T-cells. We consider a model of HIV infection
that describes T-cell and viral interactions, as well as, the production and
activation of latently infected T-cells. Upon determining equilibrium states of
the latent cell model, the local and global asymptotic behavior of solutions is
examined, and the basic reproduction number of the system is computed to be
strictly less than that of the corresponding three-component model, which omits
the effects of latent infection. In particular, this implies that a wider
variety of parameter values will lead to viral eradication as $t \to \infty$
due to the appearance of latent CD4+ T-cells. With this realization we discuss
possible alternative notions for eradication and persistence of infection other
than traditional tools. These results are further illustrated by a number of
numerical simulations.


  Mathematical modeling of biological systems is crucial to effectively and
efficiently developing treatments for medical conditions that plague humanity.
Often, systems of ordinary differential equations are a traditional tool used
to describe the spread of disease within the body. We consider the dynamics of
the Human Immunodeficiency Virus (HIV) in vivo during the initial stage of
infection. In particular, we examine the well-known three-component model and
prove the existence, uniqueness, and boundedness of solutions. Furthermore, we
prove that solutions remain biologically meaningful, i.e., are positivity
preserving, and perform a thorough, local stability analysis for the
equilibrium states of the system. Finally, we incorporate random coefficients
within the model and obtain numerical results to predict the probability of
infection given the transmission of the virus to a new individual.


  Since its first formulations almost a century ago, mathematical models for
disease spreading contributed to understand, evaluate and control the epidemic
processes.They promoted a dramatic change in how epidemiologists thought of the
propagation of infectious diseases.In the last decade, when the traditional
epidemiological models seemed to be exhausted, new types of models were
developed.These new models incorporated concepts from graph theory to describe
and model the underlying social structure.Many of these works merely produced a
more detailed extension of the previous results, but some others triggered a
completely new paradigm in the mathematical study of epidemic processes. In
this review, we will introduce the basic concepts of epidemiology, epidemic
modeling and networks, to finally provide a brief description of the most
relevant results in the field.


  In constant-payoff finite population games, when selection is weak and
population size is large, the one-third law serves as the condition for a
strategy to be advantageous. We generalize the result to the case where payoff
matrices are environment-dependent and provide a more general law. In this way
we model feedback from the environment and show its impact on the dynamics.


  After the recent advancement of sequencing strategies, mirror repeats have
been found to be present in the gene sequence of many organisms and species.
This presence of mirror repeats in most of the sequences indicates towards some
important functional role of these repeats. However, a simple and quick
strategy to search these repeats in a given sequence is not available. We in
this manuscript have proposed a simple and swift strategy named as FPCB
strategy to identify mirror repeats in a give sequence. The strategy includes
three simple steps of downloading sequencing in FASTA format (F), making its
parallel complement (PC) and finally performing a homology search with the
original sequence (B). At least twenty genes were analyzed using the proposed
study. A number and types of mirror repeats were observed. We have also tried
to give nomenclature to these repeats. We hope that the proposed FPCB strategy
will be quite helpful for the identification of mirror repeats in DNA or mRNA
sequence. Also the strategy may help in unraveling the functional role of
mirror repeats in various processes including evolution.


  Systems biology studies the structure and behavior of complex gene regulatory
networks. One of its aims is to develop a quantitative understanding of the
modular components that constitute such networks. The self-regulating gene is a
type of auto regulatory genetic modules which appears in over 40% of known
transcription factors in E. coli. In this work, using the technique of Poisson
Representation, we are able to provide exact steady state solutions for this
feedback model. By using the methods of synthetic biology (P.E.M. Purnick and
Weiss, R., Nature Reviews, Molecular Cell Biology, 2009, 10: 410-422) one can
build the system itself from modules like this.


  The rapid expansion of human activities threatens ocean-wide biodiversity
loss. Numerous marine animal populations have declined, yet it remains unclear
whether these trends are symptomatic of a chronic accumulation of global marine
extinction risk. We present the first systematic analysis of threat for a
globally-distributed lineage of 1,041 chondrichthyan fishes - sharks, rays, and
chimaeras. We estimate that one-quarter are threatened according to IUCN Red
List criteria due to overfishing (targeted and incidental). Large-bodied,
shallow-water species are at greatest risk and five out of the seven most
threatened families are rays. Overall chondrichthyan extinction risk is
substantially higher than for most other vertebrates, and only one-third of
species are considered safe. Population depletion has occurred throughout the
world's ice-free waters, but is particularly prevalent in the Indo-Pacific
Biodiversity Triangle and Mediterranean Sea. Improved management of fisheries
and trade is urgently needed to avoid extinctions and promote population
recovery.


  The physiology of voltage gated ion channels is complex and insights into
their gating mechanism is incomplete. Their function is best represented by
Markov models with relatively large number of distinct states that are
connected by thermodynamically feasible transitions. On the other hand, popular
models such as the one of Hodgkin and Huxley have empirical assumptions that
are generally unrealistic. Experimental protocols often dictate the number of
states in proposed Markov models, thus creating disagreements between various
observations on the same channel. Here we aim to propose a limit to the minimum
number of states required to model ion channels by employing a paradigm to
define stationary conductance in a class of ion-channels. A simple expression
is generated using concepts in elementary thermodynamics applied to protein
conformational transitions. Further, it matches well many published channel
current-voltage characteristics and parameters of the model are found to be
identifiable and easily determined from usual experimental protocols.


  Metagenomics enables the reconstruction of microbial genomes in complex
microbial communities without the need for culturing. Since assembly typically
results in fragmented genomes the grouping of genome fragments (contigs)
belonging to the same genome, a process referred to as binning, remains a major
informatics challenge. Here we present CONCOCT, a computer program that
combines three types of information - sequence composition, coverage across
multiple sample, and read-pair linkage - to automatically bin contigs into
genomes. We demonstrate high recall and precision rates of the program on
artificial as well as real human gut metagenome datasets.


  DNA hairpin molecules with periodic base sequences can be expected to exhibit
a regular coarse-grained free energy landscape (FEL) as function of the number
of open base pairs and applied mechanical force. Using a commonly employed
model, we first analyse for which types of sequences a particularly simple
landscape structure is predicted, where forward and backward energy barriers
between partly unfolded states are decreasing linearly with force. Stochastic
unfolding trajectories for such molecules with simple FEL are subsequently
generated by kinetic Monte Carlo simulations. Introducing probabilities that
can be sampled from these trajectories, it is shown how the parameters
characterising the FEL can be estimated. Already 300 trajectories, as typically
generated in experiments, provide faithful results for the FEL parameters.


  In mutualisms where there is exchange of resources for resources, or
resources for services, the resources are typically short lived compared with
the lives of the organisms that produce and make use of them. This fact allows
a separation of time scales, by which the numerical response of one species
with respect to the abundance of another can be derived mechanistically. These
responses can account for intra-specific competition, due to the partition of
the resources provided by mutualists, in this way connecting competition theory
and mutualism at a microscopic level. It is also possible to derive saturating
responses in the case of species that provide resources but expect a service in
return (e.g. pollination, seed dispersal) instead of food or nutrients. In both
situations, competition and saturation have the same underlying cause, which is
that the generation of resources occur at a finite velocity per individual of
the providing species, but their depletion happens much faster due to the
acceleration in growth rates that characterizes mutualism. The resulting models
can display all the basic features seen in many models of facultative and
obligate mutualisms, and they can be generalized from species pairs to larger
communities. The parameters of the numerical responses can be related with
quantities that can be in principle measured, and that can be related by
trade-offs, which can be useful for studying the evolution of mutualisms.
  Abstract Keywords: mutualism, resources, services, steady-state, functional
and numerical response


  Certain chemical reaction networks (CRNs) when modeled as a deterministic
dynamical system taken with mass-action kinetics have the property of reaction
network detailed balance (RNDB) which is achieved by imposing network-related
constraints on the reaction rate constants. Markov chains (whether arising as
models of CRNs or otherwise) have their own notion of detailed balance, imposed
by the network structure of the graph of the transition matrix of the Markov
chain. When considering Markov chains arising from chemical reaction networks
with mass-action kinetics, we will refer to this property as Markov chain
detailed balance (MCDB). Finally, we refer to the stochastic analog of RNDB as
Whittle stochastic detailed balance (WSDB). It is known that RNDB and WSDB are
equivalent. We prove that WSDB and MCDB are also intimately related but are not
equivalent. While RNDB implies MCDB, the converse is not true. The conditions
on rate constants that result in networks with MCDB but without RNDB are
stringent, and thus examples of this phenomenon are rare, a notable exception
is a network whose Markov chain is a birth and death process. We give a new
algorithm to find conditions on the rate constants that are required for MCDB.


  We report the first study of a network of connected enzyme-catalyzed
reactions, with added chemical and enzymatic processes that incorporate the
recently developed biochemical filtering steps into the functioning of this
biocatalytic cascade. New theoretical expressions are derived to allow simple,
few-parameter modeling of network components concatenated in such cascades,
both with and without filtering. The derived expressions are tested against
experimental data obtained for the realized network's responses, measured
optically, to variations of its input chemicals' concentrations with and
without filtering processes. We also describe how the present modeling approach
captures and explains several observations and features identified in earlier
studies of enzymatic processes when they were considered as potential network
components for multi-step information/signal processing systems.


  We present a numerical simulation study of the dynamics of filopodial growth
in the presence of active transport by myosin X motors. We employ both a
microscopic agent-based model, which captures the stochasticity of the growth
process, and a continuum mean-field theory which neglects fluctuations. We show
that in the absence of motors, filopodia growth is overestimated by the
continuum mean-field theory. Thus fluctuations slow down the growth, especially
when the protrusions are driven by a small number (10 or less) of F-actin
fibres, and when the force opposing growth (coming from membrane elasticity) is
large enough. We also show that, with typical parameter values for eukaryotic
cells, motors are unlikely to provide an actin transport mechanism which
enhances filopodial size significantly, unless the G-actin concentration within
the filopodium greatly exceeds that of the cytosol bulk. We explain these
observations in terms of order-of-magnitude estimates of diffusion-induced and
advection-induced growth of a bundle of Brownian ratchets.


  Viruses are biological nanosystems with a capsid of protein-made capsomer
units that encloses and protects the genetic material responsible for their
replication. Here we show how the geometrical constraints of the
capsomer-capsomer interaction in icosahedral capsids fix the form of the
shortest and universal truncated multipolar expansion of the two-body
interaction between capsomers. The structures of many of the icosahedral and
related virus capsids are located as single lowest energy states of this
potential energy surface. Our approach unveils relevant features of the natural
design of the capsids and can be of interest in fields of nanoscience and
nanotechnology where similar hollow convex structures are relevant.