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.

Conformational Mapping, Interactions, and Fluorine Impact by Combined Spectroscopic Approaches and Quantum Chemical Calculations

Noncovalent interactions and their careful variation can be crucial in understanding molecular structures, conformational topographies, and properties. Here, we examine the fluorination impact on the structure and conformational behavior of 2-(2-fluorophenyl)ethyl alcohol (2-FPEAL) by monitoring the first individual ionization-loss-stimulated Raman spectra of the jet-cooled molecule. The comparison of two different broad-range spectra and predicted equivalents discloses two distinct structures. One possesses a folded side chain (gauche) and the other an extended chain (anti) with the terminal hydrogen atom pointing opposite or toward the fluorine side, indicating the improper previous tentative assignment of the latter. These conformers resemble and differ from the nonfluorinated analogue structures. Theoretical analyses reveal interconversion pathways of 2-FPEAL conformers during expansion and the delicate balance between attractive (C-H···F and O-H···π) and repulsive interactions. These findings show the achievements of our integrated approach, suggesting its potential for overcoming future structural challenges.

Searching low-energy conformers of neutral and protonated di-, tri-, and tetra-glycine using first-principles accuracy assisted by the use of neural network potentials

In the last ten years, combinations of state-of-the-art gas-phase spectroscopies and quantum chemistry calculations have suggested several intuitive trends in the structure of small polypeptides that may not hold true. For example, the preference for the cis form of the peptide bond and multiple protonated sites was proposed by comparing experimental spectra with low-energy minima obtained from limited structural sampling using various density functional theory methods. For understanding the structures of polypeptides, extensive sampling of their configurational space with high-accuracy computational methods is required. In this work, we demonstrated the use of deep-learning neural network potential (DL-NNP) to assist in exploring the structure and energy landscape of di-, tri-, and tetra-glycine with the accuracy of high-level quantum chemistry methods, and low-energy conformers of small polypeptides can be efficiently located. We hope that the structures of these polypeptides we found and our preliminary analysis will stimulate further experimental investigations.

Indication of 310-Helix Structure in Gas-Phase Neutral Pentaalanine

We investigate the gas-phase structure of the neutral pentaalanine peptide. The IR spectrum in the 340-1820 cm^-1 frequency range is obtained by employing supersonic jet cooling, infrared multiphoton dissociation, and vacuum-ultraviolet action spectroscopy. Comparison with quantum chemical spectral calculations suggests that the molecule assumes multiple stable conformations, mainly of two structure types. In the most stable conformation theoretically found, the N-terminus forms a C5 ring and the backbone resembles that of an 3₁₀-helix with two β-turns. Additionally, the conformational preferences of pentaalanine have been evaluated using Born-Oppenheimer molecular dynamics, showing that a nonzero simulation time step causes a systematic frequency shift.

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.

Revisiting the Spectroscopy of Water Dimer in Jets

Laser spectroscopy in jets is one of the main sources of structural data from molecular aggregates. Consequently, numerous and sophisticated experimental systems have been developed to extract precise information, which is usually interpreted in the light of quantum mechanical calculations. However, even with the most sophisticated experiments, it is sometimes difficult to interpret the experimental results. We present here the example of water dimer and how after almost 70 years, the assignment of its mass-resolved IR spectrum still generates controversy that extends toward the mechanism of ionization of water aggregates.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

A rotational study of the AlaAla dipeptide

Herein, we present the first rotational study of the AlaAla dipeptide, brought into the gas phase by laser ablation. Two different structures have been unveiled in the isolated environment of a supersonic expansion by Fourier transform microwave spectroscopy. These structures have been identified through their rotational and 14N quadrupole coupling constants. The flexibility of the -NH2 and -COOH ends allows the formation of strong intramolecular interactions giving rise to five- and seven-membered ring configurations.

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.