A computational and spectroscopic study of the gas-phase conformers of adrenaline

Accurate Computation of Aqueous pKas of Biologically Relevant Organic Acids: Overcoming the Challenges Posed by Multiple Conformers, Tautomeric Equilibria, and Disparate Functional Groups with the Fully Black-Box pK-Yay Method

The use of static electronic structure calculations to compute solution-phase pKas offers a great advantage in that a macroscopic bulk property could be computed via microscopic computations involving very few molecules. There are various sources of errors in the quantum chemical calculations though. Overcoming these errors to accurately compute pKas of a plethora of acids is an active area of research in physical chemistry pursued by both computational as well as experimental chemists. We recently developed the pK-Yay method in our attempt to accurately compute aqueous pKas of strong and weak acids. The method is fully black-box, computationally inexpensive, and is very easy for even a nonexpert to use. However, the method was thus far tested on very few molecules (only 16 in all). Herein, in order to assess the future applicability of pK-Yay, we study the effect of multiple conformers, the presence of tautomers under equilibrium, and the impact of a wide variety of functional groups (derivatives of acetic acid with substituents at various positions, dicarboxylic acids, aromatic carboxylic acids, amines and amides, phenols and thiols, and fluorine bearing organic acids). Starting with more than 1000 conformers and tautomers, this study establishes that overall errors of ∼ 1.0 pKa units are routinely obtained for a majority of the molecules. Larger errors are noted in cases where multiple charges, intramolecular hydrogen bonding, and several ionizable functional groups are simultaneously present. An important conclusion to emerge from this work is that, the computed pKas are insensitive (difference <0.5) to whether we consider multiple conformers/tautomers or only choose the most stable conformer/tautomer. Further, pK-Yay captures the stereoelectronic effects arising due to differing axial vs equatorial pattern, and is useful to predict the dominant acid-base equilibrium in a system featuring several equilibria. Overall, pK-Yay may be employed in several chemical applications featuring organic molecules and biomonomers.

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.

Hydration-controlled excited-state relaxation in protonated dopamine studied by cryogenic ion spectroscopy

Ultraviolet (UV) and infrared (IR) spectra of protonated dopamine (DAH⁺) and its hydrated clusters DAH⁺(H₂O)_1-3 are measured by cryogenic ion spectroscopy. DAH⁺ monomer and hydrated clusters with up to two water molecules show a broad UV spectrum, while it turns to a sharp, well-resolved one for DAH⁺-(H₂O)₃. Excited state calculations of DAH⁺(H₂O)₃ reproduce these spectral features. The conformer-selected IR spectrum of DAH⁺(H₂O)₃ is measured by IR dip spectroscopy, and its structure is assigned with the help of quantum chemical calculations. The excited state lifetime of DAH⁺ is much shorter than 20 ps, the cross correlation of the ps lasers, revealing a fast relaxation dynamics. The minimal energy path along the NH → π proton transfer coordinate exhibits a low energy barrier in the monomer, while this path is blocked by the high energy barrier in DAH⁺(H₂O)₃. It is concluded that the excited state proton transfer in DAH⁺ is inhibited by water-insertion.

Theoretical investigation on hydrogen bond interaction between adrenaline and hydrogen sulfide

In this study, we elucidated the formation of hydrogen bond between adrenaline (AD) and hydrogen sulfide utilizing computational studies. Six potential complexes were studied including geometrical parameters, energy, vibrational frequency, topological analysis, natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM), and NMR analysis. Moreover, these calculations were examined through DFT/ωB97XD/6-311G++(d,p) level. It was found that there are no indication on formation on hydrogen bonding between the two catecholic OHs where the one formed between the amino group and the hydroxyl oxygen atom of adrenaline monomer was broken in AS₁ to form two new interactions namely SH^...N and O₁H₁^...S, while it retained in other complexes. Furthermore, the bond became stronger due to cooperativity in AS₃ and AS₆, for the presence of withdrawing effect of the phenyl ring, the H-bonds formed with the side chain oxygen atom. The adrenaline and H₂S interaction was experimentally examined via FT-IR spectrometry and thin layer chromatography for confirmation of our theoretical study. Graphical abstract.

Gas-phase endstation of electron, ion and coincidence spectroscopies for diluted samples at the FinEstBeAMS beamline of the MAX IV 1.5 GeV storage ring

Since spring 2019 an experimental setup consisting of an electron spectrometer and an ion time-of-flight mass spectrometer for diluted samples has been available for users at the FinEstBeAMS beamline of the MAX IV Laboratory in Lund, Sweden. The setup enables users to study the interaction of atoms, molecules, (molecular) microclusters and nanoparticles with short-wavelength (vacuum ultraviolet and X-ray) synchrotron radiation and to follow the electron and nuclear dynamics induced by this interaction. Test measurements of N2 and thiophene (C4H4S) molecules have demonstrated that the setup can be used for many-particle coincidence spectroscopy. The measurements of the Ar 3p photoelectron spectra by linear horizontal and vertical polarization show that angle-resolved experiments can also be performed. The possibility to compare the electron spectroscopic results of diluted samples with solid targets in the case of Co2O3 and Fe2O3 at the Co and Fe L2,3-absorption edges in the same experimental session is also demonstrated. Because the photon energy range of the FinEstBeAMS beamline extends from 4.4 eV up to 1000 eV, electron, ion and coincidence spectroscopy studies can be executed in a very broad photon energy range.

Excited-state proton transfer in protonated adrenaline revealed by cryogenic UV photodissociation spectroscopy

We report a comprehensive study of the structures and deactivation processes of protonated adrenaline through cryogenic UV photodissociation spectroscopy. Single UV and double-resonance UV-UV hole burning spectroscopies have been performed and compared to coupled-cluster SCS-CC2 calculations performed on the ground and first electronic states. Three conformers were assigned, two lowest energy gauche conformers along with a higher energy conformer with an extended structure which is indeed the global minimum in solution. This demonstrates the kinetic trapping of this high energy gas phase conformer during the electrospray process. At the band origin of all conformers, the main fragmentation channel is the C_α-C_β bond cleavage, triggered by an excited state proton transfer to the catechol ring. Internal conversion leading to the water loss channel competes with the direct dissociation and tends to prevail with the increase of excess energy brought by the UV laser. Picosecond time-resolved pump-probe spectroscopy was performed to measure the excited state lifetimes of the three conformers of AdH⁺, which decay with the increase of excess energy in the ππ* state, from 2 ns at the band origin down to few hundreds of picoseconds 0.5 eV to the blue. Finally, about 0.8 eV above the band origin, the πσ* state is directly reached, leading to the opening of the H-loss channel.

Comparative antiradical activity and molecular Docking/Dynamics analysis of octopamine and norepinephrine: the role of OH groups

Octopamine is a neurotransmitter in invertebrates and a phenol analog of norepinephrine. The crystallographic and spectral (UV-visUV, and NMR) characteristics of octopamine were investigated experimentally and theoretically by applying appropriate level of theory, B3LYP-D3BJ/6-311++G(d,p), which reproduced well the experimental bond lengths and angles. The intramolecular interactions governing the stability of conformers were described by NBO and QTAIM analyses. The antiradical potencies of octopamine and norepinephrine towards DPPH and ABTS⁺ were examined with special emphasis on the preferred mechanism and effect of catechol moiety. Several techniques were used to distinguish Hydrogen Atom Transfer (HAT) and Proton Coupled Electron Transfer (PCET) mechanisms for reaction with DPPH. The calculated rate constants of the reactions with both radicals showed that Sequential Proton Loss Electron Transfer (SPLET) mechanism was dominant both thermodynamically and kinetically, with values of thermodynamic functions and rate constants clearly proving the importance of the second hydroxyl group in structure. The Molecular Docking and afterward Molecular Dynamics calculations of formed complexes between octopamine/norepinephrine with β1- and β2- adrenergic receptors examined in details the interactions that lead to the formation of stable complexes. The number of strong interactions of amino acids with norepinephrine was higher, but the absence of hydroxyl group in octopamine did not lead to a significant change in the type of interactions and stability. The formed complexes showed higher flexibility of amino acids, similar compactness of structure as proteins and increased interatomic distances of the backbone when compared to pure proteins.

On the applicability of functional-group symmetry-adapted perturbation theory and other partitioning models for chiral recognition - the case of popular drug molecules interacting with chiral phases

The applicability of symmetry-adapted perturbation theory (SAPT) and functional-group SAPT (F-SAPT) to study chiral recognition is investigated on an example of three popular chiral drug molecules: ibuprofen, norepinephrine, and baclofen, interacting with phenethylamine or proline - two molecules that are often used as chiral phases in chromatography. The comparison of the F-SAPT with the interacting quantum atoms (IQA) approach is also provided. Accurate estimation of energetic differences of the non-covalent intermolecular complexes composed of two chiral molecules is a necessary prerequisite for the possibility of a prediction of the chiral recognition. The SAPT method with interacting molecules described on the density functional theory level provides accurate total interaction energies, while the F-SAPT approach is the most useful in determining which functional groups are responsible for strengthening or weakening of the interaction between chiral molecules. The largest difference in the interaction energies has been calculated for the baclofen-phenethylamine and norepinephrine-phenethylamine pairs, while the smallest for the ibuprofen-proline and baclofen-proline ones. In most cases, the intermolecular interaction is found to be composed of a strong directional hydrogen bond, which was stabilized by two or more weaker non-covalent interactions between groups (in accordance with the phenomological three-point rule), but in several cases more subtle factors are responsible for larger stability of one diastereoisomer, like the stabilization of the conformation involving two noninteracting functional groups attached to a chiral atom through intramolecular attraction. Additionally, the simulated IR spectra were analyzed for all pairs of diastereoisomeric complexes and the red- and blue-shifts of characteristic bond vibrations were discussed in the context of inter-group interactions.

Theoretical insights into the hydrogen bonding interaction in the complexation of epinephrine with uracil

The present study is aimed at probing the hydrogen bonding interaction between epinephrine and uracil by means of density functional theory calculations concerning their complexation's geometries, interaction energies, and vibrational frequencies. Geometry optimization was carried out giving 19 stable geometries of epinephrine-uracil complex with interaction energies in a range of - 21.51 to - 62.37 kJ mol^-1 using the basis set superposition error (BSSE) correction. The analysis of structure and vibration shows that the hydrogen bonding elongates the length of corresponding bond O(N)-H and decreases the symmetric stretching vibrational frequency, which indicates red-shifted H-bonding interactions in all the geometries. Additionally, the analysis with theories of natural bond orbital (NBO), atoms in molecules (AIM), and the reduced density gradient (RDG) of hydrogen bonding properties and characteristics of the 19 geometries suggests that the hydrogen bonding in all the optimized structures of epinephrine-uracil complex is kind of a closed-shell interaction and mainly electrostatic dominant.

Theoretical Coupling and Stability of Boronic Acid Adducts with Catecholamines

Background:

Catecholamines combined with boric/boronic acids are attractive chemical agents in drug design because some of their adducts have shown interesting biological activity. Scant information exists about their stability.

Objective:

The aim of the present theoretical study was to explore the role of boron in molecules that combine catecholamines and boric/boronic acids, with a particular interest in examining stability.

Method:

The methodology was based on the US GAMESS program using DFT with the B3LYP exchange-correlation functional and the 6-31G (d,p) split-valence basis set.

Results:

According to the current findings, the boron-containing compounds (BCCs) exhibit weaker bonding to the hydroxyls on the ethylamine moiety than to those in the aromatic ring. The strongest binding site of a hydroxyl group was often found to be in meta-position (relative to ethylamine moiety) for boron-free compounds and in para-position for BCCs. Nonetheless, the methyl substituent in the amino group was able to induce changes in this pattern. We analyzed feasible boronsubstituted structures and assessed the relative strength of the respective C-B bonds, which allowed for the identification of the favorable points for reaction and stability.

Conclusion:

It is feasible to form adducts by bonding on the amine and catechol sides of catecholamines. The presence of boron stabilizes the adducts in para-position. Since some of these BCCs are promising therapeutic agents, understanding the mechanisms of reaction is relevant for drug design.

Rotational Spectrum and Conformational Analysis of N-Methyl-2-Aminoethanol: Insights into the Shape of Adrenergic Neurotransmitters

We describe an experimental and quantum chemical study for the accurate determination of the conformational space of small molecular systems governed by intramolecular non-covalent interactions. The model systems investigated belong to the biological relevant aminoalcohol's family, and include 2-amino-1-phenylethanol, 2-methylamino-1-phenylethanol, noradrenaline, adrenaline 2-aminoethanol, and N-methyl-2-aminoethanol. For the latter molecule, the rotational spectrum in the 6–18 and 59.6–74.4 GHz ranges was recorded in the isolated conditions of a free jet expansion. Based on the analysis of the rotational spectra, two different conformational species and 11 isotopologues were observed and their spectroscopic constants, including ¹⁴N-nuclear hyperfine coupling constants and methyl internal rotation barriers, were determined. From the experimental data a structural determination was performed, which was also used to benchmark accurate quantum chemical calculations on the whole conformational space. Atom in molecules and non-covalent interactions theories allowed the characterization of the position of the intramolecular non-covalent interactions and the energies involved, highlighting the subtle balance responsible of the stabilization of all the molecular systems.

A conformational study of protonated noradrenaline by UV–UV and IR dip double resonance laser spectroscopy combined with an electrospray and a cold ion trap method

In protonated noradrenaline, 3 folded and 2 extended conformers were identified under the ultra-cold condition.

IR spectrum of the protonated neurotransmitter 2-phenylethylamine: dispersion and anharmonicity of the NH3(+)-π interaction

The structure and dynamics of the highly flexible side chain of (protonated) phenylethylamino neurotransmitters are essential for their function. The geometric, vibrational, and energetic properties of the protonated neutrotransmitter 2-phenylethylamine (H(+)PEA) are characterized in the N-H stretch range by infrared photodissociation (IRPD) spectroscopy of cold ions using rare gas tagging (Rg = Ne and Ar) and anharmonic calculations at the B3LYP-D3/(aug-)cc-pVTZ level including dispersion corrections. A single folded gauche conformer (G) protonated at the basic amino group and stabilized by an intramolecular NH(+)-π interaction is observed. The dispersion-corrected density functional theory calculations reveal the important effects of dispersion on the cation-π interaction and the large vibrational anharmonicity of the NH3(+) group involved in the NH(+)-π hydrogen bond. They allow for assigning overtone and combination bands and explain anomalous intensities observed in previous IR multiple-photon dissociation spectra. Comparison with neutral PEA reveals the large effects of protonation on the geometric and electronic structure.

Gas phase ultraviolet and infrared spectroscopy on a partial peptide of β2-adrenoceptor SIVSF-NH2 by a laser desorption supersonic jet technique

Laser desorption supersonic jet laser spectroscopy has been applied to a penta-peptide, Ser-Ile-Val-Ser-Phe-NH₂ (SIVSF-NH₂), which is a partial sequence of a binding site in a β₂-adrenaline receptor protein.

Assessing the suitability of the multilevel strategy for the conformational analysis of small ligands

Predicting the conformational preferences of flexible compounds is a challenging problem in drug design, where the recognition between ligand and receptor is affected by the ability of the interacting partners to adopt a favorable conformation for the binding. To explore the conformational space of flexible ligands and to obtain the relative free energy of the conformation wells, we have recently reported a multilevel computational strategy that relies on the predominant-state approximation-where the conformational space is partitioned into distinct conformational wells-and combines a low-level method for sampling the conformational minima and high-level ab initio calculations for estimating their relative stability. In this study, we assess the performance of the multilevel strategy for predicting the conformational preferences of a series of structurally related phenylethylamines and streptomycin in aqueous solution. The charged nature of these compounds and the chemical complexity of streptomycin make them a challenging test for the multilevel approach. Furthermore, we explore the suitability of using a molecular mechanics approach as a source of approximate ensembles in the first stage of the multilevel strategy. The results support the reliability of the multilevel approach for obtaining an accurate conformational ensemble of small (bio)organic molecules in aqueous solution.

Competing intramolecular vs. intermolecular hydrogen bonds in solution

A hydrogen bond for a local-minimum-energy structure can be identified according to the definition of the International Union of Pure and Applied Chemistry (IUPAC recommendation 2011) or by finding a special bond critical point on the density map of the structure in the framework of the atoms-in-molecules theory. Nonetheless, a given structural conformation may be simply favored by electrostatic interactions. The present review surveys the in-solution competition of the conformations with intramolecular vs. intermolecular hydrogen bonds for different types of small organic molecules. In their most stable gas-phase structure, an intramolecular hydrogen bond is possible. In a protic solution, the intramolecular hydrogen bond may disrupt in favor of two solute-solvent intermolecular hydrogen bonds. The balance of the increased internal energy and the stabilizing effect of the solute-solvent interactions regulates the new conformer composition in the liquid phase. The review additionally considers the solvent effects on the stability of simple dimeric systems as revealed from molecular dynamics simulations or on the basis of the calculated potential of mean force curves. Finally, studies of the solvent effects on the type of the intermolecular hydrogen bond (neutral or ionic) in acid-base complexes have been surveyed.

Seven Conformers of Neutral Dopamine Revealed in the Gas Phase

The rotational spectrum of neutral dopamine has been investigated for the first time using a combination of Fourier transform microwave spectroscopy with laser ablation. The parameters extracted from the analysis of the spectrum unequivocally identify the existence of seven conformers of dopamine. (14)N nuclear quadrupole coupling interactions have been used to determine the orientation of the amino group probing the existence of stabilizing N-H···π interactions for all observed conformers.

Gas-phase spectroscopy of synephrine by laser desorption supersonic jet technique

In our previous work, we found that synephrine has six conformers in the gas phase, while adrenaline, which is a catecholamine and has the same side chain as synephrine, has been reported to have only two conformers. To determine the conformational geometries of synephrine, we measured resonance enhanced multiphoton ionization, ultraviolet-ultraviolet hole burning, and infrared dip spectra by utilizing the laser desorption supersonic jet technique. By comparing the observed infrared spectra with theoretical ones, we assigned geometries except for the orientations of the phenolic OH group. Comparison between the determined structures of synephrine and those of 2-methylaminno-1-phenylethanol, which has the same side chain as synephrine but no phenol OH group, leads to the conclusion that the phenolic OH group in synephrine does not affect the conformational flexibility of the side chain. In the case of adrenaline, which is expected to have 12 conformers if there are no interactions between the catecholic OH groups and the side chain, some interactions possibly exist between them because only two conformations are observed. By estimation of the dipole-dipole interaction energy between partial dipole moments of the catecholic OH groups and the side chain, it was concluded that the dipole-dipole interaction stabilizes specific conformers which are actually observed.

Thermodynamic functions of molecular conformations of (2-fluoro-2-phenyl-1-ethyl)ammonium ion and (2-hydroxy-2-phenyl-1-ethyl)ammonium ion as models for protonated noradrenaline and adrenaline: first-principles computational study of conformations and thermodynamic functions for the noradrenaline and adrenaline models

This paper reports the structural and thermodynamic consequences of substitution of the OH group by the isoelectronic F-atom in the case of the adrenaline family of molecules. The conformational landscapes were explored for the two enantiomeric forms of N-protonated-beta-fluoro-beta-phenyl-ethylamine, also called (2-fluoro-2-phenyl-1-ethyl)-ammonium ion (Model 1) and that of N-protonated-beta-hydroxy-beta-phenyl-ethylamine, also referred to as (2-hydroxy-2-phenyl-1-ethyl)-ammonium (Model 2) models of noradrenaline and adrenaline molecules. These full conformational studies were carried out by first principles of quantum mechanical computations at the B3LYP/6-31G(d,p) and G3MP2B3 levels of theory, using the Gaussian03 program. Also, frequency calculations of the stable structures were performed at the B3LYP/6-31G(d,p), and G3MP2B3 levels of theory. The thermodynamic functions (U, H, S, and G) of the various stable conformations of the title compounds were calculated at these levels of theory for the R and S stereoisomers. Relative values of the thermodynamic functions have been calculated with respect of the chosen reference conformers in which all relevant dihedral angles assumed anti orientation for the Model 1 and Model 2. Through the combination of both point and axis chirality, the enantiomeric and diastereomeric relationships of the six structures for each molecule investigated were established. Intramolecular hydrogen bonding interactions have been studied by the atoms in molecules (AIM) analysis of the electron density. The aromaticity of phenyl group has been determined by a selective hydrogenation protocol. The pattern of the extent of aromacity, due intramolecular interactions, varies very little between the two models studied.

The effect of solvation on biomolecular conformation: 2-amino-1-phenylethanol

Small molecule neurotransmitters form one the most important classes of pharmaceutical molecules. While the behavior of these molecules in their neutral forms in the gas phase is well understood, their behavior in more biologically relevant scenarios (protonated and in aqueous solution) has received comparatively little attention. Here we address this problem by using molecular mechanics simulations to build up a detailed picture of the conformational behavior of 2-amino-1-phenylethanol, a noradrenaline analogue, in aqueous solution in both its neutral and protonated forms. For the sake of comparison, equivalent simulations are also performed on the gas-phase molecules and gas-phase hydrated clusters. These calculations reveal the important role that water has to play in determining the conformational preferences and dynamic behavior of the molecules. Water molecules are found to bridge between the various functional groups within the molecule, significantly affecting their relative stabilities in comparison to the gas-phase values. The reorganization of these solvation structures also provides a mechanism for conformational interconversion. The role of the solvent in mediating interactions between the various functional groups within the molecule suggests that in noradrenaline the catechol groups will be able to interact, albeit indirectly, with the other functional groups, thereby influencing the behavior of the molecule.