Aluminosilicate surfaces as promoters for peptide bond formation: an assessment of Bernal's hypothesis by ab initio methods

Amino Acid Polymerization on Silica Surfaces

The polymerization of unactivated amino acids (AAs) is an important topic because of its applications in various fields including industrial medicinal chemistry and prebiotic chemistry. Silica as a promoter for this reaction, is of great interest owing to its large abundance and low cost. The amide/peptide bond synthesis on silica has been largely demonstrated but suffers from a lack of knowledge regarding its reaction mechanism, the key parameters, and surface features that influence AA adsorption and reactivity, the selectivity of the reaction product, the role of water in the reaction, etc. The present review addresses these problems by summarizing experimental and modeling results from the literature and attempts to rationalize some apparent divergences in published results. After briefly presenting the main types of silica surface sites and other relevant macroscopic features, we discuss the different deposition procedures of AAs, whose importance is often neglected. We address the possible AA adsorption mechanisms including covalent grafting and H-bonding and show that they are highly dependent on silanol types and density. We then consider how the adsorption mechanisms determine the occurrence and outcome of AA condensation (formation of cyclic dimers or of long linear chains), and outline some recent results that suggest significant polymerization selectivity in systems containing several AAs, as well as the formation of specific elements of secondary structure in the growing polypeptide chains.

Enantiodetermining processes in the synthesis of alanine, serine, and isovaline

In this study, quantum chemical calculations were used to explore the synthesis of three chiral α-amino acids, specifically alanine, serine, and isovaline, from reactants found in interstellar space. Our focus is on the crucial step in the synthesis pathway that determines the chirality of the amino acids. The results indicate that in the case of alanine, the determination of enantiomer is primarily influenced by the direction of the collision of molecules or functional groups, which leads to the formation of a chirality center in a crucial intermediate. However, contrary to chemical expectations, the enantiodetermining/enantioselection step for serine and isovaline synthesis occurs prior to the creation of a chirality center.

The Formation of RNA Pre-Polymers in the Presence of Different Prebiotic Mineral Surfaces Studied by Molecular Dynamics Simulations

We used all-atom Molecular Dynamics (MD) computer simulations to study the formation of pre-polymers between the four nucleotides in RNA (AMP, UMP, CMP, GMP) in the presence of different substrates that could have been present in a prebiotic environment. Pre-polymers are C3'-C5' hydrogen-bonded nucleotides that have been suggested to be the precursors of phosphodiester-bonded RNA polymers. We simulated wet-dry cycles by successively removing water molecules from the simulations, from ~60 to 3 water molecules per nucleotide. The nine substrates in this study include three clay minerals, one mica, one phosphate mineral, one silica, and two metal oxides. The substrates differ in their surface charge and ability to form hydrogen bonds with the nucleotides. From the MD simulations, we quantify the interactions between different nucleotides, and between nucleotides and substrates. For comparison, we included graphite as an inert substrate, which is not charged and cannot form hydrogen bonds. We also simulated the dehydration of a nucleotide-only system, which mimics the drying of small droplets. The number of hydrogen bonds between nucleotides and nucleotides and substrates was found to increase significantly when water molecules were removed from the systems. The largest number of C3'-C5' hydrogen bonds between nucleotides occurred in the graphite and nucleotide-only systems. While the surface of the substrates led to an organization and periodic arrangement of the nucleotides, none of the substrates was found to be a catalyst for pre-polymer formation, neither at full hydration, nor when dehydrated. While confinement and dehydration seem to be the main drivers for hydrogen bond formation, substrate interactions reduced the interactions between nucleotides in all cases. Our findings suggest that small supersaturated water droplets that could have been produced by geysers or springs on the primitive Earth may play an important role in non-enzymatic RNA polymerization.

Tracing the Primordial Chemical Life of Glycine: A Review from Quantum Chemical Simulations

Glycine (Gly), NH₂CH₂COOH, is the simplest amino acid. Although it has not been directly detected in the interstellar gas-phase medium, it has been identified in comets and meteorites, and its synthesis in these environments has been simulated in terrestrial laboratory experiments. Likewise, condensation of Gly to form peptides in scenarios resembling those present in a primordial Earth has been demonstrated experimentally. Thus, Gly is a paradigmatic system for biomolecular building blocks to investigate how they can be synthesized in astrophysical environments, transported and delivered by fragments of asteroids (meteorites, once they land on Earth) and comets (interplanetary dust particles that land on Earth) to the primitive Earth, and there react to form biopolymers as a step towards the emergence of life. Quantum chemical investigations addressing these Gly-related events have been performed, providing fundamental atomic-scale information and quantitative energetic data. However, they are spread in the literature and difficult to harmonize in a consistent way due to different computational chemistry methodologies and model systems. This review aims to collect the work done so far to characterize, at a quantum mechanical level, the chemical life of Gly, i.e., from its synthesis in the interstellar medium up to its polymerization on Earth.

Canonical, deprotonated, or zwitterionic? II. A computational study on amino acid interaction with the TiO2(110) rutile surface: comparison with the anatase (101) surface

The adsorption of 11 amino acids (Gly, Leu, Met, Phe, Ser, Cys, Glu, Gln, Arg, Lys, and His) on the TiO₂(110) rutile surface is investigated adopting a theoretical approach, using the PBE-D2* functional as implemented in the periodic VASP code. The adsorption of the amino acids is considered in their canonical, deprotonated and zwitterionic forms. For all cases, the most stable adsorption mode adopts a bidentate (O,O) binding with surface undercoordinated Ti atoms, in agreement with previous experimental and computational studies using glycine as a test case. Such a binding mode is possible due to the surface morphology, because the Ti-Ti distances match very well with the carboxylic O-O distance. The most stable adsorption states are the deprotonated and the zwitterionic ones, the canonical one lying significantly above in energy. The relative stability between the deprotonated and the zwitterionic states results in a delicate trade-off among dative interactions (O, N, and S atoms of the amino acids with Ti atoms of the surface), H-bond interactions, dispersive forces and, to a lesser extent, steric hindrance of the amino acidic lateral chains. Finally, the difference in the amino acid adsorption between the (110) rutile and the (101) anatase surfaces is discussed both from the energetic and surface morphological standpoints, highlighting the larger reactivity of the rutile polymorph in adsorbing and deprotonating the amino acids compared with the anatase one.

A quantum mechanical study of dehydration vs. decarbonylation of formamide catalysed by amorphous silica surfaces

Formamide is abundant in the interstellar medium and was also present during the formation of the Solar system through the accretion process of interstellar dust. Under the physicochemical conditions of primordial Earth, formamide could have undergone decomposition, either via dehydration (HCN + H2O) or via decarbonylation (CO + NH3). The first reactive channel provides HCN, which is an essential molecular building block for the formation of RNA/DNA bases, crucial for the emergence of life on Earth. In this work, we studied, at the CCSD(T)/cc-pVTZ level, the two competitive routes of formamide decomposition, i.e. dehydration and decarbonylation, either in liquid formamide (by using the polarization continuum model technique) or at the interface between liquid formamide and amorphous silica. Amorphous silica was adopted as a convenient model of the crystalline silica phases ubiquitously present in the primordial (and actual) Earth's crust, and also due to its relevance in catalysis, adsorption and chromatography. Results show that: (i) silica surface sites catalyse both decomposition channels by reducing the activation barriers by about 100 kJ mol-1 with respect to the reactions in homogeneous medium, and (ii) the dehydration channel, giving rise to HCN, is strongly favoured from a kinetic standpoint over decarbonylation, the latter being, instead, slightly favoured from a thermodynamic point of view.

Experimental and theoretical elucidation of catalytic pathways in TiO2-initiated prebiotic polymerization

The tendency of glycine to form polymer chains on a rutile(110) surface under wet/dry conditions (dry-wet cycles at high temperature) is studied through a conjunction of surface sensitive experimental techniques and sequential periodic multilevel calculations that mimics the experimental procedures with models of decreasing complexity and increasing accuracy. X-ray photoemission spectroscopy (XPS) and thermal desorption spectroscopy (TDS) experimentally confirmed that the dry-wet cycles lead to Gly polymerization on the oxide support. This was supported by all the theoretical characterizations. First, classical reactive molecular dynamics (MD) simulations based on the ReaxFF approach were used to reproduce the adsorption of the experimental glycine solution droplets sprayed onto an oxide support and to identify the most probable arrangement of the molecules that triggered the polymerization mechanisms. Then, quantum chemistry density functional tight binding (DF-TB) MDs and static density functional theory (DFT) calculations were carried out to further explore favorable configurations and to evaluate the energy barriers of the most promising reaction pathways for the peptide bond-formation reactions. The results confirmed the fundamental role played by the substrate to thermodynamically and kinetically favor the process and disclosed its main function as an immobilizing agent: the molecules accommodated in the surface channels close to each other were the ones starting the key events of the dimerization process and the most favorable mechanism was the one where a water molecule acted as a proton exchange mediator in the condensation process.

Role of Mineral Surfaces in Prebiotic Chemical Evolution. In Silico Quantum Mechanical Studies

There is a consensus that the interaction of organic molecules with the surfaces of naturally-occurring minerals might have played a crucial role in chemical evolution and complexification in a prebiotic era. The hurdle of an overly diluted primordial soup occurring in the free ocean may have been overcome by the adsorption and concentration of relevant molecules on the surface of abundant minerals at the sea shore. Specific organic⁻mineral interactions could, at the same time, organize adsorbed molecules in well-defined orientations and activate them toward chemical reactions, bringing to an increase in chemical complexity. As experimental approaches cannot easily provide details at atomic resolution, the role of in silico computer simulations may fill that gap by providing structures and reactive energy profiles at the organic⁻mineral interface regions. Accordingly, numerous computational studies devoted to prebiotic chemical evolution induced by organic⁻mineral interactions have been proposed. The present article aims at reviewing recent in silico works, mainly focusing on prebiotic processes occurring on the mineral surfaces of clays, iron sulfides, titanium dioxide, and silica and silicates simulated through quantum mechanical methods based on the density functional theory (DFT). The DFT is the most accurate way in which chemists may address the behavior of the molecular world through large models mimicking chemical complexity. A perspective on possible future scenarios of research using in silico techniques is finally proposed.

When the Surface Matters: Prebiotic Peptide-Bond Formation on the TiO2 (101) Anatase Surface through Periodic DFT-D2 Simulations

The mechanism of the peptide-bond formation between two glycine (Gly) molecules has been investigated by means of PBE-D2* and PBE0-D2* periodic simulations on the TiO₂ (101) anatase surface. This is a process of great relevance both in fundamental prebiotic chemistry, as the reaction univocally belongs to one of the different organizational events that ultimately led to the emergence of life on Earth, as well as from an industrial perspective, since formation of amides is a key reaction for pharmaceutical companies. The efficiency of the surface catalytic sites is demonstrated by comparing the reactions in the gas phase and on the surface. At variance with the uncatalyzed gas-phase reaction, which involves a concerted nucleophilic attack and dehydration step, on the surface these two steps occur along a stepwise mechanism. The presence of surface Lewis and Brönsted sites exerts some catalytic effect by lowering the free energy barrier for the peptide-bond formation by about 6 kcal mol^-1 compared to the gas-phase reaction. Moreover, the co-presence of molecules acting as proton-transfer assistants (i.e., H₂ O and Gly) provide a more significant kinetic energy barrier decrease. The reaction on the surface is also favorable from a thermodynamic standpoint, involving very large and negative reaction energies. This is due to the fact that the anatase surface also acts as a dehydration agent during the condensation reaction, since the outermost coordinatively unsaturated Ti atoms strongly anchor the released water molecules. Our theoretical results provide a comprehensive atomistic interpretation of the experimental results of Martra et al. (Angew. Chem. Int. Ed. 2014, 53, 4671), in which polyglycine formation was obtained by successive feedings of Gly vapor on TiO₂ surfaces in dry conditions and are, therefore, relevant in a prebiotic context envisaging dry and wet cycles occurring, at mineral surfaces, in a small pool.

Thermal Condensation of Glycine and Alanine on Metal Ferrite Surface: Primitive Peptide Bond Formation Scenario

The amino acid condensation reaction on a heterogeneous mineral surface has been regarded as one of the important pathways for peptide bond formation. Keeping this in view, we have studied the oligomerization of the simple amino acids, glycine and alanine, on nickel ferrite (NiFe₂O₄), cobalt ferrite (CoFe₂O₄), copper ferrite (CuFe₂O₄), zinc ferrite (ZnFe₂O₄), and manganese ferrite (MnFe₂O₄) nanoparticles surfaces, in the temperature range from 50-120 °C for 1-35 days, without applying any wetting/drying cycles. Among the metal ferrites tested for their catalytic activity, NiFe₂O₄ produced the highest yield of products by oligomerizing glycine to the trimer level and alanine to the dimer level, whereas MnFe₂O₄ was the least efficient catalyst, producing the lowest yield of products, as well as shorter oligomers of amino acids under the same set of experimental conditions. It produced primarily diketopiperazine (Ala) with a trace amount of alanine dimer from alanine condensation, while glycine was oligomerized to the dimer level. The trend in product formation is in accordance with the surface area of the minerals used. A temperature as low as 50 °C can even favor peptide bond formation in the present study, which is important in the sense that the condensation process is highly feasible without any sort of localized heat that may originate from volcanoes or hydrothermal vents. However, at a high temperature of 120 °C, anhydrides of glycine and alanine formation are favored, while the optimum temperature for the highest yield of product formation was found to be 90 °C.

Abiogenic Syntheses of Lipoamino Acids and Lipopeptides and their Prebiotic Significance

Researchers have formed peptide bonds under a variety of presumed prebiotic conditions. Here it is proposed that these same conditions would have also formed amide bonds between fatty acids and amino acids, producing phosphate-free amphipathic lipoamino acids and lipopeptides. These compounds are known to form vesicles and are ubiquitous in living organisms. They could represent molecules that provided protection by membranes as well as possibilities for proto-life metabolism . It is here demonstrated that when a fatty acid is heated with various amino acids, optimally in the presence of suitable salts or minerals, lipoamino acids are formed. Magnesium and potassium carbonates as well as iron (II) sulfide are found to be particularly useful in these reactions. In this manner N-lauroylglycine, N-lauroylalanine, N-stearoylalanine and several other lipoamino acids have been synthesized. Similarly, when glycylglycine was heated with lauric acid in the presence of magnesium carbonate, the lipopeptide N-lauroylglycylglycine was formed. Such compounds are proposed to have been critical precursors to the development of life.

Attachment of ribonucleotides on α-alumina as a function of pH, ionic strength, and surface loading

The interactions between nucleic acids and mineral surfaces have been the focus of many studies in environmental sciences, in biomedicine, as well as in origin of life studies for the prebiotic formation of biopolymers. However, few studies have focused on a wide range of environmental conditions and the likely modes of attachment. Here we investigated the adsorption of ribonucleotides onto α-alumina surfaces over a wide range of pH, ionic strength, and ligand-to-solid ratio, by both an experimental and a theoretical approach. The adsorption of ribonucleotides is strongly affected by pH, with a maximum adsorption at pH values around 5. Alumina adsorbs high amounts of nucleotides >2 μmol/m(2). We used the extended triple-layer model (ETLM) to predict the speciation of the surface complexes formed as well as the stoichiometry and equilibrium constants. We propose the formation of two surface species: a monodentate inner-sphere complex, dominant at pH <7, and a bidentate outer-sphere complex, dominant at higher pH. Both complexes would involve interactions between the negatively charged phosphate group and the positively charged surface of alumina. Our results provide a better understanding of how nucleic acids attach to mineral surfaces under varying environmental conditions. Moreover, the predicted configuration of nucleotide surface species, bound via the phosphate group, could have implications for the abiotic formation of nucleic acids in the context of the origin of life.

Miller experiments in atomistic computer simulations

The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.

Peptide formation mechanism on montmorillonite under thermal conditions

The oligomerization of amino acids is an essential process in the chemical evolution of proteins, which are precursors to life on Earth. Although some researchers have observed peptide formation on clay mineral surfaces, the mechanism of peptide bond formation on the clay mineral surface has not been clarified. In this study, the thermal behavior of glycine (Gly) adsorbed on montmorillonite was observed during heating experiments conducted at 150 °C for 336 h under dry, wet, and dry-wet conditions to clarify the mechanism. Approximately 13.9 % of the Gly monomers became peptides on montmorillonite under dry conditions, with diketopiperazine (cyclic dimer) being the main product. On the other hand, peptides were not synthesized in the absence of montmorillonite. Results of IR analysis showed that the Gly monomer was mainly adsorbed via hydrogen bonding between the positively charged amino groups and negatively charged surface sites (i.e., Lewis base sites) on the montmorillonite surface, indicating that the Lewis base site acts as a catalyst for peptide formation. In contrast, peptides were not detected on montmorillonite heated under wet conditions, since excess water shifted the equilibrium towards hydrolysis of the peptides. The presence of water is likely to control thermodynamic peptide production, and clay minerals, especially those with electrophilic defect sites, seem to act as a kinetic catalyst for the peptide formation reaction.

Structural and energetic factors controlling the enantioselectivity of dinucleotide formation under prebiotic conditions

Recently, it has been reported that the montmorillonite-catalyzed oligomerization of activated nucleotides exhibits remarkable enantioselectivity. In the current paper we investigate the structures and intrinsic energies of homochiral and heterochiral cyclic dinucleotides by means of accurate quantum chemical calculations in gas-phase and in bulk water. The steric effect of the clay is represented with geometrical constraints. Our computations reveal that the heterochiral dimer geometries are systematically less stable than their homochiral counterparts due to steric clashes inside the sugar-phosphate ring geometry. Thus we suggest that the homochiral selectivity observed in the cyclic dinucleotide formation in confined spaces may arise from the energetic destabilization of the heterochiral ring geometries as compared to their homochiral analogues. In the present paper we provide the first model of the 3D structure of d,l cyclic dinucleotides, which until now has eluded experimental observation.

Pathways for the formation and evolution of peptides in prebiotic environments

α-Amino acids are easily accessible through abiotic processes and were likely present before the emergence of life. However, the role they could have played in the process remains uncertain. Chemical pathways that could have brought about features of self-organization in a peptide world are considered in this review and discussed in relation with their possible contribution to the origin of life. An overall scheme is proposed with an emphasis on possibilities that may have led to dynamically stable far from equilibrium states. This analysis defines new lines of investigation towards a better understanding of the contribution of the systems chemistry of amino acids and peptides to the emergence of life.

Glycine peptide bond formation catalyzed by faujasite

The catalysis of peptide bond formation between two glycine molecules on H-FAU zeolite was computationally studied by the M08-HX density functional. Two reaction pathways, the concerted and the stepwise mechanism, starting from three differently adsorbed reactants, amino-bound, carboxyl-bound, and hydroxyl-bound, are studied. Adsorption energies, activation energies, and reaction energies, as well as the corresponding intrinsic rate constants were calculated. A comparison of the computed energetics of the various reaction paths for glycine indicates that the catalyzed reaction proceeds preferentially via the concerted reaction mechanism of the hydroxyl-bound configuration. This involves an eight-membered ring of the transition structure instead of the four-membered ring of the others. The step from the amino-bound configuration to glycylglycine is the rate-determining step of the concerted mechanism. It has an estimated activation energy of 51.2 kcal mol(-1). Although the catalytic reaction can also occur via the stepwise reaction mechanism, this path is not favored.

Possible origin of life between mica sheets

The mica hypothesis is a new hypothesis about how life might have originated. The mica hypothesis provides simple solutions to many basic questions about the origins of life. In the mica hypothesis, the spaces between mica sheets functioned as the earliest cells. These 'cells' between mica sheets are filled with potassium ions, and they provide an environment in which: polymer entropy is low; cyclic wetting and drying can occur; molecules can evolve in isolated spaces and also migrate and ligate to form larger molecules. The mica hypothesis also proposes that mechanical energy (work) is a major energy source that could have been used on many length scales to form covalent bonds, to alter polymer conformations, and to bleb daughter cells off protocells. The mica hypothesis is consistent with many other origins hypotheses, including the RNA, lipid, and metabolic 'worlds'. Therefore the mica hypothesis has the potential to unify origins hypotheses, such that different molecular components and systems could simultaneously evolve in the spaces between mica sheets.

Formation versus hydrolysis of the peptide bond from a quantum-mechanical viewpoint: The role of mineral surfaces and implications for the origin of life

The condensation (polymerization by water elimination) of molecular building blocks to yield the first active biopolymers (e.g. of amino acids to form peptides) during primitive Earth is an intriguing question that nowadays still remains open since these processes are thermodynamically disfavoured in highly dilute water solutions. In the present contribution, formation and hydrolysis of glycine oligopeptides occurring on a cluster model of sanidine feldspar (001) surface have been simulated by quantum mechanical methods. Results indicate that the catalytic interplay between Lewis and Brønsted sites both present at the sanidine surface, in cooperation with the London forces acting between the biomolecules and the inorganic surface, plays a crucial role to: i) favour the condensation of glycine to yield oligopeptides as reaction products; ii) inhibit the hydrolysis of the newly formed oligopeptides. Both facts suggest that mineral surfaces may have helped in catalyzing, stabilizing and protecting from hydration the oligopeptides formed in the prebiotic era.

Neutral vs zwitterionic glycine forms at the water/silica interface: structure, energies, and vibrational features from B3LYP periodic simulations

B3LYP periodic calculations with a triple-xi-polarized Gaussian basis set have been used to study adsorption of glycine on a hydroxylated silica surface (2.2 OH/nm2) model derived from the (001) surface of edingtonite. The simulation envisages glycine adsorbed either as a gas-phase molecule or when microsolvated by up to five H20 molecules. Both neutral and zwitterionic forms of glycine have been considered and their structural, energetic, and spectroscopic vibrational features compared internally and with experiments. As a gas phase glycine sticks in its neutral form at the silica surface, the zwitterion being highly unstable and with transition-state character. When glycine is microsolvated at the silica interface, two H20 molecules render the zwitterion population comparable to that of the neutral form whereas with four H2O molecules the neutral glycine population is wiped out in favor of the zwitterion. With four H20 molecules the most stable structure shows no direct contact between glycine and the silica surface, H20 acting as a mediator via H-bond interactions. The B3LYP energies and structural data were also supported by comparing the scaled harmonic vibrational features with literature FTIR data of glycine adsorbed on an amorphous silica surface either from the gas phase or in water solution.

Adsorption and polymerization of amino acids on mineral surfaces: a review

The present paper offers a review of recent (post-1980) work on amino acid adsorption and thermal reactivity on oxide and sulfide minerals. This review is performed in the general frame of evaluating Bernal's hypothesis of prebiotic polymerization in the adsorbed state, but written from a surface scientist's point of view. After a general discussion of the thermodynamics of the problem and exactly what effects surfaces should have to make adsorbed-state polymerization a viable scenario, we examine some practical difficulties in experimental design and their bearing on the conclusions that can be drawn from extant works, including the relevance of the various available characterization techniques. We then present the state of the art concerning the mechanisms of the interactions of amino acids with mineral surfaces, including results from prebiotic chemistry-oriented studies, but also from several different fields of application, and discuss the likely consequences for adsorption selectivities. Finally, we briefly summarize the data concerning thermally activated amide bond formation of adsorbed amino acids without activating agents. The reality of the phenomenon is established beyond any doubt, but our understanding of its mechanism and therefore of its prebiotic potential is very fragmentary. The review concludes with a discussion of future work needed to fill the most conspicuous gaps in our knowledge of amino acids/mineral surfaces systems and their reactivity.