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.

Using Density Functional Theory for Testing the Robustness of Mobile-Proton Molecular Dynamics Simulations on Electrosprayed Ions: Structural Implications for Gaseous Proteins

Current experiments only provide low-resolution information on gaseous protein ions generated by electrospray ionization (ESI). Molecular dynamics (MD) simulations can yield complementary insights. Unfortunately, conventional MD does not capture the mobile nature of protons in gaseous proteins. Mobile-proton MD (MPMD) overcomes this limitation. Earlier MPMD data at 300 K indicated that protein ions generated by "native" ESI retain solution-like structures with a hydrophobic core and zwitterionic exterior [Bakhtiari, M.; Konermann, L. J. Phys. Chem. B 2019, 123, 1784-1796]. MPMD redistributes protons using electrostatic and proton affinity calculations. The robustness of this approach has never been scrutinized. Here, we close this gap by benchmarking MPMD against density functional theory (DFT) at the B3LYP/6-31G* level, which is well suited for predicting proton affinities. The computational cost of DFT necessitated the use of small peptides. The MPMD energetic ranking of proton configurations was found to be consistent with DFT single-point energies, implying that MPMD can reliably identify favorable protonation sites. Peptide MPMD runs converged to DFT-optimized structures only when applying 300-500 K temperature cycling, which was necessary to prevent trapping in local minima. Temperature cycling MPMD was then applied to gaseous protein ions. Native ubiquitin converted to slightly expanded structures with a zwitterionic core and a nonpolar exterior. Our data suggest that such inside-out protein structures are intrinsically preferred in the gas phase, and that they form in ESI experiments after moderate collisional excitation. This is in contrast to native ESI (with minimal collisional excitation, simulated by MPMD at 300 K), where kinetic trapping promotes the survival of solution-like structures. In summary, this work validates the MPMD approach for simulations on gaseous peptides and proteins.

Understanding the Role of Intramolecular Ion-Pair Interactions in Conformational Stability Using an Ab Initio Thermodynamic Cycle

Intramolecular ion-pair interactions yield shape and functionality to many molecules. With proper orientation, these interactions overcome steric factors and are responsible for the compact structures of several peptides. In this study, we present a thermodynamic cycle based on isoelectronic and alchemical mutation to estimate the intramolecular ion-pair interaction energy. We determine these energies for 26 benchmark molecules with common ion-pair combinations and compare them with results obtained using intramolecular symmetry-adapted perturbation theory. For systems with long linkers, the ion-pair energies evaluated using both approaches deviate by less than 2.5% in the vacuum phase. The thermodynamic cycle based on density functional theory facilitates calculations of salt-bridge interactions in model tripeptides with continuum/microsolvation modeling and four large peptides: 1EJG (crambin), 1BDK (bradykinin), 1L2Y (a mini-protein with a tryptophan cage), and 1SCO (a toxin from the scorpion venom).

Infrared multiple photon dissociation (IRMPD) spectroscopy and its potential for the clinical laboratory

Infrared multiple photon dissociation (IRMPD) spectroscopy is a powerful tool used to probe the vibrational modes-and, by extension, the structure-of an ion within an ion trap mass spectrometer. Compared to traditional FTIR spectroscopy, IRMPD spectroscopy has advantages including its sensitivity and its relative ability to handle complex mixtures. While IRMPD has historically been a technique for fundamental analyses, it is increasingly being applied in a more analytical fashion. Notable recent demonstrations pertinent to the clinical laboratory and adjacent interests include analysis of modified amino acids/residues and carbohydrates, structural elucidation (including isomeric differentiation) of metabolites, identification of novel illicit drugs, and structural studies of various biomolecules and pharmaceuticals. Improvements in analysis time, coupling to commercial instruments, and integration with separations methods are all drivers toward the realization of these analytical applications. Additional improvements in these areas, along with advances in benchtop tunable IR sources and increased cross-discipline collaboration, will continue to drive innovation and widespread adoption. The goal of this tutorial article is to briefly present the fundamentals and instrumentation of IRMPD spectroscopy, as an overview of the utility of this technique for helping to answer questions relevant to clinical analysis, and to highlight limitations to widespread adoption, as well as promising directions in which the field may be heading.

IRMPD Spectroscopy of [PC (4:0/4:0) + M]+ (M = H, Na, K) and Corresponding CID Fragment Ions

Glycerophospholipids (GPs) are highly abundant in eukaryotic cells and take part in numerous fundamental physiological processes such as molecular signaling. The GP composition of samples is often analyzed using mass spectrometry (MS), but identification of some structural features, for example, differentiation of stereospecific numbering (sn) isomers by well-established tandem MS (MS²) methods, is challenging. In particular, the formation of 1,3-dioxolane over 1,3-dioxane intermediates proposed to be responsible for the sn-selectivity of these tandem MS strategies has not been validated by spectroscopic methods. In this work, we present infrared multiple photon dissociation (IRMPD) spectra of phosphatidylcholine (PC) ions [PC 4:0/4:0 + H/Na/K]⁺ and [PC 4:0/4:0 + Na/K - 183]⁺ fragments generated by electrospray ionization (ESI)-MS and collision-induced dissociation (CID), respectively. IRMPD spectra of protonated, sodiated, and potassiated PC 4:0/4:0 differ in the phosphate- and ester-related bands, which are increasingly shifted to lower wavenumbers with higher adduct masses. Comparison of calculated and experimental IR spectra indicates the presence of multiple, two and one isomer(s) for [PC 4:0/4:0 + H]⁺, [PC 4:0/4:0 + Na]⁺, and [PC 4:0/4:0 + K]⁺, respectively. Isomers exhibiting pronounced sn-1 ester-ion interactions are computationally predicted to be energetically preferred for all species and are in line with experimental results. IRMPD spectra of [PC 4:0/4:0 + Na/K - 183]⁺ are presented and shed the first light on the fragment ion structures, rationalizing MS-based lipidomics strategies that aim to characterize the sn-isomerism of GPs.

First Direct Evidence of Interpartner Hydride/Deuteride Exchanges for Stored Sodiated Arginine/Fructose-6-phosphate Complex Anions within Salt-Solvated Structures

Mass spectrometric investigations of noncovalent binding between low molecular weight compounds revealed the existence of gas-phase (GP) noncovalent complex (NCC) ions involving zwitterionic structures. ESI MS is used to prove the formation of stable sodiated NCC anions between fructose (F6P) and arginine (R) moieties. Theoretical calculations indicate a folded solvated salt (i.e., sodiated carboxylate interacting with phosphate) rather than a charge-solvated form. Under standard CID conditions, [(F6P+R-H+Na)-H]^- competitively forms two major product ions (PIs) through partner splitting [(R-H+Na) loss] and charge-induced cross-ring cleavage while preserving the noncovalent interactions (noncovalent product ions (NCPIs)). MS/MS experiments combined with in-solution proton/deuteron exchanges (HDXs) demonstrated an unexpected labeling of PIs, i.e., a correlated D-enrichment/D-depletion. An increase in activation time up to 3000 ms favors such processes when limited to two H/D exchanges. These results are rationalized by interpartner hydride/deuteride exchanges (⟨HDX⟩) through stepwise isomerization/dissociation of sodiated NCC-d11 anions. In addition, the D-enrichment/D-depletion discrepancy is further explained by back HDX with residual water in LTQ (selective for the isotopologue NCPIs as shown by PI relaxation experiments). Each isotopologue leads to only one back HDX unlike multiple HDXs generally observed in GP. This behavior shows that NCPIs are zwitterions with charges solvated by a single water molecule, thus generating a back HDX through a relay mechanism, which quenches the charges and prevents further back HDX. By estimating back HDX impact on D-depletion, the interpartner ⟨HDX⟩ complementarity was thus illustrated. This is the first description of interpartner ⟨HDX⟩ and selective back HDX validating salt-solvated structures.

Binding Motifs in the Naked Complexes of Target Amino Acids with an Excerpt of Antitumor Active Biomolecule: An Ion Vibrational Spectroscopy Assay

The structures of proton-bound complexes of 5,7-dimethoxy-4H-chromen-4-one (1) and basic amino acids (AAs), namely, histidine (His) and lysine (Lys), have been examined by means of mass spectrometry coupled with IR ion spectroscopy and quantum chemical calculations. This selection of systems is based on the fact that 1 represents a portion of glabrescione B, a natural small molecule of promising antitumor activity, while His and Lys are protein residues lining the cavity of the alleged receptor binding site. These species are thus a model of the bioactive adduct, although clearly the isolated state of the present study bears little resemblance to the complex biological environment. A common feature of [1+AA+H]⁺ complexes is the presence of a protonated AA bound to neutral 1, in spite of the fact that the gas-phase basicity of 1 is comparable to those of Lys and His. The carbonyl group of 1 acts as a powerful hydrogen-bond acceptor. Within [1+AA+H]⁺ the side-chain substituents (imidazole group for His and terminal amino group for Lys) present comparable basic properties to those of the α-amino group, taking part to a cooperative hydrogen-bond network. Structural assignment, relying on the comparative analysis of the infrared multiple photon dissociation (IRMPD) spectrum and calculated IR spectra for the candidate geometries, derives from an examination over two frequency ranges: 900-1800 and 2900-3700 cm^-1 . Information gained from the latter one proved especially valuable, for example, pointing to the contribution of species characterized by an unperturbed carboxylic OH or imidazole NH stretching mode.

Augmenting Basin-Hopping With Techniques From Unsupervised Machine Learning: Applications in Spectroscopy and Ion Mobility

Evolutionary algorithms such as the basin-hopping (BH) algorithm have proven to be useful for difficult non-linear optimization problems with multiple modalities and variables. Applications of these algorithms range from characterization of molecular states in statistical physics and molecular biology to geometric packing problems. A key feature of BH is the fact that one can generate a coarse-grained mapping of a potential energy surface (PES) in terms of local minima. These results can then be utilized to gain insights into molecular dynamics and thermodynamic properties. Here we describe how one can employ concepts from unsupervised machine learning to augment BH PES searches to more efficiently identify local minima and the transition states connecting them. Specifically, we introduce the concepts of similarity indices, hierarchical clustering, and multidimensional scaling to the BH methodology. These same machine learning techniques can be used as tools for interpreting and rationalizing experimental results from spectroscopic and ion mobility investigations (e.g., spectral assignment, dynamic collision cross sections). We exemplify this in two case studies: (1) assigning the infrared multiple photon dissociation spectrum of the protonated serine dimer and (2) determining the temperature-dependent collision cross-section of protonated alanine tripeptide.