Proton Transfer from the Inactive Gas-Phase Nicotine Structure to the Bioactive Aqueous-Phase Structure

• Published in: Journal of the American Chemical Society 2010-12, Vol.132 (51), p.18067-18077
• Main Authors: Gaigeot, Marie-Pierre, Cimas, Alvaro, Seydou, Mahamadou, Kim, Ju-Young, Lee, Sungyul, Schermann, Jean-Pierre
• Format: Article
• Language: English
• Subjects: Entropy; Molecular Structure; Nicotine - chemistry; Nitrates - chemistry; Protons; Water - chemistry
• ISSN: 0002-7863; 1520-5126
• DOI: 10.1021/ja103759v

The role of water in the structural change of nicotine from its inactive form in the gas phase to its bioactive form in aqueous solution has been investigated by two complementary theoretical approaches, i.e., geometry optimizations and molecular dynamics. Structures of the lowest-energy nicotineH+−(H2O) n complexes protonated either on the pyridine (inactive form) or pyrrolidine (active form) ring have been calculated, as well as the free-energy barriers for the proton-transfer tautomerization between the two cycles. These structures show chains of 2−4 water molecules bridging the two protonation sites. The room-temperature free-energy barrier to tautomerization along the minimum-energy path from the pyridine to the pyrrolidine cycle drops rapidly when the number of water molecules increases from 0 to 4, but still remains rather high (16 kJ/mol with four water molecules), indicating that the proton transfer is a rather difficult and rare event. We compare results obtained through this explicit water molecule approach to those obtained by means of continuum methods. Car−Parrinello molecular dynamics (CPMD) simulations of the proton-transfer process in bulk with explicit water molecules have been conducted at room temperature. No spontaneous proton transfers have been observed during the dynamics, and biased CPMD simulations have therefore been performed in order to measure the free-energy profile of the proton transfer in the aqueous phase and to reveal the proton-transfer mechanism through water bridges. The MD bias involves pulling the proton from the pyridine ring to the surrounding bulk. Dynamics show that this triggers the tautomerization toward the pyrrolidine ring, proceeding without energy barrier. The proton transfer is extremely fast, and protonation of the pyrrolidine ring was achieved within 0.5 ps. CPMD simulations confirmed the pivotal role played by the water molecules that bridge the two protonation sites of nicotine within the bulk of the surrounding water.