Recognition in nucleic acids and the anticodon families

Errors and alternatives in prebiotic replication and catalysis

The work on nonenzymatic nucleic acid replication performed by Leslie Orgel and co-workers over the last four decades, now extended by work on artificial selection of RNA aptamers and ribozymes, is generating some pessimism concerning the 'naked gene' theories of the origin of life. It is suggested here that the low probability of finding RNA aptamers and ribozymes within pools of random sequences is not as disquieting as the poor gain in efficiency obtained with increases in information content. As acknowledged by Orgel and many other authors, primitive RNA replication and catalysis must have occurred within already complex and dynamic environments. I, thus, propose to pay attention to a number of possibilities that bridge the gap between 'naked gene' theories, on one side, and metabolic theories in which complex systems self-propagate by growth and fragmentation, on the other side. For instance, one can de-emphasize nucleotide-by-nucleotide replication leading to long informational polymers, and view instead long random polymers as storage devices, from which shorter oligomers are excised. Catalytic tasks would be mainly performed by complexes associating two or more oligomers belonging to the same or to different chemical families. It is proposed that the problems of stability, binding affinity, reactivity, and specificity could be easier to handle by heterogeneous complexes of short oligomers than by long, single-stranded polymers. Finally, I point out that replication errors in a primitive replication context should include incorporations of alternative nucleotides with interesting, chemically reactive groups. In this way, an RNA sequence could be at the same time an inert sequence when copied without error, and a ribozyme, when a chemically reactive nucleotide is inadvertently introduced during replication.

Physico-chemical constraints connected with the coding properties of the genetic system

New insights on the origin of the genetic code, based on the analysis of the physico-chemical properties of its molecular constituents (RNA and amino acids), are reported in this paper. We point out a symmetry in the genetic code table and show that it can be explained by the nature of the anticodon-codon interaction. The importance of the strength of this interaction is examined and a correlation is found between the free-energy change (DeltaG(0)) of anticodon-codon association and the volume of the corresponding amino acids. This correlation is investigated in conjunction with the well-known one linking the hydrophobicity of the anticodons with that of the amino acids. We show that they can be considerated separately and that the energy vs. volume correlation may be explained by the process implicating the peptide bond formation between two successive amino acids during translation. This interpretation is supported by a statistical pattern of bases (purines or pyrimidines), observed in present coding genes, and by considerations involving the availability of the different kinds of amino acids. Finally, we try to explain the hydrophobicity correlation when reconstructing the events at the time of the so-called "RNA World". The whole of our investigation shows that the genetic code might be sufficiently robust to exist without the participation of pre-existing proteins, and that this robustness is a consequence of the physico-chemical properties of the four bases of the genetic system.

[Prediction of secondary structures of nucleic acids: algorithmic and physical aspects]

Prediction of secondary structures in nucleic acids requires both an adequate physical model and powerful calculation algorithms. In our approach, we cut the molecules in sections of which the contributions to the global energy are context-dependent but roughly additive. The structure of minimum energy is obtained by a tree search under constraints of binary incompatibilities. Our algorithm of the "incompatibility islets" is shown to be more powerful than the "bit parallel forward checking" algorithm, well known in Artificial Intelligence. Recurrent algorithms, proposed by other authors are even more rapid, but often miss the correct structures, for they demand a strict additivity of the energetic contributions, physically unjustified. New strategies, required to deal with molecules of more than 200 nucleotides are discussed. Our physical model has been improved by considering the special case of internal loops beginning with a G-A opposition. A bonus of 1.5 kcal. is attributed to such a feature, at each side of an internal loop. To illustrate our programs, we give the computed schemes for the 3' termini of the small subunit ribosomal RNA.

Computer building and folding of fictitious transfer-RNA sequences

In order to evaluate the common occurrence with which polynucleotides may adopt the cloverleaf configuration, 1150 random sequences were computer built and folded into their most stable secondary structure. Various constraints modulated the generation of the sequences: i) the base-pairing pattern, ii) the nucleotide composition, iii) the presence of assigned bases (modified or not) at certain sites, and iv) the chain length. In many cases, artificial tRNAs appear to require a more complex organization than a cloverleaf pairing scheme to achieve, as do natural molecules, the corresponding secondary structure. Moreover, the preferred foldings of sequences from 50 to 90 nucleotide long without an imposed pairing pattern usually contain two rather than three hairpin-loops. Implications concerning the emergence and the evolution of the protein-synthesis apparatus are discussed.

Prediction of pairing schemes in RNA molecules-loop contributions and energy of wobble and non-wobble pairs

Previously published models for predicting pairing schemes in RNA molecules, when applied to tRNA, give the clover leaf structure in only half the cases. We made a systematic investigation of the predictability of the clover leaf structure under various assumptions concerning the energetic contributions of single and double-stranded regions. We tested 21 different models and variants on a set of 100 tRNA sequences and many other variants on a smaller set of sequences. In our models we allowed not only G.C, A.U and G.U pairing, but also every other pair. Under conditions which are much less restrictive than those of previous attempts, we can nevertheless reach 90 per cent predictability for the clover leaf structure of tRNA. A most surprising and far-reaching result is that we can assign to C.G and C.C pairs binding energies quite close to the energies of G.U pairs, and still predict the clover leaf. The following ranking for non-complementary pairs was obtained : G.U, G.G and C.C, U.U, C.A, A.A and G.A, U.C. The main practical innovation which made possible the improvements in predictability are: i) not counting the stacking of base pairs separated by a bulge loop; ii) making the terminal C.C's in stems more stable than the terminal A.U's by merely -- 0.7 kcal; iii) replacing the distinction between G.C and A.U-closed loops by a distinction based on the presence of loop-favoring residues; iv) carefully adjusting the energetic balance between the various kinds of loops; v) narrowing the gap between the GC/GC and the GC/AU contributions; vi) using observations on nearest-neighbours in tRNA sequences to refine the contributions of G.U pairs.

Fluorescent tRNA derivatives and ribosome recognition

The use of fluorescent derivatives of tRNAPhe (yeast) in studies on tRNA conformation and on tRNA-ribosome recognition is described. Evidence is presented which indicates that under physiological conditions with respect to ionic strength and Mg2+ concentration, tRNAPhe exists in at least two conformations. The functional significance of this behavior is discussed on the basis of aminoacylation experiments. The investigation of the ribosome complexes of tRNAPhe labeled in the anticodon and D-loops has provided evidence suggesting that the presence of the codon, although not appreciably altering the apparent association constant, leads to qualitatively different complexes in which the tRNA appears to be rigidly bound to the codon even in the P-tRNA to the ribosome occurs in several steps, which take place only in the presence of the proper codon. One or more of these steps may represent codon-induced conformational changes of the tRNA molecule, which constitute the molecular basis of the highly specific binding of the tRNA to the ribosome.

Stacking of Crick Wobble pair and Watson-Crick pair: stability rules of G-U pairs at ends of helical stems in tRNAs and the relation to codon-anticodon Wobble interaction

The occurrence of the noncomplementary G-U base pair at the end of a helix is found to be governed by stacking interactions. As a rule, a G-U pair with G on the 5'-side of a Watson-Crick base pair exhibits strikingly greater stacking overlap with the Watson-Crick base pair than a G-U pair on the 3'-side of a Watson-Crick base pair. The former arrangement is expected to be more stable and indeed is observed 29 times out of 32 in the known transfer RNA molecules. In accordance with this rule, the major wobble base pairs G-U or I-U in codon-anticodon interactions have G or I on the 5'-side of the anticodon. Similarly, in initiator tRNAs, this rule is obeyed where now the G is the first letter of the codon (5'-side). In the situation where U is in the wobble position of the anticodon, it is usually substituted at C(5) andmay also have a 2-thio group and it can read one to four codons depending on its modifications. A G at the wobble position of the anticodon can recognize the two codons ending with U or C and modification of G (unless it is I) does not change its reading properties.

Allosteric mechanism for codon-dependent tRNA selection on ribosomes

We suggest that the interaction between a codon and its cognate tRNA induces conformational changes in the tRNA. We further suggest that sites on the ribosome preferentially bind tRNA in those conformations which require proper matching of codon and anticodon. According to this model, the codon functions as an allosteric effector which influences the conformation at various sites in the tRNA. This is made possible by the ribosome, which we suggest traps tRNA molecules in those conformation states that maximize the energy difference between cognate and noncognate codon-anticodon interactions. Studies of the interactions between tRNA molecules and their cognate codons in the absence of the ribosome have suggested that triplet-triplet interaction between codon and anticodon is far too weak to account for the specificity of the tRNA selection mechanism during protein synthesis. In contrast, we suggest that such affinity measurements do not adequately describe the interaction between a codon and its cognate tRNA. Thus, such experiments can not detect conformational changes in the tRNA, and, in particular, those stabilized by the ribosome.

Kinetic amplification of enzyme discrimination

The dependence of the accuracy of enzymatic systems on the mechanism of the catalyzed reaction is investigated, using a probabilistic approach. Certain mechanisms of reaction, involving a delay in one of the steps act as kinetic amplifiers of molecular discriminations. The relationship between our scheme for a delayed reaction [1] and Hopfield's scheme [i] is discussed.