  Understanding under what conditions populations, whether they be plants,
animals, or viral particles, persist is an issue of theoretical and practical
importance in population biology. Both biotic interactions and environmental
fluctuations are key factors that can facilitate or disrupt persistence. One
approach to examining the interplay between these deterministic and stochastic
forces is the construction and analysis of stochastic difference equations
$X_{t+1}=F(X_t,\xi_{t+1})$ where $X_t \in \R^k$ represents the state of the
populations and $\xi_1,\xi_2,...$ is a sequence of random variables
representing environmental stochasticity. In the analysis of these stochastic
models, many theoretical population biologists are interested in whether the
models are bounded and persistent. Here, boundedness asserts that
asymptotically $X_t$ tends to remain in compact sets. In contrast, persistence
requires that $X_t$ tends to be "repelled" by some "extinction set" $S_0\subset
\R^k$. Here, results on both of these proprieties are reviewed for single
species, multiple species, and structured population models. The results are
illustrated with applications to stochastic versions of the Hassell and Ricker
single species models, Ricker, Beverton-Holt, lottery models of competition,
and lottery models of rock-paper-scissor games. A variety of conjectures and
suggestions for future research are presented.


  The distances over which biological molecules and their complexes can
function range from a few nanometres, in the case of folded structures, to
millimetres, for example during chromosome organization. Describing phenomena
that cover such diverse length, and also time scales, requires models that
capture the underlying physics for the particular length scale of interest.
Theoretical ideas, in particular, concepts from polymer physics, have guided
the development of coarse-grained models to study folding of DNA, RNA, and
proteins. More recently, such models and their variants have been applied to
the functions of biological nanomachines. Simulations using coarse-grained
models are now poised to address a wide range of problems in biology.


  An overwhelming majority of humans are right-handed. Numerous explanations
for individual handedness have been proposed, but this population-level
handedness remains puzzling. Here we use a minimal mathematical model to
explain this population-level hand preference as an evolved balance between
cooperative and competitive pressures in human evolutionary history. We use
selection of elite athletes as a test-bed for our evolutionary model and
account for the surprising distribution of handedness in many professional
sports. Our model predicts strong lateralization in social species with limited
combative interaction, and elucidates the rarity of compelling evidence for
"pawedness" in the animal world.


  Key insights in molecular biology, such as enzyme kinetics, protein allostery
and gene regulation emerged from quantitative analysis based on time-scale
separation, allowing internal complexity to be eliminated and resulting in the
well-known formulas of Michaelis-Menten, Monod-Wyman-Changeux and
Ackers-Johnson-Shea. In systems biology, steady-state analysis has yielded
eliminations that reveal emergent properties of multi-component networks. Here
we show that these analyses of nonlinear biochemical systems are consequences
of the same linear framework, consisting of a labelled, directed graph on which
a Laplacian dynamics is defined, whose steady states can be algorithmically
calculated. Analyses previously considered distinct are revealed as identical,
while new methods of analysis become feasible.


  Complete genome sequences contain valuable information about natural
selection, but extracting this information for short, widely scattered
noncoding elements remains a challenging problem. Here we introduce a new
computational method for addressing this problem called Inference of Natural
Selection from Interspersed Genomically coHerent elemenTs (INSIGHT). INSIGHT
uses a generative probabilistic model to contrast patterns of polymorphism and
divergence in the elements of interest with those in flanking neutral sites,
pooling weak information from many short elements in a manner that accounts for
variation among loci in mutation rates and genealogical backgrounds. The method
is able to disentangle the contributions of weak negative, strong negative, and
positive selection based on their distinct effects on patterns of polymorphism
and divergence. Information about divergence is obtained from multiple outgroup
genomes using a full phylogenetic model. The model is efficiently fitted to
genome-wide data by decomposing the maximum likelihood estimation procedure
into three straightforward stages. The key selection-related parameters are
estimated by expectation maximization. Using simulations, we show that INSIGHT
can accurately estimate several parameters of interest even in complex
demographic scenarios. We apply our methods to noncoding RNAs, promoter
regions, and transcription factor binding sites in the human genome, and find
clear evidence of natural selection. We also present a detailed analysis of
particular nucleotide positions within GATA2 binding sites and primary
micro-RNA transcripts.


  We introduce and numerically solve a mathematical model of the experimentally
established mechanisms responsible for the symmetry breaking transition
observed in the chiral crystallization experiments reported by I. Weissbuch, L.
Addadi, L. Leiserowitz and M. Lahav, J. Am. Chem. Soc. 110 (1988), 561-567. The
mathematical model is based on five basic processes: (1) The formation of
achiral glycine clusters in solution, (2) The nucleation of oriented glycine
crystals at the air/water interface in the presence of hydrophobic amino acids,
(3) A kinetic orienting effect which inhibits crystal growth, (4) The
enantioselective occlusion of the amino acids from solution, and (5) The growth
of oriented host glycine crystals at the interface. We translate these
processes into differential rate equations. We first study the model with the
orienting process (2) without (3) and then combine both allowing us to make
detailed comparisons of both orienting effects which actually act in unison in
the experiment. Numerical results indicate that the model can yield a high
percentage orientation of the mixed crystals at the interface and the
consequent resolution of the initially racemic mixture of amino acids in
solution. The model thus leads to separation of enantiomeric territories, the
generation and amplification of optical activity by enantioselective occlusion
of chiral additives through chiral surfaces of glycine crystals.


  In this work two archaea microorganisms (Haloferax volcanii and Natrialba
magadii) used as biocatalyst at a microbial fuel cell (MFC) anode were
evaluated. Both archaea are able to grow at high salt concentrations. By
increasing the media conductivity, the internal resistance was diminished,
improving the MFCs performance. Without any added redox mediator, maximum power
(Pmax) and current at Pmax were 11.87 / 4.57 / 0.12 {\mu}W cm-2 and 49.67 /
22.03 / 0.59 {\mu}A cm-2 for H. volcanii, N. magadii and E. coli, respectively.
When neutral red was used as redox mediator, Pmax was 50.98 and 5.39 {\mu}W
cm-2 for H. volcanii and N. magadii respectively. In this paper an archaea MFC
is described and compared with other MFC systems; the high salt concentration
assayed here, comparable with that used in Pt-catalyzed alkaline hydrogen fuel
cells will open new options when MFC scaling-up is the objective, necessary for
practical applications.


  The wormlike chain (WLC) model of DNA bending accurately reproduces
single-molecule force-extension profiles of long (kilobase) chains. These
bending statistics over large scales do not, however, establish a unique
microscopic model for elasticity at the 1-10 bp scale, which holds particular
interest in biological contexts. Here we examine a class of microscopic models
which allow for disruption of base pairing (i.e., a `melt' or `kink',
generically an `excitation') and consequently enhanced local flexibility. We
first analyze the effect on the excitation free energy of integrating out the
spatial degrees of freedom in a wormlike chain. Based on this analysis, we
present a formulation of these models that ensures consistency with the
well-established thermodynamics of melting in long chains. Using a new method
to calculate cyclization statistics of short chains from enhanced-sampling
Monte Carlo simulations, we compute J-factors of a meltable wormlike chain
(MWLC) over a broad range of chain lengths, including very short molecules (30
bp) that have not yet been explored experimentally. For chains longer than
about 120 bp, including most molecules studied to date in the laboratory, we
find that melting excitations have little impact on cyclization kinetics.
Strong signatures of melting, which might be resolved within typical
experimental scatter, emerge only for shorter chains.


  We suggest that the thermodynamic stability parameters (nearest neighbor
stacking and hydrogen bonding free energies) of double-stranded DNA molecules
can be inferred reliably from time series of the size fluctuations (breathing)
of local denaturation zones (bubbles). On the basis of the reconstructed bubble
size distribution, this is achieved through stochastic optimization of the free
energies in terms of Simulated Annealing. In particular, it is shown that even
noisy time series allow the identification of the stability parameters at
remarkable accuracy. This method will be useful to obtain the DNA stacking and
hydrogen bonding free energies from single bubble breathing assays rather than
equilibrium data.


  We analyze a class of partial differential equations that arise as "backwards
Kolmogorov operators" in infinite population limits of the Wright-Fisher models
in population genetics and in mathematical finance. These are degenerate
elliptic operators defined on manifolds with corners. The classical example is
the Kimura diffusion operator, which acts on functions defined on the simplex
in R^n. We introduce anisotropic Holder spaces, and prove existence, uniqueness
and regularity results for the heat and resolvent equations defined by this
class of operators. This suffices to prove that the C^0-graph closure generates
a strongly continuous semigroup, and that the associated Martingale problem has
a unique solution. We give a detailed description of the nullspace of the
forward Kolmogorov operator.


  Further to the proposal and application of a stochastic methodology and the
resulting first exit time distribution function to life table data we introduce
a theoretical framework for the estimation of the maximum deterioration age and
to explore on how "vitality," according to Halley and Strehler and Mildvan,
changes during the human lifetime. The mortality deceleration or mortality
leveling-off is also explored. The effect of the deterioration over time is
estimated as the expectation that an individual will survive from the
deterioration caused in his organism by the deterioration mechanism. A method
is proposed and the appropriate software was developed for the estimation of
life expectancy. Several applications follow. The method was applied to the
Halley life table data of Breslau. Extrapolations are done showing a gradual
improvement of vitality mechanisms during last centuries.


  Many cytoskeletal biopolymers are "active," consuming energy in large
quantities. In this Letter, we identify a fundamental difference between active
polymers and passive, equilibrium polymers: for equal mean lengths, active
polymers can reorganize faster than equilibrium polymers. We show that
equilibrium polymers are intrinsically limited to linear scaling between mean
lifetime and mean length, MFPT ~ <L>, by analogy to 1-d Potts models. By
contrast, we present a simple active-polymer model that improves upon this
scaling, such that MFPT ~ <L>^{1/2}. Since to be biologically useful,
structural biopolymers must typically be many monomers long, yet respond
dynamically to the needs of the cell, the difference in reorganization kinetics
may help to justify active polymers' greater energy cost. PACS numbers:
87.10.Ed, 87.16.ad, 87.16.Ln


  We have established an RNA Mapping Database (RMDB) to enable a new generation
of structural, thermodynamic, and kinetic studies from quantitative
single-nucleotide-resolution RNA structure mapping (freely available at
http://rmdb.stanford.edu). Chemical and enzymatic mapping is a rapid, robust,
and widespread approach to RNA characterization. Since its recent coupling with
high-throughput sequencing techniques, accelerated software pipelines, and
large-scale mutagenesis, the volume of mapping data has greatly increased, and
there is a critical need for a database to enable sharing, visualization, and
meta-analyses of these data. Through its on-line front-end, the RMDB allows
users to explore single-nucleotide-resolution chemical accessibility data in
heat-map, bar-graph, and colored secondary structure graphics; to leverage
these data to generate secondary structure hypotheses; and to download the data
in standardized and computer-friendly files, including the RDAT and
community-consensus SNRNASM formats. At the time of writing, the database
houses 38 entries, describing 2659 RNA sequences and comprising 355,084 data
points, and is growing rapidly.


  We show that in a particular model of catalytic reaction systems, known as
the binary polymer model, there is a mathematical invariance between two
versions of the model: (1) random catalysis and (2) template-based catalysis.
In particular, we derive an analytical calculation that allows us to accurately
predict the (observed) required level of catalysis in one version of the model
from that in the other version, for a given probability of having
self-sustaining autocatalytic sets exist in instances of both model versions.
This provides a tractable connection between two models that have been
investigated in theoretical origin-of-life studies.


  Three-dimensional RNA models fitted into crystallographic density maps
exhibit pervasive conformational ambiguities, geometric errors and steric
clashes. To address these problems, we present enumerative real-space
refinement assisted by electron density under Rosetta (ERRASER), coupled to
Python-based hierarchical environment for integrated 'xtallography' (PHENIX)
diffraction-based refinement. On 24 data sets, ERRASER automatically corrects
the majority of MolProbity-assessed errors, improves the average Rfree factor,
resolves functionally important discrepancies in noncanonical structure and
refines low-resolution models to better match higher-resolution models.


  We have studied the evolution of strategies in spatial public goods games
where both individual (peer) and institutional (pool) punishments are present
beside unconditional defector and cooperator strategies. The evolution of
strategy distribution is governed by imitation based on random sequential
comparison of neighbors' payoff for a fixed level of noise. Using numerical
simulations we have evaluated the strategy frequencies and phase diagrams when
varying the synergy factor, punishment cost, and fine. Our attention is focused
on two extreme cases describing all the relevant behaviors in such a complex
system. According to our numerical data peer punishers prevail and control the
system behavior in a large segments of parameters while pool punishers can only
survive in the limit of weak peer punishment when a rich variety of solutions
is observed. Paradoxically, the two types of punishment may extinguish each
other's impact resulting in the triumph of defectors. The technical
difficulties and suggested methods are briefly discussed.


  Previous work has shown that species interacting in an ecosystem and actors
transacting in an economic context may have notable similarities in behavior.
However, the specific mechanism that may underlie similarities in nature and
human systems has not been analyzed. Building on stochastic food-web models, we
propose a parsimonious bipartite-cooperation model that reproduces the key
features of mutualistic networks - degree distribution, nestedness and
modularity -- for both ecological networks and socio-economic networks. Our
analysis uses two diverse networks. Mutually-beneficial interactions between
plants and their pollinators, and cooperative economic exchanges between
designers and their contractors. We find that these mutualistic networks share
a key hierarchical ordering of their members, along with an exponential
constraint in the number and type of partners they can cooperate with. We use
our model to show that slight changes in the interaction constraints can
produce either extremely nested or random structures, revealing that these
constraints play a key role in the evolution of mutualistic networks. This
could also encourage a new systematic approach to study the functional and
structural properties of networks. The surprising correspondence across
mutualistic networks suggests their broadly representativeness and their
potential role in the productive organization of exchange systems, both
ecological and social.


  Successful animal systems often manage risk through synchronous behavior that
spontaneously arises without leadership. In critical human systems facing risk,
such as financial markets or military operations, our understanding of the
benefits associated to synchronicity is nascent but promising. Building on
previous work illuminating commonalities between ecological and human systems,
we compare the activity patterns of individual financial traders with the
simultaneous activity of other traders---an individual and spontaneous
characteristic we call synchronous trading. Additionally, we examine the
association of synchronous trading with individual performance and
communication patterns. Analyzing empirical data on day traders'
second-to-second trading and instant messaging, we find that the higher the
traders' synchronous trading, the less likely they lose money at the end of the
day. We also find that the daily instant messaging patterns of traders are
closely associated with their level of synchronous trading. This suggests that
synchronicity and vanguard technology may help cope with risky decisions in
complex systems and furnish new prospects for achieving collective and
individual goals.


  Phylogenetic networks are a generalization of phylogenetic trees that are
used in biology to represent reticulate or non-treelike evolution. Recently,
several algorithms have been developed which aim to construct phylogenetic
networks from biological data using {\em triplets}, i.e. binary phylogenetic
trees on 3-element subsets of a given set of species. However, a fundamental
problem with this approach is that the triplets displayed by a phylogenetic
network do not necessary uniquely determine or {\em encode} the network. Here
we propose an alternative approach to encoding and constructing phylogenetic
networks, which uses phylogenetic networks on 3-element subsets of a set, or
{\em trinets}, rather than triplets. More specifically, we show that for a
special, well-studied type of phylogenetic network called a 1-nested network,
the trinets displayed by a 1-nested network always encode the network. We also
present an efficient algorithm for deciding whether a {\em dense} set of
trinets (i.e. one that contains a trinet on every 3-element subset of a set)
can be displayed by a 1-nested network or not and, if so, constructs that
network. In addition, we discuss some potential new directions that this new
approach opens up for constructing and comparing phylogenetic networks.


  The tertiary structures of functional RNA molecules remain difficult to
decipher. A new generation of automated RNA structure prediction methods may
help address these challenges but have not yet been experimentally validated.
Here we apply four prediction tools to a remarkable class of double glycine
riboswitches that exhibit ligand-binding cooperativity. A novel method
(BPPalign), RMdetect, JAR3D, and Rosetta 3D modeling give consistent
predictions for a new stem P0 and kink-turn motif. These elements structure the
linker between the RNAs' double aptamers. Chemical mapping on the F. nucleatum
riboswitch with SHAPE, DMS, and CMCT probing, mutate-and-map studies, and
mutation/rescue experiments all provide strong evidence for the structured
linker. Under solution conditions that separate two glycine binding
transitions, disrupting this helix-junction-helix structure gives 120-fold and
6- to 30-fold poorer association constants for the two transitions,
corresponding to an overall energetic impact of 4.3 \pm 0.5 kcal/mol. Prior
biochemical and crystallography studies from several labs did not include this
critical element due to over-truncation of the RNA. We argue that several
further undiscovered elements are likely to exist in the flanking regions of
this and other RNA switches, and automated prediction tools can now play a
powerful role in their detection and dissection.


  The selection of the most suitable evolutionary model to analyze the given
molecular data is usually left to biologist's choice. In his famous book, J
Felsenstein suggested that certain linear equations satisfied by the expected
probabilities of patterns observed at the leaves of a phylogenetic tree could
be used for model selection. It remained open the question regarding whether
these equations were enough for characterizing the evolutionary model.
  Here we prove that, for equivariant models of evolution, the space of
distributions satisfying these linear equations coincides with the space of
distributions arising from mixtures of trees on a set of taxa. In other words,
we prove that an alignment is produced from a mixture of phylogenetic trees
under an equivariant evolutionary model if and only if its distribution of
column patterns satisfies the linear equations mentioned above. Moreover, for
each equivariant model and for any number of taxa, we provide a set of linearly
independent equations defining this space of phylogenetic mixtures. This is a
powerful tool that has already been successfully used in model selection. We
also use the results obtained to study identifiability issues for phylogenetic
mixtures.


  The Indus Valley Civilization (IVC) was one of the first great civilizations
in prehistory. This bronze age civilization flourished from the end of the
fourth millennium BC. It disintegrated during the second millennium BC; despite
much research effort, this decline is not well understood. Less research has
been devoted to the emergence of the IVC, which shows continuous cultural
precursors since at least the seventh millennium BC. To understand the decline,
we believe it is necessary to investigate the rise of the IVC, i.e., the
establishment of agriculture and livestock, dense populations and technological
developments 7000--3000 BC. Although much archaeological information is
available, our capability to investigate the system is hindered by poorly
resolved chronology, and by a lack of field work in the intermediate areas
between the Indus valley and Mesopotamia. We thus employ a complementary
numerical simulation to develop a consistent picture of technology,
agropastoralism and population developments in the IVC domain. Results from
this Global Land Use and technological Evolution Simulator show that there is
(1) fair agreement between the simulated timing of the agricultural transition
and radiocarbon dates from early agricultural sites, but the transition is
simulated first in India then Pakistan; (2) an independent agropastoralism
developing on the Indian subcontinent; and (3) a positive relationship between
archeological artifact richness and simulated population density which remains
to be quantified.


  A nearest-neighbor-interchange (NNI) walk is a sequence of unrooted
phylogenetic trees, T_0, T_1, T_2,... where each consecutive pair of trees
differ by a single NNI move. We give tight bounds on the length of the shortest
NNI-walks that visit all trees in an subtree-prune-and-regraft (SPR)
neighborhood of a given tree. For any unrooted, binary tree, T, on n leaves,
the shortest walk takes {\theta}(n^2) additional steps than the number of trees
in the SPR neighborhood. This answers Bryant's Second Combinatorial Conjecture
from the Phylogenetics Challenges List, the Isaac Newton Institute, 2011, and
the Penny Ante Problem List, 2009.


  We propose a driven optical cavity quantum electrodynamics (QED) set up aimed
at directly probing energy transport dynamics in photosynthetic biomolecules.
We show that detailed information concerning energy transfer paths and
delocalization of exciton states can be inferred (and exciton energies
estimated) from the statistical properties of the emitted photons. This
approach provides us with a novel spectroscopic tool for the interrogation of
biological systems in terms of quantum optical phenomena which have been
usually studied for atomic or solid-state systems, e.g. trapped atoms and
semiconductor quantum dots.


  Ordinary differential equations are widely-used in the field of systems
biology and chemical engineering to model chemical reaction networks. Numerous
techniques have been developed to estimate parameters like rate constants,
initial conditions or steady state concentrations from time-resolved data. In
contrast to this countable set of parameters, the estimation of entire courses
of network components corresponds to an innumerable set of parameters. The
approach presented in this work is able to deal with course estimation for
extrinsic system inputs or intrinsic reactants, both not being constrained by
the reaction network itself. Our method is based on variational calculus which
is carried out analytically to derive an augmented system of differential
equations including the unconstrained components as ordinary state variables.
Finally, conventional parameter estimation is applied to the augmented system
resulting in a combined estimation of courses and parameters. The combined
estimation approach takes the uncertainty in input courses correctly into
account. This leads to precise parameter estimates and correct confidence
intervals. In particular this implies that small motifs of large reaction
networks can be analysed independently of the rest. By the use of variational
methods, elements from control theory and statistics are combined allowing for
future transfer of methods between the two fields.


  Rift Valley fever is a vector-borne zoonotic disease which causes high
morbidity and mortality in livestock. In the event Rift Valley fever virus is
introduced to the United States or other non-endemic areas, understanding the
potential patterns of spread and the areas at risk based on disease vectors and
hosts will be vital for developing mitigation strategies. Presented here is a
general network-based mathematical model of Rift Valley fever. Given a lack of
empirical data on disease vector species and their vector competence, this
discrete time epidemic model uses stochastic parameters following several PERT
distributions to model the dynamic interactions between hosts and likely North
American mosquito vectors in dispersed geographic areas. Spatial effects and
climate factors are also addressed in the model. The model is applied to a
large directed asymmetric network of 3,621 nodes based on actual farms to
examine a hypothetical introduction to some counties of Texas, an important
ranching area in the United States of America (U.S.A.). The nodes of the
networks represent livestock farms, livestock markets, and feedlots, and the
links represent cattle movements and mosquito diffusion between different
nodes. Cattle and mosquito (Aedes and Culex) populations are treated with
different contact networks to assess virus propagation. Rift Valley fever virus
spread is assessed under various initial infection conditions (infected
mosquito eggs, adults or cattle). A surprising trend is fewer initial
infectious organisms result in a longer delay before a larger and more
prolonged outbreak. The delay is likely caused by a lack of herd immunity while
the infections expands geographically before becoming an epidemic involving
many dispersed farms and animals almost simultaneously.


  Introduction: Gap junctions play an important role in the cell proliferation
in mammalian cells as well as carcinogenesis. However, there are controversial
issues about their role in cancer pathogenesis. This study was designed to
evaluate genotoxicity and cytotoxicity of Carbenoxolone (CBX) as a prototype of
inter-cellular gap junction blocker in MCF7 and BT20 human breast cancer cells.
Methods: The MCF7and BT20 human breast cancer cell lines were cultivated, and
treated at designated confluency with different doses of CBX. Cellular
cytotoxicity was examined using standard colorimetric assay associated with
cell viability tests. Gene expression evaluation was carried out using real
time polymerase chain reaction (PCR). Results: MCF7 and BT20 cells were
significantly affected by CBX in a dose dependent manner in cell viability
assays. Despite varying expression of genes, down regulation of pro- and
anti-apoptotic genes was observed in these cells. Conclusion: Based upon this
investigation, it can be concluded that CBX could affect both low and high
proliferative types of breast cancer cell lines and disproportionate down
regulation of both pre- and anti-apoptotic genes may be related to interacting
biomolecules, perhaps via gap junctions.


  The outbreak of enterohemorrhagic Escherichia coli (EHEC) in May 2011 warns
the potential threats of the world vegetables trade network (VTN) in spreading
fatal infectious diseases. The heterogeneous weight distribution and
multi-scale activity of intermediary networks affects the diffusion,
proliferation and extinction of epidemics. Here, we constructed a dual-weighted
VTN with 118 major countries and territories from FAO 2008 statistic data about
global vegetation production and trade, and develop a reaction-diffusion model
to simulate the epidemic behaviors in through VTN. We found an emerged
asymmetric threshold of epidemic on VTN, in which local proliferation within
nodes plays a more critical role than global diffusion in spreading of
EHEC-like diseases, i.e. sufficient local proliferation is the precondition for
global diffusion. We also found that a strong modularity on VTN structure,
which restricts the spreading of EHEC-like diseases; however, within the
communities, the diffusion is quick and easy. There is, moreover, a critical
"epidemic stem", in which a serial of positive feedback loop for amplifying the
proliferation and diffusion pathogens has been identified from entire VTN.
Surprisingly, statistical analysis shows a well consistency between theoretical
composition of stem and actual pattern of EHEC. The results provide a chance to
design gradient control strategies for controlling disease global diffusion.
Our analysis provided the first inspect of global epidemics mediated by trade
networks for improved control and immunity strategies in the future.


  Biomolecular logic systems processing biochemical input signals and producing
"digital" outputs in the form of YES/NO were developed for analysis of
physiological conditions characteristic of liver injury, soft tissue injury and
abdominal trauma. Injury biomarkers were used as input signals for activating
the logic systems. Their normal physiological concentrations were defined as
logic-0 level, while their pathologically elevated concentrations were defined
as logic-1 values. Since the input concentrations applied as logic 0 and 1
values were not sufficiently different, the output signals being at low and
high values (0, 1 outputs) were separated with a short gap making their
discrimination difficult. Coupled enzymatic reactions functioning as a
biomolecular signal processing system with a built-in filter property were
developed. The filter process involves a partial back-conversion of the
optical-output-signal-yielding product, but only at its low concentrations,
thus allowing the proper discrimination between 0 and 1 output values.


  Populations of mobile and communicating agents describe a vast array of
technological and natural systems, ranging from sensor networks to animal
groups. Here, we investigate how a group-level agreement may emerge in the
continuously evolving network defined by the local interactions of the moving
individuals. We adopt a general scheme of motion in two dimensions and we let
the individuals interact through the minimal naming game, a prototypical scheme
to investigate social consensus. We distinguish different regimes of
convergence determined by the emission range of the agents and by their
mobility, and we identify the corresponding scaling behaviors of the consensus
time. In the same way, we rationalize also the behavior of the maximum memory
used during the convergence process, which determines the minimum
cognitive/storage capacity needed by the individuals. Overall, we believe that
the simple and general model presented in this paper can represent a helpful
reference for a better understanding of the behavior of populations of mobile
agents.


  Understanding relationship between noisy dynamics and biological network
architecture is a fundamentally important question, particularly in order to
elucidate how cells encode and process information. We analytically and
numerically investigate general network architectural conditions that are
necessary to generate stochastic amplified and coherent oscillations. We
enumerate all possible topologies of coupled negative feedbacks in the
underlying biochemical networks with three components, negative feedback loops,
and mass action kinetics. Using the linear noise approximation to analytically
obtain the time-dependent solution of the master equation and derive the
algebraic expression of power spectra, we find that (a) all networks with
coupled negative feedbacks are capable of generating stochastic amplified and
coherent oscillations; (b) networks with a single negative feedback are better
stochastic amplified and coherent oscillators than those with multiple coupled
negative feedbacks; (c) multiple timescale difference among the kinetic rate
constants is required for stochastic amplified and coherent oscillations.


  Understanding how an extinction event affects ecosystem is fundamental to
biodiversity conservation. For this reason, food web response to species loss
has been investigated in several ways in the last years. Several studies
focused on secondary extinction due to biodiversity loss in a bottom-up
perspective using in-silico extinction experiments in which a single species is
removed at each step and the number of secondary extinctions is recorded. In
these binary simulations a species goes secondarily extinct if it loses all its
resource species, that is, when the energy intake is zero. This pure
topological statement represents the best case scenario. In fact a consumer
species could go extinct losing a certain fraction of the energy intake and the
response of quantitative food webs to node loss could be very different with
respect to simple binary predictions. The goal of this paper is to analyze how
patterns of secondary extinctions change when higher species sensitivity are
included in the analyses. In particular, we explored how food web secondary
extinction, triggered by the removal of most connected nodes, varies as a
function of the energy intake threshold assumed as the minimum needed for
species persistence. As we will show, a very low increase of energy intake
threshold stimulates a disproportionate growth of secondary extinction.


  In the study of disease spreading on empirical complex networks in SIR model,
initially infected nodes can be ranked according to some measure of their
epidemic impact. The highest ranked nodes, also referred to as
"superspreaders", are associated to dominant epidemic risks and therefore
deserve special attention. In simulations on studied empirical complex
networks, it is shown that the ranking depends on the dynamical regime of the
disease spreading. A possible mechanism leading to this dependence is
illustrated in an analytically tractable example. In systems where the
allocation of resources to counter disease spreading to individual nodes is
based on their ranking, the dynamical regime of disease spreading is frequently
not known before the outbreak of the disease. Therefore, we introduce a
quantity called epidemic centrality as an average over all relevant regimes of
disease spreading as a basis of the ranking. A recently introduced concept of
phase diagram of epidemic spreading is used as a framework in which several
types of averaging are studied. The epidemic centrality is compared to
structural properties of nodes such as node degree, k-cores and betweenness.
There is a growing trend of epidemic centrality with degree and k-cores values,
but the variation of epidemic centrality is much smaller than the variation of
degree or k-cores value. It is found that the epidemic centrality of the
structurally peripheral nodes is of the same order of magnitude as the epidemic
centrality of the structurally central nodes. The implications of these
findings for the distributions of resources to counter disease spreading are
discussed.


  We consider a chemical reaction network governed by mass action kinetics and
composed of N different species which can reversibly form heterodimers. A fast
iterative algorithm is introduced to compute the equilibrium concentrations of
such networks. We show that the convergence is guaranteed by the Banach fixed
point theorem. As a practical example, of relevance for a quantitative analysis
of microarray data, we consider a reaction network formed by N~10^6 mutually
hybridizing different mRNA sequences. We show that, despite the large number of
species involved, the convergence to equilibrium is very rapid for most
species. The origin of slow convergence for some specific subnetworks is
discussed. This provides some insights for improving the performance of the
algorithm.


  The colonization of unoccupied territory by invading species, known as range
expansion, is a spatially heterogeneous non-equilibrium growth process. We
introduce a two-species Eden growth model to analyze the interplay between
uni-directional (irreversible) mutations and selection at the expanding front.
While the evolutionary dynamics leads to coalescence of both wild-type and
mutant clusters, the non-homogeneous advance of the colony results in a rough
front. We show that roughening and domain dynamics are strongly coupled,
resulting in qualitatively altered bulk and front properties. For beneficial
mutations the front is quickly taken over by mutants and growth proceeds
Eden-like. In contrast, if mutants grow slower than wild-types, there is an
antagonism between selection pressure against mutants and growth by merging of
mutant domains with an ensuing absorbing state phase transition to an
all-mutant front. We find that surface roughening has a marked effect on the
critical properties of the absorbing state phase transition. While reference
models, which keep the expanding front flat, exhibit directed percolation
critical behavior, the exponents of the two-species Eden model strongly deviate
from it. In turn, the mutation-selection process induces an increased surface
roughness with exponents distinct from that of the classical Eden model.


  The architecture of bipartite networks linking two classes of constituents is
affected by the interactions within each class. For the bipartite networks
representing the mutualistic relationship between pollinating animals and
plants, it has been known that their degree distributions are broad but often
deviate from power-law form, more significantly for plants than animals. Here
we consider a model for the evolution of the mutualistic networks and find that
their topology is strongly dependent on the asymmetry and non-linearity of the
preferential selection of mutualistic partners. Real-world mutualistic networks
analyzed in the framework of the model show that a new animal species
determines its partners not only by their attractiveness but also as a result
of the competition with pre-existing animals, which leads to the
stretched-exponential degree distributions of plant species.


  The vast majority of mutations are deleterious, and are eliminated by
purifying selection. Yet in finite asexual populations, purifying selection
cannot completely prevent the accumulation of deleterious mutations due to
Muller's ratchet: once lost by stochastic drift, the most-fit class of
genotypes is lost forever. If deleterious mutations are weakly selected,
Muller's ratchet turns into a mutational "meltdown" leading to a rapid
degradation of population fitness. Evidently, the long term stability of an
asexual population requires an influx of beneficial mutations that continuously
compensate for the accumulation of the weakly deleterious ones. Here we propose
that the stable evolutionary state of a population in a static environment is a
dynamic mutation-selection balance, where accumulation of deleterious mutations
is on average offset by the influx of beneficial mutations. We argue that this
state exists for any population size N and mutation rate $U$. Assuming that
beneficial and deleterious mutations have the same fitness effect s, we
calculate the fraction of beneficial mutations, \epsilon, that maintains the
balanced state. We find that a surprisingly low \epsilon suffices to maintain
stability, even in small populations in the face of high mutation rates and
weak selection. This may explain the maintenance of mitochondria and other
asexual genomes, and has implications for the expected statistics of genetic
diversity in these populations.


  We present Molecular Dynamics simulations of Chymotrypsin inhibitor II and
PEO in presence of compatible solutes. Our results indicate that the the native
compact state of the studied macromolecules is stabilized in presence of
hydroxyectoine. We are able to explain the corresponding mechanism by a
variation of the solvent properties for higher hydroxyectoine concentrations.
Our results are validated by detailed free energy calculations.


  The process of evolutionary diversification unfolds in a vast genotypic space
of potential outcomes. During the past century there have been remarkable
advances in the development of theory for this diversification, and the
theory's success rests, in part, on the scope of its applicability. A great
deal of this theory focuses on a relatively small subset of the space of
potential genotypes, chosen largely based on historical or contemporary
patterns, and then predicts the evolutionary dynamics within this pre-defined
set. To what extent can such an approach be pushed to a broader perspective
that accounts for the potential open-endedness of evolutionary diversification?
There have been a number of significant theoretical developments along these
lines but the question of how far such theory can be pushed has not been
addressed. Here a theorem is proven demonstrating that, because of the digital
nature of inheritance, there are inherent limits on the kinds of questions that
can be answered using such an approach. In particular, even in extremely simple
evolutionary systems a complete theory accounting for the potential
open-endedness of evolution is unattainable unless evolution is progressive.
The theorem is closely related to G\"odel's Incompleteness Theorem and to the
Halting Problem from computability theory.


  Oscillation has an important role in bio-dynamical systems such as circadian
rhythms and eukaryotic cell cycle. John Tyson et. al. in Nature Review Mol Cell
Biol 2008 examined a limited number of network topologies consisting of three
nodes and four or fewer edges and identified the network design principles of
biochemical oscillations. Tsai et. al. in Science 2008 studied three different
network motifs, namely a negative feedback loop, coupled negative feedback
loops, and coupled positive and negative feedback loops, and found that the
interconnected positive and negative feedback loops are capable of generating
frequency-tunable oscillations. We enumerate 249 topologically unique network
architectures consisting of three nodes and at least three cyclic inhibitory
edges, and identify network architectural commonalities among three functional
groups: (1) most frequency-tunable yet less robust oscillators, (2) least
frequency-tunable and least robust oscillators, and (3) less frequency-tunable
yet most robust oscillators. We find that Frequency-tunable networks cannot
simultaneously express high robustness, indicating a tradeoff between frequency
tunability and robustness.


  The transient fluctuation of the prosperity of firms in a network economy is
investigated with an abstract stochastic model. The model describes the profit
which firms make when they sell materials to a firm which produces a product
and the fixed cost expense to the firms to produce those materials and product.
The formulae for this model are parallel to those for population dynamics. The
swinging changes in the fluctuation in the transient state from the initial
growth to the final steady state are the consequence of a topology-dependent
time trial competition between the profitable interactions and expense. The
firm in a sparse random network economy is more likely to go bankrupt than
expected from the value of the limit of the fluctuation in the steady state,
and there is a risk of failing to reach by far the less fluctuating steady
state.


  In the present work we analyze the problem of adaptation and evolution of RNA
virus populations, by defining the basic stochastic model as a multivariate
branching process in close relation with the branching process advanced by
Demetrius, Schuster and Sigmund ("Polynucleotide evolution and branching
processes", Bull. Math. Biol. 46 (1985) 239-262), in their study of
polynucleotide evolution. We show that in the absence of beneficial forces the
model is exactly solvable. As a result it is possible to prove several key
results directly related to known typical properties of these systems like (i)
proof, in the context of the theory of branching processes, of the lethal
mutagenesis criterion proposed by Bull, Sanju\'an and Wilke ("Theory of lethal
mutagenesis for viruses", J. Virology 18 (2007) 2930-2939); (ii) a new proposal
for the notion of relaxation time with a quantitative prescription for its
evaluation and (iii) the quantitative description of the evolution of the
expected values in four distinct regimes: transient, "stationary" equilibrium,
extinction threshold and lethal mutagenesis. Moreover, new insights on the
dynamics of evolving virus populations can be foreseen.


  Large molecules such as proteins and nucleic acids are crucial for life, yet
their primordial origin remains a major puzzle. The production of large
molecules, as we know it today, requires good catalysts, and the only good
catalysts we know that can accomplish this task consist of large molecules.
Thus the origin of large molecules is a chicken and egg problem in chemistry.
Here we present a mechanism, based on autocatalytic sets (ACSs), that is a
possible solution to this problem. We discuss a mathematical model describing
the population dynamics of molecules in a stylized but prebiotically plausible
chemistry. Large molecules can be produced in this chemistry by the coalescing
of smaller ones, with the smallest molecules, the `food set', being buffered.
Some of the reactions can be catalyzed by molecules within the chemistry with
varying catalytic strengths. Normally the concentrations of large molecules in
such a scenario are very small, diminishing exponentially with their size.
ACSs, if present in the catalytic network, can focus the resources of the
system into a sparse set of molecules. ACSs can produce a bistability in the
population dynamics and, in particular, steady states wherein the ACS molecules
dominate the population. However to reach these steady states from initial
conditions that contain only the food set typically requires very large
catalytic strengths, growing exponentially with the size of the catalyst
molecule. We present a solution to this problem by studying `nested ACSs', a
structure in which a small ACS is connected to a larger one and reinforces it.
We show that when the network contains a cascade of nested ACSs with the
catalytic strengths of molecules increasing gradually with their size (e.g., as
a power law), a sparse subset of molecules including some very large molecules
can come to dominate the system.


  Recently, research that focuses on the rigorous understanding of the relation
between simulation and/or exact models on graphs and approximate counterparts
has gained lots of momentum. This includes revisiting the performance of
classic pairwise models with closures at the level of pairs and/or triples as
well as effective-degree-type models and those based on the probability
generating function formalism. In this paper, for a fully connected graph and
the simple $SIS$ (susceptible-infected-susceptible) epidemic model, a novel
closure is introduced. This is done via using the equations for the moments of
the distribution describing the number of infecteds at all times combined with
the empirical observations that this is well described/approximated by a
binomial distribution with time dependent parameters. This assumption allows us
to express higher order moments in terms of lower order ones and this leads to
a new closure. The significant feature of the new closure is that the
difference of the exact system, given by the Kolmogorov equations, from the
solution of the newly defined approximate system is of order $1/N^2$. This is
in contrast with the $\mathcal{O}(1/N)$ difference corresponding to the
approximate system obtained via the classic triple closure.


  Spatial extent is a complicating factor in mathematical biology. The
possibility that an action at point A cannot immediately affect what happens at
point B creates the opportunity for spatial nonuniformity. This nonuniformity
must change our understanding of evolutionary dynamics, as the same organism in
different places can have different expected evolutionary outcomes. Since
organism origins and fates are both determined locally, we must consider
heterogeneity explicitly to determine its effects. We use simulations of
spatially extended host--pathogen and predator--prey ecosystems to reveal the
limitations of standard mathematical treatments of spatial heterogeneity. Our
model ecosystem generates heterogeneity dynamically; an adaptive network of
hosts on which pathogens are transmitted arises as an emergent phenomenon. The
structure and dynamics of this network differ in significant ways from those of
related models studied in the adaptive-network field. We use a new technique,
organism swapping, to test the efficacy of both simple approximations and more
elaborate moment-closure methods, and a new measure to reveal the timescale
dependence of invasive-strain behavior. Our results demonstrate the failure not
only of the most straightforward ("mean field") approximation, which smooths
over heterogeneity entirely, but also of the standard correction ("pair
approximation") to the mean field treatment. In spatial contexts, invasive
pathogen varieties can prosper initially but perish in the medium term,
implying that the concepts of reproductive fitness and the Evolutionary Stable
Strategy have to be modified for such systems.


  In a recent paper, Klaere et al. modeled the impact of substitutions on
arbitrary branches of a phylogenetic tree on an alignment site by the so-called
One Step Mutation (OSM) matrix. By utilizing the concept of the OSM matrix for
the four-state nucleotide alphabet, Nguyen et al. presented an efficient
procedure to compute the minimal number of substitutions needed to translate
one alignment site into another.The present paper delivers a proof for this
computation.Moreover, we provide several mathematical insights into the
generalization of the OSM matrix to multistate alphabets.The construction of
the OSM matrix is only possible if the matrices representing the substitution
types acting on the character states and the identity matrix form a commutative
group with respect to matrix multiplication. We illustrate a means to establish
such a group for the twenty-state amino acid alphabet and critically discuss
its biological usefulness.


  In this paper we present a model describing Susceptible-Infected-Susceptible
(SIS) type epidemics spreading on a dynamic contact network with random link
activation and deletion where link ac- tivation can be locally constrained. We
use and adapt a improved effective degree compartmental modelling framework
recently proposed by Lindquist et al. [J. Lindquist et al., J. Math Biol. 62,
2, 143 (2010)] and Marceau et al. [V. Marceau et al., Phys. Rev. E 82, 036116
(2010)]. The resulting set of ordinary differential equations (ODEs) is solved
numerically and results are compared to those obtained using individual-based
stochastic network simulation. We show that the ODEs display excellent
agreement with simulation for the evolution of both the disease and the
network, and is able to accurately capture the epidemic threshold for a wide
range of parameters. We also present an analytical R0 calculation for the
dynamic network model and show that depending on the relative timescales of the
network evolution and disease transmission two limiting cases are recovered:
(i) the static network case when network evolution is slow and (ii) homogeneous
random mixing when the network evolution is rapid. We also use our threshold
calculation to highlight the dangers of relying on local stability analysis
when predicting epidemic outbreaks on evolving networks.


  In this paper we study the susceptible-infectious (SI) epidemiological model
using dynamical graphs. Dynamical structures have been recently applied in many
areas including complex systems. Dynamical structures include the mutual
interaction between the structure topology and the characteristics of its
members. Dynamical graphs applied to epidemics consider generally that the
nodes are individuals and the links represent different classes of
relationships between individuals with the potential to transmit the disease.
The main aim in this article is to study the evolution of the SI
epidemiological model and the creation of subgraphs due to the dynamic behavior
of the individuals trying to avoid the contagious of the disease. The proposed
dynamical graph model uses a single parameter which reflects the probability of
rewire that represent actions to avoid the disease. This parameter includes
also information regarding the infectivity of the disease. The numerical
simulations using Monte Carlo method show that the dynamical behavior of
individuals affects the evolution of the subgraphs. Furthermore, it is shown
that the connectivity degree of the graphs can change the arise of subgraphs
and the asymptotic state of the infectious diseases.


  Since its first publication in 2003, the Gene Set Enrichment Analysis (GSEA)
method, based on the Kolmogorov-Smirnov statistic, has been heavily used,
modified, and also questioned. Recently a simplified approach, using a one
sample t test score to assess enrichment and ignoring gene-gene correlations
was proposed by Irizarry et al. 2009 as a serious contender. The argument
criticizes GSEA's nonparametric nature and its use of an empirical null
distribution as unnecessary and hard to compute. We refute these claims by
careful consideration of the assumptions of the simplified method and its
results, including a comparison with GSEA's on a large benchmark set of 50
datasets. Our results provide strong empirical evidence that gene-gene
correlations cannot be ignored due to the significant variance inflation they
produced on the enrichment scores and should be taken into account when
estimating gene set enrichment significance. In addition, we discuss the
challenges that the complex correlation structure and multi-modality of gene
sets pose more generally for gene set enrichment methods.


  Imposing a minimum principle in the framework of the so called crystal basis
model of the genetic code, we determine the structure of the minimum set of
anticodons which allows the translational-transcription for animal
mitochondrial code. The results are in very good agreement with the observed
anticodons.


  For a population with any given number of types, we construct a new
multivariate Moran process with frequency-dependent selection and establish,
analytically, a correspondence to equilibrium Lotka-Volterra phenomenology.
This correspondence, on the one hand, allows us to infer the phenomenology of
our Moran process based on much simpler Lokta-Volterra phenomenology, and on
the other, allows us to study Lotka-Volterra dynamics within the finite
populations of a Moran process. Applications to community ecology, population
genetics, and evolutionary game theory are discussed.


  To understand the mechanism of trimethylamine N-oxide (TMAO) induced
stabilization of folded protein states, we systematically investigated the
action of TMAO on several model dipeptides (Leucine, L2, Serine, S2, Glutamine,
Q2, Lysine, K2, and Glycine, G2) in order to elucidate the effect of
residue-specific TMAO interactions on small fragments of solvent-exposed
conformations of the denatured states of proteins. We find that TMAO
preferentially hydrogen bonds with the exposed dipeptide backbone, but
generally not with nonpolar or polar side chains. However, interactions with
the positively charged Lys are substantially greater than with the backbone.
The dipeptide G2, is a useful model of pure amide backbone, interacts with TMAO
by forming a hydrogen bond between the amide nitrogen and the oxygen in TMAO.
In contrast, TMAO is depleted from the protein backbone in the hexapeptide G6,
which shows that the length of the polypeptide chain is relevant in aqueous
TMAO solutions. These simulations lead to the hypothesis that TMAO-induced
stabilization of proteins and peptides is a consequence of depletion of the
solute from the protein surface provided intramolecular interactions are more
favorable than those between TMAO and the backbone. To test our hypothesis we
performed additional simulations of the action of TMAO on an intrinsically
disordered A{\beta}16-22 (KLVFFAE) monomer. In the absence of TMAO
A{\beta}16-22 is a disordered random coil. However, in aqueous TMAO solution
A{\beta}16-22 monomer samples compact conformations. A transition from random
coil to {\alpha}-helical secondary structure is observed at high TMAO
concentrations. Our work highlights the potential similarities between the
action of TMAO on long polypeptide chains and entropic stabilization of
proteins in a crowded environment due to excluded volume interactions. In this
sense TMAO is a nano-crowding particle.


  The choice of the summary statistics used in Bayesian inference and in
particular in ABC algorithms has bearings on the validation of the resulting
inference. Those statistics are nonetheless customarily used in ABC algorithms
without consistency checks. We derive necessary and sufficient conditions on
summary statistics for the corresponding Bayes factor to be convergent, namely
to asymptotically select the true model. Those conditions, which amount to the
expectations of the summary statistics to asymptotically differ under both
models, are quite natural and can be exploited in ABC settings to infer whether
or not a choice of summary statistics is appropriate, via a Monte Carlo
validation.


  We consider networks with two types of nodes. The v-nodes, called centers,
are hyperconnected and interact one to another via many u-nodes, called
satellites. This centralized architecture, widespread in gene networks, realize
a bow-tie scheme and possesses interesting properties. Namely, this
organization creates feedback loops that are capable to generate any prescribed
patterning dynamics, chaotic or periodic, and create a number of equilibrium
states. We show that activation or silencing of a node can sharply switch the
network attractor, even if the activated or silenced node is weakly connected.
We distinguish between two dynamically different situations, "power of center"
(PC) when satellite response is fast and "satellite power" (SP) when center
response is fast. Using a simple network example we show that a centralized
network is more robust with respect to time dependent perturbations, in the PC
relative to the SP case. In theoretical molecular biology, this class of models
can be used to reveal a non-trivial relation between the architecture of
protein-DNA and protein-protein interaction networks and controllability of
space-time dynamics of cellular processes.


  We investigate an analytically tractable toy model for thermally induced
polymorphic dynamics of cooperatively rearranging biofilaments - like
microtubules. The proposed 4 -block model, which can be seen as a
coarse-grained approximation of the full polymorphic tube model, permits a
complete analytical treatment of all thermodynamic properties including
correlation functions and angular fourier mode distributions. Due to its
mathematical tractability the model straightforwardly leads to some physical
insights in recently discussed phenomena like the "length dependent persistence
length". We show that a polymorphic filament can disguise itself as a classical
worm like chain on small and on large scales and yet display distinct anomalous
tell-tale features indicating an inner switching dynamics on intermediate
length scales.


  The evolutionary trajectory of a protein through sequence space is
constrained by function and three-dimensional (3D) structure. Residues in
spatial proximity tend to co-evolve, yet attempts to invert the evolutionary
record to identify these constraints and use them to computationally fold
proteins have so far been unsuccessful. Here, we show that co-variation of
residue pairs, observed in a large protein family, provides sufficient
information to determine 3D protein structure. Using a data-constrained maximum
entropy model of the multiple sequence alignment, we identify pairs of
statistically coupled residue positions which are expected to be close in the
protein fold, termed contacts inferred from evolutionary information (EICs). To
assess the amount of information about the protein fold contained in these
coupled pairs, we evaluate the accuracy of predicted 3D structures for proteins
of 50-260 residues, from 15 diverse protein families, including a G-protein
coupled receptor. These structure predictions are de novo, i.e., they do not
use homology modeling or sequence-similar fragments from known structures. The
resulting low C{\alpha}-RMSD error range of 2.7-5.1{\AA}, over at least 75% of
the protein, indicates the potential for predicting essentially correct 3D
structures for the thousands of protein families that have no known structure,
provided they include a sufficiently large number of divergent sample
sequences. With the current enormous growth in sequence information based on
new sequencing technology, this opens the door to a comprehensive survey of
protein 3D structures, including many not currently accessible to the
experimental methods of structural genomics. This advance has potential
applications in many biological contexts, such as synthetic biology,
identification of functional sites in proteins and interpretation of the
functional impact of genetic variants.


  Using the shape space of codons and tRNAs we give a physical description of
the genetic code evolution on the basis of the codon capture and ambiguous
intermediate scenarios in a consistent manner. In the lowest dimensional
version of our description, a physical quantity, codon level is introduced. In
terms of the codon levels two scenarios are typically classified into two
different routes of the evolutional process. In the case of the ambiguous
intermediate scenario we perform an evolutional simulation implemented cost
selection of amino acids and confirm a rapid transition of the code change.
Such rapidness reduces uncomfortableness of the non-unique translation of the
code at intermediate state that is the weakness of the scenario. In the case of
the codon capture scenario the survival against mutations under the mutational
pressure minimizing GC content in genomes is simulated and it is demonstrated
that cells which experience only neutral mutations survive.


  Structure fluctuations and conformational changes accompany all biological
processes involving macromolecules. The paper presents a classification of
protein residues based on the normalized equilibrium fluctuations of the
residue centers of mass in proteins and a statistical analysis of conformation
changes in the side-chains upon binding. Normal mode analysis and an elastic
network model were applied to a set of protein complexes to calculate the
residue fluctuations and develop the residue classification. Comparison with a
classification based on normalized B-factors suggests that the B-factors may
underestimate protein flexibility in solvent. Our classification shows that
protein loops and disordered fragments are enriched with highly fluctuating
residues and depleted with weakly fluctuating residues. To calculate the
dihedral angles distribution functions, the configuration space was divided
into cells by a cubic grid. The effect of protein association on the
distribution functions depends on the amino acid type and a grid step in the
dihedral angles space. The changes in the dihedral angles increase from the
near-backbone dihedral angle to the most distant one, for most residues. On
average, one fifth of the interface residues change the rotamer state upon
binding, whereas the rest of the interface residues undergo local readjustments
within the same rotamer.


  The similarity in the three-dimensional structures of homologous proteins
imposes strong constraints on their sequence variability. It has long been
suggested that the resulting correlations among amino acid compositions at
different sequence positions can be exploited to infer spatial contacts within
the tertiary protein structure. Crucial to this inference is the ability to
disentangle direct and indirect correlations, as accomplished by the recently
introduced Direct Coupling Analysis (DCA) (Weigt et al. (2009) Proc Natl Acad
Sci 106:67). Here we develop a computationally efficient implementation of DCA,
which allows us to evaluate the accuracy of contact prediction by DCA for a
large number of protein domains, based purely on sequence information. DCA is
shown to yield a large number of correctly predicted contacts, recapitulating
the global structure of the contact map for the majority of the protein domains
examined. Furthermore, our analysis captures clear signals beyond intra- domain
residue contacts, arising, e.g., from alternative protein conformations,
ligand- mediated residue couplings, and inter-domain interactions in protein
oligomers. Our findings suggest that contacts predicted by DCA can be used as a
reliable guide to facilitate computational predictions of alternative protein
conformations, protein complex formation, and even the de novo prediction of
protein domain structures, provided the existence of a large number of
homologous sequences which are being rapidly made available due to advances in
genome sequencing.


  Programmed cell death - apoptosis is one of the most studied biological
phenomenon of recent years. Apoptotic regulatory network contains several
significant control points, including probably the most important one - Bcl--2
apoptotic switch. There are two proposed hypotheses regarding its internal
working - the indirect activation and direct activation models. Since these
hypotheses form extreme poles of full continuum of intermediate models, we have
constructed more general model with these two models as extreme cases.
  By studying relationship between model parameters and steady-state response
ultrasensitivity we have found optimal interaction pattern which reproduces
behavior of Bcl-2 apoptotic switch. Our results show, that stimulus-response
ultrasensitivity is negatively related to spontaneous activation of Bcl-2
effectors - subgroup of Bcl-2 proteins. We found that ultrasensitivity requires
effector's activation, mediated by another subgroup of Bcl-2 proteins -
activators. We have shown that the auto-activation of effectors forms
ultrasensitivity enhancing feedback loop, only if mediated by monomers, but not
by oligomers. Robustness analysis revealed that interaction pattern proposed by
direct activation hypothesis is able to conserve stimulus-response dependence
and preserve ultrasensitivity despite large changes of its internal parameters.
This ability is strongly reduced as for the intermediate to indirect side of
the models.
  Computer simulation of the more general model presented here suggest, that
stimulus-response ultrasensitivity is an emergent property of the direct
activation model, that cannot originate within model of indirect activation.
Introduction of indirect-model-specific interactions does not provide better
explanation of Bcl-2 functioning compared to direct model.


  We outline the global control architecture of genomes. A theory of genomic
control information is presented. The concept of a developmental control
network called a cene (for control gene) is introduced. We distinguish
parts-genes from control genes or cenes. Cenes are interpreted and executed by
the cell and, thereby, direct cell actions including communication, growth,
division, differentiation and multi-cellular development. The cenome is the
global developmental control network in the genome. The cenome is also a cene
that consists of interlinked sub-cenes that guide the ontogeny of the organism.
The complexity of organisms is linked to the complexity of the cenome. The
relevance to ontogeny and evolution is mentioned. We introduce the concept of a
universal cell and a universal genome.


  We calculate both the exponential and pre-factor contributions in a WKB
approximation of the master equation for a stochastic SIR model with highly
oscillatory dynamics. Fixing the basic parameters of the model we investigate
how the outbreak distribution changes with the population size. We show that
this distribution rapidly becomes highly non-Gaussian, acquiring large tails
indicating the presence of rare, but large outbreaks, as the population is made
smaller. The analytic results are found to be in excellent agreement with
simulations until the systems become so small that the dynamics are dominated
by fade-out of the disease.


  It is widely accepted that population genetics theory is the cornerstone of
evolutionary analyses. Empirical tests of the theory, however, are challenging
because of the complex relationships between space, dispersal, and evolution.
Critically, we lack quantitative validation of the spatial models of population
genetics. Here we combine analytics, on and off-lattice simulations, and
experiments with bacteria to perform quantitative tests of the theory. We study
two bacterial species, the gut microbe Escherichia coli and the opportunistic
pathogen Pseudomonas aeruginosa, and show that spatio-genetic patterns in
colony biofilms of both species are accurately described by an extension of the
one-dimensional stepping-stone model. We use one empirical measure, genetic
diversity at the colony periphery, to parameterize our models and show that we
can then accurately predict another key variable: the degree of short-range
cell migration along an edge. Moreover, the model allows us to estimate other
key parameters including effective population size (density) at the expansion
frontier. While our experimental system is a simplification of natural
microbial community, we argue it is a proof of principle that the spatial
models of population genetics can quantitatively capture organismal evolution.


  Ecologists and conservation biologists need to identify the relative
importance of species to make sound management decisions and effectively
allocate scarce resources. We introduce a new method, termed environ
centrality, to determine the relative importance of a species in an ecosystem
network with respect to ecosystem energy--matter exchange. We demonstrate the
uniqueness of environ centrality by comparing it to other common centrality
metrics and then show its ecological significance. Specifically, we tested two
hypotheses on a set of 50 empirically-based ecosystem network models. The first
concerned the distribution of centrality in the community. We hypothesized that
the functional importance of species would tend to be concentrated into a few
dominant species followed by a group of species with lower, more even
importance as is often seen in dominance--diversity curves. Second, we tested
the systems ecology hypothesis that indirect relationships homogenize the
functional importance of species in ecosystems. Our results support both
hypotheses and highlight the importance of detritus and nutrient recyclers such
as fungi and bacteria in generating the energy--matter flow in ecosystems. Our
homogenization results suggest that indirect effects are important in part
because they tend to even the importance of species in ecosystems. A core
contribution of this work is that it creates a formal, mathematical method to
quantify the importance species play in generating ecosystem activity by
integrating direct, indirect, and boundary effects in ecological systems.


  Blebs are balloon-shaped membrane protrusions that form during many
physiological processes. Using computer simulation of a particle-based model
for self-assembled lipid bilayers coupled to an elastic meshwork, we
investigated the phase behavior and kinetics of blebbing. We found that blebs
form for large values of the ratio between the areas of the bilayer and the
cytoskeleton. We also found that blebbing can be induced when the cytoskeleton
is subject to a localized ablation or a uniform compression. The results
obtained are qualitatively in agreement with the experimental evidence and the
model opens up the possibility to study the kinetics of bleb formation in
detail.


  Recent studies have shown that adaptive networks driven by simple local rules
can organize into "critical" global steady states, providing another framework
for self-organized criticality (SOC). We focus on the important convergence to
criticality and show that noise and time-scale optimality are reached at finite
values. This is in sharp contrast to the previously believed optimal zero noise
and infinite time scale separation case. Furthermore, we discover a noise
induced phase transition for the breakdown of SOC. We also investigate each of
the three new effects separately by developing models. These models reveal
three generically low-dimensional dynamical behaviors: time-scale resonance
(TR), a new simplified version of stochastic resonance - which we call steady
state stochastic resonance (SSR) - as well as noise-induced phase transitions.


  Thermal fluctuations of microtubules (MTs) and other cytoskeletal filaments
govern to a great extent the complex rheological properties of the cytoskeleton
in eukaryotic cells. In recent years, much effort has been put into capturing
the dynamics of these fluctuations by means of analytical and numerical models.
These attempts have been very successful for, but also remain limited to,
isotropic polymers. To correctly interpret experimental work on (strongly)
anisotropic semiflexible polymers, there is a need for a numerical modelling
tool that accurately captures the dynamics of polymers with anisotropic
material properties. In the current study, we present a finite element (FE)
framework for simulating the thermal dynamics of a single anisotropic
semiflexible polymer. First, we demonstrate the accuracy of our framework by
comparison of the simulated mean square displacement (MSD) of the end-to-end
distance with analytical predictions based on the worm-like chain model. Then,
we implement a transversely isotropic material model, characteristic for
biopolymers such as MTs, and study the persistence length for various ratios
between the longitudinal shear modulus and corresponding Young's modulus.
Finally, we put our findings in context by addressing a recent experimental
work on grafted transversely isotropic MTs. In that research, a simplified
static mechanical model was used to deduce a very high level of MT anisotropy
to explain the observation that the persistence length of grafted MTs increases
as contour length increases. We show, by means of our FE framework, that the
anisotropic properties cannot account for the reported length-dependent
persistence length.


This paper has been withdrawn by the author due to a missing figure


  In many cellular signaling pathways, key components form clusters at the cell
membrane. Although much work has focused on the mechanisms behind such cluster
formation, the implications for downstream signaling remain poorly understood.
Here, motivated by recent experiments, we study via particle-based simulation a
covalent modification network in which the activating component is either
clustered or randomly distributed on the membrane. We find that while
clustering reduces the response of a single-modification network, clustering
can enhance the response of a double-modification network. The reduction is a
bulk effect: a cluster presents a smaller effective target to a substrate
molecule in the bulk. The enhancement, on the other hand, is a local effect: a
cluster promotes the rapid rebinding and second activation of singly active
substrate molecules. As such, the enhancement relies upon frequent collisions
on a short timescale, which leads to a diffusion coefficient at which the
enhancement is optimal. We complement simulation with analytic results at both
the mean-field and first-passage distribution levels. Our results emphasize the
importance of spatially resolved models, showing that significant effects of
spatial correlations persist even in spatially averaged quantities such as
response curves.


  Determining the interaction partners among protein/domain families poses hard
computational problems, in particular in the presence of paralogous proteins.
Available approaches aim to identify interaction partners among protein/domain
families through maximizing the similarity between trimmed versions of their
phylogenetic trees. Since maximization of any natural similarity score is
computationally difficult, many approaches employ heuristics to maximize the
distance matrices corresponding to the tree topologies in question. In this
paper we devise an efficient deterministic algorithm which directly maximizes
the similarity between two leaf labeled trees with edge lengths, obtaining a
score-optimal alignment of the two trees in question.
  Our algorithm is significantly faster than those methods based on distance
matrix comparison: 1 minute on a single processor vs. 730 hours on a
supercomputer. Furthermore we have advantages over the current state-of-the-art
heuristic search approach in terms of precision as well as a recently suggested
overall performance measure for mirrortree approaches, while incurring only
acceptable losses in recall.
  A C implementation of the method demonstrated in this paper is available at
http://compbio.cs.sfu.ca/mirrort.htm


  We present a general computational theory of cancer and its developmental
dynamics. The theory is based on a theory of the architecture and function of
developmental control networks which guide the formation of multicellular
organisms. Cancer networks are special cases of developmental control networks.
Cancer results from transformations of normal developmental networks. Our
theory generates a natural classification of all possible cancers based on
their network architecture. Each cancer network has a unique topology and
semantics and developmental dynamics that result in distinct clinical tumor
phenotypes. We apply this new theory with a series of proof of concept cases
for all the basic cancer types. These cases have been computationally modeled,
their behavior simulated and mathematically described using a multicellular
systems biology approach. There are fascinating correspondences between the
dynamic developmental phenotype of computationally modeled {\em in silico}
cancers and natural {\em in vivo} cancers. The theory lays the foundation for a
new research paradigm for understanding and investigating cancer. The theory of
cancer networks implies that new diagnostic methods and new treatments to cure
cancer will become possible.


  S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key
role in methylation reactions including those required for maturation of viral
mRNA, is an important drug target in the discovery of antiviral agents. While
targeting the active site is a straightforward strategy of enzyme inhibition,
evidences of allosteric modulation of active site in many enzymes underscore
the molecular origin of signal transduction. Information of co-evolving
sequences in SAHH family and the key residues for functional dynamics that can
be identified using native topology of the enzyme provide glimpses into how the
allosteric signaling network, dispersed over the molecular structure,
coordinates intra- and inter-subunit conformational dynamics. To study the link
between the allosteric communication and functional dynamics of SAHHs, we
performed Brownian dynamics simulations by building a coarse-grained model
based on the holo and ligand-bound structures. The simulations of
ligand-induced transition revealed that the signal of intra-subunit closure
dynamics is transmitted to form inter-subunit contacts, which in turn invoke a
precise alignment of active site, followed by the dimer-dimer rotation that
compacts the whole tetrameric structure. Further analyses of SAHH dynamics
associated with ligand binding provided evidence of both induced fit and
population shift mechanisms, and also showed that the transition state ensemble
is akin to the ligand-bound state. Besides the formation of enzyme-ligand
contacts at the active site, the allosteric couplings from the residues distal
to the active site is vital to the enzymatic function.


  The Hill coefficient is often used as a direct measure of the cooperativity
of binding processes. It is an essential tool for probing properties of
reactions in many biochemical systems. Here we analyze existing experimental
data and demonstrate that the Hill coefficient characterizing the binding of
transcription factors to their cognate sites can in fact be larger than one --
the standard indication of cooperativity -- even in the absence of any standard
cooperative binding mechanism. By studying the problem analytically, we
demonstrate that this effect occurs due to the disordered binding energy of the
transcription factor to the DNA molecule and the steric interactions between
the different copies of the transcription factor. We show that the enhanced
Hill coefficient implies a significant reduction in the number of copies of the
transcription factors which is needed to occupy a cognate site and, in many
cases, can explain existing estimates for numbers of the transcription factors
in cells. The mechanism is general and should be applicable to other biological
recognition processes.


  Background: A classical problem in metabolic design is to maximize the
production of desired compound in a given chemical reaction network by
appropriately directing the mass flow through the network. Computationally,
this problem is addressed as a linear optimization problem over the "flux
cone". The prior construction of the flux cone is computationally expensive and
no polynomial-time algorithms are known. Results: Here we show that the output
maximization problem in chemical reaction networks is NP-complete. This
statement remains true even if all reactions are monomolecular or bimolecular
and if only a single molecular species is used as influx. As a corollary we
show, furthermore, that the detection of autocatalytic species, i.e., types
that can only be produced from the influx material when they are present in the
initial reaction mixture, is an NP-complete computational problem. Conclusions:
Hardness results on combinatorial problems and optimization problems are
important to guide the development of computational tools for the analysis of
metabolic networks in particular and chemical reaction networks in general. Our
results indicate that efficient heuristics and approximate algorithms need to
be employed for the analysis of large chemical networks since even conceptually
simple flow problems are provably intractable.


  Traditional Quartet Puzzling algorithms use maximum likelihood methods to
reconstruct quartet trees, and a puzzling algorithm to combine these quartets
into a tree for the full collection of $n$ taxa. We propose a variation of
Quartet Puzzling in which the quartet trees are reconstructed using
biologically symmetric invariants. We find that under certain conditions,
invariant based quartet puzzling outperforms Quartet Puzzling using maximum
likelihood.


  Motivation: Proteins are known to undergo conformational changes in the
course of their functions. The changes in conformation are often attributable
to a small fraction of residues within the protein. Therefore identification of
these variable regions is important for an understanding of protein function.
  Results: We propose a novel method for identification of local conformational
changes in proteins. In our method, backbone conformations are encoded into a
sequence of letters from a 16-letter alphabet (called D2 codes) to perform
structural comparison. Since we do not use clustering analysis to encode local
structures, the D2 codes not only provides a intuitively understandable
description of protein structures, but also covers wide varieties of
distortions. This paper shows that the D2 codes are better correlated with
changes in the dihedral angles than a structural alphabet and a secondary
structure description. In the case of the N37S mutant of HIV-1 protease, local
conformational changes were captured by the D2 coding method more accurately
than other methods. The D2 coding also provided a reliable representation of
the difference between NMR models of an HIV-1 protease mutant.


  We show that the effect of stress on the reproductive fitness of noisy cell
populations can be modelled as first-passage time problem, and demonstrate that
even relatively short-lived fluctuations in gene expression can ensure
long-term survival of a drug-resistant population. We examine how this effect
contributes to the development of drug-resistant cancer cells, and demonstrate
that permanent immunity can arise independently of mutations.


  Predicting protein interactions is one of the more interesting challenges of
the post-genomic era. Many algorithms address this problem as a binary
classification problem: given two proteins represented as two vectors of
features, predict if they interact or not. Importantly however, computational
predictions are only one component of a larger framework for identifying PPI.
The most promising candidate pairs can be validated experimentally by testing
if they physical bind to each other. Since these experiments are more costly
and error prone, the computational predictions serve as a filter, aimed to
produce a small number of highly promising candidates. Here we propose to
address this problem as a ranking problem: given a network with known
interactions, rank all unknown pairs based on the likelihood of their
interactions.
  In this paper we propose a ranking algorithm that trains multiple
inter-connected models using a passive aggressive on-line approach. We show
good results predicting protein-protein interactions for post synaptic density
PPI network. We compare the precision of the ranking algorithm with local
classifiers and classic support vector machine. Though the ranking algorithm
outperforms the classic SVM classification, its performance is inferior
compared to the local supervised method.


  The application of the quasi-steady-state approximation to the
Michaelis-Menten reaction embedded in large open chemical reaction networks is
a popular model reduction technique in deterministic and stochastic simulations
of biochemical reactions inside cells. It is frequently assumed that the
predictions of the reduced master equations obtained using the stochastic
quasi-steady-state approach are in very good agreement with the predictions of
the full master equations, provided the conditions for the validity of the
deterministic quasi-steady-state approximation are fulfilled. We here use the
linear-noise approximation to show that this assumption is not generally
justified for the Michaelis-Menten reaction with substrate input, the simplest
example of an open embedded enzyme reaction. The reduced master equation
approach is found to considerably overestimate the size of intrinsic noise at
low copy numbers of molecules. A simple formula is obtained for the relative
error between the predictions of the reduced and full master equations for the
variance of the substrate concentration fluctuations. The maximum error is
reached when modeling moderately or highly efficient enzymes, in which case the
error is approximately 30%. The theoretical predictions are validated by
stochastic simulations using experimental parameter values for enzymes involved
in proteolysis, gluconeogenesis and fermentation.


  The family of TGFb ligands is large and its members are involved in many
different signaling processes. These signaling processes strongly differ in
type with TGFb ligands eliciting both sustained or transient responses. Members
of the TGFb family can also act as morphogen and cellular responses would then
be expected to provide a direct read-out of the extracellular ligand
concentration. We were interested to define the set of minimal modifications
that are required to change the type of signal processing in the TGFb signaling
network. To define the key aspects for signaling plasticity we focused on the
core of the TGFb signaling network. With the help of a parameter screen we
identified ranges of kinetic parameters and protein concentrations that give
rise to transient, sustained, or oscillatory responses to constant stimuli, as
well as those parameter ranges that enable a proportional response to
time-varying ligand concentrations (as expected in the read-out of morphogens).
A combination of a strong negative feedback and fast shuttling to the nucleus
biases signaling to a transient rather than a sustained response, while
oscillations were obtained if ligand binding to the receptor is weak and the
turn-over of the I-Smad is fast. A proportional read-out required inefficient
receptor activation in addition to a low affinity of receptor-ligand binding.
We find that targeted modification of single parameters suffices to alter the
response type. The architecture of the TGFb pathway enables the observed
signaling plasticity. The observed range of signaling outputs to TGFb ligand in
different cell types and under different conditions can be explained with
differences in cellular protein concentrations and with changes in effective
rate constants due to cross-talk with other signaling pathways.


  This paper summarizes previous work linking protein aggregation to the
heterozygosity of organisms. It also cites the literature showing a correlation
between species' morphological complexity and the lengths of their proteins.
These two findings are combined to form a theory that may potentially explain
the ploidy levels of organisms. Organisms can employ heterozygosity to inhibit
protein aggregation. Hence, complex organisms tend to be diploid because they
tend to synthesize long, aggregation-prone proteins. On the other hand, simple
organisms tend to be haploid because they synthesize short proteins that are
less prone to aggregation. The theory may also explain ecological trends
associated with organisms of different ploidy level. Two mathematical models
are also developed that may explain: 1) how protein aggregation results in
truncation selection that maintains numerous polymorphisms in natural
populations, and 2) the relationship between protein turnover, metabolic
efficiency, and heterosis.


  We present the first sub-quadratic time algorithm that with high probability
correctly reconstructs phylogenetic trees for short sequences generated by a
Markov model of evolution. Due to rapid expansion in sequence databases, such
very fast algorithms are becoming necessary. Other fast heuristics have been
developed for building trees from very large alignments (Price et al, and Brown
et al), but they lack theoretical performance guarantees. Our new algorithm
runs in $O(n^{1+\gamma(g)}\log^2n)$ time, where $\gamma$ is an increasing
function of an upper bound on the branch lengths in the phylogeny, the upper
bound $g$ must be below$1/2-\sqrt{1/8} \approx 0.15$, and $\gamma(g)<1$ for all
$g$. For phylogenies with very short branches, the running time of our
algorithm is close to linear. For example, if all branch lengths correspond to
a mutation probability of less than 0.02, the running time of our algorithm is
roughly $O(n^{1.2}\log^2n)$. Via a prototype and a sequence of large-scale
experiments, we show that many large phylogenies can be reconstructed fast,
without compromising reconstruction accuracy.


  We describe a system of stochastic differential equations (SDEs) which model
the interaction between processive molecular motors, such as kinesin and
dynein, and the biomolecular cargo they tow as part of microtubule-based
intracellular transport. We show that the classical experimental environment
fits within a parameter regime which is qualitatively distinct from conditions
one expects to find in living cells. Through an asymptotic analysis of our
system of SDEs, we develop a means for applying in vitro observations of the
nonlinear response by motors to forces induced on the attached cargo to make
analytical predictions for two parameter regimes that have thus far eluded
direct experimental observation: 1) highly viscous in vivo transport and 2)
dynamics when multiple identical motors are attached to the cargo and
microtubule.


  The general time reversible model (GTR) is presently the most popular model
used in phylogentic studies. However, GTR has an undesirable mathematical
property that is potentially of significant concern. It is the purpose of this
article to give examples that demonstrate why this deficit may pose a problem
for phylogenetic analysis and interpretation.


  We model recruitment in adaptive social networks in the presence of birth and
death processes. Recruitment is characterized by nodes changing their status to
that of the recruiting class as a result of contact with recruiting nodes. Only
a susceptible subset of nodes can be recruited. The recruiting individuals may
adapt their connections in order to improve recruitment capabilities, thus
changing the network structure adaptively. We derive a mean field theory to
predict the dependence of the growth threshold of the recruiting class on the
adaptation parameter. Furthermore, we investigate the effect of adaptation on
the recruitment level, as well as on network topology. The theoretical
predictions are compared with direct simulations of the full system. We
identify two parameter regimes with qualitatively different bifurcation
diagrams depending on whether nodes become susceptible frequently (multiple
times in their lifetime) or rarely (much less than once per lifetime).


  This paper reports on a simulated evolution project, which had the goal of
simulating the refractive components of the PM eye of Australian netting spider
diopis subrufus on a desktop computer. The model for the simulation is the
anatomy of the eye described by Blest and Land [Ble & Lan 1977]. The evolution
simulation was able to produce hundreds of eyes with equivalent optical
qualities to the measured eyes for the phenotype of the netting spider. These
artificially evolved eyes began to occur in the computer simulation between
8X106 and 35X106 cycles after Demonstration of Darwinian Theory of Evolution;
arXiv 1006.0480 simulated the evolution for the ctenid spider PM eye,
cupiennius sale. This paper follows the previous paper, but the netting spider
eye is more complex than the ctenid PM eye so the simulated evolution equations
are more complex.


  Membrane tubes are important elements for living cells to organize many
functions. Experiments have found that membrane tube can be extracted from
giant lipid vesicles by a group of kinesin. How these motors cooperate in
extracting the fluid-like membrane tube is still unclear. In this paper, we
propose a new cooperation mechanism called two-track-dumbbell model, in which
kinesin is regarded as a dumbbell with an end (tail domain) tightly bound onto
the fluid-like membrane and the other end (head domain) stepping on or
unbinding from the microtubule. Taking account of the elasticity of kinesin
molecule and the exclude volume effect of both the head domain and the tail
domain of kinesin, which are not considered in previous models, we simulate the
growth process of the membrane tube pulled by kinesin motors. Our results
indicate that motors along a single microtubule protofilament can generate
enough force to extract membrane tubes from vesicles, and the average number of
motors pulling the tube is about 8~9. These results are quite different from
previous studies (Ref. \cite{camp.08}), and further experimental tests are
necessary to elucidate the cooperation mechanism.


  When a collection of phenotypically diverse organisms compete with each other
for limited resources, with competition being strongest amongst the most
similar, the population can evolve into tightly localised clusters. This
process can be thought of as a simple model of the emergence of species. Past
studies have neglected the effects of demographic noise and studied the
population on a macroscopic scale, where species formation is found to depend
upon the shape of the curve describing the decline of competition strength with
phenotypic distance. In the following, we will show how including the effects
of demographic noise leads to a radically different conclusion. Two situations
are identified: a weak-noise regime in which the population exhibits patterns
of fluctuation around the macroscopic description, and a strong-noise regime
where species appear spontaneously even in the case that all organisms have
equal fitness.