Acetylcholine is the first discovered neurotransmitter that has received a great attention regarding its capability of binding to several cellular targets. The chemical composition of acetylcholine, including a positively charged trimethylammonium and a carbonyl group, as well as its conformational flexibility was pointed out as the key factors in the stabilization of its interactions. Here, the possibilities offered by a Raman scattering-based multiconformatioal analysis to access the most stable conformers of acetylcholine, is discussed. To control the validity of this protocol, acetylcholine and one of its closely structured analogues, acetylthiocholine, were simultaneously analyzed. Solution Raman spectra revealed distinct and well resolved strong markers for each molecule. Density functional theory calculations were consistent with the fact that the energy order of the low energy conformers is considerably affected by the acyloxy oxygen→sulfur atom substitution. Raman spectra were calculated on the basis of the thermal average of the spectra arising from the low energy conformers. It has been evidenced that the carbonyl and trimethylammonium groups are the most favorable hydration sites in aqueous environment. Taking into account the large gap between the carbonyl bond-stretch and aliphatic bending bands, Raman spectra also allowed separation of the HOH bending vibrations arising from the bound and bulk water molecules.