Next Article in Journal
Structural, Vibrational and Electrochemical Analysis and Antibacterial Potential of Isomeric Chalcones Derived from Natural Acetophenone
Previous Article in Journal
Social Media Rumor Refuter Feature Analysis and Crowd Identification Based on XGBoost and NLP
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Transformation by Catalysis of Prebiotic Chemical Systems to Useful Biochemicals: A Perspective Based on IR Spectroscopy of the Primary Chemicals II. Catalysis and the Building of RNA

1
Department of Chemical Engineering, University of Lund, Box 124, SE 221 00 Lund, Sweden
2
Technologie DMI, 980 Rue Robert Brossard, Quebec, PQ 14X 1C9, Canada, [email protected]
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2020, 10(14), 4712; https://doi.org/10.3390/app10144712
Submission received: 24 February 2020 / Revised: 8 June 2020 / Accepted: 12 June 2020 / Published: 8 July 2020
(This article belongs to the Section Chemical and Molecular Sciences)

Abstract

:
Having found that carbonyl sulfide (COS), works well as a catalyst in the transformation of amino acids to polypeptides, we have now tested COS as a catalyst also for the formation of substances that might be thought of as partners in the building of RNA. The model used was selective energy transfer (SET). This model implies that a certain number of vibrational quanta are donated from the catalyst system and a corresponding number of quanta of the reactant accept the energy thus transferred. In this way, we found that carbonyl sulfide, COS, was a perfect catalyst for combining, first, five molecules of formaldehyde to form one molecule of ribose, and next, five molecules of hydrogen cyanide, HCN, to form one molecule of adenine, one of the nucleobases of RNA. However, beyond this, we found that COS was a perfect catalyst for precisely all reactions, needed to build the RNA, ribonucleic acid.

1. Introduction

In our previous paper [1] we have treated reactions involving amino acids and the formation of peptide bonds by catalysts such as carbonyl sulfide (COS) or others, like Fe(CN)64− and cyanamide (H2N-CN), finding that a resonance effect described by the selective energy transfer (SET) model works quite satisfactorily. In this paper, we will concentrate on the formation of RNA, ribonucleic acid, which is built from units such as PO43−, ribose and four kinds of nucleobases, i.e., adenine, cytosine, guanine, and uracil (Figure 1). The nucleotides, thus formed, are combined to a long chain with the phosphate ions and the sugar ribose, as coupling units. The reactions that result in these polymeric systems are greatly dependent on effective catalysts and we will also here apply the theory of “selective energy transfer” (SET), described in previous papers, e.g., [2,3,4]. The basic concept of SET is that the catalyst does have a vibration frequency that is equal to—or nearly equal to—that of the vibration in the reacting molecule that is supposed to lead to reaction. This means that a state of resonance will be established between the two systems: catalyst and reactant. When the systems are coupled in such a way, the ’energy transfer’ is carried out by transfer of one or several vibrational quanta, from catalyst to reactant, and the rate of the reaction is set by the rate of energy transfer [5].

The SET Model Exemplified

This model of catalysis was originally initiated from observations that the activation energies of catalytic reactions varied in a stepwise fashion [6]. We will in the following exemplify the function of SET by considering the formation of the primary constituents of RNA, i.e., the formation of ribose and the four nucleobases. By the action of the catalyst, the units that are building the various parts of RNA are accepting quanta from the COS molecule. In this way they reach – stepwise – their ability to react.
Furthermore, SET will be used to illustrate the binding of these nucleobases, as well as phosphate groups, towards the ribose structure In Figure 1, one can see how a chain of ribose molecules is formed by binding to two groups of phosphates using two different sites on each ribose molecule.
On the other side of this chain, four different nucleobases are bonded, one for each of the ribose molecules. By variation of the internal positions of the nucleobases a long chain is shaped, with a multitude of properties defined by this ordering and the physical parameters of the four nucleobases.
We will in the following exemplify the function of SET by considering the formation of the primary constituents of RNA, i.e., the formation of ribose and the four nucleobases. Furthermore, SET will be used to illustrate the binding of these nucleobases, as well as phosphate groups, towards the ribose structure.

2. The Formation of the Constituents of RNA

The existence of the relatively complicated molecules indicated in Figure 1 must not be taken for granted. In the conditions of an early atmosphere, one must ask from which simpler molecules they were formed. Further, is this formation possible to describe by the SET model?
In this context, we will introduce carbonyl sulfide, COS, as the catalyst. This gaseous substance is characterized by one very strong IR absorption band, at 2079 cm−1, [7], and consequently (e.g., p. 20, Formula I. 48 in [8]), also by a strong emission band of the same frequency.
Hence, in order for the SET model to work properly, the frequency of the accepting (reacting) vibration should be a factor x of 2079 cm−1. Here x could be 0.5, 1, 2, 3 etc., but also 1.5, 1.3333 or 1.25, etc.
As an example, if x = 2, or in other words, ν (reactant) = 2 × 2079 cm−1, it means that one quantum of the reacting molecule needs two quanta from the catalyst to achieve a reaction.
Similarly, if x = 1.3333, ν (reactant) = 4/3 × 2079 cm−1, or ν reactant/ω catalyst = 4:3. This means that 4 quanta of the catalyst match 3 quanta that will produce a reaction. To specify the relations between quanta of catalyst and reactant, one can form the expression 2073/ν (reactant) that we call the resonance condition (RC). The use of this quantity to describe which vibration of the reactant that is nearest a resonance situation for a given catalyst is exemplified, e.g., in Chapter 2.1 to follow.
SET has shown itself useful in many other cases. As an example, we can state that it is possible to calculate the so called “isokinetic temperature” [9].
In a series of hydrodechlorinations of chlorobenzene over Ni but with varying support, we observed that the Arrhenius lines were crossing each other at one temperature, the isokinetic temperature, 669 ± 2 K. The theory of SET made it possible to calculate the exact value to Tiso = 669.2 K in perfect agreement with the experimental value, mentioned above [10]. This result strongly points to the usefulness of SET.

A SET Approach to the Synthesis of Ribose

As shown in Figure 1, the basis of the RNA structure is the ribose link. We will now discuss the potential route of synthesis of this substance, starting from the simplest example of a compound containing carbon, oxygen, and hydrogen, namely, HCHO, formaldehyde. This simple molecule must expose a similarly simple IR spectrum. The Table 1 gives the data from the investigations of Herzberg, (table 76 in [10]), characterizing the 4 × 3 − 6 = 6 normal vibrations as shown in Figure 2. Furthermore, it is in (table 42 of [10]), that one finds the exact value, 2079 cm−1, of the ν3 vibration of the catalyst COS.
The first point of our discussion is to find, if there exists, among the available data, a pronounced resonance condition, RC, i.e., 2079/ν cm−1. With ‘pronounced’ we mean a ν−value that gives a small value of Δ (the difference between 2079 and ν cm−1), such as is explained as follows: full resonance occurs when the frequencies of catalyst and reactant are exactly equal. The relative difference between the ratio ω/ν and the nearest ratio of small integers, we designate as Δ%. Here ω is the frequency of the catalyst and ν is the frequency of the reactant, e.g., in the row for ν 2 of Table 1, the closest ratio of RC is 1.200 = 6:5.
Of the three items in the Table 1 that show Δ−values less than 1%, ν1 is a fully symmetric vibration which does not relates to a strong distortion of the molecule, and the existence of the one named 2 ν 6 has been strongly doubted (foot note 96 in [10]). However, the ν 2 vibration has a very strong IR absorption/emission and, furthermore, it is mainly a C=O vibration where 5 quanta of the 2079 cm−1 vibration contribute energy to 6 quanta of the 1743.6 cm−1 vibration. This means, that at the extreme points of the vibration, i.e., when the adsorption of quanta from the catalyst is at its maximum, the 2p–2p interaction is considerably weakened and the bond between the carbon and the oxygen is close to a single bond.
We designate this state as C*--O*, in Scheme 1, where the asterisk indicates (almost) half-filled 2p orbitals. In Scheme 1 we now illustrate the reaction steps that will follow from the excitation of formaldehyde molecules via the catalyst COS.
(1) indicates, in relation to what is said above, that energy is transferred from the catalyst COS to the formaldehyde molecule. At the point of reaction, when enough vibrational quanta have been accepted by the C=O bond of the H2CO reactant, the double bond between C and O is weakened and self-supported, half-filled 2p orbitals (designated as *) are formed.
(2) illustrates the meeting of two equally excited molecules. The first happening is that one neutral hydrogen atom of the right-hand species is transferred to the (almost) half-filled 2p orbital of the oxygen of the left-hand species. This results, pro primo, in the formation of an OH group (to the left) and, pro secundo, in the creation of a free valency at the carbon from which the H-atom has been snatched away. Thus, we get one C-atom with a half-filled 2p orbital and another C-atom with an (almost) half-filled 2p orbital. If the molecules concerned are positioned one above the other, this results in a 2p, C-C bond, binding the two molecules together as indicated also in the left-hand formula of (c).
(3) HO―CH(H)―C*(H)―O* + (H)H>C*―O* => HO―CH(H)―C(OH)H―C*(H)―O* indicates that the excitation of the right-hand part of the product from (2) has not lost its excitation. Thus, when meeting anew with a COS-excited H2CO molecule the same procedure, as described above, can occur. This produces a molecule containing three carbon atoms as shown in (4).
In (4) and (5), two other excited formaldehydes, one in each row, are added in a similar way. This results in a molecule with five carbon atoms in which the π–orbital of the last added CO group is still split into C*-O*. This excitation, however, cannot last forever, and the five-carbon molecule relaxes to an aldehyde as shown in (6).
The linear molecule of Row 1f enters an equilibrium with a five-membered ring (Figure 3) and this process puts an end to further growth via carbon additions.

3. The SET Synthesis of Adenine from Hydrogen Cyanide (HCN)

Although adenine can be synthesized under laboratory conditions from HCN, e.g., [11,12,13], it might be interesting to the possibilities for a SET synthesis under prebiotic conditions when HCN and COS came streaming out together from an undersea vent. One notes then that the molecular vibrations of HCN as given by Herzberg (table 42 and figure 61 in reference [10]) are: ν1 = 2089 cm−1, ν2 = 712 cm−1 and ν3 = 3312 cm−1. We reproduce here in Figure 4 (the data from figure 61 of Reference [10]).
As said above, the ν3 vibration of COS is 2079 cm−1 and that will give a good possibility of resonance between the two vibrations now described (ν3 of COS and ν1 of HCN).
Therefore, the more energy that is transferred from COS to HCN, the larger will be the C-N distance in HCN. Because of this, the weaker will be the 2p–2p interaction in the HCN molecule; Scheme 2.
In the Scheme 2. is shown excitation of the HCN molecule: Δ E is the energy of one vibration quantum, given from the catalyst COS; n is the number of such quanta to reach reaction. Then the triple bond between C and N is destroyed and the resulting “free” p orbitals are designated by a point.
On the other hand, there will appear “freer” 2p orbitals ready to join with a 2p orbital of a neighboring molecule, whereby a single bond will be created between the two molecules. In Scheme 3 we illustrate this event with four of the five HCN molecules that are supposed to build the adenine molecule. The bond between N and C in the HCN molecule, activated by resonance with COS, is depicted as a double bond, what might possibly be an overstatement. Anyway, the appearance of a C-N single bond between two molecules makes it imperative that the two linear molecules attach in an angle approaching 120 degrees (weakly indicated in Scheme 3).
In Scheme 4 the gap that results between the first and the fourth HCN molecule is illustrated. This gap will be filled by the carbon atom from the fifth HCN molecule. It is difficult to tell in what form the HCN molecule molecule appears, but most probably it is the isomeric form HNC, with two half-used p-orbitals available on the carbon atom (Scheme 4). The important point is that the carbon atom will find itself surrounded by three neighbouring atoms, more than willing to form new bonds with the intruder. This results, for one sake, in a weakening of the original CN triple bond so that the nitrogen atom is expelled and can bind to the carbon atom of the starting HCN unit. The two H atoms, that are set free in this maneuver, add to the N atom forming the top NH2 group.
Secondly, the result of the penetration of the last HCN molecule (or possibly CN ion or CNH) is, of course, that three new bonds are established what gives a closed structure; (Scheme 4 and Figure 5).
It is interesting to note that Jung and Choe [13] find, by CBS-QB3 calculations, that the formation in gas phase of oligomers (HCN)n where n = 2 up to n = 5 (adenine) requires very high activation energies. For each step of that reaction “a few hundred kJ“ are used [13] and consequently a catalyst is strongly needed. These authors used a proton as a catalyst, that strongly decreased the activation energy. Nevertheless only “some oligomerization stages” were predicted [13] to occur in an atmosphere like that around the moon Titan (of the planet Saturn). No formation of adenine could, however, be foreseen. Near the surface of the planet Earth things might have been better. From hot volcano eruptions lots of hydrogen cyanide streamed out, but also lots of warm, already vibrationally excited carbonyl sulfide (COS). As this substance has shown good capacity as a catalyst [1], we allow ourselves to consider it as such also in this context.
If it has entered this structure like a cork to a bottle or it is squeezed between the outer parts of that structure like a nail in the gap of a pair of tongs is difficult to tell. The creation of a C-C connection establishes the two ring structures coupled together as is also shown in Figure 5.

3.1. Splitting One Oxygen Atom From CO2

RNA is constituted from four different nitrogen bases, Figure 1. Besides Adenine, which is treated above, they are uracil (U), guanine (G), and cytosine (C). All of the latter three molecules have one common trait that is lacking for adenine: they all contain oxygen. As the reactions discussed here were taking place under non-oxidative conditions, one might ask about the origin of that oxygen. Two possibilities can be considered: CO and CO2, both of volcanic origin. In this paper we concentrate on CO2.
Thus, given the existence of CO2 in the prebiotic atmosphere, how can it split off one separate oxygen atom, to take part in the construction of the nitrogen bases?
Nitrogen atoms appear, of course, also in the substances under study. Thus, one can assume that the atmosphere contained at least CO, CO2 and ammonia. There were also other simple substances, e.g., HCN that could have acted to build the nitrogen bases.

3.1.1. Carbon Dioxide Degradation

Now, one can observe that there appears in the CO2 infrared spectrum a non-fundamental vibration composed by one ν 6 (stretching) and one ν 10 (bending of the three-atom line) [10]. This composite vibration appears at 2076.5 cm−1 and is designated as a medium strong IR absorber, (Table 56 in Reference [10]).
The frequency of this composite vibration is very close to that of the gaseous catalyst COS (2079 cm−1), with a very strong emission intensity. This equality of frequencies implies that a pronounced state of resonance could appear, resulting in an activation of the oxygen atoms. This atom can probably combine with the two hydrogens equally activated vibrating with e-symmetry in an ammonia molecule, resulting in the formation of water and a peptide group, discussion below.

3.1.2. Ammonia Vibration Modes

The catalyst for the reaction between NH3 and CO2 has been suggested to be gaseous COS as the critical vibrations are agreeing with the concept of SET. Among the vibrations of N-H stretch reported by Herzberg [10] we find one species at 4176 cm−1, meaning that it can obtain reaction energy from two quanta of COS (2 × 2079 = 4158 cm−1). This composite vibration is assigned as ν2 (a) + 2 ν4 (e), [10], i.e., the degenerate vibration appears twice that one with singular character. It seems to us that this composite vibration causes the expulsion of two hydrogen atoms, Figure 6, that reacts with the oxygen atom, activated as indicated above. One might say that both the H as well as the O atoms are in a statu nascendi, and hence have the capacity of reacting with each other.

3.2. The Formation of the Oxygen-Containing Nucleobases Uracil, Guanine, and Cytosine from a SET Point of View

Starting from the peptide group formed as described in Figure 6, (or, in one case, the two peptide groups) and the empirical formula for the molecule in question we try to understand the building of the nitrogen bases.
One notes that in all cases, Scheme 3 and Scheme 4 and Table 2, Table 3 and Table 4, catalysis is a strongly working factor. In adenine, as described above, the entire structure is built up from five HCN molecules, where the vibration of the carbon-nitrogen bond is activated by the gaseous catalyst COS.
For the three remaining nitrogen bases, the structure is started by the formation of a peptide group created by the action of the same catalyst, COS, on the splitting of a CO2 molecule mixed with ammonia as described above, Figure 6. The results are collected in Table 2, Table 3 and Table 4.
As the peptide unit covers two positions in the ring system of guanine (Table 3), only four HCN molecules are needed to complete the structure. As in the case of adenine, the last HCN molecule (or rather the isomer form CNH) is split into a carbon, entering a double bond, and the remaining N-H unit that takes part in the formation of an amino group. In both cases (adenine and guanine) the amino group is positioned on the larger ring of the molecule. In addition, in both cases, the carbon atom is squeezed by the outer part of the structure.
The remaining two nitrogen bases form only one ring, but differs between themselves in that uracil contains two peptide groups. From the empirical formula of uracil it follows that an additional C2H2 group is needed to complete the structure (Table 3). We suggest that this group originates from an ethene molecule. In addition, in this case catalysis is of decisive importance. One finds that (in Herzberg’s Table 92; Reference [10]) there is tabulated a strong fundamental vibration, ν9CH(b2u), at 3105.5 cm−1. The (inverted) resonance condition is therefore 3105.5/2079 = 1.4937; this means that the difference from 1.5000 is 0.0063 or 4.2 permillage. This is a very reasonable value for resonance.
This means that two quanta of ν9 accepts energy from three quanta of the catalyst. Now it turns out [10], that ν9 has a form as shown in Figure 7, that implies that at high excitation of their vibration, the two "upper" H-atoms are distancing themselves (and together) from the compact rest of the molecule, i.e., C2H2. This results, after bond breaking at a suitable excitation level, in a H2 molecule and a remaining rest of H-C*=C*-H, with two free valencies, that can form bonds with other -excited- atoms (Table 2 and Table 4).
Cytosine has the most complex structure of the four substances discussed here. We suggest (Table 4), that HCN, NH3 and C2H4 are involved, in a 1:1:1 ratio and that all three substances are activated by the catalyst COS as described above.
Most reasonably ammonia, NH3, will decompose into =NH and H2 (compare Figure 6) and C2H4 will decompose into C2H2 and H2 (compare Table 2 and Table 4).
Figure 8 shows the structure of uracil, and in Table 2 is shown its gradual growth.
Figure 9 shows the structure of guanine, in Table 3 is shown its gradual growth.
Figure 10 shows the structure of cytosine, and in Table 4 is shown its gradual growth.

4. Additions to the Ribose Structure

As seen from Figure 1, ribose is the carrying tool for different molecules, nucleobases, where the combining part of the ribose molecules are phosphate groups.
These chains, presenting a collection of biophysical properties are, as is well known, of great biological importance. Therefore, it is of essential importance to know how the mechanics of keeping such a chain together works.
In addition, here, we find that carbonyl sulfide catalysis works to establish good bonding conditions, whether it is the ribose-phosphate bonding, or it is the ribose nucleobase combination that is concerned.
In the following, we will discuss testing methods that might indicate how this chain formation takes place.

4.1. Phosphorylation of Carbohydrates, e.g., Glucose and Ribose

To evaluate the phosphorylation of carbohydrates, e.g., glucose or ribose, in the presence of COS, we have collected some relevant IR-data. We have concentrated here on the ν 3 vibration of COS, testing if ring vibrations of ribose are as stimulated by the catalyst as the vibrations (e.g., ρw (NH2)) described in the previous investigation [1].
This activity might to a part depend on the strong IR absorption, and consequently very strong IR emission band of COS at 2079 cm−1 [9] and on the relatively intense P-O stretching vibrations of the phosphate groups at about 1045 cm−1. As cited in Reference [4], these data will give a 2:1 ratio for the COS/PO4 vibrational quanta, 2079/1045 = 1.9895 which is as close to resonance as one can demand. Here 1045 cm−1 is the mean of 1044 cm−1 [14] and 1046 cm−1 [15], compare Table 5.

4.1.1. Glucose

The above-mentioned observation is of great interest as the vibrationally excited PO43− is a key substance for the phosphorylation of, e.g., glucose [4]. One notes that the reaction, just mentioned, also involves an activation—via the COS vibration—of the C(6)-OH vibration of glucose at about 1048 cm−1, [4], with an equally good resonance state as obtained for the above-mentioned activation of the PO4 group. The relevant vibrational data [4] are presented once more here in Table 5, where [16,17,18] refer to data for glucose and [14,15] refer to data for phosphate.
Thus, e.g., an activated PO4H2− might approach an equally activated C(6)-OH group, resulting in the formation of a H2O molecule and the binding of a phosphate group to the sugar molecule.

4.1.2. Ribose

If the above reasonings hold also for a pentose substance they will indicate the starting point for the formation of the ATP molecule (adenosine triphosphate). Actually, one finds (Table 6, References [19,20,21]) a corresponding C(5)-OH vibration of ribose, with a mean value, ca 1040 cm−1 that fits rather well for an interaction with the phosphate system (1045 cm−1) shown in Table 6.

4.2. The Attachment of Adenine to Ribose

Adenine (Figure 5) is one of the important building blocks of ATP or adenosine triphosphate. We use here ATP as a small-scale example of adenine bonding. If we accept that the C(5) position of the ribose is the one first occupied by the PO4 group it is tempting to assign the opposite position of the ribose molecule, C(1), for anchoring the adenine molecule. This view is supported by well known structure data, Figure 11. We will use ATP as a model example for describing the bonding of nucleobases, all with an N-H group, towards the C(1) position of ribose.
As the picture of ATP shows (Figure 11), the connection between the adenine N(9)-H and the ribose C(1)-OH keeps the different groups together. How does this extremely important bond come about? It is known that the C-O + C-C vibration is found at about 1040 cm−1 (table 6, in [19,20,21]).
But what is known about vibrations of the N(9)-H group that might be in close resonance with the vibrations of the C(1)-OH moiety? To follow the treatment discussed above we plot a series of FTIR data [22] as shown in Table 7.
Table 7 shows a series of FTIR data on anhydrous adenine.
From the data given in Table 7, we can safely assume that Row 5 represents a condition where the ‘extra terms’ give the least disturbance when calculating the vibration frequency of δ(N9-H). We, therefore, set forth to find a correction that will give us a ‘true’ value of this parameter, and thereby understand the transfer of energy from catalyst to the intended reactant.
As above (Table 7; Row 5) we can state that (1.521–1.500)/1.5 = 14 o/oo.
Hence, to reach an expected, somewhat smaller, value of the resonance condition, one has to increase the denominator with a term that we call Δ”.
Δ” = 0.014 × 1367 cm−1 = 19 cm−1.
Thus, the true value of δ (N9-H) = 1367 + 19 = 1386 cm−1.
We can now check the calculations by forming the expression for resonance condition, R.C. = 2079/1386 = 1.500 = 3:2. This gives full support to our treatment.
If we express the value of a vibrational quantum in cm−1 units, this means that two catalyst quanta = 2 × 2079 cm−1 = 4158 cm−1, and three reactant quanta = 3 × 1386 cm−1 = 4158 cm−1.
Thus, two vibrational quanta from the catalyst are donated and three quanta are accepted by the reactant for the reaction to go.
Furthermore, one can now check if the conditions for resonance-and consequently energy exchange-between the adenine molecule and ribose is at hand. Using the mean value of the C-OH vibration of ribose as 1040 cm−1 (Table 6), we can detect that the ratio between adenine and ribose vibrational quanta is 1386 cm−1/1040 cm−1 = 1.3327, i.e., 4:3.
This means, that for the reaction of binding adenine to ribose to occur, the energy of excitation of the reaction partners should be equal: 3 × 1386 cm−1 = 4158 cm−1, and 4 × 1040 cm−1 = 4160 cm−1.
This result suggests that this type of bonding should hold also for the other three nucleobases. Thus, the deductions made above strengthen our reasoning in the introduction. It is also of interest to note that δ(N9-H) = 1386 cm−1 is in good correspondence to the same data (1382 cm−1), measured by the Inelastic Neutron Scattering (INS) method [22].

4.3. Total Synthesis

When all four nucleobases are present, they can form RNA by combining with ribose molecules which, in their turn, are kept together by PO43− units as described in detail for adenine in Section 3.1 and Section 3.2. A fraction of the result is shown in Figure 1. It must be emphasized that the coupling between the remaining three nucleobases and ribose is like that described for adenine and ribose (Section 4.2):
Uracil, with two N-H groups (Figure 8) uses that group that gives the largest steric availability. Figure 1 shows that the choice of bonding site fulfills the requirements.
For guanine, with a structure like adenine, it is the N9-H group (compare Figure 9 and Figure 1) that reacts with the C1-OH group.
Cytosine with only one N-H group reacts with that one, cf. Figure 1 and Figure 10.
One may note that in many, not to say most, of the reactions described here, there occurs the formation and removal of one molecule of H2O, i.e., a “condensation”. The fact that COS is the catalyst in all these reactions is remarkable (see below).

5. Conclusions

We have found that it seems possible to build (or suggest a scheme for building) all the nucleobases that make up the RNA chain from a few simple molecules, such as HCN, CO2, NH3, C2H4, providing that a certain order between the reactants is kept up, and that one uses the concept of SET (selective energy transfer) for catalysis of the reaction steps involved.
We have also found that for all reactions considered there is a ratio in small numbers; 1:2, 1:1, 3:2, 4:3, etc. between the vibration quanta of catalyst and reactant, well in accordance with the rules of the SET model.
What is extraordinary and definitely interesting from the view-point of “the origin of life”, is that in all cases the catalyst is one and the same substance, the carbonyl sulfide, COS; or rather SCO to emphasize its relation to carbon dioxide (Table 8). One finds reason to use COS as a catalyst not only for the synthesis of the ribose and the four nucleobases, not only for the splitting of carbon dioxide and the splitting of ammonia and ethene, but also for the conversion of amino acids to chains of peptide groups as described in Reference [1]. Further, this gas, COS, came from the—at the time—not so uncommon volcanic eruptions! We are tempted to cite the expression [23], coined by the Nobel Laureate Christian de Duve: COS must be “the dust of life”. (This, of course, is our later-days interpretation).
But that phrase, “the dust of life”, also stresses the weak point of our results; a gaseous substance would have difficulties to stay inside the walls of a living cell. It is now well known, however, that small, solid particles of montmorillonite (a type of petrified clay) might be trapped in such protocells in which RNA might operate [24,25]. Especially interesting for us is the concept “everything first”, coined in this context by Sutherland [26,27], indicating that if all the reactions, needed to build the cell and its contents, were served by one and the same catalyst, the various components must be formed simultaneously, side by side irrespective of the catalyst being a gas or a petrified clay. Sutherland’s concept also implied that no special part of the total cell was assigned to lead the run, having to see the others slowly following up.
Therefore, hopefully, we will be able to discuss similarities and differences between gas systems and solid-state catalysts in the next part of this breath-taking tale on the power of SET.

Author Contributions

The activities for producing the manuscript relate as follows: conceptualization, R.L. and A.M.; literature search and evaluation, A.M.; calculations, R.L.; methodology, R.L. and I.O.; writing—original draft, A.M. and R.L.; data collection and evaluation, A.M. and R.L.; drawings and manuscript corrections, I.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful for the support for this work provided by the University of Lund and the Department of Chemical Engineering of LU. This support was aimed at covering the costs of publishing the results of the study in full Open Access style. One of us (R. L.) will strongly thank his daughter, Karin Lilja, B.Sc., for good guidance in computer use. Likewise, we, the authors, want to express our gratitude to Mr. Farhan Zaman of Ottawa, Canada, for his computer related help. We are also indebted to Mr. Lars Fembro for his help with the figures of vibration patterns.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Larsson, R.; Malek, A. The transformation by catalysis of prebiotic systems to useful biochemicals: A perspective based on IR spectroscopy of the primary chemicals. I. The Synthesis of Peptides by the Condensation of Amino Acids. Appl. Sci. 2019, 10, 928. [Google Scholar] [CrossRef] [Green Version]
  2. Larsson, R. A Model of Selective Energy Transfer at the Active Site of the Catalyst. J. Mol. Catal. 1989, 55, 70–83. [Google Scholar] [CrossRef]
  3. Larsson, R. Concluding remarks on the theory of selective energy transfer and exemplification on a zeolite kinetics study. Monatsh. Chem. 2013, 144, 21–28. [Google Scholar] [CrossRef]
  4. Larsson, R. A SET Approach to the Interplay of Catalysts and Reactants. Catalysts 2018, 8, 97. [Google Scholar] [CrossRef] [Green Version]
  5. Kittel, C.; Knight, W.D.; Ruderman, M.A. Mechanics; Berkeley Physics Course; McGraw-Hill: New York, NY, USA, 1966; Volume 1, Chapter 7. [Google Scholar]
  6. Larsson, R. On the Stepwise Change of the Activation Energy of Catalytic Reactions. Z. Phys. Chemie, Leipzig 1987, 268, 721–732. [Google Scholar] [CrossRef]
  7. NIST Chemistry Webbook. Available online: http://webbook.nist.gov/chemistry (accessed on 14 November 2014).
  8. Herzberg, G. Spectra of Diatomic Molecules, 1st ed.; Van Nostrand Co., Inc.: New York, NY, USA, 1950. [Google Scholar]
  9. Keane, M.A.; Larsson, R. Isokinetic behavior in gas phase catalytic hydrodechlorination of chlorobenzene over supported nickel. J. Mol. Catal. 2007, 268, 87–94. [Google Scholar] [CrossRef]
  10. Herzberg, G. Molecular Spectra and Molecular Structure, II. Infrared and Raman Spectra of Polyatomic Molecules; D van Nostrand Company, Inc.: Princeton, NJ, USA, 1962. [Google Scholar]
  11. Oró, J. Mechanism of Synthesis of Adenine from Hydrogen Cyanide under Possible Primitive Earth Conditions. Nature 1961, 191, 1193–1194. [Google Scholar] [CrossRef]
  12. Voet, A.B.; Schwartz, A.W. Prebiotic adenine synthesis from HCN. Evidence for a newly discovered major pathway. Bioorg. Chem. 1983, 12, 8–17. [Google Scholar] [CrossRef]
  13. Jung, S.H.; Choe, J.C. Mechanisms of prebiotic adenine synthesis from HCN by oligomerization in the gas phase. Astrobiology 2013, 13, 465–475. [Google Scholar] [CrossRef]
  14. Frost, R.L.; Xi, Y.; Scholtz, R.; Lopez, A.; Belatti, F.M. Vibrational spectroscopic characterization of the phosphate mineral hureaulite–(Mn,Fe)5(PO4)2(HPO4)2*4 (H2O). Vib. Spectrosc. 2013, 66, 69–75. [Google Scholar] [CrossRef] [Green Version]
  15. Mejoo, S.; Maneeprakorn, W.; Winotai, P. Phase and thermal stability of nanocrystalline hydroxyapatite prepared via microwave heating. Thermochim. Acta 2006, 447, 115–120. [Google Scholar] [CrossRef]
  16. Kovacs, A.; Ivanov, A.Y. Vibrational Analysis of D-Glucose Trapped in Ar-Matrix. J. Phys. Chem. B 2009, 113, 2151–2159. [Google Scholar] [CrossRef] [PubMed]
  17. Ibrahim, M.; Alaam, M.; El-Haes, H.; Jalbout, A.F.; de Leon, A. Analysis of the structure and vibrational spectra of glucose and fructose. Ecletica Quimica 2006, 31, 15–21. [Google Scholar] [CrossRef]
  18. Vasko, P.D.; Blackwell, J.; Koenig, J.L. Infrared and Raman Spectroscopy of Carbohydrates. Part II: Normal Coordinate Analysis of D-Glucose. Carbohydr. Res. 1972, 23, 407–416. [Google Scholar] [CrossRef]
  19. Carmona, P.; Molina, M. Raman and infrared spectra of D-ribose and D-ribose 5- phosphate. J. Raman Spectrosc. 1990, 21, 395–400. [Google Scholar] [CrossRef]
  20. Mathlouthi, M.; Seuvre, A.M.; Koenig, J.L.F.T.-I.R. Laser-Raman spectra of D-Ribose and 2-Deoxy-D-erythro-Pentose (2-Deoxy-D-Ribose). Carbohydr. Res. 1983, 122, 31–47. [Google Scholar] [CrossRef]
  21. Tipson, R.S.; Isbell, H.S. Infrared Absorption Spectra in the Study of Mutarotational Equilibria of Monosacharides. J. Res. NBS—A Phys. Chem. 1962, 66A, 31–58. [Google Scholar] [CrossRef]
  22. Lopes, R.P.; Valero, R.; Tomkinson, J.; Marques, M.P.M.; de Carvalino, L.A.E.B. Applying vibrational spectroscopy to the study of nucleobases—Adenine as a case-study. New J. Chem. 2013, 37, 2691–2699. [Google Scholar] [CrossRef] [Green Version]
  23. De Duve, C. Vital Dust: Life as a Cosmic Imperative; Basic Books: New York, NY, USA, 1995. [Google Scholar]
  24. Luisi, P.L.; Walde, P.; Oberholzer, T. Enzymatic RNA synthesis in self-reproducing vesicles: An approach to the construction of a minimal synthetic cell. Ber. Bunsenges. Phys. Chem. 1994, 98, 1160–1165. [Google Scholar] [CrossRef]
  25. Adamala, K.; Szostak, J.W. Nonenzymatic Template-Directed RNA Synthesis Inside Model Protocells. Science 2013, 242, 1098–1100. [Google Scholar] [CrossRef] [Green Version]
  26. Powner, M.W.; Gerland, B.; Sutherland, J.D. Synthesis of activated pyrimidine ribonucleotides in prebiotically plausible conditions. Nature 2009, 459, 239–242. [Google Scholar] [CrossRef] [PubMed]
  27. Patel, B.H.; Percivalle, C.; Ritson, D.J.; Duffy, C.D.; Sutherland, J.D. Common origins of RNA, protein, and lipid precursors in a cyanosulfidic protometabolism. Nat. Chem. 2015, 7, 301–307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. A fraction of the long chain of RNA, ribonucleic acid, showing the ‘back-bone’ of phosphate groups and ribose. The ribose binds to the phosphate in two positions (C (3) and C (5)) and to one out of four nucleobases in the C (1) position.
Figure 1. A fraction of the long chain of RNA, ribonucleic acid, showing the ‘back-bone’ of phosphate groups and ribose. The ribose binds to the phosphate in two positions (C (3) and C (5)) and to one out of four nucleobases in the C (1) position.
Applsci 10 04712 g001
Figure 2. Molecular vibrations of, e.g., H2CO.
Figure 2. Molecular vibrations of, e.g., H2CO.
Applsci 10 04712 g002
Scheme 1. Catalytic synthesis of ribose reaction steps 1–6.
Scheme 1. Catalytic synthesis of ribose reaction steps 1–6.
Applsci 10 04712 sch001
Figure 3. Ribose (β-D-ribofuranose). Numbering code: first carbon on the left of the oxygen atom is C1. Last position in the five-ring is C4. The carbon in the side branch is C5.
Figure 3. Ribose (β-D-ribofuranose). Numbering code: first carbon on the left of the oxygen atom is C1. Last position in the five-ring is C4. The carbon in the side branch is C5.
Applsci 10 04712 g003
Figure 4. Vibrational pattern of hydrogen cyanide (HCN). Note that ν1 implies an increase of the C-N distance.
Figure 4. Vibrational pattern of hydrogen cyanide (HCN). Note that ν1 implies an increase of the C-N distance.
Applsci 10 04712 g004
Scheme 2. Excitation of the HCN molecule.
Scheme 2. Excitation of the HCN molecule.
Applsci 10 04712 sch002
Scheme 3. An excited molecule is attacking another exited HC-N or an already formed polymer unit, HnCnNn.
Scheme 3. An excited molecule is attacking another exited HC-N or an already formed polymer unit, HnCnNn.
Applsci 10 04712 sch003
Scheme 4. A carbon atom from the fifth hydrogen cyanide (HCN) unit in the grip of the already formed H4C4N4 structure.
Scheme 4. A carbon atom from the fifth hydrogen cyanide (HCN) unit in the grip of the already formed H4C4N4 structure.
Applsci 10 04712 sch004
Figure 5. Structure of adenine. The numbering code is as follows: the N atom to the left of the NH2 carrying C-atom is no 1, and following, anti clockwise around the six-membered ring so that the above-mentioned C-atom takes position no 6. Thereafter, from positions 5 and 4, around the five-membered ring one gets positions 7, 8, and 9. This means that the N-H group is number 9.
Figure 5. Structure of adenine. The numbering code is as follows: the N atom to the left of the NH2 carrying C-atom is no 1, and following, anti clockwise around the six-membered ring so that the above-mentioned C-atom takes position no 6. Thereafter, from positions 5 and 4, around the five-membered ring one gets positions 7, 8, and 9. This means that the N-H group is number 9.
Applsci 10 04712 g005
Figure 6. The catalyzed reaction between ammonia and carbon dioxide.
Figure 6. The catalyzed reaction between ammonia and carbon dioxide.
Applsci 10 04712 g006
Figure 7. The ν9 vibration pattern of an X2Y4 planar molecule—let it be C2H4 with the hydrogen atoms in periphery —from part of a scheme by Herzberg (figure 57 in [10] here ν9 is one out of twelwe symmetry coordinates.).
Figure 7. The ν9 vibration pattern of an X2Y4 planar molecule—let it be C2H4 with the hydrogen atoms in periphery —from part of a scheme by Herzberg (figure 57 in [10] here ν9 is one out of twelwe symmetry coordinates.).
Applsci 10 04712 g007
Figure 8. The structure of uracil.
Figure 8. The structure of uracil.
Applsci 10 04712 g008
Figure 9. Structure of Guanine.
Figure 9. Structure of Guanine.
Applsci 10 04712 g009
Figure 10. Structure of cytosine.
Figure 10. Structure of cytosine.
Applsci 10 04712 g010
Figure 11. Adenosine triphosphate (ATP) is a combination of one adenine, one ribose and three phosphate groups.
Figure 11. Adenosine triphosphate (ATP) is a combination of one adenine, one ribose and three phosphate groups.
Applsci 10 04712 g011
Table 1. Treatment of observed vibrations of H2CO.
Table 1. Treatment of observed vibrations of H2CO.
Assignmentν cm−1R.C. = 2079/νΔ%ω/ν
ν 611671.7813.17:4
ν 512801.6244.35:3
ν 315031.3835.04:3
ν 21743.61.1920.676:5
2 ν 62081-0.11:1
ν 127800.7480.273: 4
ν 428740.7237.23:4
2ν 329730.6996.92:3
Table 2. The gradual growth of uracil.
Table 2. The gradual growth of uracil.
UracilEmpirical FormulaC4H4N2O2
Substance addedResulting StructureResulting empirical formulaRemaining formula
Two peptide groups Applsci 10 04712 i001C2H2N2O2C2H2
Applsci 10 04712 i002 Applsci 10 04712 i003C4H4N2O2H2
Table 3. The gradual growth of guanine.
Table 3. The gradual growth of guanine.
GuanineEmpirical FormulaC5H5N5O
Added substanceResulting StructureResulting empirical formulaRemaining formula
One peptide group Applsci 10 04712 i004CHNOC4H4N4
3 HCN Applsci 10 04712 i005C4H4N4OHCN
HNC Applsci 10 04712 i006--
- Applsci 10 04712 i007C5H5N5O-
Nothing added or remaining.
Table 4. The gradual growth of cytosine.
Table 4. The gradual growth of cytosine.
CytosineEmpirical FormulaC4H5N3O
Added substanceResulting StructureResulting empirical formulaRemaining formula
One peptide group Applsci 10 04712 i008CHNOC3H4N2
HCN Applsci 10 04712 i009C2H2N2OC2H3N
NH3 =NH + H2 Applsci 10 04712 i010C2H3N3OC2H2 + H2
Applsci 10 04712 i011 Applsci 10 04712 i012C4H5N3O2 H2
Table 5. IR data (cm−1) for glucose C(6)O-H and solid-state phosphates [4]. The quoted assignments are from respective authors.
Table 5. IR data (cm−1) for glucose C(6)O-H and solid-state phosphates [4]. The quoted assignments are from respective authors.
AssignmentGlucosePhosphateReferences
PO43−/HPO42− asymmetric stretching-1044[14]
PO43−asymmetric stretching.-1046[15]
C-O-H 1047-[16]
ν CO + ν CC1050-[17]
CO (79%) + CC (19%) PED1047-[18]
Mean value10481045-
PED = Percentage Electron Density.
Table 6. Vibrational data for phosphates and ribose.
Table 6. Vibrational data for phosphates and ribose.
RowAssignmentRibosePhosphateR.C.Δ ‰Reference
1PO4 3−/HPO4 2− asymmetric stretching-10441.9914.5[14]
2PO4 3− asymmetric stretching-10461.9886.0[15]
3νC-O + ν(C–C),etc.1044 *-1.9829.0[19]
4δ(O-C-H)1035-2.0094.5[20]
5Not assigned1041-1.9971.5[21]
Mean value10401045---
All data in cm−1. All assignments relate to the respective authors. R.C. = “resonance condition” means the ratio between the COS catalyst frequency (2079 cm 1) and the tabulated frequency. Δ = difference (in per-millage) to nearest even digits. Mean value of Δ ‰ = 5 per-millage. * Mean value of ν CO + ν CC (1078 cm−1) and νC(1)O + νC(1)C (1010 cm−1) = 1044 cm−1. [Reference 20; Table 2 of that reference]. Because of structural information, we have omitted data related to C(2).
Table 7. Series of FTIR data on anhydrous adenine, numbers 1–5, from an investigation with detailed analysis of the spectral assignments. The “resonance condition” is calculated as was done in Table 1.
Table 7. Series of FTIR data on anhydrous adenine, numbers 1–5, from an investigation with detailed analysis of the spectral assignments. The “resonance condition” is calculated as was done in Table 1.
NumberAssignment *IR Vibration Frequency, cm−1Resonance Condition
(R. C.)
Ideal ValueDeviation Δ (o/oo)
1ν (N9-H)30160.689--
2δ (N9-H) + 7 extra terms16031.2971.33327
3δ (N9-H)
+ 10 extra terms
15081.3791.33335
4δ (N9-H)
+ 8 extra terms
14501.4341.50044
5δ (N9-H)
+ 9 extra terms
13671.5211.50014
* The “extra terms” in column 2 (mostly stretching vibration but also some bending ones) are given [23], but with no quantitative estimations.
Table 8. The influence of the catalyst COS (ω = 2079 cm−1) on the reactions mentioned here, indicating the type of interaction.
Table 8. The influence of the catalyst COS (ω = 2079 cm−1) on the reactions mentioned here, indicating the type of interaction.
ReactantCritical
Vibration
cm−1
Route of
Attack
ω/νDeviation
o/oo
ProductReference
H2COν2 = 1743.6Double bond splitting6:56Ribose
C5H5O5
[10]
Table 2
HCNν1 = 2089Triple bond splitting1:15AdenineScheme 3 and Scheme 4
CO2ν 6 + ν 10 = 2076.5Vibration energy transfer1:11CO + O[10]
Figure 6
NH3ν 2(a) + 2 ν 4(e) = 4176Vibration energy transfer1:24NH + H2[10]
Figure 6
C2H4ν 9CH(b2u) = 3105.5Double bond splitting2:34C2H2 + H2Table 2
(in [10])
PO43−
glucose
1045
1048
Molecular Interaction2:1
2:1
1048/1045 = 1.0048
5 o/oo
Phospho-glucose-[4]
Table 5
PO43−
ribose
1045
1040
Molecular interaction2:1
2 = 1
1045/1040 = 1.0048
5 o/oo
Phospho-ribose[19,20]
Table 6
Ribose
adenine
1040
1386
Molecular interaction2:1
3:2
2/1.5 = 4:3
1386/1040 = 1.327
5 o/oo
Adenine-ribose[22]
Table 7
Amino-
acetate
(Cu2+)
COS
1058
2079
Molecular
interaction
1:2
1:1
2 × 1058 = 2116
2116/2079 = 1.0178
9 o/oo
-NH-CO-(Condensation of amino acids)[1]

Share and Cite

MDPI and ACS Style

Larsson, R.; Malek, A.; Odenbrand, I. The Transformation by Catalysis of Prebiotic Chemical Systems to Useful Biochemicals: A Perspective Based on IR Spectroscopy of the Primary Chemicals II. Catalysis and the Building of RNA. Appl. Sci. 2020, 10, 4712. https://doi.org/10.3390/app10144712

AMA Style

Larsson R, Malek A, Odenbrand I. The Transformation by Catalysis of Prebiotic Chemical Systems to Useful Biochemicals: A Perspective Based on IR Spectroscopy of the Primary Chemicals II. Catalysis and the Building of RNA. Applied Sciences. 2020; 10(14):4712. https://doi.org/10.3390/app10144712

Chicago/Turabian Style

Larsson, Ragnar, Abdul Malek, and Ingemar Odenbrand. 2020. "The Transformation by Catalysis of Prebiotic Chemical Systems to Useful Biochemicals: A Perspective Based on IR Spectroscopy of the Primary Chemicals II. Catalysis and the Building of RNA" Applied Sciences 10, no. 14: 4712. https://doi.org/10.3390/app10144712

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop