title: Folding thermodynamics of peptides

abstract: 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.


title: Finite size effects on calorimetric cooperativity of two-state proteins

abstract: 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.


title: Off-lattice simulation of the solid phase DNA amplification

abstract: 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.


title: Transcriptional Regulation by the Numbers 1: Models

abstract: 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.


title: Transcriptional Regulation by the Numbers 2: Applications

abstract: 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.


title: Subtree power analysis finds optimal species for comparative genomics

abstract: 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.


title: Surfaces and orientations: much to fret about?

abstract: 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.


title: Robust formation of morphogen gradients

abstract: We discuss the formation of graded morphogen profiles in a cell layer by nonlinear transport phenomena, important for patterning developing organisms. We focus on a process termed transcytosis, where morphogen transport results from binding of ligands to receptors on the cell surface, incorporation into the cell and subsequent externalization. Starting from a microscopic model, we derive effective transport equations. We show that, in contrast to morphogen transport by extracellular diffusion, transcytosis leads to robust ligand profiles which are insensitive to the rate of ligand production.


title: The Mystery of Two Straight Lines in Bacterial Genome Statistics. Release 2007

abstract: 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.


title: Correlations in Systems of Complex Directed Macromolecules

abstract: 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.


title: A double-deletion method to quantifying incremental binding energies in proteins from experiment. Example of a destabilizing hydrogen bonding pair

abstract: 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.


title: A one-dimensional Lattice Boltzmann method for modeling the dynamic pole-to-pole oscillations of Min proteins for determining the position of the midcell division plane

abstract: 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.


title: A New DNA Sequences Vector Space on a Genetic Code Galois Field

abstract: 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.


title: The Logic Backbone of a Transcription Network

abstract: 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.


title: Multiple sequence alignment based on set covers

abstract: We introduce a new heuristic for the multiple alignment of a set of sequences. The heuristic is based on a set cover of the residue alphabet of the sequences, and also on the determination of a significant set of blocks comprising subsequences of the sequences to be aligned. These blocks are obtained with the aid of a new data structure, called a suffix-set tree, which is constructed from the input sequences with the guidance of the residue-alphabet set cover and generalizes the well-known suffix tree of the sequence set. We provide performance results on selected BAliBASE amino-acid sequences and compare them with those yielded by some prominent approaches.


title: Ion-Phosphate Mode in the DNA Low-Frequency Spectra

abstract: 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.


title: The dynamics of the min proteins of Escherichia coli under the constant external fields

abstract: 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.


title: The species-area relationship and evolution

abstract: 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.


title: Conformational transformations of DNA macromolecule in heteronomous conformation

abstract: 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.


title: Force generation in small ensembles of Brownian motors

abstract: 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.


title: An RNA-centered view of eukaryotic cells

abstract: 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.


title: Forcing reversibility in the no strand-bias substitution model allows for the theoretical and practical identifiability of its 5 parameters from pairwise DNA sequence comparisons

abstract: 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.


title: Tradeoff between short-term and long-term adaptation in a changing environment

abstract: 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.


title: Love before Sex

abstract: 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.


title: Global optimization of proteins using a dynamical lattice model: Ground states and energy landscapes

abstract: 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.


title: Mechanical Stretching of Proteins: Calmodulin and Titin

abstract: 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.


title: Gene Algebra from a Genetic Code Algebraic Structure

abstract: 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.


title: The Genetic Code Boolean Lattice

abstract: 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.


title: Counting attractors in synchronously updated random Boolean networks

abstract: 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.


title: Computational Models of Adult Neurogenesis

abstract: Experimental results in recent years have shown that adult neurogenesis is a significant phenomenon in the mammalian brain. Little is known, however, about the functional role played by the generation and destruction of neurons in the context of and adult brain. Here we propose two models where new projection neurons are incorporated. We show that in both models, using incorporation and removal of neurons as a computational tool, it is possible to achieve a higher computational efficiency that in purely static, synapse-learning driven networks. We also discuss the implication for understanding the role of adult neurogenesis in specific brain areas.


title: Divergence and Shannon information in genomes

abstract: 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.


title: Equilibrium transitions in finite populations of players

abstract: 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.


title: DNA condensation and redissolution: Interaction between overcharged DNA molecules

abstract: 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.


title: Energetic model of tumor growth

abstract: A macroscopic model of the tumor Gompertzian growth is proposed. This approach is based on the energetic balance among the different cell activities, described by methods of statistical mechanics and related to the growth inhibitor factors. The model is successfully applied to the multicellular tumor spheroid data.


title: Stochastic stability in three-player games

abstract: 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.


title: The complex channel networks of bone structure

abstract: Bone structure in mammals involves a complex network of channels (Havers and Volkmann channels) required to nourish the bone marrow cells. This work describes how three-dimensional reconstructions of such systems can be obtained and represented in terms of complex networks. Three important findings are reported: (i) the fact that the channel branching density resembles a power law implies the existence of distribution hubs; (ii) the conditional node degree density indicates a clear tendency of connection between nodes with degrees 2 and 4; and (iii) the application of the recently introduced concept of hierarchical clustering coefficient allows the identification of typical scales of channel redistribution. A series of important biological insights is drawn and discussed

categories: q-bio.TO cond-mat.dis-nn physics.bio-ph q-bio.QM


title: Is abundant A-to-I RNA editing primate-specific?

abstract: 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.


title: Exact law of live nature

abstract: 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.


title: The efficiency of multi-target drugs: the network approach might help drug design

abstract: 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.


title: A protein structural alphabet and its substitution matrix CLESUM

abstract: 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.


title: Modeling Genetic Networks from Clonal Analysis

abstract: 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.


title: Vibrational energy relaxation in proteins

abstract: 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.


title: Particle Statistics and Population Dynamics

abstract: 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.


title: Chemotactic Collapse and Mesenchymal Morphogenesis

abstract: We study the effect of chemotactic signaling among mesenchymal cells. We show that the particular physiology of the mesenchymal cells allows one-dimensional collapse in contrast to the case of bacteria, and that the mesenchymal morphogenesis represents thus a more complex type of pattern formation than those found in bacterial colonies. We finally compare our theoretical predictions with recent in vitro experiments.


title: Elasticities and stabilities: lipid membranes vs cell membranes

abstract: A cell membrane can be simply regarded as composite material consisting of lipid bilayer, membrane cytoskeleton beneath lipid bilayer, and proteins embedded in lipid bilayer and linked with membrane cytoskeleton if one only concerns its mechanical properties. In this Chapter, above all, the authors give a brief introduction to some important work on mechanical properties of lipid bilayers following Helfrich's seminal work on spontaneous curvature energy of lipid bilayers. Next, the entropy of a polymer confined in a curved surface and the free energy of membrane cytoskeleton are obtained by scaling analysis. It is found that the free energy of cell membranes has the form of the in-plane strain energy plus Helfrich's curvature energy. The equations to describe equilibrium shapes and in-plane strains of cell membranes by osmotic pressures are obtained by taking the first order variation of the total free energy containing the elastic free energy, the surface tension energy and the term induced by osmotic pressure. The stability of spherical cell membrane is discussed and the critical pressure is found to be much larger than that of spherical lipid bilayer without membrane cytoskeleton. Lastly, the authors try to extend the present static mechanical model of cell membranes to the cell structure dynamics by proposing a group of coupling equations involving tensegrity architecture of cytoskeleton, fluid dynamics of cytoplasm and elasticities of cell membranes.


title: Nonlinear Protein Degradation and the Function of Genetic Circuits

abstract: 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.


title: An approximate threshold condition for non-autonomous system: an application to a vector-borne infection

abstract: 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.


title: Binding properties and evolution of homodimers in protein-protein interaction networks

abstract: We demonstrate that Protein-Protein Interaction (PPI) networks in several eucaryotic organisms contain significantly more self-interacting proteins than expected if such homodimers randomly appeared in the course of the evolution. We also show that on average homodimers have twice as many interaction partners than non-self-interacting proteins. More specifically the likelihood of a protein to physically interact with itself was found to be proportional to the total number of its binding partners. These properties of dimers are are in agreement with a phenomenological model in which individual proteins differ from each other by the degree of their ``stickiness'' or general propensity towards interaction with other proteins including oneself. A duplication of self-interacting proteins creates a pair of paralogous proteins interacting with each other. We show that such pairs occur more frequently than could be explained by pure chance alone. Similar to homodimers, proteins involved in heterodimers with their paralogs on average have twice as many interacting partners than the rest of the network. The likelihood of a pair of paralogous proteins to interact with each other was also shown to decrease with their sequence similarity. This all points to the conclusion that most of interactions between paralogs are inherited from ancestral homodimeric proteins, rather than established de novo after the duplication. We finally discuss possible implications of our empirical observations from functional and evolutionary standpoints.


title: Recoverable One-dimensional Encoding of Three-dimensional Protein Structures

abstract: Protein one-dimensional (1D) structures such as secondary structure and contact number provide intuitive pictures to understand how the native three-dimensional (3D) structure of a protein is encoded in the amino acid sequence. However, it has not been clear whether a given set of 1D structures contains sufficient information for recovering the underlying 3D structure. Here we show that the 3D structure of a protein can be recovered from a set of three types of 1D structures, namely, secondary structure, contact number and residue-wise contact order which is introduced here for the first time. Using simulated annealing molecular dynamics simulations, the structures satisfying the given native 1D structural restraints were sought for 16 proteins of various structural classes and of sizes ranging from 56 to 146 residues. By selecting the structures best satisfying the restraints, all the proteins showed a coordinate RMS deviation of less than 4\AA{} from the native structure, and for most of them, the deviation was even less than 2\AA{}. The present result opens a new possibility to protein structure prediction and our understanding of the sequence-structure relationship.


title: The posterior-Viterbi: a new decoding algorithm for hidden Markov models

abstract: Background: Hidden Markov models (HMM) are powerful machine learning tools successfully applied to problems of computational Molecular Biology. In a predictive task, the HMM is endowed with a decoding algorithm in order to assign the most probable state path, and in turn the class labeling, to an unknown sequence. The Viterbi and the posterior decoding algorithms are the most common. The former is very efficient when one path dominates, while the latter, even though does not guarantee to preserve the automaton grammar, is more effective when several concurring paths have similar probabilities. A third good alternative is 1-best, which was shown to perform equal or better than Viterbi. Results: In this paper we introduce the posterior-Viterbi (PV) a new decoding which combines the posterior and Viterbi algorithms. PV is a two step process: first the posterior probability of each state is computed and then the best posterior allowed path through the model is evaluated by a Viterbi algorithm. Conclusions: We show that PV decoding performs better than other algorithms first on toy models and then on the computational biological problem of the prediction of the topology of beta-barrel membrane proteins.


title: Pattern formation in a stochastic model of cancer growth

abstract: We investigate noise-induced pattern formation in a model of cancer growth based on Michaelis-Menten kinetics, subject to additive and multiplicative noises. We analyse stability properties of the system and discuss the role of diffusion and noises in the system's dynamics. We find that random dichotomous fluctuations in the immune response intensity along with Gaussian environmental noise lead to emergence of a spatial pattern of two phases, in which cancer cells, or, respectively, immune cells predominate.


title: Base pair interactions and hybridization isotherms of matched and mismatched oligonucleotide probes on microarrays

abstract: The lack of specificity in microarray experiments due to non-specific hybridization raises a serious problem for the analysis of microarray data because the residual chemical background intensity is not related to the expression degree of the gene of interest. We analyzed the concentration dependence of the signal intensity of perfect match (PM) and mismatch (MM) probes in terms using a microscopic binding model using a combination of mean hybridization isotherms and single base related affinity terms. The signal intensities of the PM and MM probes and their difference are assessed with regard to their sensitivity, specificity and resolution for gene expression measures. The presented theory implies the refinement of existing algorithms of probe level analysis to correct microarray data for non-specific background intensities.


title: A schematic age-structured compartment model of the impact of antiretroviral therapy on HIV incidence and prevalence

abstract: A simple deterministic model is proposed to represent the basic aspects concerning the effects of different antiretroviral treatment schedules on HIV incidence and prevalence of affected populations. The model mimics current treatment guidelines applied in Brazil. However, the model does not intend to fit the data with any acceptable degree of accuracy since uncertainties on the values of the parameters and on the precise effect of the treatment put some limits on the practical implications of our model from which only orders of magnitude and some qualitative effects can be deduced. So, this paper intends to provide a conceptual and mechanistic understanding of the possible long term effects of treatment on the dynamics of HIV transmission. According to the model, the effect of the treatment depends on the level of sexual activity of the subpopulations considered, being more pronounced on the subpopulations with the highest sexual activity levels. Also, inefficient treatment can be prejudicial depending on the level of sexual activity and on the capacity to provide adequate treatment coverages to the population affected.


title: A Solvable Sequence Evolution Model and Genomic Correlations

abstract: We study a minimal model for genome evolution whose elementary processes are single site mutation, duplication and deletion of sequence regions and insertion of random segments. These processes are found to generate long-range correlations in the composition of letters as long as the sequence length is growing, i.e., the combined rates of duplications and insertions are higher than the deletion rate. For constant sequence length, on the other hand, all initial correlations decay exponentially. These results are obtained analytically and by simulations. They are compared with the long-range correlations observed in genomic DNA, and the implications for genome evolution are discussed.


title: Information Capacity of Biological Macromoleculae Reloaded

abstract: Information capacity of a symbol sequence is a measure of the unexpectedness of a continuation of given string of symbols. Continuation of a string is determined through the maximum entropy of the reconstructed frequency dictionary; the capacity, in turn, is determined through the calculation of mutual entropy of a real frequency dictionary of a sequence with respect to the reconstructed one. The capacity does not depend on the length of strings in a dictionary. The capacity calculated for various genomes exhibits a multi-minima pattern reflecting an order observed within a sequence.


title: A New Simulated Annealing Algorithm for the Multiple Sequence Alignment Problem: The approach of Polymers in a Random Media

abstract: We proposed a probabilistic algorithm to solve the Multiple Sequence Alignment problem. The algorithm is a Simulated Annealing (SA) that exploits the representation of the Multiple Alignment between $D$ sequences as a directed polymer in $D$ dimensions. Within this representation we can easily track the evolution in the configuration space of the alignment through local moves of low computational cost. At variance with other probabilistic algorithms proposed to solve this problem, our approach allows for the creation and deletion of gaps without extra computational cost. The algorithm was tested aligning proteins from the kinases family. When D=3 the results are consistent with those obtained using a complete algorithm. For $D>3$ where the complete algorithm fails, we show that our algorithm still converges to reasonable alignments. Moreover, we study the space of solutions obtained and show that depending on the number of sequences aligned the solutions are organized in different ways, suggesting a possible source of errors for progressive algorithms.


title: Thermosynthesis as energy source for the RNA World: a new model for the origin of life

abstract: The thermosynthesis concept, biological free energy gain from thermal cycling, is combined with the concept of the RNA World. The resulting overall origin of life model gives new explanations for the emergence of the genetic code and the ribosome. The first protein named pF1 obtains the energy to support the RNA world by a thermal variation of F1 ATP synthase's binding change mechanism. This pF1 is the single translation product during the emergence of the genetic machinery. During thermal cycling pF1 condenses many substrates with broad specificity, yielding NTPs and randomly constituted protein and RNA libraries that contain (self)-replicating RNA. The smallness of pF1 permits the emergence of the genetic machinery by selection of RNA that increases the fraction of pF1s in the protein library: (1) a progenitor of rRNA that concatenates amino acids bound to (2) a chain of 'positional tRNAs' linked by mutual recognition, yielding a pF1 (or its main motif); this positional tRNA set gradually evolves to a set of regular tRNAs functioning according to the genetic code, with concomitant emergence of (3) an mRNA coding for pF1.


title: Single Stranded DNA Translocation Through a Fluctuating Nanopore

abstract: We investigate the translocation of a single stranded DNA (ssDNA) through a pore, which fluctuates between two conformations, by using coupled master equations (ME). The probability density function (PDF) of the first passage times (FPT) of the translocation process is calculated, displaying a triple, double or mono-peaked behavior, depending on the system parameters. An analytical expression for the mean first passage time (MFPT) of the translocation process is derived, and provides an extensive characterization of the translocation process.


title: Predicting Residue-wise Contact Orders of Native Protein Structure from Amino Acid Sequence

abstract: Residue-wise contact order (RWCO) is a new kind of one-dimensional protein structures which represents the extent of long-range contacts. We have recently shown that a set of three types of one-dimensional structures (secondary structure, contact number, and RWCO) contains sufficient information for reconstructing the three-dimensional structure of proteins. Currently, there exist prediction methods for secondary structure and contact number from amino acid sequence, but none exists for RWCO. Also, the properties of amino acids that affect RWCO is not clearly understood. Here, we present a linear regression-based method to predict RWCO from amino acid sequence, and analyze the regression parameters to identify the properties that correlates with the RWCO. The present method achieves the significant correlation of 0.59 between the native and predicted RWCOs on average. An unusual feature of the RWCO prediction is the remarkably large optimal half window size of 26 residues. The regression parameters for the central and near-central residues of the local sequence segment highly correlate with those of the contact number prediction, and hence with hydrophobicity.


title: Chemical organization theory: towards a theory of constructive dynamical systems

abstract: Complex dynamical networks consisting of many components that interact and produce each other are difficult to understand, especially, when new components may appear. In this paper we outline a theory to deal with such systems. The theory consists of two parts. The first part introduces the concept of a chemical organization as a closed and mass-maintaining set of components. This concept allows to map a complex (reaction) network to the set of organizations, providing a new view on the system's structure. The second part connects dynamics with the set of organizations, which allows to map a movement of the system in state space to a movement in the set of organizations.


title: Active percolation analysis of pyramidal neurons of somatosensory cortex: A comparison of wildtype and p21H-ras^Val12 transgenic mice

abstract: This article describes the investigation of morphological variations among two set of neuronal cells, namely a control group of wild type rat cells and a group of cells of a trangenic line. Special attention is given to sigular points in the neuronal structure, namely the branching points and extremities of the dendritic processes. The characterization of the spatial distribution of such points is obtained by using a recently reported morphological technique based on forced percolation and window-size compensation, which is particularly suited to the analysis of scattered points presenting several coexisting densities. Different dispersions were identified in our statistical analysis, suggesting that the transgenic line of neurons is characterized by a more pronounced morphological variation. A classification scheme based on a canonical discriminant function was also considered in order to identify the morphological differences.


title: Identification and Measurement of Neighbor Dependent Nucleotide Substitution Processes

abstract: The presence of neighbor dependencies generated a specific pattern of dinucleotide frequencies in all organisms. Especially, the CpG-methylation-deamination process is the predominant substitution process in vertebrates and needs to be incorporated into a more realistic model for nucleotide substitutions. Based on a general framework of nucleotide substitutions we develop a method that is able to identify the most relevant neighbor dependent substitution processes, measure their strength, and judge their importance to be included into the modeling. Starting from a model for neighbor independent nucleotide substitution we successively add neighbor dependent substitution processes in the order of their ability to increase the likelihood of the model describing given data. The analysis of neighbor dependent nucleotide substitutions in human, zebrafish and fruit fly is presented. A web server to perform the presented analysis is publicly available.


title: Informational Way to Protein Alphabet: Entropic Classification of Amino Acids

abstract: What are proteins made from, as the working parts of the living cells protein machines? To answer this question, we need a technology to disassemble proteins onto elementary func-tional details and to prepare lumped description of such details. This lumped description might have a multiple material realization (in amino acids). Our hypothesis is that informational approach to this problem is possible. We propose a way of hierarchical classification that makes the primary structure of protein maximally non-random. The first steps of the suggested research program are realized: the method and the analysis of optimal informational protein binary alphabet. The general method is used to answer several specific questions, for example: (i) Is there a syntactic difference between Globular and Membrane proteins? (ii) Are proteins random sequences of amino acids (a long discussion)? For these questions, the answers are as follows: (i) There exists significant syntactic difference between Globular and Membrane proteins, and this difference is described; (ii) Amino acid sequences in proteins are definitely not random.


title: Substantial regional variation in substitution rates in the human genome: importance of GC content, gene density and telomere-specific effects

abstract: This study presents the first global, 1 Mbp level analysis of patterns of nucleotide substitutions along the human lineage. The study is based on the analysis of a large amount of repetitive elements deposited into the human genome since the mammalian radiation, yielding a number of results that would have been difficult to obtain using the more conventional comparative method of analysis. This analysis revealed substantial and consistent variability of rates of substitution, with the variability ranging up to 2-fold among different regions. The rates of substitutions of C or G nucleotides with A or T nucleotides vary much more sharply than the reverse rates suggesting that much of that variation is due to differences in mutation rates rather than in the probabilities of fixation of C/G vs. A/T nucleotides across the genome. For all types of substitution we observe substantially more hotspots than coldspots, with hotspots showing substantial clustering over tens of Mbp's. Our analysis revealed that GC-content of surrounding sequences is the best predictor of the rates of substitution. The pattern of substitution appears very different near telomeres compared to the rest of the genome and cannot be explained by the genome-wide correlations of the substitution rates with GC content or exon density. The telomere pattern of substitution is consistent with natural selection or biased gene conversion acting to increase the GC-content of the sequences that are within 10-15 Mbp away from the telomere.


title: Spike timing precision and neural error correction: local behavior

abstract: The effects of spike timing precision and dynamical behavior on error correction in spiking neurons were investigated. Stationary discharges -- phase locked, quasiperiodic, or chaotic -- were induced in a simulated neuron by presenting pacemaker presynaptic spike trains across a model of a prototypical inhibitory synapse. Reduced timing precision was modeled by jittering presynaptic spike times. Aftereffects of errors -- in this communication, missed presynaptic spikes -- were determined by comparing postsynaptic spike times between simulations identical except for the presence or absence of errors. Results show that the effects of an error vary greatly depending on the ongoing dynamical behavior. In the case of phase lockings, a high degree of presynaptic spike timing precision can provide significantly faster error recovery. For non-locked behaviors, isolated missed spikes can have little or no discernible aftereffects (or even serve to paradoxically reduce uncertainty in postsynaptic spike timing), regardless of presynaptic imprecision. This suggests two possible categories of error correction: high-precision locking with rapid recovery and low-precision non-locked with error immunity.


title: Stochastic models in population biology and their deterministic analogs

abstract: In this paper we introduce a class of stochastic population models based on "patch dynamics". The size of the patch may be varied, and this allows one to quantify the departures of these stochastic models from various mean field theories, which are generally valid as the patch size becomes very large. These models may be used to formulate a broad range of biological processes in both spatial and non-spatial contexts. Here, we concentrate on two-species competition. We present both a mathematical analysis of the patch model, in which we derive the precise form of the competition mean field equations (and their first order corrections in the non-spatial case), and simulation results. These mean field equations differ, in some important ways, from those which are normally written down on phenomenological grounds. Our general conclusion is that mean field theory is more robust for spatial models than for a single isolated patch. This is due to the dilution of stochastic effects in a spatial setting resulting from repeated rescue events mediated by inter-patch diffusion. However, discrete effects due to modest patch sizes lead to striking deviations from mean field theory even in a spatial setting.


title: Predator-prey cycles from resonant amplification of demographic stochasticity

abstract: In this paper we present the simplest individual level model of predator-prey dynamics and show, via direct calculation, that it exhibits cycling behavior. The deterministic analogue of our model, recovered when the number of individuals is infinitely large, is the Volterra system (with density-dependent prey reproduction) which is well-known to fail to predict cycles. This difference in behavior can be traced to a resonant amplification of demographic fluctuations which disappears only when the number of individuals is strictly infinite. Our results indicate that additional biological mechanisms, such as predator satiation, may not be necessary to explain observed predator-prey cycles in real (finite) populations.


title: Competing Polymerization of Actin Skeleton explains Relation between Network Polarity and Cell Movements

abstract: Based on experimental observations it is known that various biological cells exhibit a persistent random walk during migration on flat substrates. The persistent random walk is characterized by `stop-and-go' movements : unidirectional motions over distances of the order of several cell diameter are separated by localized short time erratic movements. Using computer simulations the reasons for this phenomena had been unveiled and shown to be attributed to two antagonistic nucleation processes during the polymerization of the cell's actin cytoskeleton : the (ordinary) spontaneous nucleation and the dendritic nucleation processes. Whereas spontaneous nucleations generate actin filaments growing in different directions and hence create motions in random directions, dendritic nucleations provide a unidirectional growth. Since dendritic growth exhibits stochastic fluctuations, spontaneous nucleation may eventually compete or even dominate, which results in a reorientation of filament growth and hence a new direction of cell motion. The event of reorientation takes place at instants of vanishing polarity of the actin skeleton.


title: Dynamics of competing species in a model of adaptive radiations and macroevolution

abstract: We present a simple model of adaptive radiations in evolution based on species competition. Competition is found to promote species divergence and branching, and to dampen the net species production. In the model simulations, high taxonomic diversification and branching take place during the beginning of the radiation. The results show striking similarities with empirical data and highlight the mechanism of competition as an important driving factor for accelerated evolutionary transformation.


title: Synaptic Plasticity with Discrete state synapses

abstract: Experimental observations on synaptic plasticity at individual glutamatergic synapses from the CA3 Shaffer collateral pathway onto CA1 pyramidal cells in the hippocampus suggest that the transitions in synaptic strength occur among discrete levels at individual synapses (~\cite{Peter} and S. S.-H. Wang, unpublished data used with the authors' permission). This happens for both long term potentiation (LTP) and long term depression (LTD) induction protocols. O'Connor, Wittenberg, and Wang have argued that three states would account for their observations on individual synapses in the CA3-CA1 pathway. We develop a quantitative model of this three state system with transitions among the states determined by a competition between kinases and phosphatases shown by O'Connor et al. to be determinant of LTP and LTD, respectively. Specific predictions for various plasticity protocols are given by coupling this description of discrete synaptic AMPA conductance changes to a model of postsynaptic membrane potential and associated intracellular calcium fluxes to yield the transition rates among the states. We then present various LTP and LTD induction protocols to the model system and report the resulting whole cell changes in AMPA conductance. We also examine the effect of our discrete state synaptic plasticity model on the synchronization of realistic oscillating neurons. We show that one-to-one synchronization is enhanced by the plasticity we discuss here and the presynaptic and postsynaptic oscillations are in phase. Synaptic strength saturates naturally in this model and does not require artificial upper or lower cutoffs, in contrast to earlier models of plasticity.


title: Are scale-free regulatory networks larger than random ones?

abstract: Network of packages with regulatory interactions (dependences and conflicts) from Debian GNU/Linux operating system is compiled and used as analogy of a gene regulatory network. Using a trace-back algorithm we assembly networks from the potential pool of packages for both scale-free and exponential topology from real and a null model data, respectively. We calculate the maximum number of packages that can be functionally installed in the system (i.e., the active network size). We show that scale-free regulatory networks allow a larger active network size than random ones. Small genomes with scale-free regulatory topology could allow much more functionality than large genomes with an exponential one, with implications on its dynamics, robustness and evolution.


title: Evolutionary trajectories in rugged fitness landscapes

abstract: We consider the evolutionary trajectories traced out by an infinite population undergoing mutation-selection dynamics in static, uncorrelated random fitness landscapes. Starting from the population that consists of a single genotype, the most populated genotype \textit{jumps} from a local fitness maximum to another and eventually reaches the global maximum. We use a strong selection limit, which reduces the dynamics beyond the first time step to the competition between independent mutant subpopulations, to study the dynamics of this model and of a simpler one-dimensional model which ignores the geometry of the sequence space. We find that the fit genotypes that appear along a trajectory are a subset of suitably defined fitness \textit{records}, and exploit several results from the record theory for non-identically distributed random variables. The genotypes that contribute to the trajectory are those records that are not \textit{bypassed} by superior records arising further away from the initial population. Several conjectures concerning the statistics of bypassing are extracted from numerical simulations. In particular, for the one-dimensional model, we propose a simple relation between the bypassing probability and the dynamic exponent which describes the scaling of the typical evolution time with genome size. The latter can be determined exactly in terms of the extremal properties of the fitness distribution.


title: Uncorrelated two-state single molecule trajectories from reducible kinetic schemes

abstract: Trajectories of on-off events are the output of many single molecule experiments. Usually, one describes the underlying mechanism that generates the trajectory using a kinetic scheme, and by analyzing the trajectory aims at deducing this scheme. In a previous work [O. Flomenbom, J. Klafter, and A. Szabo, submitted (2004)], we showed that when successive events along a trajectory are uncorrelated, all the information in the trajectory is contained in two basic functions, which are the waiting time probability functions (PDFs) of the on state and of the off state. The kinetic schemes that lead to such uncorrelated trajectories were termed reducible. Here we discuss the reasons that lead to reducible schemes. In particular, the topology of reducible schemes is characterized and proven.


title: Molecular Phylogenetic Analyses and Real Life Data

abstract: In molecular phylogeny, relationships among organisms are reconstructed using DNA or protein sequences and are displayed as trees. A linear increase in the number of sequences results in an exponential increase of possible trees. Thus, inferring trees from molecular data was shown to be NP-hard. This causes problems, if large data sets are used. This review gives an introduction to molecular phylogenetic methods and to the problems biologists are facing in molecular phylogenetic analyses.


title: Heisenberg uncertainty principle for thermal response of the microtubules excited by ultra-short laser pulses

abstract: In this paper the heat signaling in microtubules (MT) is investigated. It is argued that for the description of the heat signaling phenomena in MT, the hyperbolic heat transport (HHT) equation must be used. It is shown that HHT is the Klein-Gordon (K-G) equation. The general solution for the K-G equation for MT is obtained. For the undistorted signal propagation in MT the Heisenberg uncertainty principle is formulated and discussed. Key words: Microtubules; Heat signaling; Klein-Gordon equation; Heisenberg principle.


title: Filament depolymerization by motor molecules

abstract: Motor proteins that specifically interact with the ends of cytoskeletal filaments can induce filament depolymerization. A phenomenological description of this process is presented. We show that under certain conditions motors dynamically accumulate at the filament ends. We compare simulations of two microscopic models to the phenomenological description. The depolymerization rate can exhibit maxima and dynamic instabilities as a function of the bulk motor density for processive depolymerization. We discuss our results in relation to experimental studies of Kin-13 family motor proteins.


title: Global divergence of microbial genome sequences mediated by propagating fronts

abstract: We model the competition between recombination and point mutation in microbial genomes, and present evidence for two distinct phases, one uniform, the other genetically diverse. Depending on the specifics of homologous recombination, we find that global sequence divergence can be mediated by fronts propagating along the genome, whose characteristic signature on genome structure is elucidated, and apparently observed in closely-related {\it Bacillus} strains. Front propagation provides an emergent, generic mechanism for microbial "speciation", and suggests a classification of microorganisms on the basis of their propensity to support propagating fronts.


title: Evolution of Variance in Offspring Number: the Effects of Population Size and Migration

abstract: It was shown by Gillespie (1974) that if two genotypes produce the same average number of offspring on but have a different variance associated within each generation, the genotype with a lower variance will have a higher effective fitness. Specifically, the effective fitness is W(effective)=w-var/N, where w is the mean fitness, var is the variance in offspring number, and N is the total population size. The model also predicts that if a strategy has a higher arithmetic mean fitness and a higher variance than the competitor, the outcome of selection will depend on the population size (with larger population sizes favoring the high variance, high mean genotype). This suggests that for metapopulations with large numbers of (relatively) small demes, a strategy with lower variance and lower mean may be favored if the migration rate is low while higher migration rates (consistent with a larger effective population size) favor the opposite strategy. Individual based simulation confirms that this is indeed the case for an island model of migration, though the effect of migration differs greatly depending on whether migration precedes or follows selection. It is noted in the appendix that while Gillespie 1974 does seem to be heuristically accurate, it is not clear that the definition of effective fitness follows from his derivation.


title: Power-law scaling in dimension-to-biomass relationship of fish schools

abstract: Motivated by the finding that there is some biological universality in the relationship between school geometry and school biomass of various pelagic fishes in various conditions, I here establish a scaling law for school dimensions: the school diameter increases as a power-law function of school biomass. The power-law exponent is extracted through the data collapse, and is close to 3/5. This value of the exponent implies that the mean packing density decreases as the school biomass increases, and the packing structure displays a mass-fractal dimension of 5/3. By exploiting an analogy between school geometry and polymer chain statistics, I examine the behavioral algorithm governing the swollen conformation of large-sized schools of pelagics, and I explain the value of the exponent.

categories: q-bio.PE cond-mat.other cond-mat.stat-mech q-bio.OT


title: A Novel Lie Algebra of the Genetic Code over the Galois Field of Four DNA Bases

abstract: By starting from the four DNA bases order in the Boolean lattice, a novel Lie Algebra of the genetic code is proposed. Here, the principal partitions of the genetic code table were obtained as equivalent classes of quotient subspaces of the genetic code vector space over the Galois field of the four DNA bases. The new algebraic structure shows strong connections among algebraic relationships, codon assignment and physicochemical properties of amino acids. Moreover, a distance function defined between codons in the Lie algebra was demonstrated to have a linear behavior respect to physical variables such as the mean of amino acids energy in proteins. It was also noticed that the distance between wild type and mutant codons approach smaller values in mutational variants of four genes, i.e, human phenylalanine hydroxylase, human beta-globin, HIV-1 protease and HIV-1 reverse transcriptase. These results strongly suggest that deterministic rules in genetic code origin must be involved.


title: Robustness and fragility of Boolean models for genetic regulatory networks

abstract: Interactions between genes and gene products give rise to complex circuits that enable cells to process information and respond to external signals. Theoretical studies often describe these interactions using continuous, stochastic, or logical approaches. We propose a new modeling framework for gene regulatory networks, that combines the intuitive appeal of a qualitative description of gene states with a high flexibility in incorporating stochasticity in the duration of cellular processes. We apply our methods to the regulatory network of the segment polarity genes, thus gaining novel insights into the development of gene expression patterns. For example, we show that very short synthesis and decay times can perturb the wild type pattern. On the other hand, separation of timescales between pre- and posttranslational processes and a minimal prepattern ensure convergence to the wild type expression pattern regardless of fluctuations.


title: Comment on: "Characterization of subthreshold voltage fluctuations in neuronal membranes" by M. Rudolph and A. Destexhe

abstract: In two recent papers, Rudolph and Destexhe (Neural Comp. {\bf 15}, 2577-2618, 2003; Neural Comp. in press, 2005) studied a leaky integrator model (i.e. an RC-circuit) driven by correlated (``colored'') Gaussian conductance noise and Gaussian current noise. In the first paper they derived an expression for the stationary probability density of the membrane voltage; in the second paper this expression was modified to cover a larger parameter regime. Here we show by standard analysis of solvable limit cases (white-noise limit of additive and multiplicative noise sources; only slow multiplicative noise; only additive noise) and by numerical simulations that their first result does not hold for the general colored-noise case and uncover the errors made in the derivation of a Fokker-Planck equation for the probability density. Furthermore, we demonstrate analytically (including an exact integral expression for the time-dependent mean value of the voltage) and by comparison to simulation results, that the extended expression for the probability density works much better but still does not solve exactly the full colored-noise problem. We also show that at stronger synaptic input the stationary mean value of the linear voltage model may diverge and give an exact condition relating the system parameters for which this takes place.


title: Functionality Encoded In Topology? Discovering Macroscopic Regulatory Modules from Large-Scale Protein-DNA Interaction Networks

abstract: The promise of discovering a functional blueprint of a cellular system from large-scale and high-throughput sequence and experimental data is predicated on the belief that the same top-down investigative approach that proved successful in other biological problems (e.g. DNA sequencing) will be as effective when it comes to inferring more complex intracellular processes. The results in this paper address this fundamental issue in the specific context of transcription regulatory networks. Although simple recurring regulatory motifs have been identified in the past, due to the size and complexity of the connectivity structure, the subdivision of such networks into larger, and possibly inter-connected, regulatory modules is still under investigation. Specifically, it is unclear whether functionally well-characterized transcriptional sub-networks can be identified by solely analyzing the connectivity structure of the overall network topology. In this paper, we show that transcriptional regulatory networks can be systematically partitioned into communities whose members are consistently functionally related. We applied the partitioning method to the transcriptional regulatory networks of the yeast Saccharomyces cerevisiae; the resulting communities of gene and transcriptional regulators can be associated to distinct functional units, such as amino acid metabolism, cell cycle regulation, protein biosynthesis and localization, DNA replication and maintenance, lipid catabolism, stress response and so on. Moreover, the observation of inter-community connectivity patterns provides a valuable tool for elucidating the inter-dependency between the discovered regulatory modules.


title: The Dynamic Evolution of the Power Exponent in a Universal Growth Model of Tumors

abstract: We have previously reported that a universal growth law, as proposed by West and collaborators for all living organisms, appears to be able to describe also the growth of tumors in vivo. In contrast to the assumption of a fixed power exponent p (assumed by West et al. to be equal to 3/4), we show in this paper the dynamic evolution of p from 2/3 to 1, using experimental data from the cancer literature and in analogy with results obtained by applying scaling laws to the study of fragmentation of solids. The dynamic behaviour of p is related to the evolution of the fractal topology of neoplastic vascular systems and might be applied for diagnostic purposes to mark the emergence of a functionally sufficient (or effective) neo-angiogenetic structure.


title: The Effects of EGF-Receptor Density on Multiscale Tumor Growth Patterns

abstract: We studied the effects of epidermal growth factor receptor (EGFR) density on tumor growth dynamics, both on the sub- and the multi-cellular level using our previously developed model. This algorithm simulates the growth of a brain tumor using a multi-scale two-dimensional agent-based approach with an integrated transforming growth factor alpha (TGFalpha) induced EGFR-gene-protein interaction network. The results confirm that increasing cell receptor density correlates with an acceleration of the tumor system's spatio-temporal expansion dynamics. This multicellular behavior cannot be explained solely on the basis of spatial sub-cellular dynamics, which remain qualitatively similar amongst the three glioma cell lines investigated here in silico. Rather, we find that cells with higher EGFR density show an early increase in the phenotypic switching activity between proliferative and migratory traits, linked to a higher level of initial auto-stimulation by the PLCgamma-mediated TGFalpha-EGFR autocrine network. This indicates a more active protein level interaction in these chemotactically acting tumor systems and supports the role of post-translational regulation for the implemented EGFR pathway. Implications of these results for experimental cancer research are discussed.


title: Autoassociative Memory Retrieval and Spontaneous Activity Bumps in Small-World Networks of Integrate-and-Fire Neurons

abstract: Qualitatively, some real networks in the brain could be characterized as 'small worlds', in the sense that the structure of their connections is intermediate between the extremes of an orderly geometric arrangement and of a geometry-independent random mesh. Small worlds can be defined more precisely in terms of their mean path length and clustering coefficient; but is such a precise description useful to better understand how the type of connectivity affects memory retrieval? We have simulated an autoassociative memory network of integrate-and-fire units, positioned on a ring, with the network connectivity varied parametrically between ordered and random. We find that the network retrieves when the connectivity is close to random, and displays the characteristic behavior of ordered nets (localized 'bumps' of activity) when the connectivity is close to ordered. Recent analytical work shows that these two behaviours can coexist in a network of simple threshold-linear units, leading to localized retrieval states. We find that they tend to be mutually exclusive behaviours, however, with our integrate-and-fire units. Moreover, the transition between the two occurs for values of the connectivity parameter which are not simply related to the notion of small worlds.


title: Telomere loss limits the rate of human epithelial tumor formation

abstract: Most human carcinomas exhibit telomere abnormalities early in the carcinogenesis process suggesting that crisis caused by telomere shortening may be a necessary event leading to human carcinomas. Epidemiological records of the age at which each patient in a population develops carcinoma are known as age-incidence data; these provide a quantitative measure of human tumor initiation and dynamics. If crisis brought on by telomere shortening is necessary for most human carcinomas, it may also be the rate limiting step. To test this, we compared a mathematical model in which telomere loss is the rate limiting step during carcinogenesis with age-incidence data compiled by the Surveillance, Epidemiology and End Results (SEER) program. We found that this model adequately explains the age-incidence data. The model also implies that two distinct paths exist for carcinoma to develop in prostate, breast, and ovary tissues. We conclude that a single step, crisis brought on by telomere shortening, limits the rate of formation of human carcinomas.


title: Cliques and duplication-divergence network growth

abstract: A population of complete subgraphs or cliques in a network evolving via duplication-divergence is considered. We find that a number of cliques of each size scales linearly with the size of the network. We also derive a clique population distribution that is in perfect agreement with both the simulation results and the clique statistic of the protein-protein binding network of the fruit fly. In addition, we show that such features as fat-tail degree distribution, various rates of average degree growth and non-averaging, revealed recently for only the particular case of a completely asymmetric divergence, are present in a general case of arbitrary divergence.


title: What can one learn from two-state single molecule trajectories?

abstract: A time trajectory of an observable that fluctuates between two values (say, on and off), stemming from some unknown multi-substate kinetic scheme, is the output of many single molecule experiments. Here we show that when all successive waiting times along the trajectory are uncorrelated the on and the off waiting time probability density functions (PDFs) contain all the information. By relating the lack of correlation in the trajectory to the topology of kinetic schemes, we can immediately specify those kinetic schemes that are equally consistent with experiment, which means that it is impossible to differentiate between them by any sophisticated analyses of the trajectory. Correlated trajectories, however, contain additional information about the underlying kinetic scheme, and we consider the strategy that one should use to extract it. An example is given on correlations in the activity of individual lipase molecules.


title: Some Properties of the Speciation Model for Food-Web Structure - Mechanisms for Degree Distributions and Intervality

abstract: We present a mathematical analysis of the speciation model for food-web structure, which had in previous work been shown to yield a good description of empirical data of food-web topology. The degree distributions of the network are derived. Properties of the speciation model are compared to those of other models that successfully describe empirical data. It is argued that the speciation model unifies the underlying ideas of previous theories. In particular, it offers a mechanistic explanation for the success of the niche model of Williams and Martinez and the frequent observation of intervality in empirical food webs.


title: Dissociation in a polymerization model of homochirality

abstract: A fully self-contained model of homochirality is presented that contains the effects of both polymerization and dissociation. The dissociation fragments are assumed to replenish the substrate from which new monomers can grow and undergo new polymerization. The mean length of isotactic polymers is found to grow slowly with the normalized total number of corresponding building blocks. Alternatively, if one assumes that the dissociation fragments themselves can polymerize further, then this corresponds to a strong source of short polymers, and an unrealistically short average length of only 3. By contrast, without dissociation, isotactic polymers becomes infinitely long.


title: Statistical analysis of simple repeats in the human genome

abstract: The human genome contains repetitive DNA at different level of sequence length, number and dispersion. Highly repetitive DNA is particularly rich in homo-- and di--nucleotide repeats, while middle repetitive DNA is rich of families of interspersed, mobile elements hundreds of base pairs (bp) long, among which the Alu families. A link between homo- and di-polymeric tracts and mobile elements has been recently highlighted. In particular, the mobility of Alu repeats, which form 10% of the human genome, has been correlated with the length of poly(A) tracts located at one end of the Alu. These tracts have a rigid and non-bendable structure and have an inhibitory effect on nucleosomes, which normally compact the DNA. We performed a statistical analysis of the genome-wide distribution of lengths and inter--tract separations of poly(X) and poly(XY) tracts in the human genome. Our study shows that in humans the length distributions of these sequences reflect the dynamics of their expansion and DNA replication. By means of general tools from linguistics, we show that the latter play the role of highly-significant content-bearing terms in the DNA text. Furthermore, we find that such tracts are positioned in a non-random fashion, with an apparent periodicity of 150 bases. This allows us to extend the link between repetitive, highly mobile elements such as Alus and low-complexity words in human DNA. More precisely, we show that Alus are sources of poly(X) tracts, which in turn affect in a subtle way the combination and diversification of gene expression and the fixation of multigene families.


title: Sympatric Speciation in a Simple Food Web

abstract: Observations of the evolution of species groups in nature, such as well recognized Galapagos finches, have motivated much theoretical research aimed at understanding the processes associated with such radiations. The Penna model is one such model and has been widely used to study aging. In this paper we use the basic Penna model to investigate the process of sympatric speciation in a simple food web model. Initially our web consists of a primary food source and a single herbivore species that feeds on this resource. Subsequently we introduce a predator that feeds on the herbivore. In both instances we directly manipulate the food source, that is, its size distribution, and monitore the changes in the populations structures. Sympatric speciation is obtained for the consumer species in both webs, and our results confirm that the speciation velocity depends on how far up,in the food chain, the focus population is feeding. Simulations are done with three different sexual imprinting-like mechanisms, in order to discuss adaptation by natural selection.


title: Determination of barrier heights and prefactors from protein folding rate data

abstract: We determine both barrier heights and prefactors for protein folding by applying constraints determined from experimental rate measurements to a Kramers theory for folding rate. The theoretical values are required to match the experimental values at two conditions of temperature and denaturant that induce the same stability. Several expressions for the prefactor in the Kramers rate equation are examined: a random energy approximation, a correlated energy approximation, and an approximation using a single Arrhenius activation energy. Barriers and prefactors are generally found to be large as a result of implementing this recipe, i.e. the folding landscape is cooperative and smooth. Interestingly, a prefactor with a single Arrhenius activation energy admits no formal solution.


title: Short-range spectroscopic ruler based on a single-molecule optical switch

abstract: We demonstrate a novel all-optical switch consisting of two molecules: a primary fluorophore that can be switched between a fluorescent and a dark state by light of different wavelengths, and a secondary chromophore that facilitates switching. The interaction between the two molecules exhibits a distance dependence much steeper than that of Forster resonance energy transfer. This enables the switch to act as a ruler with the capability to probe distances difficult to access by other spectroscopic methods, thus presenting a new tool for the study of biomolecules at the single-molecule level.

categories: q-bio.QM physics.bio-ph physics.chem-ph q-bio.BM


title: Accelerating, hyper-accelerating, and decelerating probabilistic networks

abstract: Many growing networks possess accelerating statistics where the number of links added with each new node is an increasing function of network size so the total number of links increases faster than linearly with network size. In particular, biological networks can display a quadratic growth in regulator number with genome size even while remaining sparsely connected. These features are mutually incompatible in standard treatments of network theory which typically require that every new network node possesses at least one connection. To model sparsely connected networks, we generalize existing approaches and add each new node with a probabilistic number of links to generate either accelerating, hyper-accelerating, or even decelerating network statistics in different regimes. Under preferential attachment for example, slowly accelerating networks display stationary scale-free statistics relatively independent of network size while more rapidly accelerating networks display a transition from scale-free to exponential statistics with network growth. Such transitions explain, for instance, the evolutionary record of single-celled organisms which display strict size and complexity limits.


title: Resolution exchange simulation

abstract: We extend replica exchange simulation in two ways, and apply our approaches to biomolecules. The first generalization permits exchange simulation between models of differing resolution -- i.e., between detailed and coarse-grained models. Such ``resolution exchange'' can be applied to molecular systems or spin systems. The second extension is to ``pseudo-exchange'' simulations, which require little CPU usage for most levels of the exchange ladder and also substantially reduces the need for overlap between levels. Pseudo exchanges can be used in either replica or resolution exchange simulations. We perform efficient, converged simulations of a 50-atom peptide to illustrate the new approaches.


title: A stochastic approach to multi-gene expression dynamics

abstract: In the last years, tens of thousands gene expression profiles for cells of several organisms have been monitored. Gene expression is a complex transcriptional process where mRNA molecules are translated into proteins, which control most of the cell functions. In this process, the correlation among genes is crucial to determine the specific functions of genes. Here, we propose a novel multi-dimensional stochastic approach to deal with the gene correlation phenomena. Interestingly, our stochastic framework suggests that the study of the gene correlation requires only one theoretical assumption -Markov property- and the experimental transition probability, which characterizes the gene correlation system. Finally, a gene expression experiment is proposed for future applications of the model.