A model of DNA sequence evolution

Statistical studies of gene populations on the purine/pyrimidine alphabet have shown that the mean occurrence probability of thei-motif YRY(N) _i YRY (R=purine, Y=pyrimidine, N=R or Y) is not uniform by varyingi in the range [1,99], but presents a maximum ati=6 in the following populations: protein coding genes of eukaryotes, prokaryotes, chloroplasts and mitrochondria, and also viral introns, ribosomal RNA genes and transfer RNA genes (Arquès and Michel, 1987b,J. theor. Biol. 128, 457–461). From the “universality” of this observation, we suggested that the oligonucleotide YRY(N)₆ is a primitive one and that it has a central function in DNA sequence evolution (Arquès and Michel, 1987b,J. theor. Biol. 128, 457–461). Following this idea, we introduce a concept of a model of DNA sequence evolution which will be validated according to a shema presented in three parts. In the first part, using the last version of the gene database, the YRY(N)₆YRY preferential occurrence (maximum ati=6) is confirmed for the populations mentioned above and is extended to some newly analysed populations: chloroplast introns, chloroplast 5′ regions, mitochondrial 5′ regions and small nuclear RNA genes. On the other hand, the YRY(N)₆YRY preferential occurrence and periodicities are used in order to classify 18 gene populations. In the second part, we will demonstrate that several statistical features characterizing different gene populations (in particular the YRY(N)₆YRY preferential occurrence and the periodicities) can be retrieved from a simple Markov model based on the mixing of the two oligonucleotides YRY(N)₆ and YRY(N)₃ and based on the percentages of RYR and YRY in the unspecified trinucleotides (N)₃ of YRY(N)₆ and YRY(N)₃. Several properties are identified and prove in particular that the oligonucleotide mixing is an independent process and that several different features are functions of a unique parameter. In the third part, the return of the model to the reality shows a strong correlation between reality and simulation concerning the presence of large alternating purine/pyrimidine stretches and of periodicities. It also contributes to a greater understanding of biological reality, e.g. the presence or the absence of large alternating purine/pyrimidine stretches can be explained as being a simple consequence of the mixing of two particular oligonucleotides. Finally, we believe that such an approach is the first step toward a unified model of DNA sequence evolution allowing the molecular understanding of both the origin of life and the actual biological reality.     

A purine-pyrimidine motif verifying an identical presence in almost all gene taxonomic groups

A statistical parameter identifies, with a high degree of significance, a motif which is present in protein-coding sequences of eukaryotes, prokaryotes, chloroplasts, mitochondria, viral introns, ribosomal RNA genes, and transfer RNA genes. The random probability of occurrence of such a situation is 10(-12). This motif has the following properties: (i) its significant presence in almost all present-day genes explains why it can be considered as primitive oligonucleotide, (ii) its nucleotide order is: YRY (N)6YRY, R being a purine base, Y a pyrimidine one and N any base, (iii) its length and its terminal trinucleotides YRY suggest a primordial function related to the spatial structure of the DNA sequences. This motif is found in some viral protein-coding genes, but not in eukaryotic introns.     

On the prevalence of certain codons (“RNY”) in genes for proteins

J.C. Shepherd notes that codons of the type RNY (R = purine, N = any nucleotide base, Y = pyrimidine) predominate over RNR in the genes for proteins. He has hypothesized that RNY codons are the relics of “a primitive code” composed of repeating RNY triplets. He found that RNY codons predominated in fourfold RNN codon sets (family boxes). These family boxes code for valine, threonine, alanine, and glycine. We argue that the proposed “comma-less” code composed of RNY never existed, and that, in any case, survival of such a code would have long since been erased by mutations. The excess of RNY codons in family boxes is probably attributable to preference for the corresponding tRNAs.     

Significance of the purine-pyrimidine motif present in most gene groups

The probability of the sequence YRY(Ni)YRY occurring most frequently with the same i value in seven out of nine gene classes is reassessed and found to be about 1.3 X 10(-8), more than 4000 times greater than the value calculated by Arquès & Michel (1987), but still much too small for chance to be a reasonable explanation for the observation. Even if the sequence YRYNNNNNNYRY were very frequent in the most primitive genes, it would not have survived in a recognizable form to the present day if it were selectively neutral. However, if it is selectively favoured one would expect it to exist regardless of whether it was present in primitive genes.     

Hypothetical stages in the mutational evolution of the present-day genetic code from its RNY,comma-free ancestor

The primitive comma-free genetic code may have had 16 triplets of the form RNY, where R = purine, N = purine or pyrimidine, and Y = pyrimidine, specifying eight (present-day) amino acids. Calculations reveal that in this primitive code all transition changes (A?G, C?U) are either silent or missense i.e. result in the same or another one of these particular eight amino acids. There are no single transitions to non-RNY codons. Single transversions in the primitive codons can, individually, generate new (present-day) codons for four or eight amino acids. Present-day glutamine, tryptophan and stop (UGA, UAA, UAG) codons cannot be so derived., by single transversions, from any of the eight primitive codons. The modern initiation codons, AUG and GUG, can however be generated by both C → G and U → G single transversions in primitive codons. Overall, a total of 32 modern sense codons, not represented in the primitive RNY code, can be derived from this code by single transversions. Many modern codons, including all those not generated by single transversions in the primitive code, can also be produced by either of the two types of frameshift possible in runs of U- or C-rich primitive codons. Present-day stop codons are generated by +1 (-2) type frameshifts in U-rich primitive runs; AUG and GUG initiation codons are produced by the other type, +2 (-1), frameshifts in U-rich runs.     

The effects of mispair and nonpair correction in hybrid DNA on base ratios (G + C content) and total amounts of DNA

Base ratios and total DNA amounts can vary substantially between and within higher taxa and genera, and even within species. Gene conversion is one of several mechanisms that could cause such changes. For base substitutions, disparity in conversion direction is accompanied by an equivalent disparity in base ratio at the heterozygous site. Disparity in the direction of gene conversion at meiosis is common and can be extreme. For transitions (which give purine [R]/pyrimidine [Y] mispairs) and for transversions giving unlike R/R and Y/Y mispairs in hybrid DNA, this disparity could give slow but systematic changes in G + C percentage. For transversions giving like R/R and Y/Y mispairs, it could change AT/TA and CG/GC ratios. From the extent of correction direction disparity, one can deduce properties of repair enzymes, such as the ability (1) to excise preferentially the purine from one mispair and the pyrimidine from the other for two different R/Y mispairs from a single heterozygous site and (2) to excise one base preferentially from unlike R/R or Y/Y mispairs. Frame-shifts usually show strong disparity in conversion direction, with preferential cutting of the nonlooped or the looped-out strand of the nonpair in heterozygous h-DNA. The opposite directions of disparity for frame-shifts and their intragenic suppressors as Ascobolus suggest that repair enzymes have a strong, systematic bias as to which strand is cut. The conversion spectra of mutations induced with different mutagens suggest that the nonlooped strand is preferentially cut, so that base additions generally convert to mutant and deletions generally convert to wild-type forms. Especially in nonfunctional or noncoding DNA, this could cause a general increase in DNA amounts. Conversion disparity, selection, mutation, and other processes interact, affecting rates of change in base ratios and total DNA.      

A study of the purine/pyrimidine codon occurrence with a reduced centered variable and an evaluation compared to the frequency statistic

With the three-letter alphabet [R,Y,N] (R = purine, Y = pyrimidine, N = R or Y), there are 26 codons (NNN being excluded): RNN,...,NNY (six codons at two unspecified bases N), RRN,...,NYY (12 codons at one unspecified base N), RRR,...,YYY (eight specified codons). A statistical methodology that uses the codon frequency and a reduced centered variable leads to similar results for a codon occurrence study, regardless of gene function and regardless of a particular protein coding gene taxonomic population. Therefore, this variable can be considered a new codon usage index, whose use removes certain nonsignificant results found with the frequency statistic. This methodology identifies the common and rare codons (i.e., the codons having the highest and lowest occurrence) and leads to a model of codon evolution at three successive states: RNN, then RNY, and finally RYY. Some biological relations between this model and the YRY(N)6YRY preferential occurrence are also presented.     

Novel third-letter bias in Escherichia coli codons revealed by rigorous treatment of coding constraints

A novel bias in codon third-letter usage was found in Escherichia coli genes with low fractions of "optimal codons", by comparing intact sequences with control random sequences. Third-letter usage has been found to be biased according to preference in codon usage and to doublet preference from the following first letter. The present study examines third-letter usage in the context of the nucleotide sequence when these preferences are considered. In order to exclude any influence by these factors, the random sequences were generated such that the amino acid sequence, codon usage, and the doublet frequency in each gene were all preserved. Comparison of intact sequences with these randomly generated sequences reveals that third letters of codons show a strong preference for the purine/pyrimidine pattern of the next codons: purine (R) is preferred to pyrimidine (Y) at the third site when followed by an R-Y-R codon, and pyrimidine is preferred when followed by an R-R-Y, an R-Y-Y or a Y-R-Y codon. This bias is probably related to interactions of tRNA molecules in the ribosome.     

Nuclear magnetic resonance structural studies of intramolecular purine.purine.pyrimidine DNA triplexes in solution. Base triple pairing alignments and strand direction

Recently, P.A. Beal and P.B. Dervan, expanding on earlier observations by others, have established the formation of purine.purine.pyrimidine triple helices stabilized by G.GC, A.AT and T.AT base triples where the purine-rich third strand was positioned in the major groove of the Watson-Crick duplex and anti-parallel to its purine strand. The present nuclear magnetic resonance (n.m.r.) study characterizes the base triple pairing alignments and strand direction in a 31-mer deoxyoligonucleotide that intramolecularly folds to generate a 7-mer (R/Y-)n.(R+)n(Y-)n triplex with the strands linked by two T5 loops and stabilized by potential T.AT and G.GC base triples. (R and Y stand for purine and pyrimidine, respectively, while the signs establish the strand direction.) This intramolecular triplex gives well-resolved exchangeable and non-exchangeable proton spectra with Li+ as counterion in aqueous solution. These studies establish that the T1 to C7 pyrimidine and the G8 to A14 purine strands are anti-parallel to each other and align through Watson-Crick A.T and G.C pair formation. The T15 to G21 purine-rich third strand is positioned in the major groove of this duplex and pairs through Hoogsteen alignment with the purine strand to generate T.AT and G.GC triples. Several lines of evidence establish that the thymidine and guanosine bases in the T15 to G21 purine-rich third strand adopt anti glycosidic torsion angles under conditions where this strand is aligned anti-parallel to the G8 to A14 purine strand. We have also recorded imino proton n.m.r. spectra for an (R-)n.(R+)n(Y-)n triplex stabilized by G.GC and A.AT triples through intramolecular folding of a related 31-mer deoxyoligonucleotide with Li+ as counterion. The intramolecular purine.purine.pyrimidine triplexes containing unprotonated G.GC, A.AT and T.AT triples are stable at basic pH in contrast to pyrimidine.purine.pyrimidine triplexes containing protonated C+.GC and T.AT triples, which are only stable at acidic pH.     

Amplification of the sequences displaying the pattern RNY in the RNA world: the translation --> translation/replication hypothesis

Based on previous considerations published in J. theor. Biol., new analyses of the organization of the genetic system are reported in this paper. We show that theoretical considerations about the order observed in the genetic code table support the idea of a primitive self-aminoacylation process achieved by primordial tRNAs. The physico-chemical constraints connected with this process may explain why a primitive genetic system predominantly uses sequences with the codonic pattern RNN (R=purine; Y=pyrimidine; N=any of the four bases) to polymerize the amino acids into peptides through translation. These considerations lead us to propose the Translation --> Translation/Replication hypothesis, which may explain why only RNA sequences with the pattern RNY, instead of less restrictive RNN, are susceptible to amplification. Using these ideas, supported by properties of symmetry, features of the genetic code may be connected with the replication of specific RNA sequences in the RNA world.     

A simulation of the genetic periodicities modulo 2 and 3 with processes of nucleotide insertions and deletions.

Recently, a new genetic process termed RNA editing has been identified showing insertions and deletions of nucleotides in particular RNA molecules. On the other hand, there are a few non-random statistical properties in genes: in particular, the periodicity modulo 3 (P3) associated with an open reading frame, the periodicity modulo 2 (P2) associated with alternating purine/pyrimidine stretches, the YRY(N)6YRY preferential occurrence (R = purine = adenine or guanine, Y = pyrimidine = cytosine or thymine, N = R or Y) representing a "code" of the DNA helix pitch, etc. The problem investigated here is whether a process of the type RNA editing can lead to the non-random statistical properties commonly observed in genes. This paper will show in particular that: The process of insertions and deletions of mononucleotides in the initial sequence [YRY(N)3]* [series of YRY(N)3] can lead to the periodicity modulo 2 (P2). The process of insertions and deletions of trinucleotides in the initial sequence [YRY(N)6]* [series of YRY(N)6] can lead to the periodicity modulo 3 (P3) and the YRY(N)6YRY preferential occurrence. Furthermore, these two processes lead to a strong correlation with the reality, namely the mononucleotide insertion/deletion process, with the 5' eukaryotic regions and the trinucleotide insertion/deletion process, with the eukaryotic protein coding genes.     

Molecular dynamics and mechanics calculations on a DNA duplex with A(+)-C,G-T and T-C mispairs

Despite major advances in characterizing purine(R)-purine(R), purine(R)-pyrimidine(Y) and pyrimidine(Y)-pyrimidine(Y) mismatches in DNA, there have not been any structural studies on a synthetic DNA duplex containing several different mispairs. Here, using NMR restrained molecular mechanics and dynamics simulations we have structurally characterized a 12 nucleotide long antiparallel DNA duplex with three different mispairs, namely A+-C, G-T and T-C. Our results show that the overall conformation of the antiparallel DNA duplex is B-DNA-like with slight structural distortions at or near the mispairs' sites. All these mispairs are properly stacked with their flanking base pairs. Each mispair is stabilized by two hydrogen bonds and the decreasing order of the hydrogen-bonding interactions is G-T>T-C>A+-C. G-T mispair has smaller configurational space while the structure is slightly bent at A+-C mispair's site. Overall, this study is the first ever structural characterization of a DNA duplex with three different mismatched base pairs and throws light upon the local conformations of the three mispairs present in the DNA duplex.     

Universal full-length nucleosome mapping sequence probe

For the computational sequence-directed mapping of the nucleosomes, the knowledge of the nucleosome positioning motifs – 10–11 base long sequences – and respective matrices of bendability, is not sufficient, since there is no justified way to fuse these motifs in one continuous nucleosome DNA sequence. Discovery of the strong nucleosome (SN) DNA sequences, with visible sequence periodicity allows derivation of the full-length nucleosome DNA bendability pattern as matrix or consensus sequence. The SN sequences of three species (A. thaliana, C. elegans, and H. sapiens) are aligned (512 sequences for each species), and long (115 dinucleotides) matrices of bendability derived for the species. The matrices have strong common property – alternation of runs of purine–purine (RR) and pyrimidine–pyrimidine (YY) dinucleotides, with average period 10.4 bases. On this basis the universal [R,Y] consensus of the nucleosome DNA sequence is derived, with exactly defined positions of respective penta- and hexamers RRRRR, RRRRRR, YYYYY, and YYYYYY.     

Molecular mechanics calculations on a triple stranded DNA involving C+.G-T and T.A+-C mismatched base triples

We have carried out molecular modeling of a triple stranded pyrimidine(Y). purine(R): pyrimidine(Y) (where ':' refers to Watson-Crick and '.' to Hoogsteen bonding) DNA, formed by a homopurine (d-TGAGGAAAGAAGGT) and homo-pyrimidine (d-CTCCTTTCTTCC). Molecular mechanics calculations using NMR constraints have provided a detailed three dimensional structure of the triplex. The entire stretches of purine and the pyrimidine nucleotides have a conformation close to B-DNA. The three strands are held by the canonical C+.G:C and T.A:T hydrogen bonds. The structure also contains two mismatch C+.G-T and T.A+-C base triples which have been characterized for the first time. In the A+-C base-pair of the T.A+-C triple, both hydrogen donors are situated on the purine (A+(1N) and A+(6N)). We observe a unique hydrogen bonding interaction scheme in case of C+.G-T where one acceptor, G(60), is bonded to three donors (C+(3NH), C+(4NH2) and T(3NH)). Though the C+.G-T base triple is less stable than C+.G:C, it is significantly more stable than T.A:T. On the other hand, T.A+-C is as stable as the T.A:T base triad.     

Properties of Oligonucleotides Containing the Transition Bases P and K

Abstract

The properties of the degenerate nucleosides dP and dK, in templates and primers were determined. dP was copied as either pyrimidine, dK as either purine. In primers, an equimolar mixture of the two nucleosides functioned as a universal base equivalent in both sequencing and the polymerase chain reactions. Cloning of DNA containing dP or dK produced transition mutations in vivo.     

Single-Stranded DNA and RNA Targeted Triplex-Formation: UV,CD and Molecular Modeling Studies of Foldback Triplexes Containing Different RNA, 2′-OMe-RNA and DNA Strand Combinations

Abstract

We studied the influence of different 2′-OMe-RNA and DNA strand combinations on single strand targeted foldback triplex formation in the Py.Pu:Py motif using ultraviolet (UV) and circular dichroism (CD) spectroscopy, and molecular modeling. The study of eight combinations of triplexes (D D:D, R* D:D, D D:R*, R* D:R*, D R:D, R* R:D, DR:R*, and R*-R:R*; where the first, middle, and last letters stand for the Hoogsteen Pyrimidine, Watson-Crick [WC] purine and WC pyrimidine strands, respectively, and D, R and R* stand for DNA, RNA and 2′-OMe-RNA strands, respectively) indicate more stable foldback triplex formation with a DNA purine strand than with an RNA purine strand. Of the four possible WC duplexes with RNA/DNA combinations, the duplex with a DNA purine strand and a 2′-O-Me-RNA pyrimidine strand forms the most thermally stable triplex, although its thermal stability is the lowest of all four duplexes. Irrespective of the duplex combination, a 2′-OMe-RNA Hoogsteen pyrimidine strand forms a stable foldback triplex over a DNA Hoogsteen pyrimidine strand confirming the earlier reports with conventional and circular triplexes. The CD studies suggest a B-type conformation for an all DNA homo-foldback triplex (D.D.D), while hetero-foldback triplex spectra suggest intermediate conformation to both Atype and B-type structures. A novel molecular modeling study has been carried out to understand the stereochemical feasibility of all the combinations of foldback triplexes using a geometric approach. The new approach allows use of different combinations of chain geometries depending on the nature of the chain (RNA vs. DNA).     

Synthesis and Biological Evaluation of Pyrimidine and Purine α-L-2′,3′-Dideoxy Nucleosides

Abstract

A series of α-L-2′,3′-dideoxy nucleosides was prepared as potential antiviral agents. The pyrimidine nucleosides were prepared by standard Vorbrüggen coupling reactions. The purine analogues were prepared by enzymatic transfer of the dideoxy sugar from a pyrimidine to a purine base. These compounds were inactive against HIV-1, HBV, HSV-1 and -2, VZV, and HCMV.     

Sense in antisense?

A correspondence between open reading frames in sense and antisense strands is expected from the hypothesis that the prototypic triplet code was of general form RNY, where R is a purine base, N is any base, and Y is a pyrimidine. A deficit of stop codons in the antisense strand (and thus long open reading frames) is predicted for organisms with high G + C percentages; however, two bacteria (Azotobacter vinelandii, Rhodobacter capsulatum) have larger average antisense strand open reading frames than predicted from (G + C)%. The similar Codon frequencies found in sense and antisense strands can be attributed to the wide distribution of inverted repeats (stem-loop potential) in natural DNA sequences.     

Some remarks on Rosen's quantum-mechanical approach to genetics

It has been suggested by Robert Rosen (Bull. Math. Biophysics,22, 227–255, 1960) that multiple alleles or pseudoalleles correspond to multiple cites of degenerate states of some quantum mechanical observable which acts as a source of primary genetic information. It is pointed out here that if the quantum mechanical states are determined by the different sequences of the purine and pyrimidine bases in the DNA molecule, the expected number of pseudoalleles would be much too large. The expected number is considerably reduced if we assume that a quantum mechanical state determines the coupling between a molecule of transfer RNA and the corresponding amino acid.