Molecular Mechanisms of Phosphoester Bond Formation in Water Using Tight-Binding Ab Initio Molecular Dynamics

Phosphate groups are ubiquitous in biomolecules and are usually incorporated through phosphoester bonds between alcohol groups and orthophosphate. The formation of this bond is exceptionally difficult, with associated barriers of 30-45 kcal/mol in the absence of catalysts. In abiotic conditions, polymerizing nucleic acids without enzymes remains very challenging and is still a partly unsolved problem that severely questions the RNA World hypothesis for the origins of life. Offering a solution to this problem would involve a detailed knowledge of the reaction energetics and mechanisms, yet these remain not fully understood at a molecular level, especially because of the very slow reaction rates that represent a significant challenge for the experiments. The number of involved reaction coordinates and the possible role of the solvent in assisting the reaction are challenging for computational studies. Here, we use extensive ab initio molecular dynamics simulations using semiempirical tight-binding methods and enhanced sampling to address these issues. We first show that the choice of the tight-binding method is greatly limited by the instability of the water liquid phase for most DFTB generations and parameter sets that are widely available. We then focus on a model reaction involving methanol and orthophosphate, for which the two protonation states (mono- and dianionic) that are dominant around neutral pH are considered. We compare different proton coordinates that enable (or not) the participation of solvent water molecules. Our simulations suggest that in all cases, a dissociative associative mechanism, with an intermediate metaphosphate, is favored. The main difference between the two phosphate species is that reaction with the monoanion is assisted by the substrate, while that with the dianion involves solvent water molecules. Our results are in agreement with early experimental measurements, but the reaction barriers are underestimated in our framework. We believe that our approach provides an interesting perspective on how to sample the reaction phase space efficiently, but it calls for future studies using more accurate descriptions of chemical reactivity.