Search

Motivated by the operation of myogenic (self-oscillatory) insect flight muscle, we study a model consisting of a large number of identical oscillatory contractile elements joined in a chain, whose end is attached to a damped mass-spring oscillator. When the inertial load is small, the serial coupling favors an antisynchronous state in which the extension of one oscillator is compensated by the contraction of another, in order to preserve the total length. However, a sufficiently massive load can sychronize the oscillators and can even induce oscillation in situations where isolated elements would be stable. The system has a complex phase diagram displaying quiescent, synchronous and antisynchrononous phases, as well as an unsual asynchronous phase in which the total length of the chain oscillates at a different frequency from the individual active elements., Comment: 5 pages, 4 figures, To appear in Phys. Rev. Lett

Most organisms grow according to simple laws, which in principle can be derived from energy conservation and scaling arguments, critically dependent on the relation between the metabolic rate B of energy flow and the organism mass m. Although this relation is generally recognized to be of the form B(m) = mp, the specific value of the exponent p is the object of an ongoing debate, with many mechanisms being postulated to support different predictions. We propose that multicellular tumor spheroids provide an ideal experimental model system for testing these allometric growth theories, especially under controlled conditions of malnourishment and applied mechanical stress.

Motor enzymes catalyze chemical reactions, like the hydrolysis of ATP, and in the process they also perform work. Recent studies indicate that motor enzymes perform work with specific intermediate steps in their catalyzed reactions, challenging the classic view (in Brownian motor models) that work can only be performed within biochemical states. An alternative class of models (chemical motor models) has emerged in which motors perform work with biochemical transitions, but many of these models lack a solid physicochemical foundation. In this paper, I develop a self consistent framework for chemical motor models. This novel framework accommodates multiple pathways for free energy transfer, predicts rich behaviors from the simplest multi motor systems, and provides important new insights into muscle and motor function., Comment: 26 pages, 6 figures

We present a theoretical model for retinocollicular map development, which can account for intriguing behaviors observed in gain-of-function experiments in knock-in mice by Brown et al., including bifurcation in heterozygous Isl2/EphA3 knock-ins. The model is based on known chemical labels, axonal repulsion/competition, stochasticity and uses Markov chain description. Our model suggests that the map in heterozygotes is single-valued in the temporal region of retina due to reduced gradient of ephrin in the corresponding region of SC. The remaining map is double-valued since the gradient of ephrin is high there. We predict therefore that if gradient of ephrin is reduced by a genetic manipulation, the single-valued region of the map should occupy a larger portion of temporal retina, i.e. the point of transition between single- and doulble-valued maps should move to a more nasal position in Isl2-EphA3 heterozygotes. We also discuss the importance of inhomogeneous EphA gradient and mapping in Isl2/EphB knock-ins., Comment: 19 pages, 11 color figures

We develop a propagating front analysis, in terms of a local probability of zipping, for the helicase activity of opening up a double stranded DNA (dsDNA). In a fixed-distance ensemble (conjugate to the fixed-force ensemble) the front separates the zipped and unzipped phases of a dsDNA and a drive acts locally around the front. Bounds from variational analysis and numerical estimates for the speed of a helicase are obtained. Different types of helicase behaviours can be distinguished by the nature of the drive., Comment: 5 pages, 5 eps figures; replaced by the published version

Survival in groups of termites is known to be socially facilitated since individual longevity is increased together with the size of the group. Here we report on experimental evidence that survival as a function of group size in {\textit Cornitermes cumulans} (Kollar) does not increase asysmptotically but it peaks in a group size range at which individual longevity is maximal. Strinkly, this same effect in other social insects was noted more than fifty years ago, but it was regarded as ``aberrant'' since it has been generally assumed that individual survival should increase with the size of the group. We also provide evidence that individual activity --measured as mobility-- is maximal at approximately the same group-size range as the observed survival

A mechanism of memories, especially biological memories, is studied in terms of quantum fluids. Two-dimensional flows in central potentials $V_a(\rho)=-a^2g_a\rho^{2(a-1)}$ ($a\not=0$ and $\rho=\sqrt{x^2+y^2}$) have zero-energy eigenstates that degenerate infinitely for all $a$. It is shown that stable standing waves constructed from the zero-energy flows are confined in various types of polygons which can be the minimum units of memory systems. Vortex patterns awoken in the units by stimuli correspond to the memories of the stimuli. This memory system is not a system for preserving memories as usual but that for awaking memories. The system has interesting properties; (i) the absolute economy as for the energy consumption, (ii) the infinite variety for a huge number of memories, (iii) the perfect recovery of the system from any disturbances by stimuli, and (iv) the large flexibility in the construction of the system. A process for thinking is also proposed in terms of this memory system., Comment: 8 pages

We present a study on filamentous actin solutions containing heavy meromyosin subfragments of myosin II motor molecules. We focus on the viscoelastic phase behavior and internal dynamics of such networks during ATP depletion. Upon simultaneously using micro-rheology and fluorescence microscopy as complementary experimental tools, we find a sol-gel transition accompanied by a sudden onset of directed filament motion. We interpret the sol-gel transition in terms of myosin II enzymology, and suggest a "zipping" mechanism to explain the filament motion in the vicinity of the sol-gel transition., Comment: 4 pages, 3 figures

The SIRM (Stochastic Inclined Rods Model) proposed by H. Matsuura and M. Nakano can explain the muscle's motion perfectly, but the intermolecular potential between myosin head and G-actin is too simple and only repulsive potential is considered. In this paper we study the SIRM with different complex potential and discuss the effect of the spring on the system. The calculation results show that the spring, the effective radius of the G-actin and the intermolecular potential play key roles in the motion. The sliding speed is about $4.7\times10^{-6}m/s$ calculated from the model which well agrees with the experimental data., Comment: 9 pages, 6 figures

Cu,Zn superoxide dismutase (SOD1) has been implicated in the familial form of the neurodegenerative disease Amyotrophic Lateral Sclerosis (ALS). It has been suggested that mutant mediated SOD1 misfolding/aggregation is an integral part of the pathology of ALS. We study the folding thermodynamics and kinetics of SOD1 using a hybrid molecular dynamics approach. We reproduce the experimentally observed SOD1 folding thermodynamics and find that the residues which contribute the most to SOD1 thermal stability are also crucial for apparent two-state folding kinetics. Surprisingly, we find that these residues are located on the surface of the protein and not in the hydrophobic core. Mutations in some of the identified residues are found in patients with the disease. We argue that the identified residues may play an important role in aggregation. To further characterize the folding of SOD1, we study the role of cysteine residues in folding and find that non-native disulfide bond formation may significantly alter SOD1 folding dynamics and aggregation propensity., Comment: 16 pages, 5 figures

Based on the Stochastic Inclined Rods Model (SIRM) proposed by H. Matsuura and M. Nakano, we study the microscopic motion of actin myosin system including the motion of the G-actin. Our model is composed of an inclined spring (rod), a myosin head, a myosin filament and G-actins. We discuss the stochastic resonance between the myosin head and random noise. The results of calculation show that the model can convert the random motion to one directional motion, and the myosin head works as a resonator of random noise which absorbs the energy through the stochastic resonance. The intermolecular potential between the myosin head and G-actin and the inclined rod play a key role for the muscle's motion. The energy consumed by the motor is directly supplied from the surroundings (i.e., the thermal motions of water molecules)., Comment: 11 pages, 8 figures

The logistic differential equation is used to analyze cancer cell population, in the presence of a correlated Gaussian white noise. We study the steady state properties of tumor cell growth and discuss the effects of the correlated noise. It is found that the degree of correlation of the noise can cause tumor cell extinction., Comment: 3 pages, 4 figures

We study the steady state properties of an insect (spruce budworm) outbreak model in the presence of Gaussian white noise. Based on the corresponding Fokker-Planck equation the steady state solution of the probability distribution function and its extrema have been investigated. It was found that fluctuations of the insect birth rate reduces the population of the insects while fluctuations of predation rate and the noise correlation can prevent the population of the insects from going into extinction. Noise in the model can induce a phase transition., Comment: 8 pages, 4 figures

We study the steady state properties of a genotype selection model in presence of correlated Gaussian white noise. The effect of the noise on the genotype selection model is discussed. It is found that correlated noise can break the balance of gene selection and induce the phase transition which can makes us select one type gene haploid from a gene group., Comment: 8 pages, 4 figures

A new method of absolute heat capacity measurement was developed in the Department of Biological Systems Physics at the E. Andronikashvili Institute of Physics, Georgian Academy of Sciences. This is the method of absolute heat capacity measurement of the differential calorimetric ampoule sample and reference cell temperature difference registration while alternate transferring of heat pulses to the cells occurs. The new method makes possible to determine both the heat capacities of the sample and reference cells and the difference of thereof. The method of cell temperature difference compensation is used in order to decrease in the time of equilibrium establishment. The proposed method of absolute heat capacity measurement was implemented in the model of the virtual differential pulse calorimeter. The model allows to perform calculations of the design of the differential ampoule - the main part of any calorimeter., Comment: 6 pages, 4 figures

The morphology and time-evolution of tumors are expected to depend heavily on the detailed balance of the overall physics of the cell assembly (e.g., the kinetic pressure, the cell-cell interaction, and the external trapping by the tissue) and the biological processes of mitosis, necrosis, etc. Here, for the first time, we include such a detailed balance in a theoretical model by exploiting the {\it ab initio} mathematical framework of atomic Bose-Einstein Condensation (BEC) with feeding and dissipation, to study tumor evolution dynamics. We show that, with a characteristic length scaling, the Gross-Pitaevskii (GP) equation, which describes the many-body atomic BEC characteristics, indeed explains the detailed features of a prevascular tumor culture data. The agreement suggests the prevascular that the prevascular carcinoma may be a natural analog to BEC and predicts an intercellular wave connecting the cells., Comment: 9 pages, 4 figures

Extensive Monte Carlo folding simulations for four proteins of various structural classes are carried out, using a single atomistic potential. In all cases, collapse occurs at a very early stage, and proteins fold into their native-like conformations at appropriate temperatures. The results demonstrate that the folding mechanism is controlled not only by thermodynamic factors but also by kinetic factors: The way a protein folds into its native structure, is also determined by the convergence point of early folding trajectories, which cannot be obtained by the free energy surface., Comment: 11 pages, 4 figures

We consider a model for strand separation in negatively supercoiled dna, and study the effect of supercoiling on the localization of denaturation, using notions from statistical mechanics. The model is based on on mean field Ising system, coupled with a random field one dimensional ISing model. The effect of the word in AT and Gc associated with the DNA on localization of denaturation is studied under the hypothesis of small perimeter., Comment: 20 pages

A computer model was developed for estimation of the kinetics of microbial inactivation by pulsed electric field. The model is based on the electroporation theory of individual membrane damage, where spherical cell geometry and distribution of cell sizes are assumed. The variation of microbial cell sizes was assumed to follow a statistical probability distribution of the Gaussian type. Surviving kinetics was approximated by Weibull equation. The dependencies of two Weibull parameters (shape \textit{n} and time $\tau $, respectively) versus electric field intensity E and width of cell diameters distribution was studied., Comment: 4 pages, 4 figures, RevTex4

We present a theory to investigate electro-kinetic behavior, namely, electrorotation and dielectrophoresis under alternating current (AC) applied fields for a pair of touching inhomogeneous colloidal particles and biological cells. These inhomogeneous particles are treated as graded ones with physically motivated model dielectric and conductivity profiles. The mutual polarization interaction between the particles yields a change in their respective dipole moments, and hence in the AC electrokinetic spectra. The multipolar interactions between polarized particles are accurately captured by the multiple images method. In the point-dipole limit, our theory reproduces the known results. We find that the multipolar interactions as well as the spatial fluctuations inside the particles can affect the AC electrokinetic spectra significantly., Comment: Revised version with minor changes: References added and discussion extended

We examine the dependence of the physical quantities of the rotatory molecular motor, such as the rotation velocity and the proton translocation rate, on the chemical reaction rate using the model based only on diffusion process. A peculiar behavior of proton translocation is found and the energy transduction efficiency of the motor protein is enhanced by this behavior. We give a natural explanation that this behavior is universal when certain inequalities between chemical reaction rates hold. That may give a clue to examine whether the motion of the molecular motor is dominated by diffusion process or not., Comment: 12 pages, 8 figures

Understanding of the evolutionary origins of protein structures represents a key component of the understanding of molecular evolution as a whole. Here we seek to elucidate how the features of an underlying protein structural "space" might impact protein structural evolution. We approach this question using lattice polymers as a completely characterized model of this space. We develop a measure of structural comparison of lattice structures that is analgous to the one used to understand structural similarities between real proteins. We use this measure of structural relatedness to create a graph of lattice structures and compare this graph (in which nodes are lattice structures and edges are defined using structural similarity) to the graph obtained for real protein structures. We find that the graph obtained from all compact lattice structures exhibits a distribution of structural neighbors per node consistent with a random graph. We also find that subgraphs of 3500 nodes chosen either at random or according to physical constraints also represent random graphs. We develop a divergent evolution model based on the lattice space which produces graphs that, within certain parameter regimes, recapitulate the scale-free behavior observed in similar graphs of real protein structures., Comment: 27 pages, 7 figures

We introduce a phase field approach for diffusion inside and outside a closed cell with damping and with source terms at the interface. The method is compared to exact solutions (where possible) and the more traditional finite element method. It is shown to be very accurate, easy to implement and computationally inexpensive. We apply our method to a recently introduced model for chemotaxis by Rappel et al. [Biophys. J. 83, 1361 (2002)]., Comment: 4 pages, 4 figures

Experimental observations suggest that proteins follow different pathways under different environmental conditions. We perform molecular dynamics simulations of a model of the SH3 domain over a broad range of temperatures, and identify distinct pathways in the folding transition. We determine the kinetic partition temperature --the temperature for which the SH3 domain undergoes a rapid folding transition with minimal kinetic barriers-- and observe that below this temperature the model protein may undergo a folding transition via multiple folding pathways. The folding kinetics is characterized by slow and fast pathways and the presence of only one or two intermediates. Our findings suggest the hypothesis that the SH3 domain, a protein for which only two-state folding kinetics was observed in previous experiments, may exhibit intermediates states under extreme experimental conditions, such as very low temperatures. A very recent report (Viguera et al., Proc. Natl. Acad. Sci. USA, 100:5730--5735, 2003) of an intermediate in the folding transition of the Bergerac mutant of the alpha-spectrin SH3 domain protein supports this hypothesis., Comment: 16 pages, 4 figures To be published in the "Journal of Molecular Biology"

An In Silico model to relate the properties of proteins to the structure, sequence, function and evolutionary history of proteins is shown. The derived ideal sequences for amino acid residues in proteins can then be considered as attractors for structures of actual protein families and their functions. Hopfield networks1 can then act as GCMs (General Content Addressable Memories) when they are trained to recall unique pre-determined states when presented with information associated with that state. The tertiary structure of proteins is determined by structurally enforced interactions between residues and they help to elucidate common motifs exhibited in proteins. The acceptance (or rejection) of substitution on the protein backbone conforms to a nonlinear response and can lead to emergent behavior., Comment: 19 pages, 14 figures

The hair cells of the vertebrate inner ear convert mechanical stimuli to electrical signals. Two adaptation mechanisms are known to modify the ionic current flowing through the transduction channels of the hair bundles: a rapid process involves calcium ions binding to the channels; and a slower adaptation is associated with the movement of myosin motors. We present a mathematical model of the hair cell which demonstrates that the combination of these two mechanisms can produce `self-tuned critical oscillations', i.e. maintain the hair bundle at the threshold of an oscillatory instability. The characteristic frequency depends on the geometry of the bundle and on the calcium dynamics, but is independent of channel kinetics. Poised on the verge of vibrating, the hair bundle acts as an active amplifier. However, if the hair cell is sufficiently perturbed, other dynamical regimes can occur. These include slow relaxation oscillations which resemble the hair bundle motion observed in some experimental preparations., Comment: 13 pages, 6 figures,REVTeX 4, To appear in Biophysical Journal

The ensemble of RNA secondary structures of uniform sequences is studied analytically. We calculate the partition function for very long sequences and discuss how the cross-over length, beyond which asymptotic scaling laws apply, depends on thermodynamic parameters. For realistic choices of parameters this length can be much longer than natural RNA molecules. This has to be taken into account when applying asymptotic theory to interpret experiments or numerical results., Comment: 10 pages, 13 figures, published in Phys. Rev. E

The influx of calcium ions into the dendritic spines through the N-metyl-D-aspartate (NMDA) channels is believed to be the primary trigger for various forms of synaptic plasticity. In this paper, the authors calculate analytically the mean values of the calcium transients elicited by a spiking neuron undergoing a simple model of ionic currents and back-propagating action potentials. The relative variability of these transients, due to the stochastic nature of synaptic transmission, is further considered using a simple Markov model of NMDA receptos. One finds that both the mean value and the variability depend on the timing between pre- and postsynaptic action-potentials. These results could have implications on the expected form of synaptic-plasticity curve and can form a basis for a unified theory of spike time-dependent, and rate based plasticity., Comment: 14 pages, 10 figures. A few changes in section IV and addition of a new figure

We investigate the stability of negative image equilibria in mean synaptic weight dynamics governed by spike-timing dependent plasticity (STDP). The neural architecture of the model is based on the electrosensory lateral line lobe (ELL) of mormyrid electric fish, which forms a negative image of the reafferent signal from the fish's own electric discharge to optimize detection of external electric fields. We derive a necessary and sufficient condition for stability, for arbitrary postsynaptic potential functions and arbitrary learning rules. We then apply the general result to several examples of biological interest., Comment: 13 pages, revtex4; uses packages: graphicx, subfigure; 9 figures, 16 subfigures

Qualitative attributes of the region between order and disorder are examined to explore models of genetic and protein networks. Results show how the connectivity of vertices and the strength of their connections are related and how their stability is related to their geometry. It has been possible to relate the interaction energies and chemical potentials of general lattices with the coordination numbers of such lattice models. The results seem to agree with some of those obtained by means of other methods, e.g., those of Barabasi, Strognatz, Dorogovtsev, and Kauffman (BSDK). This can yield new perspectives to the cause, treatment and remedies of disease other than the present mode of drug discovery prevalent: the docking of a single ligand on to a target molecule. In order to utilize such results more efficiently, the pathway of biological processes need be elucidated. A method is available for determining biochemical-reaction or metabolic pathways through its systematic synthesis. It is based on a rigorous graph-theoretic method for identifying pathways of catalytic reactions. It synthesizes networks of metabolic pathways using a highly exacting combinatorial method. It generates not only all feasible, independent reaction networks but also those combinations of independent pathways. This method can determine the mechanisms of complex chemical reactions and is applicable to biochemical reactions. It is important to combine this result with the mechanisms existent in gene regulatory networks. Training Genetic Regulatory Networks for feasible biochemical reaction networks or pathways and incorporating such knowledge into DNA would be a superb technique for vanquishing complex diseases., Comment: 13 pages, 7 figures

It has been argued that the sum of usage probabilities for codons, belonging to quartets, that have as third nucleotide C or A, is independent of the biological species for vertebrates. The comparison between the theoretical correlation matrix derived from these sum rules and the experimentally computed matrix for 26 species shows a satisfactory agreement. The Shannon entropy, weakly depending on the biological species, gives further support. Suppression of codons containing the dinucleotides CG or AU is put in evidence.

Localization properties of residue fluctuations in globular proteins are studied theoretically by using the Gaussian network model. Participation ratio for each residue fluctuation mode is calculated. It is found that the relationship between participation ratio and frequency is similar for all globular proteins, indicating a universal behavior in spite of their different size, shape, and architecture., Comment: 4 pages, 3 figures. To appear in Phys. Rev. E

The immunological synapse is a patterned collection of different types of receptors and ligands that forms in the intercellular junction between T Cells and antigen presenting cells (APCs) during recognition. The synapse is implicated in information transfer between cells, and is characterized by different spatial patterns of receptors at different stages in the life cycle of T cells. We obtain a minimalist model that captures this experimentally observed phenomenology. A functional RG analysis provides further insights., Comment: 6 pages, 3 figures, submitted for publication

We address theoretically aggregation of DNA segments by multivalent polyamines such as spermine and spermidine. In experiments, the aggregation occurs above a certain threshold concentration of multivalent ions. We demonstrate that the dependence of this threshold on the concentration of DNA has a simple form. When the DNA concentration c_DNA is smaller than the monovalent salt concentration, the threshold multivalent ion concentration depends linearly on c_DNA, having the form alpha c_DNA + beta. The coefficients alpha and beta are related to the density profile of multivalent counterions around isolated DNA chains, at the onset of their aggregation. This analysis agrees extremely well with recent detailed measurements on DNA aggregation in the presence of spermine. From the fit to the experimental data, the number of condensed multivalent counterions per DNA chain can be deduced. A few other conclusions can then be reached: i) the number of condensed spermine ions at the onset of aggregation decreases with the addition of monovalent salt; ii) the Poisson-Boltzmann theory over-estimates the number of condensed multivalent ions at high monovalent salt concentrations; iii) our analysis of the data indicates that the DNA charge is not over-compensated by spermine at the onset of aggregation., Comment: 12 pages, 8 figures. Biophysical Journal 2003, in press

We investigate a map-based model of paced cardiac muscle in the presence of closed-loop feedback control. The model relates the duration of an action potential to the preceding diastolic interval as well as the preceding action potential duration and thus has some amount of `memory.' We find that the domain of control depends on this memory, independently of the specific functional form of the map. The memory-dependent domain of control can encompass large feedback gains, thus providing the first possible explanation of the recent experimental results of Hall et al. [Phys. Rev. Lett, 88, 198102 (2002)] on controlling alternans in small pieces of rapidly-paced cardiac muscle., Comment: 10 pages, 2 figures, submitted to Phys. Rev. Lett

We investigate the isometric transient response of muscle using a quantitative stochastic model of the actomyosin cycle based on the swinging lever-arm hypothesis. We first consider a single pair of filaments, and show that when values of parameters such as the lever-arm displacement and the crossbridge elasticity are chosen to provide effective energy transduction, the T2 curve (the tension recovered immediately after a step displacement) displays a region of negative slope. If filament compliance and the discrete nature of the binding sites are taken into account, the negative slope is diminished, but not eliminated. This implies that there is an instability in the dynamics of individual half-sarcomeres. However, when the symmetric nature of whole sarcomeres is taken into account, filament rearrangement becomes important during the transient: as tension is recovered, some half-sarcomeres lengthen while others shorten. This leads to a flat T2 curve, as observed experimentally. In addition, we investigate the isotonic transient response and show that for a range of parameter values the model displays damped oscillations, as recently observed in experiments on single muscle fibers. We conclude that it is essential to consider the collective dynamics of many sarcomeres, rather than the dynamics of a single pair of filaments, when interpreting the transient response of muscle., Comment: 11 pages, 11 figures, Submitted to Biophysical Journal

A natural similarity in body dimensions of terrestrial animals noticed by ancient philosophers remains the main key to the problem of mammalian skeletal evolution with body mass explored in theoretical and experimental biology and tested by comparative zoologists. We discuss the long-standing problem of mammalian bone allometry commonly studied in terms of the so-called ''geometric'', ''elastic'', and ''static stress'' similarities by McMahon (1973, 1975a, 1975b). We revise the fundamental assumptions underlying these similarities and give new physical insights into geometric-shape and elastic-force constraints imposed on spatial evolution of mammalian long bones., Comment: 15 pages, 1 figure

Extracellular local field potentials (LFP) are usually modeled as arising from a set of current sources embedded in a homogeneous extracellular medium. Although this formalism can successfully model several properties of LFPs, it does not account for their frequency-dependent attenuation with distance, a property essential to correctly model extracellular spikes. Here we derive expressions for the extracellular potential that include this frequency-dependent attenuation. We first show that, if the extracellular conductivity is non-homogeneous, there is induction of non-homogeneous charge densities which may result in a low-pass filter. We next derive a simplified model consisting of a punctual (or spherical) current source with spherically-symmetric conductivity/permittivity gradients around the source. We analyze the effect of different radial profiles of conductivity and permittivity on the frequency-filtering behavior of this model. We show that this simple model generally displays low-pass filtering behavior, in which fast electrical events (such as Na$^+$-mediated action potentials) attenuate very steeply with distance, while slower (K$^+$-mediated) events propagate over larger distances in extracellular space, in qualitative agreement with experimental observations. This simple model can be used to obtain frequency-dependent extracellular field potentials without taking into account explicitly the complex folding of extracellular space., Comment: text (LaTeX), 6 figs. (ps)

Why is the genetic code the way it is? The most successful theory states that the codon assignments minimise the effects of errors arising in primordial living systems. Here a transversion is reported that leaves invariant degeneracy in the genetic code. A model was established to address the question: How should sets of synonymous codons be chosen to minimise the deleterious effects of non synonymous substitutions? It is shown that the arithmetic theory of quadratic forms provides the framework for degeneracy in the genetic code., Comment: 8 pages, 2 figures, modified title, added introduction and references, added paragraph with further concept, extended conclusion, acknowledgements and figure legends

A sound stimulus entering the inner ear excites a deformation of the basilar membrane which travels along the cochlea towards the apex. It is well established that this wave-like disturbance is amplified by an active system. Recently, it has been proposed that the active system consists of a set of self-tuned critical oscillators which automatically operate at an oscillatory instability. Here, we show how the concepts of a traveling wave and of self-tuned critical oscillators can be combined to describe the nonlinear wave in the cochlea., Comment: 5 pages, 2 figures

We introduce an approach for performing "very long" computer simulations of the dynamics of simplified, folded proteins. Using an alpha-carbon protein model and a fine grid to mimic continuum computations at increased speed, we perform unbiased simulations which exhibit many large-scale conformational transitions at low cost. In the case of the 72-residue N-terminal domain of calmodulin, the approach yields structural transitions between the calcium-free and calcium-bound structures at a rate of roughly one per day on a single Intel processor. Stable intermediates can be clearly characterized. The model employs Go-like interactions to stabilize two (or more) experimentally-determined structures. The approach is trivially parallelizable and readily generalizes to more complex potentials at minimal cost.

The growth of an actin network against an obstacle that stimulates branching locally is studied using several variants of a kinetic rate model based on the orientation-dependent number density of filaments. The model emphasizes the effects of branching and capping on the density of free filament ends. The variants differ in their treatment of side vs. end branching and dimensionality, and assume that new branches are generated by existing branches (autocatalytic behavior) or independently of existing branches (nucleation behavior). In autocatalytic models, the network growth velocity is rigorously independent of the opposing force exerted by the obstacle, and the network density is proportional to the force. The dependence of the growth velocity on the branching and capping rates is evaluated by a numerical solution of the rate equations. In side-branching models, the growth velocity drops gradually to zero with decreasing branching rate, while in end-branching models the drop is abrupt. As the capping rate goes to zero, it is found that the behavior of the velocity is sensitive to the thickness of the branching region. Experiments are proposed for using these results to shed light on the nature of the branching process., Comment: 6 figures

Protein representation and potential function are essential ingredients for studying proteins folding and protein prediction. We introduce a novel geometric representation of contact interactions using the edge simplices from alpha shape of protein structure. This representation can eliminate implausible neighbors not in physical contact, and can avoid spurious contact between two residues when a third residue is between them. We develop statistical alpha contact potential. A studentized bootstrap method is then introduced for assessing the 95% confidence intervals for each of the 210 parameters. We found with confidence that there is significant long range propensity (>30 residues apart) for hydrophobic interactions. We test alpha contact potential for native structure discrimination using several decoy sets, and found it often has comparable performance with atom-based potentials requiring more parameters. We also show that alpha contact potential has better performance than potential defined by cut-off distance between geometric centers of side chains. Clustering of alpha contact potentials reveals natural grouping of residues. To explore the relationship between shape representation and physicochemical representation, we test the minimum alphabet size for structure discrimination. We found that there is no significant difference in discrimination when alphabet size varies from 7 to 20, if geometry is represented accurately by alpha simplicial edges. This result suggests that the geometry of packing plays an important role, but the specific residue types are often interchangeable., Comment: 18 pages, 7 figures, and 6 tables. Accepted by Proteins

In this work, we revisit the basic reproduction rate $\mathcal{R}_{0}$ definition for analysis of epidemic-non-epidemic phases describing the dynamics of the discrete stochastic version of the epidemic $SIR$ model based on the Master Equation formalism. One shows that it is a very precise and efficient way to determine the epidemic threshold; using its most primitive concept, we can find exact results., Comment: Submitted to PRL

An important issue in theoretical epidemiology is the epidemic threshold phenomenon, which specify the conditions for an epidemic to grow or die out. In standard (mean-field-like) compartmental models the concept of the basic reproductive number, R0, has been systematically employed as a predictor for epidemic spread and as an analytical tool to study the threshold conditions. Despite the importance of this quantity, there are no general formulation of R0 when one considers the spread of a disease in a generic finite population, involving, for instance, arbitrary topology of inter-individual interactions and heterogeneous mixing of susceptible and immune individuals. The goal of this work is to study this concept in a generalized stochastic system described in terms of global and local variables. In particular, the dependence of R0 on the space of parameters that define the model is investigated; it is found that near of the ``classical'' epidemic threshold transition the uncertainty about the strength of the epidemic process still is significantly large. The forecasting attributes of R0 for a discrete finite system is discussed and generalized; in particular, it is shown that, for a discrete finite system, the pretentious predictive power of R0 is significantly reduced., Comment: To appear in "Journal of Biological Physics" - March 2003

DNA denaturation, wetting in two dimensions, depinning of a flux line, and other problems map onto a phase transition with effective long range interaction. It yields giant non-universal critical indexes, arbitrarily large macroscopic correlation length and fluctuations at a finite distance from the critical temperature. In the vicinity of this region the Gibbs distribution is invalid, and thermodynamics must be calculated from the first principles. There are no fluctuations above the critical temperature., Comment: 9 pages, no figures, an improved presentation of cond-mat/0212175

Biological approximations, which are universal for diverse species, are well known. With no other experimental data, their invariance to transformations from one species to another yields exact conservation (with respect to biological diversity and evolutionary history) laws, which are inconsistent with known physics and unique for self-organized live systems. The laws predict two and only two universal ways of biological diversity and evolution; their singularities; a new kind of rapid (compared to lifespan) adaptation and reversible mortality, which may be directed. Predictions agree with experimental data, and call for new concepts, insights, and microscopic theory., Comment: 15 pages, no figures, a generalization of cond-mat/0212174

Interictal spikes and sharp waves in human EEG are characteristic signatures of epilepsy. These potentials originate as a result of synchronous, pathological discharge of many neurons. The reliable detection of such potentials has been the long standing problem in EEG analysis, especially after long-term monitoring became common in investigation of epileptic patients. The traditional definition of a spike is based on its amplitude, duration, sharpness, and emergence from its background. However, spike detection systems built solely around this definition are not reliable due to the presence of numerous transients and artifacts. We use wavelet transform to analyze the properties of EEG manifestations of epilepsy. We demonstrate that the behavior of wavelet transform of epileptic spikes across scales can constitute the foundation of a relatively simple yet effective detection algorithm., Comment: 4 pages, 3 figures

We study the middle cerebral artery blood flow velocity (MCAfv) in humans using transcranial Doppler ultrasonography (TCD). Scaling properties of time series of the axial flow velocity averaged over a cardiac beat interval may be characterized by two exponents. The short time scaling exponent (STSE) determines the statistical properties of fluctuations of blood flow velocities in short-time intervals while the Hurst exponent describes the long-term fractal properties. In many migraineurs the value of the STSE is significantly reduced and may approach that of the Hurst exponent. This change in dynamical properties reflects the significant loss of short-term adaptability and the overall hyperexcitability of the underlying cerebral blood flow control system. We call this effect fractal rigidity., Comment: 4 pages, 6 figures

The interaction energies between the different types of bases of a single strand of DNA molecule have been calculated. Using these original values of energies the harmonic behavior of a number of base patterns of DNA has been studied. In view of the great interest aroused by the discovery of the role of the telomere segment of the DNA in the replication process and its possible link with the aging process, we have investigated, with simple models, the harmonic behavior of the telomeric pattern of bases as well as the thermodynamic response in the biological system. With these results a conclusion on the probable role of the telomeric pattern on aging has also been drawn. Here the calculated values of harmonic frequencies of the telomeric pattern of bases and of other possible patterns show that the telomeric pattern is associated with the highest vibrational frequency among all patterns of base combinations at the replication end of DNA. This seems to suggest that due to the existing telomeric pattern being closest to the frequencies of the electromagnetic radiation coming from sunlight, resonance of the telomeric frequency with such radiation may be responsible for damage to the reproductive ability of the cells and consequent aging and other problems. On the other hand in the last part of this work we have calculated the thermal vibrational amplitudes of the telomeric pattern and other possible patterns which show that the amplitude for the telomeric pattern is the least, and this suggests that the telomeric pattern is more mechanically and thermally stable than other possible patterns in the biological environment against damage from ordinary heat and mechanical effects., Comment: 10 pages, 3 figures, physics theories in DNA