1
|
Xu R, Jiang Z, Yang Q, Bloino J, Biczysko M. Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers. J Comput Chem 2024; 45:1846-1869. [PMID: 38682874 DOI: 10.1002/jcc.27377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 04/01/2024] [Accepted: 04/02/2024] [Indexed: 05/01/2024]
Abstract
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.
Collapse
Affiliation(s)
- Ruiqin Xu
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| | | | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Science, Prague, Czechia
| | - Julien Bloino
- Classe di Scienze, Scuola Normale Superiore, Pisa, Italy
| | - Malgorzata Biczysko
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
2
|
Moore CC, Staroverov VN, Konermann L. Using Density Functional Theory for Testing the Robustness of Mobile-Proton Molecular Dynamics Simulations on Electrosprayed Ions: Structural Implications for Gaseous Proteins. J Phys Chem B 2023; 127:4061-4071. [PMID: 37116098 DOI: 10.1021/acs.jpcb.3c01581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/30/2023]
Abstract
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.
Collapse
Affiliation(s)
- Conrad C Moore
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Viktor N Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Lars Konermann
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
3
|
Chakraborty S, Mandal K, Ramakrishnan R. Understanding the Role of Intramolecular Ion-Pair Interactions in Conformational Stability Using an Ab Initio Thermodynamic Cycle. J Phys Chem B 2023; 127:648-660. [PMID: 36638237 DOI: 10.1021/acs.jpcb.2c06803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
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).
Collapse
Affiliation(s)
| | - Kalyaneswar Mandal
- Tata Institute of Fundamental Research Hyderabad, Hyderabad500046, India
| | | |
Collapse
|
4
|
Carlo MJ, Patrick AL. Infrared multiple photon dissociation (IRMPD) spectroscopy and its potential for the clinical laboratory. J Mass Spectrom Adv Clin Lab 2022; 23:14-25. [PMID: 34993503 PMCID: PMC8713122 DOI: 10.1016/j.jmsacl.2021.12.004] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 12/09/2021] [Accepted: 12/09/2021] [Indexed: 11/29/2022] Open
Abstract
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.
Collapse
Key Words
- 2-AEP, 2-aminoethylphosphonic acid
- 2P1EA, 2-phenyl-1-ethanolamine
- CIVP, cryogenic ion vibrational predissociation spectroscopy
- CLIO, Centre Laser Infrarouge d’Orsay
- DFT, density functional theory
- FA, fluoroamphetamine
- FEL, free electron laser
- FELIX, Free Electron Laser for Infrared eXperiments
- FMA, fluoromethamphetamine
- FTICR, Fourier transform ion cyclotron resonance
- GC–MS, gas chromatography-mass spectrometry
- GSNO, S- nitro glutathione
- GlcNAc, n-Acetylglucosamine
- IR, infrared
- IR2MS3, infrared-infrared double-resonance multi-stage mass spectrometry
- IRMPD, infrared multiple photon dissociation (IRMPD)
- IRMPD-MS, infrared multiple photon dissociation spectroscopy mass spectrometry
- IRPD, infrared predissociation spectroscopy
- IVR, intramolecular vibrational redistribution
- Infrared multiple photon dissociation spectroscopy
- LC, liquid chromatography
- LC-MS, liquid chromatography-mass spectrometry
- LC-MS/MS, liquid chromatography-tandem mass spectrometry
- MDA, methylenedioxyamphetamine
- MDMA, methylenedioxymethamphetamine
- MMC, methylmethcathinone
- MS/MS, tandem mass spectrometry
- MSn, multi-stage mass spectrometry
- Mass spectrometry
- Metabolites
- NANT, N-acetyl-N-nitrosotryptophan
- OPO/A, optical parametric oscillator/amplifier
- PTM, post-translational modification
- Pharmaceuticals
- Post-translational modifications
- SNOCys, S-nitrosocysteine
- UV, ultraviolet
- UV-IR, ultraviolet-infrared
- Vibrational spectroscopy
- cw, continuous wave
- α-PVP, alpha-pyrrolidinovalerophenone
Collapse
Affiliation(s)
- Matthew J. Carlo
- Department of Chemistry, Mississippi State University, Mississippi State, MS 39762, USA
| | - Amanda L. Patrick
- Department of Chemistry, Mississippi State University, Mississippi State, MS 39762, USA
| |
Collapse
|
5
|
Becher S, Berden G, Martens J, Oomens J, Heiles S. IRMPD Spectroscopy of [PC (4:0/4:0) + M] + (M = H, Na, K) and Corresponding CID Fragment Ions. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2021; 32:2874-2884. [PMID: 34723538 DOI: 10.1021/jasms.1c00277] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
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 (MS2) 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.
Collapse
Affiliation(s)
- Simon Becher
- Institute of Inorganic and Analytical Chemistry, Justus Liebig University Giessen, 35392 Giessen, Germany
| | - Giel Berden
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, 6525 ED Nijmegen, The Netherlands
| | - Jonathan Martens
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, 6525 ED Nijmegen, The Netherlands
| | - Jos Oomens
- FELIX Laboratory, Institute for Molecules and Materials, Radboud University, 6525 ED Nijmegen, The Netherlands
| | - Sven Heiles
- Institute of Inorganic and Analytical Chemistry, Justus Liebig University Giessen, 35392 Giessen, Germany
| |
Collapse
|
6
|
Darii E, Gimbert Y, Alves S, Damont A, Perret A, Woods AS, Fenaille F, Tabet JC. First Direct Evidence of Interpartner Hydride/Deuteride Exchanges for Stored Sodiated Arginine/Fructose-6-phosphate Complex Anions within Salt-Solvated Structures. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2021; 32:1424-1440. [PMID: 33929837 DOI: 10.1021/jasms.1c00040] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ekaterina Darii
- Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, 91057 Evry, France
| | - Yves Gimbert
- Département de Chimie Moléculaire, UMR CNRS 5250, Université Grenoble Alpes, 38058 Grenoble, France
- Sorbonne Université, Faculté des Sciences et de l'Ingénierie, Institut Parisien de Chimie Moléculaire (IPCM), F-75005 Paris, France
| | - Sandra Alves
- Sorbonne Université, Faculté des Sciences et de l'Ingénierie, Institut Parisien de Chimie Moléculaire (IPCM), F-75005 Paris, France
| | - Annelaure Damont
- Université Paris-Saclay, CEA, INRAE, Département Médicaments et Technologies pour la Santé (DMTS), MetaboHUB, F-91191 Gif sur Yvette, France
| | - Alain Perret
- Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, 91057 Evry, France
| | - Amina S Woods
- NIDA IRP, NIH Structural Biology Unit Cellular Neurobiology Branch, 333 Cassell Drive, Baltimore, Maryland 21224, United States
- The Johns Hopkins University School of Medicine, Pharmacology and Molecular Sciences, Baltimore, Maryland 21205, United States
| | - François Fenaille
- Université Paris-Saclay, CEA, INRAE, Département Médicaments et Technologies pour la Santé (DMTS), MetaboHUB, F-91191 Gif sur Yvette, France
| | - Jean-Claude Tabet
- Sorbonne Université, Faculté des Sciences et de l'Ingénierie, Institut Parisien de Chimie Moléculaire (IPCM), F-75005 Paris, France
- Université Paris-Saclay, CEA, INRAE, Département Médicaments et Technologies pour la Santé (DMTS), MetaboHUB, F-91191 Gif sur Yvette, France
| |
Collapse
|
7
|
Medina FE, Prejanò M. Water Molecules Allow the Intramolecular Activation of the Thiamine Di-Phosphate Cofactor in Human Transketolase: Mechanistic Insights into a Famous Proposal. ACS Catal 2021. [DOI: 10.1021/acscatal.0c05414] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Fabiola E. Medina
- Departamento de Ciencias Químicas, Facultad de Ciencias Exactas, Universidad Andres Bello, Autopista Concepción-Talcahuano, Talcahuano 7100, Chile
| | - Mario Prejanò
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm SE-10691, Sweden
| |
Collapse
|
8
|
Chiavarino B, Sinha RK, Crestoni ME, Corinti D, Filippi A, Fraschetti C, Scuderi D, Maitre P, Fornarini S. Binding Motifs in the Naked Complexes of Target Amino Acids with an Excerpt of Antitumor Active Biomolecule: An Ion Vibrational Spectroscopy Assay. Chemistry 2021; 27:2348-2360. [PMID: 33175428 DOI: 10.1002/chem.202003555] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 09/30/2020] [Indexed: 11/06/2022]
Abstract
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.
Collapse
Affiliation(s)
- Barbara Chiavarino
- Dipartimento di Chimica e Tecnologie del Farmaco, Università degli studi di Roma La Sapienza, P.le A. Moro 5, 00185, Roma, Italy
| | - Rajeev K Sinha
- Department of Atomic and Molecular Physics, Manipal University, Manipal, 576104, Karnataka, India
| | - Maria Elisa Crestoni
- Dipartimento di Chimica e Tecnologie del Farmaco, Università degli studi di Roma La Sapienza, P.le A. Moro 5, 00185, Roma, Italy
| | - Davide Corinti
- Dipartimento di Chimica e Tecnologie del Farmaco, Università degli studi di Roma La Sapienza, P.le A. Moro 5, 00185, Roma, Italy
| | - Antonello Filippi
- Dipartimento di Chimica e Tecnologie del Farmaco, Università degli studi di Roma La Sapienza, P.le A. Moro 5, 00185, Roma, Italy
| | - Caterina Fraschetti
- Dipartimento di Chimica e Tecnologie del Farmaco, Università degli studi di Roma La Sapienza, P.le A. Moro 5, 00185, Roma, Italy
| | - Debora Scuderi
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405, Orsay, France
| | - Philippe Maitre
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405, Orsay, France
| | - Simonetta Fornarini
- Dipartimento di Chimica e Tecnologie del Farmaco, Università degli studi di Roma La Sapienza, P.le A. Moro 5, 00185, Roma, Italy
| |
Collapse
|
9
|
Zhou C, Ieritano C, Hopkins WS. Augmenting Basin-Hopping With Techniques From Unsupervised Machine Learning: Applications in Spectroscopy and Ion Mobility. Front Chem 2019; 7:519. [PMID: 31440497 PMCID: PMC6693329 DOI: 10.3389/fchem.2019.00519] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Accepted: 07/08/2019] [Indexed: 12/17/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ce Zhou
- Department of Chemistry, University of Waterloo, Waterloo, ON, Canada
| | | | | |
Collapse
|