Compositional Assemblies Behave Similarly to Quasispecies Model

Total Page:16

File Type:pdf, Size:1020Kb

Compositional Assemblies Behave Similarly to Quasispecies Model Compositional Assemblies Behave Similarly to Quasispecies Model Renan Gross, Omer Markovitch and Doron Lancet Department of Molecular Genetics, Weizmann Institute of Science, Israel Group meeting, 01/10/2013 Introduction Quasispecies are a cloud of genotypes that appear in a population at mutation-selection balance. (J.J Bull et al, PLoS 2005) - Theoretical model (equations and assumptions), with experimental support by RNA viruses. - Usually applied when mutation rates are high. - GARD composomes replicate with relatively low fidelity (high mutation rate). Could they show similar dynamic behavior? 2 Quasispecies model • Basically a population model • n different genotypes / identities 3 Quasispecies model • Basically a population model • n different genotypes / identities Reactor 4 Quasispecies model • Basically a population model • n different genotypes / identities • Their relative concentration in the environment is denoted by Reactor the fraction xi. 푛 ෍ 푥푖 = 1 푖=1 5 Quasispecies model • Basically a population model • n different genotypes / identities • Their relative concentration in the environment is denoted by Reactor the fraction xi. 푛 • Goal: To find out how xi behave as a function of time. ෍ 푥푖 = 1 푖=1 6 Quasispecies model • Each one replicates at a certain rate – how many offspring it has per unit time. • Some replicate faster than others. • This is called the replication rate, denoted Ai 푛푢푚푏푒푟 표푓 표푓푓푠푝푟푖푛푔 • 퐴 = 7 푖 푡푖푚푒 Quasispecies model • However, replication is not exact. Sometimes, the offspring is of another genotype. • The chance that a genotype j replicates into genotype i is denoted Qij. Q11 = 0.5 Q21 = 0.2 Q31 = 0.2 Q41 = 0.1 8 Quasispecies model • We can put everything in a matrix, called the transition matrix, Q. From ퟎ. ퟓ 0.1 0 0.7 To 0.2 ퟎ. ퟓ 0 0.2 0.2 0.2 ퟏ 0 0.1 0.2 0 ퟎ. ퟏ • The main diagonal is faithful self replication. 9 Quasispecies model • How can we find out Q and A? • Assume that we are dealing with RNA sequences of length v: – There is a single genotype with highest replication rate: the master sequence 10 Quasispecies model • How can we find out Q and A? • Assume that we are dealing with RNA sequences of length v: – There is a single genotype with highest q 1 replication rate: the master sequence 1 – Single digit replication: 0 ≤ q ≤ 1 1-q 0 • What is the chance for error-less replication? 11 Quasispecies model • How can we find out Q and A? • Assume that we are dealing with RNA sequences of length v: – There is a single genotype with highest q 1 replication rate: the master sequence 1 – Single digit replication: 0 ≤ q ≤ 1 1-q 0 • What is the chance for error-less replication? – All mutations of the master sequence replicate slower according to the Hamming distance. • Effectively, many genotypes are grouped together. • Q and A are built only from q and v. 12 Quasispecies model 13 Quasispecies model • q = 1 14 Quasispecies model • q = 1 exact replication 15 Quasispecies model • q = 1 exact replication • q = 0 16 Quasispecies model • q = 1 exact replication • q = 0 exact complimentary replication 17 Quasispecies model • q = 1 exact replication • q = 0 exact complimentary replication • q = 0.5 18 Quasispecies model • q = 1 exact replication • q = 0 exact complimentary replication all data is lost. • q = 0.5 19 Quasispecies model • q = 1 exact replication • q = 0 exact complimentary replication all data is lost. • q = 0.5 • Starting at q = 1, lowering it results in loss of the master sequence as the most frequent genotype. – (requires lack of back mutation) ERROR CATASTROPHE • RNA viruses may be fought by bringing them to error catastrophe. 20 Quasispecies model • Under constant population assumptions, the transition matrix Q and the replication rates A are all that are needed in order to find out the concentration dynamics. 21 Quasispecies model • Under constant population assumptions, the transition matrix Q and the replication rates A are all that are needed in order to find out the concentration dynamics. • The Eigen equation (after Manfred Eigen): 푑푥 푖 = 퐴 푄 − 퐸෨ 푡 푥 + ෍ 퐴 푄 푥 푑푡 푖 푖푖 푖 푖 푖푗 푗 푗≠푖 22 Quasispecies model • Under constant population assumptions, the transition matrix Q and the replication rates A are all that are needed in order to find out the concentration dynamics. • The Eigen equation (after Manfred Eigen): 푑푥 푖 = 퐴 푄 − 퐸෨ 푡 푥 + ෍ 퐴 푄 푥 푑푡 푖 푖푖 푖 푖 푖푗 푗 푗≠푖 • Where E is the “average excess rate” 푛 퐸෨ 푡 = ෍ 퐴푖푥푖 푖=1 23 Quasispecies model • Q and A can be combined into one matrix W, which tells how much of each genotype is produced per unit time. 푾 = 푸 ⋅ 푑푖푎푔(푨) 24 Examples: 1 0.001 0.001 0.01 0.1 2 0.01 0.01 푾 = 0.001 0.1 3 0.01 0.001 0.001 0.001 4 • Initial conditions: x1 = 1, all others are 0. 25 Examples: 1 0.001 0.001 0.01 0.1 2 0.01 0.01 푾 = 0.001 0.1 3 0.01 0.001 0.001 0.001 4 • Initial conditions: x1 = 1, all others are 0. 26 Examples: 1 0.001 0.001 1 0.1 2 0.01 1 푾 = 0.001 0.1 3 1 0.001 0.001 0.001 4 • Initial conditions: x1 = 1, all others are 0. 27 Examples: 1 0.001 0.001 1 0.1 2 0.01 1 푾 = 0.001 0.1 3 1 0.001 0.001 0.001 4 • Initial conditions: x1 = 1, all others are 0. 28 • The most frequent genotype is not necessarily the one with highest Ai. • The steady state population distribution is called the quasispecies. – That means, a vertical slice 29 GARD model (Graded Autocatalysis Replication Domain) RNA world Lipid world DNA / RNA / Polymers Sequence covalent bonds 30 GARD model (Graded Autocatalysis Replication Domain) RNA world Lipid world Assemblies / Clusters / DNA / RNA / Polymers Vesicles / Membranes Sequence Composition covalent bonds non-covalent bonds Segre and Lancet, EMBO Reports 1 (2000) 31 GARD model (Graded Autocatalysis Replication Domain) Synthetic chemistry Kinetic model Catalytic network (b) of rate-enhancement values Segre, Ben-Eli and Lancet, Proc. Natl. Acad. Sci. 97 (2000) GARD model (Graded Autocatalysis Replication Domain) Synthetic chemistry Kinetic model Catalytic network (b) of rate-enhancement values Rate enhancement dn NG n i k N k n 1 b j i 1...N Molecular repertoire f i b i ij G dt j1 N Segre, Ben-Eli and Lancet, Proc. Natl. Acad. Sci. 97 (2000) GARD model (Graded Autocatalysis Replication Domain) Synthetic chemistry Kinetic model Homeostatic growth Catalytic network (b) of rate-enhancement values b Fission / Split Rate enhancement dn NG n i k N k n 1 b j i 1...N Molecular repertoire f i b i ij G dt j1 N Segre, Ben-Eli and Lancet, Proc. Natl. Acad. Sci. 97 (2000) GARD model (Graded Autocatalysis Replication Domain) Following a single lineage. 35 GARD model (Graded Autocatalysis Replication Domain) Following a single lineage. Similarityng=30; split=1.5; seed=361 ‘carpet’ 1000 1 800 0.8 600 0.6 400 0.4 Generation 200 0.2 Compositional Similarity 0 200 400 600 800 1000 Generation 36 GARD model (Graded Autocatalysis Replication Domain) Following a single lineage. Similarityng=30; split=1.5; seed=361 ‘carpet’ 1000 1 800 0.8 600 0.6 400 0.4 Generation 200 0.2 Compositional Similarity 0 200 400 600 800 1000 Generation Composome (compositional genome): a faithfully replicating composition/assembly. Compotype (composome type): a collection of similar composomes. Molecular Compotype: the center of mass of the compotype cloud, treated as a molecular assembly. 37 GARD model (Graded Autocatalysis Replication Domain) • Idea: if we ignore the stochasticity inherent in the model, then errorless-replication occurs according to beta matrix eigenvectors. dn NG n i k N k n 1 b j i 1...N f i b i ij G dt j1 N 38 Moment Algebraux • Multiplying a matrix by a vector gives another vector • 푨 ⋅ 푥Ԧ = 푦Ԧ 푎 푎 푥 푎 푥 + 푎 푥 푦 • 11 12 ⋅ 1 = 11 1 12 2 = 1 푎21 푎22 푥2 푎21푥1 + 푎22푥2 푦2 • An eigenvector is a vector 푥Ԧ such that: • 푨 ⋅ 푥Ԧ = 휆푥Ԧ – l is called the eigenvalue. – l may be complex, and so may the eigenvector. 39 Moment Algebraux • The Perron-Frobenius Theorem: – A matrix with strictly positive entries contains a maximal real eigenvalue. – Its eigenvector is real and non-negative. In fact, it’s the only one with this property. • As molecular assemblies must contain real non- negative number of molecules, this looks interesting. 40 GARD and Quasispecies • Do GARD population dynamics behave like the quasispecies model? • What we would like to do: – For each assembly, experimentally find out the transition frequencies and replication rates – In other words: find Q and A. • Problem: 199 – N = 100, n = 100 There are possible G max 100 assemblies ( ~4 ⋅ 1058) 41 GARD and Quasispecies • Solution: group some assemblies together and treat them as one genotype. • We decided to group together by distances from the eigenvector. N space (simplified to 2d) 2 3 G 2 Molecule 1 Eigenvector Molecule 1 42 GARD and Quasispecies • Problem: how do we sample such a large space? • 30000 assemblies are randomly generated. 43 GARD and Quasispecies • Problem: how do we sample such a large space? • 30000 assemblies are randomly generated. – By filling up assemblies until they reach nmax. • Assemblies generated this way are far from the target. 43 GARD and Quasispecies • Problem: how do we sample such a large space? • 30000 assemblies are randomly generated.
Recommended publications
  • Prebiological Evolution and the Metabolic Origins of Life
    Prebiological Evolution and the Andrew J. Pratt* Metabolic Origins of Life University of Canterbury Keywords Abiogenesis, origin of life, metabolism, hydrothermal, iron Abstract The chemoton model of cells posits three subsystems: metabolism, compartmentalization, and information. A specific model for the prebiological evolution of a reproducing system with rudimentary versions of these three interdependent subsystems is presented. This is based on the initial emergence and reproduction of autocatalytic networks in hydrothermal microcompartments containing iron sulfide. The driving force for life was catalysis of the dissipation of the intrinsic redox gradient of the planet. The codependence of life on iron and phosphate provides chemical constraints on the ordering of prebiological evolution. The initial protometabolism was based on positive feedback loops associated with in situ carbon fixation in which the initial protometabolites modified the catalytic capacity and mobility of metal-based catalysts, especially iron-sulfur centers. A number of selection mechanisms, including catalytic efficiency and specificity, hydrolytic stability, and selective solubilization, are proposed as key determinants for autocatalytic reproduction exploited in protometabolic evolution. This evolutionary process led from autocatalytic networks within preexisting compartments to discrete, reproducing, mobile vesicular protocells with the capacity to use soluble sugar phosphates and hence the opportunity to develop nucleic acids. Fidelity of information transfer in the reproduction of these increasingly complex autocatalytic networks is a key selection pressure in prebiological evolution that eventually leads to the selection of nucleic acids as a digital information subsystem and hence the emergence of fully functional chemotons capable of Darwinian evolution. 1 Introduction: Chemoton Subsystems and Evolutionary Pathways Living cells are autocatalytic entities that harness redox energy via the selective catalysis of biochemical transformations.
    [Show full text]
  • Lack of Evolvability in Self-Sustaining Autocatalytic Networks Constraints Metabolism-first Scenarios for the Origin of Life
    Lack of evolvability in self-sustaining autocatalytic networks constraints metabolism-first scenarios for the origin of life Vera Vasasa,b, Eörs Szathmárya,b,c,1, and Mauro Santosa,d aCollegium Budapest, Institute for Advanced Study, H-1014 Budapest, Hungary; bInstitute of Biology, Eötvös University, H-1117 Budapest, Hungary; cParmenides Center for the Study of Thinking, Munich D-80333, Germany; and dDepartament de Genètica i de Microbiologia, Grup de Biologia Evolutiva, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Spain Edited by Gerald F. Joyce, The Scripps Research Institute, La Jolla, CA, and approved December 3, 2009 (received for review November 3, 2009) A basic property of life is its capacity to experience Darwinian reactions that may deplete certain reactants (12) and dynamical evolution. The replicator concept is at the core of genetics-first aspects of autocatalytic cycles if they are assumed to coexist in theories of the origin of life, which suggest that self-replicating abstract space (13). Even if we ignore such hurdles, the key oligonucleotides or their similar ancestors may have been the first question still remains: Was a network of chemical reactions able to “living” systems and may have led to the evolution of an RNA increase in complexity and eventually undergo Darwinian selec- world. But problems with the nonenzymatic synthesis of biopol- tion as assumed by their advocates? A basic condition for any ymers and the origin of template replication have spurred the nascent Darwinian process in a population of self-reproducing alternative metabolism-first scenario, where self-reproducing systems is that they must have a sort of hereditary transmission and evolving proto-metabolic networks are assumed to have pre- which requires, in turn, becoming familiar with a lesser-known and dated self-replicating genes.
    [Show full text]
  • The Origin of Life: What We Do and Don't Know
    The Origin of Life: What We Do and Don’t Know Jim Cleaves Blue Marble Space Institute for Science Image courtesy of the American Museum of Natural History How Do We Think Life Began on Earth? A View from 38 Years Ago… “It must be admitted from the beginning that we do not know how life began. It is generally believed that a variety of processes led to the formation of simple organic compounds on the primitive Earth. These compounds combined together to give more and more complex structures until one was formed that could be called living. No one should be satisfied with an explanation as general as this.” S.L. Miller and L.E. Orgel, The Origins of Life on Earth, 1974 NSCORT Exobiology Modern Historical Background Oparin postulated the self-organization of environmentally supplied compounds produced in an environment different from the modern one (e.g. with a different atmosphere); this was an outgrowth of 19th c. thinking on the problem. A.I. Oparin H.C. Urey S.L. Miller 1894-1980 1930-2007 1893-1981 Urey extended Oparin’s model to ideas corroborated by mid-20th c. conceptions of the origin of the solar system, backed up discoveries in isotope geochemistry. Miller’s electric discharge experiment provided an experimental validation of the possible link between the two concepts. The Origin of Biochemistry (as a Scientific Field) First amino acid discovered in 1806 (asparagine, interestingly threonine not until 1936!). Proteins were not identified as such until 1838; the role of amino acids in forming them was not understood until the 1860’s.
    [Show full text]
  • Lack of Evolvability in Self-Sustaining Autocatalytic Networks Constraints Metabolism-first Scenarios for the Origin of Life
    Lack of evolvability in self-sustaining autocatalytic networks constraints metabolism-first scenarios for the origin of life Vera Vasasa,b, Eörs Szathmárya,b,c,1, and Mauro Santosa,d aCollegium Budapest, Institute for Advanced Study, H-1014 Budapest, Hungary; bInstitute of Biology, Eötvös University, H-1117 Budapest, Hungary; cParmenides Center for the Study of Thinking, Munich D-80333, Germany; and dDepartament de Genètica i de Microbiologia, Grup de Biologia Evolutiva, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Spain Edited by Gerald F. Joyce, The Scripps Research Institute, La Jolla, CA, and approved December 3, 2009 (received for review November 3, 2009) A basic property of life is its capacity to experience Darwinian reactions that may deplete certain reactants (12) and dynamical evolution. The replicator concept is at the core of genetics-first aspects of autocatalytic cycles if they are assumed to coexist in theories of the origin of life, which suggest that self-replicating abstract space (13). Even if we ignore such hurdles, the key oligonucleotides or their similar ancestors may have been the first question still remains: Was a network of chemical reactions able to “living” systems and may have led to the evolution of an RNA increase in complexity and eventually undergo Darwinian selec- world. But problems with the nonenzymatic synthesis of biopol- tion as assumed by their advocates? A basic condition for any ymers and the origin of template replication have spurred the nascent Darwinian process in a population of self-reproducing alternative metabolism-first scenario, where self-reproducing systems is that they must have a sort of hereditary transmission and evolving proto-metabolic networks are assumed to have pre- which requires, in turn, becoming familiar with a lesser-known and dated self-replicating genes.
    [Show full text]
  • Micelles As Containers for Protocells
    Beitr¨age des Instituts f¨urUmweltsystemforschung der Universit¨atOsnabr¨uck Herausgeber: Prof. Dr. Michael Matthies Beitrag Nr. 33 Micelles as Containers for Protocells Harold Fellermann Dezember 2005 ISSN-Nr. 1433-3805 Beitr¨age des Instituts f¨urUmweltsystemforschung der Universit¨atOsnabr¨uck ISSN 1433-3805 Herausgeber Prof. Dr. Michael Matthies Universit¨at Osnabr¨uck Institut f¨urUmweltsystemforschung Artilleriestr. 34 D-49069 Osnabr¨uck Tel. 0541/969-2575 Fax. 0541/969-2599 E-Mail: [email protected] http://www.usf.uni-osnabrueck.de c USF – Institut f¨urUmweltsystemforschung, Universit¨atOsnabr¨uck Abstract We extended the mesoscopic simulation method Dissipative Particle Dynamics to incorporate chemical reactions on the basis of rate equations. With the extended model, we analyzed micellar systems in the context of early life and artificial proto- cells. We studied whether micelles have the potential to serve as embodiments for protocells, i. e. whether they are able to grow in size and divide into two daughter cells. It was found that the ability of micelles to grow depends on the pathway new surfactants are provided: micelles grow only if new surfactants are provided faster than the monomers dissociate from the assembly. We compared the growth scenario of two distinct scenarios that envision micelles as protocells: the origin of life theory “Lipid World” and the currently explored artificial protocell “Los Alamos bug”. We found that relaxation in micellar kinetics is likely to jeopardize the pathway proposed by the Lipid World theory, but not the “Los Alamos bug”. We finally developed a toy model that delineates the whole life-cycle of the “Los Alamos bug”.
    [Show full text]
  • PROSPECTS of a COMPUTATIONAL ORIGIN of LIFE ENDEAVOR 1. Introduction for Half a Century, Since the Breakthrough Work on Prebioti
    PROSPECTS OF A COMPUTATIONAL ORIGIN OF LIFE ENDEAVOR BARAK SHENHAV and DORON LANCET* Department of Molecular Genetics and the Crown Human Genome Center, the Weizmann Institute of Science, Rehovot 76100, Israel (* author for correspondence, e-mail: [email protected], phone: 972 8 9343683 or 972 8 9344121, fax: 972 8 9344487) (Received 22 January 2003; accepted in revised form 1 April 2003) Abstract. While the last century brought an exquisite understanding of the molecular basis of life, very little is known about the detailed chemical mechanisms that afforded the emergence of life on early earth. There is a broad agreement that the problem lies in the realm of chemistry, and likely resides in the formation and mutual interactions of carbon-based molecules in aqueous medium. Yet, present-day experimental approaches can only capture the synthesis and behavior of a few molecule types at a time. On the other hand, experimental simulations of prebiotic syntheses, as well as chemical analyses of carbonaceous meteorites, suggest that the early prebiotic hydrosphere contained many thousands of different compounds. The present paper explores the idea that given the limitations of test-tube approaches with regards to such a ‘random chemistry’ scenario, an alternative mode of analysis should be pursued. It is argued that as computational tools for the reconstruction of molecular interactions improve rapidly, it may soon become possible to perform adequate computer- based simulations of prebiotic evolution. We thus propose to launch a computational origin of life endeavor (http://ool.weizmann.ac.il/CORE), involving computer simulations of realistic complex prebiotic chemical networks.
    [Show full text]
  • States of Origin: Influences on Research Into the Origins of Life
    COPYRIGHT AND USE OF THIS THESIS This thesis must be used in accordance with the provisions of the Copyright Act 1968. Reproduction of material protected by copyright may be an infringement of copyright and copyright owners may be entitled to take legal action against persons who infringe their copyright. Section 51 (2) of the Copyright Act permits an authorized officer of a university library or archives to provide a copy (by communication or otherwise) of an unpublished thesis kept in the library or archives, to a person who satisfies the authorized officer that he or she requires the reproduction for the purposes of research or study. The Copyright Act grants the creator of a work a number of moral rights, specifically the right of attribution, the right against false attribution and the right of integrity. You may infringe the author’s moral rights if you: - fail to acknowledge the author of this thesis if you quote sections from the work - attribute this thesis to another author - subject this thesis to derogatory treatment which may prejudice the author’s reputation For further information contact the University’s Director of Copyright Services sydney.edu.au/copyright Influences on Research into the Origins of Life. Idan Ben-Barak Unit for the History and Philosophy of Science Faculty of Science The University of Sydney A thesis submitted to the University of Sydney as fulfilment of the requirements for the degree of Doctor of Philosophy 2014 Declaration I hereby declare that this submission is my own work and that, to the best of my knowledge and belief, it contains no material previously published or written by another person, nor material which to a substantial extent has been accepted for the award of any other degree or diploma of a University or other institute of higher learning.
    [Show full text]
  • Origin of Species Before Origin of Life: the Role of Speciation in Chemical Evolution
    life Review Origin of Species before Origin of Life: The Role of Speciation in Chemical Evolution Tony Z. Jia 1,2,* , Melina Caudan 1 and Irena Mamajanov 1,* 1 Earth-Life Science Institute, Tokyo Institute of Technology, 2-12-1-IE-1 Ookayama, Meguro-ku, Tokyo 152-8550, Japan; [email protected] 2 Blue Marble Space Institute of Science, 1001 4th Ave., Suite 3201, Seattle, WA 98154, USA * Correspondence: [email protected] (T.Z.J.); [email protected] (I.M.) Abstract: Speciation, an evolutionary process by which new species form, is ultimately responsible for the incredible biodiversity that we observe on Earth every day. Such biodiversity is one of the critical features which contributes to the survivability of biospheres and modern life. While speciation and biodiversity have been amply studied in organismic evolution and modern life, it has not yet been applied to a great extent to understanding the evolutionary dynamics of primitive life. In particular, one unanswered question is at what point in the history of life did speciation as a phenomenon emerge in the first place. Here, we discuss the mechanisms by which speciation could have occurred before the origins of life in the context of chemical evolution. Specifically, we discuss that primitive compartments formed before the emergence of the last universal common ancestor (LUCA) could have provided a mechanism by which primitive chemical systems underwent speciation. In particular, we introduce a variety of primitive compartment structures, and associated functions, that may have plausibly been present on early Earth, followed by examples of both discriminate and indiscriminate speciation affected by primitive modes of compartmentalization.
    [Show full text]
  • Prebiotic Network Evolution: Six Key Parameters
    Molecular BioSystems Prebiotic Network Evolution: Six Key Parameters Journal: Molecular BioSystems Manuscript ID MB-ART-09-2015-000593.R1 Article Type: Review Article Date Submitted by the Author: 02-Oct-2015 Complete List of Authors: Nghe, Philippe; CNRS ESPCI ParisTech, Laboratoire de Biochimie Hordijk, Wim; SmartAnalytiX.com, Kauffman, Stuart; Institute for Systems Biology, Walker, Sara; Arizona State University, School of Earth and Space Exploration and Beyond Center for Fundamental Concepts in Science Schmidt, Francis; University of Missouri, Biochemistry Kemble, Harry; CNRS ESPCI ParisTech, Yeates, Jessica; Portland State University, Chemistry Lehman, Niles; Portland State University, Department of Chemistry Page 1 of 36 Molecular BioSystems Prebiotic Network Evolution: Six Key Parameters Philippe Nghe 1†, Wim Hordijk 2†, Stuart A. Kauffman 3, Sara I. Walker 4, Francis J. Schmidt 5, Harry Kemble 1, Jessica A.M. Yeates 6, and Niles Lehman 6* 1 Laboratoire de Biochimie, CNRS - ESPCI ParisTech, France 2 SmartAnalytiX.com, Lausanne, Switzerland 3 Institute for Systems Biology, Seattle, WA 98109 USA 4 School of Earth and Space Exploration and Beyond Center for Fundamental Concepts in Science, Arizona State University, Tempe, AZ 85287 USA 5 Department of Biochemistry, University of Missouri, 117 Schweitzer Hall, Columbia, MO 65211, USA 6 Department of Chemistry, Portland State University, PO Box 751, Portland, OR 97207 USA †These authors contributed equally to this work. *To whom correspondence should be addressed at: Niles Lehman Department of Chemistry Portland State University PO Box 751 Portland, OR 97207 +1-503-725-8769 [email protected] 1 Molecular BioSystems Page 2 of 36 Abstract The origins of life likely required the cooperation among a set of molecular species interacting in a network.
    [Show full text]
  • Protobiotic Systems Chemistry Analyzed by Molecular Dynamics
    life Review Protobiotic Systems Chemistry Analyzed by Molecular Dynamics Amit Kahana and Doron Lancet * Dept. Molecular Genetics, The Weizmann Institute of Science, Rehovot 7610010, Israel; [email protected] * Correspondence: [email protected] Received: 16 April 2019; Accepted: 8 May 2019; Published: 10 May 2019 Abstract: Systems chemistry has been a key component of origin of life research, invoking models of life’s inception based on evolving molecular networks. One such model is the graded autocatalysis replication domain (GARD) formalism embodied in a lipid world scenario, which offers rigorous computer simulation based on defined chemical kinetics equations. GARD suggests that the first pre-RNA life-like entities could have been homeostatically-growing assemblies of amphiphiles, undergoing compositional replication and mutations, as well as rudimentary selection and evolution. Recent progress in molecular dynamics has provided an experimental tool to study complex biological phenomena such as protein folding, ligand-receptor interactions, and micellar formation, growth, and fission. The detailed molecular definition of GARD and its inter-molecular catalytic interactions make it highly compatible with molecular dynamics analyses. We present a roadmap for simulating GARD’s kinetic and thermodynamic behavior using various molecular dynamics methodologies. We review different approaches for testing the validity of the GARD model by following micellar accretion and fission events and examining compositional changes over time. Near-future computational advances could provide empirical delineation for further system complexification, from simple compositional non-covalent assemblies towards more life-like protocellular entities with covalent chemistry that underlies metabolism and genetic encoding. Keywords: systems chemistry; systems protobiology; molecular dynamics; gard; lipid world; micelle; origin of life 1.
    [Show full text]
  • Predicting Species Emergence in Simulated Complex Pre-Biotic Networks
    RESEARCH ARTICLE Predicting species emergence in simulated complex pre-biotic networks Omer Markovitch¤a¤b, Natalio Krasnogor* Interdisciplinary Computing and Complex Bio-Systems research group, School of Computing Science, Newcastle University, Newcastle upon Tyne, United-Kingdom ¤a Current Address: Center for Systems Chemistry, Stratingh Institute, University of Groningen, Groningen, The Netherlands ¤b Current Address: Blue Marble Space Institute of Science, Seattle, Washington, United States of America a1111111111 * [email protected] a1111111111 a1111111111 a1111111111 Abstract a1111111111 An intriguing question in evolution is what would happen if one could ªreplayº life's tape. Here, we explore the following hypothesis: when replaying the tape, the details (ªdecora- tionsº) of the outcomes would vary but certain ªinvariantsº might emerge across different OPEN ACCESS life-tapes sharing similar initial conditions. We use large-scale simulations of an in silico Citation: Markovitch O, Krasnogor N (2018) model of pre-biotic evolution called GARD (Graded Autocatalysis Replication Domain) to Predicting species emergence in simulated test this hypothesis. GARD models the temporal evolution of molecular assemblies, gov- complex pre-biotic networks. PLoS ONE 13(2): erned by a rates matrix (i.e. network) that biases different molecules' likelihood of joining or e0192871. https://doi.org/10.1371/journal. leaving a dynamically growing and splitting assembly. Previous studies have shown the pone.0192871 emergence of so called compotypes, i.e., species capable of replication and selection Editor: Ricard V. SoleÂ, Santa Fe Institute, SPAIN response. Here, we apply networks' science to ascertain the degree to which invariants Received: November 29, 2017 emerge across different life-tapes under GARD dynamics and whether one can predict Accepted: January 31, 2018 these invariant from the chemistry specification alone (i.e.
    [Show full text]
  • The Non-Genetic Dawn of Life a Review of the Contemporary Standing of the Metabolism-First Approach to the Origin of Life Enigma
    The non-genetic dawn of life a review of the contemporary standing of the metabolism-first approach to the origin of life enigma F.A.Bastiaan von Meijenfeldt Email: [email protected] 5 May 2013 A master thesis submitted to Universiteit Utrecht, Environmental Biology Examiner: Prof. dr. Paulien Hogeweg Theoretical Biology and Bioinformatics Department of Biology Faculty of Science Universiteit Utrecht Email: [email protected] Second reviewer: Dr. Berend Snel Theoretical Biology and Bioinformatics Department of Biology Faculty of Science Universiteit Utrecht Email: [email protected] The non-genetic dawn of life a review of the contemporary standing of the metabolism-first approach to the origin of life enigma F.A.Bastiaan von Meijenfeldt 5 May 2013 SUMMARY: the origin of life on early Earth is unresolved still. Two approaches to the problem dominate the debate. Replicator-first proponents search for the very first replicator molecule. A much envisioned molecule in this regard is RNA. Metabolism-first adherents deem it unlikely that a replicator molecule emerged de novo on early Earth, and instead search for a network of small molecules that could have seeded modern life. Three metabolism-first scenarios are discussed. The elaborate lipid world scenario proclaims the Darwinian selection of early metabolic networks. It is shown that even though persisting entities may have come into existence, heredity in ensemble replicators is probably too inaccurate to allow for their selection, and even if accurate enough inherently attractor based, posing severe problems for their subsequent evolution. The problems are exemplary for scenarios that pose the evolution of metabolic networks.
    [Show full text]