Cited by 122DNA and Nucleic Acid ChemistrySpectroscopy and Quantum Chemical StudiesPhotochemistry and Electron Transfer Studies

The DOI matches its Crossref record; 67 of the 67 references it cites were found in a registry or an index; 24 of the 46 sources its claims cite were checked for support.

An open quantum systems approach to proton tunnelling in DNA

Open access
Slocombe, Louie, Sacchi, Marco, Al-Khalili, Jim·2022·Communications Physics

Abstract

One of the most important topics in molecular biology is the genetic stability of DNA. One threat to this stability is proton transfer along the hydrogen bonds of DNA that could lead to tautomerisation, hence creating point mutations. We present a theoretical analysis of the hydrogen bonds between the Guanine-Cytosine (G-C) nucleotide, which includes an accurate model of the structure of the base pairs, the quantum dynamics of the hydrogen bond proton, and the influence of the decoherent and dissipative cellular environment. We determine that the quantum tunnelling contribution to the proton transfer rate is several orders of magnitude larger than the classical over-the-barrier hopping. Due to the significance of the quantum tunnelling even at biological temperatures, we find that the canonical and tautomeric forms of G-C inter-convert over timescales far shorter than biological ones and hence thermal equilibrium is rapidly reached. Furthermore, we find a large tautomeric occupation probability of 1.73 × 10 −4 , suggesting that such proton transfer may well play a far more important role in DNA mutation than has hitherto been suggested. Our results could have far-reaching consequences for current models of genetic mutations.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper models proton transfer in the Guanine–Cytosine hydrogen bonds using an open quantum systems approach with a Caldeira–Leggett thermal bath. It finds that quantum tunnelling dominates over classical over-the-barrier hopping, giving rapid canonical–tautomeric interconversion and a tautomeric occupation probability of \(1.73 \times 10^{-4}\) at 300 K.

Key findings

  • The quantum tunnelling contribution is much larger than classical hopping: the tunnelling factor \(\kappa\) is \(\sim 10^5\), with quantum forward/reverse rates \(k_f = 7.61 \times 10^5\ \mathrm{s^{-1}}\) and \(k_r = 1.69 \times 10^{13}\ \mathrm{s^{-1}}\), versus classical rates \(k_f = 7.596\ \mathrm{s^{-1}}\) and \(k_r = 1.692 \times 10^8\ \mathrm{s^{-1}}\).
  • The tautomeric occupation probability at 300 K is \(1.73 \times 10^{-4}\), more than four orders of magnitude larger than classical/semiclassical predictions; the classical thermal equilibrium constant from \(\Delta E = 0.435\ \mathrm{eV}\) is \(K_{\mathrm{eq}} = 5.4 \times 10^{-8}\).
  • The canonical and tautomeric forms interconvert on timescales far shorter than biological ones, so thermal equilibrium is rapidly reached; quantum half-lives are \(9.11 \times 10^{-7}\ \mathrm{s}\) for the reactant and \(4.09 \times 10^{-14}\ \mathrm{s}\) for the product.
  • The kinetic isotope effect at biological temperatures is \(\mathrm{KIE} = 30\), and the model uses a forward barrier of \(0.705\ \mathrm{eV}\), reverse barrier of \(0.270\ \mathrm{eV}\), and first eigenstate energy \(E_0 = 0.049\ \mathrm{eV}\).

Methods

  • The proton is treated as a quantum system in an asymmetric double-well potential derived from an ab initio potential energy surface for G–C tautomerisation, represented by a back-to-back double Morse potential.
  • Dynamics use the Caldeira–Leggett quantum Brownian motion model, with the proton coupled to an Ohmic thermal bath of harmonic oscillators; the Wigner–Moyal–Caldeira–Leggett equation is solved numerically with a low-temperature correction replacing \(k_B T\) by \(E_\omega = \frac{\hbar\omega}{2}\coth(\frac{\hbar\omega}{2k_B T})\).
  • Numerical integration uses the method of lines, finite differences, and VCABM5 in Julia, with proton mass \(1836\ \mathrm{a.u.}\), bath temperature \(298.15\ \mathrm{K}\), and friction \(\gamma = 3900\ \mathrm{cm^{-1}}\); rates are obtained from probability flux, and the KIE is computed by doubling the mass.

Limitations

  • The paper states that the master equation is rigorously valid only in the weak-coupling and high-temperature limits; at biological temperatures the high-temperature limit is invalid, so a low-temperature correction is applied.
  • The model drops third-order potential terms, assumes a single effective reaction coordinate with a single proton mass at the transition state, and neglects collective classical rearrangement of the remaining base atoms.
  • The paper notes that experimental tautomerisation rates for G–C double proton transfer are not available, that KIE measurements in biological assays are problematic, and that the KIE calculation neglects isotopic changes to the energy landscape; it also suggests added dimensions could introduce corner-cutting effects.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s42005-022-00881-8

Cites

67 of 67 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Verification

Claims the paper makes

29 Own interpretation · 25 Own result · 15 No basis given · 10 No full text of the sources · 7 Supported · 7 Partly supported · 2 Supported so far · 2 Sources not found · accuracy not yet measured

Show all 61 claims

  • Probably not a substantial claim. 11
  • It may state two separate claims as one. 27
36 of 36
  • The authors state this as an assumption.

    The two protons in G-C undergo asynchronous transfer, with the rate-determining process being the transfer of a single proton across the barrier, and the other proton then moving in the opposite direction to preserve electroneutrality.

    “We assume, as previously observed 6,14 , that the two protons undergo asynchronous transfer, suggesting that the ratedetermining process is the transfer of a single proton across the barrier, with the other proton then moving in the opposite direction to preserve electroneutrality.”
  • The authors state this as an assumption.

    The bath is dominated by the thermal fluctuations of the surrounding water molecules.

    “Assuming that the bath is dominated by the thermal fluctuations of the surrounding water molecules, we can estimate the value of γ.”
  • The paper relies on this without support.

    Experimental measurements of tautomerisation rates for G-C double proton transfer reaction are still not available.

    “Experimental measurements of tautomerisation rates for G-C double proton transfer reaction are still not available.”
  • The paper relies on this without support.

    DNA repair mechanisms are highly efficient.

    “the well-known presence of highly efficient DNA repair mechanisms”
  • The authors state this as an assumption.

    Heisenberg's position-momentum uncertainty principle must be upheld.

    “We note additionally that Heisenberg's position-momentum uncertainty principle must also be upheld”
  • The authors state this as an assumption.

    Changes to the transfer energy landscape due to isotopic substitutions (such as zero-point energy and free energy contributions) are neglected; only mass effects on the reaction rate are considered.

    “We neglect any changes to the transfer energy landscape due to isotopic substitutions, such as zero-point energy and free energy contributions, and consider only the effect on the reaction rate due to modifications of the mass terms in Eq. ( 1).”
  • The authors state this as an assumption.

    The third-order potential terms in Eq. (1) are dropped to simplify numerical calculations.

    “To simplify the numerical calculations, we drop the third-order potential terms in Eq. ( 1).”
  • The authors state this as an assumption.

    The mass of the proton is, in a very good approximation, the effective mass of the tunnelling particle in the asynchronous concerted process.

    “it is acceptable to assume that the mass of the proton mass is in, a very good approximation, the effective mass of the tunnelling particle in the asynchronous concerted process 5 .”
  • The authors state this as an assumption.

    The collective classical rearrangement of the rest of the atoms of both bases can be neglected.

    “Furthermore, the collective classical rearrangement of the rest of the atoms of both bases can be neglected.”
9 of 9

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

1 of 1

Citations

Cites

Show all 63 references

Research works printed but not linked: 1

  1. Wales, D., Saykally, R., Zewail, A. & King, D. Energy Landscapes: Applications to Clusters, Biomolecules and Glasses (Cambridge University Press, 2003). https://books.google.co.uk/books?id=YQrB6s3LALEC.

Cited by

Sorted by citations

Metadata

Authors
Slocombe, L., Sacchi, M., & Al-Khalili, J.
Affiliations
University of Surrey (GB)
Corresponding author
  • Slocombe, Louie
  • Sacchi, Marco
  • Al-Khalili, Jim
Venue
Communications Physics
Published
2022-05-05
Volume
5
Issue
1
ISSN
2399-3650 (ISSN-L)
OpenAlex
W3204256293
Cited by
122
References
63
Citations by year
2026: 30 · 2025: 31 · 2024: 24 · 2023: 30 · 2022: 7
FWCI
6.35
Type
journal article
Language
en
Open access
Locations
  • Communications Physics (publishedVersion, cc-by)
  • arXiv (Cornell University) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • arXiv (Cornell University) (cc-by)
Funders
Leverhulme Trust (DS- 2017-079)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › DNA and Nucleic Acid Chemistry
  • Physical Sciences › Physics and Astronomy › Atomic and Molecular Physics, and Optics › Spectroscopy and Quantum Chemical Studies
  • Physical Sciences › Chemistry › Physical and Theoretical Chemistry › Photochemistry and Electron Transfer Studies
Keywords
proton tunneling, DNA, hydrogen bonding, guanine-cytosine base pair, DNA mutations, open quantum systems, decoherence, energy dissipation, quantum tunneling, point mutations, genetic stability

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-19, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-19 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review