Mutations and Lethality in Simulated Prebiotic Networks

The Graded Autocatalysis Replication Domain (GARD) model describes an origin of life scenario which involves non-covalent compositional assemblies, made of monomeric mutually catalytic molecules. GARD constitutes an alternative to informational biopolymers as a mechanism of primordial inheritance. In the present work, we examined the effect of mutations, one of the most fundamental mechanisms for evolution, in the context of the networks of mutual interaction within GARD prebiotic assemblies. We performed a systematic analysis analogous to single and double gene deletions within GARD. While most deletions have only a small effect on both growth rate and molecular composition of the assemblies, ~10% of the deletions caused lethality, or sometimes showed enhanced fitness. Analysis of 14 different network properties on 2,000 different GARD networks indicated that lethality usually takes place when the deleted node has a high molecular count, or when it is a catalyst for such node. A correlation was also found between lethality and node degree centrality, similar to what is seen in real biological networks. Addressing double knockout mutations, our results demonstrate the occurrence of both synthetic lethality and extragenic suppression within GARD networks, and convey an attempt to correlate synthetic lethality to network node-pair properties. The analyses presented help establish GARD as a workable alternative prebiotic scenario, suggesting that life may have begun with large molecular networks of low fidelity, that later underwent evolutionary compaction and fidelity augmentation.     

Compositional Inheritance: Comparison of Self-assembly and Catalysis

Genetic inheritance in modern cells is due to template-directed replication of nucleic acids. However, the difficulty of prebiotic synthesis of long information-carrying polymers like RNA raises the question of whether some other form of heredity is possible without polymers. As an alternative, the lipid world theory has been proposed, which considers non-covalent assemblies of lipids, such as micelles and vesicles. Assemblies store information in the form of a non-random molecular composition, and this information is passed on when the assemblies divide, i.e. the assemblies show compositional inheritance. Here, we vary several important assumptions of previous lipid world models and show that compositional inheritance is relevant more generally than the context in which it was originally proposed. Our models assume that interaction occurs between nearest neighbour molecules only, and account for spatial segregation of molecules of different types within the assembly. We also draw a distinction between a self-assembly model, in which the composition is determined by mutually favourable interaction energies between the molecules, and a catalytic model, in which the composition is determined by mutually favourable catalysis. We show that compositional inheritance occurs in both models, although the self-assembly case seems more relevant if the molecules are simple lipids. In the case where the assemblies are composed of just two types of molecules, there is a strong analogy with the classic two-allele Moran model from population genetics. This highlights the parallel between compositional inheritance and genetic inheritance.     

Adduct size limits efficient and error-free bypass across bulky N2-guanine DNA lesions by human DNA polymerase eta

The N2 position of guanine (G) is one of the major sites for DNA modification by various carcinogens. Eight oligonucleotides with varying adduct bulk at guanine N2 were analyzed for catalytic efficiency and fidelity with human DNA polymerase (pol) eta, which is involved in translesion synthesis (TLS). Pol eta effectively bypassed N2-methyl(Me)G, N2-ethyl(Et)G, N2-isobutyl(Ib)G, N2-benzyl(Bz)G, and N2-CH2(2-naphthyl)G but was severely blocked at N2-CH2(9-anthracenyl)G (N2-AnthG) and N2-CH2(6-benzo[a]pyrenyl)G (N2-BPG). Steady-state kinetic analysis showed proportional decreases of kcat/Km in dCTP insertion opposite N2-AnthG and N2-BPG (73 and 320-fold) and also kcat/Km in next-base extension from a C paired with each adduct (15 and 51-fold relative to G). Frequencies of dATP misinsertion and extension beyond mispairs were also proportionally increased (70 and 450-fold; 12 and 44-fold) with N2-AnthG and N2-BPG, indicating the effect of adduct bulk on blocking and misincorporation in TLS by pol eta. N2-AnthG and N2-BPG also greatly decreased the pre-steady-state kinetic burst rate (25 and 125-fold) compared to unmodified G. N2-AnthG decreased dCTP binding affinity (2.6-fold) and increased DNA substrate binding affinity. These results and the small kinetic thio effects (S(p)-dCTPalphaS) suggest that the early steps, possibly conformational change, are interfered with by the bulky adducts. In contrast, human pol delta bypassed adducts effectively up to N2-EtG but was strongly blocked by N2-IbG and larger adducts. We conclude that TLS DNA polymerases may be required for the efficient bypass of pol delta-blocking N2-G adducts bulkier than N2-EtG in human cells, and the bulk size can be a major factor for efficient and error-free bypass at these adducts by TLS DNA polymerases.     

Graded Autocatalysis Replication Domain (GARD): Kinetic Analysis of Self-Replication in Mutually Catalytic Sets

A Graded Autocatalysis Replication Domain (GARD) model is proposed, which provides a rigorous kinetic analysis of simple chemical sets that manifest mutual catalysis. It is shown that catalytic closure can sustain self-replication up to a critical dilution rate, _c, related to the graded extent of mutual catalysis. We explore the behavior of vesicles containing GARD species whose mutual catalysis is governed by a previously published statistical distribution. In the population thus generated, some GARD vesicles display a significantly higher replication efficiency than most others. GARD thus represents a simple model for primordial chemical selection of mutually catalytic sets.     

A formalisation of the neural assembly concept 1. Constraints on neural assembly size

Hebb proposed the concept of a neural assembly distributed across cortical tissue as a model for representation of information in the cerebral cortex. Later developments of the concept highlight the need for overlapping membership between independent assemblies, and the spread of activity throughout the assembly once it is activated above a critical level (ignition). Formalisation of the neural assembly concept, especially in relation to quantitative data from the real cortex, is at a very early stage. We consider two constraints on neural assembly size: (1) if a neural assembly is too small the fraction of its neurons that need to be active to ignite the whole assembly becomes unrealistically large; (2) if assemblies in a block of cortical tissue become too large then the block becomes ‘unsafe’, that is, unwanted spread from an active assembly to overlapping ones becomes inevitable. We consider variations in three parameters: neuronal firing threshold; connection density; and the total number of assemblies stored in the block of cortical tissue. Given biologically plausible values for these parameters we estimate maximum assembly size compatible with ignitability of individual assemblies, low probability of unwanted spread to overlapping assemblies, and safe operation of the block as a whole. Received: 7 March 1997 / Accepted in revised form: 1 July 1997     

Kinetic evidence for inefficient and error-prone bypass across bulky N2-guanine DNA adducts by human DNA polymerase iota

DNA polymerase (pol) iota has been proposed to be involved in translesion synthesis past minor groove DNA adducts via Hoogsteen base pairing. The N2 position of G, located in minor groove side of duplex DNA, is a major site for DNA modification by various carcinogens. Oligonucleotides with varying adduct size at G N2 were analyzed for bypass ability and fidelity with human pol iota. Pol iota effectively bypassed N2-methyl (Me)G and N2-ethyl(Et)G, partially bypassed N2-isobutyl(Ib)G and N2-benzylG, and was blocked at N2-CH2(2-naphthyl)G (N2-NaphG), N2-CH2(9-anthracenyl)G (N2-AnthG), and N2-CH2(6-benzo[a]pyrenyl)G. Steady-state kinetic analysis showed decreases of kcat/Km for dCTP insertion opposite N2-G adducts according to size, with a maximal decrease opposite N2-AnthG (61-fold). dTTP misinsertion frequency opposite template G was increased 3-11-fold opposite adducts (highest with N2-NaphG), indicating the additive effect of bulk (or possibly hydrophobicity) on T misincorporation. N2-IbG, N2-NaphG, and N2-AnthG also decreased the pre-steady-state kinetic burst rate compared with unmodified G. High kinetic thio effects (S(p)-2'-deoxycytidine 5'-O-(1-thiotriphosphate)) opposite N2-EtG and N2-AnthG (but not G) suggest that the chemistry step is largely interfered with by adducts. Severe inhibition of polymerization opposite N2,N2-diMeG compared with N2-EtG by pol eta but not by pol iota is consistent with Hoogsteen base pairing by pol iota. Thus, polymerization by pol iota is severely inhibited by a bulky group at G N2 despite an advantageous mode of Hoogsteen base pairing; pol iota may play a limited role in translesion synthesis on bulky N2-G adducts in cells.     

Coevolution of compositional protocells and their environment

The coevolution of environment and living organisms is well known in nature. Here, it is suggested that similar processes can take place before the onset of life, where protocellular entities, rather than full-fledged living systems, coevolve along with their surroundings. Specifically, it is suggested that the chemical composition of the environment may have governed the chemical repertoire generated within molecular assemblies, compositional protocells, while compounds generated within these protocells altered the chemical composition of the environment. We present an extension of the graded autocatalysis replication domain (GARD) model--the environment exchange polymer GARD (EE-GARD) model. In the new model, molecules, which are formed in a protocellular assembly, may be exported to the environment that surrounds the protocell. Computer simulations of the model using an infinite-sized environment showed that EE-GARD assemblies may assume several distinct quasi-stationary compositions (composomes), similar to the observations in previous variants of the GARD model. A statistical analysis suggested that the repertoire of composomes manifested by the assemblies is independent of time. In simulations with a finite environment, this was not the case. Composomes, which were frequent in the early stages of the simulation disappeared, while others emerged. The change in the frequencies of composomes was found to be correlated with changes induced on the environment by the assembly. The EE-GARD model is the first GARD model to portray a possible time evolution of the composomes repertoire.     

Comparison of phenotypic and genetic clone delineation in quaking aspen, <Emphasis Type="Italic">Populus tremuloides</Emphasis>

Key message

Clonal delineation at nuclear microsatellites and phenotypic traits showed high correspondence and revealed an important role of both sexual and clonal reproduction for stand genetic structure.

Abstract

Quaking aspen (Populus tremuloides Michx.) grows throughout the northern and central portions of North America. Reproduction occurs both sexually via seeds and clonally from root suckers. Clonal delineation using morphological/phenological traits, and more recently, highly variable nuclear microsatellites have shown considerable variation in the size of clonal assemblies, and the relative importance of sexual versus clonal reproduction across the species range. In order to provide reliable estimates of genet size (N/G; ramets per sampled genet) and genotypic diversity (G/N; genets/ramets), and to compare genetic and phenotypic clone delineation, we characterized 181 sampled stems (ramets) at seven nuclear microsatellites, and morphological and phenological traits from six clones (genet size ≥11). Genotypic diversity was moderate (G/N = 0.18) and within the range reported in other studies across North America. Multivariate statistics revealed a high correspondence between genetic and phenotypic clone delineation, both with and without predefined genetic groups (94.2 %, 81.7 %). Moderate average genet size (5.6 ramets per genet) and the occurrence of genetically distinct single-ramet genets surrounded by larger genets suggested intermediate levels of sexual reproduction contributing to the genetic structure of this stand. Significant differences among genets were found for phenological and morphological traits such as bark thickness and leaf shape. However, most clones showed no significant differences in diameter growth which was likely caused by poor drainage in this high clay soil that inhibited the expression of genetic differences in growth.

     

LC-MS/MS identification and yeast polymerase eta bypass of a novel gamma-irradiation-induced intrastrand cross-link lesion G[8-5]C

Reactive oxygen species can give rise to intrastrand cross-link lesions, where two neighboring nucleobases are covalently bonded. Here, we employed LC-MS/MS and demonstrated for the first time that gamma irradiation of a synthetic duplex oligodeoxyribonucleotide can give rise to an intrastrand cross-link lesion G[8-5]C, where the C8 carbon atom of guanine and the C5 carbon atom of its 3'-neighboring cytosine are covalently bonded. We also carried out in vitro replication studies of a substrate containing a site-specifically incorporated G[8-5]C, and our results showed that yeast Saccharomyces cerevisiae DNA polymerase eta (pol eta) was able to replicate past the cross-link lesion. Steady-state kinetic analyses for nucleotide incorporation by pol eta showed that the 3'-cytosine moiety of the cross-link did not significantly affect either the efficiency or the fidelity of nucleotide incorporation. The 5' guanine portion of the cross-link lesion, however, markedly reduced both the efficiency and the fidelity of nucleotide incorporation; the insertion of dGMP or dAMP was slightly favored over the insertion of the correct nucleotide, dCMP, which was in turn favored over the insertion of dTMP. The above results support that the oxidative cross-link lesion, if not repaired, can be mutagenic.     

Getting from an RNA world to modern cells just got a little easier

Our understanding of the early steps in the evolution of life is hampered by a Catch-22: Darwinian selection leading to longer genomes requires as prerequisite increased replicative fidelity. Yet a genome at capacity cannot increase in size; it will be catastrophically mutated out of existence if fidelity has not already increased. Traditionally the problem has been considered for genotypes but can be down-sized if multiple genotypes specify the same phenotype. Kun and colleagues put empirical meat on theoretical bone by analysing ribozyme mutagenesis data, concluding that modest replication fidelities could permit a primordial genome with up to 100 genes.     

Structural basis of cellulosome efficiency explored by small angle X-ray scattering

Cellulose, the main structural component of plant cell walls, is the most abundant carbohydrate polymer in nature. To break down plant cell walls, anaerobic microorganisms have evolved a large extracellular enzyme complex termed cellulosome. This megadalton catalytic machinery organizes an enzymatic assembly, tenaciously bound to a scaffolding protein via specialized intermodular "cohesin-dockerin" interactions that serve to enhance synergistic activity among the different catalytic subunits. Here, we report the solution structure properties of cellulosome-like assemblies analyzed by small angle x-ray scattering and molecular dynamics. The atomic models, generated by our strategy for the free chimeric scaffoldin and for binary and ternary complexes, reveal the existence of various conformations due to intrinsic structural flexibility with no, or only coincidental, inter-cohesin interactions. These results provide primary evidence concerning the mechanisms by which these protein assemblies attain their remarkable synergy. The data suggest that the motional freedom of the scaffoldin allows precise positioning of the complexed enzymes according to the topography of the substrate, whereas short-scale motions permitted by residual flexibility of the enzyme linkers allow "fine-tuning" of individual catalytic domains.     

Biomolecular analysis of elastic fibre molecules

Elastic fibres are macromolecular extracellular matrix assemblies that endow dynamic connective tissues such as arteries, lungs and skin with the property of elastic recoil. Here, we describe how we have purified elastic fibre molecules and then analysed them using a range of biochemical and biomolecular approaches. Such approaches have provided powerful insights into the complex hierarchical processes of extracellular matrix assembly. We outline molecular interaction and kinetics assays using Biacore, biophysical approaches such as multi-angle laser light scattering and analytical ultracentrifugation which provide information on molecular and macromolecular shape and mass in solution, the visualisation of molecules and assemblies using microscopy approaches such as atomic force microscopy and environmental scanning electron microscopy, and compositional analysis of macromolecular complexes using mass spectrometry. Data from these in vitro analytical approaches can be combined to develop powerful new models of elastic fibre assembly.     

Loss of DNA polymerase beta stacking interactions with templating purines, but not pyrimidines, alters catalytic efficiency and fidelity.

Structures of DNA polymerases bound with DNA reveal that the 5'-trajectory of the template strand is dramatically altered as it exits the polymerase active site. This distortion provides the polymerase access to the nascent base pair to interrogate proper Watson-Crick geometry. Upon binding a correct deoxynucleoside triphosphate, alpha-helix N of DNA polymerase beta is observed to form one face of the binding pocket for the new base pair. Asp-276 and Lys-280 stack with the bases of the incoming nucleotide and template, respectively. To determine the role of Lys-280, site-directed mutants were constructed at this position, and the proteins were expressed and purified, and their catalytic efficiency and fidelity were assessed. The catalytic efficiency for single-nucleotide gap filling with the glycine mutant (K280G) was strongly diminished relative to wild type for templating purines (>15-fold) due to a decreased binding affinity for the incoming nucleotide. In contrast, catalytic efficiency was hardly affected by glycine substitution for templating pyrimidines (<4-fold). The fidelity of the glycine mutant was identical to the wild type enzyme for misinsertion opposite a template thymidine, whereas the fidelity of misinsertion opposite a template guanine was modestly altered. The nature of the Lys-280 side-chain substitution for thymidine triphosphate insertion (templating adenine) indicates that Lys-280 "stabilizes" templating purines through van der Waals interactions.     

Biochemical basis of genotoxicity of heterocyclic arylamine food mutagens: Human DNA polymerase eta selectively produces a two-base deletion in copying the N2-guanyl adduct of 2-amino-3-methylimidazo[4,5-f]quinoline but not the C8 adduct at the NarI G3 site

Heterocyclic arylamines are highly mutagenic and cause tumors in animal models. The mutagenicity is attributed to the C8- and N2-G adducts, the latter of which accumulates due to slower repair. The C8- and N 2-G adducts derived from 2-amino-3-methylimidazo[4,5-f]quinoline (IQ) were placed at the G1 and G3 sites of the NarI sequence, in which the G3 site is an established hot spot for frameshift mutation with the model arylamine derivative 2-acetylaminofluorene but G1 is not. Human DNA polymerase (pol) eta extended primers beyond template G-IQ adducts better than did pol kappa and much better than pol iota or delta. In 1-base incorporation studies, pol eta inserted C and A, pol iota inserted T, and pol kappa inserted G. Steady-state kinetic parameters were measured for these dNTPs opposite the C8- and N 2-IQ adducts at both sites, being most favorable for pol eta. Mass spectrometry of pol eta extension products revealed a single major product in each of four cases; with the G1 and G3 C8-IQ adducts, incorporation was largely error-free. With the G3 N 2-IQ adduct, a -2 deletion occurred at the site of the adduct. With the G1 N 2-IQ adduct, the product was error-free at the site opposite the base and then stalled. Thus, the pol eta products yielded frame-shifts with the N 2 but not the C8 IQ adducts. We show a role for pol eta and the complexity of different chemical adducts of IQ, DNA position, and DNA polymerases.     

Effects of accessory proteins on the bypass of a cis-syn thymine-thymine dimer by Saccharomyces cerevisiae DNA polymerase eta

Among several hypotheses to explain how translesion synthesis (TLS) by DNA polymerase eta (pol eta) suppresses ultraviolet light-induced mutagenesis in vivo despite the fact that pol eta copies DNA with low fidelity, here we test whether replication accessory proteins enhance the fidelity of TLS by pol eta. We first show that the single-stranded DNA binding protein RPA, the sliding clamp PCNA, and the clamp loader RFC slightly increase the processivity of yeast pol eta and its ability to recycle to new template primers. However, these increases are small, and they are similar when copying an undamaged template and a template containing a cis-syn TT dimer. Consequently, the accessory proteins do not strongly stimulate the already robust TT dimer bypass efficiency of pol eta. We then perform a comprehensive analysis of yeast pol eta fidelity. We show that it is much less accurate than other yeast DNA polymerases and that the accessory proteins have little effect on fidelity when copying undamaged templates or when bypassing a TT dimer. Thus, although accessory proteins clearly participate in pol eta functions in vivo, they do not appear to help suppress UV mutagenesis by improving pol eta bypass fidelity per se.     

A structural solution for the DNA polymerase lambda-dependent repair of DNA gaps with minimal homology

Human DNA polymerase lambda (Pol lambda) is a family X member with low frameshift fidelity that has been suggested to perform gap-filling DNA synthesis during base excision repair and during repair of broken ends with limited homology. Here, we present a 2.1 A crystal structure of the catalytic core of Pol lambda in complex with DNA containing a two nucleotide gap. Pol lambda makes limited contacts with the template strand at the polymerase active site, and superimposition with Pol beta in a ternary complex suggests a shift in the position of the DNA at the active site that is reminiscent of a deletion intermediate. Surprisingly, Pol lambda can adopt a closed conformation, even in the absence of dNTP binding. These observations have implications for the catalytic mechanism and putative DNA repair functions of Pol lambda.