Yu
Wu
abc,
Zhongjie
Zhu
*b,
Jinrong
Yang
d,
Jie
Wang
b,
Te
Ji
b,
Huachun
Zhu
b,
Weiwei
Peng
b,
Min
Chen
ab and
Hongwei
Zhao
*ab
aShanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China. E-mail: zhaohw@sari.ac.cn
bShanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China. E-mail: zhuzj@sari.ac.cn
cUniversity of Chinese Academy of Sciences, Beijing 100049, China
dEast China Normal University, Shanghai 200241, China
First published on 12th July 2024
L-Glutamic acid (L-Glu) is a basic unit of proteins and also serves as an important neurotransmitter in the central nervous system. Its structural properties are critical for biological functions and selective receptor recognition. Although this molecule has been extensively studied, the low frequency vibrational behavior that is closely related to conformational changes and the intermolecular interactions between L-Glu and its receptors are still unclear. In this study, we acquired the fingerprint spectrum of L-Glu by using air plasma terahertz (THz) time-domain spectroscopy in the 0.5–18 THz range. The low frequency vibrational characteristics of L-Glu were investigated through density functional theory (DFT) calculations. The THz responses of the ligand binding domain of the NMDAR–L-Glu complex were studied by the ONIOM method, with a focus on discussing the normal modes and interactions of ligand L-Glu and water molecules. The results illustrate that THz spectroscopy exhibits a sensitive response to the influence of L-Glu on the structure of the NMDAR. The water molecules in proteins have various strong vibration modes in the THz band, showing specificity, diversity and complexity of vibrational behavior. There is potential for influencing and regulating the structural stability of the NMDAR–L-Glu complex through water molecules.
During the process of the protein–ligand complex exerting biological functions, water molecules usually play an essential role, often participating in the change of configuration and conformation as part of the structure and affecting biological activity. The interactions between water molecules and biomolecules are complex and subtle.8–10 Water molecules have the dual ability to act as both hydrogen bond donors and acceptors and can form rich hydrogen bond networks.11,12 They are sometimes regarded as “biomolecules” with a basic structure and function.13 Water molecules located in the internal cavities and cracks of proteins can affect the shape and flexibility of proteins as well as stabilize the structure between the protein and the ligand, improve spatial complementarity,11,14 regulate interactions, and increase affinity.8 Vance et al. studied the structures of the ionotropic glutamate receptor NMDAR in complex with four ligands including L-glutamate and its analogues. Compared with other agonists, the time of receptor inactivation caused by L-Glu is significantly slower, which is due to the unique conformation induced by L-Glu at the backside of the ligand binding site. An atomic-resolution crystallographic study shows that in the complex structure of the NMDAR with L-Glu, there are numerous water molecules embedded, with one water molecule bridging the ligand L-Glu and residues through hydrogen bonding, which is critical in the ligand binding domain.15 However, due to the complexity of the microenvironment in the L-Glu–NMDAR complex, the characteristics of the hydrogen bond network established by water molecules, the impact of variations in water molecules on receptor conformation, and the specificity of interactions involving water molecules remain uncertain.
In protein–ligand complexes, weak interactions such as hydrogen bonding and van der Waals forces play an important role in receptor–ligand recognition and maintaining the structure of biological molecules. Recent studies have shown that these weak interactions are widely distributed in the THz far-infrared region and exhibit unique properties. The low frequency vibrations are dominated by intermolecular interactions, which mainly manifest as collective vibrations involving more atoms, and multiple vibrations are mixed together, which are highly sensitive to intramolecular and intermolecular structures, thus providing characteristic fingerprint information on molecular conformational states.16–18 Terahertz time-domain spectroscopy (THz-TDS) exhibits a highly sensitive response to the conformational changes of molecules, enabling effective detection of the physical and chemical properties of materials and facilitating non-destructive testing of biomolecules in the THz band.19–24 Ruggiero et al. used THz-TDS to distinguish the unique fingerprint spectra of two crystal types of L-Glu in the 0.5–4.0 THz region and evaluated the quality of the crystal intermolecular force model.25 Li et al. found that L-Glu and its structural analog exhibit different vibrational characteristics through THz-TDS and DFT calculations in the 0.5–3.0 THz region, which can be attributed to their distinct hydrogen bonding network structures.26 Hand et al.27 suggested that some of the capabilities of biomolecules depend on their sampling of different conformational states, which are determined by the vibrational dynamics of macromolecules. Broad band THz spectroscopy provides a wealth of low-frequency vibrational information about molecules, and many biomolecules exhibit rich spectra near 6 THz, corresponding to the densest transitions between states. Hence, richer information in the entire THz range is helpful for a comprehensive understanding of biomolecules. The spectral resolution and fingerprint characteristics of THz spectroscopy make it useful for label-free detection in biochemical research.
THz photons have low energy and usually do not cause ionization damage to biomolecules.28 Therefore, THz electromagnetic waves are considered to be relatively safe for living organisms.29 However, some recent research results suggest that the resonant coupling between THz photons and molecules may produce some effects, even unexpected ones, indicating that the biological effects of THz radiation cannot be ignored.30–33 Song et al. revealed the potential mechanism of THz photons released by ion oscillation in the potassium ion channel and their resonant and coherent coupling with other channel oscillations, demonstrating the macroscopic coherent state of multiple ion channels.34 Wu et al. found that THz photons facilitate stretching vibrations on the purine plane, weakening hydrogen bonds near base pairs, thereby accelerating DNA unwinding in an efficient and non-thermal way.35 Zhu et al. successfully induced one-dimensional confined water to transform from a normal permeation state to a superpermeation state by using specific frequency THz waves resonating with water in a confined space and proved that the effect is non-thermal, suggesting the potential application of THz technology in biomedicine.36 Li et al. showed that the emission of THz waves resonating with the –COO− and –CO functional groups in voltage-gated calcium channels can significantly enhance the selectivity and conductivity of Ca2+ and found that an electromagnetic wave with a frequency of 4.0 THz can expand the exit of the risperidone ligand in the dopamine D2 receptor through conformational modulation and accelerate risperidone dissociation, thereby reducing the toxic side effects of drugs.31,37 As a possible neural information processing mechanism, the generation and signal transmission of biophotons in the nervous system have been supported by a series of studies.38–41 However, because of biological complexity, the mechanisms by which THz electromagnetic waves interact with biomolecules remain poorly understood. The inherent correlations between low frequency vibrations and interactions associated with conformational changes in protein–ligand complexes are still elusive. As an important amino acid and neurotransmitter, despite extensive work on L-Glu, little is known about the intermolecular interactions between this ligand and its receptors.
In this paper, L-Glu and its receptor NMDAR are studied as a model to investigate the signal response in the 0.5–18 THz band. The rich absorption peaks and unique normal modes of L-Glu reveal the characteristics of the interaction between L-Glu and THz electromagnetic waves. In protein–ligand complexes, water molecules dominate THz absorption in multiple frequency bands and are associated with protein and ligand binding. This study may provide a reference for the future use of THz waves to regulate the interaction between the ligand and the receptor.
![]() | ||
Fig. 1 Experimental spectrum (solid red line) and the calculated spectrum (blue dotted line) of L-Glu in the frequency range of 0.5–18 THz. The absorption in the grey area is for reference only. |
Modes with no infrared activity have no absorption peaks in the THz spectrum, but that does not imply their nonexistence. Therefore, both infrared and non-infrared active modes are considered. Fig. 2a shows the intramolecular vibrational contribution of L-Glu in the range of 0.5–18 THz, in which red dots indicate infrared activity modes and blue dots indicate non-infrared activity modes. It can be seen from the figure that the intramolecular vibrational contribution increases with the increase of frequency. This indicates that the molecular vibrational distribution with frequency in the studied region of 0.5–18 THz has certain regularity,53–55 which is related to the resonance response between different interactions and THz electromagnetic waves. The resonance frequency is positively correlated with the interaction strength. Therefore, the intermolecular interaction response is significant at low frequencies, with a relatively high contribution.
![]() | ||
Fig. 2 (a) Changes in the intramolecular vibrational contribution and (b) intermolecular vibrational contribution of normal modes of L-Glu in the range of 0.5–18 THz with frequency. |
In order to further clarify the relationship between molecular vibration and THz frequency, we used the function P = Ae−λF to fit the relationship between the intermolecular vibrational contribution and frequency, where P represents the intermolecular vibrational contribution, A is a constant, e is a natural constant, F is the frequency, and λ is a constant, which is related to the decay rate of the function. We found that the λ value of all normal modes (green dots in Fig. 2b) is 0.094, and P = 79.13e−0.094F is used for fitting. The λ value of infrared active normal modes (red circle in Fig. 2b) is 0.096, and P = 80.85e−0.096F is used for fitting. The λ value of normal modes without infrared active modes (blue circle in Fig. 2b) is 0.092, and P = 76.43e−0.092F is used for fitting. We used 1/λ to represent the upper limit of the frequency effectively dominated by the intermolecular interaction as a range indicator for evaluating the THz wave response of the intermolecular interaction, which is about 10.63 for all modes, further confirming the uniqueness of the matter's vibration in the THz range. The results show that the low frequency normal modes of THz are complex, but have certain regularity. The intermolecular vibration of biomolecules is dominated by low frequency collective vibration, which involves a large number of atoms and a series of intermolecular interactions. The hydrogen bond plays a major role in the weak intermolecular interaction, and the frequency of vibration led by the hydrogen bond has a large distribution in the THz range.56 In combination with the characteristic frequency distribution of the covalent bond in the mid-infrared band,57 it explains, to a certain extent, the phenomenon that the contribution of intermolecular vibration decreases with the increase of frequency.
There are two obvious differences between the spectra of NMDAR–L-Glu and NMDAR at 5.97 THz and 10.50 THz. Their corresponding normal modes are shown in Fig. 3c and d, respectively. The vibrations are collective and mixed with each other, showing complexity. Among them, the vibrations of water molecules are more significant, and the vibrations related to water are highlighted by the green arrow. Considering the special role of water, we focus on the response of water molecules to THz waves in the NMDAR–L-Glu complex. In the NMDAR, the normal mode at 5.97 THz frequency is mainly the rocking vibration of w2 (left panel in Fig. 3c). When L-Glu acts on the receptor, the vibration is significantly inhibited (right panel in Fig. 3c), and the spectrum shows that NMDAR–L-Glu lacks the 5.97 THz absorption peak compared to NMDAR. This may be related to the interaction between γ-COOH in ligand L-Glu and w2. In the normal mode of the NMDAR at a frequency of 10.50 THz, all five water molecules are inactive (left panel in Fig. 3d). When L-Glu acts on the receptor, w1 exhibits an active vibration (right panel in Fig. 3d), mainly a vibration in the opposite direction between the –OH of w1 and the –NH2 of ligand L-Glu. And the absorption peak of NMDAR–L-Glu at 10.50 THz is much higher than that of NMDAR. This is most likely related to the interaction between –NH2 in L-Glu and w1 next to it. It has been suggested that low frequency vibrations and intermolecular interactions may be involved in the biological structure change.60,61 Chou's study has shown that the vibrational signal in IgG antibodies can be transmitted from one location to another via low frequency resonant channels, thus being able to induce relevant conformational changes required for some important biological functions.62 According to González-Jiménez et al., the origin of DNA hydrogen bond breaks in living organisms is believed to be the propagation of low frequency vibrational modes along its length in the form of phonon-like modes that can expand and contract the space between bases.63,64 The use of THz electromagnetic waves with specific frequencies may affect the vibrational behavior of substances and further achieve the effect of regulating biological functions.65,66 In the NMDAR–L-Glu complex system, the state of water molecules can be regarded as that in a certain confined space, different water molecules are present in diverse micro-environments and assume distinct roles. This phenomenon is related to the surrounding molecules, such as charge density, hydrophilicity, steric hindrance and three-dimensional architecture. Owing to environmental variations, the vibrational behavior of one-dimensional confined water differs from that of bulk water in the theoretical calculations conducted by Zhu et al. It is this fact that the absorption of 1.39 THz electromagnetic waves by the confined water can be realized, and then the resonance coupling causes the superpermeation of water, while the response of bulk water to this band is limited.36 Simulations conducted by Okumura et al. demonstrated distinct differences in the absorption spectra between the helical structure and the intermolecular β-sheet structure within amyloid-β amyloid fibrils. Upon resonance with an irradiated laser, numerous stable helical structures are formed subsequent to the destruction of intermolecular β-sheet structures.67 Zhang et al.'s simulation results showed that, under an appropriate THz stimulation, water confined on both sides of an asymmetric wettability membrane channel with different collective vibration modes has dissimilar absorbability to frequency-specific THz stimulation, resulting in a water density gradient that drives the ultrahigh water flux and can effectively penetrate into a bulk water system.68
Fig. 4 and Fig. S5† are the vibrations of water molecules corresponding to specific frequencies. For the same frequency, water molecules at different locations respond differently to THz waves. At a frequency of 5.66 THz, the response mode of w1 is shown in Fig. 4a, where one H atom amplitude is small and the other H atom amplitude is large, and the responses of four water molecules w2–w5 are not obvious. The modes of w2 and w3 at 9.34 THz are shown in Fig. 4b and they are significantly different. The other three water molecules have no obvious response. The modes of w4 and w5 at 11.50 THz are mainly attributed to wagging vibrations, but the directions of vibrations are roughly opposite, as shown by the red and blue arrow heads in Fig. 4c. In addition, an interesting phenomenon is that the same water molecules exhibit distinct response modes at different frequencies, demonstrating a diversity of behaviors. The w5 molecule is inactive at 5.66 THz and 9.34 THz, but active at 11.50 THz. The results show that the response of water molecules to THz waves is spatial and frequency specific. This could be due to the unique environment of the bound water, which leads to distinct responses of the different bound water in the molecule. Therefore, the vibrational behavior of water molecules can serve as valuable indicators for discerning diverse environments and unveiling distinct interactions.
![]() | ||
Fig. 4 Schematic diagram of normal modes of different water molecules at 5.66 THz (a), 9.34 THz (b) and 11.50 THz (c). |
Since the water molecules in the NMDAR–L-Glu complex exhibit significant THz responses, we further studied the influence of the presence or absence of water molecules on the structure of the NMDAR–L-Glu complex. As shown in Fig. 5a, L-Glu is deeply embedded in surrounding molecules, and the green areas represent the locations of interactions between L-Glu and surrounding molecules. In the absence of water molecules (Fig. 5b), the interactions between L-Glu and its surroundings are weakened, as indicated by the arrows, which are largely due to the influence of w1 and w2, while the contributions of other water molecules are relatively small. The interaction change indicated by the red arrow could be related to the change of the w1 normal mode at 10.50 THz (Fig. 3d), while the interaction change indicated by the blue arrow could be related to the change of the w2 normal mode at 5.97 THz (Fig. 3c). The results of further structural optimization are shown in Fig. 5c. The complex shrinks as a whole, and the gaps created by the absence of water molecules are filled by the surrounding molecules. The γ-COOH of L-Glu undergoes torsion, which may be caused by the disappearance of the pulling effect of w2. The interactions between L-Glu and surrounding molecules are adjusted. The calculations show that the binding energy of L-Glu to the surrounding molecules in the initial structure is −27.0 kcal mol−1. In the absence of water molecules, it decreases to −25.0 kcal mol−1. After structural optimization, it increases to −26.2 kcal mol−1, but is still lower than the initial state.
![]() | ||
Fig. 5 The conformation and binding energy of the NMDAR–L-Glu complex in the presence of water molecules (a), in the absence of water molecules (b) and after further structural optimization (c). |
These results suggest that the lack of water molecules leads to a tendency for ligands to detach from or bind less tightly to receptors, indicating that the presence of water molecules can promote the binding of the receptor and the ligand. Water molecules increase the affinity between the protein and the ligand through intermolecular interactions and regulate the structure between the protein and the ligand. This affects the shape and flexibility of protein binding sites, further causing changes in protein configuration and conformation and affects biological activity. Yang et al. studied the influence of water molecules in CDK2 active pockets on the CDK2-ligand binding free energy, and the results showed that the presence of water molecules at the active site helps to increase the value of protein–ligand binding free energy.69 A study by Schiebel et al. showed that water molecules have a key influence on the ligand binding and dissociation process.70 Calculations of Rohani et al. suggest that three specific water molecules near the phylloquinone molecule in the photosystem I protein complex are particularly important and may provide structural factors that increase the strength of hydrogen bonds.59 Arinaminpathy et al. demonstrated a correlation between the migration rate of water within the binding pockets of the glutamate receptor and both the degree of domain closure and the presence of ligands, which has implications for drug design targeting proteins.71,72
Considering that L-Glu is an acidic amino acid, its isoelectric point is less than 7, and it exists as a negative ion in neutral solution, we also investigated two negatively charged L-Glu, and preliminary studies show similar results to neutral L-Glu (see ESI Fig. S5–S11†). Therefore, we conclude that: (1) there are some differences in the THz spectra between the NMDAR and the NMDAR–L-Glu complex, which reflects the ability of the THz spectra to detect non-local conformational changes of proteins. (2) The water molecules in the NMDAR–L-Glu complex have a strong response in the THz band. Water molecules at different locations respond differently to THz waves. Water molecules exhibit distinct response modes at varying frequencies. This shows the specificity of the responses of water molecules to THz waves in a particular confined space. (3) The binding energy of L-Glu to the surrounding molecules decreased after the removal of water, which revealed that water molecules can promote the stability of the structure to a certain extent. Considering water molecules have a small volume, high activity, high flexibility and strong responses to THz waves, it is suggested that water molecules could be used as probes to reveal the structure and interaction of receptor proteins. Numerous studies have demonstrated the crucial and intricate role of water in biological structures; however, the detection of alterations in the molecular structure of these biological systems in solution remains a formidable challenge due to the pronounced absorption of THz waves by water.73,74 González-Jiménez et al., employing femtosecond optical Kerr-effect spectroscopy for the analysis of low frequency vibrational spectra of G-quadruplexes, discovered the presence of strongly underdamped delocalized phonon-like modes that hold potential implications for DNA biology at the atomic level.75 Li et al. present a time-domain terahertz optoacoustics method for the analysis of water-rich samples. Changing the THz photoacoustic signal of water by adjusting the temperature can improve the sensitivity of its analysis. This could provide a powerful tool for exploring molecular interactions and biochemical processes in aqueous solutions.76 The development of these techniques will be helpful to the study of low frequency vibration dynamics of biomolecules under water or physiological conditions.
Footnote |
† Electronic supplementary information (ESI) available: PXRD patterns, structure diagram, FTIR spectrum, and normal modes. See DOI: https://doi.org/10.1039/d4an00697f |
This journal is © The Royal Society of Chemistry 2024 |