Enol tautomers of Watson-Crick base pair models are metastable because of nuclear quantum effects

A review on point mutations via proton transfer in DNA base pairs in the absence and presence of electric fields

This is a comprehensive review that focuses on spontaneous mutations that may occur during DNA replication, the fundamental process responsible for transferring genetic information. In 1963, Löwdin postulated that these mutations are primarily a result of proton transfer reactions within the hydrogen-bonded DNA base pairs. The single and double proton transfer reactions within the base pairs result in zwitterions and rare tautomers, respectively. For persistent mutations, these products must be generated at high rates and should be thermodynamically stable. This review covers the proton transfer reactions studied experimentally and computationally. The review also examines the influence of externally applied electric fields on the thermodynamics and kinetics of proton transfer reactions within DNA base pairs, and their biological implications.

Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine

We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.

Multiscale simulations reveal the role of PcrA helicase in protecting against spontaneous point mutations in DNA

Proton transfer across hydrogen bonds in DNA can produce non-canonical nucleobase dimers and is a possible source of single-point mutations when these forms mismatch under replication. Previous computational studies have revealed this process to be energetically feasible for the guanine-cytosine (GC) base pair, but the tautomeric product (G[Formula: see text]C[Formula: see text]) is short-lived. In this work we reveal, for the first time, the direct effect of the replisome enzymes on proton transfer, rectifying the shortcomings of existing models. Multi-scale quantum mechanical/molecular dynamics (QM/MM) simulations reveal the effect of the bacterial PcrA Helicase on the double proton transfer in the GC base pair. It is shown that the local protein environment drastically increases the activation and reaction energies for the double proton transfer, modifying the tautomeric equilibrium. We propose a regime in which the proton transfer is dominated by tunnelling, taking place instantaneously and without atomic rearrangement of the local environment. In this paradigm, we can reconcile the metastable nature of the tautomer and show that ensemble averaging methods obscure detail in the reaction profile. Our results highlight the importance of explicit environmental models and suggest that asparagine N624 serves a secondary function of reducing spontaneous mutations in PcrA Helicase.

Nuclear Quantum Effects Made Accessible: Local Density Fitting in Multicomponent Methods

The simulation of nuclear quantum effects (NQEs) is crucial for an accurate description of systems and processes involving light nuclei, such as hydrogen atoms. Within the last years, the importance of those effects has been highlighted for a vast range of systems with tremendous implications in chemistry, biology, physics, and materials sciences. However, while electronic structure theory methods have become routine tools for quantum chemical investigations, there is still a lack of approaches to address NQEs that are computationally accessible and straightforward to use. To address this, we present the first combination of the nuclear-electronic orbital Hartree-Fock approach with both local and density fitting approximations (LDF-NEO-HF). This results in a low-order scaling approach that enables the inclusion of NQEs for large systems within a fraction of a day and for small to medium size systems in minutes. Moreover, we demonstrate the qualitative accuracy and robustness of our approach to retrieve NQEs for three real-use cases motivated by chemical, biological, and materials science applications.

Measuring the Enthalpy of an Individual Hydrogen Bond in a DNA Duplex with Nucleobase Isotope Editing and Variable-Temperature Infrared Spectroscopy

The level of interest in probing the strength of noncovalent interactions in DNA duplexes is high, as these weak forces dictate the range of suprastructures the double helix adopts under different conditions, in turn directly impacting the biological functions and industrial applications of duplexes that require making and breaking them to access the genetic code. However, few experimental tools can measure these weak forces embedded within large biological suprastructures in the native solution environment. Here, we develop experimental methods for detecting the presence of a single noncovalent interaction [a hydrogen bond (H-bond)] within a large DNA duplex in solution and measure its formation enthalpy (ΔH_f). We report that introduction of a H-bond into the TC2═O group from the noncanonical nucleobase 2-aminopurine produces an expected decrease ∼10 ± 0.76 cm^-1 (from ∼1720 cm^-1 in Watson-Crick to ∼1710 cm^-1 in 2-aminopurine), which correlates with an enthalpy of ∼0.93 ± 0.066 kcal/mol for this interaction.

Tautomeric Equilibrium in 1-Benzamidoisoquinoline Derivatives

In this study, the tautomeric equilibrium of a sequence of 1-benzamidoisoquinoline derivatives was investigated with the tools of NMR spectroscopy and computational chemistry. The equilibrium between different tautomers in these systems could be controlled via the substitution effect, and the relative content of the amide form varied from 74% for the strong electron-donating NMe2 substituent to 38% for the strong electron-accepting NO2 group in the phenyl ring. In contrast to the previously investigated 2-phenacylquinoline derivatives, the most stable and thus most abundant tautomer in the 1-benzamidoisoquinoline series except the two most electron-accepting substituents was an amide. The intramolecular hydrogen bond present in the enol tautomer competed with the intermolecular hydrogen bonds created with the solvent molecules and thus was not a sufficient factor to favor this tautomer in the mixture. Although routinely computational studies of tautomeric equilibrium are performed within the continuum solvent models, it is proven here that the inclusion of the explicit solvent is mandatory in order to reproduce the experimental tendencies observed for this type of system, facilitating strong intermolecular hydrogen bonds.

Proton transfer during DNA strand separation as a source of mutagenic guanine-cytosine tautomers

Proton transfer between the DNA bases can lead to mutagenic Guanine-Cytosine tautomers. Over the past several decades, a heated debate has emerged over the biological impact of tautomeric forms. Here, we determine that the energy required for generating tautomers radically changes during the separation of double-stranded DNA. Density Functional Theory calculations indicate that the double proton transfer in Guanine-Cytosine follows a sequential, step-like mechanism where the reaction barrier increases quasi-linearly with strand separation. These results point to increased stability of the tautomer when the DNA strands unzip as they enter the helicase, effectively trapping the tautomer population. In addition, molecular dynamics simulations indicate that the relevant strand separation time is two orders of magnitude quicker than previously thought. Our results demonstrate that the unwinding of DNA by the helicase could simultaneously slow the formation but significantly enhance the stability of tautomeric base pairs and provide a feasible pathway for spontaneous DNA mutations.

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Molecular mechanics (MM) is a powerful tool to study the properties of molecular systems in the fields of biology and materials science. With the development of ab initio force field and the application of ab initio potential energy surface, the nuclear quantum effect (NQE) is becoming increasingly important for the robustness of the simulation. However, the state-of-the-art path-integral molecular dynamics simulation, which incorporates NQE in MM, is still too expensive to conduct for most biological and material systems. In this work, we analyze the locality of NQE, using both analytical and numerical approaches, and conclude that NQE is an extremely localized phenomenon in nonreactive molecular systems. Therefore, we can use localized machine learning (ML) models to predict quantum force corrections both accurately and efficiently. Using liquid water as example, we show that the ML facilitated centroid MD can reproduce the NQEs in both the thermodynamical and the dynamical properties, with a minimal increase in computational time compared to classical molecular dynamics. This simple approach thus largely decreases the computational cost of quantum simulations, making it really accessible to the studies of large-scale molecular systems.

Quantum mechanical effects in acid–base chemistry

Acid–base chemistry has immense importance for explaining and predicting the chemical products formed by an acid and a base when mixed together. However, the traditional chemistry theories used to describe acid–base reactions do not take into account the effect arising from the quantum mechanical nature of the acidic hydrogen shuttling potential and its dependence on the acid base distance. Here, infrared and NMR spectroscopies, in combination with first principles simulations, are performed to demonstrate that quantum mechanical effects, including electronic and nuclear quantum effects, play an essential role in defining the acid–base chemistry when 1-methylimidazole and acetic acid are mixed together. In particular, it is observed that the acid and the base interact to form a complex containing a strong hydrogen bond, in which the acidic hydrogen atom is neither close to the acid nor to the base, but delocalized between them. In addition, the delocalization of the acidic hydrogen atom in the complex leads to characteristic IR and NMR signatures. The presence of a hydrogen delocalized state in this simple system challenges the conventional knowledge of acid–base chemistry and opens up new avenues for designing materials in which specific properties produced by the hydrogen delocalized state can be harvested.

Acid-based theories do not consider the quantum mechanical nature of the acidic hydrogen shuttling potential. Here, it is demonstrated that this particularity is needed to explain the formation acid-base complex with a delocalized acidic hydrogen.

Nuclear Quantum Effects in Liquid Water at Near Classical Computational Cost Using the Adaptive Quantum Thermal Bath

We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQEs) in molecular simulations: the adaptive quantum thermal bath (adQTB). In this method, zero-point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to that of classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.

Investigating Tunneling-Controlled Chemical Reactions through Ab Initio Ring Polymer Molecular Dynamics

We use the ab initio ring polymer molecular dynamics (RPMD) approach to investigate tunneling-controlled reactions in methylhydroxycarbene. Nuclear tunneling effects enable molecules to overcome the barriers which cannot be overcome classically. Under low-temperature conditions, intrinsic quantum tunneling effects can facilitate the chemical reaction in a pathway that is favored neither thermodynamically nor kinetically. This behavior is referred to as the tunneling-controlled chemical reaction and is regarded as the third paradigm of chemical reaction controls. In this work, we use the ab initio RPMD approach to incorporate the tunneling effects in our quantum dynamics simulations and investigate the reaction kinetics of two competitive reaction pathways at various temperatures. The reaction rate constants obtained here agree extremely well with the experimentally measured rates. We demonstrate the feasibility of using ab initio RPMD rate calculations in a realistic molecular system and provide an interesting and important example for future investigations of reaction mechanisms dominated by quantum tunneling effects.

Capturing the Catalytic Proton of Dihydrofolate Reductase: Implications for General Acid-Base Catalysis

Acid–base catalysis, which involves one or more proton transfer reactions, is a chemical mechanism commonly employed by many enzymes. The molecular basis for catalysis is often derived from structures determined at the optimal pH for enzyme activity. However, direct observation of protons from experimental structures is quite difficult; thus, a complete mechanistic description for most enzymes remains lacking. Dihydrofolate reductase (DHFR) exemplifies general acid–base catalysis, requiring hydride transfer and protonation of its substrate, DHF, to form the product, tetrahydrofolate (THF). Previous X-ray and neutron crystal structures coupled with theoretical calculations have proposed that solvent mediates the protonation step. However, visualization of a proton transfer has been elusive. Based on a 2.1 Å resolution neutron structure of a pseudo-Michaelis complex of E. coli DHFR determined at acidic pH, we report the direct observation of the catalytic proton and its parent solvent molecule. Comparison of X-ray and neutron structures elucidated at acidic and neutral pH reveals dampened dynamics at acidic pH, even for the regulatory Met20 loop. Guided by the structures and calculations, we propose a mechanism where dynamics are crucial for solvent entry and protonation of substrate. This mechanism invokes the release of a sole proton from a hydronium (H₃O⁺) ion, its pathway through a narrow channel that sterically hinders the passage of water, and the ultimate protonation of DHF at the N5 atom.

Quantum and classical effects in DNA point mutations: Watson-Crick tautomerism in AT and GC base pairs

Proton transfer along the hydrogen bonds of DNA can lead to the creation of short-lived, but biologically relevant point mutations that can further lead to gene mutation and, potentially, cancer. In this work, the energy landscape of the canonical A-T and G-C base pairs (standard, amino-keto) to tautomeric A*-T* and G*-C* (non-standard, imino-enol) Watson-Crick DNA base pairs is modelled with density functional theory and machine-learning nudge-elastic band methods. We calculate the energy barriers and tunnelling rates of hydrogen transfer between and within each base monomer (A, T, G and C). We show that the role of tunnelling in A-T tautomerisation is statistically unlikely due to the presence of a small reverse reaction barrier. On the contrary, the thermal populations of the G*-C* point mutation could be non-trivial and propagate through the replisome. For the direct intramolecular transfer, the reaction is hindered by a substantial energy barrier. However, our calculations indicate that tautomeric bases in their monomeric form have remarkably long lifetimes.

Ring polymer quantization of the photon field in polariton chemistry

We use the ring polymer (RP) representation to quantize the radiation field inside an optical cavity to investigate polariton quantum dynamics. Using a charge transfer model coupled to an optical cavity, we demonstrate that the RP quantization of the photon field provides accurate rate constants of the polariton mediated electron transfer reaction compared to Fermi's golden rule. Because RP quantization uses extended phase space to describe the photon field, it significantly reduces the computational costs compared to the commonly used Fock state description of the radiation field. Compared to the other quasi-classical descriptions of the photon field, such as the classical Wigner based mean-field Ehrenfest model, the RP representation provides a much more accurate description of the polaritonic quantum dynamics because it alleviates the potential quantum distribution leakage problem associated with the photonic degrees of freedom (DOF). This work demonstrates the possibility of using the ring polymer description to treat the quantized radiation field in polariton chemistry, offering an accurate and efficient approach for future investigations in cavity quantum electrodynamics.

Dynamical strengthening of covalent and non-covalent molecular interactions by nuclear quantum effects at finite temperature

Nuclear quantum effects (NQE) tend to generate delocalized molecular dynamics due to the inclusion of the zero point energy and its coupling with the anharmonicities in interatomic interactions. Here, we present evidence that NQE often enhance electronic interactions and, in turn, can result in dynamical molecular stabilization at finite temperature. The underlying physical mechanism promoted by NQE depends on the particular interaction under consideration. First, the effective reduction of interatomic distances between functional groups within a molecule can enhance the n → π* interaction by increasing the overlap between molecular orbitals or by strengthening electrostatic interactions between neighboring charge densities. Second, NQE can localize methyl rotors by temporarily changing molecular bond orders and leading to the emergence of localized transient rotor states. Third, for noncovalent van der Waals interactions the strengthening comes from the increase of the polarizability given the expanded average interatomic distances induced by NQE. The implications of these boosted interactions include counterintuitive hydroxyl-hydroxyl bonding, hindered methyl rotor dynamics, and molecular stiffening which generates smoother free-energy surfaces. Our findings yield new insights into the versatile role of nuclear quantum fluctuations in molecules and materials.

Quantum chemical accuracy from density functional approximations via machine learning

Kohn-Sham density functional theory (DFT) is a standard tool in most branches of chemistry, but accuracies for many molecules are limited to 2-3 kcal ⋅ mol-1 with presently-available functionals. Ab initio methods, such as coupled-cluster, routinely produce much higher accuracy, but computational costs limit their application to small molecules. In this paper, we leverage machine learning to calculate coupled-cluster energies from DFT densities, reaching quantum chemical accuracy (errors below 1 kcal ⋅ mol-1) on test data. Moreover, density-based Δ-learning (learning only the correction to a standard DFT calculation, termed Δ-DFT ) significantly reduces the amount of training data required, particularly when molecular symmetries are included. The robustness of Δ-DFT  is highlighted by correcting "on the fly" DFT-based molecular dynamics (MD) simulations of resorcinol (C6H4(OH)2) to obtain MD trajectories with coupled-cluster accuracy. We conclude, therefore, that Δ-DFT  facilitates running gas-phase MD simulations with quantum chemical accuracy, even for strained geometries and conformer changes where standard DFT fails.

Quantum effects and 1H NMR chemical shifts of a bifurcated short hydrogen bond

The monoprotonated compound N,N',N''-tris(p-tolyl)azacalix[3](2,6)pyridine (TAPH) contains an intramolecular hydrogen bond that is formed from three N atoms in its cavity. Constrained by the macrocyclic molecular structure, the separations between the N atoms in this bifurcated hydrogen bond are about 2.6 Å, considerably shorter than those typically observed for hydrogen bonded systems in the condensed phases. As such, TAPH exhibits significantly elongated N-H lengths in its hydrogen bond and a downfield ¹H NMR chemical shift of 22.1 ppm. In this work, we carry out ab initio molecular dynamics and ab initio path integral molecular dynamics simulations of TAPH in the acetonitrile solution to reveal the geometry and proton sharing conditions of the bifurcated short hydrogen bond and uncover how the interplay of electronic and nuclear quantum effects gives rise to its far downfield ¹H chemical shift. Taking a linear short hydrogen bond as a reference, we demonstrate the distinct features of competing quantum effects and electronic shielding effects in the bifurcated hydrogen bond of TAPH. We further use the degree of deshielding on the proton as a measure of the hydrogen bonding interactions and evaluate the strength of the bifurcated short hydrogen bond as compared to its linear counterpart.

A New Bias Site for Epigenetic Modifications: How Non-Canonical GC Base Pairs Favor Mechanochemical Cleavage of DNA

Properties of non-canonical GC base pairs and their relations with mechanochemical cleavage of DNA are analyzed. A hypothesis of the involvement of the transient GC wobble base pairs both in the mechanisms of the mechanochemical cleavage of DNA and epigenetic mechanisms involving of 5-methylcytosine, is proposed. The hypothesis explains the increase in the frequency of the breaks of the sugar-phosphate backbone of DNA after cytosines, the asymmetric character of these breaks, and an increase in break frequency in CpG after cytosine methylation. As an alternative hypothesis, probable implication of GC⁺ Hoogsteen base pairs is considered, which now exemplify the best-studied non-canonical GC base pairs in the DNA double helix. Also see the video abstract here https://youtu.be/EUunVWL0ptw.

Quantum kinetic energy and isotope fractionation in aqueous ionic solutions

At room temperature, the quantum contribution to the kinetic energy of a water molecule exceeds the classical contribution by an order of magnitude. The quantum kinetic energy (QKE) of a water molecule is modulated by its local chemical environment and leads to uneven partitioning of isotopes between different phases in thermal equilibrium, which would not occur if the nuclei behaved classically. In this work, we use ab initio path integral simulations to show that QKEs of the water molecules and the equilibrium isotope fractionation ratios of the oxygen and hydrogen isotopes are sensitive probes of the hydrogen bonding structures in aqueous ionic solutions. In particular, we demonstrate how the QKE of water molecules in path integral simulations can be decomposed into translational, rotational and vibrational degrees of freedom, and use them to determine the impact of solvation on different molecular motions. By analyzing the QKEs and isotope fractionation ratios, we show how the addition of the Na⁺, Cl^- and HPO₄^2- ions perturbs the competition between quantum effects in liquid water and impacts their local solvation structures.

Tautomerism of Guanine Analogues

Tautomerism of nucleic acid (NA) bases is a crucial factor for the maintenance and translation of genetic information in organisms. Only canonical tautomers of NA bases can form hydrogen-bonded complexes with their natural counterparts. On the other hand, rare tautomers of nucleobases have been proposed to be involved in processes catalysed by NA enzymes. Isocytosine, which can be considered as a structural fragment of guanine, is known to have two stable tautomers both in solution and solid states. The tautomer equilibrium of isocytosine contrasts with the remarkable stability of the canonical tautomer of guanine. This paper investigates the factors contributing to the stability of the canonical tautomer of guanine by a combination of NMR experiments and theoretical calculations. The electronic effects of substituents on the stability of the rare tautomers of isocytosine and guanine derivatives are studied by density functional theory (DFT) calculations. Selected derivatives are studied by variable-temperature NMR spectroscopy. Rare tautomers can be stabilised in solution by intermolecular hydrogen-bonding interactions with suitable partners. These intermolecular interactions give rise to characteristic signals in proton NMR spectra, which make it possible to undoubtedly confirm the presence of a rare tautomer.

A combined on-the-fly/interpolation procedure for evaluating energy values needed in molecular simulations

We propose an algorithm for molecular dynamics or Monte Carlo simulations that uses an interpolation procedure to estimate potential energy values from energies and gradients evaluated previously at points of a simplicial mesh. We chose an interpolation procedure that is exact for harmonic systems and considered two possible mesh types: Delaunay triangulation and an alternative anisotropic triangulation designed to improve performance in anharmonic systems. The mesh is generated and updated on the fly during the simulation. The procedure is tested on two-dimensional quartic oscillators and on the path integral Monte Carlo evaluation of the HCN/DCN equilibrium isotope effect.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Modeling the action of environment on proton tunneling in the adenine-thymine base pair

Effect of environment coupling on the quantum-biological phenomenon of proton tunneling in the hydrogen bonds of the adenine-thymine base pair in DNA was modeled within the framework of quantum statistics and perturbation theory. A number of important thermodynamic indicators including partition function, free energy, and entropy were then calculated and examined. The proton was then assumed to be subject to an attraction represented by a double-well potential energy surface with a small asymmetry, which was considered as the perturbation introduced to the system. The action of environment manifested itself in the form of a global minimum in the free energy curve, as an implicit implication of the tendency of the system toward randomness and disorder, at which no spontaneous change such as quantum tunneling will accordingly occur. Furthermore, assuming the free energy to be in a close neighborhood of its minimum truly explained the smallness of the contribution of environment coupling to the tunneling probability reported in the literature based on the fact that the closer the free energy to its minimum, the less the transition probability to this point.

Unraveling the structural and chemical features of biological short hydrogen bonds

Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.

The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.

From classical to quantum and back: Hamiltonian adaptive resolution path integral, ring polymer, and centroid molecular dynamics

Path integral-based methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, including path-integral molecular dynamics, which allows for the calculation of quantum statistical properties, and ring-polymer and centroid molecular dynamics, which allow the calculation of approximate quantum dynamical properties. To this end, we derive a new integration algorithm that also makes use of multiple time-stepping. The scheme is validated via adaptive classical-path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.

Proton transfer in guanine-cytosine base pair analogues studied by NMR spectroscopy and PIMD simulations

It has been hypothesised that proton tunnelling between paired nucleobases significantly enhances the formation of rare tautomeric forms and hence leads to errors in DNA replication. Here, we study nuclear quantum effects (NQEs) using deuterium isotope-induced changes of nitrogen NMR chemical shifts in a model base pair consisting of two tautomers of isocytosine, which form hydrogen-bonded dimers in the same way as the guanine-cytosine base pair. Isotope effects in NMR are consequences of NQEs, because ro-vibrational averaging of different isotopologues gives rise to different magnetic shielding of the nuclei. The experimental deuterium-induced chemical shift changes are compared with those calculated by a combination of path integral molecular dynamics (PIMD) simulations with DFT calculations of nuclear shielding. These calculations can directly link the observable isotope-induced shifts with NQEs. A comparison of the deuterium-induced changes of 15N chemical shifts with those predicted by PIMD simulations shows that inter-base proton transfer reactions do not take place in this system. We demonstrate, however, that NMR isotope shifts provide a unique possibility to study NQEs and to evaluate the accuracy of the computational methods used for modelling quantum effects in molecules. Calculations based on the PBE functional from the general-gradient-approximation family provided significantly worse predictions of deuterium isotope shifts than those with the hybrid B3LYP functional.

Quantum mechanical studies on the existence of AC mismatches through prototopic tautomerization pathway in adenine and cytosine recognition

The gas phase studies on the electronic stabilities and thermochemical properties of metastable AC mismatches have been performed and the results are analyzed to explain the features of experimentally available AC mismatches. The hydrogen bonding patterns observed in these mismatches are relevant to the formation of stable AC mismatches. In these AC pairs, the H migration mechanism to generate other tautomeric forms is not observed, which shows the compatibility of H bonding capacity of the sites involved in H bonds. The presence of hydrogen bond type –N[Formula: see text]H–N– may contribute to better AC pairing, hence cisA1-1cisC2, transA2-cisC1 and transA1-1cisC2 are found to be more favorable pairs compared to other AC pairs.

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

Proton Transfer and Structure-Specific Fluorescence in Hydrogen Bond-Rich Protein Structures

Protein structures which form fibrils have recently been shown to absorb light at energies in the near UV range and to exhibit a structure-specific fluorescence in the visible range even in the absence of aromatic amino acids. However, the molecular origin of this phenomenon has so far remained elusive. Here, we combine ab initio molecular dynamics simulations and fluorescence spectroscopy to demonstrate that these intrinsically fluorescent protein fibrils are permissive to proton transfer across hydrogen bonds which can lower electron excitation energies and thereby decrease the likelihood of energy dissipation associated with conventional hydrogen bonds. The importance of proton transfer on the intrinsic fluorescence observed in protein fibrils is signified by large reductions in the fluorescence intensity upon either fully protonating, or deprotonating, the fibrils at pH = 0 or 14, respectively. Thus, our results point to the existence of a structure-specific fluorophore that does not require the presence of aromatic residues or multiple bond conjugation that characterize conventional fluorescent systems. The phenomenon may have a wide range of implications in biological systems and in the design of self-assembled functional materials.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Stability of Complex Biomolecular Structures: van der Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects

Biomolecules are complex systems stabilized by a delicate balance of weak interactions, making it important to assess all energetic contributions in an accurate manner. However, it is a priori unclear which contributions make more of an impact. Here, we examine stacked polyglutamine (polyQ) strands, a peptide repeat often found in amyloid aggregates. We investigate the role of hydrogen bond (HB) cooperativity, van der Waals (vdW) dispersion interactions, and quantum contributions to free energies, including anharmonicities through density functional theory and ab initio path integral simulations. Of these various factors, we find that the largest impact on structural stabilization comes from vdW interactions. HB cooperativity is the second largest contribution as the size of the stacked chain grows. Competing nuclear quantum effects make the net quantum contribution small but very sensitive to anharmonicities, vdW, and the number of HBs. Our results suggest that a reliable treatment of these systems can only be attained by considering all of these components.

Modeling quantum nuclei with perturbed path integral molecular dynamics

Here we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics for modeling quantum nuclear effects.

The quantum nature of nuclear motions plays a vital role in the structure, stability, and thermodynamics of molecules and materials. The standard approach to model nuclear quantum fluctuations in chemical and biological systems is to use path-integral molecular dynamics. Unfortunately, conventional path-integral simulations can have an exceedingly large computational cost due to the need to employ an excessive number of coupled classical subsystems (beads) for quantitative accuracy. Here, we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics and derive an improved non-empirical partition function and estimators to calculate converged quantum observables. Our perturbed path-integral (PPI) method requires the same ingredients as the conventional approach, but increases the accuracy and efficiency of path integral simulations by an order of magnitude. Results are presented for the thermodynamics of fundamental model systems, an empirical water model containing 256 water molecules within periodic boundary conditions, and ab initio simulations of nitrogen and benzene molecules. For all of these examples, PPI simulations with 4 to 8 classical beads recover the nuclear quantum contribution to the total energy and heat capacity at room temperature within a 3% accuracy, paving the way toward seamless modeling of nuclear quantum effects in realistic molecules and materials.

Modelling proton tunnelling in the adenine–thymine base pair

Environment-assisted quantum tunnelling is unlikely to play a significant role in adenine–thymine tautomerisation.

Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site

Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

Direct observation of multiple tautomers of oxythiamine and their recognition by the thiamine pyrophosphate riboswitch

Structural diversification of canonical nucleic acid bases and nucleotide analogues by tautomerism has been proposed to be a powerful on/off switching mechanism allowing regulation of many biological processes mediated by RNA enzymes and aptamers. Despite the suspected biological importance of tautomerism, attempts to observe minor tautomeric forms in nucleic acid or hybrid nucleic acid-ligand complexes have met with challenges due to the lack of sensitive methods. Here, a combination of spectroscopic, biochemical, and computational tools probed tautomerism in the context of an RNA aptamer-ligand complex; studies involved a model ligand, oxythiamine pyrophosphate (OxyTPP), bound to the thiamine pyrophosphate (TPP) riboswitch (an RNA aptamer) as well as its unbound nonphosphorylated form, oxythiamine (OxyT). OxyTPP, similarly to canonical heteroaromatic nucleic acid bases, has a pyrimidine ring that forms hydrogen bonding interactions with the riboswitch. Tautomerism was established using two-dimensional infrared (2D IR) spectroscopy, variable temperature FTIR and NMR spectroscopies, binding isotope effects (BIEs), and computational methods. All three possible tautomers of OxyT, including the minor enol tautomer, were directly identified, and their distributions were quantitated. In the bound form, BIE data suggested that OxyTPP existed as a 4'-keto tautomer that was likely protonated at the N1'-position. These results also provide a mechanistic framework for understanding the activation of riboswitch in response to deamination of the active form of vitamin B1 (or TPP). The combination of methods reported here revealing the fine details of tautomerism can be applied to other systems where the importance of tautomerism is suspected.