Isomers and conformational barriers of gas-phase nicotine, nornicotine, and their protonated forms

The Gas Phase Protonation Sites of Six Naturally Occurring Nicotinoids

The gas phase protonation sites of six naturally occurring nicotinoids, namely nicotine (NIC), nornicotine (NOR), anabasine (ANB), anatabine (ANT), cotinine (COT), and myosmine (MYO), consisting of a common Pyridine and differing non-Pyridine rings, have been determined for the first time at the physiological temperature from cryogenic ion trap infrared spectroscopy and electronic structure calculations. The protonation site on either of these two rings is related to the nicotinoid's biological activity. At room temperature, NIC is a mixture of Pyridine and Pyrrolidine (non-Pyridine) protomers, whereas NOR, ANB, ANT, and COT are pure Pyridine protomers and finally MYO is mostly a Pyroline (non-Pyridine) protomer. The nearly planar structure of MYO-H+, induced by the presence of a conjugated π system and confirmed from calculations and the UV absorption spectra, breaks from the trends observed for NIC, NOR, and ANB, since its structure is drastically different from the structures of the other nicotinoids.

Affinity of Nicotinoids to a Model Nicotinic Acetylcholine Receptor (nAChR) Binding Pocket in the Human Brain

The binding affinity of nicotinoids to the binding residues of the α4β2 variant of the nicotinic acetylcholine receptor (nAChR) was identified as a strong predictor of the nicotinoid's addictive character. Using ab initio calculations for model binding pockets of increasing size composed of 3, 6, and 14 amino acids (3AA, 6AA, and 14AA) that are derived from the crystal structure, the differences in binding affinity of 6 nicotinoids, namely, nicotine (NIC), nornicotine (NOR), anabasine (ANB), anatabine (ANT), myosmine (MYO), and cotinine (COT) were correlated to their previously reported doses required for increases in intracranial self-stimulation (ICSS) thresholds, a metric for their addictive function. By employing the many-body decomposition, the differences in the binding affinities of the various nicotinoids could be attributed mainly to the proton exchange energy between the pyridine and non-pyridine rings of the nicotinoids and the interactions between them and a handful of proximal amino acids, namely Trp156, Trpβ57, Tyr100, and Tyr204. Interactions between the guest nicotinoid and the amino acids of the binding pocket were found to be mainly classical in nature, except for those between the nicotinoid and Trp156. The larger pockets were found to model binding structures more accurately and predicted the addictive character of all nicotinoids, while smaller models, which are more computationally feasible, would only predict the addictive character of nicotinoids that are similar to nicotine. The present study identifies the binding affinity of the guest nicotinoid to the host binding pocket as a strong descriptor of the nicotinoid's addiction potential, and as such it can be employed as a fast-screening technique for the potential addiction of nicotine analogs.

Conformer Weighting and Differently Sized Cluster Weighting for Nicotine and Its Phosphorus Derivatives

Weighting methods applied to systems with many conformers have been broadly employed to calculate thermodynamic properties, structural characteristics, and populations. To better understand and test the sensitivity of conventional weighting methods, the conformational distributions of nicotine and its phosphorus-substituted derivatives are investigated. The weighting schemes used for this are all based on Boltzmann statistics. Classical Boltzmann factors based on the electronic energy and the Gibbs free energy are calculated at different quantum chemical levels of theory and compared to cluster weights obtained by the quantum cluster equilibrium method. Furthermore, the influence of the modified rigid-rotor-harmonic-oscillator (mRRHO) approximation on the cluster weights is investigated. The substitution of the nitrogen atom in the methylpyrrolidine ring by a phosphorus atom results in more monomer conformers and clusters being populated. The conformational distribution of the monomers remained stable at different levels of theory and weighting methods. However, going to dimers and trimers, we observe a significant influence of the level of theory, weighting method, and mRRHO cutoff on the populations of these clusters. We show that mRRHO cutoff values of 50 and 100 cm^-1 yield similar results, which is why 50 cm^-1 is recommended as a robust choice. Furthermore, we observe that the global minimum for ΔE₀ and ΔG varies in a few cases and that the global minimum is not always the dominantly occupied structure.

Rapid Gas-Phase Autoxidation of Nicotine in the Atmosphere

Nicotine is the most abundant alkaloid chemical in smoke emission. In this work, we investigated the gas-phase oxidation mechanism of nicotine initiated by its reactions with the OH radical and ozone. Both initiation reactions start dominantly by hydrogen atom abstractions from the C1, C3, and -CH₃ groups of the methylpyrrolidinyl group and form radicals nicotinyl-1, nicotinyl-3, and nicotinyl-6, respectively. The nicotinyl radicals would recombine rapidly with O₂, forming RO₂ with rapid intramolecular hydrogen-atom transfers (HATs) with rate coefficients from 4 s^-1 to greater than 10⁴ s^-1. The rapid HATs in peroxy radicals suggest rapid autoxidation of nicotine in the gas phase. Formation of HCNO and HC(O)NH₂, being observed in previous studies, arises likely from secondary reactions or photolysis of intermediate products.

Structure of Gas Phase Monohydrated Nicotine: Implications for Nicotine's Native Structure in the Acetylcholine Binding Protein

We report a joint experimental-theoretical study of the never reported before structure and infrared spectra of gas phase monohydrated nicotine (NIC) and nornicotine (NOR) and use them to assign their protonation sites. NIC's biological activity is strongly affected by its protonation site, namely, the pyrrolidine (Pyrro-NICH⁺, anticipated active form) and pyridine (Pyri-NICH⁺) forms; however, these have yet to be directly experimentally determined in either the nicotinic acetylcholine receptor (nAChR, no water present) or the acetylcholine-binding protein (AChBP, a single water molecule is present) but can only be inferred to be Pyrro-NICH⁺ from the intermolecular distance to the neighboring residues (i.e., tryptophan). Our temperature-controlled double ion trap infrared spectroscopic experiments assisted by the collisional stripping method and high-level theoretical calculations yield the protonation ratio of Pyri:Pyrro = 8:2 at 240 K for the gas phase NICH⁺···(H₂O) complex, which resembles the molecular cluster present in the AChBP. Therefore, a single water molecule in the gas phase enhances this ratio in NICH⁺ relative to the 3:2 for the nonhydrated gas phase NICH⁺ in a trend that contrasts with the almost exclusive presence of Pyrro-NICH⁺ in aqueous solution. In contrast, the Pyri-NORH⁺ protomer is exclusively observed, a fact that may correlate with its weaker biological activity.

Gas phase protonated nicotine is a mixture of pyridine- and pyrrolidine-protonated conformers: implications for its native structure in the nicotinic acetylcholine receptor

The infrared (IR) spectra of gas phase protonated nicotine has been measured in the never-before probed N-H "fingerprint region" (3200-3500 cm^-1). The protonated molecules generated by an electrospray source are thermalized in the first ion trap with water vapor and He gas at a pre-determined temperature prior to being probed by IR spectroscopy in the second ion trap at 4 K. The IR spectra exhibit two N-H stretching bands which are assigned to the pyridine and pyrrolidine protomers with the aid of high-level electronic structure calculations. This finding is in sharp contrast to previous spectroscopic studies that suggested a single population of the pyridine protomer. The relative populations of the two protomers vary by changing the temperature of the thermalizing trap from 180-300 K. The relative conformer populations at 240 K and 300 K are well reproduced by the theoretical calculations, unequivocally determining that gas phase nicotine is a 3 : 2 mixture of both pyridine and pyrrolidine protomers at room temperature. The thermalizing anhydrous vapor does not result in any population change. It rather demonstrates the catalytic role of water in achieving equilibrium between the two protomers. The combination of IR spectroscopy and electronic structure calculations establish the small energy difference between the pyridine and pyrrolidine protomers in nicotine. One of the gas phase nicotine pyrrolidine protomers has the closest conformational resemblance among all low-lying energy isomers with the X-ray structure of nicotine in the nicotinic acetylcholine receptor (nAChR).

Early Chemistry of Nicotine Degradation in Heat-Not-Burn Smoking Devices and Conventional Cigarettes: Implications for Users and Second- and Third-Hand Smokers

Nicotine exposure results in health risks not only for smokers but also for second- and third-hand smokers. Unraveling nicotine's degradation mechanism and the harmful chemicals that are produced under different conditions is vital to assess exposure risks. We performed a theoretical study to describe the early chemistry of nicotine degradation by investigating two important reactions that nicotine can undergo: hydrogen abstraction by hydroxyl radicals and unimolecular dissociation. The former contributes to the control of the degradation mechanism below 800 K due to a non-Arrhenius kinetics, which implies an enhancement of reactivity as temperature decreases. The latter becomes important at higher temperatures due to its larger activation energy. This change in the degradation mechanism is expected to affect the composition of vapors inhaled by smokers and room occupants. Conventional cigarettes, which operate at temperatures higher than 1000 K, are more prone to yield harmful pyridinyl radicals via nicotine dissociation, while nicotine in electronic cigarettes and vaporizers, with operating temperatures below 600 K, will be more likely degraded by hydroxyl radicals, resulting in a vapor with a different composition. Although low-temperature nicotine delivery devices have been claimed to be less harmful due to their nonburning operating conditions, the non-Arrhenius kinetics that we observed for the degradation mechanism below 873 K suggests that nicotine degradation may be more rapidly initiated as temperature is reduced, indicating that these devices may be more harmful than it is commonly assumed.

Computational studies of biologically active alkaloids of plant origin: an overview

Computational studies nowadays constitute a crucial source of information for drug development, because they provide information on many molecular properties and also enable predictions of the properties of not-yet-synthesized compounds. Alkaloids are a vast group of natural products exhibiting a variety of biological activities, many of which are interesting for drug development. On the other hand, computational studies of biologically active alkaloids have so far mostly focused on few particularly relevant or “popular” molecules, such as quinine, caffeine, or cocaine, with only few works on the other molecules. The present work offers an overview of existing computational studies on alkaloid molecules, from the earliest ones to the most recent, and considering all the theoretical approaches with which studies have been performed (both quantum mechanics and molecular dynamics). The considered studies are grouped according to their objectives and outcomes, such as conformational analysis of alkaloid molecules, effects of selected solvents on their properties, docking studies aimed at better understanding of the interactions between alkaloid molecules and biological targets, studies focusing on structure activity relationships, and computational studies performed to confirm experimental results. It is concluded that it would be important that computational studies on many other alkaloid molecules are performed and their results made available, covering their different classes as well as the variety of their biological activities, to attain better understanding of the properties not only of individual molecules, but also of groups of related molecules and of the overall alkaloids family.

Conformer Specific Ultraviolet and Infrared Detection of Nicotine in the Vapor Phase

A gas-phase electronic spectrum of nicotine in a supersonic expansion has been recorded using two-color resonant two-photon ionization spectroscopy. Efficient photoionization was achievable only via the pyridine chromophore owing to poor Franck-Condon overlap in the N-methylpyrrolidine moiety. Two conformers of nicotine have been characterized and assigned by infrared-ultraviolet (IR-UV) ion depletion and IR-UV hole-burning spectroscopy, in combination with quantum chemical techniques. Trans-A with nitrogen atoms further apart is more stable by 2 kJ mol^-1 and the most populated conformer in the supersonic jet, owing this stability to a stronger inter-ring CH···N hydrogen bond than the trans-B counterpart.

Selecting and identifying gas-phase protonation isomers of nicotineH+ using combined laser, ion mobility and mass spectrometry techniques

Protonation isomers of gas-phase nicotineH⁺ are separated and assigned using a combination of FAIMS and UV photodissociation action spectroscopy.

Conformational properties of chiral tobacco alkaloids by DFT calculations and vibrational circular dichroism: (-)-S-anabasine

A thorough DFT and MM study of the conformational landscape, molecular and electronic structures of (-)-S-anabasine is reported aimed to reveal the mechanism controlling its conformational preference. Although the conformational flexibility and diversity of this system is quite extensive, only two structures are populated both in gas-phase and solution (CCl4 and DMSO). NBO-aided electronic structure analyses performed for the eight conformers representing minima in the potential energy surface of (-)-S-anabasine indicate that both steric and electrostatic factors are determinant in the conformational distribution of the sample in gas phase. Nonetheless, hyperconjugative effects are the key force tipping the balance in the conformational equilibrium between the two main rotamers. Increasing the polarity of the medium (using the IEF-PCM formalism) barely affect the conformational energy profile, although a slight increase in the theoretical population of those structures more affected by electrostatic interactions is predicted. The validity of the theoretical models and calculated conformers populations are endorsed by the accurate reproduction of the IR and VCD spectra (recorded in pure liquid and in CCl4 solution) of the sample (that have been firstly recorded and assigned in the present work) which are consistent with the occurrence of a 2:1 conformational ratio.

Vibrational circular dichroism and theoretical study of the conformational equilibrium in (-)-S-nicotine

We report an extensive study of the molecular and electronic structure of (-)-S-nicotine, to deduce the phenomenon that controls its conformational equilibrium and to solve its solution-state conformer population. Density functional theory, ab initio, and molecular mechanics calculations were used together with vibrational circular dichroism (VCD) and Fourier transform infrared spectroscopies. Calculations and experiments in solution show that the structure and the conformational energy profile of (-)-S-nicotine are not strongly dependent on the medium, thus suggesting that the conformational equilibrium is dominated by hyperconjugative interactions rather than repulsive electronic effects. The analysis of the first recorded VCD spectra of (-)-S-nicotine confirmed the presence of two main conformers at room temperature. Our results provide further evidence of the hypersensitivity of vibrational optical activity spectroscopies to the three-dimensional structure of chiral samples and prove their suitability for the elucidation of solution-state conformer distribution.