Electronic Energies Are Not Enough: An Ion Mobility-Aided, Quantum Chemical Benchmark Analysis of H+GPGG Conformers

A Direct Measurement of the Absolute Energies of Protonated Gly-Pro-Gly-Gly Conformations Using Ion Mobility Spectrometry and Guided Ion Beam Tandem Mass Spectrometry

Threshold collision-induced dissociation of conformations selected by ion mobility spectrometry is described as a method to determine their absolute relative energies. Here, the method is demonstrated with the cis and trans conformers of the protonated tetrapeptide, Gly-Pro-Gly-Gly, which differ in the orientation of the peptide bond at the proline residue. The trans conformation was found to be more stable than the cis conformation by 4.5 ± 2.5 kJ mol-1 at 0 K. Heating the molecule anneals it to the cis conformer, indicating it has a lower Gibbs energy at higher temperatures. These results are compared to theoretical values calculated here and from the literature. This experimental analysis is the first quantitative measurement of the relative stability of the conformers of a protonated peptide in the gas phase, which has far-reaching impacts on the selection of theoretical methods to describe the energetics of these systems.

Tautomerization of H+KPGG: Entropic Consequences of Strong Hydrogen-Bond Networks in Peptides

Ion mobility spectrometry-mass spectrometry and quantum chemical calculations are used to determine the structures and stabilities of the singly protonated peptide H+KPGG. The two peaks making up the IMS distribution are shown to be tautomers differing by the location of the extra proton on either the lysine side chain or the N-terminus. The lysine-protonated tautomer is strongly preferred entropically while being disfavored in terms of the electronic energy and enthalpy. This relationship is shown, through comparison of all low-lying conformers of both tautomers, to be related to the strong hydrogen-bond network of the N-terminally protonated tautomer. A general relationship is demonstrated wherein stronger cross-peptide hydrogen-bond networks result in entropically disfavored conformers. Further effects of the H+KPGG hydrogen-bond network are probed by computationally examining singly and doubly methylated analogues. These results demonstrate the importance of the entropic consequences of hydrogen bonds to peptide stability as well as techniques for perturbing the hydrogen-bond network and folding preferences of peptides via minimal chemical modification.

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.

Interaction-Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions

The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG⁺G⁺TG⁺G⁺g^-)-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as "preparation energy" or "reorganization energy". We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme's efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 3₁₀-polyglycine and H⁺GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction-deletion in 2-nitrobenzene-1,3-diol.

Untangling Hydrogen Bond Networks with Ion Mobility Spectrometry and Quantum Chemical Calculations: A Case Study on H+XPGG

Ion mobility spectrometry-mass spectrometry and quantum chemical calculations are used to determine the structures and stabilities of singly protonated XaaProGlyGly peptides: H⁺DPGG, H⁺NPGG, H⁺EPGG, and H⁺QPGG. The IMS distributions are similar, suggesting the peptides adopt closely related structures in the gas phase. Quantum chemical calculations show that all conformers seen in the experimental spectrum correspond to the cis configuration about the Xaa-Pro peptide bond, significantly different from the behavior seen previously for H⁺GPGG. Density functional theory and quantum theory of atoms in molecules (QTAIM) investigations uncover a silent drama as a minor conformer not observed in the H⁺DPGG spectrum becomes the preferred conformer in H⁺QPGG, with both conformers being coincident in collision cross section. Investigation of the highly coupled hydrogen bond network, replete with CH···O interactions and bifurcated hydrogen bonds, reveals the cause of this effect as well as the absence of trans conformers from the spectra. A series of generalized observations are provided to aid in enzyme and ligand design using these coupled hydrogen bond motifs.

A Trip to the Density Functional Theory Zoo: Warnings and Recommendations for the User

This account is written for general users of density functional theory (DFT) methods as well as experimental researchers who are new to the field and would like to conduct such calculations. Its main emphasis lies on how to find a way through the confusing ‘zoo’ of DFT by addressing common misconceptions and highlighting those modern methods that should ideally be used in calculations of energetic properties and geometries. A particular focus is on highly popular methods and the important fact that popularity does not imply accuracy. In this context, we present a new analysis of the openly available data published in Swart and co-workers’ famous annual ‘DFT poll’ (http://www.marcelswart.eu/dft-poll/) to demonstrate the existing communication gap between the DFT user and developer communities. We show that despite considerable methodological advances in the field, the perception of some parts of the user community regarding their favourite approaches has changed little. It is hoped that this account makes a contribution towards changing this status and that users are inspired to adjust their current computational protocols to accommodate strategies that are based on proven robustness, accuracy, and efficiency rather than popularity.