Good vibrations: probing biomolecular structure and interactions through spectroscopy in the gas phase

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.

Comprehensive Analysis of Methyl-β-D-ribofuranoside: A Multifaceted Spectroscopic and Theoretical Approach

This study presents a comprehensive analysis of the vibrational spectra of methyl-β-D-ribofuranoside. Employing a combination of inelastic neutron scattering, Raman, and infrared spectroscopy allows for the observation of all modes regardless of the selection rules. The experimental techniques were complemented by density functional theory computational methods using both gas-phase (Gaussian) and solid-state (CRYSTAL, CASTEP) approaches to provide an unambiguous assignment of the defining vibrational features. Two distinct structures of the molecule were identified in the unit cell, differentiated mainly by the orientation of the furanose ring O-H bonds. The low-energy region of the spectrum (<400 cm-1) is dominated by lattice vibrations and functional group rotation, while the midenergy region is dominated by out-of-plane bending motions of the furanose ring (400-900 cm-1) and by C-H bending in the methyl and methylene groups (1400-1600 cm-1). The high-energy region (>2800 cm-1) encompasses the C-H and O-H stretching modes and offers convincing evidence of at least one H-bonding interaction between the two structures of methyl-β-D-ribofuranoside.

Isotope effects on radical pair performance in cryptochrome: A new hypothesis for the evolution of animal migration: The quantum biology of migration

Mechanisms occurring at the atomic level are now known to drive processes essential for life, as revealed by quantum effects on biochemical reactions. Some macroscopic characteristics of organisms may thus show an atomic imprint, which may be transferred across organisms and affect their evolution. This possibility is considered here for the first time, with the aim of elucidating the appearance of an animal innovation with an unclear evolutionary origin: migratory behaviour. This trait may be mediated by a radical pair (RP) mechanism in the retinal flavoprotein cryptochrome, providing essential magnetic orientation for migration. Isotopes may affect the performance of quantum processes through their nuclear spin. Here, we consider a simple model and then apply the standard open quantum system approach to the spin dynamics of cryptochrome RP. We changed the spin quantum number (I) and g-factor of hydrogen and nitrogen isotopes to investigate their effect on RP's yield and magnetic sensitivity. Strong differences arose between isotopes with I = 1 and I = 1/2 in their contribution to cryptochrome magnetic sensitivity, particularly regarding Earth's magnetic field strengths (25-65 µT). In most cases, isotopic substitution improved RP's magnetic sensitivity. Migratory behaviour may thus have been favoured in animals with certain isotopic compositions of cryptochrome.

Exploring the interaction sites in glucose and galactose using phenol as a probe

Sugars, together with amino acids and nucleobases, are the fundamental building blocks of a cell. They are involved in many fundamental processes and they especially play relevant roles as part of the immune system. The latter is connected to their ability to establish a collection of intermolecular interactions, depending on the position of their hydroxyl groups. Here we explore how the position of the OH in C4, the anomeric conformation and the nature substituent affect the interaction with phenol, which serves as a probe of the preferred site for the interaction. Using mass-resolved excitation spectroscopy and density functional calculations, we unravel the structure of the dimers and compare their conformation with those found for similar systems. The main conclusion is that the hydroxymethyl group has a very strong influence, guiding the whole aggregation process and that the position of the substituent in C4 has a stronger influence on the final structure of the dimer than the anomeric conformation.

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.

Setting up the HyDRA blind challenge for the microhydration of organic molecules

The procedure leading to the first HyDRA blind challenge for the prediction of water donor stretching vibrations in monohydrates of organic molecules is described. A training set of 10 monohydrates...

Rapid IRMPD (InfraRed multiple photon dissociation) analysis for glycomics

Infrared vibrational spectroscopy in the gas phase has emerged as a powerful tool to determine complex molecular structures with high precision. Among the different approaches IRMPD (InfraRed multiple photon dissociation), which requires the use of an intense pulsed tuneable laser in the InfraRed (IR) domain, has been broadly applied to the study of complex (bio)molecules. Recently, it also emerged as a highly relevant approach for analytical purposes especially in the field of glycomics in which structural analysis is still a tremendous challenge. This opens the perspective to develop new analytical tools allowing for the determination of molecular structures with atomic precision, and to address advanced questions in the field. However, IRMPD experiments require non commercial equipment or/and long acquisition time which limits the data output. Here we show that it is possible to improve the IRMPD performances by optimizing the combination between a linear ion trap mass spectrometer and a high repetition tuneable laser. Two orders of magnitude are gained with this approach compared to the usual experiments ultimately leading to a completely resolved spectrum acquired in less than one minute. These results open the way to many new applications in glycomics with the possibility to include IRMPD in complex analytical workflows.

Exploring the Molecular Conformation Space by Soft Molecule-Surface Collision

Biomolecules function by adopting multiple conformations. Such dynamics are governed by the conformation landscape whose study requires characterization of the ground and excited conformation states. Here, the conformational landscape of a molecule is sampled by exciting an initial gas-phase molecular conformer into diverse conformation states, using soft molecule-surface collision (0.5-5.0 eV). The resulting ground and excited molecular conformations, adsorbed on the surface, are imaged at the single-molecule level. This technique permits the exploration of oligosaccharide conformations, until now, limited by the high flexibility of oligosaccharides and ensemble-averaged analytical methods. As a model for cellulose, cellohexaose chains are observed in two conformational extremes, the typical "extended" chain and the atypical "coiled" chain-the latter identified as the gas-phase conformer preserved on the surface. Observing conformations between these two extremes reveals the physical properties of cellohexaose, behaving as a rigid ribbon that becomes flexible when twisted. The conformation space of any molecule that can be electrosprayed can now be explored.

Guanidinium/ammonium competition and proton transfer in the interaction of the amino acid arginine with the tetracarboxylic 18-crown-6 ionophore

The binding of arginine by the 18-crown-6 tetracarboxylic ionophore relies on extensive host–guest redistribution of electronic charge and proton transfer.

Machine learning molecular dynamics for the simulation of infrared spectra

Machine learning has emerged as an invaluable tool in many research areas. In the present work, we harness this power to predict highly accurate molecular infrared spectra with unprecedented computational efficiency. To account for vibrational anharmonic and dynamical effects - typically neglected by conventional quantum chemistry approaches - we base our machine learning strategy on ab initio molecular dynamics simulations. While these simulations are usually extremely time consuming even for small molecules, we overcome these limitations by leveraging the power of a variety of machine learning techniques, not only accelerating simulations by several orders of magnitude, but also greatly extending the size of systems that can be treated. To this end, we develop a molecular dipole moment model based on environment dependent neural network charges and combine it with the neural network potential approach of Behler and Parrinello. Contrary to the prevalent big data philosophy, we are able to obtain very accurate machine learning models for the prediction of infrared spectra based on only a few hundreds of electronic structure reference points. This is made possible through the use of molecular forces during neural network potential training and the introduction of a fully automated sampling scheme. We demonstrate the power of our machine learning approach by applying it to model the infrared spectra of a methanol molecule, n-alkanes containing up to 200 atoms and the protonated alanine tripeptide, which at the same time represents the first application of machine learning techniques to simulate the dynamics of a peptide. In all of these case studies we find an excellent agreement between the infrared spectra predicted via machine learning models and the respective theoretical and experimental spectra.

Isolated complexes of the amino acid arginine with polyether and polyamine macrocycles, the role of proton transfer

Protonated arginine interacts with 12-crown-4 through the guanidinium side group. In the complex with the N-substituted analog cyclen, the dominant conformation is the result of the proton transfer from the carboxylic acid group of the amino acid to the macrocycle.

Isolated alkali cation complexes of the antibiotic ionophore nonactin: correlation with crystalline structures

The stabilization of the nonactin–Na⁺ complex in a S₄ or C₂ conformation constitutes a challenging benchmark for experimental and modelling methods.

Retention of Native Protein Structures in the Absence of Solvent: A Coupled Ion Mobility and Spectroscopic Study

Can the structures of small to medium‐sized proteins be conserved after transfer from the solution phase to the gas phase? A large number of studies have been devoted to this topic, however the answer has not been unambiguously determined to date. A clarification of this problem is important since it would allow very sensitive native mass spectrometry techniques to be used to address problems relevant to structural biology. A combination of ion‐mobility mass spectrometry with infrared spectroscopy was used to investigate the secondary and tertiary structure of proteins carefully transferred from solution to the gas phase. The two proteins investigated are myoglobin and β‐lactoglobulin, which are prototypical examples of helical and β‐sheet proteins, respectively. The results show that for low charge states under gentle conditions, aspects of the native secondary and tertiary structure can be conserved.

Furanosic forms of sugars: conformational equilibrium of methyl β-d-ribofuranoside

The rotational spectrum of the isolated ribofuranoside unit, the biologically active sugar form, revealed two structures with a similar puckering to the RNA units in crystals.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Distinguishing isobaric phosphated and sulfated carbohydrates by coupling of mass spectrometry with gas phase vibrational spectroscopy

An original application of the coupling of mass spectrometry with vibrational spectroscopy, used for the first time to discriminate isobaric bioactive saccharides with sulfate and phosphate functional modifications, is presented. Whereas their nominal masses and fragmentation patterns are undifferentiated by sole mass spectrometry, their distinctive OH stretching modes at 3595 cm(-1) and 3666 cm(-1), respectively, provide a reliable spectroscopic diagnostic for distinguishing their sulfate or phosphate functionalization. A detailed analysis of the 6-sulfated and 6-phosphated d-glucosamine conformations is presented, together with theoretical scaled harmonic spectra and anharmonic spectra (VPT2 and DFT-based molecular dynamics simulations). Strong anharmonic effects are observed in the case of the phosphated species, resulting in a dramatic enhancement of its phosphate diagnostic mode.

Isolated neutral peptides

This chapter examines the structural characterisation of isolated neutral amino-acids and peptides. After a presentation of the experimental and theoretical state-of-the-art in the field, a review of the major structures and shaping interactions is presented. Special focus is made on conformationally-resolved studies which enable one to go beyond simple structural characterisation; probing flexibility and excited-state photophysics are given as examples of promising future directions.

Intra-residue interactions in proteins: interplay between serine or cysteine side chains and backbone conformations, revealed by laser spectroscopy of isolated model peptides

Intra-residue interactions play an important role in proteins by influencing local folding of the backbone. Taking advantage of the capability of gas phase experiments to provide relevant information on the intrinsic H-bonding pattern of isolated peptide chains, the intra-residue interactions of serine and cysteine residues, i.e., OH/SH···OC(i) C6 and NH(i···)O/S C5 interactions in Ser/Cys residues, are probed by laser spectroscopy of isolated peptides. The strength of these local side chain-main chain interactions, elegantly documented from their IR spectral features for well-defined conformations of the main chain, demonstrates that a subtle competition exists between the two types of intra-residue bond: the C6 H-bond is the major interaction with Ser, in contrast to Cys where C5 interaction takes over. The restricted number of conformers observed in the gas phase experiment with Ser compared to Cys (where both extended and folded forms are observed) also suggests a significant mediation role of these intra-residue interactions on the competition between the several main chain folding patterns.

High energy conformers of M(+)(APE)(H2O)(0-1)Ar(0-1) clusters revealed by combined IR-PD and DFT-MD anharmonic vibrational spectroscopy

IR-PD vibrational spectroscopy and DFT-based molecular dynamics simulations are combined in order to unravel the structures of M(+)(APE)(H2O)0-1 ionic clusters (M = Na, K), where APE (2-amino-1-phenyl ethanol) is commonly used as an analogue for the noradrenaline neurotransmitter. The strength of the synergy between experiments and simulations presented here is that DFT-MD provides anharmonic vibrational spectra that unambiguously help assign the ionic clusters structures. Depending on the interacting cation, we have found that the lowest energy conformers of K(+)(APE)(H2O)0-1 clusters are formed, while the lowest energy conformers of Na(+)(APE)(H2O)0-1 clusters can only be observed through water loss channel (i.e. without argon tagged to the clusters). Trapping of higher energy conformers is observed when the argon loss channel is recorded in the experiment. This has been rationalized by transition state energies. The dynamical anharmonic vibrational spectra unambiguously provide the prominent OH stretch due to the OH···NH2 H-bond, within 10 cm(-1) of the experiment, hence reproducing the 240-300 cm(-1) red-shift (depending on the interacting cation) from bare neutral APE. When this H-bond is not present, the dynamical anharmonic spectra provide the water O-H stretches as well as the rotational motion of the water molecule at finite temperature, as observed in the experiment.

Many-body van der Waals interactions in molecules and condensed matter

This work reviews the increasing evidence that many-body van der Waals (vdW) or dispersion interactions play a crucial role in the structure, stability and function of a wide variety of systems in biology, chemistry and physics. Starting with the exact expression for the electron correlation energy provided by the adiabatic connection fluctuation-dissipation theorem, we derive both pairwise and many-body interatomic methods for computing the long-range dispersion energy by considering a model system of coupled quantum harmonic oscillators within the random-phase approximation. By coupling this approach to density functional theory, the resulting many-body dispersion (MBD) method provides an accurate and efficient scheme for computing the frequency-dependent polarizability and many-body vdW energy in molecules and materials with a finite electronic gap. A select collection of applications are presented that ascertain the fundamental importance of these non-bonded interactions across the spectrum of intermolecular (the S22 and S66 benchmark databases), intramolecular (conformational energies of alanine tetrapeptide) and supramolecular (binding energy of the 'buckyball catcher') complexes, as well as molecular crystals (cohesive energies in oligoacenes). These applications demonstrate that electrodynamic response screening and beyond-pairwise many-body vdW interactions--both captured at the MBD level of theory--play a quantitative, and sometimes even qualitative, role in describing the properties considered herein. This work is then concluded with an in-depth discussion of the challenges that remain in the future development of reliable (accurate and efficient) methods for treating many-body vdW interactions in complex materials and provides a roadmap for navigating many of the research avenues that are yet to be explored.

Photofragment ionization-loss stimulated Raman spectroscopy of a hydrated neurotransmitter: 2-phenylethylamine–water

The first use of photofragment ionization-loss stimulated Raman spectroscopy, in a broad range, for structural studies of clusters is reported.

Carbohydrates

Although carbohydrates represent one of the most important families of biomolecules, they remain under-studied in comparison to the other biomolecular families (peptides, nucleobases). Beyond their best-known function of energy source in living systems, they act as mediator of molecular recognition processes, carrying molecular information in the so-called "sugar code," just to name one of their countless functions. Owing to their high conformational flexibility, they encode extremely rich information conveyed via the non-covalent hydrogen bonds within the carbohydrate and with other biomolecular assemblies, such as peptide subunits of proteins. Over the last decade there has been tremendous progress in the study of the conformational preferences of neutral oligosaccharides, and of the interactions between carbohydrates and various molecular partners (water, aromatic models, and peptide models), using vibrational spectroscopy as a sensitive probe. In parallel, other spectroscopic techniques have recently become available to the study of carbohydrates in the gas phase (microwave spectroscopy, IRMPD on charged species).

'Naked' and hydrated conformers of the conserved core pentasaccharide of N-linked glycoproteins and its building blocks

N-glycosylation of eukaryotic proteins is widespread and vital to survival. The pentasaccharide unit -Man3GlcNAc2- lies at the protein-junction core of all oligosaccharides attached to asparagine side chains during this process. Although its absolute conservation implies an indispensable role, associated perhaps with its structure, its unbiased conformation and the potential modulating role of solvation are unknown; both have now been explored through a combination of synthesis, laser spectroscopy, and computation. The proximal -GlcNAc-GlcNAc- unit acts as a rigid rod, while the central, and unusual, -Man-β-1,4-GlcNAc- linkage is more flexible and is modulated by the distal Man-α-1,3- and Man-α-1,6- branching units. Solvation stiffens the 'rod' but leaves the distal residues flexible, through a β-Man pivot, ensuring anchored projection from the protein shell while allowing flexible interaction of the distal portion of N-glycosylation with bulk water and biomolecular assemblies.