Takakazu
Seki‡
a,
Chun-Chieh
Yu‡
a,
Xiaoqing
Yu
a,
Tatsuhiko
Ohto
b,
Shumei
Sun
acd,
Konrad
Meister
ae,
Ellen H. G.
Backus
ad,
Mischa
Bonn
*a and
Yuki
Nagata
*a
aMax Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany. E-mail: bonn@mpip-mainz.mpg.de; nagata@mpip-mainz.mpg.de
bGraduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
cDepartment of Physics, Applied Optics Beijing Area Major Laboratory, Beijing Normal University, Beijing 100875, China
dDepartment of Physical Chemistry, University of Vienna, Währinger Strasse 42, 1090 Vienna, Austria
eDepartment of Natural Sciences, University of Alaska Southeast, Juneau, AK 99801, USA
First published on 23rd April 2020
The structure of interfacial water determines atmospheric chemistry, wetting properties of materials, and protein folding. The challenge of investigating the properties of specific interfacial water molecules has frequently been confronted using surface-specific sum-frequency generation (SFG) vibrational spectroscopy using the O–H stretch mode. While perfectly suited for the water–air interface, for complex interfaces, a potential complication arises from the contribution of hydroxyl or amine groups of non-water species present at the surface, such as surface hydroxyls on minerals, or O–H and N–H groups contained in proteins. Here, we present a protocol to extract the hydrogen bond strength selectively of interfacial water, through the water bending mode. The bending mode vibrational frequency distribution provides a new avenue for unveiling the hydrogen bonding structure of interfacial water at complex aqueous interfaces. We demonstrate this method for the water–CaF2 and water–protein interfaces. For the former, we show that this method can indeed single out water O–H groups from surface hydroxyls, and that with increasing pH, the hydrogen-bonded network of interfacial water strengthens. Furthermore, we unveil enhanced hydrogen bonding of water, compared to bulk water, at the interface with human serum albumin proteins, a prototypical bio-interface.
Among them, vibrational SFG spectroscopy has been demonstrated to be a powerful technique to probe the interfacial water molecules at soft (biomolecule–water)2–7 and buried (solid–water) interfaces.8–17 SFG spectroscopy is active only in the centrosymmetry-broken region and therefore probes those molecular responses with a non-zero net transition dipole moment along the surface normal direction. As such, the molecules in the bulk and the molecular responses with their transition dipole moment being parallel to the surface are SFG-inactive. However, capturing the net contribution of the interfacial water in the O–H stretch spectra is often highly challenging at the complex aqueous interfaces, because the contribution of the water O–H stretch mode is spectrally indistinguishable from that of the O–H groups chemisorbed onto the solid surface as well as the O–H and N–H groups of proteins. A small number of these O–H/N–H groups can overwhelm the SFG signal, because the SFG signal of water arises from a very thin (sub-nm) water layer at the interfaces. These comparable contributions from the water O–H groups and the other hydroxyl groups impede a clear molecular view of the interfacial water structure. For instance, a weakly hydrogen-bonded O–H stretch peak in the SFG spectra at the water-silica interface has been conflictingly assigned to originating from interfacial water18,19 and silanol groups on the silica surface.15 As such, an unambiguous determination of the interfacial water contribution is essential for unveiling the microscopic structure of water at interfaces. A second complication is the distortion of the SFG spectrum of the O–H stretching mode due to intra- and intermolecular vibrational coupling effects.20
The net interfacial water contribution may be measured through the water bending mode, because the H–O–H bending mode originates solely from the H2O water molecules and has a spectrally distinguishable frequency from other O–H bending modes. For example, the H–O–H bending mode frequency is 1590–1650 cm−1, which is significantly distinct from the bending mode frequency of hydroxyl groups attached on solid surfaces (300–900 cm−1).21 Since the first publication of the sum-frequency generation spectrum in the frequency region of the bending mode from the water–air interface,22 there have been debates about the surface- versus bulk-contribution to the signal.23,24 We have recently demonstrated that the surface contribution to the overall spectrum can be readily extracted, paving the way for using the bending mode as a reporter for determining the interfacial water structure.24 The water bending mode may also be well-suited to study other complex interfaces such as protein interfaces. Although the strong amide I mode of proteins (1620–1700 cm−1) masks the H–O–H bending mode of water, the H–O–D bending mode of the isotopically diluted water (1400–1450 cm−1) can provide the information on the interfacial water.
A drawback to probe the water bending mode with SFG spectroscopy lies in its weak signal. The challenge in probing the bending mode SFG signal is further increased, when the bending mode signal is weakened by using the isotopically diluted water; if the bending mode lineshape is distorted by the intermolecular vibrational coupling of the bending modes, one has to obtain the characteristic frequency of the bending mode through the isotopic dilution of water. However, knowledge of the vibrational coupling on the bending mode in the SFG spectra is completely lacking. Thus, to establish a protocol to probe the interfacial water contribution through the bending mode SFG, one should first obtain the information on whether the bending mode SFG response is affected by vibrational coupling, by performing experiments with isotopic dilution.
Here, we demonstrate that the lineshape of the SFG bending mode signal is not changed upon the isotopic dilution from the SFG measurement at the water–air interface. We further compute the SFG spectra using the ab initio molecular dynamics (AIMD) simulations and confirm that the intermolecular coupling of the bending modes does not alter the lineshape of the bending mode. Subsequently, we obtain the net interfacial water contribution through the water bending mode at two aqueous interfaces. One is the water–CaF2 interface, where we decode the contribution of the interfacial water and find that the hydrogen bond strength is increased with increasing pH. The other is the water–air interface in the presence of the human serum albumin (HSA) protein, where we use the H–O–D bending mode and discuss the strength of the water–HSA protein binding interaction. The successful determination of the interfacial hydrogen bond network through the SFG water bending mode manifests that probing the water bending mode in the SFG spectra is a promising tool for understanding the interfacial water structure.
Fig. 1(b) and (c) show that ωbend is independent of the H2O concentration, while Γbend shows a weak broadening with increasing H2O concentration. Fig. 1(d) shows that Abend is linearly proportional to the H2O concentration. Here, the H2O concentration was inferred using the equilibrium relation:28–30
[HOD]2/[H2O][D2O] = 3.86. | (1) |
The dotted lines in Fig. 1(b–d) show the response expected in the absence of vibrational coupling. Clearly, we find no evidence for substantial vibrational coupling between bending modes of interfacial water.
SFG spectral calculations using AIMD trajectories at the B3LYP-D3(0) and revPBE0-D3(0) levels of theory confirm this conclusion. The simulated imaginary part spectra of χ(2)R, Im(χ(2)R), with various cross-correlation cutoff distances rt are displayed in Fig. 2. An increase in rt means that the simulated spectra contain the effects of the intermolecular vibrational coupling between the distant water molecules. These data show that Im(χ(2)R) is insensitive to the cross-correlation cutoff rt, confirming the notion that the effects of the vibrational coupling on the bending mode spectrum are negligible and thus the SFG lineshape of the water bending mode is not distorted by the isotopic dilution. This ensures that the characteristic frequency of the water bending mode can be obtained with the arbitrary H2O–D2O concentrations. This is in stark contrast to the O–H stretch mode, for which the peak is strongly distorted through the vibrational coupling of the O–H stretch modes20,31–33 and thus the characteristic frequency of the O–H stretch mode can be obtained only for substantial isotopically diluted water (e.g. 5% O–H in D2O).
This finding is important for several aspects. First, since the characteristic frequency of the bending mode can be obtained even from the pure H2O, we do not have to reduce the bending mode signal further by isotopic dilution to obtain the characteristic frequency. This is a key for probing the SFG bending mode, as the SFG bending mode signal of water is weak and thus the detection of the bending mode signal is often challenging. Second, not only the H–O–H bending mode but also the H–O–D bending mode can be used for probing the interfacial water contribution; because the H–O–D bending mode appears only when investigating the H2O–D2O mixtures, the H–O–D bending mode can be obtained through the systematic change of the H2O–D2O mixture composition. If the lineshape is changed through the concentration of H–O–D, the H–O–D bending mode contribution is not accessible. Fortunately, since the bending mode SFG lineshape is insensitive to the intermolecular vibrational couplings, one can safely assume that the H–O–D bending lineshape is unchanged. Investigating the H–O–D bending mode is particularly relevant for bio-interfaces, since the H–O–H bending mode overlaps spectrally with the amide mode of proteins. Below, we demonstrate how measuring the water bending mode SFG feature with the tuned H2O–D2O concentrations can reveal the hydrogen-bonding strength of interfacial water at water–CaF2 and water–protein interfaces.
We measured the SFG intensity (I(2)) spectra for the H–O–H bending mode of the H2O–CaF2 interface. These are displayed in Fig. 3(a). The intensity spectra show a peak-like feature at pH 2, while a dip feature was observed around 1650–1750 cm−1 at pH 12 due to the destructive interference between the non-resonant contribution and the bending mode signal. The peak-to-dip transition with pH increasing from 2 to 12 indicates that the sign of the water bending mode changes. This is attributed to flipping of the water orientation due to the change of the sign of the CaF2 surface charge. Note that the SFG signal at D2O–CaF2 does not show any resonant feature. The fit of the intensity spectra (see the ESI†) indicates that the peak frequency of the H–O–H bending mode band at pH = 2 is 1633 ± 2 cm−1, while it is 1662 ± 3 cm−1 at pH = 12 in the Im(χ(2)) spectra. Since a higher frequency of the bending mode represents a stronger hydrogen bond,24,34,35 the variation of the bending mode frequency indicates that the hydrogen bond strength of the interfacial water at the water–CaF2 interface enhances as the pH increases.
To relate the H–O–H bending mode with the O–H stretch mode, we measured the O–H stretch modes of HDO at the isotopically diluted water–CaF2 interface. The spectra are displayed in Fig. 3(b). Note that, to correlate the vibrational coupling-absent H–O–H bending mode spectra with the O–H stretch spectra which are largely modulated by the vibrational couplings,20,36,37 we used isotopically diluted water for the O–H stretch mode to suppress the vibrational couplings. The following empirical relation between the bending mode and stretch mode frequencies;24,35,38
ωgstr − ωstr = a(ωgbend − ωbend)2 + b(ωgbend − ωbend), | (2) |
Eqn (2) indicates the expected O–H stretch peak frequency, concluded from the bend frequency of interfacial water at pH = 2, is 3503 ± 9 cm−1, which agrees with the O–H stretch mode frequency obtained from the fit of the intensity spectra of 3426 ± 3 cm−1. Since the Ca–O–H groups are expected to be absent at pH 2,10,11 the O–H stretch signature can be reasonably considered to arise solely from the interfacial water. Because the O–H stretch mode is much stronger than the non-resonant contribution unlike the H–O–H bending mode, the non-resonant contribution is nearly zero for the O–H stretch spectra, while the non-resonant contribution is certainly large for the H–O–H bending mode spectra (see the ESI†). Subsequently, we focus on the case of pH 12. The H–O–H bending mode frequency of 1662 cm−1 suggests an expected O–H stretch mode frequency of 3320 ± 12 cm−1. This expected O–H stretch frequency is in good agreement with the peak frequency of 3363 ± 15 cm−1 for the 3000–3500 cm−1 O–H stretch band. This suggests that the 3000–3500 cm−1 O–H band arises from the interfacial water, not the Ca–O–H groups, consistent with the previous studies.10,11 The results also demonstrate that the 3650 cm−1 peak does not originate from water, but indeed from the Ca–O–H groups. As such, the comparison of the O–H stretch SFG spectra with the bending mode spectra can offer an unambiguous assignment of the vibrational signatures at complex aqueous interfaces.
Here, we demonstrate that the water contribution can be extracted from the bending mode SFG spectra. Probing the H–O–H bending mode in the presence of protein is challenging, because the H–O–H bending mode peak (1590–1650 cm−1) is overwhelmed by the strong amide I mode peak (1620–1700 cm−1). To avoid the frequency overlap with the amide I mode, we probe the H–O–D bending mode (1450–1500 cm−1) instead of the H–O–H bending mode. Although the H–O–D bending mode frequency overlaps with the frequency of the amide II mode (1480–1570 cm−1),41 the amide II mode signal is small in the SFG spectra, due to the small Raman cross-section of the amide II mode.42 We would like to note that, to properly disentangle the H–O–D bending mode contribution from the potentially overlapped amide II mode, and to verify the procedure, one needs to carefully compare the amplitude of the bending mode against the H–O–D concentration.
We measured the SFG intensity (I(2)) at the isotopically diluted water–air interface in the presence of HSA interfaces by varying the H2O/D2O concentration. Fig. 4(a) displays the variation of the SFG spectra. In the absence of HDO – i.e. for pure H2O and D2O-, the SFG spectra are indistinguishable and show a signal at 1410 cm−1 and a relatively weak signal between 1500 and 1550 cm−1, attributed to the H–C–H bending mode or the symmetric stretch mode of R-COO−, and the amide II mode of HSA, respectively. These spectral features closely resemble the SFG spectra obtained at pH 6.3 for the water–air interface with bovine serum albumin.43 To extract the H–O–D bending mode contribution from the spectra for the H2O/D2O mixtures, we fit the spectra (see the ESI†). The obtained imaginary parts of the response are displayed in Fig. 4(b), with the H–O–D bending contributions highlighted, while the amplitude of the H–O–D contribution is plotted in Fig. 4(c). The variation of the amplitude shows excellent agreement with the H–O–D concentration inferred from eqn (1), confirming the peak assignment. The center frequency for the H–O–D contribution is 1506 ± 3 cm−1. This frequency is substantially higher than the frequency of the H–O–D bending mode in the bulk (∼1460 cm−1),44 indicating that the interfacial water molecules near the HSA protein are more strongly hydrogen-bonded than water in the bulk.
![]() | ||
Fig. 4 SFG responses at the water–air interface in the presence of HSA. (a) SFG intensity spectra, I(2)(ω), of the water–air interface in the presence of HSA with 1/0, 3/1, 1/1, 1/3, and 0/1 H2O/D2O mixtures. The black lines are the fit. The spectra are offset by 0.05 for clarity. (b) Imaginary part spectra of the HSA and H–O–D water contributions obtained from the fits. The H-O-D contribution is highlighted in blue. The spectra are offset by 0.25 for clarity. (c) The amplitude of the H–O–D contribution vs. H2O fraction. Theoretical prediction denoted by the black broken line is obtained from eqn (1). |
For example, 1,2-dipalmitoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (DPPG) is an essential lipid for a cell membrane, and the DPPG–water interfaces have been studied with O–H stretch SFG to understand the water orientations.45–47 However, a DPPG lipid has multiple O–H groups in its headgroup, prohibiting the extraction of the net interfacial water contribution. Having established here that the characteristic frequency of the bending mode can be obtained without any isotopic dilution, it is apparent from our previous bending SFG data at the water–DPPG interface24 that the hydrogen bond strength of water near DPPG is similar to that in the bulk. We observed that the characteristic frequency of the interfacial H–O–H bending mode is 1653 cm−1, i.e., very similar to 1650 cm−1 in the bulk liquid. The conclusion that water at a DPPG interface is hydrogen-bonded with bulk-like strength also affects the interpretation of the vibrational dynamics of the O–H stretch mode near DPPG, reported to be much slower than the bulk water vibrational dynamics.55
We point out that the bending mode SFG probe will be powerful for unveiling the interfacial structure of water at the interface of water–alcohol mixtures. This system is important for dilution processes as well as general physical chemistry, because upon mixing with water, the structure of the water and alcohol mixture is drastically changed, resulting in the variation of the vapor pressure. The microscopic structures of the water–alcohol interfaces have been studied with SFG in the O–H stretch region,48–52 but the net contribution of the interfacial water has not been accessed. The different contributions to the O–H stretch spectrum can now be disentangled by probing the water bending mode in the SFG measurements.
The water bending mode signal is useful not only because its frequency carries information on the hydrogen bond strength of the interfacial water molecules, but also because its amplitude is sensitive to the chemical compositions. Theoretical studies have suggested that, at the water–TiO2 interface, chemisorbed and physisorbed conformations of the interfacial water molecules17 are easily distinguished through the bending mode amplitude.56 As such, for several cases, the water bending mode can be a very sensitive probe for the microscopic structure of the interfacial water molecules.
In summary, our combined simulations and experiments showed that the SFG feature of the water bending mode of interfacial water at the water–air interface is not distorted upon the isotopic dilution of water (Fig. 1) and thus the vibrational coupling of the bending modes of the interfacial water is virtually absent (Fig. 2). Thus, the characteristic frequency of the bending mode can be obtained for an arbitrary mixture of the H2O–D2O concentration, unlike the O–H stretch mode, where the characteristic frequency can be obtained only for the isotopically diluted water.36 Since the use of the isotopic diluted water reduces the SFG signal drastically, our findings ensure that one can use pure H2O for gaining the maximum SFG amplitude or can use the H–O–D bending mode for avoiding frequency overlap with the amide I mode. We have applied this new approach for exploring the CaF2–water interface, and found that an increase in pH enhances the hydrogen bond strength of interfacial water (Fig. 3). Even for a complex protein–water interface, the bending mode SFG could successfully probe the interfacial water contribution in the H–O–D bending mode region and revealed that the interfacial water near the HSA protein shows a stronger hydrogen bond strength than bulk water (Fig. 4). Altogether, our work manifests that the water bending mode is a promising tool for understanding interfacial water structures at various complex interfaces.
![]() | (3) |
Footnotes |
† Electronic supplementary information (ESI) available: Details of the SFG measurements, the fitting procedures, and the AIMD protocols. See DOI: 10.1039/d0cp01269f |
‡ These authors contributed equally. |
This journal is © the Owner Societies 2020 |