Peptide bond formation activated by the interplay of Lewis and Brønsted catalysts

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.

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.

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.

Water-assisted peptide bond formation between two double amino acid molecules in the gas phase

The gas phase mechanism of the peptide bond formation between two double amino acid (DAA) molecules described by the (NH₂)₂C(COOH)₂ formula is investigated in the presence of a water molecule. Formations of trans and cis DAA-DAA dipeptide products along both concerted and stepwise mechanisms have been studied at the CCSD(T)/aug-cc-pVDZ//MP2/aug-cc-pVDZ level. The results indicate that the activation energy barriers estimated for the water-assisted mechanisms are significantly reduced in comparison to the corresponding uncatalyzed reactions. The trans DAA-DAA isomer is expected to dominate in the final product due to its larger stability compared to the cis DAA-DAA product.

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.

Computational study of peptide bond formation in the gas phase through ion-molecule reactions

A computational study of peptide bond formation from gas-phase ion-molecule reactions has been carried out. We have considered the reaction between protonated glycine and neutral glycine, as well as the reaction between two neutral glycine molecules for comparison purposes. Two different mechanisms, concerted and stepwise, were studied. Both mechanisms show significant energy barriers for the neutral reaction. The energy requirements for peptide bond formation are considerably reduced upon protonation of one of the glycine molecules. For the reaction between neutral glycine and N-protonated glycine the lowest energy barrier is observed for the concerted mechanism. For the reaction between neutral glycine and protonated glycine at carbonyl oxygen, the preferred mechanism is the stepwise one, with a relatively small energy barrier (23 kJ mol(-1) at 0 K) and leading to the lowest-lying protonated glycylglycine isomer. In the case that the reaction could be initiated by protonated glycine at hydroxyl oxygen the process would be barrier-free and clearly exothermic. In that case peptide bond formation could take place even under interstellar conditions if glycine is present in space.

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.

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.

Clean, reusable and low cost heterogeneous catalyst for amide synthesis

We have developed a heterogeneous silica catalyst that can effectively catalyse amide synthesis from acid and amine, without production of toxic by-products and with the advantage of being readily available, low cost, environmentally benign and reusable.

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.

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.

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

The role in prebiotic chemistry that Brønsted and Lewis sites, both present at the surface of common aluminosilicates, may have played in favoring the peptide bond formation has been addressed by ab initio methods within a cluster approach. B3LYP/6-31+G(d,p) free energy potential energy surfaces have been fully characterized for the model reaction glycine + NH3 --> 2-NH2 acetamide (mimicking the true 2 Gly --> GlyGly one) occurring on (i) a Lewis site, (ii) a Brønsted site, and (iii) a combined action of Lewis/Brønsted sites. Compared to the gas-phase (gp) activation free energy of 50 kcal/mol, the Lewis site alone reduces the gp barrier to 41 kcal/mol, whereas the activation by the Brønsted site dramatically reduces the barrier to about 18 kcal/mol. Nevertheless, formation of the prereactant complex in this latter case will rarely occur, since water will easily displace the glycine molecule interacting with the Brønsted site. However, if a realistic feldspar surface with neighboring Brønsted and Lewis sites is considered, the proper prereactant complex is highly stabilized by a simultaneous interaction with the Lewis and the Brønsted sites, in such a way that the Lewis site strongly attaches the glycine molecule to the surface whereas the Brønsted site efficiently catalyzes the condensation reaction, showing that the interplay between Lewis/Brønsted sites is an important issue. The free energy barrier computed for the realistic feldspar surface model is 26 kcal/mol. The role of dispersive interactions on the free energy barrier and the stabilization of the final product, not accounted for by the B3LYP functional, have been estimated and shown to be substantial. Speculations about further elongation of the formed dipeptide have been put forward on the basis of the relatively strong interaction energy of the formed GlyGly dipeptide with the aluminosilicate surface.

Is the Peptide Bond Formation Activated by Cu2+ Interactions? Insights from Density Functional Calculations

The catalytic role that Cu(2+) cations play in the peptide bond formation has been addressed by means of density functional calculations. First, the Cu(2+)-(glycine)2 --> Cu(2+)-(glycylglycine) + H2O reaction was investigated since mass spectrometry low collision activated dissociation (CAD) spectra of Cu(2+)-(glycine)2 led to the elimination of a water molecule, which suggested that an intracomplex peptide bond formation might have occurred. Results show that this intracomplex condensation is associated to a very high free energy barrier (97 kcal mol(-1)) and reaction free energy (66 kcal mol(-1)) because of the loss of metal coordination during the reaction. Second, on the basis of the salt-induced peptide formation theory, the condensation reaction between two glycines was studied in aqueous solution using discrete water molecules and the conductor polarized continuum model (CPCM) continuous method. It is found that the synergy between the interaction of glycines with Cu(2+) and the presence of water molecules acting as proton-transfer helpers significantly lower the activation barrier (from 55 kcal/mol for the uncatalyzed system to 20 kcal/mol for the Cu(2+) solvated system) which largely favors the formation of the peptide bond.

Interaction of glycine with isolated hydroxyl groups at the silica surface: first principles B3LYP periodic simulation

The adsorption of a glycine molecule on a model silica surface terminated by an isolated hydroxyl group has been studied ab initio using a double-zeta polarized Gaussian basis set, the hybrid B3LYP functional, and a full periodic treatment of the silica surface/glycine system. The hydroxylated silica surface has been simulated using either a 2D slab or a single polymer strand cut out from the (001) surface of an all-silica edingtonite. A number of B3LYP-optimized structures have been found by docking glycine on the silica surface exploiting all possible hydrogen bond patterns. Whereas glycine is generally adsorbed in its neutral form, two structures show glycine adsorbed as a zwitterion, the surface playing the role of a "solid solvent" whereas intrastrand hydrogen bond cooperativity stabilizes the zwitterions. The adsorbed zwitterionic structures are no longer formed at a lower glycine coverage as simulated by enlarging the unit cell so as to break intrastrand hydrogen bonds, showing the importance of H-bond cooperativity in stabilizing the zwitterionic forms. Each structure has been characterized by computing its harmonic vibrational spectrum at the Gamma point, which also allowed us to calculate the free energy of adsorption. The experimental infrared features of chemical-vapor-deposited glycine on a silica surface are in agreement with those computed for glycine adsorbed in its neutral form and engaging three hydrogen bonds with the surface silanols, two of them involving the C=O bond and one originating from the glycine OH group. The NH(2) group plays only a minor role as a weak hydrogen bond donor.

Imidazo[1,2-a]pyrazine-3,6-diones derived from alpha-amino acids: a theoretical mechanistic study of their formation via pyrolysis and silica-catalyzed process

Imidazo[1,2-a]pyrazine-3,6-diones are unusual compounds composed of three alpha-amino acid fragments. These bicyclic amidines (BCAs) form under high temperatures or with the use of strong dehydrating reagents. We gave insight into the mechanisms of BCA formation via gas-phase pyrolytic and silica-catalyzed reactions of glycine (Gly) and alpha-aminoisobutyric acid (AIB) with related diketopiperazines (DKPs), using quantum chemical calculations. The entire process requires four steps: (1) O-acylation of DKP with free or silica-bonded amino acid, (2) acyl transfer from the oxygen to the nitrogen atom, (3) intramolecular condensation of the N-acyl DKP into a cyclol, and (4) elimination of water. To study step (1) at silica surface (modeled by H7Si8O12-OH cluster), we employed two-level ONIOM calculations (AM1:UFF, B3LYP/3-21G:UFF and B3LYP/6-31G(d):UFF); all gas-phase reactions were studied at the AM1, B3LYP/3-21G and B3LYP/6-31G(d) levels. The catalytic effect of silica was observed for both Gly and AIB: the activation energy in the O-acylation at the surface was lower by more than 9 kcal mol(-1) as compared to the gas-phase process. Contrary to the exothermic O-acylation, the gas-phase transfer reaction (step 2) was exothermic in both cases, but more favorable for Gly. The cyclocondensation of N-acylated DKPs into BCAs (steps 3 and 4) is endothermic for Gly and exothermic for AIB.

Does Silica Surface Catalyse Peptide Bond Formation? New Insights from First-Principles Calculations

The role that silica surface could have played in prebiotic chemistry as a catalyst for peptide bond formation has been addressed at the B3LYP/6-31+G(d,p) level for a model reaction involving glycine and ammonia on a silica cluster mimicking an isolated terminal silanol group present at the silica surface. Hydrogen-bond complexation between glycine and the silanol is followed by the formation of the mixed surface anhydride Si(surf)-O-C(=O)-R, which has been suggested in the literature to activate the C=O bond towards nucleophilic attack by a second glycine molecule, here simulated by the simpler NH3 molecule. However, B3LYP/6-31+G(d,p) calculations show that formation of the surface mixed anhydride Si(surf)-O-C(=O)-R is disfavoured (delta(r)G298 approximately 6 kcal mol(-1)), and that the surface bond only moderately lowers the free-energy barrier of the nucleophilic attack responsible for peptide bond formation (deltaG298(double dagger) approximately 48 kcal mol(-1)) in comparison with the uncatalysed reaction (deltaG298(double dagger) approximately 52 kcal mol(-1)). A further decrease of the free-energy barrier of peptide bond formation (deltaG298(double dagger) approximately 41 kcal mol(-1)) is achieved by a single water molecule close to the reaction centre acting as a proton-transfer helper in the activated complex. A possible role of strained silica surface defects on the formation of the surface mixed anhydride Si(surf)-O-C(=O)-R has also been addressed.