💻 bioinformatics

Glycine molecule radical: Predicted properties and dipeptide formation

Using density functional theory calculations, this study proposes a novel, catalyst-free radical chemistry pathway for the formation of glycine dipeptides under prebiotic conditions, characterized by low electronic energy barriers and supported by detailed analyses of intermediate radical properties.

Original authors: Synak, J., Blazewicz, J.

Published 2026-07-25

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Synak, J., Blazewicz, J.

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Imagine you are a detective trying to solve the ultimate cold case: How did life begin? Scientists have a strong hunch that before complex cells existed, there was a "RNA World" where simple genetic molecules did all the heavy lifting. But there's a missing piece of the puzzle: how did those molecules eventually start building proteins? Proteins are the workhorses of life, made of long chains of tiny building blocks called amino acids, held together by "peptide bonds." The problem is that in the messy, chaotic environment of early Earth, snapping these blocks together usually requires a sophisticated factory machine (an enzyme) that didn't exist yet. So, the big question remains: How did the first protein chains form without any tools?

To understand the story in this paper, you need to know a few basic concepts. First, think of a radical as a molecule that has lost a piece of itself, leaving it with an "open hand" (an unpaired electron) that makes it desperate to grab onto something else. These are like hyperactive kids at a party who can't sit still. Second, energy barriers are like hills you have to push a ball over to get it to the other side. If the hill is too high, the ball never makes it; if it's low enough, the ball rolls right over. Finally, simulations are like super-advanced video games where scientists use math to predict how atoms will behave, rather than mixing chemicals in a test tube. This paper uses these simulations to see if a wild, tool-free idea could actually work.

The Radical Recipe for Life's First Chains

This paper, written by Jaroslaw Synak and Jacek Blazewicz, dives into the world of computer simulations to test a daring idea: Could the first peptide bonds have formed using nothing but "radical chemistry"? Instead of waiting for a complex biological machine to appear, the authors suggest that simple, hyperactive glycine molecules (the simplest amino acid) might have bumped into each other, broke apart in a specific way, and reassembled into a chain all on their own.

The researchers focused on a specific character in this story: the glycine acyl radical. Imagine a glycine molecule as a small Lego brick. In this scenario, the brick gets hit by a stray hydrogen atom, loses a piece of its "arm" (the amino group), and becomes a radical. This new, unstable version is desperate to connect. The authors' simulations show that this radical can approach another normal glycine molecule. Because of the way their electrical charges are arranged—like opposite poles of a magnet—they are naturally drawn together.

When they meet, the connection happens surprisingly fast. The radical grabs the other molecule, a tiny hydrogen atom pops off and escapes, and a new bond is formed, creating a two-block chain called diglycine. The best part? The "hill" (energy barrier) they had to climb to make this happen was relatively small, at 19.34 kcal/mol. In the world of chemistry, that's a low enough hill that this reaction could have happened naturally in the wild, without any fancy catalysts to push it over.

The Cast of Characters: Stable vs. Unstable

To make sure this recipe was solid, the authors had to investigate the ingredients. They looked at different ways a glycine molecule could break apart to become a radical.

The "Simple" Radical: If you just break a carbon-hydrogen bond in glycine, you get a "simple" radical. The paper found this one is actually quite stable and interesting. It has a special "resonance" structure, which is like a molecule that can't decide which shape to be, so it exists as a blend of three different shapes at once. This makes the molecule flat and stiff, almost like a rigid ruler. It can even flip between two shapes (called cis and trans), but it takes a bit of energy (18.70 kcal/mol) to force that flip.
The "Wrong" Radicals: The authors also checked what happens if you break the bond on the oxygen side or the nitrogen side. Unfortunately, these versions are troublemakers. If you break the oxygen bond, the molecule instantly falls apart into carbon dioxide and a different radical, like a house of cards collapsing. If you break a nitrogen bond, it quickly rearranges itself and falls apart too. The paper explicitly rules these out as candidates for building long chains; they are dead ends.

The Two-Step Dance to the Acyl Radical

So, how do we get the star of the show, the glycine acyl radical? The authors suggest a two-step dance.

Step 1: A glycine molecule meets a hydrogen atom and loses its amino group. This is a bit of a bumpy ride, requiring an energy push of about 12.92 kcal/mol, but it releases a lot of energy (31.54 kcal/mol) once it happens.
Step 2: The ammonia that was released in step 1 can reattach, and a free hydrogen atom attacks the remaining part, kicking off a water molecule. This second step has a high energy hill (45.21 kcal/mol), but if you look at the whole journey from the very start, the total energy cost is manageable.

The authors note that while this two-step process seems complicated, it's actually more efficient than trying to smash a hydrogen atom directly onto the oxygen part of the molecule, which would require a massive 29.29 kcal/mol energy push.

What This Means (and What It Doesn't)

The main takeaway from this study is a "proof of concept" simulation. It suggests that a pathway exists where glycine radicals could link up to form a dipeptide (a two-amino-acid chain) with an energy barrier of 19.34 kcal/mol. This is low enough to be plausible in a prebiotic world.

However, it is important to keep the excitement in check. The authors are very clear that this is a simulation, not a physical experiment where they mixed chemicals in a lab and watched them react. They used a specific mathematical method (B3LYP) which is great for small molecules but has known limits for larger, more complex systems. The energy numbers they provide are approximations based on current computer models, not absolute, unchangeable facts.

The paper doesn't claim to have solved the mystery of life's origin. Instead, it offers a new, tool-free recipe that looks promising. It suggests that if radicals were present in the early Earth's environment, they might have been the spark that started the chain reaction leading to proteins. The authors admit that the "simple" glycine radical is interesting, but the acyl radical is the much better candidate for actually building the chain.

In short, this paper paints a picture of early Earth not as a place waiting for complex machines, but as a chaotic playground where hyperactive, broken pieces of molecules could accidentally bump into each other, stick together, and start the long journey toward life. It's a suggestion, a simulation, and a fascinating new angle on an ancient mystery.

Technical Summary: Glycine Molecule Radical: Predicted Properties and Dipeptide Formation

Problem Statement
The transition from the RNA World to protein-based life remains a significant unsolved challenge in origin-of-life research. A primary obstacle is identifying a plausible prebiotic pathway for peptide bond formation that does not rely on the sophisticated biological machinery found in modern systems. While empirical pathways involving thioesters, surface chemisorption, or specific catalysts have been demonstrated, the authors focus on the formation of diglycine (the simplest dipeptide) as a minimal test case. The study aims to predict a prebiotically plausible reaction pathway utilizing radical chemistry without external catalysts, addressing the gap in understanding how the first peptide bonds could form in a relatively simple environment.

Methodology
The authors employed Density Functional Theory (DFT) calculations using the Gaussian 16 software suite. The specific computational setup utilized the B3LYP functional with unrestricted electrons and the 6-31G(d) basis set. This combination was selected to provide sufficient precision for small molecular systems while maintaining reasonable computational efficiency to explore various scenarios.

The investigation followed a three-stage theoretical and computational approach:

Theoretical Screening: Reagents were examined using resonance hybrid models and organic chemistry principles to identify reasonable interactions.
Energetic Feasibility: Total energies of reactants and products were calculated to determine if the energy changes were reasonable.
Transition State Analysis: The most complex stage involved guessing transition states and utilizing the Gaussian "relaxed scan" function to optimize inter-atomic distances and angles. Intrinsic Reaction Coordinate (IRC) calculations were performed to map the energy changes throughout the reaction processes.

Key Contributions and Results

Peptide Bond Formation Mechanism:
The study predicts a novel pathway for diglycine formation involving the reaction between a glycine acyl radical and a neutral glycine molecule. The mechanism is driven by electrostatic attraction between the partially positive acyl carbon and the partially negative nitrogen of the amino group.
- Energetics: The reaction possesses an electronic energy barrier of 19.34 kcal/mol, which the authors consider relatively low and accessible without sophisticated catalysts.
- Thermodynamics: The process is endergonic, with products having 8.13 kcal/mol more energy than reactants. However, the mechanism involves the immediate rejection of a hydrogen atom, which facilitates the reaction progress.
Glycine Radical Properties:
The authors thoroughly analyzed the "simple" glycine radical (formed by breaking a C-H bond).
- Structure: Using resonance hybrid models and Mulliken charge analysis, the radical is shown to possess a flat geometry with $sp^2$ hybridization character for the carbon and nitrogen atoms. This results in partial double-bond character for both C-C and C-N bonds.
- Isomerism: The restricted rotation around the C-C bond leads to cis-trans (E/Z) isomerism. The energy barrier to switch from the Z-isomer to the E-isomer is 18.70 kcal/mol.
- Synthesis: The formation of the (E)-isomer via C-H bond cleavage by atomic hydrogen is highly favorable, releasing 26.49 kcal/mol with a low activation barrier of 3.19 kcal/mol.
Alternative Radical Pathways:
The study investigated radicals formed by removing hydrogen from the OH or NH groups.
- OH-radical: Unstable; rapidly decays into a methylamine radical and carbon dioxide.
- NH-radical: Less unstable but undergoes a decay process involving hydrogen bond formation and subsequent bond breaking, resulting in an exergonic reaction (net energy balance of 14.89 kcal/mol).
Glycine Acyl Radical Synthesis:
The glycine acyl radical, identified as the crucial reactant for peptide bond formation, is predicted to form via a two-step process:
1. Reaction of glycine with an atomic hydrogen atom, losing the amino group (activation energy ~12.92 kcal/mol).
2. Reattachment of ammonia and reaction of the released hydrogen with the OH group to form water.
  While the second step has a high barrier (~45.21 kcal/mol), the authors note that relative to the initial reactants, the effective barrier is lower (13.67 kcal/mol). The intermediate acetic acid radical is stabilized by resonance.

Significance and Claims
The paper claims to present a series of reactions that could lead to dipeptide formation under prebiotic conditions without the need for complex catalysts. The primary significance lies in the prediction of a pathway with an electronic energy barrier (<20 kcal/mol) that is accessible in early Earth environments.

The authors explicitly state that their results serve as a "possible starting point for further peptide chain extension" and a "good starting point for anyone researching these radical pathways." They maintain a modest stance regarding their claims, acknowledging that:

The energy values are approximations due to the limitations of current quantum chemistry models (specifically B3LYP) for larger systems.
The presence of the necessary radicals in prebiotic environments must be corroborated by other means.
The acyl radical appears to be a more viable candidate for peptide formation than the simple glycine radical, though the latter was analyzed for its interesting properties and potential role in other pathways.

Ultimately, the study suggests that radical chemistry alone, without external catalysis, offers a plausible mechanism for the emergence of the first peptide bonds, representing a potential step toward protein-based biology.

The Radical Recipe for Life's First Chains

The Cast of Characters: Stable vs. Unstable

The Two-Step Dance to the Acyl Radical

What This Means (and What It Doesn't)

Technical Summary: Glycine Molecule Radical: Predicted Properties and Dipeptide Formation

More like this