  Prion diseases are invariably fatal and highly infectious neurodegenerative
diseases that affect a wide variety of mammalian species such as sheep and
goats, cattle, deer, elks, humans and mice etc., but rabbits have a low
susceptibility to be infected by prion diseases with respect to other species.
The stability of rabbit prion protein is due to its highly ordered
{\beta}2-{\alpha}2 loop (PLoS One 5(10) e13273 (2010); Journal of Biological
Chemistry 285(41) 31682-31693 (2010)) and a hydrophobic staple helix-capping
motif (PNAS 107(46) 19808-19813 (2010); PLoS One 8 (5) e63047 (2013)). The
{\beta}2-{\alpha}2 loop and the tail of Helix 3 it interacts with have been a
focus in prion protein structure studies. For this loop we found a salt bridge
linkage ASP177-ARG163 (O-N) (Journal of Theoretical Biology 342 (7 February
2014) 70-82 (2014)). Some scientists said on the 2FJ3.pdb NMR file of the
rabbit prion protein, the distance of ASP177-ARG163 (O-N) gives the salt bridge
of about 10 {\AA} which is nearly null in terms of energy and such a salt
bridge is not observed in their work. But, from the 3O79.pdb X-ray file of the
rabbit prion protein, we can clearly observe this salt bridge. This article
analyses the NMR and X-ray structures and gives an answer to the above
question: the salt bridge presents at pH 6.5 in the X-ray structure is simply
gone at pH 4.5 in the NMR structure is simply due to the different pH values
that impact electrostatics at the salt bridge and hence also impact the
structures. Moreover, some molecular dynamics simulation results of the X-ray
structure are reported in this article to reveal the secrets of the structural
stability of rabbit prion protein.


  Being able to design genetic regulatory networks (GRNs) to achieve a desired
cellular function is one of the main goals of synthetic biology. However,
determining minimal GRNs that produce desired time-series behaviors is
non-trivial. In this paper, we propose a 'top-down' approach to evolving small
GRNs and then use these to recursively boot-strap the identification of larger,
more complex, modular GRNs. We start with relatively dense GRNs and then use
differential evolution (DE) to evolve interaction coefficients. When the target
dynamical behavior is found embedded in a dense GRN, we narrow the focus of the
search and begin aggressively pruning out excess interactions at the end of
each generation. We first show that the method can quickly rediscover known
small GRNs for a toggle switch and an oscillatory circuit. Next we include
these GRNs as non-evolvable subnetworks in the subsequent evolution of more
complex, modular GRNs. Successful solutions found in canonical DE where we
truncated small interactions to zero, with or without an interaction penalty
term, invariably contained many excess interactions. In contrast, by
incorporating aggressive pruning and the penalty term, the DE was able to find
minimal or nearly minimal GRNs in all test problems.


  Systems biology is an emerging interdisciplinary area of research that
focuses on study of complex interactions in a biological system, such as gene
regulatory networks. The discovery of gene regulatory networks leads to a wide
range of applications, such as pathways related to a disease that can unveil in
what way the disease acts and provide novel tentative drug targets. In
addition, the development of biological models from discovered networks or
pathways can help to predict the responses to disease and can be much useful
for the novel drug development and treatments. The inference of regulatory
networks from biological data is still in its infancy stage. This paper
proposes a recurrent neural network (RNN) based gene regulatory network (GRN)
model hybridized with generalized extended Kalman filter for weight update in
backpropagation through time training algorithm. The RNN is a complex neural
network that gives a better settlement between the biological closeness and
mathematical flexibility to model GRN. The RNN is able to capture complex,
non-linear and dynamic relationship among variables. Gene expression data are
inherently noisy and Kalman filter performs well for estimation even in noisy
data. Hence, non-linear version of Kalman filter, i.e., generalized extended
Kalman filter has been applied for weight update during network training. The
developed model has been applied on DNA SOS repair network, IRMA network, and
two synthetic networks from DREAM Challenge. We compared our results with other
state-of-the-art techniques that show superiority of our model. Further, 5%
Gaussian noise has been added in the dataset and result of the proposed model
shows negligible effect of noise on the results.


  Molecular dynamics simulations have the potential to provide atomic-level
detail and insight to important questions in chemical physics that cannot be
observed in typical experiments. However, simply generating a long trajectory
is insufficient, as researchers must be able to transform the data in a
simulation trajectory into specific scientific insights. Although this analysis
step has often been taken for granted, it deserves further attention as
large-scale simulations become increasingly routine. In this perspective, we
discuss the application of Markov models to the analysis of large-scale
biomolecular simulations. We draw attention to recent improvements in the
construction of these models as well as several important open issues. In
addition, we highlight recent theoretical advances that pave the way for a new
generation of models of molecular kinetics.


  Natural selection at one site shapes patterns of genetic variation at linked
sites. Quantifying the effects of 'linked selection' on levels of genetic
diversity is key to making reliable inference about demography, building a null
model in scans for targets of adaptation, and learning about the dynamics of
natural selection. Here, we introduce the first method that jointly infers
parameters of distinct modes of linked selection, notably background selection
and selective sweeps, from genome-wide diversity data, functional annotations
and genetic maps. The central idea is to calculate the probability that a
neutral site is polymorphic given local annotations, substitution patterns, and
recombination rates. Information is then combined across sites and samples
using composite likelihood in order to estimate genome-wide parameters of
distinct modes of selection. In addition to parameter estimation, this approach
yields a map of the expected neutral diversity levels along the genome. To
illustrate the utility of our approach, we apply it to genome-wide resequencing
data from 125 lines in Drosophila melanogaster and reliably predict diversity
levels at the 1Mb scale. Our results corroborate estimates of a high fraction
of beneficial substitutions in proteins and untranslated regions (UTR). They
allow us to distinguish between the contribution of sweeps and other modes of
selection around amino acid substitutions and to uncover evidence for pervasive
sweeps in untranslated regions (UTRs). Our inference further suggests a
substantial effect of linked selection from non-classic sweeps. More generally,
we demonstrate that linked selection has had a larger effect in reducing
diversity levels and increasing their variance in D. melanogaster than
previously appreciated.


  Recombination plays a fundamental role in meiosis, ensuring the proper
segregation of chromosomes and contributing to genetic diversity by generating
novel combinations of alleles. Using data derived from directUtoUconsumer
genetic testing, we investigated patterns of recombination in over 4,200
families. Our analysis revealed a number of sex differences in the distribution
of recombination. We find the fraction of male events occurring within hotspots
to be 4.6% higher than for females. We confirm that the recombination rate
increases with maternal age, while hotspot usage decreases, with no such
effects observed in males. Finally, we show that the placement of female
recombination events becomes increasingly deregulated with maternal age, with
an increasing fraction of events appearing to escape crossover interference.


  Reconstruction of family trees, or pedigree reconstruction, for a group of
individuals is a fundamental problem in genetics. The problem is known to be
NP-hard even for datasets known to only contain siblings. Some recent methods
have been developed to accurately and efficiently reconstruct pedigrees. These
methods, however, still consider relatively simple pedigrees, for example, they
are not able to handle half-sibling situations where a pair of individuals only
share one parent. In this work, we propose an efficient method, IPED2, based on
our previous work, which specifically targets reconstruction of complicated
pedigrees that include half-siblings. We note that the presence of
half-siblings makes the reconstruction problem significantly more challenging
which is why previous methods exclude the possibility of half-siblings. We
proposed a novel model as well as an efficient graph algorithm and experiments
show that our algorithm achieves relatively accurate reconstruction. To our
knowledge, this is the first method that is able to handle pedigree
reconstruction based on genotype data only when half-sibling exists in any
generation of the pedigree.


  Identifying control strategies for biological networks is paramount for
practical applications that involve reprogramming a cell's fate, such as
disease therapeutics and stem cell reprogramming. Here we develop a novel
network control framework that integrates the structural and functional
information available for intracellular networks to predict control targets.
Formulated in a logical dynamic scheme, our approach drives any initial state
to the target state with 100% effectiveness and needs to be applied only
transiently for the network to reach and stay in the desired state. We
illustrate our method's potential to find intervention targets for cancer
treatment and cell differentiation by applying it to a leukemia signaling
network and to the network controlling the differentiation of helper T cells.
We find that the predicted control targets are effective in a broad dynamic
framework. Moreover, several of the predicted interventions are supported by
experiments.


  Biomolecules exhibit conformational fluctuations near equilibrium states,
inducing uncertainty in various biological properties in a dynamic way. We have
developed a general method to quantify the uncertainty of target properties
induced by conformational fluctuations. Using a generalized polynomial chaos
(gPC) expansion, we construct a surrogate model of the target property with
respect to varying conformational states. We also propose a method to increase
the sparsity of the gPC expansion by defining a set of conformational "active
space" random variables. With the increased sparsity, we employ the compressive
sensing method to accurately construct the surrogate model. We demonstrate the
performance of the surrogate model by evaluating fluctuation-induced
uncertainty in solvent-accessible surface area for the bovine trypsin inhibitor
protein system and show that the new approach offers more accurate statistical
information than standard Monte Carlo approaches. Further more, the constructed
surrogate model also enables us to directly evaluate the target property under
various conformational states, yielding a more accurate response surface than
standard sparse grid collocation methods. In particular, the new method
provides higher accuracy in high-dimensional systems, such as biomolecules,
where sparse grid performance is limited by the accuracy of the computed
quantity of interest. Our new framework is generalizable and can be used to
investigate the uncertainty of a wide variety of target properties in
biomolecular systems.


  Tumors initiate when a population of proliferating cells accumulates a
certain number and type of genetic and/or epigenetic alterations. The
population dynamics of such sequential acquisition of (epi)genetic alterations
has been the topic of much investigation. The phenomenon of stochastic
tunneling, where an intermediate mutant in a sequence does not reach fixation
in a population before generating a double mutant, has been studied using a
variety of computational and mathematical methods. However, the field still
lacks a comprehensive analytical description since theoretical predictions of
fixation times are only available for cases in which the second mutant is
advantageous. Here, we study stochastic tunneling in a Moran model. Analyzing
the deterministic dynamics of large populations we systematically identify the
parameter regimes captured by existing approaches. Our analysis also reveals
fitness landscapes and mutation rates for which finite populations are found in
long-lived metastable states. These are landscapes in which the final mutant is
not the most advantageous in the sequence, and resulting metastable states are
a consequence of a mutation-selection balance. The escape from these states is
driven by intrinsic noise, and their location affects the probability of
tunneling. Existing methods no longer apply. In these regimes it is the escape
from the metastable states that is the key bottleneck; fixation is no longer
limited by the emergence of a successful mutant lineage. We used the so-called
Wentzel-Kramers-Brillouin method to compute fixation times in these parameter
regimes, successfully validated by stochastic simulations. Our work fills a gap
left by previous approaches and provides a more comprehensive description of
the acquisition of multiple mutations in populations of somatic cells.


  We outline a framework for describing photo-activated biological reactions as
generalized quantum measurements of external fields, for which the biological
system takes on the role of a quantum meter. By using general arguments
regarding the Hamiltonian that describes the measurement interaction, we
identify the cases where it is essential for a complex chemical or biological
system to exhibit non-equilibrium quantum coherent dynamics in order to achieve
the requisite functionality. We illustrate the analysis by considering
measurement of the solar radiation field in photosynthesis and measurement of
the earth's magnetic field in avian magnetoreception.


  Protein conformational transitions, which are essential for function, may be
driven either by entropy or enthalpy when molecular systems comprising solute
and solvent molecules are the focus. Revealing thermodynamic origin of a given
molecular process is an important but difficult task, and general principles
governing protein conformational distributions remain elusive. Here we
demonstrate that when protein molecules are taken as thermodynamic systems and
solvents being treated as the environment, conformational entropy is an
excellent proxy for free energy and is sufficient to explain protein
conformational distributions. Specifically, by defining each unique combination
of side chain torsional state as a conformer, the population distribution (or
free energy) on an arbitrarily given order parameter is approximately a linear
function of conformational entropy. Additionally, span of various microscopic
potential energy terms is observed to be highly correlated with both
conformational entropy and free energy. Presently widely utilized free energy
proxies, including minimum potential energy, average potential energy terms by
themselves or in combination with vibrational entropy\cite, are found to
correlate with free energy rather poorly. Therefore, our findings provide a
fundamentally new theoretical base for development of significantly more
reliable and efficient next generation computational tools, where the number of
available conformers,rather than poential energy of microscopic configurations,
is the central focus. We anticipate that many related research fields,
including structure based drug design and discovery, protein design, docking
and prediction of general intermolecular interactions involving proteins, are
expected to benefit greatly.


  Motivated by tumor growth and spatial population genetics, we study the
interplay between evolutionary and spatial dynamics at the surfaces of
three-dimensional, spherical range expansions. We consider range expansion
radii that grow with an arbitrary power-law in time:
$R(t)=R_0(1+t/t^*)^{\Theta}$, where $\Theta$ is a growth exponent, $R_0$ is the
initial radius, and $t^*$ is a characteristic time for the growth, to be
affected by the inflating geometry. We vary the parameters $t^*$ and $\Theta$
to capture a variety of possible growth regimes. Guided by recent results for
two-dimensional inflating range expansions, we identify key dimensionless
parameters that describe the survival probability of a mutant cell with a small
selective advantage arising at the population frontier. Using analytical
techniques, we calculate this probability for arbitrary $\Theta$. We compare
our results to simulations of linearly inflating expansions ($\Theta=1$
spherical Fisher-Kolmogorov-Petrovsky-Piscunov waves) and treadmilling
populations ($\Theta=0$, with cells in the interior removed by apoptosis or a
similar process). We find that mutations at linearly inflating fronts have
survival probabilities enhanced by factors of 100 or more relative to mutations
at treadmilling population frontiers. We also discuss the special properties of
"marginally inflating" $(\Theta=1/2)$ expansions.


  Cancer poses danger because of its unregulated growth, development of
resistant subclones, and metastatic spread to vital organs. Although the major
transitions in cancer development are increasingly well understood, we lack
quantitative theory for how chemoprevention is predicted to affect survival. We
employ master equations and probability generating functions, the latter well
known in statistical physics, to derive the dynamics of tumor growth as a
mean-field approximation. We also study numerically the associated stochastic
birth-death process. Our findings predict exponential tumor growth when a
cancer is in its early stages of development and hyper-exponential growth
thereafter. Numerical simulations are in general agreement with our analytical
approach. We evaluate how constant, low impact treatments affect both
neoplastic growth and the frequency of chemoresistant clones. We show that
therapeutic outcomes are highly predictable for treatments starting either
sufficiently early or late in terms of initial tumor size and the initial
number of chemoresistant cells, whereas stochastic dynamics dominate therapies
starting at intermediate neoplasm sizes, with high outcome sensitivity both in
terms of tumor control and the emergence of resistant subclones. The outcome of
chemoprevention can be understood in terms of both minimal physiological
impacts resulting in long-term control and either preventing or slowing the
emergence of resistant subclones. We argue that our model and results can also
be applied to the management of early, clinically detected cancers after tumor
excision.


  Clustering of molecules on biological membranes is a widely observed
phenomenon. In some cases, such as the clustering of Ras proteins on the
membranes of mammalian cells, proper cell signaling is critically dependent on
the maintenance of these clusters. Yet, the mechanism by which clusters form
and are maintained in these systems remains unclear. Recently, it has been
discovered that activated Ras promotes further Ras activation. Here we show
using particle-based simulation that this positive feedback is sufficient to
produce persistent clusters of active Ras molecules at the nanometer scale via
a dynamic nucleation mechanism. Furthermore, we find that our cluster
statistics are consistent with experimental observations of the Ras system.
Interestingly, we show that our model does not support a Turing regime of
macroscopic reaction-diffusion patterning, and therefore that the clustering we
observe is a purely stochastic effect, arising from the coupling of positive
feedback with the discrete nature of individual molecules. These results
underscore the importance of stochastic and dynamic properties of reaction
diffusion systems for biological behavior.


  Heterochromatin protein 1 (HP1) participates in establishing and maintaining
heterochromatin via its histone modification dependent chromatin interactions.
In recent papers HP1 binding to nucleosomal arrays was measured in vitro and
interpreted in terms of nearest-neighbor cooperative binding. This mode of
chromatin interactions could lead to spreading of HP1 along the nucleosome
chain. Here, we reanalyzed previous data by representing the nucleosome chain
as a one-dimensional binding lattice, and show how the experimental HP1 binding
isotherms can be explained by a simpler model without cooperative interactions
between neighboring HP1 dimers. Based on these calculations and spatial models
of dinucleosomes and nucleosome chains, we propose that binding stoichiometry
is dependent of the nucleosome repeat length (NRL) rather than protein
interactions between HP1 dimers. According to our calculations, more open
nucleosome arrays with long DNA linkers are characterized by a larger number of
binding sites in comparison to chains with short NRL. Furthermore, we
demonstrate by Monte Carlo simulations that the NRL dependent folding of the
nucleosome chain can induce allosteric changes of HP1 binding sites. Thus, HP1
chromatin interactions can be modulated by the change of binding stoichiometry
and type of binding to condensed (methylated) and non-condensed (unmethylated)
nucleosome arrays in the absence of direct interactions between HP1 dimers.


  Life in salty habitats is a function of tolerance to those chemicals at all
critical phases of a plant's life history. The ability to withstand salt as an
established plant may require different mechanisms and plant traits than those
needed to germinate in salty soils. Seeds establishing themselves in high salt
content may respond differently depending on the light conditions and seed
germinability at the time of salty water imbibition. S. faberi seed (and S.
viridis and S. pumila) plants were discovered thriving along the seacoasts of
Southern Japan. These plants possess the ability to after-ripen, germinate,
emerge and establish themselves, grow and reproduce in the salty soils and
salt-laden atmospheres present in these windy habitats. The objectives of this
paper are to determine the effect of salt (NaCl) in water imbibed by S. faberi
seed during after-ripening and germination, as well temperature and light.
Observations made also provide insights on the possible relationship between
salt and drought tolerance. Seed germination of all phenotypes inhibited by two
percent or more of NaCl. The effects of lesser amounts of NaCl on each of the
three phenotypes was highly dependent on the specific temperature and light
conditions. The three test phenotypes provided a good range to detect responses
to salinity, allowing the observation of both stimulatory and inhibitory
responses.


  We study the conditions for persistent cooperation in an off-lattice model of
mobile agents playing the Prisoner's Dilemma game with pure, unconditional
strategies. Each agent has an exclusion radius rP, which accounts for the
population viscosity, and an interaction radius rint, which defines the
instantaneous contact network for the game dynamics. We show that, differently
from the rP=0 case, the model with finite-sized agents presents a coexistence
phase with both cooperators and defectors, besides the two absorbing phases, in
which either cooperators or defectors dominate. We provide, in addition, a
geometric interpretation of the transitions between phases. In analogy with
lattice models, the geometric percolation of the contact network (i.e.,
irrespective of the strategy) enhances cooperation. More importantly, we show
that the percolation of defectors is an essential condition for their survival.
Differently from compact clusters of cooperators, isolated groups of defectors
will eventually become extinct if not percolating, independently of their size.


  The observation that some subgraphs, called motifs, appear more often in real
networks than in their randomized counterparts has attracted much attention in
the scientific community. In the prevalent approach the detection of motifs is
based on comparing subgraph counts in a network with their counterparts in the
configuration model with the same degree distribution as the network. In this
short note we derive conservation laws that relate motif counts in the
configuration model.


  We study an abstract model for the co-evolution between mutating viruses and
the adaptive immune system. In sequence space, these two populations are
localized around transiently dominant strains. Delocalization or error
thresholds exhibit a novel interdependence because immune response is
conditional on the viral attack. An evolutionary chase is induced by stochastic
fluctuations and can occur via periodic or intermittent cycles. Using
simulations and stochastic analysis, we show how the transition between these
two dynamic regimes depends on mutation rate, immune response, and population
size.


  We study fixation in large, but finite, populations with two types, and
dynamics governed by birth-death processes. By considering a restricted class
of such processes, we derive a continuous approximation for the probability of
fixation that is valid beyond the weak-selection (WS) limit. From the
continuous approximations, we then obtain asymptotic approximations for
evolutionary dynamics with at most one equilibrium, in the selection-driven
regime, that does not preclude a weak-selection regime. As an application, we
study the fixation pattern when the infinite population limit has an interior
Evolutionary Stable Strategy (ESS): (i) we show that the fixation pattern for
the Hawk and Dove game satisfies what we term the one-half law: if the
Evolutionary Stable Strategy (ESS) is outside a small interval around
$\sfrac{1}{2}$, the fixation is of dominance type; (ii) we also show that,
outside of the weak-selection regime, the long-term dynamics of large
populations can have very little resemblance to the infinite population case;
in addition, we also present results for the case of two equilibria. Finally,
we present continuous restatements valid for large populations of two classical
concepts naturally defined in the discrete case: (i) the definition of an
$\textsc{ESS}_N$ strategy; (ii) the definition of a risk-dominant strategy. We
then present two applications of these restatements: (i) we obtain an
asymptotic definition valid in the quasi-neutral regime that recovers both the
one-third law under linear fitness and the generalised one-third law for
$d$-player games; (ii) we extend the ideas behind the (generalised) one-third
law outside the quasi-neutral regime and, as a generalisation, we introduce the
concept of critical-frequency; (iii) we recover the classification of
risk-dominant strategies for $d$-player games.


  This communication reports the increase in fluorescence resonance energy
transfer (FRET) efficiency between two laser dyes in presence of
Deoxyribonucleic acid (DNA). Two types of molecular logic gates have been
designed where DNA acts as input signal and fluorescence intensity of different
bands are taken as output signal. Use of these logic gates as DNA sensor has
been demonstrated


  The applications of Fluorescence resonance energy transfer (FRET) have
expanded tremendously in the last 25 years, and the technique has become a
staple technique in many biological and biophysical fields. FRET can be used as
spectroscopic ruler in various areas such as structural elucidation of
biological molecules and their interactions, in vitro assays, in vivo
monitoring in cellular research, nucleic acid analysis, signal transduction,
light harvesting, and metallic nanomaterials etc. Based on the mechanism of
FRET a variety of novel chemical sensors and Biosensors have been developed.
This review highlights the recent applications of sensitive and selective
ratiometric FRET based sensors.


  We introduce a statistical method that can reconstruct nonlinear genetic
models (i.e., including epistasis, or gene-gene interactions) from
phenotype-genotype (GWAS) data. The computational and data resource
requirements are similar to those necessary for reconstruction of linear
genetic models (or identification of gene-trait associations), assuming a
condition of generalized sparsity, which limits the total number of gene-gene
interactions. An example of a sparse nonlinear model is one in which a typical
locus interacts with several or even many others, but only a small subset of
all possible interactions exist. It seems plausible that most genetic
architectures fall in this category. Our method uses a generalization of
compressed sensing (L1-penalized regression) applied to nonlinear functions of
the sensing matrix. We give theoretical arguments suggesting that the method is
nearly optimal in performance, and demonstrate its effectiveness on broad
classes of nonlinear genetic models using both real and simulated human
genomes.


  Insect-borne plant diseases recur commonly in wild plants and in agricultural
crops, and are responsible for severe losses in terms of produce yield and
monetary return. Mathematical models of insect-borne plant diseases are
therefore an essential tool to help predicting the progression of an epidemic
disease and aid in decision making when control strategies are to be
implemented in the field. While retaining a generalized applicability of the
proposed model to plant epidemics vectored by insects, we specifically
investigated the epidemics of Flavescence dor\'ee phytoplasma (FD) in grapevine
plant Vitis vinifera specifically transmitted by the leafhopper Scaphoideus
titanus. The epidemiological model accounted for life-cycle stage of S.
titanus, FD pathogen cycle within S. titanus and V. vinifera, vineyard setting,
and agronomic practices. The model was comprehensively tested against
biological S. titanus life cycle and FD epidemics data collected in various
research sites in Piemonte, Italy, over multiple years. The work presented here
represents a unique suite of governing equations tested on existing independent
data and sets the basis for further modelling advances and possible
applications to investigate effectiveness of real-case epidemics control
strategies and scenarios.


  Many recent exciting discoveries have revealed the versatility of RNAs and
their importance in a variety of cellular functions which are strongly coupled
to RNA structures. To understand the functions of RNAs, some structure
prediction models have been developed in recent years. In this review, the
progress in computational models for RNA structure prediction is introduced and
the distinguishing features of many outstanding algorithms are discussed,
emphasizing three dimensional (3D) structure prediction. A promising
coarse-grained model for predicting RNA 3D structure, stability and salt effect
is also introduced briefly. Finally, we discuss the major challenges in the RNA
3D structure modeling.


  We investigate protein-protein interactions in solution by small-angle X-ray
scattering (SAXS) and theoretical modeling. The structure factor for solutions
of bovine pancreatic trypsin inhibitor (BPTI), myoglobin (Mb), and intestinal
fatty acid-binding protein (IFABP) is determined from SAXS measurements at
multiple concentrations, from Monte Carlo simulations with a coarse-grained
structure-based interaction model, and from analytic approximate solutions of
two idealized colloidal interaction models without adjustable parameters. By
combining these approaches, we find that the structure factor is essentially
determined by hard-core and screened electrostatic interactions. Other soft
short-ranged interactions (van der Waals and solvation-related) are either
individually insignificant or tend to cancel out. The structure factor is also
not significantly affected by charge fluctuations. For Mb and IFABP, with small
net charge and relatively symmetric charge distribution, the structure factor
is well described by a hard-sphere model. For BPTI, with larger net charge,
screened electrostatic repulsion is also important, but the asymmetry of the
charge distribution reduces the repulsion from that predicted by a charged
hard-sphere model with the same net charge. Such charge asymmetry may also
amplify the effect of shape asymmetry on the protein-protein potential of mean
force.


  Many population genetic models have been developed for the purpose of
inferring population size and growth rates from random samples of genetic data.
We examine two popular approaches to this problem, the coalescent and the
birth-death-sampling model, in the context of estimating population size and
birth rates in a population growing exponentially according to the birth-death
branching process. For sequences sampled at a single time, we found the
coalescent and the birth-death-sampling model gave virtually indistinguishable
results in terms of the growth rates and fraction of the population sampled,
even when sampling from a small population. For sequences sampled at multiple
time points, we find that the birth-death model estimators are subject to large
bias if the sampling process is misspecified. Since birth-death-sampling models
incorporate a model of the sampling process, we show how much of the
statistical power of birth-death-sampling models arises from the sequence of
sample times and not from the genealogical tree. This motivates the development
of a new coalescent estimator, which is augmented with a model of the known
sampling process and is potentially more precise than the coalescent that does
not use sample time information.


  We discuss a curious example for the collective mechanical behavior of
coupled non-linear monomer units entrapped in a circular filament. Within a
simple model we elucidate how multistability of monomer units and exponentially
large degeneracy of the filament's ground state emerge as a collective feature
of the closed filament. Surprisingly, increasing the monomer frustration, i.e.,
the bending prestrain within the circular filament, leads to a conformational
softening of the system. The phenomenon, that we term polymorphic crunching, is
discussed and applied to a possible scenario for membrane tube deformation by
switchable dynamin or FtsZ filaments. We find an important role of cooperative
inter-unit interaction for efficient ring induced membrane fission.


  Gray wolves were reintroduced to Yellowstone National Park (YNP) in 1995. The
population initially flourished, but since 2003 the population has experience
significant reductions due to factors that may include disease-induced
mortality, illegal hunting, park control pro- grams, vehicle induced deaths and
intra-species aggression. Despite facing similar conditions, and interference
competition with the wolves, the coyote population at YNP has persisted. In
this paper we introduce an epidemiological framework that incorporates natural,
human-caused and disease-induced mortality as well as interference competition
between two species of predators. The outcomes generated by this theoretical
framework are used to explore the impact of competition and death-induced
mechanisms on predators coexistence. It is the hope that these results on the
competitive dynamics of carnivores in Yellowstone National Park will provide
park management insights that result in policies that keep the reintroduction
of wolves successful.


  Rock is wrapped by paper, paper is cut by scissors, and scissors are crushed
by rock. This simple game is popular among children and adults to decide on
trivial disputes that have no obvious winner, but cyclic dominance is also at
the heart of predator-prey interactions, the mating strategy of side-blotched
lizards, the overgrowth of marine sessile organisms, and the competition in
microbial populations. Cyclical interactions also emerge spontaneously in
evolutionary games entailing volunteering, reward, punishment, and in fact are
common when the competing strategies are three or more regardless of the
particularities of the game. Here we review recent advances on the
rock-paper-scissors and related evolutionary games, focusing in particular on
pattern formation, the impact of mobility, and the spontaneous emergence of
cyclic dominance. We also review mean-field and zero-dimensional
rock-paper-scissors models and the application of the complex Ginzburg-Landau
equation, and we highlight the importance and usefulness of statistical physics
for the successful study of large-scale ecological systems. Directions for
future research, related for example to dynamical effects of coevolutionary
rules and invasion reversals due to multi-point interactions, are outlined as
well.


  The stability analysis of socioeconomic systems has been centered on
answering whether small perturbations when a system is in a given quantitative
state will push the system permanently to a different quantitative state.
However, typically the quantitative state of socioeconomic systems is subject
to constant change. Therefore, a key stability question that has been
under-investigated is how strong the conditions of a system itself can change
before the system moves to a qualitatively different behavior, i.e., how
structurally stable the systems is. Here, we introduce a framework to
investigate the structural stability of socioeconomic systems formed by the
network of interactions among agents competing for resources. We measure the
structural stability of the system as the range of conditions in the
distribution and availability of resources compatible with the qualitative
behavior in which all the constituent agents can be self-sustained across time.
To illustrate our framework, we study an empirical representation of the global
socioeconomic system formed by countries sharing and competing for
multinational companies used as proxy for resources. We demonstrate that the
structural stability of the system is inversely associated with the level of
competition and the level of heterogeneity in the distribution of resources.
Importantly, we show that the qualitative behavior of the observed global
socioeconomic system is highly sensitive to changes in the distribution of
resources. We believe this work provides a methodological basis to develop
sustainable strategies for socioeconomic systems subject to constantly changing
conditions.


  School environments are thought to play an important role in the community
spread of airborne infections (e.g., influenza) because of the high mixing
rates of school children. The closure of schools has therefore been proposed as
efficient mitigation strategy, with however high social and economic costs:
alternative, less disruptive interventions are highly desirable. The recent
availability of high-resolution contact networks in school environments
provides an opportunity to design micro-interventions and compare the outcomes
of alternative mitigation measures. We consider mitigation measures that
involve the targeted closure of school classes or grades based on readily
available information such as the number of symptomatic infectious children in
a class. We focus on the case of a primary school for which we have
high-resolution data on the close-range interactions of children and teachers.
We simulate the spread of an influenza-like illness in this population by using
an SEIR model with asymptomatics and compare the outcomes of different
mitigation strategies. We find that targeted class closure affords strong
mitigation effects: closing a class for a fixed period of time -equal to the
sum of the average infectious and latent durations- whenever two infectious
individuals are detected in that class decreases the attack rate by almost 70%
and strongly decreases the probability of a severe outbreak. The closure of all
classes of the same grade mitigates the spread almost as much as closing the
whole school. Targeted class closure strategies based on readily available
information on symptomatic subjects and on limited information on mixing
patterns, such as the grade structure of the school, can be almost as effective
as whole-school closure, at a much lower cost. This may inform public health
policies for the management and mitigation of influenza-like outbreaks in the
community.


  Inspired by real scenarios in Northern Patagonia, we analyze a mathematical
model of a simple trophic web with two herbivores and one predator. The studied
situations represent a common practice in the steppes of Argentine Patagonia,
where livestock are raised in a semi-wild state, either on the open range or
enclosed, coexisting with competitors and predators. In the present work, the
competing herbivores represent sheep and guanacos, while the predator is
associated with the puma. The proposed model combines the concepts of
metapopulations and patches dynamics, and includes an explicit hierarchical
competition between species, which affects their prospect to colonize an empty
patch when having to compete with other species. We perform numerical
simulations of spatially extended metapopulations assemblages of the system,
which allow us to incorporate the effects of habitat heterogeneity and
destruction. The numerical results are compared with those obtained from mean
field calculations. We find that the model provides a good theoretical
framework in several situations, including the control of the wild populations
that the ranchers exert to different extent. Furthermore, the present
formulation incorporates new terms in previously analyzed models, that help to
reveal the important effects due to the heterogeneous nature of the system.


  Biochemists know that the law of mass action is not exact and not very useful
because we cannot transfer it (with unchanged parameters) from one condition to
another. I argue that exact equations require calibrated multiscale analysis to
deal with ions. Exact theories in biochemistry must use mathematics of
interactions because biological ionic solutions, derived from seawater, are
complex (not simple) fluids. The activity of one ion depends on every other
ion. Mathematics of conservative interactions is well understood but friction
is another matter. Mathematicians now have now an energetic variational
calculus dealing with friction. Complex fluids need variational methods because
everything interacts with everything else. Mathematics designed to handle
interactions is needed to produce exact equations. If interactions are not
addressed with variational mathematics, they are bewildering. The mathematics
must include the global properties of the electric field. Flow of charge in one
place changes the flow everywhere by Kirchoff and Maxwell laws. Charge changes
physical nature as it flows through a circuit. It is ions in salt water; it is
electrons in a vacuum tube; it is quasi-particles in a semiconductor; and it is
nothing much in a vacuum capacitor (i.e., displacement current). Charge is
abstract. The physical nature of charge and current is strikingly diverse; yet,
the flow of current is exactly the same in every element in a series circuit.
The global nature of electric flow prevents the law of mass action from being
exact. The law of mass action (with rate constants that are constant) does not
know about charge. The law of mass action is about mass conservation. I believe
the law of mass action must be modified to be consistent with the Kirchoff
current law if biochemistry is to be an exact science.


  The canonical Wnt signaling pathway, mediated by $\beta$-catenin, is
crucially involved in development, adult stem cell tissue maintenance and a
host of diseases including cancer. We undertake analysis of different
mathematical models of Wnt from the literature, and compare them to a new
mechanistic model of Wnt signaling that targets spatial localization of key
molecules. Using Bayesian methods we infer parameters for each of the models to
mammalian Wnt signaling data and find that all models can fit this time course.
We are able to overcome this lack of data by appealing to algebraic methods
(concepts from chemical reaction network theory and matroid theory) to analyze
the models without recourse to specific parameter values. These approaches
provide insight into Wnt signaling: The new model (unlike any other
investigated) permits a bistable switch in the system via control of shuttling
and degradation parameters, corresponding to stem-like vs committed cell states
in the differentiation hierarchy. Our analysis also identifies groups of
variables that must be measured to fully characterize and discriminate between
competing models, and thus serves as a guide for performing minimal experiments
for model comparison.


  To bridge the gap between the sequences and 3-dimensional (3D) structures of
RNAs, some computational models have been proposed for predicting RNA 3D
structures. However, the existed models seldom consider the conditions
departing from the room/body temperature and high salt (1M NaCl), and thus
generally hardly predict the thermodynamics and salt effect. In this study, we
propose a coarse-grained model with implicit salt for RNAs to predict 3D
structures, stability and salt effect. Combined with Monte Carlo simulated
annealing algorithm and a coarse-grained force field, the model folds 46 tested
RNAs (less than or equal to 45 nt) including pseudoknots into their native-like
structures from their sequences, with an overall mean RMSD of 3.5 {\AA} and an
overall minimum RMSD of 1.9 {\AA} from the experimental structures. For 30 RNA
hairpins, the present model also gives the reliable predictions for the
stability and salt effect with the mean deviation ~ 1.0 degrees Celsius of
melting temperatures, as compared with the extensive experimental data. In
addition, the model could provide the ensemble of possible 3D structures for a
short RNA at a given temperature/salt condition.


  Cellular response to a perturbation is the result of a dynamic system of
biological variables linked in a complex network. A major challenge in drug and
disease studies is identifying the key factors of a biological network that are
essential in determining the cell's fate.
  Here our goal is the identification of perturbed pathways from
high-throughput gene expression data. We develop a three-level hierarchical
model, where (i) the first level captures the relationship between gene
expression and biological pathways using confirmatory factor analysis, (ii) the
second level models the behavior within an underlying network of pathways
induced by an unknown perturbation using a conditional autoregressive model,
and (iii) the third level is a spike-and-slab prior on the perturbations. We
then identify perturbations through posterior-based variable selection.
  We illustrate our approach using gene transcription drug perturbation
profiles from the DREAM7 drug sensitivity predication challenge data set. Our
proposed method identified regulatory pathways that are known to play a
causative role and that were not readily resolved using gene set enrichment
analysis or exploratory factor models. Simulation results are presented
assessing the performance of this model relative to a network-free variant and
its robustness to inaccuracies in biological databases.


  Empirical data on the dynamics of human face-to-face interactions across a
variety of social venues have recently revealed a number of context-independent
structural and temporal properties of human contact networks. This universality
suggests that some basic mechanisms may be responsible for the unfolding of
human interactions in the physical space. Here we discuss a simple model that
reproduces the empirical distributions for the individual, group and collective
dynamics of face-to-face contact networks. The model describes agents that move
randomly in a two-dimensional space and tend to stop when meeting "attractive"
peers, and reproduces accurately the empirical distributions.


  We study how the temporal behaviours of p53 and MDM2 are affected by stress
inducing bioactive molecules NO (Nitric Oxide) in the p53-MDM2-NO regulatory
network. We also study synchronization among a group of identical stress
systems arranged in a three dimensional array with nearest neighbour diffusive
coupling. The role of NO and effect of noise are investigated. In the single
system study, we have found three distinct types of temporal behaviour of p53,
namely, oscillation death, damped oscillation and sustain oscillation,
depending on the amount of stress induced by the NO concentration, indicating
how p53 responds to the incoming stress. The correlation among the coupled
systems increases as the value of coupling constant (\epsilon) is increased
(\gamma increases) and becomes constant after certain value of \epsilon. The
permutation entropy spectra H(\epsilon) for p53 and MDM2 as a function of
\epsilon are found to be different due to direct and indirect interaction of NO
with the respective proteins. \gamma versus \epsilon for p53 and MDM2 are found
to be similar in deterministic approach, but different in stochastic approach
and the separation between \gamma of the respective proteins as a function of
\epsilon decreases as system size increases. The role of NO is found to be
twofold: stress induced by it is prominent at small and large values of
\epsilon but synchrony inducing by it dominates in moderate range of \epsilon.
Excess stress induce apoptosis to the system.


  Research investigating the origins of life usually focuses on exploring
possible life-bearing chemistries in the pre-biotic Earth, or else on synthetic
approaches. Little work has been done exploring fundamental issues concerning
the spontaneous emergence of life using only concepts (such as information and
evolution) that are divorced from any particular chemistry. Here, I advocate
studying the probability of spontaneous molecular self-replication as a
function of the information contained in the replicator, and the environmental
conditions that might enable this emergence. I show that (under certain
simplifying assumptions) the probability to discover a self-replicator by
chance depends exponentially on the rate of formation of the monomers. If the
rate at which monomers are formed is somewhat similar to the rate at which they
would occur in a self-replicating polymer, the likelihood to discover such a
replicator by chance is increased by many orders of magnitude. I document such
an increase in searches for a self-replicator within the digital life system
avida


  In a (two-type) Wright-Fisher diffusion with directional selection and
two-way mutation, let $x$ denote today's frequency of the beneficial type, and
given $x$, let $h(x)$ be the probability that, among all individuals of today's
population, the individual whose progeny will eventually take over in the
population is of the beneficial type. Fearnhead [Fearnhead, P., 2002. The
common ancestor at a nonneutral locus. J. Appl. Probab. 39, 38-54] and Taylor
[Taylor, J. E., 2007. The common ancestor process for a Wright-Fisher
diffusion. Electron. J. Probab. 12, 808-847] obtained a series representation
for $h(x)$. We develop a construction that contains elements of both the
ancestral selection graph and the lookdown construction and includes pruning of
certain lines upon mutation. Besides being interesting in its own right, this
construction allows a transparent derivation of the series coefficients of
$h(x)$ and gives them a probabilistic meaning.


  The presence of phenomena analogous to phase transition in Statistical
Mechanics, has been suggested in the evolution of a polygenic trait under
stabilizing selection, mutation and genetic drift.
  By using numerical simulations of a model system, we analyze the evolution of
a population of $N$ diploid hermaphrodites in random mating regime. The
population evolves under the effect of drift, selective pressure in form of
viability on an additive polygenic trait, and mutation. The analysis allows to
determine a phase diagram in the plane of mutation rate and strength of
selection. The involved pattern of phase transitions is characterized by a line
of critical points for weak selective pressure (smaller than a threshold),
whereas discontinuous phase transitions, characterized by metastable
hysteresis, are observed for strong selective pressure. A finite size scaling
analysis suggests the analogy between our system and the mean field Ising model
for selective pressure approaching the threshold from weaker values. In this
framework, the mutation rate, which allows the system to explore the accessible
microscopic states, is the parameter controlling the transition from large
heterozygosity (disordered phase) to small heterozygosity (ordered one).


  Most phylogenetic analyses result in a sample of trees, but summarizing and
visualizing these samples can be challenging. Consensus trees often provide
limited information about a sample, and so methods such as consensus networks,
clustering and multidimensional scaling have been developed and applied to tree
samples. This paper describes a stochastic algorithm for constructing a
principal geodesic or line through treespace which is analogous to the first
principal component in standard Principal Components Analysis. A principal
geodesic summarizes the most variable features of a sample of trees, in terms
of both tree topology and branch lengths, and it can be visualized as an
animation of smoothly changing trees. The algorithm performs a stochastic
search through parameter space for a geodesic which minimises the sum of
squared projected distances of the data points. This procedure aims to identify
the globally optimal principal geodesic, though convergence to locally optimal
geodesics is possible. The methodology is illustrated by constructing principal
geodesics for experimental and simulated data sets, demonstrating the insight
into samples of trees that can be gained and how the method improves on a
previously published approach. A java package called GeoPhytter for
constructing and visualising principal geodesics is freely available from
www.ncl.ac.uk/~ntmwn/geophytter.


  Variation in the evolutionary process across the sites of nucleotide sequence
alignments is well established, and is an increasingly pervasive feature of
datasets composed of gene regions sampled from multiple loci and/or different
genomes. Inference of phylogeny from these data demands that we adequately
model the underlying process heterogeneity; failure to do so can lead to biased
estimates of phylogeny and other parameters. Traditionally, process
heterogeneity has been accommodated by first assigning sites to data subsets
based on relevant prior information (reflecting codon positions in
protein-coding DNA, stem and loop regions of ribosomal DNA, etc.), and then
estimating the phylogeny and other model parameters under the resulting mixed
model. Here, we consider an alternative approach for accommodating process
heterogeneity that is similar in spirit to this conventional mixed-model
approach. However, rather than treating the partitioning scheme as a fixed
assumption of the analysis, we treat the process partition as a random variable
using a Dirichlet process prior model, where the phylogeny is estimated by
integrating over all possible process partitions for the specified data
subsets. We apply this method to simulated and empirical datasets, and compare
our results to those estimated previously using conventional mixed-model
selection criteria based on Bayes factors. We find that estimation under the
Dirichlet process prior model discovers novel process partitions that may more
effectively balance error variance and estimation bias, while rendering
phylogenetic inference more robust to process heterogeneity by virtue of
integrating estimates over all possible partition schemes. (Keywords: Bayesian
phylogenetic inference; Dirichlet process prior; Markov chain Monte Carlo;
maximum likelihood; partitioned analyses; process heterogeneity.)


  Dengue is a vector borne disease transmitted to humans by {\it{Aedes
Aegypti}} mosquitoes carrying Dengue virus of different serotypes. Primarily an
urban epidemic, Dengue exhibits complex spatial and temporal dynamics,
influenced by many biological, human and environmental factors. In this work,
we study the Dengue spread for a single serotype including human mobility. We
model the Dengue spreading by using PDE reaction-diffusion and stochastic
Cellular Automata (CA) with human and vector dynamics and analyze the spatial
and temporal spreading of the disease using parameters from field studies.
Mosquito density data from Ahmedabad city serves as a proxy for climate data to
our model. We predict the dynamics of Dengue incidence and compare it to the
reported data on the prevalence of the disease from 2006-2012. We find that for
certain infection rates, CA model closely reproduces observed peaks and
intensity. We have used statistical model of human mobility with exponential
step length distribution to study mobility effects on Dengue spreading within
the city. We find an interesting result that although human mobility makes the
infection spread faster, there is an apparent early suppression of the epidemic
compared to immobile humans. The primary reason for decline is that mobility
causes secondary and tertiary waves of infected individuals who recover in
short time span and act as a barricade for the primary wave.


  The process of recombination in population genetics, in its deterministic
limit, leads to a nonlinear ODE in the Banach space of finite measures on a
locally compact product space. It has an embedding into a larger family of
nonlinear ODEs that permits a systematic analysis with lattice-theoretic
methods for general partitions of finite sets. We discuss this type of system,
reduce it to an equivalent finite-dimensional nonlinear problem, and establish
a connection with an ancestral partitioning process, backward in time. We solve
the finite-dimensional problem recursively for generic sets of parameters and
briefly discuss the singular cases, and how to extend the solution to this
situation.


  The resurgence of vector-borne diseases is an increasing public health
concern, and there is a need for a better understanding of their dynamics. For
a number of diseases, e.g. dengue and chikungunya, this resurgence occurs
mostly in urban environments, which are naturally very heterogeneous,
particularly due to population circulation. In this scenario, there is an
increasing interest in both multi-patch and multi-group models for such
diseases. In this work, we study the dynamics of a vector borne disease within
a class of multi-group models that extends the classical Bailey-Dietz model.
This class includes many of the proposed models in the literature, and it can
accommodate various functional forms of the infection force. For such models,
the vector-host/host-vector contact network topology gives rise to a bipartite
graph which has different properties from the ones usually found in directly
transmitted diseases. Under the assumption that the contact network is strongly
connected, we can define the basic reproductive number $\mathcal{R}_0$ and show
that this system has only two equilibria: the so called disease free
equilibrium (DFE); and a unique interior equilibrium---usually termed the
endemic equilibrium (EE)---that exists if, and only if, $\mathcal{R}_0>1$. We
also show that, if $\mathcal{R}_0\leq1$, then the DFE equilibrium is globally
asymptotically stable, while when $\mathcal{R}_0>1$, we have that the EE is
globally asymptotically stable.


  Dynamical properties of a long polymer ring in a melt of unknotted and
unconcatenated rings are calculated. We re-examine and generalize the well
known model of a ring confined to a lattice of topological obstacles in the
light of the recently developed Flory theory of unentangled rings which maps
every ring on an annealed branched polymer and establishes that the backbone
associated with each ring follows self-avoiding rather than Gaussian random
walk statistics. We find the scaling of ring relaxation time and diffusion
coefficient with ring length, as well as time dependence of stress relaxation
modulus, zero shear viscosity and mean square averaged displacements of both
individual monomers and ring's mass center. Our results agree within error bars
with all available experimental and simulations data of the ring melt, although
the quality of the data so far is insufficient to make a definitive judgment
for or against the annealed tree theory. In the end we review briefly the
relation between our findings and experimental data on chromatin dynamics.


  In the framework of evolutionary games with institutional reciprocity,
limited incentives are at disposal for rewarding cooperators and punishing
defectors. In the simplest case, it can be assumed that, depending on their
strategies, all players receive equal incentives from the common pool. The
question arises, however, what is the optimal distribution of institutional
incentives? How should we best reward and punish individuals for cooperation to
thrive? We study this problem for the public goods game on a scale-free
network. We show that if the synergetic effects of group interactions are weak,
the level of cooperation in the population can be maximized simply by adopting
the simplest "equal distribution" scheme. If synergetic effects are strong,
however, it is best to reward high-degree nodes more than low-degree nodes.
These distribution schemes for institutional rewards are independent of payoff
normalization. For institutional punishment, however, the same optimization
problem is more complex, and its solution depends on whether absolute or
degree-normalized payoffs are used. We find that degree-normalized payoffs
require high-degree nodes be punished more lenient than low-degree nodes.
Conversely, if absolute payoffs count, then high-degree nodes should be
punished stronger than low-degree nodes.


  A recent paper by James et al. finds that mutualistic interactions decrease
the biodiversity of model ecosystems. However, this result can be reverted if
we consider ecological trade-offs and choose parameters suitable for sparse
mutualistic networks instead of fully connected networks.


  Empirical temporal networks display strong heterogeneities in their dynamics,
which profoundly affect processes taking place on these networks, such as rumor
and epidemic spreading. Despite the recent wealth of data on temporal networks,
little work has been devoted to the understanding of how such heterogeneities
can emerge from microscopic mechanisms at the level of nodes and links. Here we
show that long-term memory effects are present in the creation and
disappearance of links in empirical networks. We thus consider a simple
generative modeling framework for temporal networks able to incorporate these
memory mechanisms. This allows us to study separately the role of each of these
mechanisms in the emergence of heterogeneous network dynamics. In particular,
we show analytically and numerically how heterogeneous distributions of contact
durations, of inter-contact durations and of numbers of contacts per link
emerge. We also study the individual effect of heterogeneities on dynamical
processes, such as the paradigmatic Susceptible-Infected epidemic spreading
model. Our results confirm in particular the crucial role of the distributions
of inter-contact durations and of the numbers of contacts per link.


  In this paper, we study the structure and dynamical properties of protein
contact networks with respect to other biological networks, together with
simulated archetypal models acting as probes. We consider both classical
topological descriptors, such as the modularity and statistics of the shortest
paths, and different interpretations in terms of diffusion provided by the
discrete heat kernel, which is elaborated from the normalized graph Laplacians.
A principal component analysis shows high discrimination among the network
types, either by considering the topological and heat kernel based vector
characterizations. Furthermore, a canonical correlation analysis demonstrates
the strong agreement among those two characterizations, providing thus an
important justification in terms of interpretability for the heat kernel.
Finally, and most importantly, the focused analysis of the heat kernel provides
a way to yield insights on the fact that proteins have to satisfy specific
structural design constraints that the other considered networks do not need to
obey. Notably, the heat trace decay of an ensemble of varying-size proteins
denotes subdiffusion, a peculiar property of proteins.


  Discrete-state, continuous-time Markov models are widely used in the modeling
of biochemical reaction networks. Their complexity often precludes analytic
solution, and we rely on stochastic simulation algorithms to estimate system
statistics. The Gillespie algorithm is exact, but computationally costly as it
simulates every single reaction. As such, approximate stochastic simulation
algorithms such as the tau-leap algorithm are often used. Potentially
computationally more efficient, the system statistics generated suffer from
significant bias unless tau is relatively small, in which case the
computational time can be comparable to that of the Gillespie algorithm. The
multi-level method (Anderson and Higham, Multiscale Model. Simul. 2012) tackles
this problem. A base estimator is computed using many (cheap) sample paths at
low accuracy. The bias inherent in this estimator is then reduced using a
number of corrections. Each correction term is estimated using a collection of
paired sample paths where one path of each pair is generated at a higher
accuracy compared to the other (and so more expensive). By sharing random
variables between these paired paths the variance of each correction estimator
can be reduced. This renders the multi-level method very efficient as only a
relatively small number of paired paths are required to calculate each
correction term. In the original multi-level method, each sample path is
simulated using the tau-leap algorithm with a fixed value of $\tau$. This
approach can result in poor performance when the reaction activity of a system
changes substantially over the timescale of interest. By introducing a novel,
adaptive time-stepping approach where $\tau$ is chosen according to the
stochastic behaviour of each sample path we extend the applicability of the
multi-level method to such cases. We demonstrate the efficiency of our method
using a number of examples.


  We study genetic networks that produce many species of non-coding RNA
molecules that are present at a moderate density, as typically exists in the
cell. The associations of the many species of these RNA are modeled physically,
taking into account the equilibrium constants between bound and unbound states.
By including the pair-wise binding of the many RNA species, the network becomes
highly interconnected and shows different properties than the usual type of
genetic network. It shows much more robustness to mutation, and also rapid
evolutionary adaptation in an environment that oscillates in time. This
provides a possible explanation for the weak evolutionary constraints seen in
much of the non-coding RNA that has been studied.


  The regulating mechanism of miRNA on p53 dynamics in p53-MDM2-miRNA model
network incorporating reactive oxygen species (ROS) is studied. The study shows
that miRNA drives p53 dynamics at various states, namely, stabilized states and
oscillating states (damped and sustain oscillation). We found the co-existence
of these states within certain range of the concentartion level of miRNA in the
system. This co-existence in p53 dynamics is the signature of the system's
survival at various states, normal, activated and apoptosis driven by a
constant concentration of miRNA.


  We study the temporal and the synchronous behaviours in p53-Mdm2 regulatory
network due to the interaction of its complex network components with the
nitric oxide molecule. In single cell process, increase in nitric oxide
concentration gives rise the transition to various p53 temporal behaviours,
namely fixed point oscillation, damped oscillation and sustain oscillation
indicating stability, weakly activated and strongly activated states. The noise
in stochastic system is found to help to reach these states much faster as
compared to deterministic case which is evident from permutation entropy
dynamics. In coupled system with nitric oxide as diffusively coupling molecule,
we found nitric oxide as strong coupling molecule within a certain range of
coupling strength \epsilon beyond which it become weak synchronizing agent. We
study these effects by using correlation like synchronization indicator \gamma
obtained from permutation entropies of the coupled system, and found five
important regimes in (\epsilon-\gamma) phase diagram, indicating
desynchronized, transition, strongly synchronized, moderately synchronized and
weakly synchronized regimes respectively. We claim that there is the
competition between the toxicity and the synchronizing role of nitric oxide
that lead the cell in different stressed conditions.


  Recent studies have consistently inferred high rates of adaptive molecular
evolution between Drosophila species. At the same time, the Drosophila genome
evolves under different rates of recombination, which results in partial
genetic linkage between alleles at neighboring genomic loci. Here we analyze
how linkage correlations affect adaptive evolution. We develop a new inference
method for adaptation that takes into account the effect on an allele at a
focal site caused by neighboring deleterious alleles (background selection) and
by neighboring adaptive substitutions (hitchhiking). Using complete genome
sequence data and fine-scale recombination maps, we infer a highly
heterogeneous scenario of adaptation in Drosophila. In high-recombining
regions, about 50% of all amino acid substitutions are adaptive, together with
about 20% of all substitutions in proximal intergenic regions. In
low-recombining regions, only a small fraction of the amino acid substitutions
are adaptive, while hitchhiking accounts for the majority of these changes.
Hitchhiking of deleterious alleles generates a substantial collateral cost of
adaptation, leading to a fitness decline of about 30/2N per gene and per
million years in the lowest-recombining regions. Our results show how
recombination shapes rate and efficacy of the adaptive dynamics in eukaryotic
genomes.


  Recent tumor genome sequencing confirmed that one tumor often consists of
multiple cell subpopulations (clones) which bear different, but related,
genetic profiles such as mutation and copy number variation profiles. Thus far,
one tumor has been viewed as a whole entity in cancer functional studies. With
the advances of genome sequencing and computational analysis, we are able to
quantify and computationally dissect clones from tumors, and then conduct
clone-based analysis. Emerging technologies such as single-cell genome
sequencing and RNA-Seq could profile tumor clones. Thus, we should reconsider
how to conduct cancer systems biology studies in the genome sequencing era. We
will outline new directions for conducting cancer systems biology by
considering that genome sequencing technology can be used for dissecting,
quantifying and genetically characterizing clones from tumors. Topics discussed
in Part 1 of this review include computationally quantifying of tumor
subpopulations; clone-based network modeling, cancer hallmark-based networks
and their high-order rewiring principles and the principles of cell survival
networks of fast-growing clones.


  A tumor often consists of multiple cell subpopulations (clones). Current
chemo-treatments often target one clone of a tumor. Although the drug kills
that clone, other clones overtake it and the tumor reoccurs. Genome sequencing
and computational analysis allows to computational dissection of clones from
tumors, while singe-cell genome sequencing including RNA-Seq allows to
profiling of these clones. This opens a new window for treating a tumor as a
system in which clones are evolving. Future cancer systems biology studies
should consider a tumor as an evolving system with multiple clones. Therefore,
topics discussed in Part 2 of this review include evolutionary dynamics of
clonal networks, early-warning signals for formation of fast-growing clones,
dissecting tumor heterogeneity, and modeling of clone-clone-stroma interactions
for drug resistance. The ultimate goal of the future systems biology analysis
is to obtain a whole-system understanding of a tumor and therefore provides a
more efficient and personalized management strategies for cancer patients.


  GeneNetMiner is standalone software which parses the sentences of iHOP and
captures regulatory relations. The regulatory relations are either gene gene
regulations or gene biological processes relations. Capturing of gene
biological process relations is a unique feature for the tools of this kind.
These relations can be used to build up gene regulatory networks for specific
biological processes, diseases, or phenotypes. Users are able to search genes
and biological processes to find the regulatory relationships between them.
Each regulatory relationship has been assigned a confidence score, which
indicates the probability of the true relation. Furthermore, it reports the
sentence containing the queried terms, which allows users to manually checking
whether the relation is true if they wish. GeneNetMiner is able to accurately
capture the regulatory relationships between genes from literature.


  The aftermath of influenza infection is determined by a complex set of
host-pathogen interactions, where genomic variability on both viral and host
sides influences the final outcome. Although there exists large body of
literature describing influenza virus variability, only a very small fraction
covers the issue of host variance. The goal of this review is to explore the
variability of host genes responsible for host-pathogen interactions, paying
particular attention to genes responsible for the presence of sialylated
glycans in the host endothelial membrane, mucus, genes used by viral immune
escape mechanisms, and genes particularly expressed after vaccination, since
they are more likely to have a direct influence on the infection outcome.


  The reconstruction of a species tree from genomic data faces a double hurdle.
First, the (gene) tree describing the evolution of each gene may differ from
the species tree, for instance, due to incomplete lineage sorting. Second, the
aligned genetic sequences at the leaves of each gene tree provide merely an
imperfect estimate of the topology of the gene tree. In this note, we
demonstrate formally that a basic statistical problem arises if one tries to
avoid accounting for these two processes and analyses the genetic data directly
via a concatenation approach. More precisely, we show that, under the
multi-species coalescent with a standard site substitution model, maximum
likelihood estimation on sequence data that has been concatenated across genes
and performed under the incorrect assumption that all sites have evolved
independently and identically on a fixed tree is a statistically inconsistent
estimator of the species tree. Our results provide a formal justification of
simulation results described of Kubatko and Degnan (2007) and others, and
complements recent theoretical results by DeGorgio and Degnan (2010) and
Chifman and Kubtako (2014).


  The primary goal in cluster analysis is to discover natural groupings of
objects. The field of cluster analysis is crowded with diverse methods that
make special assumptions about data and address different scientific aims.
Despite its shortcomings in accuracy, hierarchical clustering is the dominant
clustering method in bioinformatics. Biologists find the trees constructed by
hierarchical clustering visually appealing and in tune with their evolutionary
perspective. Hierarchical clustering operates on multiple scales
simultaneously. This is essential, for instance, in transcriptome data where
one may be interested in making qualitative inferences about how lower-order
relationships like gene modules lead to higher-order relationships like
pathways or biological processes. The recently developed method of convex
clustering preserves the visual appeal of hierarchical clustering while
ameliorating its propensity to make false inferences in the presence of
outliers and noise. The current paper exploits the proximal distance principle
to construct a novel algorithm for solving the convex clustering problem. The
solution paths generated by convex clustering reveal relationships between
clusters that are hidden by static methods such as k-means clustering. Our
convex clustering software separates parameters, accommodates missing data, and
supports prior information on relationships. The software is implemented on ATI
and nVidia graphics processing units (GPUs) for maximal speed. Several
biological examples illustrate the strengths of convex clustering and the
ability of the proximal distance algorithm to handle high-dimensional problems.


  How to induce differentiated cells into pluripotent cells has elicited
researchers' interests for a long time since pluripotent stem cells are able to
offer remarkable potential in numerous subfields of biological research.
However, the nature of cell reprogramming, especially the mechanisms still
remain elusive for the sake of most protocols of inducing pluripotent stem
cells were discovered by screening but not from the knowledge of gene
regulation networks. Generally there are two hypotheses to elucidate the
mechanism termed as elite model and stochastic model which regard reprogramming
process a deterministic process or a stochastic process, respectively. However,
the difference between these two models cannot yet be discriminated
experimentally. Here we used a general mathematical model to elucidate the
nature of cell reprogramming which can fit both hypotheses. We investigated
this process from a novel perspective, the timing. We calculated the time of
reprogramming in a general way and find that noise would play a significant
role if the stochastic hypothesis holds. Thus the two hypotheses may be
discriminated experimentally by counting the time of reprogramming in different
magnitudes of noise. Because our approach is general, our results should
facilitate broad studies of rational design of cell reprogramming protocols.


  The evolution of the allelic proportion $x$ of a biallelic locus subject to
the forces of mutation and drift is investigated in a diffusion model, assuming
small scaled mutation rates. The overall scaled mutation rate is parametrized
with $\theta=(\mu_1+\mu_0)N$ and the ratio of mutation rates with
$\alpha=\mu_1/(\mu_1+\mu_0)=1-\beta$. The equilibrium density of this process
is beta with parameters $\alpha\theta$ and $\beta\theta$. Away from
equilibrium, the transition density can be expanded into a series of modified
Jacobi polynomials. If the scaled mutation rates are small, i.e., $\theta \ll
1$, it may be assumed that polymorphism derives from mutations at the
boundaries. A model, where the interior dynamics conform to the pure drift
diffusion model and the mutations are entering from the boundaries is derived.
In equilibrium, the density of the proportion of polymorphic alleles, \ie\ $x$
within the polymorphic region $[1/N,1-1/N]$, is
$\alpha\beta\theta(\tfrac1x+\tfrac1{1-x})=\tfrac{\alpha\beta\theta}{x(1-x)}$,
while the mutation bias $\alpha$ influences the proportion of monomorphic
alleles at 0 and 1. Analogous to the expansion with modified Jacobi
polynomials, a series expansion of the transition density is derived, which is
connected to Kimura's well known solution of the pure drift model using
Gegenbauer polynomials. Two temporal and two spatial regions are separated. The
eigenvectors representing the spatial component within the polymorphic region
depend neither on the on the scaled mutation rate $\theta$ nor on the mutation
bias $\alpha$. Therefore parameter changes, e.g., growing or shrinking
populations or changes in the mutation bias, can be modeled relatively easily,
without the change of the eigenfunctions necessary for the series expansion
with Jacobi polynomials.


  We study two important roles of miRNA as stress inducer and synchronizing
agent in regulating diffusively coupled genetic oscillators within stochastic
formalism and competition between them. We identify optimal value of coupling
constant at which the rate of synchronization is strongest, below this value
sychronizing activity dominates stress activity, and above it stress activity
destroys synchronizing pattern. The concentration of miRNA maintained in
individual oscillator modulates the synchronization activity achieved by
diffusing miRNA.


  Linear mixed models (LMMs) have emerged as the method of choice for
confounded genome-wide association studies. However, the performance of LMMs in
non-randomly ascertained case-control studies deteriorates with increasing
sample size. We propose a framework called LEAP (Liability Estimator As a
Phenotype, https://github.com/omerwe/LEAP) that tests for association with
estimated latent values corresponding to severity of phenotype, and demonstrate
that this can lead to a substantial power increase.


  Protein binding and function often involves conformational changes. Advanced
NMR experiments indicate that these conformational changes can occur in the
absence of ligand molecules (or with bound ligands), and that the ligands may
'select' protein conformations for binding (or unbinding). In this review, we
argue that this conformational selection requires transition times for ligand
binding and unbinding that are small compared to the dwell times of proteins in
different conformations, which is plausible for small ligand molecules. Such a
separation of timescales leads to a decoupling and temporal ordering of
binding/unbinding events and conformational changes. We propose that
conformational-selection and induced-change processes (such as induced fit) are
two sides of the same coin, because the temporal ordering is reversed in
binding and unbinding direction. Conformational-selection processes can be
characterized by a conformational excitation that occurs prior to a binding or
unbinding event, while induced-change processes exhibit a characteristic
conformational relaxation that occurs after a binding or unbinding event. We
discuss how the ordering of events can be determined from relaxation rates and
effective on- and off-rates determined in mixing experiments, and from the
conformational exchange rates measured in advanced NMR or single-molecule FRET
experiments. For larger ligand molecules such as peptides, conformational
changes and binding events can be intricately coupled and exhibit aspects of
conformational-selection and induced-change processes in both binding and
unbinding direction.


  This paper reviews strategies for solving problems encountered when analyzing
large genomic data sets and describes the implementation of those strategies in
R by packages from the Bioconductor project. We treat the scalable processing,
summarization and visualization of big genomic data. The general ideas are well
established and include restrictive queries, compression, iteration and
parallel computing. We demonstrate the strategies by applying Bioconductor
packages to the detection and analysis of genetic variants from a whole genome
sequencing experiment.


  We introduce two models for multi-type random trees motivated by studies of
trait dependence in the evolution of species. Our discrete time model, the
multi-type ERM tree, is a generalization of Markov propagation models on a
random tree generated by a binary search or `equal rates Markov' mechanism. Our
continuous time model, the multi-type Yule tree with mutations, is a multi-type
generalization of the tree generated by a pure birth or Yule process. We study
type dependent topological properties of these two random tree models. We
derive asymptotic results that allow one to infer model parameters from data on
types at the leaves and at branch-points that are one step away from the
leaves.


  We introduce an evolutionary metacommunity of multitrophic food webs on
several habitats coupled by migration. In contrast to previous studies that
focus either on evolutionary or on spatial aspects, we include both and
investigate the interplay between them. Locally, the species emerge, interact
and go extinct according to the rules of the well-known evolutionary food web
model proposed by Loeuille and Loreau in 2005. Additionally, species are able
to migrate between the habitats. With random migration, we are able to
reproduce common trends in diversity-dispersal relationships: Regional
diversity decreases with increasing migration rates, whereas local diversity
can increase in case of a low level of dispersal. Moreover, we find that the
total biomasses in the different patches become similar even when species
composition remains different. With adaptive migration, we observe species
compositions that differ considerably between patches and contain species that
are descendant from ancestors on both patches. This result indicates that the
combination of spatial aspects and evolutionary processes affects the structure
of food webs in different ways than each of them alone.


  A semiclassical approximation approach based on the Maslov complex germ
method is considered in detail for the 1D nonlocal
Fisher-Kolmogorov-Petrovskii-Piskunov equation under the supposition of weak
diffusion. In terms of the semiclassical formalism developed, the original
nonlinear equation is reduced to an associated linear partial differential
equation and some algebraic equations for the coefficients of the linear
equation with a given accuracy of the asymptotic parameter. The solutions of
the nonlinear equation are constructed from the solutions of both the linear
equation and the algebraic equations. The solutions of the linear problem are
found with the use of symmetry operators. A countable family of the leading
terms of the semiclassical asymptotics is constructed in explicit form.
  The semiclassical asymptotics are valid by construction in a finite time
interval. We construct asymptotics which are different from the semiclassical
ones and can describe evolution of the solutions of the
Fisher-Kolmogorov-Petrovskii-Piskunov equation at large times. In the example
considered, an initial unimodal distribution becomes multimodal, which can be
treated as an example of a space structure.


  Individual cancer cells carry a bewildering number of distinct genomic
alterations i.e., copy number variations and mutations, making it a challenge
to uncover genomic-driven mechanisms governing tumorigenesis. Here we performed
exome-sequencing on several breast cancer cell lines which represent two
subtypes, luminal and basal. We integrated this sequencing data, and functional
RNAi screening data (i.e., for identifying genes which are essential for cell
proliferation and survival), onto a human signaling network. Two
subtype-specific networks were identified, which potentially represent
core-signaling mechanisms underlying tumorigenesis. Within both networks, we
found that genes were differentially affected in different cell lines; i.e., in
some cell lines a gene was identified through RNAi screening whereas in others
it was genomically altered. Interestingly, we found that highly connected
network genes could be used to correctly classify breast tumors into subtypes
based on genomic alterations. Further, the networks effectively predicted
subtype-specific drug targets, which were experimentally validated.


  In recent years, cancer genome sequencing and other high-throughput studies
of cancer genomes have generated many notable discoveries. In this review,
Novel genomic alteration mechanisms, such as chromothripsis (chromosomal
crisis) and kataegis (mutation storms), and their implications for cancer are
discussed. Genomic alterations spur cancer genome evolution. Thus, the
relationship between cancer clonal evolution and cancer stems cells is
commented. The key question in cancer biology concerns how these genomic
alterations support cancer development and metastasis in the context of
biological functioning. Thus far, efforts such as pathway analysis have
improved the understanding of the functional contributions of genetic mutations
and DNA copy number variations to cancer development, progression and
metastasis. However, the known pathways correspond to a small fraction,
plausibly 5-10%, of somatic mutations and genes with an altered copy number. To
develop a comprehensive understanding of the function of these genomic
alterations in cancer, an integrative network framework is proposed and
discussed. Finally, the challenges and the directions of studying cancer omic
data using an integrative network approach are commented.


  Variability is ubiquitous in nature and a fundamental feature of complex
systems. Few studies, however, have investigated variance itself as a trait
under genetic control. By focusing primarily on trait means and ignoring the
effect of alternative alleles on trait variability, we may be missing an
important axis of genetic variation contributing to phenotypic differences
among individuals. To study genetic effects on individual-to-individual
phenotypic variability (or intragenotypic variability), we used a panel of
Drosophila inbred lines and focused on locomotor handedness, in an assay
optimized to measure variability. We discovered that some lines had
consistently high levels of intragenotypic variability among individuals while
others had low levels. We demonstrate that the degree of variability is itself
heritable. Using a genome-wide association study (GWAS) for the degree of
intragenotypic variability as the phenotype across lines, we identified several
genes expressed in the brain that affect variability in handedness without
affecting the mean. One of these genes, Ten-a, implicated a neuropil in the
central complex of the fly brain as influencing the magnitude of behavioral
variability, a brain region involved in sensory integration and locomotor
coordination. We have validated these results using genetic deficiencies, null
alleles, and inducible RNAi transgenes. This study reveals the constellation of
phenotypes that can arise from a single genotype and it shows that different
genetic backgrounds differ dramatically in their propensity for phenotypic
variability. Because traditional mean-focused GWASs ignore the contribution of
variability to overall phenotypic variation, current methods may miss important
links between genotype and phenotype.


  Phylogenetic methods typically rely on an appropriate model of how data
evolved in order to infer an accurate phylogenetic tree. For molecular data,
standard statistical methods have provided an effective strategy for extracting
phylogenetic information from aligned sequence data when each site (character)
is subject to a common process. However, for other types of data (e.g.
morphological data), characters can be too ambiguous, homoplastic or saturated
to develop models that are effective at capturing the underlying process of
change. To address this, we examine the properties of a classic but neglected
method for inferring splits in an underlying tree, namely, maximum
compatibility. By adopting a simple and extreme model in which each character
either fits perfectly on some tree, or is entirely random (but it is not known
which class any character belongs to) we are able to derive exact and explicit
formulae regarding the performance of maximum compatibility. We show that this
method is able to identify a set of non-trivial homoplasy-free characters, when
the number $n$ of taxa is large, even when the number of random characters is
large. By contrast, we show that a method that makes more uniform use of all
the data --- maximum parsimony --- can provably estimate trees in which {\em
none} of the original homoplasy-free characters support splits.


  Brazil holds approximately 1/3 of population living infected with AIDS
(acquired immunodeficiency syndrome) in Central and South Americas, and it was
also the first developing country to implement a large-scale control and
intervention program against AIDS epidemic. In this scenario, we investigate
the temporal evolution and current status of the AIDS epidemic in Brazil.
Specifically, we analyze records of annual absolute frequency of cases for more
than 5000 cities for the first 33 years of the infection in Brazil. We found
that (i) the annual absolute frequencies exhibit a logistic-type growth with an
exponential regime in the first few years of the AIDS spreading; (ii) the
actual reproduction number decaying as a power law; (iii) the distribution of
the annual absolute frequencies among cities decays with a power law behavior;
(iv) the annual absolute frequencies and the number of inhabitants have an
allometric relationship; (v) the temporal evolution of the annual absolute
frequencies have different profile depending on the average annual absolute
frequencies in the cities. These findings yield a general quantitative
description of the AIDS infection dynamics in Brazil since the beginning. They
also provide clues about the effectiveness of treatment and control programs
against the infection, that has had a different impact depending on the number
of inhabitants of cities. In this framework, our results give insights into the
overall dynamics of AIDS epidemic, which may contribute to select empirically
accurate models.


  The networks of predator-prey interactions in ecological systems are
remarkably complex, but nevertheless surprisingly stable in terms of long term
persistence of the system as a whole. In order to understand the mechanism
driving the complexity and stability of such food webs, we developed an
eco-evolutionary model in which new species emerge as modifications of existing
ones and dynamic ecological interactions determine which species are viable.
The food-web structure thereby emerges from the dynamical interplay between
speciation and trophic interactions. The proposed model is less abstract than
earlier evolutionary food web models in the sense that all three evolving
traits have a clear biological meaning, namely the average body mass of the
individuals, the preferred prey body mass, and the width of their potential
prey body mass spectrum. We observed networks with a wide range of sizes and
structures and high similarity to natural food webs. The model networks exhibit
a continuous species turnover, but massive extinction waves that affect more
than $50 \%$ of the network are not observed.


  In this paper we show the existence of three dimensional rigid, and thus
unfoldable, lattice conformations. The structure we found has 450+ bonds, and
we provide a computer assisted proof of the existence of such structures. The
existence of such rigid structures illustrates why protein folding problems are
hard also in dimension three. The existence of two rigid two dimensional
structures was shown earlier. This work answers question 8 in \cite{2D} in the
affirmative: rigid self avoiding walks exist in three dimensional lattice
configurations.


  Intracellular transport of vesicular cargos, organelles, and other
macromolecules is an essential process to move large items through a crowded,
and inhomogeneous cellular environment. In an effort to dissect the fundamental
effects of crowding and an increasingly complex cellular environment on the
transport of individual motor proteins, we have performed in vitro
reconstitution experiments with single kinesin-1 motors walking on microtubules
in the presence of crowding agents and transient microtubule-associated
proteins that more closely emulate the cellular environment. Macromolecular
crowding due to inert polymers caused enhanced run lengths of motors, but
displayed an increased tendency for non-specific motor association and
diffusion, most likely due to depletion interactions. We found that transiently
bound associated proteins slowed forward motion, but did not drastically affect
the association times, in opposition to previously reported obstacle properties
of stably associated microtubule-associated proteins, such as the neuronal
protein tau. Such studies of the transport properties of molecular motors in
increasingly complex reconstituted environments are important to illuminate the
fundamental biophysical principles underlying the essential process of
intracellular cargo transport.


  Recent experiments have shown that stochastic effects exerted at the level of
translation contribute a substantial portion of the variation in abundance of
proteins expressed at moderate to high levels. This study analyzes
translational noise arising from fluctuations in residue-specific elongation
rates. The resulting variation has multiplicative components that lead
individual protein abundances in a population to exhibit approximately
log-normal behavior. The high variability inherent in the process leads to
parameter variation that has the features of a type of noise in biological
systems that has been characterized as extrinsic. Elongation rate variation
offers an accounting for a major component of extrinsic noise, and the analysis
provided here highlights a probability distribution that is a natural extension
of the Poisson and has broad applicability to many types of multiplicative
noise processes.


  We propose and analyse an optimal control problem where the control system is
a mathematical model for tuberculosis that considers reinfection. The control
functions represent the fraction of early latent and persistent latent
individuals that are treated. Our aim is to study how these control measures
should be implemented, for a certain time period, in order to reduce the number
of active infected individuals, while minimizing the interventions
implementation costs. The optimal intervention is compared along different
epidemiological scenarios, by varying the transmission coefficient. The impact
of variation of the risk of reinfection, as a result of acquired immunity to a
previous infection for treated individuals on the optimal controls and
associated solutions, is analysed. A cost-effectiveness analysis is done, to
compare the application of each one of the control measures, separately or in
combination.