Identification of Serine Conformers by Matrix-Isolation IR Spectroscopy Aided by Near-Infrared Laser-Induced Conformational Change, 2D Correlation Analysis, and Quantum Mechanical Anharmonic Computations

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

NIR and UV induced transformations of indazole-3-carboxylic acid isolated in low temperature matrices

The molecular structure and NIR and UV induced phototransformations of indazole-3-carboxylic acid were studied in low temperature argon and nitrogen matrices by FTIR spectroscopy and B3LYP/6-311++G(2d,2p) calculations. Eleven minima of IC were located on the S₀ potential energy surface. The three most stable structures were detected experimentally in both matrices after deposition. Upon NIR irradiation a dominant process was rotamerization within the carboxylic group leading to changes in the population of the trans 1HIC1 and cis 1HIC2 forms. In turn, at UV irradiation at λ = 260 nm two new tautomers (2HIC2 and 2HIC3) were generated indicating that the hydrogen atom transfer in pyrazole ring took place. Based on the obtained kinetic curves, differences in the phototransformation rates in different matrices were indicated.

Protonation of serine: conformers, proton affinities and gas-phase basicities at the "gold standard" and beyond

The potential energy surfaces (PESs) of serine and its protonated counterparts are investigated to determine the structures of the minima. A total of 95 neutral serine, 15 N-(amino-) and 46 O-(carbonyl-)protonated serine conformers are found. Their relative energies, geometries and harmonic vibrational frequencies are determined at the MP2/aug-cc-pVDZ level of theory. To obtain highly accurate thermodynamic values, further computations are performed: the ten conformers with the lowest relative energies from each molecule type (neutral, N- and O-protonated) are further optimized using the explicitly correlated CCSD(T)-F12a/cc-pVDZ-F12 method (for neutral serine, harmonic vibrational frequencies were also computed). In addition, auxiliary corrections were determined: basis-set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation and second-order Douglas-Kroll relativistic effects along with zero-point energy contributions. Two important thermodynamic parameters (at 298.15 K), proton affinity (PA)/gas-phase basicity (GB) are calculated considering the two different protonation sites: 218.05 ± 0.2/209.86 ± 0.6 kcal mol^-1 and 205.87 ± 0.2/196.36 ± 0.3 kcal mol^-1 for the amino and carbonyl sites, respectively. The uncertainty of the determined values was approximated based on various sources including auxiliary corrections, basis-set effects, harmonic vibrational frequencies.

Coupling and decoupling CH stretching vibration of methylene and methine in serine conformers

In the molecules of the early Earth, as a building block of proteins, serine has enormous chemical and biological significance. The vibrational spectroscopy of CH bonds plays an important role in probing biomolecules. Whether the CH stretching vibration bands can be accurately assigned will affect the accuracy of the detection results. In this study, we employed the MP2/cc-pVTZ method to calculate the Raman spectra of 85 serine conformers and the corresponding species with deuterium in the CH stretching region from 2800 cm^-1 to 3050 cm^-1 and then recorded the movement of each atom and the dihedral angles, CH bond lengths, and Raman shifts before and after deuterium for each conformer. We directly observed that the stretching vibration of two CH bonds in the methylene group decoupled to vibrate independently in some conformers, and the stretching vibrations of methylene and methine could be strongly coupled in some conformers. Those results are inconsistent with the traditional understanding, which is generally believed that the CH stretching vibrations are mutually coupled in a single methyl or methylene group to generate symmetric and antisymmetric stretching vibrations, while for different methyl, methylene or methine groups, the CH stretching vibrations cannot be mutually coupled. Through the statistical analysis between several factors, we found that the level of local coupling in serine methylene was correlated with the bond length difference between two CH bonds. Our work provides a new understanding of the vibrational modes of hydrocarbon bonds and the coupling between different hydrocarbon groups.

IR-induced and tunneling reactions in cryogenic matrices: the (incomplete) story of a successful endeavor

In this article, IR-induced and tunneling-driven reactions observed in cryogenic matrices are described in a historical perspective, the entangling of the two types of processes being highlighted. The story of this still ongoing fascinating scientific endeavor is presented here following closely our own involvement in the field for more than 30 years, and thus focuses mostly on our work. It is, because of this reason, also an incomplete story. Nevertheless, it considers a large range of examples, from very selective IR-induced conformational isomerizations to IR-induced bond-breaking/bond-forming reactions and successful observations of rare heavy atom tunneling processes. As a whole, this article provides a rather general overview of the major progress achieved in the field.

Sum-frequency vibrational spectroscopy of methanol at interfaces due to Fermi resonance

We present a theoretical method of studying sum-frequency vibrational spectroscopy for the CH3 group of methanol at interfaces due to Fermi resonance, which provides a novel and untraditional point of view with respect to traditional approaches.

Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data

Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)^MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol^-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH₂, Ac-Ala-NH₂, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)^MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Monomers of Glycine and Serine Have a Limited Ability to Hydrate in the Atmosphere

The role of atmospheric aerosols on climate change is one of the biggest uncertainties in most global climate models. Organic aerosols have been identified as potential cloud condensation nuclei (CCN), and amino acids are organic molecules that could serve as CCN. Amino acids make up a significant portion of the total organic material in the atmosphere, and herein we present a systematic study of hydration for two of the most common atmospheric amino acids, glycine and serine. We compute DLPNO/CCSD(T)//M08-HX/MG3S thermodynamic properties and atmospheric concentrations of Gly(H₂O)_n and Ser(H₂O)_n, where n = 1-5. We predict that serine-water clusters have higher concentrations at n = 1 and 5, while glycine-water clusters have higher concentrations at n = 2-4. However, both glycine and serine are inferred to exist primarily in their nonhydrated monomer forms in the absence of other species such as sulfuric acid.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

A Combined Systematic-Stochastic Algorithm for the Conformational Search in Flexible Acyclic Molecules

We propose an algorithm that is a combination of systematic variation of the torsions and Monte Carlo (or stochastic) search. It starts with a trial geometry in internal coordinates and with a set of preconditioned torsional angles, i.e., torsional angles at which minima are expected according to the chemical knowledge. Firstly, the optimization of those preconditioned geometries is carried out at a low electronic structure level, generating an initial set of conformers. Secondly, random points in the torsional space are generated outside the “area of influence” of the previously optimized minima (i.e., outside a hypercube about each minima). These random points are used to build the trial structure, which is optimized by an electronic structure software. The optimized structure may correspond to a new conformer (which would be stored) or to an already existing one. Initial torsional angles (and also final ones if a new conformer is found) are stored to prevent visiting the same region of the torsional space twice. The stochastic search can be repeated as many times as desired. Finally, the low-level geometries are recovered and used as the starting point for the high-level optimizations. The algorithm has been employed in the calculation of multi-structural quasi harmonic and multi-structural torsional anharmonic partition functions for a series of alcohols ranging from n-propanol to n-heptanol. It was also tested for the amino acid L-serine.

Conformational isomerizations triggered by vibrational excitation of second stretching overtones

Vibrational excitation using frequency-tunable IR laser light has been developed as a powerful tool for selective manipulation of molecular conformations. In this methodology, vibrational excitation has been typically applied to the first stretching overtones (∼80 kJ mol^-1) but also to the fundamental modes (∼40 kJ mol^-1). Here, we demonstrate that selective conformational isomerizations are also achieved using excitation to second stretching overtones (∼120 kJ mol^-1). The extremely weak absorptions of the second stretching overtones of molecules isolated in low-temperature matrices were measured for the first time; here using three prototype molecules: hydroxyacetone (HA), glycolic acid (GAc) and glycolamide (GAm). Benchmarking of computed anharmonic IR spectra showed that the B3LYP/SNSD method provides the best agreement with experimental frequencies of the ν(OH), 2ν(OH) and 3ν(OH) modes for the studied molecules in argon matrices. Selective irradiation at the 3ν(OH) frequencies (9850-10 500 cm^-1) of HA, GAc and GAm monomers in argon matrices at 15 K successfully triggers their conformational isomerization. These results open the door to extend control over conformations separated by higher barriers and to induce other transformations not energetically accessible by excitation to the fundamental or first stretching overtone modes.

The role of dispersion and anharmonic corrections in conformational analysis of flexible molecules: the allyl group rotamerization of matrix isolated safrole

Conformational changes of the monomeric safrole (5-(2-propenyl)-1,3-benzodioxole) isolated in low temperature xenon matrices were induced thermally or using narrow-band UV radiation. The rotation of the allyl group taking place in the studied matrices was followed by FTIR spectroscopy. Safrole represents a challenging example of a flexible molecule highlighting the importance of dispersion interactions and anharmonic effects in the structural, spectroscopic and energetic analysis. Structures of the safrole conformers, their energetics and infrared spectra have been calculated using various computational methods ranging from density functional theory (DFT) to coupled cluster (CC). The best theoretical results were obtained by integrating CCSD(T) energies including complete basis set extrapolation and core-valence corrections with B2PLYP-D3 equilibrium structures and hybrid B2PLYP-D3/B3LYP-D3 anharmonic computations for IR spectra and thermodynamics.

Nunes C, Tarczay G, Fausto R. Selective conformational control by excitation of NH imino vibrational antennas

We provide experimental evidence for the occurrence of selective and reversible conformational control over the SH group by vibrational excitation of remote NH groups. Using an imino group that acts as a molecular antenna has no precedents.

Nunes C, Fausto R, Reva I. Conformational control over an aldehyde fragment by selective vibrational excitation of interchangeable remote antennas

We apply vibrational antennas (OH or NH₂ group) to achieve unprecedented conformational control over the heavy aldehyde fragment in 2-formyl-2H-azirine, using selective vibrational excitations of the OH or NH₂ stretching overtones and combination modes.

Experimental Evidence of Long-Range Intramolecular Vibrational Energy Redistribution through Eight Covalent Bonds: NIR Irradiation Induced Conformational Transformation of E-Glutaconic Acid

Long-range intramolecular vibrational energy redistribution (IVR) driven conformational changes were investigated in a matrix-isolated open-chain, asymmetrical dicarboxylic acid, E-glutaconic acid. Although the analysis was challenging due to the presence of multiple backbone conformers and short lifetimes of the prepared higher energy cis conformers, it was shown that the selective excitation of the O-H stretching overtone of one of the carboxylic groups can induce the conformational change (trans to cis) of the other carboxylic group, located at the other end of the E-glutaconic acid molecule. This is a direct proof that the IVR process can act through eight covalent bonds in a flexible molecule before the excess energy completely dissipates into the matrix. The lifetime of the prepared higher energy conformers (averaged over the different backbones) was measured to be 12 s.

Theoretical Investigations of the Chiral Transition of α-Amino Acid Confined in Various Sized Armchair Boron-Nitride Nanotubes

We computationally study the chiral transition process of the α-Ala molecule under confined different sizes of armchair SWBNNTs to explore the confinement effect. We find that the influence of a confinement environment (in armchair SWBNNTs) on the α-Ala molecule would lead to different reaction pathways. Meanwhile, the preferred reaction pathway is also different in various sizes of armchair SWBNNTs, and their energy barriers for the rate-limiting step decrease rapidly with the decreasing of the diameters of the nanotubes. It is obvious that significant decrease of the chiral transition energy barrier occurs compared with the isolated α-Ala molecule chirality conversion mechanism, by ∼15.6 kcal mol^-1, highlighting the improvement in the activity the enantiomers of α-Ala molecule. We concluded that the confinement environment has a significant impact at the nanoscale on the enantiomer transformation process of the chiral molecule.

Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases

The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.