Hydration of neurotransmitters: a computational and spectroscopic study of a noradrenaline analogue, 2-amino-1-phenyl-ethanol

Probing the selectivity of Li+ and Na+ cations on noradrenaline at the molecular level

Although several mechanisms concerning the biological function of lithium salts, drugs having tranquilizing abilities, have been proposed so far, the key mechanism for its selectivity and subsequent interaction with neurotransmitters has not been established yet. We report ultraviolet (UV) and infrared (IR) spectra under ultra-cold conditions of Li+ and Na+ complexes of noradrenaline (NAd, norepinephrine), a neurotransmitter responsible for the body's response to stress or danger, in an effort to provide a molecular level understanding of the conformational changes of NAd due to its interactions with these two cations. A detailed analysis of the IR spectra, aided by quantum chemical calculations, reveals that the Li+-noradrenaline (NAd-Li+) complex forms only an extended structure, while the NAd-Na+ and protonated (NAd-H+) complexes form both folded and extended structures. This conformational preference of the NAd-Li+ complex is further explained by considering specific conformational distributions in solution. Our results clearly discern the unique structural motifs that NAd adopts when interacting with Li+ compared with other abundant cations in the human body (Na+) and can form the basis of a molecular level understanding of the selectivity of Li+ in biological systems.

High energy conformers of M(+)(APE)(H2O)(0-1)Ar(0-1) clusters revealed by combined IR-PD and DFT-MD anharmonic vibrational spectroscopy

IR-PD vibrational spectroscopy and DFT-based molecular dynamics simulations are combined in order to unravel the structures of M(+)(APE)(H2O)0-1 ionic clusters (M = Na, K), where APE (2-amino-1-phenyl ethanol) is commonly used as an analogue for the noradrenaline neurotransmitter. The strength of the synergy between experiments and simulations presented here is that DFT-MD provides anharmonic vibrational spectra that unambiguously help assign the ionic clusters structures. Depending on the interacting cation, we have found that the lowest energy conformers of K(+)(APE)(H2O)0-1 clusters are formed, while the lowest energy conformers of Na(+)(APE)(H2O)0-1 clusters can only be observed through water loss channel (i.e. without argon tagged to the clusters). Trapping of higher energy conformers is observed when the argon loss channel is recorded in the experiment. This has been rationalized by transition state energies. The dynamical anharmonic vibrational spectra unambiguously provide the prominent OH stretch due to the OH···NH2 H-bond, within 10 cm(-1) of the experiment, hence reproducing the 240-300 cm(-1) red-shift (depending on the interacting cation) from bare neutral APE. When this H-bond is not present, the dynamical anharmonic spectra provide the water O-H stretches as well as the rotational motion of the water molecule at finite temperature, as observed in the experiment.

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.

Structure, electronic circular dichroism and Raman optical activity in the gas phase and in solution: a computational and experimental investigation

A computational (ab initio and molecular dynamics) and experimental exploration of the relative importance of molecular conformation and explicit solvent effects on the electronic circular dichroism (ECD) of chiral molecules, is presented. The exploration includes an assessment of the validity of angular correlation (sector) rules linking ECD to molecular conformation. It is based upon studies of 1-(R) phenylethanol (including its Raman optical activity spectrum), the corresponding 'benchmark' base, 1-(R)-phenylethylamine and its protonated cation; their hydrated clusters in the gas phase; and their non-polar and aqueous solutions. Emphasis is placed on the influence of specific, hydrogen bonded interactions with the aqueous solvent. The theoretical validity of the (otherwise empirical) sector rule in the neutral molecules and in their specifically hydrated clusters has been established--but with a reversal of the 'historical' sign convention. Protonation of the amine leads to a breakdown of the conventional sector rule but the change in its ECD intensity can still be related to the side chain dihedral angular dependence of its rotatory strength, computed ab initio for its explicitly hydrated clusters. Comparisons between ECD spectra measured in aqueous and in hydrocarbon solutions and the results of molecular dynamics calculations for aqueous solutions at 300 K, identify solvent induced structural change as the principal determinant of their relative ECD spectral intensities. Further links connecting the structures and conformations of chiral molecules and their explicitly solvated clusters in the gas phase, to their structures and conformational populations in solution can be expected through measurement, ab initio computation and analysis of their vibrational, ROA spectra.

Chirality recognition between neutral molecules in the gas phase

Noncovalent interactions are particularly intriguing when they involve chiral molecules, because the interactions change in a subtle way upon replacing one of the partners by its mirror image. The resulting phenomena involving chirality recognition are relevant in the biosphere, in organic synthesis, and in polymer design. They may be classified according to the permanent or transient chirality of the interacting partners, leading to chirality discrimination, chirality induction, and chirality synchronization processes. For small molecules, high-level quantum chemical calculations for such processes are feasible. To provide reliable connections between theory and experiment, such phenomena are best studied in vacuum isolation at low temperature, using rotational, vibrational, electronic, and photoionization spectroscopy. We review these techniques and the results which have become available in recent years, with special emphasis on dimers of permanently chiral molecules and on the influence of conformational flexibility. Analogies between the microscopic mechanisms and macroscopic phenomena and between intra- and intermolecular cases are drawn.

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.

A solvent induced mechanism for conformational change

Molecular dynamics simulations have been used to investigate the dynamic behaviour of two small molecule neurotransmitter analogues in aqueous solution, leading to the elucidation of a mechanism for conformational change which is driven by the presence of the solvent molecules.

Quantum Simulation of a Hydrated Noradrenaline Analog with the Torsional Path Integral Method

An extended version of the torsional path integral Monte Carlo (TPIMC) method is presented and shown to be useful for studying the conformation of flexible molecules in solvated clusters. The new technique is applied to the hydrated clusters of the 2-amino-1-phenyl-ethanol (APE) molecule. APE + nH2O clusters with n = 0-4 are studied at 100 and 300 K using both classical and quantum simulations. Only at the lower temperature is the hydration number n found to impact the conformational distribution of the APE molecule. This is shown to be a result of the temperature-dependent balance between the internal energy and entropy contributions to the relative conformer free energies. Furthermore, at 100 K, large quantum effects are observed in the calculated conformer populations. A particularly large quantum shift of 30% of the total population is calculated for the APE + 2H2O cluster, which is explained in terms of the relative zero point energy of the lowest-energy hydrated structures for this cluster. Finally, qualitative agreement is found between the reported calculations and recent spectroscopy experiments on the hydrated clusters of APE, including an entropically driven preference for the formation of AG-type hydrated structures and the formation of a water "droplet" in the APE + 4H2O cluster.

First-principles quantum chemistry in the life sciences

The area of computational quantum chemistry, which applies the principles of quantum mechanics to molecular and condensed systems, has developed drastically over the last decades, due to both increased computer power and the efficient implementation of quantum chemical methods in readily available computer programs. Because of this, accurate computational techniques can now be applied to much larger systems than before, bringing the area of biochemistry within the scope of electronic-structure quantum chemical methods. The rapid pace of progress of quantum chemistry makes it a very exciting research field; calculations that are too computationally expensive today may be feasible in a few months' time! This article reviews the current application of 'first-principles' quantum chemistry in biochemical and life sciences research, and discusses its future potential. The current capability of first-principles quantum chemistry is illustrated in a brief examination of computational studies on neurotransmitters, helical peptides, and DNA complexes.