  The integration of empirical data in computational frameworks to model the
spread of infectious diseases poses challenges that are becoming pressing with
the increasing availability of high-resolution information on human mobility
and contacts. This deluge of data has the potential to revolutionize the
computational efforts aimed at simulating scenarios and designing containment
strategies. However, the integration of detailed data sources yields models
that are less transparent and general. Hence, given a specific disease model,
it is crucial to assess which representations of the raw data strike the best
balance between simplicity and detail. We consider high-resolution data on the
face-to-face interactions of individuals in a hospital ward, obtained by using
wearable proximity sensors. We simulate the spread of a disease in this
community by using an SEIR model on top of different mathematical
representations of the contact patterns. We show that a contact matrix that
only contains average contact durations fails to reproduce the size of the
epidemic obtained with the high-resolution contact data and also to identify
the most at-risk classes. We introduce a contact matrix of probability
distributions that takes into account the heterogeneity of contact durations
between (and within) classes of individuals, and we show that this
representation yields a good approximation of the epidemic spreading properties
obtained by using the high-resolution data. Our results mark a step towards the
definition of synopses of high-resolution dynamic contact networks, providing a
compact representation of contact patterns that can correctly inform
computational models designed to discover risk groups and evaluate containment
policies. We show that this novel kind of representation can preserve in
simulation quantitative features of the epidemics that are crucial for their
study and management.


  Efforts to catalogue the structure of metabolic networks have generated
highly detailed, genome-scale atlases of biochemical reactions in the cell.
Unfortunately, these atlases fall short of capturing the kinetic details of
metabolic reactions, instead offering only \textit{topological} information
from which to make predictions. As a result, studies frequently consider the
extent to which the topological structure of a metabolic network determines its
dynamic behavior, irrespective of kinetic details. Here, we study a class of
metabolic networks known as non-autocatalytic metabolic cycles, and
analytically prove an open conjecture regarding the stability of their
steady-states. Importantly, our results are invariant to the choice of kinetic
parameters, rate laws, equilibrium fluxes, and metabolite concentrations.
Unexpectedly, our proof exposes an elementary but apparently open problem of
locating the roots of a sum of two polynomials S = P+Q, when the roots of the
summand polynomials P and Q are known. We derive two new results named the
Stubborn Roots Theorems, which provide sufficient conditions under which the
roots of S remain qualitatively identical to the roots of P. Our work
illustrates how complementary feedback, from classical fields like dynamical
systems to biology and vice versa, can expose fundamental and potentially
overlooked questions.


  Stochasticity in gene expression gives rise to fluctuations in protein levels
across a population of genetically identical cells. Such fluctuations can lead
to phenotypic variation in clonal populations, hence there is considerable
interest in quantifying noise in gene expression using stochastic models.
However, obtaining exact analytical results for protein distributions has been
an intractable task for all but the simplest models. Here, we invoke the
partitioning property of Poisson processes to develop a mapping that
significantly simplifies the analysis of stochastic models of gene expression.
The mapping leads to exact protein distributions using results for mRNA
distributions in models with promoter-based regulation. Using this approach, we
derive exact analytical results for steady-state and time-dependent
distributions for the basic 2-stage model of gene expression. Furthermore, we
show how the mapping leads to exact protein distributions for extensions of the
basic model that include the effects of post-transcriptional and
post-translational regulation. The approach developed in this work is widely
applicable and can contribute to a quantitative understanding of stochasticity
in gene expression and its regulation.


  Savannas are characterized by a discontinuous tree layer superimposed on a
continuous layer of grass. Identifying the mechanisms that facilitate this
tree-grass coexistence has remained a persistent challenge in ecology and is
known as the "savanna problem". In this work, we propose a model that combines
a previous savanna model (Calabrese et al., 2010), which includes competitive
interactions among trees and dispersal, with the Drossel-Schwabl forest fire
model, therefore representing fire in a spatially explicit manner. The model is
used to explore how the pattern of fire-spread, coupled with an explicit,
fire-vulnerable tree life stage, affects tree density and spatial pattern. Tree
density depends strongly on both fire frequency and tree-tree competition
although the fire frequency, which induces indirect interactions between trees
and between trees and grass, appears to be the crucial factor controlling the
tree-extinction transition in which the savanna becomes grassland. Depending on
parameters, adult trees may arrange in different regular or clumped patterns,
the later of two different types (compact or open). Cluster-size distributions
have fat tails but clean power-law behavior is only attained in specific cases.


  Multisite phosphorylation networks are encountered in many intracellular
processes like signal transduction, cell-cycle control or nuclear signal
integration. In this contribution networks describing the phosphorylation and
dephosphorylation of a protein at $n$ sites in a sequential distributive
mechanism are considered. Multistationarity (i.e.\ the existence of at least
two positive steady state solutions of the associated polynomial dynamical
system) has been analyzed and established in several contributions. It is, for
example, known that there exist values for he rate constants where
multistationarity occurs. However, nothing else is known about these rate
constants.
  Here we present a sign condition that is necessary and sufficient for
multistationarity in $n$-site sequential, distributive phosphorylation. We
express this sign condition in terms of linear systems and show that solutions
of these systems define rate constants where multistationarity is possible. We
then present, for $n\geq 2$, a collection of {\em feasible} linear systems and
hence give a new and independent proof that multistationarity is possible for
$n\geq 2$. Moreover, our results allow to explicitly obtain values for the rate
constants where multistationarity is possible. Hence we believe that, for the
first time, a systematic exploration of the region in parameter space where
multistationarity occurs has become possible.One consequence of our work is
that, for any pair of steady states, the ratio of the steady state
concentrations of kinase-substrate complexes equals that of
phosphatase-substrate complexes.


  Protein rigidity and flexibility can be analyzed accurately and efficiently
using the program FIRST. Previous studies using FIRST were designed to analyze
the rigidity and flexibility of proteins using a single static (snapshot)
structure. It is however well known that proteins can undergo spontaneous
sub-molecular unfolding and refolding, or conformational dynamics, even under
conditions that strongly favour a well-defined native structure. These (local)
unfolding events result in a large number of conformers that differ from each
other very slightly. In this context, proteins are better represented as a
thermodynamic ensemble of `native-like' structures, and not just as a single
static low-energy structure.
  Working with this notion, we introduce a novel FIRST-based approach for
predicting rigidity/flexibility of the protein ensemble by (i) averaging the
hydrogen bonding strengths from the entire ensemble and (ii) by refining the
mathematical model of hydrogen bonds. Furthermore, we combine our
FIRST-ensemble rigidity predictions with the ensemble solvent accessibility
data of the backbone amides and propose a novel computational method which uses
both rigidity and solvent accessibility for predicting hydrogen-deuterium
exchange (HDX). To validate our predictions, we report a novel site specific
HDX experiment which characterizes the native structural ensemble of
Acylphosphatase from hyperthermophile Sulfolobus solfataricus (Sso AcP).
  The sub-structural conformational dynamics that is observed by HDX data, is
closely matched with the FIRST-ensemble rigidity predictions, which could not
be attained using the traditional single `snapshot' rigidity analysis.
Moreover, the computational predictions of regions that are protected from HDX
and those that undergo exchange are in very good agreement with the
experimental HDX profile of Sso AcP.


  A perturbation framework is developed to analyze metastable behavior in
stochastic processes with random internal and external states. The process is
assumed to be under weak noise conditions, and the case where the deterministic
limit is bistable is considered. A general analytical approximation is derived
for the stationary probability density and the mean switching time between
metastable states, which includes the pre exponential factor. The results are
illustrated with a model of gene expression that displays bistable switching.
In this model, the external state represents the number of protein molecules
produced by a hypothetical gene. Once produced, a protein is eventually
degraded. The internal state represents the activated or unactivated state of
the gene; in the activated state the gene produces protein more rapidly than
the unactivated state. The gene is activated by a dimer of the protein it
produces so that the activation rate depends on the current protein level. This
is a well studied model, and several model reductions and diffusion
approximation methods are available to analyze its behavior. However, it is
unclear if these methods accurately approximate long-time metastable behavior
(i.e., mean switching time between metastable states of the bistable system).
Diffusion approximations are generally known to fail in this regard.


  Dysregulated microRNA (miRNA) expression is a well-established feature of
human cancer. However, the role of specific miRNAs in determining cancer
outcomes remains unclear. Using Level 3 expression data from the Cancer Genome
Atlas (TCGA), we identified 61 miRNAs that are associated with overall survival
in 469 ovarian cancers profiled by microarray (p<0.01). We also identified 12
miRNAs that are associated with survival when miRNAs were profiled in the same
specimens using Next Generation Sequencing (miRNA-Seq) (p<0.01). Surprisingly,
only 1 miRNA transcript is associated with ovarian cancer survival in both
datasets. Our analyses indicate that this discrepancy is due to the fact that
miRNA levels reported by the two platforms correlate poorly, even after
correcting for potential issues inherent to signal detection algorithms.
Further investigation is warranted.


  Motivation Gene transcription requires the orchestrated binding of various
proteins to the promoter of a gene. The binding times and binding order of
proteins allow to draw conclusions about the proteins' exact function in the
recruitment process. Time-resolved ChIP experiments are being used to analyze
the order of protein binding for these processes. However, these ChIP signals
do not represent the exact protein binding patterns. Results We show that for
promoter complexes that follow sequential recruitment dynamics the ChIP signal
can be understood as a convoluted signal and propose the application of
deconvolution methods to recover the protein binding patterns from experimental
ChIP time-series. We analyze the suitability of four deconvolution methods: two
non-blind deconvolution methods, Wiener deconvolution and Lucy-Richardson
deconvolution, and two blind deconvolution methods, blind Lucy-Richardson
deconvolution and binary blind deconvolution. We apply these methods to infer
the protein binding pattern from ChIP time-series for the pS2 gene.


  Although Vitamin D is believed to be involved in a large number of immune
responses our understanding of these processes at cellular level remained at
infancy. We develop and solve a coarse grained kinetic network model to
quantify the effects of variation of vitamin D on human immunity. The system of
equations accounts for known inter-relation between active and inactive vitamin
D, antigen presenting cells, effector T cells, regulatory T cells and pathogen.
Both time dependent and steady state solutions are obtained. The time dependent
solution of the system of equations reveals that the immune response is rather
strongly regulated in presence of vitamin D. We found quantitatively that lower
than optimum levels of concentration of active vitamin D correspond to weak
regulation where, once a pathogen/antigen enters the body, the nature of the
immune response would be less regulatory and hence more aggressive, or
inflammatory. The steady state solution of our model shows that vitamin D
enhances the tolerance level of immune system, thereby increasing resistance to
autoimmune diseases. Our model and accompanied numerical analyses reveal
another important aspect of immunity: While extremely low levels of vitamin D
could lead to increased risk of autoimmune responses, an overdose (toxic) level
would give rise to too large a tolerant response, leading to the increased risk
of tumors and cancerous cell growth.


  Recent work has revealed a new class of "zero-determinant" (ZD) strategies
for iterated, two-player games. ZD strategies allow a player to unilaterally
enforce a linear relationship between her score and her opponent's score, and
thus achieve an unusual degree of control over both players' long-term payoffs.
Although originally conceived in the context of classical, two-player game
theory, ZD strategies also have consequences in evolving populations of
players. Here we explore the evolutionary prospects for ZD strategies in the
Iterated Prisoner's Dilemma (IPD). Several recent studies have focused on the
evolution of "extortion strategies" - a subset of zero-determinant strategies -
and found them to be unsuccessful in populations. Nevertheless, we identify a
different subset of ZD strategies, called "generous ZD strategies", that
forgive defecting opponents, but nonetheless dominate in evolving populations.
For all but the smallest population sizes, generous ZD strategies are not only
robust to being replaced by other strategies, but they also can selectively
replace any non-cooperative ZD strategy. Generous strategies can be generalized
beyond the space of ZD strategies, and they remain robust to invasion. When
evolution occurs on the full set of all IPD strategies, selection
disproportionately favors these generous strategies. In some regimes, generous
strategies outperform even the most successful of the well-known Iterated
Prisoner's Dilemma strategies, including win-stay-lose-shift.


  Energy landscape theory describes how a full-length protein can attain its
native fold by sampling only a tiny fraction of all possible structures.
Although protein folding is now understood to be concomitant with synthesis on
the ribosome, there have been few attempts to modify energy landscape theory by
accounting for cotranslational folding. Here we provide a model for
cotranslational folding that leads to a natural definition of a nested energy
landscape. By applying concepts drawn from submanifold differential geometry,
the physics of protein folding on the ribosome can be explored in a
quantitative manner and conditions on the nested energy landscapes for a good
cotranslational folder are derived.


  The study reported here concerns a phenomenon, discovered and extensively
investigated by Sorin Comorosan, wherein enzyme initial reaction rates are
enhanced as a consequence of incorporation of solutions derived from previously
irradiated crystalline material into the reaction medium. Effective irradiation
times conform to a sharply oscillatory pattern. In most reports, the irradiated
crystalline material has been the substrate for the enzyme reaction to be
studied, but there have been exceptions. The experiments presented here serve
to confirm and extend this latter aspect of the phenomenon. It is found that
the initial reaction rates for the lactic acid dehydrogenase (LDH) conversion
of pyruvate to lactate can be stimulated by irradiation of crystalline deposits
of sodium chloride, sodium bromide, potassium chloride and diatomaceous earth.
Similarly, stimulation of the LDH conversion of lactate to pyruvate is
demonstrated for irradiated sodium chloride. There appears to be no required
chemical feature of the irradiated material other than crystalline state.


  In a recent work we proposed a coarse-grained methodology for studying the
response of peptides when simulated at different values of pH; in this work we
extend the methodology to analyze the pH-dependent behavior of coiled coils.
This protein structure presents a remarkable chain stiffness andis formed by
two or more long helical peptides that are interacting like the strands of a
rope. Chain length and rigidity are the key aspects needed to extend previous
peptide interaction potentials to this particular case; however the original
model is naturally recovered when the length or the ridigity of the simulated
chain are reduced. We apply the model and discuss results for two cases: (a)
the folding/unfolding transition of a generic coiled coil as a function of pH;
(b) behavior of a specific sequence as a function of the acidity conditions. In
this latter case results are compared with experimental data from the
literature in order to comment about the consistency of our approach.


  We propose a novel path sampling method based on the Onsager-Machlup (OM)
action by generalizing the multiscale enhanced sampling (MSES) technique
suggested by Moritsugu and coworkers (J. Chem. Phys. 133, 224105 (2010)). The
basic idea of this method is that the system we want to study (for example,
some molecular system described by molecular mechanics) is coupled to a
coarse-grained (CG) system, which can move more quickly and computed more
efficiently than the original system. We simulate this combined system
(original + CG system) using (underdamped) Langevin dynamics where different
heat baths are coupled to the two systems. When the coupling is strong enough,
the original system is guided by the CG system, and able to sample the
configuration and path space more efficiency. We need to correct the bias
caused by the coupling, however, by employing the Hamiltonian replica exchange
where we prepare many path replica with different coupling strengths. As a
result, an unbiased path ensemble for the original system can be found in the
weakest coupling path ensemble. This strategy is easily implemented because a
weight for a path calculated by the OM action is formally the same as the
Boltzmann weight if we properly define the path "Hamiltonian". We apply this
method to a model polymer with Asakura-Oosawa interaction, and compare the
results with the conventional transition path sampling method.


  In this paper, three multiscale methods for coupling of mesoscopic
(compartment-based) and microscopic (molecular-based) stochastic
reaction-diffusion simulations are investigated. Two of the three methods that
will be discussed in detail have been previously reported in the literature;
the two-regime method (TRM) and the compartment-placement method (CPM). The
third method that is introduced and analysed in this paper is the ghost cell
method (GCM). Presented is a comparison of sources of error. The convergent
properties of this error are studied as the time step $\Delta t$ (for updating
the molecular-based part of the model) approaches zero. It is found that the
error behaviour depends on another fundamental computational parameter $h$, the
compartment size in the mesoscopic part of the model. Two important limiting
cases, which appear in applications, are considered: (i) \Delta t approaches 0
and h is fixed; and (ii) \Delta t approaches 0 and h approaches 0 such that
\Delta t/h^2 is fixed. The error for previously developed approaches (the TRM
and CPM) converges to zero only in the limiting case (ii), but not in case (i).
It is shown that the error of the GCM converges in the limiting case (i). Thus
the GCM is superior to previous coupling techniques if the mesoscopic
description is much coarser than the microscopic part of the model.


  In this paper we present a sampling framework for RNA structures of fixed
topological genus. We introduce a novel, linear time, uniform sampling
algorithm for RNA structures of fixed topological genus $g$, for arbitrary
$g>0$. Furthermore we develop a linear time sampling algorithm for RNA
structures of fixed topological genus $g$ that are weighted by a simplified,
loop-based energy functional. For this process the partition function of the
energy functional has to be computed once, which has $O(n^2)$ time complexity.


  An important task of human genetics studies is to accurately predict disease
risks in individuals based on genetic markers, which allows for identifying
individuals at high disease risks, and facilitating their disease treatment and
prevention. Although hundreds of genome-wide association studies (GWAS) have
been conducted on many complex human traits in recent years, there has been
only limited success in translating these GWAS data into clinically useful risk
prediction models. The predictive capability of GWAS data is largely
bottlenecked by the available training sample size due to the presence of
numerous variants carrying only small to modest effects. Recent studies have
shown that different human traits may share common genetic bases. Therefore, an
attractive strategy to increase the training sample size and hence improve the
prediction accuracy is to integrate data of genetically correlated phenotypes.
Yet the utility of genetic correlation in risk prediction has not been explored
in the literature. In this paper, we analyzed GWAS data for bipolar and related
disorders (BARD) and schizophrenia (SZ) with a bivariate ridge regression
method, and found that jointly predicting the two phenotypes could
substantially increase prediction accuracy as measured by the AUC (area under
the receiver operating characteristic curve). We also found similar prediction
accuracy improvements when we jointly analyzed GWAS data for Crohn's disease
(CD) and ulcerative colitis (UC). The empirical observations were substantiated
through our comprehensive simulation studies, suggesting that a gain in
prediction accuracy can be obtained by combining phenotypes with relatively
high genetic correlations. Through both real data and simulation studies, we
demonstrated pleiotropy as a valuable asset that opens up a new opportunity to
improve genetic risk prediction in the future.


  ATP-hydrolysis is the basic energy source of many physiological processes,
but there is a lack of knowledge regarding its biological role other than
energy transfer and thermogenesis. Not all the energy released by
ATP-hydrolysis could be used in powering target biological processes and
functions. Partial energy dissipates into water. By validating the impact of
this dissipated energy, we found that energy released by ATP hydrolysis could
induce notable regulation of biomolecule's properties 100 nanometers away.
Namely ATP hydrolyzation is recycled in remote biochemical property modulation.


  Prion diseases are invariably fatal and highly infectious neurodegenerative
diseases that affect a wide variety of mammalian species such as sheep, goats,
mice, humans, chimpanzees, hamsters, cattle, elks, deer, minks, cats, chicken,
pigs, turtles, etc. These neurodegenerative diseases are caused by the
conversion from a soluble normal cellular protein into insoluble abnormally
folded infectious prions and the conversion is believed to involve
conformational change from a predominantly alpha-helical protein to one rich in
beta-sheet structure. Such conformational changes may be amenable to study by
molecular dynamics (MD) techniques. For rabbits, classical studies show they
have a low susceptibility to be infected, but in 2012 it was reported that
rabbit prion can be generated (though not directly) and the rabbit prion is
infectious and transmissible (Proceedings of the National Academy of Sciences
USA 109(13): 5080-5). This paper studies the NMR and X-ray molecular structures
of rabbit prion protein wild-type and mutants by MD techniques, in order to
understand the specific mechanism of rabbit prion protein and rabbit prions.


  Genes linked to X or Z chromosomes, which are hemizygous in the heterogametic
sex, are predicted to evolve at different rates than those on autosomes. This
faster-X effect can arise either as a consequence of hemizygosity, which leads
to more efficient selection for recessive beneficial mutations in the
heterogametic sex, or as a consequence of reduced effective population size of
the hemizygous chromosome, which leads to increased fixation of weakly
deleterious mutations due to random genetic drift. Empirical results to date
have suggested that, while the overall pattern across taxa is complicated, in
general systems with male-heterogamy show a faster-X effect primarily
attributable to more efficient selection while the only female-heterogamy taxon
studied to date (birds) shows a faster-Z effect primarily attributable to
increased drift. In order to test the generality of the faster-Z pattern seen
in birds, we sequenced the genome of the Lepidopteran insect Bombyx huttoni, a
close outgroup of the domesticated silkmoth Bombyx mori. We show that silkmoths
experience faster-Z evolution, but unlike in birds, the faster-Z effect appears
to be attributable to more efficient positive selection in females. These
results suggest that female-heterogamy alone is unlikely to be sufficient to
explain the reduced efficacy of selection on the bird Z chromosome. Instead, it
is likely that a combination of patterns of dosage compensation and overall
effective population size, among other factors, influence patterns of faster-Z
evolution.


  Drug resistance is a serious public health problem that threatens to thwart
our ability to treat many infectious diseases. Repeatedly, the introduction of
new drugs has been followed by the evolution of resistance. In principle there
are two ways to address this problem: (i) enhancing drug development, and (ii)
slowing drug resistance. We present data and a modeling approach based on
queueing theory that explores how interventions aimed at these two facets
affect the ability of the entire drug supply system to provide service.
Analytical and simulation-based results show that, all else equal, slowing the
evolution of drug resistance is more effective at ensuring an adequate supply
of effective drugs than is enhancing the rate at which new drugs are developed.
This lends support to the idea that evolution management is not only a
significant component of the solution to the problem of drug resistance, but
may in fact be the most important component.


  Electrical neural signalling typically takes place at the time-scale of
milliseconds, and is typically modeled using the cable equation. This is a good
approximation for processes when ionic concentrations vary little during the
time course of a simulation. During periods of intense neural signalling,
however, the local extracellular K+ concentration may increase by several
millimolars. Clearance of excess K+ likely depends partly on diffusion in the
extracellular space, partly on local uptake by- and intracellular transport
within astrocytes. This process takes place at the time scale of seconds, and
can not be modeled accurately without accounting for the spatiotemporal
variations in ion concentrations. The work presented here consists of two main
parts: First, we developed a general electrodiffusive formalism for modeling
ion concentration dynamics in a one-dimensional geometry, including both an
intra- and extracellular domain. The formalism was based on the Nernst-Planck
equations. It ensures (i) consistency between the membrane potential and ion
concentrations, (ii) global particle/charge conservation, and (iii) accounts
for diffusion and concentration dependent variations in resistivities. Second,
we applied the formalism to model how astrocytes exchange ions with the ECS,
and identified the key astrocytic mechanisms involved in K+ removal from high
concentration regions. We found that a local increase in extracellular
K\textsuperscript{+} evoked a local depolarization of the astrocyte membrane,
which at the same time (i) increased the local astrocytic uptake of
K\textsuperscript{+}, (ii) suppressed extracellular transport of K+, (iii)
increased transport of K+ within astrocytes, and (iv) facilitated astrocytic
relase of K+ in extracellular low concentration regions. In summary, these
mechanisms seem optimal for shielding the extracellular space from excess K+.


  Calcium is an ubiquitous second messenger that triggers a plethora of key
physiological responses. The events are initiated in micro- or nano-sized
compartments and determined by the complex interactions with calcium-binding
proteins and mechanisms of calcium clearance. Local calcium increases in the
vicinity of single channels represent an essentially non-linear
reaction-diffusion problem that have been analysed previously using various
linearized approximations. I revisited the problem of stationary patterns that
can be generated by the point calcium source in the presence of buffer and
obtained new explicit solutions. Main results of the analysis of the calcium
buffering are supplemented with pertinent derivations and discussion of
respective mathematical problems in Appendices. I show that for small calcium
influx the calcium gradients around established around channel lumen have
quasi-exponential form. For bigger fluxes, when the buffer is saturated, the
model predicts periodic patterns. The transition between the two regimes depend
on the capacity of buffer and its mobility. Theoretical predictions were
examined using a model one-dimensional system. For sufficiently big fluxes the
oscillatory calcium patterns were observed. Theoretical and experimental
results are discussed in terms of their possible physiological implications.


  Birth-and-death models are now a common mathematical tool to describe
branching patterns observed in real-world phylogenetic trees. Liggett and
Schinazi (2009) is one such example. The authors propose a simple
birth-and-death model that is compatible with phylogenetic trees of both
influenza and HIV, depending on the birth rate parameter. An interesting
special case of this model is the critical case where the birth rate equals the
death rate. This is a non-trivial situation and to study its asymptotic
behaviour we employed the Laplace transform. With this we correct the proof of
Liggett and Schinazi (2009) in the critical case.


  Single-molecule pulling experiments on unstructured proteins linked to
neurodegenerative diseases have measured rupture forces comparable to those for
stable folded proteins. To investigate the structural mechanisms of this
unexpected force resistance, we perform pulling simulations of the amyloid
{\beta}-peptide (A{\beta}) and {\alpha}-synuclein ({\alpha}S), starting from
simulated conformational ensembles for the free monomers. For both proteins,
the simulations yield a set of rupture events that agree well with the
experimental data. By analyzing the conformations right before rupture in each
event, we find that the mechanically resistant structures share a common
architecture, with similarities to the folds adopted by A{\beta} and {\alpha}S
in amyloid fibrils. The disease-linked Arctic mutation of A{\beta} is found to
increase the occurrence of highly force-resistant structures. Our study
suggests that the high rupture forces observed in A{\beta} and {\alpha}S
pulling experiments are caused by structures that might have a key role in
amyloid formation.


  Systemic approaches to the study of a biological cell or tissue rely
increasingly on the use of context-specific metabolic network models. The
reconstruction of such a model from high-throughput data can routinely involve
large numbers of tests under different conditions and extensive parameter
tuning, which calls for fast algorithms. We present FASTCORE, a generic
algorithm for reconstructing context-specific metabolic network models from
global genome-wide metabolic network models such as Recon X. FASTCORE takes as
input a core set of reactions that are known to be active in the context of
interest (e.g., cell or tissue), and it searches for a flux consistent
subnetwork of the global network that contains all reactions from the core set
and a minimal set of additional reactions. Our key observation is that a
minimal consistent reconstruction can be defined via a set of sparse modes of
the global network, and FASTCORE iteratively computes such a set via a series
of linear programs. Experiments on liver data demonstrate speedups of several
orders of magnitude, and significantly more compact reconstructions, over a
chief rival method. Given its simplicity and its excellent performance,
FASTCORE can form the backbone of many future metabolic network reconstruction
algorithms.


  The present study aims at integrating drug-releasing materials with
signal-processing biocomputing systems. Enzymes alanine transaminase (ALT) and
aspartate transaminase (AST)---biomarkers for liver injury---were logically
processed by a biocatalytic cascade realizing Boolean AND gate. Citrate
produced in the system was used to trigger a drug-mimicking release from
alginate microspheres. In order to differentiate low vs. high concentration
signals, the microspheres were coated with a protective shell composed of
layer-by-layer adsorbed poly(L-lysine) and alginate. The alginate core of the
microspheres was prepared from (Fe3+)-cross-linked alginate loaded with
rhodamine 6G dye mimicking a drug. Dye release from the core occurred only when
both biomarkers, ALT and AST, appeared at their high pathophysiological
concentrations jointly indicative of liver injury. The signal-triggered
response was studied at the level of a single microsphere, yielding information
on the dye release kinetics.


  The conformational vibrations of Z-DNA with counterions are studied in
framework of phenomenological model developed. The structure of left-handed
double helix with counterions neutralizing the negatively charged phosphate
groups of DNA is considered as the ion-phosphate lattice. The frequencies and
Raman intensities for the modes of Z-DNA with Na+ and Mg2+ ions are calculated,
and the low-frequency Raman spectra are built. At the spectra range about the
frequency 150 cm-1 new mode of ion-phosphate vibrations is found, which
characterizes the vibrations of Mg2+ counterions. The results of our
calculations show that the intensities of Z-DNA modes are sensitive to the
concentration of magnesium counterions. The obtained results describe well the
experimental Raman spectra of Z-DNA.


  The global dynamics of gene regulatory networks are known to show robustness
to perturbations in the form of intrinsic and extrinsic noise, as well as
mutations of individual genes. One molecular mechanism underlying this
robustness has been identified as the action of so-called microRNAs that
operate via feedforward loops. We present results of a computational study,
using the modeling framework of stochastic Boolean networks, which explores the
role that such network motifs play in stabilizing global dynamics. The paper
introduces a new measure for the stability of stochastic networks. The results
show that certain types of feedforward loops do indeed buffer the network
against stochastic effects.


  Several modifications of the famous mathematical Game of Life are introduced
by making Game of Life rules stochastic and mutual influence of cells in their
8-neighborhood on a rectangular lattice spatially non-uniform. Results are
reported of experimental investigation of evolutionary dynamics of the
introduced models. A number of new phenomena in the evolutionary dynamics of
the models and collective behavior of patterns they generate are revealed,
described and illustrated: formation of maze-like patterns as fixed points of
the models, "self-controlled growth", "eternal life" in a bounded space and
"coherent shrinkage".


  Identification of every single genome present in a microbial sample is an
important and challenging task with crucial applications. It is challenging
because there are typically millions of cells in a microbial sample, the vast
majority of which elude cultivation. The most accurate method to date is
exhaustive single cell sequencing using multiple displacement amplification,
which is simply intractable for a large number of cells. However, there is hope
for breaking this barrier as the number of different cell types with distinct
genome sequences is usually much smaller than the number of cells.
  Here, we present a novel divide and conquer method to sequence and de novo
assemble all distinct genomes present in a microbial sample with a sequencing
cost and computational complexity proportional to the number of genome types,
rather than the number of cells. The method is implemented in a tool called
Squeezambler. We evaluated Squeezambler on simulated data. The proposed divide
and conquer method successfully reduces the cost of sequencing in comparison
with the naive exhaustive approach.
  Availability: Squeezambler and datasets are available under
http://compbio.cs.wayne.edu/software/squeezambler/.


  In the past six decades, lifespan inequality has varied greatly within and
among countries even while life expectancy has continued to increase. How and
why does mortality change generate this diversity? We derive a precise link
between changes in age-specific mortality and lifespan inequality, measured as
the variance of age at death. Key to this relationship is a young-old threshold
age, below and above which mortality decline respectively decreases and
increases lifespan inequality. First, we show that shifts in the threshold's
location modified the correlation between changes in life expectancy and
lifespan inequality over the last two centuries. Second, we analyze the post
Second World War trajectories of lifespan inequality in a set of developed
countries, Japan, Canada and the United States (US), where thresholds centered
on retirement age. Our method reveals how divergence in the age-pattern of
mortality change drives international divergence in lifespan inequality. Most
strikingly, early in the 1980s, mortality increases in young US males led
lifespan inequality to remain high in the US, while in Canada the decline of
inequality continued. In general, our wider international comparisons show that
mortality change varied most at young working ages after the Second World War,
particularly for males. We conclude that if mortality continues to stagnate at
young ages, yet declines steadily at old ages, increases in lifespan inequality
will become a common feature of future demographic change.


  Transitions to absorbing states are of fundamental importance in
non-equilibrium physics as well as ecology. In ecology, absorbing states
correspond to the extinction of species. We here study the spatial population
dynamics of three cyclically interacting species. The interaction scheme
comprises both direct competition between species as in the cyclic
Lotka-Volterra model, and separated selection and reproduction processes as in
the May-Leonard model. We show that the dynamic processes leading to the
transient maintenance of biodiversity are closely linked to attractors of the
nonlinear dynamics for the overall species' concentrations. The characteristics
of these global attractors change qualitatively at certain threshold values of
the mobility, and depend on the relative strength of the different types of
competition between species. They give information about the scaling of
extinction times with the system size and thereby the stability of
biodiversity. We define an effective free energy as the negative logarithm of
the probability to find the system in a specific global state before reaching
one of the absorbing states. The global attractors then correspond to minima of
this effective energy landscape and determine the most probable values for the
species' global concentrations. As in equilibrium thermodynamics, qualitative
changes in the effective free energy landscape indicate and characterize the
underlying non-equilibrium phase transitions. We provide the complete phase
diagrams for the population dynamics, and give a comprehensive analysis of the
spatio-temporal dynamics and routes to extinction in the respective phases.


  Discrete dynamical systems defined by the iteration of a polynomial map of
the unit simplex to itself appear in the context of population genetic systems
evolving under mutation, recombination and weak selection. Although exceptional
progress has been made in finding particular solutions to these systems, our
knowledge of the general properties of the space of all possible dynamical
systems of this kind is still limited. We prove that the space of
bounded-degree polynomial maps of the unit simplex to itself is a compact and
convex subset of a Euclidean space. We provide an explicit characterization of
such a space and of its boundary. A special class of maps in the boundary, the
folding maps, which generalize the logistic map for any dimension and degree
are defined and constructed. Finally, we use numerical methods to study the
ergodic and mixing properties of maps in the neighborhood of several of these
folding maps.


  Motivation:
  Rapid technological progress in DNA sequencing has stimulated interest in
compressing the vast datasets that are now routinely produced. Relatively
little attention has been paid to compressing the quality scores that are
assigned to each sequence, even though these scores may be harder to compress
than the sequences themselves. By aggregating a set of reads into a compressed
index, we find that the majority of bases can be predicted from the sequence of
bases that are adjacent to them and hence are likely to be less informative for
variant calling or other applications. The quality scores for such bases are
aggressively compressed, leaving a relatively small number at full resolution.
Since our approach relies directly on redundancy present in the reads, it does
not need a reference sequence and is therefore applicable to data from
metagenomics and de novo experiments as well as to resequencing data.
  Results:
  We show that a conservative smoothing strategy affecting 75% of the quality
scores above Q2 leads to an overall quality score compression of 1 bit per
value with a negligible effect on variant calling. A compression of 0.68 bit
per quality value is achieved using a more aggressive smoothing strategy, again
with a very small effect on variant calling.
  Availability:
  Code to construct the BWT and LCP-array on large genomic data sets is part of
the BEETL library, available as a github respository at
http://git@github.com:BEETL/BEETL.git .


  The advent of "next-generation" DNA sequencing (NGS) technologies has meant
that collections of hundreds of millions of DNA sequences are now commonplace
in bioinformatics. Knowing the longest common prefix array (LCP) of such a
collection would facilitate the rapid computation of maximal exact matches,
shortest unique substrings and shortest absent words. CPU-efficient algorithms
for computing the LCP of a string have been described in the literature, but
require the presence in RAM of large data structures. This prevents such
methods from being feasible for NGS datasets.
  In this paper we propose the first lightweight method that simultaneously
computes, via sequential scans, the LCP and BWT of very large collections of
sequences. Computational results on collections as large as 800 million
100-mers demonstrate that our algorithm scales to the vast sequence collections
encountered in human whole genome sequencing experiments.


  In microbial ecology studies, the most commonly used ways of investigating
alpha (within-sample) diversity are either to apply count-only measures such as
Simpson's index to Operational Taxonomic Unit (OTU) groupings, or to use
classical phylogenetic diversity (PD), which is not abundance-weighted.
Although alpha diversity measures that use abundance information in a
phylogenetic framework do exist, but are not widely used within the microbial
ecology community. The performance of abundance-weighted phylogenetic diversity
measures compared to classical discrete measures has not been explored, and the
behavior of these measures under rarefaction (sub-sampling) is not yet clear.
In this paper we compare the ability of various alpha diversity measures to
distinguish between different community states in the human microbiome for
three different data sets. We also present and compare a novel one-parameter
family of alpha diversity measures, BWPD_\theta, that interpolates between
classical phylogenetic diversity (PD) and an abundance-weighted extension of
PD. Additionally, we examine the sensitivity of these phylogenetic diversity
measures to sampling, via computational experiments and by deriving a closed
form solution for the expectation of phylogenetic quadratic entropy under
re-sampling. In all three of the datasets considered, an abundance-weighted
measure is the best differentiator between community states. OTU-based
measures, on the other hand, are less effective in distinguishing community
types. In addition, abundance-weighted phylogenetic diversity measures are less
sensitive to differing sampling intensity than their unweighted counterparts.
Based on these results we encourage the use of abundance-weighted phylogenetic
diversity measures, especially for cases such as microbial ecology where
species delimitation is difficult.


  Facing the threats of infectious diseases, we take various actions to protect
ourselves, but few studies considered an evolving system with competing
strategies. In view of that, we propose an evolutionary epidemic model coupled
with human behaviors, where individuals have three strategies: vaccination,
self-protection and laissez faire, and could adjust their strategies according
to their neighbors' strategies and payoffs at the beginning of each new season
of epidemic spreading. We found a counter-intuitive phenomenon analogous to the
well-known \emph{Braess's Paradox}, namely a better condition may lead to worse
performance. Specifically speaking, increasing the successful rate of
self-protection does not necessarily reduce the epidemic size or improve the
system payoff. This phenomenon is insensitive to the network topologies, and
can be well explained by a mean-field approximation. Our study demonstrates an
important fact that a better condition for individuals may yield a worse
outcome for the society.


  The characterization of plasticity, robustness, and evolvability, an
important issue in biology, is studied in terms of phenotypic fluctuations. By
numerically evolving gene regulatory networks, the proportionality between the
phenotypic variances of epigenetic and genetic origins is confirmed. The former
is given by the variance of the phenotypic fluctuation due to noise in the
developmental process; and the latter, by the variance of the phenotypic
fluctuation due to genetic mutation. The relationship suggests a link between
robustness to noise and to mutation, since robustness can be defined by the
sharpness of the distribution of the phenotype. Next, the proportionality
between the variances is demonstrated to also hold over expressions of
different genes (phenotypic traits) when the system acquires robustness through
the evolution. Then, evolution under environmental variation is numerically
investigated and it is found that both the adaptability to a novel environment
and the robustness are made compatible when a certain degree of phenotypic
fluctuations exists due to noise. The highest adaptability is achieved at a
certain noise level at which the gene expression dynamics are near the critical
state to lose the robustness. Based on our results, we revisit Waddington's
canalization and genetic assimilation with regard to the two types of
phenotypic fluctuations.


  Twitter is a micro-blogging social media platform for short messages that can
have a long-term impact on how scientists create and publish ideas. We
investigate the usefulness of twitter in the development and distribution of
scientific knowledge. At the start of the life cycle of a scientific
publication, twitter provides a large virtual department of colleagues that can
help to rapidly generate, share and refine new ideas. As ideas become
manuscripts, twitter can be used as an informal arena for the pre-review of
works in progress. Finally, tweeting published findings can communicate
research to a broad audience of other researchers, decision makers, journalists
and the general public that can amplify the scientific and social impact of
publications. However, there are limitations, largely surrounding issues of
intellectual property and ownership, inclusiveness and misrepresentations of
science sound bites. Nevertheless, we believe twitter is a useful social media
tool that can provide a valuable contribution to scientific publishing in the
21st century.


  The puzzle associated with the cost of sex, an old problem of evolutionary
biology, is discussed here from the point of view of nonequilibrium statistical
mechanics. The results suggest, in a simplified model, that the prevalence of
sexual species in nature can be a natural and necessary consequence of the
discrete character of the nonlinear interactions between couples and their
pathogens/parasites. Mapped into a field theory, the stochastic processes
performed by the species are described by continuous fields in space and time.
The way that the model's parameters scale with subsequent iterations of the
renormalization group gives us information about the stationary emergent
properties of the complex interacting systems modeled. We see that the
combination of one aspect of the Red Queen theory with the stochastic processes
theory, including spatiotemporal interactions, provides interesting insights
into this old Darwinian dilemma.


  We present efficient algorithms for computing a maximum agreement forest
(MAF) of a pair of multifurcating (nonbinary) rooted trees. Our algorithms
match the running times of the currently best algorithms for the binary case.
The size of an MAF corresponds to the subtree prune-and-regraft (SPR) distance
of the two trees and is intimately connected to their hybridization number.
These distance measures are essential tools for understanding reticulate
evolution, such as lateral gene transfer, recombination, and hybridization.
Multifurcating trees arise naturally as a result of statistical uncertainty in
current tree construction methods.


  A tumor can be thought of as an ecosystem, which critically means that we
cannot just consider it as a collection of mutated cells but more as a complex
system of many interacting cellular and microenvironmental elements. At its
simplest, a growing tumor with increased proliferation capacity must compete
for space as a limited resource. Hypercellularity leads to a contact-inhibited
core with a competitive proliferating rim. Evolution and selection occurs, and
an individual cell's capacity to survive and propagate is determined by its
combination of traits and interaction with the environment. With heterogeneity
in phenotypes, the clone that will dominate is not always obvious as there are
both local interactions and global pressures. Several combinations of
phenotypes can coexist, changing the fitness of the whole.
  To understand some aspects of heterogeneity in a growing tumor we build an
off-lattice agent based model consisting of individual cells with assigned
trait values for proliferation and migration rates. We represent heterogeneity
in these traits with frequency distributions and combinations of traits with
density maps. How the distributions change over time is dependent on how traits
are passed on to progeny cells, which is our main inquiry. We bypass the
translation of genetics to behavior by focussing on the functional end result
of inheritance of the phenotype combined with the environmental influence of
limited space.


  We discuss how spreading processes on temporal networks are impacted by the
shape of their inter-event time distributions. Through simple mathematical
arguments and toy examples, we find that the key factor is the ordering in
which events take place, a property that tends to be affected by the bulk of
the distributions and not only by their tail, as usually considered in the
literature. We show that a detailed modeling of the temporal patterns observed
in complex networks can change dramatically the properties of a spreading
process, such as the ergodicity of a random walk process or the persistence of
an epidemic.


  Circadian oscillation provides selection advantages through synchronization
to the daylight cycle. However, a reliable clock must be designed through two
conflicting properties: entrainability to synchronize internal time with
periodic stimuli such as sunlight, and regularity to oscillate with a precise
period. These two aspects do not easily coexist because better entrainability
favors higher sensitivity, which may sacrifice the regularity. To investigate
conditions for satisfying the two properties, we analytically calculated the
optimal phase-response curve with a variational method. Our result indicates an
existence of a dead zone, i.e., a time period during which input stimuli
neither advance nor delay the clock. A dead zone appears only when input
stimuli obey the time course of actual solar radiation but a simple sine curve
cannot yield a dead zone. Our calculation demonstrates that every circadian
clock with a dead zone is optimally adapted to the daylight cycle.


  Many biological networks have to filter out useful information from a vast
excess of spurious interactions. We use computational evolution to predict
design features of networks processing ligand categorization. The important
problem of early immune response is considered as a case-study. Rounds of
evolution with different constraints uncover elaborations of the same network
motif we name adaptive sorting. Corresponding network substructures can be
identified in current models of immune recognition. Our work draws a deep
analogy between immune recognition and biochemical adaptation.


  We propose a new algorithm to do posterior sampling of Kingman's coalescent,
based upon the Particle Markov Chain Monte Carlo methodology. Specifically, the
algorithm is an instantiation of the Particle Gibbs Sampling method, which
alternately samples coalescent times conditioned on coalescent tree structures,
and tree structures conditioned on coalescent times via the conditional
Sequential Monte Carlo procedure. We implement our algorithm as a C++ package,
and demonstrate its utility via a parameter estimation task in population
genetics on both single- and multiple-locus data. The experiment results show
that the proposed algorithm performs comparable to or better than several
well-developed methods.


  In the last decade, rapid advances have been made in the field of
micro-rheology of cells and tissues. Given the complexity of living systems,
there is a need for the development of new types of nano- and micron-sized
probes, and in particular of probes with controlled interactions with the
surrounding medium. In the present paper, we evaluate the use of micron-sized
wires as potential probes of the mechanical properties of cells. The wire-based
micro-rheology technique is applied to living cells such as murine fibroblasts
and canine kidney epithelial cells. The mean-squared angular displacement
(MSAD) of wires associated to their rotational dynamics is obtained as a
function of the time using optical microscopy and image processing. It reveals
a Brownian-like diffusive regime where the MSA scale linearly with time and as
the inverse of the cube of the wire length. This scaling suggests that an
effective viscosity of the intracellular medium can be determined, and that in
the range 1 - 10 micrometers it does not depend on the length scale over which
it is measured.


  Maximum likelihood is a popular technique for isoform reconstruction. Here,
we show that isoform reconstruction using short RNA-Seq reads by maximum
likelihood is NP-hard.


  Many protein systems fold in a two-state manner. Random models, however,
rarely display two-state kinetics and thus such behavior should not be accepted
as a default. To date, many theories for the prevalence of two-state kinetics
have been presented, but none sufficiently explain the breadth of experimental
observations. A model, making a minimum of assumptions, is introduced that
suggests two-state behavior is likely for any system with an overwhelmingly
populated native state. We show two-state folding is emergent and strengthened
by increasing the occupancy population of the native state. Further, the model
exhibits a hub-like behavior, with slow interconversions between unfolded
states. Despite this, the unfolded state equilibrates quickly relative to the
folding time. This apparent paradox is readily understood through this model.
Finally, our results compare favorable with experimental measurements of
protein folding rates as a function of chain length and Keq, and provide new
insight into these results.


  Quantifying interaction strength between species is of interest in food web
studies for understanding population dynamics. Theory has run ahead of
experiment in solving equations describing ecological systems (the
Lotka-Volterra equations, for example). Modeling approaches depend upon
quantitative knowledge of interaction strengths among species, yet there are
few methods available for estimating strengths of interactions between species,
especially in the field. Competition and predation in fishes is traditionally
studied with a variety of methods, most requiring extensive sampling to
determine stock densities or extensive stomach content or isotope analysis.
Average relative weight of fish species in an ecosystem can usually be
accurately determined with much smaller sample sizes than required for accurate
determination of stock densities. Correlations of mean annual relative weights
may be used to test hypotheses and establish interaction strengths regarding
competition and predation among fishes. The trends are sufficiently suggestive
to use a strong positive correlation as supporting evidence of competition and
possibly indicative of the magnitude of competition if other evidence such as
analysis of stomach contents is also present. Likewise, a strong negative
correlation might be interpreted as relevant when forming a hypothesis of
predation and as supporting evidence of the magnitude of predation if other
evidence is also present. In systems where regular weight and length survey
data are available, this method can suggest the likely strength of competitive
and predatory relationships and may be more cost effective than stomach content
or isotope analysis.


  RNAs are negatively charged molecules residing in macromolecular crowding
cellular environments. Macromolecular confinement can influence the ion effects
in RNA folding. In this work, using the recently developed tightly bound ion
model for ion fluctuation and correlation, we investigate the confinement
effect on the ion-mediated RNA structural collapse for a simple model system.
We found that, for both Na$^+$ and Mg$^{2+}$, ion efficiencies in mediating
structural collapse/folding are significantly enhanced by the structural
confinement. Such an enhancement in the ion efficiency is attributed to the
decreased electrostatic free energy difference between the compact conformation
ensemble and the (restricted) extended conformation ensemble due to the spatial
restriction.


  A protein contact map is a binary symmetric adjacency matrix capturing the
distance relationship between atoms of a protein. Each cell (i, j) of a protein
contact map states whether the atoms (nodes) i and j are within some Euclidean
distance from each other. We examined the radius one Moore neighborhood
surrounding each cell (i, j) where j > (i + 2) in complete protein contact maps
by mutating them one at a time. We found that the particular configuration of a
neighborhood is generally (97%) optimal in the sense that no other
configuration could maintain or improve upon existing local and global
efficiencies of the nodes residing in a neighborhood. Local efficiency of a
node is directly related to its clustering measure. Global efficiency of a node
is inversely related to its distance to other nodes in the network. This
feature of the Moore neighborhood in complete protein contact maps may explain
how protein residue networks are able to form long-range links to reduce
average path length while maintaining a high level of clustering throughout the
process of small-world formation, and could suggest new approaches to protein
contact map prediction. Effectively, the problem of protein contact map
prediction is transformed to one of maximizing the number of optimal
neighborhoods. By comparison, Moore neighborhoods in protein contact maps with
randomized long-range links are less optimal.


  Plasmodium falciparum aquaglyceroporin (PfAQP) is a multifunctional channel
protein in the plasma membrane of the malarial parasite that causes the most
severe form of malaria infecting more than a million people a year. Finding a
novel way to inhibit PfAQP, I conducted 3+ microseconds in silico experiments
of an atomistic model of the PfAQP-membrane system and computed the
chemical-potential profiles of six permeants (erythritol, water, glycerol,
urea, ammonia, and ammonium) that can be efficiently transported across P.
falciparum's plasma membrane through PfAQP's conducting pore. The profiles show
that, with all the existent in vitro data being supportive, erythritol, a
permeant of PfAQP itself having a deep ditch in its permeation passageway,
strongly inhibits PfAQP's functions of transporting water, glycerol, urea,
ammonia, and ammonium (The IC50 is in the range of high nanomolars). This
suggests the possibility that erythritol, a sweetener generally considered
safe, may be the drug needed to kill the malarial parasite in vivo without
causing serious side effects.


  Identifying targets for personalized targeted therapy is the pathologist's
domain and a treasure. For decades, pathologists have had to learn, understand,
adopt and implement many new laboratory techniques as they arrived on the
scene. Pathologists successfully integrate the results of those tests into
final pathology reports that were, and still are, the basis of clinical
therapeutic decisions. The molecular methods are different but no more
difficult to comprehend in the era of "kit procedures". Pathologists have the
knowledge and expertise to identify particular gene mutations using the
appropriate molecular tests currently available. This review focuses on the
most important recent developments in KRAS mutation testing in metastatic
colorectal cancer (CRC), and shows that a pathologist is involved in 10 stages
of this procedure. Recent studies have shown that highly sensitive, simple,
reliable and rapid assays may significantly improve the identification of CRC
patients resistant to anti-EGFR therapy. Thus, direct sequencing does not seem
to be an optimal procedure of KRAS testing for clinical purposes. Twelve
currently available high-sensitivity diagnostic assays (with the CE-IVD mark)
for KRAS mutation testing are briefly described and compared. The suggested
pathology report content for somatic mutation tests is described. In
conclusion, evidence is presented that sending away paraffin blocks with tumor
tissue for KRAS mutation testing may not be in the best interest of patients.
Instead, an evidence-based approach indicates that KRAS mutation testing should
be performed in pathology departments, only with the use of CE-IVD/FDA-approved
KRAS tests, and with the obligatory, periodic participation in the KRAS EQA
scheme organized by the European Society of Pathology as an independent
international body.


  We study the evolution of cooperation in the spatial public goods game,
focusing on the effects that are brought about by the delayed distribution of
goods that accumulate in groups due to the continuous investments of
cooperators. We find that intermediate delays enhance network reciprocity
because of a decelerated invasion of defectors, who are unable to reap the same
high short-term benefits as they do in the absence of delayed distribution.
Long delays, however, introduce a risk because the large accumulated wealth
might fall into the wrong hands. Indeed, as soon as the curvature of a
cooperative cluster turns negative, the engulfed defectors can collect the
heritage of many generations of cooperators, and by doing so start a waning
moon pattern that nullifies the benefits of decelerated invasion. Accidental
meeting points of growing cooperative clusters may also act as triggers for the
waning moon effect, thus linking the success of cooperators with their
propensity to fail in a rather bizarre way. Our results highlight that
"investing into the future" is a good idea only if that future is sufficiently
near and not likely to be burdened by inflation.


  Random walks on multidimensional nonlinear landscapes are of interest in many
areas of science and engineering. In particular, properties of adaptive
trajectories on fitness landscapes determine population fates and thus play a
central role in evolutionary theory. The topography of fitness landscapes and
its effect on evolutionary dynamics have been extensively studied in the
literature. We will survey the current research knowledge in this field,
focusing on a recently developed systematic approach to characterizing path
lengths, mean first-passage times, and other statistics of the path ensemble.
This approach, based on general techniques from statistical physics, is
applicable to landscapes of arbitrary complexity and structure. It is
especially well-suited to quantifying the diversity of stochastic trajectories
and repeatability of evolutionary events. We demonstrate this methodology using
a biophysical model of protein evolution that describes how proteins maintain
stability while evolving new functions.


  The first realization of a biomolecular OR gate function with double-sigmoid
response (sigmoid in both inputs) is reported. Two chemical inputs activate the
enzymatic gate processes resulting in the output signal: chromogen oxidation,
which occurs when either one of the inputs or both are present (corresponding
to the OR binary function), and can be optically detected. High-quality gate
functioning in handling of sources of noise is enabled by "filtering" involving
pH control with an added buffer. The resulting gate response is sigmoid in both
inputs when proper system parameters are chosen, and the gate properties are
theoretically analyzed within a model devised to evaluate its noise-handling
properties.


  Non-traditional thermodynamics, applied to random behaviour associated with
turbulence, mixing and competition, is reviewed and analysed. Competitive
mixing represents a general framework for the study of generic properties of
competitive systems and can be used to model a wide class of non-equilibrium
phenomena ranging from turbulent premixed flames and invasion waves to complex
competitive systems. We demonstrate consistency of the general principles of
competition with thermodynamic description, review and analyse the related
entropy concepts and introduce the corresponding competitive H-theorem. A
competitive system can be characterised by a thermodynamic quantity -
competitive potential --- which determines the likely direction of evolution of
the system. Contested resources tend to move between systems from lower to
higher values of the competitive potential. There is, however, an important
difference between conventional thermodynamics and competitive thermodynamics.
While conventional thermodynamics is constrained by its zeroth law and is
fundamentally transitive, the transitivity of competitive thermodynamics
depends on the transitivity of the competition rules. Intransitivities are
common in the real world and are responsible for complex behaviour in
competitive systems.
  This work follows the ideas and methods that are originated in analysis of
turbulent combustion but reviews a much broader scope of issues linked to
mixing and competition, including thermodynamic characterisation of complex
competitive systems with self-organisation. The approach presented here is
interdisciplinary and is addressed to a general educated reader, while the
mathematical details can be found in the Appendices.


  We suggest a modification of a comb model to describe anomalous transport in
spiny dendrites. Geometry of the comb structure consisting of a one-dimensional
backbone and lateral branches makes it possible to describe anomalous
diffusion, where dynamics inside fingers corresponds to spines, while the
backbone describes diffusion along dendrites. The presented analysis
establishes that the fractional dynamics in spiny dendrites is controlled by
fractal geometry of the comb structure and fractional kinetics inside the
spines. Our results show that the transport along spiny dendrites is
subdiffusive and depends on the density of spines in agreement with recent
experiments.


  In a letter published in Molecular Biology Evolution [10], Chen and Zhang
argue that the variation of the mutation rate along the Escherichia coli genome
that we recently reported [3] cannot be evolutionarily optimised. To support
this claim they first attempt to calculate the selective advantage of a local
reduction in the mutation rate and conclude that it is not strong enough to be
favoured by selection. Second, they analyse the distribution of 166 mutations
from a wild-type E. coli K12 MG1655 strain and 1,346 mutations from a
repair-deficient strain, and claim to find a positive association between
transcription and mutation rate rather than the negative association that we
reported. Here we respond to this communication. Briefly, we explain how the
long-standing theory of mutation-modifier alleles supports the evolution of
local mutation rates within a genome by mechanisms acting on sufficiently large
regions of a genome, which is consistent with our original observations [3,4].
We then explain why caution must be exercised when comparing mutations from
repair deficient strains to data from wild-type strains, as different
mutational processes dominate these conditions. Finally, a reanalysis of the
data used by Zhang and Chen with an alternative expression dataset reveals that
their conclussions are unreliable.


  Understanding human mobility from a microscopic point of view may represent a
fundamental breakthrough for the development of a statistical physics for
cognitive systems and it can shed light on the applicability of macroscopic
statistical laws for social systems. Even if the complexity of individual
behaviors prevents a true microscopic approach, the introduction of mesoscopic
models allows the study of the dynamical properties for the non-stationary
states of the considered system. We propose to compute various entropy measures
of the individual mobility patterns obtained from GPS data that record the
movements of private vehicles in the Florence district, in order to point out
new features of human mobility related to the use of time and space and to
define the dynamical properties of a stochastic model that could generate
similar patterns. Moreover, we can relate the predictability properties of
human mobility to the distribution of time passed between two successive trips.
Our analysis suggests the existence of a hierarchical structure in the mobility
patterns which divides the performed activities into three different
categories, according to the time cost, with different information contents. We
show that a Markov process defined by using the individual mobility network is
not able to reproduce this hierarchy, which seems the consequence of different
strategies in the activity choice. Our results could contribute to the
development of governance policies for a sustainable mobility in modern cities.


  We use Brownian dynamics simulations to study the permeation properties of a
generic electrostatic model of a biological ion channel as a function of the
fixed charge $Q_f$ at its selectivity filter. We reconcile the
recently-discovered discrete calcium conduction bands M0 ($Q_f$=1e), M1 (3e),
M2 (5e) with the set of sodium conduction bands L0 (0.5-0.7e), L1 (1.5-2e)
thereby obtaining a completed pattern of conduction and selectivity bands v
$Q_f$ for the sodium-calcium channels family. An increase of $Q_f$ leads to an
increase of calcium selectivity: L0 (sodium selective, non-blocking channel) ->
M0 (non-selective channel) -> L1 (sodium selective channel with divalent block)
-> M1 (calcium selective channel exhibiting the anomalous mole fraction
effect). We create a consistent identification scheme where the L0 band is
identified with the eukaryotic (DEKA) sodium channel, and L1/L2 (speculatively)
with the bacterial NaChBac channel. The scheme created is able to account for
the experimentally observed mutation-induced transformations between
non-selective channels, sodium-selective channels, and calcium-selective
channels, which we interpret as transitions between different rows of the
identification table. By considering the potential energy changes during
permeation, we show explicitly that the multi-ion conduction bands of calcium
and sodium channels arise as the result of resonant barrier-less conduction.
Our results confirm the crucial influence of electrostatic interactions on
conduction and on the Ca/Na valence selectivity of calcium and sodium ion
channels. The model and results could be also applicable to biomimetic
nanopores with charged walls.


  Counting the frequencies of k-mers in read libraries is often a first step in
the analysis of high-throughput sequencing experiments. Infrequent k-mers are
assumed to be a result of sequencing errors. The frequent k-mers constitute a
reduced but error-free representation of the experiment, which can inform read
error correction or serve as the input to de novo assembly methods. Ideally,
the memory requirement for counting should be linear in the number of frequent
k-mers and not in the, typically much larger, total number of k-mers in the
read library.
  We present a novel method that balances time, space and accuracy requirements
to efficiently extract frequent k-mers even for high coverage libraries and
large genomes such as human. Our method is designed to minimize cache-misses in
a cache-efficient manner by using a Pattern-blocked Bloom filter to remove
infrequent k-mers from consideration in combination with a novel
sort-and-compact scheme, instead of a Hash, for the actual counting. While this
increases theoretical complexity, the savings in cache misses reduce the
empirical running times. A variant can resort to a counting Bloom filter for
even larger savings in memory at the expense of false negatives in addition to
the false positives common to all Bloom filter based approaches. A comparison
to the state-of-the-art shows reduced memory requirements and running times.
Note that we also provide the first competitive method to count k-mers up to
size 64.


  Much has been written about introduced rainbow trout interbreeding and
outcompeting native cutthroat trout. However, specific mechanisms have not been
thoroughly explored, and most data is limited to lotic ecosystems. Samples of
Snake River cutthroat trout (Oncorhynchus clarkii bouvieri), the
rainbow-cutthroat hybrid, the cutbow trout (Onchorhynchus mykiss x clarkii),
and rainbow trout (Oncorhynchus mykiss), were obtained from a lentic ecosystem
(Eleven Mile Reservoir, Colorado) by creel surveys conducted from May to
October, 2012. The total length and weight of each fish was measured and the
relative condition factor of each fish was computed using expected weight from
weight-length relationships from the Colorado Division of Parks and Wildlife
(CDPW). Data from the CDPW collected from 2003-2010 in the same lentic
ecosystem were used to compute relative condition factors for additional
comparison, as was independent creel survey data from 2011. Cutthroat trout
were plump: the mean relative condition factor of the cutthroat trout was
112.0% (+/- 1.0%). Cutbow hybrid trout were close to the expected weights with
a mean relative condition factor of 99.8% (+/- 0.6%). Rainbow trout were
thinner with a mean relative condition factor of 96.4% (+/- 1.4%). Comparing
mean relative condition factors of CDPW data from earlier years and plotting
the 2012 data relative to percentile curves also shows the same trend of
cutthroat trout being plumper than expected and rainbow trout being thinner
than the cutthroat trout, with the hybrid cutbow trout in between. This data
supports the hypothesis that rainbow trout do not outcompete cutthroat trout in
lentic ecosystems. Comparison with data from three other Colorado reservoirs
also shows that cutthroat trout tend to be more plump than rainbow trout and
their hybrids in sympatric lentic ecosystems.


  Modern microscopy techniques are aimed at imaging an individual molecule at
atomic resolution. Here we show that low-energy electrons with kinetic energies
of 50-250 eV offer a possibility of overcome the problem of radiation damage,
and obtaining images of individual biomolecules. Two experimental schemes for
obtaining images of individual molecules, holography and coherent diffraction
imaging, are discussed and compared. Images of individual molecules obtained by
both techniques, using low-energy electrons, are shown.


  Extensive research shows that more species-rich assemblages are generally
more productive and efficient in resource use than comparable assemblages with
fewer species. But the question of how diversity simultaneously affects the
wide variety of ecological functions that ecosystems perform remains relatively
understudied, and it presents several analytical and empirical challenges that
remain unresolved. In particular, researchers have developed several disparate
metrics to quantify multifunctionality, each characterizing different aspects
of the concept, and each with pros and cons. We compare four approaches to
characterizing multifunctionality and its dependence on biodiversity,
quantifying 1) magnitudes of multiple individual functions separately, 2) the
extent to which different species promote different functions, 3) the average
level of a suite of functions, and 4) the number of functions that
simultaneously exceed a critical threshold. We illustrate each approach using
data from the pan-European BIODEPTH experiment and the R multifunc package
developed for this purpose, evaluate the strengths and weaknesses of each
approach, and implement several methodological improvements. We conclude that a
extension of the fourth approach that systematically explores all possible
threshold values provides the most comprehensive description of
multifunctionality to date. We outline this method and recommend its use in
future research.


  This paper proposed a discrete stochastic dynamic programming (SDP) model for
sustainable ecosystem (SE) planning of the Loess Plateau in Northwestern,
China, and analyzed the ecological resource planning by the evolutionary game
model in the decision-making process. The main objective is to explore a new
approach of SE planning from a viewpoint of discrete SDP and evolutionary game
theory, with a specific application in the area of ecological resource planning
such as water management problems. In contrast to previous work, the proposed
SDP method focuses on the transition probability matrix of the ecosystem in a
statistic sense, and uses the DP algorithm to obtain the optimal ecological
resource planning strategies among multi-subsystems, then analyzes impacts of
decision between different users. Firstly, the application background and the
concept of SE planning are introduced. Then, a brief overview of existing
theory for analyzing sustainable ecosystem is presented. Furthermore, a
SDP-based mathematical model and its application to water resource planning of
central areas of Loess Plateau are presented as an example. Finally,
supplementary analysis of impacts between different users in SE planning as a
game playing is provided.


  Human populations have undergone dramatic changes in population size in the
past 100,000 years, including a severe bottleneck of non-African populations
and recent explosive population growth. There is currently great interest in
how these demographic events may have affected the burden of deleterious
mutations in individuals and the allele frequency spectrum of disease mutations
in populations. Here we use population genetic models to show that--contrary to
previous conjectures--recent human demography has likely had very little impact
on the average burden of deleterious mutations carried by individuals. This
prediction is supported by exome sequence data showing that African American
and European American individuals carry very similar burdens of damaging
mutations. We next consider whether recent population growth has increased the
importance of very rare mutations in complex traits. Our analysis predicts that
for most classes of disease variants, rare alleles are unlikely to contribute a
large fraction of the total genetic variance, and that the impact of recent
growth is likely to be modest. However, for diseases that have a direct impact
on fitness, strongly deleterious rare mutations likely do play important roles,
and the impact of very rare mutations will be far greater as a result of recent
growth. In summary, demographic history has dramatically impacted patterns of
variation in different human populations, but these changes have likely had
little impact on either genetic load or on the importance of rare variants for
most complex traits.


  Gene co-expression databases in model plants such as Arabidopsis and rice
have been extensively developed and utilised for predicting gene function and
identifying functional modules for underlying biological processes. However,
such tools are less available for non-model plants such as grapevine. We have
constructed a gene co-expression database, VTCdb
(http://vtcdb.adelaide.edu.au/Home.aspx) that offers an online platform for
transcriptional regulatory inference in the cultivated grapevine. Using a
condition-independent approach, the grapevine co-expression network was
constructed using 352 publicly available microarray datasets from a diverse
experimental series, profiling approximately 9000 genes (40% of the predicted
grapevine transcriptome). Using correlation rank transformation and
graph-clustering, we have identified modules putatively involved in several
fundamental biological processes such as photosynthesis, secondary metabolism
and stress responses. Inter-module network connections revealed a higher-level
organization of function, in which densely connected modules often participated
in related biological processes and had similar expression profiles in
grapevine. The database enables users to query genes, modules or biological
processes of interest. Querying a gene will give a ranked list of co-expressed
genes, functional annotations and information on the associated module.
Alternatively, browsing modules of interest will retrieve information regarding
Gene Ontology and Mapman enriched terms, as well as tissue- and
condition-specific patterns of gene expression within the module. Furthermore,
the database features interactive network visualisation via CytoscapeWeb. VTCdb
should aid researchers in their prioritization of gene candidates for further
study towards the understanding of biological processes related to many aspects
of grapevine development and metabolism.


  Ions in water are important in biology, from molecules to organs.
Classically, ions in water are treated as ideal noninteracting particles in a
perfect gas. Excess free energy of ion was zero. Mathematics was not available
to deal consistently with flows, or interactions with ions or boundaries.
Non-classical approaches are needed because ions in biological conditions flow
and interact. The concentration gradient of one ion can drive the flow of
another, even in a bulk solution. A variational multiscale approach is needed
to deal with interactions and flow. The recently developed energetic
variational approach to dissipative systems allows mathematically consistent
treatment of bio-ions Na, K, Ca and Cl as they interact and flow. Interactions
produce large excess free energy that dominate the properties of the high
concentration of ions in and near protein active sites, channels, and nucleic
acids: the number density of ions is often more than 10 M. Ions in such crowded
quarters interact strongly with each other as well as with the surrounding
protein. Non-ideal behavior has classically been ascribed to allosteric
interactions mediated by protein conformation changes. Ion-ion interactions
present in crowded solutions--independent of conformation changes of
proteins--are likely to change interpretations of allosteric phenomena.
Computation of all atoms is a popular alternative to the multiscale approach.
Such computations involve formidable challenges. Biological systems exist on
very different scales from atomic motion. Biological systems exist in ionic
mixtures (extracellular/intracellular solutions), and usually involve flow and
trace concentrations of messenger ions (e.g., 10-7 M Ca2+). Energetic
variational methods can deal with these characteristic properties of biological
systems while we await the maturation and calibration of all atom simulations
of ionic mixtures and divalents.


  Extreme enrichment of the human immunodeficiency virus (HIV-1) RNA genome for
the purine A parallels the mild purine-loading of the RNAs of most organisms.
This should militate against loop-loop "kissing" interactions between the
structured viral genome and structured host RNAs, which can generate segments
of double-stranded RNA sufficient to trigger intracellular alarms. However,
human T cell leukaemia virus (HTLV-1), with the potential to invade the same
host cell, shows extreme enrichment for the pyrimidine C. Assuming the low GC%
HIV and the high GC% HTLV-1 to share a common ancestor, it was postulated that
differences in GC% arose to prevent homologous recombination between these
emerging lentiviral species. Sympatrically isolated by this intracellular
reproductive barrier, prototypic HIV-1 seized the AU-rich (low GC%) high ground
(thus committing to purine A rather than purine G). Prototypic HTLV-1 forwent
this advantage and evolved an independent evolutionary strategy. Evidence
supporting this hypothesis since its elaboration in the 1990s is growing. The
conflict between the needs to encode accurately both a protein, and nucleic
acid structure, is often resolved in favour of the nucleic acid because, apart
from regulatory roles, structure is critical for recombination. However, above
a sequence difference threshold, structure (and hence recombination) is
impaired. New species can then arise.


  We apply optimal control theory to a tuberculosis model given by a system of
ordinary differential equations. Optimal control strategies are proposed to
minimize the cost of interventions, considering reinfection and post-exposure
interventions. They depend on the parameters of the model and reduce
effectively the number of active infectious and persistent latent individuals.
The time that the optimal controls are at the upper bound increase with the
transmission coefficient. A general explicit expression for the basic
reproduction number is obtained and its sensitivity with respect to the model
parameters is discussed. Numerical results show the usefulness of the
optimization strategies.


  Cooperative action of molecular motors is essential for many cellular
processes. One possible regulator of motor coordination is the
elasticity-mediated crosstalk (EMC) coupling between myosin II motors whose
origin is the tensile stress that they collectively generate in actin
filaments. Here, we use a statistical mechanical analysis to investigate the
influence of the EMC effect on the sarcomere - the basic contractile unit of
skeletal muscles. We demonstrate that the EMC effect leads to an increase in
the attachment probability of motors located near the end of the sarcomere
while, simultaneously, decreasing the attachment probability of the motors in
the central part. Such a polarized attachment probability would impair the
motors ability to cooperate efficiently. Interestingly, this undesired
phenomenon becomes significant only when the system size exceeds that of the
sarcomere in skeletal muscles, which provides an explanation for the remarkable
lack of sarcomere variability in vertebrates. Another phenomenon that we
investigate is the recently observed increase in the duty ratio of the motors
with the tension in muscle. We reveal that the celebrated Hill's equation for
muscle contraction is very closely related to this observation.


  We present the ProCS method for the rapid and accurate prediction of protein
backbone amide proton chemical shifts - sensitive probes of the geometry of key
hydrogen bonds that determine protein structure. ProCS is parameterized against
quantum mechanical (QM) calculations and reproduces high level QM results
obtained for a small protein with an RMSD of 0.25 ppm (r = 0.94). ProCS is
interfaced with the PHAISTOS protein simulation program and is used to infer
statistical protein ensembles that reflect experimentally measured amide proton
chemical shift values. Such chemical shift-based structural refinements,
starting from high-resolution X-ray structures of Protein G, ubiquitin, and SMN
Tudor Domain, result in average chemical shifts, hydrogen bond geometries, and
trans-hydrogen bond (h3JNC') spin-spin coupling constants that are in excellent
agreement with experiment. We show that the structural sensitivity of the
QM-based amide proton chemical shift predictions is needed to refine protein
structures to this agreement. The ProCS method thus offers a powerful new tool
for refining the structures of hydrogen bonding networks to high accuracy with
many potential applications such as protein flexibility in ligand binding.


  A small but growing number of people are finding interesting parallels
between ecosystems as studied by ecologists (think of a Savanna or the Amazon
rain forest or a Coral reef) and tumours1-3. The idea of viewing cancer from an
ecological perspective has many implications but fundamentally, it means that
we should not see cancer just as a group of mutated cells. A more useful
definition of cancer is to consider it a disruption in the complex balance of
many interacting cellular and microenvironmental elements in a specific organ.
This perspective means that organs undergoing carcinogenesis should be seen as
sophisticated ecosystems in homeostasis that cancer cells can disrupt. It also
makes cancer seem even more complex but may ultimately provides isights that
make it more treatable. Here we discuss how ecological principles can be used
to better understand cancer progression and treatment, using several
mathematical and computational models to illustrate our argument.


  The population dynamics that evolves in the radial symmetric geometry is
investigated. The nonlinear reaction-diffusion model, which depends on
population density, is employed as the governing equation for this system. The
approximate analytical solution to this equation has been found. It shows that
the population density evolves from initial state and propagates as the
traveling wave-like for the large time scale. One can be mentioned that, if the
distance is insufficient large, the curvature has ineluctable influence on
density profile and front speed. In comparison, the analytical solution is in
agreement with the numerical solution.


  Spreading processes represent a very efficient tool to investigate the
structural properties of networks and the relative importance of their
constituents, and have been widely used to this aim in static networks. Here we
consider simple disease spreading processes on empirical time-varying networks
of contacts between individuals, and compare the effect of several immunization
strategies on these processes. An immunization strategy is defined as the
choice of a set of nodes (individuals) who cannot catch nor transmit the
disease. This choice is performed according to a certain ranking of the nodes
of the contact network. We consider various ranking strategies, focusing in
particular on the role of the training window during which the nodes'
properties are measured in the time-varying network: longer training windows
correspond to a larger amount of information collected and could be expected to
result in better performances of the immunization strategies. We find instead
an unexpected saturation in the efficiency of strategies based on nodes'
characteristics when the length of the training window is increased, showing
that a limited amount of information on the contact patterns is sufficient to
design efficient immunization strategies. This finding is balanced by the large
variations of the contact patterns, which strongly alter the importance of
nodes from one period to the next and therefore significantly limit the
efficiency of any strategy based on an importance ranking of nodes. We also
observe that the efficiency of strategies that include an element of randomness
and are based on temporally local information do not perform as well but are
largely independent on the amount of information available.


  An ongoing debate in evolutionary biology is whether phenotypic change occurs
predominantly around the time of speciation or whether it instead accumulates
gradually over time. In this work I propose a general framework incorporating
both types of change, quantify the effects of speciational change via the
correlation between species and attribute the proportion of change to each
type. I discuss results of parameter estimation of Hominoid body size in this
light. I derive mathematical formulae related to this problem, the probability
generating functions of the number of speciation events along a randomly drawn
lineage and from the most recent common ancestor of two randomly chosen tip
species for a conditioned Yule tree. Additionally I obtain in closed form the
variance of the distance from the root to the most recent common ancestor of
two randomly chosen tip species.


  Analysis of network dynamics became a focal point to understand and predict
changes of complex systems. Here we introduce Turbine, a generic framework
enabling fast simulation of any algorithmically definable dynamics on very
large networks. Using a perturbation transmission model inspired by
communicating vessels, we define a novel centrality measure: perturbation
centrality. Hubs and inter-modular nodes proved to be highly efficient in
perturbation propagation. High perturbation centrality nodes of the Met-tRNA
synthetase protein structure network were identified as amino acids involved in
intra-protein communication by earlier studies. Changes in perturbation
centralities of yeast interactome nodes upon various stresses well
recapitulated the functional changes of stressed yeast cells. The novelty and
usefulness of perturbation centrality was validated in several other model,
biological and social networks. The Turbine software and the perturbation
centrality measure may provide a large variety of novel options to assess
signaling, drug action, environmental and social interventions. The Turbine
algorithm is available at: http://www.turbine.linkgroup.hu


  This paper considers a two-dimensional logistic model to study populations
with two genders. The growth behavior of a population is guided by two coupled
ordinary differential equations given by a non-differentiable vector field
whose parameters are the secondary sex ratio (the ratio of males to females at
time of birth), inter-, intra- and outer-gender competitions, fertility and
mortality rates and a mating function. For the case where there is no
inter-gender competition and the mortality rates are negligible with respect to
the density-dependent mortality, using geometrical techniques, we analyze the
singularities and the basin of attraction of the system, determining the
relationships between the parameters for which the system presents an
equilibrium point. In particular, we describe conditions on the secondary sex
ratio and discuss the role of the average number of female sexual partners of
each male for the conservation of a two-sex species.


  We propose the following model for speciation and extinction. Birth and
deaths occur according to spatially inhomogeneous contact rates. We assume that
the ratio of the birth rate over the death rate at a site converges to some
limit as the site goes to infinity. We show that the existence of a phase
transition depends critically on the value of this limit.