AUTHOR=Brovarets' Ol'ha O. , Oliynyk Timothy A. , Hovorun Dmytro M. 

TITLE=Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G*·C* DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey

JOURNAL=Frontiers in Chemistry

VOLUME=7

YEAR=2019

URL=https://www.frontiersin.org/journals/chemistry/articles/10.3389/fchem.2019.00597

DOI=10.3389/fchem.2019.00597

ISSN=2296-2646

ABSTRACT=<p>For the first time, in this study with the use of QM/QTAIM methods we have exhaustively investigated the tautomerization of the biologically-important conformers of the G<sup>*</sup>·C<sup>*</sup> DNA base pair—reverse Löwdin G<sup>*</sup>·C<sup>*</sup>(rWC), Hoogsteen G<sup>*</sup>′·C<sup>*</sup>(H), and reverse Hoogsteen G<sup>*</sup>′·C<sup>*</sup>(rH) DNA base pairs—<italic>via</italic> the single (SPT) or double (DPT) proton transfer along the neighboring intermolecular H-bonds. These tautomeric reactions finally lead to the formation of the novel G·<inline-formula><mml:math id="M1" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msubsup><mml:mtext>C</mml:mtext><mml:mrow><mml:mtext>O</mml:mtext><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>(rWC), <inline-formula><mml:math id="M2" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msubsup><mml:mtext>G</mml:mtext><mml:mrow><mml:mtext>N</mml:mtext><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>·</mml:mo></mml:math></inline-formula>C(rWC), G*′<sub>N2</sub>·C(rWC), <inline-formula><mml:math id="M4" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msubsup><mml:mtext>G</mml:mtext><mml:mrow><mml:mtext>N</mml:mtext><mml:mn>7</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>·</mml:mo></mml:math></inline-formula>C(H), and G*′<sub>N7</sub>·C(rH) DNA base mispairs. Gibbs free energies of activation for these reactions are within the range 3.64–31.65 kcal·mol<sup>−1</sup> in vacuum under normal conditions. All TSs are planar structures (C<sub>s</sub> symmetry) with a single exception—the essentially non-planar transition state TS<sub>G*·C*(rWC)↔G<sup>+</sup>·C<sup>−</sup>(rWC)</sub> (C<sub>1</sub> symmetry). Analysis of the kinetic parameters of the considered tautomerization reactions indicates that in reality only the reverse Hoogsteen G<sup>*</sup>′·C<sup>*</sup>(rH) base pair undergoes tautomerization. However, the population of its tautomerised state G*′<sub>N7</sub>·C(rH) amounts to an insignificant value−2.3·10<sup>−17</sup>. So, the G<sup>*</sup>·C<sup>*</sup>(rWC), G<sup>*</sup>′·C<sup>*</sup>(H), and G<sup>*</sup>′·C<sup>*</sup>(rH) base pairs possess a permanent tautomeric status, which does not depend on proton mobility along the neighboring H-bonds. The investigated tautomerization processes were analyzed in details by applying the author's unique methodology—sweeps of the main physical and chemical parameters along the intrinsic reaction coordinate (IRC). In general, the obtained data demonstrate the tautomeric mobility and diversity of the G<sup>*</sup>·C<sup>*</sup> DNA base pair.</p>