  Binary fluorescence time series obtained from single-molecule imaging
experiments can be used to infer protein binding kinetics, in particular,
association and dissociation rate constants from waiting time statistics of
fluorescence intensity changes. In many cases, rate constants inferred from
fluorescence time series exhibit nonintuitive dependence on ligand
concentration. Here we examine several possible mechanistic and technical
origins that may induce ligand dependence of rate constants. Using aggregated
Markov models, we show under the condition of detailed balance that
non-fluorescent bindings and missed events due to transient interactions,
instead of conformation fluctuations, may underly the dependence of waiting
times and thus apparent rate constants on ligand concentrations. In general,
waiting times are rational functions of ligand concentration. The shape of
concentration dependence is qualitatively affected by the number of binding
sites in the single molecule and is quantitatively tuned by model parameters.
We also show that ligand dependence can be caused by non-equilibrium conditions
which result in violations of detailed balance and require an energy source. As
to a different but significant mechanism, we examine the effect of ambient
buffers that can substantially reduce the effective concentration of ligands
that interact with the single molecules. To demonstrate the effects by these
mechanisms, we applied our results to analyze the concentration dependence in a
single-molecule experiment EGFR binding to fluorophore-labeled adaptor protein
Grb2 by Morimatsu et al. (PNAS,104:18013,2007).


  Genomic regions (or loci) displaying outstanding correlation with some
environmental variables are likely to be under selection and this is the
rationale of recent methods of identifying selected loci and retrieving
functional information about them. To be efficient, such methods need to be
able to disentangle the potential effect of environmental variables from the
confounding effect of population history. For the routine analysis of
genome-wide datasets, one also needs fast inference and model selection
algorithms. We propose a method based on an explicit spatial model which is an
instance of spatial generalized linear mixed model (SGLMM). For inference, we
make use of the INLA-SPDE theoretical and computational framework developed by
Rue et al. (2009) and Lindgren et al (2011). The method we propose allows one
to quantify the correlation between genotypes and environmental variables. It
works for the most common types of genetic markers, obtained either at the
individual or at the population level. Analyzing simulated data produced under
a geostatistical model then under an explicit model of selection, we show that
the method is efficient. We also re-analyze a dataset relative to nineteen pine
weevils (Hylobius abietis}) populations across Europe. The method proposed
appears also as a statistically sound alternative to the Mantel tests for
testing the association between genetic and environmental variables.


  The question of the value of the pitch angle of DNA is visited from the
perspective of a geometrical analysis of transcription. It is suggested that
for transcription to be possible, the pitch angle of B-DNA must be smaller than
the angle of zero-twist. At the zero-twist angle the double helix is maximally
rotated and its strain-twist coupling vanishes. A numerical estimate of the
pitch angle for B-DNA based on differential geometry is compared with numbers
obtained from existing empirical data. The crystallographic studies shows that
the pitch angle is approximately 38 deg., less than the corresponding
zero-twist angle of 41.8 deg., which is consistent with the suggested principle
for transcription.


  Methods for predicting the probability and timing of a species' extinction
are typically based on a combination of theoretical models and empirical data,
and focus on single species population dynamics. Of course, species also
interact with each other, forming more or less complex networks of
interactions. Models to assess extinction risk often lack explicit
incorporation of these interspecific interactions. We study a birth and death
process in which the death rate includes an effect from predation. This
predation rate is included via a general nonlinear expression for the
functional response of predation to prey density. We investigate the effects of
the foraging parameters (e.g. attack rate and handling time) on the mean time
to extinction. Mean time to extinction varies by orders of magnitude when we
alter the foraging parameters, even when we exclude the effects of these
parameters on the equilibrium population size. In particular we observe an
exponential dependence of the mean time to extinction on handling time. These
findings clearly show that accounting for the nature of interspecific
interactions is likely to be critically important when estimating extinction
risk.


  We consider the Wright-Fisher model for a population of $N$ individuals, each
identified with a sequence of a finite number of sites, and single-crossover
recombination between them. We trace back the ancestry of single individuals
from the present population. In the $N \to \infty$ limit without rescaling of
parameters or time, this ancestral process is described by a random tree, whose
branching events correspond to the splitting of the sequence due to
recombination. With the help of a decomposition of the trees into subtrees, we
calculate the probabilities of the topologies of the ancestral trees. At the
same time, these probabilities lead to a semi-explicit solution of the
deterministic single-crossover equation. The latter is a discrete-time
dynamical system that emerges from the Wright-Fisher model via a law of large
numbers and has been waiting for a solution for many decades.


  We describe a molecule-oriented modelling approach based on a collection of
Petri net models organized in the form of modules into a prototype database
accessible through a web interface. The JAK/STAT signalling pathway with the
extensive cross-talk of its components is selected as case study. Each Petri
net module represents the reactions of an individual protein with its specific
interaction partners. These Petri net modules are graphically displayed, can be
executed individually, and allow the automatic composition into coherent models
containing an arbitrary number of molecular species chosen ad hoc by the user.
Each module contains metadata for documentation purposes and can be extended to
a wiki-like minireview. The database can manage multiple versions of each
module. It supports the curation, documentation, version control, and update of
individual modules and the subsequent automatic composition of complex models,
without requiring mathematical skills. Modules can be (semi-) automatically
recombined according to user defined scenarios e.g. gene expression patterns in
given cell types, under certain physiological conditions, or states of disease.
Adding a localisation component to the module database would allow to simulate
models with spatial resolution in the form of coloured Petri nets. As synthetic
biology application we propose the fully automated generation of synthetic or
synthetically rewired network models by composition of metadata-guided
automatically modified modules representing altered protein binding sites.
Petri nets composed from modules can be executed as ODE system, stochastic,
hybrid, or merely qualitative models and exported in SMBL format.


  We investigate the competing effects and relative importance of intrinsic
demographic and environmental variability on the evolutionary dynamics of a
stochastic two-species Lotka-Volterra model by means of Monte Carlo simulations
on a two-dimensional lattice. Individuals are assigned inheritable predation
efficiencies; quenched randomness in the spatially varying reaction rates
serves as environmental noise. We find that environmental variability enhances
the population densities of both predators and prey while demographic
variability leads to essentially neutral optimization.


  We introduce and study the behavior of a tethered membrane of non-zero
thickness embedded in three dimensions subject to an effective self-attraction
induced by hydrophobicity arising from the tendency to minimize the area
exposed to a solvent. The phase behavior and the nature of the folded
conformations are found to be quite distinct in the small and large solvent
size regimes. We demonstrate spontaneous symmetry-breaking with the membrane
folding along a preferential axis, when the solvent molecules are small
compared to the membrane thickness. For large solvent molecule size, a local
crinkling mechanism effectively shields the membrane from the solvent, even in
relatively flat conformations. We discuss the binding/unbinding transition of a
membrane to a wall that serves to shield the membrane from the solvent.


  Background: Autocatalytic sets are often considered a necessary (but not
sufficient) condition for the origin and early evolution of life. Although the
idea of autocatalytic sets was already conceived of many years ago, only
recently have they gained more interest, following advances in creating them
experimentally in the laboratory. In our own work, we have studied
autocatalytic sets extensively from a computational and theoretical point of
view.
  Results: We present results from an initial study of the dynamics of
self-sustaining autocatalytic sets (RAFs). In particular, simulations of
molecular flow on autocatalytic sets are performed, to illustrate the kinds of
dynamics that can occur. Next, we present an extension of our (previously
introduced) algorithm for finding autocatalytic sets in general reaction
networks, which can also handle inhibition. We show that in this case detecting
autocatalytic sets is fixed parameter tractable. Finally, we formulate a
generalized version of the algorithm that can also be applied outside the
context of chemistry and origin of life, which we illustrate with a toy example
from economics.
  Conclusions: Having shown theoretically (in previous work) that autocatalytic
sets are highly likely to exist, we conclude here that also in terms of
dynamics such sets are viable and outcompete non-autocatalytic sets.
Furthermore, our dynamical results confirm arguments made earlier about how
autocatalytic subsets can enable their own growth or give rise to other such
subsets coming into existence. Finally, our algorithmic extension and
generalization show that more realistic scenarios (e.g., including inhibition)
can also be dealt with within our framework, and that it can even be applied to
areas outside of chemistry, such as economics.


  After being considered as a nuisance to be filtered out, it became recently
clear that biochemical noise plays a complex role, often fully functional, for
a genetic network. The influence of intrinsic and extrinsic noises on genetic
networks has intensively been investigated in last ten years, though
contributions on the co-presence of both are sparse. Extrinsic noise is usually
modeled as an unbounded white or colored gaussian stochastic process, even
though realistic stochastic perturbations are clearly bounded. In this paper we
consider Gillespie-like stochastic models of nonlinear networks, i.e. the
intrinsic noise, where the model jump rates are affected by colored bounded
extrinsic noises synthesized by a suitable biochemical state-dependent Langevin
system. These systems are described by a master equation, and a simulation
algorithm to analyze them is derived. This new modeling paradigm should enlarge
the class of systems amenable at modeling.
  We investigated the influence of both amplitude and autocorrelation time of a
extrinsic Sine-Wiener noise on: $(i)$ the Michaelis-Menten approximation of
noisy enzymatic reactions, which we show to be applicable also in co-presence
of both intrinsic and extrinsic noise, $(ii)$ a model of enzymatic futile cycle
and $(iii)$ a genetic toggle switch. In $(ii)$ and $(iii)$ we show that the
presence of a bounded extrinsic noise induces qualitative modifications in the
probability densities of the involved chemicals, where new modes emerge, thus
suggesting the possibile functional role of bounded noises.


  Ion channels are proteins with holes down their middle that control the flow
of ions and electric current across otherwise impermeable biological membranes.
The flow of sodium, potassium, calcium (divalent), and chloride ions have been
central issues in biology for more than a century. The flow of current is
responsible for the signals of the nervous system that propagate over long
distances (meters). The concentration of divalent calcium ions is a 'universal'
signal that controls many different systems inside cells. The concentration of
divalent calcium and other messenger ions has a role in life rather like the
role of the voltage in different wires of a computer. Ion channels also help
much larger solutes (e.g., organic acid and bases; perhaps polypeptides) to
cross membranes but much less is known about these systems. Ion channels can
select and control the movement of different types of ions because the holes in
channel proteins are a few times larger than the (crystal radii of the) ions
themselves. Biology uses ion channels as selective valves to control flow and
thus concentration of crucial chemical signals. For example, the concentration
of divalent calcium ions determines whether muscles contract or not. Ion
channels have a role in biology similar to the role of transistors in computers
and technology. Ion Channels Control Concentrations Important To Life The Way
Computers Control Voltages Important To Computers.


  The yield of chemical reactions involving intermediate radical pairs is
influenced by magnetic fields well beyond the levels expected from energy
considerations. This dependence can be traced back to the microscopic dynamics
of electron spins and constitutes the basis of the chemical compass. Here we
propose a new experimental approach based on molecular photoswitches to achieve
additional control on the chemical reaction and to allow short-time resolution
of the spin dynamics. Our proposal enables experiments to test some of the
standard assumptions of the radical pair model and improves the sensitivity of
chemical magnetometers by two orders of magnitude.


  We model sex-structured population dynamics to analyze pairwise competition
between groups differing both genetically and culturally. A sex-ratio allele is
expressed in the heterogametic sex only, so that assumptions of Fisher's
analysis do not apply. Sex-ratio evolution drives cultural evolution of a
group-associated trait governing mortality in the homogametic sex. The two-sex
dynamics under resource limitation induces a strong Allee effect that depends
on both sex ratio and cultural trait values. We describe the resulting
threshold, separating extinction from positive growth, as a function of female
and male densities. When initial conditions avoid extinction due to the Allee
effect, different sex ratios cannot coexist; in our model, greater female
allocation always invades and excludes a lesser allocation. But the culturally
transmitted trait interacts with the sex ratio to determine the ecological
consequences of successful invasion. The invading female allocation may permit
population persistence at self-regulated equilibrium. For this case, the
resident culture may be excluded, or may coexist with the invader culture. That
is, a single sex-ratio allele in females and a cultural dimorphism in male
mortality can persist; a low-mortality resident trait is maintained by
father-to-son cultural transmission. Otherwise, the successfully invading
female allocation excludes the resident allele and culture, and then drives the
population to extinction via a shortage of males. Finally, we show that the
results obtained under homogeneous mixing hold, with caveats, in a spatially
explicit model with local mating and diffusive dispersal in both sexes.


  Continuous-time Markov chains are a standard tool in phylogenetic inference.
If homogeneity is assumed, the chain is formulated by specifying
time-independent rates of substitutions between states in the chain. In
applications, there are usually extra constraints on the rates, depending on
the situation. If a model is formulated in this way, it is possible to
generalise it and allow for an inhomogeneous process, with time-dependent rates
satisfying the same constraints. It is then useful to require that there exists
a homogeneous average of this inhomogeneous process within the same model. This
leads to the definition of "Lie Markov models", which are precisely the class
of models where such an average exists. These models form Lie algebras and
hence concepts from Lie group theory are central to their derivation. In this
paper, we concentrate on applications to phylogenetics and nucleotide
evolution, and derive the complete hierarchy of Lie Markov models that respect
the grouping of nucleotides into purines and pyrimidines -- that is, models
with purine/pyrimidine symmetry. We also discuss how to handle the subtleties
of applying Lie group methods, most naturally defined over the complex field,
to the stochastic case of a Markov process, where parameter values are
restricted to be real and positive. In particular, we explore the geometric
embedding of the cone of stochastic rate matrices within the ambient space of
the associated complex Lie algebra.
  The whole list of Lie Markov models with purine/pyrimidine symmetry is
available at http://www.pagines.ma1.upc.edu/~jfernandez/LMNR.pdf.


  Chemistry is about chemical reactions. Chemistry is about electrons changing
their configurations as atoms and molecules react. Chemistry studies reactions
as if they occurred in ideal infinitely dilute solutions. But most reactions
occur in nonideal solutions. Then everything (charged) interacts with
everything else (charged) through the electric field, which is short and long
range extending to boundaries of the system. Mathematics has recently been
developed to deal with interacting systems of this sort. The variational theory
of complex fluids has spawned the theory of liquid crystals. In my view, ionic
solutions should be viewed as complex fluids. In both biology and
electrochemistry ionic solutions are mixtures highly concentrated (~10M) where
they are most important, near electrodes, nucleic acids, enzymes, and ion
channels. Calcium is always involved in biological solutions because its
concentration in a particular location is the signal that controls many
biological functions. Such interacting systems are not simple fluids, and it is
no wonder that analysis of interactions, such as the Hofmeister series, rooted
in that tradition, has not succeeded as one would hope. We present a
variational treatment of hard spheres in a frictional dielectric. The theory
automatically extends to spatially nonuniform boundary conditions and the
nonequilibrium systems and flows they produce. The theory is unavoidably
self-consistent since differential equations are derived (not assumed) from
models of (Helmholtz free) energy and dissipation of the electrolyte. The
origin of the Hofmeister series is (in my view) an inverse problem that becomes
well posed when enough data from disjoint experimental traditions are
interpreted with a self-consistent theory.


  Genetic feedback loops in cells break detailed balance and involve
bimolecular reactions; hence exact solutions revealing the nature of the
stochastic fluctuations in these loops are lacking. We here consider the master
equation for a gene regulatory feedback loop: a gene produces protein which
then binds to the promoter of the same gene and regulates its expression. The
protein degrades in its free and bound forms. This network breaks detailed
balance and involves a single bimolecular reaction step. We provide an exact
solution of the steady-state master equation for arbitrary values of the
parameters, and present simplified solutions for a number of special cases. The
full parametric dependence of the analytical non-equilibrium steady-state
probability distribution is verified by direct numerical solution of the master
equations. For the case where the degradation rate of bound and free protein is
the same, our solution is at variance with a previous claim of an exact
solution (Hornos et al, Phys. Rev. E {\bf 72}, 051907 (2005) and subsequent
studies). We show explicitly that this is due to an unphysical formulation of
the underlying master equation in those studies.


  Distance-based methods such as UPGMA (Unweighted Pair Group Method with
Arithmetic Mean) continue to play a significant role in phylogenetic research.
We use polyhedral combinatorics to analyze the natural subdivision of the
positive orthant induced by classifying the input vectors according to tree
topologies returned by the algorithm. The partition lattice informs the study
of UPGMA trees. We give a closed form for the extreme rays of UPGMA cones on n
taxa, and compute the normalized volumes of the UPGMA cones for small n.
  Keywords: phylogenetic trees, polyhedral combinatorics, partition lattice


  Cross-linked actomyosin bundles retract when severed in vivo by laser
ablation, or when isolated from the cell and micromanipulated in vitro in the
presence of ATP. We identify the time scale for contraction as a viscoelastic
time tau, where the viscosity is due to (internal) protein friction. We obtain
an estimate of the order of magnitude of the contraction time tau ~ 10-100 s,
consistent with available experimental data for circumferential microfilament
bundles and stress fibers. Our results are supported by an exactly solvable,
hydrodynamic model of a retracting bundle as a cylinder of isotropic, active
matter, from which the order of magnitude of the active stress is estimated.


  Experimental X-ray crystallography, NMR (Nuclear Magnetic Resonance)
spectroscopy, dual polarization interferometry, etc are indeed very powerful
tools to determine the 3-Dimensional structure of a protein (including the
membrane protein); theoretical mathematical and physical computational
approaches can also allow us to obtain a description of the protein 3D
structure at a submicroscopic level for some unstable, noncrystalline and
insoluble proteins. X-ray crystallography finds the X-ray final structure of a
protein, which usually need refinements using theoretical protocols in order to
produce a better structure. This means theoretical methods are also important
in determinations of protein structures. Optimization is always needed in the
computer-aided drug design, structure-based drug design, molecular dynamics,
and quantum and molecular mechanics. This paper introduces some optimization
algorithms used in these research fields and presents a new theoretical
computational method - an improved LBFGS Quasi-Newtonian mathematical
optimization method - to produce 3D structures of Prion AGAAAAGA amyloid
fibrils (which are unstable, noncrystalline and insoluble), from the potential
energy minimization point of view. Because the NMR or X-ray structure of the
hydrophobic region AGAAAAGA of prion proteins has not yet been determined, the
model constructed by this paper can be used as a reference for experimental
studies on this region, and may be useful in furthering the goals of medicinal
chemistry in this field.


  We find that the marine 87Sr/86Sr record shows a significant periodicity of
59.3 \pm 3 Myr. The 87Sr/86Sr record is 171{\deg} \pm 12{\deg}out of phase with
a 62 (\pm 3) Myr periodicity previously reported in the record of marine-animal
diversity. These periodicities are close to 58 (\pm 4) Myr cycles found for the
number of gap-bounded sedimentary carbonate packages of North America We
propose that these periodicities reflect the operation of a periodic "pulse of
the Earth" in large-scale, Earth processes. These may be linked to mantle or
plate-tectonic events, possibly uplift, which affects Earth's climate and
oceans, and so the geochemistry, sedimentation, and biodiversity of the marine
realm. Alternately, they may be linked to oscillation of the solar system
normal to the plane of the galaxy.


  There may be structural principles pertaining to the general behavior of
systems that lead to similarities in a variety of different contexts. Classic
examples include the descriptive power of fractals, the importance of surface
area to volume constraints, the universality of entropy in systems, and
mathematical rules of growth and form. Documenting such overarching principles
may represent a rejoinder to the Neodarwinian synthesis that emphasizes
adaptation and competition. Instead, these principles could indicate the
importance of constraint and structure on form and evolution. Here we document
a potential example of a phenomenon suggesting congruent behavior of very
different systems. We focus on the notion that universally there has been a
tendency for more volatile entities to disappear from systems such that the net
volatility in these systems tends to decline. We specifically focus on
origination and extinction rates in the marine animal fossil record, the
performance of stocks in the stock market, and the characters of stars and
stellar systems. We consider the evidence that each is experiencing declining
volatility, and also consider the broader significance of this.


  Recent work has shown that the hydrophobic protein surfaces in aqueous
solution sit near a drying transition. The tendency for these surfaces to expel
water from their vicinity leads to self assembly of macromolecular complexes.
In this article we show with a realistic model for a biologically pertinent
system how this phenomenon appears at the molecular level. We focus on the
association of the C-terminal domain (CA-C) of the human immunodeficiency virus
(HIV) capsid protein. By combining all-atom simulations with specialized
sampling techniques we measure the water density distribution during the
approach of two CA-C proteins as a function of separation and amino acid
sequence in the interfacial region. The simulations demonstrate that CA-C
protein-protein interactions sit at the edge of a dewetting transition and that
this mesoscopic manifestation of the underlying liquid-vapor phase transition
can be readily manipulated by biology or protein engineering to significantly
affect association behavior. While the wild type protein remains wet until
contact, we identify a set of in silico mutations, in which three hydrophilic
amino acids are replaced with nonpolar residues, that leads to dewetting prior
to association. The existence of dewetting depends on the size and relative
locations of substituted residues separated by nm length scales, indicating
long range cooperativity and a sensitivity to surface topography. These
observations identify important details which are missing from descriptions of
protein association based on buried hydrophobic surface area.


  Character evolution that affects ecological community interactions often
occurs contemporaneously with temporal changes in population size, potentially
altering the very nature of those dynamics. Such eco-evolutionary processes may
be most readily explored in systems with short generations and simple genetics.
Asexual and cyclically parthenogenetic organisms such as microalgae,
cladocerans, and rotifers, which frequently dominate freshwater plankton
communities, meet these requirements. Multiple clonal lines can coexist within
each species over extended periods, until either fixation occurs or a sexual
phase reshuffles the genetic material. When clones differ in traits affecting
interspecific interactions, within-species clonal dynamics can have major
effects on the population dynamics. We first consider a simple predator-prey
system with two prey genotypes, parameterized with data on a well-studied
experimental system, and explore how the extent of differences in defense
against predation within the prey population determine dynamic stability versus
instability of the system. We then explore how increased potential for
evolution affects the community dynamics in a more general community model with
multiple predator and multiple prey genotypes. These examples illustrate how
microevolutionary "details" that enhance or limit the potential for heritable
phenotypic change can have significant effects on contemporaneous
community-level dynamics and the persistence and coexistence of species.


  In this note we investigate the complexity of the Minimum Label Alignment
problem and we show that such a problem is APX-hard.


  Deformation of single stranded DNA in translocation process before reaching
the pore is investigated. By solving the Laplace equation in a suitable
coordinate system and with appropriate boundary conditions, an approximate
solution for the electric field inside and outside of a narrow pore is
obtained. With an analysis based on "electrohydrodynamic equivalence" we
determine the possibility of extension of a charged polymer due to the presence
of an electric field gradient in the vicinity of the pore entrance. With a
multi-scale hybrid simulation (LB-MD), it is shown that an effective
deformation before reaching the pore occurs which facilitates the process of
finding the entrance for the end monomers. We also highlight the role of long
range hydrodynamic interactions via comparison of the LB-MD results with those
obtained using a Langevin thermostat instead of the LB solver.


  A quantum mechanical model on histone modification is proposed. Along with
the methyl / acetate or other groups bound to the modified residues the torsion
angles of the nearby histone chain are supposed to participate in the quantum
transition cooperatively. The transition rate W is calculated based on the
non-radiative quantum transition theory in adiabatic approximation. By using
W's the reaction equations can be written for histone modification and the
histone modification level can be calculable from the equations, which is
decided by not only the atomic group bound to the modified residue, but also
the nearby histone chain. The theory can explain the mechanism for the
correlation between a pair of chromatin markers observed in histone
modification. The temperature dependence and the coherence-length dependence of
histone modification are deduced. Several points for checking the proposed
theory and the quantum nature of histone modification are suggested as follows:
1, The relationship between lnW and 1/T is same as usual protein folding. The
non-Arhenius temperature dependence of the histone modification level is
predicted. 2, The variation of histone modification level through point
mutation of some residues on the chain is predicted since the mutation may
change the coherence-length of the system. 3, Multi-site modification obeys the
quantum superposition law and the comparison between multi-site transition and
single modification transition gives an additional clue to the testing of the
quantum nature of histone modification.


  We study the dynamics and patterning of polar contractile filaments on the
surface of a cylindrical cell using active hydrodynamic equations that
incorporate couplings between curvature and filament orientation. Cables and
rings spontaneously emerge as steady state configurations on the cylinder, and
can be stationary or moving, helical or segments moving along helical
trajectories. Contractility induces coalescence of proximal rings. We observe
phase transitions in the steady state patterns upon changing cell diameter and
make several testable predictions. Our results are relevant to the dynamics and
patterning of a variety of active biopolymers in cylindrical cells.


  Field theory tools are applied to analytically study fluctuation and
correlation effects in spatially extended stochastic predator-prey systems. In
the mean-field rate equation approximation, the classic Lotka-Volterra model is
characterized by neutral cycles in phase space, describing undamped
oscillations for both predator and prey populations. In contrast, Monte Carlo
simulations for stochastic two-species predator-prey reaction systems on
regular lattices display complex spatio-temporal structures associated with
persistent erratic population oscillations. The Doi-Peliti path integral
representation of the master equation for stochastic particle interaction
models is utilized to arrive at a field theory action for spatial
Lotka-Volterra models in the continuum limit. In the species coexistence phase,
a perturbation expansion with respect to the nonlinear predation rate is
employed to demonstrate that spatial degrees of freedom and stochastic noise
induce instabilities toward structure formation, and to compute the fluctuation
corrections for the oscillation frequency and diffusion coefficient. The
drastic downward renormalization of the frequency and the enhanced diffusivity
are in excellent qualitative agreement with Monte Carlo simulation data.


  Developmental dynamics of multicellular organism is a process that takes
place in a multi-stable system in which each attractor state represents a cell
type and attractor transitions correspond to cell differentiation paths. This
new understanding has revived the idea of a quasi-potential landscape, first
proposed by Waddington as a metaphor. To describe development one is interested
in the "relative stabilities" of N attractors (N>2). Existing theories of state
transition between local minima on some potential landscape deal with the exit
in the transition between a pair attractor but do not offer the notion of a
global potential function that relate more than two attractors to each other.
Several ad hoc methods have been used in systems biology to compute a landscape
in non-gradient systems, such as gene regulatory networks. Here we present an
overview of the currently available methods, discuss their limitations and
propose a new decomposition of vector fields that permit the computation of a
quasi-potential function that is equivalent to the Freidlin-Wentzell potential
but is not limited to two attractors. Several examples of decomposition are
given and the significance of such a quasi-potential function is discussed.


  Many aspects of the historical relationships between populations in a species
are reflected in genetic data. Inferring these relationships from genetic data,
however, remains a challenging task. In this paper, we present a statistical
model for inferring the patterns of population splits and mixtures in multiple
populations. In this model, the sampled populations in a species are related to
their common ancestor through a graph of ancestral populations. Using
genome-wide allele frequency data and a Gaussian approximation to genetic
drift, we infer the structure of this graph. We applied this method to a set of
55 human populations and a set of 82 dog breeds and wild canids. In both
species, we show that a simple bifurcating tree does not fully describe the
data; in contrast, we infer many migration events. While some of the migration
events that we find have been detected previously, many have not. For example,
in the human data we infer that Cambodians trace approximately 16% of their
ancestry to a population ancestral to other extant East Asian populations. In
the dog data, we infer that both the boxer and basenji trace a considerable
fraction of their ancestry (9% and 25%, respectively) to wolves subsequent to
domestication, and that East Asian toy breeds (the Shih Tzu and the Pekingese)
result from admixture between modern toy breeds and "ancient" Asian breeds.
Software implementing the model described here, called TreeMix, is available at
http://treemix.googlecode.com


  Feature article on how networks that govern communication, growth, herd
behavior, and other key processes in nature and society are becoming
increasingly amenable to modeling, forecast, and control.


  Background: MicroRNAs, post-transcriptional repressors of gene expression,
play a pivotal role in gene regulatory networks. They are involved in core
cellular processes and their dysregulation is associated to a broad range of
human diseases. This paper focus on a minimal microRNA-mediated regulatory
circuit, in which a protein-coding gene (host gene) is targeted by a microRNA
located inside one of its introns. Results: Autoregulation via intronic
microRNAs is widespread in the human regulatory network, as confirmed by our
bioinformatic analysis, and can perform several regulatory tasks despite its
simple topology. Our analysis, based on analytical calculations and
simulations, indicates that this circuitry alters the dynamics of the host gene
expression, can induce complex responses implementing adaptation and Weber's
law, and efficiently filters fluctuations propagating from the upstream network
to the host gene. A fine-tuning of the circuit parameters can optimize each of
these functions. Interestingly, they are all related to gene expression
homeostasis, in agreement with the increasing evidence suggesting a role of
microRNA regulation in conferring robustness to biological processes. In
addition to model analysis, we present a list of bioinformatically predicted
candidate circuits in human for future experimental tests. Conclusions: The
results presented here suggest a potentially relevant functional role for
negative self-regulation via intronic microRNAs, in particular as a homeostatic
control mechanism of gene expression. Moreover, the map of circuit functions in
terms of experimentally measurable parameters, resulting from our analysis, can
be a useful guideline for possible applications in synthetic biology.


  Recently, several works have analysed the efficiency of photosynthetic
complexes in a transient scenario and how that efficiency is affected by
environmental noise. Here, following a quantum master equation approach, we
study the energy and excitation transport in fully connected networks both in
general and in the particular case of the Fenna-Matthew-Olson complex. The
analysis is carried out for the steady state of the system where the excitation
energy is constantly "flowing" through the system. Steady state transport
scenarios are particularly relevant if the evolution of the quantum system is
not conditioned on the arrival of individual excitations. By adding dephasing
to the system, we analyse the possibility of noise-enhancement of the quantum
transport.


  Inverted repeat (IR) sequences in DNA can form non-canonical cruciform
structures to relieve torsional stress. We use Monte Carlo simulations of a
recently developed coarse-grained model of DNA to demonstrate that the
nucleation of a cruciform can proceed through a cooperative mechanism. Firstly,
a twist-induced denaturation bubble must diffuse so that its midpoint is near
the centre of symmetry of the IR sequence. Secondly, bubble fluctuations must
be large enough to allow one of the arms to form a small number of hairpin
bonds. Once the first arm is partially formed, the second arm can rapidly grow
to a similar size. Because bubbles can twist back on themselves, they need
considerably fewer bases to resolve torsional stress than the final cruciform
state does. The initially stabilised cruciform therefore continues to grow,
which typically proceeds synchronously, reminiscent of the S-type mechanism of
cruciform formation. By using umbrella sampling techniques we calculate, for
different temperatures and superhelical densities, the free energy as a
function of the number of bonds in each cruciform along the correlated but
non-synchronous nucleation pathways we observed in direct simulations.


  Adaptive-network models are typically studied using deterministic
differential equations which approximately describe their dynamics. In
simulations, however, the discrete nature of the network gives rise to
intrinsic noise which can radically alter the system's behaviour. In this
article we develop a method to predict the effects of stochasticity in adaptive
networks by making use of a pair-based proxy model. The technique is developed
in the context of an epidemiological model of a disease spreading over an
adaptive network of infectious contact. Our analysis reveals that in this model
the structure of the network exhibits stochastic oscillations in response to
fluctuations in the disease dynamic.


  Proteins participating in a protein-protein interaction network can be
grouped into homology classes following their common ancestry. Proteins added
to the network correspond to genes added to the classes, so that the dynamics
of the two objects are intrinsically linked. Here, we first introduce a
statistical model describing the joint growth of the network and the
partitioning of nodes into classes, which is studied through a combined
mean-field and simulation approach. We then employ this unified framework to
address the specific issue of the age dependence of protein interactions,
through the definition of three different node wiring/divergence schemes.
Comparison with empirical data indicates that an age-dependent divergence move
is necessary in order to reproduce the basic topological observables together
with the age correlation between interacting nodes visible in empirical data.
We also discuss the possibility of nontrivial joint partition/topology
observables.


  Two-dimensional gel electrophoresis has been instrumental in the development
of proteomics. Although it is no longer the exclusive scheme used for
proteomics, its unique features make it a still highly valuable tool,
especially when multiple quantitative comparisons of samples must be made, and
even for large samples series. However, quantitative proteomics using 2D gels
is critically dependent on the performances of the protein detection methods
used after the electrophoretic separations. This chapter therefore examines
critically the various detection methods (radioactivity, dyes, fluorescence,
and silver) as well as the data analysis issues that must be taken into account
when quantitative comparative analysis of 2D gels is performed.


  Silver staining is used to detect proteins after electrophoretic separation
on polyacrylamide gels. It -combines excellent sensitivity (in the low nanogram
range) with the use of very simple and cheap equipment and chemicals. For its
use in proteomics, two important additional features must be considered,
compatibility with mass spectrometry and quantitative response. Both features
are discussed in this chapter, and optimized silver staining protocols are
proposed.


  In molecular systematics, evolutionary trees are reconstructed from sequences
at the tips under simple models of site substitution. A central question is how
much sequence data is required to reconstruct a tree accurately? The answer
depends on the lengths of the branches (edges) of the tree, with very short and
very long edges requiring long sequences for accurate tree inference,
particularly when these branch lengths are arranged in certain ways. For
four-taxon trees, the sequence length question was settled for the case of a
rapid speciation event in the distant past. Here, we generalize this result and
show that the same sequence length requirement holds even when the speciation
event is recent, provided that at least one of the four taxa is distantly
related to the others. However, this equivalence disappears if a molecular
clock applies, since the length of the long outgroup edge becomes largely
irrelevant in the estimation of the tree topology for a recent (but not a deep)
divergence. We also show how our results can be extended to models in which
substitution rates vary across sites, and to settings where more than four taxa
are involved.


  A diversity of decision-making systems has been observed in animal
collectives. In some species, choices depend on the differences of the numbers
of animals that have chosen each of the available options, while in other
species on the relative differences (a behavior known as Weber's law) or follow
more complex rules. We here show that this diversity of decision systems
corresponds to a single rule of decision-making in collectives. We first
obtained a decision rule based on Bayesian estimation that uses the information
provided by the behaviors of the other individuals to improve the estimation of
the structure of the world. We then tested this rule in decision experiments
using zebrafish (Danio rerio), and in existing rich datasets of argentine ants
(Linepithema humile) and sticklebacks (Gasterosteus aculeatus), showing that a
unified model across species can quantitatively explain the diversity of
decision systems. Further, these results show that the different counting
systems used by animals, including humans, can emerge from the common principle
of using social information to make good decisions.


  We study the problem of constructing phylogenetic trees for a given set of
species. The problem is formulated as that of finding a minimum Steiner tree on
$n$ points over the Boolean hypercube of dimension $d$. It is known that an
optimal tree can be found in linear time if the given dataset has a perfect
phylogeny, i.e. cost of the optimal phylogeny is exactly $d$. Moreover, if the
data has a near-perfect phylogeny, i.e. the cost of the optimal Steiner tree is
$d+q$, it is known that an exact solution can be found in running time which is
polynomial in the number of species and $d$, yet exponential in $q$. In this
work, we give a polynomial-time algorithm (in both $d$ and $q$) that finds a
phylogenetic tree of cost $d+O(q^2)$. This provides the best guarantees known -
namely, a $(1+o(1))$-approximation - for the case $\log(d) \ll q \ll \sqrt{d}$,
broadening the range of settings for which near-optimal solutions can be
efficiently found. We also discuss the motivation and reasoning for studying
such additive approximations.


  The recA/RAD51 gene family encodes a diverse set of recombinase proteins that
effect homologous recombination, DNA-repair, and genome stability. The recA
gene family is expressed in almost all species of Eubacteria, Archaea, and
Eukaryotes, and even in some viruses. To date, efforts to resolve the deep
evolutionary origins of this ancient protein family have been hindered, in
part, by the high sequence divergence between families (i.e. ~30% identity
between paralogous groups). Through (i) large taxon sampling, (ii) the use of a
phylogenetic algorithm designed for measuring highly divergent paralogs, and
(iii) novel Evolutionary Spatial Dynamics simulation and analytical tools, we
obtained a robust, parsimonious and more refined phylogenetic history of the
recA/RAD51 superfamily. Taken together, our model for the evolution of
recA/RAD51 family provides a better understanding of ancient origin of recA
proteins and multiple events leading to the diversification of recA homologs in
eukaryotes, including the discovery of additional RAD51 sub-families.


  Search for extraterrestrial life and intelligence constitutes one of the
major endeavors in science, but has yet been quantitatively modeled only rarely
and in a cursory and superficial fashion. We argue that probabilistic cellular
automata (PCA) represent the best quantitative framework for modeling
astrobiological history of the Milky Way and its Galactic Habitable Zone. The
relevant astrobiological parameters are to be modeled as the elements of the
input probability matrix for the PCA kernel. With the underlying simplicity of
the cellular automata constructs, this approach enables a quick analysis of
large and ambiguous input parameters' space. We perform a simple clustering
analysis of typical astrobiological histories and discuss the relevant boundary
conditions of practical importance for planning and guiding actual empirical
astrobiological and SETI projects. In addition to showing how the present
framework is adaptable to more complex situations and updated observational
databases from current and near-future space missions, we demonstrate how
numerical results could offer a cautious rationale for continuation of
practical SETI searches.


  Phylogenetically decisive sets of taxon sets have the property that if trees
are chosen for each of their elements, as long as these trees are compatible,
the resulting supertree is unique. This means that as long as the trees
describing the phylogenetic relationships of the (input) species sets are
compatible, they can only be combined into a common supertree in precisely one
way. This setting is sometimes also referred to as \enquote{perfect taxon
sampling}. While for rooted trees, the decision if a given set of input taxon
sets is phylogenetically decisive can be made in polynomial time, the decision
problem to determine whether a set of taxon sets is phylogenetically decisive
concerning \emph{unrooted} trees is unfortunately coNP-complete. This shows
that recognizing such sets is often difficult. In this paper, we explain
phylogenetic decisiveness and introduce a class of input taxon sets, namely
so-called \emph{fixing taxon traceable} sets, which are guaranteed to be
phylogenetically decisive and which can be recognized in polynomial time. Using
both combinatorial approaches as well as simulations, we compare properties of
fixing taxon traceability and phylogenetic decisiveness, e.g., by deriving
lower and upper bounds for the number of quadruple sets needed in the input set
for each of these properties. In particular, we correct an erroneous lower
bound concerning phylogenetic decisiveness from the literature.
  We have implemented the algorithm to determine if a given set of taxon sets
is fixing taxon traceable in \textsf{R} and made our software package
\verb+FixingTaxonTraceR+ publicly available.


  The idea of this project is to study the protein structure and sequence
relationship using the hidden markov model and artificial neural network. In
this context we have assumed two hidden markov models. In first model we have
taken protein secondary structures as hidden and protein sequences as observed.
In second model we have taken protein sequences as hidden and protein
structures as observed. The efficiencies for both the hidden markov models have
been calculated. The results show that the efficiencies of first model is
greater that the second one .These efficiencies are cross validated using
artificial neural network. This signifies the importance of protein secondary
structures as the main hidden controlling factors due to which we observe a
particular amino acid sequence. This also signifies that protein secondary
structure is more conserved in comparison to amino acid sequence.


  Lateral gene transfer (LGT) is a common mechanism of non-vertical evolution
where genetic material is transferred between two more or less distantly
related organisms. It is particularly common in bacteria where it contributes
to adaptive evolution with important medical implications. In evolutionary
studies, LGT has been shown to create widespread discordance between gene trees
as genomes become mosaics of gene histories. In particular, the Tree of Life
has been questioned as an appropriate representation of bacterial evolutionary
history. Nevertheless a common hypothesis is that prokaryotic evolution is
primarily tree-like, but that the underlying trend is obscured by LGT.
Extensive empirical work has sought to extract a common tree-like signal from
conflicting gene trees. Here we give a probabilistic perspective on the problem
of recovering the tree-like trend despite LGT. Under a model of randomly
distributed LGT, we show that the species phylogeny can be reconstructed even
in the presence of surprisingly many (almost linear number of) LGT events per
gene tree. Our results, which are optimal up to logarithmic factors, are based
on the analysis of a robust, computationally efficient reconstruction method
and provides insight into the design of such methods. Finally we show that our
results have implications for the discovery of highways of gene sharing.


  Protein synthesis-dependent, late long-term potentiation (LTP) and depression
(LTD) at glutamatergic hippocampal synapses are well characterized examples of
long-term synaptic plasticity. Persistent increased activity of the enzyme
protein kinase M (PKM) is thought essential for maintaining LTP. Additional
spatial and temporal features that govern LTP and LTD induction are embodied in
the synaptic tagging and capture (STC) and cross capture hypotheses. Only
synapses that have been "tagged" by an stimulus sufficient for LTP and learning
can "capture" PKM. A model was developed to simulate the dynamics of key
molecules required for LTP and LTD. The model concisely represents
relationships between tagging, capture, LTD, and LTP maintenance. The model
successfully simulated LTP maintained by persistent synaptic PKM, STC, LTD, and
cross capture, and makes testable predictions concerning the dynamics of PKM.
The maintenance of LTP, and consequently of at least some forms of long-term
memory, is predicted to require continual positive feedback in which PKM
enhances its own synthesis only at potentiated synapses. This feedback
underlies bistability in the activity of PKM. Second, cross capture requires
the induction of LTD to induce dendritic PKM synthesis, although this may
require tagging of a nearby synapse for LTP. The model also simulates the
effects of PKM inhibition, and makes additional predictions for the dynamics of
CaM kinases. Experiments testing the above predictions would significantly
advance the understanding of memory maintenance.


  We study a class of directed random graphs. In these graphs, the interval
[0,x] is the vertex set, and from each y\in [0,x], directed links are drawn to
points in the interval (y,x] which are chosen uniformly with density one. We
analyze the length of the longest directed path starting from the origin. In
the large x limit, we employ traveling wave techniques to extract the
asymptotic behavior of this quantity. We also study the size of a cascade tree
composed of vertices which can be reached via directed paths starting at the
origin.


  The expansion of tools against HIV transmission has brought increased
interest in epidemiological models that can predict the impact of these
interventions. The EMOD-HIV model was recently compared to eleven other
independently developed mathematical models of HIV transmission to determine
the extent to which they agree about the potential impact of expanded use of
antiretroviral therapy in South Africa. Here we describe in detail the modeling
methodology used to produce the results in this comparison, which we term
EMOD-HIV v0.7. We include a discussion of the structure and a full list of
model parameters. We also discuss the architecture of the model, and its
potential utility in comparing structural assumptions within a single modeling
framework.


  Bundles of polymer filaments are responsible for the rich and unique
mechanical behaviors of many biomaterials, including cells and extracellular
matrices. In fibrin biopolymers, whose nonlinear elastic properties are crucial
for normal blood clotting, protofibrils self-assemble and bundle to form
networks of semiflexible fibers. Here we show that the extraordinary
strain-stiffening response of fibrin networks is a direct reflection of the
hierarchical architecture of the fibrin fibers. We measure the rheology of
networks of unbundled protofibrils and find excellent agreement with an affine
model of extensible wormlike polymers. By direct comparison with these data, we
show that physiological fibrin networks composed of thick fibers can be modeled
as networks of tight protofibril bundles. We demonstrate that the tightness of
coupling between protofibrils in the fibers can be tuned by the degree of
enzymatic intermolecular crosslinking by the coagulation Factor XIII.
Furthermore, at high stress, the protofibrils contribute independently to the
network elasticity, which may reflect a decoupling of the tight bundle
structure. The hierarchical architecture of fibrin fibers can thus account for
the nonlinearity and enormous elastic resilience characteristic of blood clots.


  Many commonly used force fields for protein systems such as AMBER, CHARMM,
GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters of
main-chain torsion-energy terms. Here, we propose a new type of
amino-acid-dependent torsion-energy terms in the force fields. As an example,
we applied this approach to AMBER ff03 force field and determined new
amino-acid-dependent parameters for $\psi$ and $\psi'$ angles for each amino
acid by using our optimization method, which is one of the knowledge-based
approach. In order to test the validity of the new force-field parameters, we
then performed folding simulations of $\alpha$-helical and $\beta$-hairpin
peptides, using the optimized force field. The results showed that the new
force-field parameters gave structures more consistent with the experimental
implications than the original AMBER ff03 force field.


  We present an introduction to RNA databases. The history and technology
behind RNA databases is briefly discussed. We examine differing methods of data
collection and curation, and discuss their impact on both the scope and
accuracy of the resulting databases. Finally, we demonstrate these principals
through detailed examination of four leading RNA databases: Noncode, miRBase,
Rfam, and SILVA.


  Emerging high-throughput technologies have led to a deluge of putative
non-coding RNA (ncRNA) sequences identified in a wide variety of organisms.
Systematic characterization of these transcripts will be a tremendous
challenge. Homology detection is critical to making maximal use of functional
information gathered about ncRNAs: identifying homologous sequence allows us to
transfer information gathered in one organism to another quickly and with a
high degree of confidence. ncRNA presents a challenge for homology detection,
as the primary sequence is often poorly conserved and de novo secondary
structure prediction and search remains difficult. This protocol introduces
methods developed by the Rfam database for identifying "families" of homologous
ncRNAs starting from single "seed" sequences using manually curated sequence
alignments to build powerful statistical models of sequence and structure
conservation known as covariance models (CMs), implemented in the Infernal
software package. We provide a step-by-step iterative protocol for identifying
ncRNA homologs, then constructing an alignment and corresponding CM. We also
work through an example for the bacterial small RNA MicA, discovering a
previously unreported family of divergent MicA homologs in genus Xenorhabdus in
the process.


  We study the effects of motor-generated stresses in disordered three
dimensional fiber networks using a combination of a mean-field, effective
medium theory, scaling analysis and a computational model. We find that motor
activity controls the elasticity in an anomalous fashion close to the point of
marginal stability by coupling to critical network fluctuations. We also show
that motor stresses can stabilize initially floppy networks, extending the
range of critical behavior to a broad regime of network connectivities below
the marginal point. Away from this regime, or at high stress, motors give rise
to a linear increase in stiffness with stress. Finally, we demonstrate that our
results are captured by a simple, constitutive scaling relation highlighting
the important role of non-affine strain fluctuations as a susceptibility to
motor stress.


  Genotype-to-phenotype maps and the related fitness landscapes that include
epistatic interactions are difficult to measure because of their high
dimensional structure. Here we construct such a map using the recently
collected corpora of high-throughput sequence data from the 75 base pairs long
mutagenized E. coli lac promoter region, where each sequence is associated with
its phenotype, the induced transcriptional activity measured by a fluorescent
reporter. We find that the additive (non-epistatic) contributions of individual
mutations account for about two-thirds of the explainable phenotype variance,
while pairwise epistasis explains about 7% of the variance for the full
mutagenized sequence and about 15% for the subsequence associated with protein
binding sites. Surprisingly, there is no evidence for third order epistatic
contributions, and our inferred fitness landscape is essentially single peaked,
with a small amount of antagonistic epistasis. There is a significant selective
pressure on the wild type, which we deduce to be multi-objective optimal for
gene expression in environments with different nutrient sources. We identify
transcription factor (CRP) and RNA polymerase binding sites in the promotor
region and their interactions without difficult optimization steps. In
particular, we observe evidence for previously unexplored genetic regulatory
mechanisms, possibly kinetic in nature. We conclude with a cautionary note that
inferred properties of fitness landscapes may be severely influenced by biases
in the sequence data.


  Retinitis pigmentosa (RP) is a group of inherited degenerative eye diseases
characterized by mutations in the genetic structure of the photoreceptors that
leads to the premature death of both rod and cone photoreceptors. Defects in
particular genes encoding proteins that are involved in either the
photoreceptor structure, phototransduction cascades, or visual cycle are
expressed in the rods but ultimately affect both types of cells. RP is
"typically" manifested by a steady death of rods followed by a period of
stability in which cones survive initially and then inevitably die too. In some
RP cases, rods and cones die off simultaneously or even cone death precedes rod
death (reverse RP). The mechanisms and factors involved in the development of
the different types of RP are not well understood nor have researchers been
able to provide more than a limited number of short-term therapies. In this
work we trace the progression of RP to complete blindness through each subtype
via bifurcation theory. We show that the evolution of RP from one stage to
another often requires the failure of multiple components. Our results indicate
that a delicate balance between the availability of nutrients and the rates of
shedding and renewal of photoreceptors is needed at every stage of RP to halt
its progression. This work provides a framework for future physiological
investigations potentially leading to long-term targeted multi-facet
interventions and therapies dependent on the particular stage and subtype of RP
under consideration. The results of this mathematical model may also give
insight into the progression of many other degenerative eye diseases involving
genetic mutations or secondary photoreceptor death and potential ways to
circumvent these diseases.


  Dawkins' replicator-based conception of evolution has led to widespread
mis-application selectionism across the social sciences because it does not
address the paradox that inspired the theory of natural selection in the first
place: how do organisms accumulate change when traits acquired over their
lifetime are obliterated? This is addressed by von Neumann's concept of a
self-replicating automaton (SRA). A SRA consists of a self-assembly code that
is used in two distinct ways: (1) actively deciphered during development to
construct a self-similar replicant, and (2) passively copied to the replicant
to ensure that it can reproduce. Information that is acquired over a lifetime
is not transmitted to offspring, whereas information that is inherited during
copying is transmitted. In cultural evolution there is no mechanism for
discarding acquired change. Acquired change can accumulate orders of magnitude
faster than, and quickly overwhelm, inherited change due to differential
replication of variants in response to selection. This prohibits a selectionist
but not an evolutionary framework for culture. Recent work on the origin of
life suggests that early life evolved through a non-Darwinian process referred
to as communal exchange that does not involve a self-assembly code, and that
natural selection emerged from this more haphazard, ancestral evolutionary
process. It is proposed that communal exchange provides a more appropriate
evolutionary framework for culture than selectionism. This is supported by a
computational model of cultural evolution and a network-based program for
documenting material cultural history, and it is consistent with high levels of
human cooperation.


  In the standard approach to lattice proteins the models based on nearest
neighbor interaction are used. In this kind of models it is difficult to
explain the existence of secondary structures --- special preferred
conformations of protein chains.
  In the present paper a new lattice model of proteins is proposed which is
based on non-local cooperative interactions. In this model the energy of a
conformation of a polymer is equal to the sum of energies of conformations of
fragments of the polymer chain of the length five.
  It is shown that this quinary lattice model is able to describe at
qualitative level secondary structures of proteins: for this model all
conformations with minimal energy are combinations of lattice models of
alpha--helix and beta--strand. Moreover for lattice polymers of the length not
longer that 38 monomers we can describe all conformations with minimal energy.


  Flocking is a paradigmatic example of collective animal behaviour, where
decentralized interaction rules give rise to a globally ordered state. In the
emergence of order out of self-organization we find similarities between
biological systems, as bird flocks, and some physical systems, as ferromagnets.
In both cases, the tendency of individuals to align to their neighbours gives
rise to a polarized state. There is, however, one crucial difference: the
interaction network within an animal group is not necessarily fixed in time, as
each individual moves and may change its neighbours. Therefore, the dynamical
interaction mechanism in biological and physical system can be quite different,
not only due to the gross disparity in the complexity of the individual
entities, but also because of the potential role of inter-individual motion. To
assess the relevance of this mechanism it is necessary to gain quantitative
experimental information about how much individuals move with respect to each
other within the group. Here, by using data from field observations on
starlings, we study the diffusion properties of individual birds within a flock
and investigate the effect of diffusion on the dynamics of the interaction
network. We find that birds diffuse faster than Brownian particles
(superdiffusion) and in a strongly anisotropic way. We also find that
neighbours change in time exclusively as a consequence of diffusion, so that no
specific mechanism to keep one's neighbours seems to be enforced. Finally, we
study the diffusion properties of birds at the border of the flock. We find
that these individuals remain on the border significantly longer than what
would be expected on the basis of a purely diffusional model, suggesting that
there is a sort barrier a bird must cross to make the transition from border to
interior of the flock.


  RNA duplex stability depends strongly on ionic conditions, and inside cells
RNAs are exposed to both monovalent and multivalent ions. Despite recent
advances, we do not have general methods to quantitatively account for the
effects of monovalent and multivalent ions on RNA stability, and the
thermodynamic parameters for secondary structure prediction have only been
derived at 1M [Na$^+$]. Here, by mechanically unfolding and folding a 20 bp RNA
hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at
different monovalent and mixed monovalent/Mg$^{2+}$ salt conditions. We measure
the unfolding and folding rupture forces and apply Kramers theory to extract
accurate information about the hairpin free energy landscape under tension at a
wide range of ionic conditions. We obtain non-specific corrections for the free
energy of formation of the RNA hairpin and measure how the distance of the
transition state to the folded state changes with force and ionic strength. We
experimentally validate the Tightly Bound Ion model and obtain values for the
persistence length of ssRNA. Finally, we test the approximate rule by which the
non-specific binding affinity of divalent cations at a given concentration is
equivalent to that of monovalent cations taken at 100 fold that concentration
for small molecular constructs.


  Pro-social punishment, whereby cooperators punish defectors, is often
suggested as a mechanism that maintains cooperation in large human groups.
Importantly, models that support this idea have to date only allowed defectors
to be the target of punishment. However, recent empirical work has demonstrated
the existence of anti-social punishment in public goods games. That is,
individuals that defect have been found to also punish cooperators. Some recent
theoretical studies have found that such anti-social punishment can prevent the
evolution of pro-social punishment and cooperation. However, the evolution of
anti-social punishment in group-structured populations has not been formally
addressed. Previous work has informally argued that group-structure must favour
pro-social punishment. Here we formally investigate how two demographic
factors, group size and dispersal frequency, affect selection pressures on pro-
and anti-social punishment. Contrary to the suggestions of previous work, we
find that anti-social punishment can prevent the evolution of pro-social
punishment and cooperation under a range of group structures. Given that
anti-social punishment has now been found in all studied extant human cultures,
the claims of previous models showing the co-evolution of pro-social punishment
and cooperation in group-structured populations should be re-evaluated.


  We introduce a rigorous method to microscopically compute the observables
which characterize the thermodynamics and kinetics of rare macromolecular
transitions for which it is possible to identify a priori a slow reaction
coordinate. In order to sample the ensemble of statistically significant
reaction pathways, we define a biased molecular dynamics (MD) in which
barrier-crossing transitions are accelerated without introducing any unphysical
external force. In contrast to other biased MD methods, in the present approach
the systematic errors which are generated in order to accelerate the transition
can be analytically calculated and therefore can be corrected for. This allows
for a computationally efficient reconstruction of the free-energy profile as a
function of the reaction coordinate and for the calculation of the
corresponding diffusion coefficient. The transition path time can then be
readily evaluated within the Dominant Reaction Pathways (DRP) approach. We
illustrate and test this method by characterizing a thermally activated
transition on a two-dimensional energy surface and the folding of a small
protein fragment within a coarse-grained model.


  Widespread sharing of long, identical-by-descent (IBD) genetic segments is a
hallmark of populations that have experienced recent genetic drift. Detection
of these IBD segments has recently become feasible, enabling a wide range of
applications from phasing and imputation to demographic inference. Here, we
study the distribution of IBD sharing in the Wright-Fisher model. Specifically,
using coalescent theory, we calculate the variance of the total sharing between
random pairs of individuals. We then investigate the cohort-averaged sharing:
the average total sharing between one individual and the rest of the cohort. We
find that for large cohorts, the cohort-averaged sharing is distributed
approximately normally. Surprisingly, the variance of this distribution does
not vanish even for large cohorts, implying the existence of "hyper-sharing"
individuals. The presence of such individuals has consequences for the design
of sequencing studies, since, if they are selected for whole-genome sequencing,
a larger fraction of the cohort can be subsequently imputed. We calculate the
expected gain in power of imputation by IBD, and subsequently, in power to
detect an association, when individuals are either randomly selected or
specifically chosen to be the hyper-sharing individuals. Using our framework,
we also compute the variance of an estimator of the population size that is
based on the mean IBD sharing and the variance in the sharing between inbred
siblings. Finally, we study IBD sharing in an admixture pulse model, and show
that in the Ashkenazi Jewish population the admixture fraction is correlated
with the cohort-averaged sharing.


  Within a fully microscopic setting, we derive a variational principle for the
non-equilibrium steady states of chemical reaction networks, valid for
time-scales over which chemical potentials can be taken to be slowly varying:
at stationarity the system minimizes a global function of the reaction fluxes
with the form of a Hopfield Hamiltonian with Hebbian couplings, that is
explicitly seen to correspond to the rate of decay of entropy production over
time. Guided by this analogy, we show that reaction networks can be formally
re-cast as systems of interacting reactions that optimize the use of the
available compounds by competing for substrates, akin to agents competing for a
limited resource in an optimal allocation problem. As an illustration, we
analyze the scenario that emerges in two simple cases: that of toy (random)
reaction networks and that of a metabolic network model of the human red blood
cell.


  The integration of various types of genomic data into predictive models of
biological networks is one of the main challenges currently faced by
computational biology. Constraint-based models in particular play a key role in
the attempt to obtain a quantitative understanding of cellular metabolism at
genome scale. In essence, their goal is to frame the metabolic capabilities of
an organism based on minimal assumptions that describe the steady states of the
underlying reaction network via suitable stoichiometric constraints,
specifically mass balance and energy balance (i.e. thermodynamic feasibility).
The implementation of these requirements to generate viable configurations of
reaction fluxes and/or to test given flux profiles for thermodynamic
feasibility can however prove to be computationally intensive. We propose here
a fast and scalable stoichiometry-based method to explore the Gibbs energy
landscape of a biochemical network at steady state. The method is applied to
the problem of reconstructing the Gibbs energy landscape underlying metabolic
activity in the human red blood cell, and to that of identifying and removing
thermodynamically infeasible reaction cycles in the Escherichia coli metabolic
network (iAF1260). In the former case, we produce consistent predictions for
chemical potentials (or log-concentrations) of intracellular metabolites; in
the latter, we identify a restricted set of loops (23 in total) in the
periplasmic and cytoplasmic core as the origin of thermodynamic infeasibility
in a large sample ($10^6$) of flux configurations generated randomly and
compatibly with the prior information available on reaction reversibility.


  We analyse a large data set of genetic markers obtained from populations of
Cymodocea nodosa, a marine plant occurring from the East Mediterranean to the
Iberian-African coasts in the Atlantic Ocean. We fully develop and test a
recently introduced methodology to infer the directionality of gene flow based
on the concept of geographical segregation. Using the Jensen-Shannon
divergence, we are able to extract a directed network of gene flow describing
the evolutionary patterns of Cymodocea nodosa. In particular we recover the
genetic segregation that the marine plant underwent during its evolution. The
results are confirmed by natural evidence and are consistent with an
independent cross analysis.


  Many open problems in biology, as in the physical sciences, display nonlinear
and 'chaotic' dynamics, which, to the extent possible, cannot be reasonably
understood. Moreover, mathematical models which aim to predict/estimate unknown
aspects of a biological system cannot provide more information about the set of
biologically meaningful (e.g., 'hidden') states of the system than could be
understood by the designer of the model ab initio. Here, the case is made for
the utilization of such models to shift from a 'predictive' to a 'questioning'
nature, and a simple natural-logarithm variation of the logistic polynomial map
is presented that can invoke questions about protein trafficking in eukaryotic
cells.


  Most sensory cells use cross-membrane chemoreceptors to detect chemical
signals in the environment. The biochemical properties and spatial organization
of chemoreceptors play important roles in achieving and maintaining sensitivity
and accuracy of chemical sensing. Here we investigate the effects of receptor
cooperativity and adaptation on the limits of gradient sensing. We study a
single cell with aggregated chemoreceptor arrays on the cell surface and derive
general formula to the limits for gradient sensing from the uncertainty of
instantaneous receptor activity. In comparison to independent receptors, we
find that cooperativity by non-adaptative receptors could significantly lower
the sensing limit in a chemical concentration range determined by the
biochemical properties of ligand-receptor binding and ligand-induced receptor
activity. Cooperativity by adaptative receptors are beneficial to gradient
sensing within a broad range of background concentrations. Our results also
show that isotropic receptor aggregate layout on the cell surface represents an
optimal configuration to gradient sensing.


  DNA is subject to large deformations in a wide range of biological processes.
Two key examples illustrate how such deformations influence the readout of the
genetic information: the sequestering of eukaryotic genes by nucleosomes, and
DNA looping in transcriptional regulation in both prokaryotes and eukaryotes.
These kinds of regulatory problems are now becoming amenable to systematic
quantitative dissection with a powerful dialogue between theory and experiment.
Here we use a single-molecule experiment in conjunction with a statistical
mechanical model to test quantitative predictions for the behavior of DNA
looping at short length scales, and to determine how DNA sequence affects
looping at these lengths. We calculate and measure how such looping depends
upon four key biological parameters: the strength of the transcription factor
binding sites, the concentration of the transcription factor, and the length
and sequence of the DNA loop. Our studies lead to the surprising insight that
sequences that are thought to be especially favorable for nucleosome formation
because of high flexibility lead to no systematically detectable effect of
sequence on looping, and begin to provide a picture of the distinctions between
the short length scale mechanics of nucleosome formation and looping.


  Representations are internal models of the environment that can provide
guidance to a behaving agent, even in the absence of sensory information. It is
not clear how representations are developed and whether or not they are
necessary or even essential for intelligent behavior. We argue here that the
ability to represent relevant features of the environment is the expected
consequence of an adaptive process, give a formal definition of representation
based on information theory, and quantify it with a measure R. To measure how R
changes over time, we evolve two types of networks---an artificial neural
network and a network of hidden Markov gates---to solve a categorization task
using a genetic algorithm. We find that the capacity to represent increases
during evolutionary adaptation, and that agents form representations of their
environment during their lifetime. This ability allows the agents to act on
sensorial inputs in the context of their acquired representations and enables
complex and context-dependent behavior. We examine which concepts (features of
the environment) our networks are representing, how the representations are
logically encoded in the networks, and how they form as an agent behaves to
solve a task. We conclude that R should be able to quantify the representations
within any cognitive system, and should be predictive of an agent's long-term
adaptive success.


  Chlorosomes are likely the largest and most efficient natural
light-harvesting photosynthetic antenna systems. They are composed of large
numbers of bacteriochlorophylls organized into supramolecular aggregates. We
explore the microscopic origin of the fast excitation energy transfer in the
chlorosome using the recently-resolved structure and atomistic-detail
simulations. Despite the dynamical disorder effects on the electronic
transitions of the bacteriochlorophylls, our simulations show that the exciton
delocalizes over the entire aggregate in about 200 fs. The memory effects
associated to the dynamical disorder assist the exciton diffusion through the
aggregates and enhance the diffusion coefficients as a factor of two as
compared to the model without memory. Furthermore, exciton diffusion in the
chlorosome is found to be highly anisotropic with the preferential transfer
towards the baseplate, which is the next functional element in the
photosynthetic system.


  The problem of identifying splice sites consists of two sub-problems: finding
their boundaries, and characterizing their sequence markers. Other splicing
elements---including, enhancers and silencers---that occur in the intronic and
exonic regions play an important role in splicing activity. Existing methods
for detecting splicing elements are limited to finding either splice sites or
enhancers and silencers, even though these elements are well-known to co-occur.
We introduce SeeSite, an efficient and accurate tool for detecting splice sites
and their complementary exon splicing enhancers (ESEs).
  SeeSite has three stages: graph construction, finding dense subgraphs, and
recovering splice sites and ESEs along with their consensus. The third step
involves solving Consensus Sequence with Outliers, an NP-complete string
clustering problem. We prove that our algorithm for this problem outputs
near-optimal solutions in polynomial time. Using SeeSite we demonstrate that
ESEs are preferentially associated with weaker splice sites, and splice sites
of a certain canonical form co-occur with specific ESEs.


  Statistical fluctuations in population sizes of microbes may be quite large
depending on the nature of their underlying stochastic dynamics. For example,
the variance of the population size of a microbe undergoing a pure birth
process with unlimited resources is proportional to the square of its mean. We
refer to such large fluctuations, with the variance growing as square of the
mean, as Giant Number Fluctuations (GNF). Luria and Delbruck showed that
spontaneous mutation processes in microbial populations exhibit GNF. We explore
whether GNF can arise in other microbial ecologies. We study certain simple
ecological models evolving via stochastic processes: (i) bi-directional
mutation, (ii) lysis-lysogeny of bacteria by bacteriophage, and (iii)
horizontal gene transfer (HGT). For the case of bi-directional mutation
process, we show analytically exactly that the GNF relationship holds at large
times. For the ecological model of bacteria undergoing lysis or lysogeny under
viral infection, we show that if the viral population can be experimentally
manipulated to stay quasi-stationary, the process of lysogeny maps essentially
to one-way mutation process and hence the GNF property of the lysogens follows.
Finally, we show that even the process of HGT may map to the mutation process
at large times, and thereby exhibits GNF.


  The principles underlying the architectural landscape of chromatin beyond the
nucleosome level in living cells remains largely unknown despite its potential
to play a role in mammalian gene regulation. We investigated the 3-dimensional
folding of a 1 Mbp region of human chromosome 11 containing the {\beta}-globin
genes by integrating looping interactions of the insulator protein CTCF
determined comprehensively by chromosome conformation capture (3C) into a
polymer model of chromatin. We find that CTCF-mediated cell type specific
interactions in erythroid cells are organized to favor contacts known to occur
in vivo between the {\beta}-globin locus control region (LCR) and genes. In
these cells, the modeled {\beta}-globin domain folds into a globule with the
LCR and the active globin genes on the periphery. By contrast, in non-erythroid
cells, the globule is less compact with few but dominant CTCF interactions
driving the genes away from the LCR. This leads to a decrease in contact
frequencies that can exceed 1000-fold depending on the stiffness of the
chromatin and the exact positioning of the genes. Our findings show that an
ensemble of CTCF contacts functionally affects spatial distances between
control elements and target genes contributing to chromosomal organization
required for transcription.


  Empirical studies suggest that contact patterns follow heterogeneous
inter-event times, meaning that intervals of high activity are followed by
periods of inactivity. Combined with birth and death of individuals, these
temporal constraints affect the spread of infections in a non-trivial way and
are dependent on the particular contact dynamics. We propose a stochastic model
to generate temporal networks where vertices make instantaneous contacts
following heterogeneous inter-event times, and leave and enter the system at
fixed rates. We study how these temporal properties affect the prevalence of an
infection and estimate R0, the number of secondary infections, by modeling
simulated infections (SIR, SI and SIS) co-evolving with the network structure.
We find that heterogeneous contact patterns cause earlier and larger epidemics
on the SIR model in comparison to homogeneous scenarios. In case of SI and SIS,
the epidemics is faster in the early stages (up to 90% of prevalence) followed
by a slowdown in the asymptotic limit in case of heterogeneous patterns. In the
presence of birth and death, heterogeneous patterns always cause higher
prevalence in comparison to homogeneous scenarios with same average inter-event
times. Our results suggest that R0 may be underestimated if temporal
heterogeneities are not taken into account in the modeling of epidemics.


  We examine a two-person game we call Will-Testing in which the strategy space
for both players is a real number. It has no equilibrium. When an infinitely
large set of players plays this in all possible pairings, there is an
equilibrium for the distribution of strategies which requires all players to
use different strategies. We conjecture this solution could underlie some
phenomena (like pecking orders) observed in animals.


  Biological oscillators coordinate individual cellular components so that they
function coherently and collectively. They are typically composed of multiple
feedback loops, and period mismatch is unavoidable in biological
implementations. We investigated the advantageous effect of this period
mismatch in terms of a synchronization response to external stimuli.
Specifically, we considered two fundamental models of genetic circuits: smooth-
and relaxation oscillators. Using phase reduction and Floquet multipliers, we
numerically analyzed their entrainability under different coupling strengths
and period ratios. We found that a period mismatch induces better entrainment
in both types of oscillator; the enhancement occurs in the vicinity of the
bifurcation on their limit cycles. In the smooth oscillator, the optimal period
ratio for the enhancement coincides with the experimentally observed ratio,
which suggests biological exploitation of the period mismatch. Although the
origin of multiple feedback loops is often explained as a passive mechanism to
ensure robustness against perturbation, we study the active benefits of the
period mismatch, which include increasing the efficiency of the genetic
oscillators. Our findings show a qualitatively different perspective for both
the inherent advantages of multiple loops and their essentiality.


  The fields of structural biology and soft matter have independently sought
out fundamental principles to rationalize protein crystallization. Yet the
conceptual differences and the limited overlap between the two disciplines have
thus far prevented a comprehensive understanding of the phenomenon to emerge.
We conduct a computational study of proteins from the rubredoxin family that
bridges the two fields. Using atomistic simulations, we characterize their
crystal contacts, and accordingly parameterize patchy particle models.
Comparing the phase diagrams of these schematic models with experimental
results enables us to critically examine the assumptions behind the two
approaches. The study also reveals features of protein-protein interactions
that can be leveraged to crystallize proteins more generally.


  We present novel analytical results about ecosystem species diversity that
stem from a proposed coarse grained neutral model based on birth-death
processes. The relevance of the problem lies in the urgency for understanding
and synthesizing both theoretical results of ecological neutral theory and
empirical evidence on species diversity preservation. Neutral model of
biodiversity deals with ecosystems in the same trophic level where per-capita
vital rates are assumed to be species-independent. Close-form analytical
solutions for neutral theory are obtained within a coarse-grained model, where
the only input is the species persistence time distribution. Our results
pertain: the probability distribution function of the number of species in the
ecosystem both in transient and stationary states; the n-points connected time
correlation function; and the survival probability, definned as the
distribution of time-spans to local extinction for a species randomly sampled
from the community. Analytical predictions are also tested on empirical data
from a estuarine fish ecosystem. We find that emerging properties of the
ecosystem are very robust and do not depend on specific details of the model,
with implications on biodiversity and conservation biology.


  Whole and targeted sequencing of human genomes is a promising, increasingly
feasible tool for discovering genetic contributions to risk of complex
diseases. A key step is calling an individual's genotype from the multiple
aligned short read sequences of his DNA, each of which is subject to nucleotide
read error. Current methods are designed to call genotypes separately at each
locus from the sequence data of unrelated individuals. Here we propose
likelihood-based methods that improve calling accuracy by exploiting two
features of sequence data. The first is the linkage disequilibrium (LD) between
nearby SNPs. The second is the Mendelian pedigree information available when
related individuals are sequenced. In both cases the likelihood involves the
probabilities of read variant counts given genotypes, summed over the
unobserved genotypes. Parameters governing the prior genotype distribution and
the read error rates can be estimated either from the sequence data itself or
from external reference data. We use simulations and synthetic read data based
on the 1000 Genomes Project to evaluate the performance of the proposed
methods. An R-program to apply the methods to small families is freely
available at http://med.stanford.edu/epidemiology/PHGC/.


  We have developed a 3D off-lattice stochastic polymerization model to study
subcellular oscillation of Min proteins in the bacteria Escherichia coli, and
used it to investigate the experimental phenomenon of Min oscillation
stuttering. Stuttering was affected by the rate of immediate rebinding of MinE
released from depolymerizing filament tips (processivity), protection of
depolymerizing filament tips from MinD binding, and fragmentation of MinD
filaments due to MinE. Each of processivity, protection, and fragmentation
reduces stuttering, speeds oscillations, and reduces MinD filament lengths.
Neither processivity or tip-protection were, on their own, sufficient to
produce fast stutter-free oscillations. While filament fragmentation could, on
its own, lead to fast oscillations with infrequent stuttering; high levels of
fragmentation degraded oscillations. The infrequent stuttering observed in
standard Min oscillations are consistent with short filaments of MinD, while we
expect that mutants that exhibit higher stuttering frequencies will exhibit
longer MinD filaments. Increased stuttering rate may be a useful diagnostic to
find observable MinD polymerization in experimental conditions.


  Selective sweeps are typically associated with a local reduction of genetic
diversity around the adaptive site. However, selective sweeps can also quickly
carry neutral mutations to observable population frequencies if they arise
early in a sweep and hitchhike with the adaptive allele. We show that the
interplay between mutation and exponential amplification through hitchhiking
results in a characteristic frequency spectrum of the resulting novel haplotype
variation that depends only on the ratio of the mutation rate and the selection
coefficient of the sweep. Based on this result, we develop an estimator for the
selection coefficient driving a sweep. Since this estimator utilizes the novel
variation arising from mutations during a sweep, it does not rely on
preexisting variation and can also be applied to loci that lack recombination.
Compared with standard approaches that infer selection coefficients from the
size of dips in genetic diversity around the adaptive site, our estimator
requires much shorter sequences but sampled at high population depth in order
to capture low-frequency variants; given such data, it consistently outperforms
standard approaches. We investigate analytically and numerically how the
accuracy of our estimator is affected by the decay of the sweep pattern over
time as a consequence of random genetic drift and discuss potential effects of
recombination, soft sweeps, and demography. As an example for its use, we apply
our estimator to deep sequencing data from HIV populations.


  The accumulation of deleterious mutations is driven by rare fluctuations
which lead to the loss of all mutation free individuals, a process known as
Muller's ratchet. Even though Muller's ratchet is a paradigmatic process in
population genetics, a quantitative understanding of its rate is still lacking.
The difficulty lies in the nontrivial nature of fluctuations in the fitness
distribution which control the rate of extinction of the fittest genotype. We
address this problem using the simple but classic model of mutation selection
balance with deleterious mutations all having the same effect on fitness. We
show analytically how fluctuations among the fittest individuals propagate to
individuals of lower fitness and have a dramatically amplified effects on the
bulk of the population at a later time. If a reduction in the size of the
fittest class reduces the mean fitness only after a delay, selection opposing
this reduction is also delayed. This delayed restoring force speeds up Muller's
ratchet. We show how the delayed response can be accounted for using a path
integral formulation of the stochastic dynamics and provide an expression for
the rate of the ratchet that is accurate across a broad range of parameters.


This paper has been withdrawn by the author.


  Motivated by single molecule experiments on biopolymers we explore
equilibrium morphologies and force-extension behavior of copolymers with
hydrophobic segments using Langevin dynamics simulations. We find that the
interplay between different length scales, namely, the persistence length
$\ell_{p}$, and the disorder correlation length $p$, in addition to the
fraction of hydrophobic patches $f$ play a major role in altering the
equilibrium morphologies and mechanical response. In particular, we show a
plethora of equilibrium morphologies for this system, \textit{e.g.} core-shell,
looped (with hybridised hydrophilic-hydrophobic sections), and extended coils
as a function of these parameters. A competition of bending energy and
hybridisation energies between two types of beads determines the equilibrium
morphology. Further, mechanical properties of such polymer architectures are
crucially dependent on their native conformations, and in turn on the disorder
realisation along the chain backbone. Thus, for flexible chains, a globule to
extended coil transition is effected via a tensile force for all disorder
realisations. However, the exact nature of the force-extension curves are
different for the different disorder realisations. In contrast, we find that
force-extension behavior of semi-flexible chains with different equilibrium
configurations \textit{e.g.} core-shell, looped, \textit{etc.} reveal a cascade
of force-induced conformational transitions.


  Mitochondrial adenine nucleotide (AdN) content is regulated through the
Ca2+-activated, electroneutral ATP-Mg/Pi carrier (APC). The APC is a protein in
the mitochondrial carrier super family that localizes to the inner
mitochondrial membrane (IMM). It is known to modulate a number of processes
that depend on mitochondrial AdN content, such as gluconeogenesis, protein
synthesis, and citrulline synthesis. Despite this critical role, a kinetic
model of the underlying mechanism has not been developed and validated. Here, a
biophysical model of the APC is developed that is thermodynamically balanced
and accurately reproduces a number of reported data sets from isolated rat
liver and rat kidney mitochondria. The model is based on an ordered bi-bi
mechanism for hetero-exchange of ATP and Pi and also includes homo-exchanges of
ATP and Pi to explain both the initial rate and time course data on ATP and Pi
transport via the APC. The model invokes seven kinetic parameters regarding the
APC mechanism and three parameters related to matrix pH regulation by external
Pi. These parameters are estimated based on nineteen independent data curves;
the estimated parameters are validated using six additional data curves. The
model takes into account the effects of pH, Mg2+ and Ca2+ on ATP and Pi
transport via the APC and supports the conclusion that the pH gradient across
the IMM serves as the primary driving force for AdN uptake or efflux. Moreover,
computer simulations demonstrate that extra-matrix Ca2+ modulates the turnover
rate of the APC and not the binding affinity of ATP, as previously suggested.