Quantum Mechanics/Molecular Mechanics Simulations for Chiral-Selective Aminoacylation: Unraveling the Nature of Life

Biological phenomena are chemical reactions, which are inherently non-stopping or “flowing” in nature. Molecular dynamics (MD) is used to analyze the dynamics and energetics of interacting atoms, but it cannot handle chemical reactions involving bond formation and breaking. Quantum mechanics/molecular mechanics (QM/MM) umbrella sampling MD simulations gives us a significant clue about transition states of chemical reactions and their energy levels, which are the pivotal points in understanding the nature of life. To demonstrate the importance of this method, we present here the results of our application of it to the elucidation of the mechanism of chiral-selective aminoacylation of an RNA minihelix considered to be a primitive form of tRNA. The QM/MM MD simulation, for the first time, elucidated the “flowing” atomistic mechanisms of the reaction and indicated that the L-Ala moiety stabilizes the transition state more than D-Ala, resulting in L-Ala preference in the aminoacylation reaction in the RNA. The QM/MM method not only provides important clues to the elucidation of the origin of homochirality of biological systems, but also is expected to become an important tool that will play a critical role in the analysis of biomolecular reactions, combined with the development of artificial intelligence.