title: Virus evolution : the emergence of new ideas (and re-emergence of old ones)

abstract: Reputed intractable, the question of the origin of viruses has long been neglected. In the modern literature 'Virus evolution' has come to refer to study more akin to population genetics, such as the world-wide scrutiny on new polymorphisms appearing daily in the H5N1 avian flu virus [1], than to the fundamental interrogation: where do viruses come from? This situation is now rapidly changing, due to the coincidence of bold new ideas (and sometimes the revival of old ones), the unexpected features exhibited by recently isolated spectacular viruses [2] (see at URL: www.giantvirus.org), as well as the steady increase of genomic sequences for 'regular' viruses and cellular organisms enhancing the power of comparative genomics [3]. After being considered non-living and relegated in the wings by a majority of biologists, viruses are now pushed back on the center stage: they might have been at the origin of DNA, of the eukaryotic cell, and even of today's partition of biological organisms into 3 domains of life: bacteria, archaea and eukarya. Here, I quickly survey some of the recent discoveries and the new evolutionary thoughts they have prompted, before adding to the confusion with one interrogation of my own: what if we totally missed the true nature of (at least some) viruses?


title: On the weak solutions of the McKendrick equation: Existence of demography cycles

abstract: We develop the qualitative theory of the solutions of the McKendrick partial differential equation of population dynamics. We calculate explicitly the weak solutions of the McKendrick equation and of the Lotka renewal integral equation with time and age dependent birth rate. Mortality modulus is considered age dependent. We show the existence of demography cycles. For a population with only one reproductive age class, independently of the stability of the weak solutions and after a transient time, the temporal evolution of the number of individuals of a population is always modulated by a time periodic function. The periodicity of the cycles is equal to the age of the reproductive age class, and a population retains the memory from the initial data through the amplitude of oscillations. For a population with a continuous distribution of reproductive age classes, the amplitude of oscillation is damped. The periodicity of the damped cycles is associated with the age of the first reproductive age class. Damping increases as the dispersion of the fertility function around the age class with maximal fertility increases. In general, the period of the demography cycles is associated with the time that a species takes to reach the reproductive maturity.


title: Spreading of infectious diseases on heterogeneous populations: multi-type network approach

abstract: I study the spreading of infectious diseases on heterogeneous populations. I represent the population structure by a contact-graph where vertices represent agents and edges represent disease transmission channels among them. The population heterogeneity is taken into account by the agent's subdivision in types and the mixing matrix among them. I introduce a type-network representation for the mixing matrix allowing an intuitive understanding of the mixing patterns and the analytical calculations. Using an iterative approach I obtain recursive equations for the probability distribution of the outbreak size as a function of time. I demonstrate that the expected outbreak size and its progression in time are determined by the largest eigenvalue of the reproductive number matrix and the characteristic distance between agents on the contact-graph. Finally, I discuss the impact of intervention strategies to halt epidemic outbreaks. This work provides both a qualitative understanding and tools to obtain quantitative predictions for the spreading dynamics on heterogeneous populations.


title: Pattern-based phylogenetic distance estimation and tree reconstruction

abstract: We have developed an alignment-free method that calculates phylogenetic distances using a maximum likelihood approach for a model of sequence change on patterns that are discovered in unaligned sequences. To evaluate the phylogenetic accuracy of our method, and to conduct a comprehensive comparison of existing alignment-free methods (freely available as Python package decaf+py at http://www.bioinformatics.org.au), we have created a dataset of reference trees covering a wide range of phylogenetic distances. Amino acid sequences were evolved along the trees and input to the tested methods; from their calculated distances we infered trees whose topologies we compared to the reference trees. We find our pattern-based method statistically superior to all other tested alignment-free methods on this dataset. We also demonstrate the general advantage of alignment-free methods over an approach based on automated alignments when sequences violate the assumption of collinearity. Similarly, we compare methods on empirical data from an existing alignment benchmark set that we used to derive reference distances and trees. Our pattern-based approach yields distances that show a linear relationship to reference distances over a substantially longer range than other alignment-free methods. The pattern-based approach outperforms alignment-free methods and its phylogenetic accuracy is statistically indistinguishable from alignment-based distances.


title: Hierarchical modularity of nested bow-ties in metabolic networks

abstract: The exploration of the structural topology and the organizing principles of genome-based large-scale metabolic networks is essential for studying possible relations between structure and functionality of metabolic networks. Topological analysis of graph models has often been applied to study the structural characteristics of complex metabolic networks.In this work, metabolic networks of 75 organisms were investigated from a topological point of view. Network decomposition of three microbes (Escherichia coli, Aeropyrum pernix and Saccharomyces cerevisiae) shows that almost all of the sub-networks exhibit a highly modularized bow-tie topological pattern similar to that of the global metabolic networks. Moreover, these small bow-ties are hierarchically nested into larger ones and collectively integrated into a large metabolic network, and important features of this modularity are not observed in the random shuffled network. In addition, such a bow-tie pattern appears to be present in certain chemically isolated functional modules and spatially separated modules including carbohydrate metabolism, cytosol and mitochondrion respectively. The highly modularized bow-tie pattern is present at different levels and scales, and in different chemical and spatial modules of metabolic networks, which is likely the result of the evolutionary process rather than a random accident. Identification and analysis of such a pattern is helpful for understanding the design principles and facilitate the modelling of metabolic networks.


title: Methods of robustness analysis for Boolean models of gene control networks

abstract: As a discrete approach to genetic regulatory networks, Boolean models provide an essential qualitative description of the structure of interactions among genes and proteins. Boolean models generally assume only two possible states (expressed or not expressed) for each gene or protein in the network as well as a high level of synchronization among the various regulatory processes. In this paper, we discuss and compare two possible methods of adapting qualitative models to incorporate the continuous-time character of regulatory networks. The first method consists of introducing asynchronous updates in the Boolean model. In the second method, we adopt the approach introduced by L. Glass to obtain a set of piecewise linear differential equations which continuously describe the states of each gene or protein in the network. We apply both methods to a particular example: a Boolean model of the segment polarity gene network of Drosophila melanogaster. We analyze the dynamics of the model, and provide a theoretical characterization of the model's gene pattern prediction as a function of the timescales of the various processes.


title: Microarray Data Management. An Enterprise Information Approach: Implementations and Challenges

abstract: The extraction of information form high-throughput experiments is a key aspect of modern biology. Early in the development of microarray technology, researchers recognized that the size of the datasets and the limitations of both computational and visualization techniques restricted their ability to find the biological meaning hidden in the data. In addition, most researchers wanted to make their datasets accessible to others. This resulted in the development of new and advanced data storage, analysis, and visualization tools enabling the cross-platform validation of the experiments and the identification of previously undetected patterns. In order to reap the benefits of this microarray data, researchers have needed to implement database management systems providing integration of different experiments and data types. Moreover, it was necessary to standardize the basic data structure and experimental techniques for the standardization of microarray platforms. In this chapter, we introduce the reader to the major concepts related to the use of controlled vocabularies (ontologies), the definition of Minimum Information About a Microarray Experiment (MIAME) and provide an overview of different microarray data management strategies in use today. We summarize the main characteristics of microarray data storage and sharing strategies including warehouses, datamarts, and federations. The fundamental challenges involved in the distribution, and retrieval of microarray data are presented, along with an overview of some emerging technologies.


title: The role of synaptic facilitation in coincidence spike detection

abstract: Using a realistic model of activity dependent dynamical synapses and a standard integrate and fire neuron model we study, both analytically and numerically, the conditions in which a postsynaptic neuron efficiently detects temporal coincidences of spikes arriving at certain frequency from N different afferents. We extend a previous work that only considers synaptic depression as the most important mechanism in the transmission of information through synapses, to a more general situation including also synaptic facilitation. Our study shows that: 1) facilitation enhances the detection of correlated signals arriving from a subset of presynaptic excitatory neurons, with different degrees of correlation among this subset, and 2) the presence of facilitation allows for a better detection of firing rate changes. Finally, we also observed that facilitation determines the existence of an optimal input frequency which allows the best performance for a wide (maximum) range of the neuron firing threshold. This optimal frequency can be controlled by means of facilitation parameters.


title: Effective affinities in microarray data

abstract: In the past couple of years several studies have shown that hybridization in Affymetrix DNA microarrays can be rather well understood on the basis of simple models of physical chemistry. In the majority of the cases a Langmuir isotherm was used to fit experimental data. Although there is a general consensus about this approach, some discrepancies between different studies are evident. For instance, some authors have fitted the hybridization affinities from the microarray fluorescent intensities, while others used affinities obtained from melting experiments in solution. The former approach yields fitted affinities that at first sight are only partially consistent with solution values. In this paper we show that this discrepancy exists only superficially: a sufficiently complete model provides effective affinities which are fully consistent with those fitted to experimental data. This link provides new insight on the relevant processes underlying the functioning of DNA microarrays.


title: Physics-based analysis of Affymetrix microarray data

abstract: We analyze publicly available data on Affymetrix microarrays spike-in experiments on the human HGU133 chipset in which sequences are added in solution at known concentrations. The spike-in set contains sequences of bacterial, human and artificial origin. Our analysis is based on a recently introduced molecular-based model [E. Carlon and T. Heim, Physica A 362, 433 (2006)] which takes into account both probe-target hybridization and target-target partial hybridization in solution. The hybridization free energies are obtained from the nearest-neighbor model with experimentally determined parameters. The molecular-based model suggests a rescaling that should result in a "collapse" of the data at different concentrations into a single universal curve. We indeed find such a collapse, with the same parameters as obtained before for the older HGU95 chip set. The quality of the collapse varies according to the probe set considered. Artificial sequences, chosen by Affymetrix to be as different as possible from any other human genome sequence, generally show a much better collapse and thus a better agreement with the model than all other sequences. This suggests that the observed deviations from the predicted collapse are related to the choice of probes or have a biological origin, rather than being a problem with the proposed model.


title: Superstability of the yeast cell cycle dynamics: Ensuring causality in the presence of biochemical stochasticity

abstract: Gene regulatory dynamics is governed by molecular processes and therefore exhibits an inherent stochasticity. However, for the survival of an organism it is a strict necessity that this intrinsic noise does not prevent robust functioning of the system. It is still an open question how dynamical stability is achieved in biological systems despite the omnipresent fluctuations. In this paper we investigate the cell-cycle of the budding yeast Saccharomyces cerevisiae as an example of a well-studied organism. We study a genetic network model of eleven genes that coordinate the cell-cycle dynamics using a modeling framework which generalizes the concept of discrete threshold dynamics. By allowing for fluctuations in the transcription/translation times, we introduce noise in the model, accounting for the effects of biochemical stochasticity. We study the dynamical attractor of the cell cycle and find a remarkable robustness against fluctuations of this kind. We identify mechanisms that ensure reliability in spite of fluctuations: 'Catcher' states and persistence of activity levels contribute significantly to the stability of the yeast cell cycle despite the inherent stochasticity.


title: Chaotic itinerancy, temporal segmentation and spatio-temporal combinatorial codes

abstract: We study a deterministic dynamics with two time scales in a continuous state attractor network. To the usual (fast) relaxation dynamics towards point attractors (``patterns'') we add a slow coupling dynamics that makes the visited patterns to loose stability leading to an itinerant behavior in the form of punctuated equilibria. One finds that the transition frequency matrix between patterns shows non-trivial statistical properties in the chaotic itinerant regime. We show that mixture input patterns can be temporally segmented by the itinerant dynamics. The viability of a combinatorial spatio-temporal neural code is also demonstrated.


title: Stochastic Model of Yeast Cell Cycle Network

abstract: Biological functions in living cells are controlled by protein interaction and genetic networks. These molecular networks should be dynamically stable against various fluctuations which are inevitable in the living world. In this paper, we propose and study a stochastic model for the network regulating the cell cycle of the budding yeast. The stochasticity in the model is controlled by a temperature-like parameter $\beta$. Our simulation results show that both the biological stationary state and the biological pathway are stable for a wide range of "temperature". There is, however, a sharp transition-like behavior at $\beta_c$, below which the dynamics is dominated by noise. We also define a pseudo energy landscape for the system in which the biological pathway can be seen as a deep valley.


title: Stable Propagation of a Burst Through a One-Dimensional Homogeneous Excitatory Chain Model of Songbird Nucleus HVC

abstract: We demonstrate numerically that a brief burst consisting of two to six spikes can propagate in a stable manner through a one-dimensional homogeneous feedforward chain of non-bursting neurons with excitatory synaptic connections. Our results are obtained for two kinds of neuronal models, leaky integrate-and-fire (LIF) neurons and Hodgkin-Huxley (HH) neurons with five conductances. Over a range of parameters such as the maximum synaptic conductance, both kinds of chains are found to have multiple attractors of propagating bursts, with each attractor being distinguished by the number of spikes and total duration of the propagating burst. These results make plausible the hypothesis that sparse precisely-timed sequential bursts observed in projection neurons of nucleus HVC of a singing zebra finch are intrinsic and causally related.


title: A model of hyphal tip growth involving microtubule-based transport

abstract: We propose a simple model for mass transport within a fungal hypha and its subsequent growth. Inspired by the role of microtubule-transported vesicles, we embody the internal dynamics of mass inside a hypha with mutually excluding particles progressing stochastically along a growing one-dimensional lattice. The connection between long range transport of materials for growth, and the resulting extension of the hyphal tip has not previously been addressed in the modelling literature. We derive and analyse mean-field equations for the model and present a phase diagram of its steady state behaviour, which we compare to simulations. We discuss our results in the context of the filamentous fungus, Neurospora crassa.


title: Complex reaction noise in molecular quasispecies model

abstract: We have derived exact Langevin equations for a model of quasispecies dynamics. The inherent multiplicative reaction noise is complex and its statistical properties are specified completely. The numerical simulation of the complex Langevin equations is carried out using the Cholesky decomposition for the noise covariance matrix. This internal noise, which is due to diffusion-limited reactions, produces unavoidable spatio-temporal density fluctuations about the mean field value. In two dimensions, this noise strictly vanishes only in the perfectly mixed limit, a situation difficult to attain in practice.


title: Gene Function Classification Using Bayesian Models with Hierarchy-Based Priors

abstract: We investigate the application of hierarchical classification schemes to the annotation of gene function based on several characteristics of protein sequences including phylogenic descriptors, sequence based attributes, and predicted secondary structure. We discuss three Bayesian models and compare their performance in terms of predictive accuracy. These models are the ordinary multinomial logit (MNL) model, a hierarchical model based on a set of nested MNL models, and a MNL model with a prior that introduces correlations between the parameters for classes that are nearby in the hierarchy. We also provide a new scheme for combining different sources of information. We use these models to predict the functional class of Open Reading Frames (ORFs) from the E. coli genome. The results from all three models show substantial improvement over previous methods, which were based on the C5 algorithm. The MNL model using a prior based on the hierarchy outperforms both the non-hierarchical MNL model and the nested MNL model. In contrast to previous attempts at combining these sources of information, our approach results in a higher accuracy rate when compared to models that use each data source alone. Together, these results show that gene function can be predicted with higher accuracy than previously achieved, using Bayesian models that incorporate suitable prior information.


title: How T-cells use large deviations to recognize foreign antigens

abstract: A stochastic model for the activation of T-cells is analysed. T-cells are part of the immune system and recognize foreign antigens against a background of the body's own molecules. The model under consideration is a slight generalization of a model introduced by Van den Berg, Rand and Burroughs in 2001, and is capable of explaining how this recognition works on the basis of rare stochastic events. With the help of a refined large deviation theorem and numerical evaluation it is shown that, for a wide range of parameters, T-cells can distinguish reliably between foreign antigens and self-antigens.


title: Kinetic Accessibility of Buried DNA Sites in Nucleosomes

abstract: Using a theoretical model for spontaneous partial DNA unwrapping from histones, we study the transient exposure of protein-binding DNA sites within nucleosomes. We focus on the functional dependence of the rates for site exposure and reburial on the site position, which is measurable experimentally and pertinent to gene regulation. We find the dependence to be roughly described by a random walker model. Close inspection reveals a surprising physical effect of flexibility-assisted barrier crossing, which we characterize within a toy model, the "semiflexible Brownian rotor."


title: PIPER: An FFT-based Protein Docking Program with Pairwise Potentials

abstract: The Fast Fourier Transform (FFT) correlation approach to protein-protein docking can evaluate the energies of billions of docked conformations on a grid if the energy is described in the form of a correlation function. Here, this restriction is removed, and the approach is efficiently used with pairwise interactions potentials that substantially improve the docking results. The basic idea is approximating the interaction matrix by its eigenvectors corresponding to the few dominant eigenvalues, resulting in an energy expression written as the sum of a few correlation functions, and solving the problem by repeated FFT calculations. In addition to describing how the method is implemented, we present a novel class of structure based pairwise intermolecular potentials. The DARS (Decoys As the Reference State) potentials are extracted from structures of protein-protein complexes and use large sets of docked conformations as decoys to derive atom pair distributions in the reference state. The current version of the DARS potential works well for enzyme-inhibitor complexes. With the new FFT-based program, DARS provides much better docking results than the earlier approaches, in many cases generating 50\% more near-native docked conformations. Although the potential is far from optimal for antibody-antigen pairs, the results are still slightly better than those given by an earlier FFT method. The docking program PIPER is freely available for non-commercial applications.


title: Episodic synchronization in dynamically driven neurons

abstract: We examine the response of type II excitable neurons to trains of synaptic pulses, as a function of the pulse frequency and amplitude. We show that the resonant behavior characteristic of type II excitability, already described for harmonic inputs, is also present for pulsed inputs. With this in mind, we study the response of neurons to pulsed input trains whose frequency varies continuously in time, and observe that the receiving neuron synchronizes episodically to the input pulses, whenever the pulse frequency lies within the neuron's locking range. We propose this behavior as a mechanism of rate-code detection in neuronal populations. The results are obtained both in numerical simulations of the Morris-Lecar model and in an electronic implementation of the FitzHugh-Nagumo system, evidencing the robustness of the phenomenon.


title: Laws in Darwinian Evolutionary Theory

abstract: In the present article the recent works to formulate laws in Darwinian evolutionary dynamics are discussed. Although there is a strong consensus that general laws in biology may exist, opinions opposing such suggestion are abundant. Based on recent progress in both mathematics and biology, another attempt to address this issue is made in the present article. Specifically, three laws which form a mathematical framework for the evolutionary dynamics in biology are postulated. The second law is most quantitative and is explicitly expressed in the unique form of a stochastic differential equation. Salient features of Darwinian evolutionary dynamics are captured by this law: the probabilistic nature of evolution, ascendancy, and the adaptive landscape. Four dynamical elements are introduced in this formulation: the ascendant matrix, the transverse matrix, the Wright evolutionary potential, and the stochastic drive. The first law may be regarded as a special case of the second law. It gives the reference point to discuss the evolutionary dynamics. The third law describes the relationship between the focused level of description to its lower and higher ones, and defines the dichotomy of deterministic and stochastic drives. It is an acknowledgement of the hierarchical structure in biology. A new interpretation of Fisher's fundamental theorem of natural selection is provided in terms of the F-Theorem. The proposed laws are based on continuous representation in both time and population. Their generic nature is demonstrated through their equivalence to classical formulations. The present three laws appear to provide a coherent framework for the further development of the subject.

categories: q-bio.PE cond-mat.stat-mech cs.NE math.OC nlin.AO physics.bio-ph q-bio.QM


title: Spatial embedding of neuronal trees modeled by diffusive growth

abstract: The relative importance of the intrinsic and extrinsic factors determining the variety of geometric shapes exhibited by dendritic trees remains unclear. This question was addressed by developing a model of the growth of dendritic trees based on diffusion-limited aggregation process. The model reproduces diverse neuronal shapes (i.e., granule cells, Purkinje cells, the basal and apical dendrites of pyramidal cells, and the axonal trees of interneurons) by changing only the size of the growth area, the time span of pruning, and the spatial concentration of 'neurotrophic particles'. Moreover, the presented model shows how competition between neurons can affect the shape of the dendritic trees. The model reveals that the creation of complex (but reproducible) dendrite-like trees does not require precise guidance or an intrinsic plan of the dendrite geometry. Instead, basic environmental factors and the simple rules of diffusive growth adequately account for the spatial embedding of different types of dendrites observed in the cortex. An example demonstrating the broad applicability of the algorithm to model diverse types of tree structures is also presented. Key words: Diffusion-limited aggregation; Neuronal morphology; Dendrites; Growth model; tree; dendritic geometry.


title: Theory of Interaction of Memory Patterns in Layered Associative Networks

abstract: A synfire chain is a network that can generate repeated spike patterns with millisecond precision. Although synfire chains with only one activity propagation mode have been intensively analyzed with several neuron models, those with several stable propagation modes have not been thoroughly investigated. By using the leaky integrate-and-fire neuron model, we constructed a layered associative network embedded with memory patterns. We analyzed the network dynamics with the Fokker-Planck equation. First, we addressed the stability of one memory pattern as a propagating spike volley. We showed that memory patterns propagate as pulse packets. Second, we investigated the activity when we activated two different memory patterns. Simultaneous activation of two memory patterns with the same strength led the propagating pattern to a mixed state. In contrast, when the activations had different strengths, the pulse packet converged to a two-peak state. Finally, we studied the effect of the preceding pulse packet on the following pulse packet. The following pulse packet was modified from its original activated memory pattern, and it converged to a two-peak state, mixed state or non-spike state depending on the time interval.


title: Form follows function -- how PufX increases the efficiency of the light-harvesting complexes of Rhodobacter sphaeroides

abstract: Some species of purple bacteria as, e.g., Rhodobacter sphaeroides contain the protein PufX. Concurrently, the light harvesting complexes 1 (LH1) form dimers of open rings. In mutants without PufX, the LH1s are closed rings and photosynthesis breaks down, because the ubiquinone exchange at the reaction center is blocked. Thus, PufX is regarded essential for quinone exchange. In contrast to this view, which implicitly treats the LH1s as obstacles to photosynthesis, we propose that the primary purpose of PufX is to improve the efficiency of light harvesting by inducing the LH1 dimerization. Calculations with a dipole model, which compare the photosynthetic efficiency of various configurations of monomeric and dimeric core complexes, show that the dimer can absorb photons directly into the RC about 30% more efficient, when related to the number of bacteriochlorophylls, but that the performance of the more sophisticated dimeric LH1 antenna degrades faster with structural perturbations. The calculations predict an optimal orientation of the reaction centers relative to the LH1 dimer, which agrees well with the experimentally found configuration. For the increased required rigidity of the dimer additional modifications of the LH1 subunits are necessary, which would lead to the observed ubiquinone blockage, when PufX is missing.


title: Shorelines of islands of tractability: Algorithms for parsimony and minimum perfect phylogeny haplotyping problems

abstract: The problem Parsimony Haplotyping (PH) asks for the smallest set of haplotypes which can explain a given set of genotypes, and the problem Minimum Perfect Phylogeny Haplotyping (MPPH) asks for the smallest such set which also allows the haplotypes to be embedded in a perfect phylogeny, an evolutionary tree with biologically-motivated restrictions. For PH, we extend recent work by further mapping the interface between ``easy'' and ``hard'' instances, within the framework of (k,l)-bounded instances where the number of 2's per column and row of the input matrix is restricted. By exploring, in the same way, the tractability frontier of MPPH we provide the first concrete, positive results for this problem, and the algorithms underpinning these results offer new insights about how MPPH might be further tackled in the future. In addition, we construct for both PH and MPPH polynomial time approximation algorithms, based on properties of the columns of the input matrix. We conclude with an overview of intriguing open problems in PH and MPPH.


title: A method for representing and developing process models

abstract: Scientists investigate the dynamics of complex systems with quantitative models, employing them to synthesize knowledge, to explain observations, and to forecast future system behavior. Complete specification of systems is impossible, so models must be simplified abstractions. Thus, the art of modeling involves deciding which system elements to include and determining how they should be represented. We view modeling as search through a space of candidate models that is guided by model objectives, theoretical knowledge, and empirical data. In this contribution, we introduce a method for representing process-based models that facilitates the discovery of models that explain observed behavior. This representation casts dynamic systems as interacting sets of processes that act on entities. Using this approach, a modeler first encodes relevant ecological knowledge into a library of generic entities and processes, then instantiates these theoretical components, and finally assembles candidate models from these elements. We illustrate this methodology with a model of the Ross Sea ecosystem.


title: Nestedness in mutualistic networks emerges from a modified rule of preferential attachment

abstract: Replacement of papers q-bio/0605026 and q-bio/0607043 by a single manuscript (q-bio/0701029) follows suggestions from a magazine's referee


title: A new stimulus approach in the search for biological nonlocality

abstract: Preliminary research in the area of biophysics appears to indicate the existence of quantum entanglement and nonlocality at the biological level, both for human subjects and for neurons derived from human neural stem cells. The lack of clear and replicable findings, especially where human subjects are concerned, results in conflicting and marginal correlations between stimulated and nonstimulated subject's brainwaves. It is proposed to go beyond the use of just patterned photostimulation, to encompass either transcranial magnetic stimulation (TMS) or acupuncture, so that one can achieve a simultaneity of several different brain events in both stimulated and nonstimulated subjects, rather than just one EEG-type event.


title: Targeted Projection Pursuit for Gene Expression Data Classification and Visualisation

abstract: We present a novel method for finding low dimensional views of high dimensional data: Targeted Projection Pursuit. The method proceeds by finding projections of the data that best approximate a target view. Two versions of the method are introduced; one version based on Procrustes analysis and one based on a single layer perceptron. These versions are capable of finding orthogonal or non-orthogonal projections respectively. The method is quantitatively and qualitatively compared with other dimension reduction techniques. It is shown to find two-dimensional views that display the classification of cancers from gene expression data with a visual separation equal to, or better than, existing dimension reduction techniques.


title: Thriving at high hydrostatic pressure: the example of ammonoids (extinct cephalopods)

abstract: Ammonoids are a group of extinct mollusks belonging to the same class of the living genus Nautilus (Cephalopoda). In both Nautili and ammonoids, the (usually planospiral) shell is divided into chambers separated by septa that during the lifetime were filled with gas at atmospheric pressure. The intersection of septa with the external shell generates a curve called the suture line, which in living and most fossil Nautili is fairly uncomplicated. In contrast, suture lines of ancient ammonoid were gently curved and during the evolution of the group became highly complex, in some cases so extensively frilled to be considerable as fractal curves. Numerous theories have been put forward to explain the complexity of suture ammonoid lines. Calculations presented here lend support to the hypothesis that complex suture lines aided in counteracting the effect of the external water pressure. Additionally, it is found that complex suture lines diminished shell shrinkage caused by water pressure, and thus aided improve buoyancy. Understanding the reason for complex sutures in ammonoids does not only represent an important issue in paleobiology, but is also a challenging problem in the resistance of complex mechanical structures subjected to high pressure.


title: Aging cellular networks: chaperones as major participants

abstract: We increasingly rely on the network approach to understand the complexity of cellular functions. Chaperones (heat shock proteins) are key "networkers", which have among their functions to sequester and repair damaged protein. In order to link the network approach and chaperones with the aging process, we first summarize the properties of aging networks suggesting a "weak link theory of aging". This theory suggests that age-related random damage primarily affects the overwhelming majority of the low affinity, transient interactions (weak links) in cellular networks leading to increased noise, destabilization and diversity. These processes may be further amplified by age-specific network remodelling and by the sequestration of weakly linked cellular proteins to protein aggregates of aging cells. Chaperones are weakly linked hubs [i.e., network elements with a large number of connections] and inter-modular bridge elements of protein-protein interaction, signalling and mitochondrial networks. As aging proceeds, the increased overload of damaged proteins is an especially important element contributing to cellular disintegration and destabilization. Additionally, chaperone overload may contribute to the increase of "noise" in aging cells, which leads to an increased stochastic resonance resulting in a deficient discrimination between signals and noise. Chaperone- and other multi-target therapies, which restore the missing weak links in aging cellular networks, may emerge as important anti-aging interventions.


title: Structural Change of Myosin Motor Domain and Nucleotide Dissociation

abstract: We investigated the structural relaxation of myosin motor domain from the pre-power stroke state to the near-rigor state using molecular dynamics simulation of a coarse-grained protein model. To describe the structural change, we propose a "dual Go-model," a variant of the Go-like model that has two reference structures. The nucleotide dissociation process is also studied by introducing a coarse-grained nucleotide in the simulation. We found that the myosin structural relaxation toward the near-rigor conformation cannot be completed before the nucleotide dissociation. Moreover, the relaxation and the dissociation occurred cooperatively when the nucleotide was tightly bound to the myosin head. The result suggested that the primary role of the nucleotide is to suppress the structural relaxation.


title: Reverse-engineering of polynomial dynamical systems

abstract: Multivariate polynomial dynamical systems over finite fields have been studied in several contexts, including engineering and mathematical biology. An important problem is to construct models of such systems from a partial specification of dynamic properties, e.g., from a collection of state transition measurements. Here, we consider static models, which are directed graphs that represent the causal relationships between system variables, so-called wiring diagrams. This paper contains an algorithm which computes all possible minimal wiring diagrams for a given set of state transition measurements. The paper also contains several statistical measures for model selection. The algorithm uses primary decomposition of monomial ideals as the principal tool. An application to the reverse-engineering of a gene regulatory network is included. The algorithm and the statistical measures are implemented in Macaulay2 and are available from the authors.


title: Impact of observational incompleteness on the structural properties of protein interaction networks

abstract: The observed structure of protein interaction networks is corrupted by many false positive/negative links. This observational incompleteness is abstracted as random link removal and a specific, experimentally motivated (spoke) link rearrangement. Their impact on the structural properties of gene-duplication-and-mutation network models is studied. For the degree distribution a curve collapse is found, showing no sensitive dependence on the link removal/rearrangement strengths and disallowing a quantitative extraction of model parameters. The spoke link rearrangement process moves other structural observables, like degree correlations, cluster coefficient and motif frequencies, closer to their counterparts extracted from the yeast data. This underlines the importance to take a precise modeling of the observational incompleteness into account when network structure models are to be quantitatively compared to data.


title: Optimization over a class of tree shape statistics

abstract: Tree shape statistics quantify some aspect of the shape of a phylogenetic tree. They are commonly used to compare reconstructed trees to evolutionary models and to find evidence of tree reconstruction bias. Historically, to find a useful tree shape statistic, formulas have been invented by hand and then evaluated for utility. This article presents the first method which is capable of optimizing over a class of tree shape statistics, called Binary Recursive Tree Shape Statistics (BRTSS). After defining the BRTSS class, a set of algebraic expressions is defined which can be used in the recursions. The tree shape statistics definable using these expressions in the BRTSS is very general, and includes many of the statistics with which phylogenetic researchers are already familiar. We then present a practical genetic algorithm which is capable of performing optimization over BRTSS given any objective function. The chapter concludes with a successful application of the methods to find a new statistic which indicates a significant difference between two distributions on trees which were previously postulated to have similar properties.


title: A continuous model for cell sorting

abstract: The differential Adhesion Hypothesis (DAH) is a theory of the organization of cells within a tissue. In this study we introduce a stochastic model supporting the DAH, that can be seen as a continuous version of a discrete model of Graner and Glazier. Our approach is based on the mathematical framework of Gibbsian marked point processes. We provide a Markov chain Monte Carlo algorithm that can reproduce classical biological patterns, and we propose an estimation procedure for a parameter that quantifies the strength of adhesion between cells. This procedure is tested through simulations.


title: Collective evolution and the genetic code

abstract: A dynamical theory for the evolution of the genetic code is presented, which accounts for its universality and optimality. The central concept is that a variety of collective, but non-Darwinian, mechanisms likely to be present in early communal life generically lead to refinement and selection of innovation-sharing protocols, such as the genetic code. Our proposal is illustrated using a simplified computer model, and placed within the context of a sequence of transitions that early life may have made, prior to the emergence of vertical descent.


title: Mass Fractal Dimension of the Ribosome and Implication of its Dynamic Characteristics

abstract: Self-similar properties of the ribosome in terms of the mass fractal dimension are investigated. We find that both the 30S subunit and the 16S rRNA have fractal dimensions of 2.58 and 2.82, respectively; while the 50S subunit as well as the 23S rRNA has the mass fractal dimension close to 3, implying a compact three dimensional macromolecule. This finding supports the dynamic and active role of the 30S subunit in the protein synthesis, in contrast to the pass role of the 50S subunit.


title: Neural integrator - a sandpile model

abstract: We investigated a model for the neural integrator based on hysteretic units connected by positive feedback. Hysteresis is assumed to emerge from the intrinsic properties of the cells. We consider the recurrent networks containing either bistable or multistable neurons. We apply our analysis to the oculomotor velocity-to-position neural integrator that calculates the eye positions from the inputs that carry information about eye angular velocity. Using the analysis of the system in the parameter space we show the following. The direction of hysteresis in the neuronal response may be reversed for the system with recurrent connections compared to the case of unconnected neurons. Thus, for the NMDA receptor based bistability the firing rates after ON saccades may be higher than after OFF saccades for the same eye position. We suggest that this is an emergent property due to the presence of global recurrent feedback. The reversal of hysteresis occurs only when the size of hysteresis differs from neuron to neuron. We also relate the macroscopic leak time-constant of the integrator to the rate of microscopic spontaneous noise-driven transitions in the hysteretic units. Finally, we argue that the presence of neurons with small hysteresis may remove the threshold for integration.


title: A toy model of polymer stretching

abstract: We present an extremely simplified model of multiple-domains polymer stretching in an atomic force microscopy experiment. We portray each module as a binary set of contacts and decompose the system energy into a harmonic term (the cantilever) and long-range interactions terms inside each domain. Exact equilibrium computations and Monte Carlo simulations qualitatively reproduce the experimental saw-tooth pattern of force-extension profiles, corresponding (in our model) to first-order phase transitions. We study the influence of the coupling induced by the cantilever and the pulling speed on the relative heights of the force peaks. The results suggest that the increasing height of the critical force for subsequent unfolding events is an out-of-equilibrium effect due to a finite pulling speed. The dependence of the average unfolding force on the pulling speed is shown to reproduce the experimental logarithmic law.


title: Epidemic outbreaks on structured populations

abstract: Our chances to halt epidemic outbreaks rely on how accurately we represent the population structure underlying the disease spread. When analyzing global epidemics this force us to consider metapopulation models taking into account intra- and inter-community interactions. Recently Watts et al introduced a metapopulation model which accounts for several features observed in real outbreaks [Watts et al, PNAS 102, 11157 (2005)]. In this work I provide an analytical solution to this model, enhancing our understanding of the model and the epidemic outbreaks it represents. First, I demonstrate that depending on the intra-community expected outbreak size and the fraction of social bridges the epidemic outbreaks die out or there is a finite probability to observe a global epidemics. Second, I show that the global scenario is characterized by resurgent epidemics, their number increasing with increasing the intra-community average distance between individuals. Finally, I present empirical data for the AIDS epidemics supporting the model predictions.

categories: q-bio.PE cond-mat.stat-mech physics.bio-ph q-bio.QM


title: Thermodynamics of Neutral Protein Evolution

abstract: Naturally evolving proteins gradually accumulate mutations while continuing to fold to thermodynamically stable native structures. This process of neutral protein evolution is an important mode of genetic change, and forms the basis for the molecular clock. Here we present a mathematical theory that predicts the number of accumulated mutations, the index of dispersion, and the distribution of stabilities in an evolving protein population from knowledge of the stability effects (ddG values) for single mutations. Our theory quantitatively describes how neutral evolution leads to marginally stable proteins, and provides formulae for calculating how fluctuations in stability cause an overdispersion of the molecular clock. It also shows that the structural influences on the rate of sequence evolution that have been observed in earlier simulations can be calculated using only the single-mutation ddG values. We consider both the case when the product of the population size and mutation rate is small and the case when this product is large, and show that in the latter case proteins evolve excess mutational robustness that is manifested by extra stability and increases the rate of sequence evolution. Our basic method is to treat protein evolution as a Markov process constrained by a minimal requirement for stable folding, enabling an evolutionary description of the proteins solely in terms of the experimentally measureable ddG values. All of our theoretical predictions are confirmed by simulations with model lattice proteins. Our work provides a mathematical foundation for understanding how protein biophysics helps shape the process of evolution.


title: Coexistence versus extinction in the stochastic cyclic Lotka-Volterra model

abstract: Cyclic dominance of species has been identified as a potential mechanism to maintain biodiversity, see e.g. B. Kerr, M. A. Riley, M. W. Feldman and B. J. M. Bohannan [Nature {\bf 418}, 171 (2002)] and B. Kirkup and M. A. Riley [Nature {\bf 428}, 412 (2004)]. Through analytical methods supported by numerical simulations, we address this issue by studying the properties of a paradigmatic non-spatial three-species stochastic system, namely the `rock-paper-scissors' or cyclic Lotka-Volterra model. While the deterministic approach (rate equations) predicts the coexistence of the species resulting in regular (yet neutrally stable) oscillations of the population densities, we demonstrate that fluctuations arising in the system with a \emph{finite number of agents} drastically alter this picture and are responsible for extinction: After long enough time, two of the three species die out. As main findings we provide analytic estimates and numerical computation of the extinction probability at a given time. We also discuss the implications of our results for a broad class of competing population systems.


title: Dynamic scaling regimes of collective decision making

abstract: We investigate a social system of agents faced with a binary choice. We assume there is a correct, or beneficial, outcome of this choice. Furthermore, we assume agents are influenced by others in making their decision, and that the agents can obtain information that may guide them towards making a correct decision. The dynamic model we propose is of nonequilibrium type, converging to a final decision. We run it on random graphs and scale-free networks. On random graphs, we find two distinct regions in terms of the "finalizing time" -- the time until all agents have finalized their decisions. On scale-free networks on the other hand, there does not seem to be any such distinct scaling regions.


title: Emergence of the protein universe in organismal evolution

abstract: In this work we propose a physical model of organismal evolution, where phenotype, organism life expectancy, is directly related to genotype i.e. the stability of its proteins which can be determined exactly in the model. Simulating the model on a computer, we consistently observe the Big Bang scenario whereby exponential population growth ensues as favorable sequence-structure combinations (precursors of stable proteins) are discovered. After that, random diversity of the structural space abruptly collapses into a small set of preferred structural motifs. We observe that protein folds remain stable and abundant in the population at time scales much greater than mutation or organism lifetime, and the distribution of the lifetimes of dominant folds in a population approximately follows a power law. The separation of evolutionary time scales between discovery of new folds and generation of new sequences gives rise to emergence of protein families and superfamilies whose sizes are power-law distributed, closely matching the same distributions for real proteins. The network of structural similarities of the universe of evolved proteins has the same scale-free like character as the actual protein domain universe graph (PDUG). Further, the model predicts that ancient protein domains represent a highly connected and clustered subset of all protein domains, in complete agreement with reality. Together, these results provide a microscopic first principles picture of how protein structures and gene families evolved in the course of evolution.


title: Spontaneous Self-Assembly of Transcription Factor Based Gene Regulation Networks

abstract: We model the transcription factor based regulation network of yeast using a content-based network model that mimicks the recognition of binding motifs on the regulatory regions of the genes. We are thereby able to faithfully reproduce many of the topological features of the gene regulatory network of yeast once the parameters of the yeast genome, in particular the distribution of information coded by the "binding sequences" within the promoter regions is provided as input. The length distribution for the promoter regions is fixed by comparing the k-core analysis of the model network with that of yeast. Our results strongly point to the possibility that the observed topological features are generic to networks formed via sequence-matching between random strings obeying certain length distributions.


title: Different Strategies for Cancer Treatment: Mathematical Modeling

abstract: We formulate and analyze a mathematical model describing immune response to avascular tumor under the influence of immunotherapy and chemotherapy and their combinations as well as vaccine treatments. The effect of vaccine therapy is considered as a parametric perturbation of the model. In the case of a weak immune response, neither immunotherapy nor chemotherapy is found to cause tumor regression to a small size, which would be below the clinically detectable threshold. Numerical simulations show that the efficiency of vaccine therapy depends on both the tumor size and the condition of immune system as well as on the response of the organism to vaccination. In particular, we found that vaccine therapy becomes more effective when used without time delay from a prescribed date of vaccination after surgery and is ineffective without preliminary treatment. For a strong immune response, our model predicts the tumor remission under vaccine therapy. Our study of successive chemo/immuno, immuno/chemo and concurrent chemoimmunotherapy shows that the chemo/immuno sequence is more effective while concurrent chemoimmunotherapy is more sparing.


title: When does division of labor lead to increased system output?

abstract: This paper develops a set of simplified dynamical models with which to explore the conditions under which division of labor leads to optimized system output, as measured by the rate of production of a given product. We consider two models: In the first model, we consider the flow of some resource into a compartment, and the conversion of this resource into some product. In the second model, we consider the resource-limited growth of autoreplicating systems. In this case, we divide the replication and metabolic tasks among different agents. The general features that emerge from our models is that division of labor is favored when the resource to agent ratio is at intermediate values, and when the time cost associated with transporting intermediate products is small compared to characteristic process times. We discuss the results of this paper in the context of simulations with digital life. We also argue that division of labor in the context of our replication model suggests an evolutionary basis for the emergence of the stem-cell-based tissue architecture in complex organisms.


title: Transition States in Protein Folding Kinetics: The Structural Interpretation of Phi-values

abstract: Phi-values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is which structural information Phi-values contain about the transition state of folding. Traditionally, a Phi-value is interpreted as the 'nativeness' of a mutated residue in the transition state. However, this interpretation is often problematic because it assumes a linear relation between the nativeness of the residue and its free-energy contribution. We present here a better structural interpretation of Phi-values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free-energy contributions of helical residues. From a linear fit of our model to the experimental data, we obtain two structural parameters: the extent of helix formation in the transition state, and the nativeness of tertiary interactions in the transition state. We apply our model to all proteins with well-characterized helices for which more than 10 Phi-values are available: protein A, CI2, and protein L. The model captures nonclassical Phi-values <0 or >1 in these helices, and explains how different mutations at a given site can lead to different Phi-values.


title: Sequential recruitment and combinatorial assembling of multiprotein complexes in transcriptional activation

abstract: In human cells, estrogenic signals induce cyclical association and dissociation of specific proteins with the DNA in order to activate transcription of estrogen-responsive genes. These oscillations can be modeled by assuming a large number of sequential reactions represented by linear kinetics with random kinetic rates. Application of the model to experimental data predicts robust binding sequences in which proteins associate with the DNA at several different phases of the oscillation. Our methods circumvent the need to derive detailed kinetic graphs, and are applicable to other oscillatory biological processes involving a large number of sequential steps.


title: Spectrum of genetic diversity and networks of clonal organisms

abstract: Clonal organisms present a particular challenge in population genetics because, in addition to the possible existence of replicates of the same genotype in a given sample, some of the hypotheses and concepts underlying classical population genetics models are irreconcilable with clonality. The genetic structure and diversity of clonal populations was examined using a combination of new tools to analyze microsatellite data in the marine angiosperm Posidonia oceanica. These tools were based on examination of the frequency distribution of the genetic distance among ramets, termed the spectrum of genetic diversity (GDS), and of networks built on the basis of pairwise genetic distances among genets. The properties and topology of networks based on genetic distances showed a "small-world" topology, characterized by a high degree of connectivity among nodes, and a substantial amount of substructure, revealing organization in sub-families of closely related individuals. Keywords: genetic networks; small-world networks; genetic diversity; clonal organisms


title: On Mathematical Theory of Selection: Discrete-Time Models

abstract: Mathematical theory of selection systems is developed for a wide class of dynamical models of inhomogeneous populations with discrete time. The Price equation and its particular case, the Fisher Fundamental theorem of natural selection (FTNS), are well known general results of the theory. It is known that the Price equation being a mathematical identity is not dynamically sufficient, i.e., it does not allow one to predict changes in the mean of a trait beyond the immediate response if only the value of covariance of the trait and fitness at this moment is known. We show that the problem of dynamically insufficiency for the Price equations and for the FTNS can be correctly overcome if to study these equations subject to given initial distribution in framework of exact models of the population dynamics. The knowledge of the entire distribution at any given instant allows making the exact prediction for indefinite time and this prediction dramatically depends on the initial distribution. For these models, the current trait distribution and hence all statistical characteristics of interest, such as mean values of the fitness or any trait could be computed effectively at any time moment.


title: Analytical Models of Forest Dynamics

abstract: The hierarchical system of forest ecosystem models based on the theory of individual-based (structured) models of populations and communities is briefly described. New self-thinning models are integrated with tree stand models within a hierarchical system of models, aimed at assessing (quasi-) stable states of a forest ecosystem under given environment and management. Within the structural model of succession, the ergodic theorems in biology are proved. They assert that the areas of biocenoses in a climax state of association should be proportional to the proper times of their development in a succession line; the other characteristics of biocenoses should be also proportional to their own proper regeneration times in a climax state. A simple theoretical approach to estimating the deviation of the current state of forest ecosystem from its steady state is derived and applied to data from the Prioksko-Terrasny Reserve. The results obtained can serve as a theoretical basis in numerical processing of data on ecological monitoring of undisturbed forest lands.


title: Inhomogeneous maps: the basic theorems and some applications

abstract: Non-linear maps can possess various dynamical behaviors varying from stable steady states and cycles to chaotic oscillations. Most models assume that individuals within a given population are identical ignoring the fundamental role of variation. Here we develop a theory of inhomogeneous maps and apply the general approach to modeling heterogeneous populations with discrete evolutionary time step. We show that the behavior of the inhomogeneous maps may possess complex transition regimes, which depends both on the mean and the variance of the initial parameter distribution. The examples of inhomogeneous models are discussed.


title: Bifurcations of self-similar solutions of the Fokker-Plank Equation

abstract: A class of one-dimensional Fokker-Plank equations having a common stationary solution, which is a power function of the state of the process, was found. We prove that these equations also have generalized self-similar solutions which describe the temporary transition from one stationary state to another. The study was motivated by problems arising in mathematical modeling of genome size evolution.


title: Triangular lattice neurons may implement an advanced numeral system to precisely encode rat position over large ranges

abstract: We argue by observation of the neural data that neurons in area dMEC of rats, which fire whenever the rat is on any vertex of a regular triangular lattice that tiles 2-d space, may be using an advanced numeral system to reversibly encode rat position. We interpret measured dMEC properties within the framework of a residue number system (RNS), and describe how RNS encoding -- which breaks the non-periodic variable of rat position into a set of narrowly distributed periodic variables -- allows a small set of cells to compactly represent and efficiently update rat position with high resolution over a large range. We show that the uniquely useful properties of RNS encoding still hold when the encoded and encoding quantities are relaxed to be real numbers with built-in uncertainties, and provide a numerical and functional estimate of the range and resolution of rat positions that can be uniquely encoded in dMEC. The use of a compact, `arithmetic-friendly' numeral system to encode a metric variable, as we propose is happening in dMEC, is qualitatively different from all previously identified examples of coding in the brain. We discuss the numerous neurobiological implications and predictions of our hypothesis.


title: Spatially Distributed Stochastic Systems: equation-free and equation-assisted preconditioned computation

abstract: Spatially distributed problems are often approximately modelled in terms of partial differential equations (PDEs) for appropriate coarse-grained quantities (e.g. concentrations). The derivation of accurate such PDEs starting from finer scale, atomistic models, and using suitable averaging, is often a challenging task; approximate PDEs are typically obtained through mathematical closure procedures (e.g. mean-field approximations). In this paper, we show how such approximate macroscopic PDEs can be exploited in constructing preconditioners to accelerate stochastic simulations for spatially distributed particle-based process models. We illustrate how such preconditioning can improve the convergence of equation-free coarse-grained methods based on coarse timesteppers. Our model problem is a stochastic reaction-diffusion model capable of exhibiting Turing instabilities.


title: Model for unidirectional movement of axonemal and cytoplasmic dynein molecules

abstract: A model for the unidirectional movement of dynein is presented based on structural observations and biochemical experimental results available. In this model, the binding affinity of dynein for microtubule is independent of its nucleotide state and the change between strong and weak microtubule-binding is determined naturally by the variation of relative orientation between the stalk and microtubule as the stalk rotates following nucleotide-state transition. Thus the enigmatic communication from the ATP binding site in the globular domain to the far MT-binding site in the tip of the stalk, which is prerequisite in conventional models, is not required. Using the present model, the previous experimental results such as the effect of ATP and ADP bindings on dissociation of dynein from microtubule, the processive movement of single-headed axonemal dyneins at saturating ATP concentration, the load dependence of step size for the processive movement of two-headed cytoplasmic dyneins and the dependence of stall force on ATP concentration can be well explained.


title: Relevance of Quantum Mechanics in Circuit Implementation of Ion channels in Brain Dynamics

abstract: With an increasing amount of experimental evidence pouring in from neurobiological investigations, it is quite appropriate to study viable reductionist models which may explain some of the features of brain activities. It is now quite well known that the Hodgkin-Huxley (HH) Model has been quite successful in explaining the neural phenomena. The idea of circuit equivalents and the membrane voltages corresponding to neurons have been remarkable which is essentially a classical result. In view of some recent results which show that quantum mechanics may be important at suitable length scales inside the brain, the question which becomes quite important is to find out a proper quantum analogue of the HH scheme which will reduce to the well known HH model in a suitable limit. From the ideas of neuro-manifold and the relevance of quantum mechanics at some length scales in the ion channels, we investigate this situation in this paper by taking into consideration the Schr\"odinger equation in an arbitrary manifold with a metric, which is in some sense a special case of the heat kernel equation. The next important approach we have taken in order to bring about it's relevance in brain studies and to make connection with HH models is to find out a plausible circuit equivalents of it. What we do realize is that for a proper quantum mechanical description and it's circuit implementation of the same we need to incorporate the non commutativity inside the circuit model. It has been realized here that the metric is a dynamical entity governing space time and for considering equivalent circuits it plays a very distinct role. We have used the methods of stochastic quantization and have constructed a specific case here and see that HH model inductances gets renormalized in the quantum limit.


title: A Bio-Polymer Transistor: Electrical Amplification by Microtubules

abstract: Microtubules (MTs) are important cytoskeletal structures, engaged in a number of specific cellular activities, including vesicular traffic, cell cyto-architecture and motility, cell division, and information processing within neuronal processes. MTs have also been implicated in higher neuronal functions, including memory, and the emergence of "consciousness". How MTs handle and process electrical information, however, is heretofore unknown. Here we show new electrodynamic properties of MTs. Isolated, taxol-stabilized microtubules behave as bio-molecular transistors capable of amplifying electrical information. Electrical amplification by MTs can lead to the enhancement of dynamic information, and processivity in neurons can be conceptualized as an "ionic-based" transistor, which may impact among other known functions, neuronal computational capabilities.


title: Probabilistic sequence alignments: realistic models with efficient algorithms

abstract: Alignment algorithms usually rely on simplified models of gaps for computational efficiency. Based on an isomorphism between alignments and physical helix-coil models, we show in statistical mechanics that alignments with realistic laws for gaps can be computed with fast algorithms. Improved performances of probabilistic alignments with realistic models of gaps are illustrated. Probabilistic and optimization formulations are compared, with potential implications in many fields and perspectives for computationally efficient extensions to Markov models with realistic long-range interactions.

categories: q-bio.GN cond-mat.dis-nn cond-mat.stat-mech physics.bio-ph physics.comp-ph


title: Simulation in Biology

abstract: An invitation to simulation is made into a realizable challenge: the reader is taken from scratch to games with words, to mathematical models, to statistical analysis, to simulation of evolution. No special requirements are needed: we work with the inbuilt Visual Basic language that is attached to Microsoft(R) Excel and Word. The Tex source must be downloaded to copy the code that must be run.


title: Protein Structure Prediction: The Next Generation

abstract: Over the last 10-15 years a general understanding of the chemical reaction of protein folding has emerged from statistical mechanics. The lessons learned from protein folding kinetics based on energy landscape ideas have benefited protein structure prediction, in particular the development of coarse grained models. We survey results from blind structure prediction. We explore how second generation prediction energy functions can be developed by introducing information from an ensemble of previously simulated structures. This procedure relies on the assumption of a funnelled energy landscape keeping with the principle of minimal frustration. First generation simulated structures provide an improved input for associative memory energy functions in comparison to the experimental protein structures chosen on the basis of sequence alignment.


title: Nested Canalyzing, unate cascade, and polynomial functions

abstract: This paper focuses on the study of certain classes of Boolean functions that have appeared in several different contexts. Nested canalyzing functions have been studied recently in the context of Boolean network models of gene regulatory networks. In the same context, polynomial functions over finite fields have been used to develop network inference methods for gene regulatory networks. Finally, unate cascade functions have been studied in the design of logic circuits and binary decision diagrams. This paper shows that the class of nested canalyzing functions is equal to that of unate cascade functions. Furthermore, it provides a description of nested canalyzing functions as a certain type of Boolean polynomial function. Using the polynomial framework one can show that the class of nested canalyzing functions, or, equivalently, the class of unate cascade functions, forms an algebraic variety which makes their analysis amenable to the use of techniques from algebraic geometry and computational algebra. As a corollary of the functional equivalence derived here, a formula in the literature for the number of unate cascade functions provides such a formula for the number of nested canalyzing functions.


title: A Kinetic Model Describing the Processivity of Myosin-V

abstract: The precise details of how myosin-V coordinates the biochemical reactions and mechanical motions of its two head elements to engineer effective processive molecular motion along actin filaments remain unresolved. We compare a quantitative kinetic model of the myosin-V walk, consisting of five basic states augmented by two further states to allow for futile hydrolysis and detachments, with experimental results for run lengths, velocities, and dwell times and their dependence on bulk nucleotide concentrations and external loads in both directions. The model reveals how myosin-V can use the internal strain in the molecule to synchronise the motion of the head elements. Estimates for the rate constants in the reaction cycle and the internal strain energy are obtained by a computational comparison scheme involving an extensive exploration of the large parameter space. This scheme exploits the fact that we have obtained analytic results for our reaction network, e.g. for the velocity but also the run length, diffusion constant and fraction of backward steps. The agreement with experiment is often reasonable but some open problems are highlighted, in particular the inability of such a general model to reproduce the reported dependence of run length on ADP. The novel way that our approach explores parameter space means that any confirmed discrepancies should give new insights into the reaction network model.


title: Ab initio simulations of Cu binding sites in the N-terminal region of PrP

abstract: The prion protein (PrP) binds Cu2+ ions in the octarepeat domain of the N-terminal tail up to full occupancy at pH=7.4. Recent experiments show that the HGGG octarepeat subdomain is responsible for holding the metal bound in a square planar coordination. By using first principle ab initio molecular dynamics simulations of the Car-Parrinello type, the Cu coordination mode to the binding sites of the PrP octarepeat region is investigated. Simulations are carried out for a number of structured binding sites. Results for the complexes Cu(HGGGW)+(wat), Cu(HGGG) and the 2[Cu(HGGG)] dimer are presented. While the presence of a Trp residue and a H2O molecule does not seem to affect the nature of the Cu coordination, high stability of the bond between Cu and the amide Nitrogens of deprotonated Gly's is confirmed in the case of the Cu(HGGG) system. For the more interesting 2[Cu(HGGG)] dimer a dynamically entangled arrangement of the two monomers, with intertwined N-Cu bonds, emerges. This observation is consistent with the highly packed structure seen in experiments at full Cu occupancy.


title: Large-scale Oscillation of Structure-Related DNA Sequence Features in Human Chromosome 21

abstract: Human chromosome 21 is the only chromosome in human genome that exhibits oscillation of (G+C)-content of cycle length of hundreds kilobases (500 kb near the right telomere). We aim at establishing the existence of similar periodicity in structure-related sequence features in order to relate this (G+C)% oscillation to other biological phenomena. The following quantities are shown to oscillate with the same 500kb periodicity in human chromosome 21: binding energy calculated by two sets of dinucleotide-based thermodynamic parameters, AA/TT and AAA/TTT bi-/tri-nucleotide density, 5'-TA-3' dinucleotide density, and signal for 10/11-base periodicity of AA/TT or AAA/TTT. These intrinsic quantities are related to structural features of the double helix of DNA molecules, such as base-pair binding, untwisting/unwinding, stiffness, and a putative tendency for nucleosome formation.


title: Overlapping Probabilities of Top Ranking Gene Lists, Hypergeometric Distribution, and Stringency of Gene Selection Criterion

abstract: When the same set of genes appear in two top ranking gene lists in two different studies, it is often of interest to estimate the probability for this being a chance event. This overlapping probability is well known to follow the hypergeometric distribution. Usually, the lengths of top-ranking gene lists are assumed to be fixed, by using a pre-set criterion on, e.g., $p$-value for the t-test. We investigate how overlapping probability changes with the gene selection criterion, or simply, with the length of the top-ranking gene lists. It is concluded that overlapping probability is indeed a function of the gene list length, and its statistical significance should be quoted in the context of gene selection criterion.


title: Does Logarithm Transformation of Microarray Data Affect Ranking Order of Differentially Expressed Genes?

abstract: A common practice in microarray analysis is to transform the microarray raw data (light intensity) by a logarithmic transformation, and the justification for this transformation is to make the distribution more symmetric and Gaussian-like. Since this transformation is not universally practiced in all microarray analysis, we examined whether the discrepancy of this treatment of raw data affect the "high level" analysis result. In particular, whether the differentially expressed genes as obtained by $t$-test, regularized t-test, or logistic regression have altered rank orders due to presence or absence of the transformation. We show that as much as 20%--40% of significant genes are "discordant" (significant only in one form of the data and not in both), depending on the test being used and the threshold value for claiming significance. The t-test is more likely to be affected by logarithmic transformation than logistic regression, and regularized $t$-test more affected than t-test. On the other hand, the very top ranking genes (e.g. up to top 20--50 genes, depending on the test) are not affected by the logarithmic transformation.


title: Scaling law for the transient behavior of type-II neuron models

abstract: We study the transient regime of type-II biophysical neuron models and determine the scaling behavior of relaxation times $\tau$ near but below the repetitive firing critical current, $\tau \simeq C (I_c-I)^{-\Delta}$. For both the Hodgkin-Huxley and Morris-Lecar models we find that the critical exponent is independent of the numerical integration time step and that both systems belong to the same universality class, with $\Delta = 1/2$. For appropriately chosen parameters, the FitzHugh-Nagumo model presents the same generic transient behavior, but the critical region is significantly smaller. We propose an experiment that may reveal nontrivial critical exponents in the squid axon.


title: Perception proximale et distale \`{a} l'aide du dispositif de Lenay

abstract: We present in this article experimental results obtained with the "dispositif de Lenay" : for a localization task (distal perception) and an orientation estimation task (proximal perception of the orientation of a cylinder in a plane). In this last experiment, a "virtual" version of the Lenay device was used. Results are here used to illustrate methodological and theoretical proposals for the study of cognitive and sensori-motor processes involved in perception.


title: Heterogeneous animal group models and their group-level alignment dynamics; an equation-free approach

abstract: We study coarse-grained (group-level) alignment dynamics of individual-based animal group models for {\it heterogeneous} populations consisting of informed (on preferred directions) and uninformed individuals. The orientation of each individual is characterized by an angle, whose dynamics are nonlinearly coupled with those of all the other individuals, with an explicit dependence on the difference between the individual's orientation and the instantaneous average direction. Choosing convenient coarse-grained variables (suggested by uncertainty quantification methods) that account for rapidly developing correlations during initial transients, we perform efficient computations of coarse-grained steady states and their bifurcation analysis. We circumvent the derivation of coarse-grained governing equations, following an equation-free computational approach.


title: Mean Field Model of Genetic Regulatory Networks

abstract: In this paper, we propose a mean-field model which attempts to bridge the gap between random Boolean networks and more realistic stochastic modeling of genetic regulatory networks. The main idea of the model is to replace all regulatory interactions to any one gene with an average or effective interaction, which takes into account the repression and activation mechanisms. We find that depending on the set of regulatory parameters, the model exhibits rich nonlinear dynamics. The model also provides quantitative support to the earlier qualitative results obtained for random Boolean networks.


title: Subgraph Ensembles and Motif Discovery Using a New Heuristic for Graph Isomorphism

abstract: A new heuristic based on vertex invariants is developed to rapidly distinguish non-isomorphic graphs to a desired level of accuracy. The method is applied to sample subgraphs from an E.coli protein interaction network, and as a probe for discovery of extended motifs. The network's structure is described using statistical properties of its $N$-node subgraphs for $N\leq 14$. The Zipf plots for subgraph occurrences are robust power laws that do not change when rewiring the network while fixing the degree sequence -- although the specific subgraphs may exchange ranks. However the exponent depends on $N$. The study of larger subgraphs highlights some striking patterns for various $N$. Motifs, or connected pieces that are over-abundant in the ensemble of subgraphs, have more edges, for a given number of nodes, than antimotifs and generally display a bipartite structure or tend towards a complete graph. In contrast, antimotifs, which are under-abundant connected pieces, are mostly trees or contain at most a single, small loop. The extension to directed graphs is straightforward.


title: Bistable equilibrium points of mercury body burden

abstract: In the last century mercury levels in the global environment have tripled as a result of increased pollution from industrial, occupational, medicinal and domestic uses \cite{BaMe03}. Glutathione is known to be the main agent responsible for the excretion of mercury (we refer to \cite{Thim05}, \cite{ZalBar99} and \cite{Lyn02}). It has also been shown that mercury inhibits glutathione synthetase (an enzyme acting in the synthesization of Glutathione), therefore leading to decreased glutathione levels (we refer to \cite{Thim05}, \cite{GeGe05}, \cite{GeGe06} and \cite{RDeth04}). Mercury also interferes with the production of heme in the porphyrin pathway \cite{WoMaEc93}. Heme is needed for biological energy production and ability to detox organic toxins via the P450 enzymes \cite{Boy06}. The purpose of this paper is to show that body's response to mercury exposure is hysteretic, i.e. when this feedback of mercury on its main detoxifying agents is strong enough then mercury body burden has two points of equilibrium: one with normal abilities to detoxify and low levels of mercury and one with inhibited abilities to detoxify and high levels of mercury. Furthermore, a small increase of body's mercury burden may not be sufficient to trigger observable neurotoxic effects but it may be sufficient to act as a switch leading to an accumulation of mercury in the body through environmental exposure until its toxicity becomes manifest.


title: Deterministic and stochastic regimes of asexual evolution on rugged fitness landscapes

abstract: We study the adaptation dynamics of an initially maladapted asexual population with genotypes represented by binary sequences of length $L$. The population evolves in a maximally rugged fitness landscape with a large number of local optima. We find that whether the evolutionary trajectory is deterministic or stochastic depends on the effective mutational distance $d_{\mathrm{eff}}$ upto which the population can spread in genotype space. For $d_{\mathrm{eff}}=L$, the deterministic quasispecies theory operates while for $d_{\mathrm{eff}} < 1$, the evolution is completely stochastic. Between these two limiting cases, the dynamics are described by a local quasispecies theory below a crossover time $T_{\times}$ while above $T_{\times}$, the population gets trapped at a local fitness peak and manages to find a better peak either via stochastic tunneling or double mutations. In the stochastic regime $d_\mathrm{eff} < 1$, we identify two subregimes associated with clonal interference and uphill adaptive walks, respectively. We argue that our findings are relevant to the interepretation of evolution experiments with microbial populations.


title: Parametrization of the Hybrid Potential for Pairs of Neutral Atoms

abstract: The hybrid form is a combination of the Rydberg potential and the London inverse-sixth-power energy. It is accurate at all relevant distance scales and simple enough for use in all-atom simulations of biomolecules. One may compute the parameters of the hybrid potential for the ground state of a pair of neutral atoms from their internuclear separation, the depth and curvature of their potential at its minimum, and from their van der Waals coefficient of dispersion.


title: Conservation rules, their breakdown, and optimality in Caenorhabditis sinusoidal locomotion

abstract: Undulatory locomotion is common to nematodes as well as to limbless vertebrates, but its control is not understood in spite of the identification of hundred of genes involved in Caenorhabditis elegans locomotion. To reveal the mechanisms of nematode undulatory locomotion, we quantitatively analyzed the movement of C. elegans with genetic perturbations to neurons, muscles, and skeleton (cuticle). We also compared locomotion of different Caenorhabditis species. We constructed a theoretical model that combines mechanics and biophysics, and that is constrained by the observations of propulsion and muscular velocities, as well as wavelength and amplitude of undulations. We find that normalized wavelength is a conserved quantity among wild-type C. elegans individuals, across mutants, and across different species. The velocity of forward propulsion scales linearly with the velocity of the muscular wave and the corresponding slope is also a conserved quantity and almost optimal; the exceptions are in some mutants affecting cuticle structure. In theoretical terms, the optimality of the slope is equivalent to the exact balance between muscular and visco-elastic body reaction bending moments. We find that the amplitude and frequency of undulations are inversely correlated and provide a theoretical explanation for this fact. These experimental results are valid both for young adults and for all larval stages of wild type C. elegans. In particular, during development, the amplitude scales linearly with the wavelength, consistent with our theory. We also investigated the influence of substrate firmness on motion parameters, and found that it does not affect the above invariants. In general, our biomechanical model can explain the observed robustness of the mechanisms controlling nematode undulatory locomotion.


title: Does the dynamics of sine-Gordon solitons predict active regions of DNA?

abstract: In this work we analyze the possibility that soliton dynamics in a simple nonlinear model allows functionally relevant predictions of the behaviour of DNA. This suggestion was first put forward by Salerno [Phys. Rev. A, vol. 44, p. 5292 (1991)] by showing results indicating that sine-Gordon kinks were set in motion at certain regions of a DNA sequence that include promoters. We revisit that system and show that the observed behaviour has nothing to do with promoters; on the contrary, it originates from the bases at the boundary, which are not part of the studied genome. We explain this phenomenology in terms of an effective potential for the kink center. This is further extended to disprove recent claims that the dynamics of kinks [Lenholm and H\"ornquist, Physica D, vol. 177, p. 233 (2003)] or breathers [Bashford, J. Biol. Phys., vol. 32, p. 27 (2006)] has functional significance. We conclude that no such information can be extracted from this simple nonlinear model or its associated effective potential.

categories: q-bio.GN math-ph math.MP nlin.PS q-bio.QM


title: Is the intrinsic disorder of proteins the cause of the scale-free architecture of protein-protein interaction networks?

abstract: In protein-protein interaction networks certain topological properties appear to be recurrent: networks maps are considered scale-free. It is possible that this topology is reflected in the protein structure. In this paper we investigate the role of protein disorder in the network topology. We find that the disorder of a protein (or of its neighbors) is independent of its number of protein-protein interactions. This result suggests that protein disorder does not play a role in the scale-free architecture of protein networks.


title: Diagnosis of rotator cuff lesions by FT-Raman spectroscopy: a biochemical study

abstract: The biochemical changes on normal and degenerated tissues of rotator cuff supraspinatus tendons were probed by FT-Raman spectroscopy. The Raman spectra showed differences on the spectral regions of cysteine, amino acids, nucleic acids, carbohydrates, and lipids. These spectral differences were assigned to pathological biochemical alterations due to the degenerative process of the tendon. Principal Components Analysis was performed on the spectral data and enabled the correct classification of the spectra as normal (grade 1) and degenerated (grades 2 and 3). These findings indicate that Raman spectroscopy could be a very promising tool for the rotator cuff supraspinatus tendon diagnosis and for quantification of their degenerative degree.


title: Simple solvation potential for coarse-grained models of proteins

abstract: We formulate a simple solvation potential based on a coarsed-grain representation of amino acids with two spheres modeling the $C_\alpha$ atom and an effective side-chain centroid. The potential relies on a new method for estimating the buried area of residues, based on counting the effective number of burying neighbours in a suitable way. This latter quantity shows a good correlation with the buried area of residues computed from all atom crystallographic structures. We check the discriminatory power of the solvation potential alone to identify the native fold of a protein from a set of decoys and show the potential to be considerably selective.


title: Probabilistic Regulatory Networks: Modeling Genetic Networks

abstract: We describe here the new concept of $\epsilon$-Homomorphisms of Probabilistic Regulatory Gene Networks(PRN). The $\epsilon$-homomorphisms are special mappings between two probabilistic networks, that consider the algebraic action of the iteration of functions and the probabilistic dynamic of the two networks. It is proved here that the class of PRN, together with the homomorphisms, form a category with products and coproducts. Projections are special homomorphisms, induced by invariant subnetworks. Here, it is proved that an $\epsilon$-homomorphism for 0 <$\epsilon$< 1 produces simultaneous Markov Chains in both networks, that permit to introduce the concepts of $\epsilon$-isomorphism of Markov Chains, and similar networks.


title: Time Scales in Evolutionary Dynamics

abstract: Evolutionary game theory has traditionally assumed that all individuals in a population interact with each other between reproduction events. We show that eliminating this restriction by explicitly considering the time scales of interaction and selection leads to dramatic changes in the outcome of evolution. Examples include the selection of the inefficient strategy in the Harmony and Stag-Hunt games, and the disappearance of the coexistence state in the Snowdrift game. Our results hold for any population size and in the presence of a background of fitness.

categories: q-bio.PE math.DS nlin.AO physics.soc-ph q-bio.QM


title: Why are diploid genomes widespread and dominant mutations rare?

abstract: We have used the sexual Penna ageing model to show that the relation between dominance and recessiveness could be a force which optimizes the genome size. While the possibility of complementation of the damaged allele by its functional counterparts (recessiveness) leads to the redundancy of genetic information, the dominant effect of defective genes tends to diminish the number of alleles fulfilling the same function. Playing with the fraction of dominant loci in the genome it is possible to obtain the condition where the diploid state of the genome is optimal. If the status of each bit position as dominant or recessive mutations is changed for each individual randomly and rarely, then after a long time a stationary equilibrium of many recessive and few dominant loci is established in the sexual Penna model. This effect vanishes if the same changing distribution of dominant loci applies to all individuals.


title: A multiscale mathematical model of cancer, and its use in analyzing irradiation therapies

abstract: Background: Radiotherapy outcomes are usually predicted using the Linear Quadratic model. However, this model does not integrate complex features of tumor growth, in particular cell cycle regulation. Methods: In this paper, we propose a multiscale model of cancer growth based on the genetic and molecular features of the evolution of colorectal cancer. The model includes key genes, cellular kinetics, tissue dynamics, macroscopic tumor evolution and radiosensitivity dependence on the cell cycle phase. We investigate the role of gene-dependent cell cycle regulation in the response of tumors to therapeutic irradiation protocols. Results: Simulation results emphasize the importance of tumor tissue features and the need to consider regulating factors such as hypoxia, as well as tumor geometry and tissue dynamics, in predicting and improving radiotherapeutic efficacy. Conclusion: This model provides insight into the coupling of complex biological processes, which leads to a better understanding of oncogenesis. This will hopefully lead to improved irradiation therapy.


title: Evolution of Protein Interaction Networks by Whole Genome Duplication and Domain Shuffling

abstract: Successive whole genome duplications have recently been firmly established in all major eukaryote kingdoms. It is not clear, however, how such dramatic evolutionary process has contributed to shape the large scale topology of protein-protein interaction (PPI) networks. We propose and analytically solve a generic model of PPI network evolution under successive whole genome duplications. This demonstrates that the observed scale-free degree distributions and conserved multi-protein complexes may have concomitantly arised from i) intrinsic exponential dynamics of PPI network evolution and ii) asymmetric divergence of gene duplicates. This requirement of asymmetric divergence is in fact "spontaneously" fulfilled at the level of protein-binding domains. In addition, domain shuffling of multi-domain proteins is shown to provide a powerful combinatorial source of PPI network innovation, while preserving essential structures of the underlying single-domain interaction network. Finally, large scale features of PPI networks reflecting the "combinatorial logic" behind direct and indirect protein interactions are well reproduced numerically with only two adjusted parameters of clear biological significance.


title: The propagation of a cultural or biological trait by neutral genetic drift in a subdivided population

abstract: We study fixation probabilities and times as a consequence of neutral genetic drift in subdivided populations, motivated by a model of the cultural evolutionary process of language change that is described by the same mathematics as the biological process. We focus on the growth of fixation times with the number of subpopulations, and variation of fixation probabilities and times with initial distributions of mutants. A general formula for the fixation probability for arbitrary initial condition is derived by extending a duality relation between forwards- and backwards-time properties of the model from a panmictic to a subdivided population. From this we obtain new formulae, formally exact in the limit of extremely weak migration, for the mean fixation time from an arbitrary initial condition for Wright's island model, presenting two cases as examples. For more general models of population subdivision, formulae are introduced for an arbitrary number of mutants that are randomly located, and a single mutant whose position is known. These formulae contain parameters that typically have to be obtained numerically, a procedure we follow for two contrasting clustered models. These data suggest that variation of fixation time with the initial condition is slight, but depends strongly on the nature of subdivision. In particular, we demonstrate conditions under which the fixation time remains finite even in the limit of an infinite number of demes. In many cases - except this last where fixation in a finite time is seen - the time to fixation is shown to be in precise agreement with predictions from formulae for the asymptotic effective population size.


title: Use of the Metropolis algorithm to simulate the dynamics of protein chains

abstract: The Metropolis implementation of the Monte Carlo algorithm has been developed to study the equilibrium thermodynamics of many-body systems. Choosing small trial moves, the trajectories obtained applying this algorithm agree with those obtained by Langevin's dynamics. Applying this procedure to a simplified protein model, it is possible to show that setting a threshold of 1 degree on the movement of the dihedrals of the protein backbone in a single Monte Carlo step, the mean quantities associated with the off-equilibrium dynamics (e.g., energy, RMSD, etc.) are well reproduced, while the good description of higher moments requires smaller moves. An important result is that the time duration of a Monte Carlo step depends linearly on the temperature, something which should be accounted for when doing simulations at different temperatures.


title: Randomization and Feedback Properties of Directed Graphs Inspired by Gene Networks

abstract: Having in mind the large-scale analysis of gene regulatory networks, we review a graph decimation algorithm, called "leaf-removal", which can be used to evaluate the feedback in a random graph ensemble. In doing this, we consider the possibility of analyzing networks where the diagonal of the adjacency matrix is structured, that is, has a fixed number of nonzero entries. We test these ideas on a network model with fixed degree, using both numerical and analytical calculations. Our results are the following. First, the leaf-removal behavior for large system size enables to distinguish between different regimes of feedback. We show their relations and the connection with the onset of complexity in the graph. Second, the influence of the diagonal structure on this behavior can be relevant.