Toward accurate prediction of amino acid derivatives structure and energetics from DFT: glycine conformers and their interconversions

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.

Multifactorial interaction of selenium, iron, xylose, and glycine on cordycepin metabolism in Cordyceps militaris

Cordycepin, a nucleoside analog, is the main antioxidative and antimicrobial substance in Cordyceps militaris. To improve the metabolism of cordycepin, carbon sources, nitrogen sources, trace elements, and precursors were studied by single factor, Plackett-Burman, and central composite designs in C. militaris mycelial fermentation. Under the regulation of the multifactorial interactions of selenite, ferrous chloride, xylose, and glycine, cordycepin production was increased by 5.2-fold compared with the control. The gene expression of hexokinase, ATP phosphoribosyltransferase, adenylosuccinate synthetase, and cns1-3 in the glycolysis, pentose phosphate, and adenosine synthesis pathways were increased by 3.2-7.5 times due to multifactorial interactions, while the gene expression of histidine biosynthesis trifunctional protein and histidinol-phosphate aminotransferase in histidine synthesis pathway were decreased by 23.4%-56.2%. Increasing with cordycepin production, glucose uptake was accelerated, mycelia growth was inhibited, and the cell wall was damaged. Selenomethionine (SeMet), selenocysteine (SeCys), and selenium nanoparticles (SeNPs) were the major Se species in C. militaris mycelia. This study provides a new insight for promoting cordycepin production by regulating glycolysis, pentose phosphate, and histidine metabolism. KEY POINTS: • Cordycepin production in the CCDmax group was 5.2-fold than that of the control. • Glucose uptake of the CCDmax group was accelerated and cell wall was damaged. • The metabolic flux was concentrated to the cordycepin synthesis pathway.

Coupled Cluster Semiclassical Estimates of Experimental Reaction Rates: The Interconversion of Glycine Conformer VIp to Ip

We apply the full-dimensional Semiclassical Transition State Theory (SCTST) to estimate the rate constant of glycine molecule interconversion between the VIp and Ip conformers. We have reached an electronic structure accuracy up to the explicitly correlated Coupled Cluster method (CCSD(T)-F12b/cc-pVDZ-F12) thanks to our parallel implementation. The reaction has been experimentally investigated in the literature and is known to proceed by quantum mechanical tunneling. The SCTST rates improve over other theoretical methods, and our results align with the experimental measurements, thus confirming the accuracy of the fully coupled anharmonic semiclassical tunneling treatment, providing that the level of electronic structure theory gives a reliable estimate of the reaction barrier height and shape. The comparison with experimental half-life times supports the validity of SCTST for glycine VIp-Ip conformer conversion in the cryogenic temperature range, where this theory is usually not considered applicable due to the onset of the deep tunneling regime.

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.

Probing the Reactions of Thiourea (CH4N2S) with Metals (X = Au, Hf, Hg, Ir, Os, W, Pt, and Re) Anchored on Fullerene Surfaces (C59X)

Upon various investigations conducted in search for a nanosensor material with the best sensing performance, the need to explore these materials cannot be overemphasized as materials associated with best sensing attributes are of vast interest to researchers. Hence, there is a need to investigate the adsorption performances of various metal-doped fullerene surfaces: C₅₉Au, C₅₉Hf, C₅₉Hg, C₅₉Ir, C₅₉Os, C₅₉Pt, C₅₉Re, and C₅₉W on thiourea [SC(NH₂)₂] molecule using first-principles density functional theory computation. Comparative adsorption study has been carried out on various adsorption models of four functionals, M06-2X, M062X-D3, PBE0-D3, and ωB97XD, and two double-hybrid (DH) functionals, DSDPBEP86 and PBE0DH, as reference at Gen/def2svp/LanL2DZ. The visual study of weak interactions such as quantum theory of atoms in molecule analysis and noncovalent interaction analysis has been invoked to ascertain these results, and hence we arrived at a conclusive scientific report. In all cases, the weak adsorption observed is best described as physisorption phenomena, and CH₄N₂S@C₅₉Pt complex exhibits better sensing attributes than its studied counterparts in the interactions between thiourea molecule and transition metal-doped fullerene surfaces. Also, in the comparative adsorption study, DH density functionals show better performance in estimating the adsorption energies due to their reduced mean absolute deviation (MAD) and root-mean-square deviation (RMSD) values of (MAD = 1.0305, RMSD = 1.6277) and (MAD = 0.9965, RMSD = 1.6101) in DSDPBEP86 and PBE0DH, respectively.

DFT study on the synthesis of trifluoroacetophenone from palladium complex LnPd(Ph)CF3 (Ln = Xantphos or DtBPF) and CO

The carbonylative trifluoromethylation reaction mechanism of palladium complex L_nPd(Ph)CF₃ (L_n = Xantphos or DtBPF) and CO to synthesize trifluoroacetophenone was calculated using the density functional theory B3LYP method. In this paper, we conducted a computational study on the competition mechanism of two different products trifluorotoluene and trifluoroacetophenone. The calculation result reveals (1) CO insertion and reduction-elimination are two key steps in palladium-catalyzed reactions; (2) for the palladium complex (Xantphos)Pd(Ph)CF₃, the resulting product trifluoroacetyl has a lower activation energy and higher reactivity; and (3) for the metal palladium ligand DtBPF, the small energy difference between the two products indicates that the stereoselectivity of the product is relatively poor. The computational research results in this paper provide a good supplement and effective explanation to the experimental phenomenon of Domino et al. (Organometallics 39:688-697, 2020).

DFT study on the gold(I)-catalyzed cycloaddition and rearrangement reactions of allene-containing allylic silyl ether

The DFT calculation at the B3LYP/B3LYP-D3(BJ) level was carried out to explore the reaction mechanism of the synthesis of spirocyclo[4,5]decane skeleton by gold-catalyzed allenyl compounds. The more accurate energy under the CH₃CN solvent in the experiment is calculated by the single-point energy of the SMD model. Computational studies have shown that the reaction consists of three main steps: intramolecular cycloaddition of the end group carbon atoms of allenyl and vinyl groups, the semipinacol rearrangement process in which the four-membered ring is reconstructed into the five-membered ring, the elimination reaction releases the catalyst and obtains the product. The calculation results show that Zheng et al. reported that the gold-catalyzed synthesis reaction can easily occur under the experimental conditions due to its low activation free energy (12.07-15.49 kcal/mol). Furthermore, it was found that the MOMO(CH₂)₂ substituent has higher reactivity than the corresponding reactant of the phenyl substituent.

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.

Anharmonic DFT Study of Near-Infrared Spectra of Caffeine: Vibrational Analysis of the Second Overtones and Ternary Combinations

Anharmonic quantum chemical calculations were employed to simulate and interpret a near-infrared (NIR) spectrum of caffeine. First and second overtones, as well as binary and ternary combination bands, were obtained, accurately reproducing the lineshape of the experimental spectrum in the region of 10,000–4000 cm⁻¹ (1000–2500 nm). The calculations enabled performing a detailed analysis of NIR spectra of caffeine, including weak bands due to the second overtones and ternary combinations. A highly convoluted nature of NIR spectrum of caffeine was unveiled, with numerous overlapping bands found beneath the observed spectral lineshape. To properly reflect that intrinsic complexity, the band assignments were provided in the form of heat maps presenting the contributions to the NIR spectrum from various kinds of vibrational transitions. These contributions were also quantitatively assessed in terms of the integral intensities. It was found that the combination bands provide the decisively dominant contributions to the NIR spectrum of caffeine. The first overtones gain significant importance between 6500–5500 cm⁻¹, while the second overtones are meaningful in the higher wavenumber regions, particularly in the 10,000–7000 cm⁻¹ region. The obtained detailed band assignments enabled deep interpretation of the absorption regions of caffeine identified in the literature as meaningful for analytical applications of NIR spectroscopy focused on quantitative analysis of caffeine content in drugs and natural products.

Structural Water Stabilizes Protein Motifs in Liquid Protein Phase: The Folding Mechanism of Short β-Sheets Coupled to Phase Transition

Macromolecular associates, such as membraneless organelles or lipid-protein assemblies, provide a hydrophobic environment, i.e., a liquid protein phase (LP), where folding preferences can be drastically altered. LP as well as the associated phase change from water (W) is an intriguing phenomenon related to numerous biological processes and also possesses potential in nanotechnological applications. However, the energetic effects of a hydrophobic yet water-containing environment on protein folding are poorly understood. Here, we focus on small β-sheets, the key motifs of proteins, undergoing structural changes in liquid–liquid phase separation (LLPS) and also model the mechanism of energy-coupled unfolding, e.g., in proteases, during W → LP transition. Due to the importance of the accurate description for hydrogen bonding patterns, the employed models were studied by using quantum mechanical calculations. The results demonstrate that unfolding is energetically less favored in LP by ~0.3–0.5 kcal·mol⁻¹ per residue in which the difference further increased by the presence of explicit structural water molecules, where the folded state was preferred by ~1.2–2.3 kcal·mol⁻¹ per residue relative to that in W. Energetics at the LP/W interfaces was also addressed by theoretical isodesmic reactions. While the models predict folded state preference in LP, the unfolding from LP to W renders the process highly favorable since the unfolded end state has >1 kcal·mol⁻¹ per residue excess stabilization.

Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy

The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.