Computational Modeling Reveals Pathway for 3-Methylindole Formation

The work addresses a significant gap in food chemistry, where the abiotic formation of indole derivatives under thermal processing conditions has remained largely speculative compared to their well-documented biological metabolism in gut microbiota.

The Maillard reaction is a non-enzymatic browning process initiated by the condensation of reducing sugars and amino compounds under heat. While responsible for desirable sensory attributes in baked goods and roasted meats, this complex cascade also degrades essential amino acids and can generate detrimental compounds. Tryptophan, an essential aromatic amino acid abundant in legumes, participates vigorously in this chemistry due to its unique indole ring structure. When heated with reducing sugars like glucose, tryptophan yields products such as indole and 3-methylindole, which exhibit a dual sensory impact: pleasant floral notes at low concentrations and undesirable animal-like odors at higher levels.

To identify the key odor-active volatiles, the research team conducted gas chromatography-olfactometry-mass spectrometry (GC-O-MS) analysis on tryptophan-glucose mixtures heated at 120, 130, and 140 degrees Celsius. A total of 28 volatile compounds were identified, including aldehydes, ketones, and indoles. Among these, 3-methylindole stood out as the primary compound for mechanistic investigation, showing the highest concentration range of 16.87 to 21.2 micrograms per kilogram and a maximum odor intensity rating of 5, described as strong jasmine fragrance. Other significant contributors included 1-nonanal and indole, which contributed fruity and floral notes, respectively.

The core of the study employed ReaxFF molecular dynamics (MD) simulations to model the reaction at the atomic scale. The simulation box contained 29 tryptophan molecules and 50 glucose molecules, mimicking a molar ratio of approximately 29:50. To accelerate reactions within the limited simulation timescale, the system was heated to elevated temperatures of 2100, 2400, and 2700 Kelvin under pressures of 100 and 150 megapascals. These conditions allowed the researchers to track bond-breaking and bond-forming events, revealing the dynamic evolution of intermediates such as water, carbon dioxide, and critical alpha-diketones like glyoxal and methylglyoxal.

The simulations indicated that temperature significantly influenced the reaction profile, with higher temperatures promoting the formation of products and shortening the onset time for key intermediates. Methylglyoxal, a crucial intermediate in the proposed pathway, appeared at 218 picoseconds under 2100 Kelvin and 100 megapascals, but its appearance time shortened to 14 picoseconds at 2700 Kelvin. Conversely, increasing pressure reduced system density and overall potential energy, making the reaction less favorable. This temperature-dependent behavior validated the role of methylglyoxal as a key intermediate in the thermal degradation of the tryptophan-glucose system.

Building on the dynamic pathways identified by ReaxFF MD, the researchers used density functional theory (DFT) calculations to refine the mechanism and characterize transition states. Using the B3LYP/6-311G+(d,p) basis set, they optimized the structures of reactants, intermediates, and products, ensuring the absence of imaginary frequencies. Sixteen transition states were successfully identified and validated through intrinsic reaction coordinate analysis, which traced the minimum energy path connecting reactants to products. This high-precision optimization provided a reliable geometric configuration for calculating the energy barriers governing the dominant reaction channels.

The calculated energy profile revealed a complex, multi-stage process with a net endothermic feature. The pathway begins with the nucleophilic addition of the tryptophan amino group to the glucose carbonyl group, crossing an energy barrier of 45.14 kilocalories per mole. This is followed by an Amadori rearrangement with a barrier of 223.29 kilocalories per mole. The highest single barrier identified was 432.57 kilocalories per mole, associated with a hydrogen transfer step involving intermediate IM6. The final step involves the decarbonylation of indole-3-acetaldehyde to produce 3-methylindole, crossing a barrier of 241.05 kilocalories per mole. These energy values highlight the kinetic bottlenecks throughout the pathway, with the overall transformation requiring significant thermal input.

To experimentally validate the computational predictions, the team performed pyrolysis-gas chromatography-mass spectrometry (Py-GC/MS) analysis. Although the pyrolysis temperature of 550 degrees Celsius differs substantially from the actual Maillard reaction conditions, the technique functions as a structural fragment release tool. It detected key intermediates structurally consistent with those predicted by the simulations, including methylglyoxal and indole-3-acetaldehyde. The concentration trends of these compounds with temperature showed an initial increase followed by a decrease at the highest temperature, consistent with their behavior as reactive intermediates that are further degraded at elevated temperatures.

Differential scanning calorimetry (DSC) analysis provided macroscopic thermal evidence supporting the computational findings. The DSC curves for the tryptophan-glucose mixture exhibited a net endothermic feature, with major thermal events occurring between 83 and 156 degrees Celsius. The peak around 83 degrees Celsius was attributed to the glass transition of glucose, while the peak near 156 degrees Celsius likely corresponded to early-stage reactions between the reactants. The endothermic direction observed in the DSC analysis matched the DFT-calculated energy profile, which showed a net increase in energy from reactants to products, reinforcing the credibility of the proposed reaction pathway.

The proposed mechanism connects glucose-derived carbonyl chemistry with the conversion of the tryptophan side chain. The pathway involves dehydration, enolization, decarboxylation, hydrogen transfer, and Strecker degradation reactions. Specifically, the Amadori product undergoes 2,3-enolization to yield 1-deoxythreose, which then yields methylmalonaldehyde. Tryptophan subsequently undergoes Strecker degradation with methylglyoxal to form indole-3-acetaldehyde, which is finally decarbonylated to 3-methylindole. This integrated approach bridges multi-scale computational modeling with experimental validation, providing a coherent microscopic image and thermodynamic basis for understanding flavor formation.

The study establishes a transferable research framework for mechanistically studying volatile compound generation in complex food systems. By combining ReaxFF MD for dynamic pathway mapping, DFT for energy barrier quantification, and experimental techniques like GC-O-MS and Py-GC/MS for validation, the researchers clarified the heat-induced molecular source of 3-methylindole. This understanding is crucial for flavor engineering in foods such as fermented sauces, dairy products, and alcoholic beverages, where 3-methylindole is a key aroma compound. The findings also highlight the importance of controlling thermal processing parameters to manage the concentration-dependent sensory behavior of indole compounds.

Limitations of the study include the use of elevated temperatures in ReaxFF MD simulations to accelerate reactions, which may not fully replicate the kinetic details of real-world food processing at lower temperatures. Additionally, the Py-GC/MS analysis was conducted at 550 degrees Celsius, which serves as a structural probe rather than a direct simulation of Maillard conditions. Future research could further quantify the relative contributions of competitive pathways under different processing conditions and employ isotope labeling or time-resolved spectroscopy to verify key intermediates in real-time. The current work provides a foundational mechanistic map that can guide subsequent experimental and industrial applications in food science.

Subject of Research: Mechanistic analysis of 3-methylindole formation in the tryptophan-glucose Maillard reaction

Article Title: Unraveling the formation mechanism of 3-methylindole from tryptophan in the glucose-mediated Maillard reaction: a combined ReaxFF MD, DFT, and experimental study

Article References: Yang, Y.-Y., Wei, C.-K., Pan, X., Zheng, A.-R., Wang, W.-L., Yang, C.-M., & Liu, Y. (2026). Unraveling the formation mechanism of 3-methylindole from tryptophan in the glucose-mediated Maillard reaction: a combined ReaxFF MD, DFT, and experimental study. Food Chemistry: X, 40, Article 104719. https://doi.org/10.1016/j.fochx.2026.104719

Image Credits: AI Generated

DOI: 10.1016/j.fochx.2026.104719

Keywords: Maillard reaction, 3-methylindole, Tryptophan, Glucose, ReaxFF MD, DFT, Food chemistry, Flavor compounds, Molecular dynamics, Py-GC/MS, Unraveling, formation