To advance our understanding of why ribose takes on the furanose form, rather than the pyranose form, in RNA, Alonso and co-workers have examined the structure of uridine 1 in the gas phase.¹

1

Uridine is sensitive to temperature, and so the laser-ablation method long used by the Alonso group is ideal for examining uridine. The microwave spectrum is quite complicated due to the presence of many photofragments. Careful analysis lead to the identification of a number of lines and hyperfine structure that could be definitively assigned to uridine, leading to experimental values of the rotational constants and the diagonal elements of the ¹⁴N nuclear quadrupole coupling tensor for each nitrogen. These values are listed in Table 1.

Table 1. Experimental and calculated rotational constants (MHz), quadrupole coupling constants (MHz) and relative energy (kcal mol^-1).

In order to assign a 3-D structure to these experimental values, they examined the PES of uridine with molecular mechanics and semi-empirical methods, before reoptimizing the structure of the lowest 5 energy structures at MP2/6-311++G(d,p). Then, comparison of the resulting rotational constants and ¹⁴N nuclear quadrupole coupling constants of these computed structures (see Table 1) led to identification of the lowest energy structure (anti/C2’-endo-g+, see Figure 1) in best agreement with the experiment. Once again, the Alonso group has demonstrated the value of the synergy between experiment and computation in structure identification.

Figure 1. MP2/6-311++G(d,p) optimized structure of 1 (anti/C2’-endo-g+).

References

(1) Peña, I.; Cabezas, C.; Alonso, J. L. "The Nucleoside Uridine Isolated in the Gas Phase," Angew. Chem. Int. Ed. 2015, 54, 2991-2994, DOI: 10.1002/anie.201412460.

Inchis:

1: Inchi=1S/C9H12N2O6/c12-3-4-6(14)7(15)8(17-4)11-2-1-5(13)10-9(11)16/h1-2,4,6-8,12,14-15H,3H2,(H,10,13,16)/t4-,6-,7-,8-/m1/s1
InChiKey=DRTQHJPVMGBUCF-XVFCMESISA-N

Alonso and coworkers have again (see this post employed laser-ablation molecular-beam Fourier-transform microwave (LA-MB-MW)spectroscopy to discern the gas phase structure of an important biological compound: cytosine.¹ They identified five tautomers of cytosine 1-5. Comparison between the experimental and computational (MP2/6-311++G(d,p) microwave rotational constants and nitrogen nuclear quadrupole coupling constants led to the complete assignment of the spectra. The experimental and calculated rotational constants are listed in Table 1.

Table 1. Rotational constants (MHz) for 1-5.

The experimental and computed relative free energies are listed in Table 2. There is both not a complete match of the relative energetic ordering of the tautomers, nor is there good agreement in their magnitude. Previous computations² at CCSD(T)/cc-pVQZ//CCSD//cc-pVTZ are in somewhat better agreement with the gas-phase experiments.

Table 2. Relative free energies (kcal mol^-1) of 1-5.

References

(1) Alonso, J. L.; Vaquero, V.; Peña, I.; López, J. C.; Mata, S.; Caminati, W. "All Five Forms of Cytosine Revealed in the Gas Phase," Angew. Chem. Int. Ed. 2013, 52, 2331-2334, DOI: 10.1002/anie.201207744.

(2) Bazso, G.; Tarczay, G.; Fogarasi, G.; Szalay, P. G. "Tautomers of cytosine and their excited electronic states: a matrix isolation spectroscopic and quantum chemical study," Phys. Chem. Chem. Phys., 2011, 13, 6799-6807, DOI:10.1039/C0CP02354J.

InChIs

cytosine: InChI=1S/C4H5N3O/c5-3-1-2-6-4(8)7-3/h1-2H,(H3,5,6,7,8)
InChIKey=OPTASPLRGRRNAP-UHFFFAOYSA-N

Here’s another fine paper from the Alonso group employing laser ablation molecular beam Fourier transform microwave spectroscopy coupled with computation to discern molecular structure. In this work they examine the low-energy tautomers of guanine.1 The four lowest energy guanine tautomers are shown in Figure 1. (Unfortunately, Alonso does not include the optimized coordinates of these structures in the supporting information – we need to more vigorously police this during the review process!) These tautomers are predicted to be very close in energy (MP2/6-311++G(d,p), and so one might expect to see multiple signals in the microwave originating from all four tautomers. In fact, they discern all four, and the agreement between the computed and experimental rotational constants are excellent (Table 1), especially if one applies a scaling factor of 1.004. Once again, this group shows the power of combined experiment and computations!

Figure 1. Four lowest energy (kcal mol^-1, MP2/6-311++G(d,p)) tautomers of guanine.

Table 1. Experimental and computed rotational constants (MHz) of the four guanine tautomers.

References

(1) Alonso, J. L.; Peña, I.; López, J. C.; Vaquero, V., "Rotational Spectral Signatures of Four Tautomers of Guanine," Angew. Chem. Int. Ed. 2009, 48, 6141-6143, DOI: 10.1002/anie.200901462

InChIs

Guanine: InChI=1/C5H5N5O/c6-5-9-3-2(4(11)10-5)7-1-8-3/h1H,(H4,6,7,8,9,10,11)/f/h8,10H,6H2
InChIKey=UYTPUPDQBNUYGX-GSQBSFCVCX

Here is another benchmark of the performance of DFT in handling difficult situations, in this case the interaction between nucleic acid base pairs. Sherrill¹ has examined the 124 nucleic acid base pairs from the JSCH-2005 database² compiled by Hobza and coworkers. This database includes 36 hydrogen bonded complexes, and example of which is shown in Figure 1a, and 54 stacked complex, one example of which is shown in Figure 1b.

(a)

(b)

Figure 1. Optimized geometries (RI-MP2/cc-pVTZ) of two representative structures of base pairs: (a) hydrogen bonded pair and (c) stacked pair.

The energies of these base pairs computed with four different functionals: PBE, PBE-D (where Grimme’s empirical dispersion correction³), and the recently developed MO5-2X⁴ and MO6-2X⁵ methods which attempt to treat mid-range electron correlation. The aug-cc-pVDZ basis set was used. These DFT energies are compared with the CCSD(T) energies of Hobza. The mean unsigned error (MUE) for the 28 hydrogen bonded complexes and the 54 stacked complexes are listed in Table 1.

Table 1. Mean unsigned error (kcal mol^-1) of the four DFT
methods (relative to CCSD(T)) for the hydrogen bonded and stacked base pairs.

A few interesting trends are readily apparent. First, PBE (representing standard GGA DFT methods) poorly describes the energy of the hydrogen bonded complexes, but utterly fails to treat the stacking interaction. Inclusion of the dispersion correction (PBE-D) results in excellent energies for the HB cases and quite reasonable results for the stacked pairs. Both of Truhlar’s functionals dramatically outperform PBE, though MO5-2X is probably still not appropriate for the stacked case. MO6-2X however seems to be a very reasonable functional for dealing with base pair interactions, indicating that mid-range correlation correction is sufficient to describe these complexes, and that the long-range correlation correction included in the dispersion correction, while giving improved results, is not essential.

References

(1) Hohenstein, E. G.; Chill, S. T.; Sherrill, C. D., "Assessment of the Performance of the M05-2X and M06-2X Exchange-Correlation Functionals for Noncovalent Interactions in Biomolecules," J. Chem. Theory Comput., 2008, 4, 1996-2000, DOI: 10.1021/ct800308k

(2) Jurecka, P.; Sponer, J.; Cerny, J.; P., H., "Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs," Phys. Chem. Chem. Phys., 2006, 8, 1985-1993, DOI: 10.1039/b600027d.

(3) Grimme, S., "Semiempirical GGA-type density functional constructed with a long-range dispersion correction," J. Comput. Chem., 2006, 27, 1787-1799, DOI: 10.1002/jcc.20495

(4) Zhao, Y.; Schultz, N. E.; Truhlar, D. G., "Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions," J. Chem. Theory Comput., 2006, 2, 364-382, DOI: 10.1021/ct0502763.

(5) Zhao, Y.; Truhlar, D. G., "The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals," Theor. Chem. Acc., 2008, 120, 215-241, DOI: 10.1007/s00214-007-0310-x.