Nuclear Quantum Effects in Proton or Hydrogen Transfer

Proton or hydrogen transfers, basic chemical reactions, proceed either by thermally activated barrier crossing or via tunneling. Studies of molecules undergoing single or double proton or hydrogen transfer in the ground or excited electronic state reveal that tunneling can dominate under conditions usually considered to favor the thermal process. Moreover, the tunneling probability strongly varies for excitation of certain vibrational modes, which changes the effective barrier and/or proton transfer distance. When the reaction is fast compared to vibrational relaxation, the mode selectivity can still be maintained for molecules in solutions at 293 K. These observations point to dangers of relating the calculated minimum energy paths and the associated barriers to the experimentally obtained activation energies. The multidimensional character of the reaction coordinate is obvious; it can dramatically change for slowly and rapidly relaxing environments. We postulate that the hydrogen bond definition should be extended by specifically including the role of molecular vibrations.

Influence of local microenvironment on the double hydrogen transfer in porphycene

We performed time-resolved transient absorption and fluorescence anisotropy measurements in order to study tautomerization of porphycene in rigid polymer matrices at cryogenic temperatures. Studies were carried out in poly(methyl methacrylate) (PMMA), poly(vinyl butyral) (PVB), and poly(vinyl alcohol) (PVA). The results prove that in all studied media hydrogen tunnelling plays a significant role in the double hydrogen transfer which becomes very sensitive to properties of the environment below approx. 150 K. We also demonstrate that there exist two populations of porphycene molecules in rigid media: "hydrogen-transferring" molecules, in which tautomerization occurs on time scales below 1 ns and "frozen" molecules in which double hydrogen transfer is too slow to be monitored with nanosecond techniques. The number of "frozen" molecules increases when the sample is cooled. We explain this effect by interactions of guest molecules with a rigid host matrix which disturbs symmetry of porphycene and hinders tunnelling. Temperature dependence of the number of hydrogen-transferring molecules suggests that the factor which restores the symmetry of the double-minimum potential well in porphycene are intermolecular vibrations localized in separated regions of the amorphous polymer.

Generation of Tautomers Using Micro-p Ka's

Solutions of organic molecules containing one or more heterocycles with conjugated bonds may exist as a mixture of tautomers, but typically only a few of them are significantly populated even though the potential number grows combinatorially with the number of protonation and deprotonation sites. Generating the most stable tautomers from a given input structure is an important and challenging task, and numerous algorithms to tackle it have been proposed in the literature. This work describes a novel approach for tautomer prediction that involves the combined use of molecular mechanics, semiempirical quantum chemistry, and density functional theory. The key idea in our method is to identify the protonation and deprotonation sites using estimated micro-p K_a's for every atom in the molecule as well as in its nearest protonated and deprotonated forms. To generate tautomers in a systematic way with minimal bias, we then consider the full set of tautomers that arise from the combinatorial distribution of all such mobile protons among all protonatable sites, with efficient postprocessing to screen away high-energy species. To estimate the micro-p K_a's, we present a new method designed for the current task, but we emphasize that any alternative method can be used in conjunction with our basic algorithm. Our approach is therefore grounded in the computational prediction of physical properties in aqueous solution, in contrast to other approaches that may rely on the use of hard-coded rules of proton distribution, previously observed tautomerization patterns from a known chemical space, or human input. We present examples of the application of our algorithm to organic and drug-like molecules, with a focus on novel structures where traditional methods are expected to perform worse.