  The flexibility-rigidity index (FRI) is a newly proposed method for the
construction of atomic rigidity functions. The FRI method analyzes protein
rigidity and flexibility and is capable of predicting protein B-factors without
resorting to matrix diagonalization. A fundamental assumption used in the FRI
is that protein structures are uniquely determined by various internal and
external interactions, while the protein functions, such as stability and
flexibility, are solely determined by the structure. As such, one can predict
protein flexibility without resorting to the protein interaction Hamiltonian.
Consequently, bypassing the matrix diagonalization, the original FRI has a
computational complexity of O(N^2). This work introduces a fast FRI (fFRI)
algorithm for the flexibility analysis of large macromolecules. The proposed
fFRI further reduces the computational complexity to O(N). Additionally, we
propose anisotropic FRI (aFRI) algorithms for the analysis of protein
collective dynamics. The aFRI algorithms admit adaptive Hessian matrices, from
a completely global 3N*3N matrix to completely local 3*3 matrices. However,
these local 3*3 matrices have built in much non-local correlation. Furthermore,
we compare the accuracy and efficiency of FRI with some {established}
approaches to flexibility analysis, namely, normal mode analysis (NMA) and
Gaussian network model (GNM). The accuracy of the FRI method is tested. The
FRI, particularly the fFRI, is orders of magnitude more efficient and about 10%
more accurate overall than some of the most popular methods in the field. The
proposed fFRI is able to predict B-factors for alpha-carbons of the HIV virus
capsid (313,236 residues) in less than 30 seconds on a single processor using
only one core. Finally, we demonstrate the application of FRI and aFRI to
protein domain analysis.


  We formulate and analyze a theoretical model for the regulation of
microtubule (MT) polymerization dynamics by the signaling proteins Rac1 and
stathmin. In cells, the MT growth rate is inhibited by cytosolic stathmin,
which, in turn, is inactivated by Rac1. Growing MTs activate Rac1 at the cell
edge, which closes a positive feedback loop. We investigate both tubulin
sequestering and catastrophe promotion as mechanisms for MT growth inhibition
by stathmin. For a homogeneous stathmin concentration in the absence of Rac1,
we find a switch-like regulation of the MT mean length by stathmin. For
constitutively active Rac1 at the cell edge, stathmin is deactivated locally,
which establishes a spatial gradient of active stathmin. In this gradient, we
find a stationary bimodal MT length distributions for both mechanisms of MT
growth inhibition by stathmin. One subpopulation of the bimodal length
distribution can be identified with fast growing and long pioneering MTs in the
region near the cell edge, which have been observed experimentally. The
feedback loop is closed through Rac1 activation by MTs. For tubulin
sequestering by stathmin, this establishes a bistable switch with two stable
states: one stable state corresponds to upregulated MT mean length and bimodal
MT length distributions, i.e., pioneering MTs; the other stable state
corresponds to an interrupted feedback with short MTs. Stochastic effects as
well as external perturbations can trigger switching events. For catastrophe
promoting stathmin we do not find bistability.


  Solvation is essential for protein activities. To study internal solvation of
protein, site-directed mutagenesis is applied. Intrinsic fluorescent probe,
tryptophan, is inserted into desired position inside protein molecule for
ultrafast spectroscopic study. Here we review this unique method for protein
dynamics researches. We introduce the frontiers of protein solvation,
site-directed mutagenesis, protein stability and characteristics, and the
spectroscopic methods. Then we present time-resolved spectroscopic dynamics of
solvation dynamics inside caves of active sites. The studies are carried out on
a globular protein, staphylococcal nuclease. The solvation at internal sites of
the caves indicate clear characteristics of local environment. These solvation
behaviors correlated to the enzyme activity directly.


  The progress of the successive rounds of catalytic conversion of substrates
into product(s) by a single enzyme is characterized by the distribution of
turnover times. Establishing the most general form of dependence of this
distribution on the substrate concentration [S] is one of the fundamental
challenges in single molecule enzymology. The distribution of the times of
dwell of a molecular motor at the successive positions on its track is an
analogous quantity. We derive approximate series expansions for the
[ATP]-dependence of the first two moments of the dwell time distributions of
motors that catalyze hydrolysis of ATP to draw input energy. Comparison between
our results for motors with branched pathways and the corresponding expressions
reported earlier for linear enzymatic pathways provides deep insight into the
effects of the branches. Such insight is likely to help in discovering the most
general form of [S]-dependence of these fundamental distributions.


  Superoxide reductase (SOR) is an Fe protein that catalyzes the reduction of
superoxide to give H(2)O(2). Recently, the mutation of the Glu47 residue into
alanine (E47A) in the active site of SOR from Desulfoarculus baarsii has
allowed the stabilization of an iron-peroxo species when quickly reacted with
H(2)O(2) [Math{\'e} et al. (2002) J. Am. Chem. Soc. 124, 4966-4967]. To further
investigate this non-heme peroxo-iron species, we have carried out a
M{\"o}ssbauer study of the (57)Fe-enriched E47A SOR from D. baarsii reacted
quickly with H(2)O(2). Considering the M{\"o}ssbauer data, we conclude, in
conjunction with the other spectroscopic data available and with the results of
density functional calculations on related models, that this species
corresponds to a high-spin side-on peroxo-Fe(3+) complex. This is one of the
first examples of such a species in a biological system for which M{\"o}ssbauer
parameters are now available: delta(/Fe) = 0.54 (1) mm/s, DeltaE(Q) = -0.80 (5)
mm/s, and the asymmetry parameter eta = 0.60 (5) mm/s. The M{\"o}ssbauer and
spin Hamiltonian parameters have been evaluated on a model from the side-on
peroxo complex (model 2) issued from the oxidized iron center in SOR from
Pyrococcus furiosus, for which structural data are available in the literature
[Yeh et al. (2000) Biochemistry 39, 2499-2508]. For comparison, similar
calculations have been carried out on a model derived from 2 (model 3), where
the [CH(3)-S](1)(-) group has been replaced by the neutral [NH(3)](0) group
[Neese and Solomon (1998) J. Am. Chem. Soc. 120, 12829-12848]. Both models 2
and 3 contain a formally high-spin Fe(3+) ion (i.e., with empty minority spin
orbitals). We found, however, a significant fraction (approximately 0.6 for 2,
approximately 0.8 for 3) of spin (equivalently charge) spread over two occupied
(minority spin) orbitals. The quadrupole splitting value for 2 is found to be
negative and matches quite well the experimental value. The computed quadrupole
tensors are rhombic in the case of 2 and axial in the case of 3. This
difference originates directly from the presence of the thiolate ligand in 2. A
correlation between experimental isomer shifts for Fe(3+) mononuclear complexes
with computed electron densities at the iron nucleus has been built and used to
evaluate the isomer shift values for 2 and 3 (0.56 and 0.63 mm/s,
respectively). A significant increase of isomer shift value is found upon going
from a methylthiolate to a nitrogen ligand for the Fe(3+) ion, consistent with
covalency effects due to the presence of the axial thiolate ligand. Considering
that the isomer shift value for 3 is likely to be in the 0.61-0.65 mm/s range
[Horner et al. (2002) Eur. J. Inorg. Chem., 3278-3283], the isomer shift value
for a high-spin eta(2)-O(2) Fe(3+) complex with an axial thiolate group can be
estimated to be in the 0.54-0.58 mm/s range. The occurrence of a side-on peroxo
intermediate in SOR is discussed in relation to the recent data published for a
side-on peroxo-Fe(3+) species in another biological system [Karlsson et al.
(2003) Science 299, 1039-1042].


  The possibility of spreading by migration, colonizing new spatial domains
suitable for development and reproduction, can substantially relieve a
biological population from the risk of extinction. By means of a realistic
computational model based on empirical data, we study this phenomenon for the
human head louse, Pediculus humanus capitis. In particular, we show that a lice
colony infesting a single isolated host is prone to extinction by stochastic
population fluctuations within an interval of several months, while migration
over a relatively small group of hosts in contact with each other is enough to
insure the prevalence of the infestation for indefinitely long periods. We
characterize the interplay of the size of the host group with the host-to-host
contagion probability, which controls a transition between extinction of the
lice population and a situation where the infestation is endemic.


  The heterogeneity of reaction fluxes present in a metabolic network within a
single flux state can be exploited to construct the so-called backbone as a
reduced version of metabolism. The backbone maintains all significant fluxes
producing or consuming metabolites while displaying a substantially decreased
number of interconnections and, hence, it becomes a useful tool to extract
primary metabolic routes. Here, we disclose the metabolic backbone of
Escherichia coli using the computationally predicted fluxes which maximize the
growth rate in glucose minimal medium, and we compare it with the backbone of
Mycoplasma pneumoniae, a much simpler organism. We find that the central core
in both backbones is mainly composed of reactions in ancient pathways, still
playing at present a key role in energy metabolism. In E. coli, the analysis of
the backbone reveals that the synthesis of nucleotides and the metabolism of
lipids form smaller cores which rely critically on energy metabolism; but not
conversely. At the same time, an analysis of the dependence of this backbone on
media composition leads to the identification of pathways sensitive to
environmental changes. The metabolic backbone of an organism is thus useful to
trace simultaneously both its evolution and adaptation fingerprints.


  By incorporating segregated spatial domain and individual-based linkage into
the SIS (susceptible-infected-susceptible) model, we investigate the coupled
effects of random walk and intragroup interaction on contagion. Compared with
the situation where only local movement or individual-based linkage exists, the
coexistence of them leads to a wider spread of infectious disease. The roles of
narrowing segregated spatial domain and reducing mobility in epidemic control
are checked, these two measures are found to be conducive to curbing the spread
of infectious disease. Considering heterogeneous time scales between local
movement and global interaction, a log-log relation between the change in the
number of infected individuals and the timescale $\tau$ is found. A theoretical
analysis indicates that the evolutionary dynamics in the present model is
related to the encounter probability and the encounter time. A functional
relation between the epidemic threshold and the ratio of shortcuts, and a
functional relation between the encounter time and the timescale $\tau$ are
found.


  In a recent paper Siegfried et al. published a new sequence-based structural
RNA assay that utilizes mutational profiling to detect base pairing (MaP).
Output from MaP provides information about both pairing (via reactivities) and
contact (via correlations). Reactivities can be coupled to partition function
folding models for structural inference, while correlations can reveal pairs of
sites that may be in structural proximity. The possibility for inference of 3D
contacts via MaP suggests a novel approach to structural prediction for RNA
analogous to covariance structural prediction for proteins. We explore this
approach and show that partial correlation analysis outperforms na\"ive
correlation analysis. Our results should be applicable to a wide range of
high-throughput sequencing based RNA structural assays that are under
development.


  We consider a stochastic population model where the intrinsic or demographic
noise causes cycling between states before the population eventually goes
extinct. A master equation approach coupled with a WKB
(Wentzel-Kramers-Brillouin) approximation is used to construct the optimal path
to extinction. In addition, a probabilistic argument is used to understand the
pre-extinction dynamics and approximate the mean time to extinction. Analytical
results agree well with numerical Monte Carlo simulations. A control method is
implemented to decrease the mean time to extinction. Analytical results
quantify the effectiveness of the control and agree well with numerical
simulations.


  With the growing amount of available temporal real-world network data, an
important question is how to efficiently study these data. One can simply model
a temporal network as either a single aggregate static network, or as a series
of time-specific snapshots, each of which is an aggregate static network over
the corresponding time window. The advantage of modeling the temporal data in
these two ways is that one can use existing well established methods for static
network analysis to study the resulting aggregate network(s). Here, we develop
a novel approach for studying temporal network data more explicitly. We base
our methodology on the well established notion of graphlets (subgraphs), which
have been successfully used in numerous contexts in static network research.
Here, we take the notion of static graphlets to the next level and develop new
theory needed to allow for graphlet-based analysis of temporal networks. Our
new notion of dynamic graphlets is quite different than existing approaches for
dynamic network analysis that are based on temporal motifs (statistically
significant subgraphs). Namely, these approaches suffer from many limitations.
For example, they can only deal with subgraph structures of limited complexity.
Also, their major drawback is that their results heavily depend on the choice
of a null network model that is required to evaluate the significance of a
subgraph. However, choosing an appropriate null network model is a non-trivial
task. Our dynamic graphlet approach overcomes the limitations of the existing
temporal motif-based approaches. At the same time, when we thoroughly evaluate
the ability of our new approach to characterize the structure and function of
an entire temporal network or of individual nodes, we find that the dynamic
graphlet approach outperforms the static graphlet approach, which indicates
that accounting for temporal information helps.


  On rugged fitness landscapes where sign epistasis is common, adaptation can
often involve either individually beneficial "uphill" mutations or more complex
mutational trajectories involving fitness valleys or plateaus. The dynamics of
the evolutionary process determine the probability that evolution will take any
specific path among a variety of competing possible trajectories. Understanding
this evolutionary choice is essential if we are to understand the outcomes and
predictability of adaptation on rugged landscapes. We present a simple model to
analyze the probability that evolution will eschew immediately uphill paths in
favor of crossing fitness valleys or plateaus that lead to higher fitness but
less accessible genotypes. We calculate how this probability depends on the
population size, mutation rates, and relevant selection pressures, and compare
our analytical results to Wright-Fisher simulations. We find that the
probability of valley crossing depends nonmonotonically on population size:
intermediate size populations are most likely to follow a "greedy" strategy of
acquiring immediately beneficial mutations even if they lead to evolutionary
dead ends, while larger and smaller populations are more likely to cross
fitness valleys to reach distant advantageous genotypes. We explicitly identify
the boundaries between these different regimes in terms of the relevant
evolutionary parameters. Above a certain threshold population size, we show
that the degree of evolutionary "foresight" depends only on a single simple
combination of the relevant parameters.


  Feline infectious peritonitis (FIP) is associated with mutations in the
feline coronavirus (FCoV) genome that are thought to convert the subclinical
feline enteric coronavirus (FECV) into the lethal feline infectious peritonitis
virus (FIPV). A key feature of FIPV, not shared with FECV, is the productive
infection of macrophages. Therefore mutations in proteins that govern cell
tropism, such as the spike glycoprotein, may play an important role in FIP
progression. In a prior study, involving a limited number of samples, we have
shown an association of FIP with mutations in the protease cleavage-activation
site located between the receptor-binding and fusion domains of the FCoV spike
(S1/S2). Here, we extend these studies to investigate a larger sample set and
to obtain a more refined analysis of the mutations at this S1/S2 site. Our
larger data set more clearly shows that the mutations acquired by FIPV at S1/S2
are also accompanied by additional mutations at a second protease
cleavage-activation site located in the fusion domain (S2'), adjacent to the
viral fusion peptide. Overall, our data indicate a pattern of mutations across
the two protease recognition sites that results in substitutions, and/or
altered recognition, of critical basic/polar amino acid residues needed for
virus activation in the enteric tract. Typically, FIPVs have substitutions of
non-polar, aliphatic or aromatic residues in the protease recognition sites.
These changes likely modulate the proteolytic activation of the virus and its
ability to productively infect macrophages in vivo.


  The multi-level method for discrete state systems, first introduced by
Anderson and Higham [Multiscale Model. Simul. 10:146--179, 2012], is a highly
efficient simulation technique that can be used to elucidate statistical
characteristics of biochemical reaction networks. A single point estimator is
produced in a cost-effective manner by combining a number of estimators of
differing accuracy in a telescoping sum, and, as such, the method has the
potential to revolutionise the field of stochastic simulation. The first term
in the sum is calculated using an approximate simulation algorithm, and can be
calculated quickly but is of significant bias. Subsequent terms successively
correct this bias by combining estimators from approximate stochastic
simulations algorithms of increasing accuracy, until a desired level of
accuracy is reached.
  In this paper we present several refinements of the multi-level method which
render it easier to understand and implement, and also more efficient. Given
the substantial and complex nature of the multi-level method, the first part of
this work (Sections 2 - 5) is written as a tutorial, with the aim of providing
a practical guide to its use. The second part (Sections 6 - 8) takes on a form
akin to a research article, thereby providing the means for a deft
implementation of the technique, and concludes with a discussion of a number of
open problems.


  Within the field of phylogenetics there is great interest in distance
measures to quantify the dissimilarity of two trees. Recently, a new distance
measure has been proposed: the Maximum Parsimony (MP) distance. This is based
on the difference of the parsimony scores of a single character on both trees
under consideration, and the goal is to find the character which maximizes this
difference. Here we show that computation of MP distance on two \emph{binary}
phylogenetic trees is NP-hard. This is a highly nontrivial extension of an
earlier NP-hardness proof for two multifurcating phylogenetic trees, and it is
particularly relevant given the prominence of binary trees in the phylogenetics
literature. As a corollary to the main hardness result we show that computation
of MP distance is also hard on binary trees if the number of states available
is bounded. In fact, via a different reduction we show that it is hard even if
only two states are available. Finally, as a first response to this hardness we
give a simple Integer Linear Program (ILP) formulation which is capable of
computing the MP distance exactly for small trees (and for larger trees when
only a small number of character states are available) and which is used to
computationally verify several auxiliary results required by the hardness
proofs.


  The pursuit of highest payoffs in evolutionary social dilemmas is risky and
sometimes inferior to conformity. Choosing the most common strategy within the
interaction range is safer because it ensures that the payoff of an individual
will not be much lower than average. Herding instincts and crowd behavior in
humans and social animals also compel to conformity on their own right.
Motivated by these facts, we here study the impact of conformity on the
evolution of cooperation in social dilemmas. We show that an appropriate
fraction of conformists within the population introduces an effective surface
tension around cooperative clusters and ensures smooth interfaces between
different strategy domains. Payoff-driven players brake the symmetry in favor
of cooperation and enable an expansion of clusters past the boundaries imposed
by traditional network reciprocity. This mechanism works even under the most
testing conditions, and it is robust against variations of the interaction
network as long as degree-normalized payoffs are applied. Conformity may thus
be beneficial for the resolution of social dilemmas.


  We present RKappa, a framework for the development and analysis of rule-based
models within a mature, statistically empowered R environment. The
infrastructure allows model editing, modification, parameter sampling,
simulation, statistical analysis and visualisation without leaving the R
environment. We demonstrate its effectiveness through its application to Global
Sensitivity Analysis, exploring it in "parallel" and "concurrent"
implementations.
  The pipeline was designed for high performance computing platforms and aims
to facilitate analysis of the behaviour of large-scale systems with limited
knowledge of exact mechanisms and respectively sparse availability of parameter
values, and is illustrated here with two biological examples.
  The package is available on github: https://github.com/lptolik/R4Kappa


  Genomic data are rapidly resolving the tree of living species calibrated to
time, the timetree of life, which will provide a framework for research in
diverse fields of science. Previous analyses of taxonomically restricted
timetrees have found a decline in the rate of diversification in many groups of
organisms, often attributed to ecological interactions among species. Here we
have synthesized a global timetree of life from 2,274 studies representing
50,632 species and examined the pattern and rate of diversification as well as
the timing of speciation. We found that species diversity has been mostly
expanding overall and in many smaller groups of species, and that the rate of
diversification in eukaryotes has been mostly constant. We also identified, and
avoided, potential biases that may have influenced previous analyses of
diversification including low levels of taxon sampling, small clade size, and
the inclusion of stem branches in clade analyses. We found consistency in
time-to-speciation among plants and animals, approximately two million years,
as measured by intervals of crown and stem species times. Together, this
clock-like change at different levels suggests that speciation and
diversification are processes dominated by random events and that adaptive
change is largely a separate process.


  Fire is an indissoluble component of ecosystems, however quantifying the
effects of fire on vegetation is challenging task as fire lies outside the
typical experimental design attributes. A recent simulation study showed that
under increased fire regimes positive tree-tree interactions were recorded
(Bacelar et al., 2014). Data from experimental burning plots in an African
savanna, the Kruger National Park, were collected across unburnt and annual
burn plots. Indices of aggregation and spatial autocorrelation of the
distribution of trees between different fire regimes were explored. Results
show that the distribution of trees under fire were more clumped and exhibited
higher spatial autocorrelation than in unburnt plots: In burnt plots spatial
autocorrelation values were positive at finer scales and negative at coarser
scales potentially indicating co-existence of facilitation and competition
within the same ecosystem depending on the scale. The pattern derived here
provides inference for (a) fire acting as an increasing aggregation & spatial
autocorrelation force (b) tree survival under fire regimes is potentially
facilitated by forming patches of trees (c) scale-dependent facilitation and
competition coexisting within the same ecosystem with finer scale facilitation
and coarser scale competition.


  Epidemic outbreaks of new pathogens, or known pathogens in new populations,
cause a great deal of fear because they are hard to predict. For theoretical
models of disease spreading, on the other hand, quantities characterizing the
outbreak converge to deterministic functions of time. Our goal in this paper is
to shed some light on this apparent discrepancy. We measure the diversity of
(and, thus, the predictability of) outbreak sizes and extinction times as
functions of time given different scenarios of the amount of information
available. Under the assumption of perfect information -- i.e., knowing the
state of each individual with respect to the disease -- the predictability
decreases exponentially, or faster, with time. The decay is slowest for
intermediate values of the per-contact transmission probability. With a weaker
assumption on the information available, assuming that we know only the
fraction of currently infectious, recovered, or susceptible individuals, the
predictability also decreases exponentially most of the time. There are,
however, some peculiar regions in this scenario where the predictability
decreases. In other words, to predict its final size with a given accuracy, we
would need increasingly more information about the outbreak.


  Rafts, or functional domains, are transient nano- or mesoscopic structures in
the plasma membrane and are thought to be essential for many cellular processes
such as signal transduction, adhesion, trafficking and lipid/protein sorting.
Observations of these membrane heterogeneities have proven challenging, as they
are thought to be both small and short-lived. With a combination of
coarse-grained molecular dynamics simulations and neutron diffraction using
deuterium labeled cholesterol molecules we observe raft-like structures and
determine the ordering of the cholesterol molecules in binary
cholesterol-containing lipid membranes. From coarse-grained computer
simulations, heterogenous membranes structures were observed and characterized
as small, ordered domains. Neutron diffraction was used to study the lateral
structure of the cholesterol molecules. We find pairs of strongly bound
cholesterol molecules in the liquid-disordered phase, in accordance with the
umbrella model. Bragg peaks corresponding to ordering of the cholesterol
molecules in the raft-like structures were observed and indexed by two
different structures: a monoclinic structure of ordered cholesterol pairs of
alternating direction in equilibrium with cholesterol plaques, i.e., triclinic
cholesterol bilayers.


  Taylor's law (TL) states that the variance $V$ of a non-negative random
variable is a power function of its mean $M$, i.e. $V=a M^b$. The ubiquitous
empirical verification of TL, typically displaying sample exponents $b \simeq
2$, suggests a context-independent mechanism. However, theoretical studies of
population dynamics predict a broad range of values of $b$. Here, we explain
this apparent contradiction by using large deviations theory to derive a
generalized TL in terms of sample and populations exponents $b_{jk}$ for the
scaling of the $k$-th vs the $j$-th cumulant (conventional TL is recovered for
$b=b_{12}$), with the sample exponent found to depend predictably on the number
of observed samples. Thus, for finite numbers of observations one observes
sample exponents $b_{jk}\simeq k/j$ (thus $b\simeq2$) independently of
population exponents. Empirical analyses on two datasets support our
theoretical results.


  Superoxide reductase (SOR) catalyzes the reduction of O2*- to H2O2. Its
active site consists of a non-heme Fe2+ center in an unusual square-pyramidal
[His4 Cys] coordination. Like many SORs, the electronic absorption band
corresponding to the oxidized active site of the SOR from Desulfoarculus
baarsii exhibits a pH-dependent alkaline transition changing from ca. 644 to
560 nm as the pH increases and with an apparent pKa of 9.0. Variants in which
the conserved amino acids glutamate 47 and lysine 48 were replaced by the
neutral residues alanine (E47A) and isoleucine (K48I), respectively, exhibited
the same alkaline transition but at lower apparent pKa values of 6.7 and 7.6,
respectively. Previous work [Nivi{\`e}re, V.; Asso, M.; Weill, C. O.; Lombard,
M.; Guigliarelli, B.; Favaudon, V.; Hou{\'e}e-Levin, C. Biochemistry 2004, 43,
808-818] has shown that this alkaline transition is associated with the
protonation/deprotonation of an unidentified base, B-, which is neither E47 nor
K48. In this work, we show by resonance Raman spectroscopy that at basic pH a
high-spin Fe3+-OH species is formed at the active site. The presence of the HO-
ligand was directly associated with an absorption band maximum at 560 nm,
whereas upon protonation, the band shifts to 644 nm. With respect to our
previous work, B- can be identified with this high-spin Fe3+-OH species, which
upon protonation results in a water molecule at the active site. Implications
for the SOR catalytic cycle are proposed.


  Desulfoferrodoxin is a small protein found in sulfate-reducing bacteria that
contains two independent mononuclear iron centers, one ferric and one ferrous.
Expression of desulfoferrodoxin from Desulfoarculus baarsii has been reported
to functionally complement a superoxide dismutase deficient Escherichia coli
strain. To elucidate by which mechanism desulfoferrodoxin could substitute for
superoxide dismutase in E. coli, we have purified the recombinant protein and
studied its reactivity toward O-(2). Desulfoferrodoxin exhibited only a weak
superoxide dismutase activity (20 units mg(-1)) that could hardly account for
its antioxidant properties. UV-visible and electron paramagnetic resonance
spectroscopy studies revealed that the ferrous center of desulfoferrodoxin
could specifically and efficiently reduce O-(2), with a rate constant of 6-7 x
10(8) M(-1) s(-1). In addition, we showed that membrane and cytoplasmic E. coli
protein extracts, using NADH and NADPH as electron donors, could reduce the
O-(2) oxidized form of desulfoferrodoxin. Taken together, these results
strongly suggest that desulfoferrodoxin behaves as a superoxide reductase
enzyme and thus provide new insights into the biological mechanisms designed
for protection from oxidative stresses.


  Flavin reductase catalyzes the reduction of free flavins by NAD(P)H. As
isolated, Escherichia coli flavin reductase does not contain any flavin
prosthetic group but accommodates both the reduced pyridine nucleotide and the
flavin substrate in a ternary complex prior to oxidoreduction. The reduction of
riboflavin by NADPH catalyzed by flavin reductase has been studied by static
and rapid kinetics absorption spectroscopies. Static absorption spectroscopy
experiments revealed that, in the presence of riboflavin and reduced pyridine
nucleotide, flavin reductase stabilizes, although to a small extent, a
charge-transfer complex of NADP+ and reduced riboflavin. In addition, reduction
of riboflavin was found to be essentially irreversible. Rapid kinetics
absorption spectroscopy studies demonstrated the occurrence of two
intermediates with long-wavelength absorption during the catalytic cycle. Such
intermediate species exhibit spectroscopic properties similar to those of
charge-transfer complexes of oxidized flavin and NAD(P)H, and reduced flavin
and NAD(P)+, respectively, which have been identified as intermediates during
the reaction of flavoenzymes of the ferredoxin-NADP+ reductase family. Thus, a
minimal kinetic scheme for the reaction of flavin reductase with NADPH and
riboflavin can be proposed. After formation of the Michaelis complex of flavin
reductase with NADPH and riboflavin, a first intermediate, identified as a
charge-transfer complex of NADPH and riboflavin, is formed. It is followed by a
second charge-transfer intermediate of enzyme-bound NADP+ and reduced
riboflavin. The latter decays, yielding the Michaelis complex of flavin
reductase with NADP+ and reduced riboflavin, which then dissociates to complete
the reaction. These results support the initial hypothesis of a structural
similarity between flavin reductase and the enzymes of the ferredoxin-NADP+
reductase family and extend it at a functional level.


  Iron-peroxide intermediates are central in the reaction cycle of many
iron-containing biomolecules. We trapped iron(III)-(hydro)peroxo species in
crystals of superoxide reductase (SOR), a nonheme mononuclear iron enzyme that
scavenges superoxide radicals. X-ray diffraction data at 1.95 angstrom
resolution and Raman spectra recorded in crystallo revealed iron-(hydro)peroxo
intermediates with the (hydro)peroxo group bound end-on. The dynamic SOR active
site promotes the formation of transient hydrogen bond networks, which
presumably assist the cleavage of the iron-oxygen bond in order to release the
reaction product, hydrogen peroxide.


  Several independent observations have suggested that catastrophe transition
in microtubules is not a first-order process, as is usually assumed. Recent
{\it in vitro} observations by Gardner et al.[ M. K. Gardner et al., Cell
{\bf147}, 1092 (2011)] showed that microtubule catastrophe takes place via
multiple steps and the frequency increases with the age of the filament. Here,
we investigate, via numerical simulations and mathematical calculations, some
of the consequences of age dependence of catastrophe on the dynamics of
microtubules as a function of the aging rate, for two different models of
aging: exponential growth, but saturating asymptotically and purely linear
growth. The boundary demarcating the steady state and non-steady state regimes
in the dynamics is derived analytically in both cases. Numerical simulations,
supported by analytical calculations in the linear model, show that aging leads
to non-exponential length distributions in steady state. More importantly,
oscillations ensue in microtubule length and velocity. The regularity of
oscillations, as characterized by the negative dip in the autocorrelation
function, is reduced by increasing the frequency of rescue events. Our study
shows that age dependence of catastrophe could function as an intrinsic
mechanism to generate oscillatory dynamics in a microtubule population,
distinct from hitherto identified ones.


  Which reaction networks, when taken with mass-action kinetics, have the
capacity for multiple steady states? There is no complete answer to this
question, but over the last 40 years various criteria have been developed that
can answer this question in certain cases. This work surveys these
developments, with an emphasis on recent results that connect the capacity for
multistationarity of one network to that of another. In this latter setting, we
consider a network $N$ that is embedded in a larger network $G$, which means
that $N$ is obtained from $G$ by removing some subsets of chemical species and
reactions. This embedding relation is a significant generalization of the
subnetwork relation. For arbitrary networks, it is not true that if $N$ is
embedded in $G$, then the steady states of $N$ lift to $G$. Nonetheless, this
does hold for certain classes of networks; one such class is that of fully open
networks. This motivates the search for embedding-minimal multistationary
networks: those networks which admit multiple steady states but no proper,
embedded networks admit multiple steady states. We present results about such
minimal networks, including several new constructions of infinite families of
these networks.


  In computational biology and bioinformatics, the manner to understand
evolution processes within various related organisms paid a lot of attention
these last decades. However, accurate methodologies are still needed to
discover genes content evolution. In a previous work, two novel approaches
based on sequence similarities and genes features have been proposed. More
precisely, we proposed to use genes names, sequence similarities, or both,
insured either from NCBI or from DOGMA annotation tools. Dogma has the
advantage to be an up-to-date accurate automatic tool specifically designed for
chloroplasts, whereas NCBI possesses high quality human curated genes (together
with wrongly annotated ones). The key idea of the former proposal was to take
the best from these two tools. However, the first proposal was limited by name
variations and spelling errors on the NCBI side, leading to core trees of low
quality. In this paper, these flaws are fixed by improving the comparison of
NCBI and DOGMA results, and by relaxing constraints on gene names while adding
a stage of post-validation on gene sequences. The two stages of similarity
measures, on names and sequences, are thus proposed for sequence clustering.
This improves results that can be obtained using either NCBI or DOGMA alone.
Results obtained with this quality control test are further investigated and
compared with previously released ones, on both computational and biological
aspects, considering a set of 99 chloroplastic genomes.


  Motor-driven cargo transport is a complex phenomenon where multiple motor
proteins attached on to a cargo engage in pulling activity, often leading to
tug-of-war, displaying bidirectional motion. However, most mathematical and
computational models ignore the details of the motor-cargo interaction. Here,
we study a generic model in which N motors are elastically coupled to a cargo,
which itself is subjected to thermal noise in the cytoplasm and to an
additional external applied force. The motor-hopping rates are chosen to
satisfy detailed balance with respect to the energy of elastic stretching. With
these assumptions, an (N+1)-variable master equation is constructed for
dynamics of the motor-cargo complex. By expanding the hopping rates to linear
order in fluctuations in motor positions, we obtain a linear Fokker-Planck
equation. The deterministic equations governing the average quantities are
separated out and explicit analytical expressions are obtained for the mean
velocity and diffusion coefficient of the cargo. We also study the statistical
features of the force experienced by an individual motor and quantitatively
characterize the load-sharing among the cargo-bound motors. The mean cargo
velocity and the effective diffusion coefficient are found to be decreasing
functions of the stiffness. While increase in the number of motors N does not
increase the velocity substantially, it decreases the effective diffusion
coefficient which falls as 1/N asymptotically. We further show that the
cargo-bound motors share the force exerted on the cargo equally only in the
limit of vanishing elastic stiffness; as stiffness is increased, deviations
from equal load sharing are observed. Numerical simulations agree with our
analytical results where expected. Interestingly, we find in simulations that
the stall force of a cargo elastically coupled to motors is independent of the
stiffness.


  Because biological processes can make different loci have different
evolutionary histories, species tree estimation requires multiple loci from
across the genome. While many processes can result in discord between gene
trees and species trees, incomplete lineage sorting (ILS), modeled by the
multi-species coalescent, is considered to be a dominant cause for gene tree
heterogeneity. Coalescent-based methods have been developed to estimate species
trees, many of which operate by combining estimated gene trees, and so are
called summary methods. Because summary methods are generally fast, they have
become very popular techniques for estimating species trees from multiple loci.
However, recent studies have established that summary methods can have reduced
accuracy in the presence of gene tree estimation error, and also that many
biological datasets have substantial gene tree estimation error, so that
summary methods may not be highly accurate on biologically realistic
conditions. Mirarab et al. (Science 2014) presented the statistical binning
technique to improve gene tree estimation in multi-locus analyses, and showed
that it improved the accuracy of MP-EST, one of the most popular
coalescent-based summary methods. Statistical binning, which uses a simple
statistical test for combinability and then uses the larger sets of genes to
re-calculate gene trees, has good empirical performance, but using statistical
binning within a phylogenomics pipeline does not have the desirable property of
being statistically consistent. We show that weighting the recalculated gene
trees by the bin sizes makes statistical binning statistically consistent under
the multispecies coalescent, and maintains the good empirical performance.
Thus, "weighted statistical binning" enables highly accurate genome-scale
species tree estimation, and is also statistical consistent under the
multi-species coalescent model.


  Different pathogens spreading in the same host population often generate
complex co-circulation dynamics because of the many possible interactions
between the pathogens and the host immune system, the host life cycle, and the
space structure of the population. Here we focus on the competition between two
acute infections and we address the role of host mobility and cross-immunity in
shaping possible dominance/co-dominance regimes. Host mobility is modelled as a
network of traveling flows connecting nodes of a metapopulation, and the
two-pathogen dynamics is simulated with a stochastic mechanistic approach.
Results depict a complex scenario where, according to the relation among the
epidemiological parameters of the two pathogens, mobility can either be
non-influential for the competition dynamics or play a critical role in
selecting the dominant pathogen. The characterisation of the parameter space
can be explained in terms of the trade-off between pathogen's spreading
velocity and its ability to diffuse in a sparse environment. Variations in the
cross-immunity level induce a transition between presence and absence of
competition. The present study disentangles the role of the relevant biological
and ecological factors in the competition dynamics, and provides relevant
insights into the spatial ecology of infectious diseases.


  This paper proposes a novel approach for storing and retrieving massive DNA
sequences.. The method is based on a perceptual hash function, commonly used to
determine the similarity between digital images, that we adapted for DNA
sequences. Perceptual hash function presented here is based on a Discrete
Cosine Transform Sign Only (DCT-SO). Each nucleotide is encoded as a fixed gray
level intensity pixel and the hash is calculated from its significant frequency
characteristics. This results to a drastic data reduction between the sequence
and the perceptual hash. Unlike cryptographic hash functions, perceptual hashes
are not affected by "avalanche effect" and thus can be compared. The similarity
distance between two hashes is estimated with the Hamming Distance, which is
used to retrieve DNA sequences. Experiments that we conducted show that our
approach is relevant for storing massive DNA sequences, and retrieving them.


  The degrees of freedom of multi-compartment mathematical models for energy
metabolism of a neuron-astrocyte complex may offer a key to understand the
different ways in which the energetic needs of the brain are met. In this paper
we address the problem within a steady state framework and we use the
techniques of linear algebra to identify the degrees of freedom first in a
lumped model, then in its extension to a spatially distributed case. The
interpretation of the degrees of freedom in metabolic terms, more specifically
in terms of glucose and oxygen partitioning, is then leveraged to derive
constraints on the free parameters needed to guarantee that the model is
energetically feasible. We also demonstrate how the model can be used to
estimate the stoichiometric energy needs of the cells as well as the household
energy based on observed oxidative cerebral metabolic rate (CMR) of glucose,
and the glutamate cycling. Moreover, our analysis shows that in the lumped
model the direction of lactate dehydrogenase (LDH) in the cells can be deduced
from the glucose partitioning between the compartments. The extension of the
lumped model into a spatially distributed multi-compartment setting that
includes diffusion fluxes from capillary to tissue increases the number of
degrees of freedom, requiring the use of statistical sampling techniques. The
analysis of distributed model reveals that some of the conclusions, e.g.,
concerning the LDH activity and glucose partitioning, based on a spatially
lumped model may no longer hold.


  Terrestrial ecosystems are generally green and only a small part (<10%) of
the plant matter is consumed by herbivores annually,but the reason has been
unclear due to the lack of food web models for predicting the absolute
herbivore biomass in physical units. Here, I present a simple parameterized
mathematical food web model that can predict the biomass density of herbivores
h (kg protein/m3) and carnivores c from ecological factors such as the
nutritive values of plants np (kg protein/m3), herbivores nh, and carnivores
nc, searching efficiency (volume) of carnivores S (/day), eating efficiency
(speed) of herbivores eh (/day) and carnivores ec, respiratory decrease in
herbivore and carnivore biomasses, dh (/day) and dc, absorption efficiency of
herbivores and carnivores h (ratio) and c, and probabilities of carnivores
preying on herbivores or carnivores, Phc (ratio) and Pcc.The model predicts a
stable equilibrium with low herbivore biomass h sufficient to keep the world
green if the food web consists of the three trophic levels, plants, herbivores
and carnivores; intraguild predation of carnivores exists; np<nh,nc; S>>eh, and
Phc>Pcc >0,which are well-realized in above-ground terrestrial ecosystems where
plant-rich "green world" is common. The h and c calculated from our model
showed good agreement with those from empirical observations in forests, where
both h and c are ca. 100 mg (fresh biomass/m2), and in savannahs. The model
predicts that the nutritive values and digestibility of plants are positively
correlated with h and the intensity of herbivory, which theoretically explains
the out-door defensive effects of the anti-nutritive or quantitative defenses
(e.g., tannins, protease inhibitors) of plants, and predicts that c and c/h are
positively correlated with the relative growth rate of herbivores. The present
model introduced parameterized realities into food web theory.


  Background: For faithful chromosome segregation during cell division, correct
attachments must be established between sister chromosomes and microtubules
from opposite spindle poles through kinetochores (chromosome bi-orientation).
Incorrect attachments of kinetochore microtubules (kMTs) lead to chromosome
mis-segregation and aneuploidy, which is often associated with developmental
abnormalities such as Down syndrome and diseases including cancer. The
interaction between kinetochores and microtubules is highly dynamic with
frequent attachments and detachments. However, it remains unclear how
chromosome bi-orientation is achieved with such accuracy in such a dynamic
process.
  Results: To gain new insight into this essential process, we have developed a
simple mathematical model of kinetochore-microtubule interactions during cell
division in general, i.e. both mitosis and meiosis. Firstly, the model reveals
that the balance between attachment and detachment probabilities of kMTs is
crucial for correct chromosome bi-orientation. With the right balance,
incorrect attachments are resolved spontaneously into correct bi-oriented
conformations while an imbalance leads to persistent errors. In addition, the
model explains why errors are more commonly found in the first meiotic division
(meiosis I) than in mitosis and how a faulty conformation can evade the spindle
assembly checkpoint, which may lead to a chromosome loss.
  Conclusions: The proposed model, despite its simplicity, helps us understand
one of the primary causes of chromosomal instability - aberrant
kinetochore-microtubule interactions. The model reveals that chromosome
bi-orientation is a probabilistic self-organization, rather than a
sophisticated process of error detection and correction.


  We study the response of an autoregulated gene to a range of concentrations
of signal molecules. We show that transcriptional leakage and noise due to
translational bursting have the opposite effects. In a positively autoregulated
gene, increasing the noise converts the response from graded to binary, while
increasing the leakage converts the response from binary to graded. Our
findings support the hypothesis that, being a common phenomenon, leaky
expression may be a relatively easy way for evolutionary tuning of the type of
gene response without changing the type of regulation from positive to
negative.


  Motivation: Assigning RNA-seq reads to their transcript of origin is a
fundamental task in transcript expression estimation. Where ambiguities in
assignments exist due to transcripts sharing sequence, e.g. alternative
isoforms or alleles, the problem can be solved through probabilistic inference.
Bayesian methods have been shown to provide accurate transcript abundance
estimates compared to competing methods. However, exact Bayesian inference is
intractable and approximate methods such as Markov chain Monte Carlo (MCMC) and
Variational Bayes (VB) are typically used. While providing a high degree of
accuracy and modelling flexibility, standard implementations can be
prohibitively slow for large datasets and complex transcriptome annotations.
  Results: We propose a novel approximate inference scheme based on VB and
apply it to an existing model of transcript expression inference from RNA-seq
data. Recent advances in VB algorithmics are used to improve the convergence of
the algorithm beyond the standard Variational Bayes Expectation Maximisation
(VBEM) algorithm. We apply our algorithm to simulated and biological datasets,
demonstrating a significant increase in speed with only very small loss in
accuracy of expression level estimation. We carry out a comparative study
against seven popular alternative methods and demonstrate that our new
algorithm provides excellent accuracy and inter-replicate consistency while
remaining competitive in computation time.
  Availability: The methods were implemented in R and C++, and are available as
part of the BitSeq project at \url{https://github.com/BitSeq}. The method is
also available through the BitSeq Bioconductor package. The source code to
reproduce all simulation results can be accessed via
\url{https://github.com/BitSeq/BitSeqVB_benchmarking}.


  In evolutionary game theory, evolutionarily stable states are characterised
by the folk theorem because exact solutions to the replicator equation are
difficult to obtain. It is generally assumed that the folk theorem, which is
the fundamental theory for non-cooperative games, defines all Nash equilibria
in infinitely repeated games. Here, we prove that Nash equilibria that are not
characterised by the folk theorem do exist. By adopting specific reactive
strategies, a group of players can be better off by coordinating their actions
in repeated games. We call it a type-k equilibrium when a group of k players
coordinate their actions and they have no incentive to deviate from their
strategies simultaneously. The existence and stability of the type-k
equilibrium in general games is discussed. This study shows that the sets of
Nash equilibria and evolutionarily stable states have greater cardinality than
classic game theory has predicted in many repeated games.


  Cytoskeletal networks, which are essentially motor-filament assemblies, play
a major role in many developmental processes involving structural remodeling
and shape changes. These are achieved by nonequilibrium self-organization
processes that generate functional patterns and drive intracellular transport.
We construct a minimal physical model that incorporates the coupling between
nonlinear elastic responses of individual filaments and force-dependent motor
action. By performing stochastic simulations we show that the interplay of
motor processes, described as driving anti-correlated motion of the network
vertices, and the network connectivity, which determines the percolation
character of the structure, can indeed capture the dynamical and structural
cooperativity which gives rise to diverse patterns observed experimentally. The
buckling instability of individual filaments is found to play a key role in
localizing collapse events due to local force imbalance. Motor-driven
buckling-induced node aggregation provides a dynamic mechanism that stabilizes
the two dimensional patterns below the apparent static percolation limit.
Coordinated motor action is also shown to suppress random thermal noise on
large time scales, the two dimensional configuration that the system starts
with thus remaining planar during the structural development. By carrying out
similar simulations on a three dimensional anchored network, we find that the
myosin-driven isotropic contraction of a well-connected actin network, when
combined with mechanical anchoring that confers directionality to the
collective motion, may represent a novel mechanism of intracellular transport,
as revealed by chromosome translocation in the starfish oocyte.


  We investigated the viscoelastic response of model interphase chromosomes by
tracking the three-dimensional motion of hundreds of dispersed Brownian
particles of sizes ranging from the thickness of the chromatin fiber up to
slightly above the mesh size of the chromatin solution. In agreement with
previous computational studies on polymer solutions and melts, we found that
the large-time behaviour of diffusion coefficient and the experienced viscosity
of moving particles as functions of particle size deviate from the traditional
Stokes-Einstein relation, and agree with a recent scaling theory of diffusion
of non-sticky particles in polymer solutions. Interestingly, we found that at
short times large particles are temporary "caged" by chromatin spatial
constraints, which thus form effective domains whose size match remarkably well
recent experimental results for micro-tracers inside interphase nuclei.
Finally, by employing a known mathematical relation between the time
mean-square displacement of tracked particles and the complex shear modulus of
the surrounding solution, we calculated the elastic and viscous moduli of
interphase chromosomes.


  House of Aisin Gioro is the imperial family of the last dynasty in Chinese
history - Qing Dynasty (1644 - 1911). Aisin Gioro family originated from
Jurchen tribes and developed the Manchu people before they conquered China. By
investigating the Y chromosomal short tandem repeats (STRs) of 7 modern male
individuals who claim belonging to Aisin Gioro family (in which 3 have full
records of pedigree), we found that 3 of them (in which 2 keep full pedigree,
whose most recent common ancestor is Nurgaci) shows very close relationship (1
- 2 steps of difference in 17 STR) and the haplotype is rare. We therefore
conclude that this haplotype is the Y chromosome of the House of Aisin Gioro.
Further tests of single nucleotide polymorphisms (SNPs) indicates that they
belong to Haplogroup C3b2b1*-M401(xF5483), although their Y-STR results are
distant to the "star cluster", which also belongs to the same haplogroup. This
study forms the base for the pedigree research of the imperial family of Qing
Dynasty by means of genetics.


  Neutral theories have played a crucial and revolutionary role in fields such
as population genetics and biogeography. These theories are critical by
definition, in the sense that the overall growth rate of each single
allele/species/type vanishes. Thus each species in a neutral model sits at the
edge between invasion and extinction, allowing for the coexistence of
symmetric/neutral types. However, in finite systems, mono-dominated states are
ineludibly reached in relatively short times owing to demographic fluctuations,
thus leaving us with an unsatisfactory framework to rationalize
empirically-observed long-term coexistence. Here, we scrutinize the effect of
heterogeneity in quasi-neutral theories, in which there can be a local mild
preference for some of the competing species at some sites, even if the overall
species symmetry is maintained. As we show here, mild biases at a small
fraction of locations suffice to induce overall robust and durable species
coexistence, even in regions arbitrarily far apart from the biased locations.
This result stems from the long-range nature of the underlying critical bulk
dynamics and has a number of implications, for example, in conservation ecology
as it suggests that constructing local specific "sanctuaries" for different
competing species can result in global enhancement of biodiversity, even in
regions arbitrarily distant from the protected refuges.


  Advances in DNA nanotechnology have stimulated the search for simple motifs
that can be used to control the properties of DNA nanostructures. One such
motif, which has been used extensively in structures such as polyhedral cages,
two-dimensional arrays, and ribbons, is a bulged duplex, that is two helical
segments that connect at a bulge loop. We use a coarse-grained model of DNA to
characterize such bulged duplexes. We find that this motif can adopt structures
belonging to two main classes: one where the stacking of the helices at the
center of the system is preserved, the geometry is roughly straight and the
bulge is on one side of the duplex, and the other where the stacking at the
center is broken, thus allowing this junction to act as a hinge and increasing
flexibility. Small loops favor states where stacking at the center of the
duplex is preserved, with loop bases either flipped out or incorporated into
the duplex. Duplexes with longer loops show more of a tendency to unstack at
the bulge and adopt an open structure. The unstacking probability, however, is
highest for loops of intermediate lengths, when the rigidity of single-stranded
DNA is significant and the loop resists compression. The properties of this
basic structural motif clearly correlate with the structural behavior of
certain nano-scale objects, where the enhanced flexibility associated with
larger bulges has been used to tune the self-assembly product as well as the
detailed geometry of the resulting nanostructures.


  The automated inference of physically interpretable (bio)chemical reaction
network models from measured experimental data is a challenging problem whose
solution has significant commercial and academic ramifications. It is
demonstrated, using simulations, how sets of elementary reactions comprising
chemical reaction networks, as well as their rate coefficients, may be
accurately recovered from non-equilibrium time series concentration data, such
as that obtained from laboratory scale reactors. A variant of an evolutionary
algorithm called differential evolution in conjunction with least squares
techniques is used to search the space of reaction networks in order to infer
both the reaction network topology and its rate parameters. Properties of the
stoichiometric matrices of trial networks are used to bias the search towards
physically realisable solutions. No other information, such as chemical
characterisation of the reactive species is required, although where available
it may be used to improve the search process.


  Systems biology is an inter-disciplinary field that studies systems of
biological components at different scales, which may be molecules, cells or
entire organism. In particular, systems biology methods are applied to
understand functional deregulations within human cells (e.g., cancers). In this
context, we present several python packages linked to CellNOptR (R package),
which is used to build predictive logic models of signalling networks by
training networks (derived from literature) to signalling (phospho-proteomic)
data. The first package (cellnopt.wrapper) is a wrapper based on RPY2 that
allows a full access to CellNOptR functionalities within Python. The second one
(cellnopt.core) was designed to ease the manipulation and visualisation of data
structures used in CellNOptR, which was achieved by using Pandas, NetworkX and
matplotlib. Systems biology also makes extensive use of web resources and
services. We will give an overview and status of BioServices, which allows one
to access programmatically to web resources used in life science and how it can
be combined with CellNOptR.


  Single molecule F\"orster resonance energy transfer (smFRET) is a powerful
experimental technique for studying the properties of individual biological
molecules in solution. However, as adoption of smFRET techniques becomes more
widespread, the lack of available software, whether open source or commercial,
for data analysis, is becoming a significant issue. Here, we present pyFRET, an
open source Python package for the analysis of data from single-molecule
fluorescence experiments from freely diffusing biomolecules. The package
provides methods for the complete analysis of a smFRET dataset, from burst
selection and denoising, through data visualisation and model fitting. We
provide support for both continuous excitation and alternating laser excitation
(ALEX) data analysis. pyFRET is available as a package downloadable from the
Python Package Index (PyPI) under the open source three-clause BSD licence,
together with links to extensive documentation and tutorials, including example
usage and test data. Additional documentation including tutorials is hosted
independently on ReadTheDocs. The code is available from the free hosting site
Bitbucket. Through distribution of this software, we hope to lower the barrier
for the adoption of smFRET experiments by other research groups and we
encourage others to contribute modules for specific analysis needs.


  We construct an individual-based metapopulation model of population genetics
featuring migration, mutation, selection and genetic drift. In the case of a
single `island', the model reduces to the Moran model. Using the diffusion
approximation and timescale separation arguments, an effective one-variable
description of the model is developed. The effective description bears
similarities to the well-mixed Moran model with effective parameters which
depend on the network structure and island sizes, and is amenable to analysis.
Predictions from the reduced theory match the results from stochastic
simulations across a range of parameters. The nature of the fast-variable
elimination technique we adopt is further studied by applying it to a linear
system, where it provides a precise description of the slow-dynamics in the
limit of large timescale separation.


  We introduce a formulation for normal mode analyses of globular proteins that
significantly improves on an earlier, 1-parameter formulation (M. Tirion, PRL
77, 1905 (1996)) that characterized the slow modes associated with protein data
bank structures. Here we develop that empirical potential function which is
minimized at the outset to include two features essential to reproduce the
eigenspectra and associated density of states over all frequencies, not merely
the slow ones. First, introduction of preferred dihedral-angle configurations
via use of torsional stiffness constants eliminates anomalous dispersion
characteristics due to insufficiently bound surface sidechains. Second, we take
into account the atomic identities and the distance of separation of all
pairwise interactions. With these modifications we obtain stable, reliable
eigenmodes over a wide range of frequencies.


  Game theory ideas provide a useful framework for studying evolutionary
dynamics in a well-mixed environment. This approach, however, typically
enforces a strictly fixed overall population size, deemphasizing natural growth
processes. We study a competitive Lotka-Volterra model, with number
fluctuations, that accounts for natural population growth and encompasses
interaction scenarios typical of evolutionary games. We show that, in an
appropriate limit, the model describes standard evolutionary games with both
genetic drift and overall population size fluctuations. However, there are also
regimes where a varying population size can strongly influence the evolutionary
dynamics. We focus on the strong mutualism scenario and demonstrate that
standard evolutionary game theory fails to describe our simulation results. We
then analytically and numerically determine fixation probabilities as well as
mean fixation times using matched asymptotic expansions, taking into account
the population size degree of freedom. These results elucidate the interplay
between population dynamics and evolutionary dynamics in well-mixed systems.


  This article evaluates bias in one class of methods used to estimate archaic
admixture in modern humans. These methods study the pattern of allele sharing
among modern and archaic genomes. They are sensitive to "ghost" admixture,
which occurs when a population receives archaic DNA from sources not
acknowledged by the statistical model. The effect of ghost admixture depends on
two factors: branch-length bias and population-size bias. Branch-length bias
occurs because a given amount of admixture has a larger effect if the two
populations have been separated for a long time. Population-size bias occurs
because differences in population size distort branch lengths in the gene
genealogy. In the absence of ghost admixture, these effects are small. They
become important, however, in the presence of ghost admixture. Estimators
differ in the pattern of response. Increasing a given parameter may inflate one
estimator but deflate another. For this reason, comparisons among estimators
are informative. Using such comparisons, this article supports previous
findings that the archaic population was small and that Europeans received
little gene flow from archaic populations other than Neanderthals. It also
identifies an inconsistency in estimates of archaic admixture into Melanesia.


  Heritability is a central parameter in quantitative genetics, both from an
evolutionary and a breeding perspective. For plant traits heritability is
traditionally estimated by comparing within and between genotype variability.
This approach estimates broad-sense heritability, and does not account for
different genetic relatedness. With the availability of high-density markers
there is growing interest in marker based estimates of narrow-sense
heritability, using mixed models in which genetic relatedness is estimated from
genetic markers. Such estimates have received much attention in human genetics
but are rarely reported for plant traits. A major obstacle is that current
methodology and software assume a single phenotypic value per genotype, hence
requiring genotypic means. An alternative that we propose here, is to use mixed
models at individual plant or plot level. Using statistical arguments,
simulations and real data we investigate the feasibility of both approaches,
and how these affect genomic prediction with G-BLUP and genome-wide association
studies. Heritability estimates obtained from genotypic means had very large
standard errors and were sometimes biologically unrealistic. Mixed models at
individual plant or plot level produced more realistic estimates, and for
simulated traits standard errors were up to 13 times smaller. Genomic
prediction was also improved by using these mixed models, with up to a 49%
increase in accuracy. For GWAS on simulated traits, the use of individual plant
data gave almost no increase in power. The new methodology is applicable to any
complex trait where multiple replicates of individual genotypes can be scored.
This includes important agronomic crops, as well as bacteria and fungi.


  Collective behaviour in biological systems is often accompanied by strong
correlations. The question has therefore arisen of whether correlation is
amplified by the vicinity to some critical point in the parameters space.
Biological systems, though, are typically quite far from the thermodynamic
limit, so that the value of the control parameter at which correlation and
susceptibility peak depend on size. Hence, a system would need to readjust its
control parameter according to its size in order to be maximally correlated.
This readjustment, though, has never been observed experimentally. By gathering
three-dimensional data on swarms of midges in the field we find that swarms
tune their control parameter and size so as to maintain a scaling behaviour of
the correlation function. As a consequence, correlation length and
susceptibility scale with the system's size and swarms exhibit a near-maximal
degree of correlation at all sizes.


  Species richness varies widely across the tree of life, and there is great
interest in identifying ecological, geographic, and other factors that affect
rates of species proliferation. Recent methods for explicitly modeling the
relationships among character states, speciation rates, and extinction rates on
phylogenetic trees- BiSSE, QuaSSE, GeoSSE, and related models - have been
widely used to test hypotheses about character state-dependent diversification
rates. Here, we document the disconcerting ease with which neutral traits are
inferred to have statistically significant associations with speciation rate.
We first demonstrate this unfortunate effect for a known model assumption
violation: shifts in speciation rate associated with a character not included
in the model. We further show that for many empirical phylogenies, characters
simulated in the absence of state-dependent diversification exhibit an even
higher Type I error rate, indicating that the method is susceptible to
additional, unknown model inadequacies. For traits that evolve slowly, the root
cause appears to be a statistical framework that does not require replicated
shifts in character state and diversification. However, spurious associations
between character state and speciation rate arise even for traits that lack
phylogenetic signal, suggesting that phylogenetic pseudoreplication alone
cannot fully explain the problem. The surprising severity of this phenomenon
suggests that many trait-diversification relationships reported in the
literature may not be real. More generally, we highlight the need for
diagnosing and understanding the consequences of model inadequacy in
phylogenetic comparative methods.


  Standard game theory cannot describe microbial interactions mediated by
diffusible molecules. Nevertheless, we show that one can still model microbial
dynamics using game theory with parameters renormalized by diffusion. Contrary
to expectations, greater sharing of metabolites reduces the strength of
cooperation and leads to species extinction via a nonequilibrium phase
transition. We report analytic results for the critical diffusivity and the
length scale of species intermixing. We also show that fitness nonlinearities
suppress mutualism and favor the species producing slower nutrients.


  Spectra of real world networks exhibit properties which are different from
the random networks. One such property is the existence of a very high
degeneracy at zero eigenvalues. In this work, we provide possible reasons
behind occurrence of the zero degeneracy in various networks spectra.
Comparison of zero degeneracy in protein-protein interaction networks of six
different species and in their corresponding model networks sheds light in
understanding the evolution of complex biological systems.


  The human habitat is a host where microbial species evolve, function, and
continue to evolve. Elucidating how microbial communities respond to human
habitats is a fundamental and critical task, as establishing baselines of human
microbiome is essential in understanding its role in human disease and health.
However, current studies usually overlook a complex and interconnected
landscape of human microbiome and limit the ability in particular body habitats
with learning models of specific criterion. Therefore, these methods could not
capture the real-world underlying microbial patterns effectively. To obtain a
comprehensive view, we propose a novel ensemble clustering framework to mine
the structure of microbial community pattern on large-scale metagenomic data.
Particularly, we first build a microbial similarity network via integrating
1920 metagenomic samples from three body habitats of healthy adults. Then a
novel symmetric Nonnegative Matrix Factorization (NMF) based ensemble model is
proposed and applied onto the network to detect clustering pattern. Extensive
experiments are conducted to evaluate the effectiveness of our model on
deriving microbial community with respect to body habitat and host gender. From
clustering results, we observed that body habitat exhibits a strong bound but
non-unique microbial structural patterns. Meanwhile, human microbiome reveals
different degree of structural variations over body habitat and host gender. In
summary, our ensemble clustering framework could efficiently explore integrated
clustering results to accurately identify microbial communities, and provide a
comprehensive view for a set of microbial communities. Such trends depict an
integrated biography of microbial communities, which offer a new insight
towards uncovering pathogenic model of human microbiome.


  In computational biology, biological entities such as genes or proteins are
usually annotated with terms extracted from Gene Ontology (GO). The functional
similarity among terms of an ontology is evaluated by using Semantic Similarity
Measures (SSM). More recently, the extensive application of SSMs yielded to the
Semantic Similarity Networks (SSNs). SSNs are edge-weighted graphs where the
nodes are concepts (e.g. proteins) and each edge has an associated weight that
represents the semantic similarity among related pairs of nodes. The analysis
of SSNs may reveal biologically meaningful knowledge. For these aims, the need
for the introduction of tool able to manage and analyze SSN arises.
Consequently we developed SSN-Analyzer a web based tool able to build and
preprocess SSN. As proof of concept we demonstrate that community detection
algorithms applied to filtered (thresholded) networks, have better performances
in terms of biological relevance of the results, with respect to the use of raw
unfiltered networks.


  In shallow macrotidal ecosystems with large intertidal mudflats, the
sediment-water coupling plays a crucial role in structuring the pelagic
microbial food web functioning, since inorganic and organic matter and
microbial components (viruses and microbes) of the microphytobenthic biofilm
can be suspended toward the water column. Two experimental bioassays were
conducted in March and July 2008 to investigate the importance of biofilm input
for the pelagic Microbial and Viral Loops. Pelagic inocula (<0.6$\mu$ and
<10$\mu$ filtrates) were diluted either with \textless{}30kDa-ultrafiltered
seawater or with this ultrafiltrate enriched with the respective
size-fractionated benthic biofilm or with \textless{}30kDa-benthic compounds
(BC). The kinetics of heterotrophic nanoflagellates (HNF), bacteria and viruses
were assessed together with bacterial and viral genomic fingerprints, bacterial
enzymatic activities and viral life strategies. The experimental design allowed
us to evaluate the effect of BC modulated by those of benthic size-fractionated
microorganisms (virus+bacteria, +HNF). BC presented (1) in March, a positive
effect on viruses and bacteria weakened by pelagic HNF. Benthic microorganisms
consolidated this negative effect and sustained the viral production together
with a relatively diverse and uneven bacterial assemblage structure; (2) in
July, no direct impact on viruses but a positive effect on bacteria modulated
by HNF, which indirectly enhanced viral multiplication. Both effects were
intensified by benthic microorganisms and bacterial assemblage structure became
more even. HNF indirectly profited from BC more in March than in July. The
Microbial Loop would be stimulated by biofilm during periods of high resources
(March) and the Viral Loop during periods of depleted resources (July).


  Approximate Bayesian Computation (ABC) is a family of statistical inference
techniques, which is increasingly used in biology and other scientific fields.
Its main benefit is to be applicable to models for which the computation of the
model likelihood is intractable. The basic idea of ABC is to empirically
approximate the model likelihood by using intensive realizations of model runs.
Due to computing time limitations, ABC has thus been mainly applied to models
that are relatively quick to simulate. We here aim at briefly introducing the
field of statistical emulation of computer code outputs and to demonstrate its
potential for ABC applications. Emulation consists in replacing the costly to
simulate model by another (quick to simulate) statistical model called emulator
or meta-model. This emulator is fitted to a small number of outputs of the
original model, and is subsequently used as a surrogate during the inference
procedure. In this contribution, we first detail the principles of model
emulation, with a special reference to the ABC context in which the description
of the stochasticity of model realizations is as important as the description
of the trends linking model parameters and outputs. We then compare several
emulation strategies in an ABC context, using as case study a stochastic
ecological model of community dynamics. We finally describe a novel
emulation-based sequential ABC algorithm which is shown to decrease computing
time by a factor of two on the studied example, compared to previous sequential
ABC algorithms. Routines to perform emulation-based ABC were made available
within the R package EasyABC.


  Persistent homology has emerged as a popular technique for the topological
simplification of big data, including biomolecular data. Multidimensional
persistence bears considerable promise to bridge the gap between geometry and
topology. However, its practical and robust construction has been a challenge.
We introduce two families of multidimensional persistence, namely
pseudo-multidimensional persistence and multiscale multidimensional
persistence. The former is generated via the repeated applications of
persistent homology filtration to high dimensional data, such as results from
molecular dynamics or partial differential equations. The latter is constructed
via isotropic and anisotropic scales that create new simiplicial complexes and
associated topological spaces. The utility, robustness and efficiency of the
proposed topological methods are demonstrated via protein folding, protein
flexibility analysis, the topological denoising of cryo-electron microscopy
data, and the scale dependence of nano particles. Topological transition
between partial folded and unfolded proteins has been observed in
multidimensional persistence. The separation between noise topological
signatures and molecular topological fingerprints is achieved by the
Laplace-Beltrami flow. The multiscale multidimensional persistent homology
reveals relative local features in Betti-0 invariants and the relatively global
characteristics of Betti-1 and Betti-2 invariants.


  Allosteric communication in proteins is a central and yet unsolved problem of
structural biochemistry. Previous findings, from computational biology (Ota and
Agard, 2005), have proposed that heat diffuses in a protein through cognate
protein allosteric pathways. This work studied heat diffusion in the well-known
PDZ-2 protein, and confirmed that this protein has two cognate allosteric
pathways and that heat flows preferentially through these. Also, a new property
was also observed for protein structures - heat diffuses asymmetrically through
the structures. The underling structure of this asymmetrical heat flow was a
normal length hydrogen bond (~2.85 {\AA}) that acted as a thermal rectifier. In
contrast, thermal rectification was compromised in short hydrogen bonds (~2.60
{\AA}), giving rise to symmetrical thermal diffusion. Asymmetrical heat
diffusion was due, on a higher scale, to the local, structural organization of
residues that, in turn, was also mediated by hydrogen bonds. This
asymmetrical/symmetrical energy flow may be relevant for allosteric signal
communication directionality in proteins and for the control of heat flow in
materials science.


  Quorum sensing is the regulation of gene expression in response to changes in
cell density. To measure their cell density, bacterial populations produce and
detect diffusible molecules called autoinducers. Individual bacteria internally
represent the external concentration of autoinducers via the level of monitor
proteins. In turn, these monitor proteins typically regulate both their own
production and the production of autoinducers, thereby establishing internal
and external feedbacks. Here, we ask whether feedbacks can increase the
information available to cells about their local density. We quantify available
information as the mutual information between the abundance of a monitor
protein and the local cell density for biologically relevant models of quorum
sensing. Using variational methods, we demonstrate that feedbacks can increase
information transmission, allowing bacteria to resolve up to two additional
ranges of cell density when compared with bistable quorum-sensing systems. Our
analysis is relevant to multi-agent systems that track an external driver
implicitly via an endogenously generated signal.


  The human islet amyloid polypeptide (hIAPP) co-operates with insulin to
maintain glycemic balance. It also constitutes the amyloid plaques that
aggregate in the pancreas of type-II diabetic patients. We have performed
extensive in silico investigations to analyse the structural landscape of
monomeric hIAPP, which is presumed to be intrinsically disordered. For this we
construct from first principles a highly predictive energy function that
describes a monomeric hIAPP observed in a NMR experiment, as a local energy
minimum. We subject our theoretical model of hIAPP to repeated heating and
cooling simulations, back and forth between a high temperature regime where the
conformation resembles a random walker and a low temperature limit where no
thermal motions prevail. We find that the final low temperature conformations
display a high level of degeneracy, in a manner which is fully in line with the
presumed intrinsically disordered character of hIAPP. In particular, we
identify an isolated family of alpha-helical conformations that might cause the
transition to amyloidosis, by nucleation.


  A commonly recurring problem in structural protein studies, is the
determination of all heavy atom positions from the knowledge of the central
alpha-carbon coordinates. We employ advances in virtual reality to address the
problem. The outcome is a 3D visualisation based technique where all the heavy
backbone and side chain atoms are treated on equal footing, in terms of the
C-alpha coordinates. Each heavy atom can be visualised on the surfaces of the
different two-spheres, that are centered at the other heavy backbone and side
chain atoms. In particular, the rotamers are visible as clusters which display
strong dependence on the underlying backbone secondary structure. Our method
easily detects those atoms in a crystallographic protein structure which have
been been likely misplaced. Our approach forms a basis for the development of a
new generation, visualisation based side chain construction, validation and
refinement tools. The heavy atom positions are identified in a manner which
accounts for the secondary structure environment, leading to improved accuracy
over existing methods.


  The connection between optimal biological function and energy use, measured
for example by the rate of metabolite consumption, is a current topic of
interest in the systems biology literature which has been explored in several
different contexts. In [J. P. Barton and E. D. Sontag, Biophys. J. 104, 6
(2013)], we related the metabolic cost of enzymatic futile cycles with their
capacity to act as insulators which facilitate modular interconnections in
biochemical networks. There we analyzed a simple model system in which a signal
molecule regulates the transcription of one or more target proteins by
interacting with their promoters. In this note, we consider the case of a
protein with an active and an inactive form, and whose activation is controlled
by the signal molecule. As in the original case, higher rates of energy
consumption are required for better insulator performance.


  Model uncertainty and limited data are fundamental challenges to robust
management of human intervention in a natural system. These challenges are
acutely highlighted by concerns that many ecological systems may contain
tipping points, such as Allee population sizes. Before a collapse, we do not
know where the tipping points lie, if they exist at all. Hence, we know neither
a complete model of the system dynamics nor do we have access to data in some
large region of state-space where such a tipping point might exist. We
illustrate how a Bayesian Non-Parametric (BNP) approach using a Gaussian
Process (GP) prior provides a flexible representation of this inherent
uncertainty. We embed GPs in a Stochastic Dynamic Programming (SDP) framework
in order to make robust management predictions with both model uncertainty and
limited data. We use simulations to evaluate this approach as compared with the
standard approach of using model selection to choose from a set of candidate
models. We find that model selection erroneously favors models without tipping
points -- leading to harvest policies that guarantee extinction. The GPDP
performs nearly as well as the true model and significantly outperforms
standard approaches. We illustrate this using examples of simulated
single-species dynamics, where the standard model selection approach should be
most effective, and find that it still fails to account for uncertainty
appropriately and leads to population crashes, while management based on the
GPDP does not, since it does not underestimate the uncertainty outside of the
observed data.


  In this paper we report on an algorithm for aligning multiple protein
structures. The algorithm has been tested on a variety of inputs and it
performs well in comparison to well-known algorithms for this problem.


  The paper presents a geometrical model for protein secondary structure
analysis which uses only the positions of the $C_{\alpha}$-atoms. We construct
a space curve connecting these positions by piecewise polynomial interpolation
and describe the folding of the protein backbone by a succession of screw
motions linking the Frenet frames at consecutive $C_{\alpha}$-positions. Using
the ASTRAL subset of the SCOPe data base of protein structures, we derive
thresholds for the screw parameters of secondary structure elements and
demonstrate that the latter can be reliably assigned on the basis of a
$C_{\alpha}$-model. For this purpose we perform a comparative study with the
widely used DSSP (Define Secondary Structure of Proteins) algorithm.


  This is a set of lectures that I presented at the Les Houches 2014 Summer
School "Topological Aspects in Condensed Matter Physics". The lectures are an
introduction to physics of proteins. To physicists, and by a physicist. My
lectures at les Houches were also celebration of the anniversary of
Schroedinger's 1944 lectures, and for that reason I decided to share my title
with his book.


  A general class of stochastic gene expression models with self regulation is
considered. One or more genes randomly switch between regulatory states, each
having a different mRNA transcription rate. The gene or genes are self
regulating when the proteins they produce affect the rate of switching between
regulatory states. Under weak noise conditions, the deterministic forces are
much stronger than fluctuations from gene switching and protein synthesis.
Metastable transitions, such as bistable switching, can occur under weak noise
conditions, causing dramatic shifts in the expression of a gene. A general tool
used to describe metastability is the quasi stationary analysis (QSA). A large
deviation principle is derived so that the QSA can explicitly account for
random gene switching without using an adiabatic limit or diffusion
approximation, which are unreliable and inaccurate for metastable events.This
allows the existing asymptotic and numerical methods that have been developed
for continuous Markov processes to be used to analyze the full model.


  Somatic variants can be used as lineage markers for the phylogenetic
reconstruction of cancer evolution. Since somatic phylogenetics is complicated
by sample heterogeneity, novel specialized tree-building methods are required
for cancer phylogeny reconstruction. We present LICHeE (Lineage Inference for
Cancer Heterogeneity and Evolution), a novel method that automates the
phylogenetic inference of cancer progression from multiple somatic samples.
LICHeE uses variant allele frequencies of SSNVs obtained by deep sequencing to
reconstruct multi-sample cell lineage trees and infer the subclonal composition
of the samples. LICHeE is open-sourced and available at
http://viq854.github.io/lichee.


  Gene expression in individual cells is highly variable and sporadic, often
resulting in the synthesis of mRNAs and proteins in bursts. Bursting in gene
expression is known to impact cell-fate in diverse systems ranging from latency
in HIV-1 viral infections to cellular differentiation. It is generally assumed
that bursts are geometrically distributed and that they arrive according to a
Poisson process. On the other hand, recent single-cell experiments provide
evidence for complex burst arrival processes, highlighting the need for more
general stochastic models. To address this issue, we invoke a mapping between
general models of gene expression and systems studied in queueing theory to
derive exact analytical expressions for the moments associated with
mRNA/protein steady-state distributions. These moments are then used to derive
explicit conditions, based entirely on experimentally measurable quantities,
that determine if the burst distributions deviate from the geometric
distribution or if burst arrival deviates from a Poisson process. For
non-Poisson arrivals, we develop approaches for accurate estimation of burst
parameters.


  When birds come together to form a flock, the distribution of their
individual velocities narrows around the mean velocity of the flock. We argue
that, in a broad class of models for the joint distribution of positions and
velocities, this narrowing generates an entropic force that opposes the
cohesion of the flock. The strength of this force depends strongly on the
nature of the interactions among birds: if birds are coupled to a fixed number
of neighbors, the entropic forces are weak, while if they couple to all other
birds within a fixed distance, the entropic forces are sufficient to tear a
flock apart. Similar entropic forces should occur in other non-equilibrium
systems. For the joint distribution of protein structures and amino-acid
sequences, these forces favor the occurrence of "highly designable" structures.


  This article shows how 3D geometric morphometric data can be analyzed using
newly developed distance-based evolutionary tree inference methods, with
extensions to planar graphs. Application of these methods to 3D representations
of the skullcap (calvaria) of 13 diverse skulls in the genus Homo, ranging from
Homo erectus (ergaster) at about 1.6 mya, all the way forward to modern humans,
yields a remarkably clear phylogenetic tree. Various evolutionary hypotheses
are tested. Results of these tests include rejection of the monophyly of Homo
heidelbergensis, the Multi-Regional hypothesis, and the hypothesis that the
unusual 12,000 year old (12kya) Iwo Eleru skull represents a modern human.
Rather, by quantitative phylogenetic analyses the latter is seen to be an old
(200-400kya) lineage that probably represents a novel African species, Homo
iwoelerueensis. It diverged after the lineage leading to Neanderthals, and may
have been driven to extinction in the last 10kya by modern humans, Homo
sapiens, another African species of Homo that appeared about 100kya. Another
enigmatic skull, Qafzeh 6 from the Middle East about 90kya, appears to be a
hybrid of two thirds near, but not, anatomically modern human and one third of
an archaic lineage diverging close to classic European Neanderthals. Overall,
the tree clearly implies an accelerating rate of skullcap shape change, and by
extension, change of the underlying brain, over the last 400kya in Africa. This
acceleration may have extended right up to the origin of modern humans. Methods
of distance-based evolutionary tree inference are refined and extended, with
particular attention to diagnosing the model and achieving a better fit. This
includes power transformations of the input data which favor root Procrustes
distances.


  RNA molecules are essential cellular machines performing a wide variety of
functions for which a specific three-dimensional structure is required. Over
the last several years, experimental determination of RNA structures through
X-ray crystallography and NMR seems to have reached a plateau in the number of
structures resolved each year, but as more and more RNA sequences are being
discovered, need for structure prediction tools to complement experimental data
is strong. Theoretical approaches to RNA folding have been developed since the
late nineties when the first algorithms for secondary structure prediction
appeared. Over the last 10 years a number of prediction methods for 3D
structures have been developed, first based on bioinformatics and data-mining,
and more recently based on a coarse-grained physical representation of the
systems. In this review we are going to present the challenges of RNA structure
prediction and the main ideas behind bioinformatic approaches and physics-based
approaches. We will focus on the description of the more recent physics-based
phenomenological models and on how they are built to include the specificity of
the interactions of RNA bases, whose role is critical in folding. Through
examples from different models, we will point out the strengths of
physics-based approaches, which are able not only to predict equilibrium
structures, but also to investigate dynamical and thermodynamical behavior, and
the open challenges to include more key interactions ruling RNA folding.


  We use Monte-Carlo simulations to study the effect of variable rigidity on
plectoneme formation and localization in supercoiled dsDNA. We show that the
presence of soft sequences increases the number of plectoneme branches and that
the edges of the branches tend to be localized at these sequences. We propose
an experimental approach to test our results in vitro, and discuss the possible
role played by plectoneme localization in the search process of transcription
factors for their targets (promoter regions) on the bacterial genome.


  The effect of mechanical interactions between cells in the spreading of
bacterial populations was investigated in one-dimensional space. A
continuum-mechanics approach, comprising cell migration, proliferation, and
exclusion processes, was employed to elucidate the dynamics. The consequent
nonlinear reaction-diffusion-like equation describes the constitution dynamics
of a bacterial population. In this model, bacterial cells were treated as
rod-like particles that interact with each other through hard-core repulsion,
which introduces the exclusion effect that causes bacterial populations to
migrate quickly and at high density. The propagation of bacterial density as a
traveling wave front over extended times was also analysed. The analytical and
numerical solutions revealed that the front speed was enhanced by the exclusion
process, which depended upon the cell-packing fraction. Finally, we
qualitatively compared our theoretical results with experimental evidence.


  Motivation: Gene blocks are genes co-located on the chromosome. In many
cases, genes blocks are conserved between bacterial species, sometimes as
operons, when genes are co-transcribed. The conservation is rarely absolute:
gene loss, gain, duplication, block splitting, and block fusion are frequently
observed. An open question in bacterial molecular evolution is that of the
formation and breakup of gene blocks, for which several models have been
proposed. These models, however, are not generally applicable to all types of
gene blocks, and consequently cannot be used to broadly compare and study gene
block evolution. To address this problem we introduce an event-based method for
tracking gene block evolution in bacteria. Results: We show here that the
evolution of gene blocks in proteobacteria can be described by a small set of
events. Those include the insertion of genes into, or the splitting of genes
out of a gene block, gene loss, and gene duplication. We show how the
event-based method of gene block evolution allows us to determine the
evolutionary rate, and to trace the ancestral states of their formation. We
conclude that the event-based method can be used to help us understand the
formation of these important bacterial genomic structures. Availability: The
software is available under GPLv3 license on
http://github.com/reamdc1/gene_block_evolution.git Supplementary online
material: http://iddo-friedberg.net/operon-evolution Contact: Iddo Friedberg
i.friedberg@miamioh.edu