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

An influence of electronic structure theory method, thermodynamic and implicit solvation corrections on the organic carbonates conformational and binding energies

An impact of an electronic structure or force field method, gas-phase thermodynamic correction, and continuum solvation model on organic carbonate clusters (S)n conformational and binding energies is explored. None of the tested force field (GFN-FF, GAFF, MMFF94) and standard semiempirical methods (PM3, AM1, RM1, PM6, PM6-D3, PM6-D3H4, PM7) can reproduce reference RI-SCS-MP2 conformational energies. Tight-binding GFNn-xTB methods provide more realistic conformational energies which are accurate enough to discard the least stable conformers. The effect of thermodynamic correction is moderate and can be ignored if the gas phase conformational stability ranking is a goal. The influence of continuum solvation is stronger, especially if reinforced with the Gibbs free energy thermodynamic correction, and results in the reduced spread of conformational energies. The cluster formation binding energies strongly depend on a particular approach to vibrational thermochemistry with the difference between traditional harmonic and modified scaled rigid - harmonic oscillator approximations reaching 10 kcal mol-1.

Describing the Disulfide Bond: From the Density Functional Theory and Back through the "Lego Brick" Approach

Selected molecular species containing the disulfide bond, RSSR, have been considered, these ranging from hydrogen disulfide, H2S2 (R = H), to diphenyl disulfide with R = C6H5. The aim of this work is two-fold: (i) to investigate different computational approaches in order to derive accurate equilibrium structures at an affordable cost, (ii) to employ the results from the first goal in order to benchmark cheaper methodologies rooted in the density functional theory. Among the strategies used for the accurate geometrical determinations, the semiexperimental approach has been exploited in combination with a reduced-dimensionality VPT2 model, without however obtaining satisfactory results. Instead, the so-called "Lego brick" approach turned out to be very effective despite the flexibility of the systems investigated. Concerning the second target of this work, the focus was mainly on the S-S bond and the structural parameters related to it. Among those tested, PBE0(-D3BJ), M06-2X(-D3) and DSD-PBEP86-D3BJ have been found to be the best-performing functionals.

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.

DFT Meets Wave-Function Methods for Accurate Structures and Rotational Constants of Histidine, Tryptophan, and Proline

A new computational strategy has been applied to the conformational and spectroscopic properties in the gas phase of amino acids with very distinctive features, ranging from different tautomeric forms (histidine) to ring puckering (proline), and heteroaromatic structures with non-equivalent rings (tryptophan). The integration of modern double-hybrid functionals and wave-function composite methods has allowed us to obtain accurate results for a large panel of conformers with reasonable computer times. The remarkable agreement between computations and microwave experiments allows an unbiased interpretation of the latter in terms of stereoelectronic effects.

Accurate structures and spectroscopic parameters of α,α-dialkylated α-amino acids in the gas-phase: a joint venture of DFT and wave-function composite methods

Accurate computations of structural, conformational and spectroscopic properties in the gas phase have been performed for two α,α-dialkylated α-amino acids, namely aminoisobutyric acid and cyclopropylglycine. Thanks to the integration of modern double hybrid functionals and wave-function methods, several low-energy structures of the title molecules could be analyzed employing standard computer resources. The computed features of all the most stable conformers of the target amino acids closely match the corresponding spectroscopic parameters issued from microwave spectroscopic studies in the gas-phase. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the way for accurate investigations of other flexible bricks of life.

Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods

A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.

Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy

The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.

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

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.

Upon DFT-D3 dispersion correction and ECD spectral confirmation, only several conformers can stably coexist for three fungal cycloaspeptides (A, D, G)

Dispersion correction in theoretical determination of cyclopeptide conformations is emphasized. Whether in gas approximation or in solvation simulation, the density functional theory with London dispersion correction (DFT-D3) demonstrates that only 2-3 conformers can stably coexist for cycloaspeptides (A, D, G) at B3LYP-D3 and CAM-B3LYP-D3. Conformational rationality is confirmed by electronic circular dichroism (ECD). Whether for Cotton effect or for excitation energy, TD-B3LYP-D3 has better performances than TD-CAM-B3LYP-D3 because the former can better reproduce the experiment. A molecular orbital analysis is used to interpret ECD, where two energy bands observed in experiment originates from the ππ* transitions other than the σπ* transitions. Long-range correction and solvent effect make H-bonds shorten, and dispersion correction makes them further shorten.