scispace - formally typeset
Search or ask a question

Showing papers in "BioSystems in 2003"


Journal ArticleDOI
TL;DR: The most evolvable genetic architectures may often be those with an intermediate level of integration among characters, and in particular those where pleiotropic effects are variable and able to compensate for each other's constraints.
Abstract: Evolvability is the ability to respond to a selective challenge. This requires the capacity to produce the right kind of variation for selection to act upon. To understand evolvability we therefore need to understand the variational properties of biological organisms. Modularity is a variational property, which has been linked to evolvability. If different characters are able to vary independently, selection will be able to optimize each character separately without interference. But although modularity seems like a good design principle for an evolvable organism, it does not therefore follow that it is the only design that can achieve evolvability. In this essay I analyze the effects of modularity and, more generally, pleiotropy on evolvability. Although, pleiotropy causes interference between the adaptation of different characters, it also increases the variational potential of those characters. The most evolvable genetic architectures may often be those with an intermediate level of integration among characters, and in particular those where pleiotropic effects are variable and able to compensate for each other's constraints.

290 citations


Journal ArticleDOI
TL;DR: The use of a multiobjective EA (NSGA-II) has enabled a smaller gene subset size to correctly classify 100% or near 100% samples for three cancer samples and introduced a prediction strength threshold for determining a sample's belonging to one class or the other.
Abstract: In the area of bioinformatics, the identification of gene subsets responsible for classifying available disease samples to two or more of its variants is an important task. Such problems have been solved in the past by means of unsupervised learning methods (hierarchical clustering, self-organizing maps, k-mean clustering, etc.) and supervised learning methods (weighted voting approach, k-nearest neighbor method, support vector machine method, etc.). Such problems can also be posed as optimization problems of minimizing gene subset size to achieve reliable and accurate classification. The main difficulties in solving the resulting optimization problem are the availability of only a few samples compared to the number of genes in the samples and the exorbitantly large search space of solutions. Although there exist a few applications of evolutionary algorithms (EAs) for this task, here we treat the problem as a multiobjective optimization problem of minimizing the gene subset size and minimizing the number of misclassified samples. Moreover, for a more reliable classification, we consider multiple training sets in evaluating a classifier. Contrary to the past studies, the use of a multiobjective EA (NSGA-II) has enabled us to discover a smaller gene subset size (such as four or five) to correctly classify 100% or near 100% samples for three cancer samples (Leukemia, Lymphoma, and Colon). We have also extended the NSGA-II to obtain multiple non-dominated solutions discovering as much as 352 different three-gene combinations providing a 100% correct classification to the Leukemia data. In order to have further confidence in the identification task, we have also introduced a prediction strength threshold for determining a sample's belonging to one class or the other. All simulation results show consistent gene subset identifications on three disease samples and exhibit the flexibilities and efficacies in using a multiobjective EA for the gene subset identification task.

133 citations


Journal ArticleDOI
TL;DR: Although conflict mediation is pivotal to the emergence of individuality at the higher-level, the way in which the mediation is achieved can greatly affect the longer-term evolvability of the lineage.
Abstract: The continued well being of evolutionary individuals (units of selection and evolution) depends upon their evolvability, that is their capacity to generate and evolve adaptations at their level of organization, as well as their longer term capacity for diversifying into more complex evolutionary forms. During a transition from a lower- to higher-level individual, such as the transition between unicellular and multicellular organisms, the evolvability of the lower-level (cells) must be restricted, while the evolvability of the new higher-level unit (multicellular organism) must be enhanced. For these reasons, understanding the factors leading to an evolutionary transition should help us to understand the factors underlying the emergence of evolvability of a new evolutionary unit. Cooperation among lower-level units is fundamental to the origin of new functions in the higher-level unit. Cooperation can produce a new more complex evolutionary unit, with the requisite properties of heritable fitness variations, because cooperation trades fitness from a lower-level (the costs of cooperation) to the higher-level (the benefits for the group). For this reason, the evolution of cooperative interactions helps us to understand the origin of new and higher-levels of fitness and organization. As cooperation creates a new level of fitness, it also creates the opportunity for conflict between levels of selection, as deleterious mutants with differing effects at the two levels arise and spread. This conflict can interfere with the evolvability of the higher-level unit, since the lower and higher-levels of selection will often "disagree" on what adaptations are most beneficial to their respective interests. Mediation of this conflict is essential to the emergence of the new evolutionary unit and to its continued evolvability. As an example, we consider the transition from unicellular to multicellular organisms and study the evolution of an early-sequestered germ-line in terms of its role in mediating conflict between the two levels of selection, the cell and the cell group. We apply our theoretical framework to the evolution of germ/soma differentiation in the green algal group Volvocales. In the most complex member of the group, Volvox carteri, the potential conflicts among lower-level cells as to the "right" to reproduce the higher-level individual (i.e. the colony) have been mediated by restricting immortality and totipotency to the germ-line. However, this mediation, and the evolution of an early segregated germ-line, was achieved by suppressing mitotic and differentiation capabilities in all post-embryonic cells. By handicapping the soma in this way, individuality is ensured, but the solution has affected the long-term evolvability of this lineage. We think that although conflict mediation is pivotal to the emergence of individuality at the higher-level, the way in which the mediation is achieved can greatly affect the longer-term evolvability of the lineage.

98 citations


Journal ArticleDOI
TL;DR: The efficiency and superiority of the precluster scheme combined with thresholding is validated by comparison of the results for clustering both with and without preclustering for FDG-PET brain data of 13 healthy subjects.
Abstract: A new preprocessing clustering technique for quantification of kinetic PET data is presented. A two-stage clustering process, which combines a precluster and a classic hierarchical cluster analysis, provides data which are clustered according to a distance measure between time activity curves (TACs). The resulting clustered mean TACs can be used directly for estimation of kinetic parameters at the cluster level, or to span a vector space that is used for subsequent estimation of voxel level kinetics. The introduction of preclustering significantly reduces the overall time for clustering of multiframe kinetic data. The efficiency and superiority of the preclustering scheme combined with thresholding is validated by comparison of the results for clustering both with and without preclustering for FDG-PET brain data of 13 healthy subjects.

94 citations


Journal ArticleDOI
TL;DR: It is found that an univalent link exists between the dominant isomers of PABalphaCT, the dominant form of either acridine (proflavine) or acridan and the logic output of the system, which operates under the rules of Boolean algebra and performs as an "AND" logic gate.
Abstract: Enzyme-Based Logic Gates (ENLOGs) are key components in bio-molecular systems for information processing. This report and the previous one in this series address the characterization of two bio-molecular switching elements, namely the alpha-chymotrypsin (alphaCT) derivative p-phenylazobenzoyl-alpha-chymotrypsin (PABalphaCT) and its inhibitor (proflavine), as well as their assembly into a logic gate. The experimental output of the proposed system is expressed in terms of enzymic activity and this was translated into logic output (i.e. "1" or "0") relative to a predetermined threshold value. We have found that an univalent link exists between the dominant isomers of PABalphaCT (cis or trans), the dominant form of either acridine (proflavine) or acridan and the logic output of the system. Thus, of all possible combinations, only the trans-PABalphaCT and the acridan lead to an enzymic activity that can be defined as logic output "1". The system operates under the rules of Boolean algebra and performs as an "AND" logic gate.

83 citations


Journal ArticleDOI
TL;DR: It is suggested that unicellular organisms evolve largely through a process of metabolic change, resulting in biochemical diversity, not through extensive changes to cellular biochemistry, as in complex multicellular ones.
Abstract: The concept of evolvability covers a broad spectrum of, often contradictory, ideas. At one end of the spectrum it is equivalent to the statement that evolution is possible, at the other end are untestable post hoc explanations, such as the suggestion that current evolutionary theory cannot explain the evolution of evolvability. We examine similarities and differences in eukaryote and prokaryote evolvability, and look for explanations that are compatible with a wide range of observations. Differences in genome organisation between eukaryotes and prokaryotes meets this criterion. The single origin of replication in prokaryote chromosomes (versus multiple origins in eukaryotes) accounts for many differences because the time to replicate a prokaryote genome limits its size (and the accumulation of junk DNA). Both prokaryotes and eukaryotes appear to switch from genetic stability to genetic change in response to stress. We examine a range of stress responses, and discuss how these impact on evolvability, particularly in unicellular organisms versus complex multicellular ones. Evolvability is also limited by environmental interactions (including competition) and we describe a model that places limits on potential evolvability. Examples are given of its application to predator competition and limits to lateral gene transfer. We suggest that unicellular organisms evolve largely through a process of metabolic change, resulting in biochemical diversity. Multicellular organisms evolve largely through morphological changes, not through extensive changes to cellular biochemistry.

81 citations


Journal ArticleDOI
TL;DR: This paper used a hybrid algorithm combining particle swarm optimization with evolutionary algorithms to train HMMs for the alignment of protein sequences, which yields better alignments for a set of benchmark protein sequences than the most commonly applied HMM training methods.
Abstract: Multiple sequence alignment (MSA) is one of the basic problems in computational biology. Realistic problem instances of MSA are computationally intractable for exact algorithms. One way to tackle MSA is to use Hidden Markov Models (HMMs), which are known to be very powerful in the related problem domain of speech recognition. However, the training of HMMs is computationally hard and there is no known exact method that can guarantee optimal training within reasonable computing time. Perhaps the most powerful training method is the Baum-Welch algorithm, which is fast, but bears the problem of stagnation at local optima. In the study reported in this paper, we used a hybrid algorithm combining particle swarm optimization with evolutionary algorithms to train HMMs for the alignment of protein sequences. Our experiments show that our approach yields better alignments for a set of benchmark protein sequences than the most commonly applied HMM training methods, such as Baum-Welch and Simulated Annealing.

78 citations


Journal ArticleDOI
TL;DR: This paper employs a scale-invariant fitness landscape exhibiting significant ruggedness at all scales to exemplify the adaptive capacity of the composition model and explores and develops concepts using a simple abstract model of symbiotic composition to examine its impact on evolvability.
Abstract: Several of the major transitions in evolutionary history, such as the symbiogenic origin of eukaryotes from prokaryotes, share the feature that existing entities became the components of composite entities at a higher-level of organization. This composition of pre-adapted extant entities into a new whole is a fundamentally different source of variation from the gradual accumulation of small random variations, and it has some interesting consequences for issues of evolvability. Intuitively, the pre-adaptation of sets of features in reproductively independent specialists suggests a form of ‘divide and conquer’ decomposition of the adaptive domain. Moreover, the compositions resulting from one level may become the components for compositions at the next level, thus scaling-up the variation mechanism. In this paper, we explore and develop these concepts using a simple abstract model of symbiotic composition to examine its impact on evolvability. To exemplify the adaptive capacity of the composition model, we employ a scale-invariant fitness landscape exhibiting significant ruggedness at all scales. Whilst innovation by mutation and by conventional evolutionary algorithms becomes increasingly more difficult as evolution continues in this landscape, innovation by composition is not impeded as it discovers and assembles component entities through successive hierarchical levels. © 2002 Elsevier Science Ireland Ltd. All rights reserved.

76 citations


Journal ArticleDOI
TL;DR: An intracellular signalling model obtained by integrating several computational techniques into an agent-based paradigm is presented and the goal of a virtual laboratory based on this model and presently under development is discussed.
Abstract: The theory of behaviour-based systems (or autonomous agents) constitutes a useful approach for the modelling of intracellular signalling networks. In this sense, a cell can be seen as an adaptive autonomous agent or as a society of such agents, where each can exhibit a particular behaviour depending on its cognitive capabilities. We present an intracellular signalling model obtained by integrating several computational techniques into an agent-based paradigm. Cellulat, the model, takes into account two essential aspects of the intracellular signalling networks: (1) cognitive capacities, which are modelled as the agent abilities to interact with the surrounding medium and (2) a spatial organisation, this last obtained using a shared data structure through which the agents communicate between them. We propose a methodology for the modelling of intracellular signalling pathway using Cellulat and we discuss the goal of a virtual laboratory based on our model and presently under development.

72 citations


Journal ArticleDOI
TL;DR: The results suggest that mutation rate is a control parameter that governs a transition between two qualitatively different phases of evolution, an ordered phase characterized by punctuated equilibria of diversity, and a disordered phase of characterized by noisy fluctuations around an equilibrium diversity.
Abstract: We examine a simple form of the evolution of evolvability-the evolution of mutation rates-in a simple model system. The system is composed of many agents moving, reproducing, and dying in a two-dimensional resource-limited world. We first examine various macroscopic quantities (three types of genetic diversity, a measure of population fitness, and a measure of evolutionary activity) as a function of fixed mutation rates. The results suggest that (i) mutation rate is a control parameter that governs a transition between two qualitatively different phases of evolution, an ordered phase characterized by punctuated equilibria of diversity, and a disordered phase of characterized by noisy fluctuations around an equilibrium diversity, and (ii) the ability of evolution to create adaptive structure is maximized when the mutation rate is just below the transition between these two phases of evolution. We hypothesize that this transition occurs when the demands for evolutionary memory and evolutionary novelty are typically balanced. We next allow the mutation rate itself to evolve, and we observe that evolving mutation rates adapt to values at this transition. Furthermore, the mutation rates adapt up (or down) as the evolutionary demands for novelty (or memory) increase, thus supporting the balance hypothesis.

70 citations


Journal ArticleDOI
TL;DR: A simple model of genetic interactions, implemented in an evolutionary simulation, demonstrates that clustering of epistatically interacting genes increases the rate of adaptation, and can reorganize linkage patterns from random gene ordering into this more modular organization, thereby facilitating adaptation.
Abstract: A number of factors have been proposed that may affect the capacity for an evolutionary system to generate adaptation. One that has received little recent attention among biologists is linkage patterns, or the ordering of genes on chromosomes. In this study, a simple model of genetic interactions, implemented in an evolutionary simulation, demonstrates that clustering of epistatically interacting genes increases the rate of adaptation. Moreover, long-term evolution with inversion can reorganize linkage patterns from random gene ordering into this more modular organization, thereby facilitating adaptation. These results are consistent with a large body of biological observations and some mathematical theory. Although linkage patterns are neutral with respect to individual fitness in this model, they are subject to lineage level selection for evolvability. At least two candidate mechanisms may contribute to improved evolvability under epistatic clustering: clustering may reduce interference between selection on different traits, and it may allow the simultaneous optimization of different recombination rates for gene pairs with additive and epistatic fitness effects.

Journal ArticleDOI
TL;DR: A quantum formalism (Hilbert space probabilistic calculus) is developed for measurements performed over cognitive systems and is used for mathematical modelling of the functioning of consciousness as a self-measuring quantum-like system.
Abstract: We develop a quantum formalism (Hilbert space probabilistic calculus) for measurements performed over cognitive systems. In particular, this formalism is used for mathematical modelling of the functioning of consciousness as a self-measuring quantum-like system. By using this formalism, we could predict averages of cognitive observables. Reflecting the basic idea of neurophysiological and psychological studies on a hierarchic structure of cognitive processes, we use p-adic hierarchic trees as a mathematical model of a mental space. We also briefly discuss the general problem of the choice of an adequate mental geometry.

Journal ArticleDOI
TL;DR: A mathematical model of the underlying mechanisms of bone remodeling as it is mediated by PTH seems to indicate that the paradoxical observation that intermittent PTH administration causes net bone deposition while continuous administration causesNet bone loss may be attributed to the highly diversified dynamics which characterizes this nonlinear remodeling process.
Abstract: Bone, a major reservoir of body calcium, is under the hormonal control of the parathyroid hormone (PTH). Several aspects of its growth, turnover, and mechanism, occur in the absence of gonadal hormones. Sex steroids such as estrogen, nonetheless, play an important role in bone physiology, and are extremely essential to maintain bone balance in adults. In order to provide a basis for understanding the underlying mechanisms of bone remodeling as it is mediated by PTH, we propose here a mathematical model of the process. The nonlinear system model is then utilized to study the temporal effect of PTH as well as the action of estrogen replacement therapy on bone turnover. Analysis of the model is done on the assumption, supported by reported clinical evidence, that the process is characterized by highly diversified dynamics, which warrants the use of singular perturbation arguments. The model is shown to exhibit limit cycle behavior, which can develop into chaotic dynamics for certain ranges of the system's parametric values. Effects of estrogen and PTH administrations are then investigated by extending on the core model. Analysis of the model seems to indicate that the paradoxical observation that intermittent PTH administration causes net bone deposition while continuous administration causes net bone loss, and certain other reported phenomena may be attributed to the highly diversified dynamics which characterizes this nonlinear remodeling process.

Journal ArticleDOI
René Thomsen1
TL;DR: The introduced DockEA using the best settings found obtained the overall best docking solutions compared to the Lamarckian GA (LGA) provided with AutoDock and proved to be more robust than the LGA on the more difficult problems with a high number of flexible torsion angles.
Abstract: The docking of ligands to proteins can be formulated as a computational problem where the task is to find the most favorable energetic conformation among the large space of possible protein-ligand complexes. Stochastic search methods such as evolutionary algorithms (EAs) can be used to sample large search spaces effectively and is one of the commonly used methods for flexible ligand docking. During the last decade, several EAs using different variation operators have been introduced, such as the ones provided with the AutoDock program. In this paper we evaluate the performance of different EA settings such as choice of variation operators, population size, and usage of local search. The comparison is performed on a suite of six docking problems previously used to evaluate the performance of search algorithms provided with the AutoDock program package. The results from our investigation confirm that the choice of variation operators has an impact on the search-capabilities of EAs. The introduced DockEA using the best settings found obtained the overall best docking solutions compared to the Lamarckian GA (LGA) provided with AutoDock. Furthermore, the DockEA proved to be more robust than the LGA (in terms of reproducing the results in several runs) on the more difficult problems with a high number of flexible torsion angles.

Journal ArticleDOI
TL;DR: The proposed GA predicts which specific canonical base pairs will form hydrogen bonds and build helices, also known as stems, and shows that the Keep-Best Reproduction operator has similar benefits as in the traveling salesman problem domain.
Abstract: This paper presents a Genetic Algorithm (GA) to predict the secondary structure of RNA molecules, where the secondary structure is encoded as a permutation. More specifically, the proposed algorithm predicts which specific canonical base pairs will form hydrogen bonds and build helices, also known as stems. Since RNA is involved in both transcription and translation and also has catalytic and structural roles in the cell, determining the structure of RNA is of fundamental importance in helping to determine RNA function. We introduce a GA where a permutation is used to encode the secondary structure of RNA molecules. We discuss results on RNA sequences of lengths 76, 210, 681, and 785 nucleotides and present several improvements to our algorithm. We show that the Keep-Best Reproduction operator has similar benefits as in the traveling salesman problem domain. In addition, a comparison of several crossover operators is provided. We also compare the results of the permutation-based GA with a binary GA, demonstrating the benefits of the newly proposed representation.

Journal ArticleDOI
TL;DR: The goal of the MGS project is to develop an experimental programming language dedicated to the simulation of this kind of systems, enabling the specification of spatially localized computations on heterogeneous entities.
Abstract: The cell as a dynamical system presents the characteristics of having a dynamical structure. That is, the exact phase space of the system cannot be fixed before the evolution and integrative cell models must state the evolution of the structure jointly with the evolution of the cell state. This kind of dynamical systems is very challenging to model and simulate. New programming concepts must be developed to ease their modeling and simulation. In this context, the goal of the MGS project is to develop an experimental programming language dedicated to the simulation of this kind of systems. MGS proposes a unified view on several computational mechanisms (CHAM, Lindenmayer systems, Paun systems, cellular automata) enabling the specification of spatially localized computations on heterogeneous entities. The evolution of a dynamical structure is handled through the concept of transformation which relies on the topological organization of the system components. An example based on the modeling of spatially distributed biochemical networks is used to illustrate how these notions can be used to model the spatial and temporal organization of intracellular processes.

Journal ArticleDOI
TL;DR: The transition from a domain with mind and matter unseparated to separate mental and material domains can be viewed as a result of a general kind of symmetry breaking, which can be described formally in terms of inequivalent representations.
Abstract: Many philosophical and scientific discussions of topics of mind-matter research make implicit assumptions, in various guises, about the distinction between mind and matter. Currently predominant positions are based on either reduction or emergence, providing either monistic or dualistic scenarios. A more-involved framework of thinking, which can be traced back to Spinoza and Leibniz, combines the two scenarios, dualistic (with mind and matter separated) and monistic (with mind and matter unseparated), in one single picture. Based on such a picture, the transition from a domain with mind and matter unseparated to separate mental and material domains can be viewed as a result of a general kind of symmetry breaking, which can be described formally in terms of inequivalent representations. The possibility of whether this symmetry breaking might be connected to the emergence of temporal directions from temporally non-directed or even non-temporal levels of reality will be discussed. Correlations between mental and material aspects of reality could then be imagined as remnants of such primordial levels. Different conceivable types of inequivalent representations would lead to correlations with different characteristics.

Journal ArticleDOI
TL;DR: This work focuses on several scenarios in which the inherent molecular fluctuations are not merely a nuisance, but act constructively and bring about qualitative changes in the dynamics of the system.
Abstract: Biochemical and genetic regulatory systems that involve low concentrations of molecules are inherently noisy. This intrinsic stochasticity has received considerable interest recently, leading to new insights about the sources and consequences of noise in complex systems of genetic regulation. However, most prior work was devoted to the reduction of fluctuation and the robustness of cellular function with respect to intrinsic noise. Here, we focus on several scenarios in which the inherent molecular fluctuations are not merely a nuisance, but act constructively and bring about qualitative changes in the dynamics of the system. It will be demonstrated that in many typical situations biochemical and genetic regulatory systems may utilize intrinsic noise to their advantage.

Journal ArticleDOI
TL;DR: It is shown that if mutation replaces an amino acid, the change of hydrophobicity is generally weak, while that of size is strong, and the antisymmetrical correlation between the amino acid size and the degeneracy number is known.
Abstract: The first information system emerged on the earth as primordial version of the genetic code and genetic texts. The natural appearance of arithmetic power in such a linguistic milieu is theoretically possible and practical for producing information systems of extremely high efficiency. In this case, the arithmetic symbols should be incorporated into an alphabet, i.e. the genetic code. A number is the fundamental arithmetic symbol produced by the system of numeration. If the system of numeration were detected inside the genetic code, it would be natural to expect that its purpose is arithmetic calculation e.g., for the sake of control, safety, and precise alteration of the genetic texts. The nucleons of amino acids and the bases of nucleic acids seem most suitable for embodiments of digits. These assumptions were used for the analyzing the genetic code. The compressed, life-size, and split representation of the Escherichia coli and Euplotes octocarinatus code versions were considered simultaneously. An exact equilibration of the nucleon sums of the amino acid standard blocks and/or side chains was found repeatedly within specified sets of the genetic code. Moreover, the digital notations of the balanced sums acquired, in decimal representation, the unique form 111, 222...., 999. This form is a consequence of the criterion of divisibility by 037. The criterion could simplify some computing mechanism of a cell if any and facilitate its computational procedure. The cooperative symmetry of the genetic code demonstrates that possibly a zero was invented and used by this mechanism. Such organization of the genetic code could be explained by activities of some hypothetical molecular organelles working as natural biocomputers of digital genetic texts. It is well known that if mutation replaces an amino acid, the change of hydrophobicity is generally weak, while that of size is strong. The antisymmetrical correlation between the amino acid size and the degeneracy number is known as well. It is shown that these and some other familiar properties may be a physicochemical effect of arithmetic inside the genetic code. The "frozen accident" model, giving unlimited freedom to the mapping function, could optimally support the appearance of both arithmetic symbols and physicochemical protection inside the genetic code.

Journal ArticleDOI
TL;DR: This work proposes a method for coping with combinatorial explosion of possible routes across metabolic networks by suitably classifying metabolites as external or internal, and suggests to find such a classification of metabolites that minimizes the number of elementary flux modes (pathways).
Abstract: Metabolic pathway analysis based on the concept of elementary flux mode is a valuable tool for reconstruction of bacterial metabolisms and in predicting optimal conversion yields in biotechnology. However, pathway analysis of large and highly entangled metabolic networks meets the problem of combinatorial explosion of possible routes across the networks. Here we propose a method for coping with this problem by suitably classifying metabolites as external or internal. External metabolites are considered to have buffered concentrations while internal metabolites have to fulfil a balance condition at steady state. For many substances such as nutrients and excreted products, there are biochemical reasons to classify them as external. In addition, other substances (especially at central branching points) can operationally be considered external in order to avoid combinatorial explosion. We suggest to find such a classification of metabolites that minimizes the number of elementary flux modes (pathways). This is motivated by the objectives of finding such a description of the system that reduces as much as possible the amount of necessary data and of removing the ambiguity and arbitrariness in the classification of metabolites in an automated, systematic way. For networks of moderate size, the solution to this combinatorial minimization problem can be found by exhaustive search. To tackle also larger systems, a stochastic optimization program based on the Metropolis algorithm was developed. Both methods are applied, for illustration, to several reaction schemes including a larger network representing glutathione metabolism.

Journal ArticleDOI
TL;DR: It is proved how the DNA operations presented by Adleman and Lipton can be used for developing DNA algorithms to resolving the set cover problem and the problem of exact cover by 3-sets.
Abstract: Adleman wrote the first paper in which it is shown that deoxyribonucleic acid (DNA) strands could be employed towards calculating solutions to an instance of the NP-complete Hamiltonian path problem (HPP). Lipton also demonstrated that Adleman's techniques could be used to solve the NP-complete satisfiability (SAT) problem (the first NP-complete problem). In this paper, it is proved how the DNA operations presented by Adleman and Lipton can be used for developing DNA algorithms to resolving the set cover problem and the problem of exact cover by 3-sets.

Journal ArticleDOI
TL;DR: It is shown that appropriate amounts of the phasereset can prevent the model from falling, even for the perturbation that induces falling in the case without the phase reset, suggesting that those phase resets can improve the dynamic stability of the gait.
Abstract: The human walking movement shows transient changes in response to single short-lived external perturbations, termed “stumbling reactions.” During the stumbling reactions, the walking phase is reset. It has been considered that the reactions contribute to stabilizing the motion, but less evidence bridging between the rhythm reset and the dynamic stability of the gait has been provided. The present study tries to establish the relationship between them. To this end, we construct a simple dynamical system model of the human musculo-skeletal system interacting with the ground, whose joint kinematics during walking is constrained by a given periodic joint-angles-profile. We show first that the model can exhibit a stable limit cycle corresponding to the steady walking with no perturbations. The responses of the limit cycle oscillation are examined by applying a type of perturbations at various timings with various intensities, elucidating the stability of the model’s walking when no phase reset is performed. We then observe that modifications of the periodic joint-angles-profile within a short time interval in response to the perturbation can alter the responses of the limit cycle oscillation and induce phase reset of the model’s walking. It is shown that appropriate amounts of the phase reset can prevent the model from falling, even for the perturbation that induces falling in the case without the phase reset. This suggests that those phase resets can improve the dynamic stability of the gait. Moreover, the appropriate phase resets predicted by the model are compared with the experimentally observed phase resets during human stumbling reaction to show they share similar characteristics.

Journal ArticleDOI
TL;DR: A dynamical analysis of P systems is given that is focused on basic phenomena of biological relevance and a discrete formulation of the Belousov-Zhabotinsky (BZ) reaction is given in terms of PBE systems.
Abstract: A dynamical analysis of P systems is given that is focused on basic phenomena of biological relevance. After a short presentation of a new kind of P systems (PB systems), membrane systems with environment, called PBE systems, are introduced that are more suitable for modeling complex membrane interactions. Some types of periodicity and non-periodicity are considered for PBE systems by showing some “minimal” examples of systems that exhibit these properties. In particular, a discrete formulation of the Belousov–Zhabotinsky (BZ) reaction is given in terms of PBE systems. Some questions and open problems for future research are indicated.

Journal ArticleDOI
TL;DR: Whether the Petri net approach is capable of identifying biochemical networks that are consistent with disease susceptibility due to higher order nonlinear interactions between three DNA sequence variations is evaluated.
Abstract: Understanding how DNA sequence variations impact human health through a hierarchy of biochemical and physiological systems is expected to improve the diagnosis, prevention, and treatment of common, complex human diseases. We have previously developed a hierarchical dynamic systems approach based on Petri nets for generating biochemical network models that are consistent with genetic models of disease susceptibility. This modeling approach uses an evolutionary computation approach called grammatical evolution as a search strategy for optimal Petri net models. We have previously demonstrated that this approach routinely identifies biochemical network models that are consistent with a variety of genetic models in which disease susceptibility is determined by nonlinear interactions between two DNA sequence variations. In the present study, we evaluate whether the Petri net approach is capable of identifying biochemical networks that are consistent with disease susceptibility due to higher order nonlinear interactions between three DNA sequence variations. The results indicate that our model-building approach is capable of routinely identifying good, but not perfect, Petri net models. Ideas for improving the algorithm for this high-dimensional problem are presented.

Journal ArticleDOI
TL;DR: The aim of this paper is to introduce bio-inspired computing tissues that might constitute a key concept for the implementation of 'living' machines and the POE model that classifies bio- inspired machines along three axes.
Abstract: Biological inspiration in the design of computing machines could allow the creation of new machines with promising characteristics such as fault-tolerance, self-replication or cloning, reproduction, evolution, adaptation and learning, and growth. The aim of this paper is to introduce bio-inspired computing tissues that might constitute a key concept for the implementation of ‘living’ machines. We first present a general overview of bio-inspired systems and the POE model that classifies bio-inspired machines along three axes. The Embryonics project—inspired by some of the basic processes of molecular biology—is described by means of the BioWatch application, a fault-tolerant and self-repairable watch. The main characteristics of the Embryonics project are the multicellular organization, the cellular differentiation, and the self-repair capabilities. The BioWall is intended as a reconfigurable computing tissue, capable of interacting with its environment by means of a large number of touch-sensitive elements coupled with a color display. For illustrative purposes, a large-scale implementation of the BioWatch on the BioWall's computational tissue is presented. We conclude the paper with a description of bio-inspired computing tissues and POEtic machines.

Journal ArticleDOI
TL;DR: StarLogo and the Adventures in Modeling Curriculum provide an easily accessible entry point into complex systems modeling for students and other novice modelers, including specific applications to epidemiological and ecological systems.
Abstract: Research on complex, adaptive systems has made significant advances in recent years in the study of natural and social phenomena that exhibit random variation and selection, resulting in learning or evolution. Unfortunately, students (including K-12, undergraduate and graduate) in most biology programs have little opportunity to explore complex systems during the course of their studies. StarLogo and the Adventures in Modeling Curriculum [Adventures in Modeling: Exploring Complex, Dynamic Systems with StarLogo. Teachers College Press, New York] provide an easily accessible entry point into complex systems modeling for students and other novice modelers. These specialized tools can provide powerful insights into the dynamics of systems and create opportunities to explore challenging and meaningful domains in the biological sciences. Specific applications to epidemiological and ecological systems are explored, including the often debated topic of the evolution of reduced attack rates in predator-prey systems.

Journal ArticleDOI
TL;DR: It is shown that membrane systems with antiport carriers provide an appropriate model for distributed computing, particularly for message-passing algorithms interpreted here as membrane transport in both directions, namely when two chemicals behave as input and output messages and pass the membranes in both direction using antiport carrier.
Abstract: This paper presents fundamental distributed algorithms over membrane systems with antiport carriers. We describe distributed algorithms for collecting and dispersing information, leader election in these systems, and the mutual exclusion problem. Finally, we consider membrane systems producing correct results despite some failures at some of the components or the communication links. We show that membrane systems with antiport carriers provide an appropriate model for distributed computing, particularly for message-passing algorithms interpreted here as membrane transport in both directions, namely when two chemicals behave as input and output messages and pass the membranes in both directions using antiport carriers.

Journal ArticleDOI
TL;DR: The results of this work further confirm the growing indication that evolutionary computation can outperform backpropagation as a method of artificial neural network training and indicate the degree to which bias in the initial training and testing data can affect performance and the importance of bootstrapping.
Abstract: Artificial neural networks (ANNs) can be utilized to generate predictive models of quantitative structure-activity relationships between a set of molecular descriptors and activity. Evolutionary computation provides a means to appropriately search for the set of weights and bias terms associated with artificial neural networks that minimize selected functions of the error between the actual and desired outputs. This method is demonstrated by evolutionary training of artificial neural networks capable of predicting anti-HIV activity for a set of 1-[(2-hydroxyethoxy)methyl]-6-(phenylthio)thymine (HEPT) derivatives. The results of this work further confirm the growing indication that evolutionary computation can outperform backpropagation as a method of artificial neural network training. The results also indicate the degree to which bias in the initial training and testing data can affect performance and the importance of bootstrapping.

Journal ArticleDOI
TL;DR: Modulation of neuronal impulse pattern is examined by means of a simplified Hodgkin-Huxley type computer model which refers to experimental recordings of cold receptor discharges to elucidate different resonance behaviors.
Abstract: Modulation of neuronal impulse pattern is examined by means of a simplified Hodgkin-Huxley type computer model which refers to experimental recordings of cold receptor discharges. This model essentially consists of two potentially oscillating subsystems: a spike generator and a subthreshold oscillator. With addition of noise the model successfully mimics the major types of experimentally recorded impulse patterns and thereby elucidate different resonance behaviors. (1) There is a range of rhythmic spiking or bursting where the spike generator is strongly coupled to the subthreshold oscillator. (2) There is a pacemaker activity of more complex interactions where the spike generator has overtaken part of the control. (3) There is a situation where the two subsystems are decoupled and only resonate with the help of noise.

Journal ArticleDOI
TL;DR: It is shown that sexual reproduction can alleviate the problem of genetic linkage by recombining separate modules all of which incorporate either favorable or unfavorable mutations, and it is speculated that this effect may contribute to the taxonomic prevalence of sexual reproduction among higher organisms.
Abstract: What genotypic features explain the evolvability of organisms that have to accomplish many different tasks? The genotype of behaviorally complex organisms may be more likely to encode modular neural architectures because neural modules dedicated to distinct tasks avoid neural interference, i.e. the arrival of conflicting messages for changing the value of connection weights during learning. However, if the connection weights for the various modules are genetically inherited, this raises the problem of genetic linkage: favorable mutations may fall on one portion of the genotype encoding one neural module and unfavorable mutations on another portion encoding another module. We show that this can prevent the genotype from reaching an adaptive optimum. This effect is different from other linkage effects described in the literature and we argue that it represents a new class of genetic constraints. Using simulations we show that sexual reproduction can alleviate the problem of genetic linkage by recombining separate modules all of which incorporate either favorable or unfavorable mutations. We speculate that this effect may contribute to the taxonomic prevalence of sexual reproduction among higher organisms. In addition to sexual recombination, the problem of genetic linkage for behaviorally complex organisms may be mitigated by entrusting evolution with the task of finding appropriate modular architectures and learning with the task of finding the appropriate connection weights for these architectures.