  Nucleosomes organize the folding of DNA into chromatin and significantly
influence transcription, replication, regulation and repair. All atom molecular
dynamics simulations of a nucleosome and of its 146 basepairs of DNA free in
solution have been conducted. DNA helical parameters are extracted from each
trajectory to compare the conformation, effective force constants, persistence
length measures, and fluctuations of nucleosomal DNA to free DNA. A method for
disassembling and reconstructing the conformation and dynamics of the
nucleosome using Fourier analysis is presented. Results indicate that the
superhelical path of DNA in the nucleosome is irregular. Long length variations
in the conformation of nucleosomal DNA are identified other than those
associated with helix repeat. These variations are required to create a
proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns
of the nuclesomal superhelix. Free DNA achieves enough bend and shear in
solution to create an ideal nucleosome superhelix, but these deformations are
not organized so the conformation is essentially linear. Reconstruction of free
DNA using selected long wavelength variations in conformation can produce
either a left-handed or a right-handed superhelix. DNA is less flexible in the
nucleosome than when free in solution, however such measures are length scale
dependent.


  Many genes have been identified as driving cellular differentiation, but
because of their complex interactions, the understanding of their collective
behaviour requires mathematical modelling. Intriguingly, it has been observed
in numerous developmental contexts, and particularly hematopoiesis, that genes
regulating differentiation are initially co-expressed in progenitors despite
their antagonism, before one is upregulated and others downregulated. We
characterise conditions under which 3 classes of generic "master regulatory
networks", modelled at the molecular level after experimentally-observed
interactions (including bHLH protein dimerisation), and including an arbitrary
number of antagonistic components, can behave as a "multi-switch", directing
differentiation in an all-or-none fashion to a specific cell-type chosen among
more than 2 possible outcomes. bHLH dimerisation networks can readily display
coexistence of many antagonistic factors when competition is low (a simple
characterisation is derived). Decision-making can be forced by a transient
increase in competition, which could correspond to some unexplained
experimental observations related to Id proteins; the speed of response varies
with the initial conditions the network is subjected to, which could explain
some aspects of cell behaviour upon reprogramming. The coexistence of
antagonistic factors at low levels, early in the differentiation process or in
pluripotent stem cells, could be an intrinsic property of the interaction
between those factors, not requiring a specific regulatory system.


  We propose and test a pair potential that is accurate at all relevant
distances and simple enough for use in large-scale computer simulations. A
combination of the Rydberg potential from spectroscopy and the London
inverse-sixth-power energy, the proposed form fits spectroscopically determined
potentials better than the Morse, Varnshi, and Hulburt-Hirschfelder potentials
and much better than the Lennard-Jones and harmonic potentials. At long
distances, it goes smoothly to the correct London force appropriate for gases
and preserves van der Waals's "continuity of the gas and liquid states," which
is routinely violated by coefficients assigned to the Lennard-Jones 6-12 form.


  We derive an analytic expression for site-specific stationary distributions
of amino acids from the Structurally Constrained Neutral (SCN) model of protein
evolution with conservation of folding stability. The stationary distributions
that we obtain have a Boltzmann-like shape, and their effective temperature
parameter, measuring the limit of divergent evolutionary changes at a given
site, can be predicted from a site-specific topological property, the principal
eigenvector of the contact matrix of the native conformation of the protein.
These analytic results, obtained without free parameters, are compared with
simulations of the SCN model and with the site-specific amino acid
distributions obtained from the Protein Data Bank. These results also provide
new insights into how the topology of a protein fold influences its
designability, i.e. the number of sequences compatible with that fold. The
dependence of the effective temperature on the principal eigenvector decreases
for longer proteins, a possible consequence of the fact that selection for
thermodynamic stability becomes weaker in this case.


  Many modified genetic codes are found in specific genomes in which one or
more codons have been reassigned to a different amino acid from that in the
canonical code. We present a model that unifies four possible mechanisms for
reassignment, based on the observation that reassignment involves a gain and a
loss. The loss could be the deletion or loss of function of a tRNA or release
factor. The gain could be the gain of a new type of tRNA for the reassigned
codon, or the gain of function of an existing tRNA due to a mutation or a base
modification. In the codon disappearance mechanism, the codon disappears from
the genome during the period of reassignment. In the other mechanisms, the
codon does not disappear. In the ambiguous intermediate mechanism, the gain
precedes the loss; in the unassigned codon mechanism, the loss precedes the
gain; and in the compensatory change mechanism, the loss and gain spread
through the population simultaneously. We present simulations of the gain-loss
model and demonstrate that all four mechanisms are possible. The frequencies of
the different mechanisms are influenced by selection strengths, number of
codons undergoing reassignment, directional mutation pressure and the
possibility of selection for reduced genome size.


  We showed in this paper that similarity network can be used as an powerful
tools to study the relationship of tRNA genes. We constructed a network of 3719
tRNA gene sequences using simplest alignment and studied its topology, degree
distribution and clustering coefficient. It is found that the behavior of the
network shift from fluctuated distribution to scale-free distribution when the
similarity degree of the tRNA gene sequences increase. tRNA gene sequences with
the same anticodon identity are more self-organized than the tRNA gene
sequences with different anticodon identities and form local clusters in the
network. An interesting finding in our studied is some vertices of the local
cluster have a high connection with other local clusters, the probable reason
is given. Moreover, a network constructed by the same number of random tRNA
sequences is used to make comparisons. The relationships between properties of
the tRNA similarity network and the characters of tRNA evolutionary history are
discussed.


  A vast variety of biological, social, and economical networks shows
topologies drastically differing from random graphs; yet the quantitative
characterization remains unsatisfactory from a conceptual point of view.
Motivated from the discussion of small scale-free networks, a biased link
distribution entropy is defined, which takes an extremum for a power law
distribution. This approach is extended to the node-node link
cross-distribution, whose nondiagonal elements characterize the graph structure
beyond link distribution, cluster coefficient and average path length. From
here a simple (and computationally cheap) complexity measure can be defined.
This Offdiagonal Complexity (OdC) is proposed as a novel measure to
characterize the complexity of an undirected graph, or network. While both for
regular lattices and fully connected networks OdC is zero, it takes a
moderately low value for a random graph and shows high values for apparently
complex structures as scale-free networks and hierarchical trees. The
Offdiagonal Complexity apporach is applied to the Helicobacter pylori protein
interaction network and randomly rewired surrogates.


  Evolving genomes increase a number of their genes by gene duplications. To
escape degradation in a functionless pseudogene, any gene duplicate needs to be
guarded by negative (purifying) selection from otherwise inevitable fixation of
degenerative mutations. In the present study we focus on the evolutionary stage
at which new duplicates come under such surveillance.
  Our analyses of several genomes indicate that in about 10% gene pairs,
selection begins to guard a new gene copy very soon after a duplication event
whereas the vast majority (90%) of extra genes remain redundant and
unrecognised by selection. Such duplicates accumulate all mutations (including
degenerative) in neutral fashion and are actually destined to become
pseudogenes. We revealed this "two-stream" evolutionary pattern by the analysis
of mutations in 2nd versus 3rd codon positions but not by the routinely used
ratio of amino acid replacements (R) versus silent substitutions (S), i.e. the
'2nd vs. 3rd' metric proved to be more resolving than the traditional 'R vs. S'
one for distinguishing neutrally evolving future pseudogenes from their
functional counterparts controlled by negative selection.
  In gene databases for large genomes, hundreds of future pseudogenes are
annotated as functional genes because they do look like intact and valuable by
standard criteria, including even active transcription and translation.
Apparently, these "pseudogenes-to-be" over-cloud and mimic those very
infrequent gene duplicates with increased sequence evolution rates driven by
positive selection.


  Robust systems, like the molecular networks of living cells are often
resistant to single hits such as those caused by high-specificity drugs. Here
we show that partial weakening of the Escherichia coli and Saccharomyces
cerevisiae transcriptional regulatory networks at a small number (3-5) selected
nodes can have a greater impact than the complete elimination of a single
selected node. In both cases, the targeted nodes have the greatest possible
impact; still the results suggest that in some cases broad specificity
compounds or multitarget drug therapies may be more effective than individual
high-affinity, high-specificity ones. Multiple but partial attacks mimic well a
number of in vivo scenarios and may be useful in the efficient modification of
other complex systems.


  When conducting inter-species regression analyses, the phylogenetic
relationships between the individual species need to be taken into account. In
this paper, a procedure for conducting such analyses is proposed, which only
requires the use of a measure of relationship between pairs of species, rather
than a complete phylogeny, and which at the same time assesses the importance
to be attached to the relationships with regard to the conclusions reached. The
procedure is applied to previous data, relating testis size to mean hind tibia
length, duct length and spermathecal area in 15 species of Scathophagidae.


  Gene expression analysis by means of microarrays is based on the sequence
specific binding of mRNA to DNA oligonucleotide probes and its measurement
using fluorescent labels. The binding of RNA fragments involving other
sequences than the intended target is problematic because it adds a "chemical
background" to the signal, which is not related to the expression degree of the
target gene. The paper presents a molecular signature of specific and non
specific hybridization with potential consequences for gene expression
analysis. We analyzed the signal intensities of perfect match (PM) and mismatch
(MM) probes of GeneChip microarrays to specify the effect of specific and non
specific hybridization. We found that these events give rise to different
relations between the PM and MM intensities as function of the middle base of
the PMs, namely a triplet- (C>G=T>A>0) and a duplet-like (C=T>0>G=A) pattern of
the PM-MM log-intensity difference upon binding of specific and non specific
RNA fragments, respectively. The systematic behaviour of the intensity
difference can be rationalized on the level of base pairings of DNA/RNA
oligonucleotide duplexes in the middle of the probe sequence. Non-specific
binding is characterized by the reversal of the central Watson Crick (WC)
pairing for each PM/MM probe pair, whereas specific binding refers to the
combination of a WC and a self complementary (SC) pairing in PM and MM probes,
respectively. The intensity of complementary MM introduces a systematic source
of variation which decreases the precision of expression measures based on the
MM intensities.


  Food webs are networks describing who is eating whom in an ecological
community. By now it is clear that many aspects of food-web structure are
reproducible across diverse habitats, yet little is known about the driving
force behind this structure. Evolutionary and population dynamical mechanisms
have been considered. We propose a model for the evolutionary dynamics of
food-web topology and show that it accurately reproduces observed food-web
characteristic in the steady state. It is based on the observation that most
consumers are larger than their resource species and the hypothesis that
speciation and extinction rates decrease with increasing body mass. Results
give strong support to the evolutionary hypothesis.


  We present a new lattice model for proteins that incorporates a tube-like
anisotropy by introducing a preference for mutually parallel alignments in the
conformations. The model is demonstrated to capture many aspects of real
proteins.


  Electrostatic interactions play a central role in the assembly of
single-stranded RNA viruses. Under physiological conditions of salinity and
acidity, virus capsid assembly requires the presence of genomic material that
is oppositely charged to the core proteins. In this paper we apply basic
polymer physics and statistical mechanics methods to the self-assembly of a
synthetic virus encapsidating generic polyelectrolyte molecules. We find that
(i) the mean concentration of the encapsidated polyelectrolyte material depends
on the surface charge density, the radius of the capsid, and the linear charge
density of the polymer but neither on the salt concentration or the Kuhn
length, (ii) the total charge of the capsid interior is equal but opposite to
that of the empty capsid, a form of charge reversal. Unlike natural viruses,
synthetic viruses are predicted not to be under an osmotic swelling pressure.
The design condition that self-assembly only produces filled capsids is shown
to coincide with the condition that the capsid surface charge exceeds the
desorption threshold of polymer surface adsorption. We compare our results with
studies on the self-assembly of both synthetic and natural viruses.


  Coding information is the main source of heterogeneity (non-randomness) in
the sequences of bacterial genomes. This information can be naturally modeled
by analysing cluster structures in the "in-phase" triplet distributions of
relatively short genomic fragments (200-400bp). We found a universal 7-cluster
structure in bacterial genomic sequences and explained its properties. We show
that codon usage of bacterial genomes is a multi-linear function of their
genomic G+C-content with high accuracy. Based on the analysis of 143 completely
sequenced bacterial genomes available in Genbank in August 2004, we show that
there are four "pure" types of the 7-cluster structure observed. All 143
cluster animated 3D-scatters are collected in a database and is made available
on our web-site: http://www.ihes.fr/~zinovyev/7clusters The finding can be
readily introduced into any software for gene prediction, sequence alignment or
bacterial genomes classification.


  We investigate self-moving particles which prefer to hop with a certain
turning angle equally distributed to the right or left. We assume this turning
angle distribution to be given by a double Gaussian distribution. Based on the
model of Active Brownian particles and we calculate the diffusion coefficient
in dependence on the mean and the dispersion of the turning angles. It is shown
that bounded distribution of food in patches will be optimally consumed by the
objects if they hop preferably with a given angle and not straight forwardly.


  A novel information-theoretic method for reconstruction of interaction
networks is introduced. We prove that the method is exact for some class of
networks. Performance tests on large synthetic transcriptional regulatory
networks produce very encouraging results.


  Background: Elucidating gene regulatory networks is crucial for understanding
normal cell physiology and complex pathologic phenotypes. Existing
computational methods for the genome-wide ``reverse engineering'' of such
networks have been successful only for lower eukaryotes with simple genomes.
Here we present ARACNE, a novel algorithm, using microarray expression
profiles, specifically designed to scale up to the complexity of regulatory
networks in mammalian cells, yet general enough to address a wider range of
network deconvolution problems. This method uses an information theoretic
approach to eliminate the majority of indirect interactions inferred by
co-expression methods.
  Results: We prove that ARACNE reconstructs the network exactly
(asymptotically) if the effect of loops in the network topology is negligible,
and we show that the algorithm works well in practice, even in the presence of
numerous loops and complex topologies. We assess ARACNE's ability to
reconstruct transcriptional regulatory networks using both a realistic
synthetic dataset and a microarray dataset from human B cells. On synthetic
datasets ARACNE achieves very low error rates and outperforms established
methods, such as Relevance Networks and Bayesian Networks. Application to the
deconvolution of genetic networks in human B cells demonstrates ARACNE's
ability to infer validated transcriptional targets of the c MYC proto-oncogene.
We also study the effects of mis estimation of mutual information on network
reconstruction, and show that algorithms based on mutual information ranking
are more resilient to estimation errors.


  We implement the replica exchange molecular dynamics algorithm to study the
interactions of a model peptide (WALP-16) with an explicitly represented DPPC
membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16
into the DPPC bilayer and its folding into an a-helix with a trans-bilayer
orientation. We observe that the insertion of the peptide into the DPPC bilayer
precedes secondary structure formation. Although the peptide has some
propensity to form a partially helical structure in the interfacial region of
the DPPC/water system, this state is not a productive intermediate but rather
an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that
the observed insertion/folding pathway mirrors the potential of mean force
(PMF). Calculation of the enthalpic and entropic contributions to this PMF show
that the surface bound conformation of WALP-16 is significantly lower in energy
than other conformations, and that the insertion of WALP-16 into the bilayer
without regular secondary structure is enthalpically unfavorable by 5-10
kcal/mol/residue. The observed insertion/folding pathway disagrees with the
dominant conceptual model, which is that a surface bound helix is an obligatory
intermediate for the insertion of a-helical peptides into lipid bilayers. In
our simulations, the observed insertion/folding pathway is favored because of a
large (> 100 kcal/mol) increase in system entropy that occurs when the
unstructured WALP-16 peptide enters the lipid bilayer interior. The
insertion/folding pathway that is lowest in free energy depends sensitively on
the near cancellation of large enthalpic and entropic terms. This suggests that
intrinsic membrane peptides may have a diversity of insertion/folding behaviors
depending on the exact system of peptide and lipid under consideration.


  We build upon our previous analytical results for the Penna model of
senescence to include positive mutations. We investigate whether a small but
non-zero positive mutation rate gives qualitatively different results to the
traditional Penna model in which no positive mutations are considered. We find
that the high-lifespan tail of the distribution is radically changed in
structure, but that there is not much effect on the bulk of the population. Th
e mortality plateau that we found previously for a stochastic generalization of
the Penna model is stable to a small positive mutation rate.


  Cellular phenotypes are determined by the dynamical activity of networks of
co-regulated genes. Elucidating such networks is crucial for the understanding
of normal cell physiology as well as for the dissection of complex pathologic
phenotypes. Existing methods for such "reverse engineering" of genetic networks
from microarray expression data have been successful only in prokaryotes (E.
coli) and lower eukaryotes (S. cerevisiae) with relatively simple genomes.
Additionally, they have mostly attempted to reconstruct average properties
about the network connectivity without capturing the highly conditional nature
of the interactions. In this paper we extend the ARACNE algorithm, which we
recently introduced and successfully applied to the reconstruction of
whole-genome transcriptional networks from mammalian cells, precisely to link
the existence of specific network structures to the expression or lack thereof
of specific regulator genes. This is accomplished by analyzing thousands of
alternative network topologies generated by constraining the data set on the
presence or absence of putative regulator genes. By considering interactions
that are consistently supported across several such constraints, we identify
many transcriptional interactions that would not have been detectable by the
original method. By selecting genes that produce statistically significant
changes in network topology, we identify novel candidate regulator genes.
Further analysis shows that transcription factors, kinases, phosphatases, and
other gene families known to effect biochemical interactions, are significantly
overrepresented among the set of candidate regulator genes identified in
silico, indirectly supporting the validity of the approach.


  We enumerate the number of RNA contact structures according to their genus,
i.e. the topological character of their pseudoknots. By using a recently
proposed matrix model formulation for the RNA folding problem, we obtain exact
results for the simple case of an RNA molecule with an infinitely flexible
backbone, in which any arbitrary pair of bases is allowed. We analyze the
distribution of the genus of pseudoknots as a function of the total number of
nucleotides along the phosphate-sugar backbone.


  Analysis of data from an Affymetrix Latin Square spike-in experiment
indicates that measured fluorescence intensities of features on an
oligonucleotide microarray are related to spike-in RNA target concentrations
via a hyperbolic response function, generally identified as a Langmuir
adsorption isotherm. Furthermore the asymptotic signal at high spike-in
concentrations is almost invariably lower for a mismatch feature than for its
partner perfect match feature. We survey a number of theoretical adsorption
models of hybridization at the microarray surface and find that in general they
are unable to explain the differing saturation responses of perfect and
mismatch features. On the other hand, we find that a simple and consistent
explanation can be found in a model in which equilibrium hybridization followed
by partial dissociation of duplexes during the post-hybridization washing
phase.


  Recent genomic analyses on the cellular metabolic network show that reaction
flux across enzymes are diverse and exhibit power-law behavior in its
distribution. While one may guess that the reactions with larger fluxes are
more likely to be lethal under the blockade of its catalyzing gene products or
gene knockouts, we find, by in silico flux analysis, that the lethality rarely
has correlations with the flux level owing to the widespread backup pathways
innate in the genome-wide metabolism of \textit{Escherichia coli}. Lethal
reactions, of which the deletion generates cascading failure of following
reactions up to the biomass reaction, are identified in terms of the Boolean
network scheme as well as the flux balance analysis. The avalanche size of a
reaction, defined as the number of subsequently blocked reactions after its
removal, turns out to be a useful measure of lethality. As a means to elucidate
phenotypic robustness to a single deletion, we investigate synthetic lethality
in reaction level, where simultaneous deletion of a pair of nonlethal reactions
leads to the failure of the biomass reaction. Synthetic lethals identified via
flux balance and Boolean scheme are consistently shown to act in parallel
pathways, working in such a way that the backup machinery is compromised.


  In 1995 T.J.Penna introduced a simple model of biological aging. A modified
Penna model has been demonstrated to exhibit behaviour of real-life systems
including catastrophic senescence in salmon and a mortality plateau at advanced
ages. We present a general steady-state, analytic solution to the Penna model,
able to deal with arbitrary birth and survivability functions. This solution is
employed to solve standard variant Penna models studied by simulation.
Different Verhulst factors regulating both the birth rate and external death
rate are considered.


  We build upon the recent steady-state Penna model solution, Phys.Rev.Lett.
89, 288103 (2002), to study the population dynamics within the Penna model. We
show, that any perturbation to the population can be broken into a collection
of modes each of which decay exponentially with its respective time constant.
The long time behaviour of population is therefore likely to be dominated by
the modes with the largest time constants. We confirm our analytical approach
with simulation data.


  Exact numerical methods and stochastic simulation methods are developed to
study the force stretching single RNA issue on the secondary structure level in
equilibrium. By computing the force-extension curves on the constant force and
the constant extension ensembles, we find the two independent methods agree
with each other quite well. To show the precision of our methods in predicting
unfolding experiments, the unfolding forces of different RNA molecules under
different experimental conditions are calculated. We find that the ionic
corrections on the RNA free energies alone might not account for the apparent
differences between the theoretical calculations and the experimental data; an
ionic correction to the persistent length of single-stranded RNA should be
necessary.


  We analyze a series of publicly available controlled experiments (Latin
square) on Affymetrix high density oligonucleotide microarrays using a simple
physical model of the hybridization process. We plot for each gene the signal
intensity versus the hybridization free energy of RNA/DNA duplexes in solution,
for perfect matching and mismatching probes. Both values tend to align on a
single master curve in good agreement with Langmuir adsorption theory, provided
one takes into account the decrease of the effective target concentration due
to target-target hybridization in solution. We give an example of a deviation
from the expected thermodynamical behavior for the probe set 1091\_at due to
annotation problems, i.e. the surface-bound probe is not the exact complement
of the target RNA sequence, because of errors present in public databases at
the time when the array was designed. We show that the parametrization of the
experimental data with RNA/DNA free energy improves the quality of the fits and
enhances the stability of the fitting parameters compared to previous studies.


  A prominent feature of gene transcription regulatory networks is the presence
in large numbers of motifs, i.e, patterns of interconnection, in the networks.
One such motif is the feed forward loop (FFL) consisting of three genes X, Y
and Z. The protein product of x of X controls the synthesis of protein product
y of Y. Proteins x and y jointly regulate the synthesis of z proteins from the
gene Z. The FFLs, depending on the nature of the regulating interactions, can
be of eight different types which can again be classified into two categories:
coherent and incoherent. In this paper, we study the noise characteristics of
FFLs using the Langevin formalism and the Monte Carlo simulation technique
based on the Gillespie algorithm. We calculate the variances around the mean
protein levels in the steady states of the FFLs and find that, in the case of
coherent FFLs, the most abundant FFL, namely, the Type-1 coherent FFL, is the
least noisy. This is however not so in the case of incoherent FFLs. The results
suggest possible relationships between noise, functionality and abundance.


  Background: Diagnostic DNA arrays for detection of point mutations as markers
for cancer usually function in the presence of a large excess of wild type DNA.
This excess can give rise to false positives due to competitive hybridization
of the wild type target at the mutation spot. The analysis of the DNA array
data is typically qualitative aiming to establish the presence or absence of a
particular point mutation. Our theoretical approach yields methods for
quantifying the analysis so as to obtain the ratio of concentrations of mutated
and wild type DNA. Method: The theory is formulated in terms of the
hybridization isotherms relating the hybridization fraction at the spot to the
composition of the sample solutions at thermodynamic equilibrium. It focuses on
samples containing an excess of single stranded DNA and on DNA arrays with low
surface density of probes. The hybridization equilibrium constants can be
obtained by the nearest neighbor method. Results: Two approaches allow us to
obtain quantitative results from the DNA array data. In one the signal of the
mutation spot is compared with that of the wild type spot. The implementation
requires knowledge of the saturation intensity of the two spots. The second
approach requires comparison of the intensity of the mutation spot at two
different temperatures. In this case knowledge of the saturation signal is not
always necessary. Conclusions: DNA arrays can be used to obtain quantitative
results on the concentration ratio of mutated DNA to wild type DNA in studies
of somatic point mutations.


  An oscillation with a period of around 500 kb in guanine and cytosine content
(GC%) is observed in the DNA sequence of human chromosome 21. This oscillation
is localized in the rightmost one-eighth region of the chromosome, from 43.5 Mb
to 46.5 Mb. Five cycles of oscillation are observed in this region with six
GC-rich peaks and five GC-poor valleys. The GC-poor valleys comprise regions
with low density of CpG islands and, alternating between the two DNA strands,
low gene density regions. Consequently, the long-range oscillation of GC%
result in spacing patterns of both CpG island density, and to a lesser extent,
gene densities.


  Spatial fluctuations of guanine and cytosine base content (GC%) are studied
by spectral analysis for the complete set of human genomic DNA sequences. We
find that (i) the 1/f^alpha decay is universally observed in the power spectra
of all twenty-four chromosomes, and that (ii) the exponent alpha \approx 1
extends to about 10^7 bases, one order of magnitude longer than what has
previously been observed. We further find that (iii) almost all human
chromosomes exhibit a cross-over from alpha_1 \approx 1 (1/f^alpha_1) at lower
frequency to alpha_2 < 1 (1/f^alpha_2) at higher frequency, typically occurring
at around 30,000--100,000 bases, while (iv) the cross-over in this frequency
range is virtually absent in human chromosome 22. In addition to the universal
1/f^alpha noise in power spectra, we find (v) several lines of evidence for
chromosome-specific correlation structures, including a 500,000 bases long
oscillation in human chromosome 21. The universal 1/f^alpha spectrum in human
genome is further substantiated by a resistance to variance reduction in
guanine and cytosine content when the window size is increased.


  We study global fluctuations of the guanine and cytosine base content (GC%)
in mouse genomic DNA using spectral analyses. Power spectra S(f) of GC%
fluctuations in all nineteen autosomal and two sex chromosomes are observed to
have the universal functional form S(f) \sim 1/f^alpha (alpha \approx 1) over
several orders of magnitude in the frequency range 10^-7< f < 10^-5 cycle/base,
corresponding to long-ranging GC% correlations at distances between 100 kb and
10 Mb. S(f) for higher frequencies (f > 10^-5 cycle/base) shows a flattened
power-law function with alpha < 1 across all twenty-one chromosomes. The
substitution of about 38% interspersed repeats does not affect the functional
form of S(f), indicating that these are not predominantly responsible for the
long-ranged multi-scale GC% fluctuations in mammalian genomes. Several
biological implications of the large-scale GC% fluctuation are discussed,
including neutral evolutionary history by DNA duplication, chromosomal bands,
spatial distribution of transcription units (genes), replication timing, and
recombination hot spots.


  We present a vector field method for obtaining the spatial organization of 3D
patterns of gene expression based on gradients and lines of force obtained by
numerical integration. The convergence of these lines of force in local maxima
are centers of gene expression, providing a natural and powerful framework to
characterize the organization and dynamics of biological structures. We apply
this novel methodology to analyze the expression pattern of the Enhanced Green
Fluorescent Protein (EGFP) driven by the promoter of light chain myosin II
during zebrafish heart formation.


  The activities and interactions of proteins that govern the cellular response
to a signal generate a multitude of protein phosphorylation states and
heterogeneous protein complexes. Here, using a computational model that
accounts for 307 molecular species implied by specified interactions of four
proteins involved in signalling by the immunoreceptor Fc$\epsilon$RI, we
determine the relative importance of molecular species that can be generated
during signalling, chemical transitions among these species, and reaction paths
that lead to activation of the protein tyrosine kinase (PTK) Syk. By all of
these measures and over 2- and 10-fold ranges of model parameters--rate
constants and initial concentrations--only a small portion of the biochemical
network is active. The spectrum of active complexes, however, can be shifted
dramatically, even by a change in the concentration of a single protein, which
suggests that the network can produce qualitatively different responses under
different cellular conditions and in response to different inputs. Reduced
models that reproduce predictions of the full model for a particular set of
parameters lose their predictive capacity when parameters are varied over
2-fold ranges.


  Comparison of properties of three domains of titin, I1, I27 and I28, in a
simple geometry-based model shows that despite a high structural homology
between their native states different domains show similar but distinguishable
mechanical properties. Folding properties of the separate domains are predicted
to be diversified which reflects sensitivity of the kinetics to the details of
native structures. The Go-like model corresponding to the experimentally
resolved native structure of the I1 domain is found to provide the biggest
thermodynamic and mechanical stability compared to the other domains studied
here. We analyze elastic, thermodynamic and kinetic properties of several
structures corresponding to the I28 domain as obtained through homology-based
modeling. We discuss the ability of the models of the I28 domain to reproduce
experimental results qualitatively. A strengthening of contacts that involve
hydrophobic amino acids does not affect theoretical comparisons of the domains.
Tandem linkages of up to five identical or different domains unravel in a
serial fashion at low temperatures. We study the nature of the intermediate
state that arises in the early stages of the serial unraveling and find it to
qualitatively agree with the results of Marszalek et al.


  A model for the processive movement of dynein is presented based on
experimental observations available. In the model, the change from strong
microtubule-binding to weak binding of dynein is determined naturally by the
variation of the relative orientation between the two interacting surfaces of
the stalk tip and the microtubule as the stalk rotates from the ADP.Vi-state
orientation to the apo-state orientation. This means that the puzzling
communication from the ATP binding site in the globular head to the MT-binding
site in the tip of the stalk, which is prerequisite in the conventional model,
is not required. Using the present model, the previous experimental results,
such as (i) the step size of a dynein being an integer times of the period of
the MT lattice, (ii) the dependence of the step size on load, i.e., the step
size decreasing with the increase of load, and (iii) the stall force being
proportional to [ATP] at low [ATP] and becoming saturated at high [ATP], are
well explained.


  Gene sequences in the vicinity of splice sites are found to possess
dinucleotide periodicities, especially RR and YY, with the period close to the
pitch of nucleosome DNA. This confirms previously reported finding about
preferential positioning of splice junctions within the nucleosomes. The RR and
YY dinucleotides oscillate counterphase, i.e., their respective preferred
positions are shifted about half-period one from another, as it was observed
earlier for AA and TT dinucleotides. Species specificity of nucleosome
positioning DNA pattern is indicated by predominant use of the periodical
GG(CC) dinucleotides in human and mouse genes, as opposed to predominant AA(TT)
dinucleotides in Arabidopsis and C.elegans.
  Keywords: chromatin; gene splicing; intron; exon; dinucleotide; periodical
pattern


  We address the controversial hot question concerning the validity of the
loose coupling versus the lever-arm theories in the actomyosin dynamics by
re-interpreting and extending the phenomenological washboard potential model
proposed by some of us in a previous paper. In this new model a Brownian motion
harnessing thermal energy is assumed to co-exist with the deterministic swing
of the lever-arm, to yield an excellent fit of the set of data obtained by some
of us on the sliding of Myosin II heads on immobilized actin filaments under
various load conditions. Our theoretical arguments are complemented by accurate
numerical simulations, and the robustness of the model is tested via different
choices of parameters and potential profiles.


  DNA looping participates in transcriptional regulation, for instance, by
allowing distal binding sites to act synergistically. Here we study this
process and compare different regulatory mechanisms based on repression with
and without looping. Within a simple mathematical model for the lac operon, we
show that regulation based on DNA looping, in addition to increasing the
repression level, can reduce the fluctuations of transcription and, at the same
time, decrease the sensitivity to changes in the number of regulatory proteins.
Looping is thus able to circumvent some of the constraints inherent to
mechanisms based solely on binding to a single operator site and provides a
mechanism to regulate not only the average properties of transcription but also
its fluctuations.


  An abstract DNA-type system is defined by a set of nonlinear kinetic
equations with polynomial nonlinearities that admit soliton solutions
associated with helical geometry. The set of equations allows for two different
Lax representations: A von Neumann form and a Darboux-covariant Lax pair. We
explain why non-Kolmogorovian probability models occurring in soliton kinetics
are naturally associated with chemical reactions. The most general known
characterization of soliton kinetic equations is given and a class of explicit
soliton solutions is discussed. Switching between open and closed states is a
generic behaviour of the helices. The effect does not crucially depend on the
order of nonlinearity (i.e. types of reactions), a fact that may explain why
simplified models possess properties occuring in realistic systems. We explain
also why fluctuations based on Darboux transformations will not destroy the
dynamics but only switch between a finite number of helical structures.


  Exploiting recent developments in information theory, we propose, illustrate,
and validate a principled information-theoretic algorithm for module discovery
and resulting measure of network modularity. This measure is an order parameter
(a dimensionless number between 0 and 1). Comparison is made to other
approaches to module-discovery and to quantifying network modularity using
Monte Carlo generated Erdos-like modular networks. Finally, the Network
Information Bottleneck (NIB) algorithm is applied to a number of real world
networks, including the "social" network of coauthors at the APS March Meeting
2004.


  The evolutionary reason for the increase in gene length from archaea to
prokaryotes to eukaryotes observed in large scale genome sequencing efforts has
been unclear. We propose here that the increasing complexity of protein-protein
interactions has driven the selection of longer proteins, as longer proteins
are more able to distinguish among a larger number of distinct interactions due
to their greater average surface area. Annotated protein sequences available
from the SWISS-PROT database were analyzed for thirteen eukaryotes, eight
bacteria, and two archaea species. The number of subcellular locations to which
each protein is associated is used as a measure of the number of interactions
to which a protein participates. Two databases of yeast protein-protein
interactions were used as another measure of the number of interactions to
which each \emph{S. cerevisiae} protein participates. Protein length is shown
to correlate with both number of subcellular locations to which a protein is
associated and number of interactions as measured by yeast two-hybrid
experiments. Protein length is also shown to correlate with the probability
that the protein is encoded by an essential gene. Interestingly, average
protein length and number of subcellular locations are not significantly
different between all human proteins and protein targets of known, marketed
drugs. Increased protein length appears to be a significant mechanism by which
the increasing complexity of protein-protein interaction networks is
accommodated within the natural evolution of species. Consideration of protein
length may be a valuable tool in drug design, one that predicts different
strategies for inhibiting interactions in aberrant and normal pathways.


  We use the lac operon in Escherichia coli as a prototype system to illustrate
the current state, applicability, and limitations of modeling the dynamics of
cellular networks. We integrate three different levels of description
-molecular, cellular, and that of cell population- into a single model, which
seems to capture many experimental aspects of the system.


  A significant part of the thin layers of counter-ions adjacent to the
exterior and interior surfaces of a cell membrane form quasi-two-dimensional
(2D) layers of mobile charge. Collective charge density oscillations, known as
plasmon modes, in these 2D charged systems of counter-ions are predicted in the
present paper. This is based on a calculation of the self-consistent response
of this system to a fast electric field fluctuation. The possibility that the
membrane channels might be using these excitations to carry out fast
communication is suggested and experiments are proposed to reveal the existence
of such excitations.


  When highly active antiretroviral therapy is administered for long periods of
time to HIV-1 infected patients, most patients achieve viral loads that are
``undetectable'' by standard assay (i.e., HIV-1 RNA $ < 50$ copies/ml). Yet
despite exhibiting sustained viral loads below the level of detection, a number
of these patients experience unexplained episodes of transient viremia or viral
"blips". We propose here that transient activation of the immune system by
opportunistic infection may explain these episodes of viremia. Indeed, immune
activation by opportunistic infection may spur HIV replication, replenish viral
reservoirs and contribute to accelerated disease progression. In order to
investigate the effects of concurrent infection on chronically infected HIV
patients under treatment with highly active antiretroviral therapy (HAART), we
extend a simple dynamic model of the effects of vaccination on HIV infection
[Jones and Perelson, JAIDS 31:369-377, 2002] to include growing pathogens. We
then propose a more realistic model for immune cell expansion in the presence
of pathogen, and include this in a set of competing models that allow low
baseline viral loads in the presence of drug treatment. Programmed expansion of
immune cells upon exposure to antigen is a feature not previously included in
HIV models, and one that is especially important to consider when simulating an
immune response to opportunistic infection. Using these models we show that
viral blips with realistic duration and amplitude can be generated by
concurrent infections in HAART treated patients.


  Recently, several groups have used error-prone polymerase chain reactions to
construct mutant libraries containing up to 27 nucleotide mutations per gene on
average, and reported a striking observation: although retention of protein
function initially declines exponentially with mutations as has previously been
observed, orders of magnitude more proteins remain viable at the highest
mutation rates than this trend would predict. Mutant proteins having improved
or novel activity were isolated disproportionately from these heavily mutated
libraries, leading to the suggestion that distant regions of sequence space are
enriched in useful cooperative mutations and that optimal mutagenesis should
target these regions. If true, these claims have profound implications for
laboratory evolution and for evolutionary theory. Here, we demonstrate that
properties of the polymerase chain reaction can explain these results and,
consequently, that average protein viability indeed decreases exponentially
with mutational distance at all error rates. We show that high-error-rate
mutagenesis may be useful in certain cases, though for very different reasons
than originally proposed, and that optimal mutation rates are inherently
protocol-dependent. Our results allow optimal mutation rates to be found given
mutagenesis conditions and a protein of known mutational robustness.


  In this paper we investigate the role of native geometry on the kinetics of
protein folding based on simple lattice models and Monte Carlo simulations.
Results obtained within the scope of the Miyazawa-Jernigan indicate the
existence of two dynamical folding regimes depending on the protein chain
length. For chains larger than 80 amino acids the folding performance is
sensitive to the native state's conformation. Smaller chains, with less than 80
amino acids, fold via two-state kinetics and exhibit a significant correlation
between the contact order parameter and the logarithmic folding times. In
particular, chains with N=48 amino acids were found to belong to two broad
classes of folding, characterized by different cooperativity, depending on the
contact order parameter. Preliminary results based on the G\={o} model show
that the effect of long range contact interaction strength in the folding
kinetics is largely dependent on the native state's geometry.


  RNA editing by members of the double-stranded RNA-specific ADAR family leads
to site-specific conversion of adenosine to inosine (A-to-I) in precursor
messenger RNAs. Editing by ADARs is believed to occur in all metazoa, and is
essential for mammalian development. Currently, only a limited number of human
ADAR substrates are known, while indirect evidence suggests a substantial
fraction of all pre-mRNAs being affected. Here we describe a computational
search for ADAR editing sites in the human transcriptome, using millions of
available expressed sequences. 12,723 A-to-I editing sites were mapped in 1,637
different genes, with an estimated accuracy of 95%, raising the number of known
editing sites by two orders of magnitude. We experimentally validated our
method by verifying the occurrence of editing in 26 novel substrates. A-to-I
editing in humans primarily occurs in non-coding regions of the RNA, typically
in Alu repeats. Analysis of the large set of editing sites indicates the role
of editing in controlling dsRNA stability.


  Monte Carlo simulations show that long-range interactions play a major role
in determining the folding rates of 48-mer three-dimensional lattice polymers
modelled by the Go potential. For three target structures with different native
geometries we found a sharp increase in the folding time when the relative
contribution of the long-range interactions to the native state's energy is
decreased from ~50% towards zero. However, the dispersion of the simulated
folding times depends strongly on the native geometry and Go polymers folding
to one of the target structures exhibit folding times spanning three orders of
magnitude. We have also found that, depending on the target geometry, a strong
geometric coupling may exist between local and long-range contacts meaning
that, when this coupling exists, the formation of long-range contacts is forced
by the previous formation of local contacts. The absence of a strong geometric
coupling leads to kinetics that are more sensitive to the interaction energy
parameters; in this case the formation of local contacts is not sufficient to
promote the establishment of long-range ones when these are strongly penalized
energetically, leading to longer folding times.


  A dynamical picture of phylogenetic evolution is given in terms of Markov
models on a state space, comprising joint probability distributions for
character types of taxonomic classes. Phylogenetic branching is a process which
augments the number of taxa under consideration, and hence the rank of the
underlying joint probability state tensor. We point out the combinatorial
necessity for a second-quantised, or Fock space setting, incorporating discrete
counting labels for taxa and character types, to allow for a description in the
number basis. Rate operators describing both time evolution without branching,
and also phylogenetic branching events, are identified. A detailed development
of these ideas is given, using standard transcriptions from the microscopic
formulation of nonequilibrium reaction-diffusion or birth-death processes.
These give the relations between stochastic rate matrices, the matrix elements
of the corresponding evolution operators representing them, and the integral
kernels needed to implement these as path integrals. The `free' theory (without
branching) is solved, and the correct trilinear `interaction' terms
(representing branching events) are presented. The full model is developed in
perturbation theory via the derivation of explicit Feynman rules which
establish that the probabilities (pattern frequencies of leaf colourations)
arising as matrix elements of the time evolution operator are identical with
those computed via the standard analysis. Simple examples (phylogenetic trees
with 2 or 3 leaves), are discussed in detail. Further implications for the work
are briefly considered including the role of time reparametrisation covariance.


  This paper develops a point-mutation model describing the evolutionary
dynamics of a population of adult stem cells. Such a model may prove useful for
quantitative studies of tissue aging and the emergence of cancer. We consider
two modes of chromosome segregation: (1) Random segregation, where the daughter
chromosomes of a given parent chromosome segregate randomly into the stem cell
and its differentiating sister cell. (2) ``Immortal DNA strand''
co-segregation, for which the stem cell retains the daughter chromosomes with
the oldest parent strands. Immortal strand co-segregation is a mechanism,
originally proposed by Cairns (J. Cairns, {\it Nature} {\bf 255}, 197 (1975)),
by which stem cells preserve the integrity of their genomes. For random
segregation, we develop an ordered strand pair formulation of the dynamics,
analogous to the ordered strand pair formalism developed for quasispecies
dynamics involving semiconservative replication with imperfect lesion repair
(in this context, lesion repair is taken to mean repair of postreplication
base-pair mismatches). Interestingly, a similar formulation is possible with
immortal strand co-segregation, despite the fact that this segregation
mechanism is age-dependent. From our model we are able to mathematically show
that, when lesion repair is imperfect, then immortal strand co-segregation
leads to better preservation of the stem cell lineage than random chromosome
segregation. Furthermore, our model allows us to estimate the optimal lesion
repair efficiency for preserving an adult stem cell population for a given
period of time. For human stem cells, we obtain that mispaired bases still
present after replication and cell division should be left untouched, to avoid
potentially fixing a mutation in both DNA strands.


  Computational RNA secondary structure prediction is rather well established.
However, such prediction algorithms always depend on a large number of
experimentally measured parameters. Here, we study how sensitive structure
prediction algorithms are to changes in these parameters. We find that already
for changes corresponding to the actual experimental error to which these
parameters have been determined 30% of the structure are falsly predicted and
the ground state structure is preserved under parameter perturbation in only 5%
of all cases. We establish that base pairing probabilities calculated in a
thermal ensemble are a viable though not perfect measure for the reliability of
the prediction of individual structure elements. A new measure of stability
using parameter perturbation is proposed, and its limitations discussed.


  Finite size effects on the cooperative thermal denaturation of proteins are
considered. A dimensionless measure of cooperativity, Omega, scales as N^zeta,
where N is the number of amino acids. Surprisingly, we find that zeta is
universal with zeta = 1 + gamma, where the exponent gamma characterizes the
divergence of the susceptibility for a self-avoiding walk. Our lattice model
simulations and experimental data are consistent with the theory. Our finding
rationalizes the marginal stability of proteins and substantiates the earlier
predictions that the efficient folding of two-state proteins requires the
folding transition temperature to be close to the collapse temperature.


  We have analysed the possibility of scaling the sexual Penna ageing model.
Assuming that the number of genes expressed before the reproduction age grows
linearly with the genome size and that the mutation rate per genome and
generation is constant, we have found that the fraction of defective genes
expressed before the minimum reproduction age drops with the genome size, while
the number of defective genes eliminated by the genetic death grows with genome
size. Thus, the evolutionary costs decrease with enlarging the genome. After
rescaling the time scale according to the mutational clock, age distributions
of populations do not depend on the genome size. Nevertheless, enlarging the
genome increases the reproduction potential of populations.


  We show that the protein-protein interaction networks can be surprisingly
well described by a very simple evolution model of duplication and divergence.
The model exhibits a remarkably rich behavior depending on a single parameter,
the probability to retain a duplicated link during divergence. When this
parameter is large, the network growth is not self-averaging and an average
vertex degree increases algebraically. The lack of self-averaging results in a
great diversity of networks grown out of the same initial condition. For small
values of the link retention probability, the growth is self-averaging, the
average degree increases very slowly or tends to a constant, and a degree
distribution has a power-law tail.


  We study the size fluctuations of a local denaturation zone in a DNA molecule
in the presence of proteins that selectively bind to single-stranded DNA, based
on a (2+1)-dimensional master equation. By tuning the physical parameters we
can drive the system from undisturbed bubble fluctuations to full, binding
protein-induced denaturation. We determine the effective free energy landscape
of the DNA-bubble and explore its relaxation modes.


  There are several simple criteria of folding to a native state in model
proteins. One of them involves crossing of a threshold value of the RMSD
distance away from the native state. Another checks whether all native contacts
are established, i.e. whether the interacting amino acids come closer than some
characteristic distance. We use Go-like models of proteins and show that such
simple criteria may prompt one to declare folding even though fragments of the
resulting conformations have a wrong sense of chirality. We propose that a
better condition of folding should augment the simple criteria with the
requirement that most of the local values of the chirality should be nearly
native. The kinetic discrepancy between the simple and compound criteria can be
substantially reduced in the Go-like models by providing the Hamiltonian with a
term which favors native values of the local chirality. We study the effects of
this term as a function of its amplitude and compare it to other models such as
with the side groups and with the angle-dependent potentials.


  In this paper we present continuous age- and space-structured models and
numerical computations of Proteus mirabilis swarm-colony development. We base
the mathematical representation of the cell-cycle dynamics of Proteus mirabilis
on those developed by Esipov and Shapiro, which are the best understood aspects
of the system, and we make minimum assumptions about less-understood
mechanisms, such as precise forms of the spatial diffusion. The models in this
paper have explicit age-structure and, when solved numerically, display both
the temporal and spatial regularity seen in experiments, whereas the Esipov and
Shapiro model, when solved accurately, shows only the temporal regularity.
  The composite hyperbolic-parabolic partial differential equations used to
model Proteus mirabilis swarm-colony development are relevant to other
biological systems where the spatial dynamics depend on local physiological
structure. We use computational methods designed for such systems, with known
convergence properties, to obtain the numerical results presented in this
paper.


  We review and further develop an analytical model that describes how
thermodynamic constraints on the stability of the native state influence
protein evolution in a site-specific manner. To this end, we represent both
protein sequences and protein structures as vectors: Structures are represented
by the principal eigenvector (PE) of the protein contact matrix, a quantity
that resembles closely the effective connectivity of each site; Sequences are
represented through the ``interactivity'' of each amino acid type, using novel
parameters that are correlated with hydropathy scales. These interactivity
parameters are more strongly correlated than the other hydropathy scales that
we examine with: (1) The change upon mutations of the unfolding free energy of
proteins with two-states thermodynamics; (2) Genomic properties as the
genome-size and the genome-wide GC content; (3) The main eigenvectors of the
substitution matrices. The evolutionary average of the interactivity vector
correlates very strongly with the PE of a protein structure. Using this result,
we derive an analytic expression for site-specific distributions of amino acids
across protein families in the form of Boltzmann distributions whose ``inverse
temperature'' is a function of the PE component. We show that our predictions
are in agreement with site-specific amino acid distributions obtained from the
Protein Data Bank, and we determine the mutational model that best fits the
observed site-specific amino acid distributions. Interestingly, the optimal
model almost minimizes the rate at which deleterious mutations are eliminated
by natural selection.


  We determine stability and attractor properties of random Boolean genetic
network models with canalyzing rules for a variety of architectures. For all
power law, exponential, and flat in-degree distributions, we find that the
networks are dynamically stable. Furthermore, for architectures with few inputs
per node, the dynamics of the networks is close to critical. In addition, the
fraction of genes that are active decreases with the number of inputs per node.
These results are based upon investigating ensembles of networks using
analytical methods. Also, for different in-degree distributions, the numbers of
fixed points and cycles are calculated, with results intuitively consistent
with stability analysis; fewer inputs per node implies more cycles, and vice
versa. There are hints that genetic networks acquire broader degree
distributions with evolution, and hence our results indicate that for single
cells, the dynamics should become more stable with evolution. However, such an
effect is very likely compensated for by multicellular dynamics, because one
expects less stability when interactions among cells are included. We verify
this by simulations of a simple model for interactions among cells.


  A spin chain, describing a nucleotides sequence, is identified by the labels
of a vector state of an irreducible representation of U_q->0(sl2). A master
equation for the distribution function is written, where the intensity of the
one-spin flip is assumed to depend on the variation of the labels. The
numerically computed equilibrium distribution is nicely fitted by a Yule
distribution, which is the observed distribution of the ranked short
oligonucleotides frequency in DNA


  A simplified interaction potential for protein folding studies at the atomic
level is discussed and tested on a set of peptides with about 20 residues each.
The test set contains both alpha-helical (Trp cage, Fs) and beta-sheet (GB1p,
GB1m2, GB1m3, Betanova, LLM) peptides. The model, which is entirely
sequence-based, is able to fold these different peptides for one and the same
choice of model parameters. Furthermore, the melting behavior of the peptides
is in good quantitative agreement with experimental data. Apparent folded
populations obtained using different observables are compared, and are found to
be very different for some of the peptides (e.g., Betanova). In other cases (in
particular, GB1m2 and GB1m3), the different estimates agree reasonably well,
indicating a more two-state-like melting behavior.


  Finite size effects on the calorimetric cooperatity of the folding-unfolding
transition in two-state proteins are considered using the Go lattice models
with and without side chains. We show that for models without side chains a
dimensionless measure of calorimetric cooperativity kappa2 defined as the ratio
of the van't Hoff to calorimetric enthalpy does not depend on the number of
amino acids N. The average value of kappa2 is about 3/4 which is lower than the
experimental value kappa2=1. For models with side chains kappa2 approaches
unity as kappa2 \sim N^mu, where exponent mu=0.17. Above the critical chain
length Nc =135 these models can mimic the truly all-or-non folding-unfolding
transition.


  Recent simulations of the solid phase DNA amplification (SPA) by J.-F.
Mercier et al (Biophys. J. 85 (2003) 2075) are generalized to include two kinds
of primers and the off-lattice character of the primer distribution on the
surface. The sigmoidal character of the primer occupation by DNA, observed
experimentally, is reproduced in the simulation. We discuss an influence of two
parameters on the efficience of the amplification process: the initial density
p_0 of the occupied primers from the interfacial amplification and the ratio r
of the molecule length to the average distance between primers. The number of
cycles till the saturation decreases with p_0 roughly as p_0^{-0.26}. For
r=1.5, the number of occupied primers is reduced by a factor two, when compared
to the case of longer molecules. Below r=1.4, the effectivity of SPA is reduced
by a factor 100.


  The study of gene regulation and expression is often discussed in
quantitative terms. In particular, the expression of genes is regularly
characterized with respect to how much, how fast, when and where. Whether
discussing the level of gene expression in a bacterium or its precise location
within a developing embryo, the natural language for these experiments is that
of numbers. Such quantitative data demands quantitative models. We review a
class of models ("thermodynamic models") which exploit statistical mechanics to
compute the probability that RNA polymerase is at the appropriate promoter.
This provides a mathematically precise elaboration of the idea that activators
are agents of recruitment which increase the probability that RNA polymerase
will be found at the promoter of interest. We discuss a framework which
describes the interactions of repressors, activators, helper molecules and RNA
polymerase using the concept of effective concentrations, expressed in terms of
a function we call the "regulation factor". This analysis culminates in an
expression for the probability of RNA polymerase binding at the promoter of
interest as a function of the number of regulatory proteins in the cell. In a
companion paper [1], these ideas are applied to several case studies which
illustrate the use of the general formalism.


  With the increasing amount of experimental data on gene expression and
regulation, there is a growing need for quantitative models to describe the
data and relate them to the different contexts. The thermodynamic models
reviewed in the preceding paper provide a useful framework for the quantitative
analysis of bacterial transcription regulation. We review a number of
well-characterized bacterial promoters that are regulated by one or two species
of transcription factors, and apply the thermodynamic framework to these
promoters. We show that the framework allows one to quantify vastly different
forms of gene expression using a few parameters. As such, it provides a compact
description useful for higher-level studies, e.g., of genetic networks, without
the need to invoke the biochemical details of every component. Moreover, it can
be used to generate hypotheses on the likely mechanisms of transcriptional
control.


  Sequence comparison across multiple organisms aids in the detection of
regions under selection. However, resource limitations require a prioritization
of genomes to be sequenced. This prioritization should be grounded in two
considerations: the lineal scope encompassing the biological phenomena of
interest, and the optimal species within that scope for detecting functional
elements. We introduce a statistical framework for optimal species subset
selection, based on maximizing power to detect conserved sites. In a study of
vertebrate species, we show that the optimal species subset is not in general
the most evolutionarily diverged subset. Our results suggest that marsupials
are prime sequencing candidates.


  Single molecule FRET (fluorescence resonance energy transfer) is a powerful
technique for detecting real-time conformational changes and molecular
interactions during biological reactions. In this review, we examine different
techniques of extending observation times via immobilization and illustrate how
useful biological information can be obtained from single molecule FRET time
trajectories with or without absolute distance information.


  In special coordinates (codon position--specific nucleotide frequencies)
bacterial genomes form two straight lines in 9-dimensional space: one line for
eubacterial genomes, another for archaeal genomes. All the 348 distinct
bacterial genomes available in Genbank in April 2007, belong to these lines
with high accuracy. The main challenge now is to explain the observed high
accuracy. The new phenomenon of complementary symmetry for codon
position--specific nucleotide frequencies is observed. The results of analysis
of several codon usage models are presented. We demonstrate that the
mean--field approximation, which is also known as context--free, or complete
independence model, or Segre variety, can serve as a reasonable approximation
to the real codon usage. The first two principal components of codon usage
correlate strongly with genomic G+C content and the optimal growth temperature
respectively. The variation of codon usage along the third component is related
to the curvature of the mean-field approximation. First three eigenvalues in
codon usage PCA explain 59.1%, 7.8% and 4.7% of variation. The eubacterial and
archaeal genomes codon usage is clearly distributed along two third order
curves with genomic G+C content as a parameter.


  An ensemble of directed macromolecules on a lattice is considered, where the
constituting molecules are chosen as a random sequence of N different types.
The same type of molecules experiences a hard-core (exclusion) interaction. We
study the robustness of the macromolecules with respect to breaking and
substituting individual molecules, using a 1/N expansion. The properties depend
strongly on the density of macromolecules. In particular, the macromolecules
are robust against breaking and substituting at high densities.


  The contribution of a specific hydrogen bond in apoflavodoxin to protein
stability is investigated by combining theory, experiment and simulation.
Although hydrogen bonds are major determinants of protein structure and
function, their contribution to protein stability is still unclear and widely
debated. The best method so far devised to estimate the contribution of
side-chain interactions to protein stability is double-mutant-cycle analysis,
but the interaction energies so derived are not identical to incremental
binding energies (the energies quantifying net contributions of two interacting
groups to protein stability). Here we introduce double-deletion analysis of
isolated residue pairs as a means to precisely quantify incremental binding.
The method is exemplified by studying a surface-exposed hydrogen bond in a
model protein (Asp96/Asn128 in apoflavodoxin). Combined substitution of these
residues by alanines slightly destabilizes the protein, due to a decrease in
hydrophobic surface burial. Subtraction of this effect, however, clearly
indicates that the hydrogen-bonded groups in fact destabilize the native
conformation. In addition, Molecular Dynamics simulations and classic
double-mutant-cycle analysis explain quantitatively that, due to frustration,
the hydrogen bond must form in the native structure because, when the two
groups get approximated upon folding their binding becomes favorable. We would
like to remark two facts: that this is the first time the contribution of a
specific hydrogen bond to protein stability has been measured from experiment,
and that more hydrogen bonds need to be analyzed in order to draw general
conclusions on protein hydrogen bonds energetics. To that end, the double
deletion method should be of help.


  Determining the middle of the bacteria cell and the proper placement of the
septum is essential to the division of the bacterial cell. In E. coli, this
process depends on the proteins MinC, MinD, and MinE. Here, the Lattice
Boltzmann method (LBM) is used to study the dynamics of the oscillations of the
min proteins from pole to pole. This determines the midcell division plane at
the cellular level. The LBM is applied to the set of the deterministic reaction
diffusion equations proposed by Howard et. al. [1] to describe the dynamics of
the Min proteins. The LBM results are in good agreement with those of Howard et
al, and agree qualitatively with the experimental results. Our good results
indicate that the LBM can be an alternative computational tool for simulating
problems dealing with complex biological system which are described by the
reaction-diffusion equations.


  A new n-dimensional vector space of the DNA sequences on the Galois field of
the 64 codons (GF(64)) is proposed. In this vector space gene mutations can be
considered linear transformations or translations of the wild type gene. In
particular, the set of translations that preserve the chemical type of the
third base position in the codon is a subgroup which describes the most
frequent mutations observed in mutational variants of four genes: human
phenylalanine hydroxylase (PAH), human beta globin (HBG), HIV-1 Protease (HIVP)
and HIV-1 Reverse transcriptase (HIVRT). Furthermore, an inner pseudo-product
defined between codons tends to have a positive value when the codons code to
similar amino acids and a negative value when the codons code to amino acids
with extreme hydrophobic properties. Consequently, it is found that the inner
pseudo-product between the wild type and the mutant codons tends to have a
positive value in the mutational variants of the genes: PAH, HBG, HIVP, HIVRT.


  A great part of the effort in the study of coarse grained models of
transcription networks is directed to the analysis of their dynamical features.
In this letter, we consider the \emph{equilibrium} properties of such systems,
showing that the logic backbone underlying all dynamic descriptions has the
structure of a computational optimization problem. It involves variables, which
correspond to gene expression levels, and constraints, which describe the
effect of \emph{cis-}regulatory signal integration functions. In the simple
paradigmatic case of Boolean variables and signal integration functions, we
derive and discuss phase diagrams. Notably, the model exhibits a connectivity
transition between a regime of simple, but uncertain, gene control, to a regime
of complex combinatorial control.


  The vibrational dynamics of a DNA molecule with counterions neutralizing the
charged phosphate groups have been studied. With the help of elaborated model
the conformational vibrations of the DNA double helix with alkaline metal ions
have been described both qualitatively and quantitatively. For the complexes of
DNA with counterions Li+, Na+, K+, Rb+ and Cs+ the normal modes have been
found, and a mode characterized by the most notable ion displacements with
respect to the DNA backbone has been determined. The frequency of counterion
vibrations has been established to decrease as the ion mass increases. The
results of theoretical calculation have been showed to be in good agreement
with the experimental data of Raman spectroscopy.


  In E. coli the determination of the middle of the cell and the proper
placement of the septum is essential to the division of the cell. This step
depends on the proteins MinC, MinD, and MinE. Exposure to a constant external
field e.g., an electric field or magnetic field may cause the bacteria cell
division mechanism to change resulting in an abnormal cytokinesis. To have
insight into the effects of an external field on this process, we model the
process using a set of the deterministic reaction diffusion equations, which
incorporate the influence of an external field, min protein reactions, and
diffusion of all species. Using the numerical method, we have found some
changes in the dynamics of the oscillations of the min proteins from pole to
pole when compared that of without the external field. The results show some
interesting effects, which are qualitatively in good agreement with some
experimental results.


  Models relating to the Species-Area curve are usually defined at the species
level, and concerned only with ecological timescales. We examine an
individual-based model of co-evolution on a spatial lattice based on the
Tangled Nature model, and show that reproduction, mutation and dispersion by
diffusion in an interacting system produces power-law Species-Area Relations as
observed in ecological measurements at medium scales. We find that
co-evolutionary habitats form, allowing high diversity levels in a spatially
homogenous system, and these are maintained for exponentially increasing time
when increasing system size.


  To understand the mechanism of TATA-box conformational transformations we
model structure mobility and find the types of conformational excitations of
DNA macromolecule in heteronomous conformation. We have constructed the
two-component model for describing DNA conformational transformation with
simultaneous transitions in the furanos rings of the monomer link. Internal
component describes the change of the base pair position in the double helix.
External component describes the displacement of mass center of the monomer
link. Nonlinearity of the system is accounted with a form of potential energy
describing C3'-C2' and C2'-C3' sugars transitions in monomer link, and
interrelation between monomer conformational transition and macromolecule
deformation. The comparison of our results with experimental data allows to
confirm that the localized conformational excitations may realise in DNA
TATA-box. These excitations cause the deformation of the macromolecule
fragment.


  The motility of certain gram-negative bacteria is mediated by retraction of
type IV pili surface filaments, which are essential for infectivity. The
retraction is powered by a strong molecular motor protein, PilT, producing very
high forces that can exceed 150 pN. The molecular details of the motor
mechanism are still largely unknown, while other features have been identified,
such as the ring-shaped protein structure of the PilT motor. The surprisingly
high forces generated by the PilT system motivate a model investigation of the
generation of large forces in molecular motors. We propose a simple model,
involving a small ensemble of motor subunits interacting through the
deformations on a circular backbone with finite stiffness. The model describes
the motor subunits in terms of diffusing particles in an asymmetric,
time-dependent binding potential (flashing ratchet potential), roughly
corresponding to the ATP hydrolysis cycle. We compute force-velocity relations
in a subset of the parameter space and explore how the maximum force (stall
force) is determined by stiffness, binding strength, ensemble size, and degree
of asymmetry. We identify two qualitatively different regimes of operation
depending on the relation between ensemble size and asymmetry. In the
transition between these two regimes, the stall force depends nonlinearly on
the number of motor subunits. Compared to its constituents without
interactions, we find higher efficiency and qualitatively different
force-velocity relations. The model captures several of the qualitative
features obtained in experiments on pilus retraction forces, such as roughly
constant velocity at low applied forces and insensitivity in the stall force to
changes in the ATP concentration.


  Emerging evidence suggests that the introns and intergenic sequences of the
genomes of higher eukaryotes (the ``junk'' DNA) codes for a vast, RNA-based,
genetic regulatory network. It is believed that this network is responsible for
the variety and complexity of terrestrial life. We conjecture that this
regulatory network is more properly viewed as an RNA ``community'', composed of
a rich and largely unexplored biochemical web of RNA interactions. Viewed as an
RNA-community, we hypothesize that the RNA regulatory network of higher
eukaryotes can re-wire itself, and employ various and evolvable mutational
strategies in response to external pressures. Thus, we argue that much
evolutionary change is due to intracellular, RNA-mediated learning processes.
Successful strategies and pathways are then recorded (hard-wired) into the DNA
genome via reverse transcriptase. We present evidence which is consistent with
this viewpoint, and make specific predictions which could be used to test the
utility of our framework. If essentially correct, the RNA-community view of
eukaryotic cells could reconcile measured point mutation and gene duplication
rates with actual rates of evolutionary change. Futhermore, the RNA-community
view of eukaryotic cells suggests that agent-based modeling techniques, used in
mathematical economics, game theory, and neuroscience, will likely be as useful
in understanding the functioning of eukaryotic cells as the pathway-based
approaches of systems biology. We conclude this paper by arguing that a
sufficient amount of biological knowledge has been accumulated to initiate a
systematic program of experimental and computational studies of the origins and
macroevolution of terrestrial life.


  Because of the base pairing rules in DNA, some mutations experienced by a
portion of DNA during its evolution result in the same substitution, as we can
only observe differences in coupled nucleotides. Then, in the absence of a bias
between the two DNA strands, a model with at most 6 different parameters
instead of 12 is sufficient to study the evolutionary relationship between
homologous sequences derived from a common ancestor. On the other hand the same
symmetry reduces the number of independent observations which can be made. Such
a reduction can in some cases invalidate the calculation of the parameters. A
compromise between biologically acceptable hypotheses and tractability is
introduced and a five parameter reversible no-strand-bias condition (RNSB) is
presented. The identifiability of the parameters under this model is shown by
examples.


  We investigate the competition dynamics of two microbial or viral strains
that live in an environment that switches periodically between two states. One
of the strains is adapted to the long-term environment, but pays a short-term
cost, while the other is adapted to the short-term environment and pays a cost
in the long term. We explore the tradeoff between these alternative strategies
in extensive numerical simulations, and present a simple analytic model that
can predict the outcome of these competitions as a function of the mutation
rate and the time scale of the environmental changes. Our model is relevant for
arboviruses, which alternate between different host species on a regular basis.


  Much has been debated about the benefit of sexual over asexual reproduction
in terms of evolutionary fitness. Here we focus on the advantage that may be
brought about by the process of mating, where the choosing of mates contributes
to the increase in fitness in a constructive way. We carry out computer
simulations of such mating systems and investigate, on one hand, how mate
phenotypes contribute to offspring fitness, and, on the other hand, how
selection affects mate phenotypes. We discuss how helpful such a mechanism may
be in determining trajectories on rugged energy landscapes leading to global
optimum.


  A simple approach is proposed to investigate the protein structure. Using a
low complexity model, a simple pairwise interaction and the concept of global
optimization, we are able to calculate ground states of proteins, which are in
agreement with experimental data. All possible model structures of small
proteins are available below a certain energy threshold. The exact lowenergy
landscapes for the trp cage protein (1L2Y) is presented showing the
connectivity of all states and energy barriers.


  Mechanical unfolding of several domains of calmodulin and titin is studied
using a Go-like model with a realistic contact map and Lennard-Jones contact
interactions. It is shown that this simple model captures the experimentally
observed difference between the two proteins: titin is a spring that is tough
and strong whereas calmodulin acts like a weak spring with featureless
force-displacement curves. The difference is related to the dominance of the
alpha secondary structures in the native structure of calmodulin. The tandem
arrangements of calmodulin unwind simultaneously in each domain whereas the
domains in titin unravel in a serial fashion. The sequences of contact events
during unraveling are correlated with the contact order, i.e. with the
separation between contact making amino acids along the backbone in the native
state. Temperature is found to affect stretching in a profound way.


  The biological distinction between the base positions in the codon, the
chemical types of bases (purine and pyrimidine) and their hydrogen bond number
have been the most relevant codon properties used in the genetic code analysis.
Now, these properties have allowed us to build a Genetic Code ring isomorphic
to the ring (Z64, +,*) of the integer module 64. On the Z64-algebra of the set
of 64^N codon sequences of length N, gene mutations are described by means of
endomorphisms F: (Z64)^N->(Z64)^N. Endomorphisms and automorphisms helped us
describe the gene mutation pathways. For instance, 77.7% mutations in 749 HIV
protease gene sequences correspond to unique diagonal endomorphisms of the wild
type strain HXB2. In particular, most of the reported mutations that confer
drug resistance to the HIV protease gene correspond to diagonal automorphisms
of the wild type. What is more, in the human beta-globin gene a similar
situation appears where most of the single codon mutations correspond to
automorphisms. Hence, in the analyses of molecular evolution process on the DNA
sequence set of length N, the Z64-algebra will help us explain the quantitative
relationships between genes.


  The algebraic structures of the genetic code are most important to obtain
additional information about the semantic code and its applications. In this
paper we define two dual Boolean codon lattices of the genetic code using
hydrogen bond numbers and the chemical types of bases: purines and pyrimidines.
The Boolean lattices reflect the role of hydrophobicity in the distribution of
codon assignments to each amino acid. Particularly, the symmetric images of
codons with adenine as second base coding to hydrophilic amino acids are always
codons with uracil as second base coding to hydrophobic amino acids as they
represented in the Hasse diagrams. The Hamming distance between two codons in
the Hasse diagram reflects the different hydrophobicities between their
respective coded amino acids. Our experiments have demonstrated a small Hamming
distance to the wild type HXB2 of almost all the drug-resistant reported
mutations in HIV protease gene. The human beta-globin mutant genes have also
exhibited similar results. Our research suggests that the Hamming distance
between two genes in the molecular evolution process have a minimal value.


  Despite their apparent simplicity, random Boolean networks display a rich
variety of dynamical behaviors. Much work has been focused on the properties
and abundance of attractors. We here derive an expression for the number of
attractors in the special case of one input per node. Approximating some other
non-chaotic networks to be of this class, we apply the analytic results to
them. For this approximation, we observe a strikingly good agreement on the
numbers of attractors of various lengths. Furthermore, we find that for long
cycle lengths, there are some properties that seem strange in comparison to
real dynamical systems. However, those properties can be interesting from the
viewpoint of constraint satisfaction problems.


  Shannon information (SI) and its special case, divergence, are defined for a
DNA sequence in terms of probabilities of chemical words in the sequence and
are computed for a set of complete genomes highly diverse in length and
composition. We find the following: SI (but not divergence) is inversely
proportional to sequence length for a random sequence but is length-independent
for genomes; the genomic SI is always greater and, for shorter words and longer
sequences, hundreds to thousands times greater than the SI in a random sequence
whose length and composition match those of the genome; genomic SIs appear to
have word-length dependent universal values. The universality is inferred to be
an evolution footprint of a universal mode for genome growth.


  We discuss stochastic dynamics of finite populations of individuals playing
games. We review recent results concerning the dependence of the long-run
behavior of such systems on the number of players and the noise level. In the
case of two-player games with two symmetric Nash equilibria, when the number of
players increases, the population undergoes multiple transitions between its
equilibria.


  The effective DNA-DNA interaction force is calculated by computer simulations
with explicit tetravalent counterions and monovalent salt. For overcharged DNA
molecules, the interaction force shows a double-minimum structure. The
positions and depths of these minima are regulated by the counterion density in
the bulk. Using two-dimensional lattice sum and free energy perturbation
theories, the coexisting phases for DNA bundles are calculated. A
DNA-condensation and redissolution transition and a stable mesocrystal with an
intermediate lattice constant for high counterion concentration are obtained.


  Animal behavior and evolution can often be described by game-theoretic
models. Although in many situations, the number of players is very large, their
strategic interactions are usually decomposed into a sum of two-player games.
Only recently evolutionarily stable strategies were defined for multi-player
games and their properties analyzed (Broom et al., 1997). Here we study the
long-run behavior of stochastic dynamics of populations of randomly matched
individuals playing symmetric three-player games. We analyze stochastic
stability of equilibria in games with multiple evolutionarily stable
strategies. We also show that in some games, a population may not evolve in the
long run to an evolutionarily stable equilibrium.


  A-To-I RNA editing is common to all eukaryotes, associated with various
neurological functions. Recently, A-to-I editing was found to occur abundantly
in the human transcriptome. Here we show that the frequency of A-to-I editing
in humans is at least an order of magnitude higher as that of mouse, rat,
chicken or fly. The extraordinary frequency of RNA editing in human is
explained by the dominance of the primate-specific Alu element in the human
transcriptome, which increases the number of double-stranded RNA substrates.


  Exact law of mortality dynamics in changing populations and environment is
derived. The law is universal for all species, from single cell yeast to
humans. It includes no characteristics of animal- environment interactions
(metabolism etc) which are a must for life. Such law is unique for live systems
with their homeostatic self-adjustment to environment. Its universal dynamics
for all animals, with their drastically different biology, evolutionary
history, and complexity, is also unique for live systems-cf different
thermodynamics of e.g. liquids and glasses. The law which is valid for all
live, and only live, systems is a life specific law of nature. Mortality is an
instrument of natural selection and biological diversity. Its law which is
preserved in evolution of all species is a conservation law of mortality,
selection, evolution, biology. The law implies new kind of mortality and
adaptation which dominate in evolutionary unprecedented protected populations
and, in contrast to species specific natural selection, proceed via universal
stepwise rungs. The law demonstrates that intrinsic mortality and at least
certain aspects of aging are disposable evolutionary byproducts, and directed
genetic and/or biological changes may yield healthy and vital Methuselah
lifespan. This is consistent with experiments. Universality implies that single
cell yeast may provide a master key to the cellular mechanism of universal
mortality, aging, selection, evolution, and its regulation in all animals. One
may look for its manifestations in animal cells also, e.g., in their
replicative senescence.


  Despite considerable progress in genome- and proteome-based high-throughput
screening methods and rational drug design, the number of successful single
target drugs did not increase appreciably during the past decade. Network
models suggest that partial inhibition of a surprisingly small number of
targets can be more efficient than the complete inhibition of a single target.
This and the success stories of multi-target drugs and combinatorial therapies
led us to suggest that systematic drug design strategies should be directed
against multiple targets. We propose that the final effect of partial, but
multiple drug actions might often surpass that of complete drug action at a
single target. The future success of this novel drug design paradigm will
depend not only on a new generation of computer models to identify the correct
multiple hits and their multi-fitting, low-affinity drug candidates but also on
more efficient in vivo testing.


  By using a mixture model for the density distribution of the three pseudobond
angles formed by $C_\alpha$ atoms of four consecutive residues, the local
structural states are discretized as 17 conformational letters of a protein
structural alphabet. This coarse-graining procedure converts a 3D structure to
a 1D code sequence. A substitution matrix between these letters is constructed
based on the structural alignments of the FSSP database.


  In this report a systematic approach is used to determine the approximate
genetic network and robust dependencies underlying differentiation. The data
considered is in the form of a binary matrix and represent the expression of
the nine genes across the ninety-nine colonies. The report is divided into two
parts: the first part identifies significant pair-wise dependencies from the
given binary matrix using linear correlation and mutual information. A new
method is proposed to determine statistically significant dependencies
estimated using the mutual information measure. In the second, a Bayesian
approach is used to obtain an approximate description (equivalence class) of
network structures. The robustness of linear correlation, mutual information
and the equivalence class of networks is investigated with perturbation and
decreasing colony number. Perturbation of the data was achieved by generating
bootstrap realizations. The results are refined with biological knowledge. It
was found that certain dependencies in the network are immune to perturbation
and decreasing colony number and may represent robust features, inherent in the
differentiation program of osteoblast progenitor cells. The methods to be
discussed are generic in nature and not restricted to the experimental paradigm
addressed in this study.


  An overview of theories related to vibrational energy relaxation (VER) in
proteins is presented. VER of a selected mode in cytochrome c is studied using
two theoretical approaches. One is the equilibrium simulation approach with
quantum correction factors, and the other is the reduced model approach which
describes the protein as an ensemble of normal modes interacting through
nonlinear coupling elements. Both methods result in estimates of the VER time
(sub ps) for a CD stretching mode in the protein at room temperature. The
theoretical predictions are in accord with the experimental data of Romesberg's
group. A perspective on future directions for the detailed study of time scales
and mechanisms for VER in proteins is presented.


  We study a master equation system modelling a population dynamics problem in
a lattice. The problem is the calculation of the minimum size of a refuge that
can protect a population from hostile external conditions, the so called
critical patch size problem. We analize both cases in which the particles are
considered fermions and bosons and show using exact analitical methods that,
while the Fermi-Dirac statistics leads to certain extinction for any refuge
size, the Bose-Eistein statistics allows survival even for the minimal refuge.


  The functions of most genetic circuits require sufficient degrees of
cooperativity in the circuit components. While mechanisms of cooperativity have
been studied most extensively in the context of transcriptional initiation
control, cooperativity from other processes involved in the operation of the
circuits can also play important roles. In this study, we examine a simple
kinetic source of cooperativity stemming from the nonlinear degradation of
multimeric proteins. Ample experimental evidence suggests that protein subunits
can degrade less rapidly when associated in multimeric complexes, an effect we
refer to as cooperative stability. For dimeric transcription factors, this
effect leads to a concentration-dependence in the degradation rate because
monomers, which are predominant at low concentrations, will be more rapidly
degraded. Thus cooperative stability can effectively widen the accessible range
of protein levels in vivo. Through theoretical analysis of two exemplary
genetic circuits in bacteria, we show that such an increased range is important
for the robust operation of genetic circuits as well as their evolvability. Our
calculations demonstrate that a few-fold difference between the degradation
rate of monomers and dimers can already enhance the function of these circuits
substantially. These results suggest that cooperative stability needs to be
considered explicitly and characterized quantitatively in any systematic
experimental or theoretical study of gene circuits.


  An non-autonomous system is proposed to model the seasonal pattern of dengue
fever. We found that an approximate threshold condition for infection
persistence describes all possible behavior of the system. As far as we know,
the kind of analysis here proposed is entirely new. No precise mathematical
theorems are demonstrated but we give enough numerical evidence to support the
conclusions.