Next Article in Journal
Direct Regioselective C-H Cyanation of Purines
Previous Article in Journal
Application of Natural Deep Eutectic Solvents for Extraction of Bioactive Components from Rhodiola rosea (L.)
Previous Article in Special Issue
The Conformations of Isolated Gallic Acid: A Laser-Ablation Rotational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks

by
Vincenzo Barone
1,*,
Silvia Di Grande
1,2 and
Cristina Puzzarini
3,*
1
Scuola Normale Superiore, Piazza dei Cavalieri 7, I-50126 Pisa, Italy
2
Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
3
Rotational and Computational Spectroscopy Lab, Department of Chemistry “Giacomo Ciamician”, University of Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(2), 913; https://doi.org/10.3390/molecules28020913
Submission received: 1 January 2023 / Revised: 1 January 2023 / Accepted: 12 January 2023 / Published: 16 January 2023
(This article belongs to the Special Issue Structural Characterization of Biomolecules by Spectroscopic Methods)

Abstract

:
The interplay of high-resolution rotational spectroscopy and quantum-chemical computations plays an invaluable role in the investigation of biomolecule building blocks in the gas phase. However, quantum-chemical methods suffer from unfavorable scaling with the dimension of the system under consideration. While a complete characterization of flexible systems requires an elaborate multi-step strategy, in this work, we demonstrate that the accuracy obtained by quantum-chemical composite approaches in the prediction of rotational spectroscopy parameters can be approached by a model based on density functional theory. Glycine and serine are employed to demonstrate that, despite its limited cost, such a model is able to predict rotational constants with an accuracy of 0.3% or better, thus paving the way toward the accurate characterization of larger flexible building blocks of biomolecules.

1. Introduction

The experimental study of biomolecule building blocks in the gas phase has recently attracted increasing attention owing to the development of spectrometers coupling supersonic-jet expansion [1] with laser ablation [2]. This has allowed the application of high-resolution rotational spectroscopy in the microwave region to thermolabile molecules with high melting points. Among these, we find the main families of life bricks such as amino acids, peptides, DNA and RNA bases. The focus of the present contribution is the computational support to rotational spectroscopy investigations of building blocks of biomolecules.
For flexible molecules such as those of interest for this study, a great challenge is related to the large number of conformers and the fast relaxation of some of them to more stable counterparts due to the presence of low interconversion energy barriers. An inaccurate account of the relaxation processes can bias any direct thermochemical interpretation of the results provided by rotational spectroscopy experiments [3,4]. Quantum-chemical (QC) computations can help to tackle this challenge, especially because gas phase is their most natural playground [5,6]. Unfortunately, already for medium-sized systems, the usual dichotomy between accuracy and feasibility, which is the quest for accurate yet feasible predictions, comes into place [7]. In fact, state-of-the-art QC approaches are able to rival the experimental counterparts for small semi-rigid systems in the gas phase [7,8,9], but they are characterized by a very unfavorable scaling with the dimension of the system to be investigated. This prevents their brute-force application already to biomolecule building blocks containing more than a dozen atoms and characterized by several low-energy minima. Furthermore, the powerful local optimization techniques developed for semi-rigid systems are not effective for flexible systems, which requires the exploration of rugged potential energy surfaces (PESs) [10,11].
For the reasons mentioned above, the accurate characterization needed by rotational spectroscopy requires an integrated computational approach that employs QC models of increasing accuracy in the different steps of an exploration/exploitation strategy guided by machine learning (ML) tools. The main steps of this strategy [10,12,13,14] can be summarized as follows:
  • Unsupervised perception of the molecular system to identify hard and soft degrees of freedom [15];
  • Knowledge-based selection and constrained geometry optimizations of a limited number of conformers employing a fast semi-empirical method [11,16];
  • Exploration of the PES governed by soft degrees of freedom using the same semi-empirical method of the previous step, guided by a purposely tailored evolutionary algorithm with the aim of finding other low-lying minima [10];
  • Refinement of the most stable structures by hybrid and then double-hybrid density functionals [14];
  • Analysis of relaxation paths between pairs of adjacent energy minima [13];
  • Evaluation of accurate electronic energies for the final panel of low-energy minima [17,18,19];
  • Computation of zero point energies (ZPE) and thermal contributions to enthalpies and entropies [20,21,22,23,24,25];
  • Computation of spectroscopic parameters for the energy minima with non negligible populations [13,26].
The focus of the present contribution is on points 6 to 8. First, we will validate the accuracy of approximate approaches based on double-hybrid functionals and effective composite methods by exploiting results available for semi-rigid molecules and non-covalent complexes. Next, amino acids will be employed to illustrate the potentialities of the validated approaches. These compounds represent a particularly appealing playground because their rich conformational landscape is tuned by the competition among different types of intra-molecular non-covalent interactions involving, together with the amino and carboxylic acid moieties of the backbone, also side-chain groups. At the same time, results from MW experiments are available for several conformers of most natural α -amino acids [27,28,29,30,31,32,33,34,35,36,37,38], and provide accurate data for benchmarking theory. Among the α -amino acids well characterized experimentally, we have selected glycine and serine.
Since the current standard for the computation of rotational spectroscopy parameters of biomolecule building blocks (see, for example, [4,27,39,40,41]) is based on QC methodologies of limited accuracy and does not account for vibrational effects, together with the intrinsic interest of the studied molecules, the results of the present study will provide a reference for a more accurate and reliable approach to be used in the prediction and analysis of MW experiments.

2. Results and Discussion

2.1. The Methodologic Approach

Several studies have shown that the double-hybrid rev-DSD-PBEP86 functional [42] in conjunction with the jun-cc-pVTZ [43] basis set (this functional-basis combination is shortly denoted as rDSD in the following) represents a very robust and convenient choice for obtaining equilibrium geometries suitable for subsequent electronic energy evaluations [19,44,45,46]. For the latter, composite schemes rooted in the coupled-cluster (CC) theory [47] that account for the extrapolation to the complete basis set (CBS) limit and core–valence (CV) correlation are employed whenever aiming at high accuracy [13,25,46,48,49,50]. However, composited schemes entirely based on CC theory are computationally expensive. The key idea of the so-called ‘cheap’-scheme (ChS) family of composite approaches [17,18,19,51] is that both CBS and CV contributions can be estimated accurately by low-order (hence low-cost) many-body perturbative methods, with Møller–Plesset theory to second order (MP2) [52] offering the best compromise between accuracy and computational cost. The general expression of the electronic energy evaluated using the ChS model is:
E ChS = E ( CCSD ( T ) / TZ ) + Δ E ( MP 2 / CBS ) + Δ E ( MP 2 / CV ) ,
where the first term on the right-hand side is the energy at the fc-CCSD(T) [53] level in conjunction with a triple-zeta (TZ) basis set, where CCSD(T) denotes the CC method incorporating single, double, and a perturbative estimate of triple excitations, and fc stands for the frozen-core approximation. The second term is the contribution due to the extrapolation to the CBS limit evaluated at the MP2 level within the fc approximation and using the n 3 formula [54]:
Δ E ( MP 2 / CBS ) = n 3 E MP 2 / n Z ( n 1 ) 3 E MP 2 / ( n 1 ) Z n 3 ( n 1 ) 3 E MP 2 / ( n 1 ) Z ,
where n = 4 (quadruple-zeta basis set, QZ) and n-1 stands for TZ. The last term, Δ E ( MP 2 / CV ) , is the CV contribution computed as energy difference between MP2 calculations correlating all electrons and within the fc approximation, both in the cc-p(w)CVTZ basis set [55,56].
The first implementation of the ChS model employed the cc-pVnZ basis sets [57]. However, several benchmarks [18,19] have led to the conclusion that the best cost/performance ratio is obtained with the jun-cc-pVnZ basis sets [43], with jun-cc-pV(n+d)Z being used for atoms of the third period [58]. This model, hereafter referred to as junChS, has become our standard ChS methodology. A more recent alternative [19] replaces the conventional CCSD(T) and MP2 methods with their explicitly-correlated counterparts, namely CCSD(T)-F12 and MP2-F12 [59,60], still in conjunction with the jun-cc-pVnZ and cc-p(w)CVTZ basis sets. Hereafter, this model is referred to as junChSF12. In all CCSD(T)-F12 computations, the F12b approximation [61] is employed, while the default approximation (3C/FiX) is used for MP2-F12 [60].
In addition to energy evaluations, different ChS variants have been employed and tested for geometry optimizations [19,51,62,63]. For conventional QC methods, ChS and junChS are defined analogous to Equations (1) and (2), whereas the extrapolation to the CBS limit (but not the CV correction) can be avoided (and sometimes introduces non-negligible oscillations) for approaches exploiting explicitly-correlated methods. The models derived and employed in the present work are denoted as DZCCF12+CV, TZCCF12+CV, junCCF12+CV and augCCF12+CV, which correspond to CCSD(T)-F12 calculations in conjunction with cc-pVDZ-F12, cc-pVTZ-F12, jun-cc-pVTZ and aug-cc-pVTZ, respectively. The “+CV” denotes the incorporation of the CV contribution, which is evaluated at the MP2-F12/cc-p(w)CVTZ level.
When dealing with flexible molecules, a key point is the determination of the relative stability of the low-energy minima to determine those that are sufficiently populated for a spectroscopic characterization. To accomplish this task, one has to move from electronic energy differences to the corresponding relative enthalpies or free energies, evaluated at a temperature that depends on the experimental conditions. To this end, ZPEs and vibrational partition functions are needed. Within the harmonic approximation, this task requires the evaluation of harmonic vibrational frequencies, which can be obtained with a suitable accuracy at the rDSD level [64] exploiting analytical second derivatives [65]. The models of the ChS family can also be applied to the computation of harmonic frequencies. However, due to effective error compensation, improved results are obtained by neglecting the CV correction for explicitly-correlated approaches [66]. Focusing on models that exploit the F12 methodology, the DZCCF12, TZCCF12, junCCF12 and augCCF12 approaches offer a good compromise between accuracy and computational cost. The ZPE term can then be improved by incorporating anharmonic corrections, which are usually obtained at the B3LYP/jun-cc-pVDZ level (hereafter B3) within second-order vibrational perturbation theory (VPT2, [20,67]). Whenever required, the contribution of low-frequency motions to entropies can be computed using the quasi-harmonic (QH) approximation [24,68]. Within the QH approximation, below a given cut-off value, entropic terms are obtained from the free-rotor approximation, and a damping function is used to interpolate between free-rotor and harmonic oscillator expressions close to the cut-off frequency.
Moving to rotational spectroscopy, the leading terms are the rotational constants, which incorporate two contributions. The first, and by far the largest, one is the equilibrium rotational constant, B e , which is straightforwardly derived from the equilibrium geometry of the studied molecule. The second contribution is related to vibrational effects which are present because of the coupling between vibrations and rotation. This contribution is usually approximated by means of VPT2, which, for the vibrational gound state, leads to:
B 0 i = B e i 1 2 r α r i = B e i + Δ B v i b i .
where α r i denotes the vibrational-rotation interaction constants, with i being the inertial axis (a, b, or c) and the sum running over all the r vibrational modes. The computation of this contribution requires both second- and semi-diagonal third-energy derivatives of the energy with respect to normal modes. Therefore, computations of vibrational corrections to rotational constants are quite expensive. As a consequence, its inclusion is warranted only if the error on the computed equilibrium rotational constants is smaller than the expected vibrational contribution. A recent careful analysis of these aspects provides some general hints [69]. The magnitude of the vibrational contribution Δ B v i b is typically 0.1% to 0.7% that of the corresponding equilibrium rotational constant, with 0.5% being a very robust guess for semi-rigid molecules. As a consequence, an uncertainty of 10% on Δ B v i b would lead to an error of 0.05% on the rotational constants, which is thus more than acceptable. The overall conclusion is the suitability of global hybrid density functionals in conjunction with double-zeta basis sets (e.g., the B3 model previously defined) for the computation of the required anharmonic force fields, thus largely reducing the computational cost.
Focusing on the accuracy needed by rotational spectroscopy, errors of 1% on rotational constants (100 MHz for a constant with magnitude 10 GHz) are not at all helpful for the prediction and/or analysis of rotational spectra of flexible molecules, which are always characterized by several low-energy minima with non negligible populations. On the other hand, an accuracy of 0.01%, which would allow direct comparison with experiment, is presently reachable only for molecules containing two or three atoms. Based on these considerations, the optimal level of accuracy associated with predicted rotational constants should be close to 0.1% (10 MHz for a constant of 10 GHz). In terms of structural parameters, such an accuracy corresponds to errors smaller than 0.001 Å for typical bond lengths and 0.001 radians (0.05 degrees) for typical valence angles [69]. This target accuracy can be surely obtained by expensive composite schemes incorporating high excitation orders in the correlation treatment [70]. Nevertheless, the models of the ChS family are able to draw closer to such a high accuracy limit [51,62,69]. However, for quite large flexible molecules, even this computational level becomes too expensive and one has to rely on density functional theory (DFT).
To improve the optimized geometries obtained from double-hybrid functionals, one can resort to the the so-called linear regression approach (LRA) [71,72,73,74]. This corrects the computed bond lengths ( r c o m p ) for systematic errors by means of scaling factors ( a ) and offset values ( b ) that have been derived from the comparison of DFT and accurate semi-experimental (SE) equilibrium geometries for a large database:
r = ( 1 + a ) × r c o m p + b ,
with the a and b parameters depending not only on the functional considered, but also on the nature of the atoms involved. These are available in [73].
As mentioned above, the equilibrium rotational constants contribute more than 99% to the parameters derived from experiments and only depend on the molecular structure and isotopic composition [75]. Indeed, rotational constants are inversely proportional to the inertia tensor, which contains information on the mass distribution in the molecule. Therefore, a comparison between computed and experimental rotational constants can provide hints on the quality of the computed structure. In this respect, additional parameters of particular relevance are the nuclear quadrupole coupling constants ( χ i i , i referring to the inertia axis a, b or c), which strongly depend on the intra-molecular interactions [76]. Nuclear quadrupole coupling is the interaction between the quadrupole moment of a quadrupolar nucleus and the electric gradient at the nucleus itself, with quadrupolar nuclei being those having a nuclear spin I 1 [75]. Since 14 N is a quadrupolar nucleus and is present in almost all biomolecule building blocks, nuclear quadrupole coupling constants can be used, together with rotational constants, to assess the quality of computed structures. In turn, their calculation is important for accurate predictions of rotational spectra because nuclear quadrupole coupling determines a splitting of the rotational energy levels, thus causing a splitting of the rotational transitions, which is the so-called hypefine structure. It is also noted that vibrational effects on nuclear quadrupole coupling constants are usually smaller than the uncertainty affecting the computed equilibrium values, and thus have not been considered in this work.
As far as technical details are concerned, all DFT, MP2 and single-point energy CCSD(T) computations have been performed using the Gaussian program package [77], while the CFOUR program [78] has been employed for geometry optimizations at the CCSD(T) level. For explicitly-correlated calculations, the MOLPRO program [79] has been used.

2.2. The Validation Step

Non-covalent complexes represent a suitable benchmark for the validation of computational methodologies aiming at an accurate description of intra-molecular non-covalent interactions. In a series of recent papers [18,19,80], different ChS variants have been applied to the computation of the interaction energies of prototypical systems. These studies (partially summarized in Table 1) show that the junChS and, especially, junChSF12 approaches allow the determination of interaction energies with an average error smaller than 10 cm 1 (0.12 kJ mol 1 ) without the need for incorporating any empirical parameter. In Table 1, the non-covalent complexes of the A14 database [19], which is a based on a selection of systems from the A24 dataset [81], have been considered. As reference values, the highly accurate data (which also incorporate up to quadruple excitation in the CC expansion and relativistic effects) from [82] have been selected.
From an inspection of Table 1, we note that going from junCCF12+CV to junChSF12, which only differ for the extrapolation to the CBS limit, leads to a significant lowering of both the maximum and mean unsigned error (MAX and MUE, respectively), with this confirming the importance of the extrapolation to the CBS limit also for explicitly-correlated methods and that this contribution can be accurately incorporated by means of low-order perturbative methods (here MP2-F12). Indeed, MP2-F12 represents a reliable route for obtaining accurate results without any significant computational increase with respect to the underlying CCSD(T)-F12 step.
Despite the extrapolation to the CBS limit, all ChS models (both conventional and F12 variants) are affected by small, but not entirely negligible basis set superposition errors (BSSE), which can be taken into account by means of counterpoise (CP) corrections (for details, the reader is referred to [19]). However, these are slightly more pronounced when considering the model based on conventional methods. This is particularly significant for intra-molecular interactions where the BSSE becomes ill defined and CP corrections difficult to implement. An important outcome of [18,19] is that, when using the rDSD level for the evaluation of the reference structures, non-CP corrected geometries can be safely employed. This is particularly important in view of applying junChS and junChSF12 to the evaluation of the energetics for flexible molecules. Concerning interaction energy, the errors at the rDSD level without any CP correction nearly double with respect to junChS(F12) [19], but they remain significantly smaller than those issued by other density functionals. More importantly, the rDSD level is suitable for defining a correct stability order for different kinds of non-covalent interactions.
As mentioned above, the largest contribution to rotational constants comes from their equilibrium value, which is related to the corresponding equilibrium geometry. In this connection, ChS, junChS and junChSF12 have proven to provide very accurate results [19,63], with the former two models being rather well tested [83]. When employing explicitly-correlated methods, the computational cost of the ChS methodology can be further reduced by omitting the CBS extrapolation step. To inspect this, a validation of the junChS and CCF12+CV models (for CCF12, we considered the DZCCF12, TZCCF12 and junCCF12 variants) has been performed by comparing their results with accurate SE equilibrium geometries from the SE100 database [73]. The subset considered includes 24 molecules containing hydrogen, second- and third-row atoms: CH 4 , CO 2 , HCN, HNC, H 2 O, NH 3 , C 2 H 2 , C 2 H 4 , H 2 CO, t-HCOOH, CH 2 NH, BH 3 NH 3 , BH 2 OH, C 2 H 4 O (oxirane), C 2 H 4 NH (aziridine), C 3 H 6 (cyclopropane), H 2 O 2 , SO 2 , H 2 S, PH 3 , H 2 CS, CH 2 PH, C 2 H 4 S (tiirane) and H 2 S 2 . The results, collected in Table 2 in terms of maximum and mean unsigned errors (MAX, MUE), show that the CV contribution cannot be neglected and that the jun-cc-pVTZ basis set is competitive with its counterparts (cc-pVDZF12 and cc-pVTZF12) developed specifically for explicitly correlated models.
Even if not reported in Table 2, from [73], it is observed that—for the entire SE100 set—the rDSD model delivers remarkably accurate geometrical parameters with typical absolute errors of about 0.003 Å for bond lengths and 0.1–0.2 degrees for valence angles. Furthermore, the errors appear to be very systematic. Therefore, the rDSD level is particularly suitable for the application of the LRA mentioned above, which greatly reduces the expected errors [73]. In particular, in view of the focus of this work, the accuracy for the corresponding rotational constants becomes close to the 0.1% uncertainty-target.
The last issue touched in this validation step concerns the computation of reliable ZPEs and thermodynamic functions. As already mentioned, this requires accurate vibrational frequencies, whose leading terms are the corresponding harmonic values. For the purpose of finding the best compromise between accuracy and computational cost, we have computed the harmonic frequencies of eight molecules (hereafter referred to as the H8 set) for which either accurate experimental values are available (H 2 O [84], HCN [85], CO 2 [86], C 2 H 2 [87]) or whose computed force fields were adjusted in variational calculations for best agreement with experiment (H 2 CO [88], C 2 H 2 O [89], NH 3 [90] and PH 3 [91]). As already pointed out, CV contributions have not been considered.
As before, results are reported in terms of MAX, MUE and RMSD. From Table 3, it is quite apparent that the different basis sets tested in conjunction with the CCSD(T)-F12 ansatz provide very similar results. Therefore, the jun-cc-pVTZ basis set, which has been found to be the optimal choice for energy and geometry evaluations, can be confidently employed also for harmonic frequencies. Furthermore, the much cheaper rDSD level performs remarkably well and can be profitably used for the computation of harmonic frequencies of larger molecules and thus for the accurate evaluation of their ZPEs.

2.3. Amino Acids

The structure of isolated amino acids is ruled by both backbone ( ϕ = HNC α C’, ψ = NC α C’O(H) and ω = C α C’OH) and side-chain ( χ ) torsional angles (see Figure 1 for the specific case of serine). The non-planarity of the NH 2 moiety suggests that, instead of the customary ϕ dihedral angle, ϕ = LP-N-C α -C’ = ϕ + 120 can be better employed. In the definition above, LP is the nitrogen lone-pair perpendicular to the plane defined by the two aminic hydrogens and the C α atom.
The most stable backbone structures result from the formation of hydrogen bonds, which can be classified as follows, also according to the values of the dihedral angles defined above. Type I is a bifurcated NH 2 O=C hydrogen bond ( ϕ 180 , ψ 180 , ω 180 ), type II is a N⋯HO hydrogen bond ( ϕ 0 , ψ 0 , ω 0 ), and type III is a bifurcated NH 2 OH hydrogen bond ( ϕ 180 , ψ 0 , ω 180 ). The type of hydrogen bond in the backbone also leads to the principal denomination of the structure, as shown in Figure 2. It has to be noted that type III structures have never been observed in MW experiments of natural α -aminoacids without polar side-chains due to their easy relaxation to the corresponding type I counterparts. Additional hydrogen bonds can be established with polar side chains, and they further stabilize the structure. This is the case, for example, of the type III’ conformers (single NH⋯OH hydrogen bond, ϕ 180 , ψ 90 , ω 180 ). Typical conformers of types I, II and III’ are shown in Figure 2 for the case of serine (which will be discussed in detail in the following), where the g , g and t labels refer to the gauche and trans conformations of each χ dihedral angle (as evident from Figure 1, for serine, there are two χ angles).
While the strategy described in the Introduction is able to identify all possible conformers of a given amino acid, their experimental characterization is limited by the sensitivity of the technique employed. For the specific case of rotational spectroscopy, a conservative limit for the relative stability of detectable structures is around 900 cm 1 , which corresponds to a relative population of about 1% at room temperature (where kT/hc = 207 cm 1 ). In the same vein, an upper limit of 400 cm 1 for interconversion barriers is usually employed for discriminating fast relaxation processes in the case of flexible compounds showing several stable minima [4,39].
For the simplest amino acid, i.e., glycine, eight conformers have been characterized computationally, with only two of them (I and II) being sufficiently stable (and non-involved in fast relaxation processes) to be detected in MW experiments [92]. The limited size of glycine has allowed the evaluation of relative stabilities by means of state-of-the-art composite schemes including, together with CBS and CV contributions evaluated at the CCSD(T) level, full account of triple excitations, perturbative estimate of quadruple excitations, diagonal corrections to the Born–Oppenheimer approximation and relativistic effects [93]. While the reader is referred to [49,92,93,94] for details, here we point out the remarkable performance of the ChS and junChSF12 approaches in the evaluation of relative energies (Table 4), with MUEs smaller than 10 cm 1 with respect to the most accurate results [93]. Noticeably, the rDSD model performs well with a MUE of 15 cm 1 . It is worth noting that such an accuracy is by far better than that obtainable by the methods (B3LYP and MP2) usually employed in the interpretation of MW spectra of biomolecule building blocks (see Table 4). Similar remarks apply to spectroscopic parameters. In fact, the MUEs with respect to experimental rotational constants are 23.6 MHz for rDSD, 13.0 MHz for rDSD-LRA, and 17.1 MHz for ChS. These average deviations from experiment point out the great accuracy that can be reached by applying LRA corrections to rDSD structures. For glycine, the availability of an accurate SE equilibrium structure for conformer I permits to evaluate directly the accuracy of geometrical parameters for the rDSD and rDSD-LRA levels. For rDSD, the MUEs are 0.0019 Å for bond lengths and 0.0026 radians (0.15 degrees) for valence angles, with the former value decreasing to 0.0004 Å at the rDSD-LRA level. Moving to nuclear quadrupole coupling constants, we note the good agreement (nearly quantitative) with experiment obtained by the ChS and rDSD levels, with the latter requiring a greatly reduced computational effort with respect to the former. This outcome further confirms the good performance of the rDSD model.
Systematic investigations have revealed that, analogous with glycine, the natural amino acids containing simple non-polar side chains (alanine [29], valine [30], isoleucine [31] and leucine [32]) have two dominant conformers of type I and II. On the contrary, the conformational landscape of natural amino acids with polar side chains is much richer due to the cooperativity or competition between intra-backbone and backbone with side-chain hydrogen bonds. The simplest amino acid containing a polar side-chain is serine, which has two degrees of freedom in its CH 2 OH side-chain: χ 1 = N-C α -C β -O and χ 2 = C α C β -O-H (see Figure 1). The exploration of the conformational PES and refinement of low-lying energy minima end up with 12 structures within 900 cm 1 above the absolute energy minimum. However, five of them relax to more stable counterparts through low energy barriers. Thus, seven conformers are left as possibly detectable in MW experiments: three of type II, two of type III’, and one each for types I and I’. Contrary to the case of glycine, the most stable conformer of serine in terms of electronic energies is of type II. However, incorporation of ZPE and thermal effects reverses the situation. As a result, the most stable conformer is Ig g in terms of standard free energies.
The rotational constants of the two most stable conformers have been recently computed by geometry optimizations at the ChS level, with relative MAX and MUE of 0.6% and 0.3% with respect to the experiment, respectively [95]. It is noteworthy that very similar (but slightly improved) relative errors (0.4% and 0.3%, respectively) are obtained at the rDSD-LRA level (see Table 4). The strongly reduced computational requirements of this level have allowed us to obtain accurate results also for the other five low-energy conformers experimentally characterized, the results being provided in Table 5. It is noted that the comparison with the experiment confirm (i) an average relative deviation smaller than 0.4% for rotational constants and (ii) the good accuracy that rDSD is able to reach in the prediction of quadrupole coupling constants. This table also points out that the first four conformers (one of type I, one of type I’ and two of type II) are significantly more stable than the remaining ones, with the Ig g conformer becoming the most stable one once moving from relative electronic energies to standard free energies. At the rDSD/B3 level (rDSD referring to the level of theory used for electronic energies, B3 to the level employed for the computation of ZPE and thermal corrections), the IIgg conformer lies 44.1 cm 1 above (38.7 cm 1 employing ChS electronic energies).

3. Conclusions

In this paper, we have focused on the final steps of a general strategy that aims at the accurate structural and spectroscopic characterization of flexible molecules. The main outcome of this work is that accurate geometries, spectroscopic parameters and relative energies can be obtained by the double-hybrid rev-DSD-PBEP86 functional. In particular, when combined with the jun-cc-pVTZ basis set, this DFT model is able to provide results approaching the accuracy of those issued by the well-tested ChS approaches that exploit the gold-standard CCSD(T)(-F12) method. The performance of the ChS composite schemes for geometrical parameters (and thus equilibrium rotational constants) is matched (if not improved) by rDSD when this is coupled with the linear regression approach (LRA), which does not involve any additional computational cost.
In detail, we have demonstrated that rotational and quadrupole coupling constants—parameters that strongly depend on mass distribution and intra-molecular interactions—can be obtained at the rDSD level with an accuracy which is suitable for guiding and supporting experiments in the field of rotational spectroscopy. The results obtained for glycine and serine are in very good agreement with the available spectroscopic data and pave the way toward the accurate investigation of large biomolecule building blocks.

Author Contributions

V.B. and C.P. have equally contributed to the work. S.D.G. has produced the results reported in Table 2 and Table 3. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by MUR (PRIN Grant No. 202082CE3T), by the University of Bologna (RFO funds), and by the Italian Space Agency (ASI; ‘Life in Space’ project, No. 2019-3-U.0). The COST Action CA21101 “COSY-Confined molecular systems: from a new generation of materials to the stars” is also acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Marco Fusé (University of Brescia) and Federico Lazzari (Scuola Normale Superiore) are acknowledged for discussions and technical help. The SMART@SNS Laboratory (http://smart.sns.it (accessed on 19 December 2022)) is acknowledged for providing high-performance computing facilities.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Schols, G. Atomic and Molecular Beam Methods; Oxford University Press: Oxford, UK, 1988. [Google Scholar]
  2. Alonso, J.L.; López, J.C. Microwave spectroscopy of biomolecular building blocks. In Gas-Phase IR Spectroscopy and Structure of Biological Molecules; Springer: Berlin/Heidelberg, Germany, 2015; pp. 335–401. [Google Scholar]
  3. Godfrey, P.D.; Rodgers, F.M.; Brown, R.D. Theory versus Experiment in Jet Spectroscopy: Glycolic Acid. J. Am. Chem. Soc. 1997, 119, 2232–2239. [Google Scholar] [CrossRef]
  4. Florio, G.M.; Christie, R.A.; Jordan, K.D.; Zwier, T.S. Conformational Preferences of Jet-Cooled Melatonin: Probing trans- and cis-Amide Regions of the Potential Energy Surface. J. Am. Chem. Soc. 2002, 124, 10236–10247. [Google Scholar] [CrossRef]
  5. Helgaker, T.; Klopper, W.; Tew, D.P. Quantitative quantum chemistry. Mol. Phys. 2008, 106, 2107–2143. [Google Scholar] [CrossRef]
  6. Karton, A. A computational chemist’s guide to accurate thermochemistry for organic molecules. WIRES Comp. Mol. Sci. 2016, 6, 292–310. [Google Scholar] [CrossRef]
  7. Puzzarini, C.; Bloino, J.; Tasinato, N.; Barone, V. Accuracy and interpretability: The devil and the holy grail. New routes across old boundaries in computational spectroscopy. Chem. Rev. 2019, 119, 8131–8191. [Google Scholar] [CrossRef] [PubMed]
  8. Kesharwani, M.K.; Karton, A.; Martin, J.M. Benchmark ab initio conformational energies for the proteinogenic amino acids through explicitly correlated methods. Assessment of density functionale methods. J. Chem. Theory Comput. 2016, 12, 444–454. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, P.; Shu, C.; Ye, H.; Biczysko, M. Structural and energetic properties of amino acids and peptides benchmarked by accurate theoretical and experimental data. J. Phys. Chem. A 2021, 125, 9826–9837. [Google Scholar] [CrossRef]
  10. Mancini, G.; Fusè, M.; Lazzari, F.; Chandramouli, B.; Barone, V. Unsupervised search of low-lying conformers with spectroscopic accuracy: A two-step algorithm rooted into the island model evolutionary algorithm. J. Chem. Phys. 2020, 153, 124110. [Google Scholar] [CrossRef]
  11. Ferro-Costas, D.; Mosquera-Lois, I.; Fernandez-Ramos, A. Torsiflex: An automatic generator of torsional conformers. application to the twenty proteinogenic amino acids. J. Cheminf. 2021, 13, 100. [Google Scholar] [CrossRef]
  12. Barone, V.; Puzzarini, C.; Mancini, G. Integration of theory, simulation, artificial intelligence and virtual reality: A four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy. Phys. Chem. Chem. Phys. 2021, 23, 17079–17096. [Google Scholar] [CrossRef]
  13. León, I.; Fusè, M.; Alonso, E.R.; Mata, S.; Mancini, G.; Puzzarini, C.; Alonso, J.L.; Barone, V. Unbiased disentanglement of conformational baths with the help of microwave spectroscopy, quantum chemistry and artificial intelligence: The puzzling case of homocysteine. J. Chem. Phys. 2022, 157, 074107. [Google Scholar] [CrossRef]
  14. Mancini, G.; Fusè, M.; Lazzari, F.; Barone, V. Fast exploration of potential energy surfaces with a joint venture of quantum chemistry, evolutionary algorithms and unsupervised learning. Digit. Discov. 2022, 1, 10539–10547. [Google Scholar] [CrossRef]
  15. Lazzari, F.; Salvadori, A.; Mancini, G.; Barone, V. Molecular Perception for Visualization and Computation: The Proxima Library. J. Chem. Inf. Model. 2020, 60, 2668–2672. [Google Scholar] [CrossRef] [PubMed]
  16. Bannwarth, C.; Ehlert, S.; Grimme, S. GFN2-xTB—an accurate and broadly parametrized self-consistent tight-binding quantum chemical method with multipole electrostatics and density-dependent dispersion contributions. J. Chem. Theory Comput. 2019, 15, 1652–1671. [Google Scholar] [CrossRef] [Green Version]
  17. Puzzarini, C.; Biczysko, M.; Barone, V.; Largo, L.; Pena, I.; Cabezas, C.; Alonso, J.L. Accurate characterization of the peptide linkage in the gas phase: A joint quantum-chemistry and rotational spectroscopy study of the glycine dipeptide analogue. J. Phys. Chem. Lett. 2014, 5, 534–540. [Google Scholar] [CrossRef] [PubMed]
  18. Alessandrini, S.; Barone, V.; Puzzarini, C. Extension of the “Cheap” Composite Approach to Noncovalent Interactions: The jun-ChS Scheme. J. Chem. Theory Comput. 2020, 16, 988–1006. [Google Scholar] [CrossRef] [PubMed]
  19. Lupi, J.; Alessandrini, S.; Barone, V.; Puzzarini, C. junChS and junChS-F12 Models: Parameter-free Efficient yet Accurate Composite Schemes for Energies and Structures of Noncovalent Complexes. J. Chem. Theory Comput. 2021, 17, 6974–6992. [Google Scholar] [CrossRef] [PubMed]
  20. Barone, V. Anharmonic vibrational properties by a fully automated second order perturbative approach. J. Chem. Phys. 2005, 122, 014108. [Google Scholar] [CrossRef]
  21. Rosnik, A.M.; Polik, W.F. VPT2+K spectroscopic constants and matrix elements of the transformed vibrational Hamiltonian of a polyatomic molecule with resonances using Van Vleck perturbation theory. Mol. Phys. 2014, 112, 261–300. [Google Scholar] [CrossRef]
  22. Franke, P.R.; Stanton, J.F.; Douberly, G.E. How to VPT2: Accurate and Intuitive Simulations of CH Stretching Infrared Spectra Using VPT2+ K with Large Effective Hamiltonian Resonance Treatments. J. Phys. Chem. A 2021, 125, 1301–1324. [Google Scholar] [CrossRef]
  23. Mendolicchio, M.; Bloino, J.; Barone, V. Perturb-then-diagonalize vibrational engine exploiting curvilinear internal coordinates. J. Chem. Theory Comput. 2022, 18, 7603–7619. [Google Scholar] [CrossRef] [PubMed]
  24. Grimme, S. Supramolecular Binding Thermodynamics by Dispersion-Corrected Density Functional Theory. Chem. A Eur. J. 2012, 18, 9955–9964. [Google Scholar] [CrossRef] [PubMed]
  25. Li, X.; Spada, L.; Alessandrini, S.; Zheng, Y.; Lengsfels, K.G.; Grabow, J.U.; Feng, G.; Puzzarini, C.; Barone, V. Stacked but not Stuck: Unveiling the role of ππ* Interactions with the Help of the Benzofuran-Formaldehyde Complex. Angew. Chem. Int. Ed. Engl. 2022, 61, 264–270. [Google Scholar] [CrossRef] [PubMed]
  26. Alonso, E.R.; Fusé, M.; León, I.; Puzzarini, C.; Alonso, J.L.; Barone, V. Exploring the Maze of Cycloserine Conformers in the Gas Phase Guided by Microwave Spectroscopy and Quantum Chemistry. J. Phys. Chem. A 2021, 125, 2121–2129. [Google Scholar] [CrossRef] [PubMed]
  27. Sanz, M.E.; Blanco, S.; López, J.C.; Alonso, J.L. Rotational Probes of Six Conformers of Neutral Cysteine. Angew. Chem. Int. Ed. 2008, 47, 6216–6220. [Google Scholar] [CrossRef] [PubMed]
  28. Lovas, F.J.; Kawashima, Y.; Grabow, J.U.; Suenram, R.D.; Fraser, G.T.; Hirota, E. Microwave Spectra, Hyperfine Structure, and Electric Dipole Moments for Conformers I and II of Glycine. Astrophys. J. 1995, 455, L201–L204. [Google Scholar] [CrossRef]
  29. Blanco, S.; Lesarri, A.; López, J.C.; Alonso, J.L. The gas-phase structure of alanine. J. Am. Chem. Soc. 2004, 126, 11675–11683. [Google Scholar] [CrossRef]
  30. Lesarri, E.J.; Cocinero, J.C.; López, J.C.; Alonso, J.L. The Shape of Neutral Valine. Angew. Chem. Int. Ed. Engl. 2004, 43, 605–610. [Google Scholar] [CrossRef]
  31. Lesarri, E.J.; Sanchez, R.; Cocinero, E.J.; López, J.C.; Alonso, J.L. Coded Amino Acids in Gas Phase: The Shape of Isoleucine. J. Am. Chem. Soc. 2005, 127, 12952–12956. [Google Scholar] [CrossRef]
  32. Cocinero, J.C.; Lesarri, E.J.; Grabow, J.U.; López, J.C.; Alonso, J.L. The Shape of Leucine in the Gas Phase. Chem. Phys. Chem. 2007, 8, 599–604. [Google Scholar] [CrossRef]
  33. Blanco, S.; Sanz, M.E.; López, J.C.; Alonso, J.L. Revealing the multiple structures of serine. Proc. Natl. Acad. Sci. USA 2007, 104, 20183–20188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Alonso, J.L.; Peña, I.; López, J.C.; Vaquero, V. Rotational spectra signatures of four tautomers of guanine. Angew. Chem. Int. Ed. 2009, 48, 6141–6143. [Google Scholar] [CrossRef] [PubMed]
  35. Cabezas, C.; Varela, M.; Peña, I.; Mata, S.; López, J.C.; Alonso, J.L. The conformational locking of asparagine. Chem. Commun. 2012, 48, 5934–5936. [Google Scholar] [CrossRef] [PubMed]
  36. Sanz, M.E.; Cabezas, C.; Mata, S.; Alonso, J.L. Rotational spectrum of tryptophan. J. Chem. Phys. 2014, 140, 204318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Perez, C.; Mata, S.; Blanca, S.; Lopez, J.C.; Alonso, J.L. Jet-cooled rotational spectrum of laser-ablated phenylalanine. J. Phys. Chem. A 2011, 115, 9653–9657. [Google Scholar] [CrossRef] [PubMed]
  38. Perez, C.; Mata, S.; Cabezas, C.; Lopez, J.C.; Alonso, J.L. The rotational spectrum of tyrosine. J. Phys. Chem. A 2015, 119, 3731–3735. [Google Scholar] [CrossRef] [Green Version]
  39. Godfrey, P.D.; Brown, R.D. Proportions of Species Observed in Jet Spectroscopy-Vibrational Energy Effects: Histamine Tautomers and Conformers. J. Am. Chem. Soc. 1998, 120, 10724–10732. [Google Scholar] [CrossRef]
  40. Alonso, E.R.; León, I.; Alonso, J.L. Intra- and Intermolecular Interactions Between Non-Covalently Bonded Species; Elsevier: Amsterdam, The Netherlands, 2020; pp. 93–141. [Google Scholar]
  41. Nguyen, H.V.L.; Kleiner, I. Understanding (coupled) large amplitude motions: The interplay of microwave spectroscopy, spectral modeling, and quantum chemistry. Phys. Sci. Rev. 2020, 20200037. [Google Scholar] [CrossRef]
  42. Santra, G.; Sylvetsky, N.; Martin, J.M. Minimally empirical double-hybrid functionals trained against the GMTKN55 database: RevDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4. J. Phys. Chem. A 2019, 123, 5129–5143. [Google Scholar] [CrossRef]
  43. Papajak, E.; Zheng, J.; Xu, X.; Leverentz, H.R.; Truhlar, D.G. Perspectives on Basis Sets Beautiful: Seasonal Plantings of Diffuse Basis Functions. J. Chem. Theory Comput. 2011, 7, 3027–3034. [Google Scholar] [CrossRef]
  44. Lupi, J.; Puzzarini, C.; Cavallotti, C.; Barone, V. State-of-the-Art Quantum Chemistry Meets Variable Reaction Coordinate Transition State Theory to Solve the Puzzling Case of the H2S + Cl System. J. Chem. Theory Comput. 2020, 16, 5090–5104. [Google Scholar] [CrossRef] [PubMed]
  45. García de la Concepción, J.; Puzzarini, C.; Barone, V.; Jiménez-Serra, I.; Roncero, O. Formation of phosphorus monoxide (PO) in the interstellar medium: Insights from quantum-chemical and kinetic calculations. Astrophys. J. Lett. 2021, 922, 169. [Google Scholar] [CrossRef]
  46. Barone, V.; Lupi, J.; Salta, Z.; Tasinato, N. Development and Validation of a Parameter-Free Model Chemistry for the Computation of Reliable Reaction Rates. J. Chem. Theory Comput. 2021, 17, 4913–4928. [Google Scholar] [CrossRef] [PubMed]
  47. Shavitt, I.; Bartlett, R.J. Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  48. Tajti, A.; Szalay, P.G.; Császár, A.G.; Kállay, M.; Gauss, J.; Valeev, E.F.; Flowers, B.A.; Vázquez, J.; Stanton, J.F. HEAT: High accuracy extrapolated ab initio thermochemistry. J. Chem. Phys. 2004, 121, 11599–11613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Barone, V.; Biczysko, M.; Bloino, J.; Puzzarini, C. Accurate structure, thermochemistry and spectroscopic parameters from CC and CC/DFT schemes: The challenge of the conformational equilibrium in glycine. Phys. Chem. Chem. Phys. 2013, 15, 10094–10111. [Google Scholar] [CrossRef] [Green Version]
  50. Puzzarini, C.; Barone, V. The challenging playground of astrochemistry: An integrated rotational spectroscopy—Quantum chemistry strategy. Phys. Chem. Chem. Phys. 2020, 22, 6507–6523. [Google Scholar] [CrossRef]
  51. Puzzarini, C.; Barone, V. Extending the molecular size in accurate quantum-chemical calculations: The equilibrium structure and spectroscopic properties of uracil. Phys. Chem. Chem. Phys. 2011, 13, 7189–7197. [Google Scholar] [CrossRef] [PubMed]
  52. Møller, C.; Plesset, M.S. Note on an Approximation Treatment for Many-Electron Systems. Phys. Rev. 1934, 46, 618–622. [Google Scholar] [CrossRef] [Green Version]
  53. Raghavachari, K.; Trucks, G.W.; Pople, J.A.; Head-Gordon, M. A fifth-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 1989, 157, 479–483. [Google Scholar] [CrossRef]
  54. Helgaker, T.; Klopper, W.; Koch, H.; Noga, J. Basis-set convergence of correlated calculations on water. J. Chem. Phys. 1997, 106, 9639–9646. [Google Scholar] [CrossRef]
  55. Woon, D.E.; Dunning, T.H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. V. Core-Valence Basis Sets for Boron through Neon. J. Chem. Phys. 1995, 103, 4572. [Google Scholar] [CrossRef] [Green Version]
  56. Peterson, K.A.; Dunning, T.H., Jr. Accurate correlation consistent basis sets for molecular core-valence correlation effects: The second row atoms Al-Ar, and the first row atoms B-Ne revisited. J. Chem. Phys. 2002, 117, 10548–10560. [Google Scholar] [CrossRef] [Green Version]
  57. Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. JChPh 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  58. Dunning, T.H.; Peterson, K.A.; Wilson, A.K. Gaussian basis sets for use in correlated molecular calculations. X. The atoms aluminum through argon revisited. J. Chem. Phys. 2001, 114, 9244–9253. [Google Scholar] [CrossRef] [Green Version]
  59. Knizia, G.; Adler, T.B.; Werner, H.J. Simplified CCSD(T)-F12 methods: Theory and benchmarks. J. Chem. Phys. 2009, 130, 054104. [Google Scholar] [CrossRef] [PubMed]
  60. Werner, H.-J.; Adler, T.B.; Manby, F.R. General orbital invariant MP2-F12 theory. J. Chem. Phys. 2007, 126, 164102. [Google Scholar] [CrossRef]
  61. Feller, D.; Peterson, K. An expanded calibration study of the explicitly correlated CCSD(T)-F12b method using large basis set standard CCSD(T) atomization energies. J. Chem. Phys. 2013, 139, 084110. [Google Scholar] [CrossRef]
  62. Puzzarini, C.; Biczysko, M.; Barone, V.; Peña, I.; Cabezas, C.; Alonso, J.L. Accurate Molecular Structure and Spectroscopic Properties of Nucleobases: A Combined Computational-Microwave Investigation of 2-Thiouracil as a Case Study. Phys. Chem. Chem. Phys. 2013, 15, 16965–16975. [Google Scholar] [CrossRef] [Green Version]
  63. Melli, A.; Melosso, M.; Lengsfeld, K.G.; Bizzocchi, L.; Rivilla, V.M.; Dore, L.; Barone, V.; Grabow, J.U.; Puzzarini, C. Spectroscopic Characterization of 3-Aminoisoxazole, a Prebiotic Precursor of Ribonucleotides. Molecules 2022, 27, 3278. [Google Scholar] [CrossRef]
  64. Barone, V.; Ceselin, G.; Fusè, M.; Tasinato, N. Accuracy meets interpretability for computational spectroscopy by means of hybrid and double-hybrid functionals. Front. Chem. 2020, 8, 584203. [Google Scholar] [CrossRef]
  65. Biczysko, M.; Panek, P.; Scalmani, G.; Bloino, J.; Barone, V. Harmonic and Anharmonic Vibrational Frequency Calculations with the Double-Hybrid B2PLYP Method: Analytic Second Derivatives and Benchmark Studies. J. Chem. Theory Comput. 2010, 6, 2115–2125. [Google Scholar] [CrossRef] [PubMed]
  66. Martin, J.M.L.; Kesharwani, M.K. Assessment of CCSD(T)-F12 approximations and basis sets for harmonic vibrational frequencies. J. Chem. Theory Comput. 2014, 10, 2085–2090. [Google Scholar] [CrossRef] [PubMed]
  67. Barone, V.; Biczysko, M.; Bloino, J. Fully anharmonic IR and Raman spectra of medium-size molecular systems: Accuracy and interpretation. Phys. Chem. Chem. Phys. 2014, 16, 1759–1787. [Google Scholar] [CrossRef]
  68. Luchini, G.; Alegre-Requena, J.; Funes-Ardoiz, I.; Paton, R. GoodVibes: Automated thermochemistry for heterogeneous computational chemistry data [version 1; peer review: 2 approved with reservations]. F1000Research 2020, 9, 291. [Google Scholar] [CrossRef]
  69. Puzzarini, C.; Stanton, J.F. Connections between the accuracy of rotational constants and equilibrium molecular structures. Phys. Chem. Chem. Phys. 2023, 25. [Google Scholar] [CrossRef]
  70. Puzzarini, C.; Heckert, M.; Gauss, J. The accuracy of rotational constants predicted by high-level quantum-chemical calculations. I. Molecules containing first-row atoms. J. Chem. Phys. 2008, 128, 194108. [Google Scholar] [CrossRef] [PubMed]
  71. Piccardo, M.; Penocchio, E.; Puzzarini, C.; Biczysko, M.; Barone, V. Semi-experimental equilibrium structure determinations by employing B3LYP/SNSD anharmonic force fields: Validation and application to semirigid organic molecules. J. Phys. Chem. A 2015, 119, 2058–2082. [Google Scholar] [CrossRef]
  72. Penocchio, E.; Piccardo, M.; Barone, V. Semiexperimental Equilibrium Structures for Building Blocks of Organic and Biological Molecules: The B2PLYP Route. J. Chem. Theory Comput. 2015, 11, 4689–4707. [Google Scholar] [CrossRef]
  73. Ceselin, G.; Barone, V.; Tasinato, N. Accurate biomolecular structures by the nano-LEGO approach: Pick the bricks and build your geometry. J. Chem Theory Comput. 2021, 17, 7290–7311. [Google Scholar] [CrossRef]
  74. Melli, A.; Tonolo, F.; Barone, V.; Puzzarini, C. Extending the Applicability of the Semi-experimental Approach by Means of “Template Molecule” and “Linear Regression” Models on Top of DFT Computations. J. Phys. Chem. A 2021, 125, 9904–9916. [Google Scholar] [CrossRef]
  75. Puzzarini, C.; Stanton, J.F.; Gauss, J. Quantum-chemical calculation of spectroscopic parameters for rotational spectroscopy. Int. Rev. Phys. Chem. 2010, 29, 273–367. [Google Scholar] [CrossRef]
  76. Gordy, W.; Cook, R.L. Microwave Molecular Spectra; Wiley: Hoboken, NJ, USA, 1984. [Google Scholar]
  77. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16 Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  78. Stanton, J.F.; Gauss, J.; Harding, M.E.; Szalay, P.G. CFOUR. A Quantum Chemical Program Package. 2016. Available online: http://www.cfour.de (accessed on 19 December 2022).
  79. Werner, H.J.; Knowles, P.J.; Manby, F.R.; Black, J.A.; Doll, K.; Hesselmann, A.; Kats, D.; Köhn, A.; Korona, T.; Kreplin, D.A.; et al. The Molpro quantum chemistry package. J. Chem. Phys. 2020, 152, 144107. [Google Scholar] [CrossRef] [Green Version]
  80. Puzzarini, C.; Spada, L.; Alessandrini, S.; Barone, V. The challenge of non-covalent interactions: Theory meets experiment for reconciling accuracy and interpretation. J. Phys. Condens. Matter 2020, 32, 343002. [Google Scholar] [CrossRef]
  81. Řezáč, J.; Hobza, P. Describing Noncovalent Interactions beyond the Common Approximations: How Accurate Is the “Gold Standard,” CCSD(T) at the Complete Basis Set Limit? J. Chem. Theory Comput. 2013, 9, 2151–2155. [Google Scholar] [CrossRef] [PubMed]
  82. Řezáč, J.; Dubecký, M.; Jurečka, P.; Hobza, P. Extensions and applications of the A24 data set of accurate interaction energies. Phys. Chem. Chem. Phys. 2015, 17, 19268–19277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Barone, V.; Puzzarini, C. Gas-Phase Computational Spectroscopy: The Challenge of the Molecular Bricks of Life. Ann. Rev. Phys. Chem. 2023, 74. [Google Scholar] [CrossRef] [PubMed]
  84. Benedict, W.S.; Gailar, N.; Plyler, E.K. Rotation-Vibration Spectra of Deuterated Water Vapor. J. Chem. Phys. 1956, 24, 1139. [Google Scholar] [CrossRef]
  85. Smith, A.M.; Coy, S.L.; Klemperer, W.; Lehmann, K.K. Fourier Transform Spectra of Overtone Bands of HCN from 5400 to 15100 cm−1. J. Mol. Spectrosc. 1989, 134, 134–153. [Google Scholar] [CrossRef]
  86. Teffo, J.L.; Sulakshina, O.N.; Perevalov, V.I. Effective Hamiltonian for Rovibrational Energies and Line Intensities of Carbon Dioxide]. J. Mol. Spectrosc. 1992, 156, 48–64. [Google Scholar] [CrossRef]
  87. Strey, G.; Mills, I.M. Anharmonic Force Field of Acetylene]. J. Mol. Spectrosc. 1976, 59, 103–115. [Google Scholar] [CrossRef]
  88. Yachmenev, A.; Yurchenko, S.N.; Jensen, P.; Thiel, W. A New Spectroscopic Potential Energy Surface for Formaldehyde in its Ground Electronic State. J. Chem. Phys. 2011, 134, 244307. [Google Scholar] [CrossRef] [PubMed]
  89. East, A.L.L.; Allen, W.D.; Klippenstein, S.J. [The Anharmonic Force Field and Equilibrium Molecular Structure of Ketene]. J. Chem. Phys. 1995, 102, 8506–8532. [Google Scholar] [CrossRef]
  90. Huang, X.; Schwenke, D.W.; Lee, T.J. Rovibrational Spectra of Ammonia. I. Unprecedented Accuracy of a Potential Energy Surface Used with Nonadiabatic Corrections. J. Chem. Phys. 2011, 134, 044320. [Google Scholar] [CrossRef]
  91. Wang, D.; Shi, Q.; Zhu, Q.S. An ab Initio Quartic Force Field of PH3. J. Chem. Phys. 2000, 112, 9624. [Google Scholar] [CrossRef]
  92. Shu, C.; Jiang, Z.; Biczysko, M. Toward accurate prediction of amino acid derivatives structure and energetics from DFT: Glycine conformers and their interconversions. J. Mol. Model. 2020, 26, 129. [Google Scholar] [CrossRef]
  93. Nacsa, A.B.; Csako, G. Benchmark ab initio proton affinity of glycine. Phys. Chem. Chem. Phys. 2021, 23, 9663–9671. [Google Scholar] [CrossRef] [PubMed]
  94. Puzzarini, C.; Barone, V. Diving for accurate structures in the ocean of molecular systems with the help of spectroscopy and quantum chemistry. Acc. Chem. Res. 2018, 51, 548–556. [Google Scholar] [CrossRef] [Green Version]
  95. Sheng, M.; Silvestrini, F.; Biczysko, M.; Puzzarini, C. Structural and vibrational properties of amino acids from composite schemes and double-hybrid DFT: Hydrogen bonding in serine as a test case. J. Phys. Chem. A 2021, 125, 9099–9114. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Structure and main dihedral angles of serine. See the main text for further details.
Figure 1. Structure and main dihedral angles of serine. See the main text for further details.
Molecules 28 00913 g001
Figure 2. Representations of low-lying serine conformers of type I, II and III’. H-bonds are highlighted by dashed lines.
Figure 2. Representations of low-lying serine conformers of type I, II and III’. H-bonds are highlighted by dashed lines.
Molecules 28 00913 g002
Table 1. A14 database: junChS, junChSF12 and junCCF12+CV CP-corrected interaction energies a [19]. All the values are in kJ mol−1.
Table 1. A14 database: junChS, junChSF12 and junCCF12+CV CP-corrected interaction energies a [19]. All the values are in kJ mol−1.
junChSjunChSF12junCCF12+CVBest b
H2O···H2O−21.10−21.00−21.03−21.07
NH3···NH3−13.30−13.20−12.99−13.26
HF···HF−19.45−19.25−19.25−19.18
CH2O···CH2O−19.23−18.96−18.59−18.89
HCN···HCN−19.88−19.84−19.83−19.95
C2H4···C2H4−4.75−4.66−4.36−4.64
CH4···CH4−2.25−2.20−1.99−2.27
H2O···NH3−27.57−27.47−27.34−27.39
H2O···C2H4−10.86−10.77−10.55−10.82
C2H4···CH2O−6.94−6.83−6.56−6.84
NH3···C2H4−5.89−5.82−5.61−5.84
HF···CH4−7.13−7.04−6.85−6.95
H2O···CH4−2.78−2.77−2.68−2.85
NH3···CH4−3.24−3.23−3.11−3.26
MAX  c 0.340.110.30
MUE  c 0.110.060.19
RMSD  c 0.140.070.21
a Energies evaluated at the reference geometries from [81]. b Errors are with respect to the best values from [82]. c MAX, MUE and RMSD stand for maximum error, mean unsigned error, and root-mean-square deviation, respectively.
Table 2. Statistical analysis of the geometrical parameters, for a subset of the SE100 database, obtained at different levels of theory.
Table 2. Statistical analysis of the geometrical parameters, for a subset of the SE100 database, obtained at different levels of theory.
junChSDZCCF12DZCCF12+CVTZCCF12+CVjunCCF12+CV
2nd-row (17 molecules)
MAX(r)0.00420.00440.00270.00390.0044
MAX( θ , ϕ )1.520.610.620.790.62
MUE(r)0.00110.00220.00090.00080.0007
MUE( θ , ϕ )0.190.190.160.160.15
3rd-row (7 molecules)
MAX(r)0.00120.00460.00400.00210.0014
MAX( θ , ϕ )0.100.210.240.090.17
MUE(r)0.00050.00180.00100.00040.0005
MUE( θ , ϕ )0.050.100.090.020.06
Table 3. Harmonic frequencies (cm 1 ) for the molecules of the H8 set.
Table 3. Harmonic frequencies (cm 1 ) for the molecules of the H8 set.
DZCCF12junCCF12augCCF12rDSD
MAX  a 13.712.818.418.2
MUE  a 4.14.04.65.6
RMSD  a 5.05.26.17.1
a Statistical measures have been determined with respect to the best available values. See text for details.
Table 4. Equilibrium rotational constants and 14N-nuclear quadrupole coupling constants (MHz) of the two most stable conformers of glycine and serine. The relative energy differences ( Δ E, in cm 1 ) are also reported.
Table 4. Equilibrium rotational constants and 14N-nuclear quadrupole coupling constants (MHz) of the two most stable conformers of glycine and serine. The relative energy differences ( Δ E, in cm 1 ) are also reported.
ParameterExp. a ChS b rDSD b rDSD-LRA b MP2/cc-pVTZ b B3LYP/SNSD b
Glycine (I)A e 10,418.210,396.610,334.810,390.310,328.010,283.1
B e 3906.93901.13879.93897.63905.03831.1
C e 2934.42930.42913.52927.42926.22882.9
χ a a −1.208(9)−1.278−1.336
χ b b −0.343(8)−0.464−0.448
χ c c 1.552(10)1.7421.785
Glycine (II)A e 10,144.510,205.310,139.310,193.810,178.710,135.0
B e 4094.54095.64078.74097.24104.74043.4
C e 3024.73030.63021.33035.83041.02993.1
χ a a 1.773(2)1.8761.922
χ b b −3.194(4)−3.286−3.344
χ c c 1.421(4)1.4131.422
Δ E 201.5214.8 157.2237.8
Serine (Ig g)A e 4528.14499.44487.04510.44494.94485.2
B e 1838.81841.31822.71831.61830.81809.1
C e 1460.91460.01451.81459.01462.71433.3
χ a a −4.302(3)−4.416−4.554
χ b b 2.8236(6)2.8522.868
χ c c 1.479(5)1.5651.685
Serine (IIgg)A e 3585.93560.23559.33578.13554.83547.2
B e 2412.72410.22393.02404.72414.02342.0
C e 1754.41757.61739.71748.21757.91713.6
χ a a −3.462(2)−3.530−3.670
χ b b 2.0797(9)2.1492.134
χ c c 1.382(5)1.3811.536
Δ E −167.2−161.8 −185.5−100.8
a SE equilibrium rotational constants obtained by correcting the experimental ground-state rotational constants for computed vibrational corrections. Experimental data are from [28] for glycine and [33] for serine. Computed vibrational corrections are from [49] for glycine and [95] for serine. Experimental nuclear quadrupole coupling constants are from [28] for glycine and [33] for serine. b Glycine: ChS and MP2 results from [49]; rDSD and rDSD-LRA results from this work; B3LYP results from [71] and this work. Serine: ChS results from [95]; all the other data from this work.
Table 5. Ground-state rotational constants (A 0 , B 0 , and C 0 in MHz) and 14N-nuclear quadrupole coupling constants ( χ i i , i = a, b, c, in MHz) of the seven conformers of serine experimentally investigated. Computed relative standard free energies (in cm 1 ) are also reported.
Table 5. Ground-state rotational constants (A 0 , B 0 , and C 0 in MHz) and 14N-nuclear quadrupole coupling constants ( χ i i , i = a, b, c, in MHz) of the seven conformers of serine experimentally investigated. Computed relative standard free energies (in cm 1 ) are also reported.
Ig gIIggI’gg IItg III’ggIIg tIII’tg
Calc. a
A 0 4461.343549.333505.743630.863950.324508.133464.84
B 0 1823.012372.382305.212382.522222.911843.002304.68
C 0 1441.951734.671803.621515.281657.031462.051604.74
χ a a −4.5535−3.6696−0.9235−3.8114−0.6094−0.3660−1.0975
χ b b 2.86812.13412.55282.1268−0.67022.0569−0.6582
χ c c 1.68541.5355−1.62931.68471.2796−1.69091.7557
Δ G 0 0.044.1222.4295.6481.2522.7620.6
Exp. b
A 0 4479.0320(12)3557.20088(35)3524.38806(41)3638.05784(38)3931.7548(76)4517.473(17)3510.4015(35)
B 0 830.16170(25)2380.37208(40)2307.76826(70)2387.89651(99)2242.76701(70)1846.99360(30)2321.90829(24)
C 0 1443.79545(28)1740.92458(10)1805.20788(60)1519.18716(36)1664.53012(57)1463.79646(31)1584.38608(32)
χ a a −4.3023(27)−3.4616(19)−1.1343(35)−3.6257(57)−0.6733(67)−0.6066(55)−1.0486(55)
χ b b 2.82359(63)2.07974(93)2.5043(50)2.06213(26)−0.456(16)2.0723(82)−0.5637(53)
χ c c 1.4788(46)1.3819(47)−1.3701(50)1.5906(50)1.129(16)−1.466(30)1.612(21)
a All the computed data are at the rDSD level (rDSD-LRA for equilibrium rotational constants) except for vibrational corrections to equilibrium rotational constants, ZPE and thermal contributions to ΔG, which are all at the B3 level. b Experimental data are from [33]. Standard errors are shown in parentheses in units of the last digits.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barone, V.; Di Grande, S.; Puzzarini, C. Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks. Molecules 2023, 28, 913. https://doi.org/10.3390/molecules28020913

AMA Style

Barone V, Di Grande S, Puzzarini C. Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks. Molecules. 2023; 28(2):913. https://doi.org/10.3390/molecules28020913

Chicago/Turabian Style

Barone, Vincenzo, Silvia Di Grande, and Cristina Puzzarini. 2023. "Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks" Molecules 28, no. 2: 913. https://doi.org/10.3390/molecules28020913

APA Style

Barone, V., Di Grande, S., & Puzzarini, C. (2023). Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks. Molecules, 28(2), 913. https://doi.org/10.3390/molecules28020913

Article Metrics

Back to TopTop