Electrochemical Shell-Isolated Nanoparticle-Enhanced Raman Spectroscopy of Imidazole Ring Functionalized Monolayer on Smooth Gold Electrode

The imidazole ring (Im) of histidine side chains plays a unique role in the function of proteins through covalent bonding with metal ions and hydrogen bonding interactions with adjusted biomolecules and water. At biological interfaces, these interactions are modified because of the presence of an electric field. Self-assembled monolayers (SAMs) with the functional Im group mimic the histidine side chain at electrified interfaces. In this study, we applied in-situ shell-isolated nanoparticle-enhanced Raman spectroscopy (SHINERS) to probe the structure and hydrogen bonding of Im-functionalized SAM on smooth Au at the electrochemical interface. The self-assembly of molecules on the Au induced the proton shift from N1 atom (Tautomer-I), which is the dominant form of Im in the bulk sample, to N3 atom (Tautomer-II). The impact of electrode potential on the hydrogen bonding interaction strength of the Im ring was identified by SHINERS. Temperature-Raman measurements and density functional theory (DFT) analysis revealed the spectral marker for Im ring packing (mode near 1496–1480 cm−1) that allowed us to associate the confined and strongly hydrogen bonded interfacial Im groups with electrode polarization at −0.8 V. Reflection adsorption IR (RAIR) spectra of SAMs with and without Im revealed that the bulky ring prevented the formation of a strongly hydrogen bonded amide group network.


Introduction
Histidine plays a crucial role in the architecture and activity of enzymes, as it is often found to ligate transition metal ions at the active sites of metalloproteins [1,2]. Such a versatility comes from its side chain imidazole (Im) ring that can take part in many interactions, for example (i) cation-π Im , (ii) π-π stacking, (iii) hydrogen Im -π, (iv) coordination with metal cations through nitrogen lone electron pair, and (v) hydrogen bonding (H-bonding) [2]. This versatility makes Im particularly suitable for surface applications through the SAM chemistry [3][4][5][6][7][8]. The Im-derivatives functionalized SAMs have been used in biosensing applications, modulations of enzyme activity, purification of biologically relevant molecules, corrosion prevention, and others [5][6][7][8].
The pK a for imidazole is around 5.9. Thus, two forms of histidine side-chain functional group, Tautomer-I (T-I; N1-H, N3) and Tautomer-II (T-II; N1, N3-H), are possible under the physiological conditions. In neutral pH and room temperature, T-I is the more energetically favorable [9]. However, owing to the intermolecular Im interactions and the surface charge effect, the tautomeric equilibrium of surface-adsorbed Im derivatives might differ from that of the solution phase. Indeed, a clear transition from T-I to T-II was observed

Synthesis of Silicon Dioxide Covered Spherical Silver Nanoparticles (Ag@SiO2)
The reagents and solvents were used without further purifications: silver nitrate (AgNO3, 99%), trisodium citrate dihydrate (99%), (3-aminopropyl)triethoxysilane (99%), and sodium silicate solution (NaOH 10%, SiO2 27%) (Merck). All solutions were prepared with ultra-pure water (resistivity of 18.2 MΩ·cm) from Direct-Q 3UV (Merck, Darmstadt, Germany). Nanoparticles (NPs) were synthesized under microwave irradiation according to the previously published method [44]. Briefly, the bare silver-core NPs were prepared by the AgNO3 reduction with sodium citrate and capped with silica shell. Then, core-shell NPs were purified by a repeated centrifugation, supernatant removal, and resuspension. Figure 2 shows the Ag@SiO2 NPs images obtained by high-resolution transmission electron microscopy (HR-TEM). HR-TEM images were collected by an FEI Tecnai G2 F20 X-TWIN TEM (FEI, Netherlands) microscope with an accelerating voltage of 200 kV. The microscope was equipped with an EDX (EDAX) spectrometer and Gatan Orius CCD camera. The measurements were carried out in a bright-field regime.

Preparation and Characterization of SAM
A 120-nm Au film was deposited on clean glass slides by using a Quorum 150T magnetron-sputtering machine and a 99.99% Au target. After that, the slides were incubated in 10 −3 M IMHA ethanol solution for twenty-four hours unless indicated otherwise; then

Synthesis of Silicon Dioxide Covered Spherical Silver Nanoparticles (Ag@SiO 2 )
The reagents and solvents were used without further purifications: silver nitrate (AgNO 3 , 99%), trisodium citrate dihydrate (99%), (3-aminopropyl)triethoxysilane (99%), and sodium silicate solution (NaOH 10%, SiO 2 27%) (Merck). All solutions were prepared with ultra-pure water (resistivity of 18.2 MΩ·cm) from Direct-Q 3UV (Merck, Darmstadt, Germany). Nanoparticles (NPs) were synthesized under microwave irradiation according to the previously published method [44]. Briefly, the bare silver-core NPs were prepared by the AgNO 3 reduction with sodium citrate and capped with silica shell. Then, core-shell NPs were purified by a repeated centrifugation, supernatant removal, and resuspension. Figure 2 shows the Ag@SiO 2 NPs images obtained by high-resolution transmission electron microscopy (HR-TEM). HR-TEM images were collected by an FEI Tecnai G2 F20 X-TWIN TEM (FEI, Netherlands) microscope with an accelerating voltage of 200 kV. The microscope was equipped with an EDX (EDAX) spectrometer and Gatan Orius CCD camera. The measurements were carried out in a bright-field regime.

Synthesis of Silicon Dioxide Covered Spherical Silver Nanoparticles (Ag@SiO2)
The reagents and solvents were used without further purifications: silver nitrate (AgNO3, 99%), trisodium citrate dihydrate (99%), (3-aminopropyl)triethoxysilane (99%), and sodium silicate solution (NaOH 10%, SiO2 27%) (Merck). All solutions were prepared with ultra-pure water (resistivity of 18.2 MΩ·cm) from Direct-Q 3UV (Merck, Darmstadt, Germany). Nanoparticles (NPs) were synthesized under microwave irradiation according to the previously published method [44]. Briefly, the bare silver-core NPs were prepared by the AgNO3 reduction with sodium citrate and capped with silica shell. Then, core-shell NPs were purified by a repeated centrifugation, supernatant removal, and resuspension. Figure 2 shows the Ag@SiO2 NPs images obtained by high-resolution transmission electron microscopy (HR-TEM). HR-TEM images were collected by an FEI Tecnai G2 F20 X-TWIN TEM (FEI, Netherlands) microscope with an accelerating voltage of 200 kV. The microscope was equipped with an EDX (EDAX) spectrometer and Gatan Orius CCD camera. The measurements were carried out in a bright-field regime.

Preparation and Characterization of SAM
A 120-nm Au film was deposited on clean glass slides by using a Quorum 150T magnetron-sputtering machine and a 99.99% Au target. After that, the slides were incubated in 10 −3 M IMHA ethanol solution for twenty-four hours unless indicated otherwise; then

Preparation and Characterization of SAM
A 120-nm Au film was deposited on clean glass slides by using a Quorum 150T magnetron-sputtering machine and a 99.99% Au target. After that, the slides were incubated in 10 −3 M IMHA ethanol solution for twenty-four hours unless indicated otherwise; then For SHINERS experiments, the slides were mounted into an electrochemical cell and filled with Milli-Q water. After that, 3 µL of Ag@SiO 2 were carefully injected directly on the surface of Au and left to rest for 10 min for nanoparticles to adsorb. Then, the cell was thorough rinsed with 0.01 M phosphate buffer (pH 7) containing 0.1 M Na 2 SO 4 .
Reflection absorption IR spectroscopy (RAIRS) spectra were obtained by using a Vertex 80v FTIR spectrometer (Bruker, Germany) equipped with the LN-MCT narrow band detector and the horizontal reflection accessory. The spectral resolution was set at 4 cm −1 . Spectra were acquired by 400 scans at a grazing angle of 80 • by using p-polarized light. The sample chamber and the spectrometer were evacuated during the measurements to approximately 2 hPa pressure. The spectrum of deuterium substituted octadecanethiol (ODT-d38) SAM adsorbed on Au was used as a reference. FTIR transmission spectrum was recorded from KBr pellet-dispersed IMHA using an Alpha spectrometer (Bruker, Germany) equipped with an RT-DTGS detector. The resolution was set to 4 cm −1 ; 50 interferogram scans were co-added.
SHINERS spectra were recorded using high throughput (instrument NA = 0.22) Tornado HyperFlux spectrometer (Tornado Spectral Systems, Mississauga, ON, Canada) equipped with fiber-optic cable for excitation and collection of the Raman spectra. The 785 nm beam of the diode laser was used as the excitation source. The laser power at the sample was set to 38 mW, and the beam was focused on a 100 µm diameter spot on the sample. Spectra were recorded for 300 s by co-adding thirty 10-s scans. For temperaturecontrolled Raman measurements, the LinKam temperature control system PE95/T95 with the accuracy of 0.05 • C was used. IMHA powder was measured with 80 mW laser power and 100 s integration time, at 22, 130, and 135 • C. Spectroelectrochemical measurements were carried out in a three-electrode cell, where slides with Au film were used as a working electrode, platinum wire as a counter electrode, and KCl saturated Ag/AgCl as a reference electrode.
Raman wavenumbers were calibrated according to the polystyrene spectrum. Band frequencies were obtained by fitting experimental contours with mixed Gaussian-Lorentzian form components by using the GRAMS/AI 8.0 (Thermo Scientific, Waltham, MA, USA) software.
Theoretical modeling study of IMHA molecule was performed using Gaussian package version G09 D.01 [45]. Geometry optimization and frequency calculation were completed with the DFT method using the hybrid B3LYP functional and 6-311++G(2d,p) basis set. The polarizable continuum model (IEFPCM) was used to emulate the water environment. The frequency calculation ended with no imaginary wavenumbers, indicating the geometry at minimum energy. Frequency and intensity scaling were applied according to the previously described procedure [46].

Assignments of Raman Bands
The studied IMHA molecule comprises four functional units: (i) the surface-active thiol group (SH), (ii) the hydrocarbon chain (−(CH 2 ) 5 −), (iii) the amide group (−CO−NH−), and (iv) the terminal Im ring (Figure 3). To facilitate Raman band assignments, we have synthesized a compound similar to IMHA thiol without the terminal Im ring (fragment compound) (Figure 3). Detailed assignments of the Raman bands are provided in the Supplementary Data file ( Figure S2 and related discussion) and Table 1.  To evaluate Raman modes related to intra-and intermolecular interactions and those sensitive to the ordering of alkyl chain, we performed a temperature-dependent Raman study of powder IMHA compound ( Figure 4). Changes in temperature provoke structural or phase transitions and a variety of conformational alterations in a molecular system. In general, an increase in temperature induces disordering in the molecular arrangement. Such disordering may be reflected in changes of spectral mode bandwidths (related to vibrational energy distribution), peak position (energy of particular vibration), and relative intensities (distribution of two or more different forms of a molecule). Therefore, temperature-Raman provides additional means for scrutinizing a given compound. The medium intensity modes near 1566, 1321, and 983 cm −1 in the 22 °C spectrum are related to vibrations of the Tautomer-I (N1-H, N3) form of Im ring and are assigned to ν(C4=C5) stretching, ν(Im) breathing + δ(C5H), and β(CH) deformation, respectively [19]. At the elevated temperatures (130-135 °C), the shoulders appear at 1582 and 1345 cm −1 specific to the Tautomer-II (N1, N3-H) form of Im. The relative percentage of T-II determined by the integral intensity ratio A1582/(A1582 + A1564) is 16% for the compound at 130 °C and 32% for the one at 135 °C. It should be noted that the frequency of an intense band at 1321 cm −1 downshifts to 1305 cm −1 at 135 °C. In addition, the band notably broadens. Such spectral changes reflect a decrease in H-bonding interaction strength primarily at N3 site, because of a considerable contribution from ν(N3−C2) vibration to this mode for Tautomer-I [19].
Abbreviations: n.a., not applicable; sh, shoulder; G, gauche; T, trans; r, rocking; w, wagging; t, twisting; δ, deformation; β, in-plane deformation; γ, out-of-plane deformation; ν, stretching; Im, imidazole; T-I, Tautomer-I;  T-II, Tautomer-II. To evaluate Raman modes related to intra-and intermolecular interactions and those sensitive to the ordering of alkyl chain, we performed a temperature-dependent Raman study of powder IMHA compound ( Figure 4). Changes in temperature provoke structural or phase transitions and a variety of conformational alterations in a molecular system. In general, an increase in temperature induces disordering in the molecular arrangement. Such disordering may be reflected in changes of spectral mode bandwidths (related to vibrational energy distribution), peak position (energy of particular vibration), and relative intensities (distribution of two or more different forms of a molecule). Therefore, temperature-Raman provides additional means for scrutinizing a given compound. The medium intensity modes near 1566, 1321, and 983 cm −1 in the 22 • C spectrum are related to vibrations of the Tautomer-I (N1-H, N3) form of Im ring and are assigned to ν(C4=C5) stretching, ν(Im) breathing + δ(C5H), and β(CH) deformation, respectively [19]. At the elevated temperatures (130-135 • C), the shoulders appear at 1582 and 1345 cm −1 specific to the Tautomer-II (N1, N3-H) form of Im. The relative percentage of T-II determined by the integral intensity ratio A 1582 /(A 1582 + A 1564 ) is 16% for the compound at 130 • C and 32% for the one at 135 • C. It should be noted that the frequency of an intense band at 1321 cm −1 downshifts to 1305 cm −1 at 135 • C. In addition, the band notably broadens. Such spectral changes reflect a decrease in H-bonding interaction strength primarily at N3 site, because of a considerable contribution from ν(N3−C2) vibration to this mode for Tautomer-I [19].  The 1497 cm −1 band assigned to ν(N1-C2) + β(C2H) motion is not sensitive to Im tautomerism [19]. The band, however, shifts to higher wavenumbers with the coordination of metal cation at the nitrogen lone electron pair. Moreover, electrode polarization also induces shifting for imidazole-copper pair (9 ± 2 cm −1 V −1 ) [18]. We found a 17 cm −1 downshift of the spectral mode to 1480 cm −1 and a clear increase in full width at half-maximum (FWHM) with the temperature raised to 135 °C. These spectral changes are associated with the increased motional freedom of the Im ring. Therefore, the spectral changes near 1496-1480 cm −1 may serve as a useful spectral marker for Im ring packing.  The 1497 cm −1 band assigned to ν(N1-C2) + β(C2H) motion is not sensitive to Im tautomerism [19]. The band, however, shifts to higher wavenumbers with the coordination of metal cation at the nitrogen lone electron pair. Moreover, electrode polarization also induces shifting for imidazole-copper pair (9 ± 2 cm −1 V −1 ) [18]. We found a 17 cm −1 downshift of the spectral mode to 1480 cm −1 and a clear increase in full width at halfmaximum (FWHM) with the temperature raised to 135 • C. These spectral changes are associated with the increased motional freedom of the Im ring. Therefore, the spectral changes near 1496-1480 cm −1 may serve as a useful spectral marker for Im ring packing.
Nested within the alkane chain, the amide group acts as an intermolecular stabilizing agent, which forms the extended H-bond network between neighboring molecular chains. It has been already shown that such an interaction radically increases the desorption temperature and chemical stability of the monolayer [50]. From the temperature-Raman spectra of the powder compound, only the medium-low intensity Amide-I (Am-I) band at 1636 cm −1 can be recognized. This spectral mode is related to C=O stretching (83%) coupled with out-of-phase ν(C−N) and δ(C-C−N) and can serve as a diagnostic tool in identifying the secondary structure of peptides [47]. The 14 cm −1 frequency upshift during the solid-to-liquid transition with the temperature elevation, clearly shows the weakening of the H-bonding at the amide's C=O group. Figure 5 shows RAIRS spectra of the IMHA monolayer adsorbed on a smooth gold surface. The spectral features at 1645, 1556, and 1264 cm −1 are assigned to amide bands, Am-I, Am-II, and Am-III, respectively, whereas the ones near 1460 and 1380 cm −1 to methylene scissoring and wagging deformations, δ(CH 2 ). The RAIRS-surface selection rule allows interrogating the orientation of molecular groups because the intensity of particular spectral mode directly depends on the projection of the mode's transition dipole moment (TDM) on the surface normal [51]. While Am-I and Am-II have TDMs oriented perpendicular to each other in the amide bond plane, they also have perpendicular and parallel TDM orientations with respect to carbohydrate chain [51]. Thus, for the neatly packed IMHA monolayer, the Am-II is expected to dominate the spectrum. Indeed, at progressively longer incubation time, the Am-II intensity at 1556 cm −1 increases and the intensity of Am-I at 1645 cm −1 decreases. The integral intensity ratio Am-II/Am-I provides a qualitative measure of the molecular reorientation ( Figure 5B). The experimental data of Am-II/Am-I were fitted with a sigmoidal function:

RAIRS Analysis of the Monolayer Formation
where t m is the transition inflection point found at 162 min. At 10 s incubation time, the Am-II/Am-I ratio was 1.4, which is slightly above 0.9 calculated for IMHA molecules chaotically dispersed in KBr pellet ( Figure S3). The ratio doubled in 10 min and after 24 h it reached 11.2, so most molecules had a planar amide group oriented perpendicularly to the surface. Such an orientation is strengthened by the H-bonding interaction involving the amide group (C=O···H) and dipole-dipole interaction between the C=O groups [52]. The link between Amide-I and -II wavenumbers and the H-bonding strength at C=O and N-H groups has been clearly established [53][54][55]. We find marginal Am-I frequency decrease with the development of SAM (by 3 cm −1 during the 120 min incubation) revealing minor increase of the H-bonding strength at C=O. The frequency of the Am-II band remained constant.
In the RAIRS spectrum of 24 h-incubated fragment SAM, the Amide-II spectral band is blue-shifted by 12 cm −1 compared to the IMHA SAM (after 24 h incubation), showing the stronger H-bonding at the N-H group for fragment molecule. Interestingly, surfaceadsorption of molecules from solution phase induced spectral shift of Am-II mode by −20 for IMHA and 13 cm −1 for the fragment molecule (Table 2). Such frequency shifting is related to different molecular packing efficiency in bulk and SAM. The Im ring in IMHA molecule introduces sterical hindrances for the neighboring amides to engage in the formation of an optimal H-bonding network.
the amide group (C=O···H) and dipole-dipole interaction between the C=O groups [52] The link between Amide-I and -II wavenumbers and the H-bonding strength at C=O and N-H groups has been clearly established [53][54][55]. We find marginal Am-I frequency de crease with the development of SAM (by 3 cm −1 during the 120 min incubation) revealing minor increase of the H-bonding strength at C=O. The frequency of the Am-II band re mained constant. In the RAIRS spectrum of 24 h-incubated fragment SAM, the Amide-II spectral band is blue-shifted by 12 cm −1 compared to the IMHA SAM (after 24 h incubation), showing the stronger H-bonding at the N-H group for fragment molecule. Interestingly, surface adsorption of molecules from solution phase induced spectral shift of Am-II mode by −2 for IMHA and 13 cm −1 for the fragment molecule (Table 2). Such frequency shifting is re lated to different molecular packing efficiency in bulk and SAM. The Im ring in IMHA molecule introduces sterical hindrances for the neighboring amides to engage in the for mation of an optimal H-bonding network.    Figure 6A presents evidence of Ag@SiO 2 functionality. The Raman spectrum of nanoparticles has virtually no vibrational modes until the nanoparticles are placed on top of smooth surface-adsorbed IMHA monolayer. The intense mode at 702 cm −1 corresponds to ν(C-S) T stretching vibration of the molecules, which adopt a nearly vertical orientation with the surface at the C-S bond, whereas the corresponding gauche mode ν(C-S) G appears as a weak feature at 622 cm −1 . The predominant trans configuration agrees with the vertical orientation of molecules at the amide group in mature SAM as revealed by RAIRS. Weak Am-I bands at 1685 and 1640 cm −1 were immediately ascribed to surface molecules that are involved in strong and weak H-bonding at amide groups, respectively. The rather complex 1500-1700 cm −1 region contains no Am-II mode because this mode is typically either very weak or not detectable in Raman spectra [56]. Two strong bands at 1062 and 1019 cm −1 were found to be not sensitive to H/D exchange and were assigned to stretching vibrations of hydrocarbon chain, ν(C-C) T and ν(C-C), respectively ( Table 1). The higher frequency component was assigned to the in-phase vibration of -(CH 2 ) 5 − chain in extended alltrans conformation based on DFT calculations (1089 cm −1 band) and previously reported temperature-dependent SERS studies of SAMs [57]. The 1019-cm −1 band was assigned to the C-C stretching vibration of hydrocarbon chain connecting amide and Im groups based on DFT suggestion (1018 cm −1 band).

General Features of IMHA Monolayer SHINERS Spectrum
vibrations of hydrocarbon chain, ν(C-C)T and ν(C-C), respectively ( Table 1). The higher frequency component was assigned to the in-phase vibration of -(CH2)5− chain in extended all-trans conformation based on DFT calculations (1089 cm −1 band) and previously reported temperature-dependent SERS studies of SAMs [57]. The 1019-cm −1 band was assigned to the C-C stretching vibration of hydrocarbon chain connecting amide and Im groups based on DFT suggestion (1018 cm −1 band). The tautomeric equilibrium of Au-adsorbed IMHA could certainly differ from that of the solution phase. Indeed, ν(C4=C5) + ν(C4-C6) + β(C5H) mode at 1582 cm −1 fall close to the 1588−1583 cm −1 range typical for T-II [16]. While in the case of T-I conformer, this mode is expected to be observed at considerably lower wavenumbers (1573-1568 cm −1 ). Cartoon depicts the selected imidazole ring motions.
The tautomeric equilibrium of Au-adsorbed IMHA could certainly differ from that of the solution phase. Indeed, ν(C4=C5) + ν(C4-C6) + β(C5H) mode at 1582 cm −1 fall close to the 1588−1583 cm −1 range typical for T-II [16]. While in the case of T-I conformer, this mode is expected to be observed at considerably lower wavenumbers (1573-1568 cm −1 ). The downshift of this mode to 1575 cm −1 in D 2 O solution immediately confirms H/D exchange at nitrogen sites of Im ring in the monolayer. Another couple of bands, at 1261 for H 2 O and 1257 cm −1 for D 2 O solutions, confirm the T-II surface form of IMHA [16].

Potential-Controlled SHINERS Measurements of IMHA Monolayer
Negative electrode polarization at −0.8 V potential was applied for 30 s in the beginning of each potential-controlled measurement to desorb impurities of low surface affinity, which may originate from nanoparticles. Preliminary examination showed that pretreatment of the monolayer at negative potentials increases the reproducibility of results. Figure 6B shows SHINERS spectra from the smooth Au electrode adsorbed IMHA at −0.8 and 0.2 V electrode polarizations.
Plenty of SERS research on metal-sulfur bond has been carried out for the roughened metal adsorbed monolayers [57,58]. However, only with SHINERS approach the welldefined and smooth substrate adsorbed molecules become more accessible to Raman spectroscopy [37,59]. The Au-S stretching mode found at 257 cm −1 upshifts to 279 cm −1 (δ = 22 cm −1 ) when the potential is tuned from −0.8 to 0.2 V attesting the increase in Au-S bond strength. A similar value was found in a previous study on flat Au-adsorbed N-(6-mercapto)hexylpyridinium [37], whose high tuning rate was ascribed to force constant changes in Au-S bond strength because of electrode polarization and to some extent to the Stark effect [37,60]. As for C-S bond, the populations of gauche and trans molecular conformers were matching at −0.8 V, but with positive electrode polarization trans became strongly favored ( Figure S4). Both conformational bands redshift with the tuning rate of −5.3 ± 0.2 and −6.9 ± 0.2 cm −1 V −1 , respectively, indicating potential-induced decrease in C-S bond strength, while at the same time Au-S bond became increasingly stronger at more positive potentials [61]. Within the tested potential window, the trans/gauche ratio (I 701 /I 633 ) increases from 1 to 3 with transition midpoint potential of −0.34 ± 0.02 V. Notably, all spectral changes are reversible as the electrode potential is set to −0.8 V. Two strong bands at 1051 and 1014 cm −1 were assigned to stretching vibrations of hydrocarbon chain, ν(C-C) T and ν(C-C), respectively (Table 1). It should be noted that the relative intensity of ν(C−C) T band increases at more positive electrode potentials together with growing ν(C−S) T mode. In the high frequency region of Figure 6B, symmetric and asymmetric stretching vibrations of methylene groups ν(CH 2 ) near 2858 and 2921 cm −1 were found (E = −0.8 V). These modes will be discussed in more detail later in the manuscript.
To interpret the rather complex midrange region that is occupied by the imidazole ring, amide group vibrations, and the deformations of the methylene groups, a separate H/D exchange experiment was performed for the IMHA monolayer. The H/D exchange process was accomplished at the open circuit potential (0.17 V) by exchanging phosphate buffer solution (PBS) (with 0.1 M Na 2 SO 4 , pH 6.9) to pure D 2 O (Figure 7). Amide-I spectral mode is found to be composed of low-and high-energy components at 1643 and 1676 cm −1 , both of which downshift by 7-8 cm −1 upon the exchange. We find that vibrational modes in the 1300-1610 cm −1 were sensitive to the solvent exchange. For example, the Tautomer-II related C=C stretching mode of the Im ring ν(C4=C5) clearly downshifts from 1590 to 1571 cm −1 [16]. The well-defined feature at 1492 cm -1 assigned to ν(C2-N3) coupled with β(C2H) downshifts to 1483 cm −1 in accordance with literature data [18]. The intense band at 1333 cm −1 develops in D 2 O solutions. This band was identified as a Tautomer-II marker band for N3D protonated histidine [62]. Interestingly, the mode at 1545 cm −1 becomes significantly pronounced at 0.2 V potential ( Figure 7A) and almost disappears due to the H/D exchange. Its assignment remains not fully clear because the contribution from the Am-II could be almost certainly ruled out since the deuteration of amide's N atom would introduce Am-II' with additional spectral intensity near 1450 cm −1 [63,64]. The only possible proximal spectral modes are expected near 1535 cm −1 for deprotonated 4-methylimidazole and doubly protonated histidine, both of which in neutral pH are not likely [13,62]. Presented SHINERS data revealed that the imidazole ring in the SAM at relatively positive electrode potentials is in the N3H protonation state ( Tautomer-II).  Figure 8 shows the detailed SHINERS analysis of biased potential induced structural changes in IMHA SAM. An 11.2 cm −1 upshift in the frequency of dominant Amide-I component is detected with more positive electrode potentials ( Figure 8C). To complement that, we analyzed the relative intensities of Am-I modes of strongly (ca. 1631 cm −1 ) and loosely H-bonded (ca. 1680 cm −1 ) molecules. The intensity ratio of these expressed as I1631/I1680 monotonically decreased from 3.4 to 2.3 with the negative to positive potential excursion. Both frequency shift and relative intensity changes show that a positive bias potential disengages molecules from the more strongly H-bounded network at the C=O moiety of amide groups.   Figure 8C). To complement that, we analyzed the relative intensities of Am-I modes of strongly (ca. 1631 cm −1 ) and loosely H-bonded (ca. 1680 cm −1 ) molecules. The intensity ratio of these expressed as I 1631 /I 1680 monotonically decreased from 3.4 to 2.3 with the negative to positive potential excursion. Both frequency shift and relative intensity changes show that a positive bias potential disengages molecules from the more strongly H-bounded network at the C=O moiety of amide groups. As Figures 6 and 8 show, the positions of imidazole-tautomerism marker bands near 1262, and 1597 cm −1 remain consistent with N3-H protonation (Tautomer-II) within the tested potential window [19]. However, at −0.8 V potential, a shoulder at 1573 cm −1 associated with Tautomer-I protonation form of Im becomes visible. Positive polarization diminishes mode intensity, while its T-II counterpart at 1597 cm −1 becomes stronger. Thus, SHINERS data show the T-II being the preferred tautomeric form at each tested potential, with a minor portion of molecules adopting T-I at the most negative polarizations. Potential induces a slight shift of the 1597 cm −1 to lower wavenumbers by 3.6 cm −1 due to lengthening of the C4=C5 bond. Interestingly, the Im ring mode at 1490 cm −1 assigned to C2-N3 stretching coupled with in-plane deformation β(C2H), ν(C2-N1), and ν(C5-N1) shows the nonmonotonous dependency on potential [10,18]. This mode shifts by −4.9 cm −1 (tuning rate 13.3 cm −1 V −1 ) in the range from −0.8 to −0.4 V and then upshifts by 1.8 cm −1 ( Figure  8D). Besides, the spectral mode notably increases in FWHM by ca. 10 cm −1 of the −0.4 V As Figures 6 and 8 show, the positions of imidazole-tautomerism marker bands near 1262, and 1597 cm −1 remain consistent with N3-H protonation (Tautomer-II) within the tested potential window [19]. However, at −0.8 V potential, a shoulder at 1573 cm −1 associated with Tautomer-I protonation form of Im becomes visible. Positive polarization diminishes mode intensity, while its T-II counterpart at 1597 cm −1 becomes stronger. Thus, SHINERS data show the T-II being the preferred tautomeric form at each tested potential, with a minor portion of molecules adopting T-I at the most negative polarizations. Potential induces a slight shift of the 1597 cm −1 to lower wavenumbers by 3.6 cm −1 due to lengthening of the C4=C5 bond. Interestingly, the Im ring mode at 1490 cm −1 assigned to C2-N3 stretching coupled with in-plane deformation β(C2H), ν(C2-N1), and ν(C5-N1) shows the nonmonotonous dependency on potential [10,18]. This mode shifts by −4.9 cm −1 (tuning rate 13.3 cm −1 V −1 ) in the range from −0.8 to −0.4 V and then upshifts by 1.8 cm −1 ( Figure 8D). Besides, the spectral mode notably increases in FWHM by ca. 10 cm −1 of the −0.4 V spectrum compared to the −0.8 V one. From temperature-Raman measurements, decreasing mode frequency and increasing FWHM were linked with liberated imidazole ring motion. Thus, spectral data show that Im ring confinement in the monolayer at more positive electrode potentials is relaxed. The strong band at 1262 cm −1 (−0.8 V) due to Im ring breathing vibration coupled with C2-H in-plane deformation is very sensitive to potential perturbation; frequency of this band upshifts to 1266 cm −1 at electrode polarization 0.2 V. Consequently, the potential-difference spectrum clearly shows derivative-like feature at 1279/1254 cm −1 ( Figure 8B). Temperature-dependent Raman study of bulk IMHA compound revealed an upshift in the position of this mode at elevated temperature (135 • C) (Figure 4). This might be related to decreased H-bonding interaction strength at the Im ring site. Such spectral observations are consistent with decreased Hbonding interaction strength at both amide (Am-I band) and Im ring sites at more positive electrode polarization.

Raman Markers for H-Bonding Interaction
In order to better understand the way H-bonding interaction affects imidazole ring structure and its vibrational frequencies, DFT modeling was carried out for 5-ethyl-1Himidazole (Im-CH 2 -CH 3 ) molecule isolated in vacuum and coordinated with one and two explicit water molecules ( Figure 9). Calculations of T-II form molecule predict that 1608-cm -1 mode is mainly C4=C5 stretching motion; 1499 cm −1 mode is assigned to ν(C2-N1) + ν(C2-N3) + β(C2H) and associated with 1492 cm −1 mode in SHINERS spectra. The 1423 cm -1 mode emerges due to β(N3H) + ν(N3-C4) + ν(N3-C2) vibrations; 1382 cm −1 mode is tautomerism sensitive Im breathing mode. With increasing H-bonding coordination number, the 1499-cm −1 mode shifts to higher frequencies by 2-6 cm −1 (coordination number 1) and by 9 cm −1 (coordination number 2). The same holds for a compound in T-I form, for which an 8 cm -1 shift was found when in contact with two H 2 Os. The predicted wavenumber shift confirms the principal mode's sensitivity to H-bonding interaction at the Im group and agrees with the temperature-Raman data presented in Figure 4 where the higher wavenumbers were associated with stronger hydrogen bonding interactions for the sample at room temperature. We already ascribed the decrease in position of 1492 cm −1 mode with liberated Im ring based on temperature-Raman. However, DFT proves that such a liberation is hydrogen bonding strength-related. Thus, the sharp wavenumber drop from −0.8 to −0.4 V potentials in Figure 8D clearly indicates reduction of H-bonding at Im, which is followed by a slight strengthening in −0.4/0.2 V range.
Interestingly, we found stronger H-bonding interactions for Im coordinated with two H 2 O molecules rather than one H 2 O. In this case, the length of N1···H 2 O and N3-H···H 2 O bonds decreased by 1.9 and 2.7 pm with an introduction of the second H 2 O molecule (insets in Figure 9). A much higher, up to 30 cm −1 , frequency shift was found for 1423-cm −1 mode. However, this mode was difficult to identify in experimental SHINERS spectra of IMHA SAM because of the overlapping with methylene deformation vibrations.  Interestingly, we found stronger H-bonding interactions for Im coordinated with two H2O molecules rather than one H2O. In this case, the length of N1···H2O and N3-H···H2O bonds decreased by 1.9 and 2.7 pm with an introduction of the second H2O molecule (insets in Figure 9). A much higher, up to 30 cm −1 , frequency shift was found for 1423-cm −1 mode. However, this mode was difficult to identify in experimental SHINERS spectra of IMHA SAM because of the overlapping with methylene deformation vibrations.

Conclusions
Hydrogen bonding is fundamental in protein architecture, interactions, and molecular recognition. Imidazole group of histidine amino acid is a major player in active centers of enzymes and protein secondary structure because of its ability to accept and donate Hbonding. To study H-bonding interactions of imidazole ring at electrified interfaces, first, we synthesized alkanethiol molecule with interchain amide group and terminal imidazole functional group (IMHA), then studied IMHA SAMs adsorbed on atomically smooth Au surfaces by means of SHINERS and RAIRS techniques.
(i) From RAIRS data molecules adsorb on a surface chaotically but with time (transition midpoint was 162 min), their intrachain amide groups become nearly perpendicular to the surface, which points to the more neatly organized SAM. However, the formation of an optimal H-bond network among amide groups was precluded by the steric hindrances introduced by bulky IMHA's imidazole group, as revealed by a study of the imidazole-truncated fragment molecule. For further endeavors, mixed IMHA/fragment molecule SAMs are suggested.
(ii) Adsorption on a surface induced a tautomeric transition in IMHA molecules from Tautomer-I (N1-H, N3) to Tautomer-II (N1, N3-H). The Tautomer-II was strongly favored at electrode polarization close to open circuit potential. However, an increased portion of Tautomer-I molecules was observed at −0.8 V.
(iii) Temperature-Raman and DFT modeling identified the 1492 cm −1 mode as sensitive to imidazole ring confinement and H-bonding strength. The analysis of this mode revealed the potential-dependent behavior of the interfacial Im ring, such that the Im is

Conclusions
Hydrogen bonding is fundamental in protein architecture, interactions, and molecular recognition. Imidazole group of histidine amino acid is a major player in active centers of enzymes and protein secondary structure because of its ability to accept and donate H-bonding. To study H-bonding interactions of imidazole ring at electrified interfaces, first, we synthesized alkanethiol molecule with interchain amide group and terminal imidazole functional group (IMHA), then studied IMHA SAMs adsorbed on atomically smooth Au surfaces by means of SHINERS and RAIRS techniques.
(i) From RAIRS data molecules adsorb on a surface chaotically but with time (transition midpoint was 162 min), their intrachain amide groups become nearly perpendicular to the surface, which points to the more neatly organized SAM. However, the formation of an optimal H-bond network among amide groups was precluded by the steric hindrances introduced by bulky IMHA's imidazole group, as revealed by a study of the imidazoletruncated fragment molecule. For further endeavors, mixed IMHA/fragment molecule SAMs are suggested.
(ii) Adsorption on a surface induced a tautomeric transition in IMHA molecules from Tautomer-I (N1-H, N3) to Tautomer-II (N1, N3-H). The Tautomer-II was strongly favored at electrode polarization close to open circuit potential. However, an increased portion of Tautomer-I molecules was observed at −0.8 V.
(iii) Temperature-Raman and DFT modeling identified the 1492 cm −1 mode as sensitive to imidazole ring confinement and H-bonding strength. The analysis of this mode revealed the potential-dependent behavior of the interfacial Im ring, such that the Im is the most strongly confined and H-bonded at −0.8 V electrode polarization and the weakest confined and H-bonded at −0.4 V.
The presented data on construction, structure, and potential dependence of imidazole ring group-functionalized SAMs will possibly be valuable for the development of biosensors, molecular electronics, environmental, and pollution studies.