  On the base of advantages in gene geography and anthropophenetics the
phenogeographical method for anthropological research is initiated and
experienced using dental data. Statistical and cartographical analyses are
provided for 498 living Eurasian populations. Mapping principal components
supplied evidence for the phene pool structure in Eurasian populations and for
reconstructions of our species history on the continent. The longitudinal
variability seems to be the most important regularity revealed by principal
components analysis (PCA) and mapping proving the division of the whole area
into western and eastern main provinces. So, the most ancient scenario in the
history of Eurasian populations was developing from two perspective different
groups: western group related to ancient populations of West Asia and the
eastern one rooted by ancestry in South and/or East Asia. In spite of the
enormous territory and the revealed divergence the populations of the continent
have undergone wide scale and intensive time-space interaction. Many details in
the revealed landscapes could be backgrounded to different historical events.
The most amazing results are obtained for proving migrations and assimilation
as two essential phenomena in Eurasian history: the wide spread of the western
combination through the whole continent till the Pacific coastline and the
envision of the movement of the paradox combinations of eastern and western
markers from South or Central Asia to the east and to the west. Taking into
account that no additional eastern combinations in the total variation in Asian
groups have been found but mixed or western markers' sets and that eastern
dental characteristics are traced in Asia since Homo erectus, the assumption is
made in favour of the hetero-level assimilation in the Eastern province and of
net-like evolution of our species.


  The random-force (Larkin) model of a directed elastic string subject to
quenched random forces in the transverse directions has been a paradigm in the
statistical physics of disordered systems. In this brief note, we investigate a
modified version of the above model where the total transverse force along the
polymer contour and the related total torque, in each realization of disorder,
vanish. We discuss the merits of adding these constraints and show that they
leave the qualitative behavior in the strong stretching regime unchanged, but
they reduce the effects of the random force by significant numerical
prefactors. We also show that a transverse random force effectively makes the
filament softer to compression by inducing undulations. We calculate the
related linear compression coefficient in both the usual and the constrained
random force model.


  We consider the standard three-component differential equation model for the
growth of an HIV virion population in an infected host in the absence of drug
therapy. The dynamical properties of the model are determined by the set of
values of six parameters which vary across host populations. There may be one
or two critical points whose natures play a key role in determining the outcome
of infection and in particular whether the HIV population will persist or
become extinct. There are two cases which may arise. In the first case, there
is only one critical point P_1 at biological values and this is an
asymptotically stable node. The system ends up with zero virions and so the
host becomes HIV-free. In the second case, there are two critical points P_1
and P_2 at biological values. Here P_1 is an unstable saddle point and P_2 is
an asymptotically stable spiral point with a non-zero virion level. In this
case the HIV population persists unless parameters change. We let the parameter
values take random values from distributions based on empirical data, but
suitably truncated, and determine the probabilities of occurrence of the
various combinations of critical points. From these simulations the probability
that an HIV infection will persist, across a population, is estimated. It is
found that with conservatively estimated distributions of parameters, within
the framework of the standard 3-component model, the chances that a within host
HIV population will become extinct is between 0.6% and 6.9%. With less
conservative parameter estimates, the probability is estimated to be as high as
24%. The many factors related to the transmission and possible spontaneous
elimination of the virus are discussed.


  We introduce a systematic classification method for the analogs of phase
transitions in finite systems. This completely general analysis, which is
applicable to any physical system and extends towards the thermodynamic limit,
is based on the microcanonical entropy and its energetic derivative, the
inverse caloric temperature. Inflection points of this quantity signal
cooperative activity and thus serve as distinct indicators of transitions. We
demonstrate the power of this method through application to the long-standing
problem of liquid-solid transitions in elastic, flexible homopolymers.


  We discuss the advantages of parallelization by multithreading on graphics
processing units (GPUs) for parallel tempering Monte Carlo computer simulations
of an exemplified bead-spring model for homopolymers. Since the sampling of a
large ensemble of conformations is a prerequisite for the precise estimation of
statistical quantities such as typical indicators for conformational
transitions like the peak structure of the specific heat, the advantage of a
strong increase in performance of Monte Carlo simulations cannot be
overestimated. Employing multithreading and utilizing the massive power of the
large number of cores on GPUs, being available in modern but standard graphics
cards, we find a rapid increase in efficiency when porting parts of the code
from the central processing unit (CPU) to the GPU.


  We define a subclass of Chemical Reaction Networks called Post-Translational
Modification systems. Important biological examples of such systems include
MAPK cascades and two-component systems which are well-studied experimentally
as well as theoretically. The steady states of such a system are solutions to a
system of polynomial equations with as many variables as equations. Even for
small systems the task of finding the solutions is daunting. We develop a
mathematical framework based on the notion of a cut, which provides a linear
elimination procedure to reduce the number of variables in the system. The
steady states are parameterized algebraically by a set of "core" variables, and
the non-negative steady states correspond to non-negative values of the core
variables. Further, minimal cuts are the connected components in the species
graph and provide conservation laws. A criterion for when a set of independent
conservation laws can be derived from cuts is given.


  A discrete model for computer simulations of the clustering dynamics of
Social Amoebae is presented. This model incorporates the wavelike propagation
of extracellular signaling cAMP, the sporadic firing of cells at early stage of
aggregation, the signal relaying as a response to stimulus, the inertia and
purposeful random walk of the cell movement. A Monte Carlo simulation is run
which shows the existence of potential equilibriums of mean and variance of
aggregation time. The simulation result of this model could well reproduce many
phenomena observed by actual experiments.


  Max Planck Institute for the Physics of Complex Systems, Noethnitzer Str. 38,
01187 Dresden, Germany


  The Minimum-Recombinant Haplotype Configuration problem (MRHC) has been
highly successful in providing a sound combinatorial formulation for the
important problem of genotype phasing on pedigrees. Despite several algorithmic
advances and refinements that led to some efficient algorithms, its
applicability to real datasets has been limited by the absence of some
important characteristics of these data in its formulation, such as mutations,
genotyping errors, and missing data.
  In this work, we propose the Haplotype Configuration with Recombinations and
Errors problem (HCRE), which generalizes the original MRHC formulation by
incorporating the two most common characteristics of real data: errors and
missing genotypes (including untyped individuals). Although HCRE is
computationally hard, we propose an exact algorithm for the problem based on a
reduction to the well-known Satisfiability problem. Our reduction exploits
recent progresses in the constraint programming literature and, combined with
the use of state-of-the-art SAT solvers, provides a practical solution for the
HCRE problem. Biological soundness of the phasing model and effectiveness (on
both accuracy and performance) of the algorithm are experimentally demonstrated
under several simulated scenarios and on a real dairy cattle population.


  For population genetics models with recombination, obtaining an exact,
analytic sampling distribution has remained a challenging open problem for
several decades. Recently, a new perspective based on asymptotic series has
been introduced to make progress on this problem. Specifically, closed-form
expressions have been derived for the first few terms in an asymptotic
expansion of the two-locus sampling distribution when the recombination rate
$\rho$ is moderate to large. In this paper, a new computational technique is
developed for finding the asymptotic expansion to an arbitrary order.
Computation in this new approach can be automated easily. Furthermore, it is
proved here that only a finite number of terms in the asymptotic expansion is
needed to recover (via the method of Pad\'{e} approximants) the exact two-locus
sampling distribution as an analytic function of $\rho$; this function is exact
for all values of $\rho\in[0,\infty)$. It is also shown that the new
computational framework presented here is flexible enough to incorporate
natural selection.


  The 3'-monofunctional adduct of cisplatin and d(CTCTG*G*TCTC)2 duplex DNA in
solvent with explicit counter ions and water molecules were subjected to MD-
simulation with AMBER force field on a nanosecond time scale. In order to
simulate the closure of the bond between the Pt and 5'-guanine-N7 atoms, the
forces acting between them were gradually increased during MD. After 500-800 ps
the transformation of the mono-adduct (straight DNA with the cisplatin residue
linked to one guanine-N7) to the bus-adduct (bent DNA where Pt atom is
connected through the N7 atoms of neighboring guanines) was observed. A cavity
between palatinate guanines is formed and filled with solvent molecules. The
rapid inclination of the center base pairs initiates a slow transition of the
whole molecule from the linear to the bent conformation. After about 1000-1300
ps a stable structure was reached, which is very similar to the one described
experimentally. The attractive force between the Pt- atom and the N7 of the
second guanine plays the main role in the large conformational changes induced
by formation of the adduct-adduct. X-N-Pt-N-torsions accelerate the bending but
a torsion force constant greater than 0.2 Kcal/mol lead to the breaking of the
H-bonds within the base pairs. The present study is the first dynamical
simulation that demonstrates in real time scale such a large conformational
perturbation of DNA.


  Many experimental studies have shown that the prion AGAAAAGA palindrome
hydrophobic region (113-120) has amyloid fibril forming properties and plays an
important role in prion diseases. However, due to the unstable, noncrystalline
and insoluble nature of the amyloid fibril, to date structural information on
AGAAAAGA region (113-120) has been very limited. This region falls just within
the N-terminal unstructured region PrP (1-123) of prion proteins. Traditional
X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy
experimental methods cannot be used to get its structural information. Under
this background, this paper introduces a novel approach of the canonical dual
theory to address the 3D atomic-resolution structure of prion AGAAAAGA amyloid
fibrils. The novel and powerful canonical dual computational approach
introduced in this paper is for the molecular modeling of prion AGAAAAGA
amyloid fibrils, and that the optimal atomic-resolution structures of prion
AGAAAAGA amyloid fibils presented in this paper are useful for the drive to
find treatments for prion diseases in the field of medicinal chemistry.
Overall, this paper presents an important method and provides useful
information for treatments of prion diseases. Overall, this paper could be of
interest to the general readership of Theoretical Biology.


  An Ising--like model of proteins is used to investigate the mechanical
unfolding of the Green Fluorescent Protein along different directions. When the
protein is pulled from its ends, we recover the major and minor unfolding
pathways observed in experiments. Upon varying the pulling direction, we find
the correct order of magnitude and ranking of the unfolding forces. Exploiting
the direction dependence of the unfolding force at equilibrium, we propose a
force sensor whose luminescence depends on the applied force.


  Population dynamics of individuals undergoing birth and death and diffusing
by short or long ranged twodimensional spatial excursions (Gaussian jumps or
L\'{e}vy flights) is studied. Competitive interactions are considered in a
global case, in which birth and death rates are influenced by all individuals
in the system, and in a nonlocal but finite-range case in which interaction
affects individuals in a neighborhood (we also address the noninteracting
case). In the global case one single or few-cluster configurations are achieved
with the spatial distribution of the bugs tied to the type of diffusion. In the
L\'{e}vy case long tails appear for some properties characterizing the shape
and dynamics of clusters. Under non-local finite-range interactions periodic
patterns appear with periodicity set by the interaction range. This length acts
as a cut-off limiting the influence of the long L\'{e}vy jumps, so that spatial
configurations under the two types of diffusion become more similar. By
dividing initially everyone into different families and following their descent
it is possible to show that mixing of families and their competition is greatly
influenced by the spatial dynamics.


  The dialects of Madagascar belong to the Greater Barito East group of the
Austronesian family and it is widely accepted that the Island was colonized by
Indonesian sailors after a maritime trek which probably took place around 650
CE. The language most closely related to Malagasy dialects is Maanyan but also
Malay is strongly related especially for what concerns navigation terms. Since
the Maanyan Dayaks live along the Barito river in Kalimantan (Borneo) and they
do not possess the necessary skill for long maritime navigation, probably they
were brought as subordinates by Malay sailors.
  In a recent paper we compared 23 different Malagasy dialects in order to
determine the time and the landing area of the first colonization. In this
research we use new data and new methods to confirm that the landing took place
on the south-east coast of the Island. Furthermore, we are able to state here
that it is unlikely that there were multiple settlements and, therefore,
colonization consisted in a single founding event.
  To reach our goal we find out the internal kinship relations among all the 23
Malagasy dialects and we also find out the different kinship degrees of the 23
dialects versus Malay and Maanyan. The method used is an automated version of
the lexicostatistic approach. The data concerning Madagascar were collected by
the author at the beginning of 2010 and consist of Swadesh lists of 200 items
for 23 dialects covering all areas of the Island. The lists for Maanyan and
Malay were obtained from published datasets integrated by author's interviews.


  We have developed a coarse-grained formulation for modeling the dynamic
behavior of cells quantitatively, based on stochasticity and heterogeneity,
rather than on biochemical reactions. We treat each reaction as a
continuous-time stochastic process, while reducing each biochemical quantity to
a binary value at the level of individual cells. The system can be analytically
represented by a finite set of ordinary linear differential equations, which
provides a continuous time course prediction of each molecular state. In this
letter, we introduce our formalism and demonstrate it with several examples.


  In this paper we study noise-induced bistability in a specific circuit with
many biological implications, namely a single-step enzymatic cycle described by
Michaelis Menten equations with quasi-steady state assumption. We study the
system both with a Master Equation formalism, and with the Fokker-Planck
continuous approximation, characterizing the conditions in which the continuous
approach is a good approximation of the exact discrete model. An analysis of
the stationary distribution in both cases shows that bimodality can not occur
in such a system. We discuss which additional requirements can generate
stochastic bimodality, by coupling the system with a chemical reaction
involving enzyme production and turnover. This extended system shows a bistable
behaviour only in specific parameter windows depending on the number of
molecules involved, providing hints about which should be a feasible system
size in order that such a phenomenon could be exploited in real biological
systems.


  Obtaining a closed-form sampling distribution for the coalescent with
recombination is a challenging problem. In the case of two loci, a new
framework based on asymptotic series has recently been developed to derive
closed-form results when the recombination rate is moderate to large. In this
paper, an arbitrary number of loci is considered and combinatorial approaches
are employed to find closed-form expressions for the first couple of terms in
an asymptotic expansion of the multi-locus sampling distribution. These
expressions are universal in the sense that their functional form in terms of
the marginal one-locus distributions applies to all finite- and
infinite-alleles models of mutation.


  A variety of neurodegenerative diseases are associated with the formation of
amyloid plaques. Our incomplete understanding of this process underscores the
need to decipher the principles governing protein aggregation. Most
experimental and simulation studies have been interpreted largely from the
perspective of proteins: the role of solvent has been relatively overlooked.
  In this Account, we provide a perspective on how interactions with water
affect folding landscapes of A$\beta$ monomers, A$\beta_{16-22}$ oligomer
formation, and protofilament formation in a Sup35 peptide. Simulations show
that the formation of aggregation-prone structures (N$^*$) similar to the
structure in the fibril requires overcoming high desolvation barrier. The
mechanism of protofilament formation in a polar Sup35 peptide fragment
illustrates that water dramatically slows down self-assembly. Release of water
trapped in the pores as water wires creates protofilament with a dry interface.
Similarly, one of the main driving force for addition of a solvated monomer to
a preformed fibril is the entropy gain of released water.
  We conclude by postulating that two-step model for protein crystallization
must also hold for higher order amyloid structure formation starting from
N$^*$. Multiple N$^*$ structures with varying water content results in a number
of distinct water-laden polymorphic structures. In predominantly hydrophobic
sequences, water accelerates fibril formation. In contrast, water-stabilized
metastable intermediates dramatically slow down fibril growth rates in
hydrophilic sequences.


  Prion diseases are invariably fatal and highly infectious neurodegenerative
diseases affecting humans and animals. The neurodegenerative diseases such as
Creutzfeldt-Jakob disease, variant Creutzfeldt-Jakob diseases,
Gerstmann-Str$\ddot{a}$ussler-Scheinker syndrome, Fatal Familial Insomnia, Kuru
in humans, scrapie in sheep, bovine spongiform encephalopathy (or 'mad-cow'
disease) and chronic wasting disease in cattle belong to prion diseases. By now
there have not been some effective therapeutic approaches to treat all these
prion diseases. Dogs, rabbits and horses were reported to be resistant to prion
diseases. By the end of year 2010 all the NMR structures of dog, rabbit and
horse prion proteins (X-ray for rabbits too) had been finished to release into
protein data bank. Thus, at this moment it is very worth studying the NMR and
X-ray molecular structures of horse, dog and rabbit prion proteins to obtain
insights into their immunity prion diseases.
  The author found that dog and horse prion proteins have stable molecular
dynamical structures whether under neutral or low pH environments, but rabbit
prion protein has stable molecular dynamical structures only under neutral pH
environment. Under low pH environment, the stable $\alpha$-helical molecular
structures of rabbit prion protein collapse into $\beta$-sheet structures. This
article focuses the studies on rabbit prion protein (within its C-terminal NMR,
Homology and X-ray molecular structured region RaPrP$^\text{C}$ (120-230)),
compared with human and mouse prion proteins (HuPrP$^\text{C}$ (125-228) and
MoPrP$^\text{C}$ (124-226) respectively). The author finds that some salt
bridges contribute to the structural stability of rabbit prion protein under
neutral pH environment.


  Chemical reactions inside cells occur in compartment volumes in the range of
atto- to femtolitres. Physiological concentrations realized in such small
volumes imply low copy numbers of interacting molecules with the consequence of
considerable fluctuations in the concentrations. In contrast, rate equation
models are based on the implicit assumption of infinitely large numbers of
interacting molecules, or equivalently, that reactions occur in infinite
volumes at constant macroscopic concentrations. In this article we compute the
finite-volume corrections (or equivalently the finite copy number corrections)
to the solutions of the rate equations for chemical reaction networks composed
of arbitrarily large numbers of enzyme-catalyzed reactions which are confined
inside a small sub-cellular compartment. This is achieved by applying a
mesoscopic version of the quasi-steady state assumption to the exact
Fokker-Planck equation associated with the Poisson Representation of the
chemical master equation. The procedure yields impressively simple and compact
expressions for the finite-volume corrections. We prove that the predictions of
the rate equations will always underestimate the actual steady-state substrate
concentrations for an enzyme-reaction network confined in a small volume. In
particular we show that the finite-volume corrections increase with decreasing
sub-cellular volume, decreasing Michaelis-Menten constants and increasing
enzyme saturation. The magnitude of the corrections depends sensitively on the
topology of the network. The predictions of the theory are shown to be in
excellent agreement with stochastic simulations for two types of networks
typically associated with protein methylation and metabolism.


  An important problem in bioinformatics is the inference of gene regulatory
networks (GRN) from temporal expression profiles. In general, the main
limitations faced by GRN inference methods is the small number of samples with
huge dimensionalities and the noisy nature of the expression measurements. In
face of these limitations, alternatives are needed to get better accuracy on
the GRNs inference problem. This work addresses this problem by presenting an
alternative feature selection method that applies prior knowledge on its search
strategy, called SFFS-BA. The proposed search strategy is based on the
Sequential Floating Forward Selection (SFFS) algorithm, with the inclusion of a
scale-free (Barab\'asi-Albert) topology information in order to guide the
search process to improve inference. The proposed algorithm explores the
scale-free property by pruning the search space and using a power law as a
weight for reducing it. In this way, the search space traversed by the SFFS-BA
method combines a breadth-first search when the number of combinations is small
(<k> <= 2) with a depth-first search when the number of combinations becomes
explosive (<k> >= 3), being guided by the scale-free prior information.
Experimental results show that the SFFS-BA provides a better inference
similarities than SFS and SFFS, keeping the robustness of the SFS and SFFS
methods, thus presenting very good results.


  Principal components (PCA) and hierarchical clustering are two of the most
heavily used techniques for analyzing the differences between nucleic acid
sequence samples sampled from a given environment. However, a classical
application of these techniques to distances computed between samples can lack
transparency because there is no ready interpretation of the axes of classical
PCA plots, and it is difficult to assign any clear intuitive meaning to either
the internal nodes or the edge lengths of trees produced by distance-based
hierarchical clustering methods such as UPGMA. We show that more interesting
and interpretable results are produced by two new methods that leverage the
special structure of phylogenetic placement data. Edge principal components
analysis enables the detection of important differences between samples that
contain closely related taxa. Each principal component axis is simply a
collection of signed weights on the edges of the phylogenetic tree, and these
weights are easily visualized by a suitable thickening and coloring of the
edges. Squash clustering outputs a (rooted) clustering tree in which each
internal node corresponds to an appropriate "average" of the original samples
at the leaves below the node. Moreover, the length of an edge is a suitably
defined distance between the averaged samples associated with the two incident
nodes, rather than the less interpretable average of distances produced by
UPGMA. We present these methods and illustrate their use with data from the
microbiome of the human vagina.


  Many contemporary studies have shown that astrocytes play a significant role
in modulating both short and long form of synaptic plasticity. There are very
few experimental models which elucidate the role of astrocyte over Long-term
Potentiation (LTP). Recently, Perea & Araque (2007) demonstrated a role of
astrocytes in induction of LTP at single hippocampal synapses. They suggested a
purely pre-synaptic basis for induction of this N-methyl-D- Aspartate (NMDA)
Receptor-independent LTP. Also, the mechanisms underlying this pre-synaptic
induction were not investigated. Here, in this article, we propose a
mathematical model for astrocyte modulated LTP which successfully emulates the
experimental findings of Perea & Araque (2007). Our study suggests the role of
retrograde messengers, possibly Nitric Oxide (NO), for this pre-synaptically
modulated LTP.


  The simplified Lennard-Jones (LJ) potential minimization problem is
$f(x)=4\sum_{i=1}^N \sum_{j=1,j<i}^N (\frac{1}{\tau_{ij}^6}
-\frac{1}{\tau_{ij}^3}) {subject to} x\in \mathbb{R}^n,$ where
$\tau_{ij}=(x_{3i-2}-x_{3j-2})^2 +(x_{3i-1}-x_{3j-1})^2 +(x_{3i} -x_{3j})^2$,
$(x_{3i-2},x_{3i-1},x_{3i})$ is the coordinates of atom $i$ in $\mathbb{R}^3$,
$i,j=1,2,...,N(\geq 2 \quad \text{integer})$, $n=3N$ and $N$ is the whole
number of atoms. The nonconvexity of the objective function and the huge number
of local minima, which is growing exponentially with $N$, interest many
mathematical optimization experts. In this paper, the canonical dual theory
elegantly tackles this problem illuminated by the amyloid fibril molecular
model building.
  Keywords: Mathematical Canonical Duality Theory $\cdot$ Mathematical
Optimization $\cdot$ Lennard-Jones Potential Minimization Problem $\cdot$
Global Optimization.


  Many types of bacteria can survive under stress by switching stochastically
between two different phenotypes: the "normals" who multiply fast, but are
vulnerable to stress, and the "persisters" who hardly multiply, but are
resilient to stress. Previous theoretical studies of such bacterial populations
have focused on the \emph{fitness}: the asymptotic rate of unbounded growth of
the population. Yet for an isolated population of established (and not very
large) size, a more relevant measure may be the population \emph{extinction
risk} due to the interplay of adverse extrinsic variations and intrinsic noise
of birth, death and switching processes. Applying a WKB approximation to the
pertinent master equation of such a two-population system, we quantify the
extinction risk, and find the most likely path to extinction under both
favorable and adverse conditions. Analytical results are obtained both in the
biologically relevant regime when the switching is rare compared with the birth
and death processes, and in the opposite regime of frequent switching. We show
that rare switches are most beneficial in reducing the extinction risk.


  The first experimental data from single-particle scattering experiments from
free electron lasers (FELs) are now becoming available. The first such
experiments are being performed on relatively large objects such as viruses,
which produce relatively low-resolution, low-noise diffraction patterns in
so-called "diffract-and-destroy" experiments. We describe a very simple test on
the angular correlations of measured diffraction data to determine if the
scattering is from an icosahedral particle. If this is confirmed, the efficient
algorithm proposed can then combine diffraction data from multiple shots of
particles in random unknown orientations to generate a full 3D image of the
icosahedral particle. We demonstrate this with a simulation for the satellite
tobacco necrosis virus (STNV), the atomic coordinates of whose asymmetric unit
is given in Protein Data Bank entry 2BUK.


  We study a model for the translocation of proteins across membranes through a
nanopore using a ratcheting mechanism. When the protein enters the nanopore it
diffuses in and out of the pore according to a Brownian motion. Moreover, it is
bound by ratcheting molecules which hinder the diffusion of the protein out of
the nanopore, i.e. the Brownian motion is reflected such that no ratcheting
molecule exits the pore. New ratcheting molecules bind at rate gamma. Extending
our previous approach (Depperschmidt and Pfaffelhuber, 2010) we allow the
ratcheting molecules to dissociate (at rate delta) from the protein (Model I).
We also provide an approximate model (Model II) which assumes a Poisson
equilibrium of ratcheting molecules on one side of the current reflection
boundary. Using analytical methods and simulations we show that the speed of
both models are approximately the same. Our analytical results on Model II give
the speed of translocation by means of a solution of an ordinary differential
equation.


  We propose a simple model of network co-evolution in a game-dynamical system
of interacting agents that play repeated games with their neighbors, and adapt
their behaviors and network links based on the outcome of those games. The
adaptation is achieved through a simple reinforcement learning scheme. We show
that the collective evolution of such a system can be described by
appropriately defined replicator dynamics equations. In particular, we suggest
an appropriate factorization of the agents' strategies that results in a
coupled system of equations characterizing the evolution of both strategies and
network structure, and illustrate the framework on two simple examples.


  Many theoretical and experimental studies suggest that range expansions can
have severe consequences for the gene pool of the expanding population. Due to
strongly enhanced genetic drift at the advancing frontier, neutral and weakly
deleterious mutations can reach large frequencies in the newly colonized
regions, as if they were surfing the front of the range expansion. These
findings raise the question of how frequently beneficial mutations successfully
surf at shifting range margins, thereby promoting adaptation towards a
range-expansion phenotype. Here, we use individual-based simulations to study
the surfing statistics of recurrent beneficial mutations on wave-like range
expansions in linear habitats. We show that the rate of surfing depends on two
strongly antagonistic factors, the probability of surfing given the spatial
location of a novel mutation and the rate of occurrence of mutations at that
location. The surfing probability strongly increases towards the tip of the
wave. Novel mutations are unlikely to surf unless they enjoy a spatial head
start compared to the bulk of the population. The needed head start is shown to
be proportional to the inverse fitness of the mutant type, and only weakly
dependent on the carrying capacity. The second factor is the mutation
occurrence which strongly decreases towards the tip of the wave. Thus, most
successful mutations arise at an intermediate position in the front of the
wave. We present an analytic theory for the tradeoff between these factors that
allows to predict how frequently substitutions by beneficial mutations occur at
invasion fronts. We find that small amounts of genetic drift increase the
fixation rate of beneficial mutations at the advancing front, and thus could be
important for adaptation during species invasions.


  A major issue in biology is the understanding of the interactions between
proteins. These interactions can be described by a network, where the proteins
are modeled by nodes and the interactions by edges. The origin of these protein
networks is not well understood yet. Here we present a two-step model, which
generates clusters with the same topological properties as networks for
protein-protein interactions, namely, the same degree distribution, cluster
size distribution, clustering coefficient and shortest path length. The
biological and model networks are not scale free but exhibit small world
features. The model allows the fitting of different biological systems by
tuning a single parameter.


  We put forward the idea of establishing a novel interdisciplinary field of
research at the interface between quantum mechanics and proteomics. The new
field, called quantum proteomics, is defined as the large-scale study of the
electronic structure of the proteins that define an organism's proteome. The
electronic structure of proteins is unveiled with the aid of linear-scaling
quantum mechanical calculations. Such calculations provide information about
the energy levels of the proteins, the charges of their amino acid side chains,
their electrostatic potentials and permanent dipole moments ({\mu}). Since the
magnitude of the electric dipole moment of any protein is not null ({\mu}\neq0
Debye), the dipole moment can be employed to characterize the electronic
structure of each protein that belongs to an organism's proteome. As an
example, we investigate six proteins from the thermophilic bacterium
Methanobacterium thermoautotrophicum (Mth) whose atomic structures were
characterized by solution NMR spectroscopy.


  Next-generation sequencing (NGS) technologies allow new methodologies for
alternative splicing (AS) analysis. Current computational methods for AS from
NGS data are mainly focused on predicting splice site junctions or de novo
assembly of full-length transcripts. These methods are computationally
expensive and produce a huge number of full-length transcripts or splice
junctions, spanning the whole genome of organisms. Thus summarizing such data
into the different gene structures and AS events of the expressed genes is an
hard task.
  To face this issue in this paper we investigate the computational problem of
reconstructing from NGS data, in absence of the genome, a gene structure for
each gene that is represented by the isoform graph: we introduce such graph and
we show that it uniquely summarizes the gene transcripts. We define the
computational problem of reconstructing the isoform graph and provide some
conditions that must be met to allow such reconstruction.
  Finally, we describe an efficient algorithmic approach to solve this problem,
validating our approach with both a theoretical and an experimental analysis.


  Mutation rate variation across loci is well known to cause difficulties,
notably identifiability issues, in the reconstruction of evolutionary trees
from molecular sequences. Here we introduce a new approach for estimating
general rates-across-sites models. Our results imply, in particular, that large
phylogenies are typically identifiable under rate variation. We also derive
sequence-length requirements for high-probability reconstruction.
  Our main contribution is a novel algorithm that clusters sites according to
their mutation rate. Following this site clustering step, standard
reconstruction techniques can be used to recover the phylogeny. Our results
rely on a basic insight: that, for large trees, certain site statistics
experience concentration-of-measure phenomena.


  Discrete models have a long tradition in engineering, including finite state
machines, Boolean networks, Petri nets, and agent-based models. Of particular
importance is the question of how the model structure constrains its dynamics.
This paper discusses an algebraic framework to study such questions. The
systems discussed here are given by mappings on an affine space over a finite
field, whose coordinate functions are polynomials. They form a general class of
models which can represent many discrete model types. Assigning to such a
system its dependency graph, that is, the directed graph that indicates the
variable dependencies, provides a mapping from systems to graphs. A basic
property of this mapping is derived and used to prove that dynamical systems
with an acyclic dependency graph can only have a unique fixed point in their
phase space and no periodic orbits. This result is then applied to a published
model of in vitro virus competition.


  If in the sexual Penna ageing model conditions are applied leading to
complementary bit-strings, then marriages between brothers and sisters, or
between close cousins, may lead to more offspring than for unrelated couples.


  A ~ 62 My periodicity in fossil biodiversity has been observed in independent
studies of paleontology databases over ~0.5Gy. The period and phase of this
biodiversity cycle coincides with the oscillation of our solar system normal to
the galactic disk with an amplitude ~70 parsecs and a period ~64 My. Our Galaxy
is falling toward the Virgo cluster, forming a galactic shock at the north end
of our galaxy due to this motion, capable of accelerating particles and
exposing our galaxy's northern side to a higher flux of cosmic rays. These
high-energy particles strike the Earth's atmosphere initiating extensive air
showers, ionizing the atmosphere by producing charged secondary particles.
Secondary particles such as muons produced as a result of nuclear interactions
are able to reach the ground and enhance the biological radiation dose. Using a
Monte Carlo simulation package CORSIKA, we compute the biological dose
resulting from enhanced muon exposure from cosmic rays and discuss their
implications for terrestrial biodiversity variations.


  Microtubule dynamics is largely influenced by nucleotide hydrolysis and the
resultant tubulin configuration changes. The GTP cap model has been proposed to
interpret the stabilizing mechanism of microtubule growth from the view of
hydrolysis effects. Besides, the microtubule growth involves the closure of a
curved sheet at its growing end. The curvature conversion also helps to
stabilize the successive growth, and the curved sheet is referred to as the
conformational cap. However, there still lacks theoretical investigation on the
mechanical-chemical coupling growth process of microtubules. In this paper, we
study the growth mechanisms of microtubules by using a coarse-grained molecular
method. Firstly, the closure process involving a sheet-to-tube transition is
simulated. The results verify the stabilizing effect of the sheet structure,
and the minimum conformational cap length that can stabilize the growth is
demonstrated to be two dimers. Then, we show that the conformational cap can
function independently of the GTP cap, signifying the pivotal role of
mechanical factors. Furthermore, based on our theoretical results, we describe
a Tetris-like growth style of microtubules: the stochastic tubulin assembly is
regulated by energy and harmonized with the seam zipping such that the sheet
keeps a practically constant length during growth.


  Using Langevin dynamics simulations, we investigate the dynamics of
chaperone-assisted translocation of a flexible polymer through a nanopore. We
find that increasing the binding energy $\epsilon$ between the chaperone and
the chain and the chaperone concentration $N_c$ can greatly improve the
translocation probability. Particularly, with increasing the chaperone
concentration a maximum translocation probability is observed for weak binding.
For a fixed chaperone concentration, the histogram of translocation time $\tau$
has a transition from long-tailed distribution to Gaussian distribution with
increasing $\epsilon$. $\tau$ rapidly decreases and then almost saturates with
increasing binding energy for short chain, however, it has a minimum for longer
chains at lower chaperone concentration. We also show that $\tau$ has a minimum
as a function of the chaperone concentration. For different $\epsilon$, a
nonuniversal dependence of $\tau$ on the chain length $N$ is also observed.
These results can be interpreted by characteristic entropic effects for
flexible polymers induced by either crowding effect from high chaperone
concentration or the intersegmental binding for the high binding energy.


  Tumours are made up of a mixed population of different types of cells that
include normal struc- tures as well as ones associated with the malignancy, and
there are multiple interactions between the malignant cells and the local
microenvironment. These intercellular interactions, modulated by the
microenvironment, effect tumour progression and represent a largely under
appreciated therapeutic target. We use observations of primary tumor biology
from prostate cancer to extrapolate a math- ematical model: specifically; it
has been observed that in prostate cancer three disparate cellular outcomes
predominate: (i) the tumour remains well differentiated and clinically indolent
- in this case the local stromal cells may act to restrain the growth of the
cancer; (ii) early in its genesis the tumour acquires a highly malignant
phenotype, growing rapidly and displacing the original stromal population
(often referred to as small cell prostate cancer) - these less common
aggressive tumours are relatively independent of the local microenvironment;
and, (iii) the tumour co-opts the local stroma - taking on a classic
stromagenic phenotype where interactions with the local microenviron- ment are
critical to the cancer growth. We present an evolutionary game theoretical
construct that models the influence of tumour-stroma interactions in driving
these outcomes. We consider three characteristic and distinct cellular
populations: stromal cells, tumour cells that are self-reliant in terms of
microenvironmental factors and tumour cells that depend on the environment for
resources but can also co-opt stroma. Using evolutionary game theory we explore
a number of different sce- narios that elucidate the impact of tumour-stromal
interactions on the dynamics of prostate cancer growth and progression and how
different treatments in the metastatic setting can affect different types of
tumors.


  Cytotoxic T lymphocytes (CTLs) are immune system cells that are thought to
play an important role in controlling HIV infection. We develop a stochastic
ODE model of HIV-CTL interaction that extends current deterministic ODE models.
Based on this stochastic model, we consider the effect of CTL attack on
intrahost HIV lineages assuming CTLs attack several epitopes with equal
strength. In this setting, we introduce a limiting version of our stochastic
ODE under which we show that the coalescence of HIV lineages can be described
by a simple paintbox construction. Through numerical experiments, we show that
our results under the limiting stochastic ODE accurately reflect HIV lineages
under CTL attack when the HIV population size is on the low end of its
hypothesized range. Current techniques of HIV lineage construction depend on
the Kingman coalescent. Our results give an explicit connection between CTL
attack and HIV lineages.


  Statistical properties of environments experienced by biological signaling
systems in the real world change, which necessitate adaptive responses to
achieve high fidelity information transmission. One form of such adaptive
response is gain control. Here we argue that a certain simple mechanism of gain
control, understood well in the context of systems neuroscience, also works for
molecular signaling. The mechanism allows to transmit more than one bit (on or
off) of information about the signal independently of the signal variance. It
does not require additional molecular circuitry beyond that already present in
many molecular systems, and, in particular, it does not depend on existence of
feedback loops. The mechanism provides a potential explanation for abundance of
ultrasensitive response curves in biological regulatory networks.


  In evolution, the effects of a single deleterious mutation can sometimes be
compensated for by a second mutation which recovers the original phenotype.
Such epistatic interactions have implications for the structure of genome space
- namely, that networks of genomes encoding the same phenotype may not be
connected by single mutational moves. We use the folding of RNA sequences into
secondary structures as a model genotype-phenotype map and explore the neutral
spaces corresponding to networks of genotypes with the same phenotype. In most
of these networks, we find that it is not possible to connect all genotypes to
one another by single point mutations. Instead, a network for a phenotypic
structure with $n$ bonds typically fragments into at least $2^n$ neutral
components, often of similar size. While components of the same network
generate the same phenotype, they show important variations in their
properties, most strikingly in their evolvability and mutational robustness.
This heterogeneity implies contingency in the evolutionary process.


  A molecular understanding of how protein function is related to protein
structure will require an ability to understand large conformational changes
between multiple states. Unfortunately these states are often separated by high
free energy barriers and within a complex energy landscape. This makes it very
difficult to reliably connect, for example by all-atom molecular dynamics
calculations, the states, their energies and the pathways between them. A major
issue needed to improve sampling on the intermediate states is an order
parameter -- a reduced descriptor for the major subset of degrees of freedom --
that can be used to aid sampling for the large conformational change. We
present a novel way to combine information from molecular dynamics using
non-linear time series and dimensionality reduction, in order to quantitatively
determine an order parameter connecting two large-scale conformationally
distinct protein states. This new method suggests an implementation for
molecular dynamics calculations that may dramatically enhance sampling of
intermediate states.


  Single-molecule force spectroscopy has proven to be a powerful tool for
studying the kinetic behavior of biomolecules. Through application of an
external force, conformational states with small or transient populations can
be stabilized, allowing them to be characterized and the statistics of
individual trajectories studied to provide insight into biomolecular folding
and function. Because the observed quantity (force or extension) is not
necessarily an ideal reaction coordinate, individual observations cannot be
uniquely associated with kinetically distinct conformations. While
maximum-likelihood schemes such as hidden Markov models have solved this
problem for other classes of single-molecule experiments by using temporal
information to aid in the inference of a sequence of distinct conformational
states, these methods do not give a clear picture of how precisely the model
parameters are determined by the data due to instrument noise and finite-sample
statistics, both significant problems in force spectroscopy. We solve this
problem through a Bayesian extension that allows the experimental uncertainties
to be directly quantified, and build in detailed balance to further reduce
uncertainty through physical constraints. We illustrate the utility of this
approach in characterizing the three-state kinetic behavior of an RNA hairpin
in a stationary optical trap.


  Under the effects of strong genetic drift, it is highly probable to observe
gene fixation or loss in a population, shown by divergent probability density
functions, or infinite adaptive peaks on a landscape. It is then interesting to
ask what such infinite peaks imply, with or without combining other biological
factors (e.g. mutation and selection). We study the stochastic escape time from
the generated infinite adaptive peaks, and show that Kramers' classical escape
formula can be extended to the non-Gaussian distribution cases. The constructed
landscape provides a global description for system's middle and long term
behaviors, breaking the constraints in previous methods.


  A model for artificial evolution in native x86 Windows systems has been
developed at the end of 2010. In this text, further improvements and additional
analogies to natural microbiologic processes are presented. Several experiments
indicate the capability of the system - and raise the question of possible
countermeasures.


  In this article we propose a new method, which we name 'quartet neighbor
joining', or 'quartet-NJ', to infer an unrooted species tree on a given set of
taxa T from empirical distributions of unrooted quartet gene trees on all
four-taxon subsets of T. In particular, quartet-NJ can be used to estimate a
species tree on T from distributions of gene trees on T. The quartet-NJ
algorithm is conceptually very similar to classical neighbor joining, and its
statistical consistency under the multispecies coalescent model is proved by a
variant of the classical 'cherry picking'-theorem. In order to demonstrate the
suitability of quartet-NJ, coalescent processes on two different species trees
(on five resp. nine taxa) were simulated, and quartet-NJ was applied to the
simulated gene tree distributions. Further, quartet-NJ was applied to quartet
distributions obtained from multiple sequence alignments of 28 proteins of nine
prokaryotes.


  The distribution and heritability of many traits depends on numerous loci in
the genome. In general, the astronomical number of possible genotypes makes the
system with large numbers of loci difficult to describe. Multilocus evolution,
however, greatly simplifies in the limit of weak selection and frequent
recombination. In this limit, populations rapidly reach Quasi-Linkage
Equilibrium (QLE) in which the dynamics of the full genotype distribution,
including correlations between alleles at different loci, can be parameterized
by the allele frequencies. This review provides a simplified exposition of the
concept and mathematics of QLE which is central to the statistical description
of genotypes in sexual populations. We show how key results of Quantitative
Genetics such as the generalized Fisher's "Fundamental Theorem", along with
Wright's Adaptive Landscape, emerge within QLE from the dynamics of the
genotype distribution. We then discuss under what circumstances QLE is
applicable, and what the breakdown of QLE implies for the population structure
and the dynamics of selection. Understanding of the fundamental aspects of
multilocus evolution obtained through simplified models may be helpful in
providing conceptual and computational tools to address the challenges arising
in the studies of complex quantitative phenotypes of practical interest.


  Large populations may contain numerous simultaneously segregating
polymorphisms subject to natural selection. Since selection acts on individuals
whose fitness depends on many loci, different loci affect each other's
dynamics. This leads to stochastic fluctuations of allele frequencies above and
beyond genetic drift - an effect known as genetic draft. Since recombination
disrupts associations between alleles, draft is strong when recombination is
rare. Here, we study a facultatively outcrossing population in a regime where
the frequency of out-crossing and recombination, r, is small compared to the
characteristic scale of fitness differences \sigma. In this regime, fit
genotypes expand clonally, leading to large fluctuations in the number of
recombinant offspring genotypes. The power law tail in the distribution of the
latter makes it impossible to capture the dynamics of draft by an effective
neutral model. Instead, we find that the fixation time of a neutral allele
increases only slowly with the population size but depends sensitively on the
ratio r/\sigma. The efficacy of selection is reduced dramatically and alleles
behave "quasi-neutrally" even for Ns>> 1, provided that |s|< s_c, where s_c
depends strongly on r/\sigma, but only weakly on population size N. In
addition, the anomalous fluctuations due to draft change the spectrum of
(quasi)-neutral alleles from f(\nu)\sim 1/\nu, corresponding to drift, to
\sim1/\nu^2. Finally, draft accelerates the rate of two step adaptations
through deleterious intermediates.


  We introduce the use of fractional calculus, i.e., the use of integrals and
derivatives of non-integer (arbitrary) order, in epidemiology. The proposed
approach is illustrated with an outbreak of dengue disease, which is motivated
by the first dengue epidemic ever recorded in the Cape Verde islands off the
coast of west Africa, in 2009. Numerical simulations show that in some cases
the fractional models fit better the reality when compared with the standard
differential models. The classical results are obtained as particular cases by
considering the order of the derivatives to take an integer value.


  Formation of a pseudoknot in the conserved RNA core domain in the
ribonucleoprotein human telomerase is required for function. In vitro
experiments show that the pseudoknot (PK) is in equilibrium with an extended
hairpin (HP) structure. We use molecular simulations of a coarse-grained model,
which reproduces most of the salient features of the experimental melting
profiles of PK and HP, to show that crowding enhances the stability of PK
relative to HP in the wild type and in a mutant associated with dyskeratosis
congenita. In monodisperse suspensions, small crowding particles increase the
stability of compact structures to a greater extent than larger crowders. If
the sizes of crowders in a binary mixture are smaller than the unfolded RNA,
the increase in melting temperature due to the two components is additive. In a
ternary mixture of crowders that are larger than the unfolded RNA, which mimics
the composition of ribosome, large enzyme complexes and proteins in E. coli,
the marginal increase in stability is entirely determined by the smallest
component. We predict that crowding can restore partially telomerase activity
in mutants, which dramatically decrease the PK stability.


  The imaginary time path integral formalism is applied to a nonlinear
Hamiltonian for a short fragment of heterogeneous DNA with a stabilizing
solvent interaction term. Torsional effects are modeled by a twist angle
between neighboring base pairs stacked along the molecule backbone. The base
pair displacements are described by an ensemble of temperature dependent paths
thus incorporating those fluctuational effects which shape the multisteps
thermal denaturation. By summing over $\sim 10^7 - 10^8$ base pair paths, a
large number of double helix configurations is taken into account consistently
with the physical requirements of the model potential. The partition function
is computed as a function of the twist. It is found that the equilibrium twist
angle, peculiar of B-DNA at room temperature, yields the stablest helicoidal
geometry against thermal disruption of the base pair hydrogen bonds. This
result is corroborated by the computation of thermodynamical properties such as
fractions of open base pairs and specific heat.


  We investigate the impact of cyclic competition on pattern formation in the
rock-paper-scissors game. By separately considering random and prepared initial
conditions, we observe a critical influence of the competition rate $p$ on the
stability of spiral waves and on the emergence of biodiversity. In particular,
while increasing values of $p$ promote biodiversity, they may act detrimental
on spatial pattern formation. For random initial conditions, we observe a phase
transition from biodiversity to an absorbing phase, whereby the critical value
of mobility grows linearly with increasing values of $p$ on a log-log scale,
but then saturates as $p$ becomes large. For prepared initial conditions, we
observe the formation of single-armed spirals, but only for values of $p$ that
are below a critical value. Once above, the spirals break up and form
disordered spatial structures, mainly because of the percolation of vacant
sites. Thus, there exists a critical value of the competition rate $p_{c}$ for
stable single-armed spirals in finite populations. Importantly though, $p_{c}$
increases with increasing system size, because noise reinforces the
disintegration of ordered patterns. In addition, we also find that $p_{c}$
increases with the mobility. These phenomena are reproduced by a deterministic
model that is based on nonlinear partial differential equations. Our findings
indicate that competition is vital for the sustenance of biodiversity and
emergence of pattern formation in ecosystems governed by cyclical interactions.


  In disease mapping, the aim is to estimate the spatial pattern in disease
risk over an extended geographical region, so that areas with elevated risks
can be identified. A Bayesian hierarchical approach is typically used to
produce such maps, which models the risk surface with a set of spatially smooth
random effects. However, in complex urban settings there are likely to be
boundaries in the risk surface, which separate populations that are
geographically adjacent but have very different risk profiles. Therefore this
paper proposes an approach for detecting such risk boundaries, and tests its
effectiveness by simulation. Finally, the model is applied to lung cancer
incidence data in Greater Glasgow, Scotland, between 2001 and 2005.


  In this study, through phenomenological comparison of the velocity-force data
of processive motor proteins, including conventional kinesin, cytoplasmic
dynein and myosin V, we found that, the ratio between motor velocities of two
different ATP concentrations is almost invariant for any substall, superstall
or negative external loads. Therefore, the velocity of motor can be well
approximated by a Michaelis-Menten like formula $V=\atp k(F)L/(\atp +K_M)$,
with $L$ the step size, and $k(F)$ the external load $F$ dependent rate of one
mechanochemical cycle of motor motion in saturated ATP solution. The difference
of Michaelis-Menten constant $K_M$ for substall, superstall and negative
external load indicates, the ATP molecule affinity of motor head for these
three cases are different, though the expression of $k(F)$ as a function of $F$
might be unchanged for any external load $F$. Verifications of this
Michaelis-Menten like formula has also been done by fitting to the recent
experimental data.


  A new numerical characterization of symbolic sequences is proposed. The
partition of sequence based on Ke and Tong algorithm is a starting point.
Algorithm decomposes original sequence into set of distinct subsequences - a
patterns. The set of subsequences common for two symbolic sequences (their
intersection) is proposed as a measure of similarity between them. The new
similarity measure works well for short (of tens letters) sequences and the
very long (of hundred thousand letters) as well. When applied to nucleotide or
protein sequences may help to trace possible evolutionary of species. As an
illustration, similarity of several sets of nucleotide and amino acid sequences
is examined.


  Traditional studies about disease dynamics have focused on global stability
issues, due to their epidemiological importance. We study a classical SIR-SI
model for arboviruses in two different directions: we begin by describing an
alternative proof of previously known global stability results by using only a
Lyapunov approach. In the sequel, we take a different view and we argue that
vectors and hosts can have very distinctive intrinsic time-scales, and that
such distinctiveness extends to the disease dynamics. Under these hypothesis,
we show that two asymptotic regimes naturally appear: the fast host dynamics
and the fast vector dynamics. The former regime yields, at leading order, a SIR
model for the hosts, but with a rational incidence rate. In this case, the
vector disappears from the model, and the dynamics is similar to a directly
contagious disease. The latter yields a SI model for the vectors, with the
hosts disappearing from the model. Numerical results show the performance of
the approximation, and a rigorous proof validates the reduced models.


  Ecological systems comprise an astonishing diversity of species that
cooperate or compete with each other forming complex mutual dependencies. The
minimum requirements to maintain a large species diversity on long time scales
are in general unknown. Using lichen communities as an example, we propose a
model for the evolution of mutually excluding organisms that compete for space.
We suggest that chain-like or cyclic invasions involving three or more species
open for creation of spatially separated sub-populations that subsequently can
lead to increased diversity. In contrast to its non-spatial counterpart, our
model predicts robust co-existence of a large number of species, in accordance
with observations on lichen growth. It is demonstrated that large species
diversity can be obtained on evolutionary timescales, provided that
interactions between species have spatial constraints. In particular, a phase
transition to a sustainable state of high diversity is identified.


  Variability in the tumorigenic potential among cancer cells within a tumor
population is an unresolved fundamental issue in cancer biology. It is
important to know whether cancer cells with higher tumorigenic potential, such
as cancer stem cells, are only a small subpopulation. We attempt to address the
question of variability in tumorigenic potential based on the heterogeneity in
apoptosis resistance of cancer cells. We use stochastic differential equations
and kinetic Monte Carlo simulations to explore the mechanisms that generate
cell-to-cell variability in apoptosis resistance of cancer cells. In our model,
a simplified scheme of apoptosis signaling reactions is developed focusing on
the proapoptotic Bid-Bax reaction and its inhibition by Bcl-2 like
antiapoptotic proteins. We show how a combination of low probability Bid-Bax
reaction along with overexpressed reactant molecules allows specific killing of
cancer cells, especially under targeted therapy such as Bcl-2 inhibition. This
low probability Bid-Bax reaction protects normal cells from accidental
apoptosis but generates cell-to-cell stochastic variability in apoptotic
activation of cells equipped with overexpressed Bid and Bax molecules. We
further demonstrate that cellular variations in Bcl-2 / Bax ratio, within a
cancer cell population, can affect the intrinsic fluctuations arising from the
stochastic Bid-Bax reaction and thereby provides a mechanism for origin of
cells with higher tumorigenic potential. We discuss the implications of our
results for cancer therapy, such as, optimal strategies to minimize stochastic
fluctuations in cancer cell death.


  While seemingly straightforward in principle, the reliable estimation of rate
constants is seldom easy in practice. Numerous issues, such as the complication
of poor reaction coordinates, cause obvious approaches to yield unreliable
estimates. When a reliable order parameter is available, the reactive flux
theory of Chandler allows the rate constant to be extracted from the plateau
region of an appropriate reactive flux function. However, when applied to real
data from single-molecule experiments or molecular dynamics simulations, the
rate can sometimes be difficult to extract due to the numerical differentiation
of a noisy empirical correlation function or difficulty in locating the plateau
region at low sampling frequencies. We present a modified version of this
theory which does not require numerical derivatives, allowing rate constants to
be robustly estimated from the time-correlation function directly. We compare
these approaches using single-molecule force spectroscopy measurements of an
RNA hairpin.


  We describe how noise propagates through a network by calculating the
variance of the outputs. Using stochastic calculus and dynamical systems
theory, we study the network topologies that accentuate or alleviate the effect
of random variance in the network for both directed and undirected graphs.
Given a linear tree network, the variance in the output is a convex function of
the poles of the individual nodes. Cycles create correlations which in turn
increase the variance in the output. Feedforward and feedback have a limited
effect on noise propagation when the respective cycles is sufficiently long.
Crosstalk between the elements of different pathways helps reduce the output
noise, but makes the network slower. Next, we study the differences between
disturbances in the inputs and disturbances in the network parameters, and how
they propagate to the outputs. Finally, we show how noise correlations can
affect the steady state of the system in chemical reaction networks with
reactions of two or more reactants, each of which may be affected by
independent or correlated noise sources.


  In the many transitions from foraging to agropastoralism it is debated
whether the primary drivers are innovations in technology or increases of
population. The driver discussion traditionally separates Malthusian
(technology driven) from Boserupian (population driven) theories. I present a
numerical model of the transitions to agriculture and discuss this model in the
light of the population versus technology debate and in Boserup's analytical
framework in development theory. Although my model is based on ecological
-Neomalthusian- principles, the coevolutionary positive feedback relationship
between technology and population results in a seemingly Boserupian response:
innovation is greatest when population pressure is highest. This outcome is not
only visible in the theory-driven reduced model, but is also present in a
corresponding "real world" simulator which was tested against archaeological
data, demonstrating the relevance and validity of the coevolutionary model. The
lesson to be learned is that not all that acts Boserupian needs Boserup at its
core.


  Theoretical frameworks to estimate the tolerance of metabolic networks to
various failures are important to evaluate the robustness of biological complex
systems in systems biology. In this paper, we focus on a measure for robustness
in metabolic networks, namely, the impact degree, and propose an approximation
method to predict the probability distribution of impact degrees from metabolic
network structures using the theory of branching process. We demonstrate the
relevance of this method by testing it on real-world metabolic networks.
Although the approximation method possesses a few limitations, it may be a
powerful tool for evaluating metabolic robustness.


  We study the interplay of population growth and evolutionary dynamics using a
stochastic model based on birth and death events. In contrast to the common
assumption of an independent population size, evolution can be strongly
affected by population dynamics in general. Especially for fast reproducing
microbes which are subject to selection, both types of dynamics are often
closely intertwined. We illustrate this by considering different growth
scenarios. Depending on whether microbes die or stop to reproduce (dormancy),
qualitatively different behaviors emerge. For cooperating bacteria, a permanent
increase of costly cooperation can occur. Even if not permanent, cooperation
can still increase transiently due to demographic fluctuations. We validate our
analysis via stochastic simulations and analytic calculations. In particular,
we derive a condition for an increase in the level of cooperation.


  At the molecular level, most biological processes entail protein associations
which in turn rely on a small fraction of interfacial residues called hot
spots. Here we show that hot spots share a unifying molecular attribute: they
provide a third-body contribution to intermolecular cooperativity. Such motif,
based on the wrapping of interfacial electrostatic interactions, is essential
to maintain the integrity of the interface and can be exploited in rational
drug design since such regions may serve as blueprints to engineer small
molecules disruptive of protein-protein interfaces.


  We present new and improved fixed-parameter algorithms for computing maximum
agreement forests (MAFs) of pairs of rooted binary phylogenetic trees. The size
of such a forest for two trees corresponds to their subtree prune-and-regraft
distance and, if the agreement forest is acyclic, to their hybridization
number. These distance measures are essential tools for understanding
reticulate evolution. Our algorithm for computing maximum acyclic agreement
forests is the first depth-bounded search algorithm for this problem. Our
algorithms substantially outperform the best previous algorithms for these
problems.


  We have previously shown an analysis of our dimer model in the over-damped
regime to show directed transport in equilibrium. Here we analyze the full
model with inertial terms present to establish the same result. First we derive
the Fokker-Planck equation for the system following a Galilean transformation
to show that a uniformly translating equilibrium distribution is possible.
Then, we find out the velocity selection for the centre of mass motion using
that distribution on our model. We suggest generalization of our calculations
for soft collision potentials and indicate to interesting situation with
possibility of oscillatory non-equilibrium state within equilibrium.


  Various studies have shown an association between single nucleotide
polymorphisms (SNPs) and common disease. We hypothesize that information
encoded in the structure of SNP haploblock variation illumines molecular
pathways and cellular mechanisms involved in the regulation of host adaptation
to the environment. We developed and utilized the normalized information
content (NIC), a novel metric based on SNP haploblock variation. We found that
all SNP haploblocks with statistically low information content contained
putative transcription factor binding sites and microRNA motifs. We were able
to translate a biophysical, mathematical measure of common variants into a
deeper understanding of the life sciences through analysis of biochemical
patterns associated with SNP haploblock variation. We submit that this new
metric, NIC, may be useful in decoding the functional significance of common
variation in the human genome and in analyzing the regulation of molecular
pathways involved in host adaptation to environmental pathogens.


  We present a SIR+ASI epidemic model to describe the interaction between human
and dengue fever mosquito populations. A control strategy in the form of
vaccination, to decrease the number of infected individuals, is used. An
optimal control approach is applied in order to find the best way to fight the
disease.


  The reconstruction of phylogenies from DNA or protein sequences is a major
task of computational evolutionary biology. Common phenomena, notably
variations in mutation rates across genomes and incongruences between gene
lineage histories, often make it necessary to model molecular data as
originating from a mixture of phylogenies. Such mixed models play an
increasingly important role in practice. Using concentration of measure
techniques, we show that mixtures of large trees are typically identifiable. We
also derive sequence-length requirements for high-probability reconstruction.


  Adaptation often involves the acquisition of a large number of genomic
changes which arise as mutations in single individuals. In asexual populations,
combinations of mutations can fix only when they arise in the same lineage, but
for populations in which genetic information is exchanged, beneficial mutations
can arise in different individuals and be combined later. In large populations,
when the product of the population size N and the total beneficial mutation
rate U_b is large, many new beneficial alleles can be segregating in the
population simultaneously. We calculate the rate of adaptation, v, in several
models of such sexual populations and show that v is linear in NU_b only in
sufficiently small populations. In large populations, v increases much more
slowly as log NU_b. The prefactor of this logarithm, however, increases as the
square of the recombination rate. This acceleration of adaptation by
recombination implies a strong evolutionary advantage of sex.


  The formal version of our work has been published in BMC Bioinformatics and
can be found here: http://www.biomedcentral.com/1471-2105/13/S6/S1 Motivation:
To tackle the problem of huge memory usage associated with de Bruijn
graph-based algorithms, upon which some of the most widely used de novo genome
assemblers have been built, we released SparseAssembler1. SparseAssembler1 can
save as much as 90% memory consumption in comparison with the state-of-art
assemblers, but it requires rounds of denoising to accurately assemble genomes.
In this paper, we introduce a new general model for genome assembly that uses
only sparse k-mers. The new model replaces the idea of the de Bruijn graph from
the beginning, and achieves similar memory efficiency and much better
robustness compared with our previous SparseAssembler1. Results: We demonstrate
that the decomposition of reads of all overlapping k-mers, which is used in
existing de Bruijn graph genome assemblers, is overly cautious. We introduce a
sparse k-mer graph structure for saving sparse k-mers, which greatly reduces
memory space requirements necessary for de novo genome assembly. In contrast
with the de Bruijn graph approach, we devise a simple but powerful strategy,
i.e., finding links between the k-mers in the genome and traversing following
the links, which can be done by saving only a few k-mers. To implement the
strategy, we need to only select some k-mers that may not even be overlapping
ones, and build the links between these k-mers indicated by the reads. We can
traverse through this sparse k-mer graph to build the contigs, and ultimately
complete the genome assembly. Since the new sparse k-mers graph shares almost
all advantages of de Bruijn graph, we are able to adapt a Dijkstra-like
breadth-first search algorithm to circumvent sequencing errors and resolve
polymorphisms.


  Research in quantitative evolutionary genomics and systems biology led to the
discovery of several universal regularities connecting genomic and molecular
phenomic variables. These universals include the log-normal distribution of the
evolutionary rates of orthologous genes; the power law-like distributions of
paralogous family size and node degree in various biological networks; the
negative correlation between a gene's sequence evolution rate and expression
level; and differential scaling of functional classes of genes with genome
size. The universals of genome evolution can be accounted for by simple
mathematical models similar to those used in statistical physics, such as the
birth-death-innovation model. These models do not explicitly incorporate
selection, therefore the observed universal regularities do not appear to be
shaped by selection but rather are emergent properties of gene ensembles.
Although a complete physical theory of evolutionary biology is inconceivable,
the universals of genome evolution might qualify as 'laws of evolutionary
genomics' in the same sense 'law' is understood in modern physics.


  Experimental studies on enzyme evolution show that only a small fraction of
all possible mutation trajectories are accessible to evolution. However, these
experiments deal with individual enzymes and explore a tiny part of the fitness
landscape. We report an exhaustive analysis of fitness landscapes constructed
with an off-lattice model of protein folding where fitness is equated with
robustness to misfolding. This model mimics the essential features of the
interactions between amino acids, is consistent with the key paradigms of
protein folding and reproduces the universal distribution of evolutionary rates
among orthologous proteins. We introduce mean path divergence as a quantitative
measure of the degree to which the starting and ending points determine the
path of evolution in fitness landscapes. Global measures of landscape roughness
are good predictors of path divergence in all studied landscapes: the mean path
divergence is greater in smooth landscapes than in rough ones. The
model-derived and experimental landscapes are significantly smoother than
random landscapes and resemble additive landscapes perturbed with moderate
amounts of noise; thus, these landscapes are substantially robust to mutation.
The model landscapes show a deficit of suboptimal peaks even compared with
noisy additive landscapes with similar overall roughness. We suggest that
smoothness and the substantial deficit of peaks in the fitness landscapes of
protein evolution are fundamental consequences of the physics of protein
folding.


  Here we show that, given a set of clusters C on a set of taxa X, where |X|=n,
it is possible to determine in time f(k).poly(n) whether there exists a
level-<= k network (i.e. a network where each biconnected component has
reticulation number at most k) that represents all the clusters in C in the
softwired sense, and if so to construct such a network. This extends a
polynomial time result from "On the elusiveness of clusters" by Kelk,
Scornavacca and Van Iersel(2011). By generalizing the concept of "level-k
generator" to general networks, we then extend this fixed parameter
tractability result to the problem where k refers not to the level but to the
reticulation number of the whole network.


  Polarons, introduced by Davydov to explain energy transport in
$\alpha$-helices, correspond to electrons localised on a few lattice sites
because of their interaction with phonons. While the static polaron field
configurations have been extensively studied, their displacement is more
difficult to explain. In this paper we show that, when the next to nearest
neighbour interactions are included, for physical values of the parameters,
polarons can spontaneously move, at T=0, on bent chains that exhibit a positive
gradient in their curvature. At room temperature polarons perform a random walk
but a curvature gradient can induce a non-zero average speed similar to the one
observed at zero temperature. We also show that at zero temperature a polaron
bounces on sharply kinked junctions. We interpret these results in light of the
energy transport by transmembrane proteins.


  Positive feedback and cooperativity in the regulation of gene expression are
generally considered to be necessary for obtaining bistable expression states.
Recently, a novel mechanism of bistability termed emergent bistability has been
proposed which involves only positive feedback and no cooperativity in the
regulation. An additional positive feedback loop is effectively generated due
to the inhibition of cellular growth by the synthesized proteins. The
mechanism, demonstrated for a synthetic circuit, may be prevalent in natural
systems also as some recent experimental results appear to suggest. In this
paper, we study the effects of additive and multiplicative noise on the
dynamics governing emergent bistability. The calculational scheme employed is
based on the Langevin and Fokker-Planck formalisms. The steady state
probability distributions of protein levels and the mean first passage times
are computed for different noise strengths and system parameters. In the region
of bistability, the bimodal probability distribution is shown to be a linear
combination of a lognormal and a Gaussian distribution. The variances of the
individual distributions and the relative weights of the distributions are
further calculated for varying noise strengths and system parameters. The
experimental relevance of the model results is also pointed out.


  We consider an infinite-dimensional system of stochastic differential
equations describing the evolution of type frequencies in a large population.
The type of an individual is the number of deleterious mutations it carries,
where fitness of individuals carrying k mutations is decreased by \alpha k for
some \alpha>0. Along the individual lines of descent, new mutations accumulate
at rate \lambda per generation, and each of these mutations has a probability
\gamma per generation to disappear. While the case \gamma=0 is known as (the
Fleming-Viot version of) Muller's ratchet, the case \gamma>0 is associated with
compensatory mutations in the biological literature. We show that the system
has a unique weak solution. In the absence of random fluctuations in type
frequencies (i.e., for the so-called infinite population limit) we obtain the
solution in a closed form by analyzing a probabilistic particle system and show
that for \gamma>0, the unique equilibrium state is the Poisson distribution
with parameter \lambda/(\gamma+\alpha).


  We study a mechanism for reliable switching in biomolecular
signal-transduction cascades. Steady bistable states are created by system-size
cooperative effects in populations of proteins, in spite of the fact that the
phosphorylation-state transitions of any molecule, by means of which the switch
is implemented, are highly stochastic. The emergence of switching is a
nonequilibrium phase transition in an energetically driven, dissipative system
described by a master equation. We use operator and functional integral methods
from reaction-diffusion theory to solve for the phase structure, noise
spectrum, and escape trajectories and first-passage times of a class of minimal
models of switches, showing how all critical properties for switch behavior can
be computed within a unified framework.


  Genetic algorithms are considered as an original way to solve problems,
probably because of their generality and of their "blind" nature. But GAs are
also unusual since the features of many implementations (among all that could
be thought of) are principally led by the biological metaphor, while efficiency
measurements intervene only afterwards. We propose here to examine the
relevance of these biomimetic aspects, by pointing out some fundamental
similarities and divergences between GAs and the genome of living beings shaped
by natural selection. One of the main differences comes from the fact that GAs
rely principally on the so-called implicit parallelism, while giving to the
mutation/selection mechanism the second role. Such differences could suggest
new ways of employing GAs on complex problems, using complex codings and
starting from nearly homogeneous populations.


  The dynamics of noise-resilient Boolean networks with majority functions and
diverse topologies is investigated. A wide class of possible topological
configurations is parametrized as a stochastic blockmodel. For this class of
networks, the dynamics always undergoes a phase transition from a non-ergodic
regime, where a memory of its past states is preserved, to an ergodic regime,
where no such memory exists and every microstate is equally probable. Both the
average error on the network, as well as the critical value of noise where the
transition occurs are investigated analytically, and compared to numerical
simulations. The results for "partially dense" networks, comprised of
relatively few, but dynamically important nodes, which have a number of inputs
which greatly exceeds the average for the entire network, give very general
upper bounds on the maximum resilience against noise attainable on globally
sparse systems.


  We investigate the evolution of Boolean networks subject to a selective
pressure which favors robustness against noise, as a model of evolved genetic
regulatory systems. By mapping the evolutionary process into a statistical
ensemble and minimizing its associated free energy, we find the structural
properties which emerge as the selective pressure is increased and identify a
phase transition from a random topology to a "segregated core" structure, where
a smaller and more densely connected subset of the nodes is responsible for
most of the regulation in the network. This segregated structure is very
similar qualitatively to what is found in gene regulatory networks, where only
a much smaller subset of genes --- those responsible for transcription factors
--- is responsible for global regulation. We obtain the full phase diagram of
the evolutionary process as a function of selective pressure and the average
number of inputs per node. We compare the theoretical predictions with Monte
Carlo simulations of evolved networks and with empirical data for Saccharomyces
cerevisiae and Escherichia coli.


  We use computer simulations to investigate self-assembly in a system of model
chaperonin proteins, and in an Ising lattice gas. We discuss the mechanisms
responsible for rapid and efficient assembly in these systems, and we use
measurements of dynamical activity and assembly progress to compare their
propensities for kinetic trapping. We use the analytic solution of a simple
minimal model to illustrate the key features associated with such trapping,
paying particular attention to the number of ways that particles can misbind.
We discuss the relevance of our results for the design and control of
self-assembly in general.


  The spectral bound, s(a A + b V), of a combination of a resolvent positive
linear operator A and an operator of multiplication V, was shown by Kato to be
convex in b \in R. This is shown here, through an elementary lemma, to imply
that s(a A + b V) is also convex in a > 0, and notably, \partial s(a A + b V) /
\partial a <= s(A) when it exists. Diffusions typically have s(A) <= 0, so that
for diffusions with spatially heterogeneous growth or decay rates, greater
mixing reduces growth. Models of the evolution of dispersal in particular have
found this result when A is a Laplacian or second-order elliptic operator, or a
nonlocal diffusion operator, implying selection for reduced dispersal. These
cases are shown here to be part of a single, broadly general, `reduction'
phenomenon.


  The phenomenon of spatial clustering induced by death and reproduction in a
population of anomalously diffusing individuals is studied analytically. The
possibility of social behaviors affecting the migration strategies has been
taken into exam, in the case anomalous diffusion is produced by means of a
continuous time random walk (CTRW). In the case of independently diffusing
individuals, the dynamics appears to coincide with that of (dying and
reproducing) Brownian walkers. In the strongly social case, the dynamics
coincides with that of non-migrating individuals. In both limits, the growth
rate of the fluctuations becomes independent of the Hurst exponent of the CTRW.
The social behaviors that arise when transport in a population is induced by a
spatial distribution of random traps, have been analyzed.


  Here we present a new fixed parameter tractable algorithm to compute the
hybridization number r of two rooted binary phylogenetic trees on taxon set X
in time (6r)^r.poly(n), where n=|X|. The novelty of this approach is that it
avoids the use of Maximum Acyclic Agreement Forests (MAAFs) and instead
exploits the equivalence of the problem with a related problem from the
softwired clusters literature. This offers an alternative perspective on the
underlying combinatorial structure of the hybridization number problem.


  The extent to which ART (anti-retroviral therapy) reduces HIV transmission
has received attention in recent years. Using data on the relationship between
transmission and viral load we show that transmission saturates at high viral
loads. We fit a power-law model and an exponential converging to an asymptote.
Using data on the viral load in HIV-positive people we show that ART is likely
to reduce transmission by 91.6% (81.7%-96.2%) under the first and 99.5%
(98.5%-99.8%) under the second model.
  The role of the acute phase in HIV transmission is still debated. High levels
of transmission during the acute phase have been used to argue that failure to
identify people in the acute phase of HIV may compromise the impact of
treatment on preventing new infections and that having concurrent sexual
partners during the acute phase is an important driver of the epidemic. We show
that the acute phase probably accounts for less than 1% of overall
transmission. We also show that even if a significant proportion of infections
are transmitted during the acute phase, this will not compromise the impact of
treatment on population levels of transmission given the constraint implied by
the doubling time of the epidemic.
  This analysis leads to other relevant conclusions. First, it is likely that
discordant-couple studies significantly underestimate the risk of infection.
Second, attention should be paid to the variability in set point viral load
which determines both the infectiousness of HIV-positive people and the
variability in the susceptibility of HIV-negative people. Third, if ART drugs
are in short supply those with the highest viral load should be given priority,
others things including age, gender and opportunistic infections being equal,
but to reduce transmission ART should be offered to all those with a viral load
above about 10k/mm.3


  Microtubules have been in biophysical focus for several decades. Yet the
confusing and mutually contradicting results regarding their elasticity and
fluctuations have shed some doubts on their present understanding. In this
paper we expose the empirical evidence for the existence of discrete
GDP-tubulin fluctuations between a curved and a straight configuration at room
temperature as well as for conformational tubulin cooperativity. Guided by a
number of experimental findings, we build the case for a novel microtubule
model, with the principal result that microtubules can spontaneously form
micron size cooperative helical states with unique elastic and dynamic
features. The polymorphic dynamics of the microtubule lattice resulting from
the tubulin bistability quantitatively explains several experimental puzzles
including anomalous scaling of dynamic fluctuations of grafted microtubules,
their apparent length-stiffness relation and their remarkably curved-helical
appearance in general. We point out that tubulin dimers's multistability and
its cooperative switching could participate in important cellular processes,
and could in particular lead to efficient mechanochemical signalling along
single microtubules.


  We perform a full similarity analysis of an idealized ecosystem using
Buckingham's $\Pi$ theorem to obtain dimensionless similarity parameters given
that some (non- unique) method exists that can differentiate different
functional groups of individuals within an ecosystem. We then obtain the
relationship between the similarity parameters under the assumptions of (i)
that the ecosystem is in a dynamically balanced steady state and (ii) that
these functional groups are connected to each other by the flow of resource.
The expression that we obtain relates the level of complexity that the
ecosystem can support to intrinsic macroscopic variables such as density,
diversity and characteristic length scales for foraging or dispersal, and
extrinsic macroscopic variables such as habitat size and the rate of supply of
resource. This expression relates these macroscopic variables to each other,
generating commonly observed macroecological patterns; these broad trends
simply reflect the similarity property of ecosystems. We thus find that details
of the ecosystem function are not required to obtain these broad
macroecological patterns this may explain why they are ubiquitous. Departures
from our relationship may indicate that the ecosystem is in a state of rapid
change, i.e., abundance or diversity explosion or collapse. Our result provides
normalised variables that can be used to isolate the trend in one ecosystem
variable from another, providing a new method for isolating macroecological
patterns in data. A dimensionless control parameter for ecosystem complexity
emerges from our analysis and this will be a control parameter in dynamical
models for ecosystems based on energy flow and conservation and will order the
emergent behaviour of these models.


  Changes in hydration are central to the phenomenon of biomolecular
recognition, but it has been difficult to properly frame and answer questions
about their precise thermodynamic role. We address this problem by introducing
Grid Inhomogeneous Solvation Theory (GIST), which discretizes the equations of
Inhomogeneous Solvation Theory on a 3D grid in a volume of interest. Here, the
solvent volume is divided into small grid boxes and localized thermodynamic
entropies, energies and free energies are defined for each grid box.
Thermodynamic solvation quantities are defined in such a manner that summing
the quantities over all the grid boxes yields the desired total quantity for
the system. This approach smoothly accounts for the thermodynamics of not only
highly occupied water sites but also partly occupied and water depleted regions
of the solvent, without the need for ad hoc terms drawn from other theories.
The GIST method has the further advantage of allowing a rigorous end-states
analysis that, for example in the problem of molecular recognition, can account
for not only the thermodynamics of displacing water from the surface but also
for the thermodynamics of solvent reorganization around the bound complex. As a
preliminary application, we present GIST calculations at the 1-body level for
the host cucurbit[7]uril, a low molecular weight receptor molecule which
represents a tractable model for biomolecular recognition. One of the most
striking results is the observation of a toroidal region of water density, at
the center of the host's nonpolar cavity, which is significantly disfavored
entropically, and hence may contribute to the ability of this small receptor to
bind guest molecules with unusually high affinities.


  Much of the work on modeling the spread of viral infections utilized partial
differential equa- tions. Traveling-wave solutions to these PDEs are typically
concentrated on velocities and their dependence on the various parameters. Most
of the investigations into the dynamical interaction of virus and defective
interfering particles (DIP), which are incomplete forms of the virus that
replicate through co-infection, have followed the same lines. In this work we
present an agent based model of viral infection with consideration of DIP and
the negative feedback loop introduced by interferon production as part of the
host innate immune response. The model is based high resolution microscopic
images of plaques of dead cells we took from mammalian cells infected with
Sendai virus with low and high DIP. In order to investigate the effects of the
discrete stochastic microscopic mechanisms, which are responsible for virus
spreading, have on the macroscopic growth of viral plaques, we generate an
agent-based model of viral infection. The two main aims of this work are to (i)
investigate the effects of discrete microscopic randomness on the macroscopic
growth of viral plaques; and (ii) examine the dynamic interactions between the
full length virus, DIP and interferon, and interpret what may be the function
of DIP. We find that we can explain the qualitative differences between our
stochastic model and deterministic models in terms of the fractal geometry of
the resulting plaques, and that DIP have a delaying effect while the
interaction between interferon and DIP has a slowing effect on the growth of
viral plaques, potentially contributing to viral latency.


  It is well known that during homology recognition and strand exchange the
double stranded DNA (dsDNA) in DNA/RecA filaments is highly extended, but the
functional role of the extension has been unclear. We present an analytical
model that calculates the distribution of tension in the extended dsDNA during
strand exchange. The model suggests that the binding of additional dsDNA base
pairs to the DNA/RecA filament alters the tension in dsDNA that was already
bound to the filament, resulting in a non-linear increase in the mechanical
energy as a function of the number of bound base pairs. This collective
mechanical response may promote homology stringency and underlie unexplained
experimental results.


  The intensification and increased frequency of weather extremes is emerging
as one of the most important aspects of climate change. We use Monte Carlo
simulation to understand and predict the consequences of variations in trends
(i.e., directional change) and stochasticity (i.e., increase in variance) of
climate variables and consequent selection pressure by using simple models of
population dynamics. Higher variance of climate variables increases the
probability of weather extremes and consequent catastrophic disturbances.
Parameters of the model are selection pressure, mutation, directional and
stochastic variation of the environment. We follow the population dynamics and
the distribution of a trait that describes the adaptation of the individual to
the optimum phenotype defined by the environmental conditions. The survival
chances of a population depend quite strongly on the selection pressure and
decrease with increasing variance of the climate variable. In general, the
system is able to track the directional component of the optimum phenotype.
Intermediate levels of mutation generally increase the probability of tracking
the changing optimum and thus decrease the risk of extinction of a population.
With high mutation, the higher probability of maladaptation decreases the
survival chances of the populations, even with high variability of the optimum
phenotype.