  A mathematical theory on flocking serves the foundation for several
ubiquitous multi-agent phenomena in biology, ecology, sensor networks, economy,
as well as social behavior like language emergence and evolution. Directly
inspired by the recent fundamental works of Cucker and Smale on the
construction and analysis of a generic flocking model, we study the emergent
behavior of Cucker-Smale flocking under hierarchical leadership. The rates of
convergence towards asymptotically coherent group patterns in different
scenarios are established.
  The consistent convergence towards coherent patterns may well reveal the
advantages and necessities of having leaders and leadership in a complex
(biological, technological, economic, or social) system with sufficient
intelligence.


  Secretion and role of autotaxin and lysophosphatidic acid in adipose tissue
In obesity, adipocyte hypertrophy is often associated with recrutement of new
fat cells (adipogenesis) under the control of circulating and local regulatory
factors. Among the different lipids released in the extracellular compartment
of adipocytes, our group found the presence of lysophosphatidic acid (LPA). LPA
is a bioactive phospholipid able to regulate several cell responses via the
activation of specific G-protein coupled membrane receptors. Our group found
that LPA increases preadipocyte proliferation and inhibits adipogenesis via the
activation of LPA1 receptor subtype. Extracellular LPA-synthesis is catalyzed
by a lysophospholipase D secreted by adipocytes : autotaxin (ATX). Adipocyte
ATX expression strongly increases with adipogenesis as well as in individuals
exhibiting type 2 diabetes associated with massive obesity. A possible
contribution of ATX and LPA as paracrine regulators of adipogenesis and obesity
associated diabetes is proposed.


  We present a Perfect Sampling algorithm that can be applied to the Master
Equation of Gene Regulatory Networks (GRNs). The method recasts Gillespie's
Stochastic Simulation Algorithm (SSA) in the light of Markov Chain Monte Carlo
methods and combines it with the Dominated Coupling From The Past (DCFTP)
algorithm to provide guaranteed sampling from the stationary distribution. We
show how the DCFTP-SSA can be generically applied to genetic networks with
feedback formed by the interconnection of linear enzymatic reactions and
nonlinear Monod- and Hill-type elements. We establish rigorous bounds on the
error and convergence of the DCFTP-SSA, as compared to the standard SSA,
through a set of increasingly complex examples. Once the building blocks for
GRNs have been introduced, the algorithm is applied to study properly averaged
dynamic properties of two experimentally relevant genetic networks: the toggle
switch, a two-dimensional bistable system, and the repressilator, a
six-dimensional genetic oscillator.


  A recently proposed model of non-autocatalytic reactions in dipeptide
reactions leading to spontaneous symmetry breaking and homochirality is
examined. The model is governed by activation, polymerization, epimerization
and depolymerization of amino acids. Symmetry breaking is primarily a
consequence of the fact that the rates of reactions involving homodimers and
heterodimers are different, i.e., stereoselective, and on the fact that
epimerization can only occur on the N-terminal residue and not on the Cterminal
residue. This corresponds to an auto-inductive cyclic process that works only
in one sense. It is argued that epimerization mimics both autocatalytic
behavior as well as mutual antagonism - both of which were known to be crucial
for producing full homochirality.


  We present a generic computational framework for the simulation of viral
capsid assembly which is quantitative and specific. Starting from PDB files
containing atomic coordinates, the algorithm builds a coarse grained
description of protein oligomers based on graph rigidity. These reduced protein
descriptions are used in an extended Gillespie algorithm to investigate the
stochastic kinetics of the assembly process. The association rates are obtained
from a diffusive Smoluchowski equation for rapid coagulation, modified to
account for water shielding and protein structure. The dissociation rates are
derived by interpreting the splitting of oligomers as a process of graph
partitioning akin to the escape from a multidimensional well. This modular
framework is quantitative yet computationally tractable, with a small number of
physically motivated parameters. The methodology is illustrated using two
different viruses which are shown to follow quantitatively different assembly
pathways. We also show how in this model the quasi-stationary kinetics of
assembly can be described as a Markovian cascading process in which only a few
intermediates and a small proportion of pathways are present. The observed
pathways and intermediates can be related a posteriori to structural and
energetic properties of the capsid oligomers.


  Using computational modelling, we investigate mechanisms of signal
transduction focusing on the spindle assembly checkpoint where a single
unattached kinetochore is able to signal to prevent cell cycle progression.
This inhibitory signal switches off rapidly once spindle microtubules have
attached to all kinetochores. This requirement tightly constrains the possible
mechanisms. Here we investigate two possible mechanisms for spindle checkpoint
operation in metazoan cells, both supported by recent experiments. The first
involves the free diffusion and sequestration of cell-cycle regulators. This
mechanism is severely constrained both by experimental fluorescence recovery
data and also by the large volumes involved in open mitosis in metazoan cells.
Using a simple mathematical analysis and computer simulation, we find that this
mechanism can generate the inhibition found in experiment but likely requires a
two stage signal amplification cascade. The second mechanism involves spatial
gradients of a short-lived inhibitory signal that propagates first by diffusion
but then primarily via active transport along spindle microtubules. We propose
that both mechanisms may be operative in the metazoan spindle assembly
checkpoint, with either able to trigger anaphase onset even without support
from the other pathway.


  BACKGROUND. Signal recognition and information processing is a fundamental
cellular function, which in part involves comprehensive transcriptional
regulatory (TR) mechanisms carried out in response to complex environmental
signals in the context of the cell's own internal state. However, the network
topological basis of developing such integrated responses remains poorly
understood.
  RESULTS. By studying the TR network of the yeast Saccharomyces cerevisiae we
show that an intermediate layer of transcription factors naturally segregates
into distinct subnetworks. In these topological units transcription factors are
densely interlinked in a largely hierarchical manner and respond to external
signals by utilizing a fraction of these subnets.
  CONCLUSIONS. As transcriptional regulation represents the "slow" component of
overall information processing, the identified topology suggests a model in
which successive waves of transcriptional regulation originating from distinct
fractions of the TR network control robust integrated responses to complex
stimuli.


  Synaptic vesicle (SV) exocytosis mediating neurotransmitter release occurs
spontaneously at low intraterminal calcium concentrations and is stimulated by
a rise in intracellular calcium. Exocytosis is compensated for by the
reformation of vesicles at plasma membrane and endosomes. Although the adaptor
complex AP-3 was proposed to be involved in the formation of SVs from
endosomes, whether its function has an indirect effect on exocytosis remains
unknown. Using mocha mice, which are deficient in functional AP-3, we identify
an AP-3-dependent tetanus neurotoxin-resistant asynchronous release that can be
evoked at hippocampal mossy fiber (MF) synapses. Presynaptic targeting of the
tetanus neurotoxin-resistant vesicle soluble N-ethylmaleimide-sensitive factor
attachment protein receptor (SNARE) tetanus neurotoxin-insensitive
vesicle-associated membrane protein (TI-VAMP) is lost in mocha hippocampal MF
terminals, whereas the localization of synaptobrevin 2 is unaffected. In
addition, quantal release in mocha cultures is more frequent and more sensitive
to sucrose. We conclude that lack of AP-3 results in more constitutive
secretion and loss of an asynchronous evoked release component, suggesting an
important function of AP-3 in regulating SV exocytosis at MF terminals.


  In this work it is shown that three pairs of the factors appear to be the
key, i.e. main factors of a natural classification of protein (canonical) amino
acids within the amino acid (genetic) code. First pair: the factors of the
habit of an amino acid molecule (size and polarity). Second pair: the factors
of the association (type of the amino acid/enzyme reactivity and degree of the
hydrophobicity/hydrophilicity of an amino acid molecule). Third pair: the
factors of the dissociation (degree of the acidity-basicity, over acidic group,
COOH and degree of the basicity/acidity over the basic group, NH2). As a result
of the influence and interdependence of all six factors (measured through
correspondent valid parameters) it appears still one natural classification
into polar and non-polar amino acids, where polar amino acids possess negative
and non-polar, the positive values of hydropathy index.


  The mechanical properties of DNA play a critical role in many biological
functions. For example, DNA packing in viruses involves confining the viral
genome in a volume (the viral capsid) with dimensions that are comparable to
the DNA persistence length. Similarly, eukaryotic DNA is packed in DNA-protein
complexes (nucleosomes) in which DNA is tightly bent around protein spools. DNA
is also tightly bent by many proteins that regulate transcription, resulting in
a variation in gene expression that is amenable to quantitative analysis. In
these cases, DNA loops are formed with lengths that are comparable to or
smaller than the DNA persistence length. The aim of this review is to describe
the physical forces associated with tightly bent DNA in all of these settings
and to explore the biological consequences of such bending, as increasingly
accessible by single-molecule techniques.


  Fish farms represent a growing source of disturbance to shallow benthic
ecosystems like seagrass meadows. Despite some existing insights on the
mechanisms underlying decline, efficient tools to quantitatively predict the
response of benthic communities to fish farm effluents have not yet been
developed. We explored relationships of fish farm organic and nutrient input
rates to the sediments with population dynamics of the key seagrass species
(Posidonia oceanica) in deep meadows growing around four Mediterranean Sea
bream and Sea bass fish farms. We performed 2 annual shoot censuses on
permanent plots at increasing distance from cages. Before each census we
measured sedimentation rates adjacent to the plots using benthic sediment
traps. High shoot mortality rates were recorded near the cages, up to 20 times
greater than at control sites. Recruitment rates remained similar to
undisturbed meadows and could not compensate mortality, leading to rapid
seagrass decline within the first 100 meters from cages. Seagrass mortality
increased with total (R2= 0.47, p< 0.0002), organic matter (R2= 0.36, p=
0.001), nitrogen (R2= 0.34, p= 0.002) and phosphorus (R2= 0.58, p< 3 x 10-5)
sedimentation rates. P. oceanica decline accelerated above a phosphorus loading
threshold of 50 mg m-2 day-1. Benthic sedimentation rates seem a powerful
predictor of seagrass mortality from fish farming, integrating local
hydrodynamics, waste effluents variability and several environmental
mechanisms, fuelled by organic inputs and leading to seagrass loss. Coupling
direct measurements of benthic sedimentation rates with dynamics of key species
is proposed as an efficient way to predict and minimize fish farm impacts to
benthic communities.


  The `star paradox' in phylogenetics is the tendency for a particular resolved
tree to be sometimes strongly supported even when the data is generated by an
unresolved (`star') tree. There have been contrary claims as to whether this
phenomenon persists when very long sequences are considered. This note settles
one aspect of this debate by proving mathematically that there is always a
chance that a resolved tree could be strongly supported, even as the length of
the sequences becomes very large.


  DNA loop formation is one of several mechanisms used by organisms to regulate
genes. The free energy of forming a loop is an important factor in determining
whether the associated gene is switched on or off. In this paper we use an
elastic rod model of DNA to determine the free energy of forming short (50--100
basepair), protein mediated DNA loops. Superhelical stress in the DNA of living
cells is a critical factor determining the energetics of loop formation, and we
explicitly account for it in our calculations. The repressor protein itself is
regarded as a rigid coupler; its geometry enters the problem through the
boundary conditions it applies on the DNA. We show that a theory with these
ingredients is sufficient to explain certain features observed in modulation of
in vivo gene activity as a function of the distance between operator sites for
the lac repressor. We also use our theory to make quantitative predictions for
the dependence of looping on superhelical stress, which may be testable both in
vivo and in single-molecule experiments such as the tethered particle assay and
the magnetic bead assay.


  The ontogeny of osmoregulation was examined in the branchial cavity of
embryonic and early post-embryonic stages of the crayfish Astacus leptodactylus
maintained in freshwater, at the sub-cellular level through the detection of
the sodium-potassium adenosine triphosphatase (Na(+),K(+)-ATPase). The
embryonic rate of development was calculated according to the eye index (EI)
which was 430-450 microm at hatching. The distribution of the enzyme was
identified by immunofluorescence microscopy using a monoclonal antibody
IgGalpha5 raised against the avian alpha-subunit of the Na(+),K(+)-ATPase.
Immunoreactivity staining, indicating the presence of Na(+), K(+)-ATPase
appeared in the gills of late embryos (EI>/=400 microm), i.e. a few days before
hatching time, and steadily increased throughout the late embryonic and early
post-embryonic development. The appearance of the enzyme correlates with the
ability to osmoregulate which also occurs late in the embryonic development at
EI 410-420 microm and with tissue differentiation within the gill filaments.
These observations indicate that the physiological shift from osmoconforming
embryos to hyper-regulating late embryos and post-hatching stages in freshwater
must originate partly from the differentiation in the gill epithelia of
ionocytes which are the site of ion pumping, as suggested by the location of
Na(+),K(+)-ATPase. Only the gills were immunostained and a lack of specific
staining was noted in the lamina and the branchiostegites. Therefore,
osmoregulation through Na(+)active uptake is likely achieved in embryos at the
gill level; all the newly formed gills in embryos function in ion regulation;
other parts of the branchial chamber such as the branchiostegites and lamina do
not appear to be involved in osmoregulation.


  We study a Michaelis-Menten reaction for a single two-state enzyme molecule,
whose transition rates between the two conformations are modulated by an
harmonically oscillating external force. In particular, we obtain a range of
optimal driving frequencies for changing the conformation of the enzyme thereby
controlling the enzymatic activity (i.e. product formation). This analysis
demonstrates that it is, in principle, possible to obtain information about
particular rates within the kinetic scheme.


  We consider communication when there is no agreement about symbols and
meanings. We treat it within the framework of reinforcement learning. We apply
different reinforcement learning models in our studies and simplify the problem
as much as possible. We show that the modelling of the other agent is
insufficient in the simplest possible case, unless the intentions can also be
modelled. The model of the agent and its intentions enable quick agreements
about symbol-meaning association. We show that when both agents assume an
`intention model' about the other agent then the symbol-meaning association
process can be spoiled and symbol meaning association may become hard.


  Exploring the structural topology of genome-based large-scale metabolic
network is essential for investigating possible relations between structure and
functionality. Visualization would be helpful for obtaining immediate
information about structural organization. In this work, metabolic networks of
75 organisms were investigated from a topological point of view. A spread
bow-tie model was proposed to give a clear visualization of the bow-tie
structure for metabolic networks. The revealed topological pattern helps to
design more efficient algorithm specifically for metabolic networks. This
coarse-grained graph also visualizes the vulnerable connections in the network,
and thus could have important implication for disease studies and drug target
identifications. In addition, analysis on the reciprocal links and main cores
in the GSC part of bow-tie also reveals that the bow-tie structure of metabolic
networks has its own intrinsic and significant features which are significantly
different from those of random networks.


  We propose a growing network model that consists of two tunable mechanisms:
growth by merging modules which are represented as complete graphs and a
fitness-driven preferential attachment. Our model exhibits the three prominent
statistical properties are widely shared in real biological networks, for
example gene regulatory, protein-protein interaction, and metabolic networks.
They retain three power law relationships, such as the power laws of degree
distribution, clustering spectrum, and degree-degree correlation corresponding
to scale-free connectivity, hierarchical modularity, and disassortativity,
respectively. After making comparisons of these properties between model
networks and biological networks, we confirmed that our model has inference
potential for evolutionary processes of biological networks.


  Erikkson showed that singular value decomposition(SVD) of flattenings
determined a partition of a phylogenetic tree to be a split. In this paper,
based on his work, we develop new statistically consistent algorithms fit for
grid computing to construct a phylogenetic tree by computing SVD of flattenings
with the small fixed number of rows.


  We consider the evolution of populations under the joint action of mutation
and differential reproduction, or selection. The population is modelled as a
finite-type Markov branching process in continuous time, and the associated
genealogical tree is viewed both in the forward and the backward direction of
time. The stationary type distribution of the reversed process, the so-called
ancestral distribution, turns out as a key for the study of mutation-selection
balance. This balance can be expressed in the form of a variational principle
that quantifies the respective roles of reproduction and mutation for any
possible type distribution. It shows that the mean growth rate of the
population results from a competition for a maximal long-term growth rate, as
given by the difference between the current mean reproduction rate, and an
asymptotic decay rate related to the mutation process; this tradeoff is won by
the ancestral distribution.
  Our main application is the quasispecies model of sequence evolution with
mutation coupled to reproduction but independent across sites, and a fitness
function that is invariant under permutation of sites. Here, the variational
principle is worked out in detail and yields a simple, explicit result.


  The postprandial morphological changes of the intestinal epithelium of
Burmese pythons were examined using fasting pythons and at eight time points
after feeding. In fasting pythons, tightly packed enterocytes possess very
short microvilli and are arranged in a pseudostratified fashion. Enterocyte
width increases by 23% within 24 h postfeeding, inducing significant increases
in villus length and intestinal mass. By 6 days postfeeding, enterocyte volume
had peaked, following as much as an 80% increase. Contributing to enterocyte
hypertrophy is the cellular accumulation of lipid droplets at the tips and
edges of the villi of the proximal and middle small intestine, but which were
absent in the distal small intestine. At 3 days postfeeding, conventional and
environmental scanning electron microscopy revealed cracks and lipid extrusion
along the narrow edges of the villi and at the villus tips. Transmission
electron microscopy demonstrated the rapid postprandial lengthening of
enterocyte microvilli, increasing 4.8-fold in length within 24 h, and the
maintaining of that length through digestion. Beginning at 24 h postfeeding,
spherical particles were found embedded apically within enterocytes of the
proximal and middle small intestine. These particles possessed an annular-like
construction and were stained with the calcium-stain Alizarine red S suggesting
that they were bone in origin. Following the completion of digestion, many of
the postprandial responses were reversed, as observed by the atrophy of
enterocytes, the shortening of villi, and the retraction of the microvilli.
Further exploration of the python intestine will reveal the underlying
mechanisms of these trophic responses and the origin and fate of the engulfed
particles.


  The structural organisation of the viral genome within its protein container,
called the viral capsid, is an important aspect of virus architecture. Many
single-stranded (ss) RNA viruses organise a significant part of their genome in
a dodecahedral cage as a RNA duplex structure that mirrors the symmetry of the
capsid. Bruinsma and Rudnick have suggested a model for the structural
organisation of the RNA in these cages. It is the purpose of this paper to
further develop their approach based on results from the areas of graph theory
and DNA network engineering. We start by suggesting a scenario for pariacoto
virus, a representative of this class of viruses, that is energetically more
favorable than those derived previously. We then show that it is a
representative of a whole family of cage structures that abide to the same
construction principle, and then derive the energetically optimal configuration
for a second family of cage structures along similar lines. Finally, we give
reasons for the conjecture that these two families are more likely to occur in
nature than other scenarios.


  Finding coarse-grained, low-dimensional descriptions is an important task in
the analysis of complex, stochastic models of gene regulatory networks. This
task involves (a) identifying observables that best describe the state of these
complex systems and (b) characterizing the dynamics of the observables. In a
previous paper [13], we assumed that good observables were known a priori, and
presented an equation-free approach to approximate coarse-grained quantities
(i.e, effective drift and diffusion coefficients) that characterize the
long-time behavior of the observables. Here we use diffusion maps [9] to
extract appropriate observables ("reduction coordinates") in an automated
fashion; these involve the leading eigenvectors of a weighted Laplacian on a
graph constructed from network simulation data. We present lifting and
restriction procedures for translating between physical variables and these
data-based observables. These procedures allow us to perform equation-free
coarse-grained, computations characterizing the long-term dynamics through the
design and processing of short bursts of stochastic simulation initialized at
appropriate values of the data-based observables.


  We study an individual-based predator-prey model of biological coevolution,
using linear stability analysis and large-scale kinetic Monte Carlo
simulations. The model exhibits approximate 1/f noise in diversity and
population-size fluctuations, and it generates a sequence of quasi-steady
communities in the form of simple food webs. These communities are quite
resilient toward the loss of one or a few species, which is reflected in
different power-law exponents for the durations of communities and the
lifetimes of species. The exponent for the former is near -1, while the latter
is close to -2. Statistical characteristics of the evolving communities,
including degree (predator and prey) distributions and proportions of basal,
intermediate, and top species, compare reasonably with data for real food webs.


  We present studies for an individual based model of three interacting
populations whose individuals are mobile in a 2D-lattice. We focus on the
pattern formation in the spatial distributions of the populations. Also
relevant is the relationship between pattern formation and features of the
populations' time series. Our model displays travelling waves solutions,
clustering and uniform distributions, all related to the parameters values. We
also observed that the regeneration rate, the parameter associated to the
primary level of trophic chain, the plants, regulated the presence of
predators, as well as the type of spatial configuration.


  Classical models for competition between two species usually predict
exclusion or divergent evolution of resource exploitation. However, recent
experimental data show that coexistence is possible for very similar species
competing for the same resources \textit{without} niche partition. Motivated by
this experimental challenge to classical competition theory, we propose an
individual-based stochastic competition model, which is essentially a
modification of a deterministic Lotka-Volterra type model. The proposed model
of competition dynamics incorporates the effects of a discrete genotype, which
determines the individual's adaptation to the environment, as well as its
interaction with the other species.


  Using an example of physical interactions between proteins, we study how
perturbations propagate in interconnected networks whose equilibrium state is
governed by the law of mass action. We introduce a comprehensive matrix
formalism which predicts the response of this equilibrium to small changes in
total concentrations of individual molecules, and explain it using a heuristic
analogy to a current flow in a network of resistors. Our main conclusion is
that on average changes in free concentrations exponentially decay with the
distance from the source of perturbation. We then study how this decay is
influenced by such factors as the topology of a network, binding strength, and
correlations between concentrations of neighboring nodes. An exact analytic
expression for the decay constant is obtained for the case of uniform
interactions on the Bethe lattice. Our general findings are illustrated using a
real biological network of protein-protein interactions in baker's yeast with
experimentally determined protein concentrations.


  Significant fraction (about 98.5% in humans, 24% in microbe Rickettsia
prowazekii) of most animal genomes is non-coding DNA. Although recent studies
established functions of its certain portions, it remains genomic dark matter.
The paper unravels its unusual nature with time reversal approach. Any genome
emerged in evolutionary selection of the fittest survivors. Survivability of
modern species is extensively quantified. Accurate analysis establishes that
under specified conditions it is dominated by the same law in species from
human to single-cell yeast. Since all violators of the law perished in the
previous evolution, it presents the exact law of unanticipated universal
(rather than species specific natural) evolutionary selection of survivors. The
law implies their rapid hereditary, thus genetic, adaptation which is navigated
by operating system of non-coding DNA. Such adaptation to drastic environmental
changes was a must for survival, thus evolved, in otherwise lethal major mass
extinctions. Navigator genome allows for rapid, artificial included, biological
changes (e.g., Methuselah lifespan). Universal law establishes biological
relativity to age transformation in any species; quantifies applicability of
animal models to humans; implies certain universality in biological complexity
and reduces it to exact science problem. Evolutionary and experimental data
corroborate all above conclusions. Further theoretical study and test-stone
experiments are suggested.


  Consider a genetic locus carrying a strongly beneficial allele which has
recently fixed in a large population. As strongly beneficial alleles fix
quickly, sequence diversity at partially linked neutral loci is reduced. This
phenomenon is known as a selective sweep. The fixation of the beneficial allele
not only affects sequence diversity at single neutral loci but also the joint
allele distribution of several partially linked neutral loci. This distribution
can be studied using the ancestral recombination graph for samples of partially
linked neutral loci during the selective sweep. To approximate this graph, we
extend recent work by Schweinsberg & Durrett 2005 and Etheridge, Pfaffelhuber &
Wakolbinger 2006 using a marked Yule tree for the genealogy at a single neutral
locus linked to a strongly beneficial one. We focus on joint genealogies at two
partially linked neutral loci in the case of large selection coefficients
\alpha and recombination rates \rho = O(\alpha/\log\alpha) between loci. Our
approach leads to a full description of the genealogy with accuracy of O((\log
\alpha)^{-2}) in probability. As an application, we derive the expectation of
Lewontin's D as a measure for non-random association of alleles.


  An equation describing the evolution of phenotypic distribution is derived
using methods developed in statistical physics. The equation is solved by using
the singular perturbation method, and assuming that the number of bases in the
genetic sequence is large. Applying the equation to the mutation-selection
model by Eigen provides the critical mutation rate for the error catastrophe.
Phenotypic fluctuation of clones (individuals sharing the same gene) is
introduced into this evolution equation. With this formalism, it is found that
the critical mutation rate is sometimes increased by the phenotypic
fluctuations, i.e., noise can enhance robustness of a fitted state to mutation.
Our formalism is systematic and general, while approximations to derive more
tractable evolution equations are also discussed.


  The human genotope is the convex hull of all allele frequency vectors that
can be obtained from the genotypes present in the human population. In this
paper we take a few initial steps towards a description of this object, which
may be fundamental for future population based genetics studies. Here we use
data from the HapMap Project, restricted to two ENCODE regions, to study a
subpolytope of the human genotope. We study three different approaches for
obtaining informative low-dimensional projections of this subpolytope. The
projections are specified by projection onto few tag SNPs, principal component
analysis, and archetypal analysis. We describe the application of our geometric
approach to identifying structure in populations based on single nucleotide
polymorphisms.


  In this work it is shown that 20 canonical amino acids (AAs) within genetic
code appear to be a whole system with strict AAs positions; more exactly, with
AAs ordinal number in three variants; first variant 00-19, second 00-21 and
third 00-20. The ordinal number follows from the positions of belonging codons,
i.e. their digrams (or doublets). The reading itself is a reading in quaternary
numbering system if four bases possess the values within a specific logical
square: A = 0, C = 1, G = 2, U = 3. By this, all splittings, distinctions and
classifications of AAs appear to be in accordance to atom and nucleon number
balance as well as to the other physico-chemical properties, such as
hydrophobicity and polarity.


  Organisms are equipped with regulatory systems that display a variety of
dynamical behaviours ranging from simple stable steady states, to switching and
multistability, to oscillations. Earlier work has shown that oscillations in
protein concentrations or gene expression levels are related to the presence of
at least one negative feedback loop in the regulatory network. Here we study
the dynamics of a very general class of negative feedback loops. Our main
result is that in these systems the sequence of maxima and minima of the
concentrations is uniquely determined by the topology of the loop and the
activating/repressing nature of the interaction between pairs of variables.
This allows us to devise an algorithm to reconstruct the topology of
oscillating negative feedback loops from their time series; this method applies
even when some variables are missing from the data set, or if the time series
shows transients, like damped oscillations. We illustrate the relevance and the
limits of validity of our method with three examples: p53-Mdm2 oscillations,
circadian gene expression in cyanobacteria, and cyclic binding of cofactors at
the estrogen-sensitive pS2 promoter.


  Normal mode analysis offers an efficient way of modeling the conformational
flexibility of protein structures. Simple models defined by contact topology,
known as elastic network models, have been used to model a variety of systems,
but the validation is typically limited to individual modes for a single
protein. We use anisotropic displacement parameters from crystallography to
test the quality of prediction of both the magnitude and directionality of
conformational variance. Normal modes from four simple elastic network model
potentials and from the CHARMM forcefield are calculated for a data set of 83
diverse, ultrahigh resolution crystal structures. While all five potentials
provide good predictions of the magnitude of flexibility, the methods that
consider all atoms have a clear edge at prediction of directionality, and the
CHARMM potential produces the best agreement. The low-frequency modes from
different potentials are similar, but those computed from the CHARMM potential
show the greatest difference from the elastic network models. This was
illustrated by computing the dynamic correlation matrices from different
potentials for a PDZ domain structure. Comparison of normal mode results with
anisotropic temperature factors opens the possibility of using ultrahigh
resolution crystallographic data as a quantitative measure of molecular
flexibility. The comprehensive evaluation demonstrates the costs and benefits
of using normal mode potentials of varying complexity. Comparison of the
dynamic correlation matrices suggests that a combination of topological and
chemical potentials may help identify residues in which chemical forces make
large contributions to intramolecular coupling.


  Voltage-dependent ion channels determine the electric properties of axonal
cell membranes. They not only allow the passage of ions through the cell
membrane but also contribute to an additional charging of the cell membrane
resulting in the so-called capacitance loading. The switching of the channel
gates between an open and a closed configuration is intrinsically related to
the movement of gating charge within the cell membrane. At the beginning of an
action potential the transient gating current is opposite to the direction of
the current of sodium ions through the membrane. Therefore, the excitability is
expected to become reduced due to the influence of a gating current. Our
stochastic Hodgkin-Huxley like modeling takes into account both the channel
noise -- i.e. the fluctuations of the number of open ion channels -- and the
capacitance fluctuations that result from the dynamics of the gating charge. We
investigate the spiking dynamics of membrane patches of variable size and
analyze the statistics of the spontaneous spiking. As a main result, we find
that the gating currents yield a drastic reduction of the spontaneous spiking
rate for sufficiently large ion channel clusters. Consequently, this
demonstrates a prominent mechanism for channel noise reduction.


  A chain of connected genes with activation-repression links is analysed. It
is shown that for various promoter activity functions (parametrised by Hill
coefficient) the equations describing the concentrations of transcription
factors, are differential-difference KdV-type with perturbations. In the case
of large Hill coefficient the proteomic signal along the gene network is given
by a superposition of perturbed dark solitons of defocusing
differential-difference mKdV equation. Biological implications are discussed.


  This paper describes characteristic features of networks reconstructed from
gene expression time series data. Several null models are considered in order
to discriminate between informations embedded in the network that are related
to real data, and features that are due to the method used for network
reconstruction (time correlation).


  The sequence-dependent elasticity of double-helical DNA on a nm length scale
can be captured by the rigid base-pair model, whose strains are the relative
position and orientation of adjacent base-pairs. Corresponding elastic
potentials have been obtained from all-atom MD simulation and from
high-resolution structural data. On the scale of a hundred nm, DNA is
successfully described by a continuous worm-like chain model with homogeneous
elastic properties characterized by a set of four elastic constants, which have
been directly measured in single-molecule experiments. We present here a theory
that links these experiments on different scales, by systematically
coarse-graining the rigid base-pair model for random sequence DNA to an
effective worm-like chain description. The average helical geometry of the
molecule is exactly taken into account in our approach. We find that the
available microscopic parameters sets predict qualitatively similar mesoscopic
parameters. The thermal bending and twisting persistence lengths computed from
MD data are 42 and 48 nm, respectively. The static persistence lengths are
generally much higher, in agreement with cyclization experiments. All
microscopic parameter sets predict negative twist-stretch coupling. The
variability and anisotropy of bending stiffness in short random chains lead to
non-Gaussian bend angle distributions, but become unimportant after two helical
turns.


  Despite tremendous advances in the field of genomics, the amount and function
of the large non-coding part of the genome in higher organisms remains poorly
understood. Here we report an observation, made for 37 fully sequenced
eukaryotic genomes, which indicates that eukaryotes require a certain minimum
amount of non-coding DNA (ncDNA). This minimum increases quadratically with the
amount of DNA located in exons. Based on a simple model of the growth of
regulatory networks, we derive a theoretical prediction of the required
quantity of ncDNA and find it to be in excellent agreement with the data. The
amount of additional ncDNA (in basepairs) which eukaryotes require obeys Ndef =
1/2 (Nc / Np) (Nc - Np), where Nc is the amount of exonic DNA, and Np is a
constant of about 10Mb. This value Ndef corresponds to a few percent of the
genome in Homo sapiens and other mammals, and up to half the genome in simpler
eukaryotes. Thus our findings confirm that eukaryotic life depends on a
substantial fraction of ncDNA and also make a prediction of the size of this
fraction, which matches the data closely.


  The extinction of a single species due to demographic stochasticity is
analyzed. The discrete nature of the individual agents and the Poissonian noise
related to the birth-death processes result in local extinction of a metastable
population, as the system hits the absorbing state. The Fokker-Planck
formulation of that problem fails to capture the statistics of large deviations
from the metastable state, while approximations appropriate close to the
absorbing state become, in general, invalid as the population becomes large. To
connect these two regimes, a master equation based on a real space WKB method
is presented, and is shown to yield an excellent approximation for the decay
rate and the extreme events statistics all the way down to the absorbing state.
The details of the underlying microscopic process, smeared out in a mean field
treatment, are shown to be crucial for an exact determination of the extinction
exponent. This general scheme is shown to reproduce the known results in the
field, to yield new corollaries and to fit quite precisely the numerical
solutions. Moreover it allows for systematic improvement via a series expansion
where the small parameter is the inverse of the number of individuals in the
metastable state.


  Using self-organized polymer models, we predict mechanical unfolding and
refolding pathways of ribo-zymes, and the green fluorescent protein. In
agreement with experiments, there are between six and eight unfolding
transitions in the Tetrahymena ribozyme. Depending on the loading rate, the
number of rips in the force-ramp unfolding of the Azoarcus ribozymes is between
two and four. Force-quench refolding of the P4-P6 subdomain of the Tetrahymena
ribozyme occurs through a compact intermediate. Subsequent formation of
tertiary contacts between helices P5b-P6a and P5a/P5c-P4 leads to the native
state. The force-quench refolding pathways agree with ensemble experiments. In
the dominant unfolding route, the N-terminal a helix of GFP unravels first,
followed by disruption of the N terminus b strand. There is a third
intermediate that involves disruption of three other strands. In accord with
experiments, the force-quench refolding pathway of GFP is hierarchic, with the
rate-limiting step being the closure of the barrel.


  Processive molecular motors take more-or-less uniformly sized steps, along
spatially periodic tracks, mostly forwards but increasingly backwards under
loads. Experimentally, the major steps can be resolved clearly within the noise
but one knows biochemically that one or more mechanochemical substeps remain
hidden in each enzymatic cycle. In order to properly interpret experimental
data for back/forward step ratios, mean conditional step-to-step dwell times,
etc., a first-passage analysis has been developed that takes account of hidden
substeps in $N$-state sequential models. The explicit, general results differ
significantly from previous treatments that identify the observed steps with
complete mechanochemical cycles; e.g., the mean dwell times $\tau_+$ and
$\tau_-$ prior to forward and back steps, respectively, are normally {\it
unequal} although the dwell times $\tau_{++}$ and $\tau_{--}$ between {\it
successive} forward and back steps are equal. Illustrative (N=2)-state examples
display a wide range of behavior. The formulation extends to the case of two or
more detectable transitions in a multistate cycle with hidden substeps.


  The organization and mining of malaria genomic and post-genomic data is
highly motivated by the necessity to predict and characterize new biological
targets and new drugs. Biological targets are sought in a biological space
designed from the genomic data from Plasmodium falciparum, but using also the
millions of genomic data from other species. Drug candidates are sought in a
chemical space containing the millions of small molecules stored in public and
private chemolibraries. Data management should therefore be as reliable and
versatile as possible. In this context, we examined five aspects of the
organization and mining of malaria genomic and post-genomic data: 1) the
comparison of protein sequences including compositionally atypical malaria
sequences, 2) the high throughput reconstruction of molecular phylogenies, 3)
the representation of biological processes particularly metabolic pathways, 4)
the versatile methods to integrate genomic data, biological representations and
functional profiling obtained from X-omic experiments after drug treatments and
5) the determination and prediction of protein structures and their molecular
docking with drug candidate structures. Progresses toward a grid-enabled
chemogenomic knowledge space are discussed.


  Body mass is a major life history trait and provides a scale for all living
processes of organisms. Unfortunately body mass cannot be easily measured for
many species, because of the logistical difficulties involved in actually
catching and weighing them. This is particularly true for sea turtles which are
large vertebrates that spend most of their life at sea. Here we developed a
general linear model to predict body mass from 17 morphometric measurements
obtained from 49 leatherback turtles nesting in Awala Yalimapo beach, French
Guiana, South America. A stepwise backward analysis removing independent
parameters with P > 0.001 indicated that body mass can be estimated with 93%
accuracy by Body Mass (kg) = - 709.146 + 3.391MedianBodyCirc (cm) + 2.664SCCL
(cm), where SCCL is the standard curvilinear carapace length and MedianBodyCirc
is body circumference at half of SCCL. A sensibility test showed that this
parsimonious model is robust, as estimated body mass may change by 0.7-1.3% for
1-2 cm changes in SCCL, by 0.9-1.7% for 1-2 cm changes in circumference.
Leatherback turtles from French Guiana were larger and heavier than in all
other nesting sites studied so far, suggesting either that gravid leatherbacks
feed during the nesting season in French Guiana, and/or that this species may
exhibits site-specific growth strategies. Further studies are required to test
these hypotheses, including implementation of similar simple models for other
nesting populations, in order to better understand the life history of this
endangered species.


  We present the first combined analysis of diving behaviour and dispersal
patterns in gravid leatherback turtles during 3 consecutive nesting seasons in
French Guiana. In total 23 turtles were fitted with Argos satellite
transmitters and 16 individuals (including 6 concurrently satellite-tracked)
were equipped with an electronic time-depth recorder for single inter-nesting
intervals, i.e. between two consecutive ovi-positions. The leatherbacks
dispersed over the continental shelf, ranging from the coastal zone to the
shelf break and moved over 546.2 $\pm$ 154.1 km (mean $\pm$ SD) in waters of
French Guiana and neighbouring Surinam. They mostly performed shallow (9.4
$\pm$ 9.2 m) and short (4.4 $\pm$ 3.4 min) dives with a slight diurnal pattern.
They dived deeper as they moved away from the coast suggesting that they were
predominantly following the seabed. Inter-nesting intervals could be divided
into two phases: during the first 75% of the time turtles spent at sea, they
dived on average 47 min h-1 before showing a lower and more variable diving
effort as they came back to the shore. The extended movements of leatherbacks
and the fine analysis of dive shapes suggest that in French Guiana leatherbacks
may feed during the inter-nesting interval, probably to compensate for the
energy costs associated with reproduction. This results in this endangered
species being exposed to high risks of interactions with local fisheries
throughout the continental shelf.


  Considering vast amounts of genomic sequences of mostly unknown
functionality, in-silico prediction of functional regions is an important
enterprise. Many genomic browsers employ GC content, which was observed to be
elevated in gene-rich functional regions. This report shows that the entropy of
di- and tri-nucleotides distributions provides a superior measure of genomic
sequence functionality, and proposes an explanation on why the GC content must
be elevated (closer to 50%) in functional regions. Regions with high entropy
strongly co-localize with exons and provide genome-wide evidences of purifying
selection acting on non-coding regions, such as decreased SNPs density. The
observations suggest that functional non-coding regions are optimised for
mutation load in a way, that transition mutations have less impact on
functionality than transversions, leading to the decrease in transversions to
transitions ratio in functional regions.


  Splicing sites provide unique statistics in human genome due to their large
number and reasonably complete annotation. Analyses of the cumulative SNPs
distribution in splicing sites reveal a few interesting observations. While a
degree of the nucleotide conservation reflects on the SNPs density
monotonically, no detectable changes in the SNPs frequencies spectrum were
found. Semi-conserved nucleotide sites harbor transition mutations
predominantly. We propose that such transition preference is caused by
co-evolution of a site with corresponding binding agents. Since transitions in
humans and similarly in other organisms are almost twice as frequent as
transversions, this adaptation significantly lowers the mutation load.


  For the vast majority of naturally occurring, small, single domain proteins
folding is often described as a two-state process that lacks detectable
intermediates. This observation has often been rationalized on the basis of a
nucleation mechanism for protein folding whose basic premise is the idea that
after completion of a specific set of contacts forming the so-called folding
nucleus the native state is achieved promptly. Here we propose a methodology to
identify folding nuclei in small lattice polymers and apply it to the study of
protein molecules with chain length N=48. To investigate the extent to which
protein topology is a robust determinant of the nucleation mechanism we compare
the nucleation scenario of a native-centric model with that of a sequence
specific model sharing the same native fold. To evaluate the impact of the
sequence's finner details in the nucleation mechanism we consider the folding
of two non- homologous sequences. We conclude that in a sequence-specific model
the folding nucleus is, to some extent, formed by the most stable contacts in
the protein and that the less stable linkages in the folding nucleus are solely
determined by the fold's topology. We have also found that independently of
protein sequence the folding nucleus performs the same `topological' function.
This unifying feature of the nucleation mechanism results from the residues
forming the folding nucleus being distributed along the protein chain in a
similar and well-defined manner that is determined by the fold's topological
features.


  The folding of naturally occurring, single domain proteins is usually
well-described as a simple, single exponential process lacking significant
trapped states. Here we further explore the hypothesis that the smooth energy
landscape this implies, and the rapid kinetics it engenders, arises due to the
extraordinary thermodynamic cooperativity of protein folding. Studying
Miyazawa-Jernigan lattice polymers we find that, even under conditions where
the folding energy landscape is relatively optimized (designed sequences
folding at their temperature of maximum folding rate), the folding of
protein-like heteropolymers is accelerated when their thermodynamic
cooperativity enhanced by enhancing the non-additivity of their energy
potentials. At lower temperatures, where kinetic traps presumably play a more
significant role in defining folding rates, we observe still greater
cooperativity-induced acceleration. Consistent with these observations, we find
that the folding kinetics of our computational models more closely approximate
single-exponential behavior as their cooperativity approaches optimal levels.
These observations suggest that the rapid folding of naturally occurring
proteins is, at least in part, consequences of their remarkably cooperative
folding.


  This paper develops models describing asexual and sexual replication in
sporulating organisms. Replication via sporulation is the replication strategy
for all multicellular life, and may even be observed in unicellular life (such
as with budding yeast). We consider diploid populations replicating via one of
two possible sporulation mechanisms: (1) Asexual sporulation, whereby adult
organisms produce single-celled diploid spores that grow into adults
themselves. (2) Sexual sporulation, whereby adult organisms produce
single-celled diploid spores that divide into haploid gametes. The haploid
gametes enter a haploid "pool", where they may recombine with other haploids to
form a diploid spore that then grows into an adult. We consider a haploid
fusion rate given by second-order reaction kinetics. We work with a simplified
model where the diploid genome consists of only two chromosomes, each of which
may be rendered defective with a single point mutation of the wild-type. We
find that the asexual strategy is favored when the rate of spore production is
high compared to the characteristic growth rate from a spore to a reproducing
adult. Conversely, the sexual strategy is favored when the rate of spore
production is low compared to the characteristic growth rate from a spore to a
reproducing adult. As the characteristic growth time increases, or as the
population density increases, the critical ratio of spore production rate to
organism growth rate at which the asexual strategy overtakes the sexual one is
pushed to higher values. Therefore, the results of this model suggest that, for
complex multicellular organisms, sexual replication is favored at high
population densities, and low growth and sporulation rates.


  I introduce an axiomatic representation, called ecoset, to consider
interactions between species in ecological systems. For interspecific
competition, the exclusion conjecture (Gause) is put in a symbolic way and used
as a basic operational tool to consider more complex cases like two species
with two superposed resources (niche differentiation) and apparent competitors.
Competition between tortoises and invaders in Galapagos islands is considered
as a specific example. Our theory gives us an operative language to consider
elementary process in ecology and open a route to more complex systems.


  A big class of viruses self-assemble from a large number of identical capsid
proteins with long flexible N-terminal tails and ss RNA. We study the role of
the strong Coulomb interaction of positive N-terminal tails with ss RNA in the
kinetics of the in vitro virus self-assembly. Capsid proteins stick to
unassembled chain of ss RNA (which we call "antenna") and slide on it towards
the assembly site. We show that at excess of capsid proteins such
one-dimensional diffusion accelerates self-assembly more than ten times. On the
other hand at excess of ss RNA, antenna slows self-assembly down. Several
experiments are proposed to verify the role of ss RNA antenna.


  The Magicicada spp. life cycles with its prime periods and highly
synchronized emergence have defied reasonable scientific explanation since its
discovery. During the last decade several models and explanations for this
phenomenon appeared in the literature along with a great deal of discussion.
Despite this considerable effort, there is no final conclusion about this long
standing biological problem. Here, we construct a minimal automaton model
without predation/parasitism which reproduces some of these aspects. Our
results point towards competition between different strains with limited
dispersal threshold as the main factor leading to the emergence of prime
numbered life cycles.


  We build a simple model for feedback systems involving small RNA (sRNA)
molecules based on the iron metabolism system in the bacterium E. coli, and
compare it with the corresponding system in H. pylori which uses purely
transcriptional regulation. This reveals several unique features of sRNA based
regulation that could be exploited by cells. Firstly, we show that sRNA
regulation can maintain a smaller turnover of target mRNAs than transcriptional
regulation, without sacrificing the speed of response to external shocks.
Secondly, we propose that a single sRNA can prioritize the usage of different
target mRNAs. This suggests that sRNA regulation would be more common in more
complex systems which need to co-regulate many mRNAs efficiently.


  Genomic duplication-divergence events, which are the primary source of new
protein functions, occur stochastically at a wide range of genomic scales, from
single gene to whole genome duplications. Clearly, this fundamental
evolutionary process must have largely conditioned the emerging structure of
protein-protein interaction (PPI) networks, that control many cellular
activities. We propose and asymptotically solve a general
duplication-divergence model of PPI network evolution based on the statistical
selection of duplication-derived interactions. We also introduce a conservation
index, that formally defines the statistical evolutionary conservation of PPI
networks. Distinct conditions on microscopic parameters are then shown to
control global conservation and topology of emerging PPI networks. In
particular, conserved, non-dense networks, which are the only ones of potential
biological relevance, are also shown to be necessary scale-free.


  RNA co-transcriptional folding has long been suspected to play an active role
in helping proper native folding of ribozymes and structured regulatory motifs
in mRNA untranslated regions. Yet, the underlying mechanisms and coding
requirements for efficient co-transcriptional folding remain unclear.
Traditional approaches have intrinsic limitations to dissect RNA folding paths,
as they rely on sequence mutations or circular permutations that typically
perturb both RNA folding paths and equilibrium structures. Here, we show that
exploiting sequence symmetries instead of mutations can circumvent this problem
by essentially decoupling folding paths from equilibrium structures of designed
RNA sequences. Using bistable RNA switches with symmetrical helices conserved
under sequence reversal, we demonstrate experimentally that native and
transiently formed helices can guide efficient co-transcriptional folding into
either long-lived structure of these RNA switches. Their folding path is
controlled by the order of helix nucleations and subsequent exchanges during
transcription, and may also be redirected by transient antisense interactions.
Hence, transient intra- and intermolecular base pair interactions can
effectively regulate the folding of nascent RNA molecules into different native
structures, provided limited coding requirements, as discussed from an
information theory perspective. This constitutive coupling between RNA
synthesis and RNA folding regulation may have enabled the early emergence of
autonomous RNA-based regulation networks.


  An increasing number of proteins are being discovered with a remarkable and
somewhat surprising feature, a knot in their native structures. How the
polypeptide chain is able to knot itself during the folding process to form
these highly intricate protein topologies is not known. Here, we perform a
computational study on the 160-amino acid homodimeric protein YibK which, like
other proteins in the SpoU family of MTases, contains a deep trefoil knot in
its C-terminal region. In this study, we use a coarse-grained C-alpha-chain
representation and Langevin dynamics to study folding kinetics. We find that
specific, attractive nonnative interactions are critical for knot formation. In
the absence of these interactions, i.e. in an energetics driven entirely by
native interactions, knot formation is exceedingly unlikely. Further, we find,
in concert with recent experimental data on YibK, two parallel folding pathways
which we attribute to an early and a late formation of the trefoil knot,
respectively. For both pathways, knot formation occurs before dimerization. A
bioinformatics analysis of the SpoU family of proteins reveals further that the
critical nonnative interactions may originate from evolutionary conserved
hydrophobic segments around the knotted region.


  We study how correlations in the random fitness assignment may affect the
structure of fitness landscapes. We consider three classes of fitness models.
The first is a continuous phenotype space in which individuals are
characterized by a large number of continuously varying traits such as size,
weight, color, or concentrations of gene products which directly affect
fitness. The second is a simple model that explicitly describes
genotype-to-phenotype and phenotype-to-fitness maps allowing for neutrality at
both phenotype and fitness levels and resulting in a fitness landscape with
tunable correlation length. The third is a class of models in which particular
combinations of alleles or values of phenotypic characters are "incompatible"
in the sense that the resulting genotypes or phenotypes have reduced (or zero)
fitness. This class of models can be viewed as a generalization of the
canonical Bateson-Dobzhansky-Muller model of speciation. We also demonstrate
that the discrete NK model shares some signature properties of models with high
correlations. Throughout the paper, our focus is on the percolation threshold,
on the number, size and structure of connected clusters, and on the number of
viable genotypes.


  Owing to the complexity of higher eukaryotic cells, a complete proteome is
likely to be very difficult to achieve. However, advantage can be taken of the
cell compartmentalization to build organelle proteomes, which can moreover be
viewed as specialized tools to study specifically the biology and "physiology"
of the target organelle. Within this frame, we report here the construction of
the human mitochondrial proteome, using placenta as the source tissue. Protein
identification was carried out mainly by peptide mass fingerprinting. The
optimization steps in two-dimensional electrophoresis needed for proteome
research are discussed. However, the relative paucity of data concerning
mitochondrial proteins is still the major limiting factor in building the
corresponding proteome, which should be a useful tool for researchers working
on human mitochondria and their deficiencies.


  The solubilizing power of various nonionic and zwitterionic detergents as
membrane protein solubilizers for two-dimensional electrophoresis was
investigated. Human red blood cell ghosts and Arabidopsis thaliana leaf
membrane proteins were used as model systems. Efficient detergents could be
found in each class, i.e. with oligooxyethylene, sugar or sulfobetaine polar
heads. Among the commercially available nonionic detergents, dodecyl maltoside
and decaethylene glycol mono hexadecyl ether proved most efficient. They
complement the more classical sulfobetaine detergents to widen the scope of
useful detergents for the solubilization of membrane proteins in proteomics.


  Because of the outstanding separating capabilities of two-dimensional
electrophoresis for complete proteins, it would be advantageous to be able to
apply it to all types of proteins. Unfortunately, severe solubility problems
hamper the analysis of many classes of proteins, but especially membrane
proteins. These problems arise mainly in the extraction and isoelectric
focusing steps, and solutions are sought to improve protein solubility under
the conditions prevailing during isoelectric focusing. These solutions deal
mainly with chaotropes and new detergents, which are both able to enhance
protein solubility. The input of these compounds in proteomics analysis of
membrane proteins is discussed, as well as future directions.


  The structure of the self-cleaving hairpin ribozyme is well characterized,
and its folding has been examined in bulk and by single-molecule fluorescence,
establishing the importance of cations, especially magnesium in the stability
of the native fold. Here we describe the first all-atom folding simulations of
the hairpin ribozyme, using a version of a Go potential with separate secondary
and tertiary structure energetic contributions. The ratio of tertiary/secondary
interaction energies serves as a proxy for non-specific cation binding: a high
ratio corresponds to a high concentration, while a low one mimics low
concentration. By studying the unfolding behavior of the RNA over a range of
temperature and tertiary/secondary energies, a three-state phase diagram
emerges, with folded, unfolded (coil) and transient folding/unfolding tertiary
structure species. The thermodynamics were verified by paired folding
simulations in each region of the phase diagram. The three phase behaviors
correspond with experimentally observed states, so this simple model captures
the essential aspect of thermodynamics in RNA folding.


  The influence of thiol blocking on the resolution of basic proteins by
two-dimensional electrophoresis was investigated. Cysteine blocking greatly
increased resolution and decreased streaking, especially in the basic region of
the gels. Two strategies for cysteine blocking were found to be efficient:
classical alkylation with maleimide derivatives and mixed disulfide exchange
with an excess of a low molecular weight disulfide. The effect on resolution
was significant enough to allow correct resolution of basic proteins with
in-gel rehydration on wide gradients (e.g. 3-10 and 4-12), but anodic
cup-loading was still required for basic gradients (e.g. 6-12 or 8-12). These
results demonstrate that thiol-related problems are not solely responsible for
streaking of basic proteins on two-dimensional gels.


  The refolding from stretched initial conformations of ubiquitin (PDB ID:
1ubq) under the quenched force is studied using the Go model and the Langevin
dynamics. It is shown that the refolding decouples the collapse and folding
kinetics. The force quench refolding times scale as tau_F ~ exp(f_q*x_F/k_B*T),
where f_q is the quench force and x_F = 0.96 nm is the location of the average
transition state along the reaction coordinate given by the end-to-end
distance. This value is close to x_F = 0.8 nm obtained from the force-clamp
experiments. The mechanical and thermal unfolding pathways are studied and
compared with the experimental and all-atom simulation results in detail. The
sequencing of thermal unfolding was found to be markedly different from the
mechanical one. It is found that fixing the N-terminus of ubiquitin changes its
mechanical unfolding pathways much more drastically compared to the case when
the C-end is anchored. We obtained the distance between the native state and
the transition state x_UF=0.24 nm which is in reasonable agreement with the
experimental data.


  The proper functioning of mitochondria requires that both the mitochondrial
and the nuclear genome are functional. To investigate the importance of the
mitochondrial genome, which encodes only 13 subunits of the respiratory
complexes, the mitochondrial rRNAs and a few tRNAs, we performed a comparative
study on the 143B cell line and on its Rho-0 counterpart, i.e., devoid of
mitochondrial DNA. Quantitative differences were found, of course in the
respiratory complexes subunits, but also in the mitochondrial translation
apparatus, mainly mitochondrial ribosomal proteins, and in the ion and protein
import system, i.e., including membrane proteins. Various mitochondrial
metabolic processes were also altered, especially electron transfer proteins
and some dehydrogenases, but quite often on a few proteins for each pathway.
This study also showed variations in some hypothetical or poorly characterized
proteins, suggesting a mitochondrial localization for these proteins. Examples
include a stomatin-like protein and a protein sharing homologies with bacterial
proteins implicated in tyrosine catabolism. Proteins involved in apoptosis
control are also found modulated in Rho-0 mitochondria.


  Dendritic cells (DCs) display the unique ability to activate naive T cells
and to initiate primary T cell responses revealed in DC-T cell alloreactions.
DCs frequently operate under stress conditions. Oxidative stress enhances the
production of inflammatory cytokines by DCs. We performed a proteomic analysis
to see which major changes occur, at the protein expression level, during DC
differentiation and maturation. Comparative two-dimensional gel analysis of the
monocyte, immature DC, and mature DC stages was performed. Manganese superoxide
dismutase (Mn-SOD) reached 0.7% of the gel-displayed proteins at the mature DC
stage. This important amount of Mn-SOD is a primary antioxidant defense system
against superoxide radicals, but its product, H(2)O(2), is also deleterious for
cells. Peroxiredoxin (Prx) enzymes play an important role in eliminating such
peroxide. Prx1 expression level continuously increased during DC
differentiation and maturation, whereas Prx6 continuously decreased, and Prx2
peaked at the immature DC stage. As a consequence, DCs were more resistant than
monocytes to apoptosis induced by high amounts of oxidized low density
lipoproteins containing toxic organic peroxides and hydrogen peroxide.
Furthermore DC-stimulated T cells produced high levels of receptor activator of
nuclear factor kappaB ligand, a chemotactic and survival factor for monocytes
and DCs. This study provides insights into the original ability of DCs to
express very high levels of antioxidant enzymes such as Mn-SOD and Prx1, to
detoxify oxidized low density lipoproteins, and to induce high levels of
receptor activator of nuclear factor kappaB ligand by the T cells they activate
and further emphasizes the role that DCs might play in atherosclerosis, a
pathology recognized as a chronic inflammatory disorder.


  Sequence coverage in MS analysis of protein digestion-derived peptides is a
key issue for detailed characterization of proteins or identification at low
quantities. In gel-based proteomics studies, the sequence coverage greatly
depends on the protein detection method. It is shown here that ammoniacal
silver detection methods offer improved sequence coverage over standard silver
nitrate methods, while keeping the high sensitivity of silver staining. With
the development of 2D-PAGE-based proteomics, another burden is placed on the
detection methods used for protein detection on 2-D-gels. Besides the classical
requirements of linearity, sensitivity, and homogeneity from one protein to
another, detection methods must now take into account another aspect, namely
their compatibility with MS. This compatibility is evidenced by two different
and complementary aspects, which are (i) the absence of adducts and artefactual
modifications on the peptides obtained after protease digestion of a protein
detected and digested in - gel, and (ii) the quantitative yield of peptides
recovered after digestion and analyzed by the mass spectrometer. While this
quantitative yield is not very important per se, it is however a crucial
parameter as it strongly influences the S/N of the mass spectrum and thus the
number of peptides that can be detected from a given protein input, especially
at low protein amounts. This influences in turn the sequence coverage and thus
the detail of the analysis provided by the mass spectrometer.


  Significant phylogenetic codivergence between plant or animal hosts ($H$) and
their symbionts or parasites ($P$) indicate the importance of their
interactions on evolutionary time scales. However, valid and realistic methods
to test for codivergence are not fully developed. One of the systems where
possible codivergence has been of interest involves the large subfamily of
temperate grasses (Pooideae) and their endophytic fungi (epichloae). These
widespread symbioses often help protect host plants from herbivory and
stresses, and affect species diversity and food web structures. Here we
introduce the MRCALink (most-recent-common-ancestor link) method and use it to
investigate the possibility of grass-epichlo\"e codivergence. MRCALink applied
to ultrametric $H$ and $P$ trees identifies all corresponding nodes for
pairwise comparisons of MRCA ages. The result is compared to the space of
random $H$ and $P$ tree pairs estimated by a Monte Carlo method. Compared to
tree reconciliation the method is less dependent on tree topologies (which
often can be misleading), and it crucially improves on phylogeny-independent
methods such as {\tt ParaFit} or the Mantel test by eliminating an extreme (but
previously unrecognized) distortion of node-pair sampling. Analysis of 26 grass
species-epichlo\"e species symbioses did not reject random association of $H$
and $P$ MRCA ages. However, when five obvious host jumps were removed the
analysis significantly rejected random association and supported
grass-endophyte codivergence. Interestingly, early cladogenesis events in the
Pooideae corresponded to early cladogenesis events in epichloae, suggesting
concomitant origins of this grass subfamily and its remarkable group of
symbionts. We also applied our method to the well-known gopher-louse data set.


  The distribution of inclusion-rich domains in membranes with active two-state
inclusions is studied by simulations. Our study shows that typical size of
inclusion-rich domains ($L$) can be controlled by inclusion activities in
several ways. When there is effective attraction between state-1 inclusions, we
find: (i) Small domains with only several inclusions are observed for
inclusions with time scales ($\sim 10^{-3} {\rm s}$) and interaction energy
[$\sim \mathcal{O}({\rm k_BT})$] comparable to motor proteins. (ii) $L$ scales
as 1/3 power of the lifetime of state-1 for a wide range of parameters. (iii)
$L$ shows a switch-like dependence on state-2 lifetime $k_{12}^{-1}$. That is,
$L$ depends weakly on $k_{12}$ when $k_{12} < k_{12}^*$ but increases rapidly
with $k_{12}$ when $k_{12} > k_{12}^*$, the crossover $k_{12}^*$ occurs when
the diffusion length of a typical state-2 inclusion within its lifetime is
comparable to $L$. (iv) Inclusion-curvature coupling provides another length
scale that competes with the effects of transition rates.


  Natural proteins fold to a unique, thermodynamically dominant state. Modeling
of the folding process and prediction of the native fold of proteins are two
major unsolved problems in biophysics. Here, we show successful all-atom ab
initio folding of a representative diverse set of proteins, using a minimalist
transferable energy model that consists of two-body atom-atom interactions,
hydrogen-bonding, and a local sequence energy term that models
sequence-specific chain stiffness. Starting from a random coil, the native-like
structure was observed during replica exchange Monte Carlo (REMC) simulation
for most proteins regardless of their structural classes; the lowest energy
structure was close to native- in the range of 2-6 A root-mean-square deviation
(RMSD). Our results demonstrate that the successful all-atom folding of a
protein chain to its native state is governed by only a few crucial energetic
terms.


  Bacterial communities such as biofilms are widely recognised as being
important for survival and persistence of bacteria in harsh environments.
Mechanistic models of biofilm growth indicate that the way in which the surface
is seeded can effect the morphology of simulated biofilms. Experimental studies
indicate that genes which are important for chemotaxis also influence biofilm
formation, perhaps by influencing aggregation on a surface. Understanding
aggregation and microcolony formation could therefore help clarify factors
influencing biofilm formation as well as understanding how groups may influence
the fitness of bacteria. In this paper I develop an Individual Based Model to
examine how different behaviours involved in microcolony formation on a surface
determines patterns of group sizes, and link patterns to bacterial fitness. I
also provide a method for comparing data with model hypotheses to identify
bacterial behaviours in experimental systems.


  A new, versatile, multiphasic buffer system for high-resolution sodium
dodecyl sulfate-polyacrylamide gel electrophoresis of proteins in the relative
molecular weight range of 300 000-3000 Da is described. The system, based on
the theory of multiphasic zone electrophoresis, allows complete stacking and
destacking of proteins in the above M(r) range. The buffer system uses taurine
and chloride as trailing and leading ion, respectively, and Tris, at a pH close
to its pK(a), as the buffering counterion. Coupled with limited variation in
the acrylamide concentration, this electrophoresis system allows to tailor the
resolution in the 6-200 kDa M(r) range, with minimal difficulties in the post
electrophoretic identification processes.


  While the statistical mechanical description of DNA has a long tradition,
renewed interest in DNA melting from a physics perspective is nourished by
measurements of the fluctuation dynamics of local denaturation bubbles by
single molecule spectroscopy. The dynamical opening of DNA bubbles (DNA
breathing) is supposedly crucial for biological functioning during, for
instance, transcription initiation and DNA's interaction with selectively
single-stranded DNA binding proteins. Motivated by this, we consider the bubble
breathing dynamics in a heteropolymer DNA based on a (2+1)-variable master
equation and complementary stochastic Gillespie simulations, providing the
bubble size and the position of the bubble along the sequence as a function of
time. We utilize new experimental data that independently obtain stacking and
hydrogen bonding contributions to DNA stability. We calculate the spectrum of
relaxation times and the experimentally measurable autocorrelation function of
a fluorophore-quencher tagged base-pair, and demonstrate good agreement with
fluorescence correlation experiments. A significant dependence of opening
probability and waiting time between bubble events on the local DNA sequence is
revealed and quantified for a promoter sequence of the T7 phage. The strong
dependence on sequence, temperature and salt concentration for the breathing
dynamics of DNA found here points at a good potential for nanosensing
applications by utilizing short fluorophore-quencher dressed DNA constructs.


  A class of models of biological population and communities with a singular
equilibrium at the origin is analyzed; it is shown that these models can
possess a dynamical regime of deterministic extinction, which is crucially
important from the biological standpoint. This regime corresponds to the
presence of a family of homoclinics to the origin, so-called elliptic sector.
The complete analysis of possible topological structures in a neighborhood of
the origin, as well as asymptotics to orbits tending to this point, is given.
An algorithmic approach to analyze system behavior with parameter changes is
presented. The developed methods and algorithm are applied to existing
mathematical models of biological systems. In particular, we analyze a model of
anticancer treatment with oncolytic viruses, a parasite-host interaction model,
and a model of Chagas' disease.


  We study in detail a particular statistical method in genetic case-control
analysis, labeled "genotype-based association", in which the two test results
from assuming dominant and recessive model are combined in one optimal output.
This method differs both from the allele-based association which artificially
doubles the sample size, and the direct chi-square test on 3-by-2 contingency
table which may overestimate the degree of freedom. We conclude that the
comparative advantage (or disadvantage) of the genotype-based test over the
allele-based test mainly depends on two parameters, the allele frequency
difference delta and the Hardy-Weinberg disequilibrium coefficient difference
delta_epsilon. Six different situations, called "phases", characterized by the
two X^2 test statistics in allele-based and genotype-based test, are well
separated in the phase diagram parameterized by delta and delta_epsilon. For
two major groups of phases, a single parameter theta = tan^-1
(delta/delta_epsilon) is able to achieves an almost perfect phase separation.
We also applied the analytic result to several types of disease models. It is
shown that for dominant and additive models, genotype-based tests are favored
over allele-based tests.


  The chaperonin GroEL-GroES, a machine which helps some proteins to fold,
cycles through a number of allosteric states, the $T$ state, with high affinity
for substrate proteins (SPs), the ATP-bound $R$ state, and the
$R^{\prime\prime}$ ($GroEL-ADP-GroES$) complex. Structures are known for each
of these states. Here, we use a self-organized polymer (SOP) model for the
GroEL allosteric states and a general structure-based technique to simulate the
dynamics of allosteric transitions in two subunits of GroEL and the heptamer.
The $T \to R$ transition, in which the apical domains undergo counter-clockwise
motion, is mediated by a multiple salt-bridge switch mechanism, in which a
series of salt-bridges break and form. The initial event in the $R \to
R^{\prime\prime}$ transition, during which GroEL rotates clockwise, involves a
spectacular outside-in movement of helices K and L that results in K80-D359
salt-bridge formation. In both the transitions there is considerable
heterogeneity in the transition pathways. The transition state ensembles (TSEs)
connecting the $T$, $R$, and $R^{\prime\prime}$ states are broad with the the
TSE for the $T \to R$ transition being more plastic than the $R\to
R^{\prime\prime}$ TSE. The results suggest that GroEL functions as a
force-transmitting device in which forces of about (5-30) pN may act on the SP
during the reaction cycle.


  Protein-DNA interactions are vital for many processes in living cells,
especially transcriptional regulation and DNA modification. To further our
understanding of these important processes on the microscopic level, it is
necessary that theoretical models describe the macromolecular interaction
energetics accurately. While several methods have been proposed, there has not
been a careful comparison of how well the different methods are able to predict
biologically important quantities such as the correct DNA binding sequence,
total binding free energy, and free energy changes caused by DNA mutation. In
addition to carrying out the comparison, we present two important theoretical
models developed initially in protein folding that have not yet been tried on
protein-DNA interactions. In the process, we find that the results of these
knowledge-based potentials show a strong dependence on the interaction distance
and the derivation method. Finally, we present a knowledge-based potential that
gives comparable or superior results to the best of the other methods,
including the molecular mechanics force field AMBER99.


  Membrane and nuclear proteins of poor solubility have been separated by high
resolution two-dimensional (2-D) gel electrophoresis. Isoelectric focusing with
immobilized pH gradients leads to severe quantitative losses of proteins in the
resulting 2-D map, although the resolution is usually high. Protein solubility
could be improved by using denaturing solutions containing various detergents
and chaotropes. Best results were obtained with a denaturing solution
containing urea, thiourea, and detergents (both nonionic and zwitterionic). The
usefulness of thiourea-containing denaturing mixtures is shown for microsomal
and nuclear proteins as well as for tubulin, a protein highly prone to
aggregation.


  Molecular motors power directed transport of cargoes within cells. Even if a
single motor is sufficient to transport a cargo, motors often cooperate in
small teams. We discuss the cooperative cargo transport by several motors
theoretically and explore some of its properties. In particular we emphasize
how motor teams can drag cargoes through a viscous environment.


  In the interest of conservation, the importance of having a large habitat
available for a species is widely known. Here, we introduce a lattice-based
model for a population and look at the importance of fluctuations as well as
that of the population density, particularly with respect to Allee effects. We
examine the model analytically and by Monte Carlo simulations and find that,
while the size of the habitat is important, there exists a critical population
density below which extinction is assured. This has large consequences with
respect to conservation, especially in the design of habitats and for
populations whose density has become small. In particular, we find that the
probability of survival for small populations can be increased by a reduction
in the size of the habitat and show that there exists an optimal size
reduction.


  Changes in gene expression occurring during differentiation of human
monocytes into dendritic cells were studied at the RNA and protein levels.
These studies showed the induction of several gene classes corresponding to
various biological functions. These functions encompass antigen processing and
presentation, cytoskeleton, cell signalling and signal transduction, but also
an increase in mitochondrial function and in the protein synthesis machinery,
including some, but not all, chaperones. These changes put in perspective the
events occurring during this differentiation process. On a more technical
point, it appears that the studies carried out at the RNA and protein levels
are highly complementary.


  Protein splicing is a post-translational autocatalystic excision of internal
protein sequence (intein) with the subsequent ligation of the flanking
polypeptides (exteins). The high specificity of excision ensured by intein
makes it possible to use a phenomenon of protein splicing for the biotechnology
purposes. The aim of this work was optimization of obtaining and purification
of the recombinant human growth hormone using the protein splicing. It was
experimentally demonstrated that the use of modified intein as auto-removal
affine marker makes it possible to perform the rapid and cheap isolation of the
recombinant protein Hgh. Furthermore, this approach allows to obtain the human
growth hormone with native N-terminus, without formyl-metionine. Key words:
intein, human growth hormone, protein splicing


  In the three decades since its introduction, resource selection analysis
(RSA) has become a widespread method for analyzing spatial patterns of animal
relocations obtained from telemetry studies. Recently, mechanistic home range
models have been proposed as an alternative framework for studying patterns of
animal space use. In contrast to RSA models, mechanistic home range models are
derived from underlying mechanistic descriptions of individual movement
behavior and yield spatially-explicit predictions for patterns of animal
space-use. In addition, their mechanistic underpinning means that, unlike RSA,
mechanistic home range models can also be used to predict changes in space-use
following perturbation. In this paper, we develop a formal reconciliation
between these two methods of home range analysis, showing how differences in
the habitat preferences of individuals give rise to spatially-explicit patterns
of space-use. The resulting unified framework combines the simplicity of
resource selection analysis with the spatially-explicit and predictive
capabilities of mechanistic home range models.


  Much information is stored in amino acid composition of protein and base
composition of DNA. We simulated the evolution of amino acid frequencies and
genomic GC content by a linguistic model. It is showed that the evolution of
genetic code determines the evolution of amino acid frequencies and genomic GC
content. We explained the relationships among amino acid frequencies, genomic
GC content and protein length distribution in a unified theoretical framework.
Especially, the simulations of the evolution of amino acid frequencies and the
codon position GC content agree dramatically with the results based on the data
of all known genomes so far. Furthermore, we found that the space of average
protein length in proteome and ratio of amino acid frequencies is useful to
describe the phylogeny and evolution. Amazingly, the dots of all the species in
this space form an evolutionary flow. We believe that the amino acid gain and
loss is motivated by the established pattern of the variation of amino acid
frequencies. The linguistic mechanism is helpful to unveil the origin of the
genetic code.


  Binding and unbinding of ligands to specific sites of a macromolecule are one
of the most elementary molecular interactions inside the cell that embody the
computational processes of biological regulations. The interaction between
transcription factors and the operators of genes and that between ligands and
binding sites of allosteric enzymes are typical examples of such molecular
interactions. In order to obtain the general mathematical framework of
biological regulations, we formulate these interactions as finite Markov
processes and establish a computational theory of regulatory activities of
macromolecules based mainly on graphical analysis of their state transition
diagrams. The contribution is summarized as follows: (1) Stochastic
interpretation of Michaelis-Menten equation is given. (2) Notion of
\textit{probability flow} is introduced in relation to detailed balance. (3) A
stochastic analogy of \textit{Wegscheider condition} is given in relation to
loops in the state transition diagram. (4) A simple graphical method of
computing the regulatory activity in terms of ligands' concentrations is
obtained for Wegscheider cases.


  We use the Chowdhury ecosystem model, one of the most complex agent-based
ecological models, to test the latitude-niche breadth hypothesis, with regard
to habitat width, i.e., whether tropical species generally have narrower
habitats than high latitude ones. Application of the model has given realistic
results in previous studies on latitudinal gradients in species diversity and
Rapoport's rule. Here we show that tropical species with sufficient vagility
and time to spread into adjacent habitats, tend to have wider habitats than
high latitude ones, contradicting the latitude-niche breadth hypothesis.


  When beneficial mutations are rare, they accumulate by a series of selective
sweeps. But when they are common, many beneficial mutations will occur before
any can fix, so there will be many different mutant lineages in the population
concurrently. In an asexual population, these different mutant lineages
interfere and not all can fix simultaneously. In addition, further beneficial
mutations can accumulate in mutant lineages while these are still a minority of
the population. In this paper, we analyze the dynamics of such multiple
mutations and the interplay between multiple mutations and interference between
clones. These result in substantial variation in fitness accumulating within a
single asexual population. The amount of variation is determined by a balance
between selection, which destroys variation, and beneficial mutations, which
create more. The behavior depends in a subtle way on the population parameters:
the population size, the beneficial mutation rate, and the distribution of the
fitness increments of the potential beneficial mutations. The
mutation-selection balance leads to a continually evolving population with a
steady-state fitness variation. This variation increases logarithmically with
both population size and mutation rate and sets the rate at which the
population accumulates beneficial mutations, which thus also grows only
logarithmically with population size and mutation rate. These results imply
that mutator phenotypes are less effective in larger asexual populations. They
also have consequences for the advantages (or disadvantages) of sex via the
Fisher-Muller effect; these are discussed briefly.


  Despite their topological complexity almost all functional properties of
metabolic networks can be derived from steady-state dynamics. Indeed, many
theoretical investigations (like flux-balance analysis) rely on extracting
function from steady states. This leads to the interesting question, how
metabolic networks avoid complex dynamics and maintain a steady-state behavior.
Here, we expose metabolic network topologies to binary dynamics generated by
simple local rules. We find that the networks' response is highly specific:
Complex dynamics are systematically reduced on metabolic networks compared to
randomized networks with identical degree sequences. Already small topological
modifications substantially enhance the capacity of a network to host complex
dynamic behavior and thus reduce its regularizing potential. This exceptionally
pronounced regularization of dynamics encoded in the topology may explain, why
steady-state behavior is ubiquitous in metabolism.


  We study a classical two-state stochastic system in a sea of substrates and
products (absorbing states), which can be interpreted as a single
Michaelis-Menten catalyzing enzyme or as a channel on a cell surface. We
introduce a novel general method and use it to derive the expression for the
full counting statistics of transitions among the absorbing states. For the
evolution of the system under a periodic perturbation of the kinetic rates, the
latter contains a term with a purely geometrical (the Berry phase)
interpretation. This term gives rise to a pump current between the absorbing
states, which is due entirely to the stochastic nature of the system. We
calculate the first two cumulants of this current, and we argue that it is
observable experimentally.


  Accumulating observational evidence suggests an intimate connection between
rapidly expanding insect populations, deforestation, and global climate change.
We review the evidence, emphasizing the vulnerability of key planetary carbon
pools, especially the Earth's forests that link the micro-ecology of insect
infestation to climate. We survey current research regimes and insect control
strategies, concluding that at present they are insufficient to cope with the
problem's present regional scale and its likely future global scale. We propose
novel bioacoustic interactions between insects and trees as key drivers of
infestation population dynamics and the resulting wide-scale deforestation. The
bioacoustic mechanisms suggest new, nontoxic control interventions and
detection strategies.


  The near-native free energy landscape of protein G is investigated through
0.4 microseconds-long atomistic molecular dynamics simulations in explicit
solvent. A theoretical and computational framework is used to assess the
time-dependence of salient thermodynamical features. While the quasi-harmonic
character of the free energy is found to degrade in a few ns, the slow modes
display a very mild dependence on the trajectory duration. This property
originates from a striking self-similarity of the free energy landscape
embodied by the consistency of the principal directions of the local minima,
where the system dwells for several ns, and of the virtual jumps connecting
them.


  The linear instability of Lotka-Volterra orbits in the homogenous manifold of
a two-patch system is analyzed. The origin of these orbits instability in the
absence of prey migration is revealed to be the dependence of the angular
velocity on the azimuthal angle; in particular, the system desynchronizes at
the exit from the slow part of the trajectory. Using this insight, an analogous
model of a two coupled oscillator is presented and shown to yield the same type
of linear instability. This unables one to incorporate the linear instability
within a recently presented general framework that allows for comparison of all
known stabilization mechanisms and for simple classification of observed
oscillations.


  We analyze several florae (collections of plant species populating specific
areas) in different geographic and climatic regions. For every list of species
we produce a taxonomic classification tree and we consider its statistical
properties. We find that regardless of the geographical location, the climate
and the environment all species collections have universal statistical
properties that we show to be also robust in time. We then compare observed
data sets with simulated communities obtained by randomly sampling a large pool
of species from all over the world. We find differences in the behavior of the
statistical properties of the corresponding taxonomic trees. Our results
suggest that it is possible to distinguish quantitatively real species
assemblages from random collections and thus demonstrate the existence of
correlations between species.


  Populations evolving under the joint influence of recombination and
resampling (traditionally known as genetic drift) are investigated. First, we
summarise and adapt a deterministic approach, as valid for infinite
populations, which assumes continuous time and single crossover events. The
corresponding nonlinear system of differential equations permits a closed
solution, both in terms of the type frequencies and via linkage disequilibria
of all orders. To include stochastic effects, we then consider the
corresponding finite-population model, the Moran model with single crossovers,
and examine it both analytically and by means of simulations. Particular
emphasis is on the connection with the deterministic solution. If there is only
recombination and every pair of recombined offspring replaces their pair of
parents (i.e., there is no resampling), then the {\em expected} type
frequencies in the finite population, of arbitrary size, equal the type
frequencies in the infinite population. If resampling is included, the
stochastic process converges, in the infinite-population limit, to the
deterministic dynamics, which turns out to be a good approximation already for
populations of moderate size.