A Decapeptide Hydrated by Two Waters: Conformers Determined by Theory and Validated by Cold Ion Spectroscopy

Performance of Effective Harmonic Oscillator Approach for the Calculations of Vibrational Transition Energies of Large Molecules

The accuracy and performance of the effective harmonic oscillator approximation for the description of anharmonic vibrational structure calculations are tested for large molecular systems and compared with experimental values along with vibrational self-consistent field and second-order perturbation theories. The effective harmonic oscillator approach is an effective single-particle approximation where the variational parameters are the centroids and widths of the multidimensional Gaussian product functions posited as the vibrational wave functions. A comprehensive calculation for 849 transitions that include the fundamentals, two and three quanta overtone transitions, and several combination bands of three polyaromatic hydrocarbons and one DNA nucleobase with a total of 231 normal modes are assessed. A comparison of EHO results with the experimental values is done for the polyaromatic hydrocarbons, and a close agreement is found between the two results. It also offers anharmonic eigenstates and eigenfunctions that are nearly identical with vibrational self-consistent field theory. An extensive analysis on the resultant wave functions of the excited states is performed. The overall root-mean-square deviation (RMSD) between these two methods for 849 transitions understudy is only about 8.3 cm-1, suggesting the effective harmonic oscillator as a viable alternative for the reliable calculations of transition energies of large molecular systems.

Ion Spectroscopy Reveals Structural Difference for Proteins Microhydrated by Retention and Condensation of Water

Protein ubiquitin in its +7 charge state microhydrated by 5 and 10 water molecules has been interrogated in the gas phase by cold ion UV/IR spectroscopy. The complexes were formed either by condensing water onto the unfolded bare proteins in a temperature-controlled ion trap or by incomplete dehydration of the folded proteins. In the case of cryogenic condensation, the UV spectra of the complexes exhibit a resolved vibrational structure, which looks similar to the spectrum of bare unfolded ubiquitin. The spectra become, however, broad-band with no structure when complexes of the same size are produced by incomplete dehydration under soft conditions of electrospray ionization. We attribute this spectroscopic dissimilarity to the structural difference of the protein: condensing a few water molecules cannot refold the gas-phase structure of the bare ubiquitin, while the retained water preserves its solution-like folded motif through evaporative cooling. This assessment is firmly confirmed by IR spectroscopy, which reveals the presence of free NH and carboxylic OH stretching vibrations only in the complexes with condensed water.

Tracking local and global structural changes in a protein by cold ion spectroscopy

Characterization of native structures of proteins in the gas phase remains challenging due to the unpredictable conformational changes the molecules undergo during desolvation and ionization. We spectroscopically studied cryogenically cooled protonated protein ubiquitin and its microhydrated complexes prepared in the gas phase in a range of charge states under different ionization conditions. The UV spectra appear vibrationally resolved for the unfolded protein, but become redshifted and smooth for the native-like structures of ubiquitin. This spectroscopic change results from the H-bonding of the hydroxyl of Tyr to the amide group of Glu-51 in the compact structures; the minimum length of this bond was estimated to be ∼1.7 Å. IR spectroscopy reflects the global structural change by observing redshifts of free NH/OH-stretch vibrational transitions. Evaporative cooling of microhydrated complexes of ubiquitin keeps the protein chilly during ionization, enabling native-like conformers with up to eight protons to survive in the gas phase.

Performance of Vibrational Self-Consistent Field Theory for Accurate Potential Energy Surfaces: Fundamentals, Excited States, and Intensities

The performance of vibrational structure calculations beyond harmonic approximation in the framework of the vibrational self-consistent field method with second-order perturbation corrections (VSCF-PT2) is investigated in conjunction with very accurate potential energy surfaces (PESs) given by various coupled-cluster electronic structure theories. The quality of anharmonic calculations depends on the accuracy of the underlying multidimensional PES obtained from its functional form, which is given by the level of electronic structure theory. Two such highest levels of typical coupled-cluster electronic structure methods, CCSD and the ″gold standard″ CCSD(T), along with their variants such as CCD, CR-CCL (completely renormalized CR-CC(2,3) approach), and CCSD(TQ) are tested for the construction of accurate anharmonic potentials without any fitting or ad hoc scaling and using cc-pVTZ basis sets. The accuracy of VSCF-PT2 theory in comparison to experimental values is tested for a series of 16 molecules with 135 fundamental bands, 64 overtones, and combination bands and also for 39 intensities. It is found that CCD and CCSD bind the potential tighter than CCSD(T) and the computed VSCF-PT2 transitions are more blue-shifted showing higher deviation from the experiment. In general, VSCF-PT2 results computed at the CCSD(T) potential offer a good cost/accuracy ratio, with the mean absolute deviation and the mean absolute percentage error with the experiment being ∼16 cm^-1 and 1.38, respectively, for fundamentals. Additionally, while the CR-CCL and CCSD(TQ) methods offer similar levels of accuracies as compared to CCSD(T), the former offers a better accuracy/cost ratio than the latter and is a suitable alternative to CCSD(T).

My Trajectory in Molecular Reaction Dynamics and Spectroscopy

This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!

Microhydration of Biomolecules: Revealing the Native Structures by Cold Ion IR Spectroscopy

The native-like structures of protonated glycine and peptide Gly₃H⁺ were elucidated using cold ion IR spectroscopy of these biomolecules hydrated by a controlled number of water molecules. The complexes were generated directly from an aqueous solution using gentle electrospray ionization. Already with a single retained water molecule, GlyH⁺ exhibits the native-like structure characterized by a lack of intramolecular hydrogen bonds. We use our spectra to calibrate the available data for the same complexes, which are produced by cryogenic condensation of water onto the gas-phase glycine. In some conformers of these complexes, GlyH⁺ adopts the native-like structure, while in the others, it remains "kinetically" trapped in the intrinsic state. Upon condensation of 4-5 water molecules, the embedded amino acid fully adopts its native-like structure. Similarly, condensation of one water molecule onto the tripeptide is insufficient to fully eliminate its kinetically trapped intrinsic states.

Dual Basis Approach for Ab Initio Anharmonic Calculations of Vibrational Spectroscopy: Application to Microsolvated Biomolecules

A dual electronic basis set approach is introduced for more efficient but accurate calculations of the anharmonic vibrational spectra in the framework of the vibrational self-consistent field (VSCF) theory. In this approach, an accurate basis set is used to compute the vibrational spectra at the harmonic level. The results are used to scale the potential surface from a more modest but much more efficient basis set. The scaling is such that at the harmonic level the new, scaled potential agrees with one of the accurate basis sets. The approach is tested in the application of the microsolvated, protected amino acid Ac-Phe-OMe, using the scaled anharmonic hybrid potential in the VSCF and VSCF-PT2 algorithms. The hybrid potential method yields results that are in good accord with the experiment and very close to those obtained in calculations with the high-level, very costly potential from the large basis set. At the same time, the hybrid potential calculations are considerably less expensive. The results of the hybrid calculations are much more accurate than those computed from the potential surface corresponding to the modest basis set. The results are very encouraging for using the hybrid potential method for inexpensive yet sufficiently accurate anharmonic calculations for the spectra of large biomolecules.

Computational Insights into Compaction of Gas-Phase Protein and Protein Complex Ions in Native Ion Mobility-Mass Spectrometry

Native ion mobility-mass spectrometry (IM-MS) is a rapidly growing field for studying the composition and structure of biomolecules and biomolecular complexes using gas-phase methods. Typically, ions are formed in native IM-MS using gentle nanoelectrospray ionization conditions, which in many cases can preserve condensed-phase stoichiometry. Although much evidence shows that large-scale condensed-phase structure, such as quaternary structure and topology, can also be preserved, it is less clear to what extent smaller-scale structure is preserved in native IM-MS. This review surveys computational and experimental efforts aimed at characterizing compaction and structural rearrangements of protein and protein complex ions upon transfer to the gas phase. A brief summary of gas-phase compaction results from molecular dynamics simulations using multiple common force fields and a wide variety of protein ions is presented and compared to literature IM-MS data.

Efficiently Calculating Anharmonic Frequencies of Molecular Vibration by Molecular Dynamics Trajectory Analysis

Two efficient methods, the Eckart frame algorithm and the multiorder derivative algorithm, for vibrational frequency calculation directly based on the raw data of atomic trajectory from the state-of-the-art first-principles molecular dynamics simulation are presented. The Eckart frame approach is robust to retrieve the full set of anharmonic fundamental frequencies of any molecule from the atomic trajectory for a sufficiently long molecular dynamics simulation at a temperature close to 0 K. In addition to the fundamental vibrational frequencies, the multiorder derivative approach is universal for the calculations of vibrational frequencies based on the molecular dynamics result in a wide range of temperatures. The accuracy, efficiency, and applicability of these two methods are demonstrated through several successful examples in calculating the anharmonic fundamental vibrational frequencies of methane, ethylene, water, and cyclobutadiene.

Intrinsic structure of pentapeptide Leu-enkephalin: geometry optimization and validation by comparison of VSCF-PT2 calculations with cold ion spectroscopy

The intrinsic structure of an opioid peptide [Ala2, Leu5]-leucine enkephalin (ALE) has been investigated using first-principles based vibrational self-consistent field (VSCF) theory and cold ion spectroscopy. IR-UV double resonance spectroscopy revealed the presence of only one highly abundant conformer of the singly protonated ALE, isolated and cryogenically cooled in the gas phase. High-level quantum mechanical calculations of electronic structures in conjunction with a systematic conformational search allowed for finding a few low-energy candidate structures. In order to identify the observed structure, we computed vibrational spectra of the candidate structures and employed the theory at the semi-empirically scaled harmonic level and at the first-principles based anharmonic VSCF levels. The best match between the calculated "anharmonic" and the measured spectra appeared, indeed, for the most stable candidate. An average of two spectra calculated with different quantum mechanical potentials is proposed for the best match with experiment. The match thus validates the calculated intrinsic structure of ALE and demonstrates the predictive power of first-principles theory for solving structures of such large molecules.

Conformational investigation of the structure-activity relationship of GdFFD and its analogues on an achatin-like neuropeptide receptor of Aplysia californica involved in the feeding circuit

Proteins and peptides in nature are almost exclusively made from l-amino acids, and this is even more absolute in the metazoan. With the advent of modern bioanalytical techniques, however, previously unappreciated roles for d-amino acids in biological processes have been revealed. Over 30 d-amino acid containing peptides (DAACPs) have been discovered in animals where at least one l-residue has been isomerized to the d-form via an enzyme-catalyzed process. In Aplysia californica, GdFFD and GdYFD (the lower-case letter "d" indicates a d-amino acid residue) modulate the feeding behavior by activating the Aplysia achatin-like neuropeptide receptor (apALNR). However, little is known about how the three-dimensional conformation of DAACPs influences activity at the receptor, and the role that d-residues play in these peptide conformations. Here, we use a combination of computational modeling, drift-tube ion-mobility mass spectrometry, and receptor activation assays to create a simple model that predicts bioactivities for a series of GdFFD analogs. Our results suggest that the active conformations of GdFFD and GdYFD are similar to their lowest energy conformations in solution. Our model helps connect the predicted structures of GdFFD analogs to their activities, and highlights a steric effect on peptide activity at position 1 on the GdFFD receptor apALNR. Overall, these methods allow us to understand ligand-receptor interactions in the absence of high-resolution structural data.

Infrared ion spectroscopy: an analytical tool for the study of metabolites

Vibrational ion spectroscopy techniques coupled with mass spectrometry are applied to standard metabolites as a proof-of-principle demonstration for the structural identification of unknown metabolites. The traditional room temperature infrared multiple photon dissociation (IRMPD) spectroscopy technique is shown to differentiate chemical moieties in isobaric and isomeric variants. These results are compared to infrared spectra of cryogenically cooled analyte ions, showing enhanced spectral resolution, and thus also improved differentiation between closely related molecules, such as isomers. The cryogenic spectroscopy is effected in a recently developed mass-selective cryogenic linear ion trap, which is capable of high sensitivity and the ability to measure the IR spectra of multiple analytes simultaneously.