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.

Descriptors of water aggregation.

We rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (rOH) and HOH bond angle (θHOH), (ii) two-body: OO distance (rOO), OHO angle (θOHO), and HOOX dihedral angle (ϕHOOX), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θOOO), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (qm, m = 2, 3, 4, 5) and radius of gyration (Rg). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H2O)n, n = 3-25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that rOH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.

Switching of Protonation Sites in Hydrated Nicotine via a Grotthuss Mechanism.

The switching of the protonation sites in hydrated nicotine, probed by experimental infrared (IR) spectroscopy and theoretical ab initio calculations, is facilitated via a Grotthuss instead of a bimolecular proton transfer (vehicle) mechanism at the experimental temperature (T = 130 K) as unambiguously confirmed by experiments with deuterated water. In contrast, the bimolecular vehicle mechanism is preferred at higher temperatures (T = 300 K) as determined by theory. The Grotthuss mechanism for the concerted proton transfer results in the production of nicotine's bioactive and addictive pyrrolidine-protonated (Pyrro-H+) protomer with just 5 water molecules. Theoretical analysis suggests that the concerted proton transfer occurs via hydrogen-bonded bridges consisting of a 3 water molecule "core" that connects the pyridine protonated (Pyri-H+) with the pyrrolidine-protonated (Pyrro-H+) protomers. Additional water molecules attached as acceptors to the hydrogen-bonded "core" bridge result in lowering the reaction barrier of the concerted proton transfer down to less than 6 kcal/mol, which is consistent with the experimental observations.

Vibrational levels of a generalized Morse potential.

A Generalized Morse Potential (GMP) is an extension of the Morse Potential (MP) with an additional exponential term and an additional parameter that compensate for MP's erroneous behavior in the long range part of the interaction potential. Because of the additional term and parameter, the vibrational levels of the GMP cannot be solved analytically, unlike the case for the MP. We present several numerical approaches for solving the vibrational problem of the GMP based on Galerkin methods, namely, the Laguerre Polynomial Method (LPM), the Symmetrized LPM, and the Polynomial Expansion Method (PEM), and apply them to the vibrational levels of the homonuclear diatomic molecules B2, O2, and F2, for which high level theoretical near full configuration interaction (CI) electronic ground state potential energy surfaces and experimentally measured vibrational levels have been reported. Overall, the LPM produces vibrational states for the GMP that are converged to within spectroscopic accuracy of 0.01 cm-1 in between 1 and 2 orders of magnitude faster and with much fewer basis functions/grid points than the Colbert-Miller Discrete Variable Representation (CN-DVR) method for the three homonuclear diatomic molecules examined in this study. A Python library that fits and solves the GMP and similar potentials can be downloaded from https://gitlab.com/gds001uw/generalized-morse-solver.

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+···(H2O) 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).