1. Introduction
Given the role of peptides, and eventually proteins, in performing catalytic and structural functions, the formation of amino acid polymers from their monomers was a critical step for the evolution of life on Earth. Peptides are formed via condensation polymerization, where the amino group of one amino acid reacts with the carboxylic acid group of a second amino acid, forming an amide bond and releasing a water molecule in the process. Because water is a by-product of polymerization, uncatalyzed peptide bond formation is thermodynamically unfavorable in the aqueous environments predicted to have existed on prebiotic Earth [
1]. This is further compounded by the fact that the reverse reaction of condensation is hydrolysis, which is favored when water activity is high [
2]. Several mechanisms have been proposed to provide a prebiotically plausible route to peptide formation under early Earth conditions, including mineral catalysts [
3,
4,
5,
6], condensation reagents [
7,
8], metal ions [
9] or environmentally mediated dehydration such as wet–dry cycling [
10,
11].
In scenarios invoking dehydrating conditions, salts reduce water activity and facilitate polymerization, thus reducing the thermodynamic barrier to polymerization in aqueous environments. One scenario that has been investigated extensively is the salt-induced peptide formation (SIPF) reaction [
12,
13,
14,
15,
16,
17,
18,
19]. In laboratory experiments, the SIPF reaction medium is a salt solution containing monovalent and divalent metal ions (most commonly, NaCl and CuCl
2, respectively). The Na, Cl, and Cu ions aid in reducing water activity, thus helping to drive the reaction toward peptide formation, while the divalent Cu additionally lowers the activation energy for peptide formation. The Cu ion readily forms complexes with multiple α-amino acids, bringing them into close proximity and activating reactive moieties. The originally proposed reactive species consisted of one neutral amino acid chelated with the copper–chloride complex, while a protonated second amino acid binds end-on via its carboxyl group, and two water molecules remain bound to the complex, one above, one below (hereafter referred to as the protonated–neutral complex; see
Figure 1A) [
17]. Subsequent work has also identified a bi-bidentally bound CuCl–amino acid complex between two neutral amino acids species as a potential reactive species (hereafter referred to as the neutral–neutral species,
Figure 1B) [
20]. The formation of these complexes mediates amide bond formation between the amino acids via a nucleophilic attack by the amine group of the bidentate amino acid (amino acid 2) on the carboxyl carbon of amino acid 1. The complexation of the Cu (II) is weaker than the newly formed amide bond, and the Cu (II) ion is released back into solution. Further elongation of the peptide occurs via the same mechanism, with the nascent peptide binding end-on via its carboxyl group and a free amino acid chelating the copper–chloride complex; one amino acid can be added per polymerization cycle.
The requirements to facilitate the SIPF reaction are feasible in an early Earth environment:
- (1)
Amino Acids. Amino acids either formed in situ in prebiotic environments [
21,
22,
23] or delivered via meteorites [
24,
25,
26].
- (2)
Divalent Cations. Although other metal ions can catalyze the reaction (e.g., Cr (III), Al (III), Mg (II)), divalent copper confers the best reactivity [
17,
18]. Cu-bearing deposits (known as ‘greenstones’) found in Precambrian rock formations indicate that sufficient Cu(II) was present to facilitate the SIPF reaction in prebiotic environments [
27].
- (3)
Chloride Salts. NaCl was likely present in prebiotic environments, either leached from the Earth’s crust, which contains 130–1910 ppm of chlorine and 2.0–2.9 wt% sodium [
28], into freshwater sources, or were sourced from water entering the reaction environments from oceans or lagoons.
- (4)
Temperature. Experimental work has found the highest peptide yields were in SIPF reactions performed at temperatures between 60 and 90 °C, which is not unreasonable for a recently condensed primordial ocean or surface temperatures of the early Earth [
9].
- (5)
pH. Optimal pH conditions for the SIPF reaction should be below 3 to prevent Cu-catalyzed peptide hydrolysis and above 2 to prevent proton-catalyzed peptide hydrolysis. Acidic environments on prebiotic Earth, such as acidic hot springs etc., could provide suitable pH conditions.
An intriguing outcome of previous studies of the SIPF reaction is the reported enhanced reactivity of
l-amino acid enantiomers over their
d-amino acid counterparts. For valine, for example, Plankensteiner and co-workers (2005a) reported divaline yields that were as much as 20- to 380-fold for
ll peptides than
dd peptides formed via the SIPF reaction [
29]. Separately, Plankensteiner and co-workers (2005b) measured the chiral selectivity of peptide formation of the SIPF reaction with several individual amino acids [
15]. This work revealed maximum selectivities for
ll over
dd peptides of the following amounts depending on the reaction conditions (peptide, (ratio of
ll/
dd)): divaline (408.33); diproline (1.69); ditryptophan (2.54); dilysine (1.32); and diserine (1.38). However, these reactions were performed as independent enantiopure reactions that only contained all
l- or all
d-enantiomers. This can be problematic because different enantiopure enantiomers of a given amino acid often come from different sources, such as extraction from biological materials for the naturally occurring
l-enantiomers of amino acids versus the chemical synthesis of the
d-enantiomers. This can lead to the presence of different impurities, which could have impacted the reactions discussed above. For example, previous work found that
l-tyrosine was more soluble than
d-tyrosine, which was initially attributed to parity violating energy differences between the enantiomers [
30]. However, when these experiments were repeated by Lahav and co-workers (2006) using amino acids sourced from several different companies, the results were mixed, with some experiments showing
l-tyrosine as more soluble than, equally soluble as, or less soluble than
d-tyrosine [
31]. In addition, assays of the samples revealed that they each contained different species and abundances of contaminants, strongly implying that these impurities likely played a significant role in the observed apparent differences in solubility. Thus, it remains a possibility that the chiral selectivity preferences for
ll enantiomers observed in the SIPF reaction reported to date could have been similarly affected by differences in impurities between the two enantiomers.
Setting aside potential issues with enantiopure reactions, Plankensteiner et al. (2005b) postulated that the stereoselectivity of the SIPF reaction they observed could be the result of differences in geometry between different complexes that amplified the polarity-violating energetic differences between the two enantiomers [
15]. In order to evaluate the plausibility of this hypothesis, we use density functional theory calculations to compare the stability and structure of heterochiral and homochiral complexes of the individual amino acids alanine, valine and proline, in scenarios containing
d +
d,
l +
d, and
l + l enantiomer pairings for both the protonated–neutral and neutral–neutral complexes. These amino acids were chosen to explore the effects of smaller (alanine) and bulkier (valine) side-chains as well as the effects of having a secondary amine (proline) versus the primary amines alanine and valine. In addition, valine was also chosen because it was reported to show the largest stereoselectivity in the SIPF reaction with over a 400-fold increase in reactivity between
l + l and
d +
d SIPF reactions.
5. Conclusions
The present calculations suggest that energy differences between homochiral and heterochiral Cu–amino acid complexes, if present, are small enough to be obscured by rotations of functional groups within the models. However, for the neutral–neutral Cu–amino acid complexes, energy differences between trans and cis isomers were not masked by molecular rotations, with trans isomers generally being energetically preferred. This work suggests that smaller, less complex models should be explored to compare the stability of heterochiral and homochiral SIPF complexes in tandem with laboratory experiments. Comparisons of the structures of the lowest energy, protonated–neutral Cu–amino acid complexes were not consistent with previous work that found more ‘distorted’ complexes in amino acids that reportedly showed increased stereoselectivity in experimental investigations. Consistent trends for both torsion angle and CCM among amino acid types were not observed, indicating that distortion of the Cu–amino acid complex is not driving the reported stereoselectivity of the SIPF reaction. Among the neutral–neutral Cu–amino acid complexes, the trans isomers were found to be more distorted and have smaller O-N distances compared to cis isomers, suggesting that trans isomers may better facilitate polymerization. Based on these results, future experimental work should prioritize the performance of: (1) SIPF experiments using racemic mixtures of amino acids to confirm that the stereoselectivity of the reaction is not due to different impurities in commercial l and d amino acid sources; (2) SIPF experiments at varying pH values to determine if amino acid mixtures with higher proportions of neutral amino acids affect polymerization rates, and (3) SIPF experiments in solutions of varying polarity that would influence the proportions of cis vs. trans isomers.