Rotational spectrum of tryptophan

DFT Meets Wave-Function Methods for Accurate Structures and Rotational Constants of Histidine, Tryptophan, and Proline

A new computational strategy has been applied to the conformational and spectroscopic properties in the gas phase of amino acids with very distinctive features, ranging from different tautomeric forms (histidine) to ring puckering (proline), and heteroaromatic structures with non-equivalent rings (tryptophan). The integration of modern double-hybrid functionals and wave-function composite methods has allowed us to obtain accurate results for a large panel of conformers with reasonable computer times. The remarkable agreement between computations and microwave experiments allows an unbiased interpretation of the latter in terms of stereoelectronic effects.

Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods

A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.

Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy

The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.

Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks

The interplay of high-resolution rotational spectroscopy and quantum-chemical computations plays an invaluable role in the investigation of biomolecule building blocks in the gas phase. However, quantum-chemical methods suffer from unfavorable scaling with the dimension of the system under consideration. While a complete characterization of flexible systems requires an elaborate multi-step strategy, in this work, we demonstrate that the accuracy obtained by quantum-chemical composite approaches in the prediction of rotational spectroscopy parameters can be approached by a model based on density functional theory. Glycine and serine are employed to demonstrate that, despite its limited cost, such a model is able to predict rotational constants with an accuracy of 0.3% or better, thus paving the way toward the accurate characterization of larger flexible building blocks of biomolecules.

The Nicotinic Agonist Cytisine: The Role of the NH···N Interaction

We report a detailed structural study of cytisine, an alkaloid used to help with smoking cessation, looking forward to unveiling its role as a nicotinic agonist. High-resolution rotational spectroscopy has allowed us to characterize two different conformers exhibiting axial and equatorial arrangements of the piperidinic NH group. Unexpectedly, the axial form has been found as the predominant configuration, in contrast to that observed for related molecules, such as piperidine. This anomalous behavior has been justified in terms of an intramolecular NH···N hydrogen bond. Moreover, this interaction justifies the overstabilization of the axial conformer over the equatorial one and is crucial for the mechanism of action of cytisine over the nicotinic receptor, further rationalizing its behavior as a nicotinic agonist.

Unleashing the shape of L-DOPA at last

Herein, we report the first rotational study of neutral L-DOPA, an extensively used supramolecular synthon and an amino acid precursor of the neurotransmitters dopamine, norepinephrine (noradrenaline), and epinephrine (adrenaline) using broadband and narrowband Fourier transform microwave spectroscopies coupled with a laser ablation vaporization system. The spectroscopic parameters derived from the analysis of the rotational spectrum conclusively identify the existence of four distinct conformers of L-DOPA in the supersonic jet, further rejecting the previously reported catechol ring-induced conformational restriction. The analysis of the ¹⁴N nuclear quadrupole coupling hyperfine structure further revealed the orientation of the N-bearing functional group, proving the existence of stabilizing N-H⋯π interactions for the observed structures.

An Innovative Approach for the Generation of Species of the Interstellar Medium

The large amount of unstable species in the realm of interstellar chemistry drives an urgent need to develop efficient methods for the in situ generations of molecules that enable their spectroscopic characterizations. Such laboratory experiments are fundamental to decode the molecular universe by matching the interstellar and terrestrial spectra. We propose an approach based on laser ablation of nonvolatile solid organic precursors. The generated chemical species are cooled in a supersonic expansion and probed by high‐resolution microwave spectroscopy. We present a proof of concept through a simultaneous formation of interstellar compounds and the first generation of aminocyanoacetylene using diaminomaleonitrile as a prototypical precursor. With this micro‐laboratory, we open the door to generation of unsuspected species using precursors not typically accessible to traditional techniques such as electric discharge and pyrolysis.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Vibrational Predissociation Spectroscopy of Cold Protonated Tryptophan with Different Messenger Tags

Vibrational spectra of protonated tryptophan were recorded by predissociation of H₂ messenger tags using cryogenic ion traps. We explore the issue of messenger induced spectral changes by solvating TrpH⁺(H₂) _n with n = 1-5 to obtain single photon vibrational spectra of TrpH⁺ and of its partly deuterated isotopomer in the spectral region of 800-4400 cm^-1. Depending on the number of messenger molecules, the spectra of several conformational isomers associated with multiple H₂ binding locations along with two natural conformations of TrpH⁺ were found using the two photon MS³IR² conformational hole burning method. Most probable messenger positions were established by comparison with predictions from DFT calculations on various candidate structures. Mechanical anharmonicity effects associated with the charged amino group were modeled by Born-Oppenheimer ab initio molecular dynamics. The spectra of TrpH⁺(H₂O) _m=1,2, recorded by infrared multiphoton dissociation (IRMPD), reveal broad features in the NH stretching region of the NH₃⁺ group, indicating strong hydrogen bonding in acceptor-donor configuration with the benzene ring for the first water molecule, while the second water appears to attach to a less strongly perturbing site, yielding unique transitions associated with the free OH stretching fundamentals. We discuss the structural deformations induced by the water molecules and compare our results to recent experiments on similar hydrated cationic systems.

Conformational study of L-methionine and L-cysteine derivatives through quantum chemical calculations and 3JHH coupling constant analyses

The understanding of the conformational behavior of amino acids and their derivatives is a challenging task. Here, the conformational analysis of esterified and N-acetylated derivatives of L-methionine and L-cysteine using a combination of 1H NMR and electronic structure calculations is reported. The geometries and energies of the most stable conformers in isolated phase and taking into account the implicit solvent effects, according to the integral equation formalism polarizable continuum model (IEF-PCM), were obtained at the ωB97X-D/aug-cc-pVTZ level. The conformational preferences of the compounds in solution were also determined from experimental and theoretical 3JHH coupling constants analysis in different aprotic solvents. The results showed that the conformational stability of the esterified derivatives is not very sensitive to solvent effects, whereas the conformational equilibrium of the N-acetylated derivatives changes in the presence of solvent. According to the natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) and noncovalent interactions (NCI) methodologies, the conformational preferences for the compounds are not dictated by intramolecular hydrogen bonding, but by a joint contribution of hyperconjugative and steric effects.

Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases

The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.

Generation and structural characterization of aluminum cyanoacetylide

Combined spectroscopy measurements and theoretical calculations bring to light a first investigation of a metallic cyanoacetylide, AlC3N, using laser ablation molecular beam Fourier transform microwave spectroscopy. This molecule was synthesized in a supersonic expansion by the reaction of aluminum vapour with C3N, produced from solid aluminum rods and BrCCCN in a newly constructed ablation-heating nozzle device. A set of accurate rotational and (27)Al and (14)N nuclear quadrupole coupling constants have been determined from the analysis of the rotational spectrum and compared with those predicted in a high-level ab initio study, conducting to the assignment of the observed species to linear AlCCCN. We have searched for this species towards the carbon-rich evolved star IRC + 10216 but only an upper limit to its abundance has been obtained.